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:2005.13498 · The European Physical Journal Plus

Electron-light interaction in nonequilibrium: exact diagonalization for time-dependent Hubbard Hamiltonians

Michael Innerberger, Paul Worm, Paul Prauhart, Anna Kauch

Published Nov 1, 2020• 14 citations
arXivHTMLPDF

Abstract

We present a straightforward implementation scheme for solving the time-dependent Schrödinger equation for systems described by the Hubbard Hamiltonian with time-dependent hoppings. The computations can be performed for clusters of up to 14 sites with, in principle, general geometry. For the time evolution, we use the exponential midpoint rule, where the exponentials are computed via a Krylov subspace method, which only uses matrix-vector multiplication. The presented implementation uses standard libraries for constructing sparse matrices and for linear algebra. Therefore, the approach is easy to use on both desktop computers and computational clusters. We apply the method to calculate time evolution of double occupation and nonequilibrium spectral function of a photo-excited Mott-insulator. The results show that not only the double occupation increases due to creation of electron-hole pairs but also the Mott gap becomes partially filled.

DOI
10.1140/epjp/s13360-020-00919-2

Related papers

Time-dependent spectral properties of a photoexcited one-dimensional ionic Hubbard model: an exact diagonalization study

Junichi Okamoto

Abstract Motivated by the recent progress in time-resolved nonequilibrium spectroscopy in condensed matter, we study an optically excited one-dimensional ionic Hubbard model by exact diagonalization. The model is relevant to organic crystals, transition metal oxides, or ultracold atoms in optical lattices. We implement numerical pump-probe measurements to calculate time-dependent conductivity and single-particle spectral functions. In general, short optical excitation induces a metallic behavior imprinted as a Drude peak in conductivity or an in-gap density of states. In a Mott insulator, we find that the induced Drude peak oscillates at the pump frequency and its second harmonic. The former comes from the oscillation of currents, and the latter from the interference of single- and three-photon excited states. In a band insulator, the Drude peak oscillates only at the pump frequency, and quantities such as the double occupancy do not oscillate. The absence of the second harmonic oscillation is due to the degeneracy of multi-photon excited states. The in-gap density of states in both insulators correlates with the Drude weight and the energy absorption for weak pumping. Strong pumping leads to saturation of the in-gap density of states and to suppression of the Drude weight in the Mott regime. We have also checked that the above features are robust for insulators in the intermediate parameter range. Our study demonstrates the distinct natures of the multi-photon excited states in two different insulators.

arXiv:1909.02471•New Journal of Physics•Nov 27, 2019•11 citations

Nonequilibrium variational-cluster approach to real-time dynamics in the Fermi-Hubbard model

Felix Hofmann, Martin Eckstein, Michael Potthoff

The nonequilibrium variational-cluster approach is applied to study the real-time dynamics of the double occupancy in the one-dimensional Fermi-Hubbard model after different fast changes of hopping parameters. A simple reference system, consisting of isolated Hubbard dimers, is used to discuss different aspects of the numerical implementation of the approach in the general framework of nonequilibrium self-energy functional theory. Opposed to a direct solution of the Euler equation, its time derivative is found to serve as numerically tractable and stable conditional equation to fix the time-dependent variational parameters.

arXiv:1510.05866•Journal of Physics Conference Series•Mar 1, 2016•6 citations

Nonequilibrium variational-cluster approach to real-time dynamics in the Fermi-Hubbard model

Hofmann, F., Martin Eckstein, Potthoff, M.

The nonequilibrium variational-cluster approach is applied to study the real-time dynamics of the double occupancy in the one-dimensional Fermi-Hubbard model after different fast changes of hopping parameters. A simple reference system, consisting of isolated Hubbard dimers, is used to discuss different aspects of the numerical implementation of the approach in the general framework of nonequilibrium self-energy functional theory. Opposed to a direct solution of the Euler equation, its time derivative is found to serve as numerically tractable and stable conditional equation to fix the time-dependent variational parameters.

