Direct determination of the ionization energies of FeO and CuO with VUV radiation

Sub-picosecond photodynamics of small neutral copper oxide clusters

The ultrafast dynamics of neutral copper oxide clusters (CunOx, n < 5) are reported using femtosecond pump probe spectroscopy in the gas phase. The transient spectra recorded for each cluster demonstrate they relax on a 100s of fs timescale followed by a long-lived (>50 ps) response. Density functional theory calculations are performed to determine the lowest energy structures and spin states. Topological descripters for the excited states are calculated (time-dependent density functional theory) to relate the measured excited state dynamics to changes in the cluster's electronic structure with increasing oxidation. Strong field ionization is demonstrated here to be a soft form of ionization and able to record transient signals for clusters previously determined to be unstable to nanosecond multiphoton ionization. The relative cluster stability is further demonstrated by signal enhancement/depreciation that is recorded through the synergy from the two laser pulses. Once the oxygen atoms exceed the number of copper atoms, a weakly bound superoxide O2 unit forms, exhibiting a higher spin state. All clusters that are not in the lowest spin configuration demonstrate fragmentation.

M11pz: A Nonlocal Meta Functional with Zero Hartree-Fock Exchange and with Broad Accuracy for Chemical Energies and Structures

The accuracy of Kohn-Sham density functional theory depends strongly on the approximation to the exchange-correlation functional. In this work, we present a new exchange-correlation functional called M11pz (M11 plus rung-3.5 terms with zero Hartree-Fock exchange) that is built on the M11plus functional with the goal of using its rung-3.5 terms without a Hartree-Fock exchange term, especially to improve the accuracy for strongly correlated systems. The M11pz functional is optimized with the same local and rung-3.5 ingredients that are used in M11plus but without any percentage of Hartree-Fock exchange. The performance of M11pz is compared with eight local functionals, and M11pz is found to be in top three when the errors or ranks are averaged over eight grouped and partially overlapping databases: AME418/22, atomic and molecular energies; MGBE172, main-group bond energies; TMBE40, transition-metal bond energies; SR309, single-reference systems; MR54, multireference systems; BH192, barrier heights; NC579, noncovalent interaction energies; and MS20, molecular structures. For calculations of band gaps of solids, M11pz is the second best of the nine tested functionals that have zero Hartree-Fock exchange.

Bonding, Thermodynamics, and Dissociation Dynamics of NiO+ and NiS+ Determined by Photofragment Imaging and Theory

We use photofragment ion imaging and ab initio calculations to determine the bond strength and photodissociation dynamics of the nickel oxide (NiO⁺) and nickel sulfide (NiS⁺) cations. NiO⁺ photodissociates broadly from 20350 to 32000 cm^-1, forming ground state products Ni⁺(²D) + O(³P) below ∼29000 cm^-1. Above this energy, Ni⁺(⁴F) + O(³P) products become accessible and dominate over the ground state channel. In certain images, product spin-orbit levels are resolved, and spin-orbit propensities are determined. Image anisotropy and the results of MRCI calculations suggest NiO⁺ photodissociates via a 3 ⁴Σ^- ← X ⁴Σ^- transition above the Ni⁺(⁴F) threshold and via 3 ⁴Σ^-, 2 ⁴Σ^-, and/or 2 ⁴Π and 3 ⁴Π excited states below the ⁴F threshold. The photodissociation spectrum of NiS⁺ from 19900 to 23200 cm^-1 is highly structured, with ∼12 distinct vibronic peaks, each containing underlying substructure. Above 21600 cm^-1, the Ni⁺(²D_5/2) + S(³P) and Ni⁺(²D_3/2) + S(³P) product spin-orbit channels compete, with a branching ratio of ∼2:1. At lower energy, Ni⁺(²D_5/2) is formed exclusively, and S(³P₂) and S(³P₁) spin-orbit channels are resolved. MRCI calculations predict the ground state of NiS⁺ to be one of two nearly degenerate states, the 1 ⁴Σ^- and 1 ⁴Δ. Based on images and spectra, the ground state of NiS⁺ is assigned as ⁴Δ_7/2, with the 1 ⁴Σ_3/2^- and 1 ⁴Σ_1/2^- states 81 ± 30 and 166 ± 50 cm^-1 higher in energy, respectively. The majority of the photodissociation spectrum is assigned to transitions from the 1 ⁴Δ state to two overlapping, predissociative excited ⁴Δ states. Our D₀ measurements for NiO⁺ (D₀ = 244.6 ± 2.4 kJ/mol) and NiS⁺ (D₀ = 240.3 ± 1.4 kJ/mol) are more precise and closer to each other than previously reported values. Finally, using a recent measurement of D₀(NiS), we derive a more precise value for IE (NiS): 8.80 ± 0.02 eV (849 ± 1.7 kJ/mol).

