A priori calculation of molecular properties to chemical accuracy

The fundamental vibrational frequencies and spectroscopic constants of the C2O2H2 isomers: molecules known in simulated interstellar ice analogues

While trans-glyoxal may not be easily observable in astronomical sources through either IR or radioastronomy due to its C2h symmetry, its cis conformer along with the cyc-H2COCO epoxide isomer should be ready targets for astrochemical detection. The present quantum chemical study shows that not only are both molecular isomers strongly polar, they also have notable IR features and low isomerisation energies of 4.1 kcal mol-1 and 10.7 kcal mol-1, respectively. These three isomers along with two other C2O2H2 isomers have had their full set of fundamental vibrational frequencies and spectroscopic constants characterised herein. These isomers have previously been shown to occur in simulated astrophysical ices making them worthy targets of astronomical search. Furthermore, the hybrid quartic force field (QFF) approach utilized herein to produce the needed spectral data has a mean absolute percent error compared to the experimentally-available, gas phase fundamental vibrational frequencies of 0.6% and rotational constants to better than 0.1%. The hybrid QFF is defined from explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples level [CCSD(T)-F12b] including core electron correlation and a canonical CCSD(T) relativity correction for the harmonic (quadratic) terms in the QFF and simple CCSD(T)-F12b/cc-pVDZ energies for the cubic and quartic terms, the so-called "F12-TcCR+DZ QFF." This method is producing spectroscopically-accurate predictions for both fundamental vibrational frequencies and principal spectroscopic constants. Hence, the values computed in this work should be notably accurate and, hence, exceptionally useful to the spectroscopy and astrochemistry communities.

Quantum study of the rovibrational relaxation of HF by collision with 4He on a new potential energy surface

The HF molecule is considered the main reservoir of fluorine in the interstellar medium (ISM). Also, the interactions of this molecule with the most common atoms and molecules in the ISM have attracted great interest from the astrochemical community. Collisions between HF and helium have recently caused controversy following a study using a two-dimensional SAPT potential energy surface (PES) that exhibited large discrepancies with previous scattering calculations based on more recent ab initio potentials. To address this issue, our current work aims to develop the most precise three-dimensional PES for the HF+He system. We employ the size-consistent CCSD(T) method in conjunction with the aug-cc-pV6Z basis set. The main features of the new PES as well as the bound states of the He-HF complex are compared to the existing data. The new PES is then utilised to conduct close coupling calculations that demonstrate He-HF as a good instance of vibration-rotation near resonant energy transfer. The novel rate coefficients will be accessible via the BASECOL database, and the use of the new PES is advised when describing HF in helium droplets.

Transferable Water Potentials Using Equivariant Neural Networks

Machine learning interatomic potentials (MLIPs) have emerged as a technique that promises quantum theory accuracy for reduced cost. It has been proposed [J. Chem. Phys. 2023, 158, 084111] that MLIPs trained on solely liquid water data cannot accurately transfer to the vapor-liquid equilibrium while recovering the many-body decomposition (MBD) analysis of gas-phase water clusters. This suggests that MLIPs do not directly learn the physically correct interactions of water molecules, limiting transferability. In this work, we show that MLIPs using equivariant architecture and trained on 3200 liquid water structures reproduces liquid-phase water properties (e.g., density within 0.003 g/cm3 between 230 and 365 K), vapor-liquid equilibrium properties up to 550 K, the MBD analysis of gas-phase water cluster up to six-body interactions, and the relative energy and the vibrational density of states of ice phases. We show that potentials developed using equivariant MLIPs allow transferability for arbitrary phases of water that remain stable in nanosecond long simulations.

Quantum Chemistry and Astrochemistry: A Match Made in the Heavens

Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.

DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.

Towards understanding solvation effects on the conformational entropy of non-rigid molecules

