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.

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.

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 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.

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.

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.

Axial ligand influence on geometries, charge distributions and electronic structures of iron tetraazamacrocycle [Fe((II))TIM(X)(Y)](2+) complexes assessed by Density Functional Theory

We employed the Density Functional Theory along with small basis sets, B3LYP/LANL2DZ, for the study of FeTIM complexes with different pairs of axial ligands (CO, H(2)O, NH(3), imidazole and CH(3)CN). These calculations did not result in relevant changes of molecular quantities as bond lengths, vibrational frequencies and electronic populations supporting any significant back-donation to the carbonyl or acetonitrile axial ligands. Moreover, a back-donation mechanism to the macrocycle cannot be used to explain the observed changes in molecular properties along these complexes with CO or CH(3)CN. This work also indicates that complexes with CO show smaller binding energies and are less stable than complexes with CH(3)CN. Further, the electronic band with the largest intensity in the visible region (or close to this region) is associated to the transition from an occupied 3d orbital on iron to an empty pi* orbital located at the macrocycle. The energy of this Metal-to-Ligand Charge Transfer (MLCT) transition shows a linear relation to the total charge of the macrocycle in these complexes as given by Mulliken or Natural Population Analysis (NPA) formalisms. Finally, the macrocycle total charge seems to be influenced by the field induced by the axial ligands. (C) 2011 Elsevier Ltd. All rights reserved.

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.

A quantum chemical study on the formation of phosphorus mononitride

The chemical mechanism of the (1)PN formation was successfully studied by using the CCSD(T)/6-311++G(3df,3pd) level of theory. The (1)NH(3) + (3)PH and (4)P + NH(3) reaction paths are not energetically favorable to form the (1)PN molecule. However, the (3)NH + (3)PH, (4)N + (3)PH(3), (4)N + (3)PH, (4)P + (3)NH, and (4)P + (2)NH(2) reaction paths to form the (1)PN molecule are only energetically favorable by taking place through specific transition states to form the (1)PN molecule. The NH(3) + (3)PH, (4)N + (1)PH(3), NH(3) + (4)P, and (4)N + (2)PH(2) reactions are spin-forbidden and the probability of hopping for these reactions was estimated to be 0 by the Landau-Zener theory. This is the first detailed study on the chemical mechanism for the (1)PN formation. (C) 2009 Elsevier B.V. All rights reserved.

Propriedades ópticas de nanocristais de SiGe: efeitos de dopagem e desordem

We have used ab initio calculations to investigate the electronic structure of SiGe based nanocrystals (NC s). This work is divided in three parts. In the first one, we focus the excitonic properties of Si(core)/Ge(shell) and Ge(core)/Si(shell) nanocrystals. We also estimate the changes induced by the effect of strain the electronic structure. We show that Ge/Si (Si/Ge) NC s exhibits type II confinement in the conduction (valence) band. The estimated potential barriers for electrons and holes are 0.16 eV (0.34 eV) and 0.64 eV (0.62 eV) for Si/Ge (Ge/Si) NC s. In contradiction to the expected long recombination lifetimes in type II systems, we found that the recombination lifetime of Ge/Si NC s (τR = 13.39μs) is more than one order of magnitude faster than in Si/Ge NC s (τR = 191.84μs). In the second part, we investigate alloyed Si1−xGex NC s in which Ge atoms are randomly positioned. We show that the optical gaps and electron-hole binding energies decrease linearly with x, while the exciton exchange energy increases with x due to the increase of the spatial extent of the electron and hole wave functions. This also increases the electron-hole wave functions overlap, leading to recombination lifetimes that are very sensitive to the Ge content. Finally, we investigate the radiative transitions in Pand B-doped Si nanocrystals. Our NC sizes range between 1.4 and 1.8 nm of diameters. Using a three-levels model, we show that the radiative lifetimes and oscillator strengths of the transitions between the conduction and the impurity bands, as well as the transitions between the impurity and the valence bands are strongly affected by the impurity position. On the other hand, the direct conduction-to-valence band decay is practically unchanged due to the presence of the impurity

Some exact solutions of the Dirac equation

Exact analytic solutions are found to the Dirac equation for a combination of Lorentz scalar and vector Coulombic potentials with additional non-Coulombic parts. An appropriate linear combination of Lorentz scalar and vector non-Coulombic potentials, with the scalar part dominating, can be chosen to give exact analytic Dirac wave functions.