BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

D. Asthagiri

D. Asthagiri contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.chem-phBiological Physicscond-mat.softBiomoleculescond-mat.mtrl-sci

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2020arXiv

$^1$H-NMR Dipole-Dipole Relaxation in Fluids: Relaxation of Individual $^1$H-$^1$H Pairs versus Relaxation of Molecular Modes

The intra-molecular $^1$H-NMR dipole-dipole relaxation of molecular fluids has traditionally been interpreted within the Bloembergen-Purcell-Pound (BPP) theory of NMR intra-molecular relaxation. The BPP theory draws upon Debye's theory for describing the rotational diffusion of the $^1$H-$^1$H pair and predicts a mono-exponential decay of the $^1$H-$^1$H dipole-dipole autocorrelation function between distinct spin pairs. Using molecular dynamics (MD) simulations, we show that for both $n$-heptane and water this is not the case. In particular, the autocorrelation function of individual $^1$H-$^1$H intra-molecular pairs itself evinces a rich stretched-exponential behavior, implying a distribution in rotational correlation times. However for the high-symmetry molecule neopentane, the individual $^1$H-$^1$H intra-molecular pairs do conform to the BPP description, suggesting an important role of molecular symmetry in aiding agreement with the BPP model. The inter-molecular autocorrelation functions for $n$-heptane, water, and neopentane also do not admit a mono-exponential behavior of individual $^1$H-$^1$H inter-molecular pairs at distinct initial separations. We suggest expanding the auto-correlation function in terms of molecular modes, where the molecular modes do have an exponential relaxation behavior. With care, the resulting Fredholm integral equation of the first kind can be inverted to recover the probability distribution of the molecular modes. The advantages and limitations of this approach are noted.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Molecular Dynamics Simulations of NMR Relaxation and Diffusion of Heptane Confined in a Polymer Matrix

The mechanism behind the NMR surface relaxation and the large $T_1$/$T_2$ ratio of light hydrocarbons confined in the nano-pores of kerogen remains poorly understood, and consequently has engendered much debate. Towards bringing a molecular-scale resolution to this problem, we present molecular dynamics (MD) simulations of $^1$H NMR relaxation and diffusion of heptane in a polymer matrix, where the high-viscosity polymer is a model for kerogen and bitumen that provides an organic "surface" for heptane. We calculate the autocorrelation function $G(t)$ for $^1$H-$^1$H dipole-dipole interactions of heptane in the polymer matrix and use this to generate the NMR frequency ($f_0$) dependence of $T_1$ and $T_2$ relaxation times as a function of $ϕ_{C7}$. We find that increasing molecular confinement increases the correlation time of the heptane molecule, which decreases the surface relaxation times for heptane in the polymer matrix. For weak confinement ($ϕ_{C7} > 50$ vol%), we find that $T_{1S}/T_{2S} \simeq 1$. Under strong confinement ($ϕ_{C7} \lesssim $ 50 vol%), we find that the ratio $T_{1S}/T_{2S} \gtrsim 4$ increases with decreasing $ϕ_{C7}$, and that the dispersion relation $T_{1S} \propto f_0$ is consistent with previously reported measurements of polymers and bitumen. Such frequency dependence in bitumen has been previously attributed to paramagnetism, but our studies suggests that $^1$H-$^1$H dipole-dipole interactions enhanced by organic nano-pore confinement dominates the NMR response in saturated organic-rich shales, without the need to invoke paramagnetism.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Incorporating multi-body effects in SAFT by improving the description of the reference system. I. Mean activity correction for cluster integrals in the reference system

A system of patchy colloidal particles interacting with a solute that can associate multiple times in any direction is a useful model for patchy colloidal mixtures. Despite the simplicity of the interaction, because of the presence of multi-body correlations predicting the thermodynamics of such systems remains a challenge. Earlier Marshall and Chapman developed a multi-body formulation for such systems wherein the cluster partition function for the hard-sphere solvent molecules in a defined inner-shell (or coordination volume) of the hard-sphere solute is used as the reference within the statistical association fluid theory formalism. The multi-body contribution to these partition functions are obtained by ignoring the bulk solvent, thus limiting the applicability of the theory to low system densities. Deriving inspiration from the quasichemical theory of solutions where these partition functions occur in the guise of equilibrium constants for cluster formation, we develop a way to account for the multi-body correlations including the effect of the bulk solvent. We obtain the free energy to evacuate the inner-shell, the chemistry contribution within quasichemical theory, from simulations of the hard-sphere reference. This chemistry contribution reflects association in the reference in the presence of the bulk medium. The gas-phase partition functions are then augmented by a mean activity factor that is adjusted to reproduce the chemistry contribution. We show that the updated partition function provides a revised reference that better captures the distribution of solvent around the solute up to high system densities. Using this updated reference, we find that theory better describes both the bonding state and the excess chemical potential of the colloid in the physical system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Structure and thermodynamics of a mixture of patchy and spherical colloids: a multi-body association theory with complete reference fluid information