The absolute molecular entropy is a fundamental quantity for the accurate description of thermodynamic properties. For non-rigid molecules, a substantial part of the entropy can be attributed to a conformational contribution. Systems and properties where this is relevant, e.g., protein-ligand binding affinities or pK_a values refer usually to the liquid phase. In this work, the influence of solvation on the conformational entropy is investigated. A recently introduced state-of-the-art and automated computational protocol for the computation of conformational entropies [Pracht et al., Chem. Sci., 2021, 12, 6551-6568.] is applied in combination with fast and accurate semiempirical quantum-chemical methods and implicit solvation models for a set of 25 commercially available drug molecules and five transition metal compounds. Computed gas-phase conformational entropies are compared with values obtained in implicit n-hexane and water. It is found that implicit solvation can have a substantial effect of several cal mol^-1 K^-1 on the entropy as a result of large conformational changes in the different phases. We conclude that for flexible molecules chemical accuracy for free energies in solution can only be achieved if solvation effects on the conformational ensemble are considered.

A shortcut to the thermodynamic limit for quantum many-body calculations of metals

Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely dense k-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors.

Pairing double hybrid functionals with a tailored basis set for an accurate thermochemistry of hydrocarbons

A collection of five challenging datasets, including noncovalent interactions, reaction barriers and electronic rearrangements of medium-sized hydrocarbons, has been selected to verify the robustness of double-hybrid functionals used in conjunction with the small DH-SVPD basis set, especially developed for noncovalent interactions. The analysis is completed by other, more standard functionals, for a total of 17 models, including also empirical corrections for dispersion. The obtained results show that the chemical accuracy threshold, that is an error lower than 1.0 kcal mol⁻¹, can be obtained by pairing the nonempirical PBE-QIDH functional with the DH-SVPD basis set, as well as by other semi-empirical functionals, such as DSD-PBEP86, using larger basis sets and empirical corrections. More in general, a significant improvement can be obtained using the DH-SVPD basis set with DHs, without resorting to any empirical corrections. This choice leads to a fast computational protocol that, avoiding any empirical potential, remains on a fully quantum ground.

The pairing of the PBE-QIDH double-hybrid functional with a tailored split-valence basis set leads to a fast computational protocol for the accurate evaluation of hydrocarbon thermochemistry, without resorting to any empirical correction.

Cerium Oxides without U: The Role of Many-Electron Correlation

Electron transfer with changing occupation in the 4f subshell poses a considerable challenge for quantitative predictions in quantum chemistry. Using the example of cerium oxide, we identify the main deficiencies of common parameter-dependent one-electron approaches, such as density functional theory (DFT) with a Hubbard correction, or hybrid functionals. As a response, we present the first benchmark of ab initio many-electron theory for electron transfer energies and lattice parameters under periodic boundary conditions. We show that the direct random phase approximation clearly outperforms all DFT variations. From this foundation, we, then, systematically improve even further. Periodic second-order Møller-Plesset perturbation theory meanwhile manages to recover standard hybrid functional values. Using these approaches to eliminate parameter bias allows for highly accurate benchmarks of strongly correlated materials, the reliable assessment of various density functionals, and functional fitting via machine-learning.

Rovibrational Quantum Chemical Treatment of Inorganic and Organometallic Astrochemicals

ConspectusOur two groups have both independently and collaboratively been pushing quantum-chemical techniques to produce highly accurate predictions of anharmonic vibrational frequencies and spectroscopic constants for molecules containing atoms outside of the typical upper p block. Methodologies employ composite approaches, relying on various levels of coupled cluster theory-most often at the singles, doubles, and perturbative triples level-and quartic force field constructions of the potential portion of the intramolecular Watson Hamiltonian. Such methods are known to perform well for organic species, and we have extended this to molecules containing atoms outside of this realm.One notable atom that has received much attention in this application is magnesium. Mg is the second-most-abundant element in the Earth's mantle, and while molecules containing this element are among the confirmed astrochemicals, its further atomic abundance in the galaxy implies that many more molecules (both purely inorganic and organometallic) containing element 12 exist in astrophysical regions in chemical sizes between those of atoms and dust-sized nanocrystals. Our approach discussed herein is producing quality benchmarks and predicting novel data for magnesium-bearing molecules.The story is similar for Al and Si, which are also notably abundant in both rocky bodies and the universe at large. While Na, Sc, and Cu may not be as abundant as Mg, Al, and Si, molecules containing Na and transition metals have also previously been reported to be detected beyond the Earth. Consequently, the need to produce spectral reference data for molecules containing such atoms is growing. While several experimental groups (including, notably, the groups in Arizona, Boston, and France/Spain) have clearly led the way in detection of inorganic/organometallic molecules in space, computational support and even rational design can provide novel avenues for the detection of molecules containing atoms not typically studied in most laboratories. The application of quantum chemistry to other elements beyond carbon and its cronies at the top right of the periodic table promises a better understanding of the observable universe. It will also provide novel and fundamental chemical insights pushing the "central science" into new molecular territory.