•Apr 12, 2016•6 citations

Thermalization of a pump-excited Mott insulator

Martin Eckstein, Philipp Werner

We use nonequilibrium dynamical mean-field theory in combination with a recently implemented strong-coupling impurity solver to investigate the relaxation of a Mott insulator after a laser excitation with frequency comparable to the Hubbard gap. The time evolution of the double occupancy exhibits a crossover from a strongly damped transient at short times toward an exponential thermalization at long times. In the limit of strong interactions, the thermalization time is consistent with the exponentially small decay rate for artificially created doublons, which was measured in ultracold atomic gases. When the interaction is comparable to the bandwidth, however, the double occupancy thermalizes within a few times the inverse bandwidth along a rapid thermalization path in which the exponential tail is absent. Similar behavior can be observed in time-resolved photoemission spectroscopy. Our results show that a simple quasi-equilibrium description of the electronic state breaks down for pump-excited Mott insulators characterized by strong interactions.

arXiv:1103.0454•Physical Review B•Jul 25, 2011•132 citations

Photoexcitations in the Hubbard model: Generalized Loschmidt amplitude analysis of impact ionization in small clusters

Clemens Watzenböck, Markus Wallerberger, Laurenz Ruzicka +3 more

We study photoexcitations in small Hubbard clusters of up to 12 sites. After the electric field pulse some of these clusters show an increase of the double occupation through impact ionization. We treat the time-dependent electromagnetic field classically and calculate time evolution by exact diagonalization. As a tool for better analyzing the out-of-equilibrium dynamics, we generalize the Loschmidt amplitude. In this way we are able to resolve which many-body energy eigenstates are responsible for impact ionization and which ones show pronounced changes in the double occupation and spin energy. Our analysis reveals that the increase of spin energy is of little importance for impact ionization. We further demonstrate that, for one-dimensional chains, the optical conductivity has a characteristic peak structure originating solely from vertex corrections.

arXiv:2112.00685•Physical review. B./Physical review. B•Aug 26, 2022•5 citations

Time evolution with the density-matrix renormalization-group algorithm: A generic implementation for strongly correlated electronic systems

Gonzalo Álvarez, Luis G. G. V. Dias da Silva, E. Ponce +1 more

A detailed description of the time-step-targeting time evolution method within the density-matrix renormalization-group algorithm is presented. The focus of this publication is on the implementation of the algorithm and its generic application. The case of one-site excitations within a Hubbard model is analyzed as a test for the algorithm, using open chains and two-leg ladder geometries. The accuracy of the procedure in the case of the recently discussed holon-doublon photo excitations of Mott insulators is also analyzed. Performance and parallelization issues are discussed. In addition, the full open-source code is provided as supplementary material.

