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:1803.07364

Coupled Kohn-Sham equations for electrons and phonons

C. Y. Wang, Müller, T., Sangeeta Sharma, Eberhard K. U. Gross, J. K. Dewhurst

Published Mar 20, 2018• 0 citations
arXivHTMLPDF

Abstract

This work establishes the algebraic structure of the Kohn-Sham equations to be solved in a density formulation of electron and phonon dynamics, including the superconducting order parameter. A Bogoliubov transform is required to diagonalize both the fermionic and bosonic Kohn-Sham Hamiltonians since they both represent a non-interacting quantum field theory. The Bogoliubov transform for phonons is non-Hermitian in the general case, and the corresponding time-evolution is non-unitary. Several sufficient conditions for ensuring that the bosonic eigenvalues are real are provided and a practical method for solving the system is described. Finally, we produce a set of approximate mean-field potentials which are functionals of the electronic and phononic density matrices and depend on the electron-phonon vertex.

DOI
10.48550/arxiv.1803.07364

Related papers

Inverse Kohn–Sham Equations Derived from the Density Equation Theory

T. Kato, Katsuyuki Nobusada, Shinji Saito

We formulate the inverse Kohn–Sham equations based on the density equation theory. The present formulation relies on neither action functionals nor force balance equation, but is constructed by comparing the first-order density equation for the electron interacting real system and that for the electron non-interacting Kohn–Sham system. The derived inverse Kohn–Sham integral equation requires information of an exact ground-state electron density to calculate an exact Kohn–Sham potential for a ground state and it requires an exact time-dependent electron density and a time-derivative of the exact time-dependent first-order density matrix to calculate an exact time-dependent Kohn–Sham potential. A possible algorithm to construct the Kohn–Sham potential for a ground state is suggested as well as a solution algorithm for the time-dependent inverse Kohn–Sham problem without resorting to any approximations. These two solution algorithms, in combination, suggest a way to perform an exact time-dependent density functional calculation as a predictive method.

Journal of the Physical Society of Japan•Jan 27, 2020•4 citations

Exact factorization-based density functional theory of electron-phonon systems

Ryan Requist, César R. Proetto, Eberhard K. U. Gross

Density functional theory is generalized to incorporate electron-phonon coupling. A Kohn-Sham equation yielding the electronic density ${n}_{U}(\mathbf{r})$, a conditional probability density depending parametrically on the phonon normal mode amplitudes $U={{U}_{\mathbf{q}\ensuremath{\lambda}}}$, is coupled to the nuclear Schr\"odinger equation of the exact factorization method. The phonon modes are defined from the harmonic expansion of the nuclear Schr\"odinger equation. A nonzero Berry curvature on nuclear configuration space affects the phonon modes, showing that the potential energy surface alone is generally not sufficient to define the phonons. An orbital-dependent functional approximation for the nonadiabatic exchange-correlation energy reproduces the leading-order nonadiabatic electron-phonon-induced band structure renormalization in the Fr\"ohlich model.

arXiv:1901.07523•Physical review. B./Physical review. B•Apr 24, 2019•25 citations

A Rigorous Extension of the Kohn-Sham Equation for Strongly Correlated Electron Systems

Koichi Kusakabe

By introducing a set of auxiliary equations representing a many-body system, we have derived an extension of the Kohn-Sham scheme for the density functional theory. These equations consist of a Kohn-Sham-type equation determining single-particle orbitals and an eigen-value equation for an effective many-body problem. A variational method similar to the Kohn-Sham technique was utilized to derive effective interactions as well as effective potentials without artificial substitution of a Hubbard-type interaction and a mean-field correction in the energy functional. The second equation is described by an effective many-body Hamiltonian with both 2-body interactions and mean-field terms. Rigorous formulation of the extended Kohn-Sham equation is also given in accordance with the Hadjisavvas-Theophilou formulation. Our formulation can be interpreted as a way to define models of the strongly correlated electron systems, e.g. the Hubbard model

