We present a new approach for numerical solutions of ab initio quantum chemistry systems. The main idea of the approach, which we call canonical diagonalization, is to diagonalize directly the second-quantized Hamiltonian by a sequence of numerical canonical transformations.
We propose a new, alternative method for ab-initio calculations of the electronic structure of solids, which has been specifically adapted to treat many-body effects in a more rigorous way than many existing ab-initio methods. We start from a standard band-structure calculation for an effective one-particle Hamiltonian approximately describing the material of interest. This yields a suitable set of one-particle basis functions, from which well localized Wannier functions can be constructed using a method proposed by Marzari and Vanderbilt. Within this (minimal) basis of localized Wannier functions the matrix elements of the non-interacting (one-particle) Hamiltonian as well as the Coulomb matrix elements can be calculated. The result is a many-body Hamiltonian in second quantization with parameters determined from first principles calculations for the material of interest. The Hamiltonian is in the form of a multi-band Hamiltonian in second quantization (a kind of extended, multi-band Hubbard model) such that all the standard many-body methods can be applied. We explicitly show how this approach can be solved in the simplest many-body approximation, the mean-field Hartree-Fock approximation (HFA), which takes into account exact exchange and corrects for self-interaction effects.
The integration of time-dependent quantum mechanical wave equations is a fundamental problem in computational physics and computational chemistry. The wave-function's energy spectrum as well as its momentum spectrum impose fundamental limits on the performance of numerical algorithms for the solution of wave equations. We demonstrate how canonical transforms may be applied to negotiate these limitations and to increase the performance of numerical algorithms by up to several orders of magnitude. Our approach includes the so-called Kramers-Henneberger transform as a special case and puts forward modifications toward an improved numerical efficiency.
We present a math-physics modeling approach called canonical quantization with numerical mode decomposition for capturing the physics of how incoming photons interact with finite-sized dispersive media, which is not describable by the previous Fano-diagonalization methods. The main procedure is to (1) study a system where electromagnetic fields are coupled to nonuniformly distributed Lorentz oscillators in Hamiltonian mechanics, (2) derive a generalized Hermitian eigenvalue problem for conjugate pairs in coordinate space, (3) apply computational electromagnetics methods to find a countably finite set of time-harmonic eigenmodes that diagonalizes the Hamiltonian, and (4) perform the subsequent canonical quantization with mode decomposition. Moreover, we provide several numerical simulations that capture the physics of full quantum effects, impossible by classical Maxwell's equations, such as nonlocal dispersion cancellation of an entangled photon pair and the Hong-Ou-Mandel effect in a dispersive beam splitter.
Continuous unitary transformations are a powerful tool to extract valuable information out of quantum many-body Hamiltonians, in which the so-called flow equation transforms the Hamiltonian to a diagonal or block-diagonal form in second quantization. Yet, one of their main challenges is how to approximate the infinitely-many coupled differential equations that are produced throughout this flow. Here we show that tensor networks offer a natural and non-perturbative truncation scheme in terms of entanglement. The corresponding scheme is called “entanglement-CUT” or eCUT. It can be used to extract the low-energy physics of quantum many-body Hamiltonians, including quasiparticle energy gaps. We provide the general idea behind eCUT and explain its implementation for finite 1d systems using the formalism of matrix product operators. We also present proof-of-principle results for the spin-(1/2) 1d quantum Ising model and the 3-state quantum Potts model in a transverse field. Entanglement-CUTs can also be generalized to higher dimensions and to the thermodynamic limit.