BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Carlos Vega

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechgr-qcmath.DGphysics.chem-phcond-mat.mtrl-sci

Catalog footprint

What is connected

10works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

10 published item(s)

preprint2021arXiv

Spacetime distances: an exploration

What is the distance between two points in spacetime? This is a basic geometric question, which so far has no single, definitive answer. Unlike their Riemannian cousins, Lorentzian manifolds are not known to carry a canonical distance function. There is, however, a well-known way to `Riemannianize' a Lorentzian metric tensor, under suitable conditions, (e.g., `Wick rotation'). In a fairly different vein, a new `null distance function' was also introduced in [18]. Our goal here is to begin an exploration of the broad question of distance functions on spacetimes in general, including a concrete comparison of the aforementioned constructions. We devote special attention to the model `generalized Robertson-Walker' (GRW) setting, and also study how the `classical FLRW big bangs' manifest in terms of such distance functions, and associated metric completions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Homogeneous Ice Nucleation Rate in Water Droplets

To predict the radiative forcing of clouds it is necessary to know the rate with which ice homogeneously nucleates in supercooled water. Such rate is often measured in drops to avoid the presence of impurities. At large supercooling small (nanoscopic) drops must be used to prevent simultaneous nucleation events. The pressure inside such drops is larger than the atmospheric one by virtue of the Laplace equation. In this work, we take into account such pressure raise in order to predict the nucleation rate in droplets using the TIP4P/Ice water model. We start from a recent estimate of the maximum drop size that can be used at each supercooling avoiding simultaneous nucleation events [Espinosa et al. J. Chem. Phys., 2016]. We then evaluate the pressure inside the drops with the Laplace equation. Finally, we obtain the rate as a function of the supercooling by interpolating our previous results for 1 and 2000 bar [Espinosa et al. Phys. Rev. Lett. 2016] using the Classical Nucleation Theory expression for the rate. This requires, in turn, interpolating the ice-water interfacial free energy and chemical potential difference. The TIP4P/Ice rate curve thus obtained is in good agreement with most droplet-based experiments. In particular, we find a good agreement with measurements performed using nanoscopic drops, that are currently under debate. The successful comparison between model and experiments suggests that TIP4P/Ice is a reliable model to study the water-to-ice transition and that Classical Nucleation Theory is a good framework to understand it.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Seeding Approach to Bubble Nucleation in Superheated Lennard Jones Fluids

We investigate vapor homogeneous nucleation in a superheated Lennard Jones liquid with computer simulations. Special simulation techniques are required to address this study since the nucleation of a critical vapor bubble --one that has equal chance to grow or shrink-- in a moderately superheated liquid is a rare event. We use the Seeding method, that combines Classical Nucleation Theory with computer simulations of a liquid containing a vapor bubble to provide bubble nucleation rates in a wide temperature range. Seeding has been successfully applied to investigate the nucleation of crystals in supercooled fluids and here we apply it for the fist time --to the best of our knowledge-- to the liquid-to-vapor transition. We find that Seeding provides nucleation rates that are consistent with independent calculations not based on the assumptions of Classical Nucleation Theory. Different criteria to determine the radius of the critical bubble give different rate values. The accuracy of each criterion depends of the degree of superheating. Moreover, Seeding simulations show that the surface tension depends on pressure for a given temperature. Therefore, using Classical Nucleation Theory with the coexistence surface tension does not provide good estimates of the nucleation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Seeding Approach to nucleation in the NVT ensemble: the case of bubble cavitation in overstretched Lennard Jones fluids

Simulations are widely used to study nucleation in first order phase transitions due to the fact that they have access to the relevant length and time scales. However, simulations face the problem that nucleation is an activated process. Therefore, rare event simulation techniques are needed to promote the formation of the critical nucleus. The Seeding method, where the simulations are started with the nucleus already formed, has proven quite useful in efficiently providing estimates of the nucleation rate for a wide range of orders of magnitude. So far, Seeding has been employed in the NPT ensemble, where the nucleus either grows or redissolves. Thus, several trajectories have to be run in order to find the thermodynamic conditions that make the seeded nucleus critical. Moreover, the nucleus lifetime is short and the statistics for obtaining its properties is consequently poor. To deal with these shortcomings we extend the Seeding method to the NVT ensemble. We focus on the problem of bubble nucleation in a mestastable Lennard Jones fluid. We show that, in the NVT ensemble, it is possible to equilibrate and stabilise critical bubbles for a long time. The nucleation rate inferred from NVT-Seeding is fully consistent with that coming from NPT-Seeding. The former is quite suitable to obtain the nucleation rate along isotherms, whereas the latter is preferable if the dependence of the rate with temperature at constant pressure is required. Care should be taken with finite size effects when using NVT-Seeding. Further work is required to extend NVT seeding to other sorts of phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Water: one molecule, two surfaces, one mistake

