Water-Mediated Charge Transfer and Electron Localization in a Co3Fe2 Cyanide-Bridged Trigonal Bipyramidal Complex.

The effects of temperature and chemical environment on a pentanuclear cyanide-bridged, trigonal bipyramidal molecular paramagnet have been investigated. Using element- and oxidation state-specific near-ambient pressure X-ray photoemission spectroscopy (NAP-XPS) to probe charge transfer and second order, nonlinear vibrational spectroscopy, which is sensitive to symmetry changes based on charge (de)localization coupled with DFT, a detailed picture of environmental effects on charge-transfer-induced spin transitions is presented. The molecular cluster, Co3Fe2(tmphen)6(µ-CN)6(t-CN)6, abbrev. Co3Fe2, shows changes in electronic behavior depending on the chemical environment. NAP-XPS shows that temperature changes induce a metal-to-metal charge transfer (MMCT) in Co3Fe2 between a Co and Fe center, while cycling between ultrahigh vacuum and 2 mbar of water at constant temperature causes oxidation state changes not fully captured by the MMCT picture. Sum frequency generation vibrational spectroscopy (SFG-VS) probes the role of the cyanide ligand, which controls the electron (de)localization via the superexchange coupling. Spectral shifts and intensity changes indicate a change from a charge delocalized, Robin-Day class II/III high spin state to a charge-localized, class I low spin state consistent with DFT. In the presence of a H-bonding solvent, the complex adopts a localized electronic structure, while removal of the solvent delocalizes the charges and drives an MMCT. This change in Robin-Day classification of the complex as a function of chemical environment results in reversible switching of the dipole moment, analogous to molecular multiferroics. These results illustrate the important role of the chemical environment and solvation on underlying charge and spin transitions in this and related complexes.

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration.

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application ("Fragme∩t"), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

Fragment-Based Calculations of Enzymatic Thermochemistry Require Dielectric Boundary Conditions.

Electronic structure calculations on enzymes require hundreds of atoms to obtain converged results, but fragment-based approximations offer a cost-effective solution. We present calculations on enzyme models containing 500-600 atoms using the many-body expansion, comparing to benchmarks in which the entire enzyme-substrate complex is described at the same level of density functional theory. When the amino acid fragments contain ionic side chains, the many-body expansion oscillates under vacuum boundary conditions but rapid convergence is restored using low-dielectric boundary conditions. This implies that full-system calculations in the gas phase are inappropriate benchmarks for assessing errors in fragment-based approximations. A three-body protocol retains sub-kilocalorie per mole fidelity with respect to a supersystem calculation, as does a two-body calculation combined with a full-system correction at a low-cost level of theory. These protocols pave the way for application of high-level quantum chemistry to large systems via rigorous, ab initio treatment of many-body polarization.