Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Modulation of coherence and polarization using liquid crystal spatial light modulators

Open Access Open Access

Abstract

We propose a method for modulation of coherence and polarization of electromagnetic fields, employing two crossed zero-twisted nematic liquid crystal spatial light modulators. In contrast to a similar technique analyzed by Shirai and Wolf [J. Opt. Soc. Am A, 21, 1907, (2004)] our method provides a wide range simultaneous modulation of coherence and polarization. The dependence of the obtained results on different definitions of electromagnetic coherence is considered.

©2009 Optical Society of America

1. Introduction

When using laser illumination, it is frequently needed to destroy completely or partially the spatial coherence of an optical field. Such a situation arises, for example, in coherent imagery when it is necessary to eliminate undesirable speckle patterns [1] or in the problem of generating propagation-invariant fields when it is necessary to create a secondary source with a special structure of spatial coherence [2, 3]. The most efficient method to produce the desired changes of coherence is based on the use of a computer controlled liquid-crystal (LC) spatial light modulator (SLM) [4, 5]. In Ref. 5 it was shown that using a zero-twisted nematic LC SLM, under certain conditions one can obtain not only the changes of coherence but also the changes of polarization of the incident optical field. This technique can be efficiently employed to modulate coherence or polarization each taken separately. However it can not provide a wide range simultaneous modulation of both. To solve this problem, in the present paper we propose to use two crossed zero-twisted nematic LC SLMs with the different control signals. We also discuss the dependence of the obtained results on different definitions of the electromagnetic coherence. To ensure a thorough coverage of the total scope of the subject, we start the paper with a brief summary of the main results reported in Ref. 5.

2. Statistical properties of the electromagnetic field passed through a zero-twisted nematic LC SLM with random transmittance

We consider a stochastic stationary electromagnetic field propagating within a narrow solid angle around the z axis. According to the coherence theory in the space-frequency domain [6] such a field can be represented, at any typical point x = (x, y) in some plane z = const and at any frequency ν, by the statistical ensemble

{E(x,v)}={Ex(x,v)Ey(x,v)},

where Ex and Ey are the Cartesian components of the electric field vector E. For the sake of simplicity from now on we will omit the explicit dependence of the considered quantities on frequency ν. The second order correlation properties of the electromagnetic field can be completely characterized by the so-called cross-spectral density matrix [7, 8]

W(x1,x2)=[Ex*(x1)Ex(x2)Ex*(x1)Ey(x2)Ey*(x1)Ex(x2)Ey*(x2)Ey(x2)],

where the asterisk and the angle brackets denote the complex conjugation and the ensemble average, respectively. According to Wolf [7], the degree of coherence and the degree of polarization of the electromagnetic field are defined by the formulas

η(x1,x2)=TrW(x1,x2)[TrW(x1,x2)TrW(x1,x2)]1/2,
P(x)=(14DetW(x,x)[TrW(x,x)]2)1/2,

respectively, where Tr stands for the trace and Det denotes the determinant.

Now we consider the propagation of the electromagnetic field through a polarization-dependent device with a random transmittance given by the Jones matrix T(x). Let W in(x 1, x 2 ) be the cross-spectral density matrix of the incident field. Then the cross-spectral density matrix of the transmitted field is given by the expression

Wout(x1,x2)=T(x1)Win(x1,x2)T(x2),

where the dagger denotes the Hermitian conjugation. Hence, taking into account the definitions given by Eqs. (3) and (4), one can conclude that, in the general case, the degree of coherence and the degree of polarization change on propagation through a polarization-dependent device with a random transmittance. To evaluate these changes, one has to specify the form of the incident electromagnetic field and the type of the polarization-dependent device. Following Ref. 5, as an incident field we will consider the linearly polarized spatially coherent Gaussian field given by the cross-spectral density

Win(x1,x2)=E02exp(x12+x224ε2)[cos2θcosθsinθsinθcosθsin2θ],

where E 0 and ε are positive constants, and θ is the angle that the direction of polarization makes with the x axis. It can be readily verified that for such a field η in(x 1, x 2) = 1 and P in(x) = 1. As a polarization-dependent device in Ref. 5 it was employed a zero-twisted (or parallel-aligned) nematic LC SLM with the extraordinary axis aligned along the y direction. The transmittance of such a modulator is given by the matrix

T1(x)=[100exp[iβ1(x)]],

where

β1=πdλ(neno)

is the so-called birefringence, no and ne are ordinary and extraordinary indices of refraction, respectively, λ is the wavelength, and d is the thickness of the crystal. It is assumed that the birefringence has the form

