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:2008.09709 · Physical review. B./Physical review. B

Time-dependent approach to inelastic scattering spectroscopies in and away from equilibrium: Beyond perturbation theory

Krissia Zawadzki, Luhang Yang, Adrian E. Feiguin

Published Dec 18, 2020• 9 citations
arXivHTMLPDF

Abstract

We propose a nonperturbative numerical approach to calculate the spectrum of a many-body Hamiltonian with time and momentum resolution by exactly recreating a scattering event using the time-dependent Schr\"odinger equation. Akin to an actual inelastic scattering experiment, we explicitly account for the incident and scattered particles (e.g., photons, neutrons, electrons, etc.) in the Hamiltonian and obtain the spectrum by measuring the energy and momentum lost by the particle after interacting with the sample. We illustrate the method by calculating the spin excitations of a Mott-insulating Hubbard chain after a sudden quench with the aid of the time-dependent density matrix renormalization group method. Our formalism can be applied to different forms of spectroscopies, such as neutron and Compton scattering and electron-energy-loss spectroscopy, for instance.

DOI
10.1103/physrevb.102.235141

Related papers

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

Krissia Zawadzki, Alberto Nocera, Adrian E. Feiguin

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

SciPost Physics•Oct 17, 2023•5 citations

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

Adrian E. Feiguin, Alberto Nocera, Krissia Zawadzki

