Dynamic behavior of nematic liquid crystal mixtures with quantum dots in electric fields

  1. ,
  2. and
University Politehnica of Bucharest, Department of Physics, Splaiul Independenţei 313, 060042, Bucharest, Romania
  1. Corresponding author email
Guest Editor: A. Emelyanenko
Beilstein J. Nanotechnol. 2018, 9, 399–406. https://doi.org/10.3762/bjnano.9.39
Received 26 Jul 2017, Accepted 04 Jan 2018, Published 01 Feb 2018
Full Research Paper
cc by logo


The dynamic behavior of a mixture consisting of liquid crystalline 4-cyano-4’-pentylbiphenyl (5CB) and CdSe/ZnS quantum dots in electric fields was theoretically studied. The model was based on elastic continuum theory considering the interaction of the nematic molecules with the surrounding molecules, with the quantum dots and with the electric field. Experimental data obtained by dynamic measurements on a sample containing 0.89% (mass fraction) of CdSe/ZnS quantum dots revealed a decrease of the relaxation time compared to pure 5CB.


The expansion of liquid crystal (LC)-based devices in common life domains as well as science and engineering, requires improved technologies and new materials. Materials science provides great opportunities to these technologies by synthesizing new nanoparticles that can be mixed with liquid crystals: carbon nanotubes [1-3], graphene, magnetic nanoparticles [4,5], gold nanoparticles [6], quantum dots (QDs) [7-12] or other nanomaterials [13-17] that can be effectively used in electro-optical devices.

When added to liquid crystals, quantum dots may seriously influence their behavior under an electric field due to the anchoring forces acting on the quantum dot surface. Dynamic experiments performed in alternating electric fields proved that by adding a small amount of CdSe/ZnS quantum dots in thermotropic nematic liquid crystal with positive dielectric anisotropy, we obtain a decrease of the relaxation time.

When an external electric field higher than the Fréedericksz transition threshold [18,19] is applied transversely to a planar liquid crystal cell, the molecules change their orientation tending to align their director in parallel to the field. This reorientation is illustrated by an intensity variation of a laser beam crossing through the sample. The time period after which a constant intensity is obtained, which means that an equilibrium state is reached, is characterized by the relaxation time and can be experimentally evaluated from plotting the intensity as a function of the time. A theoretical analysis of the LC + QDs system based on elastic continuum theory allowed us to express the time evolution of the molecular distortion angle inside the cell and to calculate the relaxation time. Theoretical values and experimental data were in good agreement showing that the model offers an important tool for engineers to precisely control LC-based device parameters.

Theoretical Background

When a nematic liquid crystal with positive dielectric anisotropy is exposed to an electric field, and if the intensity of the field is higher than the Fréedericksz transition threshold, the molecular director changes its orientation, trying to align in parallel with the field direction. During this reorientation process, the refractive index of the sample also changes:


where no and ne are, respectively, the ordinary and extraordinary refractive indexes and θ is the angle between the direction of the light and ne. Between the extraordinary and ordinary rays, there is a path difference which can be calculated as a function of the cell thickness, L:


with the corresponding phase difference Δφ = (2πδ)/λ.

For a planar cell configuration we have:


where θm is the maximum deviation angle (in the middle in the cell). If small deviation angles are considered, the phase difference becomes:


and the intensity of the laser beam crossing through the sample is:


The equation of the maximum deviation angle as a function of time is obtained from the total free energy of the system by applying the Euler–Lagrange equation. Elastic continuum theory states that a liquid crystal is acting like a continuous fluid and the interaction forces between its molecules are elastic. Taking into account a strong anchoring on the glass support, the free energy density of such a system with added quantum dots is:


where fN is the liquid crystal free energy density, fE is a term combining the electric field action on the liquid crystal molecules in the presence of quantum dots and fINT reflects the interaction between quantum dots and liquid crystal molecules.

