The Use of Effective Core Potentials with Multiconfiguration Pair-Density Functional Theory.

The reliable and accurate prediction of chemical properties is a key goal in quantum chemistry. Transition-metal-containing complexes can often pose difficulties to quantum mechanical methods for multiple reasons, including many electron configurations contributing to the overall electronic description of the system and the large number of electrons significantly increasing the amount of computational resources required. Often, multiconfigurational electronic structure methods are employed for such systems, and the cost of these calculations can be reduced by the use of an effective core potential (ECP). In this work, we explore both theoretical considerations and performances of ECPs applied in the context of multiconfiguration pair-density functional theory (MC-PDFT). A mixed-basis set approach is used, using ECP basis sets for transition metals and all-electron basis sets for nonmetal atoms. We illustrate the effects that an ECP has on the key parameters used in the computation of MC-PDFT energies, and we explore how ECPs affect the prediction of physical observables for chemical systems. The dissociation curve for a metal dimer was explored, and ionization energies for transition metal-containing diatomic systems were computed and compared to experimental values. In general, we find that ECP approaches employed with MC-PDFT are able to predict ionization energies with improved accuracy compared to traditional Kohn-Sham density functional theory approaches.

Assessing the importance of multireference correlation in predicting reversed conductance decay.

In a classical electronic resistor, conductance decays as the device length increases according to Ohm's Law. While most molecular series display a comparable exponential decay in conductance with increasing molecular length, a class of single-molecule device series exists where conductance instead increases with molecular/device length, a phenomenon called reversed conductance decay. While reversals of conductance decay have been repeatedly theoretically predicted, they have been far more difficult to demonstrate experimentally. Previous studies have suggested that theoretical multi-reference(static) correlation errors may be a major cause of this discrepancy, yet most single-molecule transport methods are unable to treat multireference correlation. Using our unique multireference transport method based on non-equilibrium Green's function and multiconfigurational pair-density functional theory (NEGF-MCPDFT), we examined a previously predicted case of reversed conductance decay in systems of linear chains of phenyl rings with varying lengths and electrode designs. We compare our NEGF-MCPDFT results to those of non-multireference NEGF methods to quantify the exact role of static correlation in conductance decay reversals and clarify their relative importance to geometric and electrode design/coupling considerations.

A multiconfiguration pair-density functional theory-based approach to molecular junctions.

Due to their small size and unique properties, single-molecule electronics have long seen research interest from experimentalists and theoreticians alike. From a theoretical standpoint, modeling these systems using electronic structure theory can be difficult due to the importance of electron correlation in the determination of molecular properties, and this electron correlation can be computationally expensive to consider, particularly multiconfigurational correlation energy. In this work, we develop a new approach for the study of single-molecule electronic systems, denoted NEGF-MCPDFT, which combines multiconfiguration pair-density functional theory (MCPDFT) with the non-equilibrium Green's function formalism (NEGF). The use of MCPDFT with NEGF allows for the efficient inclusion of both static and dynamic electron correlations in the description of the junction's electronic structure. Complete active space self-consistent field wave functions are used as references in the MCPDFT calculation, and as with any active space method, effort must be made to determine the proper orbital character to include in the active space. We perform conductance and transmission calculations on a series of alkanes (predominantly single-configurational character) and benzyne (multiconfigurational character), exploring the role that active space selection has on the computed results. For the alkane junctions explored (where dynamic electron correlation dominates), the MCPDFT-NEGF results agree well with the DFT-NEGF results. For the benzyne junction (which has a significant static correlation), we see clear differences in the MCPDFT-NEGF and DFT-NEGF results and evidence that NEGF-MCPDFT is capturing additional electron correlation effects beyond those provided by the Perdew-Burke-Ernzerhof functional.

Acetyl nitrate mediated conversion of methyl ketones to diverse carboxylic acid derivatives.

The development of a novel acetyl nitrate mediated oxidative conversion of methyl ketones to carboxylic acid derivatives is described. By analogy to the haloform reaction and supported by experimental and computational investigation we propose a mechanism for this transformation.

Development and application of a 2-electron reduced density matrix approach to electron transport via molecular junctions.

Can an electronic device be constructed using only a single molecule? Since this question was first asked by Aviram and Ratner in the 1970s [Chem. Phys. Lett. 29, 277 (1974)], the field of molecular electronics has exploded with significant experimental advancements in the understanding of the charge transport properties of single molecule devices. Efforts to explain the results of these experiments and identify promising new candidate molecules for molecular devices have led to the development of numerous new theoretical methods including the current standard theoretical approach for studying single molecule charge transport, i.e., the non-equilibrium Green's function formalism (NEGF). By pairing this formalism with density functional theory (DFT), a wide variety of transport problems in molecular junctions have been successfully treated. For some systems though, the conductance and current-voltage curves predicted by common DFT functionals can be several orders of magnitude above experimental results. In addition, since density functional theory relies on approximations to the exact exchange-correlation functional, the predicted transport properties can show significant variation depending on the functional chosen. As a first step to addressing this issue, the authors have replaced density functional theory in the NEGF formalism with a 2-electron reduced density matrix (2-RDM) method, creating a new approach known as the NEGF-RDM method. 2-RDM methods provide a more accurate description of electron correlation compared to density functional theory, and they have lower computational scaling compared to wavefunction based methods of similar accuracy. Additionally, 2-RDM methods are capable of capturing static electron correlation which is untreatable by existing NEGF-DFT methods. When studying dithiol alkane chains and dithiol benzene in model junctions, the authors found that the NEGF-RDM predicts conductances and currents that are 1-2 orders of magnitude below those of B3LYP and M06 DFT functionals. This suggests that the NEGF-RDM method could be a viable alternative to NEGF-DFT for molecular junction calculations.

