Accurate Determination of the Reaction Course in HY2 ⇌ Y + YH (Y = O, S): Detailed Analysis of the Covalent- to Hydrogen-Bonding Transition

Cn (n=2-4): current status

The major aspects of the C₂, C₃ and C₄ elemental carbon clusters are surveyed. For C₂, a brief analysis of its current status is presented. Regarding C₃, the most recent results obtained in our group are reviewed with emphasis on modelling its potential energy surface which is particularly complicated due to the presence of multiple conical intersections. As for C₄, the most stable isomeric forms of both triplet and singlet spin states and their possible interconversion pathways are examined afresh by means of accurate ab initio calculations. The main strategies for modelling the ground triplet C₄ potential are also discussed. Starting from a truncated cluster expansion and a previously reported DMBE form for C₃, an approximate four-body term is calibrated from the ab initio energies. The final six-dimensional global DMBE form so obtained reproduces all known topographical aspects while providing an accurate description of the C₄ linear-rhombic isomerization pathway. It is therefore commended for both spectroscopic and reaction dynamics studies.This article is part of the theme issue 'Modern theoretical chemistry'.

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.

CBS extrapolation in electronic structure pushed to the end: a revival of minimal and sub-minimal basis sets

The complete basis set (CBS) limit is secluded in calculations of electronic structure, and hence CBS extrapolation draws immediate attention.

Assessing How Correlated Molecular Orbital Calculations Can Perform versus Kohn-Sham DFT: Barrier Heights/Isomerizations

To assess the title issue, 38 hydrogen transfer barrier heights and 38 non-hydrogen transfer barrier heights/isomerizations extracted from extensive databases have been considered, in addition to 4 2 p-isomerization reactions and 6 others for large organic molecules. All Kohn-Sham DFT calculations have employed the popular M06-2X functional, whereas the correlated molecular orbital (MO)-based ones are from single-reference MP2 and CCSD(T) methods. They have all utilized the same basis sets, with raw MO energies subsequently extrapolated to the complete basis set limit without additional cost. MP2 calculations are found to be as cost-effective as DFT ones and often slightly more, while showing a satisfactory accuracy when compared with the reference data. Although the focus is on barrier heights, the results may bear broader implications, in that one may see successes and difficulties of DFT when compared with traditional MO theories for the same data.

Benchmarking of Density Functionals for the Accurate Description of Thiol-Disulfide Exchange

A set of 92 density functionals was employed to accurately characterize thiol-disulfide exchange. The properties we have benchmarked throughout the study include the geometry of a 15 atoms model system, the potential energy surface, the activation barrier, and the energy of reaction for thiol-disulfide exchange. Reference energies were determined at the CCSD(T)/CBS//MP2/aug-cc-pVDZ level of theory, and reference geometries were calculated at the MP2/aug-cc-pVTZ level. M11-L, M06-2X, M06-HF, N12-SX, PBE1PBE, PBEh1PBE, and OHSE2PBE described better the geometry of the model system, with average deviations of 0.06 Å in bond lengths (0.06 Å in bond-breaking lengths) and 1.9° in bond angles. On the other hand, the potential energy surface and its gradient were more accurately described by the hybrid density functional BHandH, closely followed by mPW1N, mPW1K, and mPWB1K. The barrier height and energy of reaction were better reproduced by the BMK and M06-2X functionals (deviations of 0.17 and 0.07 kcal·mol(-1), respectively) for a set of 10 Pople's basis sets. MN12-SX and M11-L showed very good results for the widely used 6-311++G(2d,2p) basis set, with deviations of 0.02 and 0.05 kcal·mol(-1), respectively. We studied the effect of the split-valence, diffuse, and polarized functions in the activation barrier of thiol-disulfide exchange, for a set of 10 Pople's basis sets. While increasing the splitting and polarization may increase the activation barrier in approximately 1 kcal·mol(-1), diffuse functions generally contribute to decreasing it no more than 0.10 kcal·mol(-1). In general, 13 functionals provided energies within 1 kcal·mol(-1) of the reference value. The BB1K density functional is one of the best density functionals to characterize thiol-disulfide exchange reactions; however, several density functionals with modified Perdew-Wang exchange and about 40% Hartree-Fock exchange, such as mPW1K, mPW1N, and mPWB1K, show a good performance, too.

Quantum dynamics study on the CHIPR potential energy surface for the hydroperoxyl radical: the reactions O + OH⇋O2 + H

Quantum scattering calculations of the O((3)P)+OH((2)Π)⇌O2((3)Σg (-))+H((2)S) reactions are presented using the combined-hyperbolic-inverse-power-representation potential energy surface [A. J. C. Varandas, J. Chem. Phys. 138, 134117 (2013)], which employs a realistic, ab initio-based, description of both the valence and long-range interactions. The calculations have been performed with the ABC time-independent quantum reactive scattering computer program based on hyperspherical coordinates. The reactivity of both arrangements has been investigated, with particular attention paid to the effects of vibrational excitation. By using the J-shifting approximation, rate constants are also reported for both the title reactions.

Is HO3− multiple-minimum and floppy? Covalent to van der Waals isomerization and bond rupture of a peculiar anion

The HO₃⁻ anion is multiple-minimum and floppy: the two main isomers and isomerization barrier all lie quite below dissociation.

Accurate study of the two lowest singlet states of HN3: stationary structures and energetics at the MRCI complete basis set limit

The two lowest singlet potential energy surfaces of the hydrazoic acid are explored using high-level quantum chemistry calculations, revealing isomeric forms, transition states, and a new path for the nitrogen ring-closing mechanism in the ground state. The reaction of cyc-N3 with hydrogen is shown to proceed through a barrierless path that dissociates to N2 + NH((1)Δ) without reaching the HN3 global minimum. Several intersections between the two states are localized for both N3 isomers near the N3 + H dissociation. The energies of stationary structures and dissociation asymptotes are obtained with the multireference configuration interaction method at the complete basis set limit. A comparison with recent experimental data regarding the H-N3 bond strength and the dissociation barrier in the excited state shows that the present results are of chemical accuracy (~1 kcal mol(-1)).