A Simple Free Energy for the Isotropic-Nematic Phase Transition of Rods

. A free energy expression is proposed that describes the isotropic-nematic binodal concentrations of hard rods. A simple analytical form for this free energy was yet only available using a Gaussian trial function for the orientation distribution function (ODF), leading,however,toasignificantdeviationofthepredictedbinodals.Thenewfreeenergyproposedhereisbaseduponarationalized correctiontotheorientationalandpackingentropieswhenusingtheGaussianODF.IncombinationwithParsons-Leetheoryor scaledparticletheory,itenablesdescribingtheisotropic-nematicphasecoexistenceconcentrationsofrodsaccuratelyusingthe simpleGaussianODFforawiderangeofaspectratios.


Introduction
Liquid crystal phases are found in a wide variety of soft matter systems including solutions of amphiphilic molecules, polymer solutions, and colloidal dispersions (viruses, DNA, cellulose, boehmite, and vanadia) [1].These materials are composed of particles that can assume a preferred orientation making them suitable for applications due to special rheological and optical properties.Liquid crystal displays have been used since the 1980s and were later used in mobile phones, laptops, and flat screen televisions (see, e.g., [2] and references therein).
The basic understanding of the origin of liquid crystalline phases commenced thanks to theoretical work on the ability of hard rods to form a nematic, the simplest liquid crystal phase.In the 1940s, Onsager showed that an isotropic fluid of very thin, hard rods undergoes a transition to a nematic phase upon concentrating the rods [3,4].In the isotropic phase, the rods are oriented randomly whereas in the nematic phase they have a preferred, more parallel orientation.Onsager argued that the randomly oriented rods in the isotropic phase reserve less free volume per particle than the aligned rods in the nematic phase.Hence the entropy loss associated with the orientational ordering in the nematic phase is more than compensated by the reduction of excluded volume between pairs of rods [5].
For very thin, rigid, hard rods, the phase transition occurs at a very low volume fraction and the virial expansion may therefore be truncated at the level of the second (osmotic) virial coefficient, leading to an exact theory for infinitely thin particles [4][5][6].This rigorous demonstration of an entropydriven phase transition in hard particle dispersions received significant impetus [7].It stimulated computer simulation work that revealed highly ordered smectic phases in concentrated suspensions of hard rods [8].
In the following, a brief exposition of Onsager's theory for long rods is given in Section 2. Then the Gaussian approximation to the orientational distribution function is introduced in Section 3 that gives analytical results for the binodal concentrations.Subsequently, the new equation of state is proposed in Section 4 after which it is extended towards a description for arbitrary aspect ratio / using both Parsons theory The Helmholtz free energy  for a dispersion of  hard rods, modeled as spherocylinders with length  and diameter  in a volume , can be written in the second virial approximation (valid for  ≫ ) [5] as with Boltzmann constant  and temperature .The quantity  is the dimensionless concentration  = , where  = / is the number density of the rods and  = (/4) 2  is the second virial coefficient of a suspension of rods in the isotropic state.Note that (1) holds up to the binary collision limit [5,6] and is thus valid for long, thin hard rods only.
The quantity  is related to the orientational entropy loss per particle  or = − and  to the packing entropy loss per particle  pack = −.
The quantity [] can be computed from [5] with Ω being an infinitesimal surface element on the unit sphere.The orientational distribution function (ODF)  is normalized according to The quantity  can be computed from integration over the orientations between two rods: where Θ is the angle between the rods which depends on their orientations Ω and Ω  .It is noted here that  does not depend on the azimuthal angle  but only on the polar angle  (differing from Θ). Furthermore, the distribution function () must satisfy inversion symmetry, implying that the angles  and  −  are equivalent.
For the isotropic phase,  and  attain the values And hence the free energy in the isotropic phase follows as From ( 6) and (for the nematic phase) (1), standard thermodynamics provide the chemical potentials and pressures in each phase.Numerically solving the ODF () and simultaneously the coexistence equations yields the coexistence concentrations [5]: Note that the computation is involved since solving the coexistence concentrations includes numerical minimization of the ODF at each concentration of rods in the nematic phase; see [9] for an explanation of an efficient numerical approach.

