Researcher profile

L. R. Pratt

L. R. Pratt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
3topics
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

10 published item(s)

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.

preprint2014arXiv

Concentration dependence of the Flory-Huggins interaction parameter in aqueous solutions of capped PEO chains

The dependence on volume fraction $φ$ of the Flory-Huggins $χ_{\mathrm{wp}}\left(φ\right)$ describing the free energy of mixing of polymers in water is obtained by exploiting the connection of $χ_{\mathrm{wp}}\left(φ\right)$ to the chemical potential of the water, for which quasi-chemical theory is satisfactory. We test this theoretical approach with simulation data for aqueous solutions of capped PEO oligomers. For CH$_3$(CH$_2$-O-CH$_2$)$_m$CH$_3$ ($m$=11), $χ_{\mathrm{wp}}\left(φ\right)$ depends strongly on $φ$, consistent with experiment. These results identify coexisting water-rich and water-poor solutions at $T$ = 300K and $p$ = 1atm. Direct observation of the coexistence of these two solutions on simulation time scales supports that prediction for the system studied. This approach directly provides the osmotic pressures. The osmotic second virial coefficient for these chains is positive, reflecting repulsive interactions between the chains in the water, a good solvent for these chains.

preprint2014arXiv

Hydration of Kr(aq) in dilute and concentrated solutions

Molecular dynamics simulations of water with both multi-Kr and single Kr atomic solutes are carried out to implement quasi-chemical theory evaluation of the hydration free energy of Kr(aq). This approach obtains free energy differences reflecting Kr-Kr interactions at higher concentrations. Those differences are negative changes in hydration free energies with increasing concentrations at constant pressure. The changes are due to a slight reduction of packing contributions in the higher concentration case. The observed Kr-Kr distributions, analyzed with the extrapolation procedure of Krüger, \emph{et al.}, yield a modestly attractive osmotic second virial coefficient, $B_2\approx -60~\mathrm{cm}^3$/mol. The thermodynamic analysis interconnecting these two approaches shows that they are closely consistent with each other, providing support for both.

preprint2014arXiv

Loop-Closure and Gaussian Models of Collective Structural Characteristics of Capped PEO Oligomers in Water

Parallel-tempering MD results for a CH$_3$(CH$_2$-O-CH$_2$)$_m$CH$_3$ chain in water are exploited as a data-base for analysis of collective structural characteristics of the PEO globule with a goal of defining models permitting statistical thermodynamic analysis of dispersants of Corexit type. The chain structure factor, relevant to neutron scattering from a deuterated chain in neutral water, is considered specifically. The traditional continuum-Gaussian structure factor is inconsistent with the simple $k \rightarrow \infty$ behavior, but we consider a discrete-Gaussian model that does achieve that consistency. Shifting-and-scaling the discrete-Gaussian model helps to identify the low-$k$ to high-$k$ transition near $k \approx 2π/0.6 \mathrm{nm}$ when an empirically matched number of Gaussian links is about one-third of the total number of effective-atom sites. This short distance-scale boundary of 0.6 nm is directly verified with the $r$-space distributions, and this distance is thus identified with a natural size for coarsened monomers. The probability distribution of $R_g{}^2$ is compared with the classic predictions for both Gaussian model and freely-jointed chains. $\left\langle R_g{}^2(j)\right\rangle$, the contribution of the $j$-th chain segment to $\left\langle R_g{}^2\right\rangle$, depends on contour index about as expected for Gaussian chains despite significant quantitative discrepancies which express the swelling of these chains in water. Monomers central to the chain contour occupy the center of the chain globule. The density profiles of chain segments relative to their center of mass can show distinctive density structuring for smaller chains due close proximity of central elements to the globule center. But that density structuring washes-out for longer chains where many chain elements additively contribute to the density profiles.

preprint2014arXiv

The role of attractive interactions in the dynamics of molecules in liquids

The friction kernel (or memory function) $γ(t)$ characterizing single-molecule dynamics in strongly bound liquids exhibits two distinct relaxations with the longer time-scale relaxation associated with attractive intermolecular forces. This observation identifies differing roles of repulsive and attractive interaction in the motions of molecules in equilibrium liquids, and thus provides a basis for a renewed investigation of a van der Waals picture of the transport properties of liquids. This conclusion is supported by extracting $γ(t)$ from molecular dynamics simulation data for four common molecular liquids.

preprint2013arXiv

Case study of Rb+(aq), quasi-chemical theory of ion hydration, and the no split occupancies rule

