Search

The ab initio molecular dynamics (AIMD) method provides a computational route for the real-time simulation of reactive chemistry. An often-overlooked capability of this approach is the opportunity to examine the electronic evolution of a chemical system. For AIMD trajectories based on Hartree-Fock or density functional theory methods, the real-time evolution of single-particle molecular orbitals (MOs) can provide detailed insights into the time-dependent electronic structure of molecules. The evolving electronic Hamiltonians at each MD step pose problems for tracking and visualizing a given MO's character, ordering, and associated phase throughout an MD trajectory, however. This report presents and assesses a simple algorithm for correcting these deficiencies by exploiting similarity projections of the electronic structure between neighboring MD steps. Two aspects bring this analysis beyond a simple step-to-step projection scheme. First, the challenging case of coincidental orbital degeneracies is resolved via a quadrupole-field perturbation that nonetheless rigorously preserves energy conservation. Second, the resulting orbitals are shown to evolve adiabatically, in spite of the "preservation of character" concept that undergirds a projection of neighboring steps' MOs. The method is tested on water clusters, which exhibit considerable dynamic degeneracies, as well as a classic organic nucleophilic substitution reaction, in which the adiabatic evolution of the bonding orbitals clarifies textbook interpretations of the electronic structure during this reactive collision.

The copper hydroxide ion, CuOH+, serves as the catalytic core in several recently developed water-splitting catalysts, and an understanding of its chemistry is critical to determining viable cataly...

Chemical activation of water by a single chlorine atom was examined computationally for clusters of chlorine radicals and water in a size regime just prior to internal hydration of water/ions, Cl·(H2O)n=4-8,17. This investigation follows a recent analysis of this radical-molecule interaction [Christensen et al. J. Phys. Chem. A 2019, 123, 8657] for n = 1-4, which demonstrated that n = 4 marked a transition in which an oxidized-water structural motif became viable, albeit high in energy. Thousands of unique isomers were computed in the present analysis, which resulted in three structural classes of isomers, including intact hydrated chlorine, hydrogen-transferred (HCl)(OH·)(H2O)n-1, and charge-transferred (Cl-)(H3O+)(OH·)(H2O)n-2 configurations. The electronic structures of these classes were investigated, along with harmonic vibrational signatures that probed the degree of water-network perturbations and generated experimentally verifiable computational predictions. The main outcome of this analysis is that the charge-transferred isomers were stabilized considerably upon increased hydration-leading to an energetic crossover with the hydrogen-transferred forms-but the degree of hydration was surprisingly still not sufficient to achieve crossover between the intact chlorine-water complexes and these charge-separated configurations. Internal hydration of the ions appears to be necessary in order to achieve this separation, which will likely occur at larger cluster sizes.

The gas-phase vibrational spectrum of the isolated iodide–water cluster ion (I–·H2O), first reported in 1996, presents one of the most difficult, long-standing spectroscopic puzzles involving ion m...

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design. This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl•(H2O)n=1–4. Using an automated isomer-search procedure, dozen...

The 1,3-diaza-2,4-diborobutane (NBNB) molecule serves as the smallest model complex of an intramolecular "dihydrogen bond," which involves a nominally hydrogen-bonding interaction between amine and borane hydrogen atoms. In the present study, the role of this dihydrogen bond in influencing the inherent molecular dynamics of NBNB is investigated computationally with ab initio molecular dynamics and path integral molecular dynamics techniques, as well as vibrational spectra analysis and static quantum chemistry computations. These simulations indicate that the dihydrogen-bonding interaction impacts both the high- and low-frequency motions of the molecule, with the dominant motions involving low-frequency backbone isomerization and terminal amine rotation. Geometric isotope effects were found to be modest. The analysis also addresses the paradoxical fostering of amine rotation via a relatively strong dihydrogen bond interaction. Electrostatic and geometric factors most directly explain this effect, and although some orbital evidence was found for a small covalent component of this interaction, the dynamics and electronic structure suggest that electrostatic contributions are the controlling factors for molecular motion in NBNB.

