To analyze the many-body problem based on differential forms, we apply the Heisenberg equation of motion to the Hamiltonian including the many-body term, focusing on the simplest electron–phonon interaction. Simultaneous differential equations are obtained after applying the Heisenberg equation of motion. To confirm our proposed method, we compare the numerical solutions to the simultaneous differential equations with the analytical results based on the Green's function method. Good agreement is observed between the results, indicating that the analysis based on the Heisenberg equation of motion is suitable for analyzing the steady state. To examine the non-steady state, we investigate weakly pulsed phonon excitation. The numerical calculations suggest that phonons may be conserved if the excitation period is very short. Furthermore, based on the Pauli exclusion principle, the number of simultaneous differential equations that exactly describe electron–electron interactions appears to be finite but becomes extremely large as the number of electrons increases.
To analyze the many-body problem based on differential forms, we apply the Heisenberg equation of motion to the Hamiltonian including the many-body term, focusing on the simplest electron–phonon interaction. Simultaneous differential equations are obtained after applying the Heisenberg equation of motion. To confirm our proposed method, we compare the numerical solutions to the simultaneous differential equations with the analytical results based on the Green's function method. Good agreement is observed between the results, indicating that the analysis based on the Heisenberg equation of motion is suitable for analyzing the steady state. To examine the non-steady state, we investigate weakly pulsed phonon excitation. The numerical calculations suggest that phonons may be conserved if the excitation period is very short. Furthermore, based on the Pauli exclusion principle, the number of simultaneous differential equations that exactly describe electron–electron interactions appears to be finite but becomes extremely large as the number of electrons increases.
For the purpose of computer calculation to evaluate time-dependent quantum properties in finite temperature, we propose new numerical method expressed in the forms of simultaneous differential equations. At first we derive the equation of motion in finite temperature, which is found to be same expression as Heisenberg equation of motion except for the c-number. Based on this equation, we construct numerical method to estimate time-dependent physical properties in finite temperature precisely without using analytical procedures such as Keldysh formalism. Since our approach is so simple and is based on the simultaneous differential equations including no terms related to self-energies, computer programming can be easily performed. It is possible to estimate exact time-dependent physical properties, providing that Hamiltonian of the system is taken to be a one-electron picture. Furthermore, we refer to the application to the many body problem and it is numerically possible to calculate physical properties using Hartree Fock approximation. Our numerical method can be applied to the case even when perturbative Hamiltonians are newly introduced or Hamiltonian shows complex time-dependent behavior. In this article, at first, we derive the equation of motion in finite temperature. Secondly, for the purpose of verification and of exhibiting the usefulness, we show the derivation of gap equation of superconductivity and of sum rule of electrical conductivity and the application to the many body problem. Finally we apply this method to these two cases: the first case is most simplified resonance charge transfer neutralization of an ion and the second is the same process but impurity potential is newly introduced as perturbative Hamiltonian. Through both cases, it is found that neutralization process is not so sensitive to temperature, however, impurity potential as small as 10 meV strongly influences the neutralization of ion.
We lay down the {\em ab initio} many-body quantum theory of electrons and phonons in- and out-of-equilibrium at any temperature. We begin by addressing a fundamental issue concerning the {\em ab initio} Hamiltonian in the harmonic approximation, which we show must be determined {\em self-consistently} to avoid inconsistencies. After identifying the most suitable partitioning into a ``noninteracting'' and an ``interacting'' part we embark on the Green's function diagrammatic analysis. We single out key diagrammatic structures to carry on the expansion in terms of dressed propagators and screened interaction. The final outcome is the finite-temperature nonequilibrium extension of the Hedin equations, featuring the appearance of the time-local Ehrenfest diagram in the electronic self-energy. The Hedin equations have limited practical utility for real-time simulations of driven systems. We leverage the versatility of diagrammatic expansion to generate a closed system of integro-differential equations for the Green's functions and nuclear displacements. These are the Kadanoff-Baym equations for electrons and phonons. Another advantage of the diagrammatic derivation is the ability to use conserving approximations, which ensure the satisfaction of all fundamental conservation laws during the time evolution. As an example we show that the adiabatic Born-Oppenheimer approximation is not conserving whereas its dynamical extension is conserving provided that the electrons are treated in the Fan-Migdal approximation with a dynamically screened electron-phonon coupling. We also derive the formal solution of the Kadanoff-Baym equations in the long time limit and at the steady state. The expansion of the phononic Green's function around the quasi-phonon energies points to a possible correlation-induced splitting of the phonon dispersion in materials with no time-reversal invariance.
Usually, expected values for various physical quantities, such as the number of electrons occupying certain states or the Coulomb interaction between different states of electrons, can be expressed in terms of integrals. In contrast, our method, based on differential forms, shows that expected values can be obtained by averaging over time. To confirm the validity of our method, we prepare the two cases: one is a very simple case with no many-body interaction, and the other is the case where the many-body term is included (the simplest Anderson Hamiltonian). Regarding the simple case without inclusion of many-body term, we prove analytically that the number of electrons occupying any state derived from our method is equivalent to the analytical one evaluated from the Greenʼs function method. When the many-body term is included, our results show good numerical agreement with the analytical ones derived from the Greenʼs function method. By the two cases, the calculation of expected values based on our method is considered valid.
In this work, we lay down the ab initio many-body quantum theory of electrons and phonons in equilibrium as well as in steady-state or time-varying settings.Our focus is on the harmonic approximation, but the developed tools can readily incorporate anharmonic effects.We begin by showing the necessity of determining the ab initio Hamiltonian in a self-consistent manner to ensure the existence of an equilibrium state.We then identify the correct partitioning into a "noninteracting" and an "interacting" part to carry out diagrammatic expansions in terms of dressed propagators and screened interactions.The final outcome is the finite-temperature nonequilibrium extension of the Hedin equations, showcasing the emergence of the coupling between electrons and coherent phonons through the time-local Ehrenfest diagram.The Hedin equations have limited practical utility for real-time simulations of systems driven out of equilibrium by external fields.To overcome this limitation, we leverage the versatility of the diagrammatic approach to generate a closed system of differential equations for the dressed propagators and nuclear displacements.These are the Kadanoff-Baym equations for electrons and phonons.The formalism naturally merges with the theory of conserving approximations, which guarantee the satisfaction of the continuity equation and the conservation of total energy during time evolution.As an example, we show that the popular Born-Oppenheimer approximation is not conserving whereas its dynamical extension is conserving, provided that the electrons are treated in the Fan-Migdal approximation with a dynamically screened electronphonon coupling.We also derive the formal solution of the Kadanoff-Baym equations for nonequilibrium steady states, which is useful for studies in photovoltaics and optoelectronics.Interestingly, the expansion of the phononic Green's function around the quasiphonon energies points to a correlation-induced splitting of the phonon dispersion in materials with no time-reversal invariance.
Abstract Chapter 11 employs variational differential techniques and the Schwinger Action Principle to derive coupled-field Green’s function equations for a multi-component system, modeled as an interacting electron-hole-phonon system. The coupled Fermion Green’s function equations involve five interactions (electron-electron, hole-hole, electron-hole, electron-phonon, and hole-phonon). Starting with quantum Hamilton equations of motion for the various electron/hole creation/annihilation operators and their nonequilibrium average/expectation values, variational differentiation with respect to particle sources leads to a chain of coupled Green’s function equations involving differing species of Green’s functions. For example, the 1-electron Green’s function equation is coupled to the 2-electron Green’s function (as earlier), also to the 1-electron/1-hole Green’s function, and to the Green’s function for 1-electron propagation influenced by a nontrivial phonon field. Similar remarks apply to the 1-hole Green’s function equation, and all others. Higher order Green’s function equations are derived by further variational differentiation with respect to sources, yielding additional couplings. Chapter 11 also introduces the 1-phonon Green’s function, emphasizing the role of electron coupling in phonon propagation, leading to dynamic, nonlocal electron screening of the phonon spectrum and hybridization of the ion and electron plasmons, a Bohm-Staver phonon mode, and the Kohn anomaly. Furthermore, the single-electron Green’s function with only phonon coupling can be rewritten, as usual, coupled to the 2-electron Green’s function with an effective time-dependent electron-electron interaction potential mediated by the 1-phonon Green’s function, leading to the polaron as an electron propagating jointly with its induced lattice polarization. An alternative formulation of the coupled Green’s function equations for the electron-hole-phonon model is applied in the development of a generalized shielded potential approximation, analysing its inverse dielectric screening response function and associated hybridized collective modes. A brief discussion of the (theoretical) origin of the exciton-plasmon interaction follows.
In this thesis, an approximate method introduced to solve time-dependent many-body problems known as time-dependent many-body perturbation theory is studied. Many-body perturbation theory for interacting electrons and phonons is reviewed. In particular, the electron propagator G and an unconventional two-component phonon propagator, which satisfy coupled integral Dyson equations, are introduced. In practice, the associated integral kernels known as the electron Σ and phonon self-energies need to be approximated. The conserving approximations known as the Hartree (-Fock) and the first and second Born approximations, which respect the continuity equation between the electron density and current, are considered in this work. Time-dependent systems of interest are studied in this thesis by using the integro-differential forms of the Dyson equations referred to as the Kadanoff-Baym Equations (KBE). The Kadanoff-Baym equations are introduced for the electron and phonon propagators unconventionally as coupled firstorder integro-differential equations. It is reviewed how these equations are solved numerically by describing an integration rule, a class of practical methods and a parallel implementation of the numerical method. In addition, documentation of how the Kadanoff-Baym equations allow to solve the Bethe-Salpeter Equation (BSE) with the kernel δΣ/δG for the density response function, is provided. In two of the enclosed publications, we benchmarked observables obtained by using the Hartree and partially and fully self-consistent Born approximations against numerically exact results for the two- site, two-electron Holstein model. In this model, the two electrons which are constrained to move between two lattice sites interact indirectly via the electron-phonon coupling. It was observed that only the fully self-consistent Born approximation could cope qualitatively correctly with the competition between the delocalizing and localizing effects of the kinetic and interaction energies. For the other two approximations, spurious bifurcative symmetry breaking with an associated unbounded density response was observed. Nevertheless, also the self-consistent Born approximation was concluded to fail in describing the bipolaronic behavior of the true system. In the third publication, we benchmarked the Hartree-Fock and second Born approximations against an exact method for the few-site Hubbard and Pariser-Parr-Pople models in which the underlying lattice is inert and the electrons interact amongst themselves directly. It was found that the second Born approximation is capable of describing the so-called correlation induced doubly-excited states. This is not possible for time-local approximations such as Hartree-Fock. In addition to the qualitative results, which highlight successes of the applied simple self-energy approximations, the approximate and exact results were also compared on a more quantitative level. It is the quantitative and qualitative results combined which are used in this thesis to assess the quality of the many-body approximations used, with the aim to better understand when these methods are predictive.