In order to rigorously evaluate the energy and dipole moment of a certain configuration of molecules one needs to solve the Schrödinger equation. Repeating this for many different configurations allows one to determine the potential energy surface (PES) and the dipole moment surface (DMS). Since the early days of computer simulation it has been implicitly accepted that for empirical potentials the charges used to fit the PES should also be used to describe the DMS. This is a mistake. Partial charges are not observable magnitudes. They should be regarded as adjustable fitting parameters. Optimal values used to describe the PES are not necessarily the best to describe the DMS. One could use two fits: one for the PES, and another for the DMS. This is a common practice in the quantum chemistry community, but not used so often by the community performing computer simulations. This idea affects all types of modelling of water (with the exception of ab-initio calculations) from coarse grained to non-polarizable and polarizable models. We anticipate that an area that will benefit dramatically from having both, a good PES and a good DMS, is the modelling of water in the presence of electric fields.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Absence of superheating for ice Ih with a free surface : a new method of determining the melting point of different water models

Molecular dynamic simulations were performed for ice Ih with a free surface. The simulations were carried out at several temperatures and each run lasted more than 7ns. At high temperatures the ice melts. It is demonstrated that the melting process starts at the surface and propagates to the bulk of the ice block. Already at the temperatures below the melting point, we observe a thin liquid layer at the ice surface, but the block of ice remains stable along the run. As soon as the temperature reaches the melting point the entire ice block melts. Our results demonstrate that, unlike in the case of conventional simulations in the NpT ensemble, overheating of the ice Ih with a free surface does not occur. That allows to estimate the melting point of ice at zero pressure. We applied the method to the following models of water: SPC/E, TIP4P, TIP4P/Ew, TIP4P/Ice and TIP4P/2005, and found good agreement between the melting temperatures obtained by this procedure and the values obtained either from free energy calculations or from direct simulations of the ice/water interface.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Achronal limits, Lorentzian spheres, and splitting

In the early 80's S.-T. Yau posed the problem of establishing the rigidity of the Hawking-Penrose singularity theorems. Approaches to this problem have involved the introduction of Lorentzian Busemann functions and the study of the geometry of their level sets - the horospheres. The regularity theory in the Lorentzian case is considerably more complicated and less complete than in the Riemannian case. In this paper we introduce a broad generalization of the notion of horosphere in Lorentzian geometry and take a completely different (and highly geometric) approach to regularity. These generalized horospheres are defined in terms of 'achronal limits', and the improved regularity we obtain is based on regularity properties of achronal boundaries. We establish a splitting result for generalized horospheres, which when specialized to 'Cauchy horospheres' yields new results on the Bartnik splitting conjecture, a concrete realization of the problem posed by Yau. Our methods are also applied to spacetimes with positive cosmological constant. We obtain a rigid singularity result for future asymptotically de Sitter spacetimes related to results in [1,7].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water

We present a local order parameter based on the standard Steinhardt-Ten Wolde approach that is capable both of tracking and of driving homogeneous ice nucleation in simulations of all-atom models of water. We demonstrate that it is capable of forcing the growth of ice nuclei in supercooled liquid water simulated using the TIP4P/2005 model using overbiassed umbrella sampling Monte Carlo simulations. However, even with such an order parameter, the dynamics of ice growth in deeply supercooled liquid water in all-atom models of water are shown to be very slow, and so the computation of free energy landscapes and nucleation rates remains extremely challenging.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Path integral Monte Carlo simulations for rigid rotors and their application to water

In this work the path integral formulation for rigid rotors, proposed by Müser and Berne [Phys. Rev. Lett. {\bf 77}, 2638 (1996)], is described in detail. It is shown how this formulation can be used to perform Monte Carlo simulations of water. Our numerical results show that whereas some properties of water can be accurately reproduced using classical simulations with an empirical potential which, implicitly, includes quantum effects, other properties can only be described quantitatively when quantum effects are explicitly incorporated. In particular, quantum effects are extremely relevant when it comes to describing the equation of state of the ice phases at low temperatures, the structure of the ices at low temperatures, and the heat capacity of both liquid water and the ice phases. They also play a minor role in the relative stability of the ice phases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The stability of a crystal with diamond structure for patchy particles with tetrahedral symmetry

The phase diagram of model anisotropic particles with four attractive patches in a tetrahedral arrangement has been computed at two different values for the range of the potential, with the aim of investigating the conditions under which a diamond crystal can be formed. We find that the diamond phase is never stable for our longer-ranged potential. At low temperatures and pressures, the fluid freezes into a body-centred-cubic solid that can be viewed as two interpenetrating diamond lattices with a weak interaction between the two sublattices. Upon compression, an orientationally ordered face-centred-cubic crystal becomes more stable than the body-centred-cubic crystal, and at higher temperatures a plastic face-centered-cubic phase is stabilized by the increased entropy due to orientational disorder. A similar phase diagram is found for the shorter-ranged potential, but at low temperatures and pressures, we also find a region over which the diamond phase is thermodynamically favored over the body-centred-cubic phase. The higher vibrational entropy of the diamond structure with respect to the body-centred-cubic solid explains why it is stable even though the enthalpy of the latter phase is lower. Some preliminary studies on the growth of the diamond structure starting from a crystal seed were performed. Even though the diamond phase is never thermodynamically stable for the longer-ranged model, direct coexistence simulations of the interface between the fluid and the body-centred-cubic crystal and between the fluid and the diamond crystal show that, at sufficiently low pressures, it is quite probable that in both cases the solid grows into a diamond crystal, albeit involving some defects. These results highlight the importance of kinetic effects in the formation of diamond crystals in systems of patchy particles.

0 citations0 reviewsNo code linkedOpen access