Abstract
We propose a continuum theory of orientational phase transitions induced by an external magnetic field in a suspension of carbon nanotubes in a nematic liquid crystal. It is shown that in a magnetic field a nonuniform and two different uniform phases are possible in the suspension. The uniform phases of the suspension differ by the type of orientational coupling of nanotubes with the liquid crystal matrix (the planar type when the nanotubes are oriented along the matrix director, and the homeotropic type when the nanotubes are perpendicular to the director). The possibility of a redistribution of the nanotube concentration (segregation effect) is shown. The fields of orientational transitions between uniform and nonuniform phases of the suspension are found analytically. It is shown that, when the nanotubes are weakly coupled to the matrix, the magnetic field induces reentrant transitions (uniform planar phase–nonuniform phase–uniform homeotropic phase–nonuniform phase). These transitions can be of first or of second order depending on the carbon nanotubes segregation intensity.
Introduction
In recent years suspensions of anisometric particles in liquid crystals have become of great interest for researchers [1]. This is not only because liquid crystals (LCs) have found wide application in modern optoelectronic devices [2] but also because anisometric particles orient like a LC in a medium that is capable of spontaneous orientational ordering. Examples of such media are suspensions of ferromagnetic or ferroelectric particles, as well as of carbon nanotubes. In 1970, Brochard and de Gennes proposed to dope nematic liquid crystals (NLCs) with elongated ferromagnetic particles [3]. The magnetic susceptibility of such composite system (called ferronematic), turned out to be several orders of magnitude higher than that of the pure LC even for a low concentration of the dispersed phase (0.01 vol %). The new approach has opened the way for the creation of devices that operate on the basis of magnetic orientational transitions in the LC. To date, many theoretical and experimental studies on the properties of ferronematics have been published [4,5], which indicates the interest in this kind of composite materials.
Along with ferri or ferromagnetic particles it is also possible to use carbon nanotubes (CNTs) in order to increase the magnetoorientational response of the LC matrix [6,7]. Due to the highly elongated shape (aspect ratios of 10^{2} to 10^{3}) and anomalously high anisotropy of the diamagnetic susceptibility ( ≈ 10^{−5} to 10^{−4}) [811], CNTs are very attractive for the creation of nanocomposites based on LCs with high magnetoorientational response. From experimental data [1216], it is known that in the absence of external fields CNTs are oriented parallel to the director of the LC matrix, which corresponds to the planar type of coupling. However, homeotropic coupling is also possible [17]. Thus, for suspensions of CNTs based on LC with positive diamagnetic susceptibility anisotropies, one should expect decreasing of the threshold field of the magnetic Fréedericksz transition, which is confirmed by experiment [18,19], as well as a decrease in the electric field of the Fréedericksz transition [11,14,2022]. Along with this, there are experimental studies devoted to the investigation of LC suspensions with CNTs functionalized by ferromagnetic particles [19,2325], where an enhanced magnetooptical response is also observed in comparison with a pure LC.
The available theoretical approaches to the description of CNT suspensions in LCs are based on a generalization of the Landau–de Gennes theory [2629] and meanfield theory [30,31]. In these papers, the phase state of the binary mixture of CNTs in LC as a function of the concentration of CNTs, the coupling energy of subsystems and temperature in the absence of external fields were studied.
In the present paper we propose a continuum theory of dilute CNT suspension in LCs that makes it possible to study orientational transitions induced by the magnetic field.
Results and Discussion
Basic equations
We consider infinite plane layer of thickness L of a suspension of CNTs in NLC with planar texture and an absolutely rigid anchoring of NLC molecules with the boundaries. We set the origin of the coordinate system at the middle of the layer (see Figure 1). We use the unit vectors n and m, the socalled directors, to describe the preferential orientation of the LC molecules and CNTs, respectively. We assume soft and planar coupling of LC with the CNT surface, then in the absence of a magnetic field . In the case of a positive anisotropy of the diamagnetic susceptibility of the LC χ_{a} > 0, applying a magnetic field H = (0,0,H) will lead to the appearance of distortions in the orientational structure of the LC. This effect is known as the Fréedericksz transition [32]. CNTs are also oriented by the magnetic field, because they possess an anomalously strong diamagnetism and positive diamagnetic anisotropy [10,3335], so even small concentrations of CNTs in suspension should lead to a decreasing of the Fréedericksz transition threshold in comparison with the pure NLC. These two mechanisms of the field influence on the suspension prove to be interdependent due to the orientational interaction of CNTs with the LC matrix.
In the framework of continuum theory, the equilibrium state of the suspension corresponds to the minimum of free energy
Here, K_{11}, K_{22} and K_{33} are the Frank elastic moduli; χ_{a} and are the diamagnetic susceptibility anisotropies of LC and CNTs, respectively; μ_{0} is the permeability of vacuum; f is the volume fraction of CNTs in the suspension; W_{p} is the surface density of the coupling energy between the LC molecules and the surface of the CNTs; d is the transverse diameter of a CNT; ν is the volume of a CNT; k_{B} is the Boltzmann constant; T is the temperature. We assume a low concentration of CNTs in LC, ( = Nν/V, N is the number of CNTs in the suspension, V is the suspension volume). This allows us to neglect the interaction between the CNTs.
The term F_{1} in Equation 1 is the free energy density of elastic deformations of the LC [36]. The contributions F_{2} and F_{3} take into account the interaction energies of the LC matrix and CNTs with the magnetic field. The term F_{4} describes the interaction of elongated impurity particles (in our case CNTs) with LC molecules [37]. In the absence of a field and for W_{p} > 0, the minimum of F_{4} corresponds to , i.e., a planar coupling of the LC matrix and CNTs. The planar type of coupling was observed in most of the experimentally studied suspensions of CNTs in LCs [14,18,38]. The last term F_{5} is the contribution of the entropy of mixing of an ideal solution of CNTs in the LC matrix.
As it is known, the state of thermodynamic equilibrium corresponds to the minimum of the free energy (Equation 1), which is a functional of the two vectors n and m, and the scalar quantity f. Thus, the problem is reduced to deriving the equilibrium configurations of the directors n(r) and m(r), and the volume fraction f(r) of the CNTs.
In the considered geometry, the vectors n and m conveniently have the following form:
We choose the thickness L of the layer as the unit of length and the quantity as the unit of the magnetic field strength, which corresponds to magnetic Fréedericksz transition field in pure NLC. We introduce the following dimensionless parameters [39]:
Here, ζ is the dimensionless coordinate and h is the dimensionless magnetic field strength. Due to the coupling between the directors of LC, n, and CNT, m, (term F_{4} in Equation 1), the suspension possesses two orientation mechanisms under the influence of a magnetic field. Both mechanisms are quadrupole in nature and are caused by the diamagnetism of the LC matrix and the CNTs. The parameter γ, which is determined by the balance of the terms F_{3} and F_{2} of the free energy (Equation 1), characterizes the relative contribution of the mechanisms of magnetic field influence on the orientational structure of the suspension. For , the appearance of orientational distortions of the director field are caused mainly by the diamagnetism of the CNTs, and for , the distortions of orientational structure arise mainly from the diamagnetism of the LC matrix. The parameter σ characterizes the coupling energy of the LC and the impurity subsystem; g is the reduced volume fraction of CNTs in the suspension, and k is the ratio of Frank's constants of the LC.
The parameter κ is the square of the ratio between two characteristic lengths, i.e., layer thickness L and segregation length [3,40]. The characteristic size of the concentration redistribution region, L_{S}, can be determined from the balance of the contributions F_{1} and F_{5} in the volume density of the free energy (Equation 1), which allows us to introduce the dimensionless segregation parameter . For , the segregation effect is negligible because the characteristic scale of the segregation region of CNTs considerably exceeds the thickness of the layer; for κ ≤ 1, the segregation effect becomes significant.
We take typical values for NLCs [36] to estimate dimensionless quantities: χ_{a} ≈ 10^{−6} and K_{33} > K_{11} ≈ 10^{−12} N. For CNTs we can assume [711]: ≈ 10^{−5} to 10^{−4}, d ≈ 10^{−8} m, CNT length l ≈ 10^{−6} m and ν ≈ 10^{−22} m^{−3}. We also set T = 300 K and L = 20 μm. For different suspensions the coupling energy of CNTs with the LC matrix, W_{p}, varies over a wide range from 10^{−7} N·m^{−1} [26] to 10^{−3} N·m^{−1} [24]. Assuming a volume fraction of CNTs ≈ 10^{−3} [7], we obtain κ ≈ 1, γ ≈ 10^{−2} to 10^{−1}, σ ≈ 1 to 10^{3} and k ≈ 1.
The dimensionless free energy takes the following form after the substitution of Equation 2 into Equation 1:
Here, S is the surface area of the planes confining the layer, and the notation K(φ) = cos^{2} φ + k sin^{2} φ is introduced.
Minimizing the functional (Equation 4) with respect to the functions φ(ζ), ψ(ζ) and g(ζ), we obtain a system of equations describing the orientational equilibrium of the suspension:
Here, the quantity Q can be derived from the condition of a constant number of CNTs in the suspension,
Together with the condition of a rigid planar coupling of the LC director to the boundaries of the layer,
the system of Equations 5–8 allows us to find the equilibrium state of the CNT suspension in a magnetic field.
It follows from Equation 7 that in a suspension with nonuniform distribution of the directors n(r) and m(r), even a uniform magnetic field causes a redistribution of the CNT concentration f(r) in the layer, known as segregation effect [3] in the physics of ferronematics. This means that CNTs migrate to those regions of the layer where their magnetic energy in the field () and the orientational energy in the LC matrix () are minimal.
Orientational phases of the suspension
The system of Equations 5–8 allows for uniform solutions [g(ζ) = 1 and the angles φ and ψ are independent of coordinates]. One of them [φ(ζ) = ψ(ζ) ≡ 0] corresponds to the initial state , in which the long CNTs axes are parallel to the director of the LC. We call this state the planar phase of the suspension, since it is characterized by the planar coupling () of CNTs to the LC matrix. In this state, the directors of the LC and CNTs are directed orthogonally to the external magnetic field and parallel to the boundaries of the layer. Another uniform solution [φ(ζ) = 0 and ψ(ζ) = π/2] corresponds to the homeotropic phase where the LC director is parallel to the layer boundaries and the CNT director is oriented along the field, . This phase is characterized by the homeotropic coupling () of CNTs to the LC matrix. The angular phase [41] corresponds to the nonuniform solution [φ = φ(ζ), ψ = ψ(ζ) and g = g(ζ)]. In this state the angle between the directors n and m is different from zero and π/2. Schematic representations of the planar, angular and homeotropic phases of the suspension are shown in Figure 2.
Under the influence of a magnetic field, the initial planar phase of the suspension (Figure 2a) becomes unstable. When the field reaches the threshold value h_{c}, a transition to the angular phase occurs (Figure 2b) with distorted orientational structures of n and m. By analogy with pure LCs this transition can be called the Fréedericksz transition. The distortions of the orientational structure are small [, ] in the vicinity of h_{c}, and the distribution of CNTs in the layer is almost homogeneous, g(ζ) ≈ 1. Therefore, the solution of the system of Equations 5–8 can be found analytically. In the lowest order we obtain φ(ζ) = φ_{m}·cos(πζ) and ψ(ζ) = λ_{c}φ(ζ). Here, the value of corresponds to the orientation angle of the LC director in the middle of the layer, and λ_{c} is determined by the relation
Nonzero solutions of φ(ζ) exist for h ≥ h_{c}, where h_{c} has the meaning of the Fréedericksz transition field from the planar to the angular phase and it is found from the following equation:
The thermodynamically stable solution of Equation 12 has the form
and it is shown in Figure 3.
In the case of strong coupling between CNTs and LC matrix (), Equation 13 for the Fréedericksz transition field in the lowest order in the small parameter 1/σ gives
Here, the quantity corresponds to the field of the Fréedericksz transition in the case of absolutely rigid (σ→∞) coupling of LC and impurity subsystem, i.e., of the directors n and m. It follows from Equation 13 that with increasing parameter γ (i.e., as the volume fraction of CNTs or the anisotropy of their diamagnetic susceptibility increases), the Fréedericksz field decreases. Pure NLC with corresponds to the value γ = 0.
When the CNTs is weakly coupled to the LC matrix in lowest order in small σ, from Equation 13 we obtain
It follows from Equations 12–15 that for suspensions with planar coupling of the directors n and m, the Fréedericksz field is always smaller than for a pure NLC, which is confirmed by experiments [18,19,23].
We now determine the transition field between the nonuniform (angular) phase of the suspension and the uniform phase with the homeotropic coupling of CNTs with the LC matrix. In the homeotropic phase [φ(ζ) = 0 and ψ(ζ) = π/2] the LC director is parallel to the layer boundaries and the CNT director is oriented along the magnetic field (Figure 2c). In the vicinity of the transition field h_{r} between the angular and homeotropic phases, the deviations of the director of the LC, n, from the xaxis and of the CNT director, m, from the magnetic field direction H are small. Therefore in the lowest order in small φ(ζ) and π/2 − ψ(ζ), we obtain the following equation for determining the transition field between the angular and the homeotropic phase:
Equation 16 can be solved with respect to h_{r} and we obtain
The solution h_{r}_{+} describes the upper branch of the doublevalued curve h_{r}(σ) in Figure 3 and h_{r}_{−} corresponds to the lower branch.
In the case of weak coupling of the CNT and LC directors () the expressions in Equation 17 can be represented in the form
It is seen from Equation 18 that in fields the suspension can only be in the angular phase.
We note that the existence of the homeotropic phase is not possible for every value of the coupling energy of the LC matrix and CNTs, but only for σ ≤ σ_{m}, where the threshold value of the coupling energy (Figure 3) is determined by the relation
As can be seen from Equation 19, the region of existence of the homeotropic phase widens with increasing parameter γ. After substituting Equation 19 in Equation 16, we obtain the expression for the transition field h_{r} corresponding to the coupling energy σ_{m}:
As can be seen from Equation 19, the upper bound of σ_{m} is π^{2}/2 (γ→∞). Therefore, for suspensions with the coupling energies σ > π^{2}/2, the homeotropic phase cannot exist. For σ > π^{2}/2, the coupling of CNTs with LC is almost absolutely rigid. In this case the description of the orientational structure of the suspension becomes possible with the help of one director n (n = m). It is seen from the expression for the free energy (Equation 1) that in this case the suspension behaves like a pure LC with an effective anisotropy of the diamagnetic susceptibility, .
The diagram of the orientational state of the suspension, i.e., the threshold fields (Equation 13) and (Equation 16) for orientational transitions as functions of the coupling energy of CNTs with the LC matrix and for different values of γ is shown in Figure 3. The region bounded by the abscissa and curve h_{c} corresponds to the planar phase in which the directors n and m are parallel to the boundaries of the layer (Figure 2a). The region bounded by the ordinate axis and the curve h_{r} corresponds to the homeotropic phase in which the LC director is oriented along the layer boundaries and the CNT director is parallel to the magnetic field, (Figure 2c). Outside these regions is the angular phase with nonuniform distributions of n and m over the thickness of the layer (Figure 2b).
Figure 3 shows that if the coupling of CNTs with the LC matrix is weak (σ < σ_{m}), there is a sequence of transitions with increasing field for a given value of the coupling energy: planar phase–angular phase–homeotropic phase–angular phase. That is, there are reentrant phenomena. Such reentrant transitions occur only for weak coupling of CNTs with the LC matrix. A magnetic field directed perpendicularly to the CNTs makes the initial alignment of tubes energetically unfavorable, so that they begin to rotate in the field direction minimizing the contribution F_{3} in Equation 1. The LC molecules also tend to rotate into the field direction minimizing the contribution F_{2}. However, the orienting influence of the cell boundaries prevents a rotation of the LC director. Due to the orientational coupling between the LC and CNTs, the director rotation of the disperse phase (the term F_{4} in Equation 1) is transmitted to the LC matrix, and the first orientational transition from the initial planar phase to the nonuniform angular phase occurs at h = h_{c}. The Fréedericksz transition threshold field, h_{c}, appears to be smaller than that of the pure NLC. This result is in good agreement with experimental observations [18,19,23]. Due to induced gradients of the LC director and the increase in the energy of orientationally elastic deformations F_{1} in Equation 1, the forces of orientational elasticity tend to return the director to the initial planar state. Distortions of the orientational structure also cause the segregation of the CNTs, and the term F_{5} in Equation 1 is minimal when the impurity is uniformly distributed over the sample. Distortions of the LC director are induced in this case by the orientational coupling between n and m (contribution F_{4}), i.e., the CNTs director entrains the LC director, and F_{4} tends to decrease the angle between n and m. For this reason, the appearance of large orientational deformations (i.e., large gradients of n) is energetically less favorable than a minimization of F_{3}, leading to the occurrence of the uniform homeotropic phase with . So that the next transition to the uniform homeotropic phase occurs with an increase of the magnetic field to h = h_{r}_{−}. For h_{r}_{−} < h < h_{r}_{+}, the small energy loss (weak coupling) in the contribution F_{4} and F_{2} (the field is significantly lower than the Fréedericksz field of a pure LC) is compensated by the gains in the energy of orientational elastic deformations (F_{1} = 0), in the entropy contribution to the energy (F_{5} = min), and in the magnetic energy of CNTs (F_{3} = min). When the diamagnetism of the LC matrix (F_{2}) begins to predominate, at h ≈ h_{r}_{+}, the next transition from the homeotropic phase to the nonuniform state (the angular phase) occurs. A similar diagram was described in [39] for ferronematic liquid crystals. In the case of strong coupling (σ > σ_{m}), the initial planar phase undergoes a Fréedericksz transition to the nonuniform angular phase with increasing field. There is a “synchronous” rotation of the LC and CNT directors along the applied field, i.e., a monotonic increase in the deviations of the directors, since the energy loss in the contribution F_{4} is not small.
Tricritical phenomena
We now determine the character of the orientational transitions. In the vicinity of the Fréedericksz field, h_{c}, the deviations of the LC and CNT directors from the boundaries of the layer are small (small φ and ψ), so the free energy of the suspension (Equation 4) can be expanded in a power series over small φ(ζ) = φ_{m}·cos(πζ) and ψ(ζ) = λ_{c}φ(ζ). Here, the value corresponds to the orientation angle of the LC director in the middle of the layer, and λ_{c} is defined by Equation 11. After integration in the fourth order in φ_{m}, the free energy takes the form of the Landau expansion
where
Here, F_{c} is the free energy of the uniform planar phase.
The dependence of the LC director orientation angle in the middle of the layer, φ_{m}, on the applied magnetic field in the vicinity of h_{c} can be found by minimizing the expression in Equation 21 over φ_{m}:
As noted above, for γ > 0 we always have h_{c} < π, and therefore the coefficients of the expansion α_{c} > 0, and . Then, as can be seen from Equation 23, the Fréedericksz transition between the planar and angular phases is a secondorder transition (h ≥ h_{c}) for . If , real solutions of Equation 23 exist only for h ≤ h_{c}, and the Fréedericksz transition is a firstorder transition. The quantity is the tricritical value of the segregation parameter below which the character of the Fréedericksz transition changes from second order to first order.
The dependence of on the coupling energy of CNTs with the LC matrix for different values of γ is shown in Figure 4. For γ = 0.1 (curve 1) behaves nonmonotonically with increasing coupling energy σ and exhibits a pronounced maximum, while for γ = 0.5 (curve 2) the maximum is poorly distinguishable. For σ > 10, slightly varies for a given γ, and in the limiting case of rigid coupling of CNT and LC directors, it tends to the value of π^{2}γ^{2}/4k(1 + γ)^{2}. It is seen from Figure 4a for σ > 10, grows with increasing γ. For small coupling energies (Figure 4b), grows faster with increasing σ for γ = 0.1 than for γ = 0.5.
We now consider the transition from the angular to the homeotropic phase to which the field h_{r}_{−} corresponds (see Figure 3). Near the transition point h_{r}_{−}, deviations of the LC director from the boundaries of the layer and of the CNT director from the direction of the magnetic field are small. Therefore, the free energy (Equation 4) can be represented in analogy with the previous case in the form of a Landau expansion:
where
Here, F_{r} corresponds to the free energy of the uniform homeotropic phase.
After minimizing Equation 24 with respect to φ_{m}, we obtain the dependence of the LC director orientation angle in the middle of the layer on the magnetic field strength in the vicinity of h_{r}_{−}:
As can be seen from this formula, the considered transition is a secondorder transition if real solutions of Equation 26 exist for fields h ≤ h_{r}_{−}. This is only possible for , since the expansion coefficients α_{−} and are positive for h_{r} = h_{r}_{−} < h_{m}. For , the transition from the angular phase to the homeotropic one is a firstorder transition and the segregation parameter corresponds to the tricritical point.
Now we consider the possibility of changing the character of the orientational transition from the homeotropic to the angular phase, to which the field h_{r}_{+} corresponds (see Figure 3). The Landau expansion of the free energy (Equation 4) in the vicinity of h_{r}_{+} has the form
Here
Minimizing Equation 25 with respect to φ_{m} we obtain
It can be seen that the transition of the suspension from the homeotropic to the angular phase with increasing magnetic field must be a secondorder transition when real solutions of Equation 29 exist for h ≥ h_{r}_{+} (see Figure 3), i.e., since α_{+} > 0 and . For the transition from the homeotropic to the angular phase is a firstorder transition. The quantity corresponds to the tricritical value of the segregation parameter κ.
In Figure 5, the tricritical segregation parameters and are given as functions of the coupling energy of CNTs with the LC matrix, σ, for different values of the parameter γ. grows with increasing σ and γ (Figure 5a). The dependence of on σ and γ is analogous to that of (Figure 5b) with the exception of the small region near σ ≈ σ_{m} (right edge of the curves), where the decreases with growing σ.
Equation 12 and Equation 16 allow us to determine only the threshold fields of the secondorder transitions. To find the fields of the firstorder equilibrium transitions, it is necessary to use the conditions for the equality of the free energy of the nonuniform state (Equation 4) and the free energy of the planar phase, F_{c}, or of the homeotropic phase, F_{r}, respectively.
Conclusion
The orientational phase transitions in a suspension of CNTs in an NLC in a magnetic field perpendicular to the boundaries of the layer are studied. The possible phase states of the suspension in the magnetic field are established: a uniform phase with planar coupling of CNTs with the LC matrix, a nonuniform angular phase and a uniform phase with homeotropic coupling. Transitions between these phases occur when the magnetic field reaches certain threshold values. Analytical expressions that determine the fields of transitions between the uniform planar and nonuniform angular phases, as well as between the angular and the uniform homeotropic phases are found. It is shown that orientational transitions in the suspension occur at field strengths smaller than the Fréedericksz transition field of the pure NLC.
Furthermore, it is shown that when the CNTs are weakly coupled to the LC matrix reentrant transitions, planar phase – angular phase – homeotropic phase – angular phase, are possible under an increasing magnetic field. In the case of strong coupling the initial planar phase undergoes a Fréedericksz transition to the nonuniform angular phase with a growing magnetic field.
It is also shown that segregation takes place in the suspension of CNTs in the LC, i.e., under the influence of the magnetic field the CNTs are redistributed over the layer, so that the CNT concentration increases in those regions of the layer where the sum of their magnetic energy and the orientation energy of coupling to the LC matrix is minimum. It is found that, depending on the intensity of the segregation effects associated with the redistribution of CNTs over the layer thickness, all orientational transitions exhibit tricritical behavior, i.e., they can be transitions of both first or secondorder. The analytic expressions for the tricritical values of the segregation parameter are obtained.
The present paper is an extension to several studies of LC suspensions doped with dipole (ferromagnetic) particles [39,4144]. Here we discuss the similarity and difference in the orientational behavior of dipole (ferromagnetic) and quadrupole (diamagnetic) particles in LC suspensions.
In the present paper we consider a physical system that, unlike ferronematics, does not have a dipole response to the applied magnetic field. We study LC suspensions with anisometric diamagnetic particles. Thus, the present theory describes the behavior of quadrupolar particles embedded in an LC. We analyze another possible way to enhance the anisotropy of diamagnetic susceptibility of LCs by doping them with diamagnetic CNTs.
In the magnetized ferronematics (in which the magnetic moments of the ferroparticles are aligned in one direction) with planar coupling of impurity particles to the LC matrix, the magnetic field induces nonthreshold Fréedericksz transition for suspensions based on NLC with negative [42] or positive [39] anisotropy of the diamagnetic susceptibility. As we have shown above, in LC suspensions of CNTs this transition has threshold behavior. The next general feature of all considered systems is the possibility of existence of uniform homeotropic phases. In ferronematics with negative diamagnetic anisotropy of the LC matrix, the homeotropic phase remains stable with increasing field, while for positive diamagnetic anisotropy the transition from the homeotropic to the nonuniform angular phase occurs like for the CNT suspension. For all kinds of suspensions, the transitions between the angular and the homeotropic phases can be of first or of second order depending on the segregation intensity.
The socalled compensated ferronematics, which have equiprobable distributions of the particles parallel and antiparallel to the LC director in the absence of a field [43], show a magnetoorientational response which is similar to that of CNT suspensions in LCs. The compensated ferronematics exhibit a quadrupole response to the magnetic field applied perpendicularly to the initial alignment of the magnetic subsystems. They exhibit a threshold Fréedericksz transition from the initial compensated phase (planar phase) to the nonuniform angular phase, like the CNT suspensions. For weak coupling of magnetic particles and LC matrix the following sequence of reentrant transitions takes place: initial uniform compensated phase (planar phase) – nonuniform phase (angular phase) – uniform saturation phase (homeotropic phase) – nonuniform phase (planar phase). For strong coupling of ferroparticles with the LC matrix there is an analogous response of the LC director and magnetization to the external magnetic field, and the sequence of transitions is: uniform compensated phase (planar phase) – nonuniform phase. There are some very important differences between compensated ferronematics and the suspension of CNTs in the LC. For compensated ferronematics the Fréedericksz transition can be only a secondorder transition, while for the CNT suspension this transition can be of first or of second order, depending on the segregation parameter. The other feature is that the transition fields for compensated ferronematics depend on the segregation parameter, while for magnetized ferronematics and for LC suspensions of CNTs there is no such dependence [39,42,43].
In [44] we studied the orientational response of a magnetized ferronematic liquid crystal to magnetic and electric fields. In contrast to the abovementioned works, we considered the bistable coupling between the particles and LC matrix. It is shown that apart from magnetic impurity segregation, the firstorder orientational transitions can be due to the bistable orientational coupling.
Another approach for studying the LC suspensions of the dipolar (ferromagnetic) particles was proposed in [45,46]. This approach is based on the meanfield theory, and allows for an investigation of the influence of temperature and magnetic field on a suspension, including the phase transition from the ordered phase into the nematic or paranematic state. We plan to propose such a theory for LCs doped with diamagnetic particles.
References

