Dynamics of history-dependent epidemics in temporal networks

The structural properties of temporal networks often influence the dynamical processes that occur on these networks, e.g., bursty interaction patterns have been shown to slow down epidemics. In this paper, we investigate the effect of link lifetimes on the spread of history-dependent epidemics. We formulate an analytically tractable activity-driven temporal network model that explicitly incorporates link lifetimes. For Markovian link lifetimes, we use mean-field analysis for computing the epidemic threshold, while the effect of non-Markovian link lifetimes is studied using simulations. Furthermore, we also study the effect of negative correlation between the number of links spawned by an individual and the lifetimes of those links. Such negative correlations may arise due to the finite cognitive capacity of the individuals. Our investigations reveal that heavy-tailed link lifetimes slow down the epidemic, while negative correlations can reduce epidemic prevalence. We believe that our results help shed light on the role of link lifetimes in modulating diffusion processes on temporal networks.

Influence of a Counterion on the Ion Atmosphere of an Anion: A Molecular Dynamics Study of LiX and CsX (X = F-, Cl-, I-) in Methanol

We report molecular dynamics (MD) simulations to explore the influence of a counterion on the structure and dynamics of cationic and anionic solvation shells for various ions in methanol at 298 K. We show that the variation in ionic size of either the cation or the anion in an ion pair influences the solvation structure of the other ion as well as the diffusivity in an electrolyte solution of methanol. The extent of ionic association between the cation and its counteranion of different ionic sizes has been investigated by analyzing the radial distribution functions (RDFs) and the orientation of methanol molecules in the first solvation shell (FSS) of ions. It is shown that the methanol in the FSS of the anion as well the cation exhibit quite different radial and orientational structures as compared to methanol which lie in the FSS of either the anion or the cation but not both. We find that the coordination number (CN) of F-, Cr-, and I- ions decreases with increasing size of the anion which is contrary to the trend reported for the anions in H2O. The mean residence time (MRT) of methanol molecules in the FSS of ions has been calculated using the stable states picture (SSP) approach. It is seen that the ion-counterion interaction has a considerable influence on the MRT of methanol molecules in the FSS of ions. We also discuss the stability order of the ion-counterion using the potentials of mean force (PMFs) for ion pairs with ions of different sizes. The PMF plots reveal that the Li+-F- pair (small-small) is highly stable and the Li+-I- pair is least stable (small-large) in electrolyte solutions.

Scaling law characterizing the dynamics of the transition of HIV-1 to error catastrophe

Increasing the mutation rate, mu, of viruses above a threshold, mu(c), has been predicted to trigger a catastrophic loss of viral genetic information and is being explored as a novel intervention strategy. Here, we examine the dynamics of this transition using stochastic simulations mimicking within-host HIV-1 evolution. We find a scaling law governing the characteristic time of the transition: tau approximate to 0.6/(mu - mu(c)). The law is robust to variations in underlying evolutionary forces and presents guidelines for treatment of HIV-1 infection with mutagens. We estimate that many years of treatment would be required before HIV-1 can suffer an error catastrophe.

Simulations of impinging droplets with surfactant-dependent dynamic contact angle

An arbitrary Lagrangian-Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier-Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surface tension and surfactant-dependent dynamic contact angle. In particular, the dynamic contact angle (theta(d)) of the droplet is defined as a function of the surfactant concentration at the contact line and the equilibrium contact angle (theta(0)(e)) of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included into the model as surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effects of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when theta(0)(e) is less than 90 degrees, and increases it further when theta(0)(e) is greater than 90 degrees. Nevertheless, the presence of surfactants has no effect on the contact angle when theta(0)(e) = 90 degrees. The numerical study clearly demonstrates that the surfactant-dependent contact angle has to be considered, in addition to the Marangoni effect, in order to study the flow dynamics and the equilibrium states of surfactant droplet impingement accurately. The proposed numerical scheme guarantees the conservation of fluid mass and of the surfactant mass accurately. (C) 2015 Elsevier Inc. All rights reserved.

