Choice of the Right Supporting Electrolyte in Electrochemical Reductions: A Principal Component Analysis.

We present an analysis of a set of molecular, electrical, and electronic properties for a large number of the cations of quaternary ammonium salts usually employed as supporting electrolytes in cathodic reduction reactions. The goal of the present study is to define a measure for the quality of a supporting electrolyte in terms of the yield of the reaction considered. We performed a principal component analysis using the normalized values of the properties in order to lower the number of relevant reaction coordinates and find that the integral variance of 13 properties can well be represented by three principal components. The yield of the electrochemical hydrodimerization of acrylonitrile employing different quaternary ammonium salts as supporting electrolytes was determined in a series of experiments. We found only a very weak correlation between the yield and the values of the properties but a strong correlation between the yield and the values of the most important principal component. Very similar results are obtained for two further existing systematic experimental studies of the impact of the supporting electrolyte on the yield of cathodic reductions. For all three example reactions, a supervised regression using the two most important principal components as variables yields excellent values for the coefficients of determination. For comparison, we also applied our methodology to sets of purely structure-based features that are usually employed in cheminformatics and obtained results of almost similar quality. We therefore conjecture that our methodology in combination with a small number of experiments can be used to predict the yield of a given cathodic reduction on the basis of the properties of the supporting electrolyte.

Theoretical prediction of closed-shell paramagnetism for scandium and yttrium hydride.

Following chemical intuition, one would expect that all closed-shell molecules are diamagnetic. However, it is known that this is not the case for some second-row hydrides with low-lying unoccupied π orbitals due to an unquenching of the total angular momentum in the presence of an external magnetic field. In this article, the transition-metal hydrides ScH and YH are investigated, assuming a similar unquenching effect involving low-lying unoccupied π and δ orbitals formed from the metal d orbitals rather than the p orbitals. We are comparing results obtained with various quantum-chemical methods (HF, CCSD, CCSD(T), CCSDT) and basis sets. The obtained positive values for the magnetizabilities clearly indicate paramagnetic behavior. Vibrational effects on the magnetizability tensor are also considered, but these effects are small and do not change the overall conclusion that both ScH and YH are further examples for closed-shell paramagnetism.

Changing aromatic properties through stacking: the face-to-face dimer of Ni(II) bis(pentafluorophenyl)norcorrole.

Nuclear magnetic resonance (NMR) shielding constants have been calculated for Ni(II) bis(pentafluorophenyl)norcorrole and its face-to-face stacked dimer at the Hartree-Fock (HF), second-order Møller-Plesset perturbation theory (MP2), complete-active-space self-consistent-field (CASSCF) levels as well as at density functional theory (DFT) levels using several functionals. The calculated 1H NMR shielding constants agree rather well with the experimental ones. The shielding constants of N and Ni calculated at DFT, HF, and MP2 levels differ from those obtained in the CASSCF calculations due to near-degeneracy effects at the Ni atom. The calculated magnetically induced current densities show that the monomer is antiaromatic, sustaining a strong global paratropic ring current, and the dimer is aromatic, sustaining a strong diatropic ring current. Qualitatively the same current density is obtained at the employed levels of theory. The most accurate ring-current strengths are probably obtained at the MP2 level. The aromatic dimer has a short intermolecular distance of less than 3 Å. The intermolecular interaction changes the nature of the frontier orbitals leading to a formal double bond between the norcorrole macrocycles.

Error Control and Automatic Detection of Reference Active Spaces in Many-Body Expanded Full Configuration Interaction.

We present a wide-reaching revamp of the generalized many-body expanded full configuration interaction (MBE-FCI) method. First, we outline how to automatize the selection of reference active spaces, whereby the inherent bias introduced through a manual identification is reduced, also within the context of traditional complete active space methods. Second, we allow for the use of compact orbital clusters as expansion objects, which works to circumvent the unfavorable scaling with the number of orbitals included in the space complementary to the reference orbitals. Finally, we present a new algorithm for ensuring that many-body expansions can be efficiently terminated while conservatively accounting for resulting errors. These developments are all tested on a variety of molecular systems and different orbital representations to illustrate the abilities of our algorithm to produce correlation energies within predetermined error bounds, significantly broadening the overall applicability of the MBE-FCI method.

Cholesky Decomposition in Spin-Free Dirac-Coulomb Coupled-Cluster Calculations.

