Source author record

Biman Bagchi

Biman Bagchi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
13topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2022arXiv

A Stochastic Approach to Maxwell Velocity Distribution via Central Limit Theorem and Boltzmann's Entropy Formula

Maxwell's velocity distribution is known to be universally valid across systems and phases. Here we present a new and general derivation that uses the central limit theorem (CLT) of the probability theory. This essentially uses the idea that repeated intermolecular collisions introduce randomness in the velocity change in the individual components of the velocity vector, leading to, by the CLT, a Gaussian distribution. To complete the derivation, we next show that the mean-square velocity or the standard deviation follows exactly from Boltzmann's formula relating entropy to the density of states, thus avoiding the use of the ideal gas equation of state explicitly. We furthermore prove the Maxwell velocity distribution for a system of harmonic oscillators. This derivation provides a further insight into the origin of Boltzmann's constant in the Maxwell velocity distribution and also in the equipartition theorem. We propose that this derivation provides an approach that explains the universality of Maxwell Gaussian distribution.

preprint2022arXiv

Diffusion in the Presence of Correlated Returns in a Two-dimensional Energy Landscape and non-Monotonic Friction Dependence: Examination of Simulation Results by a Random Walk Model

Diffusion in a multidimensional energy surface with minima and barriers is a problem of importance in statistical mechanics and also has wide applications, such as protein folding. To understand it in such a system, we carry out theory and simulations of a tagged particle moving on a two-dimensional periodic potential energy surface, both in the presence and absence of noise. Langevin dynamics simulations at multiple temperatures are carried out to obtain the diffusion coefficient of a solute particle. Friction is varied from zero to large values. Diffusive motion emerges in the limit of long times, even in the absence of noise, although the trajectory is found to remain correlated over a long time. This correlation is manifested in correlated returns to the starting minima following a scattering by surrounding maxima. Noise destroys this correlation, induces chaos, and increases diffusion at small friction. Diffusion thus exhibits a non-monotonic friction dependence at the intermediate value of the damping, ultimately converging to our theoretically predicted value. The latter is obtained by using the well-established relation between diffusion and random walk. An excellent agreement is obtained between theory and simulations in the high friction limit, but not so in the intermediate regime. The rate of escape from one cell to another is obtained from the multidimensional rate theory of Langer. We find that enhanced dimensionality plays an important role. In order to quantify the effects of noise on the potential-imposed coherence on the trajectories, we calculate the Lyapunov exponent. At small friction values, the Lyapunov exponent mimics the friction dependence of the rate.

preprint2022arXiv

Examination of Boltzmann's H-Function: Dimensionality and Interaction Sensitivity Dependence, and a comment on his H-Theorem

Boltzmann's H-Theorem, formulated 150 years ago in terms of H-function that also bears his name, is one of the most celebrated theorems of science and paved the way for the development of nonequilibrium statistical mechanics. Nevertheless, quantitative studies of the H-function, denoted by H(t), in realistic systems are relatively scarce because of the difficulty of obtaining the time-dependent momentum distribution analytically. Also, the earlier attempts proceeded through the solution of Boltzmann's kinetic equation, which was hard. Here we investigate, by direct molecular dynamics simulations and analytic theory, the time dependence of H(t). We probe the sensitivity of nonequilibrium relaxation to interaction potential and dimensionality by using the H-function H(t). We evaluate H(t) for three different potentials in all three dimensions and find that it exhibits surprisingly strong sensitivity to these factors. The relaxation of H(t) is long in 1D, but short in 3D. We obtain, for the first time, a closed-form analytic expression for H(t) using the solution of the Fokker-Planck equation for the velocity space probability distribution and compare its predictions with the simulation results. Interestingly, H(t) is found to exhibit linear response when vastly different initial nonequilibrium conditions are employed. The oft-quoted relation of H-function with Clausius's entropy theorem is discussed.

preprint2020arXiv

Attainment of Herd Immunity: Mathematical Modelling of Survival Rate

We study the influence of the rate of the attainment of herd immunity (HI), in the absence of an approved vaccine, on the vulnerable population. We essentially ask the question: how hard the evolution towards the desired herd immunity could be on the life of the vulnerables? We employ mathematical modelling (chemical network theory) and cellular automata based computer simulations to study the human cost of an epidemic spread and an effective strategy to introduce HI. Implementation of different strategies to counter the spread of the disease requires a certain degree of quantitative understanding of the time dependence of the outcome. In this paper, our main objective is to gather understanding of the dependence of outcome on the rate of progress of HI. We generalize the celebrated SIR model (Susceptible-Infected-Removed) by compartmentalizing the susceptible population into two categories- (i) vulnerables and (ii) resilients, and study dynamical evolution of the disease progression. We achieve such a classification by employing different rates of recovery of vulnerables vis-a-vis resilients. We obtain the relative fatality of these two sub-categories as a function of the percentages of the vulnerable and resilient population, and the complex dependence on the rate of attainment of herd immunity. Our results quantify the adverse effects on the recovery rates of vulnerables in the course of attaining the herd immunity. We find the important result that a slower attainment of the HI is relatively less fatal. However, a slower progress towards HI could be complicated by many intervening factors.

preprint2020arXiv

Dynamical Theory and Cellular Automata Simulations of Pandemic Spread: Understanding Different Temporal Patterns of Infections

Here we propose and implement a generalized mathematical model to find the time evolution of population in infectious diseases and apply the model to study the recent COVID-19 pandemic. Our model at the core is a non-local generalization of the widely used Kermack-McKendrick(KM) model where the susceptible(S) population evolves into two other categories, namely infectives(I) and removed(R). This is the well-known SIR model in which we further divide both S and I into high and low risk categories. We first formulate a set of non-local dynamical equations for the time evolution of distinct population distributions under this categorization in an attempt to describe the general scenario of infectious disease progression. We then solve the non-linear coupled differential equations-(i) numerically by the method of propagation, and (ii) a more flexible and versatile cellular automata (CA) simulation which provides a coarse-grained description of the generalized non-local model. In order to account for multiple factors such as role of spreaders before containment, we introduce a time dependent rate which appears to be essential to explain the sudden spikes before the plateau observed in many cases (for example like China). We demonstrate how this generalized approach allows us to handle the effects of (i) time-dependence of the rate-constants of spread, (ii) different population density, (iii) the age ratio, (iv) quarantine, (v) lockdown, and (vi) social distancing. Our study allows us to make certain predictions regarding the nature of spread with respect to several external parameters, treated as control variables. Analysis of the model clearly shows that due to the strong heterogeneity in the epidemic process originating from the distribution of initial infectives, the theory must be local in character but at the same time connect to a global perspective.

preprint2020arXiv

Inherent structure analysis reveals origin of breakdown of Stokes-Einstein relation in aqueous binary mixtures