Molecular Dynamics Study of the Structure, Flexibility, and Hydrophilicity of PETIM Dendrimers: A Comparison with PAMAM Dendrimers

A new class of dendrimers, the poly(propyl ether imine) (PETIM) dendrimer, has been shown to be a novel hyperbranched polymer having potential applications as a drug delivery vehicle. Structure and dynamics of the amine terminated PETIM dendrimer and their changes with respect to the dendrimer generation are poorly understood. Since most drugs are hydrophobic in nature, the extent of hydrophobicity of the dendrimer core is related to its drug encapsulation and retention efficacy. In this study, we carry out fully atomistic molecular dynamics (MD) simulations to characterize the structure of PETIM (G2-G6) dendrimers in salt solution as a function of dendrimer generation at different protonation levels. Structural properties such as radius of gyration (R-g), radial density distribution, aspect ratio, and asphericity are calculated. In order to assess the hydrophilicity of the dendrimer, we compute the number of bound water molecules in the interior of dendrirner as well as the number of dendrimer-water hydrogen bonds. We conclude that PETIM dendrimers have relatively greater hydrophobicity and flexibility when compared with their extensively investigated PAMAM counterparts. Hence PETIM dendrimers are expected to have stronger interactions with lipid membranes as well as improved drug encapsulation and retention properties when compared with PAMAM dendrimers. We compute the root-mean-square fluctuation of dendrimers as well as their entropy to quantify the flexibility of the dendrimer. Finally we note that structural and solvation properties computed using force field parameters derived based on the CHARMM general purpose force field were in good quantitative agreement with those obtained using the generalized Amber force field (GAFF).

Nonlinear Dynamics of a Rotating Flexible Link

This paper deals with the study of the nonlinear dynamics of a rotating flexible link modeled as a one dimensional beam, undergoing large deformation and with geometric nonlinearities. The partial differential equation of motion is discretized using a finite element approach to yield four nonlinear, nonautonomous and coupled ordinary differential equations (ODEs). The equations are nondimensionalized using two characteristic velocities-the speed of sound in the material and a velocity associated with the transverse bending vibration of the beam. The method of multiple scales is used to perform a detailed study of the system. A set of four autonomous equations of the first-order are derived considering primary resonances of the external excitation and one-to-one internal resonances between the natural frequencies of the equations. Numerical simulations show that for certain ranges of values of these characteristic velocities, the slow flow equations can exhibit chaotic motions. The numerical simulations and the results are related to a rotating wind turbine blade and the approach can be used for the study of the nonlinear dynamics of a single link flexible manipulator.

Docking, molecular dynamics and QM/MM studies to delineate the mode of binding of CucurbitacinE to F-actin

CucurbitacinE (CurE) has been known to bind covalently to F-actin and inhibit depolymerization. However, the mode of binding of CurE to F-actin and the consequent changes in the F-actin dynamics have not been studied. Through quantum mechanical/molecular mechanical (QM/MM) and density function theory (DFT) simulations after the molecular dynamics (MD) simulations of the docked complex of F-actin and CurE, a detailed transition state (TS) model for the Michael reaction is proposed. The TS model shows nucleophilic attack of the sulphur of Cys257 at the beta-carbon of Michael Acceptor of CurE producing an enol intermediate that forms a covalent bond with CurE. The MD results show a clear difference between the structure of the F-actin in free form and F-actin complexed with CurE. CurE affects the conformation of the nucleotide binding pocket increasing the binding affinity between F-actin and ADP, which in turn could affect the nucleotide exchange. CurE binding also limits the correlated displacement of the relatively flexible domain 1 of F-actin causing the protein to retain a flat structure and to transform into a stable ``tense'' state. This structural transition could inhibit depolymerization of F-actin. In conclusion, CurE allosterically modulates ADP and stabilizes F-actin structure, thereby affecting nucleotide exchange and depolymerization of F-actin. (C) 2015 Elsevier Inc. All rights reserved.

