Generation of Naphthoquinone Radical Anions by Electrospray Ionization: Solution, Gas-Phase, and Computational Chemistry Studies

Radical anions are present in several chemical processes, and understanding the reactivity of these species may be described by their thermodynamic properties. Over the last years, the formation of radical ions in the gas phase has been an important issue concerning electrospray ionization mass spectrometry studies. In this work, we report on the generation of radical anions of quinonoid compounds (Q) by electrospray ionization mass spectrometry. The balance between radical anion formation and the deprotonated molecule is also analyzed by influence of the experimental parameters (gas-phase acidity, electron affinity, and reduction potential) and solvent system employed. The gas-phase parameters for formation of radical species and deprotonated species were achieved on the basis of computational thermochemistry. The solution effects on the formation of radical anion (Q(center dot-)) and dianion (Q(2-)) were evaluated on the basis of cyclic voltammetry analysis and the reduction potentials compared with calculated electron affinities. The occurrence of unexpected ions [Q + 15](-) was described as being a reaction between the solvent system and the radical anion, Q(center dot-).The gas-phase chemistry of the electrosprayed radical anions was obtained by collisional-induced dissociation and compared to the relative energy calculations. These results are important for understanding the formation and reactivity of radical anions and to establish their correlation with the reducing properties by electrospray ionization analyses.

New Relativistic Atomic Natural Orbital Basis Sets for Lanthanide Atoms with Applications to the Ce Diatom and LuF3

New basis sets of the atomic natural orbital (ANO) type have been developed for the lanthanide atoms La-Lu. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies and some excitation energies. Computed ionization energies have an accuracy better than 0.1 eV in most cases. Two molecular applications are inluded as illustration: the cerium diatom and the LuF3 molecule. In both cases it is shown that 4f orbitals are not involved in the chemical bond in contrast to an earlier claim for the latter molecule.

Predicting and characterizing new molecular species at a high-level theoretical approach: HSeCl, HClSe, and SeCl

