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.

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.

Isomers on the [H, S(2), Cl] potential energy surface: A high level investigation

This study reports a systematic state-of-the-art characterization of new sulfur-chlorine species on the [H, S(2), Cl] potential energy surface. Coupled cluster theory singles and doubles with perturbative contributions of connected triples, using the series of correlation consistent basis sets with extrapolations to the complete basis set limit (CBS), were employed to quantify the energetic quantities involved in the isomerization processes on this surface. The structures and vibrational frequencies are unique for some species and represent the most accurate investigation to date. These molecules are potentially a new route of coupling the sulfur and chlorine chemistries in the atmosphere, and conditions of high concentration of H(2)S (HS) like in volcanic eruptions might contribute to their formation. Also an assessment of the MP2/CBS approach relative to CCSD(T)/CBS provides insights on the expected performance of MP2/CBS on the characterization of polysulfides, and also of more complex systems containing disulfide bridges. (C) 2009 Elsevier B.V. All rights reserved.

The reaction of methyl radical with nitrogen atom on the triplet potential energy surface: A CCSD(T)/CBS characterization

CCSD(T)/cc-pVnZ (n = D, T, Q) calculations followed by extrapolations to the CBS limit are used to characterize stationary states of species participating in the N((4)S) + CH(3) (2A ``) reaction on the triplet PES. A mechanistic model is investigated and reaction rates are computed for every step and the overall reaction. Our best CBS estimate (1.93 x 10(10) cm(3) molecule(1) s(1)) for the overall rate constant leading to the formation of H(2)CN + H compares well with the experimental values (8.5 x 10 (11) and 1.3 x 10(10) cm(3) molecule(1) s(1)), thus reducing significantly the discrepancy of a previous theoretical result (9.1 x 10(12) cm(3) molecule(1) s(1)). (C) 2008 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.

The energy cycle and structural evolution of cyclones over southeastern South America in three case studies

In this paper, the Lorenz energy cycle over a limited area was applied for three cyclones with different origins and evolutions, where each of them was formed in an important cyclogenetic region near southeastern South America. The synoptic conditions and energetics were analyzed during each system`s life cycle and showed important relationships between their energy cycle and the evolution of their vertical structure. In the case of the weak baroclinic cyclone which formed on Brazil`s south-southeastern coast, the analysis showed that it originated through a midlevel cutoff low with contribution from barotropic instability. Its evolution would indicate potential transition to a hybrid system if the convective activity were stronger. The system that occurred in the La Plata River mouth had features of an oceanic bomb-type cyclogenesis and showed an important contribution from the available potential energy generation term through the latent heat release by the convection. Meanwhile, the system of the southern Argentina coast presented a classical baroclinic development of extratropical cyclogenesis in the energy cycle, from the wave amplification up to the final occlusion of the associated frontal system. These analyses revealed that the development of some cyclones that occur in eastern South America can present different mechanisms that are not related to the classical extratropical cyclogenesis.

Conservative force fields in non-Gaussian statistics

In this Letter, we determine the kappa-distribution function for a gas in the presence of an external field of force described by a potential U(r). In the case of a dilute gas, we show that the kappa-power law distribution including the potential energy factor term can rigorously be deduced in the framework of kinetic theory with basis on the Vlasov equation. Such a result is significant as a preliminary to the discussion on the role of long range interactions in the Kaniadakis thermostatistics and the underlying kinetic theory. (C) 2008 Elsevier B.V. All rights reserved.

Classical simulation of deposition of thiophene oligomers on TiO(2)-anatase: Relevance of long-range electrostatic interactions

We performed classical molecular dynamics simulations of the vapor-deposition of alpha-T4 oligomers on the TiO(2)-anatase (101) surface, comparing different sets of charges associated with the atoms of the model. The potential energy surfaces for alpha-T4 and TiO(2) were described by re-parametrizations of the Universal force field with charges given by the charge equilibration (QEq) scheme, or with fixed charges obtained by an ab initio method using the Hirshfeld partition. The two sets of charges lead to completely different results for the interface formation, and for the characteristics of the organic film, with a clearly defined alpha-T4 contact layer in the QEq case, and a more homogeneous molecular distribution when using Hirshfeld charges. The main reason for the discrepancy was found to be the incorrect charge assignment given by QEq to the sulfur and alpha-carbon atoms in thiophenes, and highlight the relevance of long-range interactions in the organization of molecular films. (C) 2009 Elsevier B.V. All rights reserved.

Instability of equilibrium points of some Lagrangian systems

In this work we show that, if L is a natural Lagrangian system such that the k-jet of the potential energy ensures it does not have a minimum at the equilibrium and such that its Hessian has rank at least n - 2, then there is an asymptotic trajectory to the associated equilibrium point and so the equilibrium is unstable. This applies, in particular, to analytic potentials with a saddle point and a Hessian with at most 2 null eigenvalues. The result is proven for Lagrangians in a specific form, and we show that the class of Lagrangians we are interested can be taken into this specific form by a subtle change of spatial coordinates. We also consider the extension of this results to systems subjected to gyroscopic forces. (C) 2008 Elsevier Inc. All rights reserved.

