BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Rakesh S. Singh

Rakesh S. Singh contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechcond-mat.softcond-mat.mtrl-sciphysics.chem-ph

Trust snapshot

Quick read

Trust 21 - Emerging
10works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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)

preprint2016arXiv

An Equation of State For The TIP4P/2005 Model of Water Including Negative Pressures

One of the most promising frameworks for understanding the anomalies of cold and especially supercooled water is that of two-structure thermodynamics, in which water is viewed as a non-ideal mixture of two interconvertible local structures. The non-ideality of this mixture may give rise, at very low temperatures, to a liquid-liquid phase transition (LLPT) and a liquid-liquid critical point (LLCP), at which thermodynamic response functions diverge. Various versions of the "two-structure equation of state" (TSEOS) based on this concept have shown remarkable agreement with both experimental data in real water and simulation results. However, recent experiments probing supercooled water at negative pressures reveal the inadequacy of extrapolations of equations of state developed for positive pressures, and have begun to shed additional light on the source of the anomalies of supercooled water. We have analyzed simulation results for the TIP4P/2005 model over a broad range of positive and negative pressures from ambient temperature to deep supercooling. We find that by explicitly incorporating a liquid-vapor spinodal into a two-structure equation of state, we are able to match the simulation data in TIP4P/2005 with striking accuracy. In particular, our equation of state reproduces the observed lines of minima and maxima in the density, isothermal compressibility, and isobaric heat capacity. Contrary to scenarios in which a retracing spinodal accounts for the thermodynamic anomalies of water, we find that the liquid-vapor spinodal in this model continues monotonically to lower pressures as temperature is decreased, influencing but not giving rise to the locus of density maxima and other thermodynamic anomalies. We explain the behavior of TIP4P/2005 in terms of two phenomena: the competition between two local structures and a monotonic liquid-vapor spinodal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Density and bond-orientational relaxations in supercooled water

Recent computational studies have reported evidence of a metastable liquid-liquid phase transition (LLPT) in molecular models of water under deeply supercooled conditions. A competing hypothesis suggests, however, that non-equilibrium artifacts associated with coarsening of the stable crystal phase have been mistaken for an LLPT in these models. Such artifacts are posited to arise due to a separation of time scales in which density fluctuations in the supercooled liquid relax orders of magnitude faster than those associated with bond-orientational order. Here, we use molecular simulation to investigate the relaxation of density and bond-orientational fluctuations in three molecular models of water (ST2, TIP5P and TIP4P/2005) in the vicinity of their reported LLPT. For each model, we find that density is the slowly relaxing variable under such conditions. We also observe similar behavior in the coarse-grained mW model of water. Our findings therefore challenge the key physical assumption underlying the competing hypothesis.e find that density relaxes significantly faster than bond-orientational order, as incorrectly predicted by this competing hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two-State Thermodynamics and the Possibility of a Liquid-Liquid Phase Transition in Supercooled TIP4P/2005 Water

Water shows intriguing thermodynamic and dynamic anomalies in the supercooled liquid state. One possible explanation of the origin of these anomalies lies in the existence of a metastable liquid-liquid phase transition (LLPT) between two (high and low density) forms of water. While the anomalies are observed in experiments on bulk and confined water and by computer simulation studies of water-like models, the existence of a LLPT in water is still debated. Unambiguous experimental proof of the existence of a LLPT in bulk supercooled water is hampered by fast ice nucleation which is a precursor of the hypothesized LLPT. Moreover, the hypothesized LLPT, being metastable, in principle cannot exist in the thermodynamic limit (infinite size, infinite time). Therefore, computer simulations of water models are crucial for exploring the possibility of the metastable LLPT and the nature of the anomalies. In this work, we present new simulation results in the NVT ensemble for one of the most accurate classical molecular models of water, TIP4P/2005. To describe the computed properties and explore the possibility of a LLPT we have applied two-structure thermodynamics, viewing water as a non-ideal mixture of two interconvertible local structures ("states"). The results suggest the presence of a liquid-liquid critical point and a LLPT in this model for the simulated length and time scales. We have compared the behavior of TIP4P/2005 with other popular water-like models, namely mW and ST2, and with real water, all of which are well described by two-state thermodynamics. In view of the current debate involving different studies of TIP4P/2005, we discuss consequences of metastability and finite size in observing the liquid-liquid separation. We also address the relationship between the phenomenological order parameter of two-structure thermodynamics and the microscopic nature of the low-density structure.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Orientational order as the origin of the long-range hydrophobic effect