Quasi-chemical theory applied to ion hydration combines statistical mechanical theory, electronic structure calculations, and molecular simulation, disciplines which are individually subjects for specialized professional attention. Because it combines activities which are themselves non-trivial, quasi-chemical theory is typically viewed with surprise. Nevertheless, it provides a fully-considered framework for analysis of ion hydration. Furthermore, the initial calculations are indeed simple, successful, and provide new information to long-standing experimental activities such as neutron diffraction by hydrated ions. Here we review quasi-chemical theory in the context of a challenging application, Rb+(aq).

preprint2013arXiv

Direct Numerical Test of the Statistical Mechanical Theory of Hydrophobic Interactions

This work tests the statistical mechanical theory of hydrophobic interactions, isolates consequences of excluded volume interactions, and obtains B2 for those purposes. Cavity methods that are particularly appropriate for study of hydrophobic interactions between atomic-size hard spheres in liquid water are developed and applied to test aspects of the Pratt-Chandler (PC) theory that have not been tested. Contact hydrophobic interactions between Ar-size hard-spheres in water are significantly more attractive than predicted by the PC theory. The corresponding results for the osmotic second virial coefficient are attractive (B2 <0), and more attractive with increasing temperature (Delta B2/Delta T < 0) in the temperature range 300K < T < 360K. This information has not been available previously, but is essential for development of the molecular-scale statistical mechanical theory of hydrophobic interactions, particularly for better definition of the role of attractive intermolecular interactions associated with the solutes.

preprint2013arXiv

Multi-scale theory in the molecular simulation of electrolyte solutions

This paper organizes McMillan-Mayer theory, the potential distribution approach, and quasi-chemical theory to provide theory for the thermodynamic effects associated with longer spatial scales involving longer time scales, thus helping to define a role for AIMD simulation directly on the time and space scales typical of those demanding methods. The theory treats composition fluctuations which would be accessed by larger-scale calculations, and also longer-ranged interactions that are of special interest for electrolyte solutions. The quasi-chemical organization breaks-up governing free energies into physically distinct contributions: packing, outer-shell, and chemical contributions. Here we study specifically the outer-shell contributions that express electrolyte screening. For that purpose we adopt a primitive model suggested by observation of ion-pairing in tetra-ethylammonium tetra-fluoroborate dissolved in propylene carbonate. Gaussian statistical models are shown to be effective physical models for outer-shell contributions, and they are conclusive for the free energies within the quasi-chemical formulation. With the present data-set the gaussian physical approximation obtains more accurate mean activity coefficients than does the Bennett direct evaluation of that free energy.

preprint2012arXiv

Pairing of 1-hexyl-3-methylimidazolium and tetrafluoroborate ions in n-pentanol

Molecular dynamics simulations are obtained and analyzed to study pairing of 1-hexyl-3-methylimidazolium and tetrafluoroborate ions in n-pentanol, in particular by evaluating the potential-of-mean-force between counter ions. The present molecular model and simulation accurately predicts the dissociation constant Kd in comparison to experiment, and thus the behavior and magnitudes for the ion-pair pmf at molecular distances, even though the dielectric constant of the simulated solvent differs from the experimental value by about 30%. A naive dielectric model does not capture molecule structural effects such as multiple conformations and binding geometries of the Hmim+ and BF4- ion-pairs. Mobilities identify multiple time-scale effects in the autocorrelation of the random forces on the ions, and specifically a slow, exponential time-decay of those long-ranged forces associated here with dielectric friction effects.

preprint2011arXiv

Generalizations of the Fuoss Approximation for Ion Pairing

An elementary statistical observation identifies generalizations of the Fuoss approximation for the probability distribution function that describes ion clustering in electrolyte solutions. The simplest generalization, equivalent to a Poisson distribution model for inner-shell occupancy, exploits measurable inter-ionic correlation functions, and is correct at the closest pair distances whether primitive electrolyte solutions models or molecularly detailed models are considered, and for low electrolyte concentrations in all cases. With detailed models these generalizations includes non-ionic interactions and solvation effects. These generalizations are relevant for computational analysis of bi-molecular reactive processes in solution. Comparisons with direct numerical simulation results show that the simplest generalization is accurate for a slightly supersaturated solution of tetraethylammonium tetrafluoroborate in propylene carbonate ([tea][BF$_4$]/PC), and also for a primitive model associated with the [tea][BF$_4$]/PC results. For [tea][BF$_4$]/PC, the atomically detailed results identify solvent-separated nearest-neighbor ion-pairs. This generalization is examined also for the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF$_4$]) where the simplest implementation is less accurate. In this more challenging situation an augmented maximum entropy procedure is satisfactory, and explains the more varied near-neighbor distributions observed in that case.