Molecular dynamics insight to phase transition in n-alkanes with carbon nanofillers

The present work aims to investigate the phase transition, dispersion and diffusion behavior of nanocomposites of carbon nanotube (CNT) and straight chain alkanes. These materials are potential candidates for organic phase change materials(PCMs) and have attracted flurry of research recently. Accurate experimental evaluation of the mass, thermal and transport properties of such composites is both difficult as well as economically taxing. Additionally it is crucial to understand the factors that results in modification or enhancement of their characteristic at atomic or molecular level. Classical molecular dynamics approach has been extended to elucidate the same. Bulk atomistic models have been generated and subjected to rigorous multistage equilibration. To reaffirm the approach, both canonical and constant-temperature, constant-pressure ensembles were employed to simulate the models under consideration. Explicit determination of kinetic, potential, non-bond and total energy assisted in understanding the enhanced thermal and transport property of the nanocomposites from molecular point of view. Crucial parameters including mean square displacement and simulated self diffusion coefficient precisely define the balance of the thermodynamic and hydrodynamic interactions. Radial distribution function also reflected the density variation, strength and mobility of the nanocomposites. It is expected that CNT functionalization could improve the dispersion within n-alkane matrix. This would further ameliorate the mass and thermal properties of the composite. Additionally, the determined density was in good agreement with experimental data. Thus, molecular dynamics can be utilized as a high throughput technique for theoretical investigation of nanocomposites PCMs. (C) 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License.

Dynamics of Plasmon in Graphene Oxide

Dynamic effects of plasmon such as scattering with defect boundaries and oxygen impurities in the graphene oxide are investigated. Study of plasmon dynamics helps in understanding electronic, opto-electronic and biological applications of graphene based nanostructures. Tuning or control over such applications is made possible by graphene nanostructure engineering. We have modeled defects with increased smoothing of defect edge in graphene keeping area of the defect constant. Scattering of plasmons in graphene with defects is modeled using an electromagnetic field coupled inter-atomic potential approach with finite element discretization of the atomic vibrational and electromagnetic field degrees of freedom. Our calculations show pi + sigma plasmon red shifting under sharp defect edges whereas pi plasmon show high extinction efficiency. Strong localization of electric fields near the sharp defect edges is observed. Observations on plasmons and its dynamics draws attention in designing novel optoelectronic devices and binders for bio-molecules.

Interaction of a vortex ring with a single bubble: bubble and vorticity dynamics

The interaction of a single bubble with a single vortex ring in water has been studied experimentally. Measurements of both the bubble dynamics and vorticity dynamics have been done to help understand the two-way coupled problem. The circulation strength of the vortex ring (Gamma) has been systematically varied, while keeping the bubble diameter (D-b) constant, with the bubble volume to vortex core volume ratio (V-R) also kept fixed at about 0.1. The other important parameter in the problem is a Weber number based on the vortex ring strength. (We = 0.87 rho(Gamma/2 pi a)(2)/(sigma/D-b); a = vortex core radius, sigma = surface tension), which is varied over a large range, We = 3-406. The interaction between the bubble and ring for each of the We cases broadly falls into four stages. Stage I is before capture of the bubble by the ring where the bubble is drawn into the low-pressure vortex core, while in stage II the bubble is stretched in the azimuthal direction within the ring and gradually broken up into a number of smaller bubbles. Following this, in stage III the bubble break-up is complete and the resulting smaller bubbles slowly move around the core, and finally in stage IV the bubbles escape. Apart from the effect of the ring on the bubble, the bubble is also shown to significantly affect the vortex ring, especially at low We (We similar to 3). In these low-We cases, the convection speed drops significantly compared to the base case without a bubble, while the core appears to fragment with a resultant large decrease in enstrophy by about 50 %. In the higher-We cases (We > 100), there are some differences in convection speed and enstrophy, but the effects are relatively small. The most dramatic effects of the bubble on the ring are found for thicker core rings at low We (We similar to 3) with the vortex ring almost stopping after interacting with the bubble, and the core fragmenting into two parts. The present idealized experiments exhibit many phenomena also seen in bubbly turbulent flows such as reduction in enstrophy, suppression of structures, enhancement of energy at small scales and reduction in energy at large scales. These similarities suggest that results from the present experiments can be helpful in better understanding interactions of bubbles with eddies in turbulent flows.