Positive semidefinite tensor factorizations of the two-electron integral matrix for low-scaling ab initio electronic structure.

Tensor factorization of the 2-electron integral matrix is a well-known technique for reducing the computational scaling of ab initio electronic structure methods toward that of Hartree-Fock and density functional theories. The simplest factorization that maintains the positive semidefinite character of the 2-electron integral matrix is the Cholesky factorization. In this paper, we introduce a family of positive semidefinite factorizations that generalize the Cholesky factorization. Using an implementation of the factorization within the parametric 2-RDM method [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)], we study several inorganic molecules, alkane chains, and potential energy curves and find that this generalized factorization retains the accuracy and size extensivity of the Cholesky factorization, even in the presence of multi-reference correlation. The generalized family of positive semidefinite factorizations has potential applications to low-scaling ab initio electronic structure methods that treat electron correlation with a computational cost approaching that of the Hartree-Fock method or density functional theory.

Energies and structures in biradical chemistry from the parametric two-electron reduced-density matrix method: applications to the benzene and cyclobutadiene biradicals.

The treatment of biradical chemistry presents a challenge for electronic structure theory, especially single-reference methods, as it requires the description of varying degrees and kinds of electron correlation. In this work we assess the ability of the parametric two-electron reduced-density matrix (p2-RDM) method to describe biradical chemistry through application to the benzene and cyclobutadiene biradicals. The relative energy of o- and m-benzynes predicted by the p2-RDM method is consistent with Wenthold et al.'s experimental determinations, while the more difficult relative energy prediction of the more multi-referenced p-benzyne is within 1.4 kcal mol(-1) of the experimental value [P. G. Wenthold et al., J. Am. Chem. Soc., 1998, 120, 5279], which is significantly better than traditional single-reference methods. We observe that the degree of multireference correlation in the biradicals depends upon the distance between their radical centers, with the largest radical separation displaying the largest degree of multireference correlation. In addition to relative and absolute electronic energies, we report molecular geometries, natural orbitals, and natural-orbital occupations for the benzene and cyclobutadiene biradicals.

Comparison of low-rank tensor expansions for the acceleration of quantum chemistry computations.

Low-rank spectral expansion and tensor hypercontraction are two promising techniques for reducing the size of the two-electron excitation tensor by factorizing it into products of smaller tensors. Both methods can potentially realize an O(r(4)) quantum chemistry method where r is the number of one-electron orbitals. We compare the two factorizations in this paper by applying them to the parametric 2-electron reduced density matrix method with the M functional [D. A. Mazziotti, Phys. Rev. Lett. 101, 253002 (2008)]. We study several inorganic molecules, alkane chains, and potential curves as well as reaction and dissociation energies. The low-rank spectral expansion, we find, is typically more efficient than tensor hypercontraction due to a faster convergence of the energy and a smaller constant prefactor in the energy optimization. Both factorizations are applicable to the acceleration of a wide range of wavefunction and reduced-density-matrix methods.

Relative energies and geometries of the cis- and trans-HO3 radicals from the parametric 2-electron density matrix method.

The parametric 2-electron reduced density matrix (2-RDM) method employing the M functional [Mazziotti, D. A. Phys. Rev. Lett. 2008, 101, 253002], also known as the 2-RDM(M) method, improves on the accuracy of coupled electron-pair theories including coupled cluster with single-double excitations at the computational cost of configuration interaction with single-double excitations. The cis- and trans-HO(3) isomers along with their isomerization transition state were examined using the recent extension of 2-RDM(M) to nonsinglet open-shell states [Schwerdtfeger, C. A.; Mazziotti, D. A. J. Chem. Phys. 2012, 137, 034107] and several coupled cluster methods. We report the calculated energies, geometries, natural-orbital occupation numbers, and reaction barriers for the HO(3) isomers. We find that the 2-RDM(M) method predicts that the trans isomer of HO(3) is lower in energy than the cis isomer by 1.71 kcal/mol in the correlation-consistent polarized valence quadruple-ζ (cc-pVQZ) basis set and 1.84 kcal/mol in the augmented correlation-consistent polarized valence quadruple-ζ (aug-cc-pVQZ) basis set. Results include the harmonic zero-point vibrational energies calculated in the correlation-consistent polarized valence double-ζ basis set. On the basis of the results of a geometry optimization in the augmented correlation consistent polarized valence triple-ζ basis set, the parametric 2-RDM(M) method predicts a central oxygen-oxygen bond of 1.6187 Å. We compare these energies and geometries to those predicted by three single-reference coupled cluster methods and experimental results and find that the inclusion of multireference correlation is important to describe properly the relative energies of the cis- and trans-HO(3) isomers and improve agreement with experimental geometries.