TIP 4 P 2005 Ice : Simulating water with two molecular states.

Rigid, non-polarizable water models are very efficient from a computational point of view, and some of them have a great ability in predicting experimental properties. There is, however, little room for improvement in simulating water with this strategy, whose main shortcoming is that water molecules do not change their interaction parameters in response to the local molecular landscape. In this work, we propose a novel modeling strategy that involves using two rigid non-polarizable models as states that water molecules can adopt depending on their molecular environment.

Effects of Terpenes on the Osteoarthritis Cytokine Profile by Modulation of IL-6: Double Face versus Dark Knight?

Background: Hemp seed oil and terpenes are emerging as a dietary supplement and complementary therapy for patients suffering from knee osteoarthritis (KOA). However, the mechanisms and effects induced by these molecules on inflammatory cytokines are not yet fully understood. The aim of this study was to evaluate the changes in the cytokine IL-1β, IL-1α, IL-2, IL-6, and TNF-α levels from two oral hemp seed oil-based dietary supplements, of which only one included the addition of terpenes, in a population of KOA patients.

Structural characterization of systems with competing interactions confined in narrow spherical shells.

Systems with short-range attraction and long-range repulsion can form ordered microphases in bulk and under confinement. In fact, confinement has been proven to be a good strategy to induce the formation of novel ordered microphases that might be appealing to the development of functional nanomaterials. Using Grand Canonical Monte Carlo (GCMC) simulations, we study a model colloidal system with competing interactions under confinement in narrow spherical shells at thermodynamic conditions under which the hexagonal phase is stable in bulk.

Phases of surface-confined trivalent colloidal particles.

Patchy colloids promise the design and modelling of complex materials, but the realization of equilibrium patchy particle structures remains challenging. Here, we assemble pseudo-trivalent particles and elucidate their phase behaviour when confined to a plane. We observe the honeycomb phase, as well as more complex amorphous network and triangular phases.

Estimation of bubble cavitation rates in a symmetrical Lennard-Jones mixture by NVT seeding simulations.

The liquid-vapor transition starts with the formation of a sufficiently large bubble in the metastable liquid to trigger the phase transition. Understanding this process is of fundamental and practical interest, but its study is challenging because it occurs over timescales that are too short for experiments but too long for simulations. The seeding method estimates cavitation rates by simulating a liquid in which a bubble is inserted, thus avoiding the long times needed for its formation.

Homogeneous nucleation rate of methane hydrate formation under experimental conditions from seeding simulations.

In this work, we shall estimate via computer simulations the homogeneous nucleation rate for the methane hydrate at 400 bars for a supercooling of about 35 K. The TIP4P/ICE model and a Lennard-Jones center were used for water and methane, respectively. To estimate the nucleation rate, the seeding technique was employed.

Visualizing defect dynamics by assembling the colloidal graphene lattice.

Graphene has been under intense scientific interest because of its remarkable optical, mechanical and electronic properties. Its honeycomb structure makes it an archetypical two-dimensional material exhibiting a photonic and phononic band gap with topologically protected states. Here, we assemble colloidal graphene, the analogue of atomic graphene using pseudo-trivalent patchy particles, allowing particle-scale insight into crystal growth and defect dynamics.

Self-Assembly of Optimally Packed Cylindrical Clusters inside Spherical Shells.

Systems with short-range attraction and long-range repulsion can form ordered microphases in bulk and under confinement. Using grand canonical Monte Carlo simulations, we study a colloidal system with competing interactions under confinement in narrow spherical shells at thermodynamic conditions at which the hexagonal phase of cylindrical clusters is stable in bulk. We observe spontaneous formation of different ordered structures.

Maximum in density of electrolyte solutions: Learning about ion-water interactions and testing the Madrid-2019 force field.

In this work, we studied the effect of Li, Na, K, Mg, and Ca chlorides and sulfates on the temperature of maximum density (TMD) of aqueous solutions at room pressure. Experiments at 1 molal salt concentration were carried out to determine the TMD of these solutions. We also performed molecular dynamics simulations to estimate the TMD at 1 and 2 m with the Madrid-2019 force field, which uses the TIP4P/2005 water model and scaled charges for the ions, finding an excellent agreement between experiment and simulation.

Freezing point depression of salt aqueous solutions using the Madrid-2019 model.

Salt aqueous solutions are relevant in many fields, ranging from biological systems to seawater. Thus, the availability of a force-field that is able to reproduce the thermodynamic and dynamic behavior of salt aqueous solutions would be of great interest. Unfortunately, this has been proven challenging, and most of the existing force-fields fail to reproduce much of their behavior.

Phase diagram of the NaCl-water system from computer simulations.