New molecular species HSeCl, HClSe, and SeCl were investigated at a high level of theory, CCSD(T), with a series of correlation consistent basis sets with extrapolation to the CBS limit. Account has been taken for valence-only and core-valence correlation effects, and of anharmonic effects on the vibrational frequencies. HSeCl is 43.25 kcal mol (1) more stable than HClSe. A barrier (Delta G(#)) of 47.20 kcal mol (1) separates these species. Internuclear distances are generally overestimated by 0.008 angstrom in the valence-only correlation calculations. Inclusion of anharmonicity leads to much improved vibrational frequencies. For SeCl, we estimate Delta H(f) (0 K) = 23.96 and Delta H(f) (298.15 K) = 24.64 kcal mol (1); for HSeCl, we had 4.20 and 4.97 kcal mol (1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

Exploring new molecular species on the (1)[H, Se, F] singlet potential energy surface: Energetics, structures, IR spectra, and heats of formation

Structural, energetic, and vibrational properties of new molecular species, HSeF and HFSe, the associated transition state, and dissociation fragments are investigated using a state-of-the-art theoretical approach, CCSD(T)/CBS. HSeF is a normal covalently bonded molecule 38.98 kcal mol (1) more stable than the complex HF-Se, which shows an unusual structure with a central fluorine atom and a bond angle of 101.8 degrees.A barrier (Delta G(#)) of 49.01 kcal mol (1) separates the two species. Vibrational frequencies are also quite distinct. Heats of formation are evaluated for the diatomic fragments and HSeF. Final Delta(f)H values depend on the experimental accuracy of those of Se(g) and H(2)Se. (c) 2009 Elsevier B.V. All rights reserved.

Electronic structure and spectra of a new molecular species: SI. A theoretical contribution

A high level theoretical approach is used to characterize for the first time a manifold of doublet and quartet A + S and Omega states correlating with the first two dissociation channels of an as yet experimentally unknown molecular species, SI, sulfur monoidide. A set of spectroscopic constants is determined, including vibrationally averaged spin-orbit coupling constants, vibrationally averaged dipole moments, and dissociation energies. The transition dipole moment function for the spin-forbidden transition a (4)Sigma -X (2)Pi, and the associated radiative lifetimes were also evaluated. Two possibilities to detect transitions experimentally and to derive spectroscopic constants are suggested. (C) 2011 Elsevier B. V. All rights reserved.

New Molecular Species of Potential Interest to Atmospheric Chemistry: Isomers on the [H, S(2), Br] Potential Energy Surface

This work reports a state-of-the-art theoretical characterization of four new sulfur-bromine species and five transition states on the [H, S(2), Br] potential energy surface. Our highest level theoretical approach employed the method coupled cluster singles and doubles with perturbative contributions of connected triples, CCSD(T), along with the series of correlation-consistent basis sets and with extrapolation to the complete basis set (CBS) limit in the optimization of the geometrical parameters and to quantify the energetic quantities. The structural and vibrational frequencies here reported are unique and represent the most accurate investigation to date of these species. The global minimum corresponds to a skewed structure HSSBr with a disulfide bond; this is followed by a pyramidal-like structure, SSHBr, 18.85 kcal/mol above the minimum. Much higher in energy, we found another skewed structure, HSBrS (50.29 kcal/mol), with one S-Br dative-type bond, and another pyramidal-like one, HBrSS (109.80 kcal/mol), with two S-Br dative-type bonds. The interconversion of HSSBr into SSHBr can occur via a transfer of either the hydrogen or the bromine atom but involves a very high barrier of about 43 kcal/mol. These molecules are potentially a new route of coupling the sulfur and bromine chemistry in the atmosphere, and conditions of high concentration of H(2)S like in volcanic eruptions might contribute to their formation. We note that HSSBr can act as a reservoir molecule for the reaction between the radicals HSS and Br. Also, an assessment of the methods DFT/B3LYP/CBS and MP2/CBS relative to CCSD(T)/CBS provides insights on the expected performance of these methods on the characterization of polysulfides and also of more complex systems containing disulfide bridges.

Nonempirical hyper-generalized-gradient functionals constructed from the Lieb-Oxford bound

A simple and completely general representation of the exact exchange-correlation functional of density-functional theory is derived from the universal Lieb-Oxford bound, which holds for any Coulomb-interacting system. This representation leads to an alternative point of view on popular hybrid functionals, providing a rationale for why they work and how they can be constructed. A similar representation of the exact correlation functional allows to construct fully nonempirical hyper-generalized-gradient approximations (HGGAs), radically departing from established paradigms of functional construction. Numerical tests of these HGGAs for atomic and molecular correlation energies and molecular atomization energies show that even simple HGGAs match or outperform state-of-the-art correlation functionals currently used in solid-state physics and quantum chemistry.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.

Metastable BrO(2+) and NBr(2+) molecules in the gas phase

The doubly positively charged gas-phase molecules BrO(2+) and NBr(2+) have been produced by prolonged high-current energetic oxygen (17 keV (16)O(-)) ion surface bombardment (ion beam sputtering) of rubidium bromide (RbBr) and of ammonium bromide (NH(4)Br) powdered ionic salt samples, respectively, pressed into indium foil. These novel species were observed at half-integer m/z values in positive ion mass spectra for ion flight times of roughly similar to 12 mu s through a magnetic-sector secondary ion mass spectrometer. Here we present these experimental results and combine them with a detailed theoretical investigation using high level ab initio calculations of the ground states of BrO(2+) and NBr(2+), and a manifold of excited electronic states. NBr(2+) and BrO(2+), in their ground states, are long-lived metastable gas-phase molecules with well depths of 2.73 x 10(4) cm(-1) (3.38 eV) and 1.62 x 10(4) cm(-1) (2.01 eV); their fragmentation channels into two monocations lie 2.31 x 10(3) cm(-1) (0.29 eV) and 2.14 x 10(4) cm(-1) (2.65 eV) below the ground state minimum. The calculated lifetimes for NBr(2+) (v '' < 35) and BrO(2+) (v '' < 18) are large enough to be considered stable against tunneling. For NBr(2+), we predicted R(e) = 3.051 a(0) and omega(e) = 984 cm(-1); for BrO(2+), we obtained 3.033 a(0) and 916 cm(-1), respectively. The adiabatic double ionization energies of BrO and NBr to form metastable BrO(2+) and NBr(2+) are calculated to be 30.73 and 29.08 eV, respectively. The effect of spin-orbit interactions on the low-lying (Lambda + S) states is also discussed. (C) 2011 American Institute of Physics. [doi:10.1063/1.3562121]

The electronic states of SeF: A reinterpretation of the chemiluminescent emission of the reaction of selenium with fluorine

The low-lying doublet and quartet electronic states of the species SeF correlating with the first dissociation channel are investigated theoretically at a high-level of electronic correlation treatment, namely, the complete active space self-consistent field/multireference single and double excitations configuration interaction (CASSCF/MRSDCI) using a quintuple-zeta quality basis set including a relativistic effective core potential for the selenium atom. Potential energy curves for (Lambda+S) states and the corresponding spectroscopic properties are derived that allows for an unambiguous assignment of the only spectrum known experimentally as due to a spin-forbidden X (2)Pi-a (4)Sigma(-) transition, and not a A (2)Pi-X (2)Pi transition as assumed so far. For the bound excited doublets, yet unknown experimentally, this study is the first theoretical characterization of their spectroscopic properties. Also the spin-orbit coupling constant function for the X (2)Pi state is derived as well as the spin-orbit coupling matrix element between the X (2)Pi and a (4)Sigma(-) states. Dipole moment functions and vibrationally averaged dipole moments show SeF to be a very polar species. An overview of the lowest-lying spin-orbit (Omega) states completes this description. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3426315]

Excitation energies from ground-state density-functionals by means of generator coordinates

The generator-coordinate method is a flexible and powerful reformulation of the variational principle. Here we show that by introducing a generator coordinate in the Kohn-Sham equation of density-functional theory, excitation energies can be obtained from ground-state density functionals. As a viability test, the method is applied to ground-state energies and various types of excited-state energies of atoms and ions from the He and the Li isoelectronic series. Results are compared to a variety of alternative DFT-based approaches to excited states, in particular time-dependent density-functional theory with exact and approximate potentials.

CASPT2 Study of the Potential Energy Surface of the HSO(2) System

The importance of the HSO(2) system in atmospheric and combustion chemistry has motivated several works dedicated to the study of associated structures and chemical reactions. Nevertheless controversy still exists in connection with the reaction SH + O(2) -> H + SO(2) and also related to the role of the HSOO isomers in the potential energy surface (PES). Here we report high-level ab initio calculation for the electronic ground state of the HSO(2) system. Energetic, geometric, and frequency properties for the major stationary states of the PES are reported at the same level of calculations:,CASPT2/aug-cc-pV(T+d)Z. This study introduces three new stationary points (two saddle points and one minimum). These structures allow the connection of the skewed HSOOs and the HSO(2) minima defining new reaction paths for SH + O(2) -> H + SO(2) and SH + O(2) -> OH + SO. In addition, the location of the HSOO isomers in the reaction pathways have been clarified.

Empirical analysis of the Lieb-Oxford bound in ions and molecules

Universal properties of the Coulomb interaction energy apply to all many-electron systems. Bounds on the exchange-correlation energy, in particular, are important for the construction of improved density functionals. Here we investigate one such universal property-the Lieb-Oxford lower bound-for ionic and molecular systems. In recent work [J Chem Phys 127, 054106 (2007)], we observed that for atoms and electron liquids this bound may be substantially tightened. Calculations for a few ions and molecules suggested the same tendency, but were not conclusive due to the small number of systems considered. Here we extend that analysis to many different families of ions and molecules, and find that for these, too, the bound can be empirically tightened by a similar margin as for atoms and electron liquids. Tightening the Lieb-Oxford bound will have consequences for the performance of various approximate exchange-correlation functionals. (C) 2008 Wiley Periodicals Inc.

Structural and spectroscopic properties of the diazocarbene radical (CNN) and its ions CNN(+) and CNN(-): a high-level theoretical investigation

The diazocarbene radical, CNN, and the ions CNN(+) and CNN(-) were investigated at a high level of theory. Very accurate structural parameters for the states X(3)Sigma(-) and A(3)Pi of CNN, and X(2)Pi of both CNN(+) and CNN(-) were obtained with the UCCSD(T) method using correlated-consistent basis functions with extrapolations to the complete basis set limit, with valence only and also with all electrons correlated. Harmonic and anharmonic frequencies were obtained for all species and the Renner parameter and average frequencies evaluated for the Pi states. At the UCCSD(T)/CBS(T-5) level of theory, Delta(f)H(0 K) = 138.89 kcal/mol and Delta(f)H(298 K) = 139.65 kcal/mol were obtained for diazocarbene; for the ionization potential and the electron affinity of CNN, 10.969 eV (252.95 kcal/mol), and 1.743 eV (40.19 kcal/mol), respectively, are predicted. Geometry optimization was also carried out with the CASSCF/MRCI/CBS(T-5) approach for the states X(3)Sigma(-) A(3)Pi, and a(1)Delta of CNN, and with the CASSCF/MRSDCI/aug-cc-pVTZ approach for the states b(1)Sigma(+), c(1)Pi, d(1)Sigma(-), and B(3)Sigma(-), and excitation energies (T(e)) evaluated. Vertical energies were calculated for 15 electronic states, thus improving on the accuracy of the five transitions already described, and allowing for a reliable overview of a manifold of other states, which is expected to guide future spectroscopic experiments. This study corroborates the experimental assignment for the vertical transition X (3)Sigma(-) <- E (3)Pi.

Electronic Structure and Chemical Bonding in the Lowest Electronic States of TcN

Multiconfiguration second-order perturbation theory, with the inclusion of relativistic effects and spin-orbit Coupling, was employed to investigate the nature of the ground and low-lying Lambda-S and Omega states of the TcN molecule. Spectroscopic constants, effective bond order, and potential energy curves for 13 low-lying Lambda-S states and 5 Omega states are given, The computed ground state of TcN is of Omega = 3 symmetry (R(e) = 1.605 angstrom and omega(e) = 1085 cm(-1)), originating mainly from the (3)Delta Lambda-S ground state. This result is contrasted with the nature of the ground state for other VIIB transtion-metal mononitrides, including X(3)Sigma(-) symmetry for MnN and Omega = 0(+) symmetry for ReN, derived also from a X(3)Sigma(-) state.