Simple Landau Model of the Liquid-R IIR I Rotator Phases of Alkanes

Simple Landau free energy function is presented to describe the Liquid-RII-RI phase sequence of alkanes and transitions between them. The order parameters necessary to describe these rotator phase transitions are identified. We present a mean-field description of the Liquid-RII and Liquid-RI transitions. General arguments are presented for the topology of the phase diagram in the vicinity of the Liquid-RII-RI triple point. Within this model the Liquid-RII and Liquid-RI transitions are found to be always strongly first order. Calculations based on this model agree qualitatively with experiments.


Introduction
During the last two decades much progress has taken place in the field of rotator phases.Rotator phases are among the most interesting condensed states of matter exhibited by normal alkanes ( 2 2 n n ), alcohols, and other hydrocarbon chain systems consisting of layered structures with three dimensional crystalline order of the center of mass, but no long range orientational order of the molecules about their long axes.Rotator phases have a number of unique and unusual properties which include surface crystallization, anomalous heat capacity, negative thermal compressibilities and unusually high thermal expansions.Five different rotator phases had been identified.The rotator-II ( II ) phase is usually described as composed of molecules that are untilted with respect to the layers that are packed in a hexagonal lattice.The layers are stacked in an ABCABC... trilayer stacking sequence.This phase is also referred to as the rhombohedral phase.In the rotator-I ( I ) phase the molecules are also untilted with respect to the layers and there is a rectangularly distorted hexagonal lattice.The layers are stacked in an ABAB... bilayer stacking sequence.This phase is referred to as the face-centered-orthorhombic (FCO) phase.In shorter chain lengths the phase sequence was reported to be Liquid-R II -R I -Crystal.

C H R R
A large number of experimental studies [1][2][3][4][5][6][7][8][9][10] are devoted to describe the structure and the phase transitions of the rotator phases.According to the X-ray scattering study by Sirota et al. [3], the R II -R I transition is first order with jump of the distortion order parameter and sharp peak on the heat capacity data [4].The binary mixtures of normal alkanes [11][12][13][14] also shows a first order character of the R II -R I transition.Sirota et al. [15] carried out a high pressure study on the R II -R I transition and confirmed the first order character of the R II -R I transition.Zammit et al. [7,8] studied the IL-R II and R II -R I transitions in pure and binary mixtures of alkanes.Over the IL-R II transition, they observed the single peak in both the specific heat and latent heat in the pure material, splits into two features at different temperatures.This indicates the first order character of the IL-R II transition.They also confirmed the first order character of the R II -R I transitions.The presence of the thermal hysteresis at these transitions indicates the first order character of the transition.
Theoretical studies of the R II -R I transition follow two main lines.The first approach consists of Monte Carlo and molecular dynamics simulations [16][17][18][19][20][21][22] which confirms the first order character of the Liquid-R II (IL-R II ) and R II -R I transitions.The second approach is pursued by Wurger [23] and Mukherjee [24][25][26][27][28][29][30][31][32][33][34].Wurger [23] developed a microscopic model for the pair interaction of hydrocarbon chains and discussed the detailed structure of the R I and R II phases in terms of a molecular-field approximation.In a series of paper Mukherjee [24][25][26] discussed R II -R I phase transition within Landau phenomenological approach and discussed in detail the various aspect of this transition including the elastic properties.
To the best of the author's knowledge, there is so far no detailed theoretical studies on the IL-R II -R I phase sequence and the transitions between them.The purpose of the present paper is to investigate the IL-R II -R I phase sequence and the transitions between them within Landau theory.We define a new order parameter to describe the IL-R II and IL-R I phase transitions.

Model Order Parameters
The R I phase differs from the R II phase only in the distortion of the hexagonal lattice.Following [5,35] we define the lattice distortion parameter where a and b are the major and minor axes of an ellipse draws through the six nearest neighbors.The distortion  is defined with respect to a plane whose normal is parallel to the long molecular axes.= 0  for the R II phase.Thus we take  as an order parameter for the R II -R I transition.Now to define the order parameter of the IL-R II or IL-R I transitions.The low temperature crystal phase of n-alkanes is found to exist in layered structures.In these cases the molecules stay in some layer stacking and the probability that a single molecule is present (partially) in two simultaneous layers is almost zero.When the temperature rises further the stacking breaks, i.e. the molecules start occupying positions which are shared by more than one layer stacking.We may choose to represent alternate layer sequences with suffixes 1,  etc.The bilayer stacking can be represented as 1 2 1 2 1 2 and the trilayer stacking can be represented as 1 2 3 1 2 3 1 2 3 We represent the probability of the k-th molecule of the system to occupy any of the j-th layer sequence as j .For bilayer the only feasible cases are and 2 k , one of which is 1 and the other being 0 in perfect layer ordering.Now we define a correlation factor 2,3, , j  This is clearly 0 in ordered layer phase and non-zero in other phases with lower or higher order.
The rotator phases of n-alkanes are found to exist in bilayer and trilayer structures.If we consider the correlation factor for the highest possible layer structure in the alkanes, i.e. trilayer we find that for the k-th molecule It is 0 (or almost zero) in R II and R I phases which in nature are bilayer and trilayer respectively.k  does have a finite value in the liquid phase which can be calculated.In isotropic liquid phase the probability density of the k-th molecule is constant everywhere in space.So, layer is the total number of layers.Hence  can be expressed as Thus k  is independent of k.Hence we define the correlation order parameter  as Here the average correlation factor is calculated for all molecules together.
Thus  is 0 is the isotropic liquid phase and has a value >0 in the R I and R II phases.

