Next Article in Journal
New Mononuclear Complex of Europium(III) and Benzoic Acid: From Synthesis and Crystal Structure Solution to Luminescence Emission
Next Article in Special Issue
The Genetic Algorithm: Using Biology to Compute Liquid Crystal Director Configurations
Previous Article in Journal
Green Synthesis of Hydrocalumite (CaAl-OH-LDH) from Ca(OH)2 and Al(OH)3 and the Parameters That Influence Its Formation and Speciation
Previous Article in Special Issue
Multiple Twisted Chiral Nematic Structures in Cylindrical Confinement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Annihilation of Highly-Charged Topological Defects

1
Department of Physics, Faculty of Natural Sciences and Mathematics, University of Maribor, Koroška 160, 2000 Maribor, Slovenia
2
Jožef Stefan International Postgraduate School, Jamova 39, 1000 Ljubljana, Slovenia
3
Faculty of Industrial Engineering, Šegova 112, 8000 Novo Mesto, Slovenia
4
Department of Physics, Case Western Reserve University, Cleveland, OH 44106, USA
5
Condensed Matter Physics, Jožef Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Crystals 2020, 10(8), 673; https://doi.org/10.3390/cryst10080673
Submission received: 24 June 2020 / Revised: 23 July 2020 / Accepted: 1 August 2020 / Published: 3 August 2020
(This article belongs to the Special Issue In Celebration of Noel A. Clark’s 80th Birthday)

Abstract

:
We studied numerically external stimuli enforced annihilation of a pair of daughter nematic topological defect (TD) assemblies bearing a relatively strong topological charge | m | = 3 / 2 . A Landau- de Gennes phenomenological approach in terms of tensor nematic order parameter was used in an effectively two-dimensional Cartesian coordinate system, where spatial variations along the z-axis were neglected. A pair of { m = 3 / 2 , m = 3 / 2 } was enforced by an appropriate surface anchoring field, mimicking an experimental sample realization using the atomic force microscope (AFM) scribing method. Furthermore, defects were confined within a rectangular boundary that imposes strong tangential anchoring. This setup enabled complex and counter-intuitive annihilation processes on varying relevant parameters. We present two qualitatively different annihilation paths, where we either gradually reduced the relative surface anchoring field importance or increased an external in-plane spatially homogeneous electric field E . The creation and depinning of additional defect pairs { 1 2 , 1 2 } mediated the annihilation in such a geometry. Furthermore, we illustrate the absorption of TDs by sharp edges of the confining boundary, accompanied by m = ± 1 / 4 1 / 4 winding reversal of edge singularities, and also E-driven zero-dimensional to one-dimensional defect core transformation.

Graphical Abstract

1. Introduction

Topological defects (TDs) [1] correspond to localized regions within ordered media where the relevant order parameter field is frustrated and topologically protected. In most cases, at the defect origin, the continuum field, which describes defectless equilibrium configurations, is not uniquely defined [2]. In general, TDs appear as a consequence of the universal concept of broken symmetry and are of interest in all branches of physics [1,3]. The key topological feature of a topological defect is fingerprinted in its topological charge [1,2], which is a conserved quantity. Corresponding topological charge conservation laws regulate transformations (decomposing, merging, depinning) of defect configurations.
Owing to their topological origin, the physics of TDs exhibits several universal features. Consequently, it is of interest to find physical systems in which TDs are relatively easily experimentally formed, manipulated, and observed, which allows detailed and controlled analysis of the behavior of interest. For example, in the last decades, several studies analyzed quench driven coarsening of TDs in different condensed matter systems to gain information on structural details in the early universe [3,4,5]. Condensed systems allow systematic and controlled experimental analysis, while related phenomena of the universe could be only indirectly and passively observed.
Nematic liquid crystals (NLC) [6] offer an ideal testbed to study TDs [7,8,9,10] due to their unique combination of liquid behavior, softness, experimentally-tractable characteristic scales of time and spatial responses, and optical anisotropy. Furthermore, NLC structures possessing TDs could be exploited in various (in particular electro-optical) applications [11,12]. The uniaxial nematic structure is commonly described by the mesoscopic nematic director field n ^ , which points along the local uniaxial NLC order, and the states ± n ^ are physically equivalent (the so-called head-to-tail invariance) [7]. In bulk, nematic equilibrium n ^ is homogeneously aligned along a symmetry breaking direction. Several studies have been performed on thin LC films, in which relatively simple quasi-two-dimensional (2D) behavior [13,14,15,16,17] is observed. Despite the simplicity, such systems display rich and complex behavior, in particular with respect to topological defect transformations. In 2D NLCs, only point defects exist. The role of topological charge is played by the winding number m [7] (also referred to as the Frank index). It determines the total rotation of n ^ on encircling the defect core counterclockwise. Due to the head-to-tail invariance, m can possess half-integer and whole integer numbers. Some typical TDs are shown in Figure 1. TDs bearing positive and negative m are commonly referred to as defects and antidefects, respectively, and pairs { m , m } in general tend to annihilate into a defectless configuration.
In general, high values of m are rarely realized. Namely, the free energy costs of a defect scales with m2 [16]. Consequently, if one enforces locally a defect bearing | m | > 1 / 2 , it tends to decompose [16,17,18] to elementary defects bearing | m | = 1 / 2 , where the charge conservation law is obeyed. Note that the core structure of single | m | = 1 / 2 disclination is well known. In terms of the tensor nematic order parameter it was originally determined by Schopohl and Sluckin [13]. Typical structural changes by traversing the defect core are shown in Figure 2. The orientational frustration is resolved via the order reconstruction (OR) mechanism [19]. In general, this mechanism is realized in cases where a relatively large orientational mismatch in nematic order is imposed on a scale comparable [19,20,21,22,23] to the nematic biaxial order parameter relaxation length ξ b .   The latter is relatively weakly temperature-dependent and is of order a few tens of nanometer [23]. The presence of the order reconstruction (OR) mechanism within the core was confirmed experimentally [24] in lyotropic chromonic LCs. Namely, such LCs exhibit a micron-scale nematic biaxial order parameter correlation length ξ b that provides experimental insight into the core structure, whose characteristic linear size is comparable to ξ b .
In this paper, we focus on the annihilation of surface enforced relatively highly charged {defect, antidefect} pair bearing the winding number | m | = 3 / 2 in an effectively 2D nematic cell confinement. We enforce the annihilation via different routes: (i) by changing effective elastic LC properties, or (ii) by increasing an external in-plane electric field strength. We demonstrate that annihilation proceeds via two qualitatively different channels, whose realization is rather counterintuitive.

