The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

Accurate Predictions of Water Cluster Formation, (H2O)n=2-10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Indium-cadmium-oxide films having exceptional electrical conductivity and optical transparency: Clues for optimizing transparent conductors

Materials with high electrical conductivity and optical transparency are needed for future flat panel display, solar energy, and other opto-electronic technologies. InxCd1-xO films having a simple cubic microstructure have been grown on amorphous glass substrates by a straightforward chemical vapor deposition process. The x = 0.05 film conductivity of 17,000 S/cm, carrier mobility of 70 cm2/Vs, and visible region optical transparency window considerably exceed the corresponding parameters for commercial indium-tin oxide. Ab initio electronic structure calculations reveal small conduction electron effective masses, a dramatic shift of the CdO band gap with doping, and a conduction band hybridization gap caused by extensive Cd 5s + In 5s mixing.

Epitaxially Grown Strained Pentacene Thin Film on Graphene Membrane

Organic-graphene system has emerged as a new platform for various applications such as flexible organic photovoltaics and organic light emitting diodes. Due to its important implication in charge transport, the study and reliable control of molecular packing structures at the graphene-molecule interface are of great importance for successful incorporation of graphene in related organic devices. Here, an ideal membrane of suspended graphene as a molecular assembly template is utilized to investigate thin-film epitaxial behaviors. Using transmission electron microscopy, two distinct molecular packing structures of pentacene on graphene are found. One observed packing structure is similar to the well-known bulk-phase, which adapts a face-on molecular orientation on graphene substrate. On the other hand, a rare polymorph of pentacene crystal, which shows significant strain along the c-axis, is identified. In particular, the strained film exhibits a specific molecular orientation and a strong azimuthal correlation with underlying graphene. Through ab initio electronic structure calculations, including van der Waals interactions, the unusual polymorph is attributed to the strong graphene-pentacene interaction. The observed strained organic film growth on graphene demonstrates the possibility to tune molecular packing via graphene-molecule interactions.

Modelling UV response of biomolecules: from electronic structure to ab-initio dynamics.

The interaction of organic chromophores with light initiates ultrafast processes in the timescale of femtoseconds. An atomistic understanding of the mechanism driving such photoinduced reactions opens up the door to exploit them for our benefit. This thesis studies the interactions of ultraviolet light with the DNA/RNA molecules and the amino-acid tryptophan. Using some of the most accurate electronic structure methods and sophisticated environmental modelling, the works documented herein enable quantitative comparisons with cutting-edge experimental data. The relaxation pathways undertaken by the excited molecule are revealed through static and dynamical investigations of the excited-state potential energy surface. The profound role played by the dynamic response of the environment to guide the excitation in these timescales is addressed thoroughly.

Electronic, vibrational and related properties of group IV metal oxides by ab initio calculations

We present our theoretical results for the structural, electronic, vibrational and optical properties of MO(2) (M = Sn, Zr, Hf and Ti) obtained by first-principles calculations. Relativistic effects are demonstrated to be important for a realistic description of the detailed structure of the electronic frequency-dependent dielectric function, as well as of the carrier effective masses. Based on our results, we found that the main contribution of the high values calculated for the oxides dielectric constants arises from the vibrational properties of these oxides, and the vibrational static dielectric constant values diminish with increasing pressure. (c) 2008 Elsevier B.V. All rights reserved.

Ab initio study of the electronic and optical properties of sillimanite (Al(2)SiO(5)) crystal

Ab initio calculations based on the density functional theory (DFT) are used to investigate the electronic and optical properties of sillimanite. The geometrical parameters of the unit cell, which contain 32 atoms, have been fully optimized and are in good agreement with the experimental data. The electronic structure shows that sillimanite has an indirect band gap of 5.18 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 sillimanite are discussed based on the band structure calculations. It is shown that the O-2p states and Al-3s, Si-3s states play the major role in optical transitions as initial and final states, respectively. (C) 2011 Elsevier B.V. All rights reserved.

Electronic and optical properties of grossular garnet (Ca(3)Al(2)Si(3)O(12)): An ab initio study