For a planar aligned cell (Figure 1) with thickness L, and when the voltage is just above the Fréedericksz transition threshold, the deviation angle is very small and the free energy caused by the elastic forces is:


where K1 and K3 are the splay and bent elastic constants and θz = ∂θ/∂z.


Figure 1: Molecular director orientation of a liquid crystal with a positive anisotropy in a planar cell exposed to an electric field.

Due to the presence of the QDs, the electric properties of the mixture change. Hence, instead of perpendicular dielectric permittivity, [Graphic 1], and electric anisotropy, εa, we use the effective perpendicular permittivity, [Graphic 2], and effective anisotropy, εaeff.


If a small spherical cavity is considered inside the liquid crystal, a polarisation field on the hole surface appears. According to classical theory developed by Landau and Lifshitz [20], the polarization field around a spherical cavity inside a dielectric medium is proportional to the applied field, so:


Inside an anisotropic dielectric (the liquid crystal) the spherical cavity changes into an ellipsoid with the semi-axes [Graphic 3] and [Graphic 4], where b is the cavity radius, [Graphic 5] is the parallel permittivity, [Graphic 6] is the perpendicular permittivity. In this case H becomes:


where d(z) the depolarization coefficient, which for an ellipsoid with almost spherical shape (horizontal axis equal b = R and vertical axis a = R + 0.0005R) is

[Graphic 7]

If a quantum dot is placed inside this cavity, another polarization field will occur but because of its semiconductor nature, this field is much smaller than the one produced at the surface of the dielectric and might be neglected. Thus, the effective field acting on the nematic molecules is smaller than the applied field E by EQD, i.e., E′ = EHE.

The free energy density may also be written as:


where D is the electric induction. Thus the free energy is:


By taking into account that the applied voltage on the sample is

[Graphic 8]

and considering small angle approximation for the configuration presented in Figure 2 we get:


Figure 2: Molecular director and quantum dot orientation in a planar LC cell subjected to electric field.

Then the free energy density is:


The first term does not depend on the molecular orientation so we may neglect it because it disappears when applying the Euler–Lagrange equation. Thus, we may write:






is the effective dielectric anisotropy of the liquid crystal containing quantum dots as a function of the dielectric anisotropy of the LC, (εa), of the dielectric constant of the QDs (εQD) and the volumetric fraction of the quantum dots (cv).

The intraction free energy between nematic molecules and a quantum dot can be calculated from the Rapini formula [21]. Considering an elliptical quantum dot with one of the axes just a little bit longer than the others and the configuration given in Figure 3, as well as homeotropic anchoring of liquid crystal molecules on the QD surface [12], the interaction energy for a particle can be obtained using a similar procedure as presented by Burylov [22]:


where w is the anchoring energy on the surface of the nanoparticle, [Graphic 9] is the nematic director on the QD surface and [Graphic 10] is the direction of molecular director far away from the QD surface.


Figure 3: Nematic molecular orientation around an ellipsoidal QD.

According to the manufacturers, the quantum dots are spherical particle with the radius R = 3 nm (6 nm in diameter). Considering their small anisotropy we approximate the sphere with an almost spherical ellipsoid having a short axis (b = R) and a long axis close to b. Thus, we obtain:


This equation can be solved and written in a simplified form as:


where fpct is a constant depending on the anchoring energy and the particle geometry and disappears when applying the Euler–Lagrange equations. Thus, it may be neglected and fp0 = 11.52wb2π . Considering the volumetric concentration of quantum dots, we obtain:


where f0 = (11.52/2)wπR2 is a constant depending on the geometrical parameters of the system and the LC + QD interaction energy calculated for the ellipsoid QD presented in Figure 3. When a dynamic study of the mixture is made, the total free energy density must also take into account the dissipative time-dependent term, and the total free energy becomes:


By applying the Euler–Lagrange equation and assuming that θ is small we obtain:


