Local order around rare earth ions during the devitrification of oxyfluoride glasses

Erbium L-3-edge extended x-ray absorption fine structure (EXAFS) measurements were performed on rare earth doped fluorosilicate and fluoroborate glasses and glass ceramics. The well known nucleating effects of erbium ions for the crystallization of cubic lead fluoride (based on x-ray diffraction measurements) and the fact that the rare earth ions are present in the crystalline phase (as indicated by Er3+ emission spectra) seem in contradiction with the present EXAFS analysis, which indicates a lack of medium range structural ordering around the Er3+ ions and suggests that the lead fluoride crystallization does not occur in the nearest neighbor distance of the rare earth ion. Molecular dynamics simulations of the devitrification process of a lead fluoride glass doped with Er3+ ions were performed, and results indicate that Er3+ ions lower the devitrification temperature of PbF2, in good agreement with the experimental results. The genuine role of Er3+ ions in the devitrification process of PbF2 has been investigated. Although Er3+ ions could indeed act as seeds for crystallization, as experiments suggest, molecular dynamics simulation results corroborate the experimental EXAFS observation that the devitrification does not occur at its nearest neighbor distance. (c) 2008 American Institute of Physics.

Análise conformacional dos petídeos protonectina e protonectina (1-6), isolados e em associação, em mistura TFE-água

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Estudo de nanofios de Au e dendrímeros

A primeira parte deste trabalho aborda a simulação computacional de dinâmica molecular clássica da interação de sistemas matriciais constituídos de nanofios paralelos de Au simuladas em função do tempo. Como resultados foram encontrados os tempos de colisões entre os fios da matriz. A segunda parte deste trabalho utiliza dinâmica molecular clássica para simular cinco gerações de dendrímeros PAMAM, cada qual interagindo individualmente com um nanotubo de carbono em função do tempo resultando num motor molecular. Além disso, foram calculados os espectros de absorção deste sistema e foi verificado que eles são nanomotores controlados pela luz. Para todos estes sistemas foram calculadas energias cinética, potencial, total, velocidade, propriedades termodinâmicas como variação de entropia molar, capacidade molar térmica e temperatura in situ. Estas grandezas nos forneceram valiosas informações sobre o comportamento destes sistemas.

Tailoring Thermal Transport Property of Graphene through Oxygen Functionalization

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The 2',4'-dihydroxychalcone could be explored to develop new inhibitors against the glycerol-3-phosphate dehydrogenase from Leishmania species

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Hb S-Sao Paulo: A new sickling hemoglobin with stable polymers and decreased oxygen affinity

Hb S-Sao Paulo (SP) [HBB:c.20A > T p.Glu6Val: c.196A > G p.Lys65Glu] is a new double-mutant hemoglobin that was found in heterozygosis in an 18-month-old Brazilian male with moderate anemia. It behaves like Hb S in acid electrophoresis, isoelectric focusing and solubility testing but shows different behavior in alkaline electrophoresis, cation-exchange HPLC and RP-HPLC. The variant is slightly unstable, showed reduced oxygen affinity and also appeared to form polymers more stable than the Hb S. Molecular dynamics simulation suggests that the polymerization is favored by interfacial electrostatic interactions. This provides a plausible explanation for some of the reported experimental observations. (C) 2012 Elsevier Inc. All rights reserved.

The Equilibrium Structure of Lithium Salt Solutions in Ether-Functionalized Ammonium Ionic Liquids

Molecular dynamics simulations have been performed for ionic liquids based on a ternary mixture of lithium and ammonium cations and a common anion, bis(trifluoromethylsulfonyl)imide, [Tf2N](-). We address structural changes resulting from adding Li+ in ionic liquids with increasing length of an ether-functionalized chain in the ammonium cation. The calculation of static structure factors reveals the lithium effect on charge ordering and intermediate range order in comparison with the neat ionic liquids. The charge ordering is modified in the lithium solution because the coordination of [Tf2N](-) toward Li+ is much stronger than ammonium cations. Intermediate range order is observed in neat ionic liquids based on ammonium cations with a long chain, but in the lithium solutions, there is also a nonhomogenous distribution of Li+ cations. The presence of Li+ enhances interactions between the ammonium cations due to correlations between the oxygen atom of the ether chain and the nitrogen atom of another ammonium cation.

Theoretical studies of microfluidics: "walking" drops and "melting" snowflakes.

