Search

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

Over the past half century, computational chemistry has evolved from a niche field to a ubiquitous pillar of modern chemical research. Driven by the increased demand for computational chemistry in research settings, the undergraduate curriculum has evolved alongside to ensure that students are well-equipped for modern research. Toward this end, many excellent computational chemistry exercises have been developed that aim to teach students what kinds of questions computational chemistry can answer and how to properly interpret the results. However, there has been far less attention given to the complexities of determining how reliable computational results are and how constraining computational scaling can be. We present an undergraduate lab exercise that uses ab initio methods to predict the combustion energy of methane. The exercise walks students through the process of benchmarking errors on a small system (methane), estimating the computational cost to perform the same analysis on a larger system (propane), and justifying an affordable yet accurate method for a hypothetical study of the larger system. Furthermore, students are introduced to other cost-saving measures like basis set extrapolation and additive corrections. The entire exercise is intentionally designed to require little technical knowledge of computational chemistry and to be flexibly grafted into a standard undergraduate curriculum. In order to ensure accessibility, the exercise utilizes the free open source software Psi4 (available on any operating system) and provides a detailed installation and use guide for completing this lab. This lab will provide students the understanding of how to properly judge, select, and justify different computational models where cost and accuracy compete, a highly desirable set of skills that generalize to any computational science.

Trinuclear and tetranuclear ruthenium carbonyls of the types Ru3(CO)n(NO)2, Ru3(N)(CO)n(NO), Ru3(N)2(CO)n, Ru3(N)(CO)n(NCO), Ru3(CO)n(NCO)(NO), Ru4(N)(CO)n(NO), Ru4(N)(CO)n(NCO), and Ru4(N)2(CO)n related to species observed experimentally in the chemistry of Ru3(CO)10(µ-NO)2 have been investigated using density functional theory. In all cases, the experimentally observed structures have been found to be low-energy structures. The low-energy trinuclear structures typically have a central strongly bent Ru–Ru–Ru chain with terminal CO groups and bridging nitrosyl, isocyanate, and/or nitride ligands across the end of the chain. The low-energy tetranuclear structures typically have a central Ru4N unit with terminal CO groups and a non-bonded pair of ruthenium atoms bridged by a nitrosyl or isocyanate group. [ABSTRACT FROM AUTHOR]

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima. [ABSTRACT FROM AUTHOR]

Psi4 is an ab initio electronic structure program providing methods such as Hartree-Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools. Several new features have been added with the aid of libraries providing easy access to techniques such as density fitting, Cholesky decomposition, and Laplace denominators. The build system has been completely rewritten to simplify interoperability with independent, reusable software components for quantum chemistry. Finally, a wide range of new theoretical methods and analyses have been added to the code base, including functional-group and open-shell symmetry adapted perturbation theory, density-fitted coupled cluster with frozen natural orbitals, orbital-optimized perturbation and coupled-cluster methods (e.g., OO-MP2 and OO-LCCD), density-fitted multiconfigurational self-consistent field, density cumulant functional theory, algebraic-diagrammatic construction excited states, improvements to the geometry optimizer, and the "X2C" approach to relativistic corrections, among many other improvements.

Preceding even the Hartree-Fock method, molecular integrals are the very foundation upon which quantum chemical molecular modeling depends. Discussions of molecular integrals are normally found only in advanced and technical texts or articles. The objective of the present article is to provide less experienced readers, or students in a physical/computational chemistry course, a thorough understanding of molecular integrals. Through a series of detailed Handouts, the student/reader can participate in the derivation of molecular integrals, and in turn implement them in computer code. Hartree-Fock theory is discussed in enough detail to motivate the molecular integrals and address such topics as the atomic orbital basis. An introduction to the programming language of choice, Python3, is provided, tailored toward developing the essential skills necessary for implementing molecular integrals. The article is intended to be useful not only to instructors of physical/computational chemistry, but also to any reader who has independently sought a primer on this elusive subject.

Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet Tl 2H 2 potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 M 2H 2 (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist. [ABSTRACT FROM AUTHOR]

Five structures of Ge2H2 and Ge2H2+ are investigated in this study. Optimized geometries at the CCSD(T)/cc-pwCVQZ-PP level of theory were obtained. Focal point analyses were performed on these optimized geometries to determine relative energies using the CCSD(T) method with polarized basis sets up to quintuple-zeta. Energy corrections include full T and pertubative (Q) coupled-cluster effects plus anharmonic corrections to the zero-point vibrational energy. Relative ordering in energy from lowest to highest of the five Ge2H2+ structures is butterfly, germylidene, monobridged, trans, then linear. In neutral Ge2H2, the monobridged structure lies lower in energy than the germylidene structure. Fundamental vibrational frequencies and IR intensities were computed for the minima at the CCSD(T)/cc-pwCVTZ-PP level of theory to compare with experimental research. Partial atomic charges and natural bonding orbital analyses indicated that the positive charge of Ge2H2+ is contained in the region of the Ge–Ge bond. [ABSTRACT FROM AUTHOR]

This paper investigates four of the constitutional isomers of Si2H2+, namely the butterfly, vinylidene‐like, monobridged, and trans structures. These isomer geometries were all studied using the CCSD(T) method with basis sets as large as cc‐pV5Z. Higher level methods CCSDT and CCSDT(Q) were used for final energetics. It is found that the butterfly isomer has the lowest energy, followed by vinylidene‐like and monobridged structures; the trans structure has the highest energy of the four. All structures were compared with their neutral counterparts. Partial charges are computed and it is found that all isomers share the positive charge between the silicon atoms. NBO analysis shows that the cation Si−Si bond order is reduced by 0.5 relative to the neutral versions of the same isomer in the case of the butterfly, vinylidene‐like, and monobridged; and reduced by 0.6 for the trans isomer. Vibrational frequencies, infrared intensities, and dipole moments are predicted to encourage the spectroscopic identification of Si2H2+. [ABSTRACT FROM AUTHOR]