arXiv:cond-mat/0103094•Journal of the Physical Society of Japan•Jul 15, 2001•23 citations

Time-dependent potential through an Ansatz for the Kohn-Sham orbitals

Rudolph J Magyar

Given the time-evolution of an electron charge density, the local potential in Kohn-Sham time-dependent density functional theory (KS-TDDFT) can be modeled as a sum of instantaneous and dynamic contributions by assuming a certain form of the time-dependent Kohn-Sham (TD-KS) orbitals. The instantaneous part is obtained numerically using methods from ground-state density functional theory (DFT) and the dynamic part is expressed in terms of a velocity potential that depends on the electron current density. The suggested form of the TD-KS orbitals satisfies several known constraints (orthonormality, N-representability, J-representability), and the domain of validity is shown to depend on the evolution of the instantaneous quantities. Through this decomposition, we can relate time-dependent and ground -state V-representability. The resulting potentials are shown numerically to approximate the exact time-dependent Kohn-Sham potentials for a set of 3 non-singlet two-particle systems (a Kohn-mode, a Coulomb explosion, and a double quantum well) where the exact solutions and reference densities are known or obtained through configuration interaction (CI) approaches.

arXiv:1611.02636•Nov 8, 2016•0 citations

On the Kohn-Sham Approach to Time-Dependent Problems in a Density-Functional Framework

Walter Tarantino

Predictivity of the Kohn-Sham approach to dynamical problems, when regarded as an initial value problem in a time-dependent density functional framework, is analysed for a class of models for which the argument devised in the work of Maitra et al. (Phys. Rev. A 78, 056501 (2008), arXiv:cond-mat/0710.0018) for the standard electronic many-body problem does not apply. The original argument is here extended and revised. As a result, predictivity for this class of problems seems possible only at the price of introducing extra unknown functionals in the corresponding Kohn-Sham equation. Furthermore, the same argument, when applied to original electronic problem, suggests that the Hartree-exchange-correlation potential is not unambiguously identified by the contemporary and past densities and initial states, but also requires knowledge of the divergence of the contemporary Kohn-Sham current.

arXiv:1506.04602•Jun 15, 2015•0 citations

Generalization of the Kohn-Sham system enabling representing arbitary one electron density matrices

van Dam, J J Hubertus

Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of non-interacting particles, is the work horse of the theory. The particular form of the Kohn-Sham wavefunction admits only idem-potent one electron density matrices whereas wavefunctions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept, and introducing a suitable dot-product as well as a probability density a non-interacting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans theorem are available. Finally, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single determinant wavefunctions.

