Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2209.02254 · Matter and Radiation at Extremes

Physical insights from imaginary-time density–density correlation functions

Tobias Dornheim, Zhandos A. Moldabekov, Panagiotis B. Tolias, Maximilian Böhme, Jan Vorberger

Published Aug 10, 2023• 48 citations
arXivHTMLPDF

Abstract

An accurate theoretical description of the dynamic properties of correlated quantum many-body systems, such as the dynamic structure factor S(q, ω), is important in many fields. Unfortunately, highly accurate quantum Monte Carlo methods are usually restricted to the imaginary time domain, and the analytic continuation of the imaginary-time density–density correlation function F(q, τ) to real frequencies is a notoriously hard problem. Here, it is argued that often no such analytic continuation is required because by definition, F(q, τ) contains the same physical information as does S(q, ω), only represented unfamiliarly. Specifically, it is shown how one can directly extract key information such as the temperature or quasi-particle excitation energies from the τ domain, which is highly relevant for equation-of-state measurements of matter under extreme conditions [T. Dornheim et al., Nat. Commun. 13, 7911 (2022)]. As a practical example, ab initio path-integral Monte Carlo results for the uniform electron gas (UEG) are considered, and it is shown that even nontrivial processes such as the roton feature of the UEG at low density [T. Dornheim et al., Commun. Phys. 5, 304 (2022)] are manifested straightforwardly in F(q, τ). A comprehensive overview is given of various useful properties of F(q, τ) and how it relates to the usual dynamic structure factor. In fact, working directly in the τ domain is advantageous for many reasons and opens up multiple avenues for future applications.

DOI
10.1063/5.0149638

Related papers

Analysing the dynamic structure of warm dense matter in the imaginary-time domain: theoretical models and simulations

