On the equation of state and atomic mean square displacement of crystals

We have presented a Green's function method for the calculation of the atomic mean square displacement (MSD) for an anharmonic Hamil toni an . This method effectively sums a whole class of anharmonic contributions to MSD in the perturbation expansion in the high temperature limit. Using this formalism we have calculated the MSD for a nearest neighbour fcc Lennard Jones solid. The results show an improvement over the lowest order perturbation theory results, the difference with Monte Carlo calculations at temperatures close to melting is reduced from 11% to 3%. We also calculated the MSD for the Alkali metals Nat K/ Cs where a sixth neighbour interaction potential derived from the pseudopotential theory was employed in the calculations. The MSD by this method increases by 2.5% to 3.5% over the respective perturbation theory results. The MSD was calculated for Aluminum where different pseudopotential functions and a phenomenological Morse potential were used. The results show that the pseudopotentials provide better agreement with experimental data than the Morse potential. An excellent agreement with experiment over the whole temperature range is achieved with the Harrison modified point-ion pseudopotential with Hubbard-Sham screening function. We have calculated the thermodynamic properties of solid Kr by minimizing the total energy consisting of static and vibrational components, employing different schemes: The quasiharmonic theory (QH), ).2 and).4 perturbation theory, all terms up to 0 ().4) of the improved self consistent phonon theory (ISC), the ring diagrams up to o ().4) (RING), the iteration scheme (ITER) derived from the Greens's function method and a scheme consisting of ITER plus the remaining contributions of 0 ().4) which are not included in ITER which we call E(FULL). We have calculated the lattice constant, the volume expansion, the isothermal and adiabatic bulk modulus, the specific heat at constant volume and at constant pressure, and the Gruneisen parameter from two different potential functions: Lennard-Jones and Aziz. The Aziz potential gives generally a better agreement with experimental data than the LJ potential for the QH, ).2, ).4 and E(FULL) schemes. When only a partial sum of the).4 diagrams is used in the calculations (e.g. RING and ISC) the LJ results are in better agreement with experiment. The iteration scheme brings a definitive improvement over the).2 PT for both potentials.

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.

Effects of Side-Chain and Electron Exchange Correlation on the Band Structure of Perylene Diimide Liquid Crystals: A Density Functional Study

The structural and electronic properties of perylene diimide liquid crystal PPEEB are studied using ab initio methods based on the density functional theory (I)FT). Using available experimental crystallographic data as a guide, we propose a detailed structural model for the packing of solid PPEEB. We find that due to the localized nature of the band edge wave function, theoretical approaches beyond the standard method, such as hybrid functional (PBE0), are required to correctly characterize the band structure of this material. Moreover, unlike previous assumptions, we observe the formation of hydrogen bonds between the side chains of different molecules, which leads to a dispersion of the energy levels. This result indicates that the side chains of the molecular crystal not only are responsible for its structural conformation but also can be used for tuning the electronic and optical properties of these materials.

sigma- and pi-defects at graphene nanoribbon edges: Building spin filters

The presence of certain kinds of defects at the edges of monohydrogenated zigzag graphene nanoribbons changes dramatically the charge transport properties inducing a spin-polarized conductance. Using an approach based on density functional theory and nonequilibrium Green`s function formalism to calculate the transmittance, we classify the defects in different classes depending on their distinct transport properties: (i) sigma-defects, which do not affect the transmittance close to the Fermi energy (E(F)); and (ii) pi-defects, which cause a spin polarization of the transmittance and that can be further divided into either electron or hole defects if the spin transport polarization results in larger transmittance for the up or down spin channel, respectively.

Sistemas Nanoestruturados: Heteroestruturas Quasi-Periódicas de Nitretos e Cálculos Ab initio em Polimorfos CaCO3

The physical properties and the excitations spectrum in oxides and semiconductors materials are presented in this work, whose the first part presents a study on the confinement of optical phonons in artificial systems based on III-V nitrides, grown in periodic and quasiperiodic forms. The second part of this work describes the Ab initio calculations which were carried out to obtain the optoeletronic properties of Calcium Oxide (CaO) and Calcium Carbonate (CaCO3) crystals. For periodic and quasi-periodic superlattices, we present some dynamical properties related to confined optical phonons (bulk and surface), obtained through simple theories, such as the dielectric continuous model, and using techniques such as the transfer-matrix method. The localization character of confined optical phonon modes, the magnitude of the bands in the spectrum and the power laws of these structures are presented as functions of the generation number of sequence. The ab initio calculations have been carried out using the CASTEP software (Cambridge Total Sequential Energy Package), and they were based on ultrasoft-like pseudopotentials and Density Functional Theory (DFT). Two di®erent geometry optimizations have been e®ectuated for CaO crystals and CaCO3 polymorphs, according to LDA (local density approximation) and GGA (generalized gradient approximation) approaches, determining several properties, e. g. lattice parameters, bond length, electrons density, energy band structures, electrons density of states, e®ective masses and optical properties, such as dielectric constant, absorption, re°ectivity, conductivity and refractive index. Those results were employed to investigate the confinement of excitons in spherical Si@CaCO3 and CaCO3@SiO2 quantum dots and in calcium carbonate nanoparticles, and were also employed in investigations of the photoluminescence spectra of CaCO3 crystal

Molecular and Electronic Structure of the Peptide Subunit of Geobacter sulfurreducens Conductive Phi from First Principles

The respiration of metal oxides by the bacterium Geobacter sulfurreducens requires the assembly of a small peptide (the GS pilin) into conductive filaments termed pili. We gained insights into the contribution of the GS pilin to the pilus conductivity by developing a homology model and performing molecular dynamics simulations of the pilin peptide in vacuo and in solution. The results were consistent with a predominantly helical peptide containing the conserved a-helix region required for pilin assembly but carrying a short carboxy-terminal random-coiled segment rather than the large globular head of other bacterial pilins. The electronic structure of the pain was also explored from first principles and revealed a biphasic charge distribution along the pilin and a low electronic HOMO-LUMO gap, even in a wet environment. The low electronic band gap was the result of strong electrostatic fields generated by the alignment of the peptide bond dipoles in the pilin's alpha-helix and by charges from ions in solution and amino acids in the protein. The electronic structure also revealed some level of orbital delocalization in regions of the pilin containing aromatic amino acids and in spatial regions of high resonance where the HOMO and LUMO states are, which could provide an optimal environment for the hopping of electrons under thermal fluctuations. Hence, the structural and electronic features of the pilin revealed in these studies support the notion of a pilin peptide environment optimized for electron conduction.

Low-energy electron collisions with glycine

We report cross sections for elastic electron scattering by gas phase glycine (neutral form), obtained with the Schwinger multichannel method. The present results are the first obtained with a new implementation that combines parallelization with OpenMP directives and pseudopotentials. The position of the well known pi* shape resonance ranged from 2.3 eV to 2.8 eV depending on the polarization model and conformer. For the most stable isomer, the present result (2.4 eV) is in fair agreement with electron transmission spectroscopy assignments (1.93 +/- 0.05 eV) and available calculations. Our results also point out a shape resonance around 9.5 eV in the A' symmetry that would be weakly coupled to vibrations of the hydroxyl group. Since electron attachment to a broad and lower lying sigma* orbital located on the OH bond has been suggested the underlying mechanism leading to dissociative electron attachment at low energies, we sought for a shape resonance around similar to 4 eV. Though we obtained cross sections with the target molecule at the equilibrium geometry and with stretched OH bond lengths, least-squares fits to the calculated eigenphase sums did not point out signatures of this anion state (though, in principle, it could be hidden in the large background). The low energy (similar to 1 eV) integral cross section strongly scales as the bond length is stretched, and this could indicate a virtual state pole, since dipole supported bound states are not expected at the geometries addressed here. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3687345]

Changes in a nanoparticle's spectroscopic signal mediated by the local environment

Using a first-principles theoretical model the adsorption of a methyl radical on different sized silver nanoparticles is compared to the adsorption of the same radical on model surfaces. Calculations of our structural, dynamical and electronic properties indicated that small changes in the local environment will lead to small changes in infrared (IR) wavenumbers, but in dramatic changes in the IR signal. Our calculations indicate the lower the adsorption site coordination, the higher is the signal strength, suggesting that small changes in the electronic charge distribution will result in bigger changes in the polarizability and hence in the spectroscopic signal intensity. This effect explains, among others, the signal magnification observed for nanoparticles in surface enhanced Raman spectroscopic (SERS) experiments.

Quantum Monte Carlo study of small aluminum clusters Al-n (n=2-13)

Using fixed node diffusion quantum Monte Carlo (FN-DMC) simulations and density functional theory (DFT) within the generalized gradient approximations, we calculate the total energies of the relaxed and unrelaxed neutral, cationic, and anionic aluminum clusters, Al-n (n = 1-13). From the obtained total energies, we extract the ionization potential and electron detachment energy and compare with previous theoretical and experimental results. Our results for the electronic properties from both the FN-DMC and DFT calculations are in reasonably good agreement with the available experimental data. A comparison between the FN-DMC and DFT results reveals that their differences are a few tenths of electron volt for both the ionization potential and the electron detachment energy. We also observe two distinct behaviors in the electron correlation contribution to the total energies from smaller to larger clusters, which could be assigned to the structural transition of the clusters from planar to three-dimensional occurring at n = 4 to 5.

CdSe/CdTe interface band gaps and band offsets calculated using spin-orbit and self-energy corrections

We performed ab initio calculations of the electronic structures of bulk CdSe and CdTe, and their interface band alignments on the CdSe in-plane lattice parameters. For this, we employed the LDA-1/2 self-energy correction scheme [L.G. Ferreira, M. Marques, L.K. Teles, Phys. Rev. B 78 (2008) 125116] to obtain corrected band gaps and band offsets. Our calculations include the spin-orbit effects for the bulk cases, which have shown to be of importance for the equilibrium systems and are possibly degraded in these strained semiconductors. Therefore, the SO showed reduced importance for the band alignment of this particular system. Moreover, the electronic structure calculated along the transition region across the CdSe/CdTe interface shows an interesting non-monotonic variation of the band gap in the range 0.8-1.8 eV, which may enhance the absorption of light for corresponding frequencies at the interface between these two materials in photovoltaic applications. (C) 2012 Elsevier B.V. All rights reserved.

Shape resonance spectra of lignin subunits

We report integral cross sections for elastic electron scattering by the lignin subunits phenol, guaiacol, and p-coumaryl alcohol. Our calculations employed the Schwinger multichannel method with pseudopotentials and indicate three to four pi* shape resonances for each of these systems, suggesting that low-energy electrons could efficiently transfer energy into the lignin matrix. We also discuss dissociation mechanisms based on the calculated cross sections, available experimental data, virtual orbital analysis, and the knowledge on electron interactions with biomolecules. Our results point out a physical-chemical basis for electron-driven biomass delignification. The latter would be an essential step for efficient biofuel production from lignocellulosic materials.

Electron collisions with the HCOOH...(H2O)n (n=1, 2 and 3) complexes in liquid phase: the influence on the π* shape resonance of formic acid

In this conference we report cross sections for elastic collisions of low-energy electrons with the HCOOH…(H2O)n complexes, with n = 1, 2 and 3. The scattering cross sections were computed with the Schwinger multichannel method [K. Takatsuka and V. McKoy, Phys. Rev. A 24 , 2473 (1981); Phys. Rev. A 30 , 1734 (1984)] with pseudopotentials [M. H. F. Bettega, L. G. Ferreira, and M. A. P. Lima, Phys. Rev. A 47, 1111 (1993)] in the static-exchange and static-exchange plus polarization approximations, for energies from 0.5 eV to 6 eV. We considered some diÆerent hydrogen-bonded structures for the complexes that were generated with classical Monte Carlo simulations [K. Coutinho and S. Canuto, J. Chem. Phys. 113, 9132, (2000)]. The aim of this work is to investigate the effect of the surrounding water molecules on the π* shape resonance of the solute. Previous theoretical and experimental studies carried out in the gas phase reported a π* state for HCOOH at around 1.9 eV. For the n = 1 case and for all complexes, the stabilization of the resonance was observed (it appears at lower energy compared to the value obtained in the gas phase), as reported previously for the CH2O…H2O complexes [T. C. Freitas, M. A. P. Lima, S. Canuto, and M. H. F. Bettega, Phys. Rev. A 80, 062710 (2009)]. This result indicates that the presence of the solvent may affect the processes related to the π* state, such as the molecular dissociation by electron impact. For the n = 2 case we have observed both stabilization and destabilization of the π* resonance, that is associated with the hydrogen bond donor or acceptor role of the water molecules in the complexes. For the n = 3 case, preliminary static-exchange results show the stabilization of the π* state. We propose an explanation of the stabilization/destabilization of the π* state in terms of the polarization of the solute due to the surrounding water molecules and the net charge in the solute.

Dissociative electron attachment to chloromethane.

A low energy electron may attach to a molecule, forming a metastable resonance, which may dissociate into a stable anion and a neutral radical. Chloromethane has been a good target for dissociative electron attachment studies, since it is a small molecule with a clear dissociative ‘sigma*’ shape resonance. We present potential energy curves for CH3Cl and its anion, as a function of the C-Cl distance. Due to the resonant nature of the anion, a correct description requires a treatment based on scattering calculations. In order to compute elastic cross sections and phase shifts we employed the Schwinger multichannel method, implemented with pseudopotentials of Bachelet, Hamann and Schlüter, at the static-exchange plus polarization approximation. At the equilibrium geometry, the resonance was found arround 3.3 eV, in accordance to experience. The incoming electron is captured by a ‘sigma*’ orbital located at the C-Cl bond, which will relax in the presence of this extra electron. We took this bond as the reaction coordinate, and performed several scattering calculations for a series of nuclear conformations. The phase shift obtained in each calculation was fitted by a two component function, consisting in the usual Breit-Wigner profile, which captures the resonant character, and a second order polynomial in the wave number, which accounts for the background contribution. That way, we obtained position and width of the resonance, which allowed us to build the potential energy curve. For larger distances, the anion becomes stable and usual electronic structure calculations suffice. Furthermore, the existence of a dipole-bound anion state is revealed when we employed a set of very diffuse functions. The knowledge on the behaviour of the neutral and anionic electronic states helps us in elucidating how the dissociation takes place.

Shape resonance spectrum of cytosine-guanine pairs.

Single and double strand breaks in DNA can be caused by low-energy electrons, the most abundant secondary products of the interaction of ionizing radiation to the biological matter. Attachment of these electrons to biomolecules lead to the formation of transient negative ions (TNIs) [1], often referred to as resonances, a process that may lead to significant vibrational excitation and dissociation. In the present study, we employ the parallel version [2] of the Schwinger Multichannel Method implemented with pseudopotentials [3] to obtain the shape resonance spectrum of cytosine-guanine (CG) pairs, with special attention to π* transient anion states. Recent experimental studies pointed out a quasi-continuum vibrational excitation spectrum for electron collisions against formic acid dimers [4], suggesting that electron attachment into π* valence orbitals could induce proton transfer in these dimers. In addition, our previous studies on the shape resonance spectra of the hydrogen-bonded complexes comprising formic acid and formamide units indicated interesting electron delocalization (localization) effects arising from the presence (absence) of inversion symmetry centers in the complexes [5]. In the present work, we extend the studies on hydrogen-bonded complexes to the CG pair, where localization of ¼¤ anions would be expected, based on the previous results. References [1]. B. Boudaïffa, P. Cloutier, D. Hunting, M. A. Huels, L. Sanche, Science 287, 1658 (2000). [2]. J. S. dos Santos, R. F. da Costa , M. T. do N. Varella, J. Chem. Phys. 136, 084307 (2012). [3]. M. H. F. Bettega, L. G. Ferreira, M. A. P. Lima, Phys. Rev. A 47, 1111 (1993). [4]. M. Allan, Phys. Rev. Lett. 98, 123201 (2007). [5]. T. C. Freitas, S. dA. Sanchez, M. T. do N. Varella, M. H. F. Bettega, Phys. Rev. A 84, 062714 (2011).

Shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil.

We report on the shape resonance spectra of uracil, 5-fluorouracil, and 5-chlorouracil, as obtained from fixed-nuclei elastic scattering calculations performed with the Schwinger multichannel method with pseudopotentials. Our results are in good agreement with the available electron transmission spectroscopy data, and support the existence of three π* resonances in uracil and 5-fluorouracil. As expected, the anion states are more stable in the substituted molecules than in uracil. Since the stabilization is stronger in 5-chlorouracil, the lowest π* resonance in this system becomes a bound anion state. The present results also support the existence of a low-lying σ ∗ CCl shape resonance in 5- chlorouracil. Exploratory calculations performed at selected C–Cl bond lengths suggest that the σ ∗ CCl resonance could couple to the two lowest π* states, giving rise to a very rich dissociation dynamics. These facts would be compatible with the complex branching of the dissociative electron attachment cross sections, even though we cannot discuss any details of the vibration dynamics based only on the present fixed-nuclei results.