Gaussian Orientational Distribution Function
It is quite useful for experimentalists and engineers to have analytical approximations at hand to estimate the concentrations at the isotropic-nematic (I-N) phase transition for practical purposes.Several trial functions (e.g., [4,6,10,11]) for the ODF have been proposed in order to avoid numerical minimization [9] and retain insight into the expressions.
Onsager [4] proposed a trial function that describes exact results quite reasonably but is numerically still quite demanding [11].
Odijk [6] realized that for rather aligned rods  can be approximated as a Gaussian distribution function: where the normalization factor Ñ() ∼ /4 for large .
When using the Gaussian ODF, the only parameter that needs to be determined from minimization of  is .The advantage of this Gaussian distribution function is that the quantities  and  can be calculated analytically for large values of  [5]: This leads to the following expression for the free energy in the nematic phase: Minimizing this expression with respect to  leads to Hence The coexistence equations now take the simple forms resulting in the following coexisting concentrations: The deviations from the numerical results   = 3.29 and   = 4.19 are considerable, being 5% for the isotropic and even 22% for the nematic phase coexistence concentrations.

Correction to the Free Energy
It should be realized that the Gaussian approximation to the ODF is an extremely rough approximation that only holds asymptotically for highly aligned rods.Within the Gaussian approach, the deviation between the results for the isotropicnematic coexistence concentrations (( 15) versus ( 8)) mainly arises from an ODF that is too sharply peaked.One can try to restore this by realizing that the orientational entropy loss  or per rod is overestimated due to a description of the rods that is too aligned.Further, the packing entropy gain is also not adequately described by the Gaussian ODF.Hence one might attempt to correct the orientational and packing entropy contributions with an unknown amount Δ: The value of Δ can be computed from the difference between exact and (Gaussian) approximated packing and orientational entropy: This is done at the numerically exact concentration   = 4.19.
For Δ = 0.139 the coexistence concentrations from (1) and ( 16) become The deviations from the exact Onsager results (see (8)) are now 1% for the isotropic phase and 10% for the nematic phase.This is a severe reduction of the deviation (a factor of 5 for the isotropic phase and a factor of 2 for the nematic phase) with respect to the uncorrected Gaussian approximation.

Aspect Ratio Dependence
The correction to the free energy upon using the Gaussian approximation for the ODF is now extended towards hard spherocylinders for arbitrary aspect ratio /.This means that interactions between the rods should be accounted for beyond the level of the second virial coefficient.Theoretically, at significant rod volume fractions, this is not possible without approximate descriptions.Here we use two approaches, scaled particle theory (SPT) and Parsons-Lee theory, the latter being used frequently by engineers.

Parsons-Lee
Theory.Parsons [13] derived an expression for the excess free energy of hard spherocylinders that can be considered as an extension of the Carnahan-Starling equation of state for hard spheres [14].The excess part of the Carnahan-Starling equation is multiplied by the ensemble-averaged excluded volume of a spherocylinder, normalized with its volume.This free energy expression was extended by Lee [15,16] to describe the isotropic-nematic phase transition of hard spherocylinders.The phase transition concentrations and osmotic pressure at coexistence, obtained after numerical minimization, have been shown to correspond quite reasonably with computer simulations [5].
For the isotropic phase, the Parsons result is where () is defined by and derives from the Carnahan-Starling equation of state, which in itself accurately describes the (osmotic) pressure of a collection of hard spheres in the fluid state up to a packing fraction of ≈50 vol%.Here, , the volume fraction of rods, which equals V  , where V  =  3 /6+ 2 /4, is the volume of a rod.Hence, using  = ,  and  are related via For the nematic phase, the Parsons result reads including the proposed correction Δ = 0.139.For  and , (10) from the Gaussian approximation are used and inserted into the free energy in (22).Minimization of this free energy with respect to  then leads to the following closed expression for the variational parameter: Now the isotropic-nematic coexistence concentrations can be obtained straightforwardly; the minimization of the ODF is incorporated into (23).
For the sake of completeness, we write down the analytical expressions for the osmotic pressure and chemical potential.The chemical potential  = (/) , follows from the Helmholtz energy  as  = ([/]/) , , and the osmotic pressure  = −(/) , can be found from V  =  2 ([/]/) , .The free energy of ( 19) yields the following results for the chemical potential of a spherocylinder in the isotropic phase: For the chemical potential of the rods in the nematic phase, (22) yields The result for the osmotic pressure is The results for the I-N binodals using Parsons-Lee theory for spherocylinders with a correction for the entropy are plotted in Figure 1.The new results are confronted with computer simulations on dispersions of rods to determine the rod concentrations at isotropic-nematic coexistence.These have been performed on relatively short rods first [18] but were performed for / up to 50 by Bolhuis and Frenkel [17] soon after.Their results are represented as the open symbols in Figure 1.The results are plotted in terms of the dimensionless (Onsager) concentration  = .It is clear that the simulation data are quite well described using the proposed new free energy expressions as a function of /, except for a slight overestimation of the concentration of rods in the coexisting nematic phase for long rods.

Scaled Particle Theory.
Scaled particle theory has the advantage that it also allows incorporating attractions in a direct manner [19] and also via, for instance, the contribution of nonadsorbing polymer chains [20].The SPT approach of Cotter [21] is used, deriving a SPT expression for the Helmholtz energy  of a dispersion of rods.For the isotropic phase, the result is where  and  are  involving the ODF via .Here  is the overall length-todiameter ratio / + 1.The quantity  is defined as  = /(1 − ).
For the nematic phase, the SPT result reads including the added proposed correction.Minimization of  with respect to  yields  = 36  ( − 1) The binodal isotropic-nematic concentrations can now be computed from ( 27)-(30).The scaled particle free energy of (27) yields the following results for the chemical potential of a spherocylinder in the isotropic phase: For the chemical potential in the nematic phase, (29) yields Advances in Condensed Matter Physics

SPT Gauss
Figure 2: Isotropic-nematic phase coexistence for hard spherocylinders as a function of the inverse of the aspect ratio /.Scaled particle theory results without (dashed curves) and with (solid curves) correction to the free energy are compared with Monte Carlo (MC) computer simulation data [17] and the Onsager limit.
The result for the osmotic pressure is where  and  in the above equations differ in the isotropic and nematic phase because of .It is noted that the above Parsons-Lee and SPT expressions for the free energy are semiempirical; third-and higher-order virial terms are included approximately [5].Implicitly, it is assumed that the correction needed to account for the Gauss ODF that is too sharply peaked in the Onsager limit similar when the aspect ratio / gets smaller.The results using scaled particle theory with a correction for the entropy for the I-N coexistence concentrations of spherocylinders are plotted in Figure 2.They are compared to SPT results using the Gaussian ODF without correction (dashed curves).Again we compare the results with computer simulations of [17], represented as the open symbols in Figure 2. Also when using SPT the simulation data are quite well described using the proposed corrected free energy expressions.Obviously, the free energy correction via the rationalized Δ correction now enables describing the I-N transitions quite well and is certainly a severe improvement to the (uncorrected) Gaussian approximation.

Concluding Remarks
A simple free energy correction was presented for the isotropic-nematic phase transition of hard spherocylinders.Basically, Onsager's free energy with a simple Gaussian orientational distribution function (ODF) was modified by realizing that this ODF is too sharply peaked.The proposed corrected free energy provides quite accurate predictions for the isotropic-nematic phase coexistence concentrations.
It can also be extended to include interactions beyond the pair level using Parsons-Lee or scaled particle theory.It is shown that inclusion of the proposed corrections into these simple approximate theories gives a fair description of the isotropic-nematic phase behavior for arbitrary rod length over diameter ratios.For more complicated applications of dispersions of rods beyond hard interactions, numerical approaches rapidly become more involved.Then this (still) simple, yet (more) accurate description of isotropic-nematic equilibria may become quite useful.

Figure 1 :
Figure1: Aspect ratio dependence of the isotropic-nematic phase coexistence for hard spherocylinders.Dashed dotted curves are predictions using the proposed free energy correction inserted into Parsons-Lee theory.Monte Carlo (MC) computer simulation data were taken from Bolhuis and Frenkel[17].