Electronic properties of liquid hydrogen fluoride: A sequential quantum mechanical/Born-Oppenheimer molecular dynamics approach

The electronic properties of liquid hydrogen fluoride (HF) were investigated by carrying out sequential quantum mechanics/Born-Oppenheimer molecular dynamics. The structure of the liquid is in good agreement with recent experimental information. Emphasis was placed on the analysis of polarisation effects, dynamic polarisability and electronic excitations in liquid HF. Our results indicate an increase in liquid phase of the dipole moment (similar to 0.5 D) and isotropic polarisability (5%) relative to their gas-phase values. Our best estimate for the first vertical excitation energy in liquid HF indicates a blue-shift of 0.4 +/- 0.2 eV relative to that of the gas-phase monomer (10.4 eV). (C) 2010 Elsevier B.V. All rights reserved.

Ab initio study of SO2 molecules interacting with pristine and transition metal covered fullerenes as a possible route for nanofilters

In this work, a systematic study of SO2 molecules interacting with pristine and transition metal (TM) covered C-60 is presented by means of first principles calculations. It is observed that the SO2 molecule interacts weakly with the pristine C-60 fullerene, although the resulting interaction is largely increased when the C-60 structure is covered with Fe, Mn, or Ti atoms and the SO2 Molecules are bounded through the TM atoms. The number of bounded SO2 molecules per TM atoms, in addition to the elevated binding energies per molecules, allows us to conclude that such composites can be used as a template for efficient devices to remove SO2 molecules or, alternatively, as SO2 gas sensor.

Hydrogen adsorption on boron doped graphene: an ab initio study

(i) The electronic and structural properties of boron doped graphene sheets, and (ii) the chemisorption processes of hydrogen adatoms on the boron doped graphene sheets have been examined by ab initio total energy calculations. In (i) we find that the structural deformations are very localized around the boron substitutional sites, and in accordance with previous studies (Endo et al 2001 J. Appl. Phys. 90 5670) there is an increase of the electronic density of states near the Fermi level. Our simulated scanning tunneling microscope (STM) images, for occupied states, indicate the formation of bright (triangular) spots lying on the substitutional boron (center) and nearest-neighbor carbon (edge) sites. Those STM images are attributed to the increase of the density of states within an energy interval of 0.5 eV below the Fermi level. For a boron concentration of similar to 2.4%, we find that two boron atoms lying on the opposite sites of the same hexagonal ring (B1-B2 configuration) represents the energetically most stable configuration, which is in contrast with previous theoretical findings. Having determined the energetically most stable configuration for substitutional boron atoms on graphene sheets, we next considered the hydrogen adsorption process as a function of the boron concentration, (ii). Our calculated binding energies indicate that the C-H bonds are strengthened near boron substitutional sites. Indeed, the binding energy of hydrogen adatoms forming a dimer-like structure on the boron doped B1-B2 graphene sheet is higher than the binding energy of an isolated H(2) molecule. Since the formation of the H dimer-like structure may represent the initial stage of the hydrogen clustering process on graphene sheets, we can infer that the formation of H clusters is quite likely not only on clean graphene sheets, which is in consonance with previous studies (Hornekaer et al 2006 Phys. Rev. Lett. 97 186102), but also on B1-B2 boron doped graphene sheets. However, for a low concentration of boron atoms, the formation of H dimer structures is not expected to occur near a single substitutional boron site. That is, the formation (or not) of H clusters on graphene sheets can be tuned by the concentration of substitutional boron atoms.

Tight-binding model for carbon nanotubes from ab initio calculations

