# Length Bidisperse Carbon Nanotubes Dispersions in Thermotropic Liquid Crystals

Authors: V. Popa-Nita; S. Buček
Journal: Physics Research International
 (2012)

---

## Abstract

We study nematic liquid crystal driven alignment of carbon nanotubes dispersed in them. We extend the mesoscopic model presented in (P. Van der Schoot et al. 2008, V. Popa-Nita, and S. Kralj 2010) including the effect of length bidispersity of carbon nanotubes. The free energy of the mixture is written as the sum of the Doi free energy for lyotropic nematic ordering of the two carbon nanotubes types, the Landau-de Gennes free energy for the thermotropic ordering of liquid crystal, and the coupling term between liquid crystal molecules and carbon nanotubes. The phase ordering of the mixtures is analyzed as a function of volume fraction, the strength of coupling, and the temperature.


---

## Body

## 1. Introduction

Carbon nanotubes (CNTs) either consist of a single sheet of carbon rolled up into a cylinder with a diameter of about one nanometer (single-walled nanotubes—SWNTs), or are built up of multiple carbon sheets producing rods with diameters ranging from a few to tens of even hundreds of nanometer (multiwalled nanotubes—MWNTs). These nanoparticles have a large aspect ratio which varies from hundred to many thousands. CNTs exhibit most of their remarkable properties in a single direction, that is, along the tube axis [1].Most of CNTs extraordinary properties of potential use in various applications could be realized in relatively well-aligned samples. Recently, it has been shown that liquid crystal alignment could trigger spontaneous ordering of CNTs with remarkably high degree of ordering [2–5]. CNTs orient parallel to average direction of liquid crystal (LC) alignment with an orientational order parameter between 0.6 to 0.9 [6–9]. Both, thermotropic [6, 7, 9] and lyotropic nematic LC phases [8] have been successfully applied as aligning solvents.Theoretically, equilibrium orientation of a single elongated particle immersed in a nematic LC phase in rather well explored [10–17]. Continuum theory predicts alignment of the particle’s longer axis along the nematic director n→ for different types of boundary conditions in the strong anchoring limit [10–12]. In the weak anchoring limit, a rod-like particle may orient either along or perpendicular to the nematic director depending on the boundary conditions [11].The theoretical study of the collective behaviour of anisotropic nanoparticles dispersed in isotropic solvents or in liquid crystals is based on the observation that they can be considered essentially as rigid-rod polymers with a large aspect ratio [18]. The Onsager’s theory for the electrostatic repulsion of long rigid rods has been used to investigate the phase behaviour of SWNTs dispersed into organic and aqueous solvents [19]. In an adequate solubility, when the van der Waals attractive interaction between CNTs is overcome by strong repulsive interrods potential, the ordered phase of CNTs can form at room temperature. On the contrary, the van der Waals attractive interactions between the rods are strong, and as a result, only extremely dilute solutions of SWNTs are thermodynamically stable and no liquid crystal phases form at room temperature. The liquid crystalline of CNTs with and without van der Waals interactions has been analyzed by using the density functional theory [20]. In the presence of van der Waals interaction, the nematic as well as the columnar phases occur in the temperature-packing fraction phase diagram in a wide range of very high temperatures. In the absence of van der Waals interaction, with an increase of packing fraction, the system undergoes an isotropic-nematic phase transition via a biphasic region. The isotropic-nematic packing fraction decreases with the increase of the aspect ratio of CNTs. To describe the dispersion of SWNTs in superacids, the Onsager theory for rigid rods was extended to include the length polydispersity and solvent-mediated attraction and repulsion [21]. The main conclusion of these theoretical models is that to obtain liquid crystal phases of CNTs at room temperature the strong van der Waals interaction between them must be screened out. This requires an adequate solubility of the CNTs down to the level of individual tube. Note that nanoparticles could strongly interact via LC ordering [22–25] and also topological defects [24, 26] as recent studies reveal. Furthermore, the CNTs could induce into system some degree of quenched disorder [27].In [28, 29], a phenomenological theory for predicting the alignment of length monodisperse CNTs dispersions in thermotropic nematic liquid crystals was developed. The Landau-de Gennes free energy for thermotropic ordering of the liquid crystal solvent was combined with the Doi free energy for the lyotropic nematic ordering of CNTs caused by excluded-volume interactions between them. In the first paper [28], the interaction between CNTs and liquid crystal molecules is thought to be sufficiently weak to not cause any director field deformations in the nematic host fluid. In the second paper [29], the strong anchoring limit case was analyzed for which the nematic ordering around nanotubes is apparently distorted and consequently relatively strong long range and anisotropic interactions can emerge within the system. In both cases, the phase ordering of the binary mixture as a function of volume fraction of CNTs, the strength of the coupling and the temperature were analyzed.Usually, after the acid treatment and ultrasonification (to enhance the dispersibility and stability of the CNTs suspensions), depending on the temperature water bath and time of ultrasonification, the obtained CNTs have different lengths and diameters [7, 30]. The influence of length bidispersity on the phase diagram and alignment of CNTs is the subjects of the present paper. As the first step we extended our mesoscopic model [28, 29] and consider length bidispersity of CNTs.The remainder of this paper is organised as follows. In Section2, we introduce the free energy of mixture of CNTs and LC. In Section 3, we present our results and in Section 4 we draw conclusions.


