Modeling the tetraphenylalanine-PEG hybrid amphiphile: from DFT calculations on the peptide to molecular dynamics simulations on the conjugate

The conformational properties of the hybrid amphiphile formed by the conjugation of a hydrophobic peptide with four phenylalanine (Phe) residues and hydrophilic poly(ethylene glycol), have been investigated using quantum mechanical calculations and atomistic molecular dynamics simulations. The intrinsic conformational preferences of the peptide were examined using the building-up search procedure combined with B3LYP/ 6-31G(d) geometry optimizations, which led to the identification of 78, 78, and 92 minimum energy structures for the peptides containing one, two, and four Phe residues. These peptides tend to adopt regular organizations involving turn-like motifs that define ribbon or helicallike arrangements. Furthermore, calculations indicate that backbone ... side chain interactions involving the N-H of the amide groups and the pi clouds of the aromatic rings play a crucial role in Phe-containing peptides. On the other hand,MD simulations on the complete amphiphile in aqueous solution showed that the polymer fragment rapidly unfolds maximizing the contacts with the polar solvent, even though the hydrophobic peptide reduce the number of waters of hydration with respect to an individual polymer chain of equivalent molecular weight. In spite of the small effect of the peptide in the hydrodynamic properties of the polymer, we conclude that the two counterparts of the amphiphile tend to organize as independent modules.

Molecular dynamics simulation insight into the mechanism of phenol adsorption at low coverages from aqueous solutions on microporous carbons

MD simulation studies showing the influence of porosity and carbon surface oxidation on phenol adsorption from aqueous solutions on carbons are reported. Based on a realistic model of activated carbon, three carbon structures with gradually changed microporosity were created. Next, a different number of surface oxygen groups was introduced. The pores with diameters around 0.6 nm are optimal for phenol adsorption and after the introduction of surface oxygen functionalities, adsorption of phenol decreases (in accordance with experimental data) for all studied models. This decrease is caused by a pore blocking effect due to the saturation of surface oxygen groups by highly hydrogen-bounded water molecules.

Photo-induced phase transitions in azobenzene-doped liquid crystals

The effect of irradiation (UV-visible light) on a nematic liquid crystal doped with a photoactive azobenzene derivative was investigated. The selective irradiation results in either an E implies Z or Z implies E isomerization of the azobenzene unit. The effect of the isomerization is to cause a reversible depression of the liquid crystal to isotropic (LC implies l) phase transition temperature of the doped mixture, which can be monitored optically as an isothermal phase transition. This depression also results in a biphasic liquid crystal+isotropic region which is discussed. The authors investigate the cause and magnitude of the phase depression as a function of the amount of doped 4-butyl-4'-methoxyazobenzene (photoactive unit) in 4-cyano-4'-n-pentylbiphenyl (liquid crystal unit), and as a function of the percentage conversion of E implies Z (caused by isomerization) in the azobenzene. The photostationary state of the doped mixtures achieved by Z implies E isomerization is considered and its effect upon the transition temperature of the mixture and response time of the system is discussed. They discuss the implications of the photostationary state with regards to the reversibility of the photo-induced phase transition and hence potential applications.

The effect of the spacer length on the nature of coupling in side chain liquid crystals polymers and elastomers

Side chain liquid crystal polymers and elastomers exhibit a rich phase behaviour which arises from the antagonistic influences of the entropically disordered polymer chain configuration and the long range orientational ordering of the mesogenic units. This competition arises since the natural macroscopic phase separation is inhibited by the inherent chemical connectivity of the system. At the heart of this connectivity is the spacer link and we consider here its influence on the phase behaviour. In particular we consider a series of elastomers in which the number of alkyl units in the spacer is systematically varied from 2 to 6. The lengthening of the coupling spacer is accompanied by an alternation of the sign of coupling between the polymer chain and the mesogenic unit. These results demonstrate the dominating influence of the so-called hinge effect in determining the phase behaviour. In addition to the alternation of the sign there is some decrease in the magnitude of the coupling with increasing spacer length.

Shear banding in molecular dynamics of polymer melts

