Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStructure of the first order Reduced Density Matrix in three electron systems: A Generalized Pauli Constraints assisted study
54
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We investigate the structure of the one-body Reduced Density Matrix (1RDM) of three electron systems, i.e. doublet and quadruplet spin configurations, corresponding to the smallest interacting system with an open-shell ground state. To this end, we use Configuration Interaction (CI) expansions of the exact wave function in Slater determinants built from natural orbitals in a finite dimensional Hilbert space. With the exception of maximally polarized systems, the natural orbitals of spin eigenstates are generally spin dependent, i.e. the spatial parts of the up and down natural orbitals form two different sets. A measure to quantify this spin dependence is introduced and it is shown that it varies by several orders of magnitude depending on the system. We also study the ordering issue of the spin-dependent occupation numbers which has practical implications in Reduced Density Matrix Functional Theory minimization schemes when Generalized Pauli Constraints are imposed and in the form of the CI expansion in terms of the natural orbitals. Finally, we discuss the aforementioned CI expansion when there are GPCs that are almost pinned.
rate research
Read More
Functionals of the one-body reduced density matrix (1-RDM) are routinely minimized under Colemans ensemble $N$-representability conditions. Recently, the topic of pure-state $N$-representability conditions, also known as generalized Pauli constraints, received increased attention following the discovery of a systematic way to derive them for any number of electrons and any finite dimensionality of the Hilbert space. The target of this work is to assess the potential impact of the enforcement of the pure-state conditions on the results of reduced density-matrix functional theory calculations. In particular, we examine whether the standard minimization of typical 1-RDM functionals under the ensemble $N$-representability conditions violates the pure-state conditions for prototype 3-electron systems. We also enforce the pure-state conditions, in addition to the ensemble ones, for the same systems and functionals and compare the correlation energies and optimal occupation numbers with those obtained by the enforcement of the ensemble conditions alone.
An active space variational calculation of the 2-electron reduced density matrix (2-RDM) is derived and implemented where the active orbitals are correlated within the pair approximation. The pair approximation considers only doubly occupied configurations of the wavefunction which enables the calculation of the 2-RDM at a computational cost of $mathcal{O}(r^3)$. Calculations were performed both with the pair active space configuration interaction (PASCI) method and the pair active space self consistent field (PASSCF) method. The latter includes a mixing of the active and inactive orbitals through unitary transformations. The active-space pair 2-RDM method is applied to the nitrogen molecule, the p-benzyne diradical, a newly synthesized BisCobalt complex, and the nitrogenase cofactor FeMoco. The FeMoco molecule is treated in a [120,120] active space. Fractional occupations are recovered in each of these systems, indicating the detection and recovery of strong electron correlation.
We consider necessary conditions for the one-body-reduced density matrix (1RDM) to correspond to a triplet wave-function of a two electron system. The conditions concern the occupation numbers and are different for the high spin projections, $S_z=pm 1$, and the $S_z=0$ projection. Hence, they can be used to test if an approximate 1RDM functional yields the same energies for both projections. We employ these conditions in reduced density matrix functional theory calculations for the triplet excitations of two-electron systems. In addition, we propose that these conditions can be used in the calculation of triplet states of systems with more than two electrons by restricting the active space. We assess this procedure in calculations for a few atomic and molecular systems. We show that the quality of the optimal 1RDMs improves by applying the conditions in all the cases we studied.
In this work, we simulate the electron dynamics in molecular systems with the Time-Dependent Density Matrix Renormalization Group (TD-DMRG) algorithm. We leverage the generality of the so-called tangent-space TD-DMRG formulation and design a computational framework in which the dynamics is driven by the exact non-relativistic electronic Hamiltonian. We show that, by parametrizing the wave function as a matrix product state, we can accurately simulate the dynamics of systems including up to 20 electrons and 32 orbitals. We apply the TD-DMRG algorithm to three problems that are hardly targeted by time-independent methods: the calculation of molecular (hyper)polarizabilities, the simulation of electronic absorption spectra, and the study of ultrafast ionization dynamics.
We present a matrix-product state (MPS)-based quadratically convergent density-matrix renormalization group self-consistent-field (DMRG-SCF) approach. Following a proposal by Werner and Knowles (JCP 82, 5053, (1985)), our DMRG-SCF algorithm is based on a direct minimization of an energy expression which is correct to second-order with respect to changes in the molecular orbital basis. We exploit a simultaneous optimization of the MPS wave function and molecular orbitals in order to achieve quadratic convergence. In contrast to previously reported (augmented Hessian) Newton-Raphson and super-configuration-interaction algorithms for DMRG-SCF, energy convergence beyond a quadratic scaling is possible in our ansatz. Discarding the set of redundant active-active orbital rotations, the DMRG-SCF energy converges typically within two to four cycles of the self-consistent procedure
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