NaCl aqueous solutions are ubiquitous. They can crystallize into ice, NaCl, or NaCl · 2HO depending on the temperature-concentration conditions. These crystallization transitions have important implications in geology, cryopreservation, or atmospheric science.

Structural and Dynamical Behaviour of Colloids with Competing Interactions Confined in Slit Pores.

Systems with short-range attractive and long-range repulsive interactions can form periodic modulated phases at low temperatures, such as cluster-crystal, hexagonal, lamellar and bicontinuous gyroid phases. These periodic microphases should be stable regardless of the physical origin of the interactions. However, they have not yet been experimentally observed in colloidal systems, where, in principle, the interactions can be tuned by modifying the colloidal solution.

How to design an icosahedral quasicrystal through directional bonding.

Icosahedral quasicrystals (IQCs) are materials that exhibit long-range order but lack periodicity in any direction. Although IQCs were the first reported quasicrystals, they have been experimentally observed only in metallic alloys, not in other materials. By contrast, quasicrystals with other symmetries (particularly dodecagonal) have now been found in several soft-matter systems.

Programming patchy particles to form three-dimensional dodecagonal quasicrystals.

Model patchy particles have been shown to be able to form a wide variety of structures, including symmetric clusters, complex crystals, and even two-dimensional quasicrystals. Here, we investigate whether we can design patchy particles that form three-dimensional quasicrystals, in particular targeting a quasicrystal with dodecagonal symmetry that is made up of stacks of two-dimensional quasicrystalline layers. We obtain two designs that are able to form such a dodecagonal quasicrystal in annealing simulations.

Formation and internal ordering of periodic microphases in colloidal models with competing interactions.

Theory and simulations predict that colloidal particles with short-range attractive and long-range repulsive interactions form periodic microphases if there is a proper balance between the attractive and repulsive contributions. However, the experimental identification of such structures has remained elusive to date. Using molecular dynamics simulations, we investigate the phase behaviour of a model system that stabilizes a cluster-crystal, a cylindrical and a lamellar phase at low temperatures.

How ice grows from premelting films and water droplets.

Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes.

The Young-Laplace equation for a solid-liquid interface.

The application of the Young-Laplace equation to a solid-liquid interface is considered. Computer simulations show that the pressure inside a solid cluster of hard spheres is smaller than the external pressure of the liquid (both for small and large clusters). This would suggest a negative value for the interfacial free energy.

Confinement of Colloids with Competing Interactions in Ordered Porous Materials.

In this work, we explore the possibility of promoting the formation of ordered microphases by confinement of colloids with competing interactions in ordered porous materials. For that aim, we consider three families of porous materials modeled as cubic primitive, diamond, and gyroid bicontinuous phases. The structure of the confined colloids is investigated by means of grand canonical Monte Carlo simulations in thermodynamic conditions at which either a cluster crystal or a cylindrical phase is stable in bulk.

Surface phase transitions and crystal habits of ice in the atmosphere.

With climate modeling predicting a raise of at least 2°C by year 2100, the fate of ice has become a serious concern, but we still do not understand how ice grows (or melts). In the atmosphere, crystal growth rates of basal and prism facets exhibit an enigmatic temperature dependence and crossover up to three times in a range between 0° and -40°. Here, we use large-scale computer simulations to characterize the ice surface and identify a sequence of previously unidentified phase transitions on the main facets of ice crystallites.

Micelle-directed chiral seeded growth on anisotropic gold nanocrystals.

Surfactant-assisted seeded growth of metal nanoparticles (NPs) can be engineered to produce anisotropic gold nanocrystals with high chiroptical activity through the templating effect of chiral micelles formed in the presence of dissymmetric cosurfactants. Mixed micelles adsorb on gold nanorods, forming quasihelical patterns that direct seeded growth into NPs with pronounced morphological and optical handedness. Sharp chiral wrinkles lead to chiral plasmon modes with high dissymmetry factors (~0.

Rounded Layering Transitions on the Surface of Ice.

Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers.

Programming patchy particles to form complex periodic structures.

We introduce a scheme to design patchy particles so that a given target crystal is the global free-energy minimum at sufficiently low temperature. A key feature is a torsional component to the potential that only allows binding when particles have the correct relative orientations. In all examples studied, the target crystal structures readily assembled on annealing from a low-density fluid phase, albeit with the simpler target structures assembling more rapidly.

The influence of confinement on the structure of colloidal systems with competing interactions.

Colloidal systems with competing interactions have a complex phase diagram with several periodic microphases, in which particles are arranged in lamellae, cylinders or clusters. Using grand canonical Monte Carlo simulations, we investigate how the structure of the colloidal fluid can be modified by confinement in channels with different cross-section geometries and sizes. We pay particular attention to the hexagonal cylindrical phase since it is the most susceptible to form new structures from.