Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals

o-Benzyne fragmentation and isomerization pathways: a CASPT2 study

The mechanism for the thermal fragmentation of o-benzyne to C₄H₂ + C₂H₂ and C₆H₂.

Differentiating between trivalent lanthanides and actinides

The reactions of LnCl(3) with molten boric acid result in the formation of Ln[B(4)O(6)(OH)(2)Cl] (Ln = La-Nd), Ln(4)[B(18)O(25)(OH)(13)Cl(3)] (Ln = Sm, Eu), or Ln[B(6)O(9)(OH)(3)] (Ln = Y, Eu-Lu). The reactions of AnCl(3) (An = Pu, Am, Cm) with molten boric acid under the same conditions yield Pu[B(4)O(6)(OH)(2)Cl] and Pu(2)[B(13)O(19)(OH)(5)Cl(2)(H(2)O)(3)], Am[B(9)O(13)(OH)(4)]·H(2)O, or Cm(2)[B(14)O(20)(OH)(7)(H(2)O)(2)Cl]. These compounds possess three-dimensional network structures where rare earth borate layers are joined together by BO(3) and/or BO(4) groups. There is a shift from 10-coordinate Ln(3+) and An(3+) cations with capped triangular cupola geometries for the early members of both series to 9-coordinate hula-hoop geometries for the later elements. Cm(3+) is anomalous in that it contains both 9- and 10-coordinate metal ions. Despite these materials being synthesized under identical conditions, the two series do not parallel one another. Electronic structure calculations with multireference, CASSCF, and density functional theory (DFT) methods reveal the An 5f orbitals to be localized and predominately uninvolved in bonding. For the Pu(III) borates, a Pu 6p orbital is observed with delocalized electron density on basal oxygen atoms contrasting the Am(III) and Cm(III) borates, where a basal O 2p orbital delocalizes to the An 6d orbital. The electronic structure of the Ce(III) borate is similar to the Pu(III) complexes in that the Ce 4f orbital is localized and noninteracting, but the Ce 5p orbital shows no interaction with the coordinating ligands. Natural bond orbital and natural population analyses at the DFT level illustrate distinctive larger Pu 5f atomic occupancy relative to Am and Cm 5f, as well as unique involvement and occupancy of the An 6d orbitals.

Large-scale symmetry-adapted perturbation theory computations via density fitting and Laplace transformation techniques: investigating the fundamental forces of DNA-intercalator interactions

Symmetry-adapted perturbation theory (SAPT) provides a means of probing the fundamental nature of intermolecular interactions. Low-orders of SAPT (here, SAPT0) are especially attractive since they provide qualitative (sometimes quantitative) results while remaining tractable for large systems. The application of density fitting and Laplace transformation techniques to SAPT0 can significantly reduce the expense associated with these computations and make even larger systems accessible. We present new factorizations of the SAPT0 equations with density-fitted two-electron integrals and the first application of Laplace transformations of energy denominators to SAPT. The improved scalability of the DF-SAPT0 implementation allows it to be applied to systems with more than 200 atoms and 2800 basis functions. The Laplace-transformed energy denominators are compared to analogous partial Cholesky decompositions of the energy denominator tensor. Application of our new DF-SAPT0 program to the intercalation of DNA by proflavine has allowed us to determine the nature of the proflavine-DNA interaction. Overall, the proflavine-DNA interaction contains important contributions from both electrostatics and dispersion. The energetics of the intercalator interaction are are dominated by the stacking interactions (two-thirds of the total), but contain important contributions from the intercalator-backbone interactions. It is hypothesized that the geometry of the complex will be determined by the interactions of the intercalator with the backbone, because by shifting toward one side of the backbone, the intercalator can form two long hydrogen-bonding type interactions. The long-range interactions between the intercalator and the next-nearest base pairs appear to be negligible, justifying the use of truncated DNA models in computational studies of intercalation interaction energies.

An overlap fitted chain of spheres exchange method

The "chain of spheres" (COS) algorithm, as part of the RIJCOSX SCF procedure, approximates the exchange term by performing analytic integration with respect to the coordinates of only one of the two electrons, whereas for the remaining coordinates, integration is carried out numerically. In the present work, we attempt to enhance the efficiency of the method by minimizing numerical errors in the COS procedure. The main idea is based on the work of Friesner and consists of finding a fitting matrix, Q, which leads the numerical and analytically evaluated overlap matrices to coincide. Using Q, the evaluation of exchange integrals can indeed be improved. Improved results and timings are obtained with the present default grid setup for both single point calculations and geometry optimizations. The fitting procedure results in a reduction of grid sizes necessary for achieving chemical accuracy. We demonstrate this by testing a number of grids and comparing results to the fully analytic and the earlier COS approximations. This turns out to be favourable for total and reaction energies, for which chemical accuracy can now be reached with a corresponding ~30% speedup over the original RIJCOSX procedure for single point energies. Results are slightly less favourable for the accuracy of geometry optimizations, but the procedure is still shown to yield geometries with errors well below the method inherent errors of the employed theoretical framework.

Assessing metal-metal multiple bonds in Cr-Cr, Mo-Mo, and W-W compounds and a hypothetical U-U compound: a quantum chemical study comparing DFT and multireference methods

To gain insights into the trends in metal-metal multiple bonding among the Group 6 elements, density functional theory has been employed in combination with multiconfigurational methods (CASSCF and CASPT2) to investigate a selection of bimetallic, multiply bonded compounds. For the compound [Ar-MM-Ar] (Ar=2,6-(C(6)H(5))(2)-C(6)H(3), M=Cr, Mo, W) the effect of the Ar ligand on the M(2) core has been compared with the analogous [Ph-MM-Ph] (Ph=phenyl, M=Cr, Mo, W) compounds. A set of [M(2)(dpa)(4)] (dpa=2,2'-dipyridylamide, M=Cr, Mo, W, U) compounds has also been investigated. All of the compounds studied here show important multiconfigurational behavior. For the Mo(2) and W(2) compounds, the σ(2)π(4)δ(2) configuration dominates the ground-state wavefunction, contributing at least 75%. The Cr(2) compounds show a more nuanced electronic structure, with many configurations contributing to the ground state. For the Cr, Mo, and W compounds the electronic absorption spectra have been studied, combining density functional theory and multireference methods to make absorption feature assignments. In all cases, the main features observed in the visible spectra may be assigned as charge-transfer bands. For all compounds investigated the Mayer bond order (MBO) and the effective bond order (EBO) were calculated by density functional theory and CASSCF methods, respectively. The MBO and EBO values share a similar trend toward higher values at shorter normalized metal-metal bond lengths.