In order to establish constitutive equations for a viscoelastic fluid uniform shear flow is usually required. However, in the last 10 years S. Q. Wang and co-workers have demonstrated that some entangled polymers do not flow with the uniform shear rate as usually assumed, but instead choose to separate into fast and slow flowing regions. This phenomenon, known as shear banding, causes flow instabilities and in principle invalidates all rheological measurements when it occurs. In this Letter we report the first observation of shear banding in molecular dynamics simulations of entangled polymer melts. We show that our observations are in a very good agreement with the phenomenology developed by Fielding and Olmsted. Our findings provide a simple way of validating the empirical macroscopic phenomenology of shear banding. © 2012 American Physical Society

Grid warehousing of molecular dynamics protein unfolding data

With the increasing awareness of protein folding disorders, the explosion of genomic information, and the need for efficient ways to predict protein structure, protein folding and unfolding has become a central issue in molecular sciences research. Molecular dynamics computer simulations are increasingly employed to understand the folding and unfolding of proteins. Running protein unfolding simulations is computationally expensive and finding ways to enhance performance is a grid issue on its own. However, more and more groups run such simulations and generate a myriad of data, which raises new challenges in managing and analyzing these data. Because the vast range of proteins researchers want to study and simulate, the computational effort needed to generate data, the large data volumes involved, and the different types of analyses scientists need to perform, it is desirable to provide a public repository allowing researchers to pool and share protein unfolding data. This paper describes efforts to provide a grid-enabled data warehouse for protein unfolding data. We outline the challenge and present first results in the design and implementation of the data warehouse.

Mg/Ca partitioning between aqueous solution and aragonite mineral: a molecular dynamics study

We have calculated the concentrations of Mg in the bulk and surfaces of aragonite CaCO3 in equilibrium with aqueous solution, based on molecular dynamics simulations and grand-canonical statistical mechanics. Mg is incorporated in the surfaces, in particular in the (001) terraces, rather than in the bulk of aragonite particles. However, the total Mg content in the bulk and surface of aragonite particles was found to be too small to account for the measured Mg/Ca ratios in corals. We therefore argue that most Mg in corals is either highly metastable in the aragonite lattice, or is located outside the aragonite phase of the coral skeleton, and we discuss the implications of this finding for Mg/Ca paleothermometry.

Scanning tunneling microscopy and molecular dynamics study of the Li2TiO3(001) surface

We have investigated the (001) surface structure of lithium titanate (Li2TiO3) using auger electron spectroscopy (AES), low-energy electron diffraction (LEED), and scanning tunneling microscopy (STM). Li2TiO3 is a potential fusion reactor blanket material. After annealing at 1200 K, LEED demonstrated that the Li2TiO3(001) surface was well ordered and not reconstructed. STM imaging showed that terraces are separated in height by about 0.3 nm suggesting a single termination layer. Moreover, hexagonal patterns with a periodicity of ∼0.4 nm are observed. On the basis of molecular dynamics (MD) simulations, these are interpreted as a dynamic arrangement of Li atoms.

Structure and magnetic properties of polydisperse ferrofluids: a molecular dynamics study

We study by Langevin molecular dynamics simulations systematically the influence of polydispersity in the particle size, and subsequently in the dipole moment, on the physical properties of ferrofluids. The polydispersity is in a first approximation modeled by a bidisperse system that consists of small and large particles at different ratios of their volume fractions. In the first part of our investigations the total volume fraction of the system is fixed, and the volume fraction phi(L) of the large particles is varied. The initial susceptibility chi and magnetization curve of the systems show a strong dependence on the value of phi(L). With the increase of phi(L), the magnetization M of the system has a much faster increment at weak fields, and thus leads to a larger chi. We performed a cluster analysis that indicates that this is due to the aggregation of the large particles in the systems. The average size of these clusters increases with increasing phi(L). In the second part of our investigations, we fixed the volume fraction of the large particles, and increased the volume fraction phi(S) of the small particles in order to study their influence on the chain formation of the large ones. We found that the average aggregate size formed by large particles decreases when phi(S) is increased, demonstrating a significant effect of the small particles on the structural properties of the system. A topological analysis of the structure reveals that the majority of the small particles remain nonaggregated. Only a small number of them are attracted to the ends of the chains formed by large particles.

