# Field induced memory effects in random nematics.

1. IntroductionImpact of weak disorder on statistical properties in phases exhibiting continuous symmetry breaking is of interest for different branches in physics [1-3]. Such systems are of particular interest due to memory effects which they can exhibit [2, 4, 5]. Pioneering studies have been done in magnetic systems [1,6]. Lately it has been shown that various liquid crystal (LC) phases are extremely adequate for such studies [7]. Namely, LCs are relatively easily experimentally accessible, they display rich diversity of different configurations, and there exists plenty of different means how qualitatively different types of disorder could be imposed to LCs. Consequently, LCs are ideal testing ground to investigate fundamental problems related to impact of weak disorder on phases reached via a symmetry breaking phase transition. In addition, randomly perturbed LCs are also technologically relevant, specially due to their unique combination of liquid character, softness, and optical anisotropy. Randomly perturbed LCs could be exploited for various devices, particularly related to electrooptic applications.

Uniaxial nematic phase represents the simplest LC configuration exhibiting only orientational ordering [8]. At mesoscopic level this ordering is commonly described by the nematic director field [??] and the scalar nematic order parameter S [8]. The unit pseudovector [??] points along local average uniaxial LC orientational ordering, exhibiting head-to-tail symmetry (i.e., states [+ or -][??] are physically equivalent). On the other hand S measures extent of fluctuations about [??]. Absence of thermal fluctuations in local orientational order are fingerprinted in S = 1. In bulk equilibrium nematic phase is spatially homogeneous where [??] is aligned along a single symmetry direction. In thermotropic LCs, to which we restrict our study to, nematic long range ordering (LRO) is obtained via the first-order phase transition from the isotropic phase, by decreasing temperature. The latter phase corresponds to ordinary liquid ordering, where S = 0.

There are several different ways how random-type disorder could be experimentally imposed to LCs [7]. For such purposes LCs are either immersed in various porous materials (e.g., aerogels, Vycor glasses, Russian glasses, controlled pore glasses, and millipore membranes...) or by mixing them with different nanoparticles. In the latter case aerosil nanoparticles are commonly used [9]. Namely, by changing their concentration, three qualitatively different types of disorders could be realized [10].

Resulting randomly perturbed nematic structures depend on disorder strength. According to the Imry-Ma theorem [1, 6], one expects that weak random field-type disorder stabilizes a domain-type pattern characterized by a single length [[xi].sub.d] obeying the scaling law [[xi].sub.d] ~ [D.sup.-2/(4-d)]. Here D stands for the disorder strength and d is the space dimensionality. The resulting phase should possess short range order (SRO). However, several subsequent studies suggest that the quasi long range order (QLRO) or even LRO might be established [2] instead.

In general, randomly perturbed nematic structures possessing SRO or QLRO exhibit glassy features [2]. Namely, strong enough disorder could stabilize enormous large variety of different metastable structures in free energy landscape. Recent experimental and simulation studies in different randomly perturbed nematics support this expectation [4, 9, 11-13]. For example, in LC-aerosil mixtures, it was demonstrated that LCs could memorize external electric or magnetic field induced orientation [4, 5, 14, 15]. Note that systems exposed to relatively strong disorder are known to exhibit several universalities, independent of systems' microscopic details [2]. However, in case of relatively weak disorder system, specific details are expected to emerge, which further increases diversity of potentially possible phenomena. In addition, glassy features in nematics could be apparently different with respect to magnetic analogues due to nematic head-to-tail symmetry, which is not observed in magnetism. Consequently, perturbed nematics structures in general display complex tangles of topologically stable line disclinations characterized by winding number m = -1/2 [5]. On the contrary, line defects could not be topologically stabilized in magnetic systems. Namely, local magnetization could exhibit only line defects with winding number m = 1, which could avoid singularity in orientational ordering by escaping in "third dimension" Furthermore, in LCs, thermal fluctuations play relatively important role due to their soft character. Consequently, energy barriers in free energy landscape could be relatively overcrossed easier with respect to magnetic analogues.

In the present paper we study disorder enabled memory effects in randomly perturbed LCs. Of our particular interest is external field driven memorization. Our study is numerical, where we use lattice version of modified Random Anisotropy Nematic model [16-18]. The plan of the paper is as follows. In Section 2 we present the model. In Section 3 results of our simulations are presented. In the last section we summarize results. Some technical details related to numerical method and mesoscopic interpretation of results are collected in two appendices.

