Periodic ordering of clusters in a one-dimensional lattice model

TL;DR

This paper introduces an exactly solvable one-dimensional lattice model with SALR interactions, revealing complex phase behavior including periodic clustering, reentrant melting, and pseudo-phase transitions, supported by analytical and Monte Carlo results.

Contribution

It provides an exact solution for a 1D SALR lattice model, uncovering rich phase phenomena and detailed correlation functions not previously characterized.

Findings

Existence of two homogeneous phases (gas and liquid) for J_2/J_1<1/3.

Appearance of a periodic cluster phase for J_2/J_1>1/3.

Prediction of pseudo-phase transitions with large correlation lengths.

Abstract

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 \rho that we associate with formation of living clusters, and if the repulsion is strong, another maximum for \rho = 1/2.