Ultrafast pump-probe spectroscopy of neutral FenOm clusters (n, m < 16)

Neutral iron oxide clusters (FenOm, n, m ≤ 16) are produced in a laser vaporization source using O2 gas seeded in He. The neutral clusters are ionized with a sequence of femtosecond laser pulses and detected using time-of-flight mass spectrometry. Small clusters are confirmed to be most prominent in the stoichiometric (n = m) distribution, with m = n + 1 clusters observed above n = 4. Pump-probe spectroscopy is employed to study the dynamics of an electron transfer from an oxygen orbital to an iron nonbonding orbital of iron oxide clusters that is driven by absorption of a 400 nm photon. A bifurcation of the initial wavepacket occurs, where a femtosecond component is attributed to electron relaxation assisted through internuclear vibrational relaxation and high density of states, and a slow relaxation shows the formation of a bound excited state. The lifetime and relative ratio of the two pathways depend on both the cluster size and iron oxidation state. The femtosecond lifetime decreases with increased cluster size until a saturation timescale is achieved at n > 5. The relative population of the long-lived excited state decreases with cluster size and suggests that the excited electron remains on the Fe atom for >20 ps.

Guided Ion Beam Tandem Mass Spectrometry and Theoretical Study of SO2 Activated by Os. J Phys Chem A 2020;124:6629-6644. [PMID: 32702982 DOI: 10.1021/acs.jpca.0c05757]

The kinetic-energy dependence of SO₂ activated by Os⁺ was studied by guided ion beam tandem mass spectrometry. Species observed in endothermic reactions were OsO⁺, OsO₂⁺ or OsS⁺, and OOsS⁺. The kinetic energy-dependent cross sections were modeled to yield 0 K bond dissociation energies (BDEs) of 5.01 ± 0.06 eV (Os⁺-O), 5.15 ± 0.07 eV (Os⁺-O₂), 4.50 ± 0.17 eV (Os⁺-S), and 4.22 ± 0.11 eV (Os⁺-SO). Among these BDE values, the values for OsO⁺ and OsO₂⁺ agree with literature values and those for OsS⁺ and OOsS⁺ are novel measurements. Theoretical calculations were performed at a B3LYP/def2-TZVPPD level for all products, and additional calculations were performed for OsS⁺, OsO₂⁺, and OsSO⁺ using the CCSD(T) level of theory, extrapolated to the complete basis set (CBS) limit, and def2-QZVPPD and aug-cc-pVxZ (x = T, Q, and 5) basis sets. These calculations indicate that the ground states of the products are ⁴Π_5/2 (OsO⁺), ²B₁ (OsO₂⁺), ⁴Π_5/2 (OsS⁺), and ²A″ (OOsS⁺) after including empirical spin-orbit corrections. The potential energy surfaces (PESs) for OsSO₂⁺ intermediates, transition states, and all products were also investigated at the B3LYP/def2-TZVPPD level. The PESs show that none of the reactions have barriers in excess of the product endothermicities. Cross sections for OsO⁺ formation are compared to those from previous guided ion beam studies of related systems (Os⁺ + O₂ and CO and Re⁺ + SO₂) to evaluate their relative behaviors.