Multiple Brake Orbits and Homoclinics in Riemannian Manifolds

Let (M, g) be a complete Riemannian manifold, Omega subset of Man open subset whose closure is homeomorphic to an annulus. We prove that if a,Omega is smooth and it satisfies a strong concavity assumption, then there are at least two distinct geodesics in starting orthogonally to one connected component of a,Omega and arriving orthogonally onto the other one. Using the results given in Giamb et al. (Adv Differ Equ 10:931-960, 2005), we then obtain a proof of the existence of two distinct homoclinic orbits for an autonomous Lagrangian system emanating from a nondegenerate maximum point of the potential energy, and a proof of the existence of two distinct brake orbits for a class of Hamiltonian systems. Under a further symmetry assumption, the result is improved by showing the existence of at least dim(M) pairs of geometrically distinct geodesics as above, brake orbits and homoclinic orbits. In our proof we shall use recent deformation results proved in Giamb et al. (Nonlinear Anal Ser A: Theory Methods Appl 73:290-337, 2010).

On the Multiplicity of Orthogonal Geodesics in Riemannian Manifold With Concave Boundary. Applications to Brake Orbits and Homoclinics

Let (M, g) be a complete Riemannian Manifold, Omega subset of M an open subset whose closure is diffeomorphic to an annulus. If partial derivative Omega is smooth and it satisfies a strong concavity assumption, then it is possible to prove that there are at least two geometrically distinct geodesics in (Omega) over bar = Omega boolean OR partial derivative Omega starting orthogonally to one connected component of partial derivative Omega and arriving orthogonally onto the other one. The results given in [6] allow to obtain a proof of the existence of two distinct homoclinic orbits for an autonomous Lagrangian system emanating from a nondegenerate maximum point of the potential energy, and a proof of the existence of two distinct brake orbits for a. class of Hamiltonian systems. Under a further symmetry assumption, it is possible to show the existence of at least dim(M) pairs of geometrically distinct geodesics as above, brake orbits and homoclinics.

Electronic Structure of the Ground and Low-Lying Electronic States of MoB and MoB(+)

Multiconfigurational second-order perturbation theory (CASSCF//CASPT2) and quadruple-zeta ANO-RCC basis sets were employed to investigate the ground and low-lying electronic states of MoB and MoB(+). Spectroscopic constants, potential energy curves, wavefunctions, Mulliken population analyses, and ionization energies are given. The ground state of MoB is of X(6)Pi symmetry (R(e) = 1.968 angstrom, omega(e) = 664 cm(-1), and mu = 2.7 D), giving rise to a Omega = 7/2 ground state after including spin-orbit coupling. For MoB(+), the ground state is computed to be of X(7)Sigma(+) symmetry (R(e) = 2.224 angstrom, omega(e) = 141 cm(-1), and mu = 1.2 D), with an adiabatic ionization energy of 7.19 eV and a vertical one of 7.53 eV. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem 111: 3362-3370, 2011

The low-lying electronic states of BeAs: a first principles characterization

The electronic structure and spectroscopic properties of a manifold of states of a new molecular species, BeAs, have been investigated theoretically at the complete active space self-consistent field/multireference single and double excitations configuration interaction (CASSCF/MRSDCI) approach, using the aug-cc-pV5Z-PP basis set for arsenic, which includes a relativistic effective core potential, and the cc-pV5Z set for beryllium. Potential energy curves of five quartet and eight doublet (I > + S) states correlating with the five lowest-lying dissociation limit are constructed. The effect of spin-orbit coupling is also included in the description of the ground state, and of the doublet states correlating with the second dissociation channel. Dipole moment functions and vibrationally averaged dipole moments are also evaluated. The similarities and differences between BeAs, BeP, and BeN are analyzed. Spin-orbit effects are small for the ground state close to the equilibrium distance, but avoided crossings between Omega = 1/2 states, and between Omega = 3/2 states changes significantly the I > + S curves for the lowest-lying doublets.

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.

Excited Electronic States, Transition Probabilities, and Radiative Lifetimes of CAs: A Theoretical Contribution and Challenge to Experimentalists

High-level CASSCF/MRCI calculations with a quintuple-zeta quality basis set are reported by characterizing for the first time a manifold of electronic states of the CAs radical yet to be investigated experimentally. Along with the potential energy curves and the associated spectroscopic constants, the dipole moment functions for selected electronic states as well as the transition dipole moment functions for the most relevant electronic transitions are also presented. Estimates of radiative transition probabilities and lifetimes complement this investigation, which also assesses the effect of spin-orbit interaction on the A (2)Pi state. Whenever pertinent, comparisons of similarities and differences with the isovalent CN and CP radicals are made.