Molecular dynamics simulation of organic-inorganic nanocomposites: Layering behavior and interlayer structure of organoclays

Isothermal-isobaric (NPT) molecular dynamics simulation has been performed to investigate the layering behavior and structure of nanoconfined quaternary alkylammoniums in organoclays. This work is focused on systems consisting of two clay layers and a number of alkylammoniums, and involves the use of modified Dreiding force field. The simulated basal spacings of organoclays agree satisfactorily with the experimental results in the literature. The atomic density profiles in the direction normal to the clay surface indicate that the alkyl chains within the interlayer space of montmorillonite exhibit an obvious layering behavior. The headgroups of long alkyl chains are distributed within two layers close to the clay surface, whereas the distributions of methyl and methylene groups are strongly dependent on the alkyl chain length and clay layer charge. Monolayer, bilayer, and pseudo-trilayer structures are found in organoclays modified with single long alkyl chains, which are identical to the structural models based on the measured basal spacings. A pseudo-quadrilayer structure, for the first time to our knowledge, is also identified in organoclays with double long alkyl chains. In the mixture structure of paraffin-type and multilayer, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer in pseudo-trilayer as well as next nearest layer in pseudo-quadrilayer.

Experimental tests of quantum nonlinear dynamics in atom optics

Cold atoms in optical potentials provide an ideal test bed to explore quantum nonlinear dynamics. Atoms are prepared in a magneto-optic trap or as a dilute Bose-Einstein condensate and subjected to a far detuned optical standing wave that is modulated. They exhibit a wide range of dynamics, some of which can be explained by classical theory while other aspects show the underlying quantum nature of the system. The atoms have a mixed phase space containing regions of regular motion which appear as distinct peaks in the atomic momentum distribution embedded in a sea of chaos. The action of the atoms is of the order of Planck's constant, making quantum effects significant. This tutorial presents a detailed description of experiments measuring the evolution of atoms in time-dependent optical potentials. Experimental methods are developed providing means for the observation and selective loading of regions of regular motion. The dependence of the atomic dynamics on the system parameters is explored and distinct changes in the atomic momentum distribution are observed which are explained by the applicable quantum and classical theory. The observation of a bifurcation sequence is reported and explained using classical perturbation theory. Experimental methods for the accurate control of the momentum of an ensemble of atoms are developed. They use phase space resonances and chaotic transients providing novel ensemble atomic beamsplitters. The divergence between quantum and classical nonlinear dynamics is manifest in the experimental observation of dynamical tunnelling. It involves no potential barrier. However a constant of motion other than energy still forbids classically this quantum allowed motion. Atoms coherently tunnel back and forth between their initial state of oscillatory motion and the state 180 out of phase with the initial state.

Molecular dynamics simulations:bring biomolecular structures alive on a computer

The molecular dynamics (MD) simulations play a very important role in science today. They have been used successfully in binding free-energy calculations and rational design of drugs and vaccines. MD simulations can help visualize and understand structures and dynamics at an atomistic level when combined with molecular graphics programs. The molecular and atomistic properties can be displayed on a computer in a time-dependent way, which opens a road toward a better understanding of the relationship of structure, dynamics, and function. In this chapter, the basics of MD are explained, together with a step-by-step description of setup and running an MD simulation.

Projected equations of motion approach to hybrid quantum/classical dynamics in dielectric-metal composites

We introduce a hybrid method for dielectric-metal composites that describes the dynamics of the metallic system classically whilst retaining a quantum description of the dielectric. The time-dependent dipole moment of the classical system is mimicked by the introduction of projected equations of motion (PEOM) and the coupling between the two systems is achieved through an effective dipole-dipole interaction. To benchmark this method, we model a test system (semiconducting quantum dot-metal nanoparticle hybrid). We begin by examining the energy absorption rate, showing agreement between the PEOM method and the analytical rotating wave approximation (RWA) solution. We then investigate population inversion and show that the PEOM method provides an accurate model for the interaction under ultrashort pulse excitation where the traditional RWA breaks down.

(Quantum-) critical dynamics of magnetoelectric materials