We present an implementation for the use of Cholesky decomposition (CD) of two-electron integrals within the spin-free Dirac-Coulomb (SFDC) scheme that enables to perform high-accuracy coupled-cluster (CC) calculations at costs almost comparable to those of their nonrelativistic counterparts. While for nonrelativistic CC calculations, atomic-orbital (AO)-based algorithms, due to their significantly reduced disk-space requirements, are the key to efficient large-scale computations, such algorithms are less advantageous in the SFDC case due to their increased computational cost in that case. Here, molecular-orbital (MO)-based algorithms exploiting the CD of the two-electron integrals allow us to reduce disk-space requirements and lead to computational cost in the CC step that is more or less the same as in the nonrelativistic case. The only remaining overhead in a CD-SFDC-CC calculation is due to the need to compute additional two-electron integrals, the somewhat higher cost of the Hartree-Fock calculation in the SFDC case, and additional cost in the transformation of the Cholesky vectors from the AO to the MO representation. However, these additional costs typically amount to less than 5-15% of the total wall time and are thus acceptable. We illustrate the efficiency of our CD scheme for SFDC-CC calculations on a series of illustrative calculations for the X(CO)4 molecules with X = Ni, Pd, Pt.

Ro-vibrational Spectrum of Linear Dialuminum Monoxide (Al2O) at 10 µm.

Dialuminum monoxide, Al2O, has been investigated in the laboratory at mid-IR wavelengths around 10 µm at high spectral resolution. The molecule was produced by laser ablation of an aluminum target with the addition of gaseous nitrous oxide, N2O. Subsequent adiabatic cooling of the gas in a supersonic beam expansion led to rotationally cold spectra. In total, 848 ro-vibrational transitions have been assigned to the fundamental asymmetric stretching mode ν3 and to five of its hot bands, originating from excited levels of the ν1 symmetric stretching mode and the ν2 bending mode. The measurements encompass 11 vibrational energy states (v1 v2l v3). The ro-vibrational transitions show spin statistical line intensity alternation of 7:5, which is caused by two identical aluminum nuclei of spin I = 5/2 at both ends of the centrosymmetric molecule of structure Al-O-Al. The less effective cooling of vibrational states in the supersonic beam expansion allowed measurement of transitions in excited vibrational states at energies of 1000 cm-1 and higher, while rotational levels within vibrational modes exhibited thermal population, with rotational temperatures around Trot = 115 K. Molecular parameters for 11 vibrational states were derived, including rotation and centrifugal distortion constants and l-type doubling constants for the states (v1 v2l v3) = (0 11 0) and (0 11 1) and an l-type resonance between the states (0 20 0) - (0 22 0) and (0 20 1) - (0 22 1). From the experimental results, rotational correction terms and the equilibrium bond length re were derived. The measurements were supported and guided by high-level quantum-chemical calculations that agree well with the derived experimental results.

A novel coupled-cluster singles and doubles implementation that combines the exploitation of point-group symmetry and Cholesky decomposition of the two-electron integrals.

A novel implementation of the coupled-cluster singles and doubles (CCSD) approach is presented that is specifically tailored for the treatment of large symmetric systems. It fully exploits Abelian point-group symmetry and the use of the Cholesky decomposition of the two-electron repulsion integrals. In accordance with modern CCSD algorithms, we propose two alternative strategies for the computation of the so-called particle-particle ladder term. The code is driven toward the optimal choice depending on the available hardware resources. As a large-scale application, we computed the frozen-core correlation energy of buckminsterfullerene (C60) with a polarized valence triple-zeta basis set (240 correlated electrons in 1740 orbitals).

Hyperfine-resolved spectra of HDS together with a global ro-vibrational analysis.

Despite their chemical simplicity, the spectroscopic investigation of light hydrides, such as hydrogen sulfide, is challenging due to strong hyperfine interactions and/or anomalous centrifugal-distortion effects. Several hydrides have already been detected in the interstellar medium, and the list includes H2S and some of its isotopologues. Astronomical observation of isotopic species and, in particular, those bearing deuterium is important to gain insights into the evolutionary stage of astronomical objects and to shed light on interstellar chemistry. These observations require a very accurate knowledge of the rotational spectrum, which is so far limited for mono-deuterated hydrogen sulfide, HDS. To fill this gap, high-level quantum-chemical calculations and sub-Doppler measurements have been combined for the investigation of the hyperfine structure of the rotational spectrum in the millimeter- and submillimeter-wave region. In addition to the determination of accurate hyperfine parameters, these new measurements together with the available literature data allowed us to extend the centrifugal analysis using a Watson-type Hamiltonian and a Hamiltonian-independent approach based on the Measured Active Ro-Vibrational Energy Levels (MARVEL) procedure. The present study thus permits to model the rotational spectrum of HDS from the microwave to far-infrared region with great accuracy, thereby accounting for the effect of the electric and magnetic interactions due to the deuterium and hydrogen nuclei.

Computation of NMR shieldings at the CASSCF level using gauge-including atomic orbitals and Cholesky decomposition.

We present an implementation of coupled-perturbed complete active space self-consistent field (CP-CASSCF) theory for the computation of nuclear magnetic resonance chemical shifts using gauge-including atomic orbitals and Cholesky decomposed two-electron integrals. The CP-CASSCF equations are solved using a direct algorithm where the magnetic Hessian matrix-vector product is expressed in terms of one-index transformed quantities. Numerical tests on systems with up to about 1300 basis functions provide information regarding both the computational efficiency and limitations of our implementation.

The He-H3 + complex. I. Vibration-rotation-tunneling states and transition probabilities.

With a He-H3 + interaction potential obtained from advanced electronic structure calculations, we computed the vibration-rotation-tunneling (VRT) states of this complex for total angular momenta J from 0 to 9, both for the vibrational ground state and for the twofold degenerate v2 = 1 excited state of H3 +. The potential has three equivalent global minima with depth De = 455.3 cm-1 for He in the plane of H3 +, three equatorial saddle points that separate these minima with barriers of 159.5 cm-1, and two axial saddle points with energies of 243.1 cm-1 above the minima. The dissociation energies calculated for the complexes of He with ortho-H3 + (oH3 +) and para-H3 + (pH3 +) are D0 = 234.5 and 236.3 cm-1, respectively. Wave function plots of the VRT states show that they may be characterized as weakly hindered internal rotor states, delocalized over the three minima in the potential and with considerable amplitude at the barriers. Most of them are dominated by the jk = 10 and 11 rotational ground states of oH3 + and pH3 +, with the intermolecular stretching mode excited up to v = 4 inclusive. However, we also found excited internal rotor states: 33 in He-oH3 +, and 22 and 21 in He-pH3 +. The VRT levels and wave functions were used to calculate the frequencies and line strengths of all allowed v2 = 0 → 1 rovibrational transitions in the complex. Theoretical spectra generated with these results are compared with the experimental spectra in Paper II [Salomon et al., J. Chem. Phys. 156, 144308 (2022)] and are extremely helpful in assigning these spectra. This comparison shows that the theoretical energy levels and spectra agree very well with the measured ones, which confirms the high accuracy of our ab initio He-H3 + interaction potential and of the ensuing calculations of the VRT states.

The He-H3 + complex. II. Infrared predissociation spectrum and energy term diagram.

The rotationally resolved infrared (IR) spectrum of the He-H3 + complex has been measured in a cryogenic ion trap experiment at a nominal temperature of 4 K. Predissociation of the stored complex has been invoked by excitation of the degenerate ν2 mode of the H3 + sub-unit using a pulsed optical parametric oscillator system. An assignment of the experimental spectrum became possible through one-to-one correlations with bands of the spectrum theoretically predicted in Paper I [Harding et al., J. Chem. Phys. 156, 144307 (2022)]. 19 bands have been assigned and analyzed, and the energy term diagram of the lower states of this floppy molecular complex has been derived from combination differences (CDs) in the experimental spectrum. Ground state combination differences (GSCDs) reveal a large part of the energy term diagram for the He-H3 + complex in its vibrational ground state, v = 0. Experimental and theoretical term energies agree within experimental accuracy for the rotational fine structure associated with the total angular momentum quantum number J and the parity e/f as well as for the coarse spacing of the lowest K states of the complex. This favorable comparison shows that the potential energy surface (PES) calculated in Paper I is accurate. The barriers between the three equivalent global minima in this PES are relatively low and the He-H3 + complex is extremely floppy, with nearly unhindered internal rotation of the H3 + sub-unit. The resulting Coriolis interactions couple the internal and end-over-end rotation of the complex and contribute significantly to the energy terms. They are observed both in experiment and theory and are, e.g., the origin of different rotational constants for states of e and f parity. Also in this respect, experiment and theory agree very well. Despite the assignment and analysis of many bands of the extremely rich IR spectrum of He-H3 +, higher levels of excitation, including the complex stretching mode, need further attention.

NMR chemical shift computations at second-order Møller-Plesset perturbation theory using gauge-including atomic orbitals and Cholesky-decomposed two-electron integrals.