Basis set extrapolation in pair natural orbital theories

We present the results of a benchmark study of the effect of Pair Natural Orbital (PNO) truncation errors on the performance of basis set extrapolation. We find that reliable conclusions from the application of Helgaker's extrapolation method are only obtained when using tight PNO thresholds of at least 10^-7. The use of looser thresholds introduces a significant risk of observing a false basis set convergence and underestimating the residual basis set errors. We propose an alternative extrapolation approach based on the PNO truncation level that only requires a single basis set and show that it is a viable alternative to hierarchical basis set extrapolation methods.

A Molecular Candle Where Few Molecules Shine: HeHHe. Molecules 2020;25:molecules25092183. [PMID: 32392765 PMCID: PMC7249080 DOI: 10.3390/molecules25092183]

HeHHe+ is the only potential molecule comprised of atoms present in the early universe that is also easily observable in the infrared. This molecule has been known to exist in mass spectrometry experiments for nearly half-a-century and is likely present, but as-of-yet unconfirmed, in cold plasmas. There can exist only a handful of plausible primordial molecules in the epochs before metals (elements with nuclei heavier than 4He as astronomers call them) were synthesized in the universe, and most of these are both rotationally and vibrationally dark. The current work brings HeHHe+ into the discussion as a possible (and potentially only) molecular candle for probing high-z and any metal-deprived regions due to its exceptionally bright infrared feature previously predicted to lie at 7.43 μm. Furthermore, the present study provides new insights into its possible formation mechanisms as well as marked stability, along with the decisive role of anharmonic zero-point energies. A new entrance pathway is proposed through the triplet state (3B1) of the He2H+ molecule complexed with a hydrogen atom and a subsequent 10.90 eV charge transfer/photon emission into the linear and vibrationally-bright 1Σg+ HeHHe+ form.

Structure, stability, infrared spectra, and bonding of OHm(H2O)7 (m = 0, ±1) clusters: ab initio study combining the particle swarm optimization algorithm

The various structural candidates of anionic, neutral, and cationic water clusters OHm(H2O)7 (m = 0, ±1) have been globally predicted by combining the particle swarm optimization method and quantum chemical calculations. Geometry optimization and vibrational analysis for the optimal structures were performed with the MP2/aug-cc-pVDZ method, and the energy profile was further refined at the CCSD(T)/CBS level. Special attention was paid to the relationships between configurations and energies, particularly the first solvation shell coordination number of OH- and OH. For OH-(H2O)7, OH(H2O)7, and OH+(H2O)7 clusters, the most stable species at room temperature are predicted to be the tetra-solvated multi-ring structure A6, the tri-solvated hemibond cage structure N1, and the single five-membered ring structure C2, respectively. The temperature effects on the stability of these three systems were also explored via Gibbs free energies. Furthermore, for the OH-(H2O)7 clusters, the assignments of vibrational transitions in the OH stretching region are in good agreement with the studies of small hydroxide ion-water clusters, and the IR spectra of two isomers (tetra-solvated multi-ring A6 and penta-solvated cage A3) may match future experimental observation well. By topological analysis and reduced density gradient analysis, the structural characteristics and bonding strengths of the studied clusters were investigated. This work indicates the excellent performance of the PSO search algorithm and CALYPSO on water clusters, and may further provide extensive insights into the chemical behavior such as the transport mechanism of OH- ions and OH radicals in the aqueous phase.

Accurate Quantum Chemical Free Energies at Affordable Cost

