Guided Ion Beam Studies of the Thorium Monocarbonyl Cation Bond Dissociation Energy and Theoretical Unveiling of Different Isomers of [Th,O,C]+ and Their Rearrangement Mechanism

On the ground and excited electronic states of LaCO and AcCO

High-level ab initio electronic structure analysis of correlated lanthanide- and actinide-based species is laborious to perform and consequently limited in the literature. In the present work, the ground and electronically excited states of LaCO and AcCO molecules were explored utilizing the multireference configuration interaction (MRCI), Davidson corrected MRCI (MRCI+Q), and coupled cluster singles doubles and perturbative triples [CCSD(T)] quantum chemical tools conjoined with correlation consistent triple-ζ and quadruple-ζ quality all-electron Douglas-Kroll (DK) basis sets. The full potential energy curves (PECs), dissociation energies (Des), excitation energies (Tes), bond lengths (res), harmonic vibrational frequencies (ωes), and chemical bonding patterns of low-lying electronic states of LaCO and AcCO are introduced. The ground electronic state of LaCO is a 4Σ- (1σ11π2) which is a product of the reaction between excited La(4F) versus CO(X1Σ+), whereas the ground state of AcCO is a 12Π (1σ21π1) deriving from ground state fragments Ac(2D) + CO(X1Σ+). The spin-orbit ground states of LaCO (14Σ-3/2) and AcCO (12Π1/2) bear ∼13 and 5 kcal mol-1D0 values, respectively. At the MRCI level, the spin-orbit curves, the spin-orbit mixing, and the Tes of spin-orbit states of LaCO and AcCO were also analyzed. Lastly, the density functional theory (DFT) calculations were performed applying 16 exchange-correlation functionals that span three rungs of "Jacob's ladder" of density functional approximations (DFAs) to assess DFT errors associated on the De and ionization energy (IE) of LaCO.

f-Block reactions of metal cations with carbon dioxide studied by inductively coupled plasma tandem mass spectrometry

f-Block chemistry offers an opportunity to test current knowledge of chemical reactivity. The energy dependence of lanthanide cation (Ln+ = Ce+, Pr+, Nd+-Eu+) and actinide cation (An+ = Th+, U+-Am+) oxidation reactions by CO2, was observed by inductively coupled plasma tandem mass spectrometry. This reaction is commonly spin-unallowed because the neutral reactant (CO2, 1Σ+g) and product (CO, 1Σ+) require the metal and metal oxide cations to have the same spin state. Correlation of the promotion energy (Ep) to the first state with two free d-electrons with the reaction efficiency indicates that spin conservation is not a primary factor in the reaction rate. The Ep likely influences the reaction rate by partially setting the crossing between the ground and reactive states. Comparison of Ln+ and An+ congener reactivity indicates that the 5f-orbitals play a small role in the An+ reactions.

Coupled Cluster Study of the Heats of Formation of UF6 and the Uranium Oxyhalides, UO2X2 (X = F, Cl, Br, I, and At)

The atomization enthalpies of the U(VI) species UF6 and the uranium oxyhalides UO2X2 (X = F, Cl, Br, I, and At) were calculated using a composite relativistic Feller-Peterson-Dixon (FPD) approach based on scalar relativistic DKH3-CCSD(T) with extrapolations to the CBS limit. The inherent multideterminant nature of the U atom was mitigated by utilizing the singly charged atomic cation in all calculations with correction back to the neutral asymptote via the accurate ionization energy of the U atom. The effects of SO coupling were recovered using full 4-component CCSD(T) with contributions due to the Gaunt Hamiltonian calculated using Dirac-Hartree-Fock. The final atomization enthalpy for UF6 (752.2 kcal/mol) was within 2.5 kcal/mol of the experimental value, but unfortunately the latter carries a ±2.4 kcal/mol uncertainty that is predominantly due to the experimental uncertainty in the formation enthalpy of the U atom. The analogous value for UO2F2 (607.6 kcal/mol) was in nearly exact agreement with the experiment, but the latter has a stated experimental uncertainty of ±4.3 kcal/mol. The FPD atomization enthalpy for UO2Cl2 (540.4 kcal/mol) was within the experimental error limit of ±5.5 kcal/mol. FPD atomization energies for the non-U-containing molecules (used for reaction enthalpies) H2O and HX (X = F, Cl, Br, I, and At) were within at most 0.3 kcal/mol of their experimental values where available. The FPD atomization enthalpies, together with FPD reaction enthalpies for two different reactions, were used to determine heats of formation for all species of this work, with estimated uncertainties of ±4 kcal/mol. The calculated heat of formation for UF6 (-511.0 kcal/mol) is within 2.5 kcal/mol of the accurately known (±0.45 kcal/mol) experimental value.