where θzz = ∂2θ/∂z2. Considering the deviation angle from Equation 3, we obtain the dependency of the deviation angle on the relaxation time:



[Graphic 11]


[Graphic 12]

Solving this equation leads to


where θm(∞) and θm(0) are the limit values of the deviation angle in the center of the cell and τon is:


When the voltage is switched off, we follow the same procedure after replacing f0E = 0 in Equation 22. Thus, we obtain:


where θm(0) is the maximum deviation angle at t = 0, ξ is a constant depending on LC properties and τoff is:



For the sample preparation we used CdSe/ZnS quantum dots in toluene with an average size of 6 nm and a fluorescence wavelength of 630 nm from Aldrich. The QD suspension was left in a bottle in clean vacuum environment until the toluene was completely evaporated. Then, 4-cyano-4’-pentylbiphenyl (5CB) was added and sonicated for three days, several hours each day until a homogenous pink mixture was obtained. The homogeneity was checked by using polarization microscopy. The mass concentration was 0.89%, which led to a volumetric fraction of 1.16% (the QD density is 0.87 g/cm3 and the density of 5CB is 1.02 g/cm3). The mixture was used to fill a 15 μm thick planar aligned cell from Instek, which was mounted in a holder placed in a Mettler Toledo hot stage. A reference sample was made using pure 5CB. The terminals of the holder were connected to a power source that produced a square wave signal with a frequency of 10 kHz. For voltages higher than the Fréedericksz transition threshold, the nematic director changes its orientation trying to align with the field. A He–Ne laser beam crossed the sample through the windows of the holder and the emergent signal was recorded with a Thor Lab photovoltaic cell. Two crossed polarizers at 45° were placed on both sides of the sample to obtain equal intensities for ordinary and extraordinary rays. An acquisition system connected to a computer allowed data recording every 20 ms. A schematic representation of the experimental setup is presented in Figure 4.


Figure 4: Experimental setup for the dynamic study of the LC + QDs mixture.

Results and Discussion

Fréedericksz transition threshold and relaxation time were determined for both samples at a constant temperature of 28 °C. The Fréedericksz transition threshold was determined by slowly increasing the applied voltage and while recording the laser intensity. Preliminary measurements showed no changes of the laser intensity below 0.5 V. Hence, the voltage steps were 0.1 V below this value and 0.01 V above it. The data were represented in an intensity versus voltage plot (Figure 5). Above the threshold value, the intensity of the emergent laser beam exhibits a series of maxima and minima.


Figure 5: Fréedericksz transition for 5CB and 5CB + QD. Lines are guide to the eyes.

For the dynamic measurements the value of the first intensity maximum was chosen. Thus, the applied voltage was without any doubt above the Fréedericksz transition and the deviation angle was small as considered in the theoretical model. The intensity versus time plot was recorded for each sample and the function in Equation 5 was used to fit the data. A power switch was connected to the acquisition system so that when the field was switched on the acquisition started at the same time to measure the intensity when the field is switched on. A similar procedure was used when the power source was disconnected and at the same time the acquisition started to measure intensity versus time plot when the field was switched off.

In Figure 6 we presented the light intensity versus time plot for the pure LC (Figure 6a,b) and for the mixture of LC with QDs (Figure 6c,d). Experimental data were fitted using Equation 5 with the time dependency of thetam given by Equation 25 and Equation 27, where L = 15 μm, K1 = 6.2 × 10−12 N, R = 3 nm and cv = 1.16% (volumetric fraction). Therefore, we obtain for the sample containing QDs, the experimental relaxation time τon−exp = 0.167 s when the field is switched on and τoff−exp = 0.139 s when the field is switched off.


Figure 6: Dynamic measurement of emergent light intensity: a) for LC sample when the field is on; b)for LC sample when the field is off; c) for LC + QDs sample when the field is on; for LC + QDs sample when the field is off. The dots represent experimental values and the red lines are theoretical representations.