We show by inherent structure (IS) analysis that the sharp composition dependent breakdown of the Stokes-Einstein relation correlates surprisingly well with an equally sharp non-monotonic variation in the average inherent structure (IS) energy of these mixtures. Further IS analysis reveals the existence of a unique ground state, stabilized by the optimum number of H-bonds at this composition. The surprisingly sharp turnaround behaviour observed in the effective hydrodynamic radius can be traced back to the formation of low energy equilibrium structures at specific compositions.

preprint2020arXiv

Origin of Multiple Infection Waves in a Pandemic: Effects of Inherent Susceptibility and External Infectivity Distributions

Two factors that are often ignored but could play a crucial role in the progression of an infectious disease are the distributions of inherent susceptibility ($σ_{inh}$) and external infectivity ($ι_{ext}$), in a given population. While the former is determined by the immunity of an individual towards a disease, the latter depends on the duration of exposure to the infection. We model the spatio-temporal propagation of a pandemic using a generalized SIR (Susceptible-Infected-Removed) model by introducing the susceptibility and infectivity distributions to understand their combined effects, which appear to remain inadequately addressed till date. We consider the coupling between $σ_{inh}$ and $ι_{ext}$ through a new Critical Infection Parameter (CIP) ($γ_c$). We find that the neglect of these distributions, as in the naive SIR model, results in an overestimation of the amount of infection in a population, which leads to incorrect (higher) estimates of the infections required to achieve the herd immunity threshold. Additionally, we include the effects of seeding of infection in a population by long-range migration. We solve the resulting master equations by performing Kinetic Monte Carlo Cellular Automata (KMC-CA) simulations. Importantly, our simulations can reproduce the multiple infection peak scenario of a pandemic. The latent interactions between disease migration and the distributions of susceptibility and infectivity can render the progression a character vastly different from the naive SIR model. In particular, inclusion of these additional features renders the problem a character of a living percolating system where the disease cluster survives by migrating from region to region.

preprint2020arXiv

Study of Entropy-Diffusion Relation in a Deterministic Hamiltonian System through Microscopic Analysis

Although an intimate relation between entropy and diffusion has been advocated for many years and even seems to have been verified in theory and experiments, a quantitatively reliable study, and any derivation of an algebraic relation between the two does not seem to exist. Here we explore the nature of this entropy-diffusion relation in three deterministic systems where an accurate estimate of both can be carried out. We study three deterministic model systems, (a) the motion of a single point particle with constant energy in a two-dimensional periodic potential energy landscape, (b) the same in regular Lorentz gas where a point particle with constant energy moves between collisions with hard disc scatterers and (c) motion of a point particle among the boxes with small apertures. These models, introduced by Zwanzig, exhibit diffusive motion in the limit where ergodicity is shown to exist. We then explore the diffusion-entropy relation by an accurate calculation of both diffusion and entropy for the aforementioned model systems. We estimate the self-diffusion coefficient of the particle by employing computer simulations and entropy by quadrature using Boltzmann's formula. We observe an interesting crossover in the diffusion-entropy relation in some specific regions which is attributed to the emergence of correlated returns. The crossover could herald a breakdown of the Rosenfeld-like exponential scaling between the two, as observed at low temperatures. Later, we modify the scaling relation to account for the correlated motions and present a detailed analysis of the dynamical entropy obtained via Lyapunov exponent which is rather an important quantity in the study of deterministic systems.

preprint2016arXiv

Anomalous dimensionality dependence of diffusion in a rugged energy landscape : How pathological is one dimension ?

Rugged (or, rough) energy landscape (REL) with spatially distributed maxima and minima are often employed in applications of physics, chemistry and biology (enzyme kinetics, protein folding, diffusion in disordered solids, transport in organic semiconductors, relaxation in random spin systems, in supercooled liquids and glasses). Sometimes the system needs to be modeled as a random walker in high dimensions (like in protein folding/unfolding) where dimensions could be the distances between different amino acid residues (as in unfolding of HP-36). Nevertheless, most of the theoretical studies of these phenomena still employ a one dimensional description. This is despite the prediction that in a rough (or, rugged) energy landscape (REL), diffusion in one dimension (1d) is predicted to be pathologically different from any higher dimension with the increased chance of encountering broken ergodicity (Stein and Newman, 2012). We explore the dimensionality dependent diffusion on REL by carrying out an effective medium approximation based analytical calculations and compare them with the available computer simulation results. We find that at intermediate level of ruggedness (assumed to have a Gaussian distribution), where diffusion is well-defined, the value of the effective diffusion coefficient depends on dimensionality and changes (increases) by several factors (~5-10) in going from 1d to 2d. In contrast, the changes in subsequent transitions (like 2d to 3d and 3d to 4d and so on) are far more modest, of the order of 10-20% only. When ruggedness is given by random traps with an exponential distribution of barrier heights, the mean square displacement is sub-diffusive (a well-known result), but the growth of MSD is described by different exponents in one and higher dimensions. The exponent of growth is larger in higher dimensions than in 1d.

preprint2016arXiv

Crossover behavior in the distance dependence of hydrophobic force law

Understanding about both the range and the strength of the effective force between two hydrophobic surfaces suspended in water is important in many areas of natural science but unfortunately has remained imperfect. Even the experimental observations have not been explained quantitatively. Here we find by varying distance (d) between two hydrophobic walls in computer simulations of water that the force exhibits a bi-exponential distance dependence. The long range part of the force can be fitted to an exponential force law with correlation length of 2 nm while the short range part displays a correlation length of only 0.5 nm. The crossover from shorter range to longer range force law is rather sharp. We show that the distance dependence of the tetrahedrality order parameter provides a reliable marker of the force law, and exhibits similar distance dependence.

preprint2016arXiv

Enhanced density fluctuations in water-ethanol mixtures at low ethanol concentrations: Temperature dependent studies

Since the structural transformations observed in water-ethanol binary mixtures are apparently driven by relatively weak intermolecular forces (like hydrophobicity and hydrogen bonding) that often cooperate to form self assembled structures, one expects the aggregation properties to show strong temperature dependence. We study the effect of temperature on the formation of transient ethanol clusters as well as on the dynamic density heterogeneity induced in the system due to such clustering. The dynamic heterogeneity is expected to occur on small length scales with short lifetime and both are expected to be temperature dependent. Indeed, a major finding of the work is strong temperature dependence of the extent of structural heterogeneity. Distinct signature of static and dynamic heterogeneity of ethanol molecules is also found to appear with lowering of temperature. This is attributed to the formation of transient ethanol clusters that are known to exhibit considerably small lifetime (order of a few picosecond). The transient dynamical features of dynamic heterogeneity are expected to affect those relaxation processes occurring at sub-picosecond time scales. On the other hand, strong temperature dependence of micro-structure formation can be anticipated to be due to enhanced structural order stimulated in the system with lowering of temperature. Present analyses reveal a number of interesting features which were not explored beforehand in this widely studied binary mixture.

