Source author record

N. G. Almarza

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

16 published item(s)

preprint2016arXiv

Lattice Model for water-solute mixtures

A lattice model for the study of mixtures of associating liquids is proposed. Solvent and solute are modeled by adapting the associating lattice gas (ALG) model. The nature of interaction solute/solvent is controlled by tuning the energy interactions between the patches of ALG model. We have studied three set of parameters, resulting on, hydrophilic, inert and hydrophobic interactions. Extensive Monte Carlo simulations were carried out and the behavior of pure components and the excess properties of the mixtures have been studied. The pure components: water (solvent) and solute, have quite similar phase diagrams, presenting: gas, low density liquid, and high density liquid phases. In the case of solute, the regions of coexistence are substantially reduced when compared with both the water and the standard ALG models. A numerical procedure has been developed in order to attain series of results at constant pressure from simulations of the lattice gas model in the grand canonical ensemble. The excess properties of the mixtures: volume and enthalpy as the function of the solute fraction have been studied for different interaction parameters of the model. Our model is able to reproduce qualitatively well the excess volume and enthalpy for different aqueous solutions. For the hydrophilic case, we show that the model is able to reproduce the excess volume and enthalpy of mixtures of small alcohols and amines. The inert case reproduces the behavior of large alcohols such as, propanol, butanol and pentanol. For last case (hydrophobic), the excess properties reproduce the behavior of ionic liquids in aqueous solution.

preprint2015arXiv

Effects of rigid or adaptive confinement on colloidal self-assembly. Fixed vs. fluctuating number of confined particles

The effects of confinement on colloidal self-assembly in the case of fixed number of confined particles are studied in the one dimensional lattice model solved exactly in the Grand Canonical Ensemble (GCE) in [J. Pękalski et al. J. Chem. Phys. 142, 014903 (2015)]. The model considers a pair interaction defined by a short-range attraction plus a longer-range repulsion. We consider thermodynamic states corresponding to self-assembly into clusters. Both, fixed and adaptive boundaries are studied. For fixed boundaries, there are particular states in which, for equal average densities, the number of clusters in the GCE is larger than in the Canonical Ensemble. The dependence of pressure on density has a different form when the system size changes with fixed number of particles and when the number of particles changes with fixed size of the system. In the former case the pressure has a nonmonotonic dependence on the system size. The anomalous increase of pressure for expanding system is accompanied by formation of a larger number of smaller clusters. In the case of elastic confining surfaces we observe a bistability, i.e. two significantly different system sizes occur with almost the same probability. The mechanism of the bistability in the closed system is different to that of the case of permeable walls, where the two equilibrium system sizes correspond to a different number of particles.

preprint2015arXiv

Pattern formation in binary fluid mixtures induced by short-range competing interactions

Molecular dynamics simulations and integral equation calculations of a simple equimolar mixture of diatomic molecules and monomers interacting via attractive and repulsive short-range potentials show the existence of pattern formation (microheterogeneity), mostly due to depletion forces away from the demixing region. Effective site-site potentials extracted from the pair correlation functions using an inverse Monte Carlo approach and an integral equation inversion procedure exhibit the features characteristic of a short-range attractive and long-range repulsive potential. When charges are incorporated into the model, this becomes a coarse grained representation of a room temperature ionic liquid, and as expected, intermediate range order becomes more pronounced and stable.

preprint2014arXiv

Demixing and confinement in slit pores

Using Monte Carlo simulation, we study the influence of geometric confinement on demixing for a series of symmetric non-additive hard spheres mixtures confined in slit pores. We consider both a wide range of positive non-additivities and a series of pore widths, ranging from the pure two dimensional limit to a large pore width where results are close to the bulk three dimensional case. Critical parameters are extracted by means of finite size analysis. We find that for this particular case in which demixing is induced by volume effects, phase separation is in most cases somewhat impeded by spatial confinement. However, a non-monotonous dependence of the critical pressure and density with pore size is found for small non-additivities. In this latter case, it turns out that an otherwise stable bulk mixture can be forced to demix by simple geometric confinement when the pore width decreases down to approximately one and a half molecular diameters.

preprint2014arXiv

Periodic ordering of clusters and stripes in a two-dimensional lattice model. II. Results of Monte Carlo simulation

The triangular lattice model with nearest-neighbor attraction and third-neighbor repulsion, introduced in [J. Pekalski, A. Ciach and N. G. Almarza, arXiv:1401.0801 [cond-mat.soft]] is studied by Monte Carlo simulation. Introduction of appropriate order parameters allowed us to construct a phase diagram, where different phases with patterns made of clusters, bubbles or stripes are thermodynamically stable. We observe, in particular, two distinct lamellar phases - the less ordered one with global orientational order and the more ordered one with both orientational and translational order. Our results concern spontaneous pattern formation on solid surfaces, fluid interfaces or membranes that is driven by competing interactions between adsorbing particles or molecules.

preprint2014arXiv

Three-dimensional patchy lattice model: ring formation and phase separation

