1
|
Trejos VM, Gámez F, Torres-Carbajal A, Martínez-Borquez A. Monte Carlo simulations and perturbation theory for highly correlated fluids: The Lennard-Jones core softened potential case. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2019.112201] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
2
|
Trejos VM, Pizio O, Sokołowski S. On the interdigitation of molecular brushes and solvation force upon adsorption of water in slit-like pores with grafted chains. Density functional approach. J Chem Phys 2019. [DOI: 10.1063/1.5116128] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Víctor M. Trejos
- Instituto de Ciencias Básicas e Ingeniería, Universidad Autónoma del Estado de Hidalgo, Carretera Pachuca-Tulancingo Km. 4.5, Col. Carboneras, C.P. 42184, Mineral de la Reforma, Hidalgo, Mexico
| | - Orest Pizio
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, 04510, Cd. de México, Mexico
| | - Stefan Sokołowski
- Department for the Modelling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin 20-031, Poland
| |
Collapse
|
3
|
de Haro ML, Rodríguez-Rivas Á, Yuste SB, Santos A. Structural properties of the Jagla fluid. Phys Rev E 2018; 98:012138. [PMID: 30110748 DOI: 10.1103/physreve.98.012138] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Indexed: 06/08/2023]
Abstract
The structural properties of the Jagla fluid are studied by Monte Carlo (MC) simulations, numerical solutions of integral equation theories, and the (semi-analytical) rational-function approximation (RFA) method. In the latter case, the results are obtained from the assumption (supported by our MC simulations) that the Jagla potential and a potential with a hard core plus an appropriate piecewise constant function lead to practically the same cavity function. The predictions obtained for the radial distribution function, g(r), from this approach are compared against MC simulations and integral equations for the Jagla model, and also for the limiting cases of the triangle-well potential and the ramp potential, with a general good agreement. The analytical form of the RFA in Laplace space allows us to describe the asymptotic behavior of g(r) in a clean way and compare it with MC simulations for representative states with oscillatory or monotonic decay. The RFA predictions for the Fisher-Widom and Widom lines of the Jagla fluid are obtained.
Collapse
Affiliation(s)
- Mariano López de Haro
- Instituto de Energías Renovables, Universidad Nacional Autónoma de México (U.N.A.M.), Temixco, Morelos 62580, Mexico
| | - Álvaro Rodríguez-Rivas
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| | - Santos B Yuste
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| | - Andrés Santos
- Departamento de Física and Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, Badajoz E-06006, Spain
| |
Collapse
|
4
|
Baran Ł, Sokołowski S. Effective interactions between a pair of particles modified with tethered chains. J Chem Phys 2018; 147:044903. [PMID: 28764361 DOI: 10.1063/1.4994919] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Using molecular dynamics, we evaluate the potential of mean force for two models of hybrid nanoparticles, namely, for the models with fixed and movable chain ligands. We also investigate the structure of segments of chains around nanoparticles and its change when one nanoparticle approaches the other. In the case of an isolated particle, we also employ a density functional theory to compute the segment density profiles. Moreover, to determine the structure of segments around a core, we have employed the concept of the so-called mass dipoles.
Collapse
Affiliation(s)
- Ł Baran
- Department for the Modelling of Physico-Chemical Processes, Maria Curie-Sklodowska University, Gliniana 33, Lublin, Poland
| | - S Sokołowski
- Department for the Modelling of Physico-Chemical Processes, Maria Curie-Sklodowska University, Gliniana 33, Lublin, Poland
| |
Collapse
|
5
|
Calculating free energy profiles using entropy as a reaction coordinate: Application to water nucleation. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.02.026] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
6
|
Lukšič M, Hribar-Lee B, Pizio O. Phase behaviour of a continuous shouldered well model fluid. A grand canonical Monte Carlo study. J Mol Liq 2017; 228:4-10. [PMID: 28450755 PMCID: PMC5403148 DOI: 10.1016/j.molliq.2016.10.030] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
The phase behavior of the continuous shouldered well model fluid proposed by Franzese [J. Mol. Liq. 136 (2007) 267] was examined using the Monte Carlo computer simulations in the grand canonical ensemble. The essential parts of the vapour-liquid and liquid-liquid coexistence envelopes were obtained. The Widom lines departing from coexistence envelopes were calculated using maxima of the fluctuations of the number of particles as a function of chemical potential along various isotherms. The region embracing anomalies in the properties of the model was located using the approximate criterion that involves the excess pair entropy.. The temperature of maximum density line was built by performing canonical Monte Carlo simulations. Our results are consistent with previous results from molecular dynamics constant pressure-constant temperature simulations and provide wider insight into the phase behavior of the model by using the chemical potential as the external parameter.
Collapse
Affiliation(s)
- Miha Lukšič
- Faculty of Chemistry and Chemical Technology, University of
Ljubljana, Večna pot 113, SI-1000 Ljubljana, Slovenia
| | - Barbara Hribar-Lee
- Faculty of Chemistry and Chemical Technology, University of
Ljubljana, Večna pot 113, SI-1000 Ljubljana, Slovenia
| | - Orest Pizio
- Instituto de Química, Universidad Nacional Autonoma de
México, Coyoacan, 04510, Cd. de México, México
| |
Collapse
|
7
|
Borówko M, Rżysko W, Sokołowski S, Staszewski T. Integral equations theory for two-dimensional systems involving nanoparticles. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1234079] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Małgorzata Borówko
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| | - Wojciech Rżysko
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| | - Stefan Sokołowski
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| | - Tomasz Staszewski
- Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska University, Lublin, Poland
| |
Collapse
|
8
|
Torres-Carbajal A, Castañeda-Priego R. Characterisation of the thermodynamics, structure and dynamics of a water-like model in 2- and 3-dimensions. Phys Chem Chem Phys 2016; 18:17335-40. [PMID: 27232761 DOI: 10.1039/c6cp01565d] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The physical properties of colloidal particles suspended in an aqueous environment are well-understood when the latter is considered to be a continuum and a structureless medium. However, this approach fails to explain complex phenomena, for example, the critical Casimir forces among colloids and the colloidal self-assembly near critical solvents, and the inertial contribution of the solvent molecules on the diffusion of non-spherical Brownian particles. Therefore, the role played by the solvent on the physical properties of colloidal dispersions is of paramount relevance. Recently, there has been an interest in the (non-trivial) diffusion mechanisms of a nano-colloidal particle in a solvent that undergoes a vapour-liquid transition. Nonetheless, the models typically used to incorporate the solvent details do not capture quantitatively the thermodynamic properties of real substances. It is then important to study the Brownian motion of colloids in more realistic models. To reach such goal, one first has to characterise the thermodynamic states and the microscopic features of the solvent. Hence, in this contribution, we have investigated the coexistence densities of a core-softened potential in two- and three-dimensions, whose potential parameters are able to capture some anomalies of water. We show that in the two-dimensional case, the potential model exhibits, besides the normal vapour-liquid coexistence region, additional liquid-liquid coexistence densities. We particularly focus our attention to the structural properties and the dynamical behaviour of the solvent around the liquid-liquid critical point and assess the differences with the three-dimensional case.
Collapse
Affiliation(s)
- Alexis Torres-Carbajal
- División de Ciencias e Ingenierías, Campus León, Universidad de Guanajuato, Loma del Bosque 103, Lomas del Campestre, 37150 León, Guanajuato, México.
| | | |
Collapse
|
9
|
Urbic T. Liquid-liquid phase transition in a two-dimensional system with anomalous liquid properties. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 88:062303. [PMID: 24483440 DOI: 10.1103/physreve.88.062303] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/25/2013] [Revised: 10/03/2013] [Indexed: 06/03/2023]
Abstract
The phase diagram of the two-dimensional particles interacting through a smooth version of Stell-Hemmer interaction was studied using Monte Carlo computer simulations. By evaluating the pressure-volume isotherms, we observed liquid-liquid, liquid-gas phase transitions and three stable crystal phases. The model shows the liquid-liquid critical point in stable liquid phase and is confirmed by observing properties of other thermodynamic functions such as heat capacity and isothermal compressibility, for example. The liquid-gas and the liquid-liquid critical points were estimated within the thermodynamic limit.
Collapse
Affiliation(s)
- Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, 1000, Slovenia
| |
Collapse
|
10
|
Huš M, Urbic T. Core-softened fluids as a model for water and the hydrophobic effect. J Chem Phys 2013; 139:114504. [PMID: 24070294 DOI: 10.1063/1.4821226] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An interaction model with core-softened potential in three dimensions was studied by Monte Carlo computer simulations and integral equation theory. We investigated the possibility that a fluid with a core-softened potential can reproduce anomalies found experimentally in liquid water, such as the density anomaly, the minimum in the isothermal compressibility as a function of temperature, and others. Critical points of the fluid were also determined. We provided additional arguments that the old notion, postulating that only angular-dependent interactions result in density anomaly, is incorrect. We showed that potential with two characteristic distances is sufficient for the system to exhibit water-like behavior and anomalies, including the famous density maximum. We also found that this model can properly describe the hydrophobic effect.
Collapse
Affiliation(s)
- Matej Huš
- Department of Chemistry and Chemical Engineering, Chair of Physical Chemistry, University of Ljubljana, Aškerčeva 5, SI-1000 Ljubljana, Slovenia
| | | |
Collapse
|
11
|
Prestipino S, Saija F, Giaquinta PV. Hexatic phase and water-like anomalies in a two-dimensional fluid of particles with a weakly softened core. J Chem Phys 2012; 137:104503. [DOI: 10.1063/1.4749260] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Agarwal M, Alam MP, Chakravarty C. Thermodynamic, diffusional, and structural anomalies in rigid-body water models. J Phys Chem B 2011; 115:6935-45. [PMID: 21553909 DOI: 10.1021/jp110695t] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Structural, density, entropy, and diffusivity anomalies of the TIP4P/2005 model of water are mapped out over a wide range of densities and temperatures. The locus of temperatures of maximum density (TMD) for this model is very close to the experimental TMD locus for temperatures between 250 and 275 K. Four different water models (mTIP3P, TIP4P, TIP5P, and SPC/E) are compared with the TIP4P/2005 model in terms of their anomalous behavior. For all the water models, the density regimes for anomalous behavior are bounded by a low-density limit at around 0.85-0.90 g cm(-3) and a high-density limit at about 1.10-1.15 g cm(-3). The onset temperatures of the density anomaly in the various models show a much greater variation, ranging from 202 K for mTIP3P to 289 K for TIP5P. The order maps for the various water models are qualitatively very similar with the structurally anomalous regions almost superimposable in the q(tet)-τ plane. Comparison of the phase diagrams of water models with the region of liquid-state anomalies shows that the crystalline phases are much more sensitive to the choice of water models than the liquid state anomalies; for example, SPC/E and TIP4P/2005 show qualitatively similar liquid state anomalies but very different phase diagrams. The anomalies in the liquid in all the models occur at much lower pressures than those at which the melting line changes from negative to positive slope. The results in this study demonstrate several aspects of structure-entropy-diffusivity relationships of water models that can be compared with experiment and used to develop better atomistic and coarse-grained models for water.
Collapse
Affiliation(s)
- Manish Agarwal
- Department of Chemistry, Indian Institute of Technology-Delhi, New Delhi, India
| | | | | |
Collapse
|
13
|
Ahmed A, Mausbach P, Sadus RJ. Strain-rate dependent shear viscosity of the Gaussian core model fluid. J Chem Phys 2009; 131:224511. [DOI: 10.1063/1.3273083] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Saija F, Prestipino S, Malescio G. Anomalous phase behavior of a soft-repulsive potential with a strictly monotonic force. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:031502. [PMID: 19905115 DOI: 10.1103/physreve.80.031502] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/23/2009] [Indexed: 05/28/2023]
Abstract
We study the phase behavior of a classical system of particles interacting through a strictly convex soft-repulsive potential which, at variance with the pairwise softened repulsions considered so far in the literature, lacks a region of downward or zero curvature. Nonetheless, such interaction is characterized by two length scales, owing to the presence of a range of interparticle distances where the repulsive force increases, for decreasing distance, much more slowly than in the adjacent regions. We investigate, using extensive Monte Carlo simulations combined with accurate free-energy calculations, the phase diagram of the system under consideration. We find that the model exhibits a fluid-solid coexistence line with multiple re-entrant regions, an extremely rich solid polymorphism with solid-solid transitions, and waterlike anomalies. In spite of the isotropic nature of the interparticle potential, we find that among the crystal structures in which the system can exist, there are also a number of non-Bravais lattices, such as cI16 and diamond.
Collapse
Affiliation(s)
- Franz Saija
- CNR-Istituto per i Processi Chimico-Fisici, Contrada Papardo, 98158 Messina, Italy.
| | | | | |
Collapse
|