Abstract

We calculate the trajectory of a monochromatic optical beam propagating in a planar-homeotropic hybrid nematic crystal cell submitted to weak anchoring conditions. We apply a uniform electric field perpendicular to the cell to control the trajectories for various values of the anchoring elastic energy. We have found that the anchoring energy has a strong influence on the ray penetration length and trajectory. Our calculations are consistent with a previously found Frederick’s type transition only present for weak anchoring in which electric fields above an anchoring-energy dependent critical field align completely the director field.

1. Introduction

The total internal reflection in hybrid nematic cell has been studied a time ago in a series of papers [13]. The authors report theory and experiments for an optical method based on observing the interference fringes of a beam reflected by a nematic cell that are used to know the local orientation of the molecular director inside the cell. Furthermore, they show that using this method it is possible to get information about the nematic anchoring conditions. In the present report we consider an analogous system but we use a different theoretical approach. Based on geometrical optics arguments, we derive the ray trajectories inside the cell since we are interested in using the total internal reflection as a mechanism for beam steering in optical devices.

Electrooptical devices such as lenses with variable focal distance, dynamic diffraction gratings, and tunable prisms are extremely useful in many applications, for instance, those concerning telecommunication, machine vision, displays, data storage, measurement equipment, and so on. Diffractive optical elements have many advantages over refractive elements [4, 5]; they can rapidly steer, stabilize, and increase reality by reducing costs. They are also small enough to be portable and can be adapted for micron-scaled optical systems. A nematic liquid crystal has a molecular architecture which makes it a medium with a direction-dependent refractive index [5]. By imposing an external electric field to the liquid crystal, it induces variations in the refractive index distribution and in the phase of light waves traveling through it which produces a modification on the direction of propagation of light through the medium. Liquid crystal devices can be much smaller and have less weight than conventional glass and plastic analogues. In addition, they can be integrated into other optical components and compact systems [5]. Different types of optical elements based on liquid crystals technology have been implemented [420]. In particular, liquid crystal technology has been applied to construct diffraction gratings that can electrically modulate the diffraction efficiency [21]. It has been shown [22] that, for a monochromatic beam that enters obliquely onto a nematic hybrid cell, the optical path of the beam can be several times larger than the cell’s thickness. Therefore, it is possible to use it as a dispersive media similar to a glass prism. As explained in [22], this device has the advantage that all the emerging beams of different wavelength will emerge parallel to each other although the emerging position is different. In other words, rays of different wavelength travel different distances but the emerging angles are equal due to ray’s momentum conservation. In the cited work [22], hard anchoring conditions for the nematic on the cell plates has been assumed; this means that the orientation of the molecules at the plates is not affected by the applied electric field. This assumption is valid only when the external field energy is weaker than the surface anchoring energy.

An anchoring energy is necessary to describe arbitrary anchoring conditions between the liquid crystal and the solid plates [23] which can be interpreted as an anisotropic part of the interface energy. The anchoring energy and its influence on the configuration provide important information of the liquid crystal surface physics. The physical reasons for the anchoring phenomenon and its microscopic mechanisms have been subject of interest [24]. Some years ago, Yang et al. [25] proposed that the interface energy should be understood as the total sum of energy potentials between the molecules of a liquid crystal and the substrate surface. They have found an expression for the anchoring energy, which has two terms, the first term is the same as the Rapini-Papoular expression; the second is related to the normal of the interface and results from the biaxial property of a NLC induced by the interface. However, the anchoring energy may originate from other physical reasons, such as adsorptive ions [2630]. General expressions for the nematic configuration considering all these phenomena have been recently obtained for a nematic slab under arbitrary boundary conditions [31]. In [22] strong anchoring conditions were considered; however, experimentally there is a wide range of values for this energy [32, 33], ranging from  J/m2 to  J/m2. For this reason, it is important to explore what could be the implications of using weak anchoring energies.

In this work, we calculate the orientational state of a liquid crystal submitted to weak anchoring conditions and to the action of a low frequency electric field. Then, we analyze the effect of changing the surface anchoring energy on the penetration length, the path, and the range of a beam traveling inside a liquid crystal cell and exhibit its notorious influence. This element can be used as the base for a beam steering or a multiplexor device.

2. Equilibrium Configuration

The system under study consists of a pure thermotropic nematic confined between two parallel substrates with refraction indexes and , respectively, as depicted in Figure 1. The cell thickness, , measured along the -axis, is small compared to the dimension, , of the cell plates. The director’s initial configuration is spatially homogeneous along the plane - and varies with so that at the boundaries the director where is the orientational angle defined with respect to the -axis.

In order to obtain the equilibrium configuration of the director, one has to consider the free energy of the nematic

The elastic part of the free energy is given by the Frank–Oseen expression [34]where the elastic moduli , , and correspond to transverse bending splay, twist, and longitudinal bending deformations, respectively. Here is the surface elastic constant. The term associated with this constant does not enter the Euler–Lagrange variational derivative for the bulk because it can be transformed in a surface term by using Gauss theorem. However, it can contribute to the energy and influence the equilibrium director through boundary conditions. Nevertheless for this particular planar geometry, this term is null as it is proved in [34].

is a parameter having units of energy per area and is associated with the interaction between the nematic molecules and the confining surfaces. The first integral involved is on the volume of the LC contained between the planar cell and the second and third ones should be performed over the upper and lower plates and are denoted by and , respectively.

A low frequency uniform electric field parallel to the -axis is applied. Thus, the electric contribution to the energy is

Here, is the electric displacement vector and is the component of the dielectric tensor.

The total free energy of the LC is obtained by expressing the integrals of (3) and (4) in Cartesian coordinates to obtain the free energy per unit length, and introducing the variable ,wherewhere , , and and are the low frequency dielectric constants perpendicular and parallel to the director and is the dielectric anisotropy.

To find the equilibrium condition we need to calculate the variation of the energy , as shown in [34]Note that is not an explicit function of ; however, it is a function of and . To proceed, we interchange the order of the variation and derivatives and use the chain rule for partial derivatives,After substitution of these results in (8) we obtainUpon application of the fundamental theorem of calculus we getFinally, to obtain the extrema of this functional we impose the condition that has to be valid for allNotice that the first condition stated in last expression is the usual Euler–Lagrange equation valid even for hard anchoring conditions whereas the two remaining expressions are mixed boundary differential equations for the orientation angle to be fulfilled at the frontiers of the region involved.

Substitution of (6) of the bulk energy in the first condition stated in the latter expression gives

After inserting (6) and (7) in the second and third conditions of (12) we obtain the following explicit expressions for the lower plate:and for the upper platewhere .

To obtain the distribution of the nematic orientation we solve (13) subjected to the boundary conditions (14) and (15). We find numerically the solutions for the nematic phase 5CB at 25.1°C, for which  N and  N. To do this, we use a Runge Kutta algorithm from the low plate by using as starting values for and , those satisfying (14) with , a given trial value. Thus, we calculate and at the top plate and check whether these values satisfy (15). If it is not the case we use another trial value for until we reach the correct configuration [35].

In Figure 2 we have depicted the orientation versus the dimensionless distance parameterized by the values of . The different panels corresponds to various values of the externally applied field . Figure 2(a) corresponds to and shows that for the smallest values of the nematic orientation is practically uniform. For larger values of the spatial dependence keeps almost linear while the slope increases until the configuration converges to the straight line connecting the angles and for the largest value of for which we recover hard anchoring conditions. The same trend is observed in Figure 2(b), however, the presence of the field makes that for the smallest values of and , the molecules of the nematic align completely parallel to the direction of the field, and thus in these cases. Figures 2(c) and 2(d) for and exhibit similar trends. Notice that the stronger the field, the larger the number of cases that align perfectly with the field. Note that the curves in Figure 2(d) corresponding to a larger field have a much more pronounced curvature as a consequence of the trend of the molecules to try to align with the field.

In calculating these curves, we have found that when the electric field surpasses certain threshold for a given value of , the director aligns completely with the external field and thus conforms a uniform configuration. This result is in agreement with a previously found Frederick’s transition obtained by Barbero and colleagues [34, 36]. For completeness we have plotted in Figure 3 this threshold versus the parameter . Above this curve, the director is aligned parallel to the field. A more extensive study of this transitional configuration has been done elsewhere [31] under the presence of some other effects as the presence of adsorptive ions. Indeed, our results are consistent with those found there.

3. Ray Propagation

An obliquely incident light beam with -polarization (-wave), that is, with the electric field contained in the incidence plane -, impinges the nematic with an incident angle as shown in Figure 1. The intensity of the beam is low enough so that it does not distort the nematic’s configuration. The dynamics of this optical field is described by Maxwell’s equations which contain the dielectric tensor , corresponding to an uniaxial medium which has the general form where and are the dielectric constants perpendicular and parallel to the director and is the dielectric anisotropy.

The usual method to solve Maxwell’s equations has been carried out in detail for a hybrid cell similar to the one considered here [37], and it is found that there is a regime for the incidence angle where the ray trajectory presents a caustic, that is, a geometrical place where the beam bends and remains inside the cell until it returns back towards the incidence substrate (see Figure 1). This trajectory is given by [37]In this equation the dimensionless variable has been introduced, is the ray component in the direction, and is a constant that is determined by the incident point of the beam on the cell; that is, it establishes an initial condition. The sign in (17) corresponds to a ray traveling with in the ± direction, that is, going from to and from to , respectively (see Figure 1).

As explained in [38], there are two regimes for . The first one corresponds to , with being a critical angle, where all the rays always reach the top substrate and part of the ray is transmitted to the top plate. On the other hand, the second regime corresponds to , namely, when the beam does not get the top substrate and it is reflected back to the interior of the cell as depicted in Figure 1. Besides , there is a second critical angle, , for which the beam no longer enters the liquid crystal cell and at which it is reflected back to the lower substrate. Here we will consider only angles for which the ray penetrates the cell and is reflected back.

The director’s angle at the returning point, , is given byfrom which the critical angles, and , can be obtained by substituting and , respectively. It is worth mentioning that since this formula was derived exclusively from Maxwell’s equations, it is independent of the type of anchoring assumed for the orientation at the boundaries.

Figure 4 exhibits the caustic position or penetration length as a function of the incidence angle parameterized with the surface elastic energy parameter . Similarly as in Figure 3, the different panels correspond to different field intensities. Figure 4(a) shows that in the absence of field () the range of incidence angles at which goes from zero to one is narrower for smaller values of . Figure 4(b) shows that the presence of the field widens the mentioned interval of incidence angles but overcoat for small values of . In Figures 4(c) and 4(d) the curves for different values of almost merge in one, demonstrating that the surface elastic energy is not playing an important role for large fields.

Figure 5 displays the beam trajectories versus the surface elastic energy parameter where Figures 5(a) and 5(b) correspond to and , respectively. Figure 5(a) exhibits how the trajectory range changes dramatically on the value of when there is no applied field. That is, for the smallest value of considered, the beam travels 35 times the thickness of the slab whereas for the largest it only travels five times the mentioned thickness. Conversely, Figure 5(b) shows that the trajectory range varies less than one cell thickness on the value of when the external electric field is strong. This implies that whether it is needed to design a device for performing a beam steering with large transverse displacements, it is necessary to use a coating for the cell plates which causes a weaker anchoring condition.

4. Conclusions

In this paper we have generalized our previous model for a monochromatic beam impinging obliquely onto a nematic hybrid cell, to take into account arbitrary anchoring conditions. We have found that the ray penetration range and trajectory are strongly affected by the anchoring energy. For example, we have shown that under low anchoring surface energies (small ) the trajectory range could be roughly an order of magnitude larger than those ranges associated with hard anchoring conditions, which already were several times the cell’s thickness. Thus, whether an electrooptical instrument is required for performing a beam steering with large transverse range, it would be essential to apply a coating substance with less stickiness at the cell plates in order to weaken the anchoring energy. Since in experiments there is a wide range of values for this energy [32, 33], ranging from  J/m2 to  J/m2, our results can be applied to a wide range of systems. Our results consistently confirms previous results [34, 36] for a Fredericks’ type transition consisting in a complete alignment of the director field for strong enough applied electric field which is only possible due to the use of weak anchoring. It is worth mentioning that in using this device all the emerging beams whose position can be controlled by the electric field will emerge parallel to each other so this effect is particularly useful for applications requiring to move parallel to a light beam.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper. The funds mentioned in Acknowledgments do not lead to any conflicts of interest regarding the publication of this manuscript.

Acknowledgments

This work was supported in part by DGAPA-UNAM under Grant nos. IN-110516 (Carlos I. Mendoza) and IN110012.3 (Juan Adrian Reyes).