No Arabic abstract
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.
In most nuclear many-body methods, observables are calculated using many-body wave functions explicitly. The variational two-particle reduced density matrix method is one of the few exceptions to the rule. Ground-state energies of both closed-shell and open-shell nuclear systems can indeed be evaluated by minimizing a constrained linear functional of the two-particle reduced density matrix. However, it has virtually never been used in nuclear theory, because nuclear ground states were found to be well overbound, contrary to those of atoms and molecules. Consequently, we introduced new constraints in the nuclear variational two-particle reduced density matrix method, developed recently for atomic and molecular systems. Our calculations then show that this approach can provide a proper description of nuclear systems where only valence neutrons are included. For the nuclear systems where both neutrons and protons are active, however, the energies obtained with the variational two-particle reduced density matrix method are still overbound. The possible reasons for the noticed discrepancies and solutions to this problem will be discussed.
By combining a parameterized Hermitian matrix, the realignment matrix of the bipartite density matrix $rho$ and the vectorization of its reduced density matrices, we present a family of separability criteria, which are stronger than the computable cross norm or realignment (CCNR) criterion. With linear contraction methods, the proposed criteria can be used to detect the multipartite entangled states that are biseparable under any bipartite partitions. Moreover, we show by examples that the presented multipartite separability criteria can be more efficient than the corresponding multipartite realignment criterion based on CCNR, multipartite correlation tensor criterion and multipartite covariance matrix criterion.
Transient absorption is a very powerful observable in attosecond experiments on atoms, molecules and solids and is frequently used in experiments employing phase-locked few-cycle infrared and XUV laser pulses derived from high harmonic generation. We show numerically and analytically that in non-centrosymmetric systems, such as many polyatomic molecules, which-way interference enabled by the lack of parity conservation leads to new spectral absorption features, which directly reveal the laser electric field. The extension of Attosecond Transient Absorption Spectroscopy (ATAS) to such targets hence becomes sensitive to global and local inversion symmetry. We anticipate that ATAS will find new applications in non-centrosymmetric systems, in which the carrier-to-envelope phase of the infrared pulse becomes a relevant parameter and in which the orientation of the sample and the electronic symmetry of the molecule can be addressed.
A new ATSP2K module is presented for evaluating the electron density function of any multiconfiguration Hartree-Fock or configuration interaction wave function in the non relativistic or relativistic Breit-Pauli approximation. It is first stressed that the density function is not a priori spherically symmetric in the general open shell case. Ways of building it as a spherical symmetric function are discussed, from which the radial electron density function emerges. This function is written in second quantized coupled tensorial form for exploring the atomic spherical symmetry. The calculation of its expectation value is performed using the angular momentum theory in orbital, spin, and quasispin spaces, adopting a generalized graphical technique. The natural orbitals are evaluated from the diagonalization of the density matrix.
We study various methods to generate ensembles of random density matrices of a fixed size N, obtained by partial trace of pure states on composite systems. Structured ensembles of random pure states, invariant with respect to local unitary transformations are introduced. To analyze statistical properties of quantum entanglement in bi-partite systems we analyze the distribution of Schmidt coefficients of random pure states. Such a distribution is derived in the case of a superposition of k random maximally entangled states. For another ensemble, obtained by performing selective measurements in a maximally entangled basis on a multi--partite system, we show that this distribution is given by the Fuss-Catalan law and find the average entanglement entropy. A more general class of structured ensembles proposed, containing also the case of Bures, forms an extension of the standard ensemble of structureless random pure states, described asymptotically, as N to infty, by the Marchenko-Pastur distribution.