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.

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.

A comparative study of carbon gasification with O-2 and CO2 by density functional theory calculations

A comparative study of carbon gasification with O-2 and CO2 was conducted by using density functional theory calculations. It was found that the activation energy and the number of active sites in carbon gasification reactions are significantly affected by both the capacity and manner of gas chemisorption. O-2 has a strong adsorption capacity and the dissociative chemisorption of O-2 is thermodynamically favorable on either bare carbon surface or even isolated edge sites. As a result, a large number of semiquinone and o-quinone oxygen can be formed indicating a significant increase in the number of active sites. Moreover, the weaker o-quinone C-C bonds can also drive the reaction forward at (ca. 30%) lower activation energy. Epoxy oxygen forms under relatively high O-2 pressure, and it can only increase the number of active sites, not further reduce the activation energy. CO2 has a lower adsorption capacity. Dissociative chemisorption of CO2 can only occur on two consecutive edge sites and o-quinone oxygen formed from CO2 chemisorption is negligible, let alone epoxy oxygen. Therefore, CO2-carbon reaction needs (ca 30%) higher activation energy. Furthermore, the effective active sites are also reduced by the manner Of CO2 chemisorption. A combination of the higher activation energy and the fewer active sites leads to the much lower reaction rate Of CO2-carbon.

Helium in polygonal nanopores at zero temperature: Density functional theory calculations

We investigate adsorption of helium in nanoscopic polygonal pores at zero temperature using a finite-range density functional theory. The adsorption potential is computed by means of a technique denoted as the elementary source method. We analyze a rhombic pore with Cs walls, where we show the existence of multiple interfacial configurations at some linear densities, which correspond to metastable states. Shape transitions and hysterectic loops appear in patterns which are richer and more complex than in a cylindrical tube with the same transverse area.

Helium in polygonal nanopores at zero temperature: Density functional theory calculations

We investigate adsorption of helium in nanoscopic polygonal pores at zero temperature using a finite-range density functional theory. The adsorption potential is computed by means of a technique denoted as the elementary source method. We analyze a rhombic pore with Cs walls, where we show the existence of multiple interfacial configurations at some linear densities, which correspond to metastable states. Shape transitions and hysterectic loops appear in patterns which are richer and more complex than in a cylindrical tube with the same transverse area.

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.

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

Cobalt incorporation in calcite: Thermochemistry of (Ca,Co)CO3 solid solutions from density functional theory simulations

The incorporation of cobalt in mixed metal carbonates is a possible route to the immobilization of this toxic element in the environment. However, the thermodynamics of (Ca,Co)CO3 solid solutions are still unclear due to conflicting data from experiment and from the observation of natural ocurrences. We report here the results of a computer simulation study of the mixing of calcite (CaCO3) and spherocobaltite (CoCO3), using density functional theory calculations. Our simulations suggest that previously proposed thermodynamic models, based only on the range of observed compositions, significantly overestimate the solubility between the two solids and therefore underestimate the extension of the miscibility gap under ambient conditions. The enthalpy of mixing of the disordered solid solution is strongly positive and moderately asymmetric: calcium incorporation in spherocobaltite is more endothermic than cobalt incorporation in calcite. Ordering of the impurities in (0001) layers is energetically favourable with respect to the disordered solid solution at low temperatures and intermediate compositions, but the ordered phase is still unstable to demixing. We calculate the solvus and spinodal lines in the phase diagram using a sub-regular solution model, and conclude that many Ca1-xCoxCO3 mineral solid solutions (with observed compositions of up to x=0.027, and above x=0.93) are metastable with respect to phase separation. We also calculate solid/aqueous distribution coefficients to evaluate the effect of the strong non-ideality of mixing on the equilibrium with aqueous solution, showing that the thermodynamically-driven incorporation of cobalt in calcite (and of calcium in spherocobaltite) is always very low, regardless of the Co/Ca ratio of the aqueous environment.

Intercalation processes and diffusion paths of lithium ions in spinel-type structured Li(1+x)Ti(2)O(4): Density functional theory study