Ab initio investigation of the ground and excited states of RuO+,0,- and their reaction with water

High-level quantum chemical calculations on RuO0,± elucidate the electronic structure of their low-lying electronic states. For thirty-two states, we report the electronic configurations, bond lengths, vibrational frequencies, spin-orbit splittings, and excitation energies. The electronic states of RuO can be generated from those of RuO+ by adding one electron to the σ non-bonding orbital closely resembling the 5s atomic orbital of Ru. The ground states for RuO and RuO- are clearly identified as 5Δ and 4Δ, but the two states (4Δ and 2Π) compete for RuO+. The difficulty of calculations is revealed by our small binding energies compared to the experimental values. In addition, we studied the reaction of the three species with water in their ground and selected low-lying electronic states. We found a consistent decrease of the activation energy barriers and higher exothermicity as we add electrons to the system. RuO- is found to facilitate the reaction for both kinetic and thermodynamic reasons.

Analytical interatomic bond-order potential for simulations of oxygen defects in iron

We present an analytical bond-order potential for the Fe-O system, capable of reproducing the basic properties of wüstite as well as the energetics of oxygen impurities in [Formula: see text]-iron. The potential predicts binding energies of various small oxygen-vacancy clusters in [Formula: see text]-iron in good agreement with density functional theory results, and is therefore suitable for simulations of oxygen-based defects in iron. We apply the potential in simulations of the stability and structure of Fe/FeO interfaces and FeO precipitates in iron, and observe that the shape of FeO precipitates can change due to formation of well-defined Fe/FeO interfaces. The interface with crystalline Fe also ensures that the precipitates never become fully amorphous, no matter how small they are.

Functionalization of the transition metal oxides FeO, CoO, and NiO with alkali metal atoms decreases their ionization potentials by 3-5 eV

The existence and stabilities of various neutral metal oxides of formula MON and MON₂ (M = Fe, Co, Ni; N = Li, Na) and their corresponding cations MON⁺ and MON₂⁺ are predicted using density functional theory (B3LYP) with the 6-311 + G(d) basis set. Ab initio calculations carried out at the CCSD(T)/6-311 + G(3df) level of theory reveal that the ionization potentials (IPs) of the oxides MO decrease by ca. 3-5 eV upon functionalization with N to give either MON or MON₂. The influences of the chemical constitution and local spin magnetic moment (on the transition metal atom) of the oxide or cation on its IP are presented and discussed.

Reactivity of 4Fe+(CO)n=0-2 + O2: oxidation of CO by O2 at an isolated metal atom

The kinetics of ⁴Fe⁺(CO)_n=0-2 + O₂ are measured under thermal conditions from 300-600 K using a selected-ion flow tube apparatus. Both the bare metal and n = 2 cations are inert to reaction over this temperature range, but ⁴Fe⁺(CO) reacts rapidly (k = 3.2 ± 0.8 × 10^-10 cm³ s^-1 at 300 K, 52% of the collisional rate coefficient) to form FeO⁺ + CO₂. This is an example of the oxidation of CO by O₂ occurring entirely on a single non-noble metal atom. The reaction of the bare metal reaction is known to be endothermic, such that this result is expected; however, the n = 2 reaction has highly exothermic product channels available, such that the lack of reaction is surprising in light of the n = 1 reactivity. Stationary points along all three reaction coordinates are calculated using the TPSSh hybrid functional. These surfaces show that the n = 1 reaction is an example of two-state reactivity; the reaction proceeds initially on the sextet surface over a submerged barrier to a structure with an O-O bond distance longer than that in O₂, but must cross to the quartet surface in order to proceed over a second submerged barrier to rearrange to form CO₂. The n = 2 reaction does not proceed because, on all spin surfaces, the transition state corresponding to O-O separation is at higher energy than the separated reactants. The difference between the n = 1 and n = 2 reactions is not a result of steric effects, but rather because the O₂ is more strongly bound to Fe in the entrance well of the n = 1 case, and that energy is available to overcome the rate-limiting barrier to O-O cleavage. Experimental verification of some of these details are provided by guided ion beam tandem mass spectrometry results. The kinetic energy dependence of the n = 1 reaction shows evidence for a curve crossing and yields relevant thermochemistry for competing reaction channels.