Free Energy
Thus we take  and  as two order parameters involved in the R I -R II , IL-R II and IL-R I phase transitions.
For simplicity we neglect the weak interlayer interaction between the stacking layers in different rotator phases so that the problem becomes two dimensional.The distortion is a two component order parameter; its components are expressed the distortion amplitude  and the azimuthal angle 2 .The multiplier 2 comes from the fact that the distortion is a symmetric traceless tensor.Since the the free energy is a scalar quantity, negative and positive  and  result in inequivalent structures, causing the Landau free energy expansion contain 3  and 3   terms, thus resulting in a first order transition.Expanding the total free energy in terms of the above mentioned order parameters yields where 0 F is free energy of the isotropic liquid phase.
The coefficients a and  are assumed to vary strongly with an external parameter.T l and where 0 is some specific chain length of the molecules.u and v are positive constants.
By lowering the temperature from the isotropic liquid phase, the R II and R I phases can appear sequentially or in partial sequence.The R II and R I phases can arise either directly from the IL phase along the curves IL-R II and IL-R I or along the curve R II -R I .In the spirit of Landau theory all the phase transitions IL-R II , IL-R I and R II -R I are first order because of the cubic invariant in the free energy expansion.If all the phase transitions involved are first order ones, then one can observe the IL-R II -R I triple point as observed in experiment [3].
 2) R II phase: 0 and a < 0. The R I phase exists for Thus it is clear from the solutions that three types of transition are possible: 1) IL-R II ; 2) IL-R I ; 3) R II -R I .
The conditions for the first order R II -R I transition can be obtained as The sufficient condition for the R I phase to be stable are The conditions for the first order IL-R I transition are given by The conditions for the first order IL-R II transition read Solving ( 10)-( 12) simultaneously will determine the various phase transition lines.The region of the R I phase shrinks and finally disappears when the IL-R II transition takes place.In experimental studies [3], alkanes C20-C27 show such a IL-R II transition via triple point.Thus the above analysis of the IL-R II -R I triple point agrees well with the experimental observations [3].Thus there is always a direct IL-R II transition is possible within the framework of our model free energy (6)  .The so-called temperature hysteresis is related to the existence of metastable states within a certain temperature range.The above analysis qualitatively agrees with experimental observations [3,4,7,8].

Conclusion
A simple model free energy has been constructed to describe the IL-R II -R I phase sequence and transitions between them.The order parameters are identified for different phase transitions.The model predicts the first order character of the IL-R II , IL-R I and R II -R I transitions and IL-R II -R I triple point in the phase diagram.The proposed interpretation of the IL-R II transition allow us to explain the various types of phase behavior observed experimentally.These results are in qualitative agreement with all experiments reported so far.

3 
For = 0  , free energy (5) describe a first order R I -R II transition for b > 0 and c > 0. In this case the minimum free energy occurs at = 0  , for b > 0 and at = π for b < 0. According to the experimental observations, in the R I phase, = 0  .H is the coupling constant. ,  , b, c and H are chosen positive.The material parameters a and  can be assumed as the experimental phase diagrams[3] one observes,

Figure 1
shows a typical phase diagram for the IL-R I , IL-R II and R II -R I phase transitions.As can be seen from the Figure 1, the R II and R I phases arise from the isotropic phase along the curves IL-R II and IL-R I or along the curve R II -R I respectively.The IL-R II and IL-R I transitions are first order because of the cubic invariant in the free energy expansion.The line of the R II -R I transition starts at the IL-R II -R I triple point as shown in Figure 1.When the temperature of the IL-R II and of R II -R I transitions coincide, a triple point appears.These three conditions determine the stability of the R I phase explicitly.The sufficient condition for the stability of the R II phase reads phase the stability conditions are

Figure 1 .
Figure 1.Possible chain length (n)-temperature (T) phase diagram in the vicinity of the Liquid-R II -R I triple point.  0 0 = n l l l  .
satisfying the above stability conditions.Of course the IL-R I transition could proceed before the IL-R II the transition temperature II to be larger than the IL-R I transi-