Based on the phenomenological model first presented by van der Schoot et al., which predicts the alignment of carbon nanotube (CNT) dispersions in thermotropic nematic liquid crystals, we present the extensive results concerning the phase diagram and the orientational properties of the mixture in this chapter.
- liquid crystal
- carbon nanotube
- phase transition
A method to obtain aligned CNTs (necessary for applications) is to disperse them into liquid crystals (LCs) [with their self-organization (long-range orientational order) and fluidity] [1, 2, 3, 4, 5, 6, 7, 8, 9]. The thermotropic [10, 11, 12, 13] and the lyotropic LCs [14, 15, 16, 17] have been used to align CNTs parallel as well as perpendicular to average direction of alignment of long axes of LC molecules called the director. As a consequence, the orientational order parameter of CNTs could have the values between and . The study of mixture composed by LCs and CNTs is also important due to the influence of CNTs on the physical properties of LCs (increased dielectric anisotropy, decreased threshold voltage, much accelerated electrooptical response) [18, 19, 20, 21, 22, 23, 24, 25].
The main hypothesis used in the theoretical study of the collective behavior of CNTs dispersed in the isotropic solvents [26, 27, 28] as well as in LCs [29, 30, 31, 32, 33, 34, 35, 36] is that they can be considered as rigid rod polymers .
Using the density functional theory, the isotropic-liquid crystal phase transition has been analyzed considering the van der Waals attractive interactions . Onsager theory of rigid rods  was used to study the phase behavior of CNTs dispersed into organic and aqueous solutions . Also, the Onsager model including length polydispersity and solvent-mediated interaction was considered to study the dispersions of CNTs in superacids . These theoretical studies lead to the conclusion that to obtain orientational order of CNTs at room temperature it is necessary that the van der Waals interactions must be screened out, i.e., the CNTs must be dispersed in a good solvent. In the case of a non-good solvent, no liquid crystalline phases of CNTs form at room temperature because only dilute solutions are thermodynamically stable.
Matsuyama [29, 30] used a mean field theory to analyze the phase behavior of a mixture composed of low-molecular-weight liquid crystal and a rigid rod-like polymer such as CNTs. The free energy is constructed on the basis of Onsager model for excluded volume interactions, the Maier-Saupe model for orientational-dependent attractive interactions, and the Flory-Huggins theory for binary mixtures.
In the previous papers [31, 32, 34, 35, 36], we have presented a phenomenological theory to describe the alignment of CNT dispersions in thermotropic nematic LCs. We combined the Landau-de Gennes [39, 40] free energy for thermotropic ordering of the LC solvent and the Doi free energy [41, 42, 43] for the lyotropic orientational order of CNTs. Because the CNT is much thinner than the elastic penetration length, the alignment of CNTs in the nematic solvent is caused by the coupling of the LC director field to the anisotropic interfacial tension of the CNTs. This is true only for very dilute solutions without large aggregates . Density functional calculations  show that CNT alignment mechanism in LC is associated with a strong interaction due to surface anchoring with a binding energy of eV for electron stacking between LC molecules and CNT.
In the present chapter, we extend this model to generalize the anisotropic interaction form between a CNT and the liquid crystal molecules including the possibility of perpendicular alignment.
The remainder of this chapter is organized as follows. In Section 2, we describe the model. In Section 3, we illustrate the results and finally give the main conclusions in Section 4.
The free energy density of the binary mixture composed of CNTs and thermotropic nematic LC contains three terms:
where represents the Doi free energy density of CNTs, represents the Landau-de Gennes free energy density of the nematic liquid crystal order, and takes into account the interaction of these two components.
2.1 Free energy of carbon nanotubes
At the mesoscopic level of description, the free energy of CNTs can be written in the following form :
where the first two terms represent the entropy of isotropic mixing of CNTs and LC components  and the third one is obtained from the Onsager theory  using the Smoluchowski Equation [41, 42]. In this form, the van der Waals attractions between CNTs are neglected. The long-ranged intermolecular attractions between CNTs were considered in . represents the volume of one LC molecule with length and diameter , and is the volume of one CNT with length and diameter . The values cited above are those used in this paper. In Eq. (2), is the Boltzmann constant, is the absolute temperature, represents the volume fraction of CNTs, defines the volume fraction of the LC, and is the aspect ratio of a CNT.
Eq. (2) defines a first-order isotropic-nematic phase transition of CNTs with the equilibrium values of the order parameters and and of the volume fractions, respectively, and . These equilibrium values of the parameters are obtained minimizing the free energy (2) with respect to and equating the chemical potentials of CNTs () and of LC ( in the two phases. The model predicts a relative variation of the volume fraction of CNTs (the Flory “chimney”) at the transition of 0.19%, while the prediction of Onsager and Flory [38, 46] for monodisperse lyotropic LC is about 1%, and the experimental results obtained by the polarized light microscopy for an aqueous dispersion of multiwall carbon nanotube indicate a biphasic region that range between 1 and 4% [47, 48]. The discrepancy is due to polydispersity, as well as the possibility of segregation which the model does not consider.
2.2 Free energy of thermotropic liquid crystal
To characterize the isotropic-nematic phase transition of thermotropic LC, we use the Landau-de Gennes  free energy:
where is the undercooling limit temperature of stability of the isotropic phase and the constant coefficients , , and depend on the material. Later, in our numerical calculations presented in Section 3, we consider the following values: K, , , and [40, 49, 50, 51], values specific for pentylcyanobiphenyl (5CB). For this nematogen, the first-order isotropic-nematic phase transition takes place at temperature with and .
2.3 The coupling free energy
The condition of weak anchoring limit of the interaction between the two components is , where the diameter of a CNT is and is the elastic penetration length ( is a typical elastic constant, and is a typical anchoring energy).
The coupling free energy has been explained in , and it has the following form:
where defines the coupling parameter. Thus, for the same anchoring energy, the coupling parameter in the case of thin rods is larger than that of the thick one. From the experimental results, we find that the value of the coupling parameter can be between . In the limit of large coupling , the orientational order of CNTs is perfect .
We note that positive values of characterize a parallel alignment of LC molecules to CNTs. In this chapter we extend our analysis considering also negative values of which characterize a perpendicular alignment. Depending on the orientational order parameters and , two nematic phases could be defined: (i) the phase with and (parallel alignment of CNT and LC molecules) and (ii) the phase with and (perpendicular alignment of CNT and LC molecules). These two phases are experimentally obtained depending on the surface treatment of the CNTs.
2.4 The total free energy density of the binary mixture
Finally, the total free energy density of the binary mixture (1) is the sum of the free energy densities of CNTs (2), thermotropic LC (3), and interaction between nematic CNTs and nematic LC (4):
We present the main results concerning the phase behavior of the CNT-LC binary mixture as a function of temperature, volume fraction of CNTs, and coupling strength .
3.1 Positive : critical point of CNT phase transition
The positive value of the interaction parameter defines the nematic phase with and (the rod and the LC molecules are parallel to each other). The isotropic phase of CNTs is transformed into a nematic phase with a small degree of order which is called paranematic phase (by analogy with the magnetism). The phase transition nematic-paranematic of the CNTs is first order for small values of the interaction parameter (so-called subcritical regime of CNTs). For large values of , the phase transition of CNTs becomes continuous. The critical point is given by the equations , giving and the effective volume fraction , while the critical value of the interaction parameter is given by the equation:
The dependence of as a function of temperature is shown in Figure 1.
In the nematic phase of the liquid crystal, the critical interaction parameter is relatively small (because has relatively large values between 0.28 and 0.8), while in the isotropic phase of the liquid crystal, the value of the order parameter is very small (), and the values of are relatively large. As a consequence there is a jump of at .
The dependence of on the lengths of CNTs for a fixed value of the temperature is shown in Figure 2 (we mention that for the value of the temperature shown () the LC is in the nematic phase). The critical value of the coupling parameter decreases when the length of CNTs increase, meaning that the continuous paranematic-nematic transition of CNTs is favored.
In Figure 3 the volume fractions of CNTs at the transition as a function of the coupling parameter for a fixed value of the temperature are shown. For negative values of (region I in the figure, perpendicular alignment), the volume fraction gap of CNTs remains practically constant, while for positive values (region II in the figure, parallel alignment), the volume fraction gap decreases with increasing until it cancels for critical value .
The CNT order parameters at the transition as a function of the coupling parameter for a fixed value of the temperature are plotted in Figure 4. For negative values of the coupling parameter (the region I in the figure), the first-order phase transition of CNTs takes place between a paranematic phase with a negative order parameter (perpendicular alignment) and a nematic phase with a positive order parameter (parallel alignment). For positive values of (region II in the figure), the transition takes place between a paranematic phase and a nematic phase, both with positive order parameters (parallel alignment). If the coupling increases the degrees of alignment of the two phases become smaller until the critical value of the coupling parameter when they are identical.
3.2 Phase diagram
The equilibrium phase diagram can be obtained by minimizing the free energy with respect to and :
and equating the chemical potentials of CNTs and LC in the two phases. The chemical potentials are given by
where the free energy density is given by Eq. (5).
3.2.1 Positive coupling constant
A positive coupling parameter corresponds to a paranematic ( very small and positive) and nematic ( large and positive) phase transition of CNTs. The corresponding phase diagram is shown in Figure 5, where the smaller case letters i, p, and n stand for isotropic, paranematic (parallel alignment), and nematic phases, respectively.
For , LC is in the nematic phase, and the paranematic-nematic phase transition of CNTs is first order. With decreasing temperature, the bi-phase region (the Flory “chimney” ) becomes smaller and cancels at the critical point (for the value , the critical temperature is given by ). On the contrary, for , LC is in the isotropic phase. The isotropic-nematic phase transition of CNTs is first order with a constant Flory “chimney” that does not depend on temperature (because , , and the problem is equivalent with the dispersion of CNTs in an isotropic solvent). The triple points are the intersection points between the CNT phase transition equilibrium curves and LC phase transition equilibrium curve (the horizontal lines).
To see in more detail the orientational order of the two components, in Figure 6, we have plotted the order parameter profiles and along the binodal lines as a function of temperature.
The transition between the isotropic and nematic phases of LC is first order (see Figure 6a), and the values of the LC order parameter depend only on temperature and are not influenced by CNT phase transition (because of , the concentration of LC component is much larger than for CNTs). But each particle contributes with to the free energy, so that the free energy of LC is significantly larger than that of CNTs. The order parameters of CNTs as functions of temperature are shown in Figure 6b. Lowering the temperature, the jump of the order parameters cancels at the critical temperature. In the isotropic phase of LC (), the jump of the CNT order parameter remains constant independent on temperature (see the discussion of the results plotted in Figure 5).
3.2.2 Negative coupling constant
The negative coupling parameter corresponds to a paranematic ( very small and negative) and nematic ( large and positive) phase transition of CNTs. The corresponding phase diagram is shown in Figure 7, where the smaller case letters i, p, and n stand for isotropic, paranematic (perpendicular alignment), and nematic phases, respectively.
For , LC is in the nematic phase, and the paranematic (perpendicular)-nematic (parallel) phase transition of CNTs is first order. Contrary to the positive coupling parameter, here with decreasing temperature, the bi-phase region (the Flory “chimney” ) becomes larger. For , the characteristic of the phase diagram remains the same; because LC is in the isotropic phase, the isotropic-nematic phase transition of CNTs is first order with a constant Flory “chimney” that does not depend on temperature.
The order parameter profiles and along the binodal lines as a function of temperature are plotted in Figure 8.
Again the transition between the isotropic and nematic phases of LC is first order (see Figure 8a), and the values of the LC order parameter depend only on temperature. The behavior of CNT orientational order with temperature is shown in Figure 8b. In the nematic phase of LC (), the CNT jump of the order parameters at transition slowly decreases with increasing temperature, while in the isotropic phase of LC (), the jump of the CNT order parameter remains constant independent on temperature.
In the present chapter, we have extended the previous mesoscopic phenomenological model [31, 32, 34, 35, 36] (used to describe the phase behavior of a binary mixture composed of CNTs and thermotropic nematic LC) to include the possibility of a perpendicular alignment of CNTs to LC molecules. The model contains the CNT free energy density (2), the Landau-de Gennes free energy (3) for thermotropic LC order, and an interaction term between the nematic CNTs and nematic LC (4). This interaction term generates the possibility of existence of two nematic phases of CNTs: (i) the mean direction of orientation of CNTs is parallel with the direction of the director of LC molecules, when and , and (ii) the two direction of orientation are perpendicular with each other, and . It is to be mentioned that the order parameter of CNTs can be negative (perpendicular alignment) only for a paranematic phase for which the degree of orientational order is very low. The model predicts in this case a first-order paranematic (perpendicular alignment)-nematic (parallel alignment) phase transition of CNTs. When the coupling constant is positive, the first-order paranematic (parallel alignment)-nematic (parallel alignment) phase transition of CNTs cancels at a critical point, when the transition becomes continuous.
The phase behavior and orientational properties of the mixture are discussed considering different values of temperature , volume fraction of CNTs , and coupling strength .
The theoretical model we have presented describes quite well (comparing with experimental results) the phase properties of the CNTs into thermotropic nematic LC. The model could be improved by considering the polydispersity of the CNT component and by considering, in more detail, the form of the interaction energy between the two components.
Lynch MD, Patrick DL. Organizing carbon nanotubes with liquid crystals. Nano Letters. 2002; 2:1197
Dierking I, Scalia G, Morales P, LeClere D, Adv D. Aligning and reoriented carbon nanotubes with nematic liquid crystals. Maternité. 2004; 16:865-869
Dierking I, San SE. Magnetically steered liquid crystal-nanotube switch. Applied Physics Letters. 2005; 87:233507
Zakri C. Carbon nanotubes and liquid crystalline phases. Liquid Crystal Today. 2007; 16:1
Lagerwall JPF, Scalia G. Carbon nanotubes in liquid crystals. Journal of Materials Chemistry. 2008; 18:2890
Rahman M, Lee W. Scientific duo of carbon nanotubes and nematic liquid crystals. Journal of Physics D: Applied Physics. 2009; 42:063001
Scalia G. Alignment of carbon nanotubes in thermotropic and lyotropic liquid crystals. ChemPhysChem. 2010; 11:333
Lagerwall JPF, Scalia G. A new era for liquid crystal research: Application of liquid crystals in soft matter nano-, bio-, and microtehnology. Current Applied Physics. 2012; 12:1387
Dierking I, Scalia G, Morales P. Liquid crystal- carbon nanotubes dispersions. Journal of Applied Physiology. 2005; 97:044309
Schymura S, Kuhnast M, Lutz V, Jagiella S, Dettiaff-Weglikowska U, Roth S, et al. Towards efficient dispersion of carbon nanotubes in thermotropic liquid crystals. Advanced Functional Materials. 2010; 20:3350
Ji Y, Huang YY, Terentjev EM. Dissolving and alignment carbon nanotubes in thermotropic liquid crystals. Langmuir. 2011; 27:13254
Lebovka NI, Lisetski LN, Goncharuk AI, Minenko SS, Ponevchinsky VV, Soskin MS. Phase transitions in smectogenic liquid crystal BBBA doped by multiwall carbon nanotubes. Phase Transitions. 2013; 86:463
Kumar S. Nanoparticles in the supramolecular order of discotic liquid crystals. Liquid Crystals. 2014; 41:353
Weiss V, Thiruvengadathan R, Regev O. Preparation and characterization of a carbon nanotube- lyotropic liquid crystal composite. Langmuir. 2006; 22:854-856
Lagerwall J, Scalia G, Haluska M, Dettlaff-Weglikowska U, Roth S, Giesselmann F. Nanotubes alignment using lyotropic liquid crystal. Advanced Materials. 2007; 19:359
Schymura S, Enz E, Roth S, Scalia G, Lagerwall JPF. Macroscopic-scale carbon nanotube alignment via self-assembly in lyotropic liquid crystals. Synthetic Metals. 2009; 159:2177
Ould-Moussa N, Blanc C, Zamora-Ledezma C, Lavrentovich OD, Smalyukh II, Islam MF, et al. Dispersion and orientation of single-walled carbon nanotubes in a chromonic liquid crystal. Liquid Crystals. 2013; 40:1628
Duran H, Gazdecki B, Yamashita A, Kyu T. Effect of carbon nanotubes on phase transition of nematic liquid crystals. Liquid Crystals. 2005; 32:815
Chen HY, Lee W, Clark NA. Faster electro-optical response characteristic of a carbon-nanotube-nematic suspension. Applied Physics Letters. 2007; 90:033510
Basu R, Iannacchione GS. Carbon nanotube dispersed liquid crystal: A nano electromechanical system. Applied Physics Letters. 2009; 95:173113
Basu R, Iannacchione GS. Orientational coupling enhancement in a carbon nanotube dispersed liquid crystal. Physical Review E. 2010; 81:051705
Sigdel KP, Iannacchione GS. Effect of carbon nanotubes on the isotropic to nematic and the nematic to smectic-A phase transitions in liquid crystal and carbon nanotubes composites. European Physical Journal E: Soft Matter and Biological Physics. 2011; 34:34
Schymura S, Scalia G. On the effect of carbon nanotubes on properties of liquid crystals. Philosophical Transactions of the Royal Society A. 2013; 371:20120261
Lee KJ, Park HG, Jeong HC, Kim DH, Seo DS, Lee JW, et al. Enhanced electro-optical behavior of a liquid crystal system via multi-walled carbon nanotube doping. Liquid Crystals. 2014; 41:25
Basu R, Iannacchione GS. Dielectric hysteresis, relaxation dynamics, and nonvolatile memory effect in carbon nanotube dispersed liquid crystal. Journal of Applied Physics. 2009; 106:124312
Somoza AM, Sagui C, Roland C. Liquid crystal phases of capped carbon nanotubes. Physical Review B. 2001; 63:081403-1
Sabba Y, Thomas EL. High-concentration dispersion of single-wall carbon nanotube. Macromolecules. 2004; 37:4815
Green MJ, Parra-Vasquez ANG, Behabtu N, Pasquali M. Modelling the phase behavior of polydispersed rigid rods with attractive interaction with application to single-walled carbon nanotubes in superacids. The Journal of Chemical Physics. 2009; 131:084901
Matsuyama A, Kato T. Phase separation in mixture of a liquid crystal and a nanocolloidal particle. Physics Review. 1999; E59:763
Matsuyama A. Theory of binary mixture of a rodlike polymer and a liquid crystal. The Journal of Chemical Physics. 2010; 132:214902
van der Schoot P, Popa-Nita V, Kralj S. Alignment of carbon nanotubes in nematic liquid crystals. The Journal of Physical Chemistry. B. 2008; 112:4512
Popa-Nita V, Kralj S. Liquid crystal-carbon nanotubes mixtures. The Journal of Chemical Physics. 2010; 132:024902
Lisetski LN, Lebovka NI, Naydenov SV, Soskin MS. Dispersions of multi-walled carbon nanotubes in liquid crystals: A physical picture of aggregation. Journal of Molecular Liquids. 2011; 164:143
Popa-Nita V, Barna V, Repnick R, Kralj S. In: Suzuki S, editor. Syntheses and Applications of Carbon Nanotubes and Their Composites. Rijeka, Croatia: InTech; 2013. Chap. 7
Popa-Nita V. Mixtures composed of liquid crystals and carbon nanotubes. The Journal of Chemical Physics. 2014; 140:164905
Popa-Nita V. The phase behavior of rigid rods in an anisotropic mean field with application to carbon nanotubes in nematic liquid crystals. The Journal of Chemical Physics. 2015; 143:094901
Green MJ, Behabtu N, Pasquali M, Adams WW. Nanotubes as polymers. Polymer. 2009; 50:4979
Onsager L. The effect of shape on the interaction of colloidal particles. Annals of the New York Academy of Sciences. 1949; 51:727
de Gennes PG, Prost J. The Physics of Liquid Crystals. Oxford: Oxford University Press; 1993
Oswald P, Pieranski P. Nematic and cholesteric liquid crystals; concepts and physical properties illustrated by experiments. In: Liquid Crystals Book Series. Boca Raton: Taylor and Francis Group/CRC Press; 2005
Doi M. Molecular dynamics and rheological properties of concentrated solutions of rodlike polymer in isotropic and liquid crystalline phases. Journal of Polymer Science Part B: Polymer Physics. 1981; 19:229
Doi M, Edwards SF. Theory of Polymer Dynamics. Oxford: Clarendon; 1989
Kuzuu N, Doi M. Constitutive equations for nematic liquid crystals under weak velocity gradient derived from a molecular kinetic equation. Journal of the Physical Society of Japan. 1983; 52:3486
Park KA, Lee SM, Lee SH, Lee YH. Anchoring a liquid crystal molecule on a single-walled carbon nanotube. Journal of Physical Chemistry C. 2007; 111:1620
Flory PJ. Principles of Polymer Chemistry. Ithaca: Cornell University Press; 1953
Flory PJ. Statistical thermodynamics of semi-flexible chain molecules. Proceedings of the Royal Society of London A. 1956; 243:73
Song W, Kinloch IA, Windle AH. Nematic liquid crystallinity of multiwall carbon nanotubes. Science. 2003; 302:1363
Song W, Windle AH. Isotropic-nematic phase transition of dispersion of multiwall carbon nanotubes. Macromolecules. 2005; 38:6181
Faetti S, Palleschi V. Measurements of the interfacial tension between nematic and isotropic phase of some cyanobiphenyls. The Journal of Chemical Physics. 1984; 81:6254
Faetti S, Palleschi V. Nematic-isotropic interface of some members of the homologous series of 4-cyano-4’-(n-alkyl)biphenyl liquid crystals. Physical Review A. 1984; 30:3241
Faetti S, Palleschi V. Molecular orientation and anchoring energy at the nematic-isotropic. Journal de Physique Lettres (France). 1984; 45:L313