2. Model and Numerical Approach

We use lattice-type simulations where the cubic simulation box of volume [L.sup.3] contains N = [(L/[a.sub.0]).sup.3] sites. Lengths L and [a.sub.0] represent the characteristic linear size of a system and the shortest distance between adjacent lattice sites, respectively. Typically, we set L/[a.sub.0] = 70. Nematic LC orientational order at each lattice site is presented by the unit vector [[??].sub.i], to which we henceforth refer to as a nematic spin. The so-called head-to-tail symmetry [+ or -] [[??].sub.i] is considered to mimic typical nematic LC behavior at the mesoscopic level. In addition, we suppose that each site might be populated by a rod-like object with probability p. Their local quenched orientations e; are determined randomly, exhibiting isotropic orientational distribution, where the head-to-tail symmetry and [absolute value of ([[??].sub.i])] = 1 are imposed. We refer to these pseudovectors as impurity spins. Note that sites occupied by impurities are not accessible to nematic spins.

Using random number generator we first set static impurity positions and their orientations. For such initial configurations we then determine fixed point configurations of nematic spins, corresponding to metastable or globally stable structures.

2.1. Hamiltonian. The interactions within a system are described by the following Hamiltonian [5,11, 19, 20]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)

The first part of (1) describes short range interaction among nematic spins, where [J.sub.ij] = J > 0 for neighboring spins (each site has 6 first neighbors) and otherwise [J.sub.ij] = 0. The factor 1/2 is added in order to avoid double counting of pair-interaction contributions. The second term takes into account presence of a global uniform field [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], which we align along the x-axis of the Cartesian coordinate system {x, y, z}. The last term takes into account short range interaction between an impurity at a site i and a nematic spin at a site j of strength W > 0. Therefore, if impurity is not present at a site i then [W.sub.ij] = 0. For such choice of interactions a nematic spin tends to be oriented parallel to a neighboring nematic or impurity spin. Furthermore, a finite external field tends to align nematic spins parallel to its direction. A resulting configuration compromises couplings represented by interaction strengths J, B, and W. Simulation details are described in Appendix A.

In simulations we use three qualitatively different histories of systems. We imposed either temperature-quenched history (TQH), field-quenched history (FQH), or temperature-annealed history (AH). The first two of these limits represent extreme histories encountered in typical experimental studies. In case of TQH (FQH) we quench system from isotropic initial distribution (homogeneously alignment along [??]) of nematic spins. In case of AH we calculate a fixed point configuration {[[??].sub.i]} at a given temperature originating from an initial fixed point nematic structure which was calculated at a slightly different temperature. Simulations were repeated [N.sub.rep] ~ 10 times for a given set of parameters (i.e., W, p, and a chosen initial condition).

2.2. Monitored Quantities. From obtained fixed point configurations we calculate the global order parameter [[Q.bar].sup.(g)] and local order parameter tensor [[Q.bar].sup.(l)]. They estimate a global and local degree of orientational ordering, respectively. We define them as follows:

[[Q.bar].sup.(g)] = 1/N [summation over (i)] 1/2 [([??].sub.i] [cross product] [[??].sub.i] - [I.bar]/3), (2a)

[[Q.bar].sup.(l)] = 1/N [summation over (i)] 1/[N.sub.nn] [summation over (j)] 1/2 ([[??].sub.i] [cross product] [[??].sub.j] - [I.bar]/3), (2b)

where the index i runs over the lattice sites occupied by nematic spins, j runs over all the neighbours (their number is [N.sub.nn]) of a site i, [I.bar] is the identity tensor, and [cross product] marks the tensor product. The average local scalar order parameter s is defined as the largest eigenvector of [[Q.bar].sup.(l)]. On the other hand the average global scalar order parameter S is defined as the largest eigenvector of [[Q.bar].sup.(g)].

We further define the orientational two-point correlation function as [11]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)

where <...> is the statistical average of the scalar products of nematic spin pairs separated by a distance r = [absolute value of ([[??].sub.ij])] = [absolute value of ([[??].sub.i] - [[??].sub.j])]. In case of short range order (SRO) or quasi long range order (QLRO), it holds G(r [right arrow] [infinity]) [right arrow] 0. In case of long range order (LRO), it follows G(r [right arrow] [infinity]) [right arrow] [G.sub.[infinity]] ~ [S.sup.2]. In general, one expects an exponential decay towards a saturated value of G(r) on increasing r for both LRO and SRO. On the other hand, for QLRO, algebraic decay of correlations is expected; that is, G(r) [varies] [r.sup.-[alpha]].

