Article Entropic Behavior of Binary Carbonaceous Mesophases

Abstract: The Maier-Saupe model for binary mixtures of uniaxial discotic nematogens, formulated in a previous study [1], is used to compute and characterize orientational entropy [2] and orientational specific heat. These thermodynamic quantities are used to determine mixture type (ideal or non-ideal) which arise due to their different intrinsic properties, determined by the molecular weight asymmetry ΔMw and the molecular interaction parameter β. These molecular properties are also used to characterize the critical concentration where the mixture behaves like a single component system and exhibits the minimum nematic to isotropic (NI) transition temperature (pseudo-pure mixture). A transition within the nematic phase takes place at this specific concentration. According to the Maier-Saupe model, in a single mesogen, entropy at NI transition is a universal value; in this work we quantify the mixing effect on this universal property. The results and analysis provide a new tool to characterize molecular interaction and molecular weight differences in mesogenic mixtures using standard calorimetric measurements.


Introduction
An ideal crystalline solid (Figure 1a) has: (i) perfect orientational order as well as (ii), perfect positional order.On the other hand, an isotropic liquid (Figure 1c) lacks any kind of ordering, either OPEN ACCESS orientational or positional.Between these two extremes lies a nematic liquid crystal that has a partially orientational order without any positional order (Figure 1b) [3].In nematic liquid crystals [4] (NLCs), the long molecular axes (u direction in Figure 2) are preferably oriented along a particular direction called the director n (Figure 2).In discotic nematic liquid crystals (DNLCs) the director n is perpendicular to the long axis of discotic molecules (Figure 2).In rod-like nematic liquid crystals (RNLCs) n is parallel to the long axis of the molecules (Figure 2).Carbonaceous mesophases (CMs), first reported by Brooks and Taylor [5], are DNLC mixtures obtained from petroleum pitches and synthetic naphthalene precursors [6].The composition, polydispersity, and molecular orientation of CMs play a significant role on the final properties of cokes [7], carbon foams, carbon/carbon composites [8], and carbon fibers [9][10][11][12].
Figure 2. Schematic of nematic liquid crystals (NLCs), the director n is the average of the molecular orientation u, and its classification into discotic and rod-like molecules.In discotic nematic liquid crystals the director n (average orientation) is perpendicular to the long axis of the molecules; however, in rod-like nematic liquid crystals the director n is parallel to the long axis of the molecules.The type and extent of NLC applications [13] are mainly determined by their thermal properties.For instance, it is known that the final structure of the fibers based on the pure liquid crystalline materials is influenced by the temperature [14][15][16][17].The results of our previous study [1] shows that the structure of the fibers based on the real CMs which are composed of DNLC species with different molecular weights and concentrations are also influenced by the thermal effect.
By the use of thermodynamics, microscopic structure as well as macroscopic energetic and/or entropic aspects is revealed, so that comprehensive understanding of materials can be achieved [18].One of the methods mostly used for this purpose is calorimetry.Calorimetric measurements are best suited for the recognition of phase transitions and the determination of their transition temperature.Specific heat, which can be measured by calorimetry, is one of the most useful quantities to investigate the thermal properties of LCs [18].By use of this quantity three fundamental thermodynamic values viz.enthalpy, entropy, and free energy can be obtained; as a result, it can give an insight into the microscopic and macroscopic aspects of the system.
Theoretically, the entropic behavior and specific heat of pure NLCs are well understood by the Maier-Saupe theory [2].For instance, this theory predicts a universal entropy value at nematic to isotropic (NI) transition for single component nematics [19].But there is no systematic study to investigate how the composition changes the entropic value of NLC mixtures at NI transition.Experimental measurements of the entropic behavior are limited to specific nematic mixtures such as the mixtures of EPPV, PBPA and PBPA, EBBA [12,[20][21][22].However, a systematic study which includes a broad range of materials and their intrinsic properties such as polydispersity and molecular interaction between the components has not been performed.Moreover, these studies are not directed to CMs which are a mixture of discotic nematogens formed by the species at a certain range of molecular weight.In this paper, we study the general thermodynamic behavior, the entropic behavior, and the characteristic features of specific heat of a mixture of two thermotropic uniaxial discotic nematogens which represents a CM that only differs in molecular weight, thus precluding phase separation.Based on our previous study [1], depending on the intrinsic properties of the system, two types of mixtures arise: ideal and non-ideal (Figure 3).The ultimate objective of the current study is to use thermal properties and to determine the type of the mixture as well as its concentration that corresponds to the minimum transition temperature, critical concentration [1].
Maier-Saupe (MS) theory is widely used to describe the thermodynamics of nematic liquid crystals [25][26].This mean field theory gives the temperature-dependence of the molecular orientation as well as macroscopic entropic aspects in mesogenic materials.It predicts the values of the experimentally measured scalar order parameters very well [25], and hence has been applied to different nematic liquid crystalline systems and can be adjusted to their mixtures [26][27][28][29].In this paper, we use the Maier-Saupe (MS) theory generalized to binary uniaxial discotic nematogens.For binary mixtures of uniaxial mesogens the mixture quadrupolar order parameter is related to that of the individual components [1]: where m i is the mole fraction of i th component.For uniaxial phase Q is given in terms of a temperature-dependent scalar order parameter s(T) and the average molecular orientation or director n: ( ) , where I is the unit tensor.For binary discotic nematogens at equilibrium we find collinear directors ( ) = n n and the mixture uniaxial scalar order parameter then is: Figure 3. Classification of the mixtures into two types; type A with the non-ideal behavior and type B with ideal behavior, based on their intrinsic properties: molecular weight asymmetry ΔMw and molecular interaction parameter β.For weakly interacting mixtures and /or small molecular weight asymmetry the NI transition line exhibits a minimum by increasing the concentration (region A); however, for sufficiently strongly interacting mixtures and /or sufficiently large molecular weight asymmetry the NI transition line is monotonic (region B) [1].
Nevertheless since the ordering states of each species are coupled and are also dependent on temperature and dilution, we expect, based on the above observations, that a binary mixture displays three states: (i) isotropic (I) : s 1 =0, s 2 =0.(ii) nematic (N 12 ) with s 1 ≥s 2 .(iii) nematic (N 21 ) with s 1 ≤s 2 .Figure 4 shows the typical phase diagram of these mixtures.Two kinds of transition take place within mixtures: (i) Nematic to Isotropic (NI) transition due to the thermal effect, and (ii) N 12 to N 21 transition within the nematic phase due to the concentration effect.For the later case s 1 = s 2 where N 21 converts to N 12 ; at this point nematic mixture behaves like a single component system, and the concentration corresponding to this transition is the critical concentration m 1c at which the NI transition temperature is a minimum.The results of our previous study [1] show that depending on the molecular weight difference and the molecular interaction, two types of uniaxial nematic mixtures arise (Figure 3): (i) non-ideal mixtures; both kinds of transitions are observed within non-ideal mixtures, and (ii) ideal mixtures; only the NI transition takes place within this type of mixtures; N 12 exists for all range of concentration; as a result, N 21 and the critical concentration do not appear in the ideal mixture.Ideal binary mixtures arise under sufficiently strong interaction and sufficiently high molecular weight differences, while non-ideal mixtures arise under weak interaction and small molecular weight differences (asymmetries).As each carbonaceous mesophase mixture leads to different carbon fiber structure, it is of a great importance to be able to determine the type of the mixture.).In N 12 (N 21 ) the higher (lower) molecular weight species has a higher molecular order parameter than the lower (higher) molecular weight component: s 1 >s 2 (s 2 >s 1 ).Depending on the type of the mixture N 21 can appear and the transition within the nematic phase between N 21 and N 12 takes place (for non-ideal mixtures) or it can disappear (for ideal mixtures).The concentration which corresponds to the transition within the nematic phase is the critical concentration [1].
In previous work [1] we used the values of the X-ray diffraction intensity and the average ordering to characterize the type of the mixture and its critical concentration.In this paper we use specific heat and entropy, to characterize the binary mixtures.Figure 5 summarizes our objective schematically.It shows the schematic of the entropic values versus ordering for different concentration.S r1(2) is the value of the entropy for a single component system.This value is universal at the NI transition.Our objective is to find out the effect of concentration, molecular weight asymmetry and the interaction between the components on the (i) entropic values and the entropic jump at transition and (ii) the rate of entropy change which is needed to calculate the values of the specific heat.
The specific objectives of this paper are: (i) To determine the entropic behavior of ideal and non-ideal mixtures.
(ii) To calculate heat capacities of the two types of mixtures.
(iii) To determine their entropic variation at transition.The organization of this paper is as follows: Section 2 presents the essential aspects of the MS binary mixture model.Orientational entropy of the binary mixture is also introduced in this section.The main parameters influencing the thermodynamics of the system are identified and the numerical solution scheme is defined.Section 3 presents the derivation of the orientational specific heat.This quantity is used to characterize the type of the mixture, ideal or non-ideal, as well as the critical concentration at which a non-ideal mixture behaves as a pure mesogen.Section 4 presents the numerical results and discussion: Section 4.1 presents entropic behavior of different types of mixtures for different concentrations.Section 4.2 discusses the specific heat of different ideal and non-ideal mixtures as a tool to characterize the type of mixture and its critical concentration.It also shows the effect of molecular weight asymmetry on specific heat values.Section 4.3 presents the effect of concentration on the entropy jump at NI transition for two types of mixtures: ideal and non-ideal.Finally, the conclusions are presented in section 5.

