1
|
Ogrin P, Urbic T. Thermodynamics perturbation theory for solvation of nonpolar solutes in rose model. Phys Rev E 2023; 108:054135. [PMID: 38115497 DOI: 10.1103/physreve.108.054135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Accepted: 11/07/2023] [Indexed: 12/21/2023]
Abstract
A simple model of water, called the rose model, is used in this work. The rose model is a very simple model that can provide insight into the anomalous properties of water. In the rose water model, the molecules are represented as two-dimensional Lennard-Jones disks with potentials for orientation-dependent pairwise interactions mimicking formations of hydrogen bonds. We have recently applied a Wertheim integral equation theory (IET) and a thermodynamic perturbation theory (TPT) to the rose model in bulk. These analytical theories offer the advantage of being computationally less intensive than computer simulations by orders of magnitudes. Here we have applied the TPT to study the transfer of a nonpolar solute into rose water, the so-called hydrophobic effect. Similarly as in our previous work for bulk water, we have found that the theory reproduces the computer simulation results quite well at higher temperatures, while the theories predict the qualitative trends at low temperatures.
Collapse
Affiliation(s)
- Peter Ogrin
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
2
|
Ogrin P, Urbic T. Angle-dependent integral equation theory improves results of thermodynamics and structure of rose water model. J Chem Phys 2023; 159:114505. [PMID: 37732557 PMCID: PMC10908565 DOI: 10.1063/5.0159438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Accepted: 08/18/2023] [Indexed: 09/22/2023] Open
Abstract
Orientation-dependent integral equation theory (ODIET) was applied to the rose water model. Structural and thermodynamic properties of water modeled with the rose model were calculated using ODIET and compared to results from orientation-averaged integral equation theory (IET) and Monte Carlo simulations. Rose water model is a simple two-dimensional water model where molecules of water are represented as Lennard-Jones disks with explicit hydrogen bonding potential in form of rose functions. Orientational dependency significantly improves IET, as the thermodynamic results obtained using ODIET are significantly more in agreement with results calculated using MC than in the case of the orientationally averaged version. At high temperatures, the agreement between the simulation and theory is quantitative; however, when temperatures lower, a slight deviation between results obtained with different methods appear. ODIET correctly predicts the radial distribution function; moreover, ODIet also enables the calculation of angular distributions. While the angular distributions obtained with ODIET are in qualitative agreement with distributions from MC simulations, the height of the peaks in angular distributions differs between methods. Using results from ODIET, the spatial distribution of water molecules was constructed, which aids in the interpretation of other structural properties. ODIET was also used to calculate fractions of molecules with different number of hydrogen bonds, which is in the agreement with the simulations. Overall, use of ODIET significantly improves the obtained results in comparison to standard IET.
Collapse
Affiliation(s)
- Peter Ogrin
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
3
|
Urbic T. Ben Naim's four-arm model with density anomaly: Theory and computer simulations. Phys Rev E 2023; 108:014136. [PMID: 37583205 DOI: 10.1103/physreve.108.014136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Accepted: 07/01/2023] [Indexed: 08/17/2023]
Abstract
Molecular dynamics, Wertheim's integral equation theory (IET), and thermodynamics perturbation theory (TPT) were used to study the thermodynamics and structure of particles interacting through angle-dependent potential. The particles are modeled as two-dimensional Lennard-Jones disks with four hydrogen bonding arms arranged symmetrically. The model was introduced by Ben-Naim and we call it the BN4 model. The BN4 model exhibits density anomaly and other anomalous properties similar to those in water and in the Mercedes-Benz (MB) model. The IET is based on the orientationally averaged version of the Ornstein-Zernike equation and correctly predicts the pair correlation function of the model at high temperatures. Both TPT and IET are in semiquantitative agreement with the simulation values of the molar volume, isothermal compressibility, thermal expansion coefficient, and heat capacity.
Collapse
Affiliation(s)
- Tomaz Urbic
- University of Ljubljana, Faculty of Chemistry and Chemical Technology, Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
4
|
Ogrin P, Urbic T. Simple rose model of water in constant electric field. Phys Rev E 2023; 107:054801. [PMID: 37329104 DOI: 10.1103/physreve.107.054801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 04/12/2023] [Indexed: 06/18/2023]
Abstract
A simple two-dimensional statistical mechanical water model, called the rose model, was used in this work. We studied how a homogeneous constant electric field affects the properties of water. The rose model is a very simple model that helps explain the anomalous properties of water. Rose water molecules are represented as two-dimensional Lennard-Jones disks with potentials for orientation-dependent pairwise interactions mimicking formations of hydrogen bonds. The original model is modified by addition of charges for interaction with the electric field. We studied what kind of influence the electric field strength has on the model's properties. To determine the structure and thermodynamics of the rose model under the influence of the electric field we used Monte Carlo simulations. Under the influence of a weak electric field the anomalous properties and phase transitions of the water do not change. On the other hand, the strong fields shift the phase transition points as well as the position of the density maximum.
Collapse
Affiliation(s)
- Peter Ogrin
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
5
|
Urbic T. The electric field changes the anomalous properties of the Mercedes Benz water model. Phys Chem Chem Phys 2023; 25:4987-4996. [PMID: 36722865 PMCID: PMC9906975 DOI: 10.1039/d2cp05670d] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The influence of a homogeneous constant electric field on water properties was assessed. We used a simple two-dimensional statistical mechanical model called the Mercedes-Benz (MB) model of water in the study. The MB water molecules are two-dimensional disks with Gaussian arms that mimic the formation of hydrogen bonds. The model is modified with added charges for interaction with the electric field. The influence of the strength of the electric field on the water's properties was studied using Monte Carlo simulations. The structure and thermodynamics of the water were determined as a function of the strength of the electric field. We observed that the properties and phase transitions of the water in the low strength electric field does not change. In contrast, the high strength electric field shifts boiling and melting points as well as the position of the density maxima. After further increasing the strength of the electric field the density anomaly disappears.
Collapse
Affiliation(s)
- Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Askerceva 5, SI-1000, Slovenia.
| |
Collapse
|
6
|
Pršlja P, Žibert T, Urbic T. Monte Carlo simulations of simple two dimensional water-alcohol mixtures. J Mol Liq 2022; 368:120692. [PMID: 37731590 PMCID: PMC10508878 DOI: 10.1016/j.molliq.2022.120692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Simple alcohols such as methanol and ethanol, are organic chemicals that can be used to store energy, which can be used as an alternative to fossil fuels. Each alcohol has at least one hydroxyl group attached to a carbon atom of an alkyl group. They can be considered as organic derivatives of water in which one of the hydrogen atoms is replaced by an alkyl group. In this work, we determined the thermodynamic and structural properties of two dimensional water-alcohol mixtures using the Monte Carlo method. We used two-dimensional Mercedes-Benz (MB) model for water and MB based models for lower alcohols. The structural and thermodynamic properties of the mixtures were studied by Monte Carlo simulations in the isothermal-isobaric ensemble. We show that 2D models display similar trends in the density maxima as in real water-alcohol mixtures. With increasing content of alcohols, the temperature of maxima increases and upon further increase starts to decrease and at high concentrations, the density maxima disappears.
Collapse
Affiliation(s)
- Paulina Pršlja
- Department of Electrical Engineering and Automation, School of Electrical Engineering, Aalto University, 02150 Espoo, Finland
| | - Taja Žibert
- Department of Physical Chemistry, Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna pot 113, SI-1000, Slovenia
| | - Tomaz Urbic
- Department of Physical Chemistry, Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna pot 113, SI-1000, Slovenia
| |
Collapse
|
7
|
Ogrin P, Urbic T. Liquid-vapour coexistence line and percolation line of rose water model. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120531] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
8
|
Ogrin P, Urbic T, Fennell CJ. Statistical-mechanical liquid theories reproduce anomalous thermodynamic properties of explicit two-dimensional water models. Phys Rev E 2022; 106:034115. [PMID: 36266898 PMCID: PMC10061499 DOI: 10.1103/physreve.106.034115] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2022] [Accepted: 08/12/2022] [Indexed: 11/07/2022]
Abstract
We have developed an analytical theory for a simple model of liquid water. We apply Wertheim's thermodynamic perturbation theory (TPT) and integral equation theory (IET) for associative liquids to the rose model, which is among the simplest models of water. The particles interact through rose potentials for orientation dependent pairwise interactions. Modifying both the shape and range of a three-petal rose function, we construct an efficient and dynamical mimic of the two-dimensional (2D) Mercedes-Benz (MB) water model. The particles in 2D MB are 2D Lennard-Jones disks with three hydrogen bonding arms arranged symmetrically, resembling the Mercedes-Benz logo. Both models qualitatively predict both the anomalous properties of pure water and the anomalous solvation thermodynamics of nonpolar molecules. The IET is based on the orientationally averaged version of the Ornstein-Zernike equation. This is one of the main approximations in the present work. IET correctly predicts the pair correlation functions at high temperatures. Both TPT and IET are in semi-quantitative agreement with the Monte Carlo values of the molar volume, isothermal compressibility, thermal expansion coefficient, and heat capacity. A major advantage of these theories is that they require orders of magnitude less computer time than the Monte Carlo simulations.
Collapse
Affiliation(s)
- Peter Ogrin
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Christopher J. Fennell
- Department of Chemistry, Oklahoma State University, Stillwater, Oklahoma 74078, United States
| |
Collapse
|
9
|
Papez P, Urbic T. Simple two-dimensional models of alcohols. Phys Rev E 2022; 105:054608. [PMID: 35706252 PMCID: PMC10040488 DOI: 10.1103/physreve.105.054608] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Accepted: 05/11/2022] [Indexed: 05/03/2023]
Abstract
Alcohols are organic compounds characterized by one or more hydroxyl groups attached to a carbon atom of an alkyl group. They can be considered as organic derivatives of water in which one of the hydrogen atoms is replaced by an alkyl group. In this work, the Mercedes-Benz model of water is used to design simple two-dimensional (2D) models of lower alcohols. The structural and thermodynamic properties of the constructed simple models are studied by conducting Monte Carlo simulations in the isothermal-isobaric ensemble. We show that 2D models display similar trends in structuring and thermodynamics as in experiments. The present work on the smallest amphiphilc organic solutes provides a simple testing ground to study the competition between polar and non-polar effects within the molecule and physical properties.
Collapse
Affiliation(s)
- Petra Papez
- Laboratory for Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, SI-1001, Slovenia and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, SI-1000, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
10
|
Ogrin P, Urbic T. Isothermal-isobaric algorithm to study the effects of rotational degrees of freedom-Benz water model. J Mol Liq 2022; 349:118152. [PMID: 37727581 PMCID: PMC10508877 DOI: 10.1016/j.molliq.2021.118152] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have developed isothermal-isobaric algorithm for non-equilibrium Monte Carlo simulations. As first we have shown that the new method correctly predict density by comparing it to the density determined in canonical Monte Carlo simulations through the virial pressure. The new method was then used to study the effect of translational and rotational degrees of freedom on the structural and thermodynamic properties of the simple Mercedes-Benz water model. By holding one of the temperatures constant and varying the other one, we investigated how the position of the density maxima changes. We have observed that upon increase of rotational temperature the fluid become more Lennard-Jones like and the density maxima disappears.
Collapse
Affiliation(s)
- Peter Ogrin
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
11
|
|
12
|
Thermodynamic perturbation theory for rotational degrees of freedom. Application to the Mercedes–Benz water model. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115671] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
13
|
Ogrin P, Urbic T. Integral equation study of the effects of rotational degrees of freedom on properties of the Mercedes–Benz water model. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.114880] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
14
|
Urbic T. Two dimensional fluid with one site-site associating point. Monte Carlo, integral equation and thermodynamic perturbation theory study. J Mol Liq 2018; 270:87-96. [PMID: 30546180 DOI: 10.1016/j.molliq.2017.12.111] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
In this paper we propose a model for the two dimensional fluid with one site-site associating point. We studied its structural and thermodynamic properties by the Monte Carlo computer simulations, the site-site integral equation theory (RISM), the Wertheim's thermodynamic perturbation theory (TPT) and the Wertheim's integral equation theory (WIET) for associative liquids. The model can have arbitrary position of the associating point from the center of particles. All particles have Lennard-Jones core while interactions between associating points are modeled as Gaussian like potential where the interaction depends only on the distance between sites. The methods were used to study the thermodynamic and structural properties as a function of the position of associating point, temperature and density. The accuracy of the analytic theories were checked by comparing the theoretical results with the corresponding Monte Carlo ones. The theories are quite accurate for cases when the associating point is on the surface and only dimers can be formed. In this case, the theories correctly predict the pair correlation functions of the model, internal energy, ratios of free and bonded particles and chemical potential. This is no longer true when associating point is away from the surface of particles and the higher clusters are formed.
Collapse
Affiliation(s)
- Tomaz Urbic
- University of Ljubljana, Department of Chemistry and Chemical Technology, Chair of Physical Chemistry, Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
15
|
Urbic T, Dill KA. Hierarchy of anomalies in the two-dimensional Mercedes-Benz model of water. Phys Rev E 2018; 98:10.1103/physreve.98.032116. [PMID: 32025599 PMCID: PMC7001678 DOI: 10.1103/physreve.98.032116] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
We investigate by Monte Carlo simulations density, diffusion, and structural anomalies of the simple two-dimensional Mercedes-Benz (MB) model of water, which is a very simple toy model for explaining the origin of water properties. MB water molecules are modeled as two-dimensional Lennard-Jones disks, with three orientation-dependent hydrogen-bonding arms, arranged as in the MB logo. The model is in a way also a variance of silica-like models. Beside the known thermodynamic anomaly for the model we also found diffusion and structural anomalies and map out the cascade of density, structural, pair entropy, and diffusivity anomalies for MB model. The orientational order parameters with three and six-fold symmetry were determined and maximum for each one observed. The anomalies occur in hierarchy order, which is a slight variation of the hierarchy order in real water. The diffusion anomaly region is the innermost in the hierarchy while for water it is the density anomaly region.
Collapse
Affiliation(s)
- Tomaz Urbic
- University of Ljubljana, Faculty of Chemistry and Chemical Technology, Večna pot 113, SI-1000 Ljubljana, Slovenia
| | - Ken A. Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-5252, USA
| |
Collapse
|
16
|
Urbic T. Modelling water with simple Mercedes-Benz models. MOLECULAR SIMULATION 2018; 45:279-294. [PMID: 31156291 PMCID: PMC6542362 DOI: 10.1080/08927022.2018.1502430] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/24/2018] [Accepted: 07/14/2018] [Indexed: 10/28/2022]
Abstract
The structures and properties of biomolecules like proteins, nucleic acids, and membranes depend on water. Water is also very important in industry. Overall, water is unusual substance with more than 70 anomalous properties. The understanding of water is advancing significantly due to theoretical and computational modeling. There are different kind of models, models with fine-scale properties and increasing structural detail with increasing computational expense and simple models which focus on global properties of water like thermodynamics, phase diagram and are less computational expensive. Simplified models give a better understanding of water in ways that complement more complex models. Here, we review a simple model, the two dimensional Mercedes-Benz (MB) model of water. We present results by Monte Carlo simulations for anomalies and phase diagram and application of various theoretical methods.
Collapse
Affiliation(s)
- Tomaz Urbic
- University of Ljubljana, Faculty of Chemistry and Chemical Technology,
Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
17
|
Urbic T. Liquid part of the phase diagram and percolation line for two-dimensional Mercedes-Benz water. Phys Rev E 2018; 96:032122. [PMID: 29346988 DOI: 10.1103/physreve.96.032122] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2017] [Indexed: 11/07/2022]
Abstract
Monte Carlo simulations and Wertheim's thermodynamic perturbation theory (TPT) are used to predict the phase diagram and percolation curve for the simple two-dimensional Mercedes-Benz (MB) model of water. The MB model of water is quite popular for explaining water properties, but the phase diagram has not been reported till now. In the MB model, water molecules are modeled as two-dimensional Lennard-Jones disks, with three orientation-dependent hydrogen-bonding arms, arranged as in the MB logo. The liquid part of the phase space is explored using grand canonical Monte Carlo simulations and two versions of Wertheim's TPT for associative fluids, which have been used before to predict the properties of the simple MB model. We find that the theory reproduces well the physical properties of hot water but is less successful at capturing the more structured hydrogen bonding that occurs in cold water. In addition to reporting the phase diagram and percolation curve of the model, it is shown that the improved TPT predicts the phase diagram rather well, while the standard one predicts a phase transition at lower temperatures. For the percolation line, both versions have problems predicting the correct position of the line at high temperatures.
Collapse
Affiliation(s)
- T Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna Pot 113, 1000 Ljubljana, Slovenia
| |
Collapse
|
18
|
Integral equation and thermodynamic perturbation theory for a two-dimensional model of chain-forming fluid. J Mol Liq 2017; 238:129-135. [PMID: 28729752 DOI: 10.1016/j.molliq.2017.04.095] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
In this paper we applied analytical theories for the two dimensional chain-forming fluid. Wertheims thermodynamic perturbation theory (TPT) and integral equation theory (IET) for associative liquids were used to study thermodynamical and structural properties of the chain-forming model. The model has polymerizing points at arbitrary position from center of the particles. Calculated analytical results were tested against corresponding results obtained by Monte Carlo computer simulations to check the accuracy of the theories. The theories are accurate for the different positions of patches of the model at all values of the temperature and density studied. The IET's pair correlation functions of the model agree well with computer simulations. Both TPT and IET are in good agreement with the Monte Carlo values of the energy, chemical potential and ratios of free, once and twice bonded particles.
Collapse
|
19
|
Urbic T. Integral equation and thermodynamic perturbation theory for a two-dimensional model of dimerising fluid. J Mol Liq 2017; 228:32-37. [PMID: 28529396 DOI: 10.1016/j.molliq.2016.09.061] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
In this paper we applied an analytical theory for the two dimensional dimerising fluid. We applied Wertheims thermodynamic perturbation theory (TPT) and integral equation theory (IET) for associative liquids to the dimerising model with arbitrary position of dimerising points from center of the particles. The theory was used to study thermodynamical and structural properties. To check the accuracy of the theories we compared theoretical results with corresponding results obtained by Monte Carlo computer simulations. The theories are accurate for the different positions of patches of the model at all values of the temperature and density studied. IET correctly predicts the pair correlation function of the model. Both TPT and IET are in good agreement with the Monte Carlo values of the energy, pressure, chemical potential, compressibility and ratios of free and bonded particles.
Collapse
Affiliation(s)
- Tomaz Urbic
- University of Ljubljana, Department of Chemistry and Chemical Technology, Chair of Physical Chemistry, Večna pot 113, SI-1000 Ljubljana, Slovenia
| |
Collapse
|
20
|
The scientific life of Vojko Vlachy. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2016.11.073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
21
|
Urbic T. Properties of the two-dimensional heterogeneous Lennard-Jones dimers: An integral equation study. J Chem Phys 2016; 145:194503. [PMID: 27875894 PMCID: PMC5116027 DOI: 10.1063/1.4967807] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2016] [Accepted: 11/02/2016] [Indexed: 11/14/2022] Open
Abstract
Structural and thermodynamic properties of a planar heterogeneous soft dumbbell fluid are examined using Monte Carlo simulations and integral equation theory. Lennard-Jones particles of different sizes are the building blocks of the dimers. The site-site integral equation theory in two dimensions is used to calculate the site-site radial distribution functions and the thermodynamic properties. Obtained results are compared to Monte Carlo simulation data. The critical parameters for selected types of dimers were also estimated and the influence of the Lennard-Jones parameters was studied. We have also tested the correctness of the site-site integral equation theory using different closures.
Collapse
Affiliation(s)
- Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Vecna pot 113, 1000 Lubljana, Slovenia
| |
Collapse
|
22
|
Kurtjak M, Urbic T. A simple water model in the presence of inert Lennard-Jones obstacles II: the hydrophobic effect. Mol Phys 2014. [DOI: 10.1080/00268976.2014.973919] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
23
|
|
24
|
Cartwright JHE, Piro O, Sánchez PA, Sintes T. Ice polyamorphism in the minimal Mercedes-Benz model of water. J Chem Phys 2012; 137:244503. [PMID: 23277941 DOI: 10.1063/1.4772202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate ice polyamorphism in the context of the two-dimensional Mercedes-Benz model of water. We find a first-order phase transition between a crystalline phase and a high-density amorphous phase. Furthermore, we find a reversible transformation between two amorphous structures of high and low density; however, we find this to be a continuous and not an abrupt transition, as the low-density amorphous phase does not show structural stability. We discuss the origin of this behavior and its implications with regard to the minimal generic modeling of polyamorphism.
Collapse
Affiliation(s)
- Julyan H E Cartwright
- Instituto Andaluz de Ciencias de la Tierra, CSIC-Universidad de Granada, E-18071 Granada, Spain.
| | | | | | | |
Collapse
|
25
|
Kalyuzhnyi YV, Hlushak SP, Cummings PT. An improved thermodynamic perturbation theory for square-well m-point model of the patchy colloids. J Chem Phys 2012; 137:244910. [DOI: 10.1063/1.4773012] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Urbic T, Holovko MF. Mercedes-Benz water molecules near hydrophobic wall: integral equation theories vs Monte Carlo simulations. J Chem Phys 2011; 135:134706. [PMID: 21992334 DOI: 10.1063/1.3644934] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Associative version of Henderson-Abraham-Barker theory is applied for the study of Mercedes-Benz model of water near hydrophobic surface. We calculated density profiles and adsorption coefficients using Percus-Yevick and soft mean spherical associative approximations. The results are compared with Monte Carlo simulation data. It is shown that at higher temperatures both approximations satisfactory reproduce the simulation data. For lower temperatures, soft mean spherical approximation gives good agreement at low and at high densities while in at mid range densities, the prediction is only qualitative. The formation of a depletion layer between water and hydrophobic surface was also demonstrated and studied.
Collapse
Affiliation(s)
- T Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Ljubljana SI-1000, Slovenia.
| | | |
Collapse
|
27
|
Urbic T, Dill KA. A statistical mechanical theory for a two-dimensional model of water. J Chem Phys 2010; 132:224507. [PMID: 20550408 DOI: 10.1063/1.3454193] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop a statistical mechanical model for the thermal and volumetric properties of waterlike fluids. Each water molecule is a two-dimensional disk with three hydrogen-bonding arms. Each water interacts with neighboring waters through a van der Waals interaction and an orientation-dependent hydrogen-bonding interaction. This model, which is largely analytical, is a variant of the Truskett and Dill (TD) treatment of the "Mercedes-Benz" (MB) model. The present model gives better predictions than TD for hydrogen-bond populations in liquid water by distinguishing strong cooperative hydrogen bonds from weaker ones. We explore properties versus temperature T and pressure p. We find that the volumetric and thermal properties follow the same trends with T as real water and are in good general agreement with Monte Carlo simulations of MB water, including the density anomaly, the minimum in the isothermal compressibility, and the decreased number of hydrogen bonds for increasing temperature. The model reproduces that pressure squeezes out water's heat capacity and leads to a negative thermal expansion coefficient at low temperatures. In terms of water structuring, the variance in hydrogen-bonding angles increases with both T and p, while the variance in water density increases with T but decreases with p. Hydrogen bonding is an energy storage mechanism that leads to water's large heat capacity (for its size) and to the fragility in its cagelike structures, which are easily melted by temperature and pressure to a more van der Waals-like liquid state.
Collapse
Affiliation(s)
- Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Askerceva 5, 1000 Ljubljana, Slovenia.
| | | |
Collapse
|
28
|
Bizjak A, Urbic T, Vlachy V, Dill KA. Theory for the three-dimensional Mercedes-Benz model of water. J Chem Phys 2010; 131:194504. [PMID: 19929057 DOI: 10.1063/1.3259970] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The two-dimensional Mercedes-Benz (MB) model of water has been widely studied, both by Monte Carlo simulations and by integral equation methods. Here, we study the three-dimensional (3D) MB model. We treat water as spheres that interact through Lennard-Jones potentials and through a tetrahedral Gaussian hydrogen bonding function. As the "right answer," we perform isothermal-isobaric Monte Carlo simulations on the 3D MB model for different pressures and temperatures. The purpose of this work is to develop and test Wertheim's Ornstein-Zernike integral equation and thermodynamic perturbation theories. The two analytical approaches are orders of magnitude more efficient than the Monte Carlo simulations. The ultimate goal is to find statistical mechanical theories that can efficiently predict the properties of orientationally complex molecules, such as water. Also, here, the 3D MB model simply serves as a useful workbench for testing such analytical approaches. For hot water, the analytical theories give accurate agreement with the computer simulations. For cold water, the agreement is not as good. Nevertheless, these approaches are qualitatively consistent with energies, volumes, heat capacities, compressibilities, and thermal expansion coefficients versus temperature and pressure. Such analytical approaches offer a promising route to a better understanding of water and also the aqueous solvation.
Collapse
Affiliation(s)
- Alan Bizjak
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Askerceva 5, 1000 Ljubljana, Slovenia
| | | | | | | |
Collapse
|
29
|
Zhou S, Solana JR. Progress in the Perturbation Approach in Fluid and Fluid-Related Theories. Chem Rev 2009; 109:2829-58. [DOI: 10.1021/cr900094p] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Shiqi Zhou
- State Key Laboratory of Powder Metallurgy, Central South University, Changsha, Hunan 410083, China, and School of Physics Science and Technology, Central South University, Changsha, Hunan 410083, China
| | - J. R. Solana
- Applied Physics Department, University of Cantabria, 39005 Santander, Spain
| |
Collapse
|