Reactivity from excited state (4)FeO(+) + CO sampled through reaction of ground state (4)FeCO(+) + N2O

The kinetics of the FeCO(+) + N2O reaction have been studied at thermal energies (300-600 K) using a variable temperature selected ion flow tube apparatus. Rate constants and product branching fractions are reported. The reaction is modestly inefficient, proceeding with a rate constant of 6.2 × 10(-11) cm(3) s(-1) at 300 K, with a small negative temperature dependence, declining to 4.4 × 10(-11) cm(3) s(-1) at 600 K. Both Fe(+) and FeO(+) products are observed, with a constant branching ratio of approximately 40:60 at all temperatures. Calculation of the stationary points along the reaction coordinate shows that only the ground state quartet surface is initially sampled resulting in N2 elimination; a submerged barrier along this portion of the surface dictates the magnitude and temperature dependence of the total rate constant. The product branching fractions are determined by the behavior of the remaining (4)OFeCO(+) fragment, and this behavior is compared to that found in the reaction of FeO(+) + CO, which initially forms (6)OFeCO(+). Thermodynamic and kinetic arguments are used to show that the spin-forbidden surface crossing in this region is efficient, proceeding with an average rate constant of greater than 10(12) s(-1).

The power of accurate energetics (or thermochemistry: what is it good for?)

The utility of measuring the energetics of ion-molecule reactions is discussed. After distinguishing between the terms of thermodynamics (macroscopic, equilibrium quantities) and energetics (microscopic and kinetically relevant quantities), the potential energy surfaces for ion-molecule reactions are reviewed and their implications discussed. Equations describing the kinetic energy dependence of ion-molecule reactions are introduced and the effects of entropy on reaction rates and branching ratios are discussed. Several case histories allow an exploration of the utility of accurate thermochemical information and probe how accurate such energetic information must be to be predictive. These case studies include decomposition of hydrated metal dications, the reaction of FeO(+) with H(2), and fragmentation of a small protonated peptide (GG). These illustrate a range of interesting systems for which accurate energetic information has been influential in understanding the observed reactivity. Comparisons with theory demonstrate that experimental information is still required for truly predictive capability.

Communication: A vibrational study of titanium dioxide cation using the vacuum ultraviolet laser pulsed field ionization-photoelectron method

We have successfully measured the vacuum ultraviolet (VUV) laser photoionization efficiency and pulsed field ionization-photoelectron (PFI-PE) spectra of cold titanium dioxide (TiO(2)) prepared by a supersonically cooled laser ablation source. The VUV-PFI-PE spectrum thus obtained exhibits long progressions of the v(2)(+)(a(1)) bending and the combination of v(1)(+)(a(1)) stretching and v(2)(+)(a(1)) bending vibrational modes of the TiO(2)(+)(X(2)B(2)) ion. The pattern of Franck-Condon factors observed indicates that the O-Ti-O bond angle of the TiO(2)(+)(X(2)B(2)) ion is significantly different from that of the TiO(2)(X(1)A(1)) neutral, whereas the change of the Ti-O bond distance is very minor upon the photoionization transition. The analysis of the PFI-PE bands has made possible the determination of the adiabatic ionization energy for TiO(2), IE(TiO(2)) = 77215.9 ± 1.2 cm(-1) (9.57355 ± 0.00015 eV), the harmonic vibrational frequencies, ω(1)(+) = 829.1 ± 2.0 cm(-1) and ω(2)(+) = 248.7 ± 0.6 cm(-1), and the anharmonic coefficients, χ(11)(+) = 5.57 ± 0.65 cm(-1), χ(22)(+) = 0.08 ± 0.06 cm(-1), and χ(12)(+) = -4.51 ± 0.30 cm(-1) for the TiO(2)(+)(X(2)B(2)) ground state.