## 2. Model

The mixture is characterized by the volume fractions of the three components:(1)Φi=Nivi∑j=13Njvjwith∑i=13Φi=1,
where Ni is the number of molecules of component i (i=1 defines the liquid crystal, i=2 the CNTs of length L02 and diameter D, and i=3 the length L03 and diameter D) and vi is the volume of a particle of component i. For latter convenience, we introduce scaled CNT lengths Lj=L0j/D (j=2,3).The degree of alignment of every component of the mixture is characterized by the scalar order parameterSi [31]. The corresponding isotropic liquid of the component i is characterized by Si=0 while a perfectly oriented nematic phase would correspond to Si=1.


### 2.1. Free Energy

The free energy per unit volume of mixture is expressed as(2)f=fCNT+fLC+fc,
where fCNT describes contribution of the two CNT components dispersed in the LC fluid, fLC represents the free energy density of nematic liquid crystal order, while fc takes into account the coupling between LC molecules and CNTs, respectively.The free energy density of CNTs is given by [28, 29]
(3)fCNT=kBT∑i=23{Φiviln⁡Φi+LiΦi26vi[(3LiΦi-1)Si2-23Si3+Si4Φiviln⁡Φi+]+LiΦi26vi[(3LiΦi-1)Si2-23Si3+Si4Φiviln⁡Φi+]}-γ23Φ2Φ3S2S3,
where kB is the Boltzmann constant and T is the temperature. The first term in the sum represents the entropic isotropic contribution due to mixing of the two CNT components neglecting their orientational degree of ordering [32]. In the limit of small volume fractions of the two CNTs components, the isotropic (Flory-Huggins) interaction terms of the form χΦiΦj can be neglected. The second term in the sum (usually called Doi free energy) that accounts for the gain in free volume and the loss of orientational entropy upon ordering of the rods for every CNT component is obtained starting from the Onsager theory [33] and using the Smoluchovsky equation. The details are presented in [34–36]. We mention that the Doi model neglects the van der Waals attractions between CNTs which are responsible for their tendency to form bundles. Minimizing the Doi free energy, a first-order isotropic-nematic phase transition is obtained which takes place at Φi(NI)=2.7/Li and Si(NI)=1/3. The last term in (3) represents the interaction energy between the different CNTs types, where the interaction parameter γ23 is given by
(4)γ23~8kBTπD3.
This result is obtained starting with the Onsager theory for the binary mixture and using the Smolukovsky equation [37]. The second term in (2) is the Landau-de Gennes free energy density [31] which describes the weakly first-order nematic-isotropic phase transition of thermotropic LC:
(5)fLC=(1-Φ2-Φ3)[a(T-T*)S12-BS13+CS14],
where (1-Φ2-Φ3) accounts for the part of the volume not taken up by LC. The quantities T*, a, B, and C are material-dependent constants, where T* denotes the spinodal temperatures limit of the isotropic phase of the pure LC. This free energy describes a weakly first-order nematic-isotropic phase transition. At T=TNI=T*+B2/(4aC), the two phases, nematic (S1(NI)=B/(2C)) and isotropic (S1=0) coexist in equilibrium.The third term in (2) represents the coupling between the liquid crystal molecules and the two CNT species. The resulting coupling term structure in both anchoring limits (weak and strong) we estimated in [28, 29]. The two limits are defined by the ratio D/de [10–12], where de=K/W stands for the surface extrapolation, W is the anchoring energy of an LC-nanotube interface, K is the average Frank nematic elastic constant [31]. For typical thermotropic LCs and nanotubes, it holds ≈10−11 N, W≈10-6N/m [38–41], and the penetration length should be vastly in excess of μm. The radius of CNTs is at the nm scale, therefore, D≪de. Consequently, only the weak-anchoring limit needs to be considered, as already concluded by Lynch and Patrick [6]. The corresponding free energy density of coupling is approximately given by [28]
(6)fc=-γ12(1-Φ2-Φ3)Φ2S1S2-γ13(1-Φ2-Φ3)Φ3S1S3.
The coupling parameters γ12 and γ13 depend only on the anchoring energy of CNTs at the LC molecules surface and the diameter of CNTs [28] (i.e., it does not depend on the length of CNT). Within this approximation, it holds
(7)γ12~γ13=γ1=4W3D.
Therefore,
(8)fc=-γ1(1-Φ2-Φ3)S1(Φ2S2+Φ3S3).