By fitting the data of the 5CB sample, we obtained a relaxation time τon−exp0 = 0.172 s when the field is switched on and τoff−exp0 = 0.157 s when the field is switched off. The relaxation time decreases when quantum dots are added. Generally, this result was reported in similar researches [9-11] but only when the field was switched on. When the field was off the results in [9] showed an increase of the relaxation time for the sample containing quantum dots. The measurements were made on similar systems both of which emplyed CdSe/ZnS in 5CB but in our case a temperature of 28 °C yielded a lower viscosity and, comsequently, a lower relaxation time especially when the field is switched off. In [8] the authors demonstrate an increase of the elastic coefficient K1 with the QD concentration, which is in good agreement with Equation 26 and Equation 28 where the effective elastic constant is proportional with the volumetric fraction of the added quantum dots.

Basu and co-workers [7] proposed a different model of a quantum dot dispersion in a nematic liquid crystal. According to this model, dots are gathering together forming long chains around which liquid crystal molecules are homogeneously aligned. Such an alignment could not explain the decrease of the relaxation time when the field is switched on because the anchoring forces on quantum dots surface will slow down the molecular movement.

A possible explanation for this behavior might be given considering a homeotropic alignment of the LC molecules on the QD surface [12]. Thus, if the applied field is high enough to exceed the Fréedericksz transition threshold some of the molecules are already parallel to the field and only those placed on the horizontal sphere circumference are parallel to the cell plates, so the response time is shorter in the quantum dots sample (Figure 7a). When the field is switched off, all the molecules are distorted and their returning to the original position is “helped” by the anchoring forces on the QD surface and a slight decrease of the relaxation time is observed (Figure 7b).


Figure 7: Schematic representation of the molecular orientation around a QD a) when the field is switched on b) when the field is switched off.


We observed an improvement in the turn-off time of the nematic liquid crystal 5CB containing CdSe/ZnS quantum dots of a nominal size of 6 nm diameter. The results are potentially useful in other cell geometries and with other materials. The theoretical model might be extended to other geometries and used to theoretically optimize the performance of QD–liquid crystal composites. The calculations presented in this paper are based on particular experimental results obtained at 28 °C, for a relatively high concentration of quantum dots. Further measurements are to be done in the future, over a wide temperature range and under consideration of other parameters, to improve this model and obtain an accurate theoretical description of these systems.


