Journal of Modern Physics
Vol.05 No.18(2014), Article ID:52632,8 pages
10.4236/jmp.2014.518204

Dynamics of a ±1/2 Defect Pair in a Confined Geometry

Lixia Lu1,2,3, Zhidong Zhang3*

1Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun, China

2University of Chinese Academy of Sciences, Beijing, China

3School of Science, Hebei University of Technology, Tianjin, China

Email: *zhidong_zhang1961@163.com

Copyright © 2014 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 4 October 2014; revised 1 November 2014; accepted 23 November 2014

ABSTRACT

This paper investigated the dynamics of a dipole of parallel wedge disclination lines in a confined geometry, based on Landau-de Gennes theory. The behavior of the pair depends on the competition between two kinds of forces: the attractive force between the two defects, aggravating the annihilation process, and the anchoring forces coming from the substrates, inhibiting the anni- hilation process. There are three states when the system is equilibrium, divided by two critical thicknesses and (existing when, is the initial distance between the two defects), both changing linearly with. When the cell gap, the two defects coalesce and annihilate. The dynamics follows the function of during the annihilation step when is sufficiently large, relative to, where is the relative distance between the pair and is the coalescence time. decreases with the decrease of or the increase of. The anni- hilation process has delicate structures: when and or and, the two defects annihilate and the system is uniaxial at equilibrium state; when and, the two defects coalesce and annihilate, and the system is not uniaxial, but biaxial in the region where the defects collide. When, the defects can be stable existence.

Keywords:

±1/2 Defect Pair, Dynamics, Stable Existence, Confined Geometry

1. Introduction

Defects are ubiquitous in nature and are important in particle physics, cosmology, and condensed matter physics [1] . This is explained by the importance of the role of defects in the course of different processes (phase transitions, plastic deformations, electronic processes, etc.) [2] . Defects in liquid crystals (LCs) affect the manifestation of a number of optical, field, hydrodynamic, and other effects [2] . They have posed important problems in optoelectric applications [3] as well as in fundamental physics. They can trap nanoparticles [4] [5] , mediate characteristic interaction between colloidal particles [5] - [7] and provide ordered templates for colloidal microassembly [8] - [11] . Therefore, location control of topological defects is important for defect-mediated colloidal assembly, and guided polymerization and crystallization in the self-organized nematic order [5] .

Defects in LCs are topological defects [2] , also called disclinations [12] . They appear spontaneously at the isotropic-nematic transition when the symmetry of the isotropic phase is broken to the symmetry of the nematic phase [13] [14] . They are moving during the equilibrium process. Therefore, the dynamics of topological defects in the ordering process has attracted the interest of many researchers and has been largely studied theoretically, numerically and experimentally over the last decades [13] - [26] . All of these studies focus on the relationship between and during the annihilation process, where is the relative distance between the defect pair and is the time for the annihilation. Most previous studies have been confined to investigating the annihilation of an isolated defect pair [15] - [26] . However, in real systems, disclinations are never isolated, but subjected to the anchoring forces coming from the substrates of the cell containing the liquid crystal. As a consequence, in a confined geometry, the substrate anchoring is expected to strongly influence the interactions between defect lines. A similar anchoring effect has been invoked to explain the annihilation dynamics of nematic point defects confined in capillary tubes [27] and in hybrid cells [13] [14] . However, the behavior of defect lines in a confined geometry remains totally unexplored.

In this paper, based on our previous study of the relaxation dynamics of a dipole of disclination lines with in a thin hybrid aligned nematic (HAN) cell [28] , we continue to study the relaxation dynamics of the dipole in a confined geometry.

2. Theoretical Basis

2.1. Basic Equations

The theoretical argument is based on Landau-de Gennes theory, in which the orientational order is described by a second-rank traceless and symmetric tensor [29]

, (1)

where and are the eigenvectors and the corresponding eigenvalues of. In the isotropic phase, vanishes. When two eigenvalues of coincide, the liquid crystal is in a uniaxial state, and can be recast in the form

, (2)