Density functional calculation of the structure and electronic properties of Cu(n)O(n) (n = 1-8) clusters

Ab initio simulations and calculations were used to study the structures and stabilities of copper oxide clusters, Cu(n)O(n) (n = 1-8). The lowest energy structures of neutral and charged copper oxide clusters were determined using primarily the B3LYP/LANL2DZ model chemistry. For n ≥ 4, the clusters are nonplanar. Selected electronic properties including atomization energies, ionization energies, electron affinities, and Bader charges were calculated and examined as a function of n.

Energy-dependent dissociation of benzylpyridinium ions in an ion-trap mass spectrometer

Benzylpyridinium ions, generated via electrospray ionization of dilute solutions of their salts in acetonitrile/water, are probed by collisional activation in an ion-trap mass spectrometer. From the breakdown diagrams obtained, phenomenological appearance energies of the fragment ions are derived. Comparison of the appearance energies with calculated reaction endothermicities shows a reasonably good correlation for this particular class of compounds. In addition, the data indirectly indicate that at threshold the dissociation of almost all of the benzylpyridinium ions under study leads to the corresponding benzylium ions, rather than the tropylium isomers. Substituent effects on the fragmentation for a series of benzylpyridinium ions demonstrate that neither mass effects nor differences in density of states seriously affect the energetics derived from the ion-trap experiments.

Vibrational spectroscopy of intermediates in benzene-to-pheno conversion by FeO+

Gas-phase FeO(+) can convert benzene to phenol under thermal conditions. Two key intermediates of this reaction are the [HO-Fe-C(6)H(5)](+) insertion intermediate and Fe(+)(C(6)H(5)OH) exit channel complex. These intermediates are selectively formed by reaction of laser ablated Fe(+) with specific organic precursors and are cooled in a supersonic expansion. Vibrational spectra of the sextet and quartet states of the intermediates in the O-H stretching region are measured by infrared multiphoton dissociation (IRMPD). For Fe(+)(C(6)H(5)OH), the O-H stretch is observed at 3598 cm(-1). Photodissociation primarily produces Fe(+) + C(6)H(5)OH; Fe(+)(C(6)H(4)) + H(2)O is also observed. IRMPD of [HO-Fe-C(6)H(5)](+) mainly produces FeOH(+) + C(6)H(5) and the O-H stretch spectrum consists of a peak at approximately 3700 cm(-1) with a shoulder at approximately 3670 cm(-1). Analysis of the experimental results is aided by comparison with hybrid density functional theory computed frequencies. Also, an improved potential energy surface for the FeO(+) + C(6)H(6) reaction is developed based on CBS-QB3 calculations for the reactants, intermediates, transition states, and products.

Chemical dynamics, molecular energetics, and kinetics at the synchrotron

Scientists at the Chemical Dynamics Beamline of the Advanced Light Source in Berkeley are continuously reinventing synchrotron investigations of physical chemistry and chemical physics with vacuum ultraviolet light. One of the unique aspects of a synchrotron for chemical physics research is the widely tunable vacuum ultraviolet light that permits threshold ionization of large molecules with minimal fragmentation. This provides novel opportunities to assess molecular energetics and reaction mechanisms, even beyond simple gas phase molecules. In this perspective, significant new directions utilizing the capabilities at the Chemical Dynamics Beamline are presented, along with an outlook for future synchrotron and free electron laser science in chemical dynamics. Among the established and emerging fields of investigations are cluster and biological molecule spectroscopy and structure, combustion flame chemistry mechanisms, radical kinetics and product isomer dynamics, aerosol heterogeneous chemistry, planetary and interstellar chemistry, and secondary neutral ion-beam desorption imaging of biological matter and materials chemistry.