The work has been funded by University Politehnica of Bucharest, through the ’Excellence Research Grants’ Program, UPB–GEX. Identifier: UPB–EXCELENŢĂ–2016 Research project title Phase modulation in nematic liquid crystals subjected to magnetic field. Contract number 24/26.09.2016.


  1. Cirtoaje, C.; Petrescu, E. Physica E 2016, 84, 244–248. doi:10.1016/j.physe.2016.06.011
    Return to citation in text: [1]
  2. Samoilov, A. N.; Minenko, S. S.; Fedoryako, A. P.; Lisetski, L. N.; Lebovka, N. I.; Soskin, M. S. Funct. Mater. 2014, 21, 190–194. doi:10.15407/fm21.02.190
    Return to citation in text: [1]
  3. Bale, S.; Liyana-Arachchi, T. P.; Hung, F. R. Mol. Simul. 2016, 42, 1242–1248. doi:10.1080/08927022.2016.1174859
    Return to citation in text: [1]
  4. Cirtoaje, C.; Petrescu, E.; Stan, C.; Creanga, D. Physica E 2016, 79, 38–43. doi:10.1016/j.physe.2015.12.006
    Return to citation in text: [1]
  5. Mănăilă-Maximean, D.; Cîrtoaje, C.; Dănilă, O.; Donescu, D. J. Magn. Magn. Mater. 2017, 438, 132–137. doi:10.1016/j.jmmm.2017.02.034
    Return to citation in text: [1]
  6. Lewandowski, W.; Fruhnert, M.; Mieczkowski, J.; Rockstuhl, C.; Górecka, E. Nat. Commun. 2015, 6, 6590. doi:10.1038/ncomms7590
    Return to citation in text: [1]
  7. Basu, R.; Iannacchione, G. S. Phys. Rev. E 2009, 80, 010701. doi:10.1103/PhysRevE.80.010701
    Return to citation in text: [1] [2]
  8. Vakulin, D. A.; Frenkel, D. A.; Gavrish, E. O.; Konshina, E. A. Mol. Cryst. Liq. Cryst. 2015, 612, 110–116. doi:10.1080/15421406.2015.1030579
    Return to citation in text: [1] [2]
  9. Konshina, E. A.; Galin, I. F.; Shcherbinin, D. P.; Gavrish, E. O. Liq. Cryst. 2014, 41, 1229–1234. doi:10.1080/02678292.2014.912764
    Return to citation in text: [1] [2] [3]
  10. Gavrish, E. O.; Galin, I. F.; Konshina, E. A.; Vakulin, D. A. Mol. Cryst. Liq. Cryst. 2015, 615, 50–56. doi:10.1080/15421406.2015.1066954
    Return to citation in text: [1] [2]
  11. Galin, I. F.; Konshina, E. A. Mol. Cryst. Liq. Cryst. 2015, 615, 57–62. doi:10.1080/15421406.2015.1066955
    Return to citation in text: [1] [2]
  12. Prodanov, M. F.; Pogorelova, N. V.; Kryshtal, A. P.; Klymchenko, A. S.; Mely, Y.; Semynozhenko, V. P.; Krivoshey, A. I.; Reznikov, Y. A.; Yarmolenko, S. N.; Goodby, J. W.; Vashchenko, V. V. Langmuir 2013, 29, 9301–9309. doi:10.1021/la401475b
    Return to citation in text: [1] [2] [3]
  13. Cîrtoaje, C.; Petrescu, E.; Stoian, V. Physica E 2015, 67, 23–27. doi:10.1016/j.physe.2014.11.004
    Return to citation in text: [1]
  14. Motoc, C.; Cirtoaje, C.; Stoica, A.; Stoian, V.; Albu, A. M. UPB Sci. Bull., Ser. A 2013, 75, 181–186.
    Return to citation in text: [1]
  15. Zakhlevnykh, A.N.; Lubnin, M. S.; Petrov, D. A. J. Magn. Magn. Mater. 2008, 431, 62–65. doi:10.1016/j.jmmm.2016.09.044
    Return to citation in text: [1]
  16. Ganea, C.; Manaila-Maximean, D. UPB Sci. Bull., Ser. A 2011, 73, 209–216.
    Return to citation in text: [1]
  17. Manaila-Maximean, D.; Rosu, C.; Yamamoto, T.; Yokoyama, H. Mol. Cryst. Liq. Cryst. 2004, 417, 215–226. doi:10.1080/15421400490478867
    Return to citation in text: [1]
  18. Fréedericksz, V.; Repiewa, A. Z. Phys. 1927, 42, 532–546. doi:10.1007/BF01397711
    Return to citation in text: [1]
  19. Fréedericksz, V.; Zolina, V. Trans. Faraday Soc. 1933, 29, 919–930. doi:10.1039/TF9332900919
    Return to citation in text: [1]
  20. Landau, L. D.; Lifshitz, E. M. Electrodynamics of Continuous Media; Pergamon Press: Oxford, United Kingdom, 1960.
    Return to citation in text: [1]
  21. Rapini, A.; Papoular, M. J. Phys. Colloques 1969, 30, 54–56. doi:10.1051/jphyscol:1969413
    Return to citation in text: [1]
  22. Burylov, S. V.; Raikher, Yu. L. J. Magn. Magn. Mater. 1993, 122, 62–65. doi:10.1016/0304-8853(93)91040-E
    Return to citation in text: [1]
Other Beilstein-Institut Open Science Activities