Generalized Path Integral Energy and Heat Capacity Estimators of Quantum Oscillators and Crystals using Harmonic Mapping

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Moustafa, Sabry G., Schultz, Andrew J.
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866914918001803264
author Moustafa, Sabry G.
Schultz, Andrew J.
author_facet Moustafa, Sabry G.
Schultz, Andrew J.
contents Imaginary-time path integral (PI) is a rigorous tool to treat nuclear quantum effects in static properties. However, with its high computational demand, it is crucial to devise precise estimators. We introduce generalized PI estimators for the energy and heat capacity that utilize coordinate mapping. While it can reduce to the standard thermodynamic and centroid virial (CVir) estimators, the formulation can also take advantage of harmonic character of quantum oscillators and crystals to construct a coordinate mapping. This yields harmonically mapped averaging (HMA) estimators, with mappings that decouple (HMAc) or couple (HMAq) the centroid and internal modes. The HMAq is constructed with normal mode coordinates (HMAq-NM) with quadratic scaling of cost or harmonic oscillator staging (HMAq-SG) coordinates with linear scaling. The estimator performance is examined for a 1D anharmonic oscillator and a 3D Lennard-Jones crystal using path integral molecular dynamics (PIMD) simulation. The HMA estimators consistently provide more precise estimates compared to CVir, with the best performance obtained by HMAq-NM, followed by HMAq-SG, and then HMAc. We also examine the effect of anharmonicity (for AO), intrinsic quantumness, and Trotter number. The HMA formulation introduced assumes the availability of forces and Hessian matrix; however, an equally efficient finite difference alternative is possible when these derivatives are inaccessible. The remarkable improvement in precision offered by HMAq estimators provides a framework for efficient PI simulation of more challenging systems, such as those based on \textit{ab initio} calculations.
format Preprint
id arxiv_https___arxiv_org_abs_2408_11214
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Generalized Path Integral Energy and Heat Capacity Estimators of Quantum Oscillators and Crystals using Harmonic Mapping
Moustafa, Sabry G.
Schultz, Andrew J.
Statistical Mechanics
Imaginary-time path integral (PI) is a rigorous tool to treat nuclear quantum effects in static properties. However, with its high computational demand, it is crucial to devise precise estimators. We introduce generalized PI estimators for the energy and heat capacity that utilize coordinate mapping. While it can reduce to the standard thermodynamic and centroid virial (CVir) estimators, the formulation can also take advantage of harmonic character of quantum oscillators and crystals to construct a coordinate mapping. This yields harmonically mapped averaging (HMA) estimators, with mappings that decouple (HMAc) or couple (HMAq) the centroid and internal modes. The HMAq is constructed with normal mode coordinates (HMAq-NM) with quadratic scaling of cost or harmonic oscillator staging (HMAq-SG) coordinates with linear scaling. The estimator performance is examined for a 1D anharmonic oscillator and a 3D Lennard-Jones crystal using path integral molecular dynamics (PIMD) simulation. The HMA estimators consistently provide more precise estimates compared to CVir, with the best performance obtained by HMAq-NM, followed by HMAq-SG, and then HMAc. We also examine the effect of anharmonicity (for AO), intrinsic quantumness, and Trotter number. The HMA formulation introduced assumes the availability of forces and Hessian matrix; however, an equally efficient finite difference alternative is possible when these derivatives are inaccessible. The remarkable improvement in precision offered by HMAq estimators provides a framework for efficient PI simulation of more challenging systems, such as those based on \textit{ab initio} calculations.
title Generalized Path Integral Energy and Heat Capacity Estimators of Quantum Oscillators and Crystals using Harmonic Mapping
topic Statistical Mechanics
url https://arxiv.org/abs/2408.11214