Sequential Bond Dissociation Energies of Th+(CO)x, x = 3-6: Guided Ion Beam Collision-Induced Dissociation and Quantum Computational Studies

Collision-induced dissociation (CID) of [Th,xC,xO]⁺, x = 3-6, with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). Products are formed exclusively by the loss of CO ligands. Analyses of the kinetic energy-dependent CID product cross sections yield bond dissociation energies (BDEs) of (CO)_x-1Th⁺-CO at 0 K as 1.09 ± 0.05, 0.82 ± 0.07, 0.63 ± 0.05, and 0.70 ± 0.05 eV, respectively. Different structures of [Th,xC,xO]⁺ were explored using various electronic structure methods, and BDEs for CO ligand loss from precursor [Th,xC,xO]⁺ complexes were computed. Both experimental and theoretical results corroborate that the structures of [Th,xC,xO]⁺, x = 3-6, formed experimentally are homoleptic thorium cation carbonyl complexes, Th⁺(CO)_x. The nonmonotonic trend in experimental BDEs is reproduced theoretically, although ambiguities in the spin states of the x = 4-6 complexes (doublet or quartet) remain. BDEs calculated at the coupled cluster with single, double, and perturbative triple excitations (CCSD(T))/cc-pVXZ//B3LYP/cc-PVXZ (X = T and Q) level and a complete basis set (CBS) extrapolation agree reasonably well with the experimental values for all complexes. Thorium oxide ketenylidene carbonyl cations, OTh⁺CCO(CO)_y, y = 1-4, were calculated to be the most stable structures of [Th,xC,xO]⁺, x = 3-6, respectively; however, these are not observed in our experiment. Potential energy profiles (PEPs) having either quartet or doublet spin calculated at the B3LYP/cc-pVQZ level suggest that the failure to observe OTh⁺CCO(CO)_y, y = 1-4, is the result of a barrier corresponding to the C-C bond formation, making the formation of OTh⁺CCO(CO)_y inaccessible kinetically under the present experimental conditions.

Interaction of Th with H0/-/+: Combined Experimental and Theoretical Thermodynamic Properties

High-level electronic structure calculations of the low-lying energy electronic states for ThH, ThH^-, and ThH⁺ are reported and compared to experimental measurements. The inclusion of spin-orbit coupling is critical to predict the ground-state ordering as inclusion of spin-orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH⁺. At the multireference spin-orbit SO-CASPT2 level, the ground states of ThH, ThH^-, and ThH⁺ are predicted to be the ²Δ_3/2, ³Φ₂, and ³Δ₁ states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH^- photoelectron spectrum has many transitions, which approximately correlate with excitations of Th⁺ and/or Th. The adiabatic ionization energy of ThH including spin-orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th⁺H^- ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH^-, and ThH⁺ using the Feller-Peterson-Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.

Experimental and computational investigation of the bond energy of thorium dicarbonyl cation and theoretical elucidation of its isomerization mechanism to the thermodynamically most stable isomer, thorium oxide ketenylidene cation, OTh+CCO

Collision-induced dissociation (CID) of [Th,2C,2O]⁺ with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). The only products observed are ThCO⁺ and Th⁺ by sequential loss of CO ligands. The experimental findings and theoretical calculations support that the structure of [Th,2C,2O]⁺ is the bent homoleptic thorium dicarbonyl cation, Th⁺(CO)₂, having quartet spin, which is both thermodynamically and kinetically stable enough in the gas phase to be observed in our GIBMS instrument. Analysis of the kinetic energy-dependent cross sections for this CID reaction yields the first experimental determination of the bond dissociation energy (BDE) of (CO)Th⁺-CO at 0 K as 1.05 ± 0.09 eV. A theoretical BDE calculated at the CCSD(T) level with cc-pVXZ (X = T and Q) basis sets and a complete basis set (CBS) extrapolation is in very good agreement with the experimental result. Although the doublet spin bent thorium oxide ketenylidene cation, OTh⁺CCO, is calculated to be the most thermodynamically stable structure, it is not observed in our experiment where [Th,2C,2O]⁺ is formed by association of Th⁺ and CO in a direct current discharge flow tube (DC/FT) ion source. Potential energy profiles of both quartet and doublet spin are constructed to elucidate the isomerization mechanism of Th⁺(CO)₂ to OTh⁺CCO. The failure to observe OTh⁺CCO is attributed to a barrier associated with C-C bond formation, which makes OTh⁺CCO kinetically inaccessible under our experimental conditions. Chemical bonding patterns in low-lying states of linear and bent Th⁺(CO)₂ and OTh⁺CCO isomers are also investigated.