A mixture of solvent particles with short-range, directional interactions and solute particles with short-range, isotropic interactions that can bond multiple times is of fundamental interest in understanding liquids and colloidal mixtures. Because of multi-body correlations predicting the structure and thermodynamics of such systems remains a challenge. Earlier Marshall and Chapman developed a theory wherein association effects due to interactions multiply the partition function for clustering of particles in a reference hard-sphere system. The multi-body effects are incorporated in the clustering process, which in their work was obtained in the absence of the bulk medium. The bulk solvent effects were then modeled approximately within a second order perturbation approach. However, their approach is inadequate at high densities and for large association strengths. Based on the idea that the clustering of solvent in a defined coordination volume around the solute is related to occupancy statistics in that defined coordination volume, we develop an approach to incorporate the complete information about hard-sphere clustering in a bulk solvent at the density of interest. The occupancy probabilities are obtained from enhanced sampling simulations but we also develop a concise parametric form to model these probabilities using the quasichemical theory of solutions. We show that incorporating the complete reference information results in an approach that can predict the bonding state and thermodynamics of the colloidal solute for a wide range of system conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermodynamics of the collapse transition of the all-backbone peptide Gly15

Simulations show Gly$_{15}$, a polypeptide lacking any side-chains, can collapse in water. We assess the hydration thermodynamics in this collapse by calculating the hydration free energy at each of the end points of the reaction coordinate, here the end-to-end distance ($r$) in the chain. To examine the role of the various conformations for a given $r$, we study the conditional distribution, $P(R_g | r)$, of the radius of gyration for a given value of $r$. $P(R_g|r)$ is found to vary more gently compared to the corresponding variation in the excess hydration free energy. Using this insight within a multistate generalization of the potential distribution theorem, we calculate a reasonable upper bound for the hydration free energy of the peptide for a given $r$. On this basis we find that peptide hydration greatly favors the expanded state of the chain, despite primitive hydrophobic effects favoring chain collapse. The net free energy of collapse is seen to be a delicate balance between opposing intra-peptide and hydration effects, with intra-peptide contributions favoring collapse by a small margin. The favorable intra-peptide interactions are primarily electrostatic in origin, and found to arise primarily from interaction between C=O dipoles, hydrogen bonding interaction between C=O and N-H groups, and favorable interaction between N-H dipoles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Molecular theory and the effects of solute attractive forces on hydrophobic interactions

The role of solute attractive forces on hydrophobic interactions is studied by coordinated development of theory and simulation results for Ar atoms in water. We present a concise derivation of the local molecular field (LMF) theory for the effects of solute attractive forces on hydrophobic interactions, a derivation that clarifies the close relation of LMF theory to the EXP approximation applied to this problem long ago. The simulation results show that change from purely repulsive atomic solute interactions to include realistic attractive interactions \emph{diminishes} the strength of hydrophobic bonds. For the Ar-Ar rdfs considered pointwise, the numerical results for the effects of solute attractive forces on hydrophobic interactions are of opposite sign and larger in magnitude than predicted by LMF theory. That comparison is discussed from the point of view of quasi-chemical theory, and it is suggested that the first reason for this difference is the incomplete evaluation within LMF theory of the hydration energy of the Ar pair. With a recent suggestion for the system-size extrapolation of the required correlation function integrals, the Ar-Ar rdfs permit evaluation of osmotic second virial coefficients $B_2$. Those $B_2$ also show that incorporation of attractive interactions leads to more positive (repulsive) values. With attractive interactions in play, $B_2$ can change from positive to negative values with increasing temperatures. This is consistent with the historical work of Watanabe, \emph{et al.,} that $B_2 \approx 0$ for intermediate cases. In all cases here, $B_2$ becomes more attractive with increasing temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Role of hydration and intramolecular interactions in the helix-coil transition and helix-helix assembly in a deca-alanine peptide

For a model deca-alanine peptide the cavity (ideal hydrophobic) contribution to hydration favors the helix state in the coil-to-helix transition and the paired helix bundle in the assembly of two helices. The energetic contributions of attractive protein-solvent interactions are separated into a short-range part arising from interactions with solvent in the first hydration shell and the remaining long-range part. In the helix-coil transition, short-range attractive protein-solvent interactions outweigh hydrophobic hydration and favor the unfolded coil states. Analysis of enthalpic effects shows that it is the favorable hydration of the peptide backbone that favors the unfolded state. Protein intramolecular interactions favor the helix state and are decisive in folding. In the pairing of two helices, the cavity contribution outweighs short-range attractive protein-water interactions. However, long-range, protein-solvent attractive interactions can either enhance or reverse this trend depending on the mutual orientation of the helices. In helix-helix assembly, change in enthalpy arising from change in attractive protein-solvent interactions favors disassembly. In helix pairing as well, favorable protein intramolecular interactions are found to be as important as hydration effects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conditional solvation of isoleucine in model extended and helical peptides: context dependence of hydrophobic hydration and the failure of the group-transfer model

The hydration thermodynamics of the GXG tripeptide relative to the reference GGG is often used to define the conditional hydration contribution of X. This quantity or the hydration thermodynamics of a small molecule analog of the side-chain or some combination of such estimates, have anchored the interpretation of seminal experiments on protein stability and folding. We show that such procedures to model protein hydration have significant limitations. We study the conditional hydration thermodynamics of the isoleucine side-chain in an extended pentapeptide and in helical deca-peptides, using as appropriate an extended penta-glycine or appropriate helical deca-peptides as reference. Hydration of butane in the gauche conformation provides a small molecule reference for the side-chain. We use the quasichemical theory to parse the hydration thermodynamics into chemical, packing, and long-range interaction contributions. The chemical contribution reflects the contribution of solvent clustering within the defined inner-shell of the solute; the chemical contribution of g-butane is substantially more negative than the conditional chemical contribution of isoleucine. The packing contribution gives the work required to create a cavity in the solvent, a quantity of interest in understanding hydrophobic hydration. The packing contribution for g-butane substantially overestimates the conditional packing of isoleucine. The net of such compensating contributions still disagrees with the conditional free energy of isoleucine but by a lesser magnitude. The excess enthalpy and entropy of hydration of g-butane model are also more negative than the corresponding conditional quantities for the side-chain. The conditional solvation of isoleucine in GGIGG also proves unsatisfactory in describing the conditional solvation of isoleucine in the helical peptides.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Examining the meaning of the peptide transfer free energy obtained from blocked (Gly)_n and cyclic-diglycine model compounds

In experiments, the free energy of transferring the peptide group from water to an osmolyte solution is obtained using the transfer free energy of (Gly)_n with the added assumption that a constant incremental change in free energy with n implies that each additional unit makes an independent contribution to the free energy. Here we test this assumption and uncover its limitations. Together with results for cyclic-diglycine, we show that, in principle, it is not possible to obtain a peptide group transfer free energy that is independent of the model system. We calculate the hydration free energy of acetyl-(Gly)_n-methyl amide (n=1..7) peptides modeled in the extended conformation in water and osmolyte solutions and find that the hydration free energy is linear in n, suggestive of independent, additive group-contributions. To probe the observed linearity further, we study the hydration of the solute bereft of water molecules in the first hydration shell. This conditioned solute arises naturally in the theoretical formulation and helps us focus on hydration effects uncluttered by the complexities of short-range solute-water interactions. We subdivide the conditioned solute into n+1 peptide groups and a methyl end group. The binding energy of each of these groups with the solvent is Gaussian distributed, but the near neighbor binding energies are themselves correlated: the (i,i+1) correlation is the strongest and tends to lower the free energy over the independent group case. We show that the observed linearity can be explained by the similarity of near neighbor correlations. Implications for group additive transfer free energy models are indicated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ion-water clusters, bulk medium effects, and ion hydration

Thermochemistry of gas-phase ion-water clusters together with estimates of the hydration free energy of the clusters and the water ligands are used to calculate the hydration free energy of the ion. Often the hydration calculations use a continuum model of the solvent. The primitive quasichemical approximation to the quasichemical theory provides a transparent framework to anchor such efforts. Here we evaluate the approximations inherent in the primitive quasichemical approach and elucidate the different roles of the bulk medium. We find that the bulk medium can stabilize configurations of the cluster that are usually not observed in the gas phase, while also simultaneously lowering the excess chemical potential of the ion. This effect is more pronounced for soft ions. Since the coordination number that minimizes the excess chemical potential of the ion is identified as the optimal or most probable coordination number, for such soft ions, the optimum cluster size and the hydration thermodynamics obtained without account of the bulk medium on the ion-water clustering reaction can be different from those observed in simulations of the aqueous ion. The ideas presented in this work are expected to be relevant to experimental studies that translate thermochemistry of ion-water clusters to the thermodynamics of the hydrated ion and to evolving theoretical approaches that combine high-level calculations on clusters with coarse-grained models of the medium.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Range of ion specific effects in the hydration of ions