Intercalation processes and corresponding diffusion paths of Li ions into spinel-type structured Li(1+x)Ti(2)O(4) (0 <= x <= 0.375) are systematically studied by means of periodic density functional theory calculations for different compositions and arrangements. An analysis of the site preference for intercalation processes is carried out, while energy barriers for the diffusion paths have been computed in detail. Our results indicate that the Li insertion is thermodynamically favorable at octahedral sites 16c in the studied composition range, and Li migration from tetrahedral sites 8a to octahedral sites 16c stabilizes the structure and becomes favorable for compositions x >= 0.25. Diffusion paths from less stable arrangements involving Li migrations between tetrahedral and octahedral sites exhibit the lowest energy barrier since the corresponding trajectories and energy profiles take place across a triangle made by three neighboring oxygen anions without structural modification. Theoretical and experimental diffusion coefficients are in reasonable agreement.

Interactions between 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid and gamma-Al2O3 (100) surface calculated by density functional theory

The main aim of this work is to investigate the 1-butyl-3-methylimidazolium tetrafluoroborate ([C4C1Im]+[BF4]-) ionic liquid (IL) adsorption on the gamma-Al2O3 (100) by density functional theory calculations to try to rationalize the adsorption as an electrostatic phenomenon. Optimized geometries and interaction energies of IL one-monolayer on the gamma-Al2O3 were obtained on high surface coverage (one cationanion pair per 94.96 nm2). A study of dispersion force was made to estimate its contribution to the adsorption. Overall, the process is ruled by electrostatic interaction between ions and surface. Adsorption of the anion [BF4]- and cation [C4C1Im]+ was also studied by Bader charge analysis and charge density difference for supported and unsupported situations. It is suggested that the IL ions have their charges maintained with significant anion cloud polarization inward to the acid aluminum sites. (c) 2012 Wiley Periodicals, Inc.

Computation of (3)J(HH) coupling constants with a combination of density functional theory and semiempirical calculations. Application to complex molecules

This work evaluates the efficiency of economic levels of theory for the prediction of (3)J(HH) spin-spin coupling constants, to be used when robust electronic structure methods are prohibitive. To that purpose, DFT methods like mPW1PW91. B3LYP and PBEPBE were used to obtain coupling constants for a test set whose coupling constants are well known. Satisfactory results were obtained in most of cases, with the mPW1PW91/6-31G(d,p)//B3LYP/6-31G(d,p) leading the set. In a second step. B3LYP was replaced by the semiempirical methods PM6 and RM1 in the geometry optimizations. Coupling constants calculated with these latter structures were at least as good as the ones obtained by pure DFT methods. This is a promising result, because some of the main objectives of computational chemistry - low computational cost and time, allied to high performance and precision - were attained together. (C) 2012 Elsevier B.V. All rights reserved.

Computing the structural and vibrational properties of polymorphic organic molecular crystals through van der Waals corrected density functional theory and the electronic properties of organic thin films through microelectrostatic calculations.

The present Thesis reports on the various research projects to which I have contributed during my PhD period, working with several research groups, and whose results have been communicated in a number of scientific publications. The main focus of my research activity was to learn, test, exploit and extend the recently developed vdW-DFT (van der Waals corrected Density Functional Theory) methods for computing the structural, vibrational and electronic properties of ordered molecular crystals from first principles. A secondary, and more recent, research activity has been the analysis with microelectrostatic methods of Molecular Dynamics (MD) simulations of disordered molecular systems. While only very unreliable methods based on empirical models were practically usable until a few years ago, accurate calculations of the crystal energy are now possible, thanks to very fast modern computers and to the excellent performance of the best vdW-DFT methods. Accurate energies are particularly important for describing organic molecular solids, since they often exhibit several alternative crystal structures (polymorphs), with very different packing arrangements but very small energy differences. Standard DFT methods do not describe the long-range electron correlations which give rise to the vdW interactions. Although weak, these interactions are extremely sensitive to the packing arrangement, and neglecting them used to be a problem. The calculations of reliable crystal structures and vibrational frequencies has been made possible only recently, thanks to development of some good representations of the vdW contribution to the energy (known as “vdW corrections”).