Systematically Convergent Correlation Consistent Basis Sets for Molecular Core−Valence Correlation Effects: The Third-Row Atoms Gallium through Krypton

Toward Reliable Dipole Moments without Single Excitations: The Role of Orbital Rotations and Dynamical Correlation

The dipole moment is a crucial molecular property linked to a molecular system's bond polarity and overall electronic structure. To that end, the electronic dipole moment, which results from the electron density of a system, is often used to assess the accuracy and reliability of new electronic structure methods. This work analyses electronic dipole moments computed with the pair coupled cluster doubles (pCCD) ansätze and its linearized coupled cluster (pCCD-LCC) corrections using the canonical Hartree-Fock and pCCD-optimized (localized) orbital bases. The accuracy of pCCD-based dipole moments is assessed against experimental and CCSD(T) reference values using relaxed and unrelaxed density matrices and different basis set sizes. Our test set comprises molecules of various bonding patterns and electronic structures, exposing pCCD-based methods to a wide range of electron correlation effects. Additionally, we investigate the performance of pCCD-in-DFT dipole moments of some model complexes. Finally, our work indicates the importance of orbital relaxation in the pCCD model and shows the limitations of the linearized couple cluster corrections in predicting electronic dipole moments of multiple-bonded systems. Most importantly, pCCD with a linearized CCD correction can reproduce the dipole moment surfaces in singly bonded molecules, which are comparable to the multireference ones.

Benchmarking Ionization Potentials from pCCD Tailored Coupled Cluster Models

The ionization potential (IP) is an important parameter providing essential insights into the reactivity of chemical systems. IPs are also crucial for designing, optimizing, and understanding the functionality of modern technological devices. We recently showed that limiting the CC ansatz to the seniority-zero sector proves insufficient in predicting reliable and accurate ionization potentials within an IP equation-of-motion coupled-cluster formalism. Specifically, the absence of dynamical correlation in the seniority-zero pair coupled cluster doubles (pCCD) model led to unacceptably significant errors of approximately 1.5 eV. In this work, we aim to explore the impact of dynamical correlation and the choice of the molecular orbital basis (canonical vs localized) in CC-type methods targeting 230 ionized states in 70 molecules, comprising small organic molecules, medium-sized organic acceptors, and nucleobases. We focus on pCCD-based approaches as well as the conventional IP-EOM-CCD and IP-EOM-CCSD. Their performance is compared to the CCSD(T) or CCSDT equivalent and experimental reference data. Our statistical analysis reveals that all investigated frozen-pair coupled cluster methods exhibit similar performance, with differences in errors typically within chemical accuracy (1 kcal/mol or 0.05 eV). Notably, the effect of the molecular orbital basis, such as canonical Hartree-Fock or natural pCCD-optimized orbitals, on the IPs is marginal if dynamical correlation is accounted for. Our study suggests that triple excitations are crucial in achieving chemical accuracy in IPs when modeling electron detachment processes with pCCD-based methods.

Weak hydrogen bonding to halogens and chirality communication in propanols: Raman and microwave spectroscopy benchmark theory

Constitutional and conformational isomers of bromopropanol are vibrationally and rotationally characterised with parallels drawn to the structural chlorine analogues. A previous microwave spectroscopic study of the chloropropanols is re-examined and all systems are explored by Raman jet spectroscopy. For bromine, the entire nuclear quadrupole coupling tensors are accurately determined and compared to their chlorine counterparts. Tensor asymmetry parameters are determined and linked with the hydrogen bond strength as indicated by the downshift of the OH-stretching frequency. The spectroscopic constants derived from the observed transitions are used as benchmarks for a large variety of electronic structure methods followed by harmonic and anharmonic rovibrational treatments. The CCSD(T) electronic structure calculations provide the best performance, in particular once anharmonic and relativistic corrections are applied or implied. Standard DFT approaches vary substantially with respect to their systematic error cancellation across the investigated species, and cost-effective compromises for the different observables are proposed.

Spectroscopic detection of the gallium methylene (GaCH2 and GaCD2) free radical in the gas phase by laser-induced fluorescence and emission spectroscopy

GaCH2, a free radical thought to play a role in the chemical vapor deposition of gallium-containing thin films and semiconductors, has been spectroscopically detected for the first time. The radical was produced in a pulsed discharge jet using a precursor mixture of trimethylgallium vapor in high pressure argon and studied by laser-induced fluorescence and wavelength resolved emission techniques. Partially rotationally resolved spectra of the hydrogenated and deuterated species were obtained, and they exhibit the nuclear statistical weight variations and subband structure expected for a 2A2-2B1 electronic transition. The measured spectroscopic quantities have been compared to our own ab initio calculations of the ground and excited state properties. The electronic spectrum of gallium methylene is similar to the corresponding spectrum of the aluminum methylene radical, which we reported in 2022.

Quantum computations in heavy noble-gas hydride cations: Reference energies and new spectroscopic data