Within the quasichemical approach, the hydration free energy of an ion is decomposed into a chemical term accounting for ion specific ion-water interactions within the coordination sphere and nonspecific contributions accounting for packing (excluded volume) and long range interactions. The change in the chemical term with a change in the radius of the coordination sphere is the compressive force exerted by the bulk solvent medium on the surface of the coordination sphere. For the Na+, K+, F-, and Cl- ions considered here this compressive force becomes equal for similarly charged ions for coordination radii of about 0.39 nm, not much larger than a water molecule. These results show that ion specific effects are short ranged and arise primarily due to differences in the local ion-water interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Regularizing binding energy distributions and thermodynamics of hydration. Application to water modeled with classical and ab initio simulations

The high-energy tail of the distribution of solute-solvent interaction energies is poorly characterized for condensed systems, but this tail region is of principal interest in determining the excess free energy of the solute. We introduce external fields centered on the solute to modulate the short-range repulsive interaction between the solute and solvent. This regularizes the binding energy distribution and makes it easy to calculate the free energy of the solute with the field. Together with the work done to apply the field in the presence and absence of the solute, we calculate the excess chemical potential of the solute. We present the formal development of this idea and apply it to study liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Water coordination structures and the excess free energy of the liquid

For a distinguished water molecule, the solute water, we assess the contribution of each coordination state to its excess chemical potential, using a molecular aufbau approach. In this approach, we define a coordination sphere, the inner-shell, and separate the excess chemical potential into packing, outer-shell, and local chemical contributions; the coordination state is defined by the number of solvent water molecules within the coordination sphere. The packing term accounts for the free energy of creating a solute-free coordination sphere in the liquid. The outer-shell term accounts for the interaction of the solute with the fluid outside the coordination sphere and it is accurately described by a Gaussian model of hydration for coordination radii greater than the minimum of the oxygen-oxygen pair correlation function. Consistent with the conventional radial cut-off used for defining hydrogen-bonds in liquid water, theory helps identify a chemically meaningful coordination radius. The local chemical contribution is recast as a sum over coordination states. The n-th term in this sum is given by the probability of observing n water molecules inside the coordination sphere in the absence of the solute water times a factor accounting for the interaction of the solute with the inner-shell solvent water molecules. Using this molecular aufbau expansion, we monitor the change in the chemical contribution due to the incremental increase in n. We find that though four water molecules are needed to fully account for the chemical term, the first added water accounts for nearly half the chemical term. Our results emphasize the need to acknowledge the intrinsic occupancy of a solute-free coordination sphere together with solute-solvent interactions in rationalizing the tetrahedral coordination of the solute water.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Molecular packing and chemical association in liquid water simulated using ab initio hybrid Monte Carlo and different exchange-correlation functionals

In the free energy of hydration of a solute, the chemical contribution is given by the free energy required to expel water molecules from the coordination sphere and the packing contribution is given by the free energy required to create the solute-free coordination sphere (the observation volume) in bulk water. With the SPC/E water model as a reference, we examine the chemical and packing contributions in the free energy of water simulated using different electron density functionals. The density is fixed at a value corresponding to that for SPC/E water at a pressure of 1 bar. The chemical contribution shows that water simulated at 300 K with BLYP is somewhat more tightly bound than water simulated at 300 K with the revPBE functional or at 350 K with the BLYP and BLYP-D functionals. The packing contribution for various radii of the observation volume is studied. In the size range where the distribution of water molecules in the observation volume is expected to be Gaussian, the packing contribution is expected to scale with the volume of the observation sphere. Water simulated at 300 K with the revPBE and at 350 K with BLYP-D or BLYP conforms to this expectation, but the results suggest an earlier onset of system size effects in the BLYP 350 K and revPBE 300 K systems than that observed for either BLYP-D 350 K or SPC/E. The implication of this observation for constant pressure simulations is indicated. For water simulated at 300 K with BLYP, in the size range where Gaussian distribution of occupation is expected, we instead find non-Gaussian behavior, and the packing contribution scales with surface area of the observation volume, suggesting the presence of heterogeneities in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Thermodynamics of water modeled using ab initio simulations

We regularize the potential distribution framework to calculate the excess free energy of liquid water simulated with the BLYP-D density functional. The calculated free energy is in fair agreement with experiments but the excess internal energy and hence also the excess entropy are not. Our work emphasizes the importance of thermodynamic characterization in assessing the quality of electron density functionals in describing liquid water and hydration phenomena.

0 citations0 reviewsNo code linkedOpen access