Dynamics of Cricket Sound Production

The clever designs of natural transducers are a great source of inspiration for man-made systems. At small length scales, there are many transducers in nature that we are now beginning to understand and learn from. Here, we present an example of such a transducer that is used by field crickets to produce their characteristic song. This transducer uses two distinct components-a file of discrete teeth and a plectrum that engages intermittently to produce a series of impulses forming the loading, and an approximately triangular membrane, called the harp, that acts as a resonator and vibrates in response to the impulse-train loading. The file-and-plectrum act as a frequency multiplier taking the low wing beat frequency as the input and converting it into an impulse-train of sufficiently high frequency close to the resonant frequency of the harp. The forced vibration response results in beats producing the characteristic sound of the cricket song. With careful measurements of the harp geometry and experimental measurements of its mechanical properties (Young's modulus determined from nanoindentation tests), we construct a finite element (FE) model of the harp and carry out modal analysis to determine its natural frequency. We fine tune the model with appropriate elastic boundary conditions to match the natural frequency of the harp of a particular species-Gryllus bimaculatus. We model impulsive loading based on a loading scheme reported in literature and predict the transient response of the harp. We show that the harp indeed produces beats and its frequency content matches closely that of the recorded song. Subsequently, we use our FE model to show that the natural design is quite robust to perturbations in the file. The characteristic song frequency produced is unaffected by variations in the spacing of file-teeth and even by larger gaps. Based on the understanding of how this natural transducer works, one can design and fabricate efficient microscale acoustic devices such as microelectromechanical systems (MEMS) loudspeakers.

Phenomenology and control of buckling dynamics in multicomponent colloidal droplets

Self-assembly of nano sized particles during natural drying causes agglomeration and shell formation at the surface of micron sized droplets. The shell undergoes sol-gel transition leading to buckling at the weakest point on the surface and produces different types of structures. Manipulation of the buckling rate with inclusion of surfactant (sodium dodecyl sulphate, SDS) and salt (anilinium hydrochloride, AHC) to the nano-sized particle dispersion (nanosilica) is reported here in an acoustically levitated single droplet. Buckling in levitated droplets is a cumulative, complicated function of acoustic streaming, chemistry, agglomeration rate, porosity, radius of curvature, and elastic energy of shell. We put forward our hypothesis on how buckling occurs and can be suppressed during natural drying of the droplets. Global precipitation of aggregates due to slow drying of surfactant-added droplets (no added salts) enhances the rigidity of the shell formed and hence reduces the buckling probability of the shell. On the contrary, adsorption of SDS aggregates on salt ions facilitates the buckling phenomenon with an addition of minute concentration of the aniline salt to the dispersion. Variation in the concentration of the added particles (SDS/AHC) also leads to starkly different morphologies and transient behaviour of buckling (buckling modes like paraboloid, ellipsoid, and buckling rates). Tuning of the buckling rate causes a transition in the final morphology from ring and bowl shapes to cocoon type of structure. (C) 2015 AIP Publishing LLC.

Catalytic pathway, substrate binding and stability in SAICAR synthetase: A structure and molecular dynamics study

The de novo purine biosynthesis is one of the highly conserved pathways among all organisms and is essential for the cell viability. A clear understanding of the enzymes in this pathway would pave way for the development of antimicrobial and anticancer drugs. Phosphoribosylaminoimidazole-succinocar boxamide (SAICAR) synthetase is one of the enzymes in this pathway that catalyzes ATP dependent ligation of carboxyaminoimidazole ribotide (CAIR) with L-aspartate (ASP). Here, we describe eight crystal structures of this enzyme, in C222(1) and H3 space groups, bound to various substrates and substrate mimics from a hyperthermophilic archaea Pyrococcus horikoshii along with molecular dynamics simulations of the structures with substrates. Complexes exhibit minimal deviation from its apo structure. The CAIR binding site displays a preference for pyrimidine nucleotides. In the ADP.TMP-ASP complex, the ASP binds at a position equivalent to that found in Saccharomyces cerevisiae structure (PDB: 2CNU) and thus, clears the ambiguity regarding ASP's position. A possible mode for the inhibition of the enzyme by CTP and UTP, observed earlier in the yeast enzyme, is clearly illustrated in the structures bound to CMP and UMP. The ADP.Mg2+.PO4.CD/MP complex having a phosphate ion between the ATP and CAIR sites strengthens one of the two probable pathways (proposed in Escherichia coli study) of catalytic mechanism and suggests the possibility of a phosphorylation taking place before the ASP's attack on CAIR. Molecular dynamic simulations of this enzyme along with its substrates at 90 degrees C reveal the relative strengths of substrate binding, possible antagonism and the role of Mg2+ ions. (C) 2015 Elsevier Inc. All rights reserved.

Dynamics of Glass Forming Liquids with Randomly Pinned Particles

It is frequently assumed that in the limit of vanishing cooling rate, the glass transition phenomenon becomes a thermodynamic transition at a temperature T-K. However, with any finite cooling rate, the system falls out of equilibrium at temperatures near T-g(> T-K), implying that the very existence of the putative thermodynamic phase transition at T-K can be questioned. Recent studies of systems with randomly pinned particles have hinted that the thermodynamic glass transition may be observed for liquids with randomly pinned particles. This expectation is based on the results of approximate calculations that suggest that the thermodynamic glass transition temperature increases with increasing concentration of pinned particles and it may be possible to equilibrate the system at temperatures near the increased transition temperature. We test the validity of this prediction through extensive molecular dynamics simulations of two model glass-forming liquids in the presence of random pinning. We find that extrapolated thermodynamic transition temperature T-K does not show any sign of increasing with increasing pinning concentration. The main effect of pinning is found to be a rapid decrease in the kinetic fragility of the system with increasing pin concentration. Implications of these observations for current theories of the glass transition are discussed.

Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man(9)GlcNAc(2)). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct beta sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

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.

