[논문 리뷰] Tangent equations of motion for nonlinear response functions

Atsushi Ono|arXiv (Cornell University)|2026. 03. 23.

Spectroscopy and Quantum Chemical Studies인용 수 0

한 줄 요약

본 논문은 Gateaux 도함수에 기반한 접선 운동방정식 TEOM을 개발하여 실시간 역학으로부터 비선 응답 함수를 직접 계산하고, 명시적 다점 상관관계 항과 불안정한 유한차분 방법을 피한다.

ABSTRACT

Nonlinear response functions, formulated as multipoint correlation functions or Volterra kernels, encode the dynamical and spectroscopic properties of physical systems and underpin a wide range of nonlinear transport and optical phenomena. However, their evaluation rapidly becomes prohibitive at high orders because of combinatorial (often factorial) scaling or severe numerical errors. Here, we establish a systematic and efficient framework to compute nonlinear response functions directly from real-time dynamics, without explicitly constructing multipoint correlators or relying on numerically unstable finite-difference methods for order-resolved extraction. Our approach is based on the Gateaux derivative with respect to the external field in function space, which yields a closed hierarchy of tangent equations of motion (TEOM). Propagating the TEOM alongside the original dynamics isolates each perturbative order with high accuracy, providing a term-by-term decomposition of physical contributions. The computational cost scales exponentially with response order in the fully general setting and reduces to polynomial complexity when all perturbation directions are identical; both regimes avoid the factorial scaling of explicit multipoint-correlator evaluations. We demonstrate the power of TEOM by computing frequency-resolved fifth-order response functions for a solid-state electron model and by obtaining nonlinear response functions up to the 49th order with controlled accuracy in a classical Duffing oscillator. We further show that our time-evolution formulation allows optical conductivities to be evaluated directly while remaining numerically stable even near zero frequency. TEOM can be incorporated seamlessly into existing real-time evolution methods, yielding a general framework for computing nonlinear response functions in quantum and classical dynamical systems.

연구 동기 및 목표

• 양자 및 고전 시스템에서 비선형 응답 함수의 중요성과 도전과제를 동기화한다.
• 실시간 역학으로부터 차수에 따라 분해된 비선 응답을 얻기 위해 Gateaux 도함수 기반의 EOM 프레임워크를 도입한다.
• 필드-함수 파생을 전개하는 접선 운동방정식(TEOM) 계층 구조를Closed하게 개발하여 상태 도함수를 전파한다.
• 다변량 비선형 응답 커널을 재구성하기 위한 주파수 해석 프로토콜을 제시한다.
• 양자 및 고전 모델에서 방법을 시연하고 계산 스케일링 및 안정성을 분석한다.

제안 방법

• 비선형 응답 함수를 n차 지연 커널과 그 주파수 영역 표현으로 정의한다.
• 외부장에 대한 Gateaux 도함수를 도입하여 방향성 함수 도함수를 얻는다.
• 상태 도함수와 원래 역학과 함께 폐쇄된 TEOM 계층 구조를 파생하고 전파한다.
• 특정 요 perturbed field(코사인 및 사인 성분)로 infinitesimal variation으로부터 다변량 커널을 재구성한다.
• infinitesimal Gateaux 도함수와 특정 섭 perturbations를 활용하는 재구성 절차를 통해 주파수 해석 커널을 추출한다.
• 확산적 구현 측면에서의 규모화, 패리티 고려, 자동 미분과의 연결성 등을 논의한다.

실험 결과

연구 질문

• RQ1 비점 상관관계를 명시적으로 구성하지 않고도 실시간 역학으로부터 비선형 응답 함수를 직접 계산할 수 있는가?
• RQ2 동적 생성기가 외부장이나 evolving 상태에 의존하는 경우에도 유효한 폐쇄 TEOM 계층 구조를 형성할 수 있는가?
• RQ3 Gateaux 도함수로부터 주파수 해석된 비대각선 비선형 커널은 어떻게 재구성될 수 있는가?
• RQ4 양자 및 고전 시스템에서 고차 응답에 대한 TEOM의 계산 스케일링 특성은 무엇인가?
• RQ5 이 프레임워크가 고체상 전자 모델 및 Duffing 진동기 같은 벤치마크 모델에서 어떤 성능을 보이는가?

주요 결과

• TEOM은 유한 차분 차감 없이 실시간 역학으로부터 비선 응답 함수의 안정적이고 차수에 따른 재구성을 제공합니다.
• TEOM 계층은 원래의 역학과 함께 혼합 방향 도함수를 전파하며 차수에 대해 지수적(팩토리얼이 아님) 스케일링이 가능하여 고차에 대한 처리 가능성이 향상됩니다.
• 다섯째 차 주파수 해석 커널이 네-대역 전자 모델에서 계산되었고, Duffing 진동기에서는 매우 높은 차수(최대 49차)까지 제어된 정확성으로 달성됩니다.
• 이 방법은 단위, Liouville 공간 및 자기 일관 평균장 등 동역학의 선형 및 비선형 제너레이터를 지원합니다.
• 시간 진화 형태식은 0 주파수 근처에서의 수치적 안정성과 함께 광학 전도도를 직접 평가할 수 있도록 허용합니다.