### 2.2. Material and Geometrical Parameters

In this paper, we will consider the following typical geometrical and material values. For CNTs, we setD=2 nm, L02=400 nm, and L03=800 nm. Considering T=300 K, the coupling constant γ23 has the value 106 J/m3, which we use in the following. For representative LC material, we chose pentylcyanobiphenyl (5CB), for which T*=307.55 K, a≈5.2·104 J/m3 K, B≈5.3·105 J/m3, C≈9.7·105 J/m3, and K~10-11 N [38–41]. This chose yields S1(NI)=0.275, TNI=308.95 K. For the LC-CNT interaction, we useγ1=4W/3D≈103 N/m2, where we have assumed W~10-6 N/m and D=2 nm.


### 2.3. Basic Equation

The equilibrium values of the order parameters are obtained by minimization of the free energy density ((2), (3), (5), and (6)) with respect to S1, S2, and S3, respectively. From the corresponding equation, we find the equilibrium values of the order parameters in the nematic phase (Sin) and paranematic phase (Sip), respectively. Once the minimization procedure has been solved, the volume fractions of the coexisting (paranematic and nematic) phases are found by solving the equilibrium conditions:
(9)μ2(S1n,S2n,S3n,Φ2n,Φ3n)=μ2(S1p,S2p,S3p,Φ2p,Φ3p),μ3(S1n,S2n,S3n,Φ2n,Φ3n)=μ3(S1p,S2p,S3p,Φ2p,Φ3p),g(S1n,S2n,S3n,Φ2n,Φ3n)=g(S1p,S2p,S3p,Φ2p,Φ3p),
where the chemical potential of the two species of CNTs μ2, μ3, and the grand potential g are defined by
(10)μ2=∂f∂Φ2;μ3=∂f∂Φ3;g=f-μ2Φ2-μ3Φ3.
There are three equilibrium equations with four variables: Φ2n, Φ3n, Φ2p, and Φ3p. We take Φ3p as freely variable and calculate the other three Φ2n, Φ3n, Φ2p from the coexisting equations (9).


## 2.1. Free Energy

