1. Introduction
It has been recently highlighted that heliconical cholesteric liquid crystals (ChOH, oblique helicoidal cholesterics) offer the opportunity to observe light-induced effects that are not possible in conventional planar cholesterics (CLC). It has been pointed out that the new phenomenology linked to the conical arrangement of the molecular director
originated from a bend elastic constant
lower than the twist elastic constant
and requires the application of a low-frequency (static) electric field
in a direction parallel to the helix axis and within a range of values between two critical fields
. If
has a conventional CLC structure with a helix axis rotated by 90° with respect to the conical configuration, for
a complete unwinding of the structure occurs [
1,
2,
3,
4,
5].
Unprecedented tuning of the Bragg resonance over the whole visible spectrum is possible by scanning a low-frequency electric field [
1,
2,
3] over fractions of V/mm, providing possible easy applications to tunable optical filtering [
6]. Additionally, light-induced optical reorientation has been demonstrated, leading to the tuning of the Bragg photonic band gap (PBG) and control of the light polarization state from linear to elliptical to circular in dependence on light intensity [
7,
8]. It has been already underlined that effective optical reorientation due to the torque action of the optical field
on the molecular director
is possible in ChOH because of the bend deformation of the conical structure, which is not present in conventional CLCs [
7], where optical reorientation would necessarily involve a twist deformation, requiring field values higher by a few orders of magnitude to occur.
The effects of such optical reorientation on the nonlinear propagation of a light-beam-inducing optical reorientation have also been theoretically investigated [
9]. In fact, light propagation through this structure may have a strong nonlinear character due to the bend deformation and the consequent optical torque originated by the interaction between the optical field and the molecular director
. This torque competes with the one originated by the static field; the result is an increase in the conical angle
and, as a consequence, of the helix pitch. The geometry of the interaction and competing actions of the static and optical fields is sketched in
Figure 1. This leads to a redshift of the Bragg resonance, producing a nonlinear transmission of the light beam strongly dependent on its wavelength as it approaches the PBG.
This is easily expected if we extend to ChOH the theory of light propagation along the helical axis of CLC using the following substitution in all the relevant mathematical expressions [
9]:
In this way, the standard expression for the reflectance of the resonant mode can be used [
10]:
where:
using the definitions given in [
9] here and in the following.
Here
is the resonant mode wavevector, which for a left-handed helix is:
Additionally, it is necessary to take into account the dependence of the conical angle and the pitch on the effective field:
This can be done analytically by considering a constant value of the intensity over the whole sample thickness. At a fixed value of the applied static field in the range of stability of the heliconical structure, the intensity dependence of light transmission spectra
can be plotted by using Equation (2) and shows a strong light-induced redshift of the whole spectrum [
9].
The constant intensity approximation (corresponding to a uniform sample layer) is certainly well fulfilled in a pump–probe configuration when the wavelength of the light beam inducing the optical reorientation is far from the Bragg resonance and the probe beam near the Bragg resonance is weak enough to induce a negligible interaction with the molecular director. On the other hand, in a single-beam configuration where optical reorientation is induced by the same beam that is traveling through the medium with a wavelength close to the Bragg resonance, this strong modulation of transmission may affect the light propagation in a more complex way if the change in intensity during propagation cannot be neglected. To take into account this effect, a more careful calculation has been recently carried out based on Ambartsumian’s layer addition modified method [
11]. This method allows calculating the components of fields reflected and transmitted by the sample and also the field inside the material, making it possible to study the features of the localization of light in the system. Under this approach, the sample is considered a stratified material, and the optical field is recalculated after traveling through each sublayer; in this way, the effects of the light field are determined more carefully. These calculations have shown that, over a typical intensity dependence on the material parameters, instabilities occur in the optical transmission as a result of the nonlinear light propagation.
We present here a detailed analysis of this phenomenon, pointing out that, over a threshold intensity, stable oscillations take place until at higher intensities a chaotic behavior is reached. We demonstrate the strong asymmetric behavior in the spectral regions around the Bragg peak, where an oscillatory behavior at longer wavelengths corresponds to a uniform transmission at shorter ones. The possibility of switching from stable to oscillatory behavior by the static electric field is also shown. The phenomenon originates from the light-induced shift of the Bragg resonance, which affects light transmission during propagation, as will be described below.
2. Method and Results
We start from the conventional dielectric permittivity and magnetic permeability tensors of a ChOH:
We remind that
,
, where
and
are the dielectric permittivity parallel and perpendicular to the director in the nematic phase. Dividing the sample thickness into a large number of sublayers small enough to consider as constant the optical parameters of each sublayer, the problem consists of the solution of the following system of matrix difference equations:
with initial conditions
and
[
12,
13].
Here
, and
are the reflectance and transmittance matrices of the system with
m and (
m − 1) sublayers, respectively;
are the reflectance and transmittance matrices of the
m-th sublayer;
is the zero matrix;
is the unit matrix; the respective matrices for the reverse light propagation are denoted by a tilde. The method allows us to find the reflection coefficient
, the transmission coefficient
, and light intensity of the wave traveling in the ChOH sample
. We consider the eigenmode with circular polarization with the same handedness as the ChOH helix, which is suffering Bragg diffraction. As usual,
,
, and
are the fields of incident, reflected, and transmitted waves, respectively, and
is the intensity of incident light. Taking
, we operate with the normalized intensity
considering in the first step (
j = 1) of the calculation:
With these parameters, we calculated the reflection, transmission coefficients as well as field intensity in each location of the liquid crystal layer following the method presented in [
12].
In this way, we carry out the calculation that allows obtaining . In the second step (j = 2), we have the new effective field allowing us to take into account the distribution of the angle and the helix pitch along the sample as determined by the local effective field. Then we have and . Again we get the reflection and transmission coefficients as well as the intensity for this new parameter, and so on by increasing the number of steps: j = 3, 4, 5, ….
The material parameters chosen for the calculation are the ones reported in Ref. [
11]:
,
,
, and
; the unperturbated helix pitch is
= 1400 nm; the layer thickness is
d = 20 µm with critical fields
and
; and we take
, where
is the refractive index of the ChOH layer. The static field applied along the helix direction is
. In our calculations all the above-mentioned sublayers have the same thickness:
d1 =
d2 =
d3 = … =
dL = 8 nm.
We remark that computer programs were compiled in Visual Basic (based on analytical expressions for the reflection and transmission matrices for an anisotropic layer [
13] and the layer addition method [
12]) and debugged by the authors, and the graphs were built using the Excel and MATLAB programs.
The basic results of our calculations are shown in
Figure 2. Spectra of the light transmittivity of the ChOH sample are reported vs. increasing values of the step index
for increasing values of the normalized beam intensity
. Since the increase in the step index corresponds to the adjustment of the material parameters to the updated value of the local light field, it is strictly correlated to the response time of the medium
; therefore, we should consider
proportional to the elapsed time. In this way, these plots give us the evolution in time of the transmission spectra. First of all, we remark that at
the transmission spectrum has the well-known form shown in
Figure 2a and the change of
has no influence on this spectrum, for this reason only the spectrum at
is reported. In
Figure 2b, data for
to
are displayed. At
, as shown in [
11], the change of transmission does not exceed 0.01; therefore, the change of
has no noticeable effect on the evolution of the transmission spectrum (see
Figure 1b). We observe that already at
oscillations in the transmission occur at the red side edge of the Bragg resonance visible as the alternate yellow and red colors highlighting a regular oscillation of transmittivity form about 0.6 to about 0.8. Increasing the intensity of the oscillations corresponds to higher modulation of the transmittivity (alternate green and red colors) and occurs for a wider range of wavelengths. They also become evident in the blue edge of the Bragg resonance, even if showing a weaker modulation and a narrower range of wavelengths. At
we observe a remarkable narrowing of the PBG (deep blue stripe), which corresponds to irregular variable transmission on the red side edge of it.
Figure 2c shows the spectra for further increase in the intensity from
to
. We also observe a broadening of the wavelength range of instabilities occurring in the blue side edge of the spectrum. A large increase in the intensity shows the loss of any regular behavior with the disappearance of the PBG.
In order to better highlight the feature of the observed phenomenon, we plot the transmission vs.
j for two wavelengths
λ1 = 494.4 nm and
λ2 = 510.4 nm symmetric with respect to the Bragg peak
λB = 502.4 nm on the blue side and on the red side of the edge of the resonance. These curves are reported in
Figure 3 for increasing values of the normalized intensity
N. We observe that at
N = 150, a stable value of the transmitted intensity is quickly achieved for both wavelengths (
Figure 3a). At
N = 500, a stable value is reached after a few damped oscillations at
λ1 while stable oscillations are present at
λ2 (
Figure 3b).
A further increase in the intensity (
N = 650) also gives rise to stable oscillations (usually of lower amplitude) at
λ1 when at
λ2 one already observes irregular oscillations (chaotic behavior) (
Figure 3c) that are present at both wavelengths for higher intensities (
Figure 3d, corresponding to
N = 700).
We also pointed out the strong difference between the oscillatory behavior and the unstable variable values shown at high intensities by calculating the Fourier transforms of these curves. We report in
Figure 4 two examples related to the behavior shown at
for
and
. In these plots, we see a clear peak in the FT spectrum, showing that oscillation is occurring with a narrow frequency range for the first case, while at higher intensities no clear oscillation frequency can be identified. This result highlights the occurrence of two different regimes in the nonlinear propagation of the light at this wavelength.
3. Discussion
The amazing self-induced oscillatory behavior occurs at intensities easily achievable by a c.w. laser with moderate power, in fact,
N = 500 (just above the threshold intensity for
λ2 = 510.4 nm) corresponds to
I = 1.9
10
5 W/cm
2. A simple model allows an understanding of the onset of the oscillations; it is illustrated in
Figure 5 where the transmission spectrum at
j = 1 is plotted. It has already been demonstrated that optical reorientation leads to a redshift of the Bragg resonance, as summarized by the above-reported Equations (5) and (6) fulfilled by the experimental data reported in [
7,
8]. Then, by considering a narrow range of wavelengths on the red side edge of the PGB (indicated by the red arrow in the figure) we notice that for these wavelengths the redshift of the Bragg resonance means a decrease in transmittivity (red curve). This results in a decrease in light intensity. The consequence is an increase in the effective static field
Eeff with consequent blueshift of the resonance leading in turn to an increase in transmittivity and light intensity inducing again a redshift of the resonance able of restarting the process (blue curve). This model also explains the strong asymmetry reported in
Figure 3, which leads to oscillatory behavior at higher intensities for wavelengths on the blue side edge of the PBG. In fact, in this case, the same process occurs around the secondary maxima of the resonant spectrum and not next to the PBG; thus the oscillatory behavior is observed at a higher intensity with a lower amplitude of modulation.
It is quite important that at the fixed impinging intensity the effect expected at a given wavelength it is determined by its distance δλ from λB and occurs at a lower intensity on the red side edge of the resonance. As a matter of fact, δλ depends on the static field ES applied to the ChOH sample, which determines the location of λB on the wavelength axis. This fact offers the opportunity of using the static field ES as a parameter to drive the transmittivity behavior for a light beam at a specific wavelength. It means that the optical transmission at a fixed wavelength and light intensity can present a stable value, a self-oscillatory behavior, or a chaotic one depending on the applied electric field.
As an example, we show in
Figure 6 the calculated transmission vs.
j at
λ1 = 494.4 nm at different intensities for two different values of the applied static field. We observe (
Figure 5a) that with
N = 150 we get a stable value for the light transmission with an applied field of
ES = 2.03 V/mm, while at
ES = 2.09 V/mm we observe the onset of stable oscillations. Alternatively, at
N = 650 (
Figure 5b) we have stable oscillations at
ES = 2.03 V/mm and chaotic behavior with
ES = 2.09 V/mm. This result shows that we have a system that can be driven to oscillatory or chaotic behavior by a static electric field and light intensity that is able to switch the system between different states. These properties need further theoretical and experimental investigation.
It is important to remark that all data reported in
Figure 2,
Figure 3 and
Figure 4 show the actual time dependence of the transmitted intensity because the step index
j can be considered as the elapsed time in units of the response time
τon required by the material to adapt its parameters to the new field values:
j~
t/τon. In this way, we are also able to foresee the time range of the observed oscillations. In fact, they are related to small changes in the tilt angle
θ determining the bend distortion of the helical structure. The maximum value of
τon is the relaxation time
τoff necessary for the system to come back to equilibrium after switching off the exciting optical field, and in general, is shorter and shorter for increasing intensities. Therefore, we can evaluate the order of magnitude of the response time from the expression used when considering a bend distortion in homeotropic nematics [
14]:
and find that with typical values for the bend elastic constant
K3 and the viscosity
for a sample thickness of 20 μm gives
τoff~40
ms. However, we must take into account that ChOH exists only when
K3 <
K2, in fact, its value has been recently measured [
15] to be 2–3 times lower than in nematics; therefore it is possible a longer response if a similar value of the viscosity is used. For this reason, we can expect these oscillations to occur in the tenth of a millisecond range. A more detailed analysis is underway on this feature to compare the variation of the oscillation frequency at different intensities and experimental data are necessary to have a deeper insight into this matter.
4. Conclusions
Here we have reported a detailed study of the nonlinear light propagation of a beam traveling along the helix direction of a heliconical cholesteric liquid crystal, using Ambartsumian’s layer addition modified method. The nonlinear light propagation is considered an effect of the optical reorientation induced by the light beam on the molecular director and the time evolution of the whole transmission spectrum is analyzed at increasing values of light intensity. We have reported how this spectrum is modified by increasing the light intensity, showing the onset of damped oscillations first followed by stable oscillations and by irregular (chaotic) behavior at higher intensities. We have demonstrated the strong asymmetry of the self-induced effects. In fact, they take place at lower intensities on the red side edge of the Bragg resonance while stable values of transmittivity are still observed on the blue side edge. A further increase in the intensity leads to oscillations of lower amplitude on the blue side while an irregular behavior of the transmission is achieved on the red side. At higher intensities, the Bragg resonance disappears and transmission becomes unstable for any light wavelength. A simple phenomenological model is proposed that explains both the onset of the oscillatory behavior and the wavelength asymmetry on the opposite sides of the Bragg resonance. In fact, since optical reorientation leads to the redshift of the Bragg peak, for wavelengths near the long-wavelength boundary of the PGB the redshift means a decrease in transmittance. This leads to a decrease in the intensity of the light. A consequence is an increase in the effective static field acting on the structure followed by a blueshift in resonance, which in turn leads to an increase in transmittance and light intensity, again causing a redshift of the resonance, which can restart the process. This model also explains the strong asymmetry, which leads to oscillatory behavior at higher intensities for wavelengths at the blue edge of the resonance. In fact, here the same process occurs around the secondary maxima of the resonance spectrum, so the oscillatory behavior occurs at higher intensities with lower amplitude.
We also point out that the static electric field is a driving parameter to switch from stable to oscillatory to irregular behavior in the transmittivity at a given light wavelength.
It is remarkable that we have a system with tunable “chaotic” behavior driven by light intensity and external static electric field that has already been shown in a previous work [
11] to correspond to a disorder of the conical structure changing both the conical angle and pitch through the sample thickness. This is potentially interesting for application development as it has been underlined in Ref. [
16] where it has been shown that a completely disordered and scattering optical medium can sometimes outperform the most carefully designed and engineered photonic device.
Finally, it is important to underline the difference between the effect of the light-induced shift of the Bragg resonance discussed in this paper as well in the previous papers related to an optical reorientation in ChOH and the changes in the resonance obtained in conventional CLC, reported in a number of articles in the past [
17,
18,
19]. Here we have the direct interaction between the optical field and the director leading to reorientation and consequent change of the helix pitch, while the changes of the Bragg resonance observed in conventional cholesterics were obtained through trans-cis isomerization, leading to a conformational change of the molecules.