β1(x)=φ0+φ(x),

where φ 0 is a constant which is much lager than π/2 and φ(x) is a zero mean random variable which is characterized by the Gaussian probability distribution density

p[φ(x)]=12πσφexp(φ2(x)2σφ2),

with variance 〈φ 2(x)〉 = σφ 2 and cross-correlation defined for two different points as

φ(x1)φ(x2)=σφ2exp(ξ22γ2),

where ξ = ∣x 1-x 2∣ and γ is a positive constant characterizing the correlation width of function φ(x).

By substituting Eqs. (6) – (9) into Eq. (5) and making use of Eqs. (A3) and (A8) from Appendix, one obtains

Wout(x1,x2)=E02exp(x12+x224ε2)
×[cos2θexp(iφ0)exp(σφ22)sinθcosθexp(iφ0)exp(σφ22)sinθcosθexp{σφ2[1exp(ξ22γ2)]}sin2θ].

Then, substituting Eq. (12) into Eqs. (3) and (4), one finds

ηout(ξ)=cos2θ+exp{σφ2[1exp(ξ22γ2)]}sin2θ,
Pout(x)={1[1exp(σφ2)]sin22θ}1/2.

To provide a physical insight into these results, let us discuss two special cases considered in Ref. 5, namely θ = 0 and θ = π/2. In the first case η out= 1 and P out = 1, i.e. no modulation of coherence or polarization occur. In the second case

ηout(ξ)=exp{σφ2[1exp(ξ22γ2)]},

and P out = 1. This time, the LC SLM provides the modulation of coherence without changing the polarization of the incident field. In the intermediate case, i.e., when 0<θ <π/2, the spectral degree of coherence and the spectral degree of polarization of the field transmitted through the zero-twisted LC SLM decrease simultaneously in accordance with Eqs. (13) and (14). In this case the particular choice θ = π/4 provides the minimum possible value of the spectral degree of polarization, namely exp(-σφ 2 / 2), and the degree of coherence

ηout(ξ)=12+12exp{σφ2[1exp(ξ22γ2)]}.

