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

Controlling optical field collapse by elliptical symmetry hybrid polarization structure

Open Access Open Access

Abstract

Engineering the diversity of the spatial structure of polarization states offers a new approach to produce controllable and designable filament patterns. Here we predict theoretically and demonstrate experimentally the novel collapsing behaviors of the optical fields by using the elliptical symmetry hybrid polarization structure in a self-focusing Kerr medium. The results reveal that the distribution of the hybrid polarization states, the spin angular momentum (SAM) gradient, and the collapsing patterns could be manipulated by changing the eccentricity, the topological charge, and the initial phase of the optical field. Owing to the synergy of the hybrid polarization states and its spatial symmetry, the collapsing behaviors are controllable, and have the robust feature insensitive to the random noise. Our idea may offer an alternative route to produce the controllable and robust multiple filamentation in other nonlinear systems, thereby facilitating the development of additional surprising applications.

© 2018 Optical Society of America

1. INTRODUCTION

Field collapse, which occurs in many physical systems including optics, fluidics, plasma, and Bose–Einstein condensates, has attracted much attention, owing to its universality, complexity, and applicability [16]. A femtosecond (fs) optical field could undergo the collapse in a self-focusing Kerr medium under the input power beyond a certain critical power. As the collapse aggravates, there will occur the higher-order nonlinearity, which combines with the intrinsic diffraction of the optical field to counterbalance the self-focusing effect, resulting ultimately in the filamentation [713]. In general, the filamention patterns are initialized from the unpredictable symmetry breaking caused by the random noise, making it to be hardly controlled.

However, it is always anticipated for the realization of controllable field collapse and sequent filamentation, due to its practical applications, such as guiding microwave radiation [14], enhancement of terahertz emission [15], generation of millijoule-level supercontinuum in solid media [16], etc. For those purposes, some methods have been proposed by controlling the input power and divergence angle [17], shaping the field profile [18,19], using the amplitude/phase mask [2023], and introducing the spatial regularization [24]. In fact, the above methods are to manipulate the spatial structure of scalar optical fields in phase and amplitude, whereas the spatial structure of polarization states was rarely involved.

Polarization, as an intrinsic and fundamental vectorial nature of light, plays an indispensable role in the light–matter interaction. We first utilized the hybrid polarized vector optical fields (VOFs) to realize the controllable and designable filamentation [25,26], in which the symmetry of light–matter nonlinear interaction is actively broken by introducing the polarization-state-dependent refractive index change. However, the diversity of the spatial geometry of polarization states is not yet applied fully. Inasmuch as the hybrid polarized VOFs in [25,26] are only a special case, in which the polarization states exhibit a cylindrical symmetry or the ellipticity (or spin angular momentum, SAM) of the polarization states is a sin/cos function of the azimuthal angle. A question is raised that whether the azimuthally inhomogeneous distribution of hybrid polarization states in the field cross section, as a new approach, can be used to effectively control the field collapse and the sequent multiple filamentation. The generalized hybrid polarized vector fields, which break the cylindrical symmetry of the polarization states of the VOFs, become an interesting approach to control the filamentation.

Here we present an idea that the actively controllable collapsing patterns are achieved by designing the geometric configurations of the initial hybrid polarization states. One of the simplest methods is to design and generate a new kind of elliptic-symmetry azimuthally variant hybrid polarized VOFs (ES-AV-HP-VOFs) to explore how the fields collapse and the sequent multiple filamentation. The eccentricity, describing the elliptic geometry of the spatial hybrid polarization states, provides an additional degree of freedom assisting in engineering the azimuthally variant gradient of the hybrid polarization states and controlling the collapsing patterns. Besides the intensity, the contribution of the polarization states to the nonlinear refractive index change strongly affects the nonlinear dynamic property of the structured optical fields.

2. THEORY

An ES-AV-HP-VOF can be written as

E=A(r)(cosδe^H+jsinδe^V),
where e^H and e^V are the orthogonally horizontal and vertical unit vectors in the Cartesian coordinate system; A(r) represents the radially variant amplitude in the cross section of the ES-AV-HP-VOF and has a top-hat-like profile with a radius of r0 excluding a polarization singularity at the center; and δ=mϕ+δ0 (m is the topological charge and δ0 is the initial phase), respectively. In particular, ϕ is not the azimuthal angle in the polar coordinate system and ϕ=arctan(a2y/b2x)=arctan(a2tanϕ/b2)=arctan(a21x2/a2/bx), which is derived from the elliptic equation x2/a2+y2/b2=1(a>b>0), where a and b are the semi-major axis and semi-minor axis of the ellipse, respectively. ϕ (tanϕ=y/x) is the azimuthal angle in the polar coordinate system attached to the cross section of the ES-AV-HP-VOF described by Eq. (1). One should point out that when ϵ=0 or a=b, ϕ degenerates into ϕ, i.e., ϕ=ϕ. The eccentricity ϵ of an ES-AV-HP-VOF is described by ϵ=1b2/a2 with a range of ϵ[0,1]. Thus the eccentricity ϵ could provide an additional degree of freedom to control the azimuthal distribution of hybrid polarization states and the field collapsing behavior.

As is well known, the polarization of the optical field associates with the optical SAM. For the ES-AV-HP-VOF described by Eq. (1), its local SAM should be σ with σ=S3, where S3 is a normalized Stokes parameter in the representation of Poincaré sphere for describing the polarization states. σ=S3=+1 and σ=S3=1 correspond to the right-handed (RH) and left-handed (LH) circular polarizations, while σ=S3=0 indicates the linear polarization, respectively. Based on the definition of the Stokes parameter, S3 of the ES-AV-HP-VOF described by Eq. (1) should be

S3(ϕ)=σ(ϕ)=sin(2δ)=sin(2mϕ+2δ0)=sin[2marctan(a2tanϕ/b2)+2δ0].
S3 or the polarization state (SAM, σ) depends on the azimuthal angle ϕ while independent of the radial coordinate. For the case of ϵ=0 or a=b, the polarization state (SAM) is a sine function of ϕ as
S3(ϕ)|ϵ=0=σ(ϕ)|ϵ=0=sin(2mϕ+2δ0).
Figure 1 shows the distribution of the polarization state (SAM) of the ES-AV-HP-VOF (m=1, δ0=0, and δ0=π/4) with different eccentricities of ϵ=0,0.3,0.6,0.8. When ϵ0, the distribution geometry of the hybrid polarization states (as shown by the red arrows) exhibits a set of concentric ellipses (as shown by the black dashed ellipses) with the same eccentricity. Meanwhile, when ϵ=0, the ES-AV-HP-VOF degenerates into the cylindric-symmetry azimuthally variant hybrid polarized VOF (CS-AV-HP-VOF) [27]. The black background corresponds to the zero SAM or the linear polarization, the green (yellow) background indicates the + () SAM or the RH (LH) circular polarization, and the blackgreen (blackyellow) gradual color background represents the 0<σ<+ (0>σ>) SAM or the RH (LH) elliptical polarization. One can be seen that the variation of the polarization state (SAM, σ) in the azimuthal direction becomes more inhomogeneous as ϵ increases. In the vicinity of ϕ=0 and ϕ=π (ϕ=π/2 and ϕ=3π/2) directions, in detail, the variation of the polarization state becomes faster (slower) as ϵ increases.

 figure: Fig. 1.

Fig. 1. Schematic distributions of the polarization states and the SAM for the ES-AV-HP-VOFs with m=1 and δ0=0 for different eccentricities (ϵ=0,0.3,0.6,0.8).

Download Full Size | PDF

To theoretically explore the nonlinear propagation behavior of an ES-AV-HP-VOF in the isotropic Kerr medium, a (2+1)-dimensional vector-version nonlinear-Schrödinger (NLS) equation should be used. Under the slowly varying amplitude approximation, the NLS equation can be divided into a pair of coupled NLS equations for two orthogonal components as follows:

ψHζ=j42ψH+jαP3PC(3|ψH|2ψH+2|ψV|2ψH+ψV2ψH*),
ψVζ=j42ψV+jαP3PC(3|ψV|2ψV+2|ψH|2ψV+ψH2ψV*),
where ψq (q=H,V) is the dimensionless q-component normalized by the total field as
ψq(ρ,ϕ;ζ)=Eq(ρ,ϕ;ζ)[|EH(ρ,ϕ;ζ)|2+|EV(ρ,ϕ;ζ)|2]ρdρdϕ.
Here P=2n0ε0c[|EH(ρ,ϕ;ζ)|2+EV(ρ,ϕ;ζ)|2]ρdρdϕ is the power and PC=αλ2/4πn0n2 is the critical power for self-focusing, respectively. ρ=r/r0 and ζ=z/Ld (Ld=πr02/λ) are the dimensionless cylindrical coordinates, n0 and n2 are the linear and nonlinear refractive indices of the Kerr medium, α is a constant dependent on the initial field shape [28], c is the speed of light in vacuum, and λ is the wavelength of light, respectively. On the right-hand side of Eqs. (4) and (5), the first term represents the contribution from the diffraction described by the transverse Laplacian 2=2/ρ2+ρ1/ρ+ρ22/ϕ2 and the second term is from the Kerr nonlinearity.

3. NUMERICAL SIMULATION

The ES-AV-HP-VOF might have some properties more useful than the CS-AV-HP-VOF, because it has an additional controllable degree of freedom (eccentricity ϵ) and the distribution diversity of hybrid polarization states. Here we focus on their collapsing behaviors in the self-focusing Kerr medium. Under P=10PC, a large number of simulations are performed by the beam propagation method [29], which is used to numerically solve Eq. (4). When the spatial random noise is added for simulations, we keep the random noise with a level of 10% in amplitude, which follows a uniform distribution within a range of 10%.

A. ES-AV-HP-VOFs (m=1,δ0=0)

First, the collapsing behaviors of the ES-AV-HP-VOFs (m=1, δ0=0) with different eccentricity ϵ in the Kerr medium are studied, as shown in Fig. 2. The initial input field is pre-focused on the input face of the Kerr medium (i.e., the focal plane is in the Kerr medium). Although all the initial input ES-AV-HP-VOFs have the top-hat-like profile in intensity excluding the central polarization singularity, the pre-focused fields incident on the input plane of the Kerr medium are different, as the simulated patterns shown in the first column of Fig. 2, due to the different ϵ.

 figure: Fig. 2.

Fig. 2. Simulated nonlinear propagation behaviors of the ES-AV-HP-VOFs (m=1,δ0=0) with different ϵ(=0,0.6,0.8) in the Kerr medium. The first and second rows show the cases without and with random noise (ϵ=0), respectively. For comparison, the ES-AV-HP-VOFs (ϵ0) without and with random noise are also shown in the third and fourth rows for the case of ϵ=0.6 and in the fifth and sixth rows for the case of ϵ=0.8. The four columns (from left to right) correspond to four propagation distances (ζ=0,0.6,1.2,1.8).

Download Full Size | PDF

Here we classify into the two cases based on the value of ϵ: ϵ=0 and ϵ0. For the case of ϵ=0, the collapsing evolutions simulated by Eqs. (4) and (5) are shown in the first and second rows of Fig. 2. In this case, the ES-AV-HP-VOF degenerates into a CS-AV-HP-VOF. The four deterministic collapsing patterns produced by the CS-AV-HP-VOF (m=1, δ0=0) are located at the azimuthal locations of ϕ=0, π/2, π, 3π/2, which is in agreement with [25]. For the case of ϵ0, including ϵ=0.6 (third and fourth rows) and ϵ=0.8 (fifth and sixth rows) in Fig. 2, (i) the larger ϵ leads to the axial-symmetry breaking behaviors by changing the eccentricity of the input field; (ii) the ES-AV-HP-VOFs (m=1, δ0=0) collapse to converge into the two (instead of four) deterministic filaments located at ϕ=π/2, 3π/2, corresponding to the linear polarization in the semi-minor axis of the ellipse; and (iii) the collapsing patterns are insensitive to the random noise.

To clearly show the evolution of azimuthal-variant polarization states during the collapsing of ES-AV-HP-VOFs (m=1, δ0=0) with ϵ=0 and ϵ0 in the isotropic Kerr medium, we simulate the Stocks parameter S3, as shown in Fig. 3. For the case of ϵ=0, S3 of the initial input field obeys S3(ϕ)=sin(2ϕ), shown by the black curve in Fig. 3(a); as the propagation distance increases, |S3| of the local polarization state has a tendency to decrease by the coupling between the two orthogonally polarized components in the Kerr medium, as shown by the other three curves in Fig. 3(a). For the case of ϵ=0.8, S3 of the initial input field obeys S3(ϕ)=sin[2marctan(tanϕ/0.36)], shown by the black curve in Fig. 3(b); however, the nonlinear coupling does not always lead to the decrease of |S3| of the local polarization state, which is different from Fig. 3(a). Interestingly, for both cases of ϵ=0 and ϵ0, so long as δ0=0 the local polarization states have no change to keep the initial linear polarization at four special azimuthal locations of ϕ(n)=nπ/2 (n=0,1,2,3). For the case of ϵ=0, the nonlinear coupling also does not result in the change of the azimuthal locations for the initial circular polarization (|S3|=1); while for the case of ϵ0, the nonlinear coupling results in the shift of the azimuthal locations for the initial circular polarization (|S3|=1). After comparing Fig. 3 with Fig. 2, we find that the final filaments could occur at those azimuthal locations, in which the initial polarization states are linearly polarized; in particular, they are more favorite to occur at the linearly polarized locations, which have the smaller azimuthal gradient (|S3/ϕ|) of S3. For instance, for the case of ϵ=0.80, the two filaments occur at ϕ=π/2 and 3π/2 instead of ϕ=0 and π, because |S3/ϕ||ϕ=π/2,3π/2<|S3/ϕ||ϕ=0,π. In contrast, for the case of ϵ=0, the four filaments occur at ϕ(n)=nπ/2 (n=0,1,2,3), because |S3/ϕ||ϕ=π/2,3π/2=|S3/ϕ||ϕ=0,π.

 figure: Fig. 3.

Fig. 3. Simulated evolution S3 of the azimuthal-variant polarization states for the ES-AV-HP-VOFs (m=1,δ0=0) with different eccentricity. (a) ϵ=0 and (b) ϵ=0.8.

Download Full Size | PDF

B. ES-AV-HP-VOFs (m=1,δ0=π/4)

We now change the initial phase δ0 (within a range of δ0[0,π/4]) of the ES-AV-HP-VOFs (m=1) with different eccentricities of ϵ=0,0.6,0.8. The collapsing behaviors simulated by Eqs. (4) and (5) are shown in Fig. 4. We also classify into the two cases: ϵ=0 and ϵ0. For the case of ϵ=0, the ES-AV-HP-VOF degenerates into the CS-AV-HP-VOF (m=1, δ0=π/4). Clearly, the four deterministic filaments produced by the CS-AV-HP-VOF (m=1, δ0=π/4) are counter-clockwise rotated by π/4 with respect to the CS-AV-HP-VOF (m=1, δ0=0) shown in Fig. 2, which correspond to the linear polarization. However, when ϵ0, new features are found: (i) the filaments do not reduce to two while still keeping to four filaments for any value of ϵ, which is different from the case (m=1, δ0=0, ϵ0) of Fig. 2; and (ii) as ϵ increases, the four filaments are changed from a square for ϵ=0 into an elongated rectangle for ϵ0.

 figure: Fig. 4.

Fig. 4. Simulated nonlinear propagation behaviors of ES-AV-HP-VOFs (m=1,δ0=π/4) with the different ϵ(=0,0.6,0.8) in the Kerr medium. The first and second rows show the cases without and with random noise (ϵ=0), respectively. For comparison, the ES-AV-HP-VOFs (ϵ0) without and with random noise are also shown in the third and fourth rows for the case ϵ=0.6 and in the fifth and sixth rows for the case ϵ=0.8. The four columns (from left to right) correspond to four propagation distances (ζ=0,0.6,1.2,1.8).

Download Full Size | PDF

Meanwhile, we also simulate the evolution of S3 (or SAM) during the collapsing of ES-AV-HP-VOFs (m=1,δ0=π/4) with ϵ=0 and ϵ0 in the Kerr medium, as shown in Figs. 5(a) and 5(b), respectively. For the case of ϵ=0, S3 of the initial input fields obeys S3(ϕ)=sin(2ϕ+π/2), shown by the black curve in Fig. 5(a). As shown by the other three curves in Fig. 5(a), the nonlinear coupling results in the decrease of |S3| of the local polarization state at any azimuthal location. In particular, in eight special azimuthal locations with the linear or circular polarizations, the nonlinear coupling does not cause the change of their polarization states. These features are the same as the case of δ0=0 in Fig. 3(a).

 figure: Fig. 5.

Fig. 5. Simulated evolution S3 of the azimuthal-variant polarization states for the ES-AV-HP-VOFs (m=1,δ0=π/4) with different eccentricity. (a) ϵ=0 and (b) ϵ=0.8.

Download Full Size | PDF

For the case of ϵ=0.80, S3 of the initial input field obeys S3(ϕ)=sin[2marctan(tanϕ/0.36)+π/2], shown by the black curve in Fig. 5(b). As shown by the other three curves, the nonlinear coupling does not always lead to the decrease of |S3| of the local polarization state, which is the same as Fig. 3(b). However, in four special azimuthal locations with the circular polarizations, the nonlinear coupling does not cause the change of their polarization states, which is different from the case of δ0=0 and ϵ0 in Fig. 3(b).

After comparing Fig. 4 with Fig. 5, we find that (i) for ϵ=0, the four filaments locate at ϕ(n)=(2n+1)π/4 (n=0,1,2,3), in which the polarization states are always linearly polarized; and (ii) for ϵ0, the filaments do not locate at ϕ(n)=(2n+1)π/4 (n=0,1,2,3), in which the initial polarization states are not linearly polarized and are changed due to the nonlinear coupling.

C. ES-AV-HP-VOFs (m=1,δ0=3π/8)

Given the ES-AV-HP-VOFs (m=1,δ0=3π/8) with different eccentricities (ϵ=0,0.6,0.8), we continue to study their collapsing behaviors based on Eqs. (4) and (5), as shown in Fig. 6. For the case of ϵ=0, the ES-AV-HP-VOF degenerates into a CS-AV-HP-VOF (m=1,δ0=3π/8), and the simulation results are shown in the first and second rows of Fig. 6. The four deterministic collapsing patterns produced by the CS-AV-HP-VOF (m=1,δ0=3π/8) are counter-clockwise rotated π/8 with respect to the case (m=1, δ0=0, ϵ=0) in Fig. 2. Four filaments are located at the azimuthal locations of ϕ=π/8, 5π/8, 9π/8, 13π/8, in which the initial polarization states are linearly polarized. However, when ϵ0, including ϵ=0.6 (third and fourth rows) and ϵ=0.8 (fifth and sixth rows), ϵ0 leads to the ES-AV-HP-VOF to collapse two (instead of four) filaments, which is the same as the case (m=1, δ0=0, ϵ0). The two filaments are located at the azimuthal locations which are very close to but slightly smaller than ϕ=5π/8,13π/8.

 figure: Fig. 6.

Fig. 6. Simulated nonlinear propagation behaviors of ES-AV-HP-VOFs (m=1 and δ0=3π/8) with the different ϵ(=0,0.6,0.8) in the Kerr medium. The first and second rows show the cases without and with random noise (ϵ=0), respectively. For comparison, the ES-AV-HP-VOFs (ϵ0) without and with random noise are also shown in the third and fourth rows for the case ϵ=0.6 and in the fifth and sixth rows for the case ϵ=0.8. The four columns (from left to right) correspond to four propagation distances (ζ=0,0.6,1.2,1.8).

Download Full Size | PDF

The evolutions of S3 (or SAM) during the collapsing of ES-AV-HP-VOFs (m=1,δ0=3π/8) with ϵ=0 and ϵ=0.80 in the isotropic Kerr medium, are also simulated, as shown in Figs. 7(a) and 7(b), respectively. For the case of ϵ=0, S3 of the initial input fields obeys S3(ϕ)=sin(2ϕ+3π/4), shown by the black curve in Fig. 7(a). As shown by the other three curves in Fig. 7(a), the nonlinear coupling results in the decrease of |S3| of the local polarization state at any azimuthal location. In particular, in eight special azimuthal locations at ϕ(n)=nπ/2+π/8 for linear polarizations and ϕ(n)=nπ/2+3π/8 (where n=0,1,2,3), the nonlinear coupling does not cause the change of their polarization states. These features are the same as the case of δ0=0 in Figs. 3(a) and 5(a).

 figure: Fig. 7.

Fig. 7. Simulated evolution S3 of the azimuthal-variant polarization states for the ES-AV-HP-VOFs (m=1 and δ0=3π/8) with different eccentricity. (a) ϵ=0 and (b) ϵ0.

Download Full Size | PDF

For the case of ϵ=0.80, S3 of the initial input field obeys S3(ϕ)=sin[2marctan(tanϕ/0.36)+3π/4], shown by the black curve in Fig. 7(b). As shown by the other three curves, the nonlinear coupling does not always lead to the decrease of |S3| of the local polarization state, which is the same as Figs. 3(b) and 5(b). However, this case has a completely different feature that the polarization states at all the azimuthal locations are changed by the nonlinear coupling.

After comparing Fig. 6 with Fig. 7, we find that (i) for ϵ=0, the four filaments locate at ϕ(n)=nπ/2+π/8 (where n=0,1,2,3), in which the polarization states are always linearly polarized; and (ii) for ϵ0, the filaments become into two located at the azimuthal locations which are very close to but slightly smaller than ϕ=5π/8,13π/8, as mentioned above. In particular, at the two azimuthal locations, the finial polarization states after the nonlinear coupling favor to tend to become into the linear polarization.

Finally, we should point out that the dashed curves in Figs. 3, 5, and 7 show the azimuthal intensity distribution after the collapsing of the ES-AV-HP-VOFs, in which the intensity peaks correspond to the locations of the filaments. The purpose is to clearly show the correspondence of filaments with the azimuthal locations and the local polarization states.

4. DISCUSSION

After summarizing all the simulation results, we obtain the following conclusions: (i) ϵ indeed provides a new degree of freedom for controlling the collapse of the ES-AV-HP-VOFs; (ii) when ϵ=0, the ES-AV-HP-VOFs with m always produce 4m filaments located at the azimuthal locations, in which the initial and final polarization states are linearly polarized; (iii) when ϵ0, the ES-AV-HP-VOFs with m always produce 4m or 2m filaments, depending on the initial phase δ0; and (iv) the filaments are more favorite to occur at azimuthal locations, in which the finial polarization states after the nonlinear coupling tend to become into the linear polarization and have the smaller azimuthal gradient (|S3/ϕ|) of S3.

To understand the physics behind the collapsing behaviors, the azimuthal self-focusing model is very beneficial. The field collapse originates from the spatially variant refractive index change Δn and requires the existence of self-focusing nucleation(s) where the induced Δn must be local maximum. As is well known, the induced refractive index change Δn is directly proportional to both the local light intensity I(ϕ) and nonlinear refraction index coefficient n2(ϕ). For the CS-AV-HP-VOFs (ϵ=0), because its polarization states constitute a cylindrical symmetrical structure, the pre-focused intensity holds the axial symmetry, as shown in the first and second rows of Figs. 2, 4, and 6. The induced refractive index change Δn is only caused by the azimuthal-variant n2, whose value depends on the local polarization state. Accordingly, the energy tends to accumulate at the local position of linear polarization, because of Δnlin>Δnell>Δncir (Δnlin=1.5Δncir). For the ES-AV-HP-VOF (ϵ0), the elliptical-symmetry structure results in the intensity redistribution from the homogeneous distribution into the two separate strong peaks. Due to the synergy of the polarization structure and intensity distribution, the induced Δn strongly affects the nonlinear dynamic behaviors of the VOFs. Interestingly, the competition between the inhomogeneities of both intensity I(ϕ) and nonlinear refraction index coefficient n2(ϕ) (originates from the azimuthal-variant polarization states) may cause the change of self-focusing nucleations, so the ES-AV-HP-VOFs (ϵ0) with different δ0 will have a variety of collapsing behaviors.

Here we should point out that the quantity and locations of filamentation are dominated by the pre-collapse process, hence in all the above simulations we take into account the Kerr nonlinearity only, which is enough to elucidate the collapsing behaviors of the ES-AV-HP-VOFs. To further test, we also simulate the multiple filamentation when taking into account the higher-order nonlinearity important to the post-collapse process. Figure 8 shows the simulation results of the CS-AV-HP-VOF (m=1,δ0=0) and the ES-AV-HP-VOFs (m=1,ϵ=0.8) with δ0=0, π/4, 3π/8. Beyond the collapse point, the sizes and the intensity of the filaments experience synchronous oscillations with the propagation distance. The filamentation patterns in the output face of the nonlinear medium are consistent with the pre-collapse points.

 figure: Fig. 8.

Fig. 8. Simulated collapsing evolution and multiple filamentation of the ES-AV-HP-VOFs. (a) CS-AV-HP-VOF with m=1 and δ0=0; (b) ES-AV-HP-VOF with m=1, ϵ=0.8, and δ0=0; (c) ES-AV-HP-VOF with m=1, ϵ=0.8, and δ0=π/4; and (d) ES-AV-HP-VOF with m=1, ϵ=0.8, and δ0=3π/8.

Download Full Size | PDF

In addition, we also carried out some simulations under the spatial random noise with a level of 30% in amplitude, we keep the random noise, which follows a uniform distribution within a range of 30%. The simulation results (not shown here) confirm that our idea may produce the robust multiple filamentation against the random noise.

5. EXPERIMENTAL VERIFICATION

We follow the same method as [27] for creating the fs ES-AV-HP-VOFs required for the experimental verification. The used light source is a Ti:sapphire regenerative amplifier fs laser system (Coherent Inc.), which provides a fundamental Gaussian mode with a central wavelength of 800 nm, a pulse duration of 35 fs, and a repetition rate of 1 kHz. An achromatic 1/2-wave plate and a broadband polarized beam splitter were used to control the laser fluence input on the sample. Another achromatic 1/2-wave plate was used to change the polarization direction of the input fs laser into the vector field generation unit, where all the elements are achromatic to suppress the pulse broadening as much as possible. The created ES-AV-HP-VOFs have a “top-hat-like” spatial profile, excluding a central singularity caused by the polarization uncertainty, with a radius of 1.5 mm, a pulse energy of 6.0 mJ, and a pulse duration of 65fs.

To experimentally confirm our theoretical prediction and numerical simulations, all the required vector fields with different m, δ0, and ϵ were created by using the experimental scheme [27]. The created ES-AV-HP-VOFs were weakly pre-focused onto the CS2 cell with a length of 30 mm, by an achromatic lens with a focal length of f=300mm. The intensity patterns of the field transmitted from the CS2 cell are imaged on a detector (Beamview, Coherent Inc.) using an achromatic lens with a focal length of f=60mm.

As shown in the first column of Fig. 9, corresponding to the CS-AV-HP-VOFs with (m=1,ϵ=0), the four filaments occur at the four azimuthal locations [ϕ(n)=(n+1)π/2δ0, where n=0,1,2,3] with the local linear polarization, as studied in [25]. As shown in the first and second rows with (m=1,δ0=0), the four filaments reduce to the two filaments located at ϕ=π/2 and ϕ=3π/2, as ϵ increases. As shown in the third and four rows with (m=1,δ0=π/4), however, the four filaments are still kept but are changed from a square into an elongated rectangle, as ϵ increases from ϵ=0 to ϵ=0.8. As shown in the fifth and sixth rows with (m=1,δ0=3π/8), as ϵ increases, the four filaments also reduce to the two filaments located at near ϕ=5π/8 and ϕ=13π/8.

 figure: Fig. 9.

Fig. 9. Measured and simulated collapsing behaviors of the ES-AV-HP-VOFs (m=1) with the different ϵ(ϵ=0,0.6,0.8) in the Kerr medium. The first and second rows show the measured and simulated collapsing patterns of the ES-AV-HP-VOFs (m=1,δ0=0), respectively. For comparison, the measured and simulated collapsing patterns of the ES-AV-HP-VOFs are also shown in the third and fourth rows for the case (m=1,δ0=π/4) and in the fifth and sixth rows for the case (m=1,δ0=3π/8), respectively. The four columns (from the left to the right) correspond to four different ϵ(ϵ=0,0.6,0.8).

Download Full Size | PDF

We also explore experimentally the collapsing behaviors of the ES-AV-HP-VOFs with m=2. As examples, Fig. 10 shows the collapsing pattens of the ES-AV-HP-VOFs (m=2,δ0=0) for different values of ϵ. The first row shows the azimuthal distributions of the SAM (S3) for ϵ=0,0.3,0.6,0.8. Clearly, the SAM (S3) exhibits the change from the homogeneous distribution to the inhomogeneous one, as ϵ increases. As shown in the second (simulation) and third (experiment) rows, the ES-AV-HP-VOFs (m=2,δ0=0) collapse into eight filaments for ϵ=0, but four filaments occupying the four vertices of a rectangle (instead of a square) for ϵ=0.8.

 figure: Fig. 10.

Fig. 10. SAM and the measured and simulated collapsing patterns of the ES-AV-HP-VOFs (m=2,δ0=0) with different eccentricity ϵ(=0,0.3,0.6,0.8).

Download Full Size | PDF

Clearly, the experimental results are in good agreement with the theoretical and numerical results for the ES-AV-HP-VOFs.

6. CONCLUSION

We have presented the theoretical prediction and performed the experimental evidence on the novel collapsing behaviors of the ES-AV-HP-VOFs in the self-focusing Kerr medium. Engineering the diversity of spatial distribution of polarization states is indeed an new approach to control the optical field collapse. The eccentricity, which describes the geometric configurations of the ES-AV-HP-VOFs, provides an additional degree of freedom assisting in controlling the variation of the SAM gradient of the hybrid polarization states and then the optical field collapse. Both theory and experiment demonstrate that the eccentricity ϵ, the topological charge m, and the initial phase δ0 of the ES-AV-HP-VOFs can indeed control the locations and number of the filaments produced in the Kerr medium. When ϵ=0, the ES-AV-HP-VOF with m degenerates into the CS-AV-HP-VOF with m, and the interaction with the Kerr medium will collapse into 4m filaments exhibiting a C4m rotation symmetry. Due to the introduction of ϵ0, the axial symmetry of the optical field is broken, the ES-AV-HP-VOF with m and ϵ0 can collapse into 4m or 2m filaments depending on the initial phase δ0, but 4m filaments do not exhibit the C4m rotation symmetry. Our idea may offer an alternative route to produce the controllable and robust multiple filamentation insensitive to the random noise, thereby facilitating the development of additional surprising applications.

Funding

National Key R&D program of China (13JCZDJC33800, 2017YFA0303700); National Natural Science Foundation of China (NSFC) (11534006, 11674184, 11774183); Natural Science Foundation of Tianjin City (16JCZDJC31300).

Acknowledgment

This work is supported by Collaborative Innovation Center of Extreme Optics.

REFERENCES

1. P. A. Robinson, “Nonlinear wave collapse and strong turbulence,” Rev. Mod. Phys. 69, 507–574 (1997). [CrossRef]  

2. C. D’Amico, A. Houard, M. Franco, B. Prade, A. Mysyrowicz, A. Couairon, and V. T. Tikhonchuk, “Conical forward THz emission from femtosecond-laser-beam filamentation in air,” Phys. Rev. Lett. 98, 235002 (2007). [CrossRef]  

3. R. R. Alfano and S. L. Shapiro, “Emission in the region 4000 to 7000 Å via four photon coupling in glass,” Phys. Rev. Lett. 24, 584–587 (1970). [CrossRef]  

4. P. Rohwetter, J. Kasparian, K. Stelmaszczyk, Z. Hao, S. Henin, N. Lascoux, W. M. Nakaema, Y. Petit, M. Queißer, R. Salamé, E. Salmon, L. Wöste, and J. P. Wolf, “Laser-induced water condensation in air,” Nat. Photonics 4, 451–456 (2010). [CrossRef]  

5. F. Belgiorno, S. L. Cacciatori, M. Clerici, V. Gorini, G. Ortenzi, L. Rizzi, E. Rubino, V. G. Sala, and D. Faccio, “Hawking radiation from ultrashort laser pulse filaments,” Phys. Rev. Lett. 105, 203901 (2010). [CrossRef]  

6. G. Fibich and B. Ilan, “Vectorial and random effects in self-focusing and in multiple filamentation,” Phys. D 157, 112–146 (2001). [CrossRef]  

7. A. Couairon and A. Mysyrowicz, “Femtosecond filamentation in transparent media,” Phys. Rep. 441, 47–189 (2007). [CrossRef]  

8. L. Bergé, S. Skupin, R. Nuter, J. Kasparian, and J. P. Wolf, “Ultrashort filaments of light in weakly ionized, optically transparent media,” Rep. Prog. Phys. 70, 1633–1713 (2007). [CrossRef]  

9. R. Y. Chiao, E. Garmire, and C. H. Townes, “Self-trapping of optical beams,” Phys. Rev. Lett. 13, 479–482 (1964). [CrossRef]  

10. P. L. Kelley, “Self-focusing of optical beams,” Phys. Rev. Lett. 15, 1005–1008 (1965). [CrossRef]  

11. S. Varma, Y. H. Chen, and H. M. Milchberg, “Trapping and destruction of long-range high-intensity optical filaments by molecular quantum wakes in air,” Phys. Rev. Lett. 101, 205001 (2008). [CrossRef]  

12. J. Kasparian and J. P. Wolf, “Physics and applications of atmospheric nonlinear optics and filamentation,” Opt. Express 16, 466–493 (2008). [CrossRef]  

13. J. Kasparian, M. Rodriguez, G. Méjean, J. Yu, E. Salmon, H. Wille, R. Bourayou, S. Frey, Y. B. André, A. Mysyrowicz, R. Sauerbrey, J. P. Wolf, and L. Wöste, “White-light filaments for atmospheric,” Science 301, 61–64 (2003). [CrossRef]  

14. Y. Ren, M. Alshershby, Z. Hao, Z. Zhao, and J. Lin, “Microwave guiding along double femtosecond filaments in air,” Phys. Rev. E 88, 013104 (2013). [CrossRef]  

15. J. Zhao, L. Guo, W. Chu, B. Zeng, H. Gao, Y. Cheng, and W. Liu, “Simple method to enhance terahertz radiation from femtosecond laser filament array with a step phase plate,” Opt. Lett. 40, 3838–3841 (2015). [CrossRef]  

16. A. Camino, Z. Hao, X. Liu, and J. Lin, “High spectral power femtosecond supercontinuum source by use of microlens array,” Opt. Lett. 39, 747–750 (2014). [CrossRef]  

17. Z. Jin, J. Zhang, M. H. Xu, X. Lu, Y. T. Li, Z. H. Wang, Z. Y. Wei, X. H. Yuan, and W. Yu, “Control of filamentation induced by femtosecond laser pulses propagating in air,” Opt. Express 13, 10424–10430 (2005). [CrossRef]  

18. A. Dubietis, G. Tamošauskas, G. Fibich, and B. Ilan, “Multiple filamentation induced by input-beam ellipticity,” Opt. Lett. 29, 1126–1128 (2004). [CrossRef]  

19. T. D. Grow and A. L. Gaeta, “Dependence of multiple filamentation on beam ellipticity,” Opt. Express 13, 4594–4599 (2005). [CrossRef]  

20. T. Pfeifer, L. Gallmann, M. J. Abel, D. M. Neumark, and S. R. Leone, “Circular phase mask for control and stabilization of single optical filaments,” Opt. Lett. 31, 2326–2328 (2006). [CrossRef]  

21. V. P. Kandidov, N. Aközbek, M. Scalora, O. G. Kosareva, A. V. Nyakk, Q. Luo, S. A. Hosseini, and S. L. Chin, “A method for spatial regularisation of a bunch of filaments in a femtosecond laser pulse,” Quantum Electron. 34, 879–880 (2004). [CrossRef]  

22. Z. Q. Hao, K. Stelmaszczyk, P. Rohwetter, W. M. Nakaema, and L. Woeste, “Femtosecond laser filament-fringes in fused silica,” Opt. Express 19, 7799–7806 (2011). [CrossRef]  

23. Z. Q. Hao, J. Zhang, T. T. Xi, X. H. Yuan, Z. Y. Zheng, X. Lu, M. Y. Yu, Y. T. Li, Z. H. Wang, W. Zhao, and Z. Y. Wei, “Optimization of multiple filamentation of femtosecond laser pulses in air using a pinhole,” Opt. Express 15, 16102–16109 (2007). [CrossRef]  

24. W. Walasik and N. M. Litchinitser, “Dynamics of large femtosecond filament arrays: possibilities, limitations, and trade-offs,” ACS Photonics 3, 640–646 (2016). [CrossRef]  

25. S. M. Li, Y. N. Li, X. L. Wang, L. J. Kong, K. Lou, C. H. Tu, Y. J. Tian, and H. T. Wang, “Taming the collapse of optical fields,” Sci. Rep. 2, 1007 (2012). [CrossRef]  

26. S. M. Li, Z. C. Ren, L. J. Kong, S. X. Qian, C. H. Tu, Y. N. Li, and H. T. Wang, “Unveiling stability of multiple filamentation caused by axial symmetry breaking of polarization,” Photonics Res. 4, B29–B34 (2016). [CrossRef]  

27. X. L. Wang, Y. N. Li, J. Chen, C. S. Guo, J. P. Ding, and H. T. Wang, “A new type of vector fields with hybrid states of polarization,” Opt. Express 18, 10786–10795 (2010). [CrossRef]  

28. G. Fibich and A. L. Gaeta, “Critical power for self-focusing in bulk media and in hollow waveguides,” Opt. Lett. 25, 335–337 (2000). [CrossRef]  

29. G. Agrawal, Nonlinear Fiber Optics (Academic, 2012).

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 (10)

Fig. 1.
Fig. 1. Schematic distributions of the polarization states and the SAM for the ES-AV-HP-VOFs with m = 1 and δ 0 = 0 for different eccentricities ( ϵ = 0 , 0.3 , 0.6 , 0.8 ).
Fig. 2.
Fig. 2. Simulated nonlinear propagation behaviors of the ES-AV-HP-VOFs ( m = 1 , δ 0 = 0 ) with different ϵ ( = 0 , 0.6 , 0.8 ) in the Kerr medium. The first and second rows show the cases without and with random noise ( ϵ = 0 ), respectively. For comparison, the ES-AV-HP-VOFs ( ϵ 0 ) without and with random noise are also shown in the third and fourth rows for the case of ϵ = 0.6 and in the fifth and sixth rows for the case of ϵ = 0.8 . The four columns (from left to right) correspond to four propagation distances ( ζ = 0 , 0.6 , 1.2 , 1.8 ).
Fig. 3.
Fig. 3. Simulated evolution S 3 of the azimuthal-variant polarization states for the ES-AV-HP-VOFs ( m = 1 , δ 0 = 0 ) with different eccentricity. (a)  ϵ = 0 and (b)  ϵ = 0.8 .
Fig. 4.
Fig. 4. Simulated nonlinear propagation behaviors of ES-AV-HP-VOFs ( m = 1 , δ 0 = π / 4 ) with the different ϵ ( = 0 , 0.6 , 0.8 ) in the Kerr medium. The first and second rows show the cases without and with random noise ( ϵ = 0 ), respectively. For comparison, the ES-AV-HP-VOFs ( ϵ 0 ) without and with random noise are also shown in the third and fourth rows for the case ϵ = 0.6 and in the fifth and sixth rows for the case ϵ = 0.8 . The four columns (from left to right) correspond to four propagation distances ( ζ = 0 , 0.6 , 1.2 , 1.8 ).
Fig. 5.
Fig. 5. Simulated evolution S 3 of the azimuthal-variant polarization states for the ES-AV-HP-VOFs ( m = 1 , δ 0 = π / 4 ) with different eccentricity. (a)  ϵ = 0 and (b)  ϵ = 0.8 .
Fig. 6.
Fig. 6. Simulated nonlinear propagation behaviors of ES-AV-HP-VOFs ( m = 1 and δ 0 = 3 π / 8 ) with the different ϵ ( = 0 , 0.6 , 0.8 ) in the Kerr medium. The first and second rows show the cases without and with random noise ( ϵ = 0 ), respectively. For comparison, the ES-AV-HP-VOFs ( ϵ 0 ) without and with random noise are also shown in the third and fourth rows for the case ϵ = 0.6 and in the fifth and sixth rows for the case ϵ = 0.8 . The four columns (from left to right) correspond to four propagation distances ( ζ = 0 , 0.6 , 1.2 , 1.8 ).
Fig. 7.
Fig. 7. Simulated evolution S 3 of the azimuthal-variant polarization states for the ES-AV-HP-VOFs ( m = 1 and δ 0 = 3 π / 8 ) with different eccentricity. (a)  ϵ = 0 and (b)  ϵ 0 .
Fig. 8.
Fig. 8. Simulated collapsing evolution and multiple filamentation of the ES-AV-HP-VOFs. (a) CS-AV-HP-VOF with m = 1 and δ 0 = 0 ; (b) ES-AV-HP-VOF with m = 1 , ϵ = 0.8 , and δ 0 = 0 ; (c) ES-AV-HP-VOF with m = 1 , ϵ = 0.8 , and δ 0 = π / 4 ; and (d) ES-AV-HP-VOF with m = 1 , ϵ = 0.8 , and δ 0 = 3 π / 8 .
Fig. 9.
Fig. 9. Measured and simulated collapsing behaviors of the ES-AV-HP-VOFs ( m = 1 ) with the different ϵ ( ϵ = 0 , 0.6 , 0.8 ) in the Kerr medium. The first and second rows show the measured and simulated collapsing patterns of the ES-AV-HP-VOFs ( m = 1 , δ 0 = 0 ), respectively. For comparison, the measured and simulated collapsing patterns of the ES-AV-HP-VOFs are also shown in the third and fourth rows for the case ( m = 1 , δ 0 = π / 4 ) and in the fifth and sixth rows for the case ( m = 1 , δ 0 = 3 π / 8 ), respectively. The four columns (from the left to the right) correspond to four different ϵ ( ϵ = 0 , 0.6 , 0.8 ) .
Fig. 10.
Fig. 10. SAM and the measured and simulated collapsing patterns of the ES-AV-HP-VOFs ( m = 2 , δ 0 = 0 ) with different eccentricity ϵ ( = 0 , 0.3 , 0.6 , 0.8 ) .

Equations (6)

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

E = A ( r ) ( cos δ e ^ H + j sin δ e ^ V ) ,
S 3 ( ϕ ) = σ ( ϕ ) = sin ( 2 δ ) = sin ( 2 m ϕ + 2 δ 0 ) = sin [ 2 m arctan ( a 2 tan ϕ / b 2 ) + 2 δ 0 ] .
S 3 ( ϕ ) | ϵ = 0 = σ ( ϕ ) | ϵ = 0 = sin ( 2 m ϕ + 2 δ 0 ) .
ψ H ζ = j 4 2 ψ H + j α P 3 P C ( 3 | ψ H | 2 ψ H + 2 | ψ V | 2 ψ H + ψ V 2 ψ H * ) ,
ψ V ζ = j 4 2 ψ V + j α P 3 P C ( 3 | ψ V | 2 ψ V + 2 | ψ H | 2 ψ V + ψ H 2 ψ V * ) ,
ψ q ( ρ , ϕ ; ζ ) = E q ( ρ , ϕ ; ζ ) [ | E H ( ρ , ϕ ; ζ ) | 2 + | E V ( ρ , ϕ ; ζ ) | 2 ] ρ d ρ d ϕ .
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.