2. Methods

We use the Landau-de Gennes mesoscopic approach in terms of the traceless and symmetric tensor nematic order parameter Q [6]. In its eigen frame, it is expressed as [20]
Q = i = 1 3 λ i e ^ i e ^ i ,
where λ i and e ^ i stand for Q eigenvalues and eigenvectors, respectively. Uniaxial states are commonly described by [6]
Q = S ( n ^ n ^ 1 3 I ) .
The nematic uniaxial order parameter S [ 1 / 2 , 1 ] quantifies the extent of fluctuations about the local nematic director n ^ , and S > 0   ( S < 0 ) reflects “prolate” (“oblate”) shape in the ellipsoidal presentation of the LC mesoscopic order. We henceforth refer to such configurations as positive uniaxial and negative uniaxial nematics.
Biaxial states could be assessed if elastic distortions are present. The degree of biaxiality is well measured by the biaxiality parameter [25]
β 2 = 1 6   tr ( Q 3 ) 2 tr ( Q 2 ) 3 [ 0 , 1 ] .
The limiting values β 2 = 0   and β 2 = 1 correspond to a uniaxial nematic state and a nematic order exhibiting the maximal biaxiality, respectively.
The eigenvalues can be parametrized as
λ 1 = 2 3   s cos ( ψ ) ,   λ 2 = 2 3   s cos ( ψ + π 3 ) ,   λ 3 = 2 3   s cos ( ψ π 3 ) ,
where
s = 3 2 tr Q 2 .
Note that the eigenvalues { λ 1 , λ 2 , λ 3 } can be expressed using two independent variables (i.e., {s, ψ }) due to the constraint i = 1 3 λ i 2 = 0 . The (s, ψ ) order parameter space is shown in Figure 2b. The isotropic order is reflected by s = 0 . Nematic configurations corresponding to   ψ = 0 ,   ψ = 2 π 3 ,   and ψ = 2 π 3 exhibit positive uniaxial states along e ^ 1 , e ^ 2 , and e ^ 3 , respectively. States ψ = π ,   ψ = π 3 ,   and ψ = π 3 correspond to the negative uniaxial states along e ^ 1 , e ^ 2 , and e ^ 3 . The remaining states are biaxial. In this parametrization it holds β 2 = s i n 2 ( 3 ψ ) . Any continuous change in nematic order is fingerprinted in a continuous change in the order parameter space.
In our study we will consider nematic structures in the Cartesian coordinate system (x,y,z), whose coordinate unit vectors are given by the triad ( e ^ x ,   e ^ y ,   e ^ y ). Of interest are quasi-two-dimensional structures, where the Q eigenvector along the z-axis remains fixed (i.e., e ^ 3 = e ^ z ). In this case, in addition to {s, ψ }, only one additional parameter (i.e. an angle   θ ) is needed to determine Q, describing the rotation of the Q eigenframe with respect to the Cartesian frame:
e ^ 1 = e ^ x s i n θ + e ^ y c o s θ ,
e ^ 2 = e ^ x c o s θ + e ^ y s i n θ .
Presentation of Q in terms of {s, ψ , θ } is useful to illustrate or guess different possible Q configuration variations along a given path in real space connecting two sites imposing different nematic order. An example is illustrated in Figure 2, where one crosses the core of | m | = 1 / 2 defect.
However, this Q representation does not provide a one-to-one mapping of each component (is not injective) [26]. For this reason, we used in simulations parametrization given by
Q = ( q 1 + q 3 ) e ^ x e ^ x + ( q 1 + q 3 ) e ^ y e ^ y + q 2 ( e ^ x e ^ y + e ^ y e ^ x ) 2 q 3 e ^ z e ^ z ,
in terms of variational fields { q 1 , q 2 , q 3 }. It holds λ 1 = q 3 + q 1 2 + q 2 2 , λ 2 = q 3 q 1 2 + q 2 2 , and λ 3 = 2 q 3 . Note that the exchange of eigenvalues [13] between λ 1 and λ 2 corresponds to the condition q 1 2 + q 2 2 = 0 . It was originally introduced to describe the elastic deformation within the core structure of | m | = 1 / 2 wedge defects [13]. The basic mechanism enabling this transformation is commonly referred to as order reconstruction [12,19].

2.1. Free Energy