Here we present a parametrized tight-binding (TB) model to calculate the band structure of single-wall carbon nanotubes (SWNTs). On the basis of ab initio calculations we fit the band structure of nanotubes of different radii with results obtained with an orthogonal TB model to third neighbors, which includes the effects of orbital hybridization by means of a reduced set of parameters. The functional form for the dependence of these parameters on the radius of the tubes can be used to interpolate appropriate TB parameters for different SWNTs and to study the effects of curvature on their electronic properties. Additionally, we have shown that the model gives an appropriate description of the optical spectra of SWNTs, which can be useful for a proper assignation of SWNTs` specific chirality from optical absorption experiments.

Electric field gradients in (111)In-doped (Hf/Zr)(3)Al(2) and (Hf/Zr)(4)Al(3) mixed compounds: ab initio calculations, perturbed angular correlation measurements and site preference

The quadrupolar hyperfine interactions of in-diffused (111)In -> (111)Cd probes in polycrystalline isostructural Zr(4)Al(3) and Hf(4)Al(3) samples containing small admixtures of the phases (Zr/Hf)(3)Al(2) were investigated. A strong preference of (111)In solutes for the contaminant (Zr/Hf)(3)Al(2) minority phases was observed. Detailed calculations of the electric field gradient (EFG) at the Cd nucleus using the full-potential augmented plane wave + local orbital formalism allowed us to assign the observed EFG fractions to the various lattice sites in the (Zr/Hf)(3)Al(2) compounds and to understand the preferential site occupation of the minority phases by the (111)In atoms. The effects of the size of the supercell and relaxation around the oversized In and Cd probe atoms were investigated in detail.

Ab initio calculation of the BCC Fe-Al-Mo (Iron-Aluminum-Molybdenum) phase diagram: Implications for the nature of the tau(2) phase

The metastable phase diagram of the BCC-based ordering equilibria in the Fe-Al-Mo system has been calculated via a truncated cluster expansion, through the combination of Full-Potential-Linear augmented Plane Wave (FP-LAPW) electronic structure calculations and of Cluster Variation Method (CVM) thermodynamic calculations in the irregular tetrahedron approximation. Four isothermal sections at 1750 K, 2000 K, 2250 K and 2500 K are calculated and correlated with recently published experimental data on the system. The results confirm that the critical temperature for the order-disorder equilibrium between Fe(3)Al-D0(3) and FeAl-B2 is increased by Mo additions, while the critical temperature for the FeAl-B2/A2 equilibrium is kept approximately invariant with increasing Mo contents. The stabilization of the Al-rich A2 phase in equilibrium with overstoichiometric B2-(Fe,Mo)Al is also consistent with the attribution of the A2 structure to the tau(2) phase, stable at high temperatures in overstoichiometric B2-FeAl. (C) 2009 Elsevier Ltd. All rights reserved.

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.

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.

Temperature dependence of molecular dynamics and supramolecular aggregation in MEH-PPV films: A solid-state NMR, X-ray and fluorescence spectroscopy study

This article presents an investigation of the temperature induced modification in the microstructure and dynamics of poly[2-methoxy-5-(2`-ethylhexyloxy)-1,4-phenylenevinylene] (MEH-PPV) cast films using Wide-Angle X-ray Scattering (WAXS), solid-state Nuclear Magnetic Resonance (NMR), and Fluorescence Spectroscopy (PL). MEH-PPV chain motions were characterized as a function of temperature by NMR. The results indicated that the solvent used to cast the films influences the activation energy of the side-chain motions. This was concluded from the comparison of the activation energy of the toluene cast film, E(a) = (54 +/- 8) kJ/mol, and chloroform cast film, E(a) = (69 +/- 5) kJ/mol, and could be attributed to the higher side-chain packing provided by chloroform, that preferentially solvates the side chain in contrast to toluene that solvates mainly the backbone. Concerning the backbone mobility, it was observed that the torsional motions in the MEH-PPV have average amplitude of similar to 10 degrees at 300 K, which was found to be independent of the solvent used to cast the films. In order to correlate the molecular dynamics processes with the changes in the microstructure of the polymer, in situ WAXS experiments as a function of temperature were performed and revealed that the interchain spacing in the MEH-PPV molecular aggregates increases as a function of temperature, particularly at temperatures where molecular relaxations occur. It was also observed that the WAXS peak associated with the bilayer spacing becomes narrower and its intensity increases whereas the peak associated with the inter-backbone planes reduces its intensity for higher temperatures. This last result Could be interpreted as a decrease in the number of aggregates and the reduction of the interchain species during the MEH-PPV relaxation processes. These WAXS results were correlated with PL spectra modifications observed upon temperature treatments. (C) 2008 Elsevier Ltd. All rights reserved.

Combined Crystallographic and Solution Molecular Dynamics Study of Allosteric Effects in Ester and Ketone p-tert-Butylcalix[4]arene Derivatives and Their Complexes with Acetonitrile, Cd(II), and Pb(II)