The free energy per unit volume of mixture is expressed as(2)f=fCNT+fLC+fc,
where fCNT describes contribution of the two CNT components dispersed in the LC fluid, fLC represents the free energy density of nematic liquid crystal order, while fc takes into account the coupling between LC molecules and CNTs, respectively.The free energy density of CNTs is given by [28, 29]
(3)fCNT=kBT∑i=23{Φiviln⁡Φi+LiΦi26vi[(3LiΦi-1)Si2-23Si3+Si4Φiviln⁡Φi+]+LiΦi26vi[(3LiΦi-1)Si2-23Si3+Si4Φiviln⁡Φi+]}-γ23Φ2Φ3S2S3,
where kB is the Boltzmann constant and T is the temperature. The first term in the sum represents the entropic isotropic contribution due to mixing of the two CNT components neglecting their orientational degree of ordering [32]. In the limit of small volume fractions of the two CNTs components, the isotropic (Flory-Huggins) interaction terms of the form χΦiΦj can be neglected. The second term in the sum (usually called Doi free energy) that accounts for the gain in free volume and the loss of orientational entropy upon ordering of the rods for every CNT component is obtained starting from the Onsager theory [33] and using the Smoluchovsky equation. The details are presented in [34–36]. We mention that the Doi model neglects the van der Waals attractions between CNTs which are responsible for their tendency to form bundles. Minimizing the Doi free energy, a first-order isotropic-nematic phase transition is obtained which takes place at Φi(NI)=2.7/Li and Si(NI)=1/3. The last term in (3) represents the interaction energy between the different CNTs types, where the interaction parameter γ23 is given by
(4)γ23~8kBTπD3.
This result is obtained starting with the Onsager theory for the binary mixture and using the Smolukovsky equation [37]. The second term in (2) is the Landau-de Gennes free energy density [31] which describes the weakly first-order nematic-isotropic phase transition of thermotropic LC:
(5)fLC=(1-Φ2-Φ3)[a(T-T*)S12-BS13+CS14],
where (1-Φ2-Φ3) accounts for the part of the volume not taken up by LC. The quantities T*, a, B, and C are material-dependent constants, where T* denotes the spinodal temperatures limit of the isotropic phase of the pure LC. This free energy describes a weakly first-order nematic-isotropic phase transition. At T=TNI=T*+B2/(4aC), the two phases, nematic (S1(NI)=B/(2C)) and isotropic (S1=0) coexist in equilibrium.The third term in (2) represents the coupling between the liquid crystal molecules and the two CNT species. The resulting coupling term structure in both anchoring limits (weak and strong) we estimated in [28, 29]. The two limits are defined by the ratio D/de [10–12], where de=K/W stands for the surface extrapolation, W is the anchoring energy of an LC-nanotube interface, K is the average Frank nematic elastic constant [31]. For typical thermotropic LCs and nanotubes, it holds ≈10−11 N, W≈10-6N/m [38–41], and the penetration length should be vastly in excess of μm. The radius of CNTs is at the nm scale, therefore, D≪de. Consequently, only the weak-anchoring limit needs to be considered, as already concluded by Lynch and Patrick [6]. The corresponding free energy density of coupling is approximately given by [28]
(6)fc=-γ12(1-Φ2-Φ3)Φ2S1S2-γ13(1-Φ2-Φ3)Φ3S1S3.
The coupling parameters γ12 and γ13 depend only on the anchoring energy of CNTs at the LC molecules surface and the diameter of CNTs [28] (i.e., it does not depend on the length of CNT). Within this approximation, it holds
(7)γ12~γ13=γ1=4W3D.
Therefore,
(8)fc=-γ1(1-Φ2-Φ3)S1(Φ2S2+Φ3S3).


## 2.2. Material and Geometrical Parameters

In this paper, we will consider the following typical geometrical and material values. For CNTs, we setD=2 nm, L02=400 nm, and L03=800 nm. Considering T=300 K, the coupling constant γ23 has the value 106 J/m3, which we use in the following. For representative LC material, we chose pentylcyanobiphenyl (5CB), for which T*=307.55 K, a≈5.2·104 J/m3 K, B≈5.3·105 J/m3, C≈9.7·105 J/m3, and K~10-11 N [38–41]. This chose yields S1(NI)=0.275, TNI=308.95 K. For the LC-CNT interaction, we useγ1=4W/3D≈103 N/m2, where we have assumed W~10-6 N/m and D=2 nm.


## 2.3. Basic Equation

The equilibrium values of the order parameters are obtained by minimization of the free energy density ((2), (3), (5), and (6)) with respect to S1, S2, and S3, respectively. From the corresponding equation, we find the equilibrium values of the order parameters in the nematic phase (Sin) and paranematic phase (Sip), respectively. Once the minimization procedure has been solved, the volume fractions of the coexisting (paranematic and nematic) phases are found by solving the equilibrium conditions:
(9)μ2(S1n,S2n,S3n,Φ2n,Φ3n)=μ2(S1p,S2p,S3p,Φ2p,Φ3p),μ3(S1n,S2n,S3n,Φ2n,Φ3n)=μ3(S1p,S2p,S3p,Φ2p,Φ3p),g(S1n,S2n,S3n,Φ2n,Φ3n)=g(S1p,S2p,S3p,Φ2p,Φ3p),
where the chemical potential of the two species of CNTs μ2, μ3, and the grand potential g are defined by
(10)μ2=∂f∂Φ2;μ3=∂f∂Φ3;g=f-μ2Φ2-μ3Φ3.
There are three equilibrium equations with four variables: Φ2n, Φ3n, Φ2p, and Φ3p. We take Φ3p as freely variable and calculate the other three Φ2n, Φ3n, Φ2p from the coexisting equations (9).


