Source author record

Lawrence R. Pratt

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

ResearcherUnclaimed source record

physics.chem-ph Biological Physics cond-mat.soft cond-mat.stat-mech physics.atm-clus q-bio

Catalog footprint

What is connected

9works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2022arXiv

Hydrated anions: From clusters to bulk solution with quasi-chemical theory

The interactions of hydrated ions with solution and interface partners are strong on a chemical energy scale. Here, we test the foremost \textit{ab initio theory} for evaluation of hydration free energies of ions, namely, \textit{quasi-chemical theory} (QCT). We focus on halide anions, but also the hydroxide anion, since they have been outstanding challenges for all theories. QCT is built by identification of inner-shell clusters, separate treatment of those clusters, then integration of those results into the broader-scale solution environment. We exploit a close comparison with mass-spectrometric measurements of ion-hydration equilibria. That theory-experiment comparison is excellent with moderate computational effort here. This agreement reinforces both theory and experiment, and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. QCT provides a favorable setting for exploitation of the polarizable continuum model (PCM) when the inner-shell material shields the ion from the outer solution environment. For the asymmetrically hydrated, and less effectively shielded, heavier halide ions, we investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement and our final results are in close agreement with a standard tabulation of hydration free energies. Comparison of anion hydration cluster structures with bulk solutions from AIMD simulations emphasize slight differences: the asymmetries of bulk solution inner-shell structures are moderated, but still present; and inner shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

preprint2015arXiv

Interfaces of Propylene Carbonate

Propylene carbonate (PC) wets graphite with a contact angle of 31 deg at ambient conditions. Molecular dynamics simulations agree with this contact angle after 40% reduction of the strength of graphite-C atom Lennard-Jones interactions with the solvent, relative to the models used initially. A simulated nano-scale PC droplet on graphite displays a pronounced layering tendency and an Aztex pyramid structure for the droplet. Extrapolation of the computed tensions of PC liquid-vapor interface estimates the critical temperature of PC accurately to about 3%. PC molecules lie flat on the PC liquid-vapor surface, and tend to project the propyl carbon toward the vapor phase. For close PC neighbors in liquid PC, an important packing motif stacks carbonate planes with the outer oxygen of one molecule snuggled into the positively charged propyl end of another molecule so that neighboring molecule dipole moments are approximately antiparallel. The calculated thermal expansion coefficient and the dielectric constants for liquid PC agree well with experiment. The distribution of PC molecule binding energies is closely Gaussian. Evaluation of the density of the coexisting vapor then permits estimation of the packing contribution to the PC chemical potential, and that contribution is about 2/3rds of the magnitude of the contributions due to attractive interactions, with opposite sign.

preprint2012arXiv

Microstructures of capped ethylene oxide oligomers in water and n-hexane

This report documents microstructural features of CH3(CH2-O-CH2)mCH3 dissolved in water and n-hexane for m = 11, 21, and 31. Probability densities for end-to-end distance, and the associated potential-of-mean-force (pmf), are more revealing of chain microstructures than are the corresponding results for the radii of gyration. For water, the pmf identifies three distinct regions: loop-closure, globule, and high-extension regions. The globule region affirms a water-swollen chain, and is not evident in the n-hexane results. Chain C-atom density profiles from the chain centroid are also different in the water and n-hexane cases. For n-hexane (but not water), the density profiles are similar for the different chain lengths when the distances are scaled by the observed <Rg^2>^{1/2}. For water (but not n-hexane) and the smaller chains considered, the carbon material exhibits a distinctive enhanced concentration, or internal condensation, at the centroid core of the structure.

preprint2010arXiv

Direct observation of a hydrophobic bond in loop-closure of a capped (-OCH2CH2-)n oligomer in water

The small r variation of the probability density P(r) for end-to-end separations of a -CH2CH3 capped (-OCH2CH2-)n oligomer in water is computed to be closely similar to the CH4 ... CH4 potential of mean force under the same circumstances. Since the aqueous solution CH4 ... CH4 potential of mean force is the natural physical definition of a primitive hydrophobic bond, the present result identifies an experimentally accessible circumstance for direct observation of a hydrophobic bond which has not been observed previously because of the low solubility of CH4 in water. The physical picture is that the soluble chain molecule carries the capping groups into aqueous solution, and permits them to find one another with reasonable frequency. Comparison with the corresponding results without the solvent shows that hydration of the solute oxygen atoms swells the chain molecule globule. This supports the view that the chain molecule globule might have a secondary effect on the hydrophobic interaction which is of first interest here. The volume of the chain molecule globule is important for comparing the probabilities with and without solvent because it characterizes the local concentration of capping groups. Study of other capping groups to enable X-ray and neutron diffraction measurements of P(r) is discussed.