Maier-Saupe Binary Mixture Model
Details of the extension of Maier-Saupe model to a binary mixture of NLCs is given in our pervious study [1].Here we briefly present the main features of the model which are necessary to obtain the entropic behavior and equilibrium phase diagram.The Helmholtz free energy per unit mole of the homogeneous mixture A is: where N A , k B , T are the Avogadro's number, Boltzmann's constant and temperature respectively, and E mix , S mix , and Z are the internal energy , the orientational entropy, and the partition function of the mixture per molecule respectively.E mix is given by the summation of three contributions: where { } W ,W ,W are the composition and molecular weight-dependent parameters [1].S mix as the entropy of mixing per molecule and is given by: and Z as the mixture partition function is factorized as: where Z i , i Φ and m i are the partition function, the partial internal potentials, and the mole fraction (concentration) of the i th species respectively.The two equations of equilibrium are obtained by minimizing the free energy A [eqn. (3)] with respect to the species order parameters ( ) s ,s .Scaling the Helmholtz energy A with the bare interaction parameter 11 U , minimizing the resulting dimensionless free energy with respect to s 1 , s 2 we find: where ϕ i and℘ i are effective mole fraction and the energetic parameter of the component "i" respectively.α is the angle between the directors of two components.The asymptotic limits of eqn.( 6) (m 1 =0,1) correspond to the pure NLC.As discussed in our pervious work [1] for the present CM case, the mixture is uniaxial and hence the relative director aperture is α=0.The thermodynamics of the mixture is defined by the dimensionless temperature T r and two effective mole fractions ( ) The two material parameters of the model are: ; which are functions of the molecular weights Mw i of the two components: where the parameters (a, b), based on the data of [30] are taken to be: a=-150 and b=0.75.Since the molecular weight dependence of 1 ℘ is not suggested by actual data, we use: where β is a constant whose sign depends on the geometrical nature of the species.For similar components, such disk/disk or rods/rod mixtures, it is positive; however, for dissimilar ones it is negative [31].As CMs are mixtures of discotic molecules the interaction parameter is positive (β>0).We choose component "1" as a representative higher molecular weight component of a CM with Mw 1 =1,400 [32] and vary the lower molecular weight of the second component, Mw 2 , so that the molecular weight asymmetry ΔMw=Mw 1 -Mw 2 > 0 changes.
The present thermodynamic model is given by the two nonlinear integral equations [eqns.(7)]; the solution vector consists of the two scalar order parameters (s 1 ,s 2 ); the two material parameters are β and ΔMw; the thermodynamic phase diagram is obtained by sweeping over temperature T r and concentration m 1 .Equations (7a&b) are solved by the Newton-Raphson method, using an eighth order Simpson integration method.Stability, accuracy, and convergence were ensured using standard methods [33].In the discussion of results we use the following nomenclature for dimensionless (reduced) entropy:

Specific Heat
The difference between dimensionless heat capacity of the nematic phase and that of the isotropic phase, ΔC pr , can be obtained by the orientational entropy which results from the ordering of the mixture [34]: As the isotropic entropy is assumed to be zero, the value given by eqn.13 is the orientational part of the specific heat.Therefore, .In this case, ΔC pr does not behave monotonically and exhibits a local extremum.

Results and Discussion
The Maier-Saupe for pure mesogens, whose interaction parameter is W, predicts the following universal values: Hence, deviations form these values will provide information on the mixture, as shown below.Figure 6a shows the dimensionless entropy, S r , as a function of temperature, for an ideal mixture with ΔMw=800 and interaction β=1 (Figure 3); each curve shows the entropic trend for a specific concentration m 1 .The following features are observed in this figure: (i) For any given concentration, entropy increases monotonically with increasing temperature (Figure 5).(ii) For any given temperature, entropy increases with decreasing concentration: component "1" has higher Mw, so it is more ordered in the nematic range.Therefore, as the concentration m 1 increases the mixture ordering increases and entropy decreases.(iii) Entropy jumps at the nematic-isotropic first order phase transition that sets in at temperature T NI (the entropy of the isotropic phase is assumed to be zero; as a result, this entropy represents the orientational entropy); the entropy jump value which is defined as (entropy of the isotropic phase) -(entropy of the nematic phase) at NI transition temperature for the almost pure mixtures (m 1 =0.001,0.999)are essentially the universal value rNI S 0.417719 = , and mixing just decreases the magnitude of jump.(iv) Depending on the relative population of two components there are three distinct concentration regions:

Entropic Behavior
. The rate of entropy change (δS r /δT r ) which corresponds to the rate of ordering change in the mixture shows a rapid decrease in region (ii).Component "2" has lower Mw; as a result, it has a lower T NI (eqn.10).Therefore, while component "1" tends to remain nematic (N 12 state), component "2" tends to transform to the isotropic (I) state and lose the ordering, though due to the interaction between the components it retains a low nematic ordering.As a result, when approaching the NI transition temperature of "2" the rate of ordering change in this component decreases and because of eqn.
(2) the overall rate of change in the mixture decreases.The appearance of the inflection point in the graphs corresponds to this change in the rate of entropy trend.This phenomenon becomes enhanced when the relative concentrations of both components are significant (region (ii) m 1 [0.2-0.3]).Figure 6b shows the entropic behavior for ΔMw=400 and interaction β=0.1 which represents a nonideal mixture (Figure 3); the concentration m 1 is shown for each curve.The following features are observed in this Figure : (i) Like the ideal case increasing the temperature decreases the ordering and increases the entropy (Figure 3) of all the mixtures.(ii) Unlike the ideal case, the effect of concentration on the ordering and entropy is different at different temperatures, increasing concentration does not necessarily decreases the entropy.The following three distinctive dilution regions arise: . Region (i) is located in N 21 state (s 2 >s 1 ) of the phase diagram (lower left region in Figure 3); in this region the lower molecular weight species "2" is the majority component.Therefore, increasing m 1 , dilutes the mixture, decreases the ordering and increases the entropy (Figure 5).The second and the third regions, both are located in N 12 state (s 1 >s 2 ); therefore, the trend is opposite.As component "1" is dominant in this region, increasing its concentration makes the mixture more ordered and decreases the entropy.(iii) there is a difference between region (ii) and region (iii).Like the ideal case, the rate of entropy change (δS r /δT r ) in the mixture shows a rapid decrease in region (ii).When approaching the NI transition temperature of "2" the rate of ordering change in this component decreases; as a result, the overall rate of change in the mixture decreases [eqn.(2)].This phenomenon becomes enhanced in the region (ii), where m 1 =[0.4-0.6] and the concentration of both components is significant.

Figure 6.
Dimensionless entropy Sr as a function of dimensionless temperature Tr of ideal (with ∆Mw=800 and β=1) (a) and a non-ideal (with ∆Mw=400 and β=0.1)(b) mixtures for different concentrations, computed using eqns.(5,12).Ideal mixtures show monotonic behavior with respect to dilution, while non-ideal mixtures show non-monotonicity due to the transition between N 12 and N 21 .The entropy behavior reflects the scalar order parameters of the species and their relative concentration.