## 3. Results

Our aim is to analyze phase behavior of immersed CNTs in nematic LC phase as a function ofΦ2, Φ3, and γ1. The form of the coupling between CNTs and the LC molecules (8) induces two important effects: (i) the isotropic phase of CNTs transforms to a paranematic phase exhibiting relatively low degree of orientational ordering and (ii) the appearance of a critical line γ1(c)(T). The reason behind this is that the coupling term fc is linearly proportional with Si.In the following, we first determine qualitatively different regimes in the(T,γ1) phase diagram on varying concentrations of CNTs. The critical line γ1(c)(T) separates two qualitatively different regimes. CNTs exhibit first order transition for γ1<γ1(c). On the other hand, CNTs display gradual variations in Si on varying Φi if γ1>γ1(c). Afterwards, we analyze phase behavior in the (Φ2,Φ3) plane for γ1<γ1(c) and γ1>γ1(c). We, henceforth, refer to the regimes γ1<γ1(c) and γ1>γ1(c) as subcritical and supercritical regime, respectively.


### 3.1. Phase Behaviour in the(T,γ1) Plane

The critical lineγ1(c)(T) is obtained by solving the equations ∂2f/∂Si2=∂3f/∂Si3=0, (i=2,3). They yield at a critical point universal values for the order parameter Si(c)=1/6 and volume fraction Φi(c)=18/7Li. The γ1(c)(T) dependence for the two species of CNTs is presented in Figure 1.The critical value of the coupling parameterγ1(c) as function of temperature calculated for both types of CNTs in the nematic LC phase.
(a)
(b)On decreasing the temperature the external field felt by the CNTs is increasing due to increasing value ofS1 and the nematic-isotropic phase transition of CNTs becomes gradual for increasingly lower values of interaction parameter γ1(c). Furthermore, a critical value γ1 decreases with increasing the length of CNTs. Therefore, the continuity of paranematic-nematic phase transition is favoured by longer CNTs.


### 3.2. Phase Behaviour in the(Φ2,Φ3) Plane

In what follows, we present the results in the two different cases: (i)γ1<γ1(c), when the two types of CNTs undergo first order parametric-nematic phase transitions and (ii) γ1>γ1(c), where the respective transitions are replaced by gradual evolution of ordering.


### 3.3. Subcritical Regime