A Combined Experimental and Computational Study on the Ionization Energies of the Cyclic and Linear C3H Isomers

For the first time, two hydrogen-deficient hydrocarbon radicals are generated in situ via laser ablation of graphite and seeding the ablated species in acetylene gas, which acts as a carrier and reactant simultaneously. By recording photoionization efficiency curves (PIE) and simulating the experimental spectrum with computed Franck-Condon (FC) factors, we can reproduce the general pattern of the PIE curve of m/z=37. We recover ionization energies of 9.15 eV and 9.76 eV for the linear and cyclic isomers, respectively. Our combined experimental and theoretical studies provide an unprecedented, versatile pathway to investigate the ionization energies of even more complex hydrocarbon radicals in situ, which are difficult to prepare by classical synthesis, in future experiments.

Degradation of Ionized OV(OCH3)3 in the Gas Phase. From the Neutral Compound All the Way down to the Quasi-terminal Fragments VO+ and VOH+

The consecutive fragmentation of ionized trimethyl vanadate(V), OV(OCH3)3 (1), is examined by experiment and theory. After an elimination of formaldehyde from the molecular ion 1+, subsequent dissociations proceed via losses of first H2 and then two molecules of formaldehyde to finally yield the VOH+ cation; these redox reactions involve the V(II)/V(IV) manifold. At elevated energies, expulsion of CH3O* from 1+ can efficiently compete to afford OV(OCH3)2+, a formal V(V) compound, from which subsequent losses of H2 and two units of CH2O lead to bare VO+, thereby exploring the V(III)/V(V) redox manifold. Experiments using complementary mass spectrometric techniques, i.e., neutralization-reionization experiments and ion/molecule reactions, in conjunction with extensive computational studies provide deep insight into the ion structures and the relative energetics of these dissociation reactions. In particular, a quantitative energetic scheme is obtained that ranges from neutral OV(OCH3)3 all the way down to the quasi-terminal fragment ions VOH+ and VO+, respectively.

A kinetic study of the reactions FeO++ O, Fe+·N2+ O, Fe+·O2+ O and FeO++ CO: implications for sporadic E layers in the upper atmosphere

These gas-phase reactions were studied by pulsed laser ablation of an iron target to produce Fe(+) in a fast flow tube, with detection of the ions by quadrupole mass spectrometry. Fe(+).N(2) and Fe(+).O(2) were produced by injecting N(2) and O(2), respectively, into the flow tube. FeO(+) was produced from Fe(+) by addition of N(2)O, or by ligand-switching from Fe(+).N(2) following the addition of atomic O. The following rate coefficients were measured: k(FeO(+) + O --> Fe(+) + O(2), 186-294 K) = (3.2 +/- 1.5) x 10(-11); k(Fe(+).N(2) + O --> FeO(+)+ N(2), 294 K) = (4.6 +/- 2.5) x 10(-10); k(Fe(+).O(2) + O --> FeO(+) + O(2), 294 K) = (6.3 +/- 2.7) x 10(-11); and k(FeO(+) + CO --> Fe(+) + CO(2), 294 K) = (1.59 +/- 0.34) x 10(-10) cm(3) molecule(-1) s(-1), where the quoted uncertainties are a combination of the 1sigma standard errors in the kinetic data and the systematic experimental errors. The surprisingly slow reaction between FeO(+) and O is examined using ab initio quantum calculations of the relevant potential energy surfaces. The importance of this reaction for controlling the lifetime of sporadic E layers is then demonstrated using a model of the upper mesosphere and lower thermosphere.