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.

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.

Vibrational properties, crystal X-ray diffraction structure and quantum chemical calculations on a divalent sulfur substituted phthalimide: 1H-Isoindole-1,3(2H)-dione, 2-[(methoxycarbonyl)thio]

Structural and conformational properties of 1H-Isoindole-1,3(2H)-dione, 2-[(methoxycarbonyl)thio] (S-phthalimido O-methyl thiocarbonate) are analyzed using a combined approach including X-ray diffraction, vibrational spectra and theoretical calculation methods. The vibrational properties have been studied by infrared and Raman spectroscopies along with quantum chemical calculations (B3LYP and B3PW91 functional in connection with the 6-311++G** and aug-cc-pVDZ basis sets). The crystal structure was determined by X-ray diffraction methods. The substance crystallizes in the monoclinic P2(1)/c space group with a = 6.795(1), b = 5.109(1), c = 30.011(3) angstrom, beta = 90.310(3)degrees and Z = 4 molecules per unit cell. The conformation adopted by the N-S-C=O group is syn (C=O double bond in synperiplanar orientation with respect to the N-S single bond). The experimental molecular structure is well reproduced by the MP2/aug-cc-pVDZ method. (C) 2010 Elsevier B.V. All rights reserved.

Cyclocreatine, an anticancer and neuroprotective agent. Spectroscopic, structural and theoretical study

The structural, spectroscopic and theoretical study of cyclocreatine (1-carboxymethyl-2-iminoimidazolidine, CyCre) has been performed prompted by the biological relevance of the molecule and its potential role as a ligand in biometalic compounds. The crystal structure of CyCre has been determined by X-ray diffraction methods. The compound crystallizes as a zwitterion in the monoclinic system, space group P2(1)/c. The crystal is further stabilized by a network of N-H center dot center dot center dot O bonds. Infrared and Raman spectra of the solid, electronic spectra of aqueous solutions at different pH values and (1)H and (13)C NMR spectra have been recorded and analyzed. Band assignments were accomplished with the help of theoretical calculations. Optimized molecular geometries, harmonic vibrational frequencies and molecular electrostatic potentials were calculated using methods based on the density functional theory. (C) 2010 Elsevier B.V. All rights reserved.

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

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.

Solvatochromic Shift of the pi -> pi* Transition in All-Trans, Cis-13, Cis-11, Cis-9, and Cis-7 Retinal Isomers Induced by Water and Methanol

The solvatochromic shift of the lowest singlet it pi -> pi* electronic transition in the all-trans, cis-13, cis-11, cis-9, and cis-7 retinal isomers were computed under the influence of water, methanol, and benzene solvents. Excitation energies were calculated in gas phase and in solution. The calculations in solution were performed considering the sequential Monte Carlo (MC) /Quantum Mechanical approach. The MC simulations were performed considering the full retinal isomer molecules and 900 water molecules, 900 methanol, or 400 benzene ones. The OPLS/AA parametrization was chosen for retinal, methanol, and benzene molecules and the SPC model was used for water one. From the MC calculations 100 independent configurations were selected, with 100 solvent molecules in thermodynamical equilibrium at T = 298.15 K. Average point-charges were obtained from those independent configurations for water, methanol, and benzene solvent. TDDFT and CASSCF//CASPT2 methodologies were used to compute the vertical excitation energy of the retinal isomers in different environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2076-2087, 2010

Electronic structure and chemical bonding in W(2) molecule

The electronic structure of the lowest-lying electronic states of W(2) were investigated at the CASPT2 level. The ground state is a X(1)Sigma(+)(g) state, followed by the a(3)Delta(u), b(3)Sigma(+)(u) and A(1)Delta(u) electronic states. Seven low-lying Omega-states were computed: (1)0(g)(+), (2)3(u), (3)2(u), (4)1(u), (5)0(u)(-), (6)1(u), and (7)2(u), with the ground state corresponding to the (1)0(g)(+)(X(1)Sigma(+)(g)) state. Comparison with the other VIB transition metal group dimers indicates a common pattern of electronic structure and spectroscopic properties. (C) 2010 Elsevier B.V. All rights reserved.

Solvatochromic and lonochromic Effects of Iron(II)bis-(1,10-phenanthroline)dicyano: a Theoretical Study

Solvatochromic and ionochromic effects of the iron(II)bis(1,10-phenanthroline)dicyano (Fe(phen)(2)(CN)(2)) complex were investigated by means of combined DFT/TDDFT calculations using the PBE and B3LYP functionals. Extended solvation models of Fe(phen)(2)(CN)(2) in acetonitrile and aqueous solution, as well as including interaction with Mg(2+), were constructed. The calculated vertical excitation energies reproduce well the observed solvatochromism in acetonitrile and aqueous solutions, the ionochromism in acetonitrile in the presence of Mg(2+), and the absence of ionochromic effect in aqueous solution. The vertical excitation energies and the nature of the transitions were reliably predicted after inclusion of geometry relaxation upon aqueous micro- and global solvation and solvent polarization effect in the TDDFT calculations. The two intense UV-vis absorption bands occurring for all systems studied are interpreted as transitions from a hybrid Fe(II)(d)/cyano N(p) orbital to a phenanthroline pi* orbital rather than a pure metal-to-ligand-charge transfer (MLCT). The solvatochromic and ionochromic blue band shifts of Fe(phen)(2)(CN)(2) were explained with preferential stabilization of the highest occupied Fe(II)(d)/cyano N(p) orbitals as a result of specific interactions with water solvent molecules or Mg(2+) ions in solution. Such interactions occur through the CN(-) groups in the complex, and they have a decisive role for the observed blue shifts of UV-vis absorption bands.

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.