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.

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.

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

Systematic investigation of a family of gradient-dependent functionals for solids

Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. On a test set of 60 solids, we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.

Low-energy electron scattering from methanol and ethanol

Measured and calculated differential cross sections for elastic (rotationally unresolved) electron scattering from two primary alcohols, methanol (CH(3)OH) and ethanol (C(2)H(5)OH), are reported. The measurements are obtained using the relative flow method with helium as the standard gas and a thin aperture as the collimating target gas source. The relative flow method is applied without the restriction imposed by the relative flow pressure conditions on helium and the unknown gas. The experimental data were taken at incident electron energies of 1, 2, 5, 10, 15, 20, 30, 50, and 100 eV and for scattering angles of 5 degrees-130 degrees. There are no previous reports of experimental electron scattering differential cross sections for CH(3)OH and C(2)H(5)OH in the literature. The calculated differential cross sections are obtained using two different implementations of the Schwinger multichannel method, one that takes all electrons into account and is adapted for parallel computers, and another that uses pseudopotentials and considers only the valence electrons. Comparison between theory and experiment shows that theory is able to describe low-energy electron scattering from these polyatomic targets quite well.

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]

DFT and electrochemical studies on nortriptyline oxidation sites

A study on the possible sites of oxidation and epoxidation of nortriptyline was performed using electrochemical and quantum chemical methods; these sites are involved in the biological responses (for example, hepatotoxicity) of nortriptyline and other similar antidepressants. Quantum chemical studies and electrochemical experiments demonstrated that the oxidation and epoxidation sites are located on the apolar region of nortriptyline, which will useful for understanding the molecule`s activity. Also, for the determination of the compound in biological fluids or in pharmaceutical formulations, we propose a useful analytical methodology using a graphite-polyurethane composite electrode, which exhibited the best performance when compared with boron-doped diamond or glassy carbon surfaces.

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters,. but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.

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.

Electronic properties of a methane-water solution

Electronic properties of a methane-water solution were investigated by a sequential quantum mechanical/molecular dynamics approach. Upon hydration methane acquires an induced dipole moment of similar to 0.5 +/- 0.2 D. This is related to polarisation effects and to weak methane-water hydrogen bond interactions. From gas phase to solution, the first vertical excitation and ionisation energies of methane are red-shifted by 0.45 +/- 0.25 and 0.87 +/- 0.40 eV, respectively. We also report results for the dynamic polarisability of methane in water. In comparison with water, no difference was found for the average monomeric dipole moment of water molecules in close interaction with methane. (c) 2011 Elsevier B.V. All rights reserved.

Large Changes of Static Electric Properties Induced by Hydrogen Bonding: An ab Initio Study of Linear HCN Oligomers

We report the partitioning of the interaction-induced static electronic dipole (hyper)polarizabilities for linear hydrogen cyanide complexes into contributions arising from various interaction energy terms. We analyzed the nonadditivities of the studied properties and used these data to predict the electric properties of an infinite chain. The interaction-induced static electric dipole properties and their nonadditivities were analyzed using an approach based on numerical differentiation of the interaction energy components estimated in an external electric field. These were obtained using the hybrid variational-perturbational interaction energy decomposition scheme, augmented with coupled-cluster calculations, with singles, doubles, and noniterative triples. Our results indicate that the interaction-induced dipole moments and polarizabilities are primarily electrostatic in nature; however, the composition of the interaction hyperpolarizabilities is much more complex. The overlap effects substantially quench the contributions due to electrostatic interactions, and therefore, the major components are due to the induction and exchange induction terms, as well as the intramolecular electron-correlation corrections. A particularly intriguing observation is that the interaction first hyperpolarizability in the studied systems not only is much larger than the corresponding sum of monomer properties, but also has the opposite sign. We show that this effect can be viewed as a direct consequence of hydrogen-bonding interactions that lead to a decrease of the hyperpolarizability of the proton acceptor and an increase of the hyperpolarizability of the proton donor. In the case of the first hyperpolarizability, we also observed the largest nonadditivity of interaction properties (nearly 17%) which further enhances the effects of pairwise interactions.

Study of the spectroscopic properties and first hyperpolarizabilities of disperse azo dyes derived from 2-amino-5-nitrothiazole

The solvatochromism and other spectroscopic and photophysical characteristics of four azo disperse dyes, derived from 2-amino-5-nitrothiazole, were evaluated and interpreted with the aid of experimental data and quantum mechanical calculations. For the non-substituted compound two conformers, E and Z, were proposed for the isolated molecules, being the second one considerably less stable. The optimization of these structures in combination with a SCRF methodology (IEFPCM, Simulating the molecules in a continuum dielectric with characteristics of methanol), suggests that the Z form is not stable in solution. This same behaviour is expected for the substituted compounds, which is corroborated by experimental data presented in previous investigations [A.E.H. Machado, L.M. Rodrigues, S. Gupta, A.M.F. Oliveira-Campos, A.M.S. Silva, J. Mol. Struct. 738 (2005) 239-245]. For the substituted compounds, two forms derived from E conformer (A and R) are possible. Quantum mechanical data suggest for the isolated molecules, that the low energy absorption hand of the E conformers involve at least two close electronic states. having the low-lying excited state a (1)(n,pi*) nature, and being the S-2 state attributed to a (1)(pi,pi*) transition. The data also suggest a small energy gap between the absorption peaks of A and B, related to the easy conversion between these forms. For the structures optimized in combination with the applied SCRF methodology, an states inversion is observed for the Substituted compounds, with a considerable diminish of the energy gap between A and B absorption peaks. The electronic spectra of these compounds are quite sensitive to changes in the solvent polarity. The positive solvatochromism is more evident in aprotic solvents, probably due to the polarization induced by the solute. These compounds do not fluoresce at 298 K, but present a small but perceptible fluorescence at 77 K, which seems to be favoured by the nature of the group in the 2 `-position of the phenyl ring. Moreover, such compounds present expressive values for first hyperpolarizability, which implies in good non-linear optics (NLO) responses and photoswitching capability. (C) 2008 Elsevier B.V. All rights reserved.

Spectroscopic and structural studies of bis[isopropoxy(thiocarbonyl)]sulfide, [(CH(3))(2)CHOC(S)](2)S

Structural and conformational properties of the molecule bis[isopropoxy(thiocarbonyl)]sulfide, [(CH(3))(2)CHOC(S)](2)S, have been studied by vibrational spectroscopy (IR and Raman) and quantum chemical calculations (HF and B3LYP with 6-31+G* basis sets). The crystal and molecular structure of the title compound was determined by X-ray diffraction methods. It crystallizes in the monoclinic C2/c space group with a = 8.4007(4), b = 13.5936(5), c = 10.3648(5) angstrom, beta = 106.024(4)degrees and Z = 4 molecules per unit cell. The molecules are sited on a crystallographic twofold axis passing through the sulphide atom and arranged in layers perpendicular to the b-axis. The solid state IR and Raman spectra of the compound give no sign of any other rotamer. (C) 2009 Elsevier B.V. All rights reserved.