In this thesis the critical dynamics of several magnetoelectric compounds at their phase transition were examined. Mostly measurements of the dielectric properties in the frequency range of below 1 Hz up to 5 GHz were employed to evaluate the critical exponents for both magnetic field and temperature-dependent measurements. Most of the materials that are part of this work show anomalous behavior, especially at very low temperatures where quantum fluctuations are of the order of or even dominate those induced thermally. This anomalous behavior manifests in different forms. In Dy2Ti2O7 we demonstrate the existence of electric dipoles on magnetic monopoles. Here the dynamics at the critical endpoint located at 0.36K and in a magnetic field of 1T parallel to the [111] direction are of special interest. At this critical endpoint the expected critical slowing down of the dynamics could not only not be observed but instead the opposite, critical speeding-up by several orders of magnitude, could be demonstrated. Furthermore, we show that the phase diagram of Dy2Ti2O7 in this field direction can be reproduced solely from the dynamical properties, for example the resonance frequency of the observed relaxation that is connected to the monopole movement. Away from this point of the phase diagram the dynamics are slowing-down with reduction of temperature as one would expect. Additional measurements on Y2Ti2O7, a structurally identical but non-magnetic material, show only slowing down with reduction of temperature and no additional features. A possible explanation for the observed critical speeding-up is a coherent movement of magnetic monopoles close to the critical field that increases the resonance frequency by reducing the damping of the process. LiCuVO4 on the other hand behaves normally at its phase transition as long as the temperature is higher than 0.4 K. In this temperature regime the dynamics show critical slowing-down analogous to classical ferroelectric materials. This analogy extends also towards higher frequencies where the permittivity displays a ‘dispersion’ minimum that is temperature-dependent but of the order of 2 GHz. Below 0.4K the observed behavior changes drastically. Here we found no longer relaxational behavior but instead an excitation with very low energy. This low energy excitation was predicted by theory and is caused by nearly gapless soliton excitations within the 1D Cu2+ chains of LiCuVO4. Finally, in TbMnO3 the dynamics of the phase transition into the multiferroic phase was observed at roughly 27 K, a much higher temperature compared to the other materials. Here the expected critical slowing-down was observed, even though in low-frequency measurements this transition into the ferroelectric phase is overshadowed by the so-called c-axis relaxation. Therefore, only frequencies above 1MHz could be used to determine the critical exponents for both temperatureand magnetic-field-dependent measurements. This was done for both the peak frequency as well as the relaxation strength. In TbMnO3 an electromagnetic soft-mode with small optical weight causes the observed fluctuations, similar to the case of multiferroic MnWO4.

Distribution of neutral prilocaine in a phospholipid bilayer: Insights from molecular dynamics simulations

In this work, we report a 20-ns constant pressure molecular dynamics simulation of prilocaine (PLC), in amine-amide local anesthetic, in a hydrated liquid crystal bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine. The partition of PLC induces the lateral expansion of the bilayer and a concomitant contraction in its thickness. PLC molecules are preferentially found in the hydrophobic acyl chains region, with a maximum probability at similar to 12 angstrom from the center of the bilayer (between the C(4) and C(5) methylene groups). A decrease in the acyl chain segmental order parameter, vertical bar S-CD vertical bar, compared to neat bilayers, is found, in good agreement with experimental H-2-NMR studies. The decrease in vertical bar S-CD vertical bar induced by PLC is attributed to a larger accessible volume per lipid in the acyl chain region. (C) 2008 Wiley Periodicals, Inc.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

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.

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.

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.

Fragility and glassy dynamics of 2Ca(NO(3))(2)center dot 3KNO(3) under pressure: Molecular dynamics simulations

Molecular dynamics simulations of the glass-forming liquid 2Ca(NO(3))(2)center dot 3KNO(3) (CKN) were performed from high temperature liquid states down to low temperature glassy states at six different pressures from 10(-4) to 5.0 GPa. The temperature dependence of the structural relaxation time indicates that the fragility of liquid CKN changes with pressure. In line with recent proposal [Scopigno , Science 302, 849 (2003)], the change on liquid fragility is followed by a proportional change of the nonergodicity factor of the corresponding glass at low temperature. (c) 2008 American Institute of Physics.

Molecular dynamics, flexible docking, virtual screening, ADMET predictions, and molecular interaction field studies to design novel potential MAO-B inhibitors

Monoamine oxidase is a flavoenzyme bound to the mitochondrial outer membranes of the cells, which is responsible for the oxidative deamination of neurotransmitter and dietary amines. It has two distinct isozymic forms, designated MAO-A and MAO-B, each displaying different substrate and inhibitor specificities. They are the well-known targets for antidepressant, Parkinson`s disease, and neuroprotective drugs. Elucidation of the x-ray crystallographic structure of MAO-B has opened the way for the molecular modeling studies. In this work we have used molecular modeling, density functional theory with correlation, virtual screening, flexible docking, molecular dynamics, ADMET predictions, and molecular interaction field studies in order to design new molecules with potential higher selectivity and enzymatic inhibitory activity over MAO-B.