The phase diagram of the bidisperse CNTs suspension in the nematic LC forT-TNI=-2.5 K for γ1=20N/m2<γ1(c) in terms of volume fractions is shown in Figure 2.Figure 2
The (Φ3,Φ2) phase diagram for a length bidisperse CNTs dispersed in nematic phase of the LC for T-TNI=-2.8K and γ1=20N/m2<γ1(c). Thick lines indicates phase boundary, while the representative thin lines connects the coexisting pairs (Φ2p,Φ2n) and (Φ3p,Φ3n). For the explication of the meaning of regions I, II, and III see the main text. The letters N and P denote nematic and paranematic ordering of CNTs, respectively.The differences of volume fractions of both CNTs types in the nematic LC as functions of volume fraction of longer CNTs in the paranematic phase are shown in Figure3. The temperature and the coupling parameter (T and γ1) are the same as in Figure 2.Figure 3
The difference of volume fraction of both CNTs types dispersed in nematic phase of LC forT-TNI=-2.8K and γ1=20N/m2<γ1(c). For the meaning of regions I, II, and III see the main text.The figures reveal two qualitatively different regimes, to which we refer asdecoupled and coupled regime, respectively. In the decoupled regime, systems exhibit monodisperse-type behavior. We label these regions by I and II in Figures 2 and 3. In the limit of very low volume fraction of the longer CNTs (Φ3→0), the system is monodisperse containing only one types of CNTs of length L02=400nm with volume fraction gap at the transition ΔΦ2=Φ2n-Φ2p=0.0143-0.0128=0.0015 (subregion I). In the limit of very low volume fraction of the shorter CNTs (Φ2→0), the system is monodisperse containing only one types of CNTs of length L03=800nm with a volume fraction gap at the transition ΔΦ3=Φ3n-Φ3p=0.007022-0.006333=0.000689 (subregion II). In these two subregions, due to very small values of Φ3 and Φ2, respectively, the coupling term between CNTs types (the γ23 term in (3)) is relatively small and the types are independent. On the contrary, in the coupled region (region III in Figures 2 and 3), the interaction term in (3) becomes important and the two CNTs types influence each other. In this region, the volume fraction of the longer CNTs increases in the nematic phase, while that of the shorter CNTs decreases.To see in more details the influence of bidispersity on the order parameters we have plotted in Figure4 the order parameters profiles for a monodisperse CNTs, while in Figure 5, we have shown the order parameters profiles for the bidisperse system for two fixed values of the volume fraction of the shorter CNTs. The values of temperatures and coupling parameter are the same as in previous figure: T-TNI=-2.8Kand γ1=20N/m2<γ1(c).Figure 4
Order parameters profiles for monodisperse CNTs forT-TNI=-2.8K and γ1=20N/m2<γ1(c).Order parameters profiles for two fixed values of the volume fractions of shorter CNTs forT-TNI=-2.8K, γ1=20N/m2<γ1(c).
(a)
(b)Figure4 shows that in the case of a monodisperse system below the critical point (γ1<γ1(c)) the paranematic-nematic phase transition is first order (discontinuity in the order parameter). The transition takes place at a value of volume fraction depending on the ratio L0i/D (Φ2(NI)=0.0133 for L2=200 and Φ3(NI)=0.0065 for L3=400). Longer CNTs exhibit nematic ordering at a lower volume fraction. In the nematic phase, the longer CNTs are more ordered than the shorter ones. Note that a value of the order parameter of the LC is almost constant on varying Φ if Figure 4 (i.e., value of S1 is dominated by temperature, not by the volume fraction of CNTs). Due to much larger free energy of the LC (there are many more molecules of LC per unit volume than there are CNTs) the CNTs types are enslaved by LC (the LC fluid dictates the behavior of CNTs).The order parameters profiles in Figure5(a) are represented for a value of Φ2 less than the corresponding Φ2(NI) for the monodisperse system with L2=200. However, the increasing volume fraction of longer CNTs induces a very weakly first-order paranematic-nematic phase transition of shorter CNTs, with a very small jump in the order parameter. This transition takes place at the same value Φ3 as the stronger paranematic-nematic first-order phase transition of longer CNTs (it is a hidden paranematic-nematic first-order phase transition which takes place inside the coexistence region in Figure 3). In the nematic phase, the longer CNTs are more ordered than the shorter ones. In Figure 5(b), we have shown the order parameter profiles for a value of Φ2 very close to ΦNI for themonodisperse system with L2=200. Now, the jump of order parameter of the shorter CNTs becomes larger that the corresponding quantity for longer CNTs, but in the nematic phase the order in longer CNTs remains larger than that of the shorter CNTs.


### 3.4. Supercritical Regime

In the case of a more realistic value of the coupling parameter between the LC molecules and CNTsγ1=103N/m2>γ1(c), we plot in Figure 6 the order parameter profiles for a monodisperse system. In Figure 7, we plot order parameter profiles for the bidisperse system for a fixed value of the volume fraction of shorter CNTs at the same temperature T-TNI=-2.8K.Figure 6
Order parameters profiles for monodisperse CNTs forT-TNI=-2.8K, γ1=103N/m2>γ1(c).Figure 7
Order parameters profiles for a fixed value of the volume fraction of shorter CNTs forT-TNI=-2.8K, γ1=103N/m2>γ1(c).Figure6 shows that in the case of monodisperse system above a critical point (γ1>γ1(c)), the alignment of the CNTs gradually increases on increasing their volume fraction. The longer CNTs are more ordered than the shorter ones and the CNTs types are enslaved by LC (the order parameter of LC depends only on temperature, not on the volume fractions of CNTs).The degree of orderingSi in bidispersed system is demonstrated in Figure 7. We display Si=Si(Φ3) dependencies at constant Φ2. The order parameter of the shorter CNTs increases but not significantly with the volume fraction of the longer CNTs (Figure 7). In relatively well-oriented state, the orientational ordering of longer CNTs is larger than that of a shorter rods.


## 3.1. Phase Behaviour in the(T,γ1) Plane