where is the uniaxial scalar order parameter, and the unit vector is the nematic director. When all eigen values of are distinct, the liquid crystal is in a biaxial state. The degree of biaxiality of can be defined as [29] - [31]

, (3)

is a convenient parameter for illustrating spatial inhomogeneities of and ranges in the interval. In all uniaxial states, , while states with maximal biaxiality correspond to.

In the reduced space defined by Schopohl and Sluckin [32] , the dynamics of can be described as

. (4)

is the free-energy density in the reduced space. is given by, and. Here, is a nematic rotational viscosity. For more details, please refer to our previous study [28] [33] [34] or see Appendix.

2.2. Numerical Methods

Time evolution of in the Equation (4) is computed using a two-dimensional finite-difference-iterative- method employed in our previous studies [28] [33] [34] . The local values of the scalar order parameter, , and the director, , can be obtained from through its largest eigenvalue and its associated eigenvector, respectively.

The reduced space is discretized into grids with the same interval of. The discretization of time steps given by is sufficient to guarantee the stability of the numerical procedure. According to the parameters of 5CB given in [35] , we have J/m3, J/m3, J/m3, J/m, Pa∙s and [32] . Then nm,. The bulk nematic-isotropic transition occurs at. The scaled temperature is set at, which guarantees the system being in the nematic state.

Consider the pair positioned at in a cell with the cell gap, which guarantees the substrates do not affect the pair. The lengths, and, of the cell along the x- and the y-axes are much larger than. The system is relaxed followed the Equation (4) from the initial condition:

,

where, [12] . is the distance between the origin

of coordinate and the observation point, is the angle between the singularity―observation point line and the x-axis, is the angle between the singularity―observation point line and the x-axis, as shown in Figure 1. is the angle between and the x-axis. When the distance between the defects is relaxed to (any value wanted), liquid crystal cells with different cell gap are obtained using the plan-parallel technique. Take them as the initial conditions of our numerical calculations. The strong anchoring conditions on the bounding plates and free boundaries in the x-direction are used.

3. Results and Discussion

In order to analyze the influence of on the behavior of the defect pair, the time dependence of the distance between the defect pair for different with and moving velocity of the defects as a function of time are shown in Figure 2.

The squares, circles and triangles represent numerical results with different cell gap, respectively, and the full lines are the corresponding fitting curves with function of [13] [21] in Figure 2(a), where is the coalescence time. decreases and the fitting error increases with the decrease of, respectively. decreases from to and the relative error increases from 3.9% to 4.3% when decreases from to. Further decrease to, decreases to, the relative error is up to 13.4%, which means that the time dependence of has deviated from the exponential law, as shown in Figure 2(a).

Figure 1. The geometry of the problem.

Figure 2. (a) The time dependence of the distance r between the defect pair in a cell for different d with; (b) Moving velocity of the defects as a function of time for different d with.

The distance between two isolated oppositely charged defects is described as [13] [22] . In a confined geometry, the substrate anchoring forces influences the interactions between the defects. The thinner the cell gap, the greater the impact of the anchoring forces. Consequently, decreases and the fitting error increases with the decrease of d. When decreases to, the anchoring forces are so strong that they make the time dependence of deviate from the exponential law.

The moving velocity of the defects increases rapidly with time because the elastic force between the two defects is [12] , but decreases with the decrease of. It decreases slightly when decreases from to, but rapidly when the cell gap decreases to and, as shown in Figure 2(b). This results in the increase of the coalescence time with the decrease of. The coalescence time slightly increases from to when d decreases from to, but soars to and when d de- creases to and, respectively.

It indicates that the anchoring forces coming from the substrates produce a negligible effect on the defects when. But, when, the anchoring forces strongly affect the behavior of the dipole: inhibiting the moving movement of the defects. The thinner the cell gap, the greater the impact of the anchoring forces. Consequently, the moving velocity of the defects decreases with the decrease of the cell gap.

In order to analyze the influence of on the behavior of the defect pair, the time dependence of the distance between the defect pair and moving velocity of the defects as a function of time for different with are shown in Figure 3.

