Adiabatic approximation in time-dependent reduced-density-matrix functional theory

With the aim of describing real-time electron dynamics, we introduce an adiabatic approximation for the equation of motion of the one-body reduced-density matrix (one-matrix). The eigenvalues of the one-matrix, which represent the occupation numbers of single-particle orbitals, are obtained from the constrained minimization of the instantaneous ground state energy functional rather than from their dynamical equations. To clarify the motivation for this minimization condition, we discuss a sequence of adiabatic energy functionals, each obeying a minimum principle. The performance of the approximation vis-a`-vis nonadiabatic effects is assessed in real-time simulations for a two-site Hubbard model. Due to the presence of Landau-Zener-type transitions, the system evolves into a nonstationary state with persistent oscillations in the observables. The amplitude and phase of the oscillations exhibit resonance behavior both with respect to the strength of the electron-electron interaction and the rate of variation of the external potential. Both types of resonances have the same origin -- the interference of dynamical and scattering phases.Comment: 18 pages, 8 figures; thoroughly revise

Time-dependent occupation numbers in reduced-density-matrix functional theory: Application to an interacting Landau-Zener model

We prove that if the two-body terms in the equation of motion for the one-body reduced density matrix are approximated by ground-state functionals, the eigenvalues of the one-body reduced density matrix (occupation numbers) remain constant in time. This deficiency is related to the inability of such an approximation to account for relative phases in the two-body reduced density matrix. We derive an exact differential equation giving the functional dependence of these phases in an interacting Landau-Zener model and study their behavior in short- and long-time regimes. The phases undergo resonances whenever the occupation numbers approach the boundaries of the interval [0,1]. In the long-time regime, the occupation numbers display correlation-induced oscillations and the memory dependence of the functionals assumes a simple form.Comment: 6 pages, revised, Fig. 2 adde

Approximate formula for the macroscopic polarization including quantum fluctuations

The many-body Berry phase formula for the macroscopic polarization is approximated by a sum of natural orbital geometric phases with fractional occupation numbers accounting for the dominant correlation effects. This reduced formula accurately reproduces the exact polarization in the Rice-Mele-Hubbard model across the band insulator-Mott insulator transition. A similar formula based on a one-body reduced Berry curvature accurately predicts the interaction-induced quenching of Thouless topological charge pumping

Model Hamiltonian for strongly-correlated systems: Systematic, self-consistent, and unique construction

An interacting lattice model describing the subspace spanned by a set of strongly-correlated bands is rigorously coupled to density functional theory to enable ab initio calculations of geometric and topological material properties. The strongly-correlated subspace is identified from the occupation number band structure as opposed to a mean-field energy band structure. The self-consistent solution of the many-body model Hamiltonian and a generalized Kohn-Sham equation exactly incorporates momentum-dependent and crystal-symmetric correlations into electronic structure calculations in a way that does not rely on a separation of energy scales. Calculations for a multiorbital Hubbard model demonstrate that the theory accurately reproduces the many-body polarization.Comment: 19 pages, 11 figure

Quantum Electrical Dipole in Triangular Systems: a Model for Spontaneous Polarity in Metal Clusters

Triangular symmetric molecules with mirror symmetry perpendicular to the 3-fold axis are forbidden to have a fixed electrical dipole moment. However, if the ground state is orbitally degenerate and lacks inversion symmetry, then a ``quantum'' dipole moment does exist. The system of 3 electrons in D_3h symmetry is our example. This system is realized in triatomic molecules like Na_3. Unlike the fixed dipole of a molecule like water, the quantum moment does not point in a fixed direction, but lies in the plane of the molecule and takes quantized values +/- mu_0 along any direction of measurement in the plane. An electric field F in the plane leads to a linear Stark splitting +/- mu_0 F}. We introduce a toy model to study the effect of Jahn-Teller distortions on the quantum dipole moment. We find that the quantum dipole property survives when the dynamic Jahn-Teller effect is included, if the distortion of the molecule is small. Linear Stark splittings are suppressed in low fields by molecular rotation, just as the linear Stark shift of water is suppressed, but will be revealed in moderately large applied fields and low temperatures. Coulomb correlations also give a partial suppression.Comment: 10 pages with 7 figures included; thoroughly revised with a new coauthor; final minor change