Lagerwall, J. P. F.; Scalia, G., Eds. Soft Condensed Matter; World Scientific Publishing Co Pte Ltd: Singapore, 2017; Vol. 7. doi:10.1142/9789814619264_fmatter
Return to citation in text: [1] 
Lagerwall, J. P. F.; Scalia, G. Curr. Appl. Phys. 2012, 12, 1387–1412. doi:10.1016/j.cap.2012.03.019
Return to citation in text: [1] 
Brochard, F.; de Gennes, P. G. J. Phys. (Paris) 1970, 31, 691–708. doi:10.1051/jphys:01970003107069100
Return to citation in text: [1] [2] [3] 
Garbovskiy, Y. A.; Glushchenko, A. V. Solid State Phys. 2010, 62, 1–74. doi:10.1016/B9780123742933.000018
Return to citation in text: [1] 
Mertelj, A.; Lisjak, D. Liq. Cryst. Rev. 2017, 5, 1–33. doi:10.1080/21680396.2017.1304835
Return to citation in text: [1] 
Lagerwall, J. P. F.; Scalia, G. J. Mater. Chem. 2008, 18, 2890–2898. doi:10.1039/b802707b
Return to citation in text: [1] 
Yadav, S. P.; Singh, S. Prog. Mater. Sci. 2016, 80, 38–76. doi:10.1016/j.pmatsci.2015.12.002
Return to citation in text: [1] [2] [3] 
Ramirez, A. P.; Haddon, R. C.; Zhou, O.; Fleming, R. M.; Zhang, J.; McClure, S. M.; Smalley, R. E. Science 1994, 265, 84–86. doi:10.1126/science.265.5168.84
Return to citation in text: [1] [2] 
Fujiwara, M.; Kawakami, K.; Tanimoto, Y. Mol. Phys. 2002, 100, 1085–1088. doi:10.1080/00268970110109501
Return to citation in text: [1] [2] 
Zaric, S.; Ostojic, G. N.; Kono, J.; Shaver, J.; Moore, V. C.; Hauge, R. H.; Smalley, R. E.; Wei, X. Nano Lett. 2004, 4, 2219–2221. doi:10.1021/nl0486012
Return to citation in text: [1] [2] [3] 
Cîrtoaje, C.; Petrescu, E. Physica E 2016, 84, 244–248. doi:10.1016/j.physe.2016.06.011
Return to citation in text: [1] [2] [3] 
Lynch, M. D.; Patrick, D. L. Nano Lett. 2002, 2, 1197–1201. doi:10.1021/nl025694j
Return to citation in text: [1] 
Dierking, I.; Scalia, G.; Morales, P.; LeClere, D. Adv. Mater. 2004, 16, 865–869. doi:10.1002/adma.200306196
Return to citation in text: [1] 
Dierking, I.; Scalia, G.; Morales, P. J. Appl. Phys. 2005, 97, 044309. doi:10.1063/1.1850606
Return to citation in text: [1] [2] [3] 
Jeong, H. S.; Youn, S. C.; Kim, Y. H.; Jung, H.T. Phys. Chem. Chem. Phys. 2013, 15, 9493–9497. doi:10.1039/c3cp00052d
Return to citation in text: [1] 
Scalia, G. ChemPhysChem 2010, 11, 333–340. doi:10.1002/cphc.200900747
Return to citation in text: [1] 
Galerne, Y. Phys. Rev. E 2016, 93, 042702. doi:10.1103/PhysRevE.93.042702
Return to citation in text: [1] 
Dierking, I.; San, S. E. Appl. Phys. Lett. 2005, 87, 233507. doi:10.1063/1.2140069
Return to citation in text: [1] [2] [3] [4] 
Tomašovičová, N.; Timko, M.; Mitróová, Z.; Koneracká, M.; Rajňak, M.; Éber, N.; TóthKatona, T.; Chaud, X.; Jadzyn, J.; Kopčanský, P. Phys. Rev. E 2013, 87, 014501. doi:10.1103/PhysRevE.87.014501
Return to citation in text: [1] [2] [3] [4] 
Chen, H.Y.; Lee, W. Appl. Phys. Lett. 2006, 88, 222105. doi:10.1063/1.2208373
Return to citation in text: [1] 
Chen, H.Y.; Lee, W. Opt. Rev. 2005, 12, 223–225. doi:10.1007/s1004300502237
Return to citation in text: [1] 
Cîrtoaje, C.; Stoian, V.; Petrescu, E.; Moţoc, C. Proc. SPIE 2015, 9517, 95171G. doi:10.1117/12.2178333
Return to citation in text: [1] 
Buluy, O.; Nepijko, S.; Reshetnyak, V.; Ouskova, E.; Zadorozhnii, V.; Leonhardt, A.; Ritschel, M.; Schönhense, G.; Reznikov, Y. Soft Matter 2011, 7, 644–649. doi:10.1039/c0sm00131g
Return to citation in text: [1] [2] [3] 
Mitróová, Z.; Tomašovičová, N.; Timko, M.; Koneracká, M.; Kováč, J.; Jadzyn, J.; Vávra, I.; Éber, N.; TóthKatona, T.; Beaugnon, E.; Chaud, X.; Kopčanský, P. New J. Chem. 2011, 35, 1260–1264. doi:10.1039/c1nj20017h
Return to citation in text: [1] [2] 
Yoo, H. J.; Lee, S. Y.; You, N.H.; Lee, D. S.; Yeo, H.; Choi, Y. M.; Goh, M.; Park, J.; Akagi, K.; Cho, J. W. Synth. Met. 2013, 181, 10–17. doi:10.1016/j.synthmet.2013.07.023
Return to citation in text: [1] 
van der Schoot, P.; PopaNita, V.; Kralj, S. J. Phys. Chem. B 2008, 112, 4512–4518. doi:10.1021/jp712173n
Return to citation in text: [1] [2] 
PopaNita, V.; Kralj, S. J. Chem. Phys. 2010, 132, 024902. doi:10.1063/1.3291078
Return to citation in text: [1] 
PopaNita, V. J. Chem. Phys. 2014, 140, 164905. doi:10.1063/1.4872241
Return to citation in text: [1] 
PopaNita, V. J. Chem. Phys. 2015, 143, 094901. doi:10.1063/1.4928655
Return to citation in text: [1] 
Matsuyama, A.; Kato, T. Phys. Rev. E 1999, 59, 763–770. doi:10.1103/PhysRevE.59.763
Return to citation in text: [1] 
Matsuyama, A. J. Chem. Phys. 2010, 132, 214902. doi:10.1063/1.3447892
Return to citation in text: [1] 
de Gennes, P. G.; Prost, J. The Physics of Liquid Crystals, 2nd ed.; Oxford Scientific: London, United Kingdom, 1993.
Return to citation in text: [1] 
Byszewski, P.; Baran, M. Europhys. Lett. 1995, 31, 363–366. doi:10.1209/02955075/31/7/004
Return to citation in text: [1] 
Tsui, F.; Jin, L.; Zhou, O. Appl. Phys. Lett. 2000, 76, 1452–1454. doi:10.1063/1.126061
Return to citation in text: [1] 
Searles, T. A.; Imanaka, Y.; Takamasu, T.; Ajiki, H.; Fagan, J. A.; Hobbie, E. K.; Kono, J. Phys. Rev. Lett. 2010, 105, 017403. doi:10.1103/PhysRevLett.105.017403
Return to citation in text: [1] 
Blinov, L. M.; Chigrinov, V. G. Electrooptic Effects in Liquid Crystal Materials; Springer: New York, NY, U.S.A., 1994. doi:10.1007/9781461226925
Return to citation in text: [1] [2] 
Burylov, S. V.; Raikher, Y. L. Mol. Cryst. Liq. Cryst. 1995, 258, 123–141. doi:10.1080/10587259508034553
Return to citation in text: [1] 
Agha, H.; Galerne, Y. Phys. Rev. E 2016, 93, 042703. doi:10.1103/PhysRevE.93.042703
Return to citation in text: [1] 
Zakhlevnykh, A. N.; Petrov, D. A. J. Mol. Liq. 2014, 198, 223–233. doi:10.1016/j.molliq.2014.06.028
Return to citation in text: [1] [2] [3] [4] [5] 
Zakhlevnykh, A. N.; Petrov, D. A. Tech. Phys. 2012, 57, 1208–1218. doi:10.1134/S1063784212090265
Return to citation in text: [1] 
Zakhlevnykh, A. N. J. Magn. Magn. Mater. 2004, 269, 238–244. doi:10.1016/S03048853(03)006127
Return to citation in text: [1] [2] 
Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2015, 393, 517–525. doi:10.1016/j.jmmm.2015.06.021
Return to citation in text: [1] [2] [3] 
Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2016, 401, 188–195. doi:10.1016/j.jmmm.2015.10.025
Return to citation in text: [1] [2] [3] 
Zakhlevnykh, A. N.; Petrov, D. A. J. Exp. Theor. Phys. 2016, 123, 687–698. doi:10.1134/S1063776116090168
Return to citation in text: [1] [2] 
Zakhlevnykh, A. N.; Lubnin, M. S.; Petrov, D. A. J. Exp. Theor. Phys. 2016, 123, 908–917. doi:10.1134/S1063776116100101
Return to citation in text: [1] 
Zakhlevnykh, A. N.; Lubnin, M. S.; Petrov, D. A. J. Magn. Magn. Mater. 2017, 431, 62–65. doi:10.1016/j.jmmm.2016.09.044
Return to citation in text: [1]
18.  Dierking, I.; San, S. E. Appl. Phys. Lett. 2005, 87, 233507. doi:10.1063/1.2140069 
19.  Tomašovičová, N.; Timko, M.; Mitróová, Z.; Koneracká, M.; Rajňak, M.; Éber, N.; TóthKatona, T.; Chaud, X.; Jadzyn, J.; Kopčanský, P. Phys. Rev. E 2013, 87, 014501. doi:10.1103/PhysRevE.87.014501 
23.  Buluy, O.; Nepijko, S.; Reshetnyak, V.; Ouskova, E.; Zadorozhnii, V.; Leonhardt, A.; Ritschel, M.; Schönhense, G.; Reznikov, Y. Soft Matter 2011, 7, 644–649. doi:10.1039/c0sm00131g 
18.  Dierking, I.; San, S. E. Appl. Phys. Lett. 2005, 87, 233507. doi:10.1063/1.2140069 
19.  Tomašovičová, N.; Timko, M.; Mitróová, Z.; Koneracká, M.; Rajňak, M.; Éber, N.; TóthKatona, T.; Chaud, X.; Jadzyn, J.; Kopčanský, P. Phys. Rev. E 2013, 87, 014501. doi:10.1103/PhysRevE.87.014501 
23.  Buluy, O.; Nepijko, S.; Reshetnyak, V.; Ouskova, E.; Zadorozhnii, V.; Leonhardt, A.; Ritschel, M.; Schönhense, G.; Reznikov, Y. Soft Matter 2011, 7, 644–649. doi:10.1039/c0sm00131g 
39.  Zakhlevnykh, A. N.; Petrov, D. A. J. Mol. Liq. 2014, 198, 223–233. doi:10.1016/j.molliq.2014.06.028 
1.  Lagerwall, J. P. F.; Scalia, G., Eds. Soft Condensed Matter; World Scientific Publishing Co Pte Ltd: Singapore, 2017; Vol. 7. doi:10.1142/9789814619264_fmatter 
6.  Lagerwall, J. P. F.; Scalia, G. J. Mater. Chem. 2008, 18, 2890–2898. doi:10.1039/b802707b 
7.  Yadav, S. P.; Singh, S. Prog. Mater. Sci. 2016, 80, 38–76. doi:10.1016/j.pmatsci.2015.12.002 
10.  Zaric, S.; Ostojic, G. N.; Kono, J.; Shaver, J.; Moore, V. C.; Hauge, R. H.; Smalley, R. E.; Wei, X. Nano Lett. 2004, 4, 2219–2221. doi:10.1021/nl0486012 
33.  Byszewski, P.; Baran, M. Europhys. Lett. 1995, 31, 363–366. doi:10.1209/02955075/31/7/004 
34.  Tsui, F.; Jin, L.; Zhou, O. Appl. Phys. Lett. 2000, 76, 1452–1454. doi:10.1063/1.126061 
35.  Searles, T. A.; Imanaka, Y.; Takamasu, T.; Ajiki, H.; Fagan, J. A.; Hobbie, E. K.; Kono, J. Phys. Rev. Lett. 2010, 105, 017403. doi:10.1103/PhysRevLett.105.017403 
45.  Zakhlevnykh, A. N.; Lubnin, M. S.; Petrov, D. A. J. Exp. Theor. Phys. 2016, 123, 908–917. doi:10.1134/S1063776116100101 
46.  Zakhlevnykh, A. N.; Lubnin, M. S.; Petrov, D. A. J. Magn. Magn. Mater. 2017, 431, 62–65. doi:10.1016/j.jmmm.2016.09.044 
4.  Garbovskiy, Y. A.; Glushchenko, A. V. Solid State Phys. 2010, 62, 1–74. doi:10.1016/B9780123742933.000018 
5.  Mertelj, A.; Lisjak, D. Liq. Cryst. Rev. 2017, 5, 1–33. doi:10.1080/21680396.2017.1304835 
36.  Blinov, L. M.; Chigrinov, V. G. Electrooptic Effects in Liquid Crystal Materials; Springer: New York, NY, U.S.A., 1994. doi:10.1007/9781461226925 
3.  Brochard, F.; de Gennes, P. G. J. Phys. (Paris) 1970, 31, 691–708. doi:10.1051/jphys:01970003107069100 
30.  Matsuyama, A.; Kato, T. Phys. Rev. E 1999, 59, 763–770. doi:10.1103/PhysRevE.59.763 
31.  Matsuyama, A. J. Chem. Phys. 2010, 132, 214902. doi:10.1063/1.3447892 
39.  Zakhlevnykh, A. N.; Petrov, D. A. J. Mol. Liq. 2014, 198, 223–233. doi:10.1016/j.molliq.2014.06.028 
42.  Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2015, 393, 517–525. doi:10.1016/j.jmmm.2015.06.021 
43.  Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2016, 401, 188–195. doi:10.1016/j.jmmm.2015.10.025 
2.  Lagerwall, J. P. F.; Scalia, G. Curr. Appl. Phys. 2012, 12, 1387–1412. doi:10.1016/j.cap.2012.03.019 
32.  de Gennes, P. G.; Prost, J. The Physics of Liquid Crystals, 2nd ed.; Oxford Scientific: London, United Kingdom, 1993. 
44.  Zakhlevnykh, A. N.; Petrov, D. A. J. Exp. Theor. Phys. 2016, 123, 687–698. doi:10.1134/S1063776116090168 
18.  Dierking, I.; San, S. E. Appl. Phys. Lett. 2005, 87, 233507. doi:10.1063/1.2140069 
19.  Tomašovičová, N.; Timko, M.; Mitróová, Z.; Koneracká, M.; Rajňak, M.; Éber, N.; TóthKatona, T.; Chaud, X.; Jadzyn, J.; Kopčanský, P. Phys. Rev. E 2013, 87, 014501. doi:10.1103/PhysRevE.87.014501 
19.  Tomašovičová, N.; Timko, M.; Mitróová, Z.; Koneracká, M.; Rajňak, M.; Éber, N.; TóthKatona, T.; Chaud, X.; Jadzyn, J.; Kopčanský, P. Phys. Rev. E 2013, 87, 014501. doi:10.1103/PhysRevE.87.014501 
23.  Buluy, O.; Nepijko, S.; Reshetnyak, V.; Ouskova, E.; Zadorozhnii, V.; Leonhardt, A.; Ritschel, M.; Schönhense, G.; Reznikov, Y. Soft Matter 2011, 7, 644–649. doi:10.1039/c0sm00131g 
24.  Mitróová, Z.; Tomašovičová, N.; Timko, M.; Koneracká, M.; Kováč, J.; Jadzyn, J.; Vávra, I.; Éber, N.; TóthKatona, T.; Beaugnon, E.; Chaud, X.; Kopčanský, P. New J. Chem. 2011, 35, 1260–1264. doi:10.1039/c1nj20017h 
25.  Yoo, H. J.; Lee, S. Y.; You, N.H.; Lee, D. S.; Yeo, H.; Choi, Y. M.; Goh, M.; Park, J.; Akagi, K.; Cho, J. W. Synth. Met. 2013, 181, 10–17. doi:10.1016/j.synthmet.2013.07.023 
39.  Zakhlevnykh, A. N.; Petrov, D. A. J. Mol. Liq. 2014, 198, 223–233. doi:10.1016/j.molliq.2014.06.028 
26.  van der Schoot, P.; PopaNita, V.; Kralj, S. J. Phys. Chem. B 2008, 112, 4512–4518. doi:10.1021/jp712173n 
27.  PopaNita, V.; Kralj, S. J. Chem. Phys. 2010, 132, 024902. doi:10.1063/1.3291078 
28.  PopaNita, V. J. Chem. Phys. 2014, 140, 164905. doi:10.1063/1.4872241 
29.  PopaNita, V. J. Chem. Phys. 2015, 143, 094901. doi:10.1063/1.4928655 
43.  Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2016, 401, 188–195. doi:10.1016/j.jmmm.2015.10.025 
12.  Lynch, M. D.; Patrick, D. L. Nano Lett. 2002, 2, 1197–1201. doi:10.1021/nl025694j 
13.  Dierking, I.; Scalia, G.; Morales, P.; LeClere, D. Adv. Mater. 2004, 16, 865–869. doi:10.1002/adma.200306196 
14.  Dierking, I.; Scalia, G.; Morales, P. J. Appl. Phys. 2005, 97, 044309. doi:10.1063/1.1850606 
15.  Jeong, H. S.; Youn, S. C.; Kim, Y. H.; Jung, H.T. Phys. Chem. Chem. Phys. 2013, 15, 9493–9497. doi:10.1039/c3cp00052d 
16.  Scalia, G. ChemPhysChem 2010, 11, 333–340. doi:10.1002/cphc.200900747 
39.  Zakhlevnykh, A. N.; Petrov, D. A. J. Mol. Liq. 2014, 198, 223–233. doi:10.1016/j.molliq.2014.06.028 
41.  Zakhlevnykh, A. N. J. Magn. Magn. Mater. 2004, 269, 238–244. doi:10.1016/S03048853(03)006127 
42.  Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2015, 393, 517–525. doi:10.1016/j.jmmm.2015.06.021 
43.  Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2016, 401, 188–195. doi:10.1016/j.jmmm.2015.10.025 
44.  Zakhlevnykh, A. N.; Petrov, D. A. J. Exp. Theor. Phys. 2016, 123, 687–698. doi:10.1134/S1063776116090168 
8.  Ramirez, A. P.; Haddon, R. C.; Zhou, O.; Fleming, R. M.; Zhang, J.; McClure, S. M.; Smalley, R. E. Science 1994, 265, 84–86. doi:10.1126/science.265.5168.84 
9.  Fujiwara, M.; Kawakami, K.; Tanimoto, Y. Mol. Phys. 2002, 100, 1085–1088. doi:10.1080/00268970110109501 
10.  Zaric, S.; Ostojic, G. N.; Kono, J.; Shaver, J.; Moore, V. C.; Hauge, R. H.; Smalley, R. E.; Wei, X. Nano Lett. 2004, 4, 2219–2221. doi:10.1021/nl0486012 
11.  Cîrtoaje, C.; Petrescu, E. Physica E 2016, 84, 244–248. doi:10.1016/j.physe.2016.06.011 
11.  Cîrtoaje, C.; Petrescu, E. Physica E 2016, 84, 244–248. doi:10.1016/j.physe.2016.06.011 
14.  Dierking, I.; Scalia, G.; Morales, P. J. Appl. Phys. 2005, 97, 044309. doi:10.1063/1.1850606 
20.  Chen, H.Y.; Lee, W. Appl. Phys. Lett. 2006, 88, 222105. doi:10.1063/1.2208373 
21.  Chen, H.Y.; Lee, W. Opt. Rev. 2005, 12, 223–225. doi:10.1007/s1004300502237 
22.  Cîrtoaje, C.; Stoian, V.; Petrescu, E.; Moţoc, C. Proc. SPIE 2015, 9517, 95171G. doi:10.1117/12.2178333 
42.  Zakhlevnykh, A. N.; Petrov, D. A. J. Magn. Magn. Mater. 2015, 393, 517–525. doi:10.1016/j.jmmm.2015.06.021 
39.  Zakhlevnykh, A. N.; Petrov, D. A. J. Mol. Liq. 2014, 198, 223–233. doi:10.1016/j.molliq.2014.06.028 
37.  Burylov, S. V.; Raikher, Y. L. Mol. Cryst. Liq. Cryst. 1995, 258, 123–141. doi:10.1080/10587259508034553 
14.  Dierking, I.; Scalia, G.; Morales, P. J. Appl. Phys. 2005, 97, 044309. doi:10.1063/1.1850606 
18.  Dierking, I.; San, S. E. Appl. Phys. Lett. 2005, 87, 233507. doi:10.1063/1.2140069 
38.  Agha, H.; Galerne, Y. Phys. Rev. E 2016, 93, 042703. doi:10.1103/PhysRevE.93.042703 
3.  Brochard, F.; de Gennes, P. G. J. Phys. (Paris) 1970, 31, 691–708. doi:10.1051/jphys:01970003107069100 
41.  Zakhlevnykh, A. N. J. Magn. Magn. Mater. 2004, 269, 238–244. doi:10.1016/S03048853(03)006127 
24.  Mitróová, Z.; Tomašovičová, N.; Timko, M.; Koneracká, M.; Kováč, J.; Jadzyn, J.; Vávra, I.; Éber, N.; TóthKatona, T.; Beaugnon, E.; Chaud, X.; Kopčanský, P. New J. Chem. 2011, 35, 1260–1264. doi:10.1039/c1nj20017h 
7.  Yadav, S. P.; Singh, S. Prog. Mater. Sci. 2016, 80, 38–76. doi:10.1016/j.pmatsci.2015.12.002 
7.  Yadav, S. P.; Singh, S. Prog. Mater. Sci. 2016, 80, 38–76. doi:10.1016/j.pmatsci.2015.12.002 
8.  Ramirez, A. P.; Haddon, R. C.; Zhou, O.; Fleming, R. M.; Zhang, J.; McClure, S. M.; Smalley, R. E. Science 1994, 265, 84–86. doi:10.1126/science.265.5168.84 
9.  Fujiwara, M.; Kawakami, K.; Tanimoto, Y. Mol. Phys. 2002, 100, 1085–1088. doi:10.1080/00268970110109501 
10.  Zaric, S.; Ostojic, G. N.; Kono, J.; Shaver, J.; Moore, V. C.; Hauge, R. H.; Smalley, R. E.; Wei, X. Nano Lett. 2004, 4, 2219–2221. doi:10.1021/nl0486012 
11.  Cîrtoaje, C.; Petrescu, E. Physica E 2016, 84, 244–248. doi:10.1016/j.physe.2016.06.011 
26.  van der Schoot, P.; PopaNita, V.; Kralj, S. J. Phys. Chem. B 2008, 112, 4512–4518. doi:10.1021/jp712173n 
3.  Brochard, F.; de Gennes, P. G. J. Phys. (Paris) 1970, 31, 691–708. doi:10.1051/jphys:01970003107069100 
40.  Zakhlevnykh, A. N.; Petrov, D. A. Tech. Phys. 2012, 57, 1208–1218. doi:10.1134/S1063784212090265 
36.  Blinov, L. M.; Chigrinov, V. G. Electrooptic Effects in Liquid Crystal Materials; Springer: New York, NY, U.S.A., 1994. doi:10.1007/9781461226925 
© 2017 Petrov et al.; licensee BeilsteinInstitut.
This is an Open Access article under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (http://www.beilsteinjournals.org/bjnano)