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.

Density functional theory investigation of 3d, 4d, and 5d 13-atom metal clusters

The knowledge of the atomic structure of clusters composed by few atoms is a basic prerequisite to obtain insights into the mechanisms that determine their chemical and physical properties as a function of diameter, shape, surface termination, as well as to understand the mechanism of bulk formation. Due to the wide use of metal systems in our modern life, the accurate determination of the properties of 3d, 4d, and 5d metal clusters poses a huge problem for nanoscience. In this work, we report a density functional theory study of the atomic structure, binding energies, effective coordination numbers, average bond lengths, and magnetic properties of the 3d, 4d, and 5d metal (30 elements) clusters containing 13 atoms, M(13). First, a set of lowest-energy local minimum structures (as supported by vibrational analysis) were obtained by combining high-temperature first- principles molecular-dynamics simulation, structure crossover, and the selection of five well-known M(13) structures. Several new lower energy configurations were identified, e. g., Pd(13), W(13), Pt(13), etc., and previous known structures were confirmed by our calculations. Furthermore, the following trends were identified: (i) compact icosahedral-like forms at the beginning of each metal series, more opened structures such as hexagonal bilayerlike and double simple-cubic layers at the middle of each metal series, and structures with an increasing effective coordination number occur for large d states occupation. (ii) For Au(13), we found that spin-orbit coupling favors the three-dimensional (3D) structures, i.e., a 3D structure is about 0.10 eV lower in energy than the lowest energy known two-dimensional configuration. (iii) The magnetic exchange interactions play an important role for particular systems such as Fe, Cr, and Mn. (iv) The analysis of the binding energy and average bond lengths show a paraboliclike shape as a function of the occupation of the d states and hence, most of the properties can be explained by the chemistry picture of occupation of the bonding and antibonding states.

Density functional theory and ab initio molecular dynamics study of NO adsorption on Pd(111) and Pt(111) surfaces

The origin of the unique geometry for nitric oxide (NO) adsorption on Pd(111) and Pt(111) surfaces as well as the effect of temperature were studied by density functional theory calculations and ab initio molecular dynamics at finite temperature. We found that at low coverage, the adsorption geometry is determined by electronic interactions, depending sensitively on the adsorption sites and coverages, and the effect of temperature on geometries is significant. At coverage of 0.25 monolayer (ML), adsorbed NO at hollow sites prefer an upright configuration, while NO adsorbed at top sites prefer a tilting configuration. With increase in the coverage up to 0.50 ML, the enhanced steric repulsion lead to the tilting of hollow NO. We found that the tilting was enhanced by the thermal effects. At coverage of 0.75 ML with p(2 x 2)-3NO(fcc+hcp+top) structure, we found that there was no preferential orientation for tilted top NO. The interplay of the orbital hybridization, thermal effects, steric repulsion, and their effects on the adsorption geometries were highlighted at the end.

Semiconductor nanoparticle modeling via density functional theory

