Abstract:The quantum harmonic oscillator is the fundamental building block to compute thermal properties of virtually any dielectric crystal at low temperatures in terms of phonons, extended further to cases with anharmonic couplings, or even disordered solids. In general, Path Integral Monte Carlo (PIMC) or Molecular Dynamics (PIMD) methods are powerful tools to determine stochastically thermodynamic quantities without systematic bias, not relying on perturbative schemes. Addressing transport properties, for instance calculating thermal conductivity from PIMC, however, is substantially more difficult. Although correlation functions of current operators can be determined by PIMC from analytic continuation on the imaginary-time axis, Bayesian methods are usually employed for the numerical inversion back to real-time response functions. This task not only strongly relies on the accuracy of the PIMC data, but also introduces noticeable dependence on the model used for the inversion. Here, we address both difficulties with care. In particular, we first devise improved estimators for current correlations which substantially reduce the variance of the PIMC data. Next, we provide a neat statistical approach to the inversion problem, blending into a fresh workflow the classical stochastic maximum entropy method together with recent notions borrowed from statistical learning theory. We test our ideas on a single harmonic oscillator and a collection of oscillators with a continuous distribution of frequencies, and provide indications of the performance of our method in the case of a particle in a double well potential. This work establishes solid grounds for an unbiased, fully quantum mechanical calculation of transport properties in solids.

Using The Thermal Gaussian Approximation For The Boltzmann Operator In Semiclassical Initial Value Time Correlation Functions

Linearized Semiclassical Initial Value Time Correlation Functions Using the Thermal Gaussian Approximation: Applications to Condensed Phase Systems.

Recent Advances in the Linearized Semiclassical Initial Value Representation/classical Wigner Model for the Thermal Correlation Function

Real Time Correlation Function in a Single Phase Space Integral--beyond the Linearized Semiclassical Initial Value Representation

Finite-Temperature Correlation Functions Obtained from Combined Real- and Imaginary-Time Propagation of Variational Thawed Gaussian Wavepackets

A Semiclassical Initial-Value Representation for Quantum Propagator and Boltzmann Operator

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

A Simple Model For The Treatment Of Imaginary Frequencies In Chemical Reaction Rates And Molecular Liquids

Semiclassical approximation of the Wigner function for the canonical ensemble

Truncated Gaussian basis approach for simulating many-body dynamics

Generalized Discrete Truncated Wigner Approximation for Nonadiabtic Quantum-Classical Dynamics

Test of the Consistency of Various Linearized Semiclassical Initial Value Time Correlation Functions in Application to Inelastic Neutron Scattering from Liquid Para-Hydrogen

Semiclassical Reaction Rate Constant of Non-Adiabatic System by Means of Initial Value Representation

Frozen Gaussian Series Representation of the Imaginary Time Propagator Theory and Numerical Tests

Time-dependent ghost-Gutzwiller non-equilibrium dynamics

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

A Semiclassical Study of the Thermal Conductivity of Low Temperature Liquids

Thermal quantum time-correlation functions from classical-like dynamics

Thermal Broadening of Phonon Spectral Function in Classical Lattice Models: Projective Truncation Approximation

Semiclassical study of diagonal and offdiagonal functions in the eigenstate thermalization hypothesis

Time correlation functions for quantum systems: validating Bayesian approaches for harmonic oscillators and beyond