We describe here a procedure to bridge the gap in the field of calixarene physicochemistry between solid-state atomic-resolution structural information and the liquid-state low-resolution thermodynamics and spectroscopic data. We use MD simulations to study the kinetics and energetics involved in the complexation of lower rim calix[4]arene derivatives (L), containing bidentate ester (1) and ketone (2) pendant groups, with acetonitrile molecule (MeCN) and Cd2+ and Pb2+ ions (M2+) in acetonitrile solution. On one hand, we found that the prior inclusion of MeCN into the calix to form a L(MeCN) adduct has only a weak effect in preorganizing the hydrophilic cavity toward metal ion binding. On the other hand, the strong ion-hydrophilic cavity interaction produces a wide open calix which enhances the binding of one MeCN molecule (allosteric effect) to stabilize the whole (M2+)1(MeCN) bifunctional complex. We reach two major conclusions: (i) the MD results for the (M2+)1(MeCN) binding are in close agreement with the ""endo"", fully encapsulated, metal complex found by X-ray diffraction and in vacuo MD calculations, and (ii) the MD structure for the more flexible 2 ligand, however, differs from the also endo solid-state molecule. In fact, it shows strong solvation effects at the calixarene lower bore by competing MeCN molecules that share the metal coordination sphere with the four C=O oxygens of an ""exo"" (M2+)2(MeCN) complex.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: Evidence from molecular dynamics simulations

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17 beta-estradiol (E-2) and the selective ER modulator raloxifene (RAL)] from the human ER alpha ligand-binding domain in monomeric and dimeric forms. E-2 dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II`). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I`). Remarkably, ER dimerization strongly suppresses Paths II and II` for E-2 dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

The effects of anticalcification treatments and hydration on the molecular dynamics of bovine pericardium collagen as revealed by (13)C solid-state NMR

This article describes a solid-state NMR (SSNMR) investigation of the influence of hydration and chemical cross-linking on the molecular dynamics of the constituents of the bovine pericardium (BP) tissues and its relation to the mechanical properties of the tissue. Samples of natural phenetylamine-diepoxide (DE)- and glutaraldehyde (GL)-fixed BP were investigated by (13)C cross-polarization SSNMR to probe the dynamics of the collagen, and the results were correlated to the mechanical properties of the tissues, probed by dynamical mechanical analysis. For samples of natural BP, the NMR results show that the higher the hydration level the more pronounced the molecular dynamics of the collagen backbone and sidechains, decreasing the tissue`s elastic modulus. In contrast, in DE- and GL-treated samples, the collagen molecules are more rigid, and the hydration seems to be less effective in increasing the collagen molecular dynamics and reducing the mechanical strength of the samples. This is mostly attributed to the presence of cross-links between the collagen plates, which renders the collagen mobility less dependent on the water absorption in chemically treated samples. Copyright (C) 2010 John Wiley & Sons, Ltd.

Molecular Basis of the Thermostability and Thermophilicity of Laminarinases: X-ray Structure of the Hyperthermostable Laminarinase from Rhodothermus marinus and Molecular Dynamics Simulations

Glycosyl hydrolases are enzymes capable of breaking the glycosidic linkage of polysaccharides and have considerable industrial and biotechnological applications. Driven by the later applications, it is frequently desirable that glycosyl hydrolases display stability and activity under extreme environment conditions, such as high temperatures and extreme pHs. Here, we present X-ray structure of the hyperthermophilic laminarinase from Rhodothermus marinus (RmLamR) determined at 1.95 angstrom resolution and molecular dynamics simulation studies aimed to comprehend the molecular basis, for the thermal stability of this class of enzymes. As most thermostable proteins, RmLamR contains a relatively large number of salt bridges, which are not randomly distributed on the structure. On the contrary, they form clusters interconnecting beta-sheets of the catalytic domain. Not all salt bridges, however, are beneficial for the protein thermostability: the existence of charge-charge interactions permeating the hydrophobic core of the enzymes actually contributes to destabilize the structure by facilitating water penetration into hydrophobic cavities, as can be seen in the case of mesophilic enzymes. Furthermore, we demonstrate that the mobility of the side-chains is perturbed differently in each class of enzymes. The side-chains of loop residues surrounding the catalytic cleft in the mesophilic laminarinase gain mobility and obstruct the active site at high temperature. By contrast, thermophilic laminarinases preserve their active site flexibility, and the active-site cleft remains accessible for recognition of polysaccharide substrates even at high temperatures. The present results provide structural insights into the role played by salt-bridges and active site flexibility on protein thermal stability and may be relevant for other classes of proteins, particularly glycosyl hydrolases.

PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009