We express the free energy F of a confined LC as a sum of nematic condensation ( f c ), elastic ( f e ) , external electric field ( f f ) , and surface ( f s ( i ) ) free energy density contributions:
F = ( f c + f e + f f ) d 3 r + i f s ( i ) d 2 r .
Here d 3 r and d 2 r stand for the volume and area measures, and the superscript (i) refers to the i-th confining substrate. The free energy densities are expressed using the standard Landau-de Gennes expansion [6] in Q . We take into account only the most essential symmetry allowed contributions to illustrate phenomena of interest [6,20,27]:
f c = a 0 ( T T * ) t r Q 2 b   t r Q 3 + c ( t r Q 2 ) 2 ,
f e = L 2 | Q | 2 ,
f f = ε 0 Δ ε 2 E . Q E ,
f s ( i ) = w ( i ) 2 tr ( Q Q s ) 2 .
Here a 0 ,   b ,   c are positive material-dependent constants, T * stands for the isotropic supercooling temperature, L is a positive representative nematic elastic constant (i.e., we use the single elastic constant approximation), E stands for an external electric field, Δ ε is the dielectric anisotropy (we limit to materials with Δ ε > 0 ), w ( i ) > 0 is the anchoring strength coefficient at the i-th confining surface, which enforces nematic order described by Q s [27].

2.2. Geometry of the Problem

We study nematic configurations inside a three-dimensional well with square cross section of length R and thickness h<<R. The geometry of the problem is depicted in Figure 3, where we use Cartesian coordinate system (x,y,z). At the lateral walls (Figure 3a) we enforce strong uniaxial tangential anchoring conditions:
Q s ( x = R / 2 , y , z ) = Q s ( x = R / 2 , y , z ) = S e q ( e ^ y e ^ y 1 3 I ) ,
Q s ( x , y = R / 2 , z ) = Q s ( x , y = R / 2 , z ) = S e q ( e ^ x e ^ x 1 3 I ) .
At the bottom (command) plate we enforce the nematic pattern consisting of a pair {defect, antidefect} of strength { m = m 1 = 3 / 2 ,   m = m 2 = 3 / 2 } placed at ( x 1 , y 1 ) and ( x 2 , y 2 ):
Q s ( x , y , z = 0 ) = S e q ( e ^ e ^ 1 3 I ) ,
e ^ = e ^ x c o s θ + e ^ y s i n θ ,   θ = i = 1 2 m i t a n 1 ( y y i x x i ) .
Typically, we set   ( x 1 = 0 , y 1 = 0.25   R )   and   ( x 2 = 0 , y 2 = 0.25   R ) . The opposite plate enforces degenerate tangential anchoring.
This setup mimics samples that could be realized experimentally using, for instance, the atomic force microscope (AFM) scribing method [28]. In a typical experimental setup, a nematic LC is confined within a thin plane-parallel cell, where at least one (command) surface imposes anchoring conditions inscribed via an AFM stylus [28], while the opposing surface imposes a degenerate tangential anchoring. Furthermore, in some simulations we switch on a spatially homogeneous electric field E = E e ^ y along the y-axis.
In simulations, we assume that the cell is thin enough so that the nematic structure is effectively two dimensional, exhibiting only variations in the (x,y) plane. Consequently, we use in numerical simulation the parametrization for the nematic tensor order parameter [22] given by Equation (7), where q 1 = q 1 ( x , y ) , q 2 = q 2 ( x , y ) , and q 3 = q 3 ( x , y ) are the variational parameters. In this parametrisation e ^ z   is always an eigenvector of   Q .
For presentational purposes, we introduce the dimensionless temperature t = ( T T * ) / ( T * * T * ) and material characteristic lengths [20]:
ξ f = L S * * ε 0 | Δ ε | E 2   ,
ξ b = 4 3 L c b 2 ( 1 + 1 t )   ,
d e = L w   .
The quantity ξ f stands for the external field nematic extrapolation length, which we express at the superheating temperature T = T * * , ξ b is the biaxial correlation length, d e is the surface extrapolation length, w stands for the anchoring strength imposed by the command surface, and S * * = 3 b / 16 c stands for the superheating bulk value of S. We henceforth utilize the relative anchoring strength of the command surface plate and of the external field E = E e ^ y   in terms of dimensionless ratios, respectively.
μ w = ξ b d e   ,     μ f = ( ξ b ξ f ) 2 ,
We minimize the free energy of the system with respect to the variational parameters q 1 , q 2 , and q 3 . The resulting Euler-Lagrange equations are solved numerically using the standard over-relaxation method. Technical and numerical details are given in [29].

3. Results

Of interest are different paths in which annihilation of a master surface anchoring enforced pair { m = 3 / 2 ,   m = 3 / 2 } annihilates in the geometry and boundary conditions described in the Methods section. The initial nematic configuration is shown in Figure 4a (biaxiality pattern) and Figure 5a (director field pattern). Each of the command surface imposed | m | = 3 / 2 defects at ( x 1 = 0 , y 1 = 0.25   R )   and   ( x 2 = 0 , y 2 = 0.25   R ) decay in three elementary | m | = 1 / 2 daughter defects. Due to boundary conditions, two additional defects m = 1 / 2 and m = 1 / 2 at the top and the bottom boundaries of Figure 4a, respectively, are introduced. In addition, at the lateral boundaries, relatively (non-singular) strongly elastically deformed regions exist, which are fingerprinted by relatively strong local biaxiality. Furthermore, sharp corners of their structure enforce | m | = 1 / 4 type surface distortions [30].
For latter convenience, we label the three positively charged daughter defects (emerging from m = 3 / 2 ) by P1, P2, P3, and the corresponding negatively charged daughter defects (emerging from m = 3 / 2 ) by N1, N2, N3 as shown in Figure 4a. Boundary condition generated defects are labeled by P4 (upper boundary) and N4 (bottom boundary). Furthermore, at some stages at lateral boundaries, elastic distortions are strong enough to trigger nucleation of creations of pairs {defect, antidefect}, see Figure 4c. We label these defect pairs by (N5,P5) and (N6,P6).
First, we enforce the annihilation by gradually decreasing the relative importance of the surface anchoring field. In practice, this could be achieved by decreasing temperature towards the second-order nematic-SmA phase transition temperature T N A , which is accompanied by divergence [6] in the nematic twist and bend elastic constants. The relative importance of elastic and surface penalties of a confined NLC is measured by the ratio R W / K , where R stands for the characteristic confinement length, and K ~ L S 2 is the representative Frank nematic elastic constant, and W measures the dimensional anchoring strength. In our simulations this ratio is measured by μ w . Therefore, on approaching T N A , the ratio μ w would decrease. The limit   μ w 1 corresponds to the strong anchoring regime, while   μ w < 1 corresponds to weak anchoring. The initial configuration, shown in Figure 4a, is evaluated for μ w = 1 , which is strong enough to prevent the annihilation of {3/2,-3/2} pair; however, weak enough to allow decomposition of | m | = 3 / 2 singularities. Such conditions are realized in typical samples prepared using AFM scribing method [17].
On decreasing   μ w , in the first stage, the closest oppositely charged daughter defects (N1,P1) begin to approach (see Supplementary movie Annihilation 1) until they mutually annihilate (See Figure 4b). Then the defect pairs (N2,P4), (P2,N4), and afterward (P3,N3) annihilate. During these processes, rearrangement of defects imposes sufficiently strong elastic distortions at the lateral boundaries to trigger the formation and depinning of two pairs of defects (P5,N5) and (P6,N6), as shown in Figure 4c. Finally, these defects are adsorbed by sharp edges of the confining boundary. Consequently, the winding number, characterizing strong distortions in corners, revert their winding number. This is evidently shown in Figure 5c,d, where one observes ± 1 / 4 1 / 4 winding number transformations of the director field profiles within the edges.
Next, we trigger annihilation by applying and gradually increasing a spatially homogeneous external electric field E = E e ^ y . The initial configuration for E = 0 is shown in Figure 6a, where we set ( x 1 = 0 , y 1 = 0.3   R )   and   ( x 2 = 0 , y 2 = 0.3   R ) . On increasing E (see Supplementary movie Annihilation 2), in the first stage, the inner (N1,P1) defect approach and finally annihilate (Figure 6c). The remaining defects annihilate in a qualitatively different way compared to the process shown in Figure 4, where lateral defect pairs (N5,P5) and (N6,P6) are strongly involved. Namely, the oppositely charged daughter and lateral defects approach each other and finally annihilate (Figure 6d). Furthermore, initially localized | m | = 1 / 2 structures of P4 and P5, clearly visible in Figure 6a–e, transform to nonlocalized order reconstruction planes. This defect core “explosion” was studied in detail in [31]. Here we summarize the key stages of the process. In this E-driven transformation, the initially spherically shaped biaxial core of a defect begins to grow as it is pushed towards a limiting wall (in our setting y = R/2 boundary for P4, and y = −R/2 boundary for M4) on increasing E. For E > 0, the ellipsoidal core-shape becomes progressively elongated along the x-axis. At some critical value of E, the core diverges (i.e., its longest linear size equals R), and a planar order reconstruction plane is formed. In this case order reconstruction takes form near the bottom and top confining plate. The corresponding nematic director profile, reflecting the Q eigen frame orientation, is shown Figure 7b. Note that on crossing the OR plane the LC molecules experience a mesoscopic shape transformation, similar to the one shown in Figure 2c along the direction 1-5. The OR plane is formed at   μ f ~ 3.5 , which for typical LCs [6] corresponds to E ~ 10 8   V/m.

4. Discussion

We studied numerically the annihilation of master surface enforced highly charged pair {defect, antidefect}, bearing topological charges { 3 / 2 , 3 / 2 } . A Landau-de Gennes mesoscopic approach in terms of the tensor nematic order parameter was used. We assumed quasi-two-dimensional order, which in practice is sensible for relatively thin nematic cells. We chose geometrical setups and boundary conditions that could be realized experimentally in future work. For example, highly charged defect structures could be enforced using, e.g., the Atomic Force Measurement scribing method [30], or plasmonic photoalignment technique [32]. In simulations we enforce a pair of | m | = 3 / 2 defects via the command surface scribed pattern. We used anchoring strength coefficients that are strong enough to prevent annihilation of the pair, however weak enough to allow decomposition [33] of each surface pattern favored | m | = 3 / 2 defect into three daughter elementary defects bearing | m | = 1 / 2 . Furthermore, we confined the pair in a rectangular confinement region. One would expect mutual annihilation of the two daughter assemblies on varying relevant control parameters. However, the chosen setup enabled a nontrivial annihilation scenario of TDs, also involving creation of additional defects, absorption of TDs within edges of the confining substrate, and one dimensional “explosion” of defect cores. We illustrate two qualitatively different annihilation channels, either by progressively weakening the relative importance of the surface anchoring field or increasing external in-plane electric field. Both processes could be realized experimentally by using relatively simple experimental setups. In our setting, the reference structures exhibit eight elementary TDs, where two additional defects are introduced by the rectangular confinement. In the anchoring weakening-driven annihilation, only four daughter defects were mutually annihilated, and the remaining two were annihilated with confinement enabled TDs. On the contrary, in E-driven annihilation only one pair of daughter defects were mutually annihilated and the remaining four defects were annihilated via TDs, which were created via creation and depinning of additional pairs {defect, antidefect}.
Note that in realistic LC samples the characteristic nematic elastic constants are different from each other. Nematic elasticity is commonly expressed in terms of splay ( K 11 ) , twist ( K 22 ) , bend ( K 33 ) , and saddle-splay ( K 24 ) elastic constants, introduced in the Frank-Oseen uniaxial description [34]. In our geometry the elastic distortions involving twist and saddle-splay deformations are absent. Therefore, only K 11 and K 33 are expected to make an impact on defect structures. In conventional LCs it holds K 33 > K 11 . In this case the defect rearrangement would be carried out via LC structural variations that prefer the splay nematic deformation. Furthermore, this elastic anisotropy would break the mirror symmetry along the y-axis observed in Figure 4a. Namely, the core structures of m = 1 / 2 and | m | = 1 / 2 defects would be different, because the relative contributions of bend and splay deformations are different in these defects.
On approaching a second order N-SmA phase transition K 33 diverges, whereas K 11 remains finite. Because any arrangement of TDs in our study involves bend elastic distortions, TDs would annihilate on increasing K 33 . This elastic anisotropy would, for sure, at least quantitatively influence the annihilation scenario. The defect rearrangement would be realized via LC rearrangements preferring the splay deformation. However, for our scenario we do not expect a qualitatively different rearrangement of TDs from those presented in Figure 4 and Figure 6.
Nevertheless, it is a demanding task to obtain a numerically “realistic” result of elastic anisotropy on rearrangement of TDs. In this paper we work in the approximation of equal elastic constants, represented by a temperature independent constant L. In reality, several different elastic constants would appear, which determine the relative importance of different symmetry-allowed elastic contributions expressed in terms of Q and its spatial derivatives. For example, in Cartesian coordinates (x,y,z) = ( x 1 , x 2 , x 3 ), the expansion up to the 2nd order in Q is commonly expressed as [35,36]
f e = f e ( 2 ) = L 1 Q i j , k Q i j , k + L 2 Q i j , j Q i k , k + L 3 Q i j , k Q i k , j + L 4 Q i j Q i j , k k + L 5 Q j k Q i k , j k .
Here we use the summation convention over repeated indices, commas indicate spatial derivatives (e.g., Q i j , k = Q i j x k ), and the superscript (2) in f e ( 2 ) indicates the 2nd order expansion contributions in the elastic free energy density. Far from defect cores, the constants L 1 - L 5 could be mapped to the more familiar Frank elastic constants. However, for the expansion up to the 2nd order in Q it holds that K 11 = K 33 [35,36], while in most LCs these constants can be significantly different. Higher order terms need to be added [35,36] to lift this degeneracy, for example f e = f e ( 2 ) + L 6 Q k l Q i j , l Q i j , k . However, according to recent studies in lyotropic chromonic LCs [24], even taking into account the above mentioned elastic anisotropy, one could not exactly reproduce experimentally the observed core structures of | m | = 1 / 2 TDs. Therefore, current theoretical modeling does not correctly describe cores of single | m | = 1 / 2 TDs, and the core structure details might influence annihilation scenarios. For this reason, we performed our study at the proof-of-principle level.
In our selected geometry pairs of defects {−1/2,1/2} are typically created as an intermediate step and they actively mediate external stimulus driven annihilation processes. To estimate the free energy costs of a pair {−1/2,1/2} creation we use the Lyuksyutov constraint [26,37] in which the cubic term in f c is neglected. This is sensible approximation for LCs exhibiting weakly first order I-N phase transition. With this in mind, it follows f c ~ a 0 ( T T * ) t r Q 2 + c ( t r Q 2 ) 2 , where f c is minimized in the nematic phase for t r Q 2 = a 0 ( T * T ) 2 c ,   corresponding to f c ~ a 0 2 ( T T * ) 2 4 c .   The main penalty of introducing a defect is “melting” of LC order within the core of TDs. Note that for a pair {−1/2,1/2}, the nematic director in the far field is spatially homogeneous because the total charge of the pair equals to zero. The free energy cost of introducing a single | m | = 1 / 2   defect of length h is roughly F ~ | f c | V a 0 2 ( T T * ) 2 4 c π ξ b 2 h , where V ~ π ξ b 2 h estimates the defect core volume. For a typical nematic LC it holds that a 0 ~ 0.05   10 6 J K m 3 , c ~   10 6 J m 3 ,     ξ b ~ 30 nm [23]. Therefore, deep in the nematic phase it costs roughly 2 F / ( k b T * )   ~ 10 3 to create a pair of defects, where we take | T T * | = 10 K, h ~ ξ b , T * ~ 330 K, and k b is the Boltzmann constant.

5. Conclusions

In summary, we studied the external stimulus-enforced annihilation of highly-charged assembly of command surface enforced TDs. In bulk, these TDs would simply mutually annihilate on varying an appropriate control parameter. However, we showed that in appropriate confinement these processes could be extremely complex and display qualitatively different realizations. In particular, we showed that when the characteristic confinement size and separation of highly charged defects are comparable, additional defect pairs are created and involved in the annihilation processes. In these events the total topological charge of the system is conserved in accordance with topological charge conservation laws. Our simulations exhibit a rich diversity of defect configurations for which only a relatively small set of control parameters was varied.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4352/10/8/673/s1, Surface anchoring and an external electric field enforced annihilations pf TDs. “Annihilation1” movie shows the β 2 ( x , y ) texture evolution of the annihilation process of command surface enforced assemblies {3/2,−3/2} TDs on decreasing the dimensionless surface anchoring field strength from μ w = 1 , to μ w = 0 , Some representative stages are shown in Figure 4. “Annihilation2” movie shows the β 2 ( x , y ) texture evolution of the annihilation process of command surface enforced assemblies {3/2,−3/2} TDs on increasing the dimensionless external electric field strength from μ f = 0 , to μ f = 3.5 , Some representative stages are shown in Figure 6.