arXiv:1512.06120•Dec 7, 2015•0 citations
This work establishes the algebraic structure of the Kohn-Sham equations to be solved in a density formulation of electron and phonon dynamics, including the superconducting order parameter. A Bogoliubov transform is required to diagonalize both the fermionic and bosonic Kohn-Sham Hamiltonians since they both represent a non-interacting quantum field theory. The Bogoliubov transform for phonons is non-Hermitian in the general case, and the corresponding time-evolution is non-unitary. Several sufficient conditions for ensuring that the bosonic eigenvalues are real are provided and a practical method for solving the system is described. Finally, we produce a set of approximate mean-field potentials which are functionals of the electronic and phononic density matrices and depend on the electron-phonon vertex.
We formulate the inverse Kohn–Sham equations based on the density equation theory. The present formulation relies on neither action functionals nor force balance equation, but is constructed by comparing the first-order density equation for the electron interacting real system and that for the electron non-interacting Kohn–Sham system. The derived inverse Kohn–Sham integral equation requires information of an exact ground-state electron density to calculate an exact Kohn–Sham potential for a ground state and it requires an exact time-dependent electron density and a time-derivative of the exact time-dependent first-order density matrix to calculate an exact time-dependent Kohn–Sham potential. A possible algorithm to construct the Kohn–Sham potential for a ground state is suggested as well as a solution algorithm for the time-dependent inverse Kohn–Sham problem without resorting to any approximations. These two solution algorithms, in combination, suggest a way to perform an exact time-dependent density functional calculation as a predictive method.
Density functional theory is generalized to incorporate electron-phonon coupling. A Kohn-Sham equation yielding the electronic density nU(r){n}_{U}(\mathbf{r})nU​(r), a conditional probability density depending parametrically on the phonon normal mode amplitudes $U={{U}_{\mathbf{q}\ensuremath{\lambda}}}$, is coupled to the nuclear Schr\"odinger equation of the exact factorization method. The phonon modes are defined from the harmonic expansion of the nuclear Schr\"odinger equation. A nonzero Berry curvature on nuclear configuration space affects the phonon modes, showing that the potential energy surface alone is generally not sufficient to define the phonons. An orbital-dependent functional approximation for the nonadiabatic exchange-correlation energy reproduces the leading-order nonadiabatic electron-phonon-induced band structure renormalization in the Fr\"ohlich model.
By introducing a set of auxiliary equations representing a many-body system, we have derived an extension of the Kohn-Sham scheme for the density functional theory. These equations consist of a Kohn-Sham-type equation determining single-particle orbitals and an eigen-value equation for an effective many-body problem. A variational method similar to the Kohn-Sham technique was utilized to derive effective interactions as well as effective potentials without artificial substitution of a Hubbard-type interaction and a mean-field correction in the energy functional. The second equation is described by an effective many-body Hamiltonian with both 2-body interactions and mean-field terms. Rigorous formulation of the extended Kohn-Sham equation is also given in accordance with the Hadjisavvas-Theophilou formulation. Our formulation can be interpreted as a way to define models of the strongly correlated electron systems, e.g. the Hubbard model
Given the time-evolution of an electron charge density, the local potential in Kohn-Sham time-dependent density functional theory (KS-TDDFT) can be modeled as a sum of instantaneous and dynamic contributions by assuming a certain form of the time-dependent Kohn-Sham (TD-KS) orbitals. The instantaneous part is obtained numerically using methods from ground-state density functional theory (DFT) and the dynamic part is expressed in terms of a velocity potential that depends on the electron current density. The suggested form of the TD-KS orbitals satisfies several known constraints (orthonormality, N-representability, J-representability), and the domain of validity is shown to depend on the evolution of the instantaneous quantities. Through this decomposition, we can relate time-dependent and ground -state V-representability. The resulting potentials are shown numerically to approximate the exact time-dependent Kohn-Sham potentials for a set of 3 non-singlet two-particle systems (a Kohn-mode, a Coulomb explosion, and a double quantum well) where the exact solutions and reference densities are known or obtained through configuration interaction (CI) approaches.
Predictivity of the Kohn-Sham approach to dynamical problems, when regarded as an initial value problem in a time-dependent density functional framework, is analysed for a class of models for which the argument devised in the work of Maitra et al. (Phys. Rev. A 78, 056501 (2008), arXiv:cond-mat/0710.0018) for the standard electronic many-body problem does not apply. The original argument is here extended and revised. As a result, predictivity for this class of problems seems possible only at the price of introducing extra unknown functionals in the corresponding Kohn-Sham equation. Furthermore, the same argument, when applied to original electronic problem, suggests that the Hartree-exchange-correlation potential is not unambiguously identified by the contemporary and past densities and initial states, but also requires knowledge of the divergence of the contemporary Kohn-Sham current.
Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of non-interacting particles, is the work horse of the theory. The particular form of the Kohn-Sham wavefunction admits only idem-potent one electron density matrices whereas wavefunctions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept, and introducing a suitable dot-product as well as a probability density a non-interacting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans theorem are available. Finally, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single determinant wavefunctions.