Free energy sampling methods allow studying the full dynamics of activated processes. Unfortunately, the affordable accuracy of the potential describing the energy and forces of the system is usually rather low. Here we introduce a new method that by combining metadynamics and free energy perturbation allows calculating accurate quantum chemical free energies for chemical reactions. To prove the effectiveness of this new approach we study the S_N2 reaction of CH₃F + Cl^- → CH₃Cl + F^- in vacuo and solvated by water. Comparisons are made with harmonic transition-state theory to show how this method could provide accurate equilibrium and rate constants for complex systems.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol^-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.

Identifying Molecular Structural Aromaticity for Hydrocarbon Classification

Determination of aromaticity in hydrocarbons may be as simple as determining the average bond length for the molecule of interest. This would greatly assist in classifying the nature of hydrocarbon chemistry, especially for large molecules such as polycyclic aromatic hydrocarbons (PAHs) where today's aromatic classification methods are prohibitively expensive. The average C-C bond lengths for a test set of known aromatic, antiaromatic, and aliphatic cyclic hydrocarbons are computed here, and they show strong delineating patterns for the structural discernment of these aromaticity classifications. Aromatic molecules have average C-C bond lengths of 1.41 Å or less with the largest molecules, PAHs, having the longest average C-C bond lengths; aliphatic species have such lengths of 1.50 Å or more; and antiaromatic species fall between the two. Consequently, a first-order guess as to the aromaticity of a system may simply arise from its geometry. Although this prediction will likely have exceptions, such simple screening can easily classify most cases, and more advanced techniques can be brought to bear on the cases that lie in the boundaries. Benchmarks for hydrocarbons are provided here, but other classes of molecular structural aromaticity likely will have to be defined on an ad hoc basis.

Straightening the Hierarchical Staircase for Basis Set Extrapolations: A Low-Cost Approach to High-Accuracy Computational Chemistry

Because the one-electron basis set limit is difficult to reach in correlated post-Hartree-Fock ab initio calculations, the low-cost route of using methods that extrapolate to the estimated basis set limit attracts immediate interest. The situation is somewhat more satisfactory at the Hartree-Fock level because numerical calculation of the energy is often affordable at nearly converged basis set levels. Still, extrapolation schemes for the Hartree-Fock energy are addressed here, although the focus is on the more slowly convergent and computationally demanding correlation energy. Because they are frequently based on the gold-standard coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)], correlated calculations are often affordable only with the smallest basis sets, and hence single-level extrapolations from one raw energy could attain maximum usefulness. This possibility is examined. Whenever possible, this review uses raw data from second-order Møller-Plesset perturbation theory, as well as CCSD, CCSD(T), and multireference configuration interaction methods. Inescapably, the emphasis is on work done by the author's research group. Certain issues in need of further research or review are pinpointed.

Characterizing Chemical Similarity with Vibrational Spectroscopy: New Insights into the Substituent Effects in Monosubstituted Benzenes

A novel approach is presented to assess chemical similarity based the local vibrational mode analysis developed by Konkoli and Cremer. The local mode frequency shifts are introduced as similarity descriptors that are sensitive to any electronic structure change. In this work, 59 different monosubstituted benzenes are compared. For a subset of 43 compounds, for which experimental data was available, the ortho-/para- and meta-directing effect in electrophilic aromatic substitution reactions could be correctly reproduced, proving the robustness of the new similarity index. For the remaining 16 compounds, the directing effect was predicted. The new approach is broadly applicable to all compounds for which either experimental or calculated vibrational frequency information is available.

Towards completing the cyclopropenylidene cycle: rovibrational analysis of cyclic N3+, CNN, HCNN+, and CNC. Phys Chem Chem Phys 2017;19:22860-22869. [PMID: 28812071 DOI: 10.1039/c7cp04257d]