We investigate the structural and thermodynamic properties of a model of particles with $2$ patches of type $A$ and $10$ patches of type $B$. Particles are placed on the sites of a face centered cubic lattice with the patches oriented along the nearest neighbor directions. The competition between the self-assembly of chains, rings and networks on the phase diagram is investigated by carrying out a systematic investigation of this class of models, using an extension of Wertheim's theory for associating fluids and Monte Carlo numerical simulations. We varied the ratio $r\equivε_{AB}/ε_{AA}$ of the interaction between patches $A$ and $B$, $ε_{AB}$, and between $A$ patches, $ε_{AA}$ ($ε_{BB}$ is set to $0$) as well as the relative position of the $A$ patches, i.e., the angle $θ$ between the (lattice) directions of the $A$ patches. We found that both $r$ and $θ$ ($60^\circ,90^\circ,$ or $120^\circ$) have a profound effect on the phase diagram. In the empty fluid regime ($r < 1/2$) the phase diagram is re-entrant with a closed miscibility loop. The region around the lower critical point exhibits unusual structural and thermodynamic behavior determined by the presence of relatively short rings. The agreement between the results of theory and simulation is excellent for $θ=120^\circ$ but deteriorates as $θ$ decreases, revealing the need for new theoretical approaches to describe the structure and thermodynamics of systems dominated by small rings.

preprint2013arXiv

Periodic ordering of clusters in a one-dimensional lattice model

A generic lattice model for systems containing particles interacting with short-range attraction long-range repulsion (SALR) potential that can be solved exactly in one dimension is introduced. We assume attraction J_1 between the first neighbors and repulsion J_2 between the third neighbors. The ground state of the model shows existence of two homogeneous phases (gas and liquid) for J_2/J_1<1/3. In addition to the homogeneous phases, the third phase with periodically distributed clusters appears for J_2/J_1>1/3. Phase diagrams obtained in the self-consistent mean-field approximation for a range of values of J_2/J_1 show very rich behavior, including reentrant melting, and coexistence of two periodic phases (one with strong and the other one with weak order) terminated at a critical point. We present exact solutions for the equation of state as well as for the correlation function for characteristic values of J_2/J_1. Based on the exact results, for J_2/J_1>1/3 we predict pseudo-phase transitions to the ordered cluster phase indicated by a rapid change of density for a very narrow range of pressure, and by a very large correlation length for thermodynamic states where the periodic phase is stable in mean field. For 1/9<J_2/J_1<1/3 the correlation function decays monotonically below certain temperature, whereas above this temperature exponentially damped oscillatory behavior is obtained. Thus, even though macroscopic phase separation is energetically favored and appears for weak repulsion at T=0, local spatial inhomogeneities appear for finite T. Monte Carlo simulations in canonical ensemble show that specific heat has a maximum for low density ρthat we associate with formation of living clusters, and if the repulsion is strong, another maximum for ρ= 1/2.

preprint2013arXiv

Phase behaviour of the confined lattice gas Lebwohl-Lasher model

The phase behaviour of the Lebwohl-Lasher lattice gas model (one of the simplest representations of a nematogenic fluid) confined in a slab is investigated by means of extensive Monte Carlo simulations. The model is known to yield a first order gas-liquid transition in both the 2D and 3D limits, that is coupled with an orientational order-disorder transition. This latter transition happens to be first order in the 3D limit and it shares some characteristic features with the continuous defect mediated Berezinskii-Kosterlitz-Thouless transition in 2D. In this work we will analyze in detail the behaviour of this system taking full advantage of the lattice nature of the model and the particular symmetry of the interaction potential, which allows for the use of efficient cluster algorithms.

preprint2012arXiv

Closed-loop liquid-vapor equilibrium in a one-component system

We report Monte Carlo simulations that show closed-loop liquid-vapor equilibrium in a pure substance. As far as we know, this is the first time that such a topology of the phase diagram has been found for one-component systems. This finding has been achieved on a two-dimensional lattice model for patchy particles that can form network fluids. We have considered related models with a slightly different patch distribution in the order to understand the features of the distribution of patches on the surface of the particles that make possible the presence of the closed-loop liquid-vapor equilibrium, and its relation with the phase diagram containing the so-called empty liquids. Finally we discuss the likelihood of finding the closed-loop liquid-vapor equilibria on related models for three dimensional models of patchy particles in the continuum, and speculate on the possible relationship between the mechanism behind the closed-loop liquid vapor equilibrium of our simple lattice model and the salt-induced reentrant condensation found in complex systems.

preprint2012arXiv

Reply to Comment on Effect of polydispersity on the ordering transition of adsorbed self-assembled rigid rods

We comment on the nature of the ordering transition of a model of equilibrium polydisperse rigid rods, on the square lattice, which is reported by Lopez et al. to exhibit random percolation criticality in the canonical ensemble, in sharp contrast to (i) our results of Ising criticality for the same model in the grand canonical ensemble [Phys. Rev. E 82, 061117 (2010)] and (ii) the absence of exponent(s) renormalization for constrained systems with logarithmic specific heat anomalies predicted on very general grounds by Fisher [M.E. Fisher, Phys. Rev. 176, 257 (1968)].

preprint2012arXiv

