Advanced
Platform for Manipulating Polarization Modes Realized with Jones Vectors in MATHEMATICA
Platform for Manipulating Polarization Modes Realized with Jones Vectors in MATHEMATICA
Journal of Astronomy and Space Sciences. 2015. Jun, 32(2): 151-159
Copyright © 2015, The Korean Space Science Society
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : March 03, 2015
  • Accepted : April 04, 2015
  • Published : June 15, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Yong-Dae Choi
Department of Microbial and Nanomaterials, Mokwon University, Daejeon 302-729, Korea
Bogyeong Kim
Department of Astronomy and Space Science, Chungnam National University, Daejeon 305-764, Korea
Hee-Joong Yun
Korean Institute of Science and Technology Information, Daejeon 305-806, Korea
heejy@reseat.re.kr
Abstract
The fundamental conception in physics of the propagation of the electromagnetic wave polarization in matter is newly understood as the cardinal keyword in free-space quantum communication technology and cosmology in astrophysics. Interactive visualization of the propagation mechanism of polarized electromagnetism in a medium with its helicity has accordingly received attention from scientists exploiting the protocol of quantum key distribution (QKD) to guarantee unconditional security in cryptography communication. We have provided a dynamic polarization platform for presenting the polarization modes of a transverse electromagnetic wave, converting the state of polarization through the arrangement of optical elements, using Jones vectors calculations in Methematica . The platform graphically simulates the mechanism of production and propagation of the polarized waves in a medium while satisfying Maxwell's equations.
Keywords
1. INTRODUCTION
Polarization is a coherent characteristic of the electromagnetic (EM) wave in a medium. In a plane wave, both the electric field vector
PPT Slide
Lager Image
and magnetic field vector
PPT Slide
Lager Image
of the electromagnetic radfield vectoriation always oscillates parallel to a fixed direction in space. Light of such character is said to be linearly polarized, and maintains a constant direction of oscillation, and does vary spatially in a regular manner, producing either elliptically polarized or circularly polarized light ( Fowles 1975 ; Jackson 1975 ; Pedrotti & Pedrotti 1987 ). Polarization characteristic have been used in radio transmission to reduce interference between channels, particularly at VHF frequencies and beyond ( Masayoshi 2007 ; Yao et al. 2007 ). Free-space communication has forced the use of circular polarization, which has the fundamental advantage of precluding disturbance from the reflections of signals ( Leitch et al. 2002 ; Elser et al. 2009 ). Free-space optical (FSO) communication utilizes a spatial diversity receiver to receive the binary signals, which are modulated by two circular polarizations. FSO communication employing a binary polarization shift keying coherent modulation scheme are utilized in atmospheric turbulence channel ( Tang et al. 2010 ) and free-space quantum key distribution (QKD) by rotation-invariant twisted photons are used to guarantee unconditional security in cryptographic communication ( Vallone 2014 ). On 17 March 2014, astronomers from the Harvard Smithsonian Center for Astrophysics announced their detection of signature patterns of polarized light in the Cosmic Microwave Background (CMB) ( Ade 2014 ; Calvin 2014 ; CfA 2014 ). The team hunted for a special type of polarization called ‘B-mode’ which represents a twisting or ‘curl’ pattern in the polarized orientations of the ancient light. This is the strongest confirmation yet of the cosmic inflation theory ( Boyle 2006 ). Gravitational waves squeeze space as they travel, and the squeezing produces a distinct pattern in the CMB. Gravitational waves have a “handedness” much like light wave and can have left- and right-handed polarization.
In the Standard Model, the weak bosons (W±, Z) mediate the weak interactions between different flavors (all quarks and leptons). Experimental results have shown that all produced and observed neutrinos have left-handed helicity, and all antineutrinos have right-handed helicity ( Aad et al. 2012 ; CfA 2014 ). The helicity of the elementary particle could be a keyword to determine the Standard Model.
The scope of application of polarization has expanded explosively this decade with the development of communication technology. Accordingly, it is fundamentally important for application of a manipulation scheme to understand the physics of polarization’s conception and the process of producing polarization modes. Presently the definition of polarization has been modified, and its nomenclature upgraded, which can be confusing to students of physics and researchers. For instance, there are similar customary nomenclature for the right circularly polarization: right circularly ( Fowles 1975 ; Jackson 1975 ), right-circularly ( Pedrotti & Pedrotti 1987 ), right-hand circularly ( Reitz et al. 1993 ), and right-handed ( Born & Wolf 1999 ; Georgi 1982 ) circularly polarization. Right circularly polarization and right-handed polarization are different types of polarization, although with similar names. Recently, the visualization of polarization propagating in matter has drawn physicist's attention for potential applications in modern physics and information technology ( Tamm 1997 ; Mooleskamp & Stokes 2015 ; Yun & Choi 2013 ).
We have provided a dynamic polarization modes platform for simulating polarization modes with Jones matrices calculations, corresponding to the physical arrangement of optical elements, in a Mathematica computing environment ( Mathematica 2015 ).
2.MATHEMATICASIMULATION FOR THE CONVERTING POLARIZATION MODES
- 2.1 Electromagnetic waves in solids
The propagating process of electromagnetic waves in solids is different from that in the vacuum, since
PPT Slide
Lager Image
and
PPT Slide
Lager Image
of the electromagnetic waves interact with the electrons in a solid ( Pedrotti & Pedrotti 1987 ; Fowles 1975 ). In particular, although the electromagnetic waves are an harmonic plane wave, the fields may pull or push the electrons in the orbital of a solid, which is responsible for inducing dipole moment
PPT Slide
Lager Image
and magnetization
PPT Slide
Lager Image
im the solid. If we assume that the medium is not a magnetic material and
PPT Slide
Lager Image
, the Helmholtz wave equation in a solid ( Fowles 1975 ; Jackson 1975 ) is
PPT Slide
Lager Image
where
PPT Slide
Lager Image
. If we suppose the EM wave will be a form of solution such as
PPT Slide
Lager Image
, we can rewrite Eq. (1) using a wave vector
PPT Slide
Lager Image
PPT Slide
Lager Image
It then becomes a vector equation for
PPT Slide
Lager Image
, which indicates that the propagation process varies with the component of electric susceptibility x . Therefore, we will write the
PPT Slide
Lager Image
solution of Eq. (2) in a solid as an inhomogeneous plane wave solution below ( Pedrotti & pedrotti 1987 ; Fowles 1975 ; Jackson 1975 )
PPT Slide
Lager Image
Aa shown above,
PPT Slide
Lager Image
field vector may be presented as the Jones vector ( Jones 1941 ) with a complex vector amplitude {A, B ± iC} oscillating in an inhomogeneous plane. The wave vector
PPT Slide
Lager Image
and form a real mutually orthogonal unit vectors
PPT Slide
Lager Image
. Here
PPT Slide
Lager Image
is a complex wave vector and N = n + iκ is the complex refraction index. For Eq. (3) to be a solution of Eq. (2) as a homogeneous plane harmonic wave, it should be,
PPT Slide
Lager Image
. Then we get the relations:
PPT Slide
Lager Image
, and
PPT Slide
Lager Image
, which result in propagation speeds that are different along the direction in the medium. Therefore, there will be a cumulative phase difference Δ ε between the two components of the
PPT Slide
Lager Image
field vector as they emerge in uniaxial crystals (Quartz, Calcite, etc.). After the wave has traveled a distance d , the phase difference is
PPT Slide
Lager Image
between E x wave and E y wave when the radiation is propagating along the
PPT Slide
Lager Image
direction. If the Δ ε =0 while the amplitude of
PPT Slide
Lager Image
is real, the vector is responsible for the linearly polarized such as Jones vector {A, B}, otherwise the amplitude is the complex vector responsible for the elliptically polarized Jones vector as {A, B±C}. Specifically, if B=0 and A=C then the wave is a circularly polarized Jones vector such as {1, i}.
A quarter wave plate is a thin birefringent crystal the thickness of which has been adjusted to produce a ± π /4 phase difference between the ordinary and extraordinary rays at the operating wavelength. We desire a matrix that will transform the element E 0x e x into E 0x e i(εxx) and E 0y e y into E 0y e i(εyy) . The general form of a matrix representing a phase retarder will transform the elements by the matrix operation as follows
PPT Slide
Lager Image
where ε x and ε y represent the advance in phase of E x - and E y - component of the incident light. As an example, consider a quarter-wave plate (QWP) which makes Δ ε = π /2. We may write the Jones matrix, M transforming the Jones vector which makes Δ ε = π /2 to produce right-handed circular polarizing light as
PPT Slide
Lager Image
This is the case of fast axis vertical (FA vertical). Similarly, we can determine the corresponding Jones matrix for a half-wave plate (HWP) or eighth-wave plate (EWP) or arbitrary phase of retarded. Jones matrices derived for various wave plates are summarized in Table 1 .
Summary of Jones vectors in the most common and Jones matrices of the wave plates and phase retarders (Pedrotti & Perotti 1987).
PPT Slide
Lager Image
†Wave plate(V/H) stand for the fast transmission axis of vertical/horizontal.
We desire now to create a new simulation presenting a means of producing polarizing modes from a Jones calculation corresponding to the physical arrangement of optical elements in the Mathematica computing environment. In a plane wave, the electric field vector
PPT Slide
Lager Image
always oscillates parallel to the fixed direction in space. Light of such character is said to be linearly polarized. If the linearly polarized light passes through a quarter-wave plate, elliptically polarized light emerges. The same can be said of the magnetic field vector
PPT Slide
Lager Image
, which maintains an orientation perpendicular to the electric field vector such that the direction of
PPT Slide
Lager Image
is everywhere the direction of wave propagation. Thus, the possibility of polarizing light is essentially due to its transverse character. Therefore, the Mathematica simulation should show the transverse character of the
PPT Slide
Lager Image
with its vectorial behaviors dynamically satisfying Maxwell's wave equations. In addition to the transverse character of the polarizing wave, the helicity of a polarized wave is a critical factor in quantum cryptography communication technology or cosmology in moderen physics.
- 2.2 Polarization modesMathematicasimulation
We have used Mathematica to implement interactive Jones matrix calculations and animations for the generation and propagation of the polarization modes in th solid state. First, we desire to confirm that the complex vector field
PPT Slide
Lager Image
of Eq. (3) with Jones vectors in the matter satisfy the Maxwell's vector equations Eqs. (6) and (7) in both numeric and graphic simulations.
PPT Slide
Lager Image
PPT Slide
Lager Image
where,
PPT Slide
Lager Image
,
PPT Slide
Lager Image
, and
PPT Slide
Lager Image
ard both complex vectors and
PPT Slide
Lager Image
,
PPT Slide
Lager Image
and
PPT Slide
Lager Image
are real vectors with
PPT Slide
Lager Image
where
PPT Slide
Lager Image
is a magnetic intensity vector in the matter. We examine the process of generating polarzation modes with the normalized Jones vectors calculations in the Mathematica simulation as below. The Mathematica input code is shown below:
PPT Slide
Lager Image
Mathematica returns calculations;
PPT Slide
Lager Image
The Mathematica calculations show the emerging polarization mode through the physical arrangement of optic elements, which promptly confirm those Jones vectors from the animating platform. For instance, in In[23]:= E0*LP∙QWV/. α ->1/4 π produces Out[23]:= {I/2, I/2}, that is the right-handed circularly polarized light (RHCP) as shown in Fig. 2 . From the arrangement of linearly polarizers in right angle no wave emerging shown as Out[29]:= {0, 0} while the Out[28]:={0.341506, 0.591506} which show the linear polarization (see Figs. 3(g) and 3(h) ). This enables us to switch on or off the polarization by the combination of arrangements of polarizers and wave plates.
- 2.3 Helicity of the elliptically polarized wave
The handedness of an elementary particle depends on the correlation between its spin and its momentum ( Goldhaver et al. 1958 ). If the spin and momentum are parallel, the particle can be said to be right-handed or have a helicity of 1. If they are anti parallel, the particle is left-handed or have a helicity of -1. We may also adopt this definition to modern optics, since the circularly polarized electromagnetic wave is just a helical motion with helicity. The helicity of the polarized electromagnetic wave is a critical factor in modern communication technology and photonics ( Aad et al. 2012 ; Goldhaver et al. 1958 ; Rubenhok et al. 2013 ). However, determination of the helicity is a perplex issue because we need to aware of the vectorial behaviors accurately unless we may observe propagating polarized wave ( Goldhaver et al. 1958 ). Therefore, it is very helpful to simulate advanced polarization modes in the medium, and further, it would be more helpful to evaluate the helicity operator in a physics state.
Here we simply estimate the helicity with a calculation of phase shift in the process in the Mathematica ( Wolfram 2015a ). We can determine the helicity based on a calculation of the phase shift in the block respectively; phase shifts of the E 1 =E0*Jv after passing through the polarizer and of E 2 =Jm∙E 1 after the wave plate respectively. For evaluation of the helicity of the polarized waves in progression, we have provided a Mathematica module, helicity [JvC_,JmC_]:=Module[{E1, E2, dp11, dp12, dp13, dp21, dp22, dp23}]. If you type in two numbers assigned to the polarizer JvC and wave plate JmC in the module, then helicity [JvC, JmC] will return helicity with a list of phase shifts. The outputs of the helicity [JvC, JmC] are shown as below:
PPT Slide
Lager Image
The helicity[1, 2] returns two helicities with two lists of three phase differences in the two divisions for the optical device arrangement of linearly polarizer (JvC=1) and quarter-wave plate (JmC=2). The helicity[2, 7] returns the helicities of the optical arrangement of left-circularly polarizer (JvC=2) and eighth-wave plate (JmC=7). It needs to be noticed that the three phase differences are all 1.5708 ( π /2) of the helicity[1, 2] which interprets the vertical component E 2z lead continuously the horizontal component E 2y with a constant phase shift ( π /2) on the yz plane of
PPT Slide
Lager Image
field vector propagating in x 1 direction perpendicular to this plane with velocity ω /k. This results in
PPT Slide
Lager Image
vector rotating in a counter colckwise direction
PPT Slide
Lager Image
around the advancing x 1 direction. That is, if we grasp our finger along the spin direction with right hand, then thumb directs the advancing direction of propagation wave, hence the wave is right-handed polarized wave regardless of viewer . This is very helpful in confirming the helicity of the polarized wave on the end block of the platform as shown in Figs. 1 and 2 . We can correctly determine the helicity of the polarized wave in the E2 block by seeing the spin direction of the wave; either the spin direction is parallel
PPT Slide
Lager Image
or anti parallel
PPT Slide
Lager Image
to the advancing x 1 direction (momentum direction). Helicity of E 2 can changed simply by the choice of JmC optical element like as In[83] and In[84] codes. The half-wave plate (HWH) shifts phase - π in helicity[2, 5] while eighth-wave plate (EWH) shifts phase - π /4 in helicity[2, 7]. Furthermore we may easily calculate the helicities of the wave train on the helicity for any combination of optical elements in Mathematica , so that we can predict exactly helicity of the prod[JvC, JmC] moduleucing polarization modes. For more information for the helicity [JvC, JmC] module refer to citation ( Yun 2015 ).
PPT Slide
Lager Image
The starting platform of the dynamic polarization modes platform with Jones vectors (jdpmp). Unpolarized EM waves 0 injected to the polarizer (LP) have not passed the polarizer yet. The simulation will run when you click the ▶ appearing while you spread the ⊕ of t1 panel of the platform.
PPT Slide
Lager Image
Dynamic polarization modes platform (jdpmp). The picture shown is a snapshot of the right-handed circular polarized (RHCP) propagation EM wave train. Unpolarized EM waves 0 injected into the linear polarizer (LP) with 45° to the horizontal transmission axis come from the vacuum and pass through the quarter-wave plate (QWPV), which result in the right-handed circular polarized EM wave train. This picture shows the E0*LP∙QWPV/.α->1/4π process in the Jones vectors manipulation with a helicity of +1.
3. DYNAMIC POLARIZATION PLATFORM WITH JONES VECTORS IN MATHEMATICA
We have provided a dynamic polarization platform with Jones vectors (jdpmp) in the Graphics3D in Mathematica , which simulates the polarizing modes dynamically while presenting the helicity of a running polarization mode. The platform manipulates three zones graphically using the Piecewise function of Mathematica depending on the polarizer and wave plate, with the Manipulate function of Mathematica . To present the transverse characteristic property satisfying Eq. (1) together with Eqs. (6) and (7), we used the Arrow function in Mathematica for drawing the vector array: Table[ Arrow [ { {x1, 0, 0}, {x1, Ey, Ez} } ] ] in the orthogonal {x1, Ey, Ez} coordinate system. In addition, the same was used for the magnetic field vector
PPT Slide
Lager Image
calculated from the relation
PPT Slide
Lager Image
of Eq. (6) also in solid. Graphics3D can draw the vector array of the
PPT Slide
Lager Image
(blue) and
PPT Slide
Lager Image
(red) fields of every point in the block as real polarization modes propagating in the advancing
PPT Slide
Lager Image
direction, in our case, advancing in the x 1 direction. This simulation scheme differs basically from that of an animation representing the envelope of the polarization propagation of
PPT Slide
Lager Image
field only using the Animate function of graphic tools ( Harrison 2015 ; Mooleskamp & Stokes 2015 ). Also, in this work, we simulated polarization process with vector fields satisfying wave equation and Maxwell's vector equations at every point during the process. Similar numerical method has been applied once to solve the Lagrange's equation in non-inertial frame ( Kim & Yun 2014 ).
The platform starts by clicking the button ▶ on pop up ⊕ of the t1 panel on the platform in Fig. 1 . We can choose an incoming angle to the polarizer and rotate the polarizer to compose an arrangement of devices for the polarization mode indent shown as Fig. 2 . While the program is running, we will change the polarizer or wave plate, with the changed animation running continuously. Even if the animation is stopped, we can change the configuration of the platform and observe the polarizing mode of the changed state.
Fig. 3 shows snapshots of three kinds of polarization modes. Figs. 3a and 3b shows a circularly polarization mode, Figs. 3c and 3d linearly polarization, Figs. 3e and 3f elliptically polarization according to the arrangement of optical devices. If we utilize linear polarizer JmLinear [ β ] instead of a wave plate, then either it passes a linearly polarized wave or block in the angle β . While making right angles between two polarizers no wave passes through the polarizer plate as shown in Figs. 3g and 3h . Thus, controlling of phase of polarized light waves is available for the phase encoding in the quantum key distribution (QKD) of the unconditional security in quantum communication technology ( Rubenhok et al. 2013 ).
PPT Slide
Lager Image
Producing different polarized modes by the different physical arrangements of polarizer and wave plate: (a) RHCP of helicity +1, (b) LHCP of helicity −1, (c) LP, (d) LP, (e) RHEP of helicity +1, (f) LHEP of helicity −1, (g) LP, (h) No wave.
Fig. 4 presents snapshots of three kinds of polarization modes viewed in the ViewPoint (200, 0, 0) against the propagation direction for confirmation of the transverse behavior and progression of the Poynting vector. It shows the orthogonal transverse of the propagating process of the waves and confirms the flow of Poynting vector. This is the visualization of the graphic version of
PPT Slide
Lager Image
satisfying Maxwell’s vector equations Eqs. (6) and (7). The traces of both
PPT Slide
Lager Image
and
PPT Slide
Lager Image
are dynamically shown in the polarization modes at the edge of the jdpmp platform.
PPT Slide
Lager Image
Snapshots of the polarized modes at the ViewPoint (200, 0, 0); (a) LP of Fig. 3c, (b) No wave of Fig.3h, (c) RHCP of Fig. 3a, (d) RHEP of Fig. 3e. It presents fields by the blue arrows and fields by the red arrows.
The platform will run faster or slower and stop and restart again by clicking the pop-up menu of the tl panel for more analytical observations while the polarization mode is running. During the platform stop, you will change to another mode only if you click a panel, then the changed mode will be presented automatically. At that time you can save the presenting mode or print out the mode status. While the platform is running, you can confirm the helicity by observing the spin direction of the field vectors along the orange trace of the helical propagation of the
PPT Slide
Lager Image
and
PPT Slide
Lager Image
vector fields at the edge of the platform. The helicities on the helicity panel are correct in the case of E1=E0*LP(α) manipulation on the jdpmp platform. The helicities of other arrangements of optical devices are calculated promptly on the helicity (JvC, JmC) module in Mathematica .
As far as we know, other than jdpmp , there is no instructive platform that simulates a transverse EM wave satisfying Maxwell's equation in the vectorial version of polarization modes ( Harrison 2015 ; Mooleskamp & Stokes 2005 ; Tamm 1997 ). The complete version of jdpmp.nb program including helicity [JvC, JmC] module is available from the site ( Yun 2015 ). If your PC has not installed Mathematica , you can run the jdpmp.cdf instead on the Wolfram CDF Player ( Wolfram 2015 ).
4. SUMMARY
We have provided a dynamic polarization mode platform, jdpmp , for simulating and producing polarization modes with the Jones calculations, corresponding to the physical arrangement of optical elements, in Mathematica . The platform animates the polarization process of the
PPT Slide
Lager Image
vector field together with
PPT Slide
Lager Image
using the Arrow vectors, so that the transverse EM wave advance in the propagation direction, satisfying Maxwell's wave equation at evey point on the advancing axes of the platform. Consequently, the vectors
PPT Slide
Lager Image
,
PPT Slide
Lager Image
,
PPT Slide
Lager Image
(in that order) form a right-hand orthogonal set. jdpmp can be accomplished graphically in Graphics3D of Mathematica based on the numerical vector calculations in Mathematica . The platform can be manipulated dynamically and interactively by advancing the polarized mode while clicking the panels of the polarizer and wave plate using the Manipulate function in Mathematica , so that the program will continuously simulate the changed mode. While the program is running the helicity of the polarized mode is displayed on the panel. The module helicity [JvC, JmC] will return the helicity of the process on Notebook of Mathematica when you're typing in a polarizer and a wave plate. The module is helpful for students or researchers to inspect the phase difference or helicity of a polarized mode for various kinds of physics arrangements. We expect the platform jdpmp will be a useful starting platform for physics students and researchers to explore polarizations in science.
Acknowledgements
This research was partially supported by the ReSEAT program funded by the Korean Ministry of Education, Science and Technology through, the National Research Foundation of Korea and the Korea Lottery Commission grants.
References
Aad G , Abbott B , Abdallah J , Khalek SA , Abdelalim AA 2012 Search for the Standard Model Higgs Boson in the Diphoton Decay Channel with 4.9 fd−1 of pp Collision Data atTeV with ATLAS Phys. Rev. Lett. http://dx.doi.org/10.1103/PhysRevLett.108.111803 108 111803 -
Ade PAR , Aikin RW , Bento SJ , Bischoff CA , Bock JJ 2014 Detection of B-Model Polarization at Degree Angular Scales by BICEP Phys. Rev. Lett. http://dx.doi.org/10.1103/PhysRevLett.112.241101 112 241101 -
Born M , Wolf E 1999 Principle of Optics 7th ed. Cambridge University Press Cambridge 795 -
Boyle LA , Steinhardt PJ , Turok N 2006 Inflationary predictions for scalar and tensor fluctuations reconsidered Phys.Rev. Lett. http://dx.doi.org/10.1103/PhysRevLett.96.111301 96 111301 -
Harrison DM Optics Applets [Internet], Available from:
Calvin W NASA Technology Views Birth of the Universe [Internet], Available from:
CfA, First Direct evidence of cosmic inflation [Internet], Available from:
Elser D , Bartley T , Heim B , Wittmann Ch , Sych D 2009 Feasibility of free space quantum key distribution with coherent polarization states New J. Phys. http://dx.doi.org/10.1088/1367-2630/11/4/045014 11 045014 -
Fowles GR 1975 Introduction to Modern Optics Holt, Reinhart and Winston
Georgi H 1982 Lie Algebra in Particle Physics The Benjamin INC 162 -
Goldhaver M , Grodzins L , Sunyar AW 1958 Helicity of Neutrinos Phys. Rev. Lett. 109 1015 -
Jackson JD 1975 Classical Electrodynamics John Wiley & Sons 274 -
Jones RC 1941 A New Calculus for the Treatment of Optical Systems J. Opt. Soc. Am. http://dx.doi.org/10.1364/JOSA.31.000488 31 488 - 493
Kim B , Yun HJ 2014 Pedagogical Mathematica Platform Visualizing the Coriolis Effects in 3-Cell Atmospheric Circulation Model J. Astron. Space Sci. http://dx.doi.org/10.5140/JASS.2014.31.1.91 31 91 - 99
Leitch EM , Kovac JM , Pryke C , Carlstrom JE , Halverson NW 2002 Measurement of Polarization with the Degree Angular Scale Interferometer Nature http://dx.doi.org/10.1038/nature01272 420 673 - 771
Masayoshi T 2007 Cutting-edge terahertz technology Nat. Photonics http://dx.doi.org/10.1038/nphoton.2007.3 1 97 -
Mathematica, Wolfram Mathematica [Internet], Available from:
Mooleskamp FE , Stokes KL Polarization of an Optical Wave through Polarizers and Wave Plates [Internet], Available from:
Pedrotti FL , Pedrotti LS 1987 Introduction to Optics Prentice-Hall
Reitz R , Milford FJ , Christy W 1993 Foundation of Electromagnetic Theory 4th ed. Addison Wesley Reading, Massachusetts
Rubenhok A , Slater JA , Chan P , Lucio-Martinez I , Tittel W 2013 Real-World Two-Photon Interference and Proof-of-Principle Quantum Key Distribution Immune to Detector Attacks Phys. Rev. Lett. http://dx.doi.org/10.1103/PhysRevLett.111.130501 111 130501 -
Tamm P 1997 A Physicist's Guide to Mathematica Academic Press San Diego 291 - 293
Tang X , Ghassemlooy Z , Rajbhandari S , Popoola WO , Lee CG Free-space optical communication employing polarization shift keying coherent modulation in atmospheric turbulence channel 7th International Symposium Newcastle, England 21-23 July 2010
Vallone G , D'Ambrosio V , Sponselli A , Slussarenko S , Lorenzo, Marrucci L 2014 Free-Space Quantum Key Distribution by Rotation-Invariant Twisted Photons Phys. Rev. Lett. http://dx.doi.org/10.1103/PhysRevLett.113.060503 113 060503 -
Wolfram, Wolfram CDF Player for Interactive Computable Document Format [Internet], Available from:
Yao XS , Yan LS , Zhang B , Willer AE , Jiang J 2007 All-optic scheme for automatic polarization division demultiplexing Opt. Express http://dx.doi.org/10.1364/OE.15.007407 15 7407 - 7414
Yun HJ , Choi YD 2013 Dynamic Polarization Mathematica Platform Realized with Poynting Vectors New Phys.: Sae Mulli http://dx.doi.org/10.3938/NPSM.63.1118 63 1118 - 1127
Yun HJ jdpmp programs [Internet], Available from: http://home.mokwon.ac.kr/-heejy/program.htm