[논문 리뷰] Thermal response of a Fermi-Pasta-Ulam chain with Andersen thermostats

Federico D’Ambrosio, Marco Baiesi|arXiv (Cornell University)|2017. 06. 21.

Advanced Thermodynamics and Statistical Mechanics참고 문헌 46인용 수 2

한 줄 요약

이 논문은 앤더슨 열조절기를 사용하여 관성 있는 비평형 시스템에 대한 변동-응답 관계를 개발하여, 작은 외부 편미분 없이도 열 감도(예: 열용량 및 열팽창계수)를 정확하게 계산할 수 있도록 한다. 이 방법은 강한 비평형 영역에서 직접 편미분 시뮬레이션보다 우수하며, 경계 열은행이 있는 페르미-파스토-울람(Fermi-Pasta-Ulam) 사슬에서 스무스하고 낮은 오차를 가진 감도 추정치를 제공한다.

ABSTRACT

The linear response to temperature variations is well characterised for equilibrium systems but a similar theory is not available, for example, for inertial heat conducting systems, whose paradigm is the Fermi-Pasta-Ulam (FPU) model driven by two different boundary temperatures. For models of inertial systems out of equilibrium, including relaxing systems, we show that Andersen thermostats are a natural tool for studying the thermal response. We derive a fluctuation-response relation that allows to predict thermal expansion coefficients or the heat capacitance in nonequilibrium regimes. Simulations of the FPU chain of oscillators suggest that estimates of susceptibilities obtained with our relation are better than those obtained via a small perturbation.

연구 동기 및 목표

• 표준 평형 변동-응답 관계가 실패하는 관성 있는 비평형 시스템에서 열 감도에 대한 선형 응답 이론을 개발하는 것.
• 예를 들어 구동된 페르미-파스토-울람(Fermi-Pasta-Ulam, FPU) 사슬처럼 열류가 존재하는 시스템에서 열용량 및 열팽창계수와 같은 열응답 계수를 신뢰할 수 있는 프레임워크 없이 계산할 수 없는 문제를 다루는 것.
• 앤더슨 열조절기가 라운지엔 열은행에 비해 수학적으로 다룰 수 있고 물리적으로 의미 있는 대안을 제공함을 보여주는 것.
• 수치적으로 제안된 변동-응답 관계가 직접적인 작은 편미분보다 더 정확하고 안정적인 감도 추정치를 제공함을 보여주는 것.

제안 방법

• 스토케스틱 궤적의 경로가중비를 바탕으로, 앤더슨 열조절기를 통해 구동되는 관성 시스템에서 열 감도에 대한 변동-응답 관계를 유도한다.
• 앤더슨 열조절기에서 속도 재설정의 스토케스틱 성격을 이용하여, 라운지엔 기반 접근법에서 흔히 발생하는 특이성을 피하는 표준 응답 함수 기법을 적용한다.
• 감도를 관측량과 엔트로피 생성률의 기능으로 표현하며, 핵심 결과는 대칭 및 비대칭 성분을 포함하는 경로가중비의 표현인 식 (12)–(14)에 의해 주어진다.
• 등온 및 비평형 정 steady 상태를 모두 사용하여 FPU 사슬 길이의 감도를 계산하기 위해 이 이론을 적용한다.
• 서로 다른 온도에서 열조절된 경계를 가진 FPU 사슬의 수치 시뮬레이션을 수행하여 유도된 공식을 통해 감도를 계산하고, 직접 편미분 결과와 비교한다.
• 정적 감도 χs_L를 t → ∞ 일 때 1/t에 대한 선형 외삽을 통해 추정하여 수렴성을 향상시키고 통계적 노이즈를 줄인다.

실험 결과

연구 질문

• RQ1 경계 열은행이 있는 FPU 사슬과 같은 관성 있는 비평형 시스템에서 열 감도에 대한 변동-응답 관계를 도출할 수 있는가?
• RQ2 유도된 변동-응답 공식은 열 감도 추정에서 직접적인 작은 편미분 방법보다 성능가능성이 어떻게 뛰어나게 되는가?
• RQ3 등온 상태와 비평형 정 steady 상태에서 시스템 길이의 감도가 경계 온도 변화에 대해 다를 수 있으며, 이는 비평형 효과를 어떻게 드러내는가?
• RQ4 비평형 상태에서 발산할 수 있는 경로가중비의 대칭 및 비대칭 성분(C+ 및 C−)이 여전히 감도를 정확하게 재구성하는 데 사용될 수 있는가?
• RQ5 정적 감도 χs_L는 비평형 정 steady 상태에서 평균 온도와 열류 방향에 따라 달라지는가?

주요 결과

• 유도된 변동-응답 관계(식 12–14)는 작은 외부 편미분에 의존하지 않고도 관성 있는 비평형 시스템에서 열 감도를 강력하고 정확하게 계산할 수 있는 방법을 제공한다.
• 수치 결과는 변동-응답 공식을 통해 추정된 감도가 직접 편미분 시뮬레이션으로부터 얻은 결과보다 더 스무스하고 노이즈가 적음을 보여준다.
• 등온 상태에서는 예상대로 대칭 및 비대칭 응답 성분(C+ 및 C−)이 동일하며, 감도는 표준 쿠보 공식 χ = 2C−와 일치한다.
• 비평형 정 steady 상태에서는 C+ 및 C−가 크게 발산하여 강한 비평형 효과를 나타내지만, 전체 공식은 여전히 정확한 감도 추정치를 제공한다.
• 정적 감도 χs_L는 온도가 증가함에 따라 증가한다: 등온 상태에서 T = 0.1일 때 χs_L = 0.524(3), T = 0.3일 때 0.601(2), T = 0.5일 때 0.625(2).
• 비평형 상태에서는 열류 방향에 따라 χs_L가 달라진다: T1 = 0.1, TN = 0.5일 때 χs_L = 0.625(2), 반대로 T1 = 0.5, TN = 0.1일 때는 0.553(3)로 나타나, 평형 평균에서는 포착되지 않는 진정한 비평형 효과를 보여준다.