The degree of coherence given by Eq. (16) is illustrated in Fig. 1 for several values of parameter γ. The width of function ηout (ξ, which according to Ref. 1 is evaluated by

Δηout=0ηout2(ξ),

may be associated with the so-called transverse coherence length of the field [6]. As it is obvious from Fig. 1, in the considered case ∆ηout= ∞ for any value of parameter γ, i.e. the transmitted field is practically coherent. Therefore, we can conclude that the technique analyzed in Ref. 5 does not provide the appropriate simultaneous modulation of coherence and polarization.

 figure: Fig. 1.

Fig. 1. Degree of coherence given by Eq. (16), σφ = 2 and γ = 1,2,3.

Download Full Size | PDF

3. Statistical properties of the electromagnetic field passed through two crossed zero-twisted nematic LC SLMs with random transmittance

To attain a wide range simultaneous modulation of coherence and polarization, we will place at the output of the first LC SLM, considered in the previous section, a second SLM (identical to the first one) with the extraordinary axis aligned along the x direction. The transmission matrix of the second SLM is

T2(x)=[exp[iβ2(x)]001],

where the birefringence β2(x) has the form

β2(x)=φ0φ(x),

with φ0 and φ(x) of the same meaning as stated in the context of Eq. (9). The transmittance of the system formed by the crossed LC SLMs is given by

T(x)=T2(x)T1(x)=exp(iφ0)[exp[(x)]00exp[(x)]].

At first sight it seems that the above system of two crossed SLMs can be replaced by a sole 90°-twisted LC SLM operating in the adiabatic limit [9]. But it is not the case since the Jones matrix for the twisted LC SLM differs from the one in Eq. (20) by the presence of a common exponential factor depending on birefringence too (see, e.g., Ref. 10).

By substituting Eqs. (6) and (20) into Eq. (5) and making use of formula (A8) from Appendix, we obtain

Wout(x1,x2)=E02exp(x12+x224ε2)exp(σφ2)
×[exp[σφ2exp(ξ22γ2)]sin2θexp[σφ2exp(ξ22γ2)]cosθsinθexp[σφ2exp(ξ22γ2)]sinθcosθexp[σφ2exp(ξ22γ2)]cos2θ].

Then, substituting Eq. (21) into Eqs. (3) and (4), we find

ηout(ξ)=exp{σφ2[1exp(ξ22γ2)]},
Pout(x)={1[1exp(4σφ2)]sin22θ}1/2.

As can be seen from Eq. (22) and (23), the output degree of coherence in this case does not depend on the direction of polarization of the incident field and the output degree of polarization reaches it minimum value of exp(-2σ2 φ ) with the choice θ = π/4. The degree of coherence given by Eq. (22), for σφ = 2 and different values of γ, is plotted in Fig. 2. As can be seen from this figure, the curve ηout(ξ) is approximately asymptotic to the ξ, axis, and hence the transverse coherence length given by Eq. (17) tends to zero with the decrease of parameter γ. This result confirms that our object is successfully achieved.

 figure: Fig. 2.

Fig. 2. Degree of coherence given by Eq. (22), σφ = 2 and γ = 1,2,3.

Download Full Size | PDF

Concluding this section, we note that the random function φ(x), which serves to change the birefringence of the LC SLM, can readily be generated with the help of the control computer. When doing this one has to calculate the needed value of the parameter σφ corresponding to the desired value of P out(x) and then to choose the needed value of the parameter γ from the desired plot of function ηout(ξ).

4. Dependence of the results on alternative definition of the degree of coherence

The above analysis of coherence modulation was done using Wolf’s definition of electromagnetic coherence given by Eq. (3). However, there are other possible definitions of the electromagnetic coherence [11–14]. It is obvious that applying these alternative definitions, one can develop the analysis done in previous sections, obtaining different results. Here, we do not set ourselves the task of recalculating the results for all known definitions of electromagnetic coherence (it is not a trivial mathematical problem). We will limit ourselves to only one of the alternative definitions, namely the definition given in Ref. 11.

According to Ref. 11, the degree of coherence of the electromagnetic field is defined as the normalized Frobenius norm of the cross-spectral density matrix W(x 1, x 2), i.e.,

μ(x1,x2)=(Tr[W(x1,x2)W(x1,x2)TrW(x1,x1)TrW(x2,x2))1/2.

Then, substituting Eq. (21) into Eq. (24), we obtain

μout(ξ)=exp(σφ2)
×{exp[2σφ2exp(ξ22γ2)](sin4θ+cos4θ)+2exp[2σφ2exp(ξ22γ2)]sin2θcos2θ}12,

In particular, for the choice θ = π/4

μout(ξ)=12exp(σφ2){exp[2σφ2exp(ξ22γ2)]+exp[2σφ2exp(ξ22γ2)]}1/2.

The behavior of the degree of coherence given by Eq. (26) in comparison with the one given by Eq. (22) is illustrated in Fig. 3. The difference between curves ηout(ξ) and μ out(ξ) is not so large to affect significantly the transverse coherence length, i.e. the main result of the previous section.

 figure: Fig. 3.

Fig. 3. Degree of coherence given by Eq. (22) (solid line) and degree of coherence given by Eq. (26) (dotted line), σφ =2 and γ = 2.

Download Full Size | PDF

5. Conclusion

We have shown that the LC-based technique analyzed in Ref. 5 does not provide the appropriate simultaneous modulation of coherence and polarization. At the same time the proposed technique solves successfully this problem. Of course this technique is more complicated than the former one, but this is the price we pay to obtain the desired result. We have also shown that our results, as well as the results obtained in Ref. 5, depend on the used definition of the degree of electromagnetic coherence. Nevertheless this dependence can not be considered as significant in practice. The proposed technique can be used in practical applications to generate the secondary electromagnetic source with the desired statistical properties.

Appendix: To derivation of Eqs. (12) and (21)

Taking into account Eq. (10), one can state the relation

exp[±iφ(x)]=12πσφexp[±iφ(x)]exp(φ2(x)2σφ2).

Then, making use of well known Fourier-transform relation

exp(πa2φ2)exp(±i2πφu)=1aexp(πu2a2),

we find

exp[±iφ(x)]=exp(σφ22).

Now we will introduce a random variable

ψ(x1,x2)=φ(x2)±φ(x1).

Employing Eq. (10), one can write the probability distribution of this function as

p[ψ(x1,x2)]=12πσψexp(ψ2(x1,x2)2σψ2),

where

σψ2=ψ2(x1,x2)=[φ(x2)±φ(x1)]2,

Calculating the square in Eq. (A6) and applying Eq. (11), we find

σψ2=2σφ2[1±exp(ξ22γ2)].

Then, by analogy with derivation of Eq. (A3), but this time for argument ψ, we find

exp{+i[φ(x2)±φ(x1)]}=exp{i[φ(x2)±φ(x1)]}=exp{σφ2[1±exp(ξ22γ2)].}

Equations (A3) and (A8) are used in the text for derivation of Eqs. (12) and (21).

Acknowledgments

The authors gratefully acknowledge the financial support from the National Council for Science and Technology of Mexico and from the Autonomous University of Puebla under projects 55108 and OSA-EXT09-G, respectively. A.S. Ostrovsky thanks also the National Institute of Astrophysics, Optics and Electronics, where he is currently in his sabbatical leave.

References and links

1. J.W. Goodman, Statistical Optics (Wiley, New York, 1985).

2. A.S. Ostrovsky, Coherent-Mode Representations in Optics, (SPIE Press, Bellingham, WA, 2006). [CrossRef]  

3. J. Turunen, A. Vasara, and A.T. Friberg, “Propagation invariance and self-imaging in variable-coherence optics,” J. Opt. Soc. Am. A 8, 282–289 (1990). [CrossRef]  

4. A.S. Ostrovsky and E. Hernández-García, “Modulation of spatial coherence of optical field by means of liquid crystal light modulator,” Revista Mexicana de Física 51, 442–446 (2005).

5. T. Shirai and E. Wolf, “Coherence and polarization of electromagnetic beams modulated by random phase screens and their changes on propagation in free space,” J. Opt. Soc. Am. A 21, 1907–1916 (2004). [CrossRef]  

6. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, UK, 1995).

7. E. Wolf, “Unified theory of coherence and polarization of random electromagnetic beams,” Phys. Lett. A 312, 263–267 (2003). [CrossRef]  

8. E. Wolf, Introduction to the Theory of Coherence and Polarization of Light (Cambridge University Press, Cambridge, UK, 2007).

9. A. Yariv and P. Yen, Optical Waves in Crystals, (Wiley, New York, 1984).

10. K. Lu and B.E.A. Saleh, “Theory and design of the liquid crystal TV as an optical spatial phase modulator,” Opt. Eng. 29, 240–246 (1990). [CrossRef]  

11. J. Tervo, T. Setälä, and A.T. Friberg, “Theory of partially coherent electromagnetic fields in the space-frequency domain,” J. Opt. Soc. Am. A 21, 2205–2215 (2004). [CrossRef]  

12. P. Réfrégier and F. Goudail, “Invariant degrees of coherence of partially polarized light,” Opt. Express 13, 6051–6060 (2005). [CrossRef]   [PubMed]  

13. A. Luis, “Degree of coherence for vectorial electromagnetic fields as the distance between correlation matrices,” J. Opt. Soc. Am. A 24, 1063–1068 (2007). [CrossRef]  

14. F. Gori, M. Santarsiero, and R. Borghi, “Maximizing Young’s fringe visibility through reversible optical transformations,“ Opt. Lett. 32, 588–590 (2007). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (3)

Fig. 1.
Fig. 1. Degree of coherence given by Eq. (16), σφ = 2 and γ = 1,2,3.
Fig. 2.
Fig. 2. Degree of coherence given by Eq. (22), σφ = 2 and γ = 1,2,3.
Fig. 3.
Fig. 3. Degree of coherence given by Eq. (22) (solid line) and degree of coherence given by Eq. (26) (dotted line), σφ =2 and γ = 2.

Equations (37)

Equations on this page are rendered with MathJax. Learn more.

{ E ( x , v ) } = { E x ( x , v ) E y ( x , v ) } ,
W ( x 1 , x 2 ) = [ E x * ( x 1 ) E x ( x 2 ) E x * ( x 1 ) E y ( x 2 ) E y * ( x 1 ) E x ( x 2 ) E y * ( x 2 ) E y ( x 2 ) ] ,
η ( x 1 , x 2 ) = Tr W ( x 1 , x 2 ) [ Tr W ( x 1 , x 2 ) Tr W ( x 1 , x 2 ) ] 1 / 2 ,
P ( x ) = ( 1 4 Det W ( x , x ) [ Tr W ( x , x ) ] 2 ) 1 / 2 ,
W out ( x 1 , x 2 ) = T ( x 1 ) W in ( x 1 , x 2 ) T ( x 2 ) ,
W in ( x 1 , x 2 ) = E 0 2 exp ( x 1 2 + x 2 2 4 ε 2 ) [ cos 2 θ cos θ sin θ sin θ cos θ sin 2 θ ] ,
T 1 ( x ) = [ 1 0 0 exp [ i β 1 ( x ) ] ] ,
β 1 = π d λ ( n e n o )
β 1 ( x ) = φ 0 + φ ( x ) ,
p [ φ ( x ) ] = 1 2 π σ φ exp ( φ 2 ( x ) 2 σ φ 2 ) ,
φ ( x 1 ) φ ( x 2 ) = σ φ 2 exp ( ξ 2 2 γ 2 ) ,
W out ( x 1 , x 2 ) = E 0 2 exp ( x 1 2 + x 2 2 4 ε 2 )
× [ cos 2 θ exp ( i φ 0 ) exp ( σ φ 2 2 ) sin θ cos θ exp ( i φ 0 ) exp ( σ φ 2 2 ) sin θ cos θ exp { σ φ 2 [ 1 exp ( ξ 2 2 γ 2 ) ] } sin 2 θ ] .
η out ( ξ ) = cos 2 θ + exp { σ φ 2 [ 1 exp ( ξ 2 2 γ 2 ) ] } sin 2 θ ,
P out ( x ) = { 1 [ 1 exp ( σ φ 2 ) ] sin 2 2 θ } 1 / 2 .
η out ( ξ ) = exp { σ φ 2 [ 1 exp ( ξ 2 2 γ 2 ) ] } ,
η out ( ξ ) = 1 2 + 1 2 exp { σ φ 2 [ 1 exp ( ξ 2 2 γ 2 ) ] } .
Δ η out = 0 η out 2 ( ξ ) ,
T 2 ( x ) = [ exp [ i β 2 ( x ) ] 0 0 1 ] ,
β 2 ( x ) = φ 0 φ ( x ) ,
T ( x ) = T 2 ( x ) T 1 ( x ) = exp ( i φ 0 ) [ exp [ ( x ) ] 0 0 exp [ ( x ) ] ] .
W out ( x 1 , x 2 ) = E 0 2 exp ( x 1 2 + x 2 2 4 ε 2 ) exp ( σ φ 2 )
× [ exp [ σ φ 2 exp ( ξ 2 2 γ 2 ) ] sin 2 θ exp [ σ φ 2 exp ( ξ 2 2 γ 2 ) ] cos θ sin θ exp [ σ φ 2 exp ( ξ 2 2 γ 2 ) ] sin θ cos θ exp [ σ φ 2 exp ( ξ 2 2 γ 2 ) ] cos 2 θ ] .
η out ( ξ ) = exp { σ φ 2 [ 1 exp ( ξ 2 2 γ 2 ) ] } ,
P out ( x ) = { 1 [ 1 exp ( 4 σ φ 2 ) ] sin 2 2 θ } 1 / 2 .
μ ( x 1 , x 2 ) = ( Tr [ W ( x 1 , x 2 ) W ( x 1 , x 2 ) Tr W ( x 1 , x 1 ) Tr W ( x 2 , x 2 ) ) 1 / 2 .
μ out ( ξ ) = exp ( σ φ 2 )
× { exp [ 2 σ φ 2 exp ( ξ 2 2 γ 2 ) ] ( sin 4 θ + cos 4 θ ) + 2 exp [ 2 σ φ 2 exp ( ξ 2 2 γ 2 ) ] sin 2 θ cos 2 θ } 1 2 ,
μ out ( ξ ) = 1 2 exp ( σ φ 2 ) { exp [ 2 σ φ 2 exp ( ξ 2 2 γ 2 ) ] + exp [ 2 σ φ 2 exp ( ξ 2 2 γ 2 ) ] } 1 / 2 .
exp [ ± i φ ( x ) ] = 1 2 π σ φ exp [ ± i φ ( x ) ] exp ( φ 2 ( x ) 2 σ φ 2 ) .
exp ( π a 2 φ 2 ) exp ( ± i 2 πφ u ) = 1 a exp ( π u 2 a 2 ) ,
exp [ ± i φ ( x ) ] = exp ( σ φ 2 2 ) .
ψ ( x 1 , x 2 ) = φ ( x 2 ) ± φ ( x 1 ) .
p [ ψ ( x 1 , x 2 ) ] = 1 2 π σ ψ exp ( ψ 2 ( x 1 , x 2 ) 2 σ ψ 2 ) ,
σ ψ 2 = ψ 2 ( x 1 , x 2 ) = [ φ ( x 2 ) ± φ ( x 1 ) ] 2 ,
σ ψ 2 = 2 σ φ 2 [ 1 ± exp ( ξ 2 2 γ 2 ) ] .
exp { + i [ φ ( x 2 ) ± φ ( x 1 ) ] } = exp { i [ φ ( x 2 ) ± φ ( x 1 ) ] } = exp { σ φ 2 [ 1 ± exp ( ξ 2 2 γ 2 ) ] . }
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.