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

Lixia Lu^{1,2,3}, Zhidong Zhang^{3*}

^{1}Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun, China.

^{2}University of Chinese Academy of Sciences, Beijing, China.

^{3}School of Science, Hebei University of Technology, Tianjin, China.

**DOI: **10.4236/jmp.2014.518204
PDF
HTML XML
4,690
Downloads
5,209
Views
Citations

This paper
investigated the dynamics of a dipole of ±1/2 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 annihilation process. There are
three states when the system is equilibrium, divided by two critical thicknesses* d _{c1}* and

Keywords

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

Share and Cite:

Lu, L. and Zhang, Z. (2014) Dynamics of a ±1/2 Defect Pair in a Confined Geometry. *Journal of Modern Physics*, **5**, 2080-2088. doi: 10.4236/jmp.2014.518204.

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/m^{3}, J/m^{3}, J/m^{3}, 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 r_{0}. 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 r_{0}, as shown in Figure 3(b). Therefore, the coalescence time soars from to when r_{0} 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 r_{0} 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 r_{0} 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 d_{c}_{1} and d_{c}_{2}, 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.

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.

Conflicts of Interest

The authors declare no conflicts of interest.

[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., Copar, S., Zumer, S. and Musevic, I. (2011) Science, 333, 62-65. http://dx.doi.org/10.1126/science.1205705 |

[7] |
Ravnik, M., Alexander, G.P., Yeomans, J.M. and Zumer, 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., Levi, G., Felidj, 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 Zumer, 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] | Bradac, Z., Kralj, S., Svetec, M. and Zumer, S. (2003) Physical Review E, 67, Article ID: 050702. |

[24] |
Svetec, M., Kralj, S., Bradac, Z. and Zumer, S. (2006) The European Physical Journal E, 20, 71-79. http://dx.doi.org/10.1140/epje/i2005-10120-9 |

[25] | Guimaraes, 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 |

Journals Menu

Contact us

+1 323-425-8868 | |

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2024 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.