The copper hydroxide ion, CuOH

Chemical activation of water by a single chlorine atom was examined computationally for clusters of chlorine radicals and water in a size regime just prior to internal hydration of water/ions, Cl

PURPOSE: Although Common Variable Immunodeficiency (CVID) is considered the most prevalent symptomatic primary antibody deficiency (PAD), there is a population with symptomatic PADs that do not meet criteria for CVID. We analyzed clinical and immunological profiles of patients with different PADs to better understand the differences and similarities between CVID and other PADs. METHODS: We extracted clinical and laboratory data of patients with PADs from electronic medical records. Patients were categorized into CVID, IgG subclass 2 deficiency (IgG2D), IgG deficiency (IgGD) and selective antibody deficiency (sAbD) based on basal immunoglobulin levels and pneumococcal vaccine responses. We compared clinical and immunological characteristics in these groups. RESULTS: All patients, regardless of PAD types, showed similar frequencies of infections, bronchiectasis, and interstitial lung disease (ILD). Hematopoietic malignancies were more frequently found in the CVID than in the IgG2D, IgGD and sAbD groups, while the latter groups trended towards an increased frequency of connective tissue diseases (CTD). Low counts of natural killer (NK) cells were associated with malignancy, autoimmunity, and ILD in CVID but not in other PAD groups. CONCLUSIONS: Higher frequency of hematopoietic malignancy in CVID than in the other PADs and association of lower NK cell counts with non-infectious complications in CVID suggest a relationship between immune alterations and the development of non-infectious manifestations in PADs.

The partial chemical activation of water by reactive radicals was examined computationally for small clusters of chlorine and water, Cl

The isomers of a hydrated Cu(I) ion with n = 1–10 water molecules were investigated by using ab initio quantum chemistry and an automated isomer-search algorithm. The electronic structure and vibrational spectra of the hundreds of resulting isomers were used to analyze the source of the observed bonding patterns. A structural evolution from dominantly two-coordinate structures (n = 1–4) toward a mixture of two- and three-coordinate structures was observed at n = 5–6, where the stability provided by expanded hydrogen-bonding was competitive with the dominantly electrostatic interaction between the water ligand and remaining binding sites of the metal ion. Further hydration (n = 7–10) led to a mixture of three- and four-coordinate structures. The metal ion was found, through spectroscopic signatures, to appreciably perturb the O–H bonds of even third-shell water molecules, which highlighted the ability of this nominally simple ion to partially activate the surrounding water network.

Acceleration of ab initio molecular dynamics (AIMD) simulations can be reliably achieved by extrapolation of electronic data from previous timesteps. Existing techniques utilize polynomial least-squares regression to fit previous steps’ Fock or density matrix elements. In this work, the recursive Burg ‘linear prediction’ technique is shown to be a viable alternative to polynomial regression, and the extrapolation-predicted Fock matrix elements were three orders of magnitude closer to converged elements. Accelerations of 1.8–3.4× were observed in test systems, and in all cases, linear prediction outperformed polynomial extrapolation. Importantly, these accelerations were achieved without reducing the MD integration timestep.

The structures, properties, and spectroscopic signatures of oxidized water clusters,(H2O)+n=6–21, are examined in this work, to provide fundamental insight into renewable energy and radiological processes. Computational quantum chemistry approaches are employed to sample cluster morphologies, yielding hundreds of low-lying isomers with low barriers to interconversion. The ion–radical pair-separation trend, however, which was observed in previous computational studies and in small-cluster spectroscopy experiments, is shown to continue in this larger cluster size regime. The source of this trend is preferential solvation of the hydronium ion by water, including effects beyond the first solvation shell. The fundamental conclusion of this work, therefore, is that the initially formed ion–radical dimer, which has served as a prototypical model of oxidized water, is a nascent species in large, oxidized water clusters and, very likely, bulk water.