arXiv:1103.5391•Physical Review E•Nov 21, 2011•20 citations
We present a straightforward implementation scheme for solving the time-dependent Schrödinger equation for systems described by the Hubbard Hamiltonian with time-dependent hoppings. The computations can be performed for clusters of up to 14 sites with, in principle, general geometry. For the time evolution, we use the exponential midpoint rule, where the exponentials are computed via a Krylov subspace method, which only uses matrix-vector multiplication. The presented implementation uses standard libraries for constructing sparse matrices and for linear algebra. Therefore, the approach is easy to use on both desktop computers and computational clusters. We apply the method to calculate time evolution of double occupation and nonequilibrium spectral function of a photo-excited Mott-insulator. The results show that not only the double occupation increases due to creation of electron-hole pairs but also the Mott gap becomes partially filled.
Abstract Motivated by the recent progress in time-resolved nonequilibrium spectroscopy in condensed matter, we study an optically excited one-dimensional ionic Hubbard model by exact diagonalization. The model is relevant to organic crystals, transition metal oxides, or ultracold atoms in optical lattices. We implement numerical pump-probe measurements to calculate time-dependent conductivity and single-particle spectral functions. In general, short optical excitation induces a metallic behavior imprinted as a Drude peak in conductivity or an in-gap density of states. In a Mott insulator, we find that the induced Drude peak oscillates at the pump frequency and its second harmonic. The former comes from the oscillation of currents, and the latter from the interference of single- and three-photon excited states. In a band insulator, the Drude peak oscillates only at the pump frequency, and quantities such as the double occupancy do not oscillate. The absence of the second harmonic oscillation is due to the degeneracy of multi-photon excited states. The in-gap density of states in both insulators correlates with the Drude weight and the energy absorption for weak pumping. Strong pumping leads to saturation of the in-gap density of states and to suppression of the Drude weight in the Mott regime. We have also checked that the above features are robust for insulators in the intermediate parameter range. Our study demonstrates the distinct natures of the multi-photon excited states in two different insulators.
The nonequilibrium variational-cluster approach is applied to study the real-time dynamics of the double occupancy in the one-dimensional Fermi-Hubbard model after different fast changes of hopping parameters. A simple reference system, consisting of isolated Hubbard dimers, is used to discuss different aspects of the numerical implementation of the approach in the general framework of nonequilibrium self-energy functional theory. Opposed to a direct solution of the Euler equation, its time derivative is found to serve as numerically tractable and stable conditional equation to fix the time-dependent variational parameters.
The nonequilibrium variational-cluster approach is applied to study the real-time dynamics of the double occupancy in the one-dimensional Fermi-Hubbard model after different fast changes of hopping parameters. A simple reference system, consisting of isolated Hubbard dimers, is used to discuss different aspects of the numerical implementation of the approach in the general framework of nonequilibrium self-energy functional theory. Opposed to a direct solution of the Euler equation, its time derivative is found to serve as numerically tractable and stable conditional equation to fix the time-dependent variational parameters.
We use nonequilibrium dynamical mean-field theory in combination with a recently implemented strong-coupling impurity solver to investigate the relaxation of a Mott insulator after a laser excitation with frequency comparable to the Hubbard gap. The time evolution of the double occupancy exhibits a crossover from a strongly damped transient at short times toward an exponential thermalization at long times. In the limit of strong interactions, the thermalization time is consistent with the exponentially small decay rate for artificially created doublons, which was measured in ultracold atomic gases. When the interaction is comparable to the bandwidth, however, the double occupancy thermalizes within a few times the inverse bandwidth along a rapid thermalization path in which the exponential tail is absent. Similar behavior can be observed in time-resolved photoemission spectroscopy. Our results show that a simple quasi-equilibrium description of the electronic state breaks down for pump-excited Mott insulators characterized by strong interactions.
We study photoexcitations in small Hubbard clusters of up to 12 sites. After the electric field pulse some of these clusters show an increase of the double occupation through impact ionization. We treat the time-dependent electromagnetic field classically and calculate time evolution by exact diagonalization. As a tool for better analyzing the out-of-equilibrium dynamics, we generalize the Loschmidt amplitude. In this way we are able to resolve which many-body energy eigenstates are responsible for impact ionization and which ones show pronounced changes in the double occupation and spin energy. Our analysis reveals that the increase of spin energy is of little importance for impact ionization. We further demonstrate that, for one-dimensional chains, the optical conductivity has a characteristic peak structure originating solely from vertex corrections.
A detailed description of the time-step-targeting time evolution method within the density-matrix renormalization-group algorithm is presented. The focus of this publication is on the implementation of the algorithm and its generic application. The case of one-site excitations within a Hubbard model is analyzed as a test for the algorithm, using open chains and two-leg ladder geometries. The accuracy of the procedure in the case of the recently discussed holon-doublon photo excitations of Mott insulators is also analyzed. Performance and parallelization issues are discussed. In addition, the full open-source code is provided as supplementary material.