In order to obtain structural details from G(r) for a finite system the correlation function is fitted using the empirical ansatz [11]:

[G.sup.(1)] (r) = a[e.sup.-kr]/[r.sup.[alpha]] + b, (4)

where a, b, k, and [alpha] are adjustable parameters. The quantity k estimates an average linear size [xi] = 1/k of relatively well-correlated regions, referred to as domains.

Note that, in case of QLRO, the decay of correlations with distance is relatively weak, and finite-size effects are expected to be important. Consequently, the power law coefficient a determining algebraic decay of correlations is estimated by the following equation [11,21, 22]:

S [varies] [L.sup.-[beta]], (5)

where [beta] ~ [alpha]/2, see (4). In case of SRO and LRO, one expects [beta] ~ 3/2 and [beta] ~ 0, respectively, and QLRO is fingerprinted by 0 < [beta] <3/2 [11, 23].

3. Results

We calculated fixed point configurations of nematic spins on varying parameters {W, p, T, B} and histories of systems in three-dimensional systems. We limited to concentrations of impurities below the percolation threshold, corresponding to the critical concentration [p.sub.c] ~ 0.3 [14]. Of particular interest is impact of external field B on degree of ordering and memory effects.

In Figure 1 we plot a typical temperature evolution of average degree of ordering s(T) in a temperature loop using AH absence of external field. We started simulations at a relatively high temperature T = [T.sub.i] > [T.sub.IN], where nematic spins exhibit essentially isotropic ordering. Here [T.sub.IN] marks the isotropic-nematic phase transition temperature of a bulk sample. In our simulations we obtain [T.sub.IN] ~ 1.21 which is in line with existing simulation results using similar modelling [11]. Then we were gradually decreasing temperature until the temperature [T.sub.f] = 0.5 was reached, corresponding to conditions deep in the nematic phase. Afterwards, we were gradually increasing temperature till T = [T.sub.i]. The plot s(T) reveals that, for example, for {W, p} = {0.25, 0.1}, the system displays thermal hysteresis, suggesting first-order phase transition. Therefore, the effective disorder imposed by impurities is below the critical value (see Appendix B), above which the first-order phase transition between isotropic (or paranematic) and nematic phase is replaced by a gradual evolution of ordering on varying temperature. Furthermore, one sees that, for {W, p} = {1, 0.1}, the hysteresis is almost totally suppressed indicating gradual evolution of nematic ordering on varying temperature [17]. Therefore, for this choice of parameters, the effective disorder is supercritical.

Next, we consider impact of B on nematic ordering. We first demonstrate impact of B on range of ordering. In Figure 2, we plot G(r) dependencies using parameter set {p, W, T} = {0.1, 3, 0.5}, for (i) B = 0, (ii) B = 0.5, and (iii) B = 0.5 and both histories FQH and TQH. The parameter set {p, W, T} was chosen in such a way that SRO is realized for (i) B = 0. In the case (ii) the external field is strong enough to replace SRO with QLRO. For still stronger fields, see (iii); LRO is realized. The range of order was obtained via finite size analysis, using (5). For example, in Figure 3(a), we show the finite size analysis for TQH. For (i) B = 0 (Figure 3(a)), we obtain [gamma] = 1.54 [+ or -] 0.05 fingerprinting SRO. For (ii) B = 0.5, (Figure 3(b)) simulations yield [gamma] = 0.34 [+ or -] 0.05, suggesting QLRO. On the other hand for (iii) B = l.5, we get [gamma] = 0.00 [+ or -] 0.02, indicating LRO (Figure 3(c)). In Table 1 we summarize simulation results for [gamma] on varying B for extreme histories TQH and FQH.

