Application of a new reverse Monte Carlo algorithm to polyatomic molecular systems. I. Liquid water

Using a new reverse Monte Carlo algorithm, we present simulations that reproduce very well several structural and thermodynamic properties of liquid water. Both Monte Carlo, molecular dynamics simulations and experimental radial distribution functions used as input are accurately reproduced using a small number of molecules and no external constraints. Ad hoc energy and hydrogen bond analysis show the physical consistency and limitations of the generated RMC configurations. (C) 2001 American Institute of Physics.

Electromagnetic form factor of the pion in the space- and time-like regions within the front-form dynamics

The pion electromagnetic form factor is calculated in the space- and time-like regions from -10 (GeV/c)2 up to 10 (GeV/c)2, within a front-form model. The dressed photon vertex where a photon decays in a quark-antiquark pair is depicted generalizing the vector meson dominance ansatz, by means of the vector meson vertex functions. An important feature of our model is the description of the on-mass-shell vertex functions in the valence sector, for the pion and the vector mesons, through the front-form wave functions obtained within a realistic quark model. The theoretical results show an excellent agreement with the data in the space-like region, while in the time-like region the description is quite encouraging. © 2003 Elsevier B.V. All rights reserved.

Estudos conformacionais por dinâmica molecular de peptídeos antimicrobianos da família dos Mastoparanos em misturas de TFE-água

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

Aspergillus niger beta-Glucosidase Has a Cellulase-like Tadpole Molecular Shape INSIGHTS INTO GLYCOSIDE HYDROLASE FAMILY 3 (GH3) beta-GLUCOSIDASE STRUCTURE AND FUNCTION

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

Combining Experimental Evidence and Molecular Dynamic Simulations To Understand the Mechanism of Action of the Antimicrobial Octapeptide Jelleine-I

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

Mechanical properties and fracture dynamics of silicene membranes

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

Prospecção de novos peptídeos com ação leishmanicida por dinâmica molecular

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

Conformational changes in a hyperthermostable glycoside hydrolase: enzymatic activity is a consequence of the loop dynamics and protonation balance

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

Estudos da interação de peptídeos antimicrobianos com modelo de membrana por simulações de dinâmica molecular

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

Computer simulation study of collective dynamics in the glass former Ca(NO3)(2)center dot 4H(2)O

Time correlation functions of current fluctuations were calculated by molecular dynamics (MD) simulations in order to investigate sound waves of high wavevectors in the glass-forming liquid Ca(NO3)(2)center dot 4H(2)O. Dispersion curves, omega(k), were obtained for longitudinal (LA) and transverse acoustic (TA) modes, and also for longitudinal optic (LO) modes. Spectra of LA modes calculated by MD simulations were modeled by a viscoelastic model within the memory function framework. The viscoelastic model is used to rationalize the change of slope taking place at k similar to 0.3 angstrom(-1) in the omega(k) curve of acoustic modes. For still larger wavevectors, mixing of acoustic and optic modes is observed. Partial time correlation functions of longitudinal mass currents were calculated separately for the ions and the water molecules. The wavevector dependence of excitation energies of the corresponding partial LA modes indicates the coexistence of a relatively stiff subsystem made of cations and anions, and a softer subsystem made of water molecules. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751548]

Cytochrome c-promoted cardiolipin oxidation generates singlet molecular oxygen

The interaction of cytochrome c (cyt c) with cardiolipin (CL) induces protein conformational changes that favor peroxidase activity. This process has been correlated with CL oxidation and the induction of cell death. Here we report evidence demonstrating the generation of singlet molecular oxygen [O-2((1)Delta(g))] by a cyt c-CL complex in a model membrane containing CL. The formation of singlet oxygen was directly evidenced by luminescence measurements at 1270 nm and by chemical trapping experiments. Singlet oxygen generation required cyt c-CL binding and occurred at pH values higher than 6, consistent with lipid-protein interactions involving fully deprotonated CL species and positively charged residues in the protein. Moreover, singlet oxygen formation was specifically observed for tetralinoleoyl CL species and was not observed with monounsaturated and saturated CL species. Our results show that there are at least two mechanisms leading to singlet oxygen formation: one with fast kinetics involving the generation of singlet oxygen directly from CL hydroperoxide decomposition and the other involving CL oxidation. The contribution of the first mechanism was clearly evidenced by the detection of labeled singlet oxygen [O-18(2)((1)Delta(g))] from liposomes supplemented with 18-oxygen-labeled CL hydroperoxides. However quantitative analysis showed that singlet oxygen yield from CL hydroperoxides was minor (<5%) and that most of the singlet oxygen is formed from the second mechanism. Based on these data and previous findings we propose a mechanism of singlet oxygen generation through reactions involving peroxyl radicals (Russell mechanism) and excited triplet carbonyl intermediates (energy transfer mechanism).

Selectivity and Mechanisms Driven by Reaction Dynamics: The Case of the Gas-Phase OH-+CH3ONO2 Reaction

Well-established statistical approaches such as transition-state theory based on high-level calculated potential energy profiles are unable to account for the selectivity observed in the gas-phase OH- + CH3ONO2 reaction. This reaction can undergo bimolecular nucleophilic displacement at either the carbon center (S(N)2@C) or the nitrogen center (S(N)2@N) as well as a proton abstraction followed by dissociation (E(CO)2) pathway. Direct dynamics simulations yield an S(N)2:E(CO)2 product ratio in close agreement with experiment and show that the lack of reactivity at the nitrogen atom is due to the highly negative electrostatic potential generated by the oxygen atoms in the ONO2 group that scatters the incoming OH-. In addition to these dynamical effects, the nonstatistical behavior of these reactions is attributed to the absence of equilibrated reactant complexes and to the large number of recrossings, which might be present in several ion-molecule gas-phase reactions.

Coarse-graining methods and polymer dynamics

This thesis studies molecular dynamics simulations on two levels of resolution: the detailed level of atomistic simulations, where the motion of explicit atoms in a many-particle system is considered, and the coarse-grained level, where the motion of superatoms composed of up to 10 atoms is modeled. While atomistic models are capable of describing material specific effects on small scales, the time and length scales they can cover are limited due to their computational costs. Polymer systems are typically characterized by effects on a broad range of length and time scales. Therefore it is often impossible to atomistically simulate processes, which determine macroscopic properties in polymer systems. Coarse-grained (CG) simulations extend the range of accessible time and length scales by three to four orders of magnitude. However, no standardized coarse-graining procedure has been established yet. Following the ideas of structure-based coarse-graining, a coarse-grained model for polystyrene is presented. Structure-based methods parameterize CG models to reproduce static properties of atomistic melts such as radial distribution functions between superatoms or other probability distributions for coarse-grained degrees of freedom. Two enhancements of the coarse-graining methodology are suggested. Correlations between local degrees of freedom are implicitly taken into account by additional potentials acting between neighboring superatoms in the polymer chain. This improves the reproduction of local chain conformations and allows the study of different tacticities of polystyrene. It also gives better control of the chain stiffness, which agrees perfectly with the atomistic model, and leads to a reproduction of experimental results for overall chain dimensions, such as the characteristic ratio, for all different tacticities. The second new aspect is the computationally cheap development of nonbonded CG potentials based on the sampling of pairs of oligomers in vacuum. Static properties of polymer melts are obtained as predictions of the CG model in contrast to other structure-based CG models, which are iteratively refined to reproduce reference melt structures. The dynamics of simulations at the two levels of resolution are compared. The time scales of dynamical processes in atomistic and coarse-grained simulations can be connected by a time scaling factor, which depends on several specific system properties as molecular weight, density, temperature, and other components in mixtures. In this thesis the influence of molecular weight in systems of oligomers and the situation in two-component mixtures is studied. For a system of small additives in a melt of long polymer chains the temperature dependence of the additive diffusion is predicted and compared to experiments.

Mode-coupling theory: generalizations, high dimensions and microscopic dynamics

This work contains several applications of the mode-coupling theory (MCT) and is separated into three parts. In the first part we investigate the liquid-glass transition of hard spheres for dimensions d→∞ analytically and numerically up to d=800 in the framework of MCT. We find that the critical packing fraction ϕc(d) scales as d²2^(-d), which is larger than the Kauzmann packing fraction ϕK(d) found by a small-cage expansion by Parisi and Zamponi [J. Stat. Mech.: Theory Exp. 2006, P03017 (2006)]. The scaling of the critical packing fraction is different from the relation ϕc(d)∼d2^(-d) found earlier by Kirkpatrick and Wolynes [Phys. Rev. A 35, 3072 (1987)]. This is due to the fact that the k dependence of the critical collective and self nonergodicity parameters fc(k;d) and fcs(k;d) was assumed to be Gaussian in the previous theories. We show that in MCT this is not the case. Instead fc(k;d) and fcs(k;d), which become identical in the limit d→∞, converge to a non-Gaussian master function on the scale k∼d^(3/2). We find that the numerically determined value for the exponent parameter λ and therefore also the critical exponents a and b depend on the dimension d, even at the largest evaluated dimension d=800. In the second part we compare the results of a molecular-dynamics simulation of liquid Lennard-Jones argon far away from the glass transition [D. Levesque, L. Verlet, and J. Kurkijärvi, Phys. Rev. A 7, 1690 (1973)] with MCT. We show that the agreement between theory and computer simulation can be improved by taking binary collisions into account [L. Sjögren, Phys. Rev. A 22, 2866 (1980)]. We find that an empiric prefactor of the memory function of the original MCT equations leads to similar results. In the third part we derive the equations for a mode-coupling theory for the spherical components of the stress tensor. Unfortunately it turns out that they are too complex to be solved numerically.