The critical lineγ1(c)(T) is obtained by solving the equations ∂2f/∂Si2=∂3f/∂Si3=0, (i=2,3). They yield at a critical point universal values for the order parameter Si(c)=1/6 and volume fraction Φi(c)=18/7Li. The γ1(c)(T) dependence for the two species of CNTs is presented in Figure 1.The critical value of the coupling parameterγ1(c) as function of temperature calculated for both types of CNTs in the nematic LC phase.
(a)
(b)On decreasing the temperature the external field felt by the CNTs is increasing due to increasing value ofS1 and the nematic-isotropic phase transition of CNTs becomes gradual for increasingly lower values of interaction parameter γ1(c). Furthermore, a critical value γ1 decreases with increasing the length of CNTs. Therefore, the continuity of paranematic-nematic phase transition is favoured by longer CNTs.


## 3.2. Phase Behaviour in the(Φ2,Φ3) Plane

In what follows, we present the results in the two different cases: (i)γ1<γ1(c), when the two types of CNTs undergo first order parametric-nematic phase transitions and (ii) γ1>γ1(c), where the respective transitions are replaced by gradual evolution of ordering.


## 3.3. Subcritical Regime

The phase diagram of the bidisperse CNTs suspension in the nematic LC forT-TNI=-2.5 K for γ1=20N/m2<γ1(c) in terms of volume fractions is shown in Figure 2.Figure 2
The (Φ3,Φ2) phase diagram for a length bidisperse CNTs dispersed in nematic phase of the LC for T-TNI=-2.8K and γ1=20N/m2<γ1(c). Thick lines indicates phase boundary, while the representative thin lines connects the coexisting pairs (Φ2p,Φ2n) and (Φ3p,Φ3n). For the explication of the meaning of regions I, II, and III see the main text. The letters N and P denote nematic and paranematic ordering of CNTs, respectively.The differences of volume fractions of both CNTs types in the nematic LC as functions of volume fraction of longer CNTs in the paranematic phase are shown in Figure3. The temperature and the coupling parameter (T and γ1) are the same as in Figure 2.Figure 3
The difference of volume fraction of both CNTs types dispersed in nematic phase of LC forT-TNI=-2.8K and γ1=20N/m2<γ1(c). For the meaning of regions I, II, and III see the main text.The figures reveal two qualitatively different regimes, to which we refer asdecoupled and coupled regime, respectively. In the decoupled regime, systems exhibit monodisperse-type behavior. We label these regions by I and II in Figures 2 and 3. In the limit of very low volume fraction of the longer CNTs (Φ3→0), the system is monodisperse containing only one types of CNTs of length L02=400nm with volume fraction gap at the transition ΔΦ2=Φ2n-Φ2p=0.0143-0.0128=0.0015 (subregion I). In the limit of very low volume fraction of the shorter CNTs (Φ2→0), the system is monodisperse containing only one types of CNTs of length L03=800nm with a volume fraction gap at the transition ΔΦ3=Φ3n-Φ3p=0.007022-0.006333=0.000689 (subregion II). In these two subregions, due to very small values of Φ3 and Φ2, respectively, the coupling term between CNTs types (the γ23 term in (3)) is relatively small and the types are independent. On the contrary, in the coupled region (region III in Figures 2 and 3), the interaction term in (3) becomes important and the two CNTs types influence each other. In this region, the volume fraction of the longer CNTs increases in the nematic phase, while that of the shorter CNTs decreases.To see in more details the influence of bidispersity on the order parameters we have plotted in Figure4 the order parameters profiles for a monodisperse CNTs, while in Figure 5, we have shown the order parameters profiles for the bidisperse system for two fixed values of the volume fraction of the shorter CNTs. The values of temperatures and coupling parameter are the same as in previous figure: T-TNI=-2.8Kand γ1=20N/m2<γ1(c).Figure 4
Order parameters profiles for monodisperse CNTs forT-TNI=-2.8K and γ1=20N/m2<γ1(c).Order parameters profiles for two fixed values of the volume fractions of shorter CNTs forT-TNI=-2.8K, γ1=20N/m2<γ1(c).
(a)
(b)Figure4 shows that in the case of a monodisperse system below the critical point (γ1<γ1(c)) the paranematic-nematic phase transition is first order (discontinuity in the order parameter). The transition takes place at a value of volume fraction depending on the ratio L0i/D (Φ2(NI)=0.0133 for L2=200 and Φ3(NI)=0.0065 for L3=400). Longer CNTs exhibit nematic ordering at a lower volume fraction. In the nematic phase, the longer CNTs are more ordered than the shorter ones. Note that a value of the order parameter of the LC is almost constant on varying Φ if Figure 4 (i.e., value of S1 is dominated by temperature, not by the volume fraction of CNTs). Due to much larger free energy of the LC (there are many more molecules of LC per unit volume than there are CNTs) the CNTs types are enslaved by LC (the LC fluid dictates the behavior of CNTs).The order parameters profiles in Figure5(a) are represented for a value of Φ2 less than the corresponding Φ2(NI) for the monodisperse system with L2=200. However, the increasing volume fraction of longer CNTs induces a very weakly first-order paranematic-nematic phase transition of shorter CNTs, with a very small jump in the order parameter. This transition takes place at the same value Φ3 as the stronger paranematic-nematic first-order phase transition of longer CNTs (it is a hidden paranematic-nematic first-order phase transition which takes place inside the coexistence region in Figure 3). In the nematic phase, the longer CNTs are more ordered than the shorter ones. In Figure 5(b), we have shown the order parameter profiles for a value of Φ2 very close to ΦNI for themonodisperse system with L2=200. Now, the jump of order parameter of the shorter CNTs becomes larger that the corresponding quantity for longer CNTs, but in the nematic phase the order in longer CNTs remains larger than that of the shorter CNTs.