Next, we consider impact of B on characteristic linear length [xi] of the system. The following behavior is expected. For relatively weak external field strengths, the impurity imposed disorder dominantly influences nematic structure. According to Larkin-Imry-Ma theorem one expects a domain pattern characterized by [[xi].sub.d] ~ [D.sup.-2/(4-d)], where D measures the disorder strength. Our mean-field estimate shown in Appendix B suggests D [varies] [square root of (p)] (see (B.6)). On the contrary, for relatively large external fields, one expects that the external field coherence length [[xi].sub.B] [varies] 1/B characterizing the impact of B in bulk LC system would be observed. Using our simulations we estimate [xi] by fitting G(r) dependences, calculated from fixed point configurations, with (4), where [xi] ~ 1/k. Characteristic dependencies [xi](B) are shown in Figure 4. One indeed observes tendency [xi] ~ [[xi].sub.B] [varies] 1/B ([xi] ~ [[xi].sub.d], roughly independent of B) for relatively strong (weak) external field B. We refer to regimes where [xi] ~ [[xi].sub.B] and [xi] ~ [[xi].sub.d] are external field dominated regime and disorder dominated regime, respectively. We estimate crossover values B = [B.sub.c], which roughly separate both regimes, via intersection of linearized [xi](B) dependencies obtained in limiting regimes B ~ 0 and B [much greater than] 10.

Different regimes are indicated in Figure 5 in the {B, T} plane for the parameter set {W, p} = {3, 0.1}. Here full symbols refer to FQH and empty symbols to TQH. Squares (triangles) indicate the line separating SRO-QLRO (QLRO-LRO) regimes. On the other hand spheres indicate crossover regime separating external field and disorder dominated regime. We determine the respective crossover fields by [B.sup.(QLRO).sub.c], [B.sub.c], and [B.sup.(LRO).sub.c], corresponding to SRO-QLRO transition, disorder-external field dominance crossover, and QLRO-LRO transition, respectively. With increasing temperature these lines are pushed towards larger values of B. Namely, due to increased thermal fluctuations, higher external fields are needed to reach regime exhibiting qualitatively stronger degree of ordering.

Finally, we analyzed external field controlled memory effects in systems. We performed the following B-loops. (i) Initial configurations were obtained for a given parameter set {p, W, T} via TQH for B = 0. (ii) Then we switched on the external field of strength B, for which we calculated a fixed point configuration. (iii) Afterwards we switched the field off and calculated the corresponding fixed point nematic structure. We refer to respective configurations as the (i) reference, (ii) field-on, and (iii) field-off structures, respectively.

In Figure 6, we plot the average degree of orientational order S obtained from [[Q.bar].sup.(g)] on increasing B for different anchoring strengths deep in the nematic phase for W = 3. For this anchoring strength the reference structures exhibit SRO, yielding S = 0 within numerical error. One sees that in field-on state the degree of ordering is substantially enhanced due to external field induced preference alignment. In field-off structures finite degree of ordering is preserved providing P < 0.2. We refer to the resulting ordering as remanent order [S.sub.r]. In Figure 6 one sees that [S.sub.r](B) displays two different regimes. In the weak field regime [S.sub.r](B) monotonously increases with B until the saturation regime is entered, in which [S.sub.r](B) variations are relatively small on varying B.

4. Conclusions

We study numerically memory effects using modified Random Anisotropy Nematic-type lattice model [16,17]. In our model impurities are randomly spatially distributed enforcing to nearby nematic spins quenched randomly chosen orientation. In simulations we varied concentration p of impurities, coupling strength W between nematic spins and impurities, temperature, and external field strength B. Furthermore, in most simulations we used two extreme histories of nematic spins, namely, field quenched history (FQH) and temperature quenched history (TQH). Using them, we gain information on memory capabilities of systems of interest. We calculated fixed point nematic configurations, revealing interplay between nematic interactions, local disordering, and global ordering fields. From fixed point configurations we calculated two-point orientational correlation function G(r) and average local and global nematic order parameter. From calculated G(r) profiles we determined characteristic length [xi] of the system. Furthermore, using finite size analysis, we determine range of orientational ordering for a given set of parameters.

Key results of our study are as follows. We determined regimes exhibiting LRO, QLRO, or SRO orders as a function of temperature and external ordering field B. The range of order was determined from values of [gamma] using finite size analysis. Note that the Lakrin-Imry-Ma [1, 6] theorem predicts SRO in absence of B. The theorem was derived for random-field type disorder. However, we used slightly different type of disorder, which is closer to random-bond type imposed frustrations. In simulations we obtained SRO only for strong enough disorder. Furthermore, our simulations confirm that range of order strongly depends on history of samples. In particular, we demonstrate that strong enough applied field B could change range of order. In the {E, T} plane we estimate critical line separating SRO-QLRO and QLRO-LRO for different histories of samples. Histories play significant role revealing glass-type properties of the systems. Furthermore, we determined impact of B on characteristic linear correlation length [xi] within samples. We distinguish between disorder dominated and external field dominated regime. In the first regime impact of B is negligible. On the other hand in the second regime [xi] is roughly proportional with 1/B. We roughly determined crossover region separating these regimes for different temperatures. We also study impact of temporarily present field B on remanent average nematic ordering [S.sub.r]. For this purpose we used parameter for which [S.sub.r] = 0 (within numerical order) in the reference state obtained in absence of B. We determined magnitude of [S.sub.r] as a function of B. Our simulations reveal that [S.sub.r] monotonously increases with B until saturation region is reached.