The long range attractive force between two hydrophobic surfaces immersed in water is observed to decrease exponentially with their separation -- this distance-dependence of effective force is known as the hydrophobic force law (HFL). We explore the microscopic origin of HFL by studying distance-dependent attraction between two parallel rods immersed in 2D Mercedes Benz model of water. This model is found to exhibit a well-defined HFL. Although the phenomenon is conventionally explained by density-dependent theories, we identify orientation, rather than density, as the relevant order parameter. The range of density variation is noticeably shorter than that of orientational heterogeneity. The latter is comparable to the observed distances of hydrophobic force. At large separation, attraction between the rods arises primarily from a destructive interference among the inwardly propagating oppositely oriented heterogeneity generated in water by the two rods. As the rods are brought closer, the interference increases leading to a decrease in heterogeneity and concomitant decrease in free energy of the system, giving rise to the effective attraction. We notice formation of hexagonal ice-like structures at the onset of attractive region which suggests that metastable free energy minimum may play a role in the origin of HFL.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Correlation between thermodynamic anomalies and pathways of ice nucleation in supercooled water

The well-known classical nucleation theory (CNT) for the free energy barrier towards formation of a nucleus of critical size of the new stable phase within the parent metastable phase fails to take into account the influence of other metastable phases having density/order intermediate between the parent metastable phase and the final stable phase. This lacuna can be more serious than capillary approximation or spherical shape assumption made in CNT. This issue is particularly significant in ice nucleation because liquid water shows rich phase diagram consisting of two (high and low density) liquid phases in supercooled state. The explanations of thermodynamic and dynamic anomalies of supercooled water often invoke the possible influence of a liquid-liquid transition between two metastable liquid phases. To investigate both the role of thermodynamic anomalies and presence of distinct metastable liquid phases in supercooled water on ice nucleation, we employ density functional theoretical approach to find nucleation free energy barrier in different regions of phase diagram. The theory makes a number of striking predictions, such as a dramatic lowering of nucleation barrier due to presence of a metastable intermediate phase and crossover in the dependence of free energy barrier on temperature near liquid-liquid critical point. These predictions can be tested by computer simulations as well as by controlled experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nucleation of a stable solid from melt in the presence of multiple metastable intermediate phases: Wetting, Ostwald step rule and vanishing polymorphs

In many systems, nucleation of a stable solid may occur in the presence of other (often more than one) metastable phases. These may be polymorphic solids or even liquid phases. In such cases, nucleation of the solid phase from the melt may be facilitated by the metastable phase because the latter can "wet" the interface between the parent and the daughter phases, even though there may be no signature of the existence of metastable phase in the thermodynamic properties of the parent liquid and the stable solid phase. Straightforward application of classical nucleation theory (CNT) is flawed here as it overestimates the nucleation barrier since surface tension is overestimated (by neglecting the metastable phases of intermediate order) while the thermodynamic free energy gap between daughter and parent phases remains unchanged. In this work we discuss a density functional theory (DFT) based statistical mechanical approach to explore and quantify such facilitation. We construct a simple order parameter dependent free energy surface that we then use in DFT to calculate (i) the order parameter profile, (ii) the overall nucleation free energy barrier and (iii) the surface tension between the parent liquid and the metastable solid and also parent liquid and stable solid phases. The theory indeed finds that the nucleation free energy barrier can decrease significantly in the presence of wetting. This approach can provide a microscopic explanation of Ostwald step rule and the well-known phenomenon of "disappearing polymorphs" that depends on temperature and other thermodynamic conditions. Theory reveals a diverse scenario for phase transformation kinetics some of which may be explored via modern nanoscopic synthetic methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solid-solid collapse transition in a two dimensional model molecular system

Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real difficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises from the interplay between different energy and length scales involved in molecular systems, often mediated through a solvent. In this work we employ Monte Carlo (MC) simulations to study the collapse transition in a model molecular system interacting via both isotropic as well as anisotropic interactions having different length and energy scales. The model we use is known as Mercedes-Benz (MB) which for a specific set of parameters sustains three solid phases: honeycomb, oblique and triangular. In order to study the temperature induced collapse transition, we start with a metastable honeycomb solid and induce transition by heating. High density oblique solid so formed has two characteristic length scales corresponding to isotropic and anisotropic parts of interaction potential. Contrary to the common believe and classical nucleation theory, interestingly, we find linear strip-like nucleating clusters having significantly different order and average coordination number than the bulk stable phase. In the early stage of growth, the cluster grows as linear strip followed by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate balance between two types of interactions which enables the dominance of stabilizing energy over the destabilizing surface energy. The nuclei of stable oblique phase are wetted by intermediate order particles which minimizes the surface free energy. We observe different pathways for pressure and temperature induced transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Strongly first order melting of a two dimensional molecular solid

