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.

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.

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.

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.

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.

Thiothionyl bromide (S=SBr2): The elusive isomer of dibromodisulfane (BrSSBr)

CCSD(T) with a series of correlation consistent basis up to quadruple-zeta is used to investigate the structures, vibrational spectra, relative stability, heats of formation, and barrier to isomerization of S=SBr2 and BrSSBr. It represents the most accurate and detailed characterization of these molecules to date. We show that the frequency mode at 302 cm(-1), detected in various studies and assigned to impurities by some authors, and to the anti-symmetric SBr stretch in BrSSBr by others, thus in fact corresponds to the anti-symmetric SBr stretch in the elusive S=SBr2 species; it thus corroborates and complements an earlier partial IR spectra study attributable to S=SBr2. Including corrections for relativistic and core-valence correlation effects, we also predict 26.33 (12.74) kcal/mol for Delta H-f (298.15 K) of S=SBr2 (BrSSBr). For the S=SBr2 -> BrSSBr reaction, our best estimates for the Gibbs free energy and the enthalpy of the reaction at 298.15 K are -13.71 and -13.44 kcal/mol, respectively. For a value of Delta G(#) equal to 23.52 kcal/mol, we estimate a TST rate constant, at 298.15 K, of 3.57 x 10(-5) s(-1). (c) 2007 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.

Investigating the Nature of Noble Gas-Copper Bonds by the Quantum Theory of Atoms in Molecules

We investigated noble gas copper bonds in linear complexes represented by the NgCuX general formula in which Ng and X stand for a noble gas (neon, argon, krypton, or xenon) and a halogen (fluorine, chlorine or bromine), respectively, by coupled cluster methods and modified cc-pVQZ basis sets. The quantum theory of atoms in molecules (QTAIM) shows a linear relation between the dissociation energy or noble gas-copper bonds and the amount of electronic charge transferred mainly from the noble gas to copper during complexation. Large changes in the QTAIM quadrupole moments of copper and noble gases resulting from this bonding and a comparison between NgCuX and NgNaCl systems indicate that these noble gas-copper bonds should be better interpreted as predominantly covalent. Finally, QTAIM atomic dipoles of noble gases in NgNaCl systems agree satisfactorily with atomic dipoles given by a simple model for these NgNa van der Waals bonds.

Measurement of transverse single-spin asymmetries for J/psi production in polarized p plus p collisions at root s=200 GeV

We report the first measurement of transverse single-spin asymmetries in J/psi production from transversely polarized p + p collisions at root s = 200 GeV with data taken by the PHENIX experiment in 2006 and 2008. The measurement was performed over the rapidity ranges 1.2 < vertical bar y vertical bar < 2.2 and vertical bar y vertical bar < 0.35 for transverse momenta up to 6 GeV/c. J/psi production at the Relativistic Heavy Ion Collider is dominated by processes involving initial-state gluons, and transverse single-spin asymmetries of the J/psi can provide access to gluon dynamics within the nucleon. Such asymmetries may also shed light on the long-standing question in QCD of the J/psi production mechanism. Asymmetries were obtained as a function of J/psi transverse momentum and Feynman-x, with a value of -0.086 +/- 0.026(stat) +/- 0.003(syst) in the forward region. This result suggests possible nonzero trigluon correlation functions in transversely polarized protons and, if well defined in this reaction, a nonzero gluon Sivers distribution function.

Measurement of the Parity-Violating Longitudinal Single-Spin Asymmetry for W(+/-) Boson Production in Polarized Proton-Proton Collisions at root s=500 GeV

We report the first measurement of the parity-violating single-spin asymmetries for midrapidity decay positrons and electrons from W(+) and W(-) boson production in longitudinally polarized proton-proton collisions at root s = 500 GeV by the STAR experiment at RHIC. The measured asymmetries, A(L)(W+) = -0.27 +/- 0.10(stat.) +/- 0.02(syst.) +/- 0.03(norm.) and A(L)(W-) = 0.14 +/- 0.19(stat.) +/- 0.02(syst.) +/- 0.01(norm.), are consistent with theory predictions, which are large and of opposite sign. These predictions are based on polarized quark and antiquark distribution functions constrained by polarized deep-inelastic scattering measurements.

Cross Section and Parity-Violating Spin Asymmetries of W(+/-) Boson Production in Polarized p plus p Collisions at root s=500 Gev

Large parity-violating longitudinal single-spin asymmetries A(L)(e+) = 0.86(-0.14)(+0.30) and Ae(L)(e-) = 0.88(-0.71)(+0.12) are observed for inclusive high transverse momentum electrons and positrons in polarized p + p collisions at a center-of-mass energy of root s = 500 GeV with the PHENIX detector at RHIC. These e(+/-) come mainly from the decay of W(+/-) and Z(0) bosons, and their asymmetries directly demonstrate parity violation in the couplings of the W(+/-) to the light quarks. The observed electron and positron yields were used to estimate W(+/-) boson production cross sections for the e(+/-) channels of sigma(pp -> W(+)X) X BR(W(+) -> e(+) nu(e)) = 144.1 +/- 21.2(stat)(-10.3)(+3.4)(syst) +/- 21.6(norm) pb, and sigma(pp -> W(-)X) X BR(W(-) -> e(-) (nu) over bar (e)) = 3.17 +/- 12.1(stat)(-8.2)(+10.1)(syst) +/- 4.8(norm) pb.

Gluon-Spin Contribution to the Proton Spin from the Double-Helicity Asymmetry in Inclusive pi(0) Production in Polarized p plus p Collisions at root s=200 GeV

The double helicity asymmetry in neutral pion production for p(T) = 1 to 12 GeV/c was measured with the PHENIX experiment to access the gluon-spin contribution, Delta G, to the proton spin. Measured asymmetries are consistent with zero, and at a theory scale of mu 2 = 4 GeV(2) a next to leading order QCD analysis gives Delta G([0.02,0.3]) = 0.2, with a constraint of -0.7 < Delta G([0.02,0.3]) < 0.5 at Delta chi(2) = 9 (similar to 3 sigma) for the sampled gluon momentum fraction (x) range, 0.02 to 0.3. The results are obtained using predictions for the measured asymmetries generated from four representative fits to polarized deep inelastic scattering data. We also consider the dependence of the Delta G constraint on the choice of the theoretical scale, a dominant uncertainty in these predictions.

Spin-wave fluctuations in ferrimagnetic Mg(x)Fe(3-x)O(4) nanoparticles

We have performed a systematic study of the magnetic properties of a series of ferrimagnetic nanoparticles of Mg(x)Fe(3-x)O(4) (0.8 <= x <= 1.5) prepared by the combustion reaction method. The magnetization data can be well fitted by Bloch's law with T(3/2). Bloch's constant B determined from the fitting procedure was found to increase with Mg content x from similar to 3.09 X 10(-5) K(-3/2) for x = 0.8 to 6.27 X 10(-5) K(-3/2) for x=1.5. The exchange integral J(AB) and the spin-wave stiffness constant D of Mg(x)Fe(3-x)O(4) nanoparticles were also determined as similar to 0.842 and 0.574 meV and 296 and 202 meV angstrom(2) for specimens with x=0.8 and 1.5, respectively. These results are discussed in terms of cation redistribution among A and B sites on these nanostructured spinel ferrites. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3359709]