Author Contributions

S.K. and P.K. proposed the research; S.K. and C.R. guided the research; E.K., P.K., and M.A. did numerical simulations; E.K. and P.K. prepared figures. All authors were involved in writing the paper. All authors have read and agreed to the published version of the manuscript.

Funding

S.K. acknowledges the support of a Slovenian Research Agency (ARRS) grant P1–0099. C.R. was supported by the National Science Foundation Condensed Matter Physics program under grant DMR-1901797.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mermin, N. The Topological Theory of Defects in Ordered Media. Rev. Mod. Phys. 1979, 51, 591–648. [Google Scholar] [CrossRef]
  2. Kurik, M.V.; Lavrentovich, O.D. Defects in Liquid Crystals: Homotopy Theory and Experimental Studies. Sov. Phys. Usp. 1988, 31, 196–224. [Google Scholar] [CrossRef]
  3. Zurek, W.H. Cosmological experiments in superfluid helium? Nature 1985, 317, 505–508. [Google Scholar] [CrossRef] [Green Version]
  4. Kibble, T.W.B. Topology of cosmic domains and strings. J. Phys. A Math. Gen. 1976, 9, 1387–1398. [Google Scholar] [CrossRef]
  5. Chesler, P.M.; García-García, A.M.; Liu, H. Defect Formation beyond Kibble-Zurek Mechanism and Holography. Phys. Rev. X 2015, 5, 021015. [Google Scholar] [CrossRef] [Green Version]
  6. De Gennes, P.G.; Prost, J. The Physics of Liquid Crystals; Clarendon: Oxford, UK, 1994. [Google Scholar]
  7. Lavrentovich, O.D. Topological defects in dispersed liquid crystals, or words and words around liquid crystal drops. Liq. Cryst. 1998, 24, 117–125. [Google Scholar] [CrossRef]
  8. Tkalec, U.; Ravnik, M.; Čopar, S.; Žumer, S.; Muševič, I. Reconfigurable Knots and Links in Chiral Nematic Colloids. Science 2011, 333, 62–65. [Google Scholar] [CrossRef] [Green Version]
  9. Martinez, A.; Ravnik, M.; Lucero, B.; Visvanathan, R.; Žumer, S.; Smalyukh, I.I. Mutually tangled colloidal knots and induced defect loops in nematic fields. Nat. Mater. 2014, 13, 258–263. [Google Scholar] [CrossRef]
  10. Foster, D.; Kind, C.; Ackerman, P.J.; Tai, J.S.B.; Dennis, M.R.; Smalyukh, I.I. Two-dimensional skyrmion bags in liquid crystals and ferromagnets. Nat. Phys. 2019, 15, 655–659. [Google Scholar] [CrossRef]
  11. Luo, C.; Majumdar, A.; Erban, R. Multistability in planar liquid crystal wells. Phys. Rev. E 2012, 85, 061702. [Google Scholar] [CrossRef] [Green Version]
  12. Palffy-Muhoray, P. The diverse world of liquid crystals. Phys. Today 2007, 60, 54–60. [Google Scholar] [CrossRef]
  13. Schopohl, N.; Sluckin, T.J. Defect Core Structure in Nematic Liquid Crystals. Phys. Rev. Lett. 1987, 59, 2582–2584. [Google Scholar] [CrossRef] [PubMed]
  14. Vitelli, V.; Turner, A.M. Anomalous coupling between topological defects and curvature. Phys. Rev. Lett. 2004, 93, 215301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Bowick, M.; Nelson, D.R.; Travesset, A. Curvature-induced defects unbinding in toroidal geometries. Phys. Rev. E 2004, 69, 041102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Stannarius, R.; Harth, K. Defect Interactions in Anisotropic Two-Dimensional Fluids. Phys. Rev. Lett. 2016, 117, 157801. [Google Scholar] [CrossRef]
  17. Kralj, S.; Murray, B.S.; Rosenblatt, C. Decomposition of Strongly Charged Topological Defects. Phys. Rev. E 2017, 95, 042702. [Google Scholar] [CrossRef] [Green Version]
  18. Svenšek, D.; Žumer, S. Instability Modes of High-Strength Disclinations in Nematics. Phys. Rev. E 2004, 70, 061707. [Google Scholar] [CrossRef] [Green Version]
  19. Palffy-Muhoray, P.; Gartland, E.C.; Kelly, J.R. A new configurational transition in inhomogeneous nematics. Liq. Cryst. 1994, 16, 713–718. [Google Scholar] [CrossRef]
  20. Bisi, F.; Gartland, E.C.; Rosso, R.; Virga, E.G. Order reconstruction in frustrated nematic twist cells. Phys. Rev. E 2003, 68, 021707. [Google Scholar] [CrossRef] [Green Version]
  21. Barberi, R.; Ciuchi, F.; Lombardo, G.; Bartolino, R.; Durand, G.E. Time resolved experimental analysis of the electric field induced biaxial order reconstruction in nematics. Phys. Rev. Lett. 2004, 93, 137801. [Google Scholar] [CrossRef]
  22. Kralj, S.; Rosso, R.; Virga, E.G. Finite-size effects on order reconstruction around nematic defects. Phys. Rev. E 2010, 81, 021702. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Lombardo, G.; Amoddeo, A.; Hamdi, R.; Ayeb, H.; Barberi, R. Biaxial surface order dynamics in calamitic nematics. Eur. Phys. J. E 2012, 35, 32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhou, S.; Shiyanovskii, S.V.; Park, H.-S.; Lavrentovich, O.D. Fine structure of the topological defect cores studied for disclinations in lyotropic chromonic liquid crystals. Nat. Commun. 2017, 8, 14974. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kaiser, P.; Wiese, W.; Hess, S. Stability and Instability of an Uniaxial Alignment Against Biaxial Distortions in the Isotropic and Nematic Phases of Liquid Crystals. J. Non-Equilib. Thermodyn. 1992, 17, 153–169. [Google Scholar] [CrossRef]
  26. Kralj, S.; Virga, E.G.; Žumer, S. Biaxial torus around nematic point defects. Phys. Rev. E 1999, 60, 1858–1866. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Nobili, M.; Durand, G. Disorientation-induced disordering at a nematic-liquid-crystal-solid interface. Phys. Rev. A 1992, 46, R6174–R6177. [Google Scholar] [CrossRef]
  28. Murray, B.S.; Pelcovits, R.A.; Rosenblatt, C. Creating arbitrary arrays of two-dimensional topological defects. Phys. Rev. E 2014, 90, 052501. [Google Scholar] [CrossRef] [Green Version]
  29. Kralj, S.; Majumdar, A. Order reconstruction patterns in nematic liquid crystal wells. Proc. R. Soc. A. 2014, 470, 20140276. [Google Scholar] [CrossRef]
  30. Senyuk, B.; Liu, O.; Yuan, Y.; Smalyukh, I.I. Edge pinning and transformation of defect lines induced by faceted colloidal rings in nematic liquid crystals. Phys. Rev. E 2016, 93, 062704. [Google Scholar] [CrossRef] [Green Version]
  31. Ambrožič, M.; Kralj, S.; Virga, E.G. Defect-enhanced nematic surface order reconstruction. Phys. Rev. E 2007, 75, 031708. [Google Scholar] [CrossRef] [Green Version]
  32. Guo, Y.; Jiang, M.; Peng, C.; Sun, K.; Yaroshchuk, O.; Lavrentovich, O.; Wei, Q.H. High-Resolution and High-Throughput Plasmonic Photopatterning of Complex Molecular Orientations in Liquid Crystals. Adv. Mater. 2016, 28, 2353. [Google Scholar] [CrossRef] [PubMed]
  33. Murray, B.S.; Kralj, S.; Rosenblatt, C. Decomposition vs. Escape of Topological Defects in a Nematic Liquid Crystal. Soft Matter. 2017, 13, 8442. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Selinger, J.V. Interpretation of saddle-splay and the Oseen-Frank free energy in liquid crystals. Liq. Cryst. Rev. 2018, 6, 129–142. [Google Scholar] [CrossRef] [Green Version]
  35. Longa, L.; Monselesan, D.; Trebin, H.R. An extension of the Landau-Ginzburg-de Gennes theory for liquid crystals. Liq. Cryst. 1987, 2, 769–796. [Google Scholar] [CrossRef] [Green Version]
  36. Kralj, S.; Žumer, S. The saddle-splay surface elasticity of nematic structures confined to a cylindrical capillary. Phys. Rev. E. 1995, 51, 366–379. [Google Scholar] [CrossRef]
  37. Lyuksyutov, I.F. Topological instability of singularities at small distances in nematics. Sov. Phys. JETP 1987, 48, 178–179. [Google Scholar]