We report on a formulation and implementation of a scheme to compute nuclear magnetic resonance (NMR) shieldings at second-order Møller-Plesset (MP2) perturbation theory using gauge-including atomic orbitals (GIAOs) to ensure gauge-origin independence and Cholesky decomposition (CD) to handle unperturbed and perturbed two-electron integrals. We investigate the accuracy of the CD for the derivatives of the two-electron integrals with respect to an external magnetic field and for the computed NMR shieldings, before we illustrate the applicability of our CD-based GIAO-MP2 scheme in calculations involving up to about 100 atoms and more than 1000 basis functions.

DLPNO-MP2 second derivatives for the computation of polarizabilities and NMR shieldings.

We present a derivation and efficient implementation of the formally complete analytic second derivatives for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory (MP2) method, applicable to electric or magnetic field-response properties but not yet to harmonic frequencies. We also discuss the occurrence and avoidance of numerical instability issues related to singular linear equation systems and near linear dependences in the projected atomic orbital domains. A series of benchmark calculations on medium-sized systems is performed to assess the effect of the local approximation on calculated nuclear magnetic resonance shieldings and the static dipole polarizabilities. Relative deviations from the resolution of the identity-based MP2 (RI-MP2) reference for both properties are below 0.5% with the default truncation thresholds. For large systems, our implementation achieves quadratic effective scaling, is more efficient than RI-MP2 starting at 280 correlated electrons, and is never more than 5-20 times slower than the equivalent Hartree-Fock property calculation. The largest calculation performed here was on the vancomycin molecule with 176 atoms, 542 correlated electrons, and 4700 basis functions and took 3.3 days on 12 central processing unit cores.

The Grignard Reaction - Unraveling a Chemical Puzzle.

More than 100 years since its discovery, the mechanism of the Grignard reaction remains unresolved. Ambiguities arise from the concomitant presence of multiple organomagnesium species and the competing mechanisms involving either nucleophilic addition or the formation of radical intermediates. To shed light on this topic, quantum-chemical calculations and ab initio molecular dynamics simulations are used to study the reaction of CH3MgCl in tetrahydrofuran with acetaldehyde and fluorenone as prototypical reagents. All organomagnesium species coexisting in solution due to the Schlenk equilibrium are found to be competent reagents for the nucleophilic pathway. The range of activation energies displayed by all of these compounds is relatively small. The most reactive species are a dinuclear Mg complex in which the substrate and the nucleophile initially bind to different Mg centers and the mononuclear dimethyl magnesium. The radical reaction, which requires the homolytic cleavage of the Mg-CH3 bond, cannot occur unless a substrate with a low-lying π*(CO) orbital coordinates the Mg center. This rationalizes why a radical mechanism is detected only in the presence of substrates with a low reduction potential. This feature, in turn, does not necessarily favor the nucleophilic addition, as shown for the reaction with fluorenone. The solvent needs to be considered as a reactant for both the nucleophilic and the radical reactions, and its dynamics is essential for representing the energy profile. The similar reactivity of several species in fast equilibrium implies that the reaction does not occur via a single process but by an ensemble of parallel reactions.

Vibrational Excitation Hindering an Ion-Molecule Reaction: The c-C_{3}H_{2}

Experiments within a cryogenic 22-pole ion trap have revealed an interesting reaction dynamic phenomenon, where rovibrational excitation of an ionic molecule slows down a reaction with a neutral partner. This is demonstrated for the low-temperature hydrogen abstraction reaction c-C_{3}H_{2}^{+}+H_{2}, where excitation of the ion into the ν_{7} antisymmetric C-H stretching mode decreased the reaction rate coefficient toward the products c-C_{3}H_{3}^{+}+H. Supported by high-level quantum-chemical calculations, this observation is explained by the reaction proceeding through a c-C_{3}H_{2}^{+}-H_{2} collision complex in the entrance channel, in which the hydrogen molecule is loosely bound to the hydrogen atom of the c-C_{3}H_{2}^{+} ion. This discovery enables high-resolution vibrational action spectroscopy for c-C_{3}H_{2}^{+} and other molecular ions with similar reaction pathways. Moreover, a detailed kinetic model relating the extent of the observed product depletion signal to the rate coefficients of inelastic collisions reveals that rotational relaxation of the vibrationally excited ions is significantly faster than the rovibrational relaxation, allowing for a large fraction of the ions to be vibrationally excited. This result provides fundamental insight into the mechanism for an important class of chemical reactions, and is capable of probing the inelastic collisional dynamics of molecular ions.

Ground and excited state first-order properties in many-body expanded full configuration interaction theory.

