1. Introduction
Optical beams with orbital angular momentum (OAM), e.g., optical vortices, have been widely used in optical communications and optical manipulation [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13]. The high-order Bessel-Gauss beam with the non-diffractive and self-reconstructing properties was used for the three-dimensional trapping of microparticles [
10]. Recently, unconventional types of optical beams have been proposed and analyzed. For example, spiral-type beams can be generated utilizing the astigmatic transforms [
12], superposition of diffraction-free beams [
13], focusing of shifted vortex beams [
14], spiral toroidal lens [
15], and refractive twisted microaxicons [
16]. A new kind of beam, generated with inseparable helical and conical phases, has been proven to have a spiral-like intensity profile at the focal plane. The beam may find application for optical guiding [
17]. Subsequently, the phase structure and interference characteristics of the helico-conical optical beam have been also analyzed [
18,
19,
20]. In addition, some experts have investigated the self-reconstruction property of the optical beam [
21]. The helical vortex structures are also connected with the depolarization of the laser beams and pulses [
22,
23]. We modified the helico-conical optical beam by adding a power exponent [
24,
25]. The proposed beams have controllable openings along the intensity trajectory, which can be dependent on the power exponent. Nathaniel Hermosa et al. proposed a method to control the intensity patterns of the helical beam, i.e., modifying the phase by boring a hole at the center of the helical phase [
26,
27]. The results show that the area of the bored hole has a great influence on the intensity patterns of the beam. The proposed method can also be introduced into the helico-conical phases. With the advancement of nanotechnology, the modulated beam can potentially be applied in the laser-induced nano-joining of nanoscale materials and the generation of light-induced helical-structured materials [
28,
29]. In isotropic polymers, 3D chiral microstructures can be achieved under the illumination of the spiral lobes and chirality generated by the helical wavefronts [
30].
In this paper, we customized the helico-conical phases to control the length of the helico-conical beams. The customized phases are obtained by eliminating the partial helico-conical phase which is restricted with a filter. We will analyze the intensity patterns of the controllable helico-conical beams in the focal region theoretically and experimentally. Based on the local spatial frequency, we discuss the dependence of the focal field intensity distributions of the controllable helico-conical beams on the filter parameter k. We also analyze the properties of the beam utilizing the energy flow density. The controllable helico-conical beams will find application for the optical guiding and light-induced helical-structured materials.
2. Controllable Helico-Conical Beams
The inseparable helical and conical phase profiles of the helico-conical beam can be written as [
17]
where the parameter
K is equal to 0 or 1,
l is the topological charge which determines the number and direction of the spiral wavefront,
θ is the azimuth angle ranging from 0 to 2
π in the polar coordinate system, and
r0 is the normalization factor of the radial coordinate
r. In this paper, the usable size of the spatial light modulator (SLM) is 1920 × 1080 pixels with a pixel pitch of 8 μm, so the value of the radius
r0 can be set as 4.32 mm and
r can range from 0 to 4.32 mm. The wavelength of the laser beam in the simulations and experiments is set as 532 nm.
Based on Equation (1), the focal-field distributions of the beams can be simulated by the Fourier transform. In this paper, to generate the controllable spiral intensity patterns, we customized the phase hologram generated with the phase function in Equation (1). The bored phase in
Figure 1c,f was obtained by customizing a partial helico-conical region along the screw dislocation of the helico-conical phase, which can be expressed by a matrix in the simulations. The bored phase profile can be written as
, where
is a filter function. The filter can be given by
Thus, the bored phase is also written as
The helico-conical phases with
l = 20,
K = 0, and
K = 1 are shown in
Figure 1a,d respectively.
Figure 1b,e show the filter
F with the parameters
l = 20,
k = 6,
K = 0, and
K = 1, respectively.
Figure 1c,f show the corresponding bored helico-conical phase with
K = 0 and
K = 1, respectively. The variable parameter
k is an integer ranging from 0 to
l. In fact, when the parameter
k is equal to 0, the bored phase profile is the whole helico-conical phase. In the simulation and experiments, we eliminated the central intensity peak by adding the blazed gratings into the holograms in order to analyze the bored intensity distribution better.
3. Intensity Distributions at the Focal Plane
We analyzed the focal-field distributions of the controllable helico-conical beams. For ease of observation, we used normalized intensity for numerical simulations. The focal-field intensity distributions can be calculated by
is the filter function in Equation (2), and is the phase function. It is convenient to evaluate Equation (4) numerically with the FFT algorithm. The intensity distribution can be calculated by .
Figure 2a shows the helico-conical phase profile with
K = 0 and
l = 20.
Figure 2b–e show the bored phase profiles with
K = 0,
l = 20, and
k = 1, 4, 7, and 10, respectively.
Figure 2f–j show the corresponding focal-field intensity distributions.
Figure 3a shows the helico-conical phase profile with
K = 1 and
l = 20.
Figure 3b–e show the bored phase profiles with
K = 1,
l = 20, and
k =1, 4, 7, and 10, respectively.
Figure 3f–j show the corresponding focal-field intensity distributions. The heads of the helico-conical beams shown in
Figure 2f and
Figure 3f can be marked with yellow arrows. The heads of the beams gradually vanish as the bored phase profile increases. It can also be seen from
Figure 2 and
Figure 3 that the spiral intensity trajectories of the controllable helico-conical beams gradually become shorter with the increasing parameter
k.
In this paper, the local spatial frequency distribution was used to analyze the spiral intensity trajectories of the controllable helico-conical beams. Generally, the approximate mapping of the local spatial frequency can be expressed by [
17]
Thus, the local spatial frequency of the controllable helico-conical beams can be described by the equation
and
, respectively [
12]. Based on the phase function in Equation (1) and the filter function
F in Equation (2), the approximate mapping of the local spatial frequency for the beams with
K = 0 and
K = 1 in polar coordinates are expressed as Equations (6) and (7), respectively.
The plots (
,
) of the controllable helico-conical beams with
K = 0,
l = 20, and
k = 0, 1, 4, 7, and 10 are shown in
Figure 4a–e, respectively. The plots (
,
) of the controllable helico-conical beams with
K = 1,
l = 20, and
k = 0, 1, 4, 7, and 10 are shown in
Figure 4f–j, respectively. In
Figure 4a, the points accumulate in a spiral, corresponding to the helico-conical beam with
K = 0 and
l = 20. With
k increasing, the partial region of the phase was bored, as is shown in
Figure 2b–e. The resulting spiral intensity trajectories of the beams became shorter. The inner and outer phases of the helico-conical phase with
K = 0 corresponded to the head and tail of the spiral intensity trajectories, respectively. Contrariwise, the inner and outer phases of the helico-conical phase with
K = 1 corresponded to the tail and head of the spiral intensity trajectories, respectively. From
Figure 4b,g, it also can be seen that an observable dislocation of the spiral’s head occurs. The spot diagrams of the local spatial frequency shown in
Figure 4a–j agree well with the results shown in
Figure 2f–j and
Figure 3f–j, respectively.
The dependence of the focal-field intensity for the controllable helico-conical beams on the parameter
k is analyzed in this paper. The length of the spiral intensity trajectories shown in
Figure 4 cannot be measured simply and directly. Thus, we used the angular dimension to describe the change of the spiral intensity trajectories indirectly. The spiral’s head was treated as the origin, and the endpoint along the spiral trajectory was treated as the tail. The azimuth between the head-to-tail connecting line and the horizontal axis (the
x-axis in
Figure 5) describes the length of the bored spiral intensity trajectories. The calculation of the azimuth is shown in
Figure 5. In the spot diagrams of the local spatial frequency for the controllable beams with
K = 0,
l = 20, and
k = 0 (see
Figure 5a), the corresponding azimuth for the whole helico-conical beam is 3
π/2. As an example, the corresponding azimuth for the controllable helico-conical beams with
K = 0,
l = 20, and
k = 12 is about 0.6712
π, as is shown in
Figure 5b.
With the method shown in
Figure 5a,b, the plots of the azimuth versus the parameter
k for the controllable helico-conical beams with
K = 0 and
K = 1 are shown in
Figure 5c,d respectively. The azimuth calculated with the local spatial spectrum corresponding to the different topological charges
l = 5, 10, 15, 20 is marked with different legends and colors. The linear fitting of the data was implemented, and the corresponding linear regression lines are marked with the corresponding colors. The results demonstrate that the azimuth scales linearly with the parameter
k and the slope of the linear regression line can be dependent on the topological charge
l. The dependence of the azimuth
φ on the parameter
k can be empirically written as
In Equation (6), the magnitude of the slope is and , respectively. The minus sign before the parameter k indicates that the angular dimension decreases as the variable parameter k increases. Thus, we can generate customizable helico-conical beams with the bored helico-conical phases.
We also analyzed the focal-field intensity distributions with the Poynting vector. The energy flow in the focal-field region of the controllable helico-conical beams was calculated. The Poynting vector can be described by the following formula [
1],
where
c,
and
represent light speed, electric field intensity vector, and magnetic field intensity vector, respectively;
,
, and
denote the unit vectors along the
x,
y, and
z directions, respectively.
ω and
k (
k = 2π/
λ) are the angular frequency and the wave number of the controllable helico-conical beam, respectively. The three components in the
x,
y, and
z directions denote the energy flow in the different directions, respectively.
The transversal energy flows (i.e.,
x-
y plane) of the controllable helico-conical beams were calculated and analyzed in this paper.
Figure 6a–j show the focal-field transverse energy flow of the controllable helico-conical with
K = 0,
l = 20, and
k = 0, 1, 4, 7, and 10;
K = 1,
l = 20, and
k = 0, 1, 4, 7, and 10, respectively. The direction and magnitude of the green arrows in the figures demonstrate the direction and size of the energy flow at the Fourier transform plane. From the figures, it can be seen that the energy flow is flowing towards the interior of the bored helico-conical beam, which demonstrated that the beam can trap particles and bind particles towards the spiral trajectories.
The controllable helico-conical beams also carry the OAM at the focal plane. The focal-field OAM density of the controllable helico-conical beams in free space can be described by
where,
E and
B denote the electric and magnetic fields, respectively;
r = (
x2 +
y2)
1/2,
Sx and
Sy are the Poynting vector along the
x and
y directions, respectively.
Figure 7a–e demonstrate the focal-field OAM density distributions of the controllable helico-conical beams with
K = 0,
l = 20, and
k = 0, 1, 4, 7, and 10, respectively.
Figure 7f–j illustrate the focal-field OAM density distributions of the controllable helico-conical beam with
K = 1,
l = 20, and
k = 0, 1, 4, 7, and 10, respectively. The OAM density distributions of the controllable beams at the focal plane are roughly consistent with the focal-field intensity distributions. With the missing part of the helico-conical phase profile increased, the OAM density distributions are totally changed and follow the bored spiral lines. The OAM density also decreased gradually. We analyzed the relationship between the normalized OAM density and the parameter
k. The plots of the normalized OAM density versus the parameter
k are shown in
Figure 8. When
K = 1, the focal-field OAM density of the controllable helico-conical beam is relatively large. However, as the parameter
k increased, the OAM density decreased. When the parameter
k is larger (e.g.,
k = 20), the whole helico-conical phase profile is missing basically and the corresponding OAM density is close to 0. Thus, the helico-conical beam can be modulated by considering the above results and can be applied to optical trapping and light-induced helical-structured materials.
In this paper, the experimental focal field intensity patterns of the controllable helico-conical beams have been analyzed.
Figure 9a shows the schematic setup for generating the controllable helico-conical beams. A collimated and expanded laser (
λ = 532 nm) beam impinged onto the SLM (reflective type), which was encoded with a computer-generated hologram. The flat convex lens
L1 (
f1 = 30 mm) and
L2 (
f2 = 200 mm) was used for the collimation and expansion of the optical beams. The lens (
L3,
f3 = 150 mm) was a Fourier-transform lens. The intensity cross-sections of the desired beams were captured by the CCD camera. During the experiments, the corresponding phase profiles in
Figure 2a–e and
Figure 3a–e were loaded on the SLM, respectively. The corresponding intensity patterns at the focal plane are shown in
Figure 9b–k, respectively. The experiment results are consistent with the simulated ones shown in
Figure 2f–j and
Figure 3f–j, respectively. The focal field intensity distributions were determined by the parameter
k.