Combining Monte Carlo simulation and density-functional theory to describe the spectral changes of Na(2) in liquid helium

Spectral changes of Na(2) in liquid helium were studied using the sequential Monte Carlo-quantum mechanics method. Configurations composed by Na(2) surrounded by explicit helium atoms sampled from the Monte Carlo simulation were submitted to time-dependent density-functional theory calculations of the electronic absorption spectrum using different functionals. Attention is given to both line shift and line broadening. The Perdew, Burke, and Ernzerhof (PBE1PBE, also known as PBE0) functional, with the PBE1PBE/6-311++G(2d,2p) basis set, gives the spectral shift, compared to gas phase, of 500 cm(-1) for the allowed X (1)Sigma(+)(g) -> B (1)Pi(u) transition, in very good agreement with the experimental value (700 cm(-1)). For comparison, cluster calculations were also performed and the first X (1)Sigma(+)(g) -> A (1)Sigma(+)(u) transition was also considered.

Electronic structure of In(2)O(3) and Sn-doped In(2)O(3) by hard x-ray photoemission spectroscopy

The valence and core levels of In(2)O(3) and Sn-doped In(2)O(3) have been studied by hard x-ray photoemission spectroscopy (hv = 6000 eV) and by conventional Al K alpha (hv = 1486.6 eV) x-ray photoemission spectroscopy. The experimental spectra are compared with density-functional theory calculations. It is shown that structure deriving from electronic levels with significant In or Sn 5s character is selectively enhanced under 6000 eV excitation. This allows us to infer that conduction band states in Sn-doped samples and states at the bottom of the valence band both contain a pronounced In 5s contribution. The In 3d core line measured at hv = 1486.6 eV for both undoped and Sn-doped In(2)O(3) display an asymmetric lineshape, and may be fitted with two components associated with screened and unscreened final states. The In 3d core line spectra excited at hv = 6000 eV for the Sn-doped samples display pronounced shoulders and demand a fit with two components. The In 3d core line spectrum for the undoped sample can also be fitted with two components, although the relative intensity of the component associated with the screened final state is low, compared to excitation at 1486.6 eV. These results are consistent with a high concentration of carriers confined close to the surface of nominally undoped In(2)O(3). This conclusion is in accord with the fact that a conduction band feature observed for undoped In(2)O(3) in Al K alpha x-ray photoemission is much weaker than expected in hard x-ray photoemission.

The role of electron localization in the atomic structure of transition-metal 13-atom clusters: the example of Co(13), Rh(13), and Hf(13)

The crystalline structure of transition-metals (TM) has been widely known for several decades, however, our knowledge on the atomic structure of TM clusters is still far from satisfactory, which compromises an atomistic understanding of the reactivity of TM clusters. For example, almost all density functional theory (DFT) calculations for TM clusters have been based on local (local density approximation-LDA) and semilocal (generalized gradient approximation-GGA) exchange-correlation functionals, however, it is well known that plain DFT fails to correct the self-interaction error, which affects the properties of several systems. To improve our basic understanding of the atomic and electronic properties of TM clusters, we report a DFT study within two nonlocal functionals, namely, the hybrid HSE (Heyd, Scuseria, and Ernzerhof) and GGA + U functionals, of the structural and electronic properties of the Co(13), Rh(13), and Hf(13) clusters. For Co(13) and Rh(13), we found that improved exchange-correlation functionals decrease the stability of open structures such as the hexagonal bilayer (HBL) and double simple-cubic (DSC) compared with the compact icosahedron (ICO) structure, however, DFT-GGA, DFT-GGA + U, and DFT-HSE yield very similar results for Hf(13). Thus, our results suggest that the DSC structure obtained by several plain DFT calculations for Rh(13) can be improved by the use of improved functionals. Using the sd hybridization analysis, we found that a strong hybridization favors compact structures, and hence, a correct description of the sd hybridization is crucial for the relative energy stability. For example, the sd hybridization decreases for HBL and DSC and increases for ICO in the case of Co(13) and Rh(13), while for Hf(13), the sd hybridization decreases for all configurations, and hence, it does not affect the relative stability among open and compact configurations.

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]

Full-dimensional analytical ab initio potential energy surface of the ground state of HOI

Extensive ab initio calculations using a complete active space second-order perturbation theory wavefunction, including scalar and spin-orbit relativistic effects with a quadruple-zeta quality basis set were used to construct an analytical potential energy surface (PES) of the ground state of the [H, O, I] system. A total of 5344 points were fit to a three-dimensional function of the internuclear distances, with a global root-mean-square error of 1.26 kcal mol(-1). The resulting PES describes accurately the main features of this system: the HOI and HIO isomers, the transition state between them, and all dissociation asymptotes. After a small adjustment, using a scaling factor on the internal coordinates of HOI, the frequencies calculated in this work agree with the experimental data available within 10 cm(-1). (C) 2011 American Institute of Physics. [doi: 10.1063/1.3615545]

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

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.