decreases with the increase of r0. It decreases from to when increases from to. The time dependence of follows the exponential law when. When increases to, decreases to 0.300 ± 0.046, hence the relative error is up to 15%, which indicates the time dependence of has deviated from the exponential law. The direct elastic interaction between the defects gives an elastic force [12] , resulting in the rapid decrease of the moving velocity with the increase of r0, as shown in Figure 3(b). Therefore, the coalescence time soars from to when r0 increases from to.

Increasing the initial distance has the same effect as reducing the cell gap: inhibiting the relative movement of the defects. The behavior of the pair depends on the competition between the attractive force, determined by the distance between the defect pair, and the anchoring forces, determined by the cell gap d. The attractive force aggravates the annihilation process, while the anchoring forces inhibit the annihilation process. If the role of one kind of force is weakened, then the impact of the other will be relatively enhanced. Therefore, the impact of the anchoring forces are relatively enhanced with the increase of, which results in the decrease of with the increase of and the time dependence of deviating from the exponential law when increases to.

The anchoring forces and the attractive force must be balanced and the defects can be stable existence under certain conditions. In order to find the conditions, the dynamic behavior of the dipole is studied in detail by varying the cell gap and the initial distance. The results are shown in Figure 4 and Figure 5.

Figure 3. (a) The time dependence of the distance r between the defect pair in a cell for different r0 with d = 20ξ. Squares, circles and triangles represent numerical results. The full lines are the corresponding fitting curves with; (b) Moving velocity of the defects as a function of time for different r0 with d = 20ξ.

When, there are two equilibrium states, as shown in Figure 4(a), Figure 4(b) and Figure 5. The critical thicknesses are and, respectively. When, the two defects annihilate because the attractive force between them plays a major role and a uniform uniaxial state is formed when the system is equilibrium. When, the defects can be stable existence because the anchoring forces are so strong that they can inhibit the relative movement of the defects. The thinner the cell gap, the greater the impact of the anchoring forces. As a consequence, the relative moving distance decreases with the decrease of.

When, there are three equilibrium states, divided by two critical thicknesses: and, as shown in Figures 4(c)-(e) and Figure 5. The two critical thicknesses both change linearly with the initial distance. The corresponding fitting equations are and, respectively. When, the two defects annihilate. The annihilation process has delicate structures. When, the two defects annihilate and a uniform uniaxial state is formed when the system is equilibrium. When, the two defects annihilate. However, the system is not uniform uniaxial, but biaxial in the region where the defects collide. Since the two defects are oppositely charged and the total charge is conserved, the defects disappear when they collide. But, the director in the region where the defects collide distorts inhomogeneously, which creates the biaxial layers. The range of the cell gap corresponding to this state increases with the decrease of. When, the dipole can be stable existence as the system is equilibrium because the anchoring forces play a major role under this condition.

4. Conclusions

The behavior of the dipole in a confined geometry depends on the competition between two kinds of forces: the attractive force and the anchoring force. The former aggravates the annihilation process, while the latter inhibits it. The system has three equilibrium states, divided by two critical thicknesses and. When,

the pair annihilates, and the dynamics follows the function of as is sufficiently large. de-

creases with the decrease of or increase of. The annihilation process has delicate structures: when and or and, the system is uniaxial at equilibrium state; while, when and, the system is biaxial in the region where the defects collide. When, the dipole can be stable existence.

As far as we know, it is the first time to discover the stable existence of the oppositely charged defects under certain conditions. This research plays a major role in the formation and control of topological defects, and has significant academic value for mediation of defects on colloidal particles in nematic liquid crystals.

Acknowledgements

We thank the Editor and the referee for their comments. Research of Z. Zhang is funded by National Natural

Figure 4. The time dependence of across the defect center parallel to the x-axis at equilibrium states: (a); (b); (c); (d); (e).

Figure 5. Phase diagram. The circles and squires represent numerical results of dc1 and dc2, respectively; the full lines are the corresponding linear fits.

