Amorphous HfO(2) and Hf(1-x)Si(x)O via a melt-and-quench scheme using ab initio molecular dynamics

We have performed ab initio molecular dynamics simulations to generate an atomic structure model of amorphous hafnium oxide (a-HfO(2)) via a melt-and-quench scheme. This structure is analyzed via bond-angle and partial pair distribution functions. These results give a Hf-O average nearest-neighbor distance of 2.2 angstrom, which should be compared to the bulk value, which ranges from 1.96 to 2.54 angstrom. We have also investigated the neutral O vacancy and a substitutional Si impurity for various sites, as well as the amorphous phase of Hf(1-x)Si(x)O(2) for x=0.25, 0375, and 0.5.

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.

Crystal and molecular structures of three 2-sulphur-substituted cyclohexanones studied by X-ray crystallography and by ab initio molecular orbital calculations

The conformational features of three 2-sulphur-substituted cyclohexanone derivatives, which differ in the number of sulphur-bound oxygen atoms, i.e. zero (I), one (II) and two (III), were investigated by single crystal X-ray crystallography and geometry optimized structures determined using Hartree-Fock method. In each of (I)-(III) an intramolecular S center dot center dot center dot O(carbonyl) interaction is found with the magnitude correlated with the oxidation state of the sulphur atom, i.e. 2.838(3) angstrom in (I) to 2.924(2) angstrom in (II) to 3.0973(18) angstrom in (III). There is an inverse relationship between the strength of this interaction and the magnitude of the carbonyl bond. The supramolecular aggregation patterns are primarily determined by C-H center dot center dot center dot O contacts and are similarly influenced by the number of oxygen atoms in the molecular structures. Thus, a supramolecular chain is found in the crystal structure of (I). With an additional oxygen atom available to participate in C-H center dot center dot center dot O interactions, as in (II), a two-dimensional array is found. Finally, a three-dimensional network is found for (III). Despite there being differences in conformations between the experimental structures and those calculated in the gas-phase, the S center dot center dot center dot O interactions persist. The presence of intermolecular C-H center dot center dot center dot O interactions involving the cyclohexanone-carbonyl group in the solid-state, disrupts the stabilising intramolecular C-H center dot center dot center dot O interaction in the energetically-favoured conformation. (I): C(12)H(13)NO(3)S, triclinic space group P (1) over bar with a = 5.392(3) angstrom b = 10.731(6) angstrom, c = 11.075(6) angstrom, alpha = 113.424(4)degrees, beta = 94.167(9)degrees, gamma = 98.444(6)degrees, V = 575.5(6) angstrom(3), Z = 2, R(1) = 0.052; (II): C(12)H(13)NO(4)S, monoclinic P2(1)/n, a = 7.3506(15) angstrom, b = 6.7814(14) angstrom, c = 23.479(5) angstrom, beta = 92.94(3)degrees, V = 1168.8(4) angstrom(3), Z = 4, R(1) = 0.046; (III): C(12)H(13)NO(5)S, monoclinic P2(1)/c, a = 5.5491(11) angstrom, b = 24.146(3) angstrom, c = 11.124(3) angstrom, beta = 114.590(10)degrees, V = 1355.3(5) angstrom(3), Z = 4, R(1) = 0.051.

Formation of atomic carbon chains from graphene nanoribbons

The formation of one-dimensional carbon chains from graphene nanoribbons is investigated using ab initio molecular dynamics. We show under what conditions it is possible to obtain a linear atomic chain via pulling of the graphene nanoribbons. The presence of dimers composed of two-coordinated carbon atoms at the edge of the ribbons is necessary for the formation of the linear chains, otherwise there is simply the full rupture of the structure. The presence of Stone-Wales defects close to these dimers may lead to the formation of longer chains. The local atomic configuration of the suspended atoms indicates the formation of single and triple bonds, which is a characteristic of polyynes.

Atomistic molecular dynamics study of interface formation: Al on poly(p-phenylene vinylene)

We model interface formation by metal deposition on the conjugated polymer poly-para-phenylene vinylene, studying direct aluminum and layered aluminum-calcium structures Al/PPV and Al/Ca/PPV. To do that we use classical molecular dynamics simulations, checked by ab initio density-functional theory calculations, for selected relevant configurations. We find that Al not only migrates easily into the film, with a strong charge transfer to the neighboring chains, but also promotes rearrangement of the polymer in the interfacial region to the hexagonal structure. On the other hand, our results indicate that a thin Ca layer is sufficient to protect the film and maintain a well-defined metal/polymer interface, and that also a thin Al capping layer may protect the whole from environmental degradation.

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

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]

Ab initio study of the EFG at the N sites in imidazole

We study the nuclear quadrupole interaction at the nitrogen sites in the molecular and crystalline phases of the imidazole compound. We use PAW which is a state-of-the-art method to calculate the electronic structure and electric field gradient at the nucleus in the framework of the density functional theory. The quadrupole frequencies at both imino and amino N sites are in excellent agreement with measurements. This is the first time that the electric field gradient at crystalline imidazole is correctly treated by an ab initio theoretical approach.

Ab Initio Calculations of Structural Evolution and Conductance of Benzene-1,4-dithiol on Gold Leads

By performing at) initio density functional theory (DFT) calculations and electronic transport simulations based on the OFT nonequilibrium Green`s functions method we investigate how the conformational changes of a benzene-1,4-dithiol molecule bonded to gold affect the molecular transport as the electrodes are separated from each other. In particular we consider the full evolution of the stretching process until the Junction breaking point and compare results obtained with a standard semilocal exchange and correlation functional to those computed with a self-interaction corrected method. We conclude that the inclusion of self-interaction corrections is fundamental for describing both the molecule conductance and its stability against conformational fluctuations.

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.

Large Changes of Static Electric Properties Induced by Hydrogen Bonding: An ab Initio Study of Linear HCN Oligomers

We report the partitioning of the interaction-induced static electronic dipole (hyper)polarizabilities for linear hydrogen cyanide complexes into contributions arising from various interaction energy terms. We analyzed the nonadditivities of the studied properties and used these data to predict the electric properties of an infinite chain. The interaction-induced static electric dipole properties and their nonadditivities were analyzed using an approach based on numerical differentiation of the interaction energy components estimated in an external electric field. These were obtained using the hybrid variational-perturbational interaction energy decomposition scheme, augmented with coupled-cluster calculations, with singles, doubles, and noniterative triples. Our results indicate that the interaction-induced dipole moments and polarizabilities are primarily electrostatic in nature; however, the composition of the interaction hyperpolarizabilities is much more complex. The overlap effects substantially quench the contributions due to electrostatic interactions, and therefore, the major components are due to the induction and exchange induction terms, as well as the intramolecular electron-correlation corrections. A particularly intriguing observation is that the interaction first hyperpolarizability in the studied systems not only is much larger than the corresponding sum of monomer properties, but also has the opposite sign. We show that this effect can be viewed as a direct consequence of hydrogen-bonding interactions that lead to a decrease of the hyperpolarizability of the proton acceptor and an increase of the hyperpolarizability of the proton donor. In the case of the first hyperpolarizability, we also observed the largest nonadditivity of interaction properties (nearly 17%) which further enhances the effects of pairwise interactions.

Ab Initio Analysis of Monomers and Dimers of Trialkylphosphine Oxides: Structural and Thermodynamic Stability

Structural and thermodynamic stabilities of monomers and dimers of trialkylphosphine oxides (TRPO) were Studied using quantum chemistry calculations. Density functional theory calculations were carried Out and the structures Of four TRPO have been determined: TMPO (methyl; R = CH(3)), TEPO (ethyl; R = CH(3)CH(2)), TBPO (n-butyl; R = CH(3)(CH(2))(3)), and TOPO (n-octyl; R = CH(3)(CH(2))(7)). TRPO homodimers were investigated considering two isomeric possibilities for each dimer. Relative binding energies and the enthalpic and entropic contributions to the Gibbs free energy were Calculated for all dimers. The formation of dimers from the individual monomeric TRPO species as a function of temperature was also analyzed. (C) 2008 Wiley Periodicals, Inc. Int J Quantum Chem 109: 250-258, 2009

The role of disulfide bridges in the 3-D structures of the antimicrobial peptides gomesin and protegrin-1: a molecular dynamics study

Some antimicrobial peptides have a broad spectrum of action against many different kinds of microorganisms. Gomesin and protegrin-1 are examples of such antimicrobial peptides, and they were studied by molecular dynamics in this research. Both have a beta-hairpin conformation stabilized by two disulfide bridges and are active against Gram-positive and Gram-negative bacteria, as well as fungi. In this study, the role of the disulfide bridge in the maintenance of the tertiary peptide structure of protegrin-1 and gomesin is analyzed by the structural characteristics of these peptides and two of their respective variants, gomy4 and proty4, in which the four cysteines are replaced by four tyrosine residues. The absence of disulfide bridges in gomy4 and proty4 is compensated by overall reinforcement of the original hydrogen bonds and extra attractive interactions between the aromatic rings of the tyrosine residues. The net effects on the variants with respect to the corresponding natural peptides are: i) maintenance of the original beta-hairpin conformation, with great structural similarities between the mutant and the corresponding natural peptide; ii) combination of positive F and. Ramachandran angles within the hairpin head region with a qualitative change to a combination of positive (F) and negative (.) angles, and iii) significant increase in structural flexibility. Experimental facts about the antimicrobial activity of the gomesin and protegrin-1 variants have also been established here, in the hope that the detailed data provided in the present study may be useful for understanding the mechanism of action of these peptides.

Study of the antimicrobial peptide indolicidin and a mutant in micelle medium by molecular dynamics simulation

The antimicrobial peptide indolicidin (IND) and the mutant CP10A in hydrated micelles were studied using molecular dynamics simulations in order to observe whether the molecular dynamics and experimental data could be sufficiently correlated and a detailed description of the interaction of the antimicrobial peptides with a model of the membrane provided by a hydrated micelle system could be obtained. In agreement with the experiments, the simulations showed that the peptides are located near the surface of the micelles. Peptide insertions agree with available experimental data, showing deeper insertion of the mutant compared with the peptide IND. Major insertion into the hydrophobic core of the micelle by all tryptophan and mutated residues of CP10A in relation to IND was observed. The charged residues of the terminus regions of both peptides present similar behavior, indicating that the major differences in the interactions with the micelles of the peptides IND and CP10A occur in the case of the hydrophobic residues.

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.