A quantitative description of the interactions between ions and water is key to characterizing the role played by ions in mediating fundamental processes that take place in aqueous environments. At the molecular level, vibrational spectroscopy provides a unique means to probe the multidimensional potential energy surface of small ion−water clusters. In this study, we combine the MB-nrg potential energy functions recently developed for ion−water interactions with perturbative corrections to vibrational self-consistent field theory and the local-monomer approximation to disentangle many-body effects on the stability and vibrational structure of the Cs+(H2O)3 cluster. Since several low-energy, thermodynamically accessible isomers exist for Cs+(H2O)3, even small changes in the description of the underlying potential energy surface can result in large differences in the relative stability of the various isomers. Our analysis demonstrates that a quantitative account for three-body energies and explicit treatment of cross-monomer vibrational couplings are required to reproduce the experimental spectrum.

This work describes an approach to accelerate ab initio Born-Oppenheimer molecular dynamics (MD) simulations by exploiting the inherent timescale separation between contributions from different atom-centered Gaussian basis sets. Several MD steps are propagated with a cost-efficient, low-level basis set, after which a dynamical correction accounts for large basis set relaxation effects in a time-reversible fashion. This multiple-timestep scheme is shown to generate valid MD trajectories, on the basis of rigorous testing for water clusters, the methanol dimer, an alanine polypeptide, protonated hydrazine, and the oxidized water dimer. This new approach generates observables that are consistent with those of target basis set trajectories, including MD-based vibrational spectra. This protocol is shown to be valid for Hartree-Fock, density functional theory, and second-order Møller-Plesset perturbation theory approaches. Recommended pairings include 6-31G as a low-level basis set for 6-31G** or 6-311G**, as well as cc-pVDZ as the subset for accurate dynamics with aug-cc-pVTZ. Demonstrated cost savings include factors of 2.6-7.3 on the systems tested and are expected to remain valid across system sizes.

Ionized water clusters serve as a model of water-splitting chemistry for energetic purposes, as well as postradiolytic events in condensed-phase systems. Structures, properties, and relative energies are presented for oxidized water clusters, (H2O)n=1-5(+), using equation-of-motion coupled-cluster theory approaches. In small clusters, an ion-radical contact pair OH···H3O+ is known to form upon ionization. The transition from n = 4 to n = 5 molecules in the cluster, however, is found to demarcate a size regime in which a propensity for the ion and radical to separate exists. This trend is consistent with recent experimental vibrational analyses. Decomposition of the cluster energetics reveals that preferential solvation of the hydronium cation by water serves as the dominant driving force for this pair separation, which should persist in larger clusters and bulk water ionization.

The dynamic, quantum structure of [Mg(H2)n=1–6]2+complexes is investigated via ab initio path integral molecular dynamics simulations. These complexes represent the strong, σ-complex regime of meta...

In combined quantum mechanical/molecular mechanical (QM/MM) free energy calculations, it is often advantageous to have a frozen geometry for the quantum mechanical (QM) region. For such multiple-environment single-system (MESS) cases, two schemes are proposed here for estimating the polarization energy: the first scheme, termed MESS-E, involves a Roothaan step extrapolation of the self-consistent field (SCF) energy; whereas the other scheme, termed MESS-H, employs a Newton-Raphson correction using an approximate inverse electronic Hessian of the QM region (which is constructed only once). Both schemes are extremely efficient, because the expensive Fock updates and SCF iterations in standard QM/MM calculations are completely avoided at each configuration. They produce reasonably accurate QM/MM polarization energies: MESS-E can predict the polarization energy within 0.25 kcal/mol in terms of the mean signed error for two of our test cases, solvated methanol and solvated β-alanine, using the M06-2X or ωB97X-D functionals; MESS-H can reproduce the polarization energy within 0.2 kcal/mol for these two cases and for the oxyluciferin-luciferase complex, if the approximate inverse electronic Hessians are constructed with sufficient accuracy.

Because of both experimental and computational challenges, protonated tryptophan has remained the last aromatic amino acid for which the intrinsic structures of low-energy conformers have not been unambiguously solved. The IR-IR-UV hole-burning spectroscopy technique has been applied to overcome the limitations of the commonly used IR-UV double resonance technique and to measure conformer-specific vibrational spectra of TrpH(+), cooled to T = 10 K. Anharmonic ab initio vibrational spectroscopy simulations unambiguously assign the dominant conformers to the two lowest-energy geometries from benchmark coupled-cluster structure computations. The match between experimental and ab initio spectra provides an unbiased validation of the calculated structures of the two experimentally observed conformers of this benchmark ion. Furthermore, the vibrational spectra provide conformer-specific signatures of the stabilizing interactions, including hydrogen bonding and an intramolecular cation-π interaction.

It's the Network Numerous reactions of small molecules and ions in the atmosphere take place in the confines of watery aerosols. Relph et al. (p. 308 ; see the Perspective by Siefermann and Abel ) explored the specific influence of a water cluster's geometry on the transformation of solvated nitrosonium (NO + ) to nitrous acid (HONO). The reaction involves (O)N–O(H) bond formation with one water molecule, concomitant with proton transfer to additional, surrounding water molecules. Vibrational spectroscopy and theoretical simulations suggest that certain arrangements of the surrounding water network are much more effective than others in accommodating this charge transfer, and thus facilitating the reaction.

Ground state structures with side-on nitrosyl ({eta}{sup 2}-NO) and isonitrosyl (ON) ligands have been observed in a variety of transition-metal complexes. In contrast, excited state structures with bent-NO ligands have been proposed for years but never directly observed. Here we use picosecond time-resolved infrared spectroscopy and density functional theory (DFT) modeling to study the photochemistry of Co(CO){sub 3}(NO), a model transition-metal-NO compound. Surprisingly, we have observed no evidence for ON and {eta}{sup 2}-NO structural isomers, but have observed two bent-NO complexes. DFT modeling of the ground and excited state potentials indicates that the bent-NO complexes correspond to triplet excited states. Photolysis of Co(CO){sub 3}(NO) with a 400-nm pump pulse leads to population of a manifold of excited states which decay to form an excited state triplet bent-NO complex within 1 ps. This structure relaxes to the ground triplet state in ca. 350 ps to form a second bent-NO structure.

A multiple-timestep ab initio molecular dynamics scheme based on varying the two-electron integral screening method used in Hartree–Fock or density functional theory calculations is presented. Although screening is motivated by numerical considerations, it is also related to separations in the length- and timescales characterizing forces in a molecular system: Loose thresholds are sufficient to describe fast motions over short distances, while tight thresholds may be employed for larger length scales and longer times, leading to a practical acceleration of ab initio molecular dynamics simulations. Standard screening approaches can lead, however, to significant discontinuities in (and inconsistencies between) the energy and gradient when the screening threshold is loose, making them inappropriate for use in dynamics. To remedy this problem, a consistent window-screening method that smooths these discontinuities is devised. Further algorithmic improvements reuse electronic-structure information within the dynamics step and enhance efficiency relative to a naı̈ve multiple-timestepping protocol. The resulting scheme is shown to realize meaningful reductions in the cost of Hartree–Fock and B3LYP simulations of a moderately large system, the protonated sarcosine/glycine dipeptide embedded in a 19-water cluster. This document is the unedited Author's version of a Submitted Work that was subsequently accepted for publication in Journal of Chemical Theory and Computation, copyright © American Chemical Society after peer review. To access the final edited and published work see http://pubs.acs.org/articlesonrequest/AOR-3PynRTiASCcZCZep2V3h.

International audience; A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Moller-Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr-2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube.

Analytic gradients of dual-basis Hartree-Fock and density functional theory energies have been derived and implemented, which provide the opportunity for capturing large basis-set gradient effects at reduced cost. Suggested pairings for gradient calculations are 6-31G/6-31G**, dual[-f,-d]/cc-pVTZ, and 6-311G*/6-311 + +G(3df,3pd). Equilibrium geometries are produced within 0.0005 A of large-basis results for the latter two pairings. Though a single, iterative SCF response equation must be solved (unlike standard SCF gradients), it may be obtained in the smaller basis set, and integral screening further reduces the cost for well-chosen subsets. Total nuclear force calculations exhibit up to 75% savings, relative to large-basis calculations.