In this work, a 2.0 nm nanoparticle (low limit synthesized system) is compared to possible simplified models: passivated clusters, small (1.3 nm) nanoparticles and sets of plane surfaces. Our density functional theory results suggest that even when geometric aspects are properly described by the simplifications considered, electronic properties might be very different, especially when edge atoms are not properly taken into account in the nanoparticle`s modeling. In addition, we propose a protocol that might help future theoretical descriptions of nanoparticles.

Analyzing the n→π* Electronic Transition of Formaldehyde in Water. A Sequential Monte Carlo/Time-Dependent Density Functional Theory

The n→π* absorption transition of formaldehyde in water is analyzed using combined and sequential classical Monte Carlo (MC) simulations and quantum mechanics (QM) calculations. MC simulations generate the liquid solute-solvent structures for subsequent QM calculations. Using time-dependent density functional theory in a localized set of gaussian basis functions (TD-DFT/6-311++G(d,p)) calculations are made on statistically relevant configurations to obtain the average solvatochromic shift. All results presented here use the electrostatic embedding of the solvent. The statistically converged average result obtained of 2300 cm-1 is compared to previous theoretical results available. Analysis is made of the effective dipole moment of the hydrogen-bonded shell and how it could be held responsible for the polarization of the solvent molecules in the outer solvation shells.

Density functional theory study of Fe(3)Ga

First-principles scalar relativistic calculations in supercells of 16 atoms are used to represent disordered B2 ordering of Fe(3)Ga in order to observe the effect of Ga-Ga pairs on the electronic structure of this alloy. From a comparison with pure bcc Fe it is observed that the energy position and occupation of e(g) and t(2g) states are largely affected by the Ga-Ga pairs and strengthened intraplane interactions takes place. The results show that a larger hybridization of the conduction band is in the source of the magnetostriction enhancement experimentally observed in Galfenol. (C) 2011 American Institute of Physics. [doi:10.1063/1.3525609]

Edge effects in bilayer graphene nanoribbons: Ab initio total-energy density functional theory calculations

We show that the ground state of zigzag bilayer graphene nanoribbons is nonmagnetic. It also possesses a finite gap, which has a nonmonotonic dependence with the width as a consequence of the competition between bulk and strongly attractive edge interactions. All results were obtained using ab initio total-energy density functional theory calculations with the inclusion of parametrized van der Waals interactions.

Approximation to density functional theory for the calculation of band gaps of semiconductors

The local-density approximation (LDA) together with the half occupation (transitionstate) is notoriously successful in the calculation of atomic ionization potentials. When it comes to extended systems, such as a semiconductor infinite system, it has been very difficult to find a way to half ionize because the hole tends to be infinitely extended (a Bloch wave). The answer to this problem lies in the LDA formalism itself. One proves that the half occupation is equivalent to introducing the hole self-energy (electrostatic and exchange correlation) into the Schrodinger equation. The argument then becomes simple: The eigenvalue minus the self-energy has to be minimized because the atom has a minimal energy. Then one simply proves that the hole is localized, not infinitely extended, because it must have maximal self-energy. Then one also arrives at an equation similar to the self- interaction correction equation, but corrected for the removal of just 1/2 electron. Applied to the calculation of band gaps and effective masses, we use the self- energy calculated in atoms and attain a precision similar to that of GW, but with the great advantage that it requires no more computational effort than standard LDA.

Spin gaps and spin-flip energies in density-functional theory

Energy gaps are crucial aspects of the electronic structure of finite and extended systems. Whereas much is known about how to define and calculate charge gaps in density-functional theory (DFT), and about the relation between these gaps and derivative discontinuities of the exchange-correlation functional, much less is known about spin gaps. In this paper we give density-functional definitions of spin-conserving gaps, spin-flip gaps and the spin stiffness in terms of many-body energies and in terms of single-particle (Kohn-Sham) energies. Our definitions are as analogous as possible to those commonly made in the charge case, but important differences between spin and charge gaps emerge already on the single-particle level because unlike the fundamental charge gap spin gaps involve excited-state energies. Kohn-Sham and many-body spin gaps are predicted to differ, and the difference is related to derivative discontinuities that are similar to, but distinct from, those usually considered in the case of charge gaps. Both ensemble DFT and time-dependent DFT (TDDFT) can be used to calculate these spin discontinuities from a suitable functional. We illustrate our findings by evaluating our definitions for the Lithium atom, for which we calculate spin gaps and spin discontinuities by making use of near-exact Kohn-Sham eigenvalues and, independently, from the single-pole approximation to TDDFT. The many-body corrections to the Kohn-Sham spin gaps are found to be negative, i.e., single-particle calculations tend to overestimate spin gaps while they underestimate charge gaps.

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.

SiO(2) in density functional theory and beyond

We present the first-principle electronic structure calculation on an amorphous material including many-body corrections within the GW approximation. We show that the inclusion of the local field effects in the exchange-correlation potential is crucial to quantitatively describe amorphous systems and defect states. We show that the mobility gap of amorphous silica coincides with the band gap of quartz, contrary to the traditional picture and the densityfunctional theory results. (C) 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Theoretical investigation of electronic and optical properties of andalusite within density functional theory

The electronic and optical properties of andalusite were studied by using quantum-mechanical calculations based on the density functional theory (DFT). The electronic structure shows that andalusite has a direct band gap of 5.01 eV. The complex dielectric function and optical constants, such as extinction coefficient, refractive index, reflectivity and energy-loss spectrum, are calculated. The optical properties of andalusite are discussed based on the band structure calculations. It is shown that the O-2p states and Al-3s states play a major role in optical transitions as initial and final states, respectively. (C) 2010 Elsevier Ltd. All rights reserved.

Orbital-polarization terms: From a phenomenological to a first-principles description of orbital magnetism in density-functional theory

Phenomenological orbital-polarizition (OP) terms have been repeatedly introduced in the single-particle equations of spin-density-functional theory, in order to improve the description of orbital magnetic moments in systems containing transition metal ions. Here we show that these ad hoc corrections can be interpreted as approximations to the exchange-correlation vector potential A(xc) of current-density functional theory (CDFT). This connection provides additional information on both approaches: phenomenological OP terms are connected to first-principles theory, leading to a rationale for their empirical success and a reassessment of their limitations and the approximations made in their derivation. Conversely, the connection of OP terms with CDFT leads to a set of simple approximations to the CDFT potential A(xc), with a number of desirable features that are absent from electron-gas-based functionals. (C) 2008 Wiley Periodicals, Inc.

Detailed analysis of the charge transfer complex N,N-dimethylaniline-SO(2) by Raman spectroscopy and density functional theory calculations

Although the amine sulfur dioxide chemistry was well characterized in the past both experimentally and theoretically, no systematic Raman spectroscopic study describes the interaction between N,N-dimethylaniline (DMA) and sulfur dioxide (SO(2)). The formation of a deep red oil by the reaction of SO(2) with DMA is an evidence of the charge transfer (CT) nature of the DMA-SO(2) interaction. The DMA -SO(2) normal Raman spectrum shows the appearance of two intense bands at 1110 and 1151 cm(-1), which are enhanced when resonance is approached. These bands are assigned to nu(s)(SO(2)) and nu(phi-N) vibrational modes, respectively, confirming the interaction between SO(2) and the amine via the nitrogen atom. The dimethyl group steric effect favors the interaction of SO(2) with the ring pi electrons, which gives rise to a pi-pi* low-energy CT electronic transition, as confirmed by time-dependent density functional theory (TDDFT) calculations. In addition, the calculated Raman DMA-SO(2) spectrum at the B3LYP/6-311++g(3df,3pd) level shows good agreement with the experimental results (vibrational wavenumbers and relative intensities), allowing a complete assignment of the vibrational modes. A better understanding of the intermolecular interactions in this model system can be extremely useful in designing new materials to absorb, detect, or even quantify SO(2). Copyright (C) 2009 John Wiley & Sons, Ltd.

FT-Raman, FTIR and density functional theory studies of a hydrogen-bonded formamide: pyridine complex

Raman and IR experiments have been carried out on formamide (FA) and pyridine (Py) mixtures at different compositions. The appearance of a new Raman band at 996 cm(-1) (nu(1) region of Py), whose intensity depends on the FA concentration, is assigned to an FA: Py adduct and this result is in excellent agreement with those of other authors who employed noisy light-based coherent Raman scattering spectroscopy (I((2)) CARS). Another band at 1587 cm(-1) (nu(8) region of Py) has been observed for the first time by using Raman and IR spectroscopies. Its intensity shows the same dependence on the FA concentration and this fact allows us to also attribute it to an FA: Py adduct. The good relationship between the Raman and IR data demonstrates the potential of the vibrational spectroscopy for this kind of study. Owing to higher absolute Raman scattering cross section, the nu(1) region of Py has been chosen for the quantitative analysis and a stoichiometry of 1 : 1 FA: Py is reported. The experimental data are very well supported by the density functional theory (OFT) calculation, which was employed for the first time to the present system. Furthermore, the actual investigation shows an excellent agreement with those reported from computational calculations for similar systems. A comparison with our previous studies confirms that: the solvent dielectric constant determines the stoichiometry of a given Lewis acid-base adduct in the infinite dilution limit. Copyright (C) 2009 John Wiley & Sons, Ltd.