preprint2001arXiv

Molecular Theory of Hydrophobic Effects: ``She is too mean to have her name repeated.''

This paper reviews the molecular theory of hydrophobic effects relevant to biomolecular structure and assembly in aqueous solution. Recent progress has resulted in simple, validated molecular statistical thermodynamic theories and clarification of confusing theories of decades ago. Current work is resolving effects of wider variations of thermodynamic state, e.g. pressure denaturation of soluble proteins, and more exotic questions such as effects of surface chemistry in treating stability of macromolecular structures in aqueous solution

preprint1998arXiv

Molecular theories and simulation of ions and polar molecules in water

Recent developments in molecular theories and simulation of ions and polar molecules in water are reviewed. The hydration of imidazole and imidazolium solutes is used to exemplify the theoretical issues. The treatment of long-ranged electrostatic interactions in simulations is discussed extensively. It is argued that the Ewald approach is an easy way to get correct hydration free energies in the thermodynamic limit from molecular calculations; and that molecular simulations with Ewald interactions and periodic boundary conditions can also be more efficient than many common alternatives. The Ewald treatment permits a conclusive extrapolation to infinite system size. The picture that emerges from testing of simple models is that the most prominent failings of the simplest theories are associated with solvent proton conformations that lead to non-gaussian fluctuations of electrostatic potentials. Thus, the most favorable cases for the second-order perturbation theories are monoatomic positive ions. For polar and anionic solutes, continuum or gaussian theories are less accurate. The appreciation of the specific deficiencies of those simple models have led to new concepts, multistate gaussian and quasi-chemical theories, that address the cases for which the simpler theories fail. It is argued that, relative to direct dielectric continuum treatments, the quasi-chemical theories provide a better theoretical organization for the computational study of the electronic structure of solution species.

preprint1998arXiv

Quasi-chemical Theories of Associated Liquids

It is shown how traditional development of theories of fluids based upon the concept of physical clustering can be adapted to an alternative local clustering definition. The alternative definition can preserve a detailed valence description of the interactions between a solution species and its near-neighbors, i.e., cooperativity and saturation of coordination for strong association. These clusters remain finite even for condensed phases. The simplest theory to which these developments lead is analogous to quasi-chemical theories of cooperative phenomena. The present quasi-chemical theories require additional consideration of packing issues because they don't impose lattice discretizations on the continuous problem. These quasi-chemical theories do not require pair decomposable interaction potential energy models. Since calculations may be required only for moderately sized clusters, we suggest that these quasi-chemical theories could be implemented with computational tools of current electronic structure theory. This can avoid an intermediate step of approximate force field generation.

preprint1997arXiv

Hydrolysis of ferric ion in water and conformational equilibrium

Reported here are results of theoretical calculations on the hexaaquoferric complex and deprotonated products to investigate the molecular mechanisms of hydrolysis of ferric ion in water. The combination of density functional electronic structure techniques and a dielectric continuum model for electrostatic solvation applied to the Fe(H$_2$O)$_6$$^{3+}$ complex yields -1020~kcal/mol (experimental values -1037~kcal/mol to -1019~kcal/mol) for the absolute free energy of the aqueous ferric ion. The predicted reaction energy change for the first hydrolysis is 2 kcal/mol (compared to 3 kcal/mol experimental). For the second hydrolysis, we found an unexpected low energy isomer of Fe(H$_2$O)$_4$(OH)$_2$$^{+}$ with five ligands in the inner sphere and one water outside. The cis and trans isomers are, respectively, slightly lower and higher in energy. Extrusion of the additional water to the outer sphere is nearly thermoneutral. The reaction free energy for the second hydrolysis is predicted in the range 16-18~kcal/mol, higher than the experimental value of 5~kcal/mol. Since the theoretical predictions for the second hydrolysis are higher than experimental values and novel structures were encountered, we argue that conformational entropy is important in these reactions.

preprint1996arXiv

Origin of entropy convergence in hydrophobic hydration and protein folding

An information theory model is used to construct a molecular explanation why hydrophobic solvation entropies measured in calorimetry of protein unfolding converge at a common temperature. The entropy convergence follows from the weak temperature dependence of occupancy fluctuations for molecular-scale volumes in water. The macroscopic expression of the contrasting entropic behavior between water and common organic solvents is the relative temperature insensitivity of the water isothermal compressibility. The information theory model provides a quantitative description of small molecule hydration and predicts a negative entropy at convergence. Interpretations of entropic contributions to protein folding should account for this result.