## 3.4. Supercritical Regime

In the case of a more realistic value of the coupling parameter between the LC molecules and CNTsγ1=103N/m2>γ1(c), we plot in Figure 6 the order parameter profiles for a monodisperse system. In Figure 7, we plot order parameter profiles for the bidisperse system for a fixed value of the volume fraction of shorter CNTs at the same temperature T-TNI=-2.8K.Figure 6
Order parameters profiles for monodisperse CNTs forT-TNI=-2.8K, γ1=103N/m2>γ1(c).Figure 7
Order parameters profiles for a fixed value of the volume fraction of shorter CNTs forT-TNI=-2.8K, γ1=103N/m2>γ1(c).Figure6 shows that in the case of monodisperse system above a critical point (γ1>γ1(c)), the alignment of the CNTs gradually increases on increasing their volume fraction. The longer CNTs are more ordered than the shorter ones and the CNTs types are enslaved by LC (the order parameter of LC depends only on temperature, not on the volume fractions of CNTs).The degree of orderingSi in bidispersed system is demonstrated in Figure 7. We display Si=Si(Φ3) dependencies at constant Φ2. The order parameter of the shorter CNTs increases but not significantly with the volume fraction of the longer CNTs (Figure 7). In relatively well-oriented state, the orientational ordering of longer CNTs is larger than that of a shorter rods.


## 4. Conclusions

In the present paper, we extend the mesoscopic model [28, 29] to include length bidispersity of CNTs dispersed in LC in the weak anchoring limit of the coupling between LC molecules and CNTs (in this limit, the coupling is dominated by the anisotropy of the surface tension not by the deformation of the director field). The main conclusions of our study can be summarized as follows.(1)
Depending on the coupling between the LC molecules and CNTs (the value of the coupling constantγ1), two different regimes can be defined: (i) if γ1<γ1(c), the nematic-isotropic phase transition of CNTs dispersed in LC is first order and (ii) γ1>γ1(c) corresponds to a supercritical transformation (the evolution of the order parameter is gradual). γ1(c) is the critical value of the coupling parameter depending on the temperature (Figure 1). In both regimes, the isotropic phase of CNTs transforms into a phase with a small degree of ordering, a paranematic phase. Above the critical point, this degree of orientational order is strongly increasing.(2)
The CNTs are enslaved by the LC (the nematic LC order parameter depends only on the temperature not on the volume fractions of the CNTs).(3)
The longer CNTs are driven into the nematic phase (Φ3n-Φ3p>0 in Figures 2 and 3).(4)
In region III from Figures2 and 3, the longer CNTs induce a larger differences in the volume fraction Φ2n-Φ2p>0 for the shorter CNTs.(5)
In the nematic phase, the longer CNTs are more ordered than the shorter ones (S3n-S2n>0 in Figures 5 and 7).Finally, we point out the mesoscopic model for length bidisperse CNTs dispersed in LC presented here is only a first step in considering the important influence of polydispersity on the ordering of CNTs in nematic fluids. Together with considering the attraction interaction between CNTs, this subject will be studied in a future work.



---

*Source: 750890-2012-08-29.xml*