Nanomechanical Mapping, Hierarchical Polymer Dynamics, and Miscibility in the Presence of Chain-End Grafted Nanoparticles

To improve the spatial distribution of nano particles in a polymeric host and to enhance the interfacial interaction with the host, the use of chain-end grafted nanoparticle has gained popularity in the field of polymeric nanocomposites. Besides changing the material properties of the host, these grafted nanoparticles strongly alter the dynamics of the polymer chain at both local and cooperative length scales (relaxations) by manipulating the enthalpic and entropic interactions. It is difficult to map the distribution of these chain-end grafted nanoparticles in the blend by conventional techniques, and herein, we attempted to characterize it by unique technique(s) like peak force quantitative nanomechanical mapping (PFQNM) through AFM (atomic force microscopy) imaging and dielectric relaxation spectroscopy (DRS). Such techniques, besides shedding light on the spatial distribution of the nanoparticles, also give critical information on the changing elasticity at smaller length scales and hierarchical polymer chain dynamics in the vicinity of the nanoparticles. The effect of one-dimensional rodlike multiwall carbon nanotubes (MWNTs), with the characteristic dimension of the order of the radius of gyration of the polymeric chain, on the phase miscibility and chain dynamics in a classical LCST mixture of polystyrene/ poly(vinyl methyl ether) (PS/PVME) was examined in detail using the above techniques. In order to tune the localization of the nanotubes, different molecular weights of PS (13, 31, and 46 kDa), synthesized using RAFT (reversible addition fragmentation chain transfer) polymerization, was grafted onto MWNTs in situ. The thermodynamic miscibility in the blends was assessed by low-amplitude isochronal temperature sweeps, the spatial distribution of MWNTs in the blends was evaluated by PFQNM, and the hierarchical polymer chain dynamics was studied by DRS. It was observed that the miscibility, concentration fluctuation, and cooperative relaxations of the PS/PVME blends are strongly governed by the spatial distribution of MWNTs in the blends. These findings should help guide theories and simulations of hierarchical chain dynamics in LCST mixtures containing rodlike nanoparticles.

Nanomechanical Mapping, Hierarchical Polymer Dynamics, and Miscibility in the Presence of Chain-End Grafted Nanoparticles

To improve the spatial distribution of nano particles in a polymeric host and to enhance the interfacial interaction with the host, the use of chain-end grafted nanoparticle has gained popularity in the field of polymeric nanocomposites. Besides changing the material properties of the host, these grafted nanoparticles strongly alter the dynamics of the polymer chain at both local and cooperative length scales (relaxations) by manipulating the enthalpic and entropic interactions. It is difficult to map the distribution of these chain-end grafted nanoparticles in the blend by conventional techniques, and herein, we attempted to characterize it by unique technique(s) like peak force quantitative nanomechanical mapping (PFQNM) through AFM (atomic force microscopy) imaging and dielectric relaxation spectroscopy (DRS). Such techniques, besides shedding light on the spatial distribution of the nanoparticles, also give critical information on the changing elasticity at smaller length scales and hierarchical polymer chain dynamics in the vicinity of the nanoparticles. The effect of one-dimensional rodlike multiwall carbon nanotubes (MWNTs), with the characteristic dimension of the order of the radius of gyration of the polymeric chain, on the phase miscibility and chain dynamics in a classical LCST mixture of polystyrene/ poly(vinyl methyl ether) (PS/PVME) was examined in detail using the above techniques. In order to tune the localization of the nanotubes, different molecular weights of PS (13, 31, and 46 kDa), synthesized using RAFT (reversible addition fragmentation chain transfer) polymerization, was grafted onto MWNTs in situ. The thermodynamic miscibility in the blends was assessed by low-amplitude isochronal temperature sweeps, the spatial distribution of MWNTs in the blends was evaluated by PFQNM, and the hierarchical polymer chain dynamics was studied by DRS. It was observed that the miscibility, concentration fluctuation, and cooperative relaxations of the PS/PVME blends are strongly governed by the spatial distribution of MWNTs in the blends. These findings should help guide theories and simulations of hierarchical chain dynamics in LCST mixtures containing rodlike nanoparticles.

