Temperature, pressure, and length scale dependence of the solvation of simple solvophobic solutes is investigated in the Jagla liquid, a simple liquid consisting of particles that interact via a spherically symmetric potential combining hard and soft core interactions. The results are compared with identical calculations for a model of a typical atomic liquid, the Lennard-Jones (LJ) potential, and with predictions for hydrophobic solvation in water using the recently developed cavity equation of state and the extended simple point charge model. We find that the Jagla liquid captures the qualitative thermodynamic behavior of hydrophobic hydration as a function of temperature and pressure for both small and large length scale solutes. In particular, for both the Jagla liquid and water, we observe temperature-dependent enthalpy and entropy of solvation for all solute sizes as well as a negative solvation entropy for sufficiently small solutes at low temperature. This feature of water-like solvation is distinct from the strictly positive and temperature independent enthalpy and entropy of cavity solvation observed in the Lennard-Jones fluid. The results suggest that a competition between two energy scales that favors low-density, open structures as temperature is decreased is an essential interaction of a liquid that models hydrophobic hydration. In addition the Jagla liquid dewets surfaces of large radii of curvature less readily than the Lennard-Jones liquid, and the so-called ``length scale crossover'' in solvation, whereby solvation free energies change from scaling with the solute volume to scaling with the solute surface area, occurs at length scales that are larger relative to the solvent size. Both features reflect a greater flexibility or elasticity in the Jagla liquid structure than that of a typical liquid, similar to water's ability to maintain its hydrogen bond network. The implications of the differences in crossover behavior between water-like and typical liquids are examined in the context of a simple thought experiment on the aggregation of solvophobic solutes that builds on ideas from Chandler and Rajamani et al. We find that water-like crossover behavior exposes a size range of solvophobic aggregates to destabilization upon cooling and pressurizing, which may thereby precipitate phenomena such as cold and pressure denaturation of proteins. Statistics of density fluctuations, void space, and pair distributions are analyzed for molecular-scale volumes. The pair distribution functions are used to provide an estimate of the size of the Jagla particle with a physical basis. The void distributions are observed to be distinct in the three liquids, with low temperature distributions in the LJ and Jagla liquids demonstrating a high degree of skewness. The void distributions observed in LJ liquid are hard sphere-like, while those of water and the Jagla liquid exhibit a higher degree of density inhomogeneity relative to a hard sphere system. The well-known Gaussian behavior of density fluctuations in molecular volumes in water is not generally observed in other liquids, as evidenced by the fact that this behavior is not consistently observed in either the LJ or the Jagla liquids. An exploratory study of the effects of explicit solvent on the sequence energy landscape of model heteropolymers has been performed. For a fixed set of configurations, the energy landscape of all possible sequences taken from a two letter alphabet consisting of only solvophilic and solvophobic monomers is characterized at different solvent temperatures. Non-trivial solvent and temperature effects are manifest in the distribution of sequences, confirming that the negation of these effects may have profound consequences on designability. / text

Kyoto University (京都大学) / 0048 / 新制・課程博士 / 博士(エネルギー科学) / 甲第19092号 / エネ博第316号 / 新制||エネ||64 / 32043 / 京都大学大学院エネルギー科学研究科エネルギー基礎科学専攻 / (主査)教授 木下 正弘, 教授 森井 孝, 教授 片平 正人 / 学位規則第4条第1項該当

abstract: Molecular dynamics simulations were used to study properties of water at the interface with nanometer-size solutes. We simulated nonpolar attractive Kihara cavities given by a Lennard-Jones potential shifted by a core radius. The dipolar response of the hydration layer to a uniform electric field substantially exceeds that of the bulk. For strongly attractive solutes, the collective dynamics of the hydration layer become slow compared to bulk water, as the solute size is increased. The statistics of electric field fluctuations at the solute center are Gaussian and tend toward the dielectric continuum limit with increasing solute size. A dipolar probe placed at the center of the solute is sensitive neither to the polarity excess nor to the slowed dynamics of the hydration layer. A point dipole was introduced close to the solute-water interface to further study the statistics of electric field fluctuations generated by the water. For small dipole magnitudes, the free energy surface is single-welled, with approximately Gaussian statistics. When the dipole is increased, the free energy surface becomes double-welled, before landing in an excited state, characterized again by a single-welled surface. The intermediate region is fairly broad and is characterized by electrostatic fluctuations significantly in excess of the prediction of linear response. We simulated a solute having the geometry of C180 fullerene, with dipoles introduced on each carbon. For small dipole moments, the solvent response follows the results seen for a single dipole; but for larger dipole magnitudes, the fluctuations of the solute-solvent energy pass through a second maximum. The juxtaposition of the two transitions leads to an approximately cubic scaling of the chemical potential with the dipole strengh. Umbrella sampling techniques were used to generate free energy surfaces of the electric potential fluctuations at the heme iron in Cytochrome B562. The results were unfortunately inconclusive, as the ionic background was not effectively represented in the finite-size system. / Dissertation/Thesis / Ph.D. Chemistry 2012

京都大学 / 0048 / 新制・課程博士 / 博士(エネルギー科学) / 甲第19092号 / エネ博第316号 / 新制||エネ||64(附属図書館) / 32043 / 京都大学大学院エネルギー科学研究科エネルギー基礎科学専攻 / (主査)教授 木下 正弘, 教授 森井 孝, 教授 片平 正人 / 学位規則第4条第1項該当 / Doctor of Energy Science / Kyoto University / DFAM

Solvation-thermodynamics models based on computational quantum mechanics, such as the conductor-like screening model (COSMO), provide a good alternative to traditional group-contribution methods for predicting thermodynamic phase behavior. Two COSMO-based thermodynamic models are COSMO-RS (real solvents) and COSMO-SAC (segment activity coefficient). The main molecule-specific input for these models is the sigma profile, or the probability distribution of a molecular surface segment having a specific charge density. Generating the sigma profiles represents the most time-consuming and computationally expensive aspect of using COSMO-based methods. A growing number of scientists and engineers are interested in the COSMO-based thermodynamic models, but are intimidated by the complexity of generating the sigma profiles. This thesis presents the first free, open-literature database of 1,513 self-consistent sigma profiles, together with two validation examples. The offer of these profiles will enable interested scientists and engineers to use the quantum-mechanics-based, COSMO methods without having to do quantum mechanics. This thesis summarizes the application experiences reported up to October 2004 to guide the use of the COSMO-based methods. Finally, this thesis also provides a FORTRAN program and a procedure to generate additional sigma profiles consistent with those presented here, as well as a FORTRAN program to generate binary phase-equilibrium predictions using the COSMO-SAC model. / Master of Science