Computational quantum chemistry has become a powerful tool with a wide range of possibilities to solve chemical-physical problems. As a result of this, the interest in the applications of computational quantum chemistry has expanded considerably, and has opened up novel research opportunities. In particular, those related to the characterization of heavy-atoms complexes, as most electronic structure calculations for such systems struggle with the problem posed by the large number of electrons present in them, and consequently, the introduction of relativistic effects. The present study performed an exhaustive assess to characterized the uncommon NgH+ (Ng = Kr, Xe, and Rn) hydride cations in order to provide accurate rovibrational data of their isotopes to assist in the laboratory characterization or even their astronomical detection. Scalar relativistic effects were included, and the ground and first electronically exited states potential curves were obtained from benchmark ab initio CCSD(T)/CBS and MRCI+Q electronic structure calculations. Next, such interaction potentials, correctly extended to long-range asymptotic regions, were employed in quantum bound state calculations and molecular spectroscopic constants were determined for the most abundance 84Kr, 132Xe, and 222Rn isotopes. Our results were discussed in comparison with available experimental and previous theoretical estimates, aiming to treat accuracy issues. The new sets provide reference data that could serve for spectroscopic characterization of such low abundance and high radioactive species.

Accuracy of Different Electronic Basis Set Families for Anharmonic Molecular Vibrations: A Comprehensive Benchmark Study

In this work, the accuracy and convergence of different electronic basis set families for the computation of anharmonic molecular vibrational spectroscopic calculations are benchmarked. A series of 39 different basis sets from different families following their hierarchy are assessed on VSCF and VSCF-PT2 algorithms with commonly used MP2 and DFT based B3LYP-D potentials for a set of molecular systems. Such an effort has been validated in a previous work ( J. Phys. Chem. A 2020, 124, 9203-9221) with split-valence basis sets for fundamentals and intensities. Here, fundamental transitions, vibrationally excited states, and intensities are compared with the experimental data to estimate the accuracy for a series of Jensen, Dunning, Calendar, Karlsruhe, and Sapporo basis set families. The convergence of basis sets are also compared with the large ANO basis set. Comprehensive statistical error analysis in terms of accuracy and precision was carried out to assess the performance of each basis set. It is observed that the improvement for the calculated harmonic and anharmonic values from the smaller basis sets to the medium (i.e., triple-ξ) is considerable. Beyond this, from medium to large basis sets, the convergence is slow and mostly posits nearly converged values. Basis sets with and without diffuse functions offer characteristically different accuracies and convergence patterns. Finally, recommendations are given on the choice of basis set chosen as black-box which can balance between accuracy and computational time, estimation of the errors, and their selections especially for large molecules.

New efficient pecS-n (n = 1, 2) basis sets for quantum chemical calculations of 31P NMR chemical shifts

The basis sets that are used in the quantum chemical calculations of ³¹P NMR chemical shifts have always been one of the most important factors of accuracy. Regardless of what high-quality approach is employed, using basis sets of insufficient flexibility in the important angular regions may give poor results and lead to misassignments of the signals in the ³¹P NMR spectra. In this work, it was found that the existing nonrelativistic basis sets for phosphorus atom of double- and triple-ζ quality, specialized for the ³¹P NMR chemical shifts calculations, are essentially undersaturated in the d-angular space that occurred to play a significant role in the overall accuracy of these calculations. This problem has been thoroughly investigated, and new pecS-n (n = 1, 2) basis sets for phosphorus chemical shifts calculations were proposed. The exponents and contraction coefficients for the pecS-n basis sets were generated with the property-energy consistent method that has been introduced in our earlier paper, and has been proven useful in the creation of efficient property-oriented basis sets. New basis sets were optimized using the GIAO-DFT method with the B97-2 functional. Extensive benchmark calculations showed that the pecS-1 and pecS-2 basis sets provide very good accuracy, characterized by the corrected mean absolute percentage errors against the experiment of about 7.03 and 4.42 ppm, respectively. In particular, the accuracy of ³¹P NMR chemical shifts calculations achieved with the pecS-2 basis set is one of the most favorable accuracies for today. We believe that our new pecS-n (n = 1, 2) basis sets for phosphorus atom will prove useful in modern large-scale quantum chemical calculations of ³¹P NMR chemical shifts.

Bromine Pentafluoride BrF5 , the Formation of [BrF6 ]- Salts, and the Stereochemical (In)activity of the Bromine Lone Pairs

BrF₅ can be prepared by treating BrF₃ with fluorine under UV light in the region of 300 to 400 nm at room temperature. It was analyzed by UV-Vis, NMR, IR and Raman spectroscopy. Its crystal structure was redetermined by X-ray diffraction, and its space group was corrected to Pnma. Quantum-chemical calculations were performed for the band assignment of the vibrational spectra. A monoclinic polymorph of BrF₅ was quantum chemically predicted and then observed as its low-temperature modification in space group P2₁ /c by single crystal X-ray diffraction. BrF₅ reacts with the alkali metal fluorides AF (A=K, Rb) to form alkali metal hexafluoridobromates(V), A[BrF₆ ] the crystal structures of which have been determined. Both compounds crystallize in the K[AsF₆ ] structure type (R 3 ‾ ${\bar 3}$ , no. 148, hR24). For the species [BrF₆ ]⁺ , BrF₅ , [BrF₆ ]^- , and [IF₆ ]^- , the chemical bonds and lone pairs on the heavy atoms were investigated by means of intrinsic bond orbital analysis.

High-Level Theoretical Study of the C1Π States of GaCl, GaBr, InCl, and InBr