Science Foundation of China grant No. 11374087 and Key Subject Construction Project of Hebei Province University. This support is greatly appreciated.

References

  1. Carbone, G., Lombardo, G., Barberi, R., Musevic, I. and Tkalec, U. (2009) Physical Review Letters, 103, Article ID: 167801.
  2. Kurik, M.V. and Lavrentovich, O.D. (1988) Soviet Physics―Uspekhi, 31, 196-224. http://iopscience.iop.org/0038-5670 http://dx.doi.org/10.1070/PU1988v031n03ABEH005710
  3. Tiribocchi, A., Gonnella, G., Marenduzzo, D. and Orlandini, E. (2010) Applied Physics Letters, 97, Article ID: 143505. http://dx.doi.org/10.1063/1.3496472
  4. Takuya, O. and Jun-Ichi, F. (2012) Nature Communications, 3, 701-707. http://dx.doi.org/10.1038/ncomms1709
  5. Takuya, O., Jun-Ichi, F., Kosuke, S. and Yamaguchi, T. (2012) Physical Review E, 86, Article ID: 030701(R).
  6. Tkalec, U., Ravnik, M., Čopar, S., Žumer, S. and Muševič, I. (2011) Science, 333, 62-65. http://dx.doi.org/10.1126/science.1205705
  7. Ravnik, M., Alexander, G.P., Yeomans, J.M. and Žumer, S. (2011) Proceedings of the National Academy of Sciences of the United States of America, 108, 5188-5192. http://dx.doi.org/10.1073/pnas.1015831108
  8. Fleury, J.B., Pires, D. and Galerne, Y. (2009) Physical Review Letters, 103, Article ID: 267801. http://dx.doi.org/10.1103/PhysRevLett.103.267801
  9. Yoon, D.K., Choi, M.C., Kim, Y.H., Kim, M.W., Lavrentovich, O.D. and Jung, H.T. (2007) Nature Materials, 6, 866-870. http://dx.doi.org/10.1038/nmat2029
  10. Milette, J., Cowling, S.J., Toader, V., Lavigne, C., Saez, I.M., Lennox , R.B., Goodby, J.W. and Reven, L. (2012) Soft Materials, 8, 173-179. http://dx.doi.org/10.1039/c1sm06604h
  11. Coursault, D., Grand, J., Zappone, B., Ayeb, H., Lévi, G., Félidj, N. and Lacaze, E. (2012) Advanced Materials, 24, 1461-1465. http://dx.doi.org/10.1002/adma.201103791
  12. Gennes, P.G. and Prost, J. (2008) The Physics of Liquid Crystals. 2nd Edition, Science Press, Beijing .
  13. Minoura, K., Kimura, Y., Ito, K., Hayakawa, R. and Miura, T. (1998) Physical Review E, 58, 643-649. http://dx.doi.org/10.1103/PhysRevE.58.643
  14. Bogi, A., Martinot-Lagarde, P., Dozov, I. and Nobili, M. (2002) Physical Review Letters, 89, Article ID: 225501. http://dx.doi.org/10.1103/PhysRevLett.89.225501
  15. Chuang, I., Turok, N. and Yurke, B. (1991) Physical Review Letters, 66, 2472-2475. http://dx.doi.org/10.1103/PhysRevLett.66.2472
  16. Pargellis, A., Turok, N. and Yurke, B. (1991) Physical Review Letters, 67, 1570-1573. http://dx.doi.org/10.1103/PhysRevLett.67.1570
  17. Chuang, I., Yurke, B., Pargellis, A.N. and Turok, N. (1993) Physical Review E, 47, 3343-3356. http://dx.doi.org/10.1103/PhysRevE.47.3343
  18. Pargellis, A.N., Mendez, J., Srinivasarao, M. and Yurke, B. (1996) Physical Review E, 53, R25-R28. http://dx.doi.org/10.1103/PhysRevE.53.R25
  19. Minoura, K., Kimura, Y., Ito, K. and Hayakawa, R. (1997) Molecular Crystals and Liquid Crystals Science and Technology. Section A: Molecular Crystals and Liquid Crystals, 302, 345-355. http://dx.doi.org/10.1080/10587259708041847
  20. Wang, W., Shiwaku, T. and Hashimoto, T. (1998) The Journal of Chemical Physics, 108, 1618-1625. http://dx.doi.org/10.1063/1.475532
  21. Svensek, D. and Žumer, S. (2002) Physical Review E, 66, Article ID: 021712. http://dx.doi.org/10.1103/PhysRevE.66.021712
  22. Tóth, G., Denniston, C. and Yeomans, J.M. (2002) Physical Review Letters, 88, Article ID: 105504. http://dx.doi.org/10.1103/PhysRevLett.88.105504
  23. Bradač, Z., Kralj, S., Svetec, M. and Žumer, S. (2003) Physical Review E, 67, Article ID: 050702.
  24. Svetec, M., Kralj, S., Bradac, Z. and Žumer, S. (2006) The European Physical Journal E, 20, 71-79. http://dx.doi.org/10.1140/epje/i2005-10120-9
  25. Guimarães, R.R., Mendes, R.S., Fernandes, P.R.G. and Mukai, H. (2013) Journal of Physics: Condensed Matter, 25, Article ID: 404203.
  26. Giomi, L., Bowick, M.J., Ma, X. and Marchetti, M.C. (2013) Physical Review Letters, 110, Article ID: 228101. http://dx.doi.org/10.1103/PhysRevLett.110.228101
  27. Peroli, G.G. and Virga, E.G. (1999) Physical Review E, 59, 3027-3032. http://dx.doi.org/10.1103/PhysRevE.59.3027
  28. Lu, L.X. and Zhang, Z.D. Chinese Physics B, in Press.
  29. Bisi, F., Gartland, E.C., Rosso, R. and Virga, E.G. (2003) Physical Review E, 68, Article ID: 021707. http://dx.doi.org/10.1103/PhysRevE.68.021707
  30. Barberi, R., Ciuchi, F., Durand, G., Iovane, M., Sikharulidze, D., Sonnet, A. and Virga, E. (2004) The European Physical Journal E: Soft Matter, 13, 61-71. http://dx.doi.org/10.1140/epje/e2004-00040-5
  31. Lombardo, G., Ayeb, H. and Barberi, R. (2008) Physical Review E, 77, Article ID: 051708.
  32. Schopohl, N. and Sluckin, T.J. (1987) Physical Review Letters, 59, 2582-2584. http://dx.doi.org/10.1103/PhysRevLett.59.2582
  33. Lu, L.X., Zhang, Z.D. and Zhou, X. (2013) Acta Physica Sinica, 62, Article ID: 226101.
  34. Zhou, X. and Zhang, Z.D. (2013) International Journal of Molecular Sciences, 14, 24135-24153. http://dx.doi.org/10.3390/ijms141224135
  35. Qian, T.Z. and Sheng, P. (1997) Physical Review E, 55, 7111-7120. http://dx.doi.org/10.1103/PhysRevE.55.7111

Appendix

The free-energy density of a NLC without external field can be expressed as

. (A1)

They are, respectively, the elastic and the bulk free-energy densities. The former is induced by the inhomogeneous order in LCs. It can be given the form

, (A2)

where and are elastic constants. The latter is a potential that depends on. It is conventionally described by an expansion in up to the fourth order

. (A3)

Usually it is assumed that, where is temperature and is the supercooling temperature., and are coefficients. The traceless and symmetry of are taken into account by the introduction of the Lagrange parameter tensor [32] :. The full free-energy functional is

. (A4)

The calculation can be simplified by dimensionless variables. Here, we follow the rescaling of Schopohl and Sluckin [32] by defining the following dimensionless quantities

,

where is the order parameter at the isotropic-nematic phase transition point and

is the characteristic length for order-parameter changes. Hence, the dynamics

of can be described as

. (4)

NOTES

*Corresponding author.