Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

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.

CH Mode Mixing Determines the Band Shape of the Carboxylate Symmetric Stretch in Apo-EDTA, Ca2+-EDTA, and Mg2+-EDTA

The infrared spectra of EDTA complexed with Ca²⁺ and Mg²⁺ contain, to date, unidentified vibrational bands. This study assigns the peaks in the linear and two-dimensional infrared spectra of EDTA, with and without either Ca²⁺ or Mg²⁺ ions. Two-dimensional infrared spectroscopy and DFT calculations reveal that, in both the presence and absence of ions, the carboxylate symmetric stretch and the terminal CH bending vibrations mix. We introduce a method to calculate participation coefficients that quantify the contribution of the carboxylate symmetric stretch, CH wag, CH twist, and CH scissor in the 1400-1550 cm^-1 region. With the help of participation coefficients, we assign the 1400-1430 cm^-1 region to the carboxylate symmetric stretch, which can mix with CH modes. We assign the 1000-1380 cm^-1 region to CH twist modes, the 1380-1430 cm^-1 region to wag modes, and the 1420-1650 cm^-1 region to scissor modes. The difference in binding geometry between the carboxylate-Ca²⁺ and carboxylate-Mg²⁺ complex manifests as new diagonal and cross-peaks between the mixed modes in the two complexes. The small Mg²⁺ ion binds EDTA tighter than the Ca²⁺ ion, which causes a redshift of the COO symmetric stretches of the sagittal carboxylates. Energy decomposition analysis further characterizes the importance of electrostatics and deformation energy in the bound complexes.

Molecular design principles of ionic liquids with a sulfonyl fluoride moiety

New family of SO₂F-functionalized ionic liquids.

Buckybowls as gas adsorbents: binding of gaseous pollutants and their electric-field induced release

The adsorption of nitric oxide and nitrogen dioxide (NOx) to the Buckybowls sumanene and corannulene was investigated. Binding energies were up to 1.8× larger than for coronene as the planar analogue, demonstrating the advantages of Buckybowls for gas adsorption. In agreement with previous reports on carbon dioxide and methane adsorption, the favorable binding energies for NOx were shown to be associated with the curvature of the Buckybowls. It is shown that applying an electric field along the bowl symmetry axis modifies the bowl curvatures and impacts adsorbate binding energies, including the potential to desorb adsorbates. As a proof of concept, it is shown that applying electric fields of different strengths and orientations selectively controls sumanene's preference to bind nitric oxide, nitrogen dioxide, and carbon dioxide, suggesting potential applications for dynamically tunable gas adsorption. Moreover, it is demonstrated that adsorbates can be desorbed by applying suitable electric field strengths, allowing cleaning of the Buckybowls for renewed usage.

Ligand-Substrate Dispersion Facilitates the Copper-Catalyzed Hydroamination of Unactivated Olefins

The current understanding of ligand effects in transition metal catalysis is mostly based on the analysis of catalyst-substrate through-bond and through-space interactions, with the latter commonly considered to be repulsive in nature. The dispersion interaction between the ligand and the substrate, a ubiquitous type of attractive noncovalent interaction, is seldom accounted for in the context of transition-metal-catalyzed transformations. Herein we report a computational model to quantitatively analyze the effects of different types of catalyst-substrate interactions on reactivity. Using this model, we show that in the copper(I) hydride (CuH)-catalyzed hydroamination of unactivated olefins, the substantially enhanced reactivity of copper catalysts based on bulky bidentate phosphine ligands originates from the attractive ligand-substrate dispersion interaction. These computational findings are validated by kinetic studies across a range of hydroamination reactions using structurally diverse phosphine ligands, revealing the critical role of bulky P-aryl groups in facilitating this process.

Polymerization of silyl ketenes using alkoxide initiators: a combined computational and experimental study

The atomic composition and chemical structure of polymers is fundamental to dictating properties and applications.

Ultrafast vibrational spectroscopy (2D-IR) of CO2 in ionic liquids: Carbon capture from carbon dioxide's point of view

The CO2ν3 asymmetric stretching mode is established as a vibrational chromophore for ultrafast two-dimensional infrared (2D-IR) spectroscopic studies of local structure and dynamics in ionic liquids, which are of interest for carbon capture applications. CO2 is dissolved in a series of 1-butyl-3-methylimidazolium-based ionic liquids ([C4C1im][X], where [X](-) is the anion from the series hexafluorophosphate (PF6 (-)), tetrafluoroborate (BF4 (-)), bis-(trifluoromethyl)sulfonylimide (Tf2N(-)), triflate (TfO(-)), trifluoroacetate (TFA(-)), dicyanamide (DCA(-)), and thiocyanate (SCN(-))). In the ionic liquids studied, the ν3 center frequency is sensitive to the local solvation environment and reports on the timescales for local structural relaxation. Density functional theory calculations predict charge transfer from the anion to the CO2 and from CO2 to the cation. The charge transfer drives geometrical distortion of CO2, which in turn changes the ν3 frequency. The observed structural relaxation timescales vary by up to an order of magnitude between ionic liquids. Shoulders in the 2D-IR spectra arise from anharmonic coupling of the ν2 and ν3 normal modes of CO2. Thermal fluctuations in the ν2 population stochastically modulate the ν3 frequency and generate dynamic cross-peaks. These timescales are attributed to the breakup of ion cages that create a well-defined local environment for CO2. The results suggest that the picosecond dynamics of CO2 are gated by local diffusion of anions and cations.

Piezoelectric Effects of Applied Electric Fields on Hydrogen-Bond Interactions: First-Principles Electronic Structure Investigation of Weak Electrostatic Interactions

The piezoelectric properties of 2-methyl-4-nitroaniline crystals were explored qualitatively and quantitatively using an electrostatically embedded many-body (EE-MB) expansion scheme for the correlation energies of a system of monomers within the crystal. The results demonstrate that hydrogen bonding is an inherently piezoelectric interaction, deforming in response to the electrostatic environment. We obtain piezo-coefficients in excellent agreement with the experimental values. This approach reduces computational cost and reproduces the total resolution of the identity (RI)-Møller-Plesset second-order perturbation theory (RI-MP2) energy for the system to within 1.3 × 10(-5)%. Furthermore, the results suggest novel ways to self-assemble piezoelectric solids and suggest that accurate treatment of hydrogen bonds requires precise electrostatic evaluation. Considering the ubiquity of hydrogen bonds across chemistry, materials, and biology, a new electromechanical view of these interactions is required.

A Linear-Scaling MP2 Method for Large Molecules by Rigorous Integral-Screening Criteria

A brief review of our linear-scaling method for atomic-orbital (AO) second-order Møller-Plesset perturbation theory (MP2) is given. The key feature of our method is the rigorous preselection of numerically significant four-center two-electron integrals based on multipole-based integral estimates (MBIE) that do not only account for the exponential coupling between the Gaussian-type basis functions forming charge distributions, but also for the 1/R coupling between the charge distributions. This coupling turns for the required integral products in AO-MP2 into at least a 1/R4 or even a 1/R6 decay behavior. Using MBIE we attain linear scaling, which is illustrated for DNA fragments with up to 1052 atoms and 10 674 basis functions as computed on a single processor. The largest molecule calculated in our present work at the scaled-opposite spin (SOS-) AO-MP2 level is an RNA system comprising 1664 atoms and 19 182 basis functions. Furthermore, we present results for the use of Cholesky-decomposed pseudo-density matrices in Laplace-based MP2, that offers the advantage of exploiting occupied/virtual blocking both with and without auxiliary basis sets.

Current status of the AMOEBA polarizable force field

Molecular force fields have been approaching a generational transition over the past several years, moving away from well-established and well-tuned, but intrinsically limited, fixed point charge models toward more intricate and expensive polarizable models that should allow more accurate description of molecular properties. The recently introduced AMOEBA force field is a leading publicly available example of this next generation of theoretical model, but to date, it has only received relatively limited validation, which we address here. We show that the AMOEBA force field is in fact a significant improvement over fixed charge models for small molecule structural and thermodynamic observables in particular, although further fine-tuning is necessary to describe solvation free energies of drug-like small molecules, dynamical properties away from ambient conditions, and possible improvements in aromatic interactions. State of the art electronic structure calculations reveal generally very good agreement with AMOEBA for demanding problems such as relative conformational energies of the alanine tetrapeptide and isomers of water sulfate complexes. AMOEBA is shown to be especially successful on protein-ligand binding and computational X-ray crystallography where polarization and accurate electrostatics are critical.

Multipole-based integral estimates for the rigorous description of distance dependence in two-electron integrals

We derive multipole-based integral estimates (MBIE) as rigorous and tight upper bounds to four-center two-electron integrals in order to account for the 1/R distance decay between the charge distributions, which is missing in the Schwarz screening commonly used in ab initio methods. Our screening criteria are valid for all angular momenta and can be formulated for any order of multipoles. We have found the expansion limited to dipoles to be sufficiently tight for estimating the integrals in Hartree-Fock and density-functional theories, while the screening effort is negligible. For, e.g., a DNA fragment with 1052 atoms and 10,674 basis functions (6-31G*) the exchange part is faster by a factor of 2.1 as compared to the Schwarz screening both within our linear exchange scheme, whereas a smaller factor of 1.3 is gained for the Coulomb part within the continuous fast multipole method. Most importantly, our new MBIE screening is perfectly suited to exploit the strong distance decay of electron-correlation effects of at least 1/R4 in atomic-orbital-based formulations of correlation methods.

Rigorous integral screening for electron correlation methods

We derive rigorous multipole-based integral estimates (MBIE) in order to account for the distance dependence occurring in atomic-orbital (AO) formulations of electron correlation theory, where our focus is on AO-MP2 theory within a Laplace scheme. We find for the exact transformed integral products an extremely early onset of a linear-scaling behavior and a very small number of significant products. To preselect the significant integral products we adapt our MBIE method as rigorous upper bound. In this way it is possible to exploit the favorable scaling behavior observed and to reduce the scaling of estimated products asymptotically to linear, without sacrificing accuracy or reliability. By separating Coulomb- and exchange-type contractions only half-transformed integrals need to be computed. Furthermore, our scheme of rigorously preselecting transformed integral products via MBIE seems to offer particularly interesting perspectives for a direct formation of half- or fully transformed integrals by using multipole expansions and auxiliary basis sets.