The simple aromatic hydrocarbon, cyclopropenylidene (c-C₃H₂), is a known, naturally-occurring molecule. The question remains as to whether its isoelectronic, cyclic, fellow aromatics of c-N₃⁺, c-CNN, HCNN⁺, and c-CNC^- are as well. Each of these are exciting objects for observation of Titan, and the rotational constants and vibrational frequencies produced here will allow for remote sensing of Titan's atmosphere or other astrophysical or terrestrial sources. None of these four aromatic species are vibrationally strong absorbers/emitters, but the two ions, HCNN⁺ and c-CNC^-, have dipole moments of greater than 3 D and 1 D, respectively, making them good targets for rotational spectroscopic observation. Each of these molecules is shown here to exhibit its own, unique vibrational properties, but the general trends put the vibrational behavior for corresponding fundamental modes within close ranges of one another, even producing nearly the same heavy atom, symmetric stretching frequencies for HCNN⁺ and c-C₃H₂ at 1600 cm^-1. The c-N₃⁺ cation is confirmed to be fairly unstable and has almost no intensity in its ν₂ fundamental. Hence, it will likely remain difficult to characterize experimentally.

Symmetry breaking and spectral considerations of the surprisingly floppy c-C3H radical and the related dipole-bound excited state of c-C3H. J Chem Phys 2017;146:224303. [PMID: 29166048 DOI: 10.1063/1.4985095]

The C₃H radical is believed to be prevalent throughout the interstellar medium and may be involved in the formation of polycyclic aromatic hydrocarbons. C₃H exists as both a linear and a cyclic isomer. The C_2v cyclopropenylidenyl radical isomer was detected in the dark molecular cloud TMC-1, and the linear propenylidenyl radical isomer has been observed in various dark molecular clouds. Even though the c-C₃H radical has been classified rotationally, the vibrational frequencies of this seemingly important interstellar molecule have never been directly observed. Established, highly accurate quartic force field methodologies are employed here to compute useful geometrical data, spectroscopic constants, and vibrational frequencies. The computed rotational constants are consistent with the experimental results. Consequently, the three a₁ (ν₁, ν₂, and ν₃) and one b₁ (ν₆) anharmonic vibrational frequencies at 3117.7 cm^-1, 1564.3 cm^-1, 1198.5 cm^-1, and 826.7 cm^-1, respectively, are reliable predictions for these, as of yet unseen, observables. Unfortunately, the two b₂ fundamentals (ν₄ and ν₅) cannot be treated adequately in the current approach due to a flat and possible double-well potential described in detail herein. The dipole-bound excited state of the anion suffers from the same issues and may not even be bound. However, the trusted fundamental vibrational frequencies described for the neutral radical should not be affected by this deformity and are the first robustly produced for c-C₃H. The insights gained here will also be applicable to other structures containing three-membered bare and exposed carbon rings that are surprisingly floppy in nature.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Rate constants for H abstraction from benzo(a)pyrene and chrysene: a theoretical study

A theoretical study of H abstraction reactions from benzo[a]pyrene and chrysene shows differences in kinetic effectiveness of various radicals and a clear distinction between zigzag and armchair edges.

Formation of Potential Interstellar Noble Gas Molecules in Gas and Adsorbed Phases

The discovery of naturally occurring ArH⁺ in various regions of the interstellar medium has shown the need for more understanding of the reactions that lead to covalently bonded noble gas molecules. The test comes with trying to predict the formation of other small noble gas molecules. Many molecules have been observed in various interstellar environments, which possess the possibility of bonding with noble gases. This work explores how both argon and neon can form bonds to ligands made of these species through quantum chemical computations. Argon and neon are chosen as they are among the most abundant atoms in the universe but are more polarizable than the more common but smaller helium atom. Reactions leading to noble gas molecules are modeled in the gas phase as well as through the adsorbed phase by catalysis with a polycyclic aromatic hydrocarbon (PAH) surface. The adsorption energy of the neutral noble gas atoms to the surface increases as the size of the PAH also increases but this is still less than 10 kcal/mol. It is proposed and supported herein that an incoming molecule can bond with the noble gas atom adsorbed onto the PAH, form a stable structure, and allow the PAH to function as the leaving group. This work shows that the noble gas molecules ArCCH⁺, ArOH⁺, ArNH⁺, and NeCCH⁺ are not only stable minima on their respective potential energy surfaces but also can be formed in either the gas phase or through PAH adsorption with known or hypothesized interstellar molecules. Most notably, NeCCH⁺ does not appear to form in the gas phase but could be catalyzed on PAH surfaces. Hence, the interstellar detection of such molecules could also serve as a probe for the observation of interstellar PAHs.