Displayed memory capabilities could be of use for various external field driven memory devices. Modelled systems mimic to some extent LCs confined to various porous matrices or mixtures of nematic LCs with appropriate nanoparticles (e.g., aerosils). In our future studies we intend to explore in more detail analogies between our model and random magnets and spin glasses.

Appendices

A. Brownian Dynamics-Type Numerical Simulation

We consider the following Hamiltonian:

[H.sup.*] = H + [summation over (i)] [[lambda].sub.i] ([[??].sub.i] x [[??].sub.i] - 1), (A.1)

where H is given by (1) and [[lambda].sub.i] are Lagrange multipliers taking into account the normalization constraint [absolute value of ([[??].sub.i])] = 1. Within a sweep the orientations of all nematic spins are updated using the following equations:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.2)

Here [[??].sup.(old).sub.i] and [[??].sup.(new).sub.i] represent orientations of nematic spins before and after the sweep, [DELTA]t measures the time step of each sweep, D is the representative LC rotational diffusion constant, [k.sub.B] is the Boltzmann constant,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.3)

is the so called residuum, and the quantity [??] is defined in terms of unit vectors [[??].sub.1] and [[??].sub.2] as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.4)

The first term on the right side of (A.2) corresponds to a mechanical elastic torque tending to rotate the nematic spin towards a local energy minimum. The second term [DELTA][[??].sup.(T).sub.i] mimics random thermal fluctuations.

In simulations we use the dimensionless time step [DELTA][t.sup.*] and dimensionless temperature [T.sup.*] defined as follows:

[DELTA][t.sup.*] = [DELTA]tD/J, [T.sup.*] = [k.sub.B]T/J. (A.5)

Note that when the residuum drops to zero (i.e., it is small enough with respect to equilibrium-like fluctuations), the system reaches a fixed point configuration reflecting either a metastable or globally stable nematic configuration.

To calculate a thermal fluctuation vector [DELTA][[??].sup.(T).sub.i], we first consider rotation of a nematic spin in its eigen frame, within which the local z-axis points along the spin direction. Therefore, only rotations of the spin about perpendicular local x- and y-axes are relevant and we set that both rotations are mutually independent. We express the corresponding rotational energy contribution as follows:

[DELTA][W.sup.(rot)] = [c.sub.1]J/2 ([[OMEGA].sub.x.sup.2] + [[OMEGA].sub.y.sup.2]). (A.6)

Here [[OMEGA].sub.x] and [[OMEGA].sub.y] are the corresponding rotation angles and [c.sub.1] is the appropriate dimensionless parameter which enables phenomenological tuning of thermal fluctuations. Assuming the canonical distribution, we express the infinitesimal probability for the range of angles within the phase-space element d[[OMEGA].sub.x]d[[OMEGA].sub.y] as follows:

dP = C exp [-[c.sub.1] J([[OMEGA].sub.x.sup.2] + [[OMEGA].sub.y.sup.2])/2[k.sub.B]T] d[[OMEGA].sub.x]d[[OMEGA].sub.y], (A.7)

where C is the normalization constant. By expressing (A.6) in the two-dimensional polar coordinate system it follows

dP = C exp [-[c.sup.1][[OMEGA].sup.2]/2[T.sup.*]] [OMEGA]d[OMEGA]d[psi]. (A.8)

Here [OMEGA] is the magnitude of the vector ([[OMEGA].sub.x], [[OMEGA].sub.y]) and the angle [PSI] defines its direction with respect to the local x-axis. Therefore, we mimic thermal fluctuations by the rotational kinetic free energy with two independent components of the angular velocity [[omega].sub.x] and [[omega].sub.y]. Within a sweep, each spin is rotated by the angle [OMEGA] about the axis that is perpendicular to the local z-axis and makes the angle [PSI} with the local x-axis. The distribution of [PSI] is uniform in the interval (0, 2[pi]) and the probability distribution function for Q is given by p([OMEGA]) = C exp[-[c.sub.1][[OMEGA].sup.2]/2[T.sup.*]][OMEGA]. After the thermal rotation of the spin in its own frame is calculated the rotated spin is transformed to the laboratory system using the following transformation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.9)

