Loren Greenman

Strong Correlation in Acene Sheets from the Active-Space Variational Two-Electron Reduced Density Matrix Method: Effects of Symmetry and Size

Polyaromatic hydrocarbons (PAHs) are a class of organic molecules with importance in several branches of science, including medicine, combustion chemistry, and materials science. The delocalized π-orbital systems in PAHs require highly accurate electronic structure methods to capture strong electron correlation. Treating correlation in PAHs has been challenging because (i) traditional wave function methods for strong correlation have not been applicable since they scale exponentially in the number of strongly correlated orbitals, and (ii) alternative methods such as the density-matrix renormalization group and variational two-electron reduced density matrix (2-RDM) methods have not been applied beyond linear acene chains. In this paper we extend the earlier results from active-space variational 2-RDM theory [Gidofalvi, G.; Mazziotti, D. A. J. Chem. Phys. 2008, 129, 134108] to the more general two-dimensional arrangement of rings--acene sheets--to study the relationship between geometry and electron correlation in PAHs. The acene-sheet calculations, if performed with conventional wave function methods, would require wave function expansions with as many as 1.5 × 10(17) configuration state functions. To measure electron correlation, we employ several RDM-based metrics: (i) natural-orbital occupation numbers, (ii) the 1-RDM von Neumann entropy, (iii) the correlation energy per carbon atom, and (iv) the squared Frobenius norm of the cumulant 2-RDM. The results confirm a trend of increasing polyradical character with increasing molecular size previously observed in linear PAHs and reveal a corresponding trend in two-dimensional (arch-shaped) PAHs. Furthermore, in PAHs of similar size they show significant variations in correlation with geometry. PAHs with the strictly linear geometry (chains) exhibit more electron correlation than PAHs with nonlinear geometries (sheets).

Strong electron correlation in the decomposition reaction of dioxetanone with implications for firefly bioluminescence

Dioxetanone, a key component of the bioluminescence of firefly luciferin, is itself a chemiluminescent molecule due to two conical intersections on its decomposition reaction surface. While recent calculations of firefly luciferin have employed four electrons in four active orbitals [(4,4)] for the dioxetanone moiety, a study of dioxetanone [F. Liu et al., J. Am. Chem. Soc. 131, 6181 (2009)] indicates that a much larger active space is required. Using a variational calculation of the two-electron reduced-density-matrix (2-RDM) [D. A. Mazziotti, Acc. Chem. Res. 39, 207 (2006)], we present the ground-state potential energy surface as a function of active spaces from (4,4) to (20,17) to determine the number of molecular orbitals required for a correct treatment of the strong electron correlation near the conical intersections. Because the 2-RDM method replaces exponentially scaling diagonalizations with polynomially scaling semidefinite optimizations, we readily computed large (18,15) and (20,17) active spaces that are inaccessible to traditional wave function methods. Convergence of the electron correlation with active-space size was measured with complementary RDM-based metrics, the von Neumann entropy of the one-electron RDM as well as the Frobenius and infinity norms of the cumulant 2-RDM. Results show that the electron correlation is not correctly described until the (14,12) active space with small variations present through the (20,17) space. Specifically, for active spaces smaller than (14,12), we demonstrate that at the first conical intersection, the electron in the σ(∗) orbital of the oxygen-oxygen bond is substantially undercorrelated with the electron of the σ orbital and overcorrelated with the electron of the carbonyl oxygens p orbital. Based on these results, we estimate that in contrast to previous treatments, an accurate calculation of the strong electron correlation in firefly luciferin requires an active space of 28 electrons in 25 orbitals, beyond the capacity of traditional multireference wave function methods.

Strong correlation in hydrogen chains and lattices using the variational two-electron reduced density matrix method