The electronic and optical properties of grossular garnet are investigated using density functional theory (DFT) within generalized gradient approximation (GGA). The calculated lattice parameters are in good agreement with the experiment data. The electronic structure shows that grossular has a direct band gap of 5.22 eV. The dielectric functions, reflective index, extinction coefficient, reflectivity and energy-loss spectrum are calculated. The optical properties of grossular are discussed based on the band structure calculations. The O 2p states and Si 3s play a major role in these optical transitions as initial and final states, respectively. The absorption spectrum is localized in the ultraviolet range between 30 and 250 nm. Finally, we concluded that pure grossular crystal does not absorb radiation in the visible range. (c) 2009 Elsevier B.V. All rights reserved.

Ab-initio determination of x-ray absorption near edge structure (xanes) spectra in an ultrasoft and norm conserving pseudopotentials scheme

X-ray absorption spectroscopy (XAS) is a powerful means of investigation of structural and electronic properties in condensed -matter physics. Analysis of the near edge part of the XAS spectrum, the so – called X-ray Absorption Near Edge Structure (XANES), can typically provide the following information on the photoexcited atom: - Oxidation state and coordination environment. - Speciation of transition metal compounds. - Conduction band DOS projected on the excited atomic species (PDOS). Analysis of XANES spectra is greatly aided by simulations; in the most common scheme the multiple scattering framework is used with the muffin tin approximation for the scattering potential and the spectral simulation is based on a hypothetical, reference structure. This approach has the advantage of requiring relatively little computing power but in many cases the assumed structure is quite different from the actual system measured and the muffin tin approximation is not adequate for low symmetry structures or highly directional bonds. It is therefore very interesting and justified to develop alternative methods. In one approach, the spectral simulation is based on atomic coordinates obtained from a DFT (Density Functional Theory) optimized structure. In another approach, which is the object of this thesis, the XANES spectrum is calculated directly based on an ab – initio DFT calculation of the atomic and electronic structure. This method takes full advantage of the real many-electron final wavefunction that can be computed with DFT algorithms that include a core-hole in the absorbing atom to compute the final cross section. To calculate the many-electron final wavefunction the Projector Augmented Wave method (PAW) is used. In this scheme, the absorption cross section is written in function of several contributions as the many-electrons function of the finale state; it is calculated starting from pseudo-wavefunction and performing a reconstruction of the real-wavefunction by using a transform operator which contains some parameters, called partial waves and projector waves. The aim of my thesis is to apply and test the PAW methodology to the calculation of the XANES cross section. I have focused on iron and silicon structures and on some biological molecules target (myoglobin and cytochrome c). Finally other inorganic and biological systems could be taken into account for future applications of this methodology, which could become an important improvement with respect to the multiscattering approach.

Spin filtering and disorder-induced magnetoresistance in carbon nanotubes: Ab initio calculations

Nitrogen-doped carbon nanotubes can provide reactive sites on the porphyrin-like defects. It is well known that many porphyrins have transition-metal atoms, and we have explored transition-metal atoms bonded to those porphyrin-like defects inN-doped carbon nanotubes. The electronic structure and transport are analyzed by means of a combination of density functional theory and recursive Green's function methods. The results determined the heme B-like defect (an iron atom bonded to four nitrogens) is the most stable and has a higher polarization current for a single defect. With randomly positioned heme B defects in nanotubes a few hundred nanometers long, the polarization reaches near 100%, meaning they are effective spin filters. A disorder-induced magnetoresistance effect is also observed in those long nanotubes, and values as high as 20 000% are calculated with nonmagnectic eletrodes.

Barrier-free substitutional doping of graphene sheets with boron atoms: Ab initio calculations

Using ab initio methods, we propose a simple and effective way to substitutionally dope graphene sheets with boron. The method consists of selectively exposing each side of the graphene sheet to different elements. We first expose one side of the membrane to boron while the other side is exposed to nitrogen. Proceeding this way, the B atoms will be spontaneously incorporated into the graphene membrane without any activation barrier. In a second step, the system should be exposed to a H-rich environment, which will remove the CN radical from the layer and form HCN, leading to a perfect substitutional doping.