In this paper, we propose an imaginary-time evolution method to obtain self-consistent static electronic states according to real-time TDDFT calculations. This method provides very good convergence accuracy and stability with respect to the ground state, as well as by a large number of computational steps due to obtain all states to evolve simultaneously. Also, we can get a diffusion equation from this method as same as the real-time TDDFT and can probably unify these time-evolution calculations in the same framework which is used in electronic structure calculations. In this study, we show that this method is effective in practical electronic state calculations. We will show several improvements and calculations with the computational process using the imaginary time evolution method in combination with the conventional DFT algorithm for practical use. The results of calculations of systems such as hydrogen atoms, hydrogen molecules and ethylene can show the computational efficiency and accuracy of this method, which probably improve traditional self-consistent electronic structure calculations for practical systems.
In this paper, we propose an imaginary-time evolution method to obtain self-consistent static electronic states according to real-time TDDFT calculations. This method provides very good convergence accuracy and stability with respect to the ground state, as well as by a large number of computational steps due to obtain all states to evolve simultaneously. Also, we can get a diffusion equation from this method as same as the real-time TDDFT and can probably unify these time-evolution calculations in the same framework which is used in electronic structure calculations. In this study, we show that this method is effective in practical electronic state calculations. We will show several improvements and calculations with the computational process using the imaginary time evolution method in combination with the conventional DFT algorithm for practical use. The results of calculations of systems such as hydrogen atoms, hydrogen molecules and ethylene can show the computational efficiency and accuracy of this method, which probably improve traditional self-consistent electronic structure calculations for practical systems.
Abstract In imaginary-time evolution, the same technique, Crank–Nicolson’s method, as for the analysis of real-time evolution is employed. For the verification using the model oscillators such as quantum well, Morse, and harmonic potentials, practical density functional theory (DFT) calculations are performed for hydrogen atoms and molecules. The convergence is confirmed to be faster under appropriate initial and boundary conditions. The results are consistent with the analytical solutions and those of other DFT calculations. Imaginary-time evolution is quite convenient to subsequent real-time real-space time-dependent DFT calculation. Although the operators are different, similar evolution is adopted. This technique is expected to unify into practical electronic state calculations for the self-consistent electronic structure and system’s response due to perturbation because calculations are in the same framework as the time-dependent equation.
Accurate descriptions of the electronic properties of molecules, solids and other materials are crucial in the development of many modern technologies, and are therefore the subject of much scientific research. In tandem with experimental research, computational simulations play a huge role. Density functional theory (DFT) is the most popular method for performing ground-state electronic structure calculations due to its accuracy in many situations along with its computational efficiency. Its time-dependent extension (TDDFT) allows the modelling of the dynamics of systems of interacting electrons. (TD)DFT is an exact theory in principle, but requires approximations in practice, specifically to the many-body effects of electron exchange and correlation (xc), which are described by some unknown universal functional. Despite much success, it often proves inadequate in many scenarios of increasing interest. Therefore, we aim to contribute to the development of improved approximate density functionals. \n \nIn this thesis, we take the approach of studying small model systems for which we can compute the exact functionals. These systems are designed to exhibit the key features present in realistic quantum systems, such as atoms and molecules, which the exact universal (TD)DFT functional is capable of describing, but are often ill-described by the commonly used approximations. By analysing the nature of these exact functionals, we gain valuable insight into their fundamental properties, together with identifying where the current approximations are lacking, thereby allowing us to suggest how these approximations may be improved.
We explore the usefulness of a quantum fluid dynamics (QFD) approach for quantitative electronic structure calculations of many-electron systems. By combining QFD and density functional theory, a single time-dependent nonlinear QFD equation can be derived. The equation is further transformed into a diffusion-type form by an imaginary-time evolution method, whose asymptotic solution reaches a global minimum and the many-body ground-state wavefunction. The time-dependent generalized pseudospectral method is extended to solve the diffusion equation in spherical coordinates, allowing optimal and nonuniform spatial discretization and accurate and efficient solution of the diffusion function in space and time. The procedure is applied to the study of electronic energies, densities and other ground-state properties of noble gas atoms (He, Ne, Ar, Kr, Xe). The results are in good agreement with other best available values. The method offers a conceptually appealing and computationally practical procedure for the treatment of many-electron systems beyond the Hartree-Fock level.
Abstract Most quantum algorithms designed to generate or probe properties of the ground state of a quantum many-body system require as input an initial state with a large overlap with the desired ground state. One approach for preparing such a ground state is Imaginary Time Evolution (ITE). Recent work by [Motta, M. , Sun, C. , Tan, A.T.K. et al (2020)] introduced an algorithm—which we will refer to as Quantum Imaginary Time Evolution (QITE)—that shows how ITE can be approximated by a sequence of unitary operators, making QITE potentially implementable on early fault-tolerant quantum computers. In this work, we provide a heuristic study of the capabilities of the QITE algorithm in approximating the ITE of lattice and molecular electronic structure Hamiltonians. We numerically study the performance of the QITE algorithm when provided with a good classical initial state for a large class of systems, some of which are of interest to industrial applications, and check if QITE is able to qualitatively replicate the ITE behavior and improve over a classical mean-field solution. The systems we consider in this work range from one- and two-dimensional lattice systems of various lattice geometries displaying short- and long-range interactions, to active spaces of molecular electronic structure Hamiltonians. In addition to the comparison of QITE and ITE, we explicitly show how imaginary time evolved fermionic Gaussian states can serve as initial states which can be efficiently computed on classical computers and efficiently implemented on quantum computers for generic spin Hamiltonians in arbitrary lattice geometries and dimensions, which can be of independent interest.
Zombie states are a recently introduced formalism to describe coupled coherent fermionic states that address the fermionic sign problem in a computationally tractable manner. Previously, it has been shown that Zombie states with fractional occupations of spin orbitals obeyed the correct fermionic creation and annihilation algebra and presented results for real-time evolution [D. V. Shalashilin, J. Chem. Phys. 148, 194109 (2018)]. In this work, we extend and build on this formalism by developing efficient algorithms for evaluating the Hamiltonian and other operators between Zombie states and address their normalization. We also show how imaginary time propagation can be used to find the ground state of a system. We also present a biasing method, for setting up a basis set of random Zombie states, that allows much smaller basis sizes to be used while still accurately describing the electronic structure Hamiltonian and its ground state and describe a technique of wave function "cleaning" that removes the contributions of configurations with the wrong number of electrons, improving the accuracy further. We also show how low-lying excited states can be calculated efficiently using a Gram-Schmidt orthogonalization procedure. The proposed algorithm of imaginary time propagation on biased random grids of Zombie states may present an alternative to the existing quantum Monte Carlo methods.
The path-integral renormalization group and direct energy minimization method of practical first-principles electronic structure calculations for multi-body systems within the framework of the real-space finite-difference scheme are introduced. These two methods can handle higher dimensional systems with consideration of the correlation effect. Furthermore, they can be easily extended to the multicomponent quantum systems which contain more than two kinds of quantum particles. The key to the present methods is employing linear combinations of nonorthogonal Slater determinants (SDs) as multi-body wavefunctions. As one of the noticeable results, the same accuracy as the variational Monte Carlo method is achieved with a few SDs. This enables us to study the entire ground state consisting of electrons and nuclei without the need to use the Born-Oppenheimer approximation. Recent activities on methodological developments aiming towards practical calculations such as the implementation of auxiliary field for Coulombic interaction, the treatment of the kinetic operator in imaginary-time evolutions, the time-saving double-grid technique for bare-Coulomb atomic potentials and the optimization scheme for minimizing the total-energy functional are also introduced. As test examples, the total energy of the hydrogen molecule, the atomic configuration of the methylene and the electronic structures of two-dimensional quantum dots are calculated, and the accuracy, availability and possibility of the present methods are demonstrated.