Molecular Dynamics Simulation of Peeling a DNA Molecule on Substrate

Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.

The Uniaxial Tensile Deformation of Ni Nanowire: Atomic-Scale Computer Simulations

The mechanical deformations of nickel nanowire subjected to uniaxial tensile strain at 300 K are simulated by using molecular dynamics with the quantum corrected Sutten-Chen many-body force field. We have used common neighbor analysis method to investigate the structural evolution of Ni nanowire during the elongation process. For the strain rate of 0.1%/ps, the elastic limit is up to about 11% strain with the yield stress of 8.6 GPa. At the elastic stage, the deformation is carried mainly through the uniform elongation of the distances between the layers (perpendicular to the Z-axis) while the atomic structure remains basically unchanged. With further strain, the slips in the {111} planes start to take place in order to accommodate the applied strain to carry the deformation partially, and subsequently the neck forms. The atomic rearrangements in the neck region result in a zigzag change in the stress-strain curve; the atomic structures beyond the region, however, have no significant changes. With the strain close to the point of the breaking, we observe the formation of a one-atom thick necklace in Ni nanowire. The strain rates have no significant effect on the deformation mechanism, but have some influence on the yield stress, the elastic limit, and the fracture strain of the nanowire.

An analysis of subgrid-resolved scale interactions with use of results from direct numerical simulations

Subgrid nonlinear interaction and energy transfer are analyzed using direct numerical simulations of isotropic turbulence. Influences of cutoff wave number at different ranges of scale on the energetics and dynamics have been investigated. It is observed that subgrid-subgrid interaction dominates the turbulent dynamics when cut-off wave number locates in the energy-containing range while resolved-subgrid interaction dominates if it is in the dissipation range; By decomposing the subgrid energy transfer and nonlinear interaction into 'forward' and 'backward' groups according to the sign of triadic interaction, we find that individually each group has very large contribution, but the net of them is much smaller, implying that tremendous cancellation happens between these two groups.

Molecular dynamics simulation on dislocation emission process

A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.

Forced Dissociation of Selectin-ligand Complexes Using Steered Molecular Dynamics Simulation

Selectin-ligand interactions are crucial to such biological processes as inflammatory cascade or tumor metastasis. How transient formation and dissociation of selectin-ligand bonds in blood flow are coupled to molecular conformation at atomic level, however, has not been well understood. In this study, steered molecular dynamics (SMD) simulations were used to elucidate the intramolecular and intermolecular conformational evolutions involved in forced dissociation of three selectin-ligand systems: the construct consisting of P-selectin lectin (Lec) and epidermal growth factor (EGF)-like domains (P-LE) interacting with synthesized sulfoglycopeptide or SGP-3, P-LE with sialyl Lewis X (sLeX), and E-LE with sLeX. SMD simulations were based on newly built-up force field parameters including carbohydrate units and sulfated tyrosine(s) using an analogy approach. The simulations demonstrated that the complex dissociation was coupled to the molecular extension. While the intramolecular unraveling in P-LESGP-3 system mainly resulted from the destroy of the two anti-parallel sheets of EGF domain and the breakage of hydrogen-bond cluster at the Lec-EGF interface, the intermolecular dissociation was mainly determined by separation of fucose (FUC) from Ca2+ ion in all three systems. Conformational changes during forced dissociations depended on pulling velocities and forces, as well as on how the force was applied. This work provides an insight into better understanding of conformational changes and adhesive functionality of selectin-ligand interactions under external forces.

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated.We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard–Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.