Molecular dynamics study on the equilibrium magnetization properties and structure of ferrofluids

We investigate in detail the initial susceptibility, magnetization curves, and microstructure of ferrofluids in various concentration and particle dipole moment ranges by means of molecular dynamics simulations. We use the Ewald summation for the long-range dipolar interactions, take explicitly into account the translational and rotational degrees of freedom, coupled to a Langevin thermostat. When the dipolar interaction energy is comparable with the thermal energy, the simulation results on the magnetization properties agree with the theoretical predictions very well. For stronger dipolar couplings, however, we find systematic deviations from the theoretical curves. We analyze in detail the observed microstructure of the fluids under different conditions. The formation of clusters is found to enhance the magnetization at weak fields and thus leads to a larger initial susceptibility. The influence of the particle aggregation is isolated by studying ferro-solids, which consist of magnetic dipoles frozen in at random locations but which are free to rotate. Due to the artificial suppression of clusters in ferrosolids the observed susceptibility is considerably lowered when compared to ferrofluids.

Mode-specific chemisorption of CH4 on Pt{110}-(1 x 2) explored by first-principles molecular dynamics

The chemisorption of CH4 on Pt{110}-(1 x 2) has been studied by vibrational analysis of the reaction pathway defined by the potential energy surface and, in time reversal, by first-principles molecular dynamics simulations of CH4 associative desorption, with the electronic structure treated explicitly using density functional theory. We find that the symmetric stretch vibration ν1 is strongly coupled to the reaction coordinate; our results therefore provide a firm theoretical basis for recently reported state-resolved reactivity measurements, which show that excitation of the ν1 normal mode is the most efficient way to enhance the reaction probability

Bond-selective energy redistribution in the chemisorption of CH3D and CD3H on Pt{1 1 0}-(1 2): a first-principles molecular dynamics study

We have investigated the chemisorption of CH3D and CD3H on Pt{11 0}-(1 2) by performing first-principles molecular dynamics simulations of the recombinative desorption of CH3D (from adsorbed methyl and deuterium) and of CD3H (from adsorbed trideuteromethyl and hydrogen). Vibrational analysis of the symmetry adapted internal coordinates of the desorbing molecules shows that excitation of the single C– D (C–H) bond in the parent molecule is strongly correlated with energy excess in the reaction coordinate. The results of the molecular dynamics simulations are consistent with observed mode- and bond-specific reactivity measurements for chemisorption of methane and its isotopomers on platinum and nickel surfaces.

Polymer brushes under shear: molecular dynamics simulations compared to experiments

Surfaces coated with polymer brushes in a good solvent are known to exhibit excellent tribological properties. We have performed coarse-grained equilibrium and nonequilibrium molecular dynamics (MD) simulations to investigate dextran polymer brushes in an aqueous environment in molecular detail. In a first step, we determined simulation parameters and units by matching experimental results for a single dextran chain. Analyzing this model when applied to a multichain system, density profiles of end-tethered polymer brushes obtained from equilibrium MD simulations compare very well with expectations based on self-consistent field theory. Simulation results were further validated against and correlated with available experimental results. The simulated compression curves (normal force as a function of surface separation) compare successfully with results obtained with a surface forces apparatus. Shear stress (friction) obtained via nonequilibrium MD is contrasted with nanoscale friction studies employing colloidal-probe lateral force microscopy. We find good agreement in the hydrodynamic regime and explain the observed leveling-off of the friction forces in the boundary regime by means of an effective polymer–wall attraction.

A numerical method for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field

A finite difference technique, based on a projection method, is developed for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field. The governing equations in this situation are derived using primitive variables and are solved using the ideas behind the GENSMAC methodology (Tome and McKee [32]; Tome et al. [34]). The resulting numerical technique is then validated by comparing the numerical solution against an analytic solution for steady three-dimensional flow between two-parallel plates subject to a strong magnetic field. The validated code is then employed to solve channel flow for which there is no analytic solution. (C) 2009 Elsevier B.V. All rights reserved.