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.

Relaxation of surface tension in the liquid-solid interfaces of Lennard-Jones liquids

We have established the surface tension relaxation time in the liquid-solid interfaces of Lennard-Jones (LJ) liquids by means of direct measurements in molecular dynamics (MD) simulations. The main result is that the relaxation time is found to be almost independent of the molecular structures and viscosity of the liquids (at seventy-fold change) used in our study and lies in such a range that in slow hydrodynamic motion the interfaces are expected to be at equilibrium. The implications of our results for the modelling of dynamic wetting processes and interpretation of dynamic contact angle data are discussed.

A SAXS study of flow alignment of thermotropic liquid crystal mixtures

We report on the capillary flow behaviour of thermotropic liquid crystal mixtures containing 4-n-octyl-4'-cyanobiphenyl (8CB) and 4-n-pentyl-4'-cyanobiphenyl (5CB). The liquid crystal mixtures are studied in the Nematic (N) and Smectic (SA) phases at room temperature. Polarised optical microscopy (POM), rheology and simultaneous X-ray diffraction (XRD)/capillary flow experiments are performed to characterise the system. Polarised optical microscopy reveals a dramatic change in optical texture when the 5CB content is increased from 20 to 30% in the mixtures. X-ray diffraction results show that the system goes through a SA-N phase transition, such that the mixtures are smectic for 10-20% 5CB and nematic for 30-90% 5CB. Smectic mixtures flow with the layers aligned along the flow direction (mesogens perpendicular to flow) while nematic mixtures flow with the mesogens aligned in the flow direction. Simultaneous XRD/shear flow experiments show that the SA-N transition is independent of the flow rate in the range 1-6 ml min-1. The correlation length of the liquid crystal order decreases with increasing 5CB content. Rheology is used to prove that the correlation length behaviour is related to a reduction in the viscosity of the mixture.

Studying reactive processes with classical dynamics: rebinding dynamics in MbNO

A new surface-crossing algorithm suitable for describing bond-breaking and bond-forming processes in molecular dynamics simulations is presented. The method is formulated for two intersecting potential energy manifolds which dissociate to different adiabatic states. During simulations, crossings are detected by monitoring an energy criterion. If fulfilled, the two manifolds are mixed over a finite number of time steps, after which the system is propagated on the second adiabat and the crossing is carried out with probability one. The algorithm is extensively tested (almost 0.5 mu s of total simulation time) for the rebinding of NO to myoglobin. The unbound surface ((FeNO)-N-...) is represented using a standard force field, whereas the bound surface (Fe-NO) is described by an ab initio potential energy surface. The rebinding is found to be nonexponential in time, in agreement with experimental studies, and can be described using two time constants. Depending on the asymptotic energy separation between the manifolds, the short rebinding timescale is between 1 and 9 ps, whereas the longer timescale is about an order of magnitude larger. NO molecules which do not rebind within 1 ns are typically found in the Xenon-4 pocket, indicating the high affinity of NO to this region in the protein.

Structure-based energetics of protein interfaces guides foot-and-mouth disease virus vaccine design

Virus capsids are primed for disassembly, yet capsid integrity is key to generating a protective immune response. Foot-and-mouth disease virus (FMDV) capsids comprise identical pentameric protein subunits held together by tenuous noncovalent interactions and are often unstable. Chemically inactivated or recombinant empty capsids, which could form the basis of future vaccines, are even less stable than live virus. Here we devised a computational method to assess the relative stability of protein-protein interfaces and used it to design improved candidate vaccines for two poorly stable, but globally important, serotypes of FMDV: O and SAT2. We used a restrained molecular dynamics strategy to rank mutations predicted to strengthen the pentamer interfaces and applied the results to produce stabilized capsids. Structural analyses and stability assays confirmed the predictions, and vaccinated animals generated improved neutralizing-antibody responses to stabilized particles compared to parental viruses and wild-type capsids.

Molecular insights of p47phox phosphorylation dynamics in the regulation of NADPH oxidase activation and superoxide production

Phagocyte superoxide production by a multicomponent NADPH oxidase is important in host defense against microbial invasion. However inappropriate NADPH oxidase activation causes inflammation. Endothelial cells express NADPH oxidase and endothelial oxidative stress due to prolonged NADPH oxidase activation predisposes many diseases. Discovering the mechanism of NADPH oxidase activation is essential for developing novel treatment of these diseases. The p47phox is a key regulatory subunit of NADPH oxidase; however, due to the lack of full protein structural information, the mechanistic insight of p47phox phosphorylation in NADPH oxidase activation remains incomplete. Based on crystal structures of three functional domains, we generated a computational structural model of the full p47phox protein. Using a combination of in silico phosphorylation, molecular dynamics simulation and protein/protein docking, we discovered that the C-terminal tail of p47phox is critical for stabilizing its autoinhibited structure. Ser-379 phosphorylation disrupts H-bonds that link the C-terminal tail to the autoinhibitory region (AIR) and the tandem Src homology 3 (SH3) domains, allowing the AIR to undergo phosphorylation to expose the SH3 pocket for p22phox binding. These findings were confirmed by site-directed mutagenesis and gene transfection of p47phox_/_ coronary microvascular cells. Compared with wild-type p47phoxcDNAtransfected cells, the single mutation of S379A completely blocked p47phox membrane translocation, binding to p22phox and endothelial O2 . production in response to acute stimulation of PKC. p47phox C-terminal tail plays a key role in stabilizing intramolecular interactions at rest. Ser-379 phosphorylation is a molecular switch which initiates p47phox conformational changes and NADPH oxidase-dependent superoxide production by cells.