Characterization and spectroscopic analysis of phenol-ethanol hydrogen bonded clusters

Different hydrogen bonded clusters involving phenol and ethanol are studied theoretically using MP2/aug-cc-pVDZ. Nine different 1: 1 clusters are obtained and analyzed according to their stability and spectroscopic properties. Different isomeric forms of ethanol are considered. Attention is also devoted to the spectral shift of the characteristic pi -> pi* transition of phenol. Using TDHF, CIS, CIS(D) and TDB3LYP in aug-cc-pVDZ basis set, all results agree that a red shift is obtained when phenol is the hydrogen donor and a blue shift is obtained in the opposite case. These results are used to rationalize the red shift observed for phenol in liquid ethanol. (C) 2010 Elsevier B.V. All rights reserved.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Self-assembling of phenothiazine compounds investigated by small-angle X-ray scattering and electron paramagnetic resonance spectroscopy

Small-angle X-ray scattering (SAXS) and electron paramagnetic resonance (EPR) have been carried out to investigate the structure of the self-aggregates of two phenothiazine drugs, chlorpromazine (CPZ) and trifluoperazine (TFP), in aqueous solution. In the SAXS studies, drug solutions of 20 and 60 mM, at pH 4.0 and 7.0, were investigated and the best data fittings were achieved assuming several different particle form factors with a homogeneous electron density distribution in respect to the water environment. Because of the limitation of scattering intensity in the q range above 0.15 angstrom(-1), precise determination of the aggregate shape was not possible and all of the tested models for ellipsoids, cylinders, or parallelepipeds fitted the experimental data equally well. The SAXS data allows inferring, however, that CPZ molecules might self-assemble in a basis set of an orthorhombic cell, remaining as nanocrystallites in solution. Such nanocrystals are composed of a small number of unit cells (up to 10, in c-direction), with CPZ aggregation numbers of 60-80. EPR spectra of 5- and 16-doxyl stearic acids bound to the aggregates were analyzed through simulation, and the dynamic and magnetic parameters were obtained. The phenothiazine concentration in EPR experiments was in the range of 5-60 mM. Critical aggregation concentration of TFP is lower than that for CPZ, consistent with a higher hydrophobicity of TFP. At acidic pH 4.0 a significant residual motion of the nitroxide relative to the aggregate is observed, and the EPR spectra and corresponding parameters are similar to those reported for aqueous surfactant micelles. However, at pH 6.5 a significant motional restriction is observed, and the nitroxide rotational correlation times correlate very well with those estimated for the whole aggregated particle from SAXS data. This implies that the aggregate is densely packed at this pH and that the nitroxide is tightly bound to it producing a strongly immobilized EPR spectrum. Besides that, at pH 6.5 the differences in motional restriction observed between 5- and 16-DSA are small, which is different from that observed for aqueous surfactant micelles.

Size effects on silver nanoparticles` properties

In this work a systematic study of the dependence of the structural, electronic, and vibrational properties on nanoparticle size is performed. Based on our total energy calculations we identified three characteristic regimes associated with the nanoparticle`s dimensions: (i) below 1.5 nm (100 atoms) where remarkable molecular aspects are observed; (ii) between 1.5 and 2.0 nm (100 and 300 atoms) where the molecular behavior is influenced by the inner core crystal properties; and (iii) above 2.0 nm (more than 300 atoms) where the crystal properties are preponderant. In all considered regimes the nanoparticle`s surface modulates its properties. This modulation decreases with the increasing of the nanoparticle`s size.

DFT study of the electronic, vibrational, and optical properties of SnO(2)

We report results on the electronic, vibrational, and optical properties of SnO(2) obtained using first-principles calculations performed within the density functional theory. All the calculated phonon frequencies, real and imaginary parts of complex dielectric function, the energy-loss spectrum, the refractive index, the extinction, and the absorption coefficients show good agreement with experimental results. Based on our calculations, the SnO(2) electron and hole effective masses were found to be strongly anisotropic. The lattice contribution to the low-frequency region of the SnO(2) dielectric function arising from optical phonons was also determined resulting the values of E > (1aSyen) (latt) (0) = 14.6 and E > (1ayen) (latt) (0) = 10.7 for directions perpendicular and parallel to the tetragonal c-axis, respectively. This is in excellent agreement with the available experimental data. After adding the electronic contribution to the lattice contribution, a total average value of E >(1)(0) = 18.2 is predicted for the static permittivity constant of SnO(2).

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.