On-the-fly, ab initio classical molecular dynamics are demonstrated with an underlying dual basis set potential energy surface. Dual-basis self-consistent field (Hartree-Fock and density functional theory) and resolution-of-the-identity second-order Moller-Plesset perturbation theory (RI-MP2) dynamics are tested for small systems, including the water dimer. The resulting dynamics are shown to be faithful representations of their single-basis analogues for individual trajectories, as well as vibrational spectra. Computational cost savings of 58% are demonstrated for SCF methods, even relative to Fock-extrapolated dynamics, and savings are further increased to 71% with RI-MP2. Notably, these timings outperform an idealized estimate of extended-Lagrangian molecular dynamics. The method is subsequently demonstrated on the vibrational absorption spectrum of two NO(+)(H₂O)₃ isomers and is shown to recover the significant width of the shared-proton bands observed experimentally.

Accurate potential energy surfaces for benzene.M complexes (M = Li+, Na+, K+, and NH4+) are obtained using coupled-cluster theory through perturbative triple excitations, CCSD(T). Our computations show that off-axis cation-pi interactions, where the cation is not directly above the aromatic ring, can be favorable and may influence molecular recognition. Even perpendicular, side-on interactions retain 18-32% of their pi-face interaction energy in the gas phase, making their bond strengths comparable to hydrogen bonds in the gas phase. Solvent effects have been explored for each complex using the polarizable continuum model.

For a system which undergoes electron or energy transfer in a polar solvent, we define the diabatic states to be the initial and final states of the system, before and after the nonequilibrium transfer process. We consider two models for the system-solvent interactions: A solvent which is linearly polarized in space and a solvent which responds linearly to the system. From these models, we derive two new schemes for obtaining diabatic states from ab initio calculations of the isolated system in the absence of solvent. These algorithms resemble standard approaches for orbital localization, namely, the Boys and Edmiston-Ruedenberg (ER) formalisms. We show that Boys localization is appropriate for describing electron transfer [Subotnik et al., J. Chem. Phys. 129, 244101 (2008)] while ER describes both electron and energy transfer. Neither the Boys nor the ER methods require definitions of donor or acceptor fragments and both are computationally inexpensive. We investigate one chemical example, the case of oligomethylphenyl-3, and we provide attachment/detachment plots whereby the ER diabatic states are seen to have localized electron-hole pairs.

This article presents a general computational approach for efficient simulations of anharmonic vibrational spectra in chemical systems. An automated local-mode vibrational approach is presented, which borrows techniques from localized molecular orbitals in electronic structure theory. This approach generates spatially localized vibrational modes, in contrast to the delocalization exhibited by canonical normal modes. The method is rigorously tested across a series of chemical systems, ranging from small molecules to large water clusters and a protonated dipeptide. It is interfaced with exact, grid-based approaches, as well as vibrational self-consistent field methods. Most significantly, this new set of reference coordinates exhibits a well-behaved spatial decay of mode couplings, which allows for a systematic, a priori truncation of mode couplings and increased computational efficiency. Convergence can typically be reached by including modes within only about 4 Å. The local nature of this truncation suggests particular promise for the ab initio simulation of anharmonic vibrational motion in large systems, where connection to experimental spectra is currently most challenging.

A time-reversible, multiple-timestep protocol is presented for ab initio molecular dynamics simulations using correlated, wavefunction-based underlying potentials. The method is motivated by the observation that electron correlation contributions to forces vary on a slower timescale than their Hartree-Fock counterparts. An efficient dynamics algorithm, involving short-timestep Hartree-Fock and long-timestep Moøller-Plesset perturbation theory, is presented and tested. Results indicate stable trajectories and relative speedups comparable to those seen in force field-based multiple-timestep schemes, with the highest efficiency improvement occurring for large systems.