The recently proposed many-body expanded full configuration interaction (MBE-FCI) method is extended to excited states and static first-order properties different from total, ground state correlation energies. Results are presented for excitation energies and (transition) dipole moments of two prototypical, heteronuclear diatomics-LiH and MgO-in augmented correlation consistent basis sets of up to quadruple-ζ quality. Given that MBE-FCI properties are evaluated without recourse to a sampled wave function and the storage of corresponding reduced density matrices, the memory overhead associated with the calculation of general first-order properties only scales with the dimension of the desired property. In combination with the demonstrated performance, the present developments are bound to admit a wide range of future applications by means of many-body expanded treatments of electron correlation.

Coupled-cluster techniques for computational chemistry: The CFOUR program package.

An up-to-date overview of the CFOUR program system is given. After providing a brief outline of the evolution of the program since its inception in 1989, a comprehensive presentation is given of its well-known capabilities for high-level coupled-cluster theory and its application to molecular properties. Subsequent to this generally well-known background information, much of the remaining content focuses on lesser-known capabilities of CFOUR, most of which have become available to the public only recently or will become available in the near future. Each of these new features is illustrated by a representative example, with additional discussion targeted to educating users as to classes of applications that are now enabled by these capabilities. Finally, some speculation about future directions is given, and the mode of distribution and support for CFOUR are outlined.

Supramolecular Packing Drives Morphological Transitions of Charged Surfactant Micelles.

The shape and size of self-assembled structures upon local organization of their molecular building blocks are hard to predict in the presence of long-range interactions. Combining small-angle X-ray/neutron scattering data, theoretical modelling, and computer simulations, sodium dodecyl sulfate (SDS), over a broad range of concentrations and ionic strengths, was investigated. Computer simulations indicate that micellar shape changes are associated with different binding of the counterions. By employing a toy model based on point charges on a surface, and comparing it to experiments and simulations, it is demonstrated that the observed morphological changes are caused by symmetry breaking of the irreducible building blocks, with the formation of transient surfactant dimers mediated by the counterions that promote the stabilization of cylindrical instead of spherical micelles. The present model is of general applicability and can be extended to all systems controlled by the presence of mobile charges.

Gas-Phase Vibrational Spectroscopy of the Hydrocarbon Cations l-C3H+, HC3H+, and c-C3H2+: Structures, Isomers, and the Influence of Ne-Tagging.

We report the first gas-phase vibrational spectra of the hydrocarbon ions C3H+ and C3H2+. The ions were produced by electron impact ionization of allene. Vibrational spectra of the mass-selected ions tagged with Ne were recorded using infrared predissociation spectroscopy in a cryogenic ion trap instrument using the intense and widely tunable radiation of a free electron laser. Comparison of high-level quantum chemical calculations and resonant depletion measurements revealed that the C3H+ ion is exclusively formed in its most stable linear isomeric form, whereas two isomers were observed for C3H2+. Bands of the energetically favored cyclic c-C3H2+ are in excellent agreement with calculated anharmonic frequencies, whereas for the linear open-shell HCCCH+ (2Πg) a detailed theoretical description of the spectrum remains challenging because of Renner-Teller and spin-orbit interactions. Good agreement between theory and experiment, however, is observed for the frequencies of the stretching modes for which an anharmonic treatment was possible. In the case of linear l-C3H+, small but non-negligible effects of the attached Ne on the ion fundamental band positions and the overall spectrum were found.

Infrared Spectroscopy of Disilicon-Carbide, Si2C: The ν3 Fundamental Band.

The ν3 antisymmetric stretching mode of disilicon-carbide, Si2C, was studied using a narrow line width infrared quantum cascade laser spectrometer operating at 8.3 µm. The Si2C molecules were produced in an Nd:YAG laser ablation source from a pure silicon sample with the addition of a few percent methane diluted in a helium buffer gas. Subsequent adiabatic expansion was used to cool the gas down to rotational temperatures of a few tens of kelvin. A total of 183 infrared transitions recorded in the spectral range between 1200 and 1220 cm-1 were assigned to the fundamental ν3 mode of Si2C. In addition, pure rotational transitions of Ka = 1 and 2 between 278 and 375 GHz were recorded using a supersonic jet spectrometer for submillimeter wavelengths. Molecular parameters for the ( v1 v2 v3) = (001) vibrationally excited state were derived and improved molecular parameters for the vibrational ground-state (000) were obtained from a global fit data analysis, which includes our new laboratory data and millimeter wavelength data from the literature. We found the rotational levels Ka = 0 and Ka = 2 in the vibrationally excited (001) state being perturbed by a Coriolis-type interaction with energetically close lying levels of the symmetric stretching and triple-excited bending mode (130). The data analysis was supported by quantum chemical calculations performed at the coupled-cluster level of theory. All experimental results were found to be in excellent agreement with the theory.