Lattice solution model for order-disorder transitions in membranes and Langmuir monolayers

Generalized van der Waals theory for phase behavior of two-dimensional nematic liquid crystals: Phase ordering and the equation of state

Liquid crystalline ordering of anisotropic particles in two dimensions is important in many physical and biological systems and their phase behavior is still a topic of interest. A generalized van der Waals theory is formulated, accounting for repulsive excluded volume and attractive van der Waals and Maier-Saupe interactions, for rectangles confined to two dimensions. The phase ordering transitions and equation of state are analyzed as a function of the model parameters (aspect ratioL/B and isotropic and anisotropic interaction parameters χ and ν). Different phase transitions are observed: continuous isotropic-nematic (high L/B and ν), first-order isotropic-nematic (intermediate L/B and small ν), and continuous isotropic-tetratic (small L/B and ν) followed by a continuous tetratic-nematic transition at higher densities. Increasing L/B decreases the pressure, and this effect is more pronounced in the nematic than in the isotropic phase. Increasing both interaction parameters decreases pressure and can lead to phase separation.

Phase Transitions in Phospholipid Monolayers: Theory Versus Experiments

The Doniach lattice gas (DLG) represents a ternary-mixture statistical model, whose components, water molecules (w), ordered-chain lipids (o), and disordered-chain lipids (d)-the latter carrying a high degenerescence ω ≫ 1-are located at each site of a two-dimensional lattice. The DLG model was introduced to describe phospholipid Langmuir films at the air-water interface and can be mapped into a spin-1 model, with the single-site states s _i = 0, +1, and -1 representing the three types of molecules in the system (w, o, and d), respectively. The model allows lipid-density fluctuations and has been analyzed at the mean-field approximation (Guidi, H. S.; Henriques, V. B. Phys. Rev. E 2014, 90, 052705) as well as at the pair approximation (de Oliveira, F. O.; Tamashiro, M. N. Phys. Rev. E 2019, 99, 012147). In this work, we focus on performing an explicit comparison of the theoretical predictions obtained for the DLG model at the pair approximation with isothermal monolayer compression experiments (Nielsen, L. K.; Bjørnholm, T.; Mouritsen, O. G. Langmuir 2007, 23, 11684) for the two most commonly studied saturated zwitterionic phospholipids, DMPC (1,2-dimyristoyl- sn-glycero-3-phosphocholine) and DPPC (1,2-dipalmitoyl- sn-glycero-3-phosphocholine). The model parameters obtained by fitting to the experimental data yield phase diagrams that are qualitatively consistent with the observed phase transitions on DMPC and DPPC monolayers, with the absence of a low-density gas phase. Quantitative agreement, however, was less significant partially because of the challenging reproducibility of Langmuir monolayer compression experiments, claimed in the literature to be influenced by kinetic effects.

Phase transitions in phospholipid monolayers: Statistical model at the pair approximation

A Langmuir film, consisting of a phospholipid monolayer at the air-water interface, was modeled as a two-dimensional lattice gas corresponding to a ternary mixture of water molecules (w), ordered-chain lipids (o), and disordered-chain lipids (d). The statistical problem is formulated in terms of a spin-1 model, in which the disordered-chain lipid states possess a high degenerescence ω≫1, and was termed Doniach lattice gas (DLG). Motivated by some open questions in the analysis of the DLG model at the mean-field approximation (MFA) [Phys. Rev. E 90, 052705 (2014)PLEEE81539-375510.1103/PhysRevE.90.052705], we have reconsidered it at the pair-approximation level by solving the model on a Cayley tree of coordination z. The attractors of the corresponding discrete-map problem are associated with the thermodynamic solutions on the Bethe lattice (the central region of an asymptotically infinite Cayley tree). To check the thermodynamic stability of the possible phases, the grand-potential density was obtained by the method proposed by Gujrati [Phys. Rev. Lett. 74, 809 (1995)PRLTAO0031-900710.1103/PhysRevLett.74.809]. In general, the previous MFA results are confirmed at the pair-approximation level, but a novel staggered phase, overlooked in the MFA analysis, was found when the condition ε_{wd}>1/2(ε_{ww}+ε_{dd}) is satisfied, where ε_{xy} represents the nearest-neighbor intermolecular interactions between single-site states x and y. Model parameters obtained by fitting to experimental data for the two most commonly studied zwitterionic phospholipids, 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), yield phase diagrams consistent with the phase transitions observed on Langmuir films of the same lipids under isothermal compression, which present a liquid-condensed to a liquid-expanded first-order transition line ending at a critical point.

Multicritical behavior of the ferromagnetic Blume-Emery-Griffiths model with repulsive biquadratic couplings

The ferromagnetic (J>0) version of the Blume-Emery-Griffiths model in the region of repulsive biquadratic couplings (K<0) is considered on a Cayley tree of coordination z, reducing the statistical problem to the analysis of a two-dimensional nonlinear discrete map. In order to investigate the effect of the coordination z on the system multicritical behavior, we study the particular case K/J=-3.5 with the inclusion of crystal fields (D≠0), but vanishing external magnetic fields (H=0), for two distinct lattice coordinations (z=4 and z=6). The thermodynamic solutions on the Bethe lattice (the central region of a large Cayley tree) are associated with the attractors of the two-dimensional map. The phase diagrams display several thermodynamic phases (paramagnetic, ferromagnetic, ferrimagnetic, and staggered quadrupolar). In some cases, there are regions of numerical costability of two different attractors of the map, associated with discontinuous phase transitions between the corresponding phases. To verify the thermodynamic stability of the phases and to locate the first-order boundaries, the analytical expression of the Gibbs free energy was obtained by the method proposed by Gujrati [Phys. Rev. Lett. 74, 809 (1995)PRLTAO0031-900710.1103/PhysRevLett.74.809]. For lower coordinations (z=4) the transition between the ferrimagnetic and the staggered quadrupolar phases is always continuous, while the transition between the ferromagnetic and the ferrimagnetic phases is discontinuous at low temperatures, turning into continuous for temperatures above a tricritical point. On the other hand, for higher coordinations (z=6), the transition between the ferromagnetic and the ferrimagnetic phases is always continuous. However, the transition between the ferrimagnetic and the staggered quadrupolar phases is continuous for higher temperatures and discontinuous for temperatures below a tricritical point, in agreement with previous results obtained in the mean-field approximation (infinity-coordination limit). In both cases, the occurrence and the thermodynamic stability of the ferrimagnetic phase is confirmed.

Critical phenomena on k-booklets

We define a "k-booklet" to be a set of k semi-infinite planes with -∞<x<∞ and y≥0, glued together at the edges (the "spine") y=0. On such booklets we study three critical phenomena: self-avoiding random walks, the Ising model, and percolation. For k=2, a booklet is equivalent to a single infinite lattice, and for k=1 to a semi-infinite lattice. In both these cases the systems show standard critical phenomena. This is not so for k≥3. Self-avoiding walks starting at y=0 show a first-order transition at a shifted critical point, with no power-behaved scaling laws. The Ising model and percolation show hybrid transitions, i.e., the scaling laws of the standard models coexist with discontinuities of the order parameter at y≈0, and the critical points are not shifted. In the case of the Ising model, ergodicity is already broken at T=T_{c}, and not only for T<T_{c} as in the standard geometry. In all three models, correlations (as measured by walk and cluster shapes) are highly anisotropic for small y.