Issue 
A&A
Volume 568, August 2014



Article Number  A8  
Number of page(s)  12  
Section  Astronomical instrumentation  
DOI  https://doi.org/10.1051/00046361/201323229  
Published online  05 August 2014 
Combined calculi for photon orbital and spin angular momenta
^{1} Zentrum für Astronomie der Universität Heidelberg, Landessternwarte, Königstuhl 12, 69117 Heidelberg, Germany
^{2} OAM Solutions, LLC, 9300 Stardust Trail, Flagstaff AZ 86004, USA
email: nick.elias@oamsolutions.com
Received: 10 December 2013
Accepted: 12 May 2014
Context. Wavelength, photon spin angular momentum (PSAM), and photon orbital angular momentum (POAM), completely describe the state of a photon or an electric field (an ensemble of photons). Wavelength relates directly to energy and linear momentum, the corresponding kinetic quantities. PSAM and POAM, themselves kinetic quantities, are colloquially known as polarization and optical vortices, respectively. Astrophysical sources emit photons that carry this information.
Aims. PSAM characteristics of an electric field (intensity) are compactly described by the Jones (Stokes/Mueller) calculus. Similarly, I created calculi to represent POAM characteristics of electric fields and intensities in an astrophysical context. Adding wavelength dependence to all of these calculi is trivial. The next logical steps are to 1) form photon total angular momentum (PTAM = POAM + PSAM) calculi; 2) prove their validity using operators and expectation values; and 3) show that instrumental PSAM can affect measured POAM values for certain types of electric fields.
Methods. I derive the PTAM calculi of electric fields and intensities by combining the POAM and PSAM calculi. I show how these quantities propagate from celestial sphere to image plane. I also form the PTAM operator (the sum of the POAM and PSAM operators), with and without instrumental PSAM, and calculate the corresponding expectation values.
Results. Apart from the vector, matrix, dot product, and direct product symbols, the PTAM and POAM calculi appear superficially identical. I provide tables with all possible forms of PTAM calculi. I prove that PTAM expectation values are correct for instruments with and without instrumental PSAM. I also show that POAM measurements of “unfactored” PTAM electric fields passing through nonzero instrumental circular PSAM can be biased.
Conclusions. The combined PTAM calculi provide insight into mathematically modeling PTAM sources and calibrating POAM and PSAMinduced measurement errors.
Key words: instrumentation: miscellaneous / methods: analytical / methods: miscellaneous / methods: observational / techniques: miscellaneous
© ESO, 2014
1. Introduction
Elias (2008) derived propagation calculi to describe astronomical photon orbital angular momentum (POAM; colloquially known as optical vortices). He employed a semiclassical/semiquantum framework where electric fields are analogous to photon wave functions and intensities are analogous to probabilities. These calculi link POAM quantities on the celestial sphere to POAM quantities at instrument backends. He tacitly assumes that the electric fields on the celestial sphere are spatially uncorrelated (the “Standard Astronomical Assumption”, or SAA). Elias (2012) used these calculi to describe POAM and torque metrics for single telescopes and interferometers.
Like most other workers in the POAM field, Elias (2008, 2012) dealt only with optical systems that ignored photon spin angular momentum (PSAM; colloquially known as polarization), in order to simplify calculations. Since POAM and PSAM are complementary properties that will eventually be measured simultaneously, combined calculi are required for modeling source and instrument behavior. Failing to take nonzero PSAM into account will yield incorrect POAM values under certain conditions.
The goals and results of this work are multiple. First, I present the simplest and most general photon total angular momentum (PTAM = POAM + PSAM) electric field forms. Second, I combine the POAM (Elias 2008) and PSAM propagation calculi to create the PTAM propagation calculi. Third, I create the POAM, PSAM, and PTAM operators and calculate the corresponding expectation values for perfect and imperfect instruments. Last, I show that POAM measurements can be biased when unfactored electric fields pass through nonzero instrumental circular PSAM.
2. Electric fields
Elias (2008, 2012) treated the electric field as a scalar quantity when he constructed the POAM state expansions (1)where is the vector in a plane (e.g., celestial sphere, image plane, etc.), t is time, the m are the POAM quantum numbers (−∞ ≤ m ≤ ∞), and the E_{m}(H;t) are the POAM states (azimuthal Fourier components of the electric field) for each radius (perpendicular to the propagation direction). An azimuthal Fourier series is performed for each H and t.
The standard way of describing the PSAM behavior of an electric field is the Jones vector (2)where and are the rightcircular and leftcircular components. The circular PSAM basis is ideal for the subsequent analyses of this paper.
A “factored” electric field, where the POAM and PSAM parts are separate factors, is the simplest and likely the most well known PTAM form (3a)The most general “unfactored” PTAM form, on the other hand, comes from independently expanding each PSAM component into POAM components, or (3b)Although POAM and PSAM appear completely intertwined, these expressions work with the PTAM calculi (Sect. 3) and lead to the correct operators and expectation values (Sects. 4 and 5). This type of electric field can be prepared in the laboratory. At present there are no known astrophysical mechanisms that generate unfactored PTAM, but I use this form anyway for the sake of mathematical completeness and in the event that such mechanisms will eventually be found.
3. PTAM calculi
Elias (2008) created POAM propagation calculi for electric fields and intensities using SAA. He also treated the electric fields as scalars, ignoring PSAM. Their timeaveraged square magnitudes are intensities, which are analogous to Stokes I.
In this section, I combine the Elias (2008) POAM calculi with the electricfield PSAM calculi of Jones (1941) and the intensity PSAM calculi of Stokes (1852) and Mueller (1948) to form the PTAM calculi. I also employ the mathematics of Schmeider (1969) and Barakat (1981), hereafter collectively SB, to more easily link Jones vectors, Jones matrices, Stokes vectors, and Mueller matrices via direct products and coherence matrices (assuming no system depolarization).
3.1. POAM correlations
According to SB, the coherence vector is the direct product of the electric field from Eq. (2) (4)where ⊗ is the direct (outer) product, and ⟨⟩ is the time average. If I substitute Eq. (3b) (instead of Eq. (2)) into Eq. (4), I obtain (5a)where (5b)is the (m,n)th POAM correlation of the coherence vector. By comparing Eqs. (4) and (5ab), I find that the individual correlations can be expanded into double sums and (6d)The have a similar form to and the identical units as the I_{m,n}(H) POAM correlations defined by Elias (2008).
Coherence vectors are used mostly by engineers because some instruments, such as radio interferometers, employ right and leftcircular feeds. Astronomers prefer Stokes vectors because their components represent the total intensity and the polarization parameters required for scientific analysis. The Stokes vector is related to the coherence vector via a simple matrix transformation (SB) (7a)where · is the dot (inner) product, and (7b)is the coherencetoStokes transformation matrix in the circular basis. If I substitute Eqs. (5ab) into Eq. (7a), I obtain (8a)where (8b)is the (m,n)th POAM correlation of the Stokes vector. By comparing Eqs. (7a) and (8ab), I find that the individual Stokes parameters can be expanded into double sums and (9d)The have a similar form to and the same units as the I_{m,n}(H) POAM correlations defined by Elias (2008). As a matter of fact, Eq. (9a) is identical to the expansion derived by Elias (2008) using scalar electric fields.
The Stokes Q, U, and V expansions are unnecessary, so I rewrite the Stokes vector as (10)This PTAM form maintains both the POAM and PSAM information while minimizing complications.
3.2. POAM rancors
Elias (2008) defined a quantity called rancor, which is the azimuthal Fourier series versus radius of the intensity (11a)where ℐ_{m}(H) is the mth POAM rancor. This quantity is interesting because it identical to the infinite sum over a subset of POAM correlations (11b)Rancors, which may be easier to determine in some cases, contain a limited amount of POAM information. As an analogy, I point out that squared visbilities and closure phases in optical interferometry can provide important physical data about astronomical sources, in spite of the fact that they contain less information than complex visibilities.
In Sect. 3.1, I combined POAM correlations with PSAM Stokes vectors. Since rancors can be expressed in terms of correlations, it follows that all intensity formulae in Sect. 3.1 can be written in terms of rancors. I will not list all possible expressions here, since those expansions are identical to those in Elias (2008, Sects. 4, 5, Appendix C) apart from the fact that scalar quantities are replaced by vectors and matrices and scalar products are replaced by dot and direct products.
3.3. Propagating POAM quantities
Elias (2008) derived scalar electricfield and intensity calculi for propagating POAM from celestial sphere to image plane and listed them in several tables. He employed system forms and SAA. In this section, I extend these expressions to combine POAM and PSAM propagation calculi, thus creating PTAM propagation calculi.
Consider the system form for propagation of the scalar electric field from celestial sphere to image plane (12a)where (ρ′cosφ′,ρ′sinφ′) is the coordinate in the image plane, (ρcosφ,ρsinφ) is the coordinate on the celestial sphere, (12b)is the diffraction function, is the Fraunhofer propagator (it can be replaced with the Fresnel propagator), is the coordinate in the pupil plane normalized by wavelength, and is the pupil function which describes the telescope aberrations, atmospheric turbulence, etc. If these scalar electric fields are changed to 2 × 1 Jones vectors, the diffraction function must become a 2 × 2 Jones matrix (12c)In principle, the scalar and matrix diffraction functions can also be functions of time, although their variability time scales are much slower than those of the electric fields.
POAM expansions of .
POAM expansions of .
If the Jones vector components are expanded into independent POAM states (unfactored form, Eq. (3b)), Eq. (12c) becomes the PTAM state expansion (13a)where (13b)is the output POAM state p, (13c)is the input POAM state m, and (13d)is the diffraction function gain between output POAM state p and the input POAM state m. I summarize all PTAM electric field expansions in Tables 1 and 2. is not a true PTAM state, which means that the input expansion is of limited use but included for the sake of completeness.
The intensity is the squared magnitude of the electric field. Using SAA and Eq. (12a), the scalar intensity becomes (14a)where is the pointspread function (PSF), and . SAA collapses one of the integrals over the celestial sphere. If I employ Eqs. (4), (7a), and (12c) as well as SAA, the scalar Eq. (14a) becomes the vector equation (14b)The point spread function is now a 4 × 4 Mueller matrix. I summarize all PTAM intensity expansions in Tables 3 and 4. and are not true PTAM quantities, which means that the input expansions are of limited use but included for the sake of completeness. Also, note that the intensity equations are cannot be derived from the electric field equations when a system has depolarization (Mueller matrices cannot be uniquely determined from Jones matrices).
Now consider the StokesI parameter the image plane (15)where is the detector operator, and the are the elements of the top row of the Mueller matrix PSF. In Sects. 3.1 and 3.2, I point out that only the StokesI parameter must be expanded in terms of POAM correlations and rancors, even though the complete derivations involve POAMlike expansions of the other Stokes parameters. Similarly, only the upperleft element of the Mueller matrix must be expanded in terms of POAM correlations or rancors. Equation (15) indicates that nonzero Stokes Q, U, and V terms could introduce measurement biases which must be calibrated when measuring the POAM of the StokesI parameter. I present a simple example in Sect. 5 using operators and expectation values.
4. Operators and expectation values
Expectation values are specific quantities that can be measured by instruments. In this section, I: 1) define the POAM, PSAM, and PTAM operators; 2) derive the corresponding expectation values; 3) show how the operators and expectation values are modified by imperfect instruments.
4.1. Perfect instrument
In the paraxial case, the scalar quantum mechanical POAM operator along the + z propagation axis is (16)where , and ħ is Planck’s constant h divided by 2π. The POAM expectation value is measured when this operator is applied to the scalar electric field (17a)where (17b)is the integrated intensity of the scalar electric field. The numerator is a quantummechanicslike product of states and matrix elements, and the denominator is the normalization. Substituting Eqs. (1) and (16) into Eq. (17a), I obtain (18a)where M_{eff} is the effective quantum number, (18b)is the probability of a photon (or an enesmble of photons) being in state m, r_{max} is the maximum radius which contains all of the flux, and I_{m,m} is the radially integrated autocorrelation of POAM state m. The expection value is simply the effective quantum number times ħ.
Similarly, the quantum mechanical PSAM operator along the + z propagation axis is (19)where is the third Pauli spin matrix. The PSAM expectation value is measured when this operator is applied to the vector electric field (20a)where (20b)is the integrated intensity of the vector electric field, and the T superscript indicates the transpose. Substituting Eqs. (2) and (19) into Eq. (20a), I obtain (21a)where v is the normalized StokesV parameter, (21b)and (21c)are the probabilities of a photon (or an ensemble of photons) being in the RCP and LCP states, and I_{R,R} and I_{L,L} are the integrated autocorrelations of the RCP and LCP states. For an unpolarized and/or linearly polarized source v = 0, which means that . For a fully circularly polarized source, v = + 1 (v = −1), p_{R,R} = 1 and p_{L,L} = 0 (p_{R,R} = 0 and p_{L,L} = 1).
The PTAM expectation value is the sum of the POAM and PSAM expectation values, or . The PTAM expectation value can be measured directly with the PTAM operator instead, but the POAM operator must first be converted to a matrix (22)where is the zeroth Pauli spin matrix (2 × 2 identity matrix). With this redefined POAM operator, the PTAM operator becomes (23a)Thus, (23b)The choice of measuring Ĵ_{Z} using separate and operators or the combined operator depends on the application.
4.2. Imperfect instrument
An instrument with nonzero instrumental PSAM, subject to the equations of Sect. 3.3, modifies the expectation values derived in Sect. 4.1. For the sake of simplicity, I assume that the circular telescope aperture is uniformly unaberrated with a nonzero instrumental PSAM, which means that and (24c)where R_{tel} is the telescope radius in units of wavelength, , and is the Bessel function of the first kind of order one. The quantity in the curly braces approaches the Dirac delta function when R_{tel} → ∞. A perfect instrument implies that , where η is a complex constant (0 < η ≤ 1). Conversely, when the instrument mixes the PSAM components.
To keep the notation consistent with Sect. 4.1, I let , , and . The PTAM expectation value for this imperfect instrument is (25a)where is the operator that includes PSAM mixing effects from the imperfect instrument, (25b)is the integrated intensity through the imperfect instrument, and (25c)is a 2 × 2 matrix. Because of the linearity of the POAM and PSAM operators, the mixed PTAM operator becomes (26a)where (26b)is the POAM operator including the imperfect instrument, (26c)is the PSAM operator including the imperfect instrument, and (26d)is another 2 × 2 matrix. When : η^{2}1, , , , , , , and Ŝ_{Z}. The η factor does not modify the expectation values because they are normalized quantities.
5. Simple example
In Sect. 1, I show that the most general unfactored PTAM electric field has PSAM states with different POAM expansions. In this section, I demonstrate how the measured POAM expectation value can be affected by source and instrumental PSAM using the simplest unfactored PTAM electric field (27)The PTAM expectation value ultimately depends on the behavior of three intensitybased quantities, and (28c)where δ_{m,n} is the Kronecker delta function. They can be rearranged to become more familiar quantities, namely the Stokes parameters I =I_{A,A} + I_{B,B}, Q = 2Re{I_{A,B}}δ_{m,n}, U = 2Im{I_{A,B}}δ_{m,n}, and V = I_{A,A} − I_{B,B}. This electric field contains linear polarization only when the PSAM states are temporally correlated have the same POAM state, or m = n.
Using the definitions in Sects. 2 and 4, the PTAM expectation value becomes (29)Note that p_{A,A} and p_{B,B} are part of both the POAM and PSAM expectation values. For a purely unpolarized and/or linearly polarized source p_{A,A} = and . Also, Ĵ_{Z} = ^{(}m + 1^{)}ħ for a purely righthanded circularly polarized source (p_{A,A} = 1 and p_{B,B} = 0) and Ĵ_{Z} = ^{(}n − 1^{)}ħ for a purely lefthanded circularly polarized source (p_{A,A} = 0 and p_{B,B} = 1). If m= n (factored electric field), those PTAM expectation values become Ĵ_{Z} = mħ (PTAM expectation value is independent of PSAM expectation value), ^{(}m + 1^{)}ħ (PTAM expectation value is POAM expectation value plus RCP PSAM expectation value), and ^{(}m − 1^{)}ħ (PTAM expectation value is POAM expectation value minus LCP PSAM expectation value), respectively.
An instrument with nonzero instrumental PSAM modifies the result of Eq. (29). Using a Jones matrix in the form of Eq. (24a), the electric field becomes (30)When (D^{A,B} = D^{B,A} ≠ 0), the instrument mixes the PSAM components. The integrated intensity of this electric field (cf. Eq. (25b)) can be rewritten in terms of the source Stokes parameters (31a)where and (31e)are instrumentdependent coefficients, and q = Q/I, u = U/I, and v= V/I are the normalized source Stokes parameters. These equations are interesting, because 1) source circular polarization contributes to I′; and 2) if m = n (unfactored electric field) the linear PSAM contributes to I′. Note that the system gain ℳ→ ℳ(q,u,v), i.e., it depends on both the instrumental and source PSAM.
Using Eqs. (30) and (31a) and the mathematics of Sect. 4.2, the POAM and PSAM expectation values including instrumental PSAM are (32a)and (32b)where and (32f)are other instrumentdependent coefficients, and is the “transitional probability”. Equation (32b) shows that instrumental PSAM changes the PSAM expectation value. Elias (2008) showed that instrumental POAM changes the POAM expectation value. These results are not unexpected and not particularly exciting. On the other hand, changes in the POAM expectation value due to instrumental and source PSAM deserve further analysis. I call this effect “PSAMModified POAM Measurement”, or SMOM.
After analysing the complete set of use cases (unpolarized, linearly polarized, circularly polarized, elliptically polarized source and instrument polarization; m ≠ n or m = n), I found that (33)where m^{I,V} = M^{I,V}/M^{I,I} is the normalized circular PSAM gain. This equation completely describes the conditions required for SMOM in this simplified example.
SMOM is possible only when m ≠ n, i.e., for unfactored PTAM electric fields. If the PTAM electric field is factored (m = n), on the other hand, (34)does not lead to SMOM. SMOM only occurs in the presence of instrumental circular PSAM, or m^{I,V} ≠ 0, because only it can mix the different PSAM states leading to modified POAM expectation values. No PSAM or partial/full linear PSAM corresponds to v = 0, which leads to the maximum for a given instrumental circular PSAM. The 0 <v< 1 cases correspond to circular source PSAM plus a combination of unpolarized and/or linear PSAM. Increasing v decreases the . The v = + 1 and v = −1 cases correspond to (35a)and (35b)respectively. Because only one PSAM component is nonzero, no SMOM is possible. These electric fields are similar to Eq. (34) because the POAM exponential can also be factored outside the PSAM vector.
6. Conclusions
I present the most general “unfactored” PTAM electric field form, where each PSAM component has its own POAM expansion. It is slightly more general than the more commonly invoked “factored” PTAM electric field form where the PSAM and POAM components are separable. I then combine the POAM and PSAM calculi to obtain the PTAM calculi. Apart from the vectors, matrices, dot products, and direct products, the PTAM and POAM calculi appear superficially identical. I derive the PTAM operator and expectation value in terms of POAM/PSAM operators and expectation values for systems with and without instrumental PSAM. Last, I prove using a simple example that POAM measurements of sources with unfactored PTAM electric fields passing through instrumental circular PSAM yield systematic POAM measurement errors.
Acknowledgments
The author acknowledges partial support from the OPTICON Fizeau Program, the Heidelberger Landessternwarte, and the Air Force Research Laboratory (clearance number 377ABW20130902). The author also acknowledges the efforts of the anonymous referee whose comments greatly improved this work.
References
 Barakat, R. 1981, Opt. Commun., 38, 159 [NASA ADS] [CrossRef] [Google Scholar]
 Elias II, N. M. 2008, A&A, 492, 883 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Elias II, N. M. 2012, A&A, 541, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jones, R. C. 1941, J. Opt. Soc. Am., 31, 488 [Google Scholar]
 Mueller, H. 1948, J. Opt. Soc. Am., 38, 661 [Google Scholar]
 Schmeider, R.W. 1969, J. Opt. Soc. Am., 59, 297 [NASA ADS] [CrossRef] [Google Scholar]
 Stokes, G. 1852, Trans. Cambridge Phil. Soc., 9, 399 [Google Scholar]
All Tables
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.