preprint2016arXiv

Exciton transfer dynamics and quantum diffusion in a lattice of two level systems: Interplay between transport and coherent population transfer dynamics

We study excitation transfer dynamics in a lattice of two level systems characterized by dynamic disorder. The diagonal and off-diagonal energy disorders arise from the coupling of system and bath. We consider both the same and the independent bath limits. In case of independent bath all diagonal and off-diagonal bath coupling elements fluctuate independently of each other and the dynamics is complicated. We obtain the time dependent population distribution by solving quantum stochastic Liouville equation (QSLE) derived by Kubo. The main result of our study is both the population transfer dynamics and the mean square displacement of the exciton behave the similar way in the same and independent bath cases in the Markovian limit. However, these two baths can give rise to markedly different behavior in the non-Markovian limit where coherent transport becomes important. There are also several additional new results as follows. (i) Exciton migration remains coherent all the time for an average, non-zero off-diagonal coupling value J for the same bath case while it becomes incoherent for independent bath case in the Markovian limit. (ii) An oscillatory behavior of the population transfer dynamics supports the coherent mode of transfer of exciton. (iii) Agreement with available analytical expression of mean squared displacement is good in the Markovian limit for independent bath case with off-diagonal fluctuation but only qualitative in the non-Markovian limit for which no complete analytical solution is available. (iv) Transition from coherent to incoherent transport is observed in the independent bath case with diagonal fluctuation when the bath is made progressively more Markovian. (v) The correlation time of the bath plays a unique role in dictating the diffusive spread that is not anticipated in a Markovian treatment.

preprint2016arXiv

How universal is the Lindemann criteria in melting of Lennard-Jones polydisperse solids?

It is commonly believed that melting occurs when mean square displacement (MSD) of a particle of crystalline solid exceeds a threshold value. This is known as the Lindemann criterion, first introduced in the year of 1910 by Lindemann. However, Chakravarty et al. demonstrated that this common wisdom is inadequate because the MSD at melting can be temperature dependent when pressure is also allowed to vary along the coexistence line of the phase diagram [Chakravarty C, Debenedetti P G and Stillinger F H 2007 J. Chem. Phys.126 204508]. We show here by extensive molecular dynamics simulation of both two and three dimensional polydisperse Lennard-Jones solids that particles on the small and large limits of size distribution exhibit substantially different Lindemann ratio at melting. Despite all the dispersion in MSD, melting is found to be first order in both the dimensions at 5-10% dispersity in size.

preprint2016arXiv

Lambda-type sharp rise in the widths of Raman and infra-red line shape near the Widom line in super-critical water above its gas-liquid critical temperature

A lambda-type divergent rise of Raman linewidth of liquid nitrogen near its critical temperature has been a subject of many discussions in the past[1-5]. Here we explore the possibility of such an anomaly in infra-red and Raman spectroscopy of super-critical water (SCW) by varying the density across the Widom line just above its critical temperature. Vibrational phase relaxation is expected to be a sensitive probe of fluid dynamics. We carry out computer simulations of two different model potentials (SPC/E and TIP4P/2005) to obtain the necessary time correlation functions. An additional feature of this work is a quantum chemical calculation of the anharmonicity parameter that largely controls frequency fluctuations. We find a sharp rise in the vibrational relaxation rate (or the line widths) for both the models as we travel across the Widom line. The rise is noticeably less sharp in water than in nitrogen. We attribute this difference to the faster relaxation rate in water. We demonstrate that the anomalous rise is due to sharp increase in the mean square frequency fluctuation and not due to any dynamical critical slowing down because the time constant of the normalized frequency-frequency time correlation functions show lack of any noticeable change as we move across the Widom line. We present an explanation of the observed results using the mode coupling theory of liquid dynamics.

preprint2016arXiv

Nanoscale heterogeneous phase separation kinetics in binary mixtures: Multistage dynamics

In order to find a measure of the dynamical features of phase separation kinetics during spinodal decomposition of a liquid binary mixture (like water and cyclohexane , water and 2,6 lutidiene or methanol and cyclohexane), we study both the initial fast exponential-like growth (the Cahn-Hilliard regime) and the subsequent cross-over to a much slower, non-exponential long time growth (the so-called scaling regime), by atomistic molecular dynamics (MD) simulation of a structure breaking binary liquid mixture. In particular, we combine our MD simulations with a coarse grained multi scale modelling (CGMSM) capable of capturing both length and time scales of phase separation kinetics within simulation box. The system is quenched from a higher temperature to two lower temperatures well below the coexistence temperature of the phase diagram. We observe a multiscale phase separation dynamics. Initially the growth is exponential up to a regime of 80-100ps having strong dependence over quench depth. Subsequently a cross-over regime appears where the dynamics slows down considerably due to coarsening through a power law phase. For deeper quench power law growth dominates over the initial exponential and the cross-over becomes transient. We find that for the present parameter values for the binary mixture, the initial rapid growth of structure formation is practically over within 200 ps which is followed by slow structural coarsening. When scaled by the respective viscosities, this time translates to 50-200 ns for water-lutidine binary mixture. The last part of dynamics may extend into ms. The dynamics of phase separation is slowest in regions that have equitable distribution of the two species, and one can observe the signatures of "up-hill diffusion" that is a trade mark of spinodal decomposition.

preprint2016arXiv

Signatures of quantum coherence in the optical line shape of an exciton in the presence of dynamic disorder

We address the effects of quantum coherences on the optical line shape of an exciton in the presence of dynamic disorder. We consider a one-dimensional excitonic system that consists of two levels placed at regular intervals. Detailed analytical calculations of line shape have been carried out by using Kubo's stochastic Liouville equation (K-QSLE). We make use of the observation that in the site representation, the Hamiltonian of our system with constant off-diagonal coupling J is a tridiagonal Toeplitz matrix (TDTM) whose eigenvalues and eigen functions are known analytically. This identification is particularly useful for long chains where the eigen values of TDTM help to understand crossover between static and fast modulation limits. We summarize the new results as follows. (i) In the slow modulation limit when the bath correlation time is large, the effects of spatial correlation are not negligible. Here the line shape is broadened and the number of peak increases beyond the ones obtained from TDTM (constant off-diagonal coupling element J and no fluctuation). (ii) However, in the fast modulation limit when the bath correlation time is small, the spatial correlation is less important (iii) Importantly, we find that the line shape can capture that quantum coherence affects in the two limits differently.

preprint2016arXiv

Water rotational jump driven large amplitude molecular motions of nitrate ions in aqueous potassium nitrate solution

Molecular dynamics simulations of aqueous potassium nitrate solution reveal a highly complex rotational dynamics of nitrate ions where, superimposed on the expected continuous Brownian motion, are large amplitude angular jumps that are coupled to and at least partly driven by similar large amplitude jump motions in water molecules which are associated with change in the hydrogen bonded water molecule. These jumps contribute significantly to rotational and translational motions of these ions. We explore the detailed mechanism of these correlated (or, coupled) jumps and introduce a new time correlation function to decompose the coupled orientational- jump dynamics of solvent and solute in the aqueous electrolytic solution. Time correlation function provides for the unequivocal determination of the time constant involved in orientational dynamics originating from making and breaking of hydrogen bonds. We discover two distinct mechanisms-both are coupled to density fluctuation but are of different types.

preprint2015arXiv

Control of human immune response function by T-cell population fluctuation and relaxation dynamics

Clinical studies have indicated that in malignant surveillances fluctuations in the population of certain effector T-cell repertoire become suppressed. Motivated by such observations and in an attempt to quantify adaptive human response to pathogens, we define an immune response function (IMRF) in terms of mean square fluctuations of T-cell concentrations. We employ a recently developed kinetic model of T-cell regulation that contains the essential immunosuppressive effects of vitamin-D. We employ Gillespie algorithm to make the first study of fluctuations along the stochastic trajectories. This fluctuation-based IMRF can differentiate responses of different individuals after pathogenic incursion both under healthy and disease conditions. We find that relative fluctuations in T-cells (and hence IMRF) are different in strongly regulated (malignant prone) and weakly regulated (autoimmune prone) regions. The cross-over from one steady state (weakly regulated) to the other (strongly regulated) is accompanied by a divergence-like growth in the fluctuation of both the effector and regulatory T-cell concentration over a wide range of pathogenic stimulation, displaying a dynamical phase transition like behavior. The growth in fluctuation in this desired immune response regime is found to arise from an intermittent fluctuation between regulatory and effector T-cells that results in a bimodal distribution of population of each, indicating bistability. The signature of intermittent behavior is further confirmed by calculating the power spectrum of the corresponding fluctuation of time correlation function. The calculated time correlation functions of fluctuations show that the slow fluctuation causes the bistabilty in healthy state. Thus, in diseases diagnosis process, such steady state response parameters can provide immense information which might become helpful to define an immune status.

preprint2015arXiv