The nature of the ordered phase of the confined self-assembled rigid rod model

We investigate the nature of the ordered phase and the orientational correlations between adjacent layers of the confined three-dimensional self-assembled rigid rod model, on the cubic lattice. We find that the ordered phase at finite temperatures becomes uniaxial in the thermodynamic limit, by contrast to the ground state (partial) order where the orientation of the uncorrelated layers is perpendicular to one of the three lattice directions. The increase of the orientational correlation between layers as the number of layers increases suggests that the unconfined model may also exhibit uniaxial ordering at finite temperatures.

preprint2012arXiv

Three-dimensional patchy lattice model for empty fluids

The phase diagram of a simple model with two patches of type A and ten patches of type B (2A10B) on the face centred cubic lattice has been calculated by simulations and theory. Assuming that there is no interaction between the B patches the behavior of the system can be described in terms of the ratio of the AB and AA interactions, r. Our results show that, similarly to what happens for related off-lattice and two-dimensional lattice models, the liquid-vapor phase equilibria exhibits reentrant behavior for some values of the interaction parameters. However, for the model studied here the liquid-vapor phase equilibria occurs for values of r lower than 1/3, a threshold value which was previously thought to be universal for 2AnB models. In addition, the theory predicts that below r = 1/3 (and above a new condensation threshold which is < 1/3) the reentrant liquid-vapor equilibria is so extreme that it exhibits a closed loop with a lower critical point, a very unusual behavior in single-component systems. An order-disorder transition is also observed at higher densities than the liquid-vapor equilibria, which shows that the liquid-vapor reentrancy occurs in an equilibrium region of the phase diagram. These findings may have implications in the understanding of the condensation of dipolar hard spheres given the analogy between that system and the 2AnB models considered here.

preprint2011arXiv

The condensation and ordering of models of empty liquids

We consider a simple model consisting of particles with four bonding sites ("patches"), two of type A and two of type B, on the square lattice, and investigate its global phase behavior by simulations and theory. We set the interaction between B patches to zero and calculate the phase diagram as the ratio between the AB and the AA interactions, ε^\star_AB, varies. In line with previous work, on three-dimensional off-lattice models, we show that the liquid-vapor phase diagram exhibits a re-entrant or "pinched" shape for the same range of ε^\star_AB, suggesting that the ratio of the energy scales - and the corresponding empty fluid regime - is independent of the dimensionality of the system and of the lattice structure. In addition, the model exhibits an order-disorder transition that is ferromagnetic in the re-entrant regime. The use of low-dimensional lattice models allows the simulation of sufficiently large systems to establish the nature of the liquid-vapor critical points and to describe the structure of the liquid phase in the empty fluid regime, where the size of the "voids" increases as the temperature decreases. We have found that the liquid-vapor critical point is in the 2D Ising universality class, with a scaling region that decreases rapidly as the temperature decreases. The results of simulations and theoretical analysis suggest that the line of order-disorder transitions intersects the condensation line at a multicritical point at zero temperature and density, for patchy particle models with a re-entrant, empty fluid, regime.

preprint2011arXiv

The criticality of self-assembled rigid rods on triangular lattices

The criticality of self-assembled rigid rods on triangular lattices is investigated using Monte Carlo simulation. We find a continuous transition between an ordered phase, where the rods are oriented along one of the three (equivalent) lattice directions, and a disordered one. We conclude that equilibrium polydispersity of the rod lengths does not affect the critical behavior, as we found that the criticality is the same as that of monodisperse rods on the same lattice, in contrast with the results of recently published work on similar models.

preprint2011arXiv

Theory and simulation of the confined Lebwohl-Lasher model

We discuss the Lebwohl-Lasher model of nematic liquid crystals in a confined geometry, using Monte Carlo simulation and mean-field theory. A film of material is sandwiched between two planar, parallel plates that couple to the adjacent spins via a surface strength $ε_s$. We consider the cases where the favoured alignments at the two walls are the same (symmetric cell) or different (asymmetric or hybrid cell). In the latter case, we demonstrate the existence of a {\it single} phase transition in the slab for all values of the cell thickness. This transition has been observed before in the regime of narrow cells, where the two structures involved correspond to different arrangements of the nematic director. By studying wider cells, we show that the transition is in fact the usual isotropic-to-nematic (capillary) transition under confinement in the case of antagonistic surface forces. We show results for a wide range of values of film thickness, and discuss the phenomenology using a mean-field model.

preprint2010arXiv

The effect of polydispersity on the ordering transition of adsorbed self-assembled rigid rods

Extensive Monte Carlo simulations were carried out to investigate the nature of the ordering transition of a model of adsorbed self-assembled rigid rods on the bonds of a square lattice [Tavares et. al., Phys. Rev E 79, 021505 (2009)]. The polydisperse rods undergo a continuous ordering transition that is found to be in the two-dimensional Ising universality class, as in models where the rods are monodisperse. This finding is in sharp contrast with the recent claim that equilibrium polydispersity changes the nature of the phase transition in this class of models [L`opez et. al., Phys. Rev E 80, 040105(R)(2009)].