Methylidyne-replaced boron nitride fullerenes and nanotubes: a wave function study

Replacement of nitrogen atoms in boron nitride fullerenes and nanotubes with methylidyne groups stabilizes the structure and provides bonding points.

The Hubbard dimer: a density functional case study of a many-body problem

This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of density functional theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green's function, demonstrating the 'failure' of the Kohn-Sham (KS) method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe ansatz local density approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for density functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H2.

The performance of low-cost commercial cloud computing as an alternative in computational chemistry

The growth of commercial cloud computing (CCC) as a viable means of computational infrastructure is largely unexplored for the purposes of quantum chemistry. In this work, the PSI4 suite of computational chemistry programs is installed on five different types of Amazon World Services CCC platforms. The performance for a set of electronically excited state single-point energies is compared between these CCC platforms and typical, "in-house" physical machines. Further considerations are made for the number of cores or virtual CPUs (vCPUs, for the CCC platforms), but no considerations are made for full parallelization of the program (even though parallelization of the BLAS library is implemented), complete high-performance computing cluster utilization, or steal time. Even with this most pessimistic view of the computations, CCC resources are shown to be more cost effective for significant numbers of typical quantum chemistry computations. Large numbers of large computations are still best utilized by more traditional means, but smaller-scale research may be more effectively undertaken through CCC services.

Quartic force fields for excited electronic states: rovibronic reference data for the 1 (2)A' and 1 (2)A″ states of the isoformyl radical, HOC

Quartic force fields (QFFs) have been shown to be an effective, accurate, and relatively compact means of computing rovibrational spectroscopic data for numerous molecules with numerous applications. However, excited states have been nearly excluded from the this approach since most accurate QFFs are based on the "gold standard" coupled cluster singles, doubles, and perturbative triples [CCSD(T)] method which is not readily extended to excited states. In this work, rovibronic spectroscopic data is provided for the isoformyl radical, a molecule of significance in combustion and astrochemistry, both through the traditional means of variational access to excited states with CCSD(T) and in the novel extension of QFFs routinely to treat electronically excited states through the standard coupled cluster excited state approach, equation of motion (EOM) CCSD. It is shown here that the new EOM-based QFF provides structural parameters and rotational constants that are quite close to those from a related CCSD(T)-based QFF for the 1 (2)A(″) excited state of HOC. The anharmonic vibrational frequency percent differences between the two QFFs are less than 0.4% for the O-H stretch, less than 1.9% for the C-O stretch, and around 3.0% for the bend. Even so, the pure excited state EOM-QFF anharmonic frequencies are still very good abinitio representations that may be applied to systems where electronically excited states are not variationally accessible. Additionally, rovibrational spectroscopic data is provided for the 1 (2)A(') ground state of HOC and for both the ground and excited state of DOC.

Electron delocalization in polyenes: a semiexperimental equilibrium structure for (3E)-1,3,5-hexatriene and theoretical structures for (3Z,5Z)-, (3E,5E)-, and (3E,5Z)-1,3,5,7-octatetraene

Electronic structure theory reveals that π-electron delocalization increases with the chain length in polyenes. To analyze quantitatively this effect a semiexperimental equilibrium structure has been determined for trans-hexatriene by the mixed estimation method. For this fit rotational constants for a number of carbon and hydrogen isotopologues as well as a high-level ab initio structure have been used. The accuracy is 0.001 Å for bond lengths and 0.1° for bond angles. For the three isomers of octatetraene, high-level ab initio calculations have given a comparably accurate structure. These structures have been used in comparison with the structure of s-trans-butadiene to show that "C═C" bonds increase in length and "C-C" bonds decrease in length as the polyene chain lengthens. These structural effects of π-electron delocalization increase toward the center of polyenes. Most likely, π-π conjugation in the molecules studied plays a large part in their planarity that, in turn, forces the hydrogen atoms of cis fragments in bay regions to be in a close contact. Their distance is indeed shorter than the sum of their van der Waals radii, and they seem to participate in a six-membered ring.