B. Phenomenological Estimates

In the appendix we present estimates on phase and structural LC behavior in presence of impurity spins, where we assume that the system exhibits a domain-type pattern characterized by a single characteristic length. Randomly distributed impurities are assumed to produce random field-type disturbances. Impurities are homogeneously distributed within the sample of volume V with the volume concentration p = [N.sub.im][v.sub.im]/V. Here [v.sub.im] determines the volume of an impurity and [N.sub.im] is the total impurity number inside V.

We estimate the nematic ordering at mesoscopic level by the uniaxial tensor order parameter [Q.bar] = S([??] [cross product] [??] - [I.bar]/3). The quantity [??] is the nematic director field and S is the uniaxial orientational order parameter. The corresponding total free energy is given by

F = [integral][integral][integral] ([f.sub.c] + [f.sub.e]) + [d.sup.3][??] + [integral][integral] [f.sub.i][d.sup.2][??]. (B.1)

The first integral is carried over the LC volume, and the second one over LC-impurity interfaces.

The condensation ([f.sub.c]), elastic ([f.sub.e]), and interfacial ([f.sub.i]) terms are expressed as follows:

[f.sub.c] = 3/2 a(T - [T.sup.*]) Tr [[Q.bar].sup.2] - 9/2 b Tr [[Q.bar].sup.3] + 9/4c [(Tr [[Q.bar].sup.2]).sup.2] (B.2a)

[f.sub.e] = [L.sub.n]/2 [[absolute value of ([nabla][Q.bar])].sup.2], (B.2b)

[f.sub.i] = -W/2 [[??].sub.i] x [Q.bar][[??].sub.i]. (B.2c)

The latter contribution is describing interaction between LC molecules and impurities. We took into account only the most essential contributions. Here a, b, c, and [T.sup.*] are material constants, [L.sub.n] is the representative nematic elastic constant, W is a positive anchoring strength, and [[??].sub.i] represents the locally impurity preferred orientation, [absolute value of ([[??].sub.i])] = 1. We assume that impurities are essentially homogeneously distributed with concentration p and that orientational probability distribution of [[??].sub.i] is spatially isotropic.

We adopt the Larkin-Imry-Ma picture in which the disorder breaks the system into domains of characteristic size [xi]. Our aim is to express it in terms of system parameters. For this purpose we considered the mean free energy of a single domain of volume [V.sub.d] ~ [[xi].sup.3], within which there are impurities possessing surface area aim. Within each domain the nematic director is regarded to be on average aligned along a single symmetry breaking direction. The corresponding free energy penalty within a domain volume is estimated by

[DELTA][F.sub.d] ~ (a(T - [T.sup.*])[s.sup.2] - b[s.sup.3] + c[s.sup.4] + [L.sub.n][S.sup.2]/2[[xi].sup.2])[V.sub.d]

-1/3 [N.sup.(d).sub.im] [W.sub.s] <[P.sub.2]> [a.sub.im]. (B.3)