Specific Heat
Figure 7a shows the value of the specific heat versus temperature for different concentration for the ideal case with ΔMw=800 and β=1.The following significant features are observed: (i) At any given temperature, increasing the concentration decreases the ordering and increases the specific heat.(ii) By sweeping over temperature, specific heat behaves differently for different concentrations.
Depending on the concentration region ((i) m 1 [0-0.1],(ii) m 1 [0.2-0.3],(iii) m 1 [0.4-1] , (see Section 4.1) two different trends are observed in the specific heat values versus temperature: in region (i) and (iii), when the population of one of the components is insignificant, and (δS r /δT r ) changes monotonically, specific heat increases monotonically with temperature; however, in region (ii), where m 1 =[0.2-0.3],specific heat increases by increasing temperature and then it shows a decrease which is followed by another increase.As mentioned in section 3, the nonmonotonicity of the specific heat is due to the direct contributions of the T r and (δS r /δT r ) to its value.According to section 4.1 (δS r /δT r ) in region (ii) shows a rapid decrease in the vicinity of T NI of component "2"; as a result, ΔC pr decreases in this range.On the other hand, increasing temperature increases ΔC pr ; as a result, a minimum appears in the ΔC pr as a function of temperature.Figure 7b shows the value of the specific heat versus temperature for different concentrations for the non-ideal case with ΔMw=400 and β=0.1 (Figure 3).The significant features observed are: (i) A transition is observed in ΔC pr trend vs. concentration.This transition takes place at the critical concentration where N 21 → N 12 .Region (i) is located in N 21 where increasing m 1 decreases the ordering and increases ΔC pr ; however, regions (ii) and (iii) are both located in N 12 where increasing m 1 increases the ordering and decreases ΔC pr .(ii) Depending on the concentration region ((i) m 1 [0-0.3],(ii) m 1 [0.4-0.6],(iii) m 1 [0.7-1], (see Section 4.1) two different trends are observed in the specific heat values versus temperature: in region (i) and (iii), when the effect of one of the components is weak, (δS r /δT r ) changes monotonically and specific heat increases monotonically with temperature; however, in region (ii), where m 1 =[0.4-0.6] and hence essentially no majority component, (δS r /δT r ) changes non-monotonically (section 4.2).As (δS r /δT r ) has a contribution to the specific heat values, ΔC pr exhibits a minimum vs. temperature in this region (see section 3).To better understand the effect of species relative population on ΔC pr values we investigate ΔC pr for three mixtures with different molecular weight differences (asymmetries) ΔMw=400 (a), 600(b), and 800(c) for m 1 =0.2 and β=0.1.The results are shown in Figure 8.Based on Figure 3 all three mixtures belong to the non-ideal type.The values of the critical concentration for these mixtures are: 0.31, 0.205 and 0.11 respectively [1].As a result, for mixtures (a) and (b) component "2" is the majority species, so they exhibit a monotonic ΔC pr ; however, in mixture (c) there is no majority component and hence ΔC pr exhibits an extremum.
The interaction parameter, as another intrinsic property changes the type of the mixture and the critical concentration (Figure 3); as a result, it also affects the ΔC pr values.For brevity the related graphs are not presented here.In partial summary, the higher Mw species always controls the overall ordering of the ideal mixture; as a result, the specific heat of ideal mixtures shows a decrease by increasing concentration at any given temperature; On the other hand, in non-ideal mixtures of weakly interacting species a transition in the trend of ΔC pr versus concentration is observed.This transition takes place within the nematic phase from N 12 to the N 12 state.For both ideal and non-ideal type, whenever there is no majority component in the mixture, the trend of specific heat versus temperature is non-monotonic because each component tends to show different ordering at different temperature range.Figure 8. Specific heat ΔC pr as a function of dimensionless temperature for m 1 =0.2 and β=0.1 for three different molecular weight asymmetries: ∆Mw=400 (a), 600 (b), and 800(c).All three cases are non-ideal mixtures; in case (c) the relative concentration of the components is significant, so it shows a non-monotonic trend with respect to temperature.

Entropy Jump at Transitions
Figure 9a shows the entropy jump at the NI transition for an ideal case, ΔMw=800 and β=1.For the pure systems (two asymptotic limits m 1 =0, 1) the value of the entropy jump is the universal value 0.417 [19].On the other hand, diluting the system disturbs the orientation and increases the entropy.As a result, there is a concentration which corresponds to the minimum entropy jump at NI transition.This concentration depends on the intrinsic properties of the mixture and is about 0.2 in the present case.Increasing dissimilarity of the components (larger ΔMw) enhances the dominancy of component "1" and shifts the minimum into lower concentrations.
Figure 9b shows the values of the entropy jump at transition for a non-ideal case, ΔMw=400 and β=0.1.In this case a local maximum is observed in the graph.The concentration corresponding to the local maximum is the critical concentration where the mixture behaves like a single component NLC (pseudo-pure mixture).Therefore, the value of the entropy for this concentration is the universal value 0.417719 and identical to that of the pure systems.Like the ideal mixture, diluting the pure (or pseudopure) mixture increases the entropy.So, the first local minimum is observed between m 1 =0 (pure system) and m 1c (pseudo-pure mixture) which is located in the N 21 region and the second one is observed between m 1 = m 1c (pseudo-pure mixture) and m 1 =1 (pure system) which is located in the N 12 region.

Conclusions
This paper uses the Maier-Saupe model for binary mixtures of discotic nematogenes to calculate the orientational entropy and orientational specific heat.These quantities are used to characterize the type of binary mixtures, ideal or non-ideal, which arise due to different intrinsic parameters, Mw asymmetry ΔMw and the interaction parameter β: (i) for ideal mixtures entropy and specific heat changes monotonically as a function of concentration at any given temperature.The entropic jump at NI transition exhibits a minimum for this kind of mixture.(ii) For non-ideal mixtures entropy and specific heat change non-monotonically as a function of concentration.A transition within the nematic phase from N 21 to N 12 is observed in this type of mixture.This transition takes place at the critical concentration and can be detected by the entropic and specific heat trend: the trend of these quantities in N 21 is opposite to their trend in N 12 .The entropic jump at the NI transition exhibits two local minima; one located in the N 21 region and the other in N 12 region.It also has a local maximum between these two regions at critical concentration, where the mixture behaves like a single component system and shows the universal entropy jump at NI transition.For both type of mixtures whenever there is a distinct majority component, ΔC pr and entropy behaves monotonically versus temperature; however, when there is no majority component ΔC pr and entropic trends are non-monotonic with temperature.
In the summary, heat capacity and transition entropies are shown to be useful tools to characterize the type of the discotic nematic liquid crystal and carbonaceous mesophase mixtures, to determine nematic ordering and to asses the degree of molecular interaction and molecular weights using standard calorimetric methods.

Figure 1 .
Figure 1.The comparison of three states of the material: a perfect crystalline solid with positional and orientational order (a), a liquid crystal with positional disorder and partial orientational order (b) and an isotropic liquid with positional and orientational disorder (c).

Figure 4 .
Figure 4. Schematic of the phase diagram of the binary mixture which includes three states: (i) isotropic (I), (ii) nematic (N 12) and, (iii) nematic (N 21 ).In N 12 (N 21 ) the higher (lower) molecular weight species has a higher molecular order parameter than the lower (higher) molecular weight component: s 1 >s 2 (s 2 >s 1 ).Depending on the type of the mixture N 21 can appear and the transition within the nematic phase between N 21 and N 12 takes place (for non-ideal mixtures) or it can disappear (for ideal mixtures).The concentration which corresponds to the transition within the nematic phase is the critical concentration[1].
pr C Δ shows the difference between C pr of the nematic phase and C pr of the isotropic phase.According to this equation two terms have direct contributions in the value of the specific heat: temperature, r T , and the rate of entropy, as a function of temperature changes.Once it increases monotonically, ΔC pr increases; however, once it decreases ΔC pr can either increases or decreases, depending the magnitude of each (i) NI transition temperature: B NI k T / W 0.22019 = (ii) scalar order parameter in the N phase at transition: s NI =0.4289 (iii) Latent heat at the NI transition:

Figure 5
Figure 5 shows the dimensionless entropy S r as a function of ordering s, for a pure mesogen and an schematic for a generic mixture.As expected, entropy decreases with increase in ordering.At the NI transition the pure component displays the universal values (s NI =0.4289, rNI S 0.417719 = − ), but for the generic mixture, the Figure shows a characteristic deviation for the universal transition values, which are discussed and quantified in detail in the rest of the paper.

Figure 5 .
Figure 5. Schematic of the effect of concentration m 1 , molecular weight asymmetry ΔMw and the interaction parameter β on the type, entropy and ordering of the mixture.It shows the general trend of the dimensionless entropy S r as a function of scalar order parameter s mix .As ordering increases the entropy decreases.It also shows that diluting the pure system changes the ordering and the entropy at NI transition.

Figure 7 .
Figure 7. Specific heat ΔC pr as a function of dimensionless temperature for ideal (with ∆Mw=800 and β=1) (a) and non-ideal (with ∆Mw=400 and β=0.1)(b) mixtures for different concentrations.Ideal mixtures show monotonic behavior of specific heat while non-ideal mixtures display non-monotonicity with respect to dilution.Both cases show non-monotonicity with respect to temperature whenever the relative concnetration of the components is significant.

Figure 9 .
Figure 9. Entropic variations at NI transition as a function of concentration for ideal (with ∆Mw=800 and β=1) (a) and non-ideal (with ∆Mw=400 and β=0.1)(b) mixtures for different concentrations.Both pure and pseudo pure mixtures have the universal entropy at NI transition.A local minimum is observed for the ideal case, where there is only N 12 state; however, two local minima (one in N 21 and the other in N 12 ) and a maximum (in N 12 →N 21 transition, which corresponds to the critical concentration) are observed for the non-ideal case.