Abstract While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response functon(Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schr\"odinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

•Oct 11, 2022•0 citations

A time-dependent scattering approach to core-level spectroscopies

Krissia Zawadzki, Alberto Nocera, Adrian E. Feiguin

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

arXiv:2002.04142•Feb 11, 2020•2 citations

Nonequilibrium DMFT approach to time-resolved Raman spectroscopy

Philipp Werner, Martin Eckstein, Naoto Tsuji

Raman spectroscopy uses light scattering to extract information on low-energy excitations of solids. The Raman process is described by diagrams which are fourth order in the light-matter interaction and in particular the resonant contribution, which involves four different space-time arguments, is difficult to evaluate. If one instead simulates explicitly the incoming (classical) light pulse, the Raman signal is given by the outgoing photon flux and can be determined from a two-point correlation function. Such a formalism can be used to compute the time-resolved Raman spectrum of nonequilibrium systems, as well as nonlinear signals which are higher order in the incoming field, such as hyper-Raman scattering. Here we explain how to implement this time-dependent formalism within the dynamical mean field theory framework. The method is illustrated with applications to the Holstein-Hubbard model in the strong electron-phonon coupling regime. We demonstrate hyper-Raman scattering in measurements with strong probe fields and frequency mixing signals in the presence of a pump field and simulate the evolution of Stokes and anti-Stokes features after photoexcitations of metallic and Mott insulating systems.

arXiv:2309.02866•Physical review. B./Physical review. B•Dec 27, 2023•2 citations

Developments in many-body theory of quantum transport and spectroscopy with non-equilibrium Green's functions and time-dependent density functional theory

Anna-Maija Uimonen

The problem of quantum dynamics in open systems has gained attention in recent decades and not the least due to the advances made in quantum transport in molecular systems. The main motivation behind quantum transport and molecular electronics is the futuristic goal to be able at some point to replace, or to complement, the silicon-based technology and to make the electronic devices faster. On a fundamental level, one has to deal with time-dependent processes where electron-electron or electron- phonon interactions are of great importance, and they can cause profound quantitative and qualitative changes on the physical and dynamical properties of electronic systems compared to the non-interacting case. Most of the studies of quantum transport have been focused on the steady-state description while neglecting the short-time dynamics. However, the dynamical effects are of great importance since fast- switching processes play a pivotal role in the operation of future devices. We studied the problem of time-dependent electron transport through the Anderson impurity model by using many-body perturbation theory (MBPT) together with Keldysh Green’s functions as well as with time-dependent density functional theory (TDDFT). These methods were compared with numerically exact time-dependent density-matrix renormalization group (tDMRG) method. We found that the many-body perturbation theory results beyond Hartree-Fock approximation were in close agreement with tDMRG results. In addition we studied the possibility of multistablity in the density and current of an interacting nanoscale junction as well as how to reversibly switch between the multiple solutions in time domain. An accurate theoretical treatment of electron correlation even in as simple model as an interacting electron gas at metallic densities still continues to be a challenge; especially description of features in the photoemission spectra due to electron correlations provides a theoretical challenge. The many-body perturbation theory yields a systematic way to study electron-electron (electron-phonon) correlations in various systems. One of the widely used approximations in MBPT is the GW approximation in which the bare interaction line is replaced with screened interaction line in the first order exchange diagram. The GW approximation gives good estimates for the band gap values close to experimental ones but especially the self-consistent GW approximation has a number of deficiencies like washing out of plasmon features and overestimation of bandwidths compared to experiment. One way to improve GW calculations is to include vertex corrections. Unfortunately, the straightforward inclusion of vertex corrections yields negative spectra in some frequency regions. We developed a diagrammatic approach to construct self-energy approximations with positive spectral properties. Our approach consists of expressing a self-energy of response diagram as a product of half-diagrams after which a minimal set of additional diagrams is identified to construct a perfect square. We applied this method to study vertex corrections in a homogeneous electron gas. In addition we analyzed the diagrammatic content of photocurrent with density functional theory. The expression for the photocurrent was obtained as an integral over the Kohn-Sham spectral function renormalized by effective potentials that depend on the exchange correlation kernel of current density functional theory. The expression for the photocurrent gives us the angular dependence of the photocurrent but it does not provide a direct access to the kinetic energy distribution of the photoelectrons.

•Jan 1, 2015•2 citations

Quantum simulation of nuclear inelastic scattering

Weijie Du, James P. Vary, Xingbo Zhao +1 more

We present a time-dependent quantum algorithm for nuclear inelastic scattering in the time-dependent basis function on qubits approach. This algorithm aims to quantum simulate a subset of the nuclear inelastic scattering problems that are of physical interest in which the internal degrees of freedom of the reaction system are excited by time-dependent external interactions. We expect that our algorithm will enable an exponential speedup in simulating the dynamics of the subset of the inelastic scattering problems, which would also be advantageous for the applications to more complicated scattering problems. For a demonstration problem, we solve for the Coulomb excitation of the deuteron, where the quantum simulations are performed with IBM Qiskit.

arXiv:2006.01369•Physical Review A•Jul 23, 2021•36 citations
We propose a nonperturbative numerical approach to calculate the spectrum of a many-body Hamiltonian with time and momentum resolution by exactly recreating a scattering event using the time-dependent Schr"odinger equation. Akin to an actual inelastic scattering experiment, we explicitly account for the incident and scattered particles (e.g., photons, neutrons, electrons, etc.) in the Hamiltonian and obtain the spectrum by measuring the energy and momentum lost by the particle after interacting with the sample. We illustrate the method by calculating the spin excitations of a Mott-insulating Hubbard chain after a sudden quench with the aid of the time-dependent density matrix renormalization group method. Our formalism can be applied to different forms of spectroscopies, such as neutron and Compton scattering and electron-energy-loss spectroscopy, for instance.
While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.
Abstract While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response functon(Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schr\"odinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.
While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.
Raman spectroscopy uses light scattering to extract information on low-energy excitations of solids. The Raman process is described by diagrams which are fourth order in the light-matter interaction and in particular the resonant contribution, which involves four different space-time arguments, is difficult to evaluate. If one instead simulates explicitly the incoming (classical) light pulse, the Raman signal is given by the outgoing photon flux and can be determined from a two-point correlation function. Such a formalism can be used to compute the time-resolved Raman spectrum of nonequilibrium systems, as well as nonlinear signals which are higher order in the incoming field, such as hyper-Raman scattering. Here we explain how to implement this time-dependent formalism within the dynamical mean field theory framework. The method is illustrated with applications to the Holstein-Hubbard model in the strong electron-phonon coupling regime. We demonstrate hyper-Raman scattering in measurements with strong probe fields and frequency mixing signals in the presence of a pump field and simulate the evolution of Stokes and anti-Stokes features after photoexcitations of metallic and Mott insulating systems.
The problem of quantum dynamics in open systems has gained attention in recent decades and not the least due to the advances made in quantum transport in molecular systems. The main motivation behind quantum transport and molecular electronics is the futuristic goal to be able at some point to replace, or to complement, the silicon-based technology and to make the electronic devices faster. On a fundamental level, one has to deal with time-dependent processes where electron-electron or electron- phonon interactions are of great importance, and they can cause profound quantitative and qualitative changes on the physical and dynamical properties of electronic systems compared to the non-interacting case. Most of the studies of quantum transport have been focused on the steady-state description while neglecting the short-time dynamics. However, the dynamical effects are of great importance since fast- switching processes play a pivotal role in the operation of future devices. We studied the problem of time-dependent electron transport through the Anderson impurity model by using many-body perturbation theory (MBPT) together with Keldysh Green’s functions as well as with time-dependent density functional theory (TDDFT). These methods were compared with numerically exact time-dependent density-matrix renormalization group (tDMRG) method. We found that the many-body perturbation theory results beyond Hartree-Fock approximation were in close agreement with tDMRG results. In addition we studied the possibility of multistablity in the density and current of an interacting nanoscale junction as well as how to reversibly switch between the multiple solutions in time domain. An accurate theoretical treatment of electron correlation even in as simple model as an interacting electron gas at metallic densities still continues to be a challenge; especially description of features in the photoemission spectra due to electron correlations provides a theoretical challenge. The many-body perturbation theory yields a systematic way to study electron-electron (electron-phonon) correlations in various systems. One of the widely used approximations in MBPT is the GW approximation in which the bare interaction line is replaced with screened interaction line in the first order exchange diagram. The GW approximation gives good estimates for the band gap values close to experimental ones but especially the self-consistent GW approximation has a number of deficiencies like washing out of plasmon features and overestimation of bandwidths compared to experiment. One way to improve GW calculations is to include vertex corrections. Unfortunately, the straightforward inclusion of vertex corrections yields negative spectra in some frequency regions. We developed a diagrammatic approach to construct self-energy approximations with positive spectral properties. Our approach consists of expressing a self-energy of response diagram as a product of half-diagrams after which a minimal set of additional diagrams is identified to construct a perfect square. We applied this method to study vertex corrections in a homogeneous electron gas. In addition we analyzed the diagrammatic content of photocurrent with density functional theory. The expression for the photocurrent was obtained as an integral over the Kohn-Sham spectral function renormalized by effective potentials that depend on the exchange correlation kernel of current density functional theory. The expression for the photocurrent gives us the angular dependence of the photocurrent but it does not provide a direct access to the kinetic energy distribution of the photoelectrons.
We present a time-dependent quantum algorithm for nuclear inelastic scattering in the time-dependent basis function on qubits approach. This algorithm aims to quantum simulate a subset of the nuclear inelastic scattering problems that are of physical interest in which the internal degrees of freedom of the reaction system are excited by time-dependent external interactions. We expect that our algorithm will enable an exponential speedup in simulating the dynamics of the subset of the inelastic scattering problems, which would also be advantageous for the applications to more complicated scattering problems. For a demonstration problem, we solve for the Coulomb excitation of the deuteron, where the quantum simulations are performed with IBM Qiskit.