Showing papers by "Martin Head-Gordon published in 2015"
TL;DR: A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided in this paper, covering approximately the last seven years, including developments in density functional theory 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.
Abstract: 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.
2,396 citations
TL;DR: The final density functional, B97M-V, is remarkably accurate for non-bonded interactions and very satisfactory for thermochemical quantities such as atomization energies, but inherits the demonstrable limitations of existing local density functionals for barrier heights.
Abstract: A meta-generalized gradient approximation density functional paired with the VV10 nonlocal correlation functional is presented. The functional form is selected from more than 1010 choices carved out of a functional space of almost 1040 possibilities. Raw data come from training a vast number of candidate functional forms on a comprehensive training set of 1095 data points and testing the resulting fits on a comprehensive primary test set of 1153 data points. Functional forms are ranked based on their ability to reproduce the data in both the training and primary test sets with minimum empiricism, and filtered based on a set of physical constraints and an often-overlooked condition of satisfactory numerical precision with medium-sized integration grids. The resulting optimal functional form has 4 linear exchange parameters, 4 linear same-spin correlation parameters, and 4 linear opposite-spin correlation parameters, for a total of 12 fitted parameters. The final density functional, B97M-V, is further assessed on a secondary test set of 212 data points, applied to several large systems including the coronene dimer and water clusters, tested for the accurate prediction of intramolecular and intermolecular geometries, verified to have a readily attainable basis set limit, and checked for grid sensitivity. Compared to existing density functionals, B97M-V is remarkably accurate for non-bonded interactions and very satisfactory for thermochemical quantities such as atomization energies, but inherits the demonstrable limitations of existing local density functionals for barrier heights.
299 citations
TL;DR: In this article, a new parameter set (P2) was derived by rescaling the previously reported characteristic energies of the Lennard-Jones potential in P1, and the accuracy of thermal correction for adsorption enthalpies determined by the rigid rotor-harmonic oscillator approximation (RRHO) was examined and shown to be improved by treating low-lying vibrational modes as free translational and rotational modes via a quasi-RRHO model.
Abstract: Quantum mechanics/molecular mechanics (QM/MM) simulations provide an efficient avenue for studying reactions catalyzed in zeolite systems; however, the accuracy of such calculations is highly dependent on the zeolite MM parameters used. Previously reported parameters (P1), which were chosen to minimize the root mean square (RMS) deviations of adsorption energies compared with full QM ωB97X-D/6-31+G** adsorption energies, are shown to overestimate binding energies compared with experimental values, particularly for larger substrates. To address this issue, a new parameter set (P2) is derived by rescaling the previously reported characteristic energies of the Lennard-Jones potential in P1. The accuracy of the thermal correction for adsorption enthalpies determined by the rigid rotor-harmonic oscillator approximation (RRHO) is examined and shown to be improved by treating low-lying vibrational modes as free translational and rotational modes via a quasi-RRHO model. With P2 and quasi-RRHO, adsorption energies...
96 citations
TL;DR: It is established that intermolecular distance, as measured by, e.g., the center-of-mass separation of two molecules, is the geometric parameter that deviates most profoundly among the various methods and makes it an ideal metric for the development and evaluation of electronic structure methods.
Abstract: Electronic structure approaches for calculating intermolecular interactions have traditionally been benchmarked almost exclusively on the basis of energy-centric metrics. Herein, we explore the idea of utilizing a metric related to geometry. On a diverse series of noncovalently interacting systems of different sizes, from the water dimer to the coronene dimer, we evaluate a variety of electronic structure approximations with respect to their abilities to reproduce coupled-cluster-level geometries. Specifically, we examine Hartree-Fock, second-order Moller-Plesset perturbation theory (MP2), attenuated MP2, scaled MP2, and a number of density functionals, many of which include empirical or nonempirical van der Waals dispersion corrections. We find a number of trends that transcend system size and interaction type. For instance, functionals incorporating VV10 nonlocal correlation tend to yield highly accurate geometries; ωB97X-V and B97M-V, in particular, stand out. We establish that intermolecular distance, as measured by, e.g., the center-of-mass separation of two molecules, is the geometric parameter that deviates most profoundly among the various methods. This property of the equilibrium intermolecular separation, coupled with its accessibility via a small series of well-defined single-point calculations, makes it an ideal metric for the development and evaluation of electronic structure methods.
87 citations
TL;DR: A fundamental difference between the polarization degrees of freedom in the fragment-blocked approaches and in constrained density schemes is proved, which proves that only the former are capable of properly prohibiting charge delocalization during polarization.
Abstract: The polarization energy in intermolecular interactions treated by self-consistent field electronic structure theory is often evaluated using a constraint that the atomic orbital (AO) to molecular orbital transformation is blocked by fragments. This approach is tied to AO basis sets, overestimates polarization energies in the overlapping regime, particularly in large AO basis sets, and lacks a useful complete basis set limit. These problems are addressed by the construction of polarization subspaces based on the responses of isolated fragments to weak electric fields. These subspaces are spanned by fragment electric-field response functions, which can capture effects up to the dipole (D), or quadrupole (DQ) level, or beyond. Schemes are presented for the creation of both non-orthogonal and orthogonal fragment subspaces, and the basis set convergence of the polarization energies computed using these spaces is assessed. Numerical calculations for the water dimer, water-Na(+), water-Mg(2+), water-F(-), and water-Cl(-) show that the non-orthogonal DQ model is very satisfactory, with small differences relative to the orthogonalized model. Additionally, we prove a fundamental difference between the polarization degrees of freedom in the fragment-blocked approaches and in constrained density schemes. Only the former are capable of properly prohibiting charge delocalization during polarization.
87 citations
TL;DR: In this paper, a series of cobalt complexes bearing redox-active pyrazines were designed and analyzed to evaluate the effects of positioning redox non-innocent ligands on hydrogen evolution catalysis.
Abstract: Mononuclear metalloenzymes in nature can function in cooperation with precisely positioned redox-active organic cofactors in order to carry out multielectron catalysis. Inspired by the finely tuned redox management of these bioinorganic systems, we present the design, synthesis, and experimental and theoretical characterization of a homologous series of cobalt complexes bearing redox-active pyrazines. These donor moieties are locked into key positions within a pentadentate ligand scaffold in order to evaluate the effects of positioning redox non-innocent ligands on hydrogen evolution catalysis. Both metal- and ligand-centered redox features are observed in organic as well as aqueous solutions over a range of pH values, and comparison with analogs bearing redox-inactive zinc(ii) allows for assignments of ligand-based redox events. Varying the geometric placement of redox non-innocent pyrazine donors on isostructural pentadentate ligand platforms results in marked effects on observed cobalt-catalyzed proton reduction activity. Electrocatalytic hydrogen evolution from weak acids in acetonitrile solution, under diffusion-limited conditions, reveals that the pyrazine donor of axial isomer 1-Co behaves as an unproductive electron sink, resulting in high overpotentials for proton reduction, whereas the equatorial pyrazine isomer complex 2-Co is significantly more active for hydrogen generation at lower voltages. Addition of a second equatorial pyrazine in complex 3-Co further minimizes overpotentials required for catalysis. The equatorial derivative 2-Co is also superior to its axial 1-Co congener for electrocatalytic and visible-light photocatalytic hydrogen generation in biologically relevant, neutral pH aqueous media. Density functional theory calculations (B3LYP-D2) indicate that the first reduction of catalyst isomers 1-Co, 2-Co, and 3-Co is largely metal-centered while the second reduction occurs at pyrazine. Taken together, the data establish that proper positioning of non-innocent pyrazine ligands on a single cobalt center is indeed critical for promoting efficient hydrogen catalysis in aqueous media, akin to optimally positioned redox-active cofactors in metalloenzymes. In a broader sense, these findings highlight the significance of electronic structure considerations in the design of effective electron-hole reservoirs for multielectron transformations.
81 citations
TL;DR: In this paper, a scaling relation between the binding energy of the catalyst to hydrogen and that to COOH, a key intermediate in the reduction of CO2 to CO; however, the M−H bond is stronger than the M-COOH bond, driving the reaction toward the competitive hydrogen evolution reaction (HER).
Abstract: Density functional theory (DFT) calculations are performed to investigate the energetics of the CO2 electrochemical reduction on metal (M) porphyrin-like motifs incorporated into graphene layers. The objective is to develop strategies that enhance CO2 reduction while suppressing the competitive hydrogen evolution reaction (HER). We find that there exists a scaling relation between the binding energy of the catalyst to hydrogen and that to COOH, a key intermediate in the reduction of CO2 to CO; however, the M–H bond is stronger than the M–COOH bond, driving the reaction toward the HER rather than the reduction of CO2 to CO. This scaling relation holds even with axial ligation to the metal cation coordinated to the porphyrin ring. When 4f lanthanide or 5f actinide elements are used as the reactive center, the scaling relation still holds but the M–COOH bond is stronger than the M–H bond, and the reaction favors the reduction of CO2 to CO. By contrast, there is no scaling relation between the binding energy ...
77 citations
TL;DR: In this article, the structure and performance of Pt-Ir catalysts for ethane and propane dehydrogenation was examined and compared with the Pt and Pt-Sn catalysts, and it was shown that at high conversions, Ir exhibited lower initial activity than Pt3Sn/Mg(Al)O but greater stability to coke deposition.
Abstract: Increased demand for light alkenes has motivated research on the catalytic dehydrogenation of the light alkanes and on understanding the role of catalyst composition on the activity, selectivity, and stability of Pt-based catalysts used for this purpose. The present study examines the structure and performance of Pt–Ir catalysts for ethane and propane dehydrogenation, and compares them with the performance of Pt and Pt–Sn catalysts. Nanoparticles of Pt, PtSn, and PtIr were prepared in a colloidal suspension and then dispersed onto calcined hydrotalcite (Mg(Al)O). After characterization to confirm formation of a bimetallic alloy, it was observed that at high conversions, Pt3Ir/Mg(Al)O exhibited lower initial activity than Pt3Sn/Mg(Al)O but greater stability to coke deposition. Intrinsic rate measurements at low feed residence time revealed the following trend in activity: Pt3Sn > Pt3Ir > Pt. DFT calculations carried out on tetrahedral clusters (Pt4, Pt3Ir, Pt3Sn) reveals that this trend in activity can be replicated and Ir is capable of alkane activation, a trait unique to this bimetallic system.
63 citations
TL;DR: It is found that a combination of weakly acidic silanol/aluminol with secondary amine groups can mimic proline catalysts and are more effective in catalyzing the selective dimerization reaction than the combination of amines with organic acids.
Abstract: A highly efficient water-tolerant, solid-base catalyst for the self-condensation of biomass-derived methyl ketones to jet-diesel fuel precursors was developed by grafting site-isolated secondary amines on silica-alumina supports. It is shown that apart from the nature and density of amine groups and the spatial separation of the acidic and basic sites, the acidity of the support material plays a critical role in defining the catalytic activity. It is also found that a combination of weakly acidic silanol/aluminol with secondary amine groups can mimic proline catalysts and are more effective in catalyzing the selective dimerization reaction than the combination of amines with organic acids. In situ FTIR measurements demonstrate that acidic groups activate methyl ketones through their carbonyl groups leading to a favorable CC bond formation step involving an enamine intermediate. DFT analysis of the reaction pathway confirms that CC bond formation is the rate-limiting step.
59 citations
TL;DR: A key benefit is that it is possible to obtain completely spin-pure exchange coupling constants with inexpensive ab initio calculations, and this procedure relies on neither spin-contaminated states nor nonunique spin-projection formulas.
Abstract: We provide a simple procedure for using inexpensive ab initio calculations to compute exchange coupling constants, J(AB), for multiradical molecules containing both an arbitrary number of radical sites and an arbitrary number of unpaired electrons For a system comprised of 2M unpaired electrons, one needs only to compute states having the Ŝ(z) quantum number M - 1 Conveniently, these are precisely the states that are accessed by the family of single spin-flip methods Building an effective Hamiltonian with these states allows one to extract all of the J(AB) constants in the molecule Unlike approaches based on density functional theory, this procedure relies on neither spin-contaminated states nor nonunique spin-projection formulas A key benefit is that it is possible to obtain completely spin-pure exchange coupling constants with inexpensive ab initio calculations A couple of examples are provided to illustrate the approach, including a 4-nickel cubane complex and a 6-chromium horseshoe complex with 18 entangled electrons
56 citations
TL;DR: A compressed exchange matrix is computed that yields the exact exchange energy, gradient, and direct inversion of the iterative subspace (DIIS) error vector that is well-behaved, that the accuracy of computed energetics is excellent, and that significant speedups can be obtained over existing integral-direct and RI-K methods.
Abstract: Construction of the exact exchange matrix, K, is typically the rate-determining step in hybrid density functional theory, and therefore, new approaches with increased efficiency are highly desirable. We present a framework with potential for greatly improved efficiency by computing a compressed exchange matrix that yields the exact exchange energy, gradient, and direct inversion of the iterative subspace (DIIS) error vector. The compressed exchange matrix is constructed with one index in the compact molecular orbital basis and the other index in the full atomic orbital basis. To illustrate the advantages, we present a practical algorithm that uses this framework in conjunction with the resolution of the identity (RI) approximation. We demonstrate that convergence using this method, referred to hereafter as occupied orbital RI-K (occ-RI-K), in combination with the DIIS algorithm is well-behaved, that the accuracy of computed energetics is excellent (identical to conventional RI-K), and that significant speedups can be obtained over existing integral-direct and RI-K methods. For a 4400 basis function C68H22 hydrogen-terminated graphene fragment, our algorithm yields a 14 × speedup over the conventional algorithm and a speedup of 3.3 × over RI-K.
TL;DR: An efficient new molecular orbital (MO) basis algorithm is reported implementing the pair atomic resolution of the identity approximation (PARI) to evaluate the exact exchange contribution (K) to self-consistent field methods, such as hybrid and range-separated hybrid density functionals.
Abstract: An efficient new molecular orbital (MO) basis algorithm is reported implementing the pair atomic resolution of the identity approximation (PARI) to evaluate the exact exchange contribution (K) to self-consistent field methods, such as hybrid and range-separated hybrid density functionals. The PARI approximation, in which atomic orbital (AO) basis function pairs are expanded using auxiliary basis functions centered only on their two respective atoms, was recently investigated by Merlot et al. [J. Comput. Chem. 2013, 34, 1486]. Our algorithm is significantly faster than quartic scaling RI-K, with an asymptotic exchange speedup for hybrid functionals of (1 + X/N), where N and X are the AO and auxiliary basis dimensions. The asymptotic speedup is 2 + 2X/N for range separated hybrids such as CAM-B3LYP, ωB97X-D, and ωB97X-V which include short- and long-range exact exchange. The observed speedup for exchange in ωB97X-V for a C68 graphene fragment in the cc-pVTZ basis is 3.4 relative to RI-K. Like conventional R...
TL;DR: An energy decomposition analysis (EDA) of intermolecular interactions is proposed for second-order Møller-Plesset perturbation theory (MP2) based on absolutely localized molecular orbitals (ALMOs) as an extension to a previous ALMO-based EDA for self-consistent field methods.
Abstract: An energy decomposition analysis (EDA) of intermolecular interactions is proposed for second-order Moller-Plesset perturbation theory (MP2) based on absolutely localized molecular orbitals (ALMOs), as an extension to a previous ALMO-based EDA for self-consistent field methods. It decomposes the canonical MP2 binding energy by dividing the double excitations that contribute to the MP2 wave function into classes based on how the excitations involve different molecules. The MP2 contribution to the binding energy is decomposed into four components: frozen interaction, polarization, charge transfer, and dispersion. Charge transfer is defined by excitations that change the number of electrons on a molecule, dispersion by intermolecular excitations that do not transfer charge, and polarization and frozen interactions by intra-molecular excitations. The final two are separated by evaluations of the frozen, isolated wave functions in the presence of the other molecules, with adjustments for orbital response. Unlike previous EDAs for electron correlation methods, this one includes components for the electrostatics, which is vital as adjustment to the electrostatic behavior of the system is in some cases the dominant effect of the treatment of electron correlation. The proposed EDA is then applied to a variety of different systems to demonstrate that all proposed components behave correctly. This includes systems with one molecule and an external electric perturbation to test the separation between polarization and frozen interactions and various bimolecular systems in the equilibrium range and beyond to test the rest of the EDA. We find that it performs well on these tests. We then apply the EDA to a halogen bonded system to investigate the nature of the halogen bond.
TL;DR: The results indicate that more sophisticated methods using complex basis functions are worth pursuing in the search for accurate and computationally feasible methods for computing resonance energies in molecular systems.
Abstract: The method of complex basis functions for computing positions and widths of molecular resonances is revisited. An open-ended and efficient implementation is described. The basis set requirements of the complex basis are investigated within the computationally inexpensive static-exchange approximation, and the results of this investigation lead to a hierarchy of basis sets for complex basis function calculations on small molecules. These basis sets are then applied in static-exchange calculations on some larger molecules with multiple low energy shape resonances: carbon tetrafluoride, benzene, pyridine, pyrimidine, pyrazine, and s-triazine. The results indicate that more sophisticated methods using complex basis functions are worth pursuing in the search for accurate and computationally feasible methods for computing resonance energies in molecular systems.
TL;DR: In situ FT-IR measurements and electronic structure calculations are reported for the reduction of CO2 catalyzed by the macrocyclic complex [CoIIN4H]2+ (N4H = 2,12-dimethyl-3,7,11,17-tetraazabicyclo-[11.3.1]-heptadeca-1(17),2, 11,13,15-pentaene) as mentioned in this paper.
Abstract: In situ FT-IR measurements and electronic structure calculations are reported for the reduction of CO2 catalyzed by the macrocyclic complex [CoIIN4H]2+ (N4H = 2,12-dimethyl-3,7,11,17-tetraazabicyclo-[11.3.1]-heptadeca-1(17),2,11,13,15-pentaene). Beginning from the [CoIIN4H]2+ resting state of the complex in wet acetonitrile solution, two different visible light sensitizers with substantially different reducing power are employed to access reduced states. Accessing reduced states of the complex with a [Ru(bpy)3]2+ sensitizer yields an infrared band at 1670 cm–1 attributed to carboxylate, which is also observed for an authentic sample of the one-electron reduced complex [CoN4H(MeCN)]+ in CO2 saturated acetonitrile solution. The results are interpreted based on calculations using the pure BP86 functional that correctly reproduces experimental geometries. Continuum solvation effects are also included. The calculations show that Co is reduced to CoI in the first reduction, which is consistent with experimental...
TL;DR: An orbital pairing theorem is introduced, with which it is demonstrated in genuine molecular systems that cRHF is capable of removing certain potential energy surface singularities that plague real-valued RHF and related methods.
Abstract: Restricted Hartree Fock using complex-valued orbitals (cRHF) is studied. We introduce an orbital pairing theorem, with which we obtain a concise connection between cRHF and real-valued RHF, and use it to uncover the close relationship between cRHF, unrestricted Hartree Fock, and generalized valence bond perfect pairing. This enables an intuition for cRHF, contrasting with the generally unintuitive nature of complex orbitals. We also describe an efficient computer implementation of cRHF and its corresponding stability analysis. By applying cRHF to the Be + H2 insertion reaction, a Woodward-Hoffmann violating reaction, and a symmetry-driven conical intersection, we demonstrate in genuine molecular systems that cRHF is capable of removing certain potential energy surface singularities that plague real-valued RHF and related methods. This complements earlier work that showed this capability in a model system. We also describe how cRHF is the preferred RHF method for certain radicaloid systems like singlet oxygen and antiaromatic molecules. For singlet O2, we show that standard methods fail even at the equilibrium geometry. An implication of this work is that, regardless of their individual efficacies, cRHF solutions to the HF equations are fairly commonplace.
TL;DR: In this article, the authors formulate the Hessian-vector product within the iterative diagonalization procedure as a finite difference of the electronic gradient with respect to orbital perturbations in the direction of the vector.
Abstract: Wavefunction stability analysis is commonly applied to converged self-consistent field (SCF) solutions to verify whether the electronic energy is a local minimum with respect to second-order variations in the orbitals. By iterative diagonalisation, the procedure calculates the lowest eigenvalue of the stability matrix or electronic Hessian. However, analytical expressions for the electronic Hessian are unavailable for most advanced post-Hartree–Fock (HF) wave function methods and even some Kohn–Sham (KS) density functionals. To address such cases, we formulate the Hessian-vector product within the iterative diagonalisation procedure as a finite difference of the electronic gradient with respect to orbital perturbations in the direction of the vector. As a model application, following the lowest eigenvalue of the orbital-optimised second-order Moller–Plesset perturbation theory (OOMP2) Hessian during H2 dissociation reveals the surprising stability of the spin-restricted solution at all separations, with a...
TL;DR: Spin-unrestricted NH-SCF is used to compute a complex potential energy surface for the carbon monoxide anion which correctly describes dissociation.
Abstract: This work describes the implementation and applications of non-Hermitian self-consistent field (NH-SCF) theory with complex basis functions for the ab initio computation of positions and widths of shape resonances in molecules. We utilize both the restricted open-shell and the previously unexplored spin-unrestricted variants to compute Siegert energies of several anionic shape resonances in small diatomic and polyatomic molecules including carbon tetrafluoride which has been the subject of several recent experimental studies. The computation of general molecular properties from a non-Hermitian wavefunction is discussed, and a density-based analysis is applied to the (2)B1 shape resonance in formaldehyde. Spin-unrestricted NH-SCF is used to compute a complex potential energy surface for the carbon monoxide anion which correctly describes dissociation.
TL;DR: In this paper, the authors applied quantum mechanics/molecular mechanics (QM/MM) models to investigate the adsorption and cracking of n-hexane on ZSM-5 and Faujasite zeolite structures.
Abstract: Quantum mechanics/molecular mechanics (QM/MM) models are applied to investigate the adsorption and cracking of n-hexane on ZSM-5 and Faujasite zeolite structures. These simulations account for the long-range electrostatic and midrange van-der-Waals interactions in the zeolite and provide energy barriers that are close to experimental data. The active acidic site was modeled by dispersion corrected density functional theory (DFT, ω B97X-D6-311/G*). The long-range interactions were calculated by molecular mechanics (MM). The adsorbed molecules under investigation are characterized by their thermodynamic properties (adsorption energy and enthalpy). The influence of the zeolite type on the thermodynamic properties is also pointed out. The results reveal that the kinetics of cracking is insensitive to differences in acid strengths. The thermodynamic data obtained are mainly influenced by the adsorption energy of n-hexane on ZSM-5 and/or Faujasite (Y) structures. The pore sizes of the zeolite types can lead to ...
TL;DR: DFT analysis of the reaction pathway indicates that, for weak Brønsted acid groups, such as Si--OH, the rate-limiting step is C--C bond formation, whereas for stronger Brénstedacid groups,such as Ti and Al, hydrolysis of iminium species produced upon C-- C bond formation is the Rate-Limiting step.
Abstract: The catalytic activity of secondary amines supported on mesoporous silica for the self-condensation of n-butanal to 2-ethylhexenal can be altered significantly by controlling the Bronsted acidity of M--OH species present on the surface of the support In this study, M--OH (M=Sn, Zr, Ti, and Al) groups were doped onto the surface of SBA-15, a mesoporous silica, prior to grafting secondary propyl amine groups on to the support surface The catalytic activity was found to depend critically on the synthesis procedure, the nature and amount of metal species introduced and the spatial separation between the acidic sites and amine groups DFT analysis of the reaction pathway indicates that, for weak Bronsted acid groups, such as Si--OH, the rate-limiting step is C--C bond formation, whereas for stronger Bronsted acid groups, such as Ti and Al, hydrolysis of iminium species produced upon C--C bond formation is the rate-limiting step Theoretical analysis shows further that the apparent activation energy decreases with increasing Bronsted acidity of the M--OH groups, consistent with experimental observation
TL;DR: This work introduces a necessary and sufficient condition for an arbitrary wavefunction to be collinear, i.e., its spin is quantized along some axis, and adapts the procedure for Generalized Hartree Fock, and uses it to study two dissociation pathways in CO2.
Abstract: We introduce a necessary and sufficient condition for an arbitrary wavefunction to be collinear, i.e., its spin is quantized along some axis. It may be used to obtain a cheap and simple computational procedure to test for collinearity in electronic structure theory calculations. We adapt the procedure for Generalized Hartree Fock (GHF), and use it to study two dissociation pathways in CO2. For these dissociation processes, the GHF wave functions transform from low-spin Unrestricted Hartree Fock (UHF) type states to noncollinear GHF states and on to high-spin UHF type states, phenomena that are succinctly illustrated by the constituents of the collinearity test. This complements earlier GHF work on this molecule.
TL;DR: Insight into the origin of observations made in plasma experiments mimicking interstellar and circumstellar conditions is sought and structural information at a high level of electronic structure theory, as well as information on the vibrational and electronic spectra of the low-energy isomers is reported.
Abstract: We seek insight into the origin of observations made in plasma experiments mimicking interstellar and circumstellar conditions. To this end theory is applied to the low-energy isomers of C4H2+, C6H2+ and C6H4+ and their formation paths from acetylene and its fragments. Ab initio molecular dynamics trajectories are performed to explore which isomers are readily accessible from acetylene and its ion fragments. Structural information at a high level of electronic structure theory [CCSD(T)/cc-pVTZ], as well as information on the vibrational [UMP2] and electronic spectra [ωB97X] of the low-energy isomers is reported.
TL;DR: A local excited-state method, based on the configuration interaction singles (CIS) wave function, for large atomic and molecular clusters, which results in formal scaling with the third power of the system size for computing the full spectrum of ALMO-CIS excited states.
Abstract: We develop a local excited-state method, based on the configuration interaction singles (CIS) wave function, for large atomic and molecular clusters. This method exploits the properties of absolutely localized molecular orbitals (ALMOs), which strictly limits the total number of excitations, and results in formal scaling with the third power of the system size for computing the full spectrum of ALMO-CIS excited states. The derivation of the equations and design of the algorithm are discussed in detail, with particular emphasis on the computational scaling. Clusters containing ∼500 atoms were used in evaluating the scaling, which agrees with the theoretical predictions, and the accuracy of the method is evaluated with respect to standard CIS. A pioneering application to the size dependence of the helium cluster spectrum is also presented for clusters of 25-231 atoms, the largest of which results in the computation of 2310 excited states per sampled cluster geometry.
TL;DR: A spin-component scaled version of attenuated MP2 is produced using the W4-11 database as a supplemental thermochemistry training set, and the resulting method reproduces the quality of MP2-V(terfc, aTZ) for noncovalent interactions and exceeds the performance of SCS-MP2/aTZ for thermochemistry.
Abstract: Attenuated second order Moller-Plesset theory (MP2) captures intermolecular binding energies at equilibrium geometries with high fidelity with respect to reference methods, yet must fail to reproduce dispersion energies at stretched geometries due to the removal of fully long-range dispersion. For this problem to be ameliorated, long-range correction using the VV10 van der Waals density functional is added to attenuated MP2, capturing short-range correlation with attenuated MP2 and long-range dispersion with VV10. Attenuated MP2 with long-range VV10 dispersion in the aug-cc-pVTZ (aTZ) basis set, MP2-V(terfc, aTZ), is parametrized for noncovalent interactions using the S66 database and tested on a variety of noncovalent databases, describing potential energy surfaces and equilibrium binding energies equally well. Further, a spin-component scaled (SCS) version, SCS-MP2-V(2terfc, aTZ), is produced using the W4-11 database as a supplemental thermochemistry training set, and the resulting method reproduces the quality of MP2-V(terfc, aTZ) for noncovalent interactions and exceeds the performance of SCS-MP2/aTZ for thermochemistry.
TL;DR: The synthesis of the 2-methyl-1-silacycloprop-2-enylidene molecule in the bimolecular gas-phase reaction of silylidyne with allene enriches the understanding toward the formation of organosilicon species in the gas phase of the interstellar medium in particular via exoergic reactions of no entrance barrier.
Abstract: The crossed molecular beam reactions of the ground-state silylidyne radical (SiH; X(2)Π) with allene (H2CCCH2; X(1)A1) and D4-allene (D2CCCD2; X(1)A1) were carried out at collision energies of 30 kJ mol(-1). Electronic structure calculations propose that the reaction of silylidyne with allene has no entrance barrier and is initiated by silylidyne addition to the π electron density of allene either to one carbon atom (C1/C2) or to both carbon atoms simultaneously via indirect (complex forming) reaction dynamics. The initially formed addition complexes isomerize via two distinct reaction pathways, both leading eventually to a cyclic SiC3H5 intermediate. The latter decomposes through a loose exit transition state via an atomic hydrogen loss perpendicularly to the plane of the decomposing complex (sideways scattering) in an overall exoergic reaction (experimentally: -19 ± 13 kJ mol(-1); computationally: -5 ± 3 kJ mol(-1)). This hydrogen loss yields the hitherto elusive 2-methyl-1-silacycloprop-2-enylidene molecule (c-SiC3H4), which can be derived from the closed-shell cyclopropenylidene molecule (c-C3H2) by replacing a hydrogen atom with a methyl group and the carbene carbon atom by the isovalent silicon atom. The synthesis of the 2-methyl-1-silacycloprop-2-enylidene molecule in the bimolecular gas-phase reaction of silylidyne with allene enriches our understanding toward the formation of organosilicon species in the gas phase of the interstellar medium in particular via exoergic reactions of no entrance barrier. This facile route to 2-methyl-1-silacycloprop-2-enylidene via a silylidyne radical reaction with allene opens up a versatile approach to form hitherto poorly characterized silicon-bearing species in extraterrestrial environments; this reaction class might represent the missing link, leading from silicon-bearing radicals via organosilicon chemistry eventually to silicon-carbon-rich interstellar grains even in cold molecular clouds where temperatures are as low as 10 K.
TL;DR: Author(s): Zuev, Dmitry; Jagau, Thomas-C; Bravaya, Ksenia B; Epifanovsky, Evgeny; Shao, Yihan; Sundstrom, Eric; Head-Gordon, Martin; Krylov, Anna I
Abstract: Author(s): Zuev, Dmitry; Jagau, Thomas-C; Bravaya, Ksenia B; Epifanovsky, Evgeny; Shao, Yihan; Sundstrom, Eric; Head-Gordon, Martin; Krylov, Anna I
TL;DR: The proposed intermolecular similarity-transformed perturbation theory represents a frozen-orbital variant of the "(2)"-type theories shown to be competitive with CCSD(T) and of similar cost if all terms are retained.
Abstract: Your correspondents develop and apply fully nonorthogonal, local-reference perturbation theories describing non-covalent interactions. Our formulations are based on a Lowdin partitioning of the similarity-transformed Hamiltonian into a zeroth-order intramonomer piece (taking local CCSD solutions as its zeroth-order eigenfunction) plus a first-order piece coupling the fragments. If considerations are limited to a single molecule, the proposed intermolecular similarity-transformed perturbation theory represents a frozen-orbital variant of the “(2)”-type theories shown to be competitive with CCSD(T) and of similar cost if all terms are retained. Different restrictions on the zeroth- and first-order amplitudes are explored in the context of large-computation tractability and elucidation of non-local effects in the space of singles and doubles. To accurately approximate CCSD intermolecular interaction energies, a quadratically growing number of variables must be included at zeroth-order.