Tobias Dornheim, Jan Vorberger, Zhandos A. Moldabekov +1 more

Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [ Nat. Commun. 13 , 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density–density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.

Philosophical Transactions of the Royal Society A Mathematical Physical and Engineering Sciences•Jul 2, 2023•34 citations

Extraction of the frequency moments of spectral densities from imaginary-time correlation function data

Tobias Dornheim, Damar Wicaksono, Juan-Esteban Suarez Cardona +5 more

We introduce an exact framework to compute the positive frequency moments ${M}^{(\ensuremath{\alpha})}(\mathbf{q})=\ensuremath{\langle}{\ensuremath{\omega}}^{\ensuremath{\alpha}}\ensuremath{\rangle}$ of different dynamic properties from imaginary-time quantum Monte Carlo data. As a practical example, we obtain the first five moments of the dynamic structure factor $S(\mathbf{q},\ensuremath{\omega})$ of the uniform electron gas at the electronic Fermi temperature based on ab initio path integral Monte Carlo simulations. We find excellent agreement with known sum rules for $\ensuremath{\alpha}=1,3$, and present results for $\ensuremath{\alpha}=2,4,5$. Our idea can be straightforwardly generalized to other dynamic properties such as the single-particle spectral function $A(\mathbf{q},\ensuremath{\omega})$, and will be useful for a number of applications, including the study of ultracold atoms, exotic warm dense matter, and condensed matter systems.

arXiv:2301.08444•Physical review. B./Physical review. B•Apr 27, 2023•30 citations

Nonlinear density response from imaginary-time correlation functions: Ab initio path integral Monte Carlo simulations of the warm dense electron gas

Tobias Dornheim, Zhandos A. Moldabekov, Jan Vorberger

The ab initio path integral Monte Carlo (PIMC) approach is one of the most successful methods in quantum many-body theory. A particular strength of this method is its straightforward access to imaginary-time correlation functions (ITCFs). For example, the well-known density–density ITCF F(q, τ) allows one to estimate the linear response of a given system for all wave vectors q from a single simulation of the unperturbed system. Moreover, it constitutes the basis for the reconstruction of the dynamic structure factor S(q, ω)—a key quantity in state-of-the-art scattering experiments. In this work, we present analogous relations between the nonlinear density response in the quadratic and cubic order of the perturbation strength and generalized ITCFs measuring correlations between up to four imaginary-time arguments. As a practical demonstration of our new approach, we carry out simulations of the warm dense electron gas and find excellent agreement with previous PIMC results that had been obtained with substantially larger computational effort. In addition, we give a relation between a cubic ITCF and the triple dynamic structure factor S(q1, ω1; q2, ω2), which evokes the enticing possibility to study dynamic three-body effects on an ab initio level.

arXiv:2106.01619•The Journal of Chemical Physics•Aug 4, 2021•55 citations

Spectral and dynamical properties of strongly correlated systems

Giuseppe Carleo

In the first part of the Thesis we mostly concentrate on spectral properties of strongly correlated systems and on their equilibrium properties. This is accomplished by the general concept of imaginary-time dynamics which we apply to a number of different problems in which different strengths of this approach emerge. In Chapter 1 we introduce the formalism that allows for a connection between the quantum and the classical worlds. The connection is established by means of the imaginary-time quantum evolution which, under certain circumstances, is shown to be equivalent to a classical stochastic process. It is further shown that exact static and spectral properties of correlated systems can be obtained when this mapping is feasible. The relationship between the imaginary-time dynamics in different frameworks such as the path-integral and the perturbative one is also underlined. In Chapter 2 we present a specific implementation of the general ideas previously presented. In particular we introduced an extension to lattice systems of the Reptation Monte Carlo algorithm [30] which benefits of a sampling scheme based on directed updates. Specific improvements over the existing methodologies consist in the unbiased evaluation of the imaginary-time path integrals for bosons and a systematic scheme to improve over the Fixed-node approximation for fermions. Applications to the Hubbard and the Heisenberg models are presented. In Chapter 3 we demonstrate the application of the imaginary-time dynamics to the exact study of spectral properties. Subject of our attention is a highly anharmonic and correlated quantum crystal such as Helium 4 at zero temperature.[33] Concerning this system, we have obtained the first ab-initio complete phonon dispersion in good agreement with neutron spectroscopy experiments. Moreover, we have also studied the density excitations of solid helium in a region of wave-vectors in between the collective (phonon) and the single-particle regimes, where the presence of residual coherence in the dynamics shows analogies between the highly anharmonic crystal and the superfluid phase. In Chapter 4 we introduce a novel method, based on the imaginary-time dynamics, to obtain unbiased estimates of fermionic properties.[34] By means of this method and of a very accurate variational state, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, we observe a discontinuous transition to a paramagnetic phase, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition scenario. In Chapter 5 the use of imaginary-time dynamics in the context of finite-temperature response functions is highlighted. As an application, we study an intriguing quantum phase featuring both glassy order and Bose-Einstein condensation. [35] We introduce and validate a model for the role of geometrical frustration in the coexistence of off-diagonal long range order with an amorphous density profile. The exact characterization of the response of the system to an external density perturbation is what allows here to establish the existence of a spin-glass phase. The differences between such a phase and the otherwise insulating "Bose glasses" are further elucidated in the Chapter. In the second part of the Thesis we focus our attention on the dynamics of closed systems out of equilibrium. This is accomplished by both non-stochastic exact methods for the dynamics and the introduction of a novel time-dependent Variational Monte Carlo scheme. In Chapter 6 exact diagonalization schemes and renormalization-based methods for one-dimensional systems are introduced. We identify key phenomenological traits resulting from the many-body correlation in closed systems driven sufficiently away from equilibrium.[31] We provide evidences that the dynamics of interacting lattice bosons away from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. In Chapter 7 we present an extension to dynamical properties of the Variational Quantum Monte Carlo method.[32] This is accomplished by introducing a general class of time-dependent variational states which is based on the mapping of the many-body dynamics onto an instantaneous ground-state problem. The application of the method to the experimentally relevant quantum quenches of interacting bosons reveals the accuracy and the reliability of the introduced numerical scheme. We indeed obtain for the first time a consistent variational description of the approach to the equilibrium of local observables and underline the origin of the metastability and glassy behavior previously identified. In the very last part we draw our conclusions and show some possible paths for stimulating future research.

•Oct 28, 2011•3 citations

Monte Carlo studies of the dynamical response of quantum many-body systems

H.-B. Schüttler, Douglas J. Scalapino

We present a method to obtain dynamical correlation functions of quantum many-body systems in the real-frequency domain from Monte Carlo data of the corresponding imaginary-time Green's functions. The method is based on a least-squares-fit procedure to solve the integral equation relating imaginary-time and real-frequency correlation functions. To demonstrate its feasibility, we have applied our method to imaginary-time Monte Carlo data for the density correlation function of a simple model of interacting spinless fermions in one dimension. We compare our results to the analytical results, available in the limits of zero and strong interaction. We find that despite the presence of noise in the input data, the real-frequency spectra obtained with our method are qualitatively correct. They reflect accurately the particle-hole and the soliton-antisoliton excitations, present in the noninteracting and in the strongly coupled fermion system, respectively. We also compare our method to other approaches that have been proposed for obtaining dynamical correlation functions via Monte Carlo simulation.

Physical review. B, Condensed matter•Oct 1, 1986•39 citations

Quantum time correlation functions from complex time Monte Carlo simulations: A maximum entropy approach

Goran Krilov, Eunji Sim, B. J. Berne

We present a way of combining real-time path integral Monte Carlo simulations with a maximum entropy numerical analytic continuation scheme in a new approach for calculating time correlation functions for finite temperature many body quantum systems. The real-time dynamics is expressed in the form of the symmetrized time correlation function, which is suitable for Monte Carlo methods, and several simulation techniques are presented for evaluating this function accurately up to moderate values of time. The symmetrized time correlation function is then analytically continued in combination with imaginary time data to obtain the real-time correlation function. We test this approach on several exactly solvable problems, including two one-dimensional systems, as well two cases of vibrational relaxation of a system coupled to a dissipative environment. The computed time correlation functions are in good agreement with exact results over several multiples of the thermal time βℏ, and exhibit a significant improvement over analytic continuation of imaginary time correlation functions. Moreover, we show how the method can be systematically improved.

The Journal of Chemical Physics•Jan 15, 2001•77 citations
An accurate theoretical description of the dynamic properties of correlated quantum many-body systems, such as the dynamic structure factor S(q, ω), is important in many fields. Unfortunately, highly accurate quantum Monte Carlo methods are usually restricted to the imaginary time domain, and the analytic continuation of the imaginary-time density–density correlation function F(q, τ) to real frequencies is a notoriously hard problem. Here, it is argued that often no such analytic continuation is required because by definition, F(q, τ) contains the same physical information as does S(q, ω), only represented unfamiliarly. Specifically, it is shown how one can directly extract key information such as the temperature or quasi-particle excitation energies from the τ domain, which is highly relevant for equation-of-state measurements of matter under extreme conditions [T. Dornheim et al., Nat. Commun. 13, 7911 (2022)]. As a practical example, ab initio path-integral Monte Carlo results for the uniform electron gas (UEG) are considered, and it is shown that even nontrivial processes such as the roton feature of the UEG at low density [T. Dornheim et al., Commun. Phys. 5, 304 (2022)] are manifested straightforwardly in F(q, τ). A comprehensive overview is given of various useful properties of F(q, τ) and how it relates to the usual dynamic structure factor. In fact, working directly in the τ domain is advantageous for many reasons and opens up multiple avenues for future applications.
Rigorous diagnostics of experiments with warm dense matter are notoriously difficult. A key method is X-ray Thomson scattering (XRTS), but the interpretation of XRTS measurements is usually based on theoretical models that entail various approximations. Recently, Dornheim et al. [ Nat. Commun. 13 , 7911 (2022)] introduced a new framework for temperature diagnostics of XRTS experiments that is based on imaginary-time correlation functions. On the one hand, switching from the frequency to the imaginary-time domain gives one direct access to a number of physical properties, which facilitates the extraction of the temperature of arbitrarily complex materials without relying on any models or approximations. On the other hand, the bulk of theoretical work in dynamic quantum many-body theory is devoted to the frequency domain, and, to the best of our knowledge, the manifestation of physics properties within the imaginary-time density–density correlation function (ITCF) remains poorly understood. In the present work, we aim to fill this gap by introducing a simple, semi-analytical model for the imaginary-time dependence of two-body correlations within the framework of imaginary-time path integrals. As a practical example, we compare our new model to extensive ab initio path integral Monte Carlo results for the ITCF of a uniform electron gas, and find excellent agreement over a broad range of wavenumbers, densities and temperatures. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.
We introduce an exact framework to compute the positive frequency moments ${M}^{(\ensuremath{\alpha})}(\mathbf{q})=\ensuremath{\langle}{\ensuremath{\omega}}^{\ensuremath{\alpha}}\ensuremath{\rangle}$ of different dynamic properties from imaginary-time quantum Monte Carlo data. As a practical example, we obtain the first five moments of the dynamic structure factor $S(\mathbf{q},\ensuremath{\omega})$ of the uniform electron gas at the electronic Fermi temperature based on ab initio path integral Monte Carlo simulations. We find excellent agreement with known sum rules for $\ensuremath{\alpha}=1,3$, and present results for $\ensuremath{\alpha}=2,4,5$. Our idea can be straightforwardly generalized to other dynamic properties such as the single-particle spectral function $A(\mathbf{q},\ensuremath{\omega})$, and will be useful for a number of applications, including the study of ultracold atoms, exotic warm dense matter, and condensed matter systems.
The ab initio path integral Monte Carlo (PIMC) approach is one of the most successful methods in quantum many-body theory. A particular strength of this method is its straightforward access to imaginary-time correlation functions (ITCFs). For example, the well-known density–density ITCF F(q, τ) allows one to estimate the linear response of a given system for all wave vectors q from a single simulation of the unperturbed system. Moreover, it constitutes the basis for the reconstruction of the dynamic structure factor S(q, ω)—a key quantity in state-of-the-art scattering experiments. In this work, we present analogous relations between the nonlinear density response in the quadratic and cubic order of the perturbation strength and generalized ITCFs measuring correlations between up to four imaginary-time arguments. As a practical demonstration of our new approach, we carry out simulations of the warm dense electron gas and find excellent agreement with previous PIMC results that had been obtained with substantially larger computational effort. In addition, we give a relation between a cubic ITCF and the triple dynamic structure factor S(q1, ω1; q2, ω2), which evokes the enticing possibility to study dynamic three-body effects on an ab initio level.
In the first part of the Thesis we mostly concentrate on spectral properties of strongly correlated systems and on their equilibrium properties. This is accomplished by the general concept of imaginary-time dynamics which we apply to a number of different problems in which different strengths of this approach emerge. In Chapter 1 we introduce the formalism that allows for a connection between the quantum and the classical worlds. The connection is established by means of the imaginary-time quantum evolution which, under certain circumstances, is shown to be equivalent to a classical stochastic process. It is further shown that exact static and spectral properties of correlated systems can be obtained when this mapping is feasible. The relationship between the imaginary-time dynamics in different frameworks such as the path-integral and the perturbative one is also underlined. In Chapter 2 we present a specific implementation of the general ideas previously presented. In particular we introduced an extension to lattice systems of the Reptation Monte Carlo algorithm [30] which benefits of a sampling scheme based on directed updates. Specific improvements over the existing methodologies consist in the unbiased evaluation of the imaginary-time path integrals for bosons and a systematic scheme to improve over the Fixed-node approximation for fermions. Applications to the Hubbard and the Heisenberg models are presented. In Chapter 3 we demonstrate the application of the imaginary-time dynamics to the exact study of spectral properties. Subject of our attention is a highly anharmonic and correlated quantum crystal such as Helium 4 at zero temperature.[33] Concerning this system, we have obtained the first ab-initio complete phonon dispersion in good agreement with neutron spectroscopy experiments. Moreover, we have also studied the density excitations of solid helium in a region of wave-vectors in between the collective (phonon) and the single-particle regimes, where the presence of residual coherence in the dynamics shows analogies between the highly anharmonic crystal and the superfluid phase. In Chapter 4 we introduce a novel method, based on the imaginary-time dynamics, to obtain unbiased estimates of fermionic properties.[34] By means of this method and of a very accurate variational state, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, we observe a discontinuous transition to a paramagnetic phase, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition scenario. In Chapter 5 the use of imaginary-time dynamics in the context of finite-temperature response functions is highlighted. As an application, we study an intriguing quantum phase featuring both glassy order and Bose-Einstein condensation. [35] We introduce and validate a model for the role of geometrical frustration in the coexistence of off-diagonal long range order with an amorphous density profile. The exact characterization of the response of the system to an external density perturbation is what allows here to establish the existence of a spin-glass phase. The differences between such a phase and the otherwise insulating "Bose glasses" are further elucidated in the Chapter. In the second part of the Thesis we focus our attention on the dynamics of closed systems out of equilibrium. This is accomplished by both non-stochastic exact methods for the dynamics and the introduction of a novel time-dependent Variational Monte Carlo scheme. In Chapter 6 exact diagonalization schemes and renormalization-based methods for one-dimensional systems are introduced. We identify key phenomenological traits resulting from the many-body correlation in closed systems driven sufficiently away from equilibrium.[31] We provide evidences that the dynamics of interacting lattice bosons away from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. In Chapter 7 we present an extension to dynamical properties of the Variational Quantum Monte Carlo method.[32] This is accomplished by introducing a general class of time-dependent variational states which is based on the mapping of the many-body dynamics onto an instantaneous ground-state problem. The application of the method to the experimentally relevant quantum quenches of interacting bosons reveals the accuracy and the reliability of the introduced numerical scheme. We indeed obtain for the first time a consistent variational description of the approach to the equilibrium of local observables and underline the origin of the metastability and glassy behavior previously identified. In the very last part we draw our conclusions and show some possible paths for stimulating future research.
We present a method to obtain dynamical correlation functions of quantum many-body systems in the real-frequency domain from Monte Carlo data of the corresponding imaginary-time Green's functions. The method is based on a least-squares-fit procedure to solve the integral equation relating imaginary-time and real-frequency correlation functions. To demonstrate its feasibility, we have applied our method to imaginary-time Monte Carlo data for the density correlation function of a simple model of interacting spinless fermions in one dimension. We compare our results to the analytical results, available in the limits of zero and strong interaction. We find that despite the presence of noise in the input data, the real-frequency spectra obtained with our method are qualitatively correct. They reflect accurately the particle-hole and the soliton-antisoliton excitations, present in the noninteracting and in the strongly coupled fermion system, respectively. We also compare our method to other approaches that have been proposed for obtaining dynamical correlation functions via Monte Carlo simulation.
We present a way of combining real-time path integral Monte Carlo simulations with a maximum entropy numerical analytic continuation scheme in a new approach for calculating time correlation functions for finite temperature many body quantum systems. The real-time dynamics is expressed in the form of the symmetrized time correlation function, which is suitable for Monte Carlo methods, and several simulation techniques are presented for evaluating this function accurately up to moderate values of time. The symmetrized time correlation function is then analytically continued in combination with imaginary time data to obtain the real-time correlation function. We test this approach on several exactly solvable problems, including two one-dimensional systems, as well two cases of vibrational relaxation of a system coupled to a dissipative environment. The computed time correlation functions are in good agreement with exact results over several multiples of the thermal time βℏ, and exhibit a significant improvement over analytic continuation of imaginary time correlation functions. Moreover, we show how the method can be systematically improved.