Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userExtending Density Matrix Embedding: A Static Two-Particle Theory
323
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We introduce Extended Density Matrix Embedding Theory (EDMET), a static quantum embedding theory explicitly self-consistent with respect to two-body environmental interactions. This overcomes the biggest practical and conceptual limitation of more traditional one-body embedding methods, namely the lack of screening and treatment of long-range correlations. This algebraic zero-temperature embedding augments the correlated cluster with a minimal number of bosons from the random phase approximation, and admits an analytic approach to build a self-consistent Coulomb-exchange-correlation kernel. For extended Hubbard models with non-local interactions, this leads to the accurate description of phase transitions, static quantities and dynamics. We also move towards {em ab initio} systems via the Parriser--Parr--Pople model of conjugated coronene derivatives, finding good agreement with experimental optical gaps.
rate research
Read More
Quantum embedding based on the (one-electron reduced) density matrix is revisited by means of the unitary Householder transformation. While being exact and equivalent to (but formally simpler than) density matrix embedding theory (DMET) in the non-interacting case, the resulting Householder transformed density matrix functional embedding theory (Ht-DMFET) preserves, by construction, the single-particle character of the bath when electron correlation is introduced. In Ht-DMFET, the projected impurity+bath clusters Hamiltonian (from which approximate local properties of the interacting lattice can be extracted) becomes an explicit functional of the density matrix. In the spirit of single-impurity DMET, we consider in this work a closed (two-electron) cluster constructed from the full-size non-interacting density matrix. When the (Householder transformed) interaction on the bath site is taken into account, per-site energies obtained for the half-filled one-dimensional Hubbard lattice match almost perfectly the exact Bethe Ansatz results in all correlation regimes. In the strongly correlated regime, the results deteriorate away from half-filling. This can be related to the electron number fluctuations in the (two-site) cluster which are not described neither in Ht-DMFET nor in regular DMET. As expected, the per-site energies dramatically improve when increasing the number of embedded impurities. Formal connections with density/density matrix functional theories have been briefly discussed and should be explored further. Work is currently in progress in this direction.
Based on recent progress on fermionic exchange symmetry we propose a way to develop new functionals for reduced density matrix functional theory. For some settings with an odd number of electrons, by assuming saturation of the inequalities stemming from the generalized Pauli principle, the many-body wave-function can be written explicitly in terms of the natural occupation numbers and natural orbitals. This leads to an expression for the two-particle density matrix and therefore for the correlation energy functional. This functional was then tested for a three-electron Hubbard model where it showed excellent performance both in the weak and strong correlation regimes.
The idea of using fragment embedding to circumvent the high computational scaling of accurate electronic structure methods while retaining high accuracy has been a long-standing goal for quantum chemists. Traditional fragment embedding methods mainly focus on systems composed of weakly correlated parts and are insufficient when division across chemical bonds is unavoidable. Recently, density matrix embedding theory (DMET) and other methods based on the Schmidt decomposition have emerged as a fresh approach to this problem. Despite their success on model systems, these methods can prove difficult for realistic systems because they rely on either a rigid, non-overlapping partition of the system or a specification of some special sites (i.e. `edge and `center sites), neither of which is well-defined in general for real molecules. In this work, we present a new Schmidt decomposition-based embedding scheme called Incremental Embedding that allows the combination of arbitrary overlapping fragments without the knowledge of edge sites. This method forms a convergent hierarchy in the sense that higher accuracy can be obtained by using fragments involving more sites. The computational scaling for the first few levels is lower than that of most correlated wave function methods. We present results for several small molecules in atom-centered Gaussian basis sets and demonstrate that Incremental Embedding converges quickly with fragment size and recovers most static correlation in small basis sets even when truncated at the second lowest level.
Quantum embedding approaches involve the self-consistent optimization of a local fragment of a strongly correlated system, entangled with the wider environment. The `energy-weighted density matrix embedding theory (EwDMET) was established recently as a way to systematically control the resolution of the fragment-environment coupling, and allow for true quantum fluctuations over this boundary to be self-consistently optimized within a fully static framework. In this work, we reformulate the algorithm to ensure that EwDMET can be considered equivalent to an optimal and rigorous truncation of the self-consistent dynamics of dynamical mean-field theory (DMFT). A practical limitation of these quantum embedding approaches is often a numerical fitting of a self-consistent object defining the quantum effects. However, we show here that in this formulation, all numerical fitting steps can be entirely circumvented, via an effective Dyson equation in the space of truncated dynamics. This provides a robust and analytic self-consistency for the method, and an ability to systematically and rigorously converge to DMFT from a static, wave function perspective. We demonstrate that this improved approach can solve the correlated dynamics and phase transitions of the Bethe lattice Hubbard model in infinite dimensions, as well as one- and two-dimensional Hubbard models where we clearly show the benefits of this rapidly convergent basis for correlation-driven fluctuations. This systematically truncated description of the effective dynamics of the problem also allows access to quantities such as Fermi liquid parameters and renormalized dynamics, and demonstrates a numerically efficient, systematic convergence to the zero-temperature dynamical mean-field theory limit.
We examine the performance of the density matrix embedding theory (DMET) recently proposed in [G. Knizia and G. K.-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)]. The core of this method is to find a proper one-body potential that generates a good trial wave function for projecting a large scale original Hamiltonian to a local subsystem with a small number of basis. The resultant ground state of the projected Hamiltonian can locally approximate the true ground state. However, the lack of the variational principle makes it difficult to judge the quality of the choice of the potential. Here we focus on the entanglement spectrum (ES) as a judging criterion; accurate evaluation of the ES guarantees that the corresponding reduced density matrix well reproduces all physical quantities on the local subsystem. We apply the DMET to the Hubbard model on the one-dimensional chain, zigzag chain, and triangular lattice and test several variants of potentials and cost functions. It turns out that ES serves as a more sensitive quantity than the energy and double occupancy to probe the quality of the DMET outcomes. A symmetric potential reproduces the ES of the phase that continues from a noninteracting limit. The Mott transition as well as symmetry-breaking transitions can be detected by the singularities in the ES. However, the details of the ES in the strongly interacting parameter region depends much on these variants, meaning that the present DMET algorithm allowing for numerous variant is insufficient to fully characterize the particular phases that require characterization by the ES.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