The variational two-electron reduced-density-matrix (2-RDM) method, scaling polynomially with the size of the system, was applied to linear chains and three-dimensional clusters of atomic hydrogen as large as H(64). In the case of the 4x4x4 hydrogen lattice of 64 hydrogen atoms, a correct description of the dissociation requires about 10(18) equally weighted determinants in the wave function, which is too large for traditional multireference methods. The correct energy in the dissociation limit was obtained from the variational 2-RDM method in contrast to Hartree-Fock and single-reference methods. Analysis of the occupation numbers demonstrates that even for 1.0 A bond distances the presence of strong electron correlation requires a multireference method. Three-dimensional systems exhibit a marked increase in electron correlation from one-dimensional systems regardless of size. The metal-to-insulator transition upon expansion of the clusters was studied using the decay of the 1-RDM off-diagonal elements. The variational 2-RDM method was shown to capture the metal-to-insulator transition and dissociation behavior accurately for all systems.

Decoherence in Attosecond Photoionization

The creation of superpositions of hole states via single-photon ionization using attosecond extreme-ultraviolet pulses is studied with the time-dependent configuration-interaction singles (TDCIS) method. Specifically, the degree of coherence between hole states in atomic xenon is investigated. We find that interchannel coupling not only affects the hole populations, but it also enhances the entanglement between the photoelectron and the remaining ion, thereby reducing the coherence within the ion. As a consequence, even if the spectral bandwidth of the ionizing pulse exceeds the energy splittings among the hole states involved, perfectly coherent hole wave packets cannot be formed. For sufficiently large spectral bandwidth, the coherence can only be increased by increasing the mean photon energy.

Highly multireferenced arynes studied with large active spaces using two-electron reduced density matrices

Using the active-space two-electron reduced density matrix (2-RDM) method, which scales polynomially with the size of the active space [G. Gidofalvi and D. A. Mazziotti, J. Chem. Phys. 129, 134108 (2008)], we were able to use active spaces as large as 24 electrons in 24 orbitals in computing the ground-state energies and properties of highly multireferenced arynes. Because the conventional complete-active-space self-consistent-field (CASSCF) method scales exponentially with the size of the active space, its application to arynes was mainly limited to active spaces of 12 electrons in 12 orbitals. For these smaller active spaces the active-space 2-RDM method accurately reproduces the results of CASSCF. However, we show that the larger active spaces are necessary for describing changes in energies and properties with aryne chain length such as the emergence of polyradical character. Furthermore, the addition of further electron correlation by multireference perturbation theory is demonstrated to be inadequate for removing the limitations of the smaller active spaces.

Energy Barriers of Vinylidene Carbene Reactions from the Anti-Hermitian Contracted Schrödinger Equation

Computational studies of carbenes must take into account the possibility of multireference correlation because the highest occupied and lowest unoccupied molecular orbitals can be nearly energetically degenerate. We apply the anti-Hermitian contracted Schrodinger equation (ACSE) [Mazziotti, D. A. Phys. Rev. Lett. 2006, 97, 143002] to compute two-electron reduced density matrices (2-RDMs) and their energies for two carbene reactions: (i) the acetylene-vinylidene rearrangement and (ii) the rearrangement of pent-1-en-4-yn-3-one to acryloylvinylidene, which then cyclizes to cyclopenta-2,4-dienone. The ACSE has some unique advantages in the treatment of carbene reactions and more general families of reactions in which the importance of multireference correlation is not known a priori: (i) the ACSE is more reliable than single-reference methods for confirming the presence or absence of multireference correlation and (ii) in the absence of multireference correlation, unlike multireference second-order perturbation theory (MRPT2), the ACSE recovers more single-reference correlation energy than similarly scaling coupled-cluster methods. Because MRPT2 does not recover as much single-reference correlation as the coupled-cluster or ACSE methods, it tends to underestimate reaction barriers within the carbene reactions. For example, in the rearrangement of pent-1-en-4-yn-3-one, the ACSE and CCSD(T) methods produce cyclization barriers of 18.9 and 18.7 kcal/mol with the 6-31G(d) basis set, whereas MRPT2 predicts this barrier to be 12.1 kcal/mol; furthermore, both the ACSE and CCSD(T) determine the energy of the transition state for acryloylvinylidene formation to be 6.6-6.7 kcal/mol above that of the carbene, and yet, MRPT2 does not predict a transition state.

Electronic excited-state energies from a linear response theory based on the ground-state two-electron reduced density matrix

Ground-state two-particle reduced density matrices (2-RDMs) are used to calculate excited-state energy spectra. Solving the Schrodinger equation for excited states dominated by single excitations from the ground-state wavefunction requires the ground-state 2- and 3-RDMs. The excited states, however, can be obtained without a knowledge of the ground-state 3-RDM by two methods: (i) cumulant expansion methods which build the 3-RDM from the 2-RDM, and (ii) double commutator methods which eliminate the 3-RDM. Previous work [Mazziotti, Phys. Rev. A 68, 052501 (2003)] examined the accuracy of excited states extracted from ground-state 2-RDMs, which were calculated by full configuration interaction or the variational 2-RDM method. In this work we employ (i) advances in semidefinite programming to treat the excited states of water and hydrogen fluoride and chains of hydrogen atoms, and (ii) the addition of partial three-particle N-representability conditions to compute more accurate ground-state 2-RDMs. With the hydrogen chains we examine the metal-to-insulator transition as measured by the band gap (the difference between the ground-state and the first excited-state energies), which is difficult for excited-state methods to capture.

Impact of multichannel and multipole effects on the Cooper minimum in the high-order-harmonic spectrum of argon

We investigate the relevance of multiple-orbital and multipole effects during high-harmonic generation (HHG). The time-dependent configuration-interaction singles (TDCIS) approach is used to study the impact of the detailed description of the residual electron-ion interaction on the HHG spectrum. We find that the shape and position of the Cooper minimum in the HHG spectrum of argon changes significantly whether or not interchannel interactions are taken into account. The HHG yield can be underestimated by up to 2 orders of magnitude in the energy regio of 30-50 eV. We show that the argument of low ionization probability is not sufficient to justify ignoring multiple-orbital contributions. Additionally, we find the HHG yield is sensitive to the nonspherical multipole character of the electron-ion interaction.

Balancing single- and multi-reference correlation in the chemiluminescent reaction of dioxetanone using the anti-Hermitian contracted Schrödinger equation

Direct computation of energies and two-electron reduced density matrices (2-RDMs) from the anti-Hermitian contracted Schrödinger equation (ACSE) [D. A. Mazziotti, Phys. Rev. Lett. 97, 143002 (2006)], it is shown, recovers both single- and multi-reference electron correlation in the chemiluminescent reaction of dioxetanone especially in the vicinity of the conical intersection where strong correlation is important. Dioxetanone, the light-producing moiety of firefly luciferin, efficiently converts chemical energy into light by accessing its excited-state surface via a conical intersection. Our previous active-space 2-RDM study of dioxetanone [L. Greenman and D. A. Mazziotti, J. Chem. Phys. 133, 164110 (2010)] concluded that correlating 16 electrons in 13 (active) orbitals is required for realistic surfaces without correlating the remaining (inactive) orbitals. In this paper we pursue two complementary goals: (i) to correlate the inactive orbitals in 2-RDMs along dioxetanones reaction coordinate and compare these results with those from multireference second-order perturbation theory (MRPT2) and (ii) to assess the size of the active space-the number of correlated electrons and orbitals-required by both MRPT2 and ACSE for accurate energies and surfaces. While MRPT2 recovers very different amounts of correlation with (4,4) and (16,13) active spaces, the ACSE obtains a similar amount of correlation energy with either active space. Nevertheless, subtle differences in excitation energies near the conical intersection suggest that the (16,13) active space is necessary to determine both energetic details and properties. Strong electron correlation is further assessed through several RDM-based metrics including (i) total and relative energies, (ii) the von Neumann entropy based on the 1-electron RDM, as well as the (iii) infinity and (iv) squared Frobenius norms based on the cumulant 2-RDM.

Robust Control Pulses Design for Electron Shuttling in Solid-State Devices

In this brief, we study robust pulse design for electron shuttling in solid-state devices. This is crucial for many practical applications of coherent quantum mechanical systems. Our objective is to design control pulses that can transport an electron along a chain of donors and that also make this process robust to parameter uncertainties. We formulate this problem here as a set of optimal control problems and derive explicit expressions for the gradients of the aggregate transfer fidelity. Numerical results for a donor chain of ionized phosphorus atoms in bulk silicon demonstrate the efficacy of our algorithm.