Time scales of the stick-slip dynamics of the peeling of an adhesive tape

The stick-slip dynamics of the peeling of an adhesive tape is characterized by bifurcations that have been experimentally well studied. In this work, we investigate the time scale in which the the stick-slips happen leading to the bifurcations. This is fundamental to understanding the triboluminescence and acoustic emissions associated with the bifurcations. We establish a relationship between the time scale of the bifurcations and the inherent mathematical structure of the peeling dynamics by studying a characteristic time quantity associated with the dynamics.

EFFECT OF PLUME DYNAMICS FOR HEAT TRANSFER ENHANCEMENT AT SOLID-LIQUID INTERFACE

The present paper analyzes the effects of plumes for heat transfer enhancement at solid-liquid interface taking both smooth and grooved surfaces. The experimental setup consists of a tank of dimensions 265 x 265 x 300 (height) containing water. The bottom surface was heated and free surface of the water was left open to the ambient. In the experiments, the bottom plate had either a smooth surface or a grooved surface. We used 90 V-grooved rough surfaces with two groove heights, 10mm and 3mm. The experiment was done with water layer depths of 90mm and 140mm, corresponding to values of aspect ratio(AR) equal to 2.9 and 1.8 respectively. Thymol blue, a pH sensitive dye, was used to visualize the flow near the heated plate. The measured heat transfer coefficients over the grooved surfaces were higher compared that over the smooth surface. The enhanced heat transport in the rough cavities cannot be ascribed to the increase in the contact area, rather it must be the local dynamics of the thermal boundary layer that changes the heat transport over the rough surface.

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.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

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.