Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userInterpretation of multiple solutions in fully iterative GF2 and GW schemes using local analysis of two-particle density matrices
149
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Due to non-linear structure, iterative Greens function methods can result in multiple different solutions even for simple molecular systems. In contrast to the wave-function methods, a detailed and careful analysis of such molecular solutions was not performed before. In this work, we use two-particle density matrices to investigate local spin and charge correlators that quantify the charge-resonance and covalent characters of these solutions. When applied within unrestricted orbital set, spin correlators elucidate the broken symmetry of the solutions, containing necessary information for building effective magnetic Hamiltonians. Based on GW and GF2 calculations of simple molecules and transition metal complexes, we construct Heisenberg Hamiltonians, four-spin-four-center corrections, as well as biquadratic spin-spin interactions. These Hamiltonian parametrizations are compared to prior wave-function calculations.
rate research
Read More
One-particle Greens function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Greens functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $langle{S^2}rangle$ and $langle{N^2}rangle$ can be challenging, because they require a solution of the computationally expensive Bethe--Salpeter equation to find two-particle Greens functions. We demonstrate that the solution of the Bethe--Salpeter equation can be complitely avoided. Applying the thermodynamic Hellmann--Feynman theorem to self-consistent one-particle Greens function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $langle{S^2}rangle$ and $langle{N^2}rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.
In this chapter, we provide a review of ground-state Kohn-Sham density-functional theory of electronic systems and some of its extensions, we present exact expressions and constraints for the exchange and correlation density functionals, and we discuss the main families of approximations for the exchange-correlation energy: semilocal approximations, single-determinant hybrid approximations, multideterminant hybrid approximations, dispersion-corrected approximations, as well as orbital-dependent exchange-correlation density functionals. The chapter aims at providing both a consistent birds-eye view of the field and a detailed description of some of the most used approximations. It is intended to be readable by chemists/physicists and applied mathematicians.
We report unphysical irregularities and discontinuities in some key experimentally-measurable quantities computed within the GW approximation of many-body perturbation theory applied to molecular systems. In particular, we show that the solution obtained with partially self-consistent GW schemes depends on the algorithm one uses to solve self-consistently the quasi-particle (QP) equation. The main observation of the present study is that each branch of the self-energy is associated with a distinct QP solution, and that each switch between solutions implies a significant discontinuity in the quasiparticle energy as a function of the internuclear distance. Moreover, we clearly observe ripple effects, i.e., a discontinuity in one of the QP energies induces (smaller) discontinuities in the other QP energies. Going from one branch to another implies a transfer of weight between two solutions of the QP equation. The case of occupied, virtual and frontier orbitals are separately discussed on distinct diatomics. In particular, we show that multisolution behavior in frontier orbitals is more likely if the HOMO-LUMO gap is small.
We introduce the Nuclear Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrodinger equation simultaneously for electrons and other quantum species. In contrast to already existing multicomponent approaches, in this work we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parametrization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate full configuration interaction energies of molecular systems with more than three nuclei and twelve particles in total, which is currently a major challenge for other multicomponent approaches. We present NEAP-DMRG results for two few-body systems, i.e., H$_2$ and H$_3^+$, and one larger system, namely BH$_3$
Describing time-dependent many-body systems where correlation effects play an important role remains a major theoretical challenge. In this paper we develop a time-dependent many-body theory that is based on the two-particle reduced density matrix (2-RDM). We develop a closed equation of motion for the 2-RDM employing a novel reconstruction functional for the three-particle reduced density matrix (3-RDM) that preserves norm, energy, and spin symmetries during time propagation. We show that approximately enforcing $N$-representability during time evolution is essential for achieving stable solutions. As a prototypical test case which features long-range Coulomb interactions we employ the one-dimensional model for lithium hydride (LiH) in strong infrared laser fields. We probe both one-particle observables such as the time-dependent dipole moment and two-particle observables such as the pair density and mean electron-electron interaction energy. Our results are in very good agreement with numerically exact solutions for the $N$-electron wavefunction obtained from the multiconfigurational time-dependent Hartree-Fock method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