Figure 1. Examples of topological defects (TDs) bearing (a) m = 1/2; (b) m = 1; (c) m = 3/2; (d) m = −1/2; (e) m = −1; (f) m = −3/2.
Figure 1. Examples of topological defects (TDs) bearing (a) m = 1/2; (b) m = 1; (c) m = 3/2; (d) m = −1/2; (e) m = −1; (f) m = −3/2.
Crystals 10 00673 g001
Figure 2. Core structure of a typical m = 1/2 disclination. (a) The biaxial spatial variation, measured via the biaxiality parameter β 2 (see Equation (3)). The core structure is characterized by a rim, where β 2 = 1 . The rim is circular in the approximation of equal elastic constants. The center of the core exhibits negative uniaxiality. (b) Characteristic changes in the ( s , ψ ) phase space on crossing the center of disclination. (c) Schematic variation of the nematic mesoscopic local profile within the core. Here numbers 1 to 5 mark points representing mesoscopic changes on traversing the defect core along a chosen direction.
Figure 2. Core structure of a typical m = 1/2 disclination. (a) The biaxial spatial variation, measured via the biaxiality parameter β 2 (see Equation (3)). The core structure is characterized by a rim, where β 2 = 1 . The rim is circular in the approximation of equal elastic constants. The center of the core exhibits negative uniaxiality. (b) Characteristic changes in the ( s , ψ ) phase space on crossing the center of disclination. (c) Schematic variation of the nematic mesoscopic local profile within the core. Here numbers 1 to 5 mark points representing mesoscopic changes on traversing the defect core along a chosen direction.
Crystals 10 00673 g002
Figure 3. The geometry of the problem. (a) Nematic liquid crystals (NLC) cell of length R and thickness h<<R. At the lateral walls we enforce strong tangential anchoring (along e ^ x or e ^ y ). At the bottom command plate a pair { m = 3 / 2 ,   m = 3 / 2 } of surface defects is enforced. The top plate imposes degenerate tangential anchoring. (b) The command plate imposed nematic pattern in the (x,y) plane.
Figure 3. The geometry of the problem. (a) Nematic liquid crystals (NLC) cell of length R and thickness h<<R. At the lateral walls we enforce strong tangential anchoring (along e ^ x or e ^ y ). At the bottom command plate a pair { m = 3 / 2 ,   m = 3 / 2 } of surface defects is enforced. The top plate imposes degenerate tangential anchoring. (b) The command plate imposed nematic pattern in the (x,y) plane.
Crystals 10 00673 g003
Figure 4. Annihilation by decreasing the surface anchoring field strength. The degree of biaxiality β 2 in ( x , y ) plane is shown. We impose via the command surface a pair of point defects with charges m { 3 / 2 , 3 / 2 } , which decomposes to elementary | m | = 1 / 2 defects. With decreasing anchoring strength, the defects annihilate. ( a )   μ w = 1.0 , ( b )   μ w = 0.6 , ( c )   μ w = 0.4 , ( d )   μ w = 0.3 , ( e )   μ w = 0.1 ,   ( f )   μ w = 0.0 . R / ξ b = 25 , t = 8 . Selected corresponding qualitatively different director profiles are shown in Figure 5.
Figure 4. Annihilation by decreasing the surface anchoring field strength. The degree of biaxiality β 2 in ( x , y ) plane is shown. We impose via the command surface a pair of point defects with charges m { 3 / 2 , 3 / 2 } , which decomposes to elementary | m | = 1 / 2 defects. With decreasing anchoring strength, the defects annihilate. ( a )   μ w = 1.0 , ( b )   μ w = 0.6 , ( c )   μ w = 0.4 , ( d )   μ w = 0.3 , ( e )   μ w = 0.1 ,   ( f )   μ w = 0.0 . R / ξ b = 25 , t = 8 . Selected corresponding qualitatively different director profiles are shown in Figure 5.
Crystals 10 00673 g004
Figure 5. Annihilation by decreasing the surface anchoring field strength. Selected qualitatively different director profiles are shown. ( a )   μ w = 1.0 , ( b )   μ w = 0.3 , ( c )   μ w = 0.1 , ( d )   μ w = 0.0 .   R / ξ b = 25 , t = 8 . TDs bearing m = 1/2 and m = −1/2 are indicated by red and blue circles, respectively.
Figure 5. Annihilation by decreasing the surface anchoring field strength. Selected qualitatively different director profiles are shown. ( a )   μ w = 1.0 , ( b )   μ w = 0.3 , ( c )   μ w = 0.1 , ( d )   μ w = 0.0 .   R / ξ b = 25 , t = 8 . TDs bearing m = 1/2 and m = −1/2 are indicated by red and blue circles, respectively.
Crystals 10 00673 g005
Figure 6. Annihilation by increasing an external in-plane electric field. 2 D biaxiality profiles β 2 ( x , y ) are shown. We gradually increase the external electric field in the y-direction. ( a )   μ f = 0 , ( b )     μ f = 0.4 , ( c )     μ f = 0.7 , ( d )     μ f = 0.85 ,     ( e )     μ f = 1 , ( f )     μ f = 3.5 .   R / ξ b = 25 , t = 8 . Selected corresponding qualitatively different director profiles are shown in Figure 7.
Figure 6. Annihilation by increasing an external in-plane electric field. 2 D biaxiality profiles β 2 ( x , y ) are shown. We gradually increase the external electric field in the y-direction. ( a )   μ f = 0 , ( b )     μ f = 0.4 , ( c )     μ f = 0.7 , ( d )     μ f = 0.85 ,     ( e )     μ f = 1 , ( f )     μ f = 3.5 .   R / ξ b = 25 , t = 8 . Selected corresponding qualitatively different director profiles are shown in Figure 7.
Crystals 10 00673 g006
Figure 7. Annihilation by increasing an external in-plane electric field. Selected qualitatively different director profiles are shown. ( a )   μ f = 0.7 , ( b )   μ f = 3.5 .   R / ξ b = 25 , t = 8 . TDs bearing m = 1/2 and m = −1/2 are indicated by red and blue circles, respectively.
Figure 7. Annihilation by increasing an external in-plane electric field. Selected qualitatively different director profiles are shown. ( a )   μ f = 0.7 , ( b )   μ f = 3.5 .   R / ξ b = 25 , t = 8 . TDs bearing m = 1/2 and m = −1/2 are indicated by red and blue circles, respectively.
Crystals 10 00673 g007

Share and Cite

MDPI and ACS Style

Klemenčič, E.; Kurioz, P.; Ambrožič, M.; Rosenblatt, C.; Kralj, S. Annihilation of Highly-Charged Topological Defects. Crystals 2020, 10, 673. https://doi.org/10.3390/cryst10080673

AMA Style

Klemenčič E, Kurioz P, Ambrožič M, Rosenblatt C, Kralj S. Annihilation of Highly-Charged Topological Defects. Crystals. 2020; 10(8):673. https://doi.org/10.3390/cryst10080673

Chicago/Turabian Style

Klemenčič, Eva, Pavlo Kurioz, Milan Ambrožič, Charles Rosenblatt, and Samo Kralj. 2020. "Annihilation of Highly-Charged Topological Defects" Crystals 10, no. 8: 673. https://doi.org/10.3390/cryst10080673

APA Style

Klemenčič, E., Kurioz, P., Ambrožič, M., Rosenblatt, C., & Kralj, S. (2020). Annihilation of Highly-Charged Topological Defects. Crystals, 10(8), 673. https://doi.org/10.3390/cryst10080673

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop