General vectorial decomposition of electromagnetic fields with application to propagation-invariant and rotating fields

doi:10.1364/OE.10.000949

« Show journal navigation

General vectorial decomposition of electromagnetic fields with application to propagation-invariant and rotating fields

Pertti Paakkonen, Jani Tervo, Pasi Vahimaa, Jari Turunen, and Franco Gori »View Author Affiliations

Optics Express, Vol. 10, Issue 18, pp. 949-959 (2002)
http://dx.doi.org/10.1364/OE.10.000949

View Full Text Article

Acrobat PDF (526 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (21)
Cited By
Figures (4)
Multimedia (2)
Metrics

Abstract

A novel decomposition of the transversal part of the electric field vector of a general non-paraxial electromagnetic field is presented, which is an extension of the radial/aximuthal decomposition and is known as γζ decomposition. Purely γ and ζ polarized fields are examined and the decomposition is applied to propagation-invariant, rotating, and self-imaging electromagnetic fields. An experimental example on the effect of state of polarization in the propagation characteristics of the field: its is shown that a simple modification of the polarization conditions of the angular spectrum converts a self-imaging field into a propagation-invariant field.

[Optical Society of America ]

1 Introduction

Self-imaging [1

W. D. Montgomery, “Self-imaging objects of infinite aperture,” J. Opt. Soc. Am. 57, 772–778 (1967). [CrossRef]

, 2

K. Patorski, “The self-imaging phenomenon and its applications,” in Progress in Optics Vol. XXVII, E. Wolf, ed. (Elsevier, Amsterdam, 1989), Chap. 1.

], propagation-invariance [3

J. Durnin, “Exact solutions for nondiffracting beams. I. The scalar theory,” J. Opt. Soc. Am. A 4, 651–654 (1987). [CrossRef]

, 4

J. Durnin, J. J. Miceli Jr., and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987). [CrossRef] [PubMed]

], and rotation [5

Y. Y. Schechner, R. Piestun, and J. Shamir, “Wave propagating with rotating intensity distributions,” Phys. Rev. E 54, R50–R53 (1996). [CrossRef]

, 6

S. Chávez-Cerda, G. S. McDonald, and G. H. S. New, “Nondiffracting Beams: travelling, standing, rotating and spiral waves,” Opt. Commun. 123, 225–233 (1996). [CrossRef]

, 7

C. Paterson and R. Smith, “Higher-order Bessel waves produced by axicon-type computer-generated holograms,” Opt. Commun. 124, 121–130 (1996). [CrossRef]

] of the intensity distribution of a scalar field are closely related phenomena that have attracted considerable recent interest. In the scalar case the angular spectrum of plane waves associated with the field must be confined on one ring, while in the case of rotating and self-imaging fields two or more rings (known as Montgomery’s rings) are involved.

In the case of propagation-invariant fields, the scalar analysis has been extended to the electromagnetic case [8

S. R. Mishra, “A vector wave analysis of a Bessel beam,” Opt. Commun. 85, 159–161 (1991). [CrossRef]

, 9

J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]

, 10

Z. Bouchal and M. Olivík, “Non-diffractive vector Bessel beams,” J. Mod. Opt. 42, 1555–1566 (1995). [CrossRef]

, 11

Z. Bouchal, R. Horák, and J. Wagner, “Propagation-invariant electromagnetic fields,” J. Mod. Opt. 43, 1905–1920 (1996). [CrossRef]

], which is now also rather well understood. Some important new features, such as the broken rotational symmetry of pure Bessel field modes in the non-paraxial domain and the fact that in certain cases two Montgomery’s rings may be involved [9

J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]

], have been uncovered. However, in the case of rotating fields the electromagnetic analysis differs more fundamentally from the scalar analysis [12

J. Tervo and J. Turunen, “Rotating scale-invariant electromagnetic fields,” Opt. Express 9, 9–15 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9. [CrossRef] [PubMed]

] because the longitudinal field component does not usually rotate when the transverse cartesian components do so. In paraxial geometries, though, the longitudinal component has an insignificant contribution to the total energy density and the electromagnetic extension of the scalar analysis is accurate [13

P. Pääkkönen, J. Lautanen, M. Honkanen, M. Kuittinen, J. Turunen, S. N. Khonina, V. V. Kotl-yar, V. A. Soifer, and A. T. Friberg, “Rotating optical fields: experimental demonstration with diffractive optics,” J. Mod. Opt. 46, 2355–2369 (1998). [CrossRef]

The purpose of this paper is extend a novel decomposition of the electric field into two orthogonal components [14

F. Gori, “Polarization basis for vortex beams,” J. Opt. Soc. Am. A 18, 1612–1617 (2001). [CrossRef]

], which we call the γζ decomposition, to arbitrary non-paraxial electromagnetic fields, and to investigate the implications of this decomposition in the theory of propagation-invariant, rotating, and self-imaging fields.

The general theory of the γζ decomposition is presented in Sect. 2 and applied in Sect. 3 to purely γ and ζ polarized electromagnetic fields. In all cases exact propagation formulas are derived, which are based on the angular spectrum representation of the electromagnetic field. In Sect. 4 the γζ decomposition is applied to propagation-invariant fields, including the paraxial case. Some of the results are illustrated experimentally in Sect. 5 and an analysis of rotating fields is provided in Sect. 6.

2 The γζ decomposition of the transverse electric field

The propagation of a time-harmonic electromagnetic field in free space is conveniently described by means of the angular spectrum representation [15

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

]. Let us express both the position vector r = (x, y, z) and the wave vector k = (k_x , k_y , k_z ) in circular cylindrical coordinates by using the relations

{\begin{matrix} x = ρ \cos ϕ \\ y = ρ \sin ϕ \end{matrix}

(1)

and

{\begin{matrix} k_{x} = α \cos ψ \\ k_{y} = α \sin ψ \\ k_{z} = β \end{matrix}

(2)

The expression for the electric field now takes the form

E (ρ, ϕ, z) = \int_{0}^{2 π} \int_{0}^{\infty} A (α, ψ) \exp [i αρ \cos (ϕ - ψ) + i βz] α d α d ψ,

(3)

where

β = {\begin{matrix} \sqrt{k^{2} - α^{2}} if α \leq k \\ i \sqrt{α^{2} - k^{2}} otherwise \end{matrix}

(4)

and k = | k | is the wave number. The angular spectrum A (α, Ψ) is obtained by Fourier-inversion at z = 0:

A (α, ψ) = \frac{1}{{(2 π)}^{2}} \int_{0}^{2 π} \int_{0}^{\infty} E (ρ, ϕ, 0) \exp [- i αρ \cos (ϕ - ψ)] ρ d ρ d ϕ .

(5)

Generally, two of the three cartesian components of E may be assumed to be independent. The third is obtained from the Maxwell’s divergence equation

\nabla \cdot E (ρ, ϕ, z) = 0 .

(6)

Let us next decompose the transversal part of the field E_xx̂ + E_yŷ into two components by the operation

{\begin{matrix} E_{γ}^{(q)} (ρ, ϕ, z) = \cos (qϕ) E_{x} (ρ, ϕ, z) + \sin (qϕ) E_{y} (ρ, ϕ, z) \\ E_{ζ}^{(q)} (ρ, ϕ, z) = - \sin (qϕ) E_{x} (ρ, ϕ, z) + \cos (qϕ) E_{y} (ρ, ϕ, z), \end{matrix}

(7)

where q is an integer. The longitudinal component of the field vector is retained as original. Equation (7) clearly defines also a new pair of unit vectors, denoted by

\hat{γ}

and

\hat{ζ}

, which are orthogonal in every point of space. These unit vectors remain invariant in the radial direction, but not in the azimuthal direction, in which they rotate with a constant rate along the ϕ coordinate.

In the following considerations we will call this basis the γζ basis. The (vectorial) components

E_{γ}^{(q)}

and E _ζ^(q) are equivalent to the basis vectors discussed recently in Ref. [14

F. Gori, “Polarization basis for vortex beams,” J. Opt. Soc. Am. A 18, 1612–1617 (2001). [CrossRef]

], although the notations used here are different. The rotation of the basis vectors is clearly dependent on the value of q, since they experience a rotation of q2π within 0 ≤ ϕ ≤ 2π. In the case q = 1 the basis vectors are rotationally symmetric and are customarily called the radial and azimuthal components denoted by E_ρ and E_ϕ , respectively. For values q > 1 the basis vectors experience several full rotations in the counterclockwise direction. On the other hand, for negative values of q the direction of rotation of the basis vectors is reversed. For example, if q = -2 the basis vectors experience a rotation of - 4π radians in the clockwise direction when 0 ≤ ϕ ≤ 2π. In the case q = 0 the basis vectors represent unmodified cartesian vectors. An example illustrating the directions of the basis vectors with values q = 1 and q = -2 is given in Fig. 1.

Fig. 1. Directions of the transversal basis vectors in the cases q = 1 and q = - 2.

Combining Eqs. (3) and (7) yields general propagation formulas for the transverse part of the electric field in terms of the γ and ζ components:

E_{γ}^{(q)} (ρ, ϕ, z) = \sum_{m = - \infty}^{\infty} \exp (i mϕ) \int_{0}^{\infty} J_{m} (αρ) \exp (i βz)

\times [a_{m}^{x} (α) \cos (qϕ) + a_{m}^{y} (α) \sin (qϕ)] dα,

(8)

and

E_{ζ}^{(q)} (ρ, ϕ, z) = \sum_{m = - \infty}^{\infty} \exp (i mϕ) \int_{0}^{\infty} J_{m} (αρ) \exp (i βz)

\times [{- a}_{m}^{x} (α) \sin (qϕ) + a_{m}^{y} (α) \cos (qϕ)] dα .

(9)

In the derivation of Eqs. (8

S. R. Mishra, “A vector wave analysis of a Bessel beam,” Opt. Commun. 85, 159–161 (1991). [CrossRef]

) and (9

J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]

) we have made use of the Jacobi–Anger expansion [16

G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists (Academic Press, New York, 2001), p. 681.

]

\exp (iϑ \cos τ) = \sum_{m = - \infty}^{\infty} i^{m} J_{m} (ϑ) \exp (i mτ),

(10)

where J_m denotes the Bessel function of the first kind and of order m, and defined the functions

a_{m}^{j}

(α) as

a_{m}^{j} (α) = i^{m} α \int_{0}^{2 π} A_{j} (α, ψ) \exp (- i mψ) d ψ,

(11)

where j = x or y.

3 Purely γ- or ζ-polarized fields

In this Section we derive general expressions for fields that have only either a γ component or a ζ component. These results represent a generalization of those given in Refs. [17

R. H. Jordan and D. G. Hall, “Highly directional surface emission from concentric-circle gratings on planar optical waveguides: the field-expansion method,” J. Opt. Soc. Am. A 12, 84–94 (1995). [CrossRef]

] and [18

J. Tervo, P. Vahimaa, and J. Turunen, “On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,” J. Mod. Opt. 49, 1537–1543 (2002). [CrossRef]

], where the field was assumed to be either radially or azimuthally polarized. In addition, we assume no radial symmetry or propagation-invariance of the field. In the following, we assume that q ≠ 0, since in the case q = 0 the field representation is clearly given by Eq. (3).

Let us first consider the case in which the ζ component of the field vanishes. By inserting this requirement into Eq. (9) and rearranging the indices in the summation, we obtain

\sum_{m = - \infty}^{\infty} \exp (i mϕ) \int_{0}^{\infty} {J_{m - q} (αρ) [i a_{m - q}^{x} (α) - a_{m - q}^{y} (α)]

+ J_{m + q} (αρ) [- i a_{m + q}^{x} (α) + a_{m + q}^{y} (α)]} \exp (i βz) d α = 0 .

(12)

Using the uniqueness of the Fourier-series representation and the fact that Eq. (12) must hold for all values of ρ and z, we arrive at the requirement that the integrand must vanish identically:

J_{m - q} (αρ) [i a_{m - q}^{x} (α) + a_{m - q}^{y} (α)] + J_{m + q} (αρ) [- i a_{m + q}^{x} (α) + a_{m + q}^{y} (α)] \equiv 0,

(13)

must hold for all m. However, because the functions J _m-q and J _m+q are linearly independent for all m ≠ 0, Eq. (13) implies that only the functions

a_{m}^{j}

(α) with m = ±q may have non-zero values. In addition, these functions are connected by the relations

{\begin{matrix} a_{- q}^{y} (α) = i {(- 1)}^{q} a_{q}^{x} (α) \\ a_{q}^{y} (α) = - i a_{q}^{x} (α) . \\ a_{- q}^{x} (α) = {(- 1)}^{q} a_{q}^{x} (α) \end{matrix}

(14)

Inserting these into Eq. (8) we obtain the propagation formula

E_{γ}^{(q)} (ρ, ϕ, z) = E_{γ}^{(q)} (ρ, z) = \int_{0}^{\infty} f_{q} (α) J_{q} (αρ) \exp (i βz) d α,

(15)

where

f_{q} (α) = 2 a_{q}^{x} (α) .

(16)

Thus, the transverse vectorial component of the field may be handled as a rotationally symmetric field, although the propagation integral differs radically from Eq. (5). However, the entire field symmetric only in the case q = 1. This is seen also by inserting Eqs. (3) and (14) into Eq. (6), which yields an propagation formula for the longitudinal component of the field:

E_{z}^{(q)} (ρ, ϕ, z) = i \cos [(q - 1) ϕ] \int_{0}^{\infty} \frac{α}{β} f_{q} (α) J_{q - 1} (αρ) \exp (i βz) d α .

(17)

The absolute value of

E_{z}^{(q)}

(ρ, ϕ, z) clearly depends on the ϕ-coordinate for all q = 1.

Derivation of a purely ζ polarized field is performed essentially identically, by assuming that E _γ(ρ, ϕ, z) ≡ 0 and then repeating similar steps as above. The results are

E_{ζ}^{(q)} (ρ, ϕ, z) = E_{ζ}^{(q)} (ρ, z) = \int_{0}^{\infty} g_{q} (α) J_{q} (αρ) \exp (i βz) d α

(18)

and

E_{z}^{(q)} (ρ, ϕ, z) = i \sin [(q - 1) ϕ] \int_{0}^{\infty} \frac{α}{β} g_{q} (α) J_{q - 1} (αρ) \exp (i βz) d α,

(19)

where

g_{q} (α) = - i f_{q} (α) .

(20)

Equations (15) and (18) state, e.g., that a field which is either radially or azimuthally polarized in one transversal plane can retain its polarization state only if it is rotationally symmetric. In addition, this field must be a superposition of J ₁ subfields with different radii. This is a generalization of the result given by Jordan and Hall [17

], who derived Eq. (18) in the case q = 1 assuming that the azimuthally polarized field is rotationally symmetric. The derivation introduced above shows that the rotationally symmetric case examined by Jordan and Hall is the only possible one, which explains the results reported recently by Lapucci and Ciofini [19

A. Lapucci and M. Ciofini, “Polarization state modifications in the propagation of high azimuthal order annular beams,” Opt. Express 9, 603–609 (2001),http: //www. opticsexpress.org/abstract.cfm?URI=OPEX-9-12-603. [CrossRef] [PubMed]

The results derived above are valid only if either the γ or the ζ component of the field vanishes. In all other cases both components are superpositions of Bessel functions with different orders. However, the linearity of Maxwell’s equations implies that the results hold also for all linear combinations of

E_{γ}^{(q)}

(ρ, ϕ, z) and

E_{ζ}^{(q)}

(ρ, ϕ, z).

4 Paraxial propagation-invariant electromagnetic fields

Probably the most natural definition of propagation-invariance is the condition that the time-averaged energy density of the field 〈w(r,t)〉 remains exactly the same in every transversal plane z = constant [9

J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]

, 18

J. Tervo, P. Vahimaa, and J. Turunen, “On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,” J. Mod. Opt. 49, 1537–1543 (2002). [CrossRef]

], i.e.,

〈 w (ρ, ϕ, z + Δ z, t) 〉 = 〈 w (ρ, ϕ, z, t) 〉

(21)

for all values of Δz. The condition (21) may also be expressed in the form

\frac{\partial}{\partial z} = 〈 w (ρ, ϕ, z, t) 〉 \equiv 0,

(22)

which is sometimes more suitable.

By comparison of Eqs. (3), (6), and (21) we see that any electromagnetic field with propagation-invariant scalar components is propagation-invariant. However, as pointed out by Turunen and Friberg [9], that condition is sufficient but not necessary when the electromagnetic nature of the field is taken into account. For example, propagation-invariant electromagnetic fields with non-propagation-invariant x- and y-components may be constructed rather easily [9

J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]

, 18

J. Tervo, P. Vahimaa, and J. Turunen, “On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,” J. Mod. Opt. 49, 1537–1543 (2002). [CrossRef]

]. These fields may be obtained by assuming that both radial and azimuthal components of the field are propagation-invariant, but not necessarily localized on the same Montgomery’s ring. Because the azimuthally polarized field does not have a z-component, the total energy density of this kind of field must remain constant upon propagation [18

J. Tervo, P. Vahimaa, and J. Turunen, “On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,” J. Mod. Opt. 49, 1537–1543 (2002). [CrossRef]

One may expect that this phenomenon may be extended to other cases in which one of the two component-fields does not have a z-component. However, as can be easily verified by inserting Eq. (6) into (3), the condition that the z-component disappears leads immediately to an azimuthally polarized angular spectrum, i.e., A_x (ρ, Ψ) cos Ψ + A_y (ρ, Ψ) sin Ψ = 0 [20

J. Tervo and J. Turunen, “Self-imaging of electromagnetic fields,” Opt. Express 9, 622–630 (2001), http: //www. opticsexpress.org/abstract.cfm?URI=OPEX-9-12-622. [CrossRef] [PubMed]

]. It is well known that, for such fields, the direction of the electric field vector changes in the radial direction, except for the purely azimuthally polarized field, and hence the orthogonality of the transverse components is possible only if both component-fields are confined to the same Montgomery’s ring [18

J. Tervo, P. Vahimaa, and J. Turunen, “On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,” J. Mod. Opt. 49, 1537–1543 (2002). [CrossRef]

, 21

J. Tervo and J. Turunen, “Generation of vectorial propagation-invariant fields by polarization-grating axicons,” Opt. Commun. 192, 13–18 (2001). [CrossRef]

Let us next assume that both γ and ζ components of the field are propagation-invariant, i.e.,

E_{γ}^{(q)} (ρ, ϕ, z + Δ z) = \exp [iξ (ρ, ϕ, Δ z)] E_{γ}^{(q)} (ρ, ϕ, z)

(23)

and

E_{ζ}^{(q)} (ρ, ϕ, z + Δ z) = \exp [iν (ρ, ϕ, Δ z)] E_{ζ}^{(q)} (ρ, ϕ, z)

(24)

where ξ and ν are as-yet arbitrary real functions. In view of Eqs. (15) and (18), the field expressions are of the form

E_{γ}^{(q)} (ρ, z) = c_{γ} J_{q} (α_{γ} ρ) \exp (i β_{γ} z)

(25)

and

E_{ζ}^{(q)} (ρ, z) = c_{ζ} J_{q} (α_{ζ} ρ) \exp (i β_{ζ} z),

(26)

where c _γ and c _ζ are arbitrary complex constants and α_Γ and α_ζ denote the constant α ₀ for the γ and ζ components, respectively.

The z-component of the field is obtained straightforwardly by using Eqs. (17) and (19) and it takes the form

E_{z} (ρ, ϕ, z) = i c_{γ} \frac{α_{γ}}{β_{γ}} \cos [(q - 1) ϕ] J_{q - 1} (α_{γ} ρ) \exp (i β_{γ} z)

+ i c_{ζ} \frac{α_{ζ}}{β_{ζ}} \sin [(q - 1) ϕ] J_{q - 1} (α_{ζ} ρ) \exp (i β_{ζ} z) .

(27)

By squaring the absolute value of the field in Eq. (27), it is immediately seen that when α_γ ≠ α_ζ the z-component is modified upon propagation in almost all cases, although both γ and ζ components are propagation-invariant. Thus, the electric energy density w_e ∝ ∥E(ρ, ϕ, z)∥² (and also the magnetic energy density) is not generally propagation-invariant, but instead self-imaging. Although there exist several special cases, in which the modifications of the electric and magnetic energy densities cancel each other [9

J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]

, 20

J. Tervo and J. Turunen, “Self-imaging of electromagnetic fields,” Opt. Express 9, 622–630 (2001), http: //www. opticsexpress.org/abstract.cfm?URI=OPEX-9-12-622. [CrossRef] [PubMed]

], in the case considered here that does not happen, as can be easily verified by a straightforward calculation.

The only known cases in which the field confined on two Montgomery’s rings may be propagation-invariant are the case in which either α_γ or α_ζ is equal to zero [9

J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]

] and the case q = 1, i.e., when E _γ and E_ζ represent radial and azimuthal components, respectively [18

J. Tervo, P. Vahimaa, and J. Turunen, “On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,” J. Mod. Opt. 49, 1537–1543 (2002). [CrossRef]

]. In the paraxial limit α_j /β_j → 0 the contribution to the energy density from the z-component of the field may be neglected and the field becomes propagation-invariant.

5 Experimental results on propagation-invariance

A simple way to demonstrate finite-aperture approximations of propagation-invariant fields is to use a ring aperture, whose far-field diffraction pattern is a Bessel field. The radius of the ring corresponds to the radius of the Montgomery’s ring in the Fourier-plane.

In the case of an aperture consisting of two concentric rings, a self-imaging field is generally produced [1

W. D. Montgomery, “Self-imaging objects of infinite aperture,” J. Opt. Soc. Am. 57, 772–778 (1967). [CrossRef]

]. However, as we saw in the previous Section, the field may be almost propagation-invariant if the transversal parts of the sub-fields confined on different rings are orthogonal. This situation may be achieved, for example, by using the setup illustrated in Fig. 2:An annular half-wave plate is used to alter the polarization state of the sub-field confined to the outer ring. A Fourier-transform of the field is then performed by using a simple refractive lens.

The orientation of the optical axis of the half-wave plate with respect to the direction of linear polarization of the incident field, denoted by θ, directly affects the polarization state of the field. If a linearly polarized input field is used, the angle θ = π/4 converts an x-polarized wave into a y-polarized wave and vice versa. On the other hand, if the optical axis is parallel (or perpendicular) to the direction of polarization, the polarization state is not modified. These cases are illustrated in Fig. 3.

Fig. 2. The principle of the experimental setup. Here a and b denote an annular half-wave plate and a ring aperture, respectively, and c is a thin lens.

We made the experiments by using a setup illustrated in Fig. 2. The used light source was a Helium–Neon laser with λ = 632.8 nm, whose phase-front was made approximately planar by using a lens-system (not shown in the Figure). The diameters of the rings in the aperture were 6 and 10 mm, whereas the (inner) diameter of the annular half-wave plate was 8 mm. The focal length of the lens was f = 1000 mm.

Fig. 3. (1.21 MB) An animation of the experimental fields within one self-imaging distance z_T ≈ 20 mm. The field in the left-hand side is obtained with parallel polarization states and the field in the right-hand side with orthogonal polarization states.

By examining Fig. 3, we notice that the field produced by orthogonally polarized rings is propagation-invariant, whereas the parallel case leads to an ordinary self-imaging field, as expected.

6 Rotating intensity distributions

The orthogonality of the γ and ζ polarized fields may be used to produce not only propagation-invariant fields, but rotating fields as well. The condition for rotation of an electromagnetic field can be defined as [12

J. Tervo and J. Turunen, “Rotating scale-invariant electromagnetic fields,” Opt. Express 9, 9–15 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9. [CrossRef] [PubMed]

]

〈 w (ρ, ϕ + η Δ z, z + Δ z, t) 〉 = 〈 w (ρ, ϕ, z, t) 〉

(28)

for all Δz. Here η is a constant that defines both the direction of propagation and the self-imaging distance. Contrary to the case of propagation-invariant fields, the energy-density distributions of electromagnetic extensions of rotating scalar fields do not generally rotate [12

J. Tervo and J. Turunen, “Rotating scale-invariant electromagnetic fields,” Opt. Express 9, 9–15 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9. [CrossRef] [PubMed]

]. However, in the paraxial domain the scalar theory predicts the behavior of the transverse components of the field accurately, and hence in many practical applications the prediction of the scalar theory coincides with the experimental results [13

It is clear that scalar fields are not the only solutions whose intensity distributions rotate in the paraxial domain, for it is possible to combine non-rotating scalar x- and y- components of the field to form rotating intensity distributions. However, the task of finding that kind of fields directly from Eq. (28) appears to be a very demanding task and hence some other method must be used. In the following, we use the γζ decomposition to find a class of such solutions.

In view of Eq. (28) it is clear that the intensity distribution of the paraxial field which fulfills the conditions

{\begin{matrix} E_{γ}^{(q)} (ρ, ϕ + η Δ z, z + Δ z) = \exp [i ξ (ρ, ϕ, Δ z)] E_{γ}^{(q)} (ρ, ϕ, z) \\ E_{ζ}^{(q)} (ρ, ϕ + η Δ z, z + Δ z) = \exp [i ξ (ρ, ϕ, Δ z)] E_{ζ}^{(q)} (ρ, ϕ, z) \end{matrix},

(29)

where ξ(ρ, ϕ, Δz) is an arbitrary real function, rotates. Inserting this requirement into Eqs. (8) and (9) yields, similarly to the case of exactly rotating fields [12

J. Tervo and J. Turunen, “Rotating scale-invariant electromagnetic fields,” Opt. Express 9, 9–15 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9. [CrossRef] [PubMed]

], the connections

a_{m}^{y} (α) = i s a_{m}^{x} (α)

(30)

and

β_{ms} = β_{0} - (m + sq) η,

(31)

where s is a sign function which can assume values s = -1 or s = +1 for each m, and β ₀ is a constant. Naturally, in the cases q = 0 or q = 1, we obtain the rotation conditions for scalar or electromagnetic fields, respectively [12

J. Tervo and J. Turunen, “Rotating scale-invariant electromagnetic fields,” Opt. Express 9, 9–15 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9. [CrossRef] [PubMed]

By inserting Eqs. (30) and (31) into Eqs. (8) and (9), we obtain the following expressions for the transversal components:

E_{γ}^{(q)} (ρ, ϕ, z) = \sum_{m \in 𝓜} \sum_{s = - 1,1} a_{ms} J_{m} (α_{ms} ρ) \exp {i [(m + sq) ϕ + β_{ms} z]}

(32)

and

E_{ζ}^{(q)} (ρ, ϕ, z) = \sum_{m \in 𝓜} \sum_{s = - 1,1} i s a_{ms} J_{m} (α_{ms} ρ) \exp {i [(m + sq) ϕ + β_{ms} z]},

(33)

where 𝓜 denotes the set of values of m leading to real values of β_ms , in view of Eqs. (4) and (31). The coefficients a_ms are defined by a_ms =

a_{m}^{x}

(α_ms)/α_ms . An example of the field satisfying Eqs. (30) and (31) is illustrated in Fig. 4. It can be clearly seen that the x-component of the field is not rotating. However, the γ-component, as well as the ζ-component (which is not shown), do rotate. This means that the contribution to the intensity form the transversal part of the field rotates.

Table 1. The parameters assumed in Fig. 4. The constant β_C must be less than k but is otherwise arbitrary.

m	β_m	n	a_mn
0	β_c /2	1	i
-2	β_c /2	-1	1
-1	β_c	1	-i
-3	β_c	-1	1

Fig. 4. (1.38 MB) An animation of the intensity distributions of the central parts of the x- and γ-components and the total electric energy density as functions of the z-coordinate, calculated with the parameters given in Table 1. The movie is best appreciated if viewed repeatedly.

Let us next consider the case q = 1, which provides also rotating z-components, of course assuming that both conditions (30) and (31) are met. Now the z-component of the electric field takes the form [12

J. Tervo and J. Turunen, “Rotating scale-invariant electromagnetic fields,” Opt. Express 9, 9–15 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9. [CrossRef] [PubMed]

]

E_{z} (ρ, ϕ, z) = \sum_{m \in 𝓜} \sum_{s = - 1,1} s a_{ms} \frac{α_{ms}}{i β_{ms}} J_{m + s} (α_{ms} ρ) \exp {i [(m + s) ϕ + β_{ms} z]} .

(34)

If we now retain only one value of m in Eqs. (32) and (33), we find that the transversal part of the field takes the form

E_{⊥} (ρ, ϕ, z) = \exp (i β_{m, 1} z) a_{m, 1} J_{m} (α_{m, 1} ρ) \exp (i m ϕ) [\begin{matrix} 1 \\ i \end{matrix}]

+ \exp (i β_{m, - 1} z) a_{m, - 1} J_{m} (α_{m, - 1} ρ) \exp (i m ϕ) [\begin{matrix} 1 \\ - i \end{matrix}],

(35)

where the column vectors are the usual Jones vectors. Thus, the left- and right-handed components are of propagation-invariant form. Hence, if this kind of a field is paraxial, the intensity distribution is approximately propagation-invariant, although the field is exactly rotating. Such a field may be easily demonstrated by using the setup discussed in Section 6, but by using a circularly polarized input field. In that case, the half-wave plate inverts the handedness of the field and hence a field of the form of Eq. (35) is produced.

7 Conclusions

In conclusion, we have introduced some new aspects into the theory of propagation-invariant, rotating, and self-imaging electromagnetic fields. Starting from a general and exact γζ decomposition of an arbitrary non-paraxial field we have considered the propagation and the properties of purely γ and ζ polarized fields, or which radially and azimuthally polarized fields are important special cases. It was then shown that if both the γ and ζ components are propagation-invariant, so is the electric energy density even though the two components are confined on Montgomery’s rings with different radii. A generalized class of paraxial rotating electromagnetic fields was also introduced, and an experimental demonstration illustrating a dramatic effect of the polarization properties of the angular spectrum into the propagation characteristics of electromagnetic fields: an apparently minor modification of the polarization properties converts a self-imaging field into a propagation-invariant field.

Acknowledgements

Pertti Pääakkönen thanks the Emil Aaltonen foundation for a generous grant. Jani Tervo thanks the Nokia Foundation and the Finnish Cultural Foundation for partial financial support.

References and links

1.	W. D. Montgomery, “Self-imaging objects of infinite aperture,” J. Opt. Soc. Am. 57, 772–778 (1967). [CrossRef]
2.	K. Patorski, “The self-imaging phenomenon and its applications,” in Progress in Optics Vol. XXVII, E. Wolf, ed. (Elsevier, Amsterdam, 1989), Chap. 1.
3.	J. Durnin, “Exact solutions for nondiffracting beams. I. The scalar theory,” J. Opt. Soc. Am. A 4, 651–654 (1987). [CrossRef]
4.	J. Durnin, J. J. Miceli Jr., and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987). [CrossRef] [PubMed]
5.	Y. Y. Schechner, R. Piestun, and J. Shamir, “Wave propagating with rotating intensity distributions,” Phys. Rev. E 54, R50–R53 (1996). [CrossRef]
6.	S. Chávez-Cerda, G. S. McDonald, and G. H. S. New, “Nondiffracting Beams: travelling, standing, rotating and spiral waves,” Opt. Commun. 123, 225–233 (1996). [CrossRef]
7.	C. Paterson and R. Smith, “Higher-order Bessel waves produced by axicon-type computer-generated holograms,” Opt. Commun. 124, 121–130 (1996). [CrossRef]
8.	S. R. Mishra, “A vector wave analysis of a Bessel beam,” Opt. Commun. 85, 159–161 (1991). [CrossRef]
9.	J. Turunen and A. T. Friberg, “Self-imaging and propagation-invariance in electromagnetic fields,” Pure Appl. Opt. 2, 51–60 (1993). [CrossRef]
10.	Z. Bouchal and M. Olivík, “Non-diffractive vector Bessel beams,” J. Mod. Opt. 42, 1555–1566 (1995). [CrossRef]
11.	Z. Bouchal, R. Horák, and J. Wagner, “Propagation-invariant electromagnetic fields,” J. Mod. Opt. 43, 1905–1920 (1996). [CrossRef]
12.	J. Tervo and J. Turunen, “Rotating scale-invariant electromagnetic fields,” Opt. Express 9, 9–15 (2001), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9. [CrossRef] [PubMed]
13.	P. Pääkkönen, J. Lautanen, M. Honkanen, M. Kuittinen, J. Turunen, S. N. Khonina, V. V. Kotl-yar, V. A. Soifer, and A. T. Friberg, “Rotating optical fields: experimental demonstration with diffractive optics,” J. Mod. Opt. 46, 2355–2369 (1998). [CrossRef]
14.	F. Gori, “Polarization basis for vortex beams,” J. Opt. Soc. Am. A 18, 1612–1617 (2001). [CrossRef]
15.	L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, 1995), Sect. 3.2.
16.	G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists (Academic Press, New York, 2001), p. 681.
17.	R. H. Jordan and D. G. Hall, “Highly directional surface emission from concentric-circle gratings on planar optical waveguides: the field-expansion method,” J. Opt. Soc. Am. A 12, 84–94 (1995). [CrossRef]
18.	J. Tervo, P. Vahimaa, and J. Turunen, “On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,” J. Mod. Opt. 49, 1537–1543 (2002). [CrossRef]
19.	A. Lapucci and M. Ciofini, “Polarization state modifications in the propagation of high azimuthal order annular beams,” Opt. Express 9, 603–609 (2001),http: //www. opticsexpress.org/abstract.cfm?URI=OPEX-9-12-603. [CrossRef] [PubMed]
20.	J. Tervo and J. Turunen, “Self-imaging of electromagnetic fields,” Opt. Express 9, 622–630 (2001), http: //www. opticsexpress.org/abstract.cfm?URI=OPEX-9-12-622. [CrossRef] [PubMed]
21.	J. Tervo and J. Turunen, “Generation of vectorial propagation-invariant fields by polarization-grating axicons,” Opt. Commun. 192, 13–18 (2001). [CrossRef]

OCIS Codes
(260.5430) Physical optics : Polarization
(350.5500) Other areas of optics : Propagation

ToC Category:
Research Papers

History
Original Manuscript: July 31, 2002
Revised Manuscript: August 30, 2002
Published: September 9, 2002

Citation
Pertti Paakkonen, Jani Tervo, Pasi Vahimaa, Jari Turunen, and Franco Gori, "General vectorial decomposition of electromagnetic fields with application to propagation-invariant and rotating fields," Opt. Express 10, 949-959 (2002)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-10-18-949

Sort: Journal | Reset

References

W. D. Montgomery, ???Self-imaging objects of in.nite aperture,??? J. Opt. Soc. Am. 57, 772???778 (1967). [CrossRef]
K. Patorski, ???The self-imaging phenomenon and its applications,??? in Progress in Optics Vol. XXVII, E. Wolf, ed. (Elsevier, Amsterdam, 1989), Chap. 1.
J. Durnin, ???Exact solutions for nondifiracting beams. I. The scalar theory,??? J. Opt. Soc. Am. A 4, 651???654 (1987). [CrossRef]
J. Durnin, J. J. Miceli, Jr., and J. H. Eberly, ???Diffraction-free beams,??? Phys. Rev. Lett. 58, 1499???1501 (1987). [CrossRef] [PubMed]
Y. Y. Schechner, R. Piestun, and J. Shamir, ???Wave propagating with rotating intensity distributions,??? Phys. Rev. E 54, R50???R53 (1996). [CrossRef]
S. Chavez-Cerda, G. S. McDonald, and G. H. S. New, ???Nondiffracting Beams: travelling, standing, rotating and spiral waves,??? Opt. Commun. 123, 225???233 (1996). [CrossRef]
C. Paterson and R. Smith, ???Higher-order Bessel waves produced by axicon-type computergenerated holograms,??? Opt. Commun. 124, 121???130 (1996). [CrossRef]
S. R. Mishra, ???A vector wave analysis of a Bessel beam,??? Opt. Commun. 85, 159???161 (1991). [CrossRef]
J. Turunen and A. T. Friberg, ???Self-imaging and propagation-invariance in electromagnetic fields,??? Pure Appl. Opt. 2, 51???60 (1993). [CrossRef]
Z. Bouchal and M. Olivýk, ???Non-diffractive vector Bessel beams,??? J. Mod. Opt. 42, 1555???1566 (1995). [CrossRef]
Z. Bouchal, R. Horak and J. Wagner, ???Propagation-invariant electromagnetic fields,??? J. Mod. Opt. 43, 1905???1920 (1996). [CrossRef]
J. Tervo and J. Turunen, ???Rotating scale-invariant electromagnetic fields,??? Opt. Express 9, 9???15 (2001), <a href="http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9">http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-1-9</a>. [CrossRef] [PubMed]
P. Paakkonen, J. Lautanen, M. Honkanen, M. Kuittinen, J. Turunen, S. N. Khonina, V. V. Kotlyar, V. A. Soifer and A. T. Friberg, ???Rotating optical fields: experimental demonstration with diffractive optics,??? J. Mod. Opt. 46, 2355???2369 (1998). [CrossRef]
F. Gori, ???Polarization basis for vortex beams,??? J. Opt. Soc. Am. A 18, 1612???1617 (2001). [CrossRef]
L. Mandel, and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, Cambridge, 1995), Sect. 3.2.
G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists (Academic Press, New York, 2001), p. 681.
R. H. Jordan and D. G. Hall, ???Highly directional surface emission from concentric-circle gratings on planar optical waveguides: the field-expansion method,??? J. Opt. Soc. Am. A 12, 84???94 (1995). [CrossRef]
J. Tervo, P. Vahimaa, and J. Turunen, ???On propagation-invariant and self-imaging intensity distributions of electromagnetic fields,??? J. Mod. Opt. 49, 1537???1543 (2002). [CrossRef]
A. Lapucci and M. Ciofini, ???Polarization state modifications in the propagation of high azimuthal order annular beams,??? Opt. Express 9, 603???609 (2001), <a href="http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-12-603">http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-12-603</a>. [CrossRef] [PubMed]
J. Tervo and J. Turunen, ???Self-imaging of electromagnetic fields,??? Opt. Express 9, 622???630 (2001), <a href="http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-12-622">http://www.opticsexpress.org/abstract.cfm?URI=OPEX-9-12-622</a>. [CrossRef] [PubMed]
J. Tervo and J. Turunen, ???Generation of vectorial propagation-invariant fields by polarizationgrating axicons,??? Opt. Commun. 192, 13???18 (2001). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper