Exploring behavioural change through an agent-oriented system dynamics model: the use of personal protective equipment among pesticide applicators in Colombia

The misuse of Personal Protective Equipment results in health risk among smallholders in developing countries, and education is often proposed to promote safer practices. However, evidence point to limited effects of education. This paper presents a System Dynamics model which allows the identification of risk-minimizing policies for behavioural change. The model is based on the IAC framework and survey data. It represents farmers' decision-making from an agent-oriented standpoint. The most successful intervention strategy was the one which intervened in the long term, targeted key stocks in the systems and was diversified. However, the results suggest that, under these conditions, no policy is able to trigger a self sustaining behavioural change. Two implementation approaches were suggested by experts. One, based on constant social control, corresponds to a change of the current model's parameters. The other, based on participation, would lead farmers to new thinking, i.e. changes in their decision-making structure.

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.

Modeling entangled dynamics: comparison between stochastic single-chain and multichain models

To test the effectiveness of stochastic single-chain models in describing the dynamics of entangled polymers, we systematically compare one such model; the slip-spring model; to a multichain model solved using stochastic molecular dynamics(MD) simulations (the Kremer-Grest model). The comparison involves investigating if the single-chain model can adequately describe both a microscopic dynamical and a macroscopic rheological quantity for a range of chain lengths. Choosing a particular chain length in the slip-spring model, the parameter values that best reproduce the mean-square displacement of a group of monomers is determined by fitting toMDdata. Using the same set of parameters we then test if the predictions of the mean-square displacements for other chain lengths agree with the MD calculations. We followed this by a comparison of the time dependent stress relaxation moduli obtained from the two models for a range of chain lengths. After identifying a limitation of the original slip-spring model in describing the static structure of the polymer chain as seen in MD, we remedy this by introducing a pairwise repulsive potential between the monomers in the chains. Poor agreement of the mean-square monomer displacements at short times can be rectified by the use of generalized Langevin equations for the dynamics and resulted in significantly improved agreement.

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.

Time-dependent orientation coupling in equilibrium polymer melts

The motion in concentrated polymer systems is described by either the Rouse or the reptation model, which both assume that the relaxation of each polymer chain is independent of the surrounding chains. This, however, is in contradiction with several experiments. In this Letter, we propose a universal description of orientation coupling in polymer melts in terms of the time-dependent coupling parameter κ(t). We use molecular dynamics simulations to show that the coupling parameter increases with time, reaching about 50% at long times, independently of the chain length or blend composition. This leads to predictions of component dynamics in mixtures of different molecular weights from the knowledge of monodisperse dynamics for unentangled melts. Finally, we demonstrate that entanglements do not play a significant role in the observed coupling. © 2010 The American Physical Society

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Segmental dynamics in entangled linear polymer melts

We present molecular dynamics (MD) and slip-springs model simulations of the chain segmental dynamics in entangled linear polymer melts. The time-dependent behavior of the segmental orientation autocorrelation functions and mean-square segmental displacements are analyzed for both flexible and semiflexible chains, with particular attention paid to the scaling relations among these dynamic quantities. Effective combination of the two simulation methods at different coarse-graining levels allows us to explore the chain dynamics for chain lengths ranging from Z ≈ 2 to 90 entanglements. For a given chain length of Z ≈ 15, the time scales accessed span for more than 10 decades, covering all of the interesting relaxation regimes. The obtained time dependence of the monomer mean square displacements, g1(t), is in good agreement with the tube theory predictions. Results on the first- and second-order segmental orientation autocorrelation functions, C1(t) and C2(t), demonstrate a clear power law relationship of C2(t) C1(t)m with m = 3, 2, and 1 in the initial, free Rouse, and entangled (constrained Rouse) regimes, respectively. The return-to-origin hypothesis, which leads to inverse proportionality between the segmental orientation autocorrelation functions and g1(t) in the entangled regime, is convincingly verified by the simulation result of C1(t) g1(t)−1 t–1/4 in the constrained Rouse regime, where for well-entangled chains both C1(t) and g1(t) are rather insensitive to the constraint release effects. However, the second-order correlation function, C2(t), shows much stronger sensitivity to the constraint release effects and experiences a protracted crossover from the free Rouse to entangled regime. This crossover region extends for at least one decade in time longer than that of C1(t). The predicted time scaling behavior of C2(t) t–1/4 is observed in slip-springs simulations only at chain length of 90 entanglements, whereas shorter chains show higher scaling exponents. The reported simulation work can be applied to understand the observations of the NMR experiments.

Significance of cross correlations in the stress relaxation of polymer melts

According to linear response theory, all relaxation functions in the linear regime can be obtained using time correlation functions calculated under equilibrium. In this paper, we demonstrate that the cross correlations make a significant contribution to the partial stress relaxation functions in polymer melts. We present two illustrations in the context of polymer rheology using (1) Brownian dynamics simulations of a single chain model for entangled polymers, the slip-spring model, and (2) molecular dynamics simulations of a multichain model. Using the single chain model, we analyze the contribution of the confining potential to the stress relaxation and the plateau modulus. Although the idea is illustrated with a particular model, it applies to any single chain model that uses a potential to confine the motion of the chains. This leads us to question some of the assumptions behind the tube theory, especially the meaning of the entanglement molecular weight obtained from the plateau modulus. To shed some light on this issue, we study the contribution of the nonbonded excluded-volume interactions to the stress relaxation using the multichain model. The proportionality of the bonded/nonbonded contributions to the total stress relaxation (after a density dependent "colloidal" relaxation time) provides some insight into the success of the tube theory in spite of using questionable assumptions. The proportionality indicates that the shape of the relaxation spectrum can indeed be reproduced using the tube theory and the problem is reduced to that of finding the correct prefactor. (c) 2007 American Institute of Physics