Mode coupling theory of electrolyte dynamics: Time dependent diffusion, dynamic structure factor and solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) dynamic structure factor of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics (BD) with implicit water molecules and molecular dynamics (MD) method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding dynamic structure factor evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of total dynamic structure factor F(k,t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behaviour allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times- the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. This solvation time correlation and ion-ion dynamic structure factor indeed exhibits highly interesting, non-trivial dynamical behaviour at intermediate to longer times that require further experimental and theoretical studies.

preprint2015arXiv

Orientational order as the origin of the long-range hydrophobic effect

The long range attractive force between two hydrophobic surfaces immersed in water is observed to decrease exponentially with their separation -- this distance-dependence of effective force is known as the hydrophobic force law (HFL). We explore the microscopic origin of HFL by studying distance-dependent attraction between two parallel rods immersed in 2D Mercedes Benz model of water. This model is found to exhibit a well-defined HFL. Although the phenomenon is conventionally explained by density-dependent theories, we identify orientation, rather than density, as the relevant order parameter. The range of density variation is noticeably shorter than that of orientational heterogeneity. The latter is comparable to the observed distances of hydrophobic force. At large separation, attraction between the rods arises primarily from a destructive interference among the inwardly propagating oppositely oriented heterogeneity generated in water by the two rods. As the rods are brought closer, the interference increases leading to a decrease in heterogeneity and concomitant decrease in free energy of the system, giving rise to the effective attraction. We notice formation of hexagonal ice-like structures at the onset of attractive region which suggests that metastable free energy minimum may play a role in the origin of HFL.

preprint2015arXiv

Relationship between Entropy and Diffusion: A statistical mechanical derivation of Rosenfeld expression for a rugged energy landscape

Diffusion-a measure of dynamics, and entropy-a measure of disorder in the system, are found to be intimately correlated in many systems, and the correlation is often strongly non-linear. We explore the origin of this complex dependence by studying diffusion of a point Brownian particle on a model potential energy surface characterized by ruggedness. If we assume that the ruggedness has a Gaussian distribution then for this model, one can obtain the excess entropy exactly for any dimension. By using the expression for the mean first passage time (MFPT), we present a statistical mechanical derivation of the well-known and well-tested scaling relation proposed by Rosenfeld between diffusion and excess entropy. In anticipation that Rosenfeld diffusion-entropy scaling (RDES) relation may continue to be valid in higher dimensions (where the mean first passage time approach is not available), we carry out an effective medium approximation (EMA) based analysis of the effective transition rate, and hence of the effective diffusion coefficient. We show that the EMA expression can be used to derive the RDES scaling relation for any dimension higher than unity. However, RDES is shown to break down in the presence of spatial correlation among the energy landscape values.

preprint2014arXiv

A Comparative Study of Protein Unfolding in Aqueous Urea and DMSO Solutions: Surface Polarity, Solvent Specificity and Sequence of Secondary Structure Melting

Elucidation of possible pathways between folded (native) and unfolded states of a protein is a challenging task, as the intermediates are often hard to detect. Here we alter the solvent environment in a controlled manner by choosing two different co-solvents of water, urea and dimethyl sulphoxide (DMSO), and study unfolding of four different proteins to understand the respective sequence of melting by computer simulation methods. We indeed find interesting differences in the sequence of melting of alpha-helices and beta-sheets in these two solvents. For example, at 8M urea solution, beta-sheet parts of a protein is found to unfold preferentially, followed by the unfolding of alpha helices. In contrast, 8M DMSO solution unfolds alpha helices first, followed by the separation of beta-sheets for majority of proteins. Sequence of unfolding events in four different alpha/beta proteins and also in chicken villin head piece (HP-36) both in urea and DMSO solution demonstrate that the unfolding pathways are determined jointly by relative exposure of polar and non-polar residues of a protein and the mode of molecular action of a solvent on that protein.

preprint2014arXiv

A stochastic chemical dynamic approach to correlate autoimmunity and optimal vitamin-D range

Motivated by several recent experimental observations that vitamin-D could interact with antigen presenting cells (APCs) and T-lymphocyte cells (T-cells) to promote and to regulate different stages of immune response, we developed a coarse grained kinetic model in an attempt to quantify the role of vitamin-D in immunomodulatory responses. Our kinetic model, developed using the ideas of chemical network theory, leads to a system of nine coupled equations that we solve both by direct and by stochastic (Gillespie) methods. Both the analyses consistently provide detail information on the dependence of immune response to the variation of critical rate parameters. We find that although vitamin-D plays a negligible role in the initial immune response, it exerts a profound influence in the long term, especially in helping the system to achieve a new, stable steady state. The study explores the role of vitamin-D in preserving an observed bistability in the phase diagram (spanned by system parameters) of immune regulation, thus allowing the response to tolerate a wide range of pathogenic stimulation which could help in resisting autoimmune diseases. We also study how vitamin-D affects the time dependent population of dendritic cells that connect between innate and adaptive immune responses. Variations in dose dependent response in anti-inflammatory and pro-inflammatory T-cell populations to vitamin-D correlate well with recent experimental results. Our kinetic model allows for an estimation of the range of optimum level of vitamin-D required for smooth functioning of the immune system and for control of both hyper-regulation and inflammation. Most importantly, the present study reveals that an overdose or toxic level of vitamin-D or any steroid analogue could give rise to too large a tolerant response, leading to an inefficacy in adaptive immune function.

preprint2014arXiv

Anomalous Power Law Decay in Solvation Dynamics of DNA: A Mode Coupling Theory Analysis of Ion Contribution

Several time domain fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Lastly, we propose several experiments to check the prediction of the present theoretical work.

preprint2014arXiv

Correlation between thermodynamic anomalies and pathways of ice nucleation in supercooled water

The well-known classical nucleation theory (CNT) for the free energy barrier towards formation of a nucleus of critical size of the new stable phase within the parent metastable phase fails to take into account the influence of other metastable phases having density/order intermediate between the parent metastable phase and the final stable phase. This lacuna can be more serious than capillary approximation or spherical shape assumption made in CNT. This issue is particularly significant in ice nucleation because liquid water shows rich phase diagram consisting of two (high and low density) liquid phases in supercooled state. The explanations of thermodynamic and dynamic anomalies of supercooled water often invoke the possible influence of a liquid-liquid transition between two metastable liquid phases. To investigate both the role of thermodynamic anomalies and presence of distinct metastable liquid phases in supercooled water on ice nucleation, we employ density functional theoretical approach to find nucleation free energy barrier in different regions of phase diagram. The theory makes a number of striking predictions, such as a dramatic lowering of nucleation barrier due to presence of a metastable intermediate phase and crossover in the dependence of free energy barrier on temperature near liquid-liquid critical point. These predictions can be tested by computer simulations as well as by controlled experiments.

preprint2014arXiv

Diffusion in a rough potential revisited

Rugged energy landscapes find wide applications in diverse fields ranging from astrophysics to protein folding. We study the dependence of diffusion coefficient $(D)$ of a Brownian particle on the distribution width $(\varepsilon)$ of randomness in a Gaussian random landscape by simulations and theoretical analysis. We first show that the elegant expression of Zwanzig [PNAS, 85, 2029 (1988)] for $D(\varepsilon)$ can be reproduced exactly by using the Rosenfeld diffusion-entropy scaling relation. Our simulations show that Zwanzig's expression overestimates $D$ in an uncorrelated Gaussian random lattice - differing by almost an order of magnitude at moderately high ruggedness. The disparity originates from the presence of "three-site traps" (TST) on the landscape -- which are formed by the presence of deep minima flanked by high barriers on either side. Using mean first passage time formalism, we derive a general expression for the effective diffusion coefficient in the presence of TST, that quantitatively reproduces the simulation results and which reduces to Zwanzig's form only in the limit of infinite spatial correlation. We construct a continuous Gaussian field with inherent correlation to establish the effect of spatial correlation on random walk. The presence of TSTs at large ruggedness $(\varepsilon \gg k_{B}T)$ give rise to an apparent breakdown of ergodicity of the type often encountered in glassy liquids.

preprint2014arXiv

Fluctuating micro-heterogeneity in water - tert-butyl alcohol mixtures and lambda-type divergence of the mean cluster size with phase transition-like multiple anomalies

Water - tert-butyl alcohol (TBA) binary mixture exhibits a large number of thermodynamic and dynamic anomalies. These anomalies are observed at surprisingly low TBA mole fraction, with $x_{\text{TBA}} \approx 0.03 - 0.07$. We demonstrate here that the origin of the anomalies lies in the local structural changes that occur due to self-aggregation of TBA molecules. We observe a percolation transition of the TBA molecules at $x_{\text{TBA}} \approx 0.05$. We note that "islands" of TBA clusters form even below this mole fraction, while a large spanning cluster emerges above that mole fraction. At this percolation threshold, we observe a lambda-type divergence in the fluctuation of the size of the largest TBA cluster, reminiscent of a critical point. Alongside, the structure of water is also perturbed, albeit weakly, by the aggregation of TBA molecules. There is a monotonic decrease in the tetrahedral order parameter of water, while the dipole moment correlation shows a weak non-linearity. Interestingly, water molecules themselves exhibit a reverse percolation transition at higher TBA concentration, $x_{\text{TBA}} \approx 0.45$, where large spanning water clusters now break-up into small clusters. This is accompanied by significant divergence of the fluctuations in the size of largest water cluster. This second transition gives rise to another set of anomalies around. Both the percolation transitions can be regarded as manifestations of Janus effect at small molecular level.

preprint2014arXiv

Hydrophobic hydration driven self-assembly of Curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting strongly with hydrophilic groups at the surface of curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer.

preprint2014arXiv

Melting/freezing transition in polydisperse Lennard-Jones system: Remarkable agreement between predictions of inherent structure, bifurcation phase diagram, Hansen-Verlet rule and Lindemann criteria

We use polydispersity in size as a control parameter to explore certain aspects of melting and freezing transitions in a system of Lennard-Jones spheres. Both analytical theory and computer simulations are employed to establish a potentially interesting relationship between observed terminal polydispersity in Lennard-Jones polydisperse spheres and prediction of the same in the integral equation based theoretical analysis of liquid-solid transition. As we increase polydispersity, solid becomes inherently unstable because of the strain built up due to the size disparity. This aspect is studied here by calculating the inherent structure (IS) calculation. With polydispersity at constant volume fraction we find initially a sharp rise of the average IS energy of the crystalline solid until transition polydispersity, followed by a cross over to a weaker dependence of IS energy on polydispersity in the amorphous state. This cross over from FCC to amorphous state predicted by IS analysis agrees remarkably well with the solid-liquid phase diagram (with extension into the metastable phase) generated by non-linear integral equation theories of freezing. Two other well-known criteria of freezing/melting transitions, the Hansen-Verlet rule of freezing and the Lindemann criterion of melting are both shown to be remarkably in good agreement with the above two estimates. Together they seem to indicate a small range of metastability in the liquid-solid transition in polydisperse solids.

preprint2013arXiv

Frequency dependence of specific heat in supercooled liquid water and emergence of correlated dynamics

Molecular origin of the well-known specific heat anomaly in supercooled liquid water is investigated here by using extensive computer simulations and theoretical analyses. A rather sharp increase in the values of isobaric specific heat with lowering temperature and the weak temperature dependence of isochoric specific heat in the same range are reproduced in simulations. We calculated the spatiotemporal correlation among temperature fluctuations and examined the frequency dependent specific heat. The latter shows a rapid growth in the low frequency regime as temperature is cooled below 270 K. In order to understand the microscopic basis of this increase, we have performed a shell wise decomposition of contributions of distant molecules to the temperature fluctuations in a central molecule. This decomposition reveals the emergence, at low temperatures, of temporally slow, spatially long ranged large temperature fluctuations. The temperature fluctuation time correlation function (TFCF) can be fitted to a William-Watts stretched exponential form with the stretching parameter close to 0.6 at low temperatures, indicating highly non-exponential relaxation. Temperature dependence of the relaxation time of the correlation function can be fitted to Vogel-Fulcher-Tamermann expression which provides a quantitative measure of the fragility of the liquid. Interestingly, we find that the rapid growth in the relaxation time of TFCF with lowering temperature undergoes a sharp crossover from a markedly fragile state to a weakly fragile state around 220 K.

preprint2013arXiv

Kinetic proofreading at single molecular level: Aminoacylation of tRNA^{Ile} and the role of water as an editor

Proofreading/editing in protein synthesis is essential for accurate translation of information from the genetic code. In this article we present a theoretical investigation of efficiency of a kinetic proofreading mechanism that employs hydrolysis of the wrong substrate as the discriminatory step in enzyme catalytic reactions. We consider aminoacylation of tRNA^{Ile} which is a crucial step in protein synthesis and for which experimental results are now available. We present an augmented kinetic scheme and then employ methods of stochastic simulation algorithm to obtain time dependent concentrations of different substances involved in the reaction and their rates of formation. We obtain the rates of product formation and ATP hydrolysis for both correct and wrong substrates (isoleucine and valine in our case), in single molecular enzyme as well as ensemble enzyme kinetics. The present theoretical scheme correctly reproduces (i) the amplitude of the discrimination factor in the overall rates between isoleucine and valine which is obtained as (1.8 \times 10^2).(4.33 \times 10^2) = 7.8 \times 10^4, (ii) the rates of ATP hydrolysis for both Ile and Val at different substrate concentrations in the aminoacylation of tRNA^{Ile}. The present study shows a non-michaelis type dependence of rate of reaction on tRNA^{Ile} concentration in case of valine. The editing in steady state is found to be independent of amino acid concentration. Interestingly, the computed ATP hydrolysis rate for valine at high substrate concentration is same as the rate of formation of Ile-tRNA^{Ile} whereas at intermediate substrate concentration the ATP hydrolysis rate is relatively low.

preprint2013arXiv

Nucleation of a stable solid from melt in the presence of multiple metastable intermediate phases: Wetting, Ostwald step rule and vanishing polymorphs

In many systems, nucleation of a stable solid may occur in the presence of other (often more than one) metastable phases. These may be polymorphic solids or even liquid phases. In such cases, nucleation of the solid phase from the melt may be facilitated by the metastable phase because the latter can "wet" the interface between the parent and the daughter phases, even though there may be no signature of the existence of metastable phase in the thermodynamic properties of the parent liquid and the stable solid phase. Straightforward application of classical nucleation theory (CNT) is flawed here as it overestimates the nucleation barrier since surface tension is overestimated (by neglecting the metastable phases of intermediate order) while the thermodynamic free energy gap between daughter and parent phases remains unchanged. In this work we discuss a density functional theory (DFT) based statistical mechanical approach to explore and quantify such facilitation. We construct a simple order parameter dependent free energy surface that we then use in DFT to calculate (i) the order parameter profile, (ii) the overall nucleation free energy barrier and (iii) the surface tension between the parent liquid and the metastable solid and also parent liquid and stable solid phases. The theory indeed finds that the nucleation free energy barrier can decrease significantly in the presence of wetting. This approach can provide a microscopic explanation of Ostwald step rule and the well-known phenomenon of "disappearing polymorphs" that depends on temperature and other thermodynamic conditions. Theory reveals a diverse scenario for phase transformation kinetics some of which may be explored via modern nanoscopic synthetic methods.

preprint2013arXiv

Solid-solid collapse transition in a two dimensional model molecular system

Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real difficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises from the interplay between different energy and length scales involved in molecular systems, often mediated through a solvent. In this work we employ Monte Carlo (MC) simulations to study the collapse transition in a model molecular system interacting via both isotropic as well as anisotropic interactions having different length and energy scales. The model we use is known as Mercedes-Benz (MB) which for a specific set of parameters sustains three solid phases: honeycomb, oblique and triangular. In order to study the temperature induced collapse transition, we start with a metastable honeycomb solid and induce transition by heating. High density oblique solid so formed has two characteristic length scales corresponding to isotropic and anisotropic parts of interaction potential. Contrary to the common believe and classical nucleation theory, interestingly, we find linear strip-like nucleating clusters having significantly different order and average coordination number than the bulk stable phase. In the early stage of growth, the cluster grows as linear strip followed by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate balance between two types of interactions which enables the dominance of stabilizing energy over the destabilizing surface energy. The nuclei of stable oblique phase are wetted by intermediate order particles which minimizes the surface free energy. We observe different pathways for pressure and temperature induced transitions.

preprint2013arXiv

Solvent Sensitivity of Protein Unfolding: Study of Chicken Villin Headpiece Subdomain in Water-Ethanol and Water-DMSO Mixtures

In the present work we study and compare unfolding of a small protein, chicken villin headpiece (HP-36) in two different aqueous binary mixtures, namely water-ethanol (EtOH) and water-dimethyl sulphoxide (DMSO). In both the binary mixtures, HP-36 is found to unfold (fully or partially, depending on the mixture) under ambient conditions, that otherwise requires temperature as high as ~600 K to denature in pure aqueous solvent. In all the cases, first step of unfolding is found to be similar, i.e. separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely Phe-7, Phe-11 and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. Subsequent unfolding steps follow different paths in different chemical environments. As both water-DMSO and water-ethanol show composition dependent anomalies, so do the details of unfolding dynamics. With an increase of co-solvent concentration different partially unfolded intermediates are found to be formed in both the cases. This is reflected in a remarkable non-monotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states is particularly attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol and methyl groups of DMSO. While in DMSO the protein gradually attains a completely unfolded state at xDMSO=0.30, unfolding in water-ethanol appears to be more complex and sensitive to solvent composition.

preprint2013arXiv

Strongly first order melting of a two dimensional molecular solid

Melting and freezing transitions in two dimensional systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems; the melting behavior in two dimensional molecular solids is still largely unexplored. In order to understand the role of multiple energy and length scales present in molecular systems on nature of melting transition, here we report computer simulation studies of melting of a two dimensional Mercedes-Benz (MB) system. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram. The computed solid-liquid phase diagram consists of isotropic liquid and two crystalline phases - honeycomb and oblique. In contradiction to the celebrated KTHNY theory, we observe strongly one step first order melting transitions for both the honeycomb and oblique solids. The defects in both solids and liquids near the transition are more complex compared to the atomic systems.

preprint2013arXiv

Vitamin D sensitivity to the immune responses and autoimmunity: A chemical network model study

Although Vitamin D is believed to be involved in a large number of immune responses our understanding of these processes at cellular level remained at infancy. We develop and solve a coarse grained kinetic network model to quantify the effects of variation of vitamin D on human immunity. The system of equations accounts for known inter-relation between active and inactive vitamin D, antigen presenting cells, effector T cells, regulatory T cells and pathogen. Both time dependent and steady state solutions are obtained. The time dependent solution of the system of equations reveals that the immune response is rather strongly regulated in presence of vitamin D. We found quantitatively that lower than optimum levels of concentration of active vitamin D correspond to weak regulation where, once a pathogen/antigen enters the body, the nature of the immune response would be less regulatory and hence more aggressive, or inflammatory. The steady state solution of our model shows that vitamin D enhances the tolerance level of immune system, thereby increasing resistance to autoimmune diseases. Our model and accompanied numerical analyses reveal another important aspect of immunity: While extremely low levels of vitamin D could lead to increased risk of autoimmune responses, an overdose (toxic) level would give rise to too large a tolerant response, leading to the increased risk of tumors and cancerous cell growth.

preprint2012arXiv

Phase transition in magnetically coupled spins on a ring (SOR) model

We have considered a new type of 'XY' model where spins are placed on concentric ring with constant spin density in every ring. The spin executes continuous rotation under a modified Shore-Zwanzig Hamiltonian (J. Chem. Phys. 63, 5445 (1975)). We have performed Monte Carlo simulation using Glauber acceptance criteria. Computations of Binder's cumulant, specific heat and magnetic susceptibility all show presence of a finite temperature order-disorder phase transition in this spin system. The system size dependence of Binder's cumulant suggests the existence of a phase transition with a transition temperature of T* = 1.2. However, we have found no signature of the occurrence of vortex in our SOR model. The absence of hysteresis rules out the possibility of first order phase transition. We have found two "stable" states for T* = 0 phase. The perfectly ordered true ground state is obtained by gradual cooling of the system, while the other is obtained by starting the simulation with a random configuration at T* = 0. This second state has higher energy than the perfectly aligned ground state.

preprint2011arXiv

Dynamic control of catalysis within biological cells

We develop a theory of enzyme catalysis within biological cells where the substrate concentration [S](t) is time dependent, in contrast to the Michaelis-Menten theory that assumes a steady state. We find that the time varying concentration can combine, in a non-linear way, with the ruggedness of the free energy landscape of enzymes (discovered both in single molecule studies and in simulations) to provide a highly efficient switch (or, bifurcation) between two catalytically active states, at a critical substrate concentration. This allows a dynamic control of product synthesis in cell.

preprint2011arXiv

Terminal polydispersity in the crystallization of polydisperse Lennard-Jones liquid

We find through computer simulations that the fractional volume change on freezing of polydisperse Lennard-Jones liquid decreases with increasing polydispersity and approaches zero near a terminal polydispersity of 0.11, independent of temperature. The transition however remains first order at terminal polydispersity. Average inherent structure (IS) energy of the crystalline phase increases nearly quadratically with polydispersity indices (δ) and marked by a crossover to nearly constant IS energy in the amorphous phase.

preprint2010arXiv

A kinetic Ising model study of dynamical correlations in confined fluids: Emergence of both fast and slow time scales

Experiments and computer simulation studies have revealed existence of rich dynamics in the orientational relaxation of molecules in confined systems such as water in reverse micelles, cyclodextrin cavities and nano-tubes. Here we introduce a novel finite length one dimensional Ising model to investigate the propagation and the annihilation of dynamical correlations in finite systems and to understand the intriguing shortening of the orientational relaxation time that has been reported for small sized reverse micelles. In our finite sized model, the two spins at the two end cells are oriented in the opposite directions, to mimic the effects of surface that in real system fixes water orientation in the opposite directions. This produces opposite polarizations to propagate inside from the surface and to produce bulk-like condition at the centre. This model can be solved analytically for short chains. For long chains we solve the model numerically with Glauber spin flip dynamics (and also with Metropolis single-spin flip Monte Carlo algorithm). We show that model nicely reproduces many of the features observed in experiments. Due to the destructive interference among correlations that propagate from the surface to the core, one of the rotational relaxation time components decays faster than the bulk. In general, the relaxation of spins is non-exponential due to the interplay between various interactions. In the limit of strong coupling between the spins or in the limit of low temperature, the nature of relaxation of the spins undergoes a qualitative change with the emergence of a homogeneous dynamics where decay is predominantly exponential, again in agreement with experiments.

preprint2010arXiv

Correlation between inherent structures and phase separation mechanism in binary mixtures

Binary mixtures are known to phase separate via both nucleation and spinodal decomposition, depending on the initial composition and extent of the quench. Here we develop an energy landscape view of phase separation and non-ideality in binary mixtures by exploring their potential energy landscape (PEL) as functions of temperature and composition, by employing a molecular model that promotes structure breaking abilities of the solute-solvent parent binary liquid at low temperatures. PEL that provides the inherent structure (IS) of a system is obtained by removing the kinetic energy (including that of intermolecular vibrations). Broad Distribution of inherent structure energy demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of inherent structure. Over broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy, we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy (<EIS>) exhibits a drop with temperature which signals onset of crystallization in one of the phases while the other remains in the liquid state. The non-ideal composition dependence of viscosity is anti-correlated with average inherent structure energy.

preprint2010arXiv

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of ATP and AMP to ADP by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of LID and NMP domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motions of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis and multi-dimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed RMSD variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder type theoretical model to explain the observed dynamic coupling between the motions of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations.

preprint2010arXiv

Gas-liquid Nucleation at Large Metastability

Nucleation at large metastability is still largely an unsolved problem, although is a problem of tremendous current interest, with wide practical value. It is well-accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this article, we present a powerful alternative formalism to treat nucleation at large supersaturation. This formalism goes over to the classical picture at small supersaturation where CNT is expected to be valid. The new theory is based on an extended set of order parameters in terms of k-th largest liquid-like clusters where k=1 is the largest cluster in the system, k=2 is the second largest cluster and so on. We derive an analytic expression for the free energy of formation of the k-th largest cluster which shows that at large metastability the barrier of growth for the few largest liquid-like clusters disappear, the nucleation becomes collective and the approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts a weaker supersaturation dependence than that of CNT at large metastability. Such a cross-over behavior has indeed been observed in recent experiments but eluded an explanation till now. In order to understand the large numerical difference between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interaction of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for a Lennard-Jones potential, both in two and three dimensions.

preprint2010arXiv

Interplay between the range of attractive potential and metastability in gas-liquid nucleation

We find an interesting interplay between the range of the attractive part of the interaction potential and the extent of metastability (as measured by supersaturation) in gas-liquid nucleation. We explore and exploit this interplay to obtain new insight into nucleation phenomena. Just like its dependence on supersaturation (S), the free energy barrier of nucleation is found to depend strongly on the range of the interaction potential. Actually, the entire free energy surface, F(n), where n is the size of the liquid-like cluster, shows this dependence. The evidences and the reasons for this strong dependence are as follows. (i) The surface tension increases dramatically on increasing the range of interaction potential. In three dimensional Lennard-Jones system, the value of the surface tension increases from 0.494 for a cut-off of 2.5 σ to 1.09 when the full range of the potential is involved. In two dimensional LJ system, the value of the line tension increases from 0.05 to 0.18, under the same variation of the potential range. (ii) The density of the gas phase at coexistence decreases while that of the liquid phase increases substantially on increasing the range of the interaction potential. (iii) As a result of the above, at a given supersaturation S, the size of the critical nucleus and the free energy barrier both increase with increase in the range of interaction potential. (iv) Surprisingly, however, we find that the functional form predicted by the classical nucleation theory (CNT) for the dependence of the free energy barrier on the size of the nucleus to remain valid except at the largest value of S studied. (v) The agreement between CNT prediction and simulated values of the barrier is supersaturation dependent and worsens with increase in the range of interaction potential, and increases above 10 kBT at the largest supersaturation that could be studied.

preprint2010arXiv

Propagating rotational jump events drive liquid-liquid transition in supercooled liquid water

We observe, at low temperature, the appearance of propagating events that originate from the inter-conversion between adjacent four and five coordinated water molecules in supercooled liquid water, resulting in the migration of the coordination number five in a string-like fashion, creating rotational jumps along the way. The length of the connected events increases with lowering temperature. Each propagation event is terminated by a 3-coordinated species, present in a small number at large supercooling, which interacts cooperatively with 5-coordinated molecules to annihilate both the species. We find that these growing dynamical correlations manifest in a divergent-like growth of a non-linear density response function, \times4(t), which is given by a four point time correlation function (FPTCF). The locus of the maximum of \times4(t), when plotted against the time t* of maximum, exhibits a sharp peak precisely at the temperature where the static response functions (specific heat, isothermal compressibility) also show similar sharp, divergent-like, peak. While the decay of population fluctuation time correlation function of both 4- and 5- coordinated species slows down dramatically and a step-like feature of the relaxation becomes evident, the lifetime itself of 5-corrdinated species remains short. These results suggest a new molecular mechanism of low temperature anomalies and of the liquid-liquid transition in terms of initiation, growth and termination of propagating jump and inter-conversion events.

preprint2010arXiv

Role of water in the enzymatic catalysis: study of ATP + AMP -> 2ADP conversion by adenylate kinase

The catalytic conversion ATP + AMP -> 2ADP by the enzyme adenylate kinase (ADK) involves the binding of one ATP molecule to the LID domain and one AMP molecule to the NMP domain. The latter is followed by a phosphate transfer, and then the release of two ADP molecules. We have computed a novel two dimensional configurational free energy surface (2DCFES), with one reaction coordinate each for the LID and the NMP domain motions, with explicit interactions with water. Our computed 2DCFES clearly reveals the existence of a stable half-open-half-closed (HOHC) intermediate state of the enzyme. Cycling of the enzyme through the HOHC state reduces the conformational free energy barrier for the reaction by about 20 kJ/mol. We find that the stability of the half-open-half-closed state (missed in all earlier studies with implicit solvent model) is largely because of the increase of specific interactions of the polar amino acid side chains with water, particularly with the arginine and the histidine residues. Free energy surface of the LID domain is rather rugged, which can conveniently slow down LID's conformational motion, thus facilitating a new substrate capture after the product release in the catalytic cycle.

preprint2009arXiv

Vibrational dynamics and boson peak in a supercooled polydisperse liquid

Vibrational density of states (VDOS) in a supercooled polydisperse liquid is computed by diagonalizing the Hessian matrix evaluated at the potential energy minima for systems with different values of polydispersity. An increase of polydispersity leads to an increase in relative population of the localized high-frequency modes. At low frequencies, the density of states show an excess compared to the Debye squared-frequency law, which has been identified with the boson peak. The height of the boson peak increases with polydispersity. The values of the participation ratio as well as the level spacing statistics demonstrate that the modes comprising the boson peak are largely delocalized. Interestingly, the intensity of the boson peak shows a rather narrow sensitivity to changes in temperature and is seen to persist even at high temperatures. Study of the difference spectrum at two different polydispersity reveals that the increase in the height of boson peak is due to a population shift from modes with frequencies above the maximum in the VDOS to that below the maximum, indicating an increase in the fraction of the unstable modes in the system. The latter is further supported by the facilitation of the observed dynamics by polydispersity. Since the strength of the liquid increases with polydispersity, the present result provides an evidence that the intensity of boson peak correlates positively with the strength of the liquid, as observed earlier in many experimental systems.

preprint2006arXiv

Structure and Dynamics of DNA-dendrimer complexation: Role of counterions, water and base pair sequence

We study sequence dependent complexation between oligonucleotides (single strand DNA) and various generation ethylene diamine (EDA) cored poly amido amide (PAMAM) dendrimers through atomistic molecular dynamics simulations accompanied by free energy calculations and inherent structure determination. Simulations reveal formation of a stable complex and provide a detailed molecular level understanding of the structure and dynamics of such a complexation. The reaction free energy surface in the initial stage is found to be funnel-like with a significant barrier arising in the late stage due to the occurrence of misfolded states of DNA. Complexation shows surprisingly strong sensitivity to the ssDNA sequence which is found to arise from a competition between enthalpic versus entropic rigidity of ssDNA.

preprint2004arXiv

Rotational friction on small globular proteins: Combined dielectric and hydrodynamic effect

Rotational friction on proteins and macromolecules is known to derive contributions from at least two distinct sources -- hydrodynamic (due to viscosity) and dielectric friction (due to polar interactions). In the existing theoretical approaches, the effect of the latter is taken into account in an {\it ad hoc} manner, by increasing the size of the protein with the addition of a hydration layer. Here we calculate the rotational dielectric friction on a protein ($ζ_{DF}$) by using a generalized arbitrary charge distribution model (where the charges are obtained from quantum chemical calculation) and the hydrodynamic friction with stick boundary condition, ($ζ_{hyd}^{stick}$) by using the sophisticated theoretical technique known as tri-axial ellipsoidal method, formulated by Harding [S. E. Harding, Comp. Biol. Med. {\bf 12}, 75 (1982)]. The calculation of hydrodynamic friction is done with only the dry volume of the protein (no hydration layer). We find that the total friction obtained by summing up $ζ_{DF}$ and $ζ_{hyd}^{stick}$ gives reasonable agreement with the experimental results, i.e., $ζ_{exp} \approx ζ_{DF} + ζ_{hyd}^{stick}$.