A local density fitting scheme is considered in which atomic orbital (AO) products are approximated using only auxiliary AOs located on one of the nuclei in that product. The possibility of variational collapse to an unphysical "attractive electron" state that can affect such density fitting [P. Merlot, T. Kjærgaard, T. Helgaker, R. Lindh, F. Aquilante, S. Reine, and T. B. Pedersen, J. Comput. Chem. 34, 1486 (2013)] is alleviated by including atom-wise semidiagonal integrals exactly. Our approach leads to a significant decrease in the computational cost of density fitting for Hartree-Fock theory while still producing results with errors 2-5 times smaller than standard, nonlocal density fitting. Our method allows for large Hartree-Fock and density functional theory computations with exact exchange to be carried out efficiently on large molecules, which we demonstrate by benchmarking our method on 200 of the most widely used prescription drug molecules. Our new fitting scheme leads to smooth and artifact-free potential energy surfaces and the possibility of relatively simple analytic gradients.

The potential energy surface for the interaction between benzene and hydroxyl radical is studied in detail using quantum mechanical methods, with a particular focus on the hydrogen abstraction pathway. Geometric parameters are optimized using a variety of density functional methods as well as perturbation theory. Energies are refined using coupled cluster singles and doubles with perturbative triples [CCSD(T)] extrapolated to the complete basis set limit. At our most reliable level of theory, complexation energies are found to be (with zero-point corrected energies in parentheses) 3.7 (2.8) kcal/mol for the benzene-hydroxyl radical complex and 2.9 (-1.7) kcal/mol for the phenyl radical-water complex. The barrier to H abstraction lies 6.5 (4.2) kcal/mol above the infinitely separated benzene and hydroxyl radical monomers.

Explicitly correlated ab initio methods have been used to compute full quartic force fields for the three chain minima for HOOOOH, which are found to lie within 1 kcal mol(-1). The CCSD(T)-F12 method with the cc-pVTZ-F12 basis set was used to compute equilibrium structures, anharmonic vibrational frequencies, and rotational constants for HOOH, HOOOH, and three chain isomers of HOOOOH, with the two former force fields being used as benchmarks for the latter three. The full quartic force fields were computed in such a way as to yield fundamental frequencies for all isotopologues at once. The present research confirms the recent experimental identification of HOOOH and provides reliable force fields in support of future experimental work on the enigmatic bonding paradigms involved in the HOOOOH chain.