Melting and freezing transitions in two dimensional systems are known to show highly unusual characteristics. Most of the earlier studies considered atomic systems; the melting behavior in two dimensional molecular solids is still largely unexplored. In order to understand the role of multiple energy and length scales present in molecular systems on nature of melting transition, here we report computer simulation studies of melting of a two dimensional Mercedes-Benz (MB) system. We find that the interplay between the strength of isotropic and anisotropic interactions can give rise to rich phase diagram. The computed solid-liquid phase diagram consists of isotropic liquid and two crystalline phases - honeycomb and oblique. In contradiction to the celebrated KTHNY theory, we observe strongly one step first order melting transitions for both the honeycomb and oblique solids. The defects in both solids and liquids near the transition are more complex compared to the atomic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gas-liquid Nucleation at Large Metastability

Nucleation at large metastability is still largely an unsolved problem, although is a problem of tremendous current interest, with wide practical value. It is well-accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this article, we present a powerful alternative formalism to treat nucleation at large supersaturation. This formalism goes over to the classical picture at small supersaturation where CNT is expected to be valid. The new theory is based on an extended set of order parameters in terms of k-th largest liquid-like clusters where k=1 is the largest cluster in the system, k=2 is the second largest cluster and so on. We derive an analytic expression for the free energy of formation of the k-th largest cluster which shows that at large metastability the barrier of growth for the few largest liquid-like clusters disappear, the nucleation becomes collective and the approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts a weaker supersaturation dependence than that of CNT at large metastability. Such a cross-over behavior has indeed been observed in recent experiments but eluded an explanation till now. In order to understand the large numerical difference between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interaction of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for a Lennard-Jones potential, both in two and three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interplay between the range of attractive potential and metastability in gas-liquid nucleation

We find an interesting interplay between the range of the attractive part of the interaction potential and the extent of metastability (as measured by supersaturation) in gas-liquid nucleation. We explore and exploit this interplay to obtain new insight into nucleation phenomena. Just like its dependence on supersaturation (S), the free energy barrier of nucleation is found to depend strongly on the range of the interaction potential. Actually, the entire free energy surface, F(n), where n is the size of the liquid-like cluster, shows this dependence. The evidences and the reasons for this strong dependence are as follows. (i) The surface tension increases dramatically on increasing the range of interaction potential. In three dimensional Lennard-Jones system, the value of the surface tension increases from 0.494 for a cut-off of 2.5 σ to 1.09 when the full range of the potential is involved. In two dimensional LJ system, the value of the line tension increases from 0.05 to 0.18, under the same variation of the potential range. (ii) The density of the gas phase at coexistence decreases while that of the liquid phase increases substantially on increasing the range of the interaction potential. (iii) As a result of the above, at a given supersaturation S, the size of the critical nucleus and the free energy barrier both increase with increase in the range of interaction potential. (iv) Surprisingly, however, we find that the functional form predicted by the classical nucleation theory (CNT) for the dependence of the free energy barrier on the size of the nucleus to remain valid except at the largest value of S studied. (v) The agreement between CNT prediction and simulated values of the barrier is supersaturation dependent and worsens with increase in the range of interaction potential, and increases above 10 kBT at the largest supersaturation that could be studied.

0 citations0 reviewsNo code linkedOpen access