Here [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is the mean value of the second Legendre polynomial of ([MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]) within the domain. The penultimate term evaluates approximately the effect of the gradient term. The final term in (B.3) comes from evaluating [f.sub.i] over the domain. The final term can be expressed in terms of [xi] as follows. We assume the isotropic orientational distribution of orientations [[??].sub.i]. In an infinitely large domain one expects <[P.sub.2]> = 0. In a finite volume the cancellation is only partial. According to the central limit theorem one expects <[P.sub.2]> ~ [([N.sup.(d).sub.im]).sup.1/3]. The number of impurities and the characteristic impurity separation [l.sub.im] are related by ~ [N.sup.(d).sub.im] ~ [([xi]/[l.sub.im]).sup.3], where [l.sub.im] ~ [([v.sub.im]/p).sup.1/3]. An expression for [xi] is obtained by balancing the [xi]-dependent elastic and interface contributions in [DELTA][F.sub.d]. It follows

[xi] ~ 9[L.sup.2.sub.n] [S.sup.2]/4p[W.sup.2] [v.sub.im]/[a.sup.2.sub.im]. (B.4)

On the other hand, the Larkin-Imry-Ma prediction for [xi] is given by

[xi] ~ [D.sup.-2/(4-d)], (B.5)

where D measures the disorder strength and d is the space dimensionality. For d = 3, it follows

D ~ 2[square root of (p)] W[a.sub.im]/3[L.sub.n]s[square root of ([v.sub.im])]. (B.6)

We next consider the phase behavior. We neglect temperature dependence of [xi] and introduce the following new-scaled variables: the scaled order parameter q = s/[s.sup.0], where [s.sub.0] = s(T = [T.sub.NI]; p = 0); the dimensionless temperature [tau] = (T - [T.sup.*])/([T.sub.NI] - [T.sup.*]), where TOT = T* + [b.sup.2]/(4ac), denotes the bulk isotropic-nematic phase transition; [[xi].sup.(0).sub.n] = [square root of ([L.sub.n]/(a([T.sub.NI] - [T.sup.*])))] is the nematic order parameter correlation length calculated at T = [T.sub.NI]; [d.sup.(0).sub.e] = [L.sub.n][s.sup.2.sub.0]/W[s.sub.0] is the nematic surface extrapolation length likewise calculated at T = [T.sub.NI]. Using this scaling we rewrite the free energy in non-dimensional form as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (B.7)

where [sigma] = ([square root of (p)]/3)([a.sub.im]/[square root of ([v.sub.im])])([[xi].sup.(0).sub.2]/[d.sup.(0).sub/e] [[xi].sup.3/2]) and [[tau].sub.eff] = [tau] + [[xi].sup.(0)2.sub.n]/2[[xi].sup.2] play roles of dimensionless field and effective dimensionless temperature, respectively. In the limit of strong enough anchoring, where [xi] ~ [l.sub.im] ~ [([v.sub.im]/p).sup.1/3], and for spherical impurities of radius r, one obtains

[sigma] = p[[xi].sup.(0)2.sub.n]/[d.sup.(0).sub.e]r. (B-8)

It is now possible to link this model to the theory of nematics in a constant external field. According to it the 1st order phase transition exists at critical temperature corresponding to [[tau].sub.eff] = 1 + [sigma] as long as [sigma] < [[sigma].sub.c] [equivalent to] 0.5. For [sigma] > [[sigma].sub.c] one obtains gradual evolution of ordering on varying temperature.

http://dx.doi.org/10.1155/2014/423693

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Amid Ranjkesh acknowledges the support of Slovenian Research Agency (ARRS) for young researcher grant.

References

[1] Y. Imry and S. Ma, "Random-field instability of the ordered state of continuous symmetry," Physical Review Letters, vol. 35, no. 21, pp. 1399-1401, 1975.

[2] D. E. Feldman, "Quasi-long-range order in nematics confined in random porous media," Physical Review Letters, vol. 85, p. 48-86, 2000.

[3] L. Radzihovsky and J. Toner, "Nematic-to-smectic-a transition in aerogel," Physical Review Letters, vol. 79, no. 21, pp. 4214-4217, 1997.

[4] M. Buscaglia, T Bellini, C. Chiccoli et al., "Memory effects in nematics with quenched disorder" Physical Review E, vol. 74, no. 1, Article ID 011706, 2006.

[5] M. Rotunno, M. Buscaglia, C. Chiccoli et al., "Nematics with quenched disorder: pinning out the origin of memory," Physical Review Letters, vol. 94, no. 9, Article ID 097802, 2005.

[6] A. I. Larkin, "Effect of inhomogeneities on the structure of the mixed state of superconductors," Soviet Physics--JETP, vol. 31, no. 4, p. 784, 1970.

[7] G. P. Crawford and S. Zumer, Liquid Crystals in Complex Geometries Formed by Polymer and Porous Networks, Taylor & Francis, London, UK, 1996.

[8] P. G. de Gennes and J. Prost, The Physics of Liquid Crystals, Oxford University Press, Oxford, UK, 1993.

[9] T. Bellini, M. Buscaglia, C. Chiccoli, F. Mantegazza, P. Pasini, and C. Zannoni, "Nematics with quenched disorder: how long will it take to heal?" Physical Review Letters, vol. 88, no. 24, Article ID 245506, 2002.

[10] T. Jin and D. Finotello, "Controlling disorder in liquid crystal aerosil dispersions," Physical Review E, vol. 69, no. 4, Article ID 041704, 12 pages, 2004.

[11] A. Ranjkesh, M. AmbroZic, S. Kralj, and T J. Sluckin, "Computational studies of history dependence in nematic liquid crystals in random environments," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 89, no. 2, Article ID 022504, 2014.

[12] V. Popa-Nita and S. Romano, "Nematic-smectic--a phase transition in porous media," Chemical Physics, vol. 264, no. 1, pp. 91-99, 2001.

[13] T Bellini, M. Buscaglia, C. Chiccoli, F. Mantegazza, P. Pasini, and C. Zannoni, "Nematics with quenched disorder: what is left when long range order is disrupted?" Physical Review Letters, vol. 85, no. 5, pp. 1008-1011, 2000.

[14] M. Cvetko, M. AmbroCic, and S. Kralj, "Competition between local disordering and global ordering fields in nematic liquid crystals," Beilstein Journal of Organic Chemistry, vol. 6, article 2, 2010.

[15] S. Relaix, R. L. Leheny, L. Reven, and M. Sutton, "Memory effect in composites of liquid crystal and silica aerosil," Physical Review E, vol. 84, no. 6, Article ID 061705, 2011.

[16] R. Harris, M. Plischke, and M. J. Zuckermann, "New model for amorphous magnetism," Physical Review Letters, vol. 31, no. 3, pp. 160-162, 1973.

[17] D. J. Cleaver, S. Kralj, T. J. Sluckin, and M. P. Allen, Liquid Crystals in Complex Geometries Formed by Polymer and Porous Networks, edited byG. P. Crawford and S. Zumer, Oxford University Press, London, UK, 1996.

[18] V Popa-Nita, "Statics and kinetics at the nematic-isotropic interface in porous media," European Physical Journal B, vol. 12, no. 1, pp. 83-90, 1999.

[19] P. A. Lebwohl and G. Lasher, "Nematic-liquid-crystal order--a Monte Carlo calculation," Physical Review A, vol. 6, no. 1, pp. 426-429, 1972.

[20] J. Chakrabarti, "Simulation evidence of critical behavior of isotropic-nematic phase transition in a porous medium," Physical Review Letters, vol. 81, p. 385, 1998.

[21] R. Fisch, "Quasi-long-range order in random-anisotropy Heisenberg models," Physical Review B, vol. 58, no. 9, pp. 5684-5691, 1998.

[22] D. P. Landau, "Monte Carlo studies of finite size effects at first and second order phase transitions," in Finite Size Scaling and Numerical Simulation of Statistical Systems, V. Privman, Ed., chapter 5, World Scientific, Singapore, 1990.

[23] D. R. Denholm, Simulation studies of random-anisotropy magnets [Ph.D. thesis], University of Southampton, Southampton, UK, 1995.

Amid Ranjkesh, (1) Milan Ambrozic, (1) Samo Kralj, (1,2) and T. J. Sluckin (3)

(1) Faculty of Natural Sciences and Mathematics, University of Maribor, Koroska 160, 2000 Maribor, Slovenia

(2) Condensed Matter Physics Department, Jozef Stefan Institute, Jamova 39,1000 Ljubljana, Slovenia

(3) Division of Mathematical Sciences, University of Southampton, Southampton SO171BJ, UK

Correspondence should be addressed to Amid Ranjkesh; amidranjkesh@gmail.com

Received 4 August 2014; Accepted 15 October 2014; Published 30 November 2014

Academic Editor: Ales Iglic

TABLE 1: Finite size analysis data obtained from (A.2). Values of [gamma] reveal range of order within nematic structures. We used different histories of samples (FQH and TQH) and different values of the external field B. Values [gamma] = 0, 0 < [gamma] < 1.5, and [gamma] = 1.5 reveal LRO, QLRO, and SRO, respectively. W = 3, p = 0.1, and T = 0.5. B 0 0.5 TQH [gamma] = 1.54 [+ or -] 0.05 [gamma] = 0.34 [+ or -] 0.05 FQH [gamma] = 1.51 [+ or -] 0.08 [gamma] = 0.23 [+ or -] 0.08 B 1 2 TQH [gamma] = 0.08 [+ or -] 0.07 [gamma] = 0.0 [+ or -] 0.02 FQH [gamma] = 0.03 [+ or -] 0.04 [gamma] = 0.0 [+ or -] 0.01

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Ranjkesh, Amid; Ambrozic, Milan; Kralj, Samo; Sluckin, T.J. |

Publication: | Advances in Condensed Matter Physics |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 5799 |

Previous Article: | Laser cleaning using Q-switched Nd:YAG laser of low carbon steel alloys. |

Next Article: | Effect of UV laser conditioning on the structure of KDP crystal. |

Topics: |