Electronic structure density matrices

Keywords Grand-canonical ensemble Electronic structure Density matrix Descriptors... 

Back in the 1960s, hopes for future progress in electronic structure theory were associated with correlated wave function techniques and the tantalizing possibility of variational calculations based on the two-electron reduced density matrix [294]. DFT was not on the quantum chemistry agenda at that time. The progress of wave function techniques has been remarkable, as documented elsewhere in this volume. In contrast, the density matrix approach has not yet materialized into a competitive computational method, despite many persistent efforts [295]. Meanwhile, approximate DFT has become the most widely used method of quantum chemistry, offering an unprecedented accuracy-to-cost ratio. 

The authors are honored to dedicate this article in memory of Isaiah Shavitt whose remarkable contributions transformed electronic structure theory. Observables that depend on pairwise interactions can be directly computed with the two-electron reduced density matrix (2-RDM) without the A-electron wave function [10, 36]. Integration of the A-electron density matrix over all electrons save two yields the 2-RDM... 

Analogically to the representation of the wave-function in structural terms, there is a way to separate (hyper)polarizabilities into the individual contributions from individual atoms. A method for such separation was developed by Bredas [15, 16] and is called the real-space finite-field method. The approach can be easily implemented for a post-Hartree-Fock method in the r-electron approximation due to the simplicity of e calculation of the one-electron reduced density matrix (RDMl) elements. In our calculations we are using a simple munerical-derivative two-points formula for RDMl matrix elements (Z ) [88] (see also [48]) ... 

Iditional importance is that the vibrational modes are dependent upon the reciprocal e vector k. As with calculations of the electronic structure of periodic lattices these cal-ions are usually performed by selecting a suitable set of points from within the Brillouin. For periodic solids it is necessary to take this periodicity into account the effect on the id-derivative matrix is that each element x] needs to be multiplied by the phase factor k-r y). A phonon dispersion curve indicates how the phonon frequencies vary over tlie luin zone, an example being shown in Figure 5.37. The phonon density of states is ariation in the number of frequencies as a function of frequency. A purely transverse ition is one where the displacement of the atoms is perpendicular to the direction of on of the wave in a pmely longitudinal vibration tlie atomic displacements are in the ition of the wave motion. Such motions can be observed in simple systems (e.g. those contain just one or two atoms per unit cell) but for general three-dimensional lattices of the vibrations are a mixture of transverse and longitudinal motions, the exceptions... 

We have extended the linear combination of Gaussian-type orbitals local-density functional approach to calculate the total energies and electronic structures of helical chain polymers[35]. This method was originally developed for molecular systems[36-40], and extended to two-dimensionally periodic sys-tems[41,42] and chain polymers[34j. The one-electron wavefunctions here are constructed from a linear combination of Bloch functions c>>, which are in turn constructed from a linear combination of nuclear-centered Gaussian-type orbitals Xylr) (in ihis case, products of Gaussians and the real solid spherical harmonics). The one-electron density matrix is given by... 

In the DC-biased structures considered here, the dynamics are dominated by electronic states in the conduction band [1]. A simplified version of the theory assumes that the excitation occurs only at zone center. This reduces the problem to an n-level system (where n is approximately equal to the number of wells in the structure), which can be solved using conventional first-order perturbation theory and wave-packet methods. A more advanced version of the theory includes all of the hole states and electron states subsumed by the bandwidth of the excitation laser, as well as the perpendicular k states. In this case, a density-matrix picture must be used, which requires a solution of the time-dependent Liouville equation. Substituting the Hamiltonian into the Liouville equation leads to a modified version of the optical Bloch equations [13,15]. These equations can be solved readily, if the k states are not coupled (i.e., in the absence of Coulomb interactions). 

Millam, J. M., Scuseria, G. E., 1997, Linear Scaling Conjugate Gradient Density Matrix Search as an Alternative to Diagonalization for First Principles Electronic Structure Calculations , J. Chem. Phys., 106, 5569. 

Yang W, Lee T (1995) A density-matrix divide-and-conquer approach for electronic-structure calculations of large molecules. J Chem Phys 103(13) 5674—5678... 

Li, X. P., R. W. Nunes, and D. Vanderbildt. 1993. Density Matrix Electronic Structure Method with Linear System-Size Scaling. Phys. Rev. B 48,10891. 

Perhaps the most fruitful of these studies was the radiolysis of HCo(C0)4 in a Kr matrix (61,62). Free radicals detected in the irradiated material corresponded to processes of H-Co fission, electron capture, H-atom additions and clustering. Initial examination at 77 K or lower temperatures revealed the presence of two radicals, Co(C0)4 and HCo(C0)4 , having similar geometries (IV and V) and electronic structures. Both have practically all of the unpaired spin-density confined to nuclei located on the three-fold axis, in Co 3dz2, C 2s or H Is orbitals. Under certain conditions, a radical product of hydrogen-atom addition, H2Co(C0)3, was observed this species is believed to have a distorted trigonal bipyramidal structure in which the H-atoms occupy apical positions. 

The identification of unknown chemical compounds isolated in inert gas matrices is nowadays facilitated by comparison of the measured IR spectra with those computed at reliable levels of ab initio or density functional theory (DFT). Furthermore, the observed reactivity of matrix isolated species can in some instances be explained with the help of computed reaction energies and barriers for intramolecular rearrangements. Hence, electronic structure methods developed into a useful tool for the matrix isolation community. In this chapter, we will give an overview of the various theoretical methods and their limitations when employed in carbene chemistry. For a more detailed qualitative description of the merits and drawbacks of commonly used electronic structure methods, especially for open-shell systems, the reader is referred to the introductory guide of Bally and Borden.29... 

Since the early days of quantum mechanics, the wave function theory has proven to be very successful in describing many different quantum processes and phenomena. However, in many problems of quantum chemistry and solid-state physics, where the dimensionality of the systems studied is relatively high, ab initio calculations of the structure of atoms, molecules, clusters, and crystals, and their interactions are very often prohibitive. Hence, alternative formulations based on the direct use of the probability density, gathered under what is generally known as the density matrix theory [1], were also developed since the very beginning of the new mechanics. The independent electron approximation or Thomas-Fermi model, and the Hartree and Hartree-Fock approaches are former statistical models developed in that direction [2]. These models can be considered direct predecessors of the more recent density functional theory (DFT) [3], whose principles were established by Hohenberg,... 

In this paper we present the first application of the ZORA (Zeroth Order Regular Approximation of the Dirac Fock equation) formalism in Ab Initio electronic structure calculations. The ZORA method, which has been tested previously in the context of Density Functional Theory, has been implemented in the GAMESS-UK package. As was shown earlier we can split off a scalar part from the two component ZORA Hamiltonian. In the present work only the one component part is considered. We introduce a separate internal basis to represent the extra matrix elements, needed for the ZORA corrections. This leads to different options for the computation of the Coulomb matrix in this internal basis. The performance of this Hamiltonian and the effect of the different Coulomb matrix alternatives is tested in calculations on the radon en xenon atoms and the AuH molecule. In the atomic cases we compare with numerical Dirac Fock and numerical ZORA methods and with non relativistic and full Dirac basis set calculations. It is shown that ZORA recovers the bulk of the relativistic effect and that ZORA and Dirac Fock perform equally well in medium size basis set calculations. For AuH we have calculated the equilibrium bond length with the non relativistic Hartree Fock and ZORA methods and compare with the Dirac Fock result and the experimental value. Again the ZORA and Dirac Fock errors are of the same order of magnitude. 

Crosslinked low-density polyethylene foams with a closedcell structure were investigated using differential scanning calorimetry, scanning electron microscopy, density, and thermal expansion measurements. At room temperature, the coefficient of thermal expansion decreased as the density increased. This was attributed to the influence of gas expansion within the cells. At a given material density, the expansion increased as the cell size became smaller. At higher temperatures, the relationship between thermal expansion and density was more complex, due to physical transitions in the matrix polymer. Materials with high density and thick cell walls were concluded to be the best for low expansion applications. 16 refs. 

Z. Zhao, B. J. Braams, H. Fukuda, M. L. Overton, and J. K. Percus, The reduced density matrix method for electronic structure calculations and the role of three-index representabihty conditions. J. Chem. Phys. 120, 2095 (2004). 

By replacing the wavefunction with a density matrix, the electronic structure problem is reduced in size to that for a two- or three-electron system. Rather than solve the Schrodinger equation to determine the wavefunction, the lower bound method is invoked to determine the density matrix this requires adjusting parameters so that the energy content of the density matrix is minimized. More precisely, the lower bound method requires finding a solution to the energy problem,... 

The RDMs for atoms and molecules have a special structure from the spin of the electrons. To each spatial orbital, we associate a spin of either a or f. Because the two spins are orthogonal upon integration of the N-particle density matrix, only RDM blocks where the net spin of the upper indices equals the net spin of the lower indices do not vanish. Hence a p-RDM is block diagonal with (p -f 1) nonzero blocks. Specifically, the 1-RDM has two nonzero blocks, an a-block and a -block ... 

Abstract The density matrix renormalization group (DMRG) is an electronic structure... 

Raghu, C., Anusooya Pati, Y., Ramasesha, S. Structural and electronic instabilities in polyacenes density-matrix renormalization group study of a long-range interacting model. Phys. Rev. B 2002, 65(15), 155204. 

---