Longitudinal double-spin asymmetry and cross section for inclusive neutral pion production at midrapidity in polarized proton collisions at s=200 GeV

We report a measurement of the longitudinal double-spin asymmetry A(LL) and the differential cross section for inclusive pi(0) production at midrapidity in polarized proton collisions at s=200 GeV. The cross section was measured over a transverse momentum range of 1 < p(T)< 17 GeV/c and found to be in good agreement with a next-to-leading order perturbative QCD calculation. The longitudinal double-spin asymmetry was measured in the range of 3.7 < p(T)< 11 GeV/c and excludes a maximal positive gluon polarization in the proton. The mean transverse momentum fraction of pi(0)'s in their parent jets was found to be around 0.7 for electromagnetically triggered events.

Longitudinal spin transfer to Lambda and Lambda hyperons in polarized proton-proton collisions at s=200 GeV

The longitudinal spin transfer, D(LL), from high energy polarized protons to Lambda and Lambda hyperons has been measured for the first time in proton-proton collisions at s=200 GeV with the STAR detector at the Relativistic Heavy Ion Collider. The measurements cover pseudorapidity, eta, in the range |eta|< 1.2 and transverse momenta, p(T), up to 4 GeV/c. The longitudinal spin transfer is found to be D(LL)=-0.03 +/- 0.13(stat)+/- 0.04(syst) for inclusive Lambda and D(LL)=-0.12 +/- 0.08(stat)+/- 0.03(syst) for inclusive Lambda hyperons with <>=0.5 and << p(T)>>=3.7 GeV/c. The dependence on eta and p(T) is presented.