The present PhD thesis summarizes two examples of research in microfluidics. Both times water was the subject of interest, once in the liquid state (droplets adsorbed on chemically functionalized surfaces), the other time in the solid state (ice snowflakes and their fractal behaviour). The first problem deals with a slipping nano-droplet of water adsorbed on a surface with photo-switchable wettability characteristics. Main focus was on identifying the underlying driving forces and mechanical principles at the molecular level of detail. Molecular Dynamics simulation was employed as investigative tool owing to its record of successfully describing the microscopic behaviour of liquids at interfaces. To reproduce the specialized surface on which a water droplet can effectively “walk”, a new implicit surface potential was developed. Applying this new method the experimentally observed droplet slippage could be reproduced successfully. Next the movement of the droplet was analyzed at various conditions emphasizing on the behaviour of the water molecules in contact with the surface. The main objective was to identify driving forces and molecular mechanisms underlying the slippage process. The second part of this thesis is concerned with theoretical studies of snowflake melting. In the present work snowflakes are represented by filled von Koch-like fractals of mesoscopic beads. A new algorithm has been developed from scratch to simulate the thermal collapse of fractal structures based on Monte Carlo and Random Walk Simulations (MCRWS). The developed method was applied and compared to Molecular Dynamics simulations regarding the melting of ice snowflake crystals and new parameters were derived from this comparison. Bigger snow-fractals were then studied looking at the time evolution at different temperatures again making use of the developed MCRWS method. This was accompanied by an in-depth analysis of fractal properties (border length and gyration radius) in order to shed light on the dynamics of the melting process.

Ab-initio study of hydrogen technology materials for hydrogen storage and proton conduction

This dissertation deals with two specific aspects of a potential hydrogen-based energy economy, namely the problems of energy storage and energy conversion. In order to contribute to the solution of these problems, the structural and dynamical properties of two promising materials for hydrogen storage (lithium imide/amide) and proton conduction (poly[vinyl phosphonic acid]) are modeled on an atomistic scale by means of first principles molecular dynamics simulation methods.rnrnrnIn the case of the hydrogen storage system lithium amide/imide (LiNH_2/Li_2NH), the focus was on the interplay of structural features and nuclear quantum effects. For these calculations, Path-Integral Molecular Dynamics (PIMD) simulations were used. The structures of these materials at room temperature were elucidated; in collaboration with an experimental group, a very good agreement between calculated and experimental solid-state 1H-NMR chemical shifts was observed. Specifically, the structure of Li_2NH features a disordered arrangement of the Li lattice, which was not reported in previous studies. In addition, a persistent precession of the NH bonds was observed in our simulations. We provide evidence that this precession is the consequence of a toroid-shaped effective potential, in which the protons in the material are immersed. This potential is essentially flat along the torus azimuthal angle, which might lead to important quantum delocalization effects of the protons over the torus.rnrnOn the energy conversion side, the dynamics of protons in a proton conducting polymer (poly[vinyl phosphonic acid], PVPA) was studied by means of a steered ab-initio Molecular Dynamics approach applied on a simplified polymer model. The focus was put on understanding the microscopic proton transport mechanism in polymer membranes, and on characterizing the relevance of the local environment. This covers particularly the effect of water molecules, which participate in the hydrogen bonding network in the material. The results indicate that these water molecules are essential for the effectiveness of proton conduction. A water-mediated Grotthuss mechanism is identified as the main contributor to proton conduction, which agrees with the experimentally observed decay on conductivity for the same material in the absence of water molecules.rnrnThe gain in understanding the microscopic processes and structures present in this materials can help the development of new materials with improved properties, thus contributing to the solution of problems in the implementation of fuel cells.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.

Comparisons of diffusive and viscous contributions to transport coefficients of light gases in single-walled carbon nanotubes

We examine here the relative importance of different contributions to transport of light gases in single walled carbon nanotubes, using methane and hydrogen as examples. Transport coefficients at 298 K are determined using molecular dynamics simulation with atomistic models of the nanotube wall, from which the diffusive and viscous contributions are resolved using a recent approach that provides an explicit expression for the latter. We also exploit an exact theory for the transport of Lennard-Jones fluids at low density considering diffuse reflection at the tube wall, thereby permitting the estimation of Maxwell coefficients for the wall reflection. It is found that reflection from the carbon nanotube wall is nearly specular, as a result of which slip flow dominates, and the viscous contribution is small in comparison, even for a tube as large as 8.1 nm in diameter. The reflection coefficient for hydrogen is 3-6 times as large as that for methane in tubes of 1.36 nm diameter, indicating less specular reflection for hydrogen and greater sensitivity to atomic detail of the surface. This reconciles results showing that transport coefficients for hydrogen and methane, obtained in simulation, are comparable in tubes of this size. With increase in adsorbate density, the reflection coefficient increases, suggesting that adsorbate interactions near the wall serve to roughen the local potential energy landscape perceived by fluid molecules.

Interfacial interactions and structure of polyurethane intercalated nanocomposite

Understanding the interfacial interactions between the nanofiller and polymer matrix is important to improve the design and manufacture of polymer nanocomposites. This paper reports a molecular dynamic Study on the interfacial interactions and structure of a clay-based polyurethane intercalated nanocomposite. The results show that the intercalation of surfactant (i.e. dioctadecyldlmethyl ammonium) and polyurethane (PU) into the nanoconfined gallery of clay leads to the multilayer structure for both surfactant and PU, and the absence of phase separation for PU chains. Such structural characteristics are attributed to the result of competitive interactions among the surfactant, PU and the clay surface, including van der Waals, electrostatic and hydrogen bonding.

Quantum effects on adsorption and diffusion of hydrogen and deuterium in microporous materials

Monte Carlo and molecular dynamics simulations and neutron scattering experiments are used to study the adsorption and diffusion of hydrogen and deuterium in zeolite Rho in the temperature range of 30-150 K. In the molecular simulations, quantum effects are incorporated via the Feynman-Hibbs variational approach. We suggest a new set of potential parameters for hydrogen, which can be used when Feynman-Hibbs variational approach is used for quantum corrections. The dynamic properties obtained from molecular dynamics simulations are in excellent agreement with the experimental results and show significant quantum effects on the transport at very low temperature. The molecular dynamics simulation results show that the quantum effect is very sensitive to pore dimensions and under suitable conditions can lead to a reverse kinetic molecular sieving with deuterium diffusing faster than hydrogen.

An Adaptation of the Hoshen-Kopelman Cluster Counting Algorithm for Honeycomb Networks

We develop a simplified implementation of the Hoshen-Kopelman cluster counting algorithm adapted for honeycomb networks. In our implementation of the algorithm we assume that all nodes in the network are occupied and links between nodes can be intact or broken. The algorithm counts how many clusters there are in the network and determines which nodes belong to each cluster. The network information is stored into two sets of data. The first one is related to the connectivity of the nodes and the second one to the state of links. The algorithm finds all clusters in only one scan across the network and thereafter cluster relabeling operates on a vector whose size is much smaller than the size of the network. Counting the number of clusters of each size, the algorithm determines the cluster size probability distribution from which the mean cluster size parameter can be estimated. Although our implementation of the Hoshen-Kopelman algorithm works only for networks with a honeycomb (hexagonal) structure, it can be easily changed to be applied for networks with arbitrary connectivity between the nodes (triangular, square, etc.). The proposed adaptation of the Hoshen-Kopelman cluster counting algorithm is applied to studying the thermal degradation of a graphene-like honeycomb membrane by means of Molecular Dynamics simulation with a Langevin thermostat. ACM Computing Classification System (1998): F.2.2, I.5.3.

Theoretical study of chloroperoxidase catalyzed chlorination of beta-cyclopentanedione and role of water in the chlorination mechanism

Chloroperoxidase (CPO) is a potential biocatalyst for use in asymmetric synthesis. The mechanisms of CPO catalysis are therefore of interest. The halogenation reaction, one of several chemical reactions that CPO catalyzes, is not fully understood and is the subject of this dissertation. The mechanism by which CPO catalyzes halogenation is disputed. It has been postulated that halogenation of substrates occurs at the active site. Alternatively, it has been proposed that hypochlorous acid, produced at the active site via oxidation of chloride, is released prior to reaction, so that halogenation occurs in solution. The free-solution mechanism is supported by the observation that halogenation of most substrates often occurs non-stereospecifically. On the other hand, the enzyme-bound mechanism is supported by the observation that some large substrates undergo halogenation stereospecifically. The major purpose of this research is to compare chlorination of the substrate β-cyclopentanedione in the two environments. One study was of the reaction with limited hydration because such a level of hydration is typical of the active site. For this work, a purely quantum mechanical approach was used. To model the aqueous environment, the limited hydration environment approach is not appropriate. Instead, reaction precursor conformations were obtained from a solvated molecular dynamics simulation, and reaction of potentially reactive molecular encounters was modeled with a hybrid quantum mechanical/molecular mechanical approach. Extensive work developing parameters for small molecules was pre-requisite for the molecular dynamics simulation. It is observed that a limited and optimized (active-site-like) hydration environment leads to a lower energetic barrier than the fully solvated model representative of the aqueous environment at room temperature, suggesting that the stable water network near the active site is likely to facilitate the chlorination mechanism. The influence of the solvent environment on the reaction barrier is critical. It is observed that stabilization of the catalytic water by other solvent molecules lowers the barrier for keto-enol tautomerization. Placement of water molecules is more important than the number of water molecules in such studies. The fully-solvated model demonstrates that reaction proceeds when the instantaneous dynamical water environment is close to optimal for stabilizing the transition state.