For the first singlet excited states C¹Π of IIIA group monohalides GaCl, GaBr, InCl, and InBr, a very shallow potential well may exist in the flat potential energy curve, which poses a challenge to the theoretical accuracy. In this study, high-level theoretical calculations have been performed through the Feller-Peterson-Dixon composite approach to study the C¹Π states, and the obtained spectroscopic constants agree well with the known experimental ones. It is found that the C¹Π states are sensitive to the size of basis functions, spin-orbit coupling, and strong correlations mainly due to triple excitations. The final results show that the C¹Π states of GaCl and InCl are quasi-bound with one (v' = 0) and four (v' = 0-3) vibrational levels, respectively, being consistent with the experimental findings, whereas the C¹Π states of GaBr and InBr are repulsive. Our conclusions deny the existence of higher vibrational levels v' = 1 for GaCl, v' > 3 for InCl, and v' ≥ 0 for InBr in previous experimental and theoretical studies of C¹Π.

Atomic Mean-Field Approach within Exact Two-Component Theory Based on the Dirac-Coulomb-Breit Hamiltonian

An extension of the exact two-component theory with atomic mean-field integrals (the X2CAMF scheme) to the treatment of the Breit term together with efficient implementation using an atomic Dirac-Coulomb-Breit Hartree-Fock program is reported. The accuracy of the X2CAMF scheme for treating the contributions from the Breit term to the molecular properties is demonstrated using benchmark calculations of equilibrium bond lengths, harmonic frequencies, and dipole moments for molecules containing elements across the periodic table. Calculations of the properties for molecules containing period four elements aiming at high accuracy as well as for Th- and U-containing molecules are also presented and compared with experimental results to demonstrate the usefulness of the X2CAMF scheme in combination with accurate treatments of electron correlation by the coupled-cluster (CC) methods. The combination of CC methods and the X2CAMF scheme shows potential to extend the accuracy of CC calculations to heavy elements, e.g., to computational heavy-element thermochemistry.

Observation of Selectively Populated Monohalide Excited States from the Reactions of Group 3 Metal (Sc, Y, and La) Monomers and Dimers with Halogen-Containing Molecules

The chemiluminescent reactions of the group 3 metals Sc and Y with F₂, Cl₂, Br₂, ClF, ICl (Sc), IBr (Y), and SF₆ and La with F₂, SF₆, Cl₂, and ClF have been studied at low pressures (6 × 10^-6 to 4 × 10^-4 Torr) using a beam-gas arrangement and extended to the 10^-3 Torr multiple collision pressure range. Contrary to previous reports, the observed chemiluminescent spectra are primarily attributed to emission from the metal monohalides. Extensive pressure and temperature dependence studies and high-level correlated molecular orbital theory calculations of the bond dissociation energies support this conclusion and the attribution of the chemiluminescence. Evidence for the "selective" production of a monohalide excited electronic state is obtained for several of the Sc and Y reactions. All reactions producing the metal monofluorides are first order with respect to the oxidant, while reactions producing the monochlorides and monobromides are found to be "faster than first order" with respect to the oxidant. This difference is associated with the metal halide bond dissociation energies and the metal halide product internal density of states. Analysis of the temperature dependence for six representative reactions indicates that the "selective" excited-state formation of the metal monohalides proceeds via a direct mechanism with negligible activation energy. We compare and contrast the present results with previous experiments and interpretations which have assigned the selective emission from these systems to the group 3 dihalides produced in a two-step reaction sequence analogous to an electron jump process. The current results suggest a distinctly different interpretation of the observed processes in these systems. The observed selectivity observed in these studies is remarkable given the significant number of known and potential excited states in the scandium and yttrium halides as well as their different electronic configurations.

Energies, structures, and harmonic frequencies of small water clusters from the direct random phase approximation

The binding energies, structures, and vibrational frequencies of water clusters up to 20 molecules are computed at the direct random phase approximation (RPA) level of theory and compared to theoretical benchmarks. Binding energies of the WATER27 set, which includes neutral and positively and negatively charged clusters, are predicted to be too low in the complete basis set limit by an average of 7 kcal/mol (9%) and are worse than the results from the best density functional theory methods or from the Møller-Plesset theory. The RPA shows significant basis set size dependence for binding energies. The order of the relative energies of the water hexamer and dodecamer isomers is predicted correctly by the RPA. The mean absolute deviation for angles and distances for neutral clusters up to the water hexamer are 0.2° and 0.6 pm, respectively, using quintuple-ζ basis sets. The relative energetic order of the hexamer isomers is preserved upon optimization. Vibrational frequencies for these systems are underestimated by several tens of wavenumbers for large basis sets, and deviations increase with the basis set size. Overall, the direct RPA method yields accurate structural parameters but systematically underestimates binding energies and shows strong basis set size dependence.

Bond Dissociation Energies in Heavy Element Chalcogen and Halogen Small Molecules

Thermodynamic properties including bond dissociation energies (BDEs), heats of formation, and gas-phase acidities for the hydrides and dimers of chalcogens and halogens, H₂Y, HX, Y₂, and X₂ for Y = Se, Te, and At and X = Br, I, and At, have been predicted using the Feller-Peterson-Dixon composite-correlated molecular orbital theory approach. A full four-component CCSD(T) approach was used to calculate the spin-orbit effects on thermodynamic properties, except for Se₂, where the AoC-DHF value was used due to strong multireference effects in Se₂ for the SO calculations. The calculated results show that the At₂ BDE is quite small, 19.5 kcal/mol, with much of the low bond energy due to spin-orbit effects. H₂Po is not predicted to be stable to dehydrogenation to Po + H₂ in terms of the free energy at 298 K. In the gas phase, HAt is predicted to be a stronger acid than H₂SO₄. The current results provide insights into potential difficulties in the actual experimental observation of such species for heavy elements.

Exploring the excited states of the GeH+ radical cation including spin-orbit interaction: A revisited study

In this paper, we have carried out high-level ab initio calculations on the electronic states of GeH⁺ with the configuration interaction method. The spin-orbit coupling (SOC), core-valence correlation (CV), scalar relativistic effects and Davidson correction (+Q) are included. The potential energy curves (PECs) of 13 Λ-S correlated with the four lowest dissociation limits and 32 Ω electronic states generated from those Λ-S states are obtained. Our results indicate that the first ³Σ^- and second ³Π states are adiabatically correlated with the dissociation limit Ge(³P_g) + H⁺(¹S), which is different from the previously reported Ge⁺(⁴P_g) + H(²S_g). From the computed PECs, the spectroscopic constants of the Λ-S and Ω states are determined, which are in good agreement with previous experiments. The dipole moments (DMs) for Λ-S electronic states are also investigated. With the help of spin-orbit coupling matrix involving the 1³Σ^- and 2³Π states, the intricate couplings related with the crossing states are revealed, and the weak predissociation for ν' ≥0 vibrational levels of 1³Σ^- and the perturbations for vibrational levels of 2¹Σ⁺ (ν' ≥1) and 1¹Δ (ν' ≥0) states are analyzed. Finally, the transition properties of five transitions are predicted, including the Franck-Condon Factors (FCFs), transition dipole moments (TDMs), and the spontaneous radiative lifetimes of lower vibrational states. This study will improve our comprehension on the detailed electronic structure and spectroscopy of GeH⁺ radical cation.

Quantum chemical benchmark databases of gold-standard dimer interaction energies

Advances in computational chemistry create an ongoing need for larger and higher-quality datasets that characterize noncovalent molecular interactions. We present three benchmark collections of quantum mechanical data, covering approximately 3,700 distinct types of interacting molecule pairs. The first collection, which we refer to as DES370K, contains interaction energies for more than 370,000 dimer geometries. These were computed using the coupled-cluster method with single, double, and perturbative triple excitations [CCSD(T)], which is widely regarded as the gold-standard method in electronic structure theory. Our second benchmark collection, a core representative subset of DES370K called DES15K, is intended for more computationally demanding applications of the data. Finally, DES5M, our third collection, comprises interaction energies for nearly 5,000,000 dimer geometries; these were calculated using SNS-MP2, a machine learning approach that provides results with accuracy comparable to that of our coupled-cluster training data. These datasets may prove useful in the development of density functionals, empirically corrected wavefunction-based approaches, semi-empirical methods, force fields, and models trained using machine learning methods.

Spectroscopic properties and spin-orbit coupling of electronic excited states of the germanium dimer

The electronic structure and spectroscopic properties of the germanium dimer have been investigated by high-level ab initio calculations with consideration of spin-orbit coupling (SOC). The potential energy curves (PECs) of 19 Λ-S states, as well as those of 52 Ω states generated from the Λ-S states, are calculated by the multireference configuration interaction plus Davidson correction (MRCI + Q) method. The properties of the F3Σ+u1-X3Σ-g1 and H3Σ-u1-X3Σ-g1 transitions as well as the interactions of the F3Σ+u and H3Σ-u states with other excited states induced by SOC are investigated based on the calculated SO matrix and the PECs of the Ω states. Our results indicate that the previously observed spectra of Ge2 in the range 20 500-22 000 cm-1 should be assigned as the transition between the Ω = 1g component of the X3Σ-g state and Ω = 1u of the F3Σ+u state. Moreover, owing to the strong SOC with the repulsive 25Πu state, the H3Σ-u state is predissociative, leading to the Ge(3P2) + Ge(3P1) channel at vibrational levels higher than v' = 6. Our theoretical study would provide comprehensive information and shed light on understanding the spectroscopy and dynamics of the electronic excited states of Ge2.

Probing the Potential Energy Surfaces of BrCN- by Dissociative Electron Attachment

State coupling certainly determines the topologic features of the molecular potential energy surface (PES) and potentially diversifies chemical reaction pathways. Here we report the new PESs of BrCN^- in the low-lying electronic states that are distinctly different from the previous predictions in the short Br-CN bond region but validated by the high-resolution ion velocity imaging measurements of low-energy dissociative electron attachment (DEA) to BrCN. Besides the vibrating CN^- ions produced in the fast Br-CN bond stretching motions, we confirm that the ro-vibrating CN^- ions with a nearly isotropic angular distribution are produced by receiving a torque in the combinational motion of Br-CN bond bending and stretching. The latter process is closely related to the potential well of BrCN^- at the first excited state A²Π_3/2 that arises from the Π-Σ state couplings. Our findings not only suggest that the PESs of other anionic cyanogen halides are in dire need of reexamination but also show that ion velocity imaging of the DEA process is a powerful experimental method for evaluating the theoretical PESs of molecular anions.

A dual attack on the peroxide bond. The common principle of peroxidatic cysteine or selenocysteine residues

The (seleno)cysteine residues in some protein families react with hydroperoxides with rate constants far beyond those of fully dissociated low molecular weight thiol or selenol compounds. In case of the glutathione peroxidases, we could demonstrate that high rate constants are achieved by a proton transfer from the chalcogenol to a residue of the active site [Orian et al. Free Radic. Biol. Med. 87 (2015)]. We extended this study to three more protein families (OxyR, GAPDH and Prx). According to DFT calculations, a proton transfer from the active site chalcogenol to a residue within the active site is a prerequisite for both, creating a chalcogenolate that attacks one oxygen of the hydroperoxide substrate and combining the delocalized proton with the remaining OH or OR, respectively, to create an ideal leaving group. The “parking postions” of the delocalized proton differ between the protein families. It is the ring nitrogen of tryptophan in GPx, a histidine in GAPDH and OxyR and a threonine in Prx. The basic principle, however, is common to all four families of proteins. We, thus, conclude that the principle outlined in this investigation offers a convincing explanation for how a cysteine residue can become peroxidatic.

•

In some protein families, (seleno)cysteine residues react with hydroperoxides with very high rate constants.

•

In GPx, DFT models of the oxidation of the catalytic site support a two-step mechanism for the H₂O₂ reduction.

•

This mechanism is here found to operate in other thiol-based enzymes, i.e. OxyR, GAPDH and Prx.

The spectroscopic properties of the low-lying excited states and laser cooling scheme of SrBr molecule

The potential energy curves and the transition dipole moments for seven electronic states of SrBr molecule are obtained via the multi-reference configuration interaction method and the all-electron basis sets. The Davidson and relativistic corrections are also included. Based on the obtained potential energy curves, the rotational and vibrational energy levels of each electronic state are determined by solving the nuclear motion equation of the molecule. The spectroscopic parameters are fitted from the obtained energy levels by using Dunham expression. Moreover, the spin-orbit coupling splits of the A²Π state are considered to construct the optical laser cooling scheme. The Frank-Condon factors, radiation lifetimes, radiation widths between the ground electronic state and ²Π_{1/2, 3/2}/B²Σ⁺ states are calculated. Then, the feasibility of laser cooling is explored and the optical scheme is proposed. The results demonstrate that the SrBr molecule is a promising candidate for laser cooling.

Relatives of cyanomethylene: replacement of the divalent carbon by B-, N+, Al-, Si, P+, Ga-, Ge, and As. Phys Chem Chem Phys 2019;21:26438-26452. [PMID: 31774089 DOI: 10.1039/c9cp05777c]

The lowest lying singlet and triplet states of HBCN-, HCCN, HNCN+, HAlCN-, HSiCN, HPCN+, HGaCN-, HGeCN, and HAsCN+ were studied using the CCSDT(Q)/CBS//CCSD(T)/aug-cc-pVQZ level of theory. Periodic trends in geometries, singlet-triplet gaps, and barriers to linearity were established and analyzed. The first row increasingly favors the triplet state, with a singlet-triplet gap (ΔEST = Esinglet - Etriplet) of 3.5 kcal mol-1, 11.9 kcal mol-1, and 22.6 kcal mol-1, respectively, for HBCN-, HCCN, and HNCN+. The second row increasing favors the singlet state, with singlet-triplet gaps of -20.4 kcal mol-1 (HAlCN-), -26.6 kcal mol-1 (HSiCN), and -26.8 kcal mol-1 (HPCN+). The third row also favors the singlet state, with singlet-triplet gaps of -26.8 kcal mol-1 (HGaCN-), -33.5 kcal mol-1 (HGeCN), and -33.1 kcal mol-1 (HAsCN+). The HXCN species have larger absolute singlet-triplet energy gaps compared to their parent species XH2 except for the case of X = N+. The effect of the substitution of hydrogen with a cyano group was analyzed with isodesmic bond separation analysis and NBO.

Generation and structural characterization of Ge carbides GeCn (n = 4, 5, 6) by laser ablation, broadband rotational spectroscopy, and quantum chemistry

Following the recent discovery of T-shaped GeC₂, rotational spectra of three larger Ge carbides, linear GeC₄, GeC₅, and GeC₆ have been observed using chirped pulse and cavity Fourier transform microwave spectroscopy and a laser ablation molecule source, guided by new high-level quantum chemical calculations of their molecular structure. Like their isovalent Si-bearing counterparts, Ge carbides with an even number of carbon atoms beyond GeC₂ are predicted to possess ¹Σ ground electronic states, while odd-numbered carbon chains are generally ³Σ; all are predicted to be highly polar. For the three new molecules detected in this work, rotational lines of four of the five naturally occurring Ge isotopic variants have been observed between 6 and 22 GHz. Combining these measurements with ab initio force fields, the Ge-C bond lengths have been determined to high precision: the derived values of 1.776 Å for GeC₄, 1.818 Å for GeC₅, and 1.782 Å for GeC₆ indicate a double bond between these two atoms. Somewhat surprisingly, the spectrum of GeC₅ very closely resembles that of a ¹Σ molecule, implying a spin-spin coupling constant λ in excess of 770 GHz for this radical, a likely consequence of the large spin-orbit constant of atomic Ge (∼1000 cm^-1). A systematic comparison between the production of SiC_n and GeC_n chains by laser ablation has also been undertaken. The present work suggests that other large metal-bearing molecules may be amenable to detection by similar means.

Beyond chemical accuracy in the heavy p-block: The first ionization potentials and electron affinities of Ga-Kr, In-Xe, and Tl-Rn

A relativistic coupled-cluster version of the Feller-Peterson-Dixon composite method has been used to accurately calculate the first ionization potentials (IPs) and electron affinities (EAs) of the post-d, p-block elements Ga-Rn. Complete basis set extrapolations including outer-core correlation at the CCSD(T) level of theory were combined with contributions from higher order electron correlation up to CCSDTQ, quantum electrodynamic effects (Lamb shift), and spin-orbit (SO) coupling including the Gaunt contribution. Several methods for including SO were investigated, in which all involved the four-component (4c) Dirac-Coulomb (DC) Hamiltonian. The treatment of SO coupling was the contribution that limited the final accuracy of the present results. In the cases where 4c-DC-CCSD(T) could be reliably used for the SO contributions, the final composite IPs and EAs agreed with the available experimental values to within an unsigned average error of just 0.16 and 0.20 kcal/mol, respectively. In all cases, the final IPs and EAs were within 1 kcal/mol of the available experimental values, except for the EAs of the group 13 elements (Ga, In, and Tl), where the currently accepted experimental values appear to be too large by as much as 4 kcal/mol. The values predicted in this work, which have estimated uncertainties of ±0.5 kcal/mol, are 5.25 (Ga), 7.69 (In), and 7.39 (Tl) kcal/mol. For the EAs of Po and At, which do not have experimental values, the current calculations predict values of 34.2 and 55.8 kcal/mol with estimated uncertainties of ±0.6 and ±0.3 kcal/mol, respectively.

Ab initio structure and vibration-rotation dynamics of the formyl and isoformyl cations, HCO+/HOC. J Chem Phys 2019;150:154307. [PMID: 31005073 DOI: 10.1063/1.5089718]

Accurate structure and potential energy surface of the formyl and isoformyl cation system, HCO⁺/HOC⁺, in its ground electronic state X̃ ¹Σ⁺ have been determined from ab initio calculations using the coupled-cluster approach in conjunction with the correlation-consistent basis sets up to septuple-zeta quality. Both the isomers are confirmed to be linear at equilibrium, with the total energy minimum of HOC⁺ lying 14 120 cm^-1 above that of HCO⁺ and the HCO⁺ → HOC⁺ isomerization energy barrier being 26 870 cm^-1 (in the Born-Oppenheimer approximation). The equilibrium structural parameters for HCO⁺ are estimated to be r_e(HC) = 1.0919 Å and r_e(CO) = 1.1058 Å, whereas those for HOC⁺ are estimated to be r_e(HO) = 0.9899 Å and r_e(CO) = 1.1544 Å. The vibration-rotation energy levels were predicted for various isotopologues using a variational approach and compared with the experimental data. For the spectroscopically well characterized formyl cation, the observed vibration-rotation energies and the effective rotational constants are reproduced to within about 2.3 cm^-1 and 1.7 MHz, respectively. The role of the core-electron correlation, higher-order valence-electron correlation, scalar relativistic, and adiabatic effects in determining the structure and vibration-rotation dynamics of both the isomers is discussed.

Titanocene Selenide Sulfides Revisited: Formation, Stabilities, and NMR Spectroscopic Properties

[TiCp₂S₅] (phase A), [TiCp₂Se₅] (phase F), and five solid solutions of mixed titanocene selenide sulfides [TiCp₂Se_xS_5−x] (Cp = C₅H₅⁻) with the initial Se:S ranging from 1:4 to 4:1 (phases B–E) were prepared by reduction of elemental sulfur or selenium or their mixtures by lithium triethylhydridoborate in thf followed by the treatment with titanocene dichloride [TiCp₂Cl₂]. Their ⁷⁷Se and ¹³C NMR spectra were recorded from the CS₂ solution. The definite assignment of the ⁷⁷Se NMR spectra was based on the PBE0/def2-TZVPP calculations of the ⁷⁷Se chemical shifts and is supported by ¹³C NMR spectra of the samples. The following complexes in varying ratios were identified in the CS₂ solutions of the phases B–E: [TiCp₂Se₅] (5₁), [TiCp₂Se₄S] (4₁), [TiCp₂Se₃S₂] (3₁), [TiCp₂SSe₃S] (3₆), [TiCp₂SSe₂S₂] (2₅), [TiCp₂SSeS₃] (1₂), and [TiCp₂S₅] (0₁). The disorder scheme in the chalcogen atom positions of the phases B–E observed upon crystal structure determinations is consistent with the spectral assignment. The enthalpies of formation calculated for all twenty [TiCp₂Se_xS_5−x] (x = 0–5) at DLPNO-CCSD(T)/CBS level including corrections for core-valence correlation and scalar relativistic, as well as spin-orbit coupling contributions indicated that within a given chemical composition, the isomers of most favourable enthalpy of formation were those, which were observed by ⁷⁷Se and ¹³C NMR spectroscopy.

A comparison between hydrogen and halogen bonding: the hypohalous acid–water dimers, HOX⋯H2O (X = F, Cl, Br)

Hypohalous acids (HOX) are a class of molecules that play a key role in the atmospheric seasonal depletion of ozone and have the ability to form both hydrogen and halogen bonds.

Atomic Charges in Highly Ionic Diatomic Molecules

Atomic charges were investigated as functions of detectable atomic and molecular constants at equilibrium structures. It was found based upon the variation idea that atomic charges in highly ionic molecules can be expressed as a function of molecular dipole moments, polarizabilities of free cations, and polarizabilities of free neutral atoms of the corresponding anions. The function can be given in the form of classical Rittner's relationship (J. Chem. Phys. 1951, 19, 1030). For the ground states of alkali halide molecules, the predicted atomic charges are close to an elementary charge e and the predicted dipole moments are in good agreement with the observed values; for spin-restricted high-ionic systems such as the lowest ⁹Σ electronic states of BN, AlN, GaN, BP, AlP, GaP, BAs, AlAs, and GaAs molecules, the predicted atomic charges are also near 1e and in good agreement with the results of natural population analysis at MRCI/cc-pvqz and HF/6-311+G(3df) levels. Polarizabilities for the lowest quintet states of B^-, Al^-, Ga^-, N⁺, P⁺, and As⁺ ions were also obtained based upon high-level ab initio computations. Atomic charges from other related methods are also investigated for comparison. The results demonstrate that high-quality atomic charges can be obtained with detectable variables, such as molecular dipole moment, vibrational frequency, as well as polarizabilities of the related free atoms and ions.

Theoretical Study of the Reactions of H Atoms with CH3I and CH2I2

High level ab initio methods have been used to provide reliable kinetic data for the H + CH₃I and H + CH₂I₂ gas-phase reactions. The (H, I)-abstraction and I-substitution reaction pathways were identified. The structures were determined on the potential energy surface at the MP2/aug-cc-pVTZ level of theory. The energetics was then refined using the coupled cluster theory. For the iodinated species, the spin-orbit coupling was calculated using the MRCI approach. The core valence and the scalar relativistic corrections were considered. Thermal rate constants were reported using the canonical transition-state theory (TST) and compared to computed values with the canonical variational transition-state theory (CVT) using the zero curvature tunneling (ZCT) and the small curvature tunneling (SCT) corrections over a wide temperature range (250-2500 K) to show the importance of quantum tunneling effects at low temperatures. They are given by the following expressions for the overall reactions using the CVT/SCT method: k_H+CH3I( T) = 1.07 × 10^-17 × T^2.13 exp(2.68 (kJ mol^-1)/ RT) and k_H+CH2I2( T) = 5.73 × 10^-21 × T^2.97 exp(3.15 (kJ mol^-1)/ RT). The I-abstraction is predicted to be the major pathway for both H + CH₃I and H + CH₂I₂ reactions. The obtained kinetic parameters for the H + CH₃I reaction are in excellent agreement with their experimental counterparts over the temperature range 300-750 K. On the basis of our calculated reaction enthalpies, a new evaluation of the standard enthalpy of formation at 298 K of CH₂I and CHI₂ has been provided. Obtained values are Δ_f H°_298K (CH₂I) = 219.5 kJ mol^-1 and Δ_f H°_298K(CHI₂) = 296.3 kJ mol^-1.

State-of-the-art ab initio potential energy curve for the xenon atom pair and related spectroscopic and thermophysical properties

A new ab initio interatomic potential energy curve for two ground-state xenon atoms is presented. It is based on supermolecular calculations at the coupled-cluster level with single, double, and perturbative triple excitations [CCSD(T)] employing basis sets up to sextuple-zeta quality, which were developed as part of this work. In addition, corrections were determined for higher coupled-cluster levels up to CCSDTQ as well as for scalar and spin-orbit relativistic effects at the CCSD(T) level. A physically motivated analytical function was fitted to the calculated interaction energies and used to compute the vibrational spectrum of the dimer, the second virial coefficient, and the dilute gas transport properties. The agreement with the best available experimental data for the investigated properties is excellent; the new potential function is superior not only to previous ab initio potentials but also to the most popular empirical ones.

Reactivity of Hydrogen Peroxide with Br and I Atoms

The reaction mechanisms of Br and I atoms with H₂O₂ have been investigated using DFT and high-level ab initio calculations. The H-abstraction and OH-abstraction channels were highlighted. The geometries of the stationary points were optimized at the B3LYP/aug-cc-pVTZ level of theory, and the energetics were recalculated with the coupled cluster theory. Spin-orbit coupling for each halogenated species was also explicitly computed by employing the MRCI level of theory. Thermochemistry for HOBr and HOI has been revised and updated standard enthalpies of formation at 298 K for HOBr and HOI are the following: Δ_fH°_298K(HOBr) = (-66.2 ± 4.6) kJ mol^-1 and Δ_fH°_298K(HOI) = (-66.8 ± 4.7) kJ mol^-1. The rate constants have been estimated using transition state theory (TST), canonical variational transition state theory (CVT), and CVT with small curvature tunneling (CVT/SCT) over a wide temperature range (250-2500 K). For the direct abstraction mechanism, the overall rate constant at 300 K was predicted to be 2.58 × 10^-16 and 7.42 × 10^-25 cm³ molecule^-1s^-1 for the Br + H₂O₂ and I + H₂O₂ reactions, respectively. The modified Arrhenius parameters have been estimated for the overall reactions: k_Br+H2O2(T) = 4.80 × 10^-26 T^4.31 exp(-5.51 (kJ mol^-1)/RT) and k_I+H2O2(T) = 3.41 × 10^-23 T^3.29 exp(-56.32 (kJ mol^-1)/RT).

Accurate virial coefficients of gaseous krypton from state-of-the-art ab initio potential and polarizability of the krypton dimer

We have developed a new krypton-krypton interaction-induced isotropic dipole polarizability curve based on high-level ab initio methods. The determination was carried out using the coupled-cluster singles and doubles plus perturbative triples method with very large basis sets up to augmented correlation-consistent sextuple zeta as well as the corrections for core-core and core-valence correlation and relativistic effects. The analytical function of polarizability and our recently constructed reference interatomic potential [J. M. Waldrop et al., J. Chem. Phys. 142, 204307 (2015)] were used to predict the thermophysical and electromagnetic properties of krypton gas. The second pressure, acoustic, and dielectric virial coefficients were computed for the temperature range of 116 K-5000 K using classical statistical mechanics supplemented with high-order quantum corrections. The virial coefficients calculated were compared with the generally less precise available experimental data as well as with values computed from other potentials in the literature {in particular, the recent highly accurate potential of Jäger et al. [J. Chem. Phys. 144, 114304 (2016)]}. The detailed examination in this work suggests that the present theoretical prediction can be applied as reference values in disciplines involving thermophysical and electromagnetic properties of krypton gas.

The Lewis superacid Al[N(C6F5)2]3 and its higher homolog Ga[N(C6F5)2]3 - structural features, theoretical investigation and reactions of a metal amide with higher fluoride ion affinity than SbF5

Herein we present the synthesis of the two Lewis acids Al[N(C6F5)2]3 (ALTA) and Ga[N(C6F5)2]3 (GATA) via salt elimination reactions. The metal complexes were characterized by NMR-spectroscopic methods and X-ray diffraction analysis revealing the stabilization of the highly Lewis acidic metal centers by secondary metal-fluorine contacts. The Lewis acidic properties of Al[N(C6F5)2]3 and Ga[N(C6F5)2]3 are demonstrated by reactions with Lewis bases resulting in the formation of metallates accompanied by crucial structural changes. The two metallates [Cs(Tol)3]+[FAl(N(C6F5)2)3]- and [AsPh4]+[ClGa(N(C6F5)2)3]- contain interesting weakly coordinating anions. The reaction of Al[N(C6F5)2]3 with trityl fluoride yielded [CPh3]+[FAl(N(C6F5)2)3]- which could find application in the activation of metallocene polymerization catalysts. The qualitative Lewis acidity of Al[N(C6F5)2]3 and Ga[N(C6F5)2]3 was investigated by means of competition experiments for chloride ions in solution. DFT calculations yielded fluoride ion affinities in the gas phase (FIA) of 555 kJ mol-1 for Al[N(C6F5)2]3 and 472 kJ mol-1 for Ga[N(C6F5)2]3. Thus, Al[N(C6F5)2]3 can be considered a Lewis superacid with a fluoride affinity higher than SbF5 (493 kJ mol-1) whereas the FIA of the corresponding gallium complex is slightly below the threshold to Lewis superacidity.

NMR shielding constants in group 15 trifluorides

By combining large basis and complete basis set (CBS) extrapolations of the coupled-cluster equilibrium geometry results with rovibrational and relativistic corrections, we demonstrate that it is possible to achieve near-quantitative accuracy for the NMR shielding constants in three group 15 trifluorides – NF₃, PF₃ and AsF₃.

Potential Energy Curves, Transition Dipole Moments, and Franck-Condon Factors of the 12 Low-Lying States of BrO- Anion

This work investigates the spectroscopic parameters, vibrational levels, and transition probabilities of 12 low-lying states, which are generated from the first dissociation limit, Br(²P_u) + O^-(²P_u), of the BrO^- anion. The 12 states are X¹Σ⁺, 2¹Σ⁺, 1¹Σ^-, 1¹Π, 2¹Π, 1¹Δ, a³Π, 1³Σ⁺, 2³Σ⁺, 1³Σ^-, 2³Π, and 1³Δ. The potential energy curves are calculated with the complete active-space self-consistent field method, which is followed by the internally contracted multireference configuration interaction approach with Davidson modification. The dissociation energy D₀ of X¹Σ⁺ state is determined to be approximately 26876.44 cm^-1, which agrees well with the experimental one of 26494.50 cm^-1. Of these 12 states, the 2¹Σ⁺, 1¹Σ^-, 2¹Π, 1¹Δ, 1³Σ⁺, 2³Σ⁺, 2³Π, and 1³Δ states are very weakly bound states, whose well depths are only several-hundred cm^-1. The a³Π, 2³Π, and 1³Δ states are inverted and account for the spin-orbit coupling effect. No states are repulsive regardless of whether the spin-orbit coupling effect is included. The spectroscopic parameters and vibrational levels are determined. The transition dipole moments of 12-pair electronic states are calculated. Franck-Condon factors of a number of transitions of more than 20-pair electronic states are evaluated. The electronic transitions are discussed. The spin-orbit coupling effect on the spectroscopic parameters and vibrational properties is profound for all the states except for X¹Σ⁺, a³Π, and 1¹Π. The spectroscopic parameters and transition probabilities obtained in this paper can provide some powerful guidelines for observing these states in a proper spectroscopy experiment, in particular the states that have very shallow potential wells.