Next Article in Journal
Fuzzy Logic Applied to Sustainable Development Goals and Human Trafficking
Previous Article in Journal
Generation of Julia and Mandelbrot Sets via Fixed Points
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computer-Aided Laser-Fiber Output Beam 3D Spatial and Angular Design

by
Rocío Pérez de Prado
1,*,
Sebastián García-Galán
1,
José Enrique Muñoz-Expósito
1 and
Adam Marchewka
2
1
Telecommunication Engineering Department, University of Jaén, Science and Technology Campus, 23700 Linares (Jaén), Spain
2
Institute of Telecommunications and Informatics, University of Technology and Life Sciences, Prof. S. Kaliskiego 7, 85-791 Bydgoszcz, Poland
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(1), 83; https://doi.org/10.3390/sym12010083
Submission received: 15 November 2019 / Revised: 21 December 2019 / Accepted: 30 December 2019 / Published: 2 January 2020

Abstract

:
Multiple laser beams and single-mode optical fibers output can be approximated by assuming that the emitted light has a symmetrical Gaussian intensity profile, which corresponds to the transverse electromagnetic mode (TEM00), which is designated as a Gaussian beam. Current free-accessible design tools are limited to the spatial analysis of the beams, in general, and to the intensity, in particular, and to the graphical visualization in 2D with very limited options. In this work, a novel a computer-aided laser-fiber output beam TEM00 designer, CATEM00, is presented based on the 3D representations typically provided by camera beam profilers, and on the fundamentals of the wave theory of light, including diverse flexibility capabilities for graphical manipulation and parameter comprehension both in terms of spatial behavior and in angular confinement. It must be highlighted that not only is the spatial limitation design of light impact relevant in TEM00 applications but, also, the angle with which the light reaches the target. Hence, the availability of capabilities of phase design in TEM00 following the paraxial limitations is highly convenient. Results and discussion in terms of intensity, power, divergence and wave fronts are presented considering a set of study cases, showing its coherence with Gaussian beam theory.

1. Introduction

Most laser applications consist of the impact of light on objects or materials on a given limited surface with a striking angle and spatial distribution of power strictly bounded. The objects to impact can be both light detection equipment for information systems, such as cameras, or industrial or biological materials to be drilled, melted, or heated [1]. The quality of the performance of the light impact highly depends on the features of the beam of the light source and the way these features change across the beam and along the optical axis as it propagates in free space. Particularly, the transverse electromagnetic mode (TEM00) describes the observed output of many lasers, which are designed and tuned to emit this mode, and it also approximates the light of the output of a single-mode optical fiber [2]. The ideal TEM00 beam presents an intensity profile that corresponds to a Gaussian expression and, thus, it is also known as a Gaussian beam [3]. In TEM00 mode, the light beam presents an intensity transverse section that is maximum in the center and that it characterized by its revolution symmetry around the optical axis. Analogously, the angular deviation of the beam reaches its lowest value in the optical axis and it is increased in its surroundings, ruled by the paraxial approximation, also presenting symmetry with respect to this axis. Hence, both the space and angular parameters of the light beam, i.e., intensity, the power distribution and divergence, and wave fronts, respectively, of the TEM00 profile are featured by its symmetry around the optical axis, which is convenient and demanded in most related applications.
TEM00 mode with Gaussian beam applications are included within those of paraxial waves. Paraxial waves are of relevance in applications where it is necessary to achieve the propagation of light in free space not only with a spatial confinement, but also with an angular confinement. Spatial confinement refers to the characteristic of the wave of being limited to a specific region of space while the angular confinement corresponds to the condition of the normal directions to the wave fronts that imposes small angles with the optical axis or axial axis z. In this way, the beam power is concentrated in a small cone around the optical axis. There are different parameters describing the light beam both in terms of spatial and angular behavior, joined to different concepts of intensity, power, divergence and wave fronts. In diverse application areas ranging from medicine, e.g., ophthalmic laser surgery [4], to industry, e.g., fusion and material drilling [5], and military applications, e.g., weapons and signals detectors, the correct understanding of the beam parameters is critical for its design and use. Furthermore, the precise understanding of the light source parameters associated to the performance of the laser source is relevant to detect the need for tuning and profiling of the beam equipment, which is typically related to the measurement of symmetry and smoothness [1].
There are diverse computer-aided tools for TEM00 beam features that have empowered researchers and engineers to design lasers to the highest quality. Particularly, in regard to TEM00 beam computer-aided open-source designers, there exist many beam software packages offering analytical and visual results of diverse beams featuring parameters such as spot size, beam waist, and divergence on the basis of main characteristics of the light source performance, such as frequency and Rayleigh distance. This is the case of [6,7,8], which are the major open-accessible software in the field to the best of the authors’ knowledge. Nevertheless, these tools present different deficiencies that reduce the manageability, intuition, and versatility. Mainly, these deficiencies arise from the limitation to spatial analysis of the beam, typically to intensity features, and to the static visualization of its profile, including non- or very reduced flexibility for graphical design in 3D and 2D. In this sense, [6] just offers the possibility to visualize in 2D the Gaussian beam waist in different positions of the optical axis. Additionally, in this regard, [7] allows the visualization of the Gaussian beam only in the spatial dimension, and not in the angular one. Moreover, the majority of available tools simply consist of analytical calculators for power and intensity study and no graphical user interfaces for both visual and quantitative analyses are found. In this sense, for example, [8] only offers the analytical end value of the power of the Gaussian beam as a function of position and distance.
In this work, a computer-aided laser-fiber output beam for 3D spatial and angular design tool, CATEM00, is developed in MATLAB [9], integrated in a simplified graphical user interface (GUI), following recent published applications in optics and electronics [10,11]. Additionally, CATEM00 is made accessible to the scientific community in [12] as well as in the Supplementary Files attached to this paper. The contribution of this work is highlighted as follows:
  • Spatial and angular beam TEM00 design. The previous works that can be currently found in the same area mainly focus on the analysis of the intensity of the Gaussian beam. The proposal in this work offers the possibility to consider a design based on spatial characteristics, such as intensity, power, divergence, as well as on angular characteristics. In particular, the beam can be configured based on wave fronts, and, therefore, it performs not only considering the intensity or power that has to impact on a target zone, but also with the angle with which it will impact on said target zone. The relevance of simultaneous design in both spatial and angular features allows a complete characterization of the beam, which is necessary in high-precision applications.
  • Intuitive and flexible GUI. CATEM00 allows the design of TEM00 beams based on the complex wave theory intuitively, both in the comprehension and introduction of the parameters that govern said theory, and in the visualization and interpretation of the results. Thus, the tool allows a direct and simplified design based on the key parameters in both the spatial and angular dimensions, such as the Rayleigh distance and the separation to the optical axis. In addition, the graphical user interface offers the necessary flexibility in determining the area around the optical axis where the paraxial wave approximation is considered. Only through this option an adapted beam design can be made for each application in terms of angular characteristics.
  • Visual 2D and 3D design: Unlike other Gaussian beam design tools, the proposal of this work allows simultaneous design in both 3D and 2D following scaled intensity contour schemes for the measured parameters of intensity, power and divergence. In this way, the design result is visually assimilated with the outputs of beam profiling tools for light sources tuning in the operational phase, as obtained, for example, in the case of intensity, with cameras of beam profiling like Spiron SP907 and BeamMic software from Ophir, as shown in Figure 1 [13]. In addition, on existing tuning tools, the visualizations can be rotated at the discretion of the designer and be arranged with the choice of azimuth, elevation, and convenient distance in each application.
  • Quantitative results extracted from visual representations of the beam and supporting analysis. The quantitative value at each point of representation of each parameter in design can be extracted directly by clicking on it. In addition, based on the chosen parameters, for the design of the beam in intensity, its characterization in far and near-field is provided, and also for power design, the analysis of power concentrated in a fixed section in the space is displayed. These diagrams are intended to facilitate the choice of input data based on their variation in ranges.
  • Freely-available tool. As pointed out above, the proposal is accessible in [12] as well as in the Supplementary Files attached to this paper, with the purpose of facilitating engineers and researchers in the area tools to improve and facilitate the design and comprehension of the spatial and angular design of TEM00 mode.
Finally, it must be highlighted here that quality of any laser-fiber output defined by a TEM00 is directly related to the symmetry of its spatial and angular features, whose guarantee and knowledge is the end purpose of this work. This comes from the fact that the expected symmetry to be found in the impact target in commented critical applications such as ophthalmic laser surgery, strictly depends on the symmetry of the light source beam of impact and thus, designing and testing the beam according to symmetry final outcome is extremely relevant. Hence, with this work, authors aim to make a further step in the symmetry research in the optical field.
The rest of the paper is structured as follows: In Section 2, the main principles and definitions related to the TEM00 mode are presented and the description of the fundamentals of both spatial and angular dimension for its design are introduced. Section 3 is devoted to the explanation of the proposal, including both the tool performance essentials and the GUI main characteristics and properties. Next, in Section 4 the results and associated discussion are provided. Specifically, Section 4 differentiates the results and discussion in terms of intensity, power, divergence, and wave fronts, and multiple examples are analyzed to validate the performance of the tool according to the wave theory for paraxial waves, in general, and Gaussian beams, in particular. Finally, in Section 5 the main conclusions of the work are drawn.

2. Materials and Methods

2.1. Transverse Electromagnetic Mode TEM00 Principles and Definitions

The complex amplitude of the Gaussian beam is obtained from the wave theory of light as a particular case of the complex amplitude of paraxial waves [14]:
U ( r ) = A ( r ) exp ( j k z )
where k represents the wave number, z corresponds to the optical axis-z and A(r) satisfies:
T 2 A j 2 k A z = 0
with the following particular form:
A ( r ) = A 1 q ( z ) exp [ j k ρ 2 2 q ( z ) ] ,   q ( z ) = z ξ
with ρ the distance to the optical axis, q(z) the so-called q parameter and A1 a constant value. This expression can be regarded as a space-shifted version of the complex envelope for parabolic waves. Specifically, the shift corresponds to ξ = j z 0 , and so the complex envelope of the Gaussian beam can be defined as:
A ( r ) = A 1 q ( z ) exp [ j k ρ 2 2 q ( z ) ] ,   q ( z ) = z + j z 0
where z0 denotes the Rayleigh range. Separating the real and imaginary parts of 1/q(z):
1 q ( z ) = 1 R ( z ) j λ π W 2 ( z )
the key parameters of the Gaussian beam can be deduced, being λ the wavelength. On the one hand, it can be inferred that the width or radius of the Gaussian beam corresponds to:
W ( z ) = W 0 [ 1 + ( z z 0 ) 2 ] 1 / 2
where W0 describes the beam waist or radius of the light source:
W 0 = ( λ z 0 π ) 1 / 2
Twice the value of the beam waist is known as the spot size or beam initial diameter, and it represents the diameter of the light source. The width of the Gauss beam determines the intensity spreading or spatial confinement from the optical axis.
Figure 2 schematically represents the propagation of the TEM00 mode as a function of distance to the light source in the optical axis-z in both spatial and angular dimension. As it can be observed, in the spatial dimension, the form of the intensity corresponds to a Gaussian curve and, thus, the TEM00 laser beams are known as Gaussian beams. Additionally, it can be appreciated that the boundaries of the beam are located around the optical axis, and expands as the beam gets further along in it. The beam waist W0 and the beam width W(z) are also depicted in Figure 2. As described by Equation (6), the beam diameter of the TEM00 mode is defined as the separation from the center of the beam, where the maximum intensity is always achieved, to the point where the power percentage corresponds to 86.47% approximately. In addition, it can be appreciated in Figure 2 that the beam waist corresponds to the beam width at the light origin.
On the other hand, in regard to the angular definition of the TEM00 beam, the curvature radius of the Gaussian beam it is given by the expression:
R ( z ) = z [ 1 + ( z 0 z ) 2 ]
which describes the angular confinement of the beam. In Figure 2, a representation of the wave fronts is presented, as an introduction to their definition, studied in the next section. As shown, the curvature radius R(z) represents the radius of the obtained fronts following the paraxial approximation. Hence, as it can be appreciated, the radio of curvature and, so, the angular confinement of the TEM00 mode, decreases from the light source up to the Rayleigh distance z0. From the Rayleigh distance, the angular confinement increases. In this way, the wave front evolves from confined plane wave fronts (with infinite value for R(z)) to limited spherical wave fronts (with R(z) = z), that is, confined plane waves in the far-field (with infinite value for R(z)).
In order to easily visualize that the complex amplitude of a Gaussian beam follows a Gauss form it is convenient to substitute q(z) in terms of its module and phase in the amplitude term of Equation (4), and q(z) in terms of W(z) and R(z), i.e., in terms of Equation (5), in the phase term of Equation (4). Hence, the following expression for the complex envelope is obtained:
U ( r ) = A 0 W 0 W ( z ) exp ( ρ 2 W 2 ( z ) ) exp ( j k z j k ρ 2 2 R ( z ) + j ζ ( z ) )
with A0 a constant value. From this expression all the properties of the Gaussian beam can be deduced, as it will be summarized as follows.

2.2. Spatial Dimension

First, the beam intensity expression can be obtained from the complex amplitude [14]:
I ( r ) = | U ( r ) | 2
and substituting Equation (9) in Equation (10):
I ( ρ , z ) = I 0 [ W 0 W ( z ) ] 2 exp ( 2 ρ 2 W 2 ( z ) )
with I 0 = | A 0 | 2 , where it can be observed that the beam intensity of a Gaussian beam follows a Gaussian expression.
Second, the power is defined as the integral of the optical intensity in a transverse circular plane at a given distance z of infinite radius:
P = 0 I ( ρ , z ) 2 π ρ d ρ
where substituting Equation (11) gives us the following expression for the power of the Gaussian beam:
P = 1 2 I 0 ( π W 0 2 )
Considering this expression, Equation (11) can also be expressed in terms of the power of the Gaussian beam:
I ( ρ , z ) = 2 P π W 2 ( z ) exp ( 2 ρ 2 W 2 ( z ) )
Also, from the power in a circular cross-section of radius ρ0, the power fraction on this selected section can be calculated as:
1 P 0 ρ 0 I ( ρ , z ) 2 π ρ d ρ = 1 exp ( 2 ρ 2 0 W 2 ( z ) )
Third, the divergence or spatial confinement of the beam is given by the beam width described in Equation (6). Specifically, the divergence in far-field corresponds to:
W ( z ) W 0 z 0 z = 2 π λ 2 W 0 = θ 0 z    : z > > z 0
where θ 0 corresponds to the divergence parameter. Particularly, it is convenient to study the divergence of the beam considering the revolution symmetry of W(z) following Equations (6) and (16) with respect to the optical axis. Moreover, in relation to the divergence of the beam it is relevant to point out the concept of depth of focus, which represents the capacity of the beam to focus light and it corresponds to twice the value of the Rayleigh distance:
2 z 0 = 2 π W 2 0 λ

2.3. Angular Dimension

From the complex amplitude of the Gaussian beam, Equation (9), the phase of the Gaussian beam can be easily deduced [14]:
φ ( ρ , z ) = k z ζ ( z ) + k ρ 2 2 R ( z )
where it can be inferred that in the optical axis, the phase corresponds to that of plane waves with a slight modification. Wave fronts are the surfaces of constant phase that satisfy:
k [ z + ρ 2 2 R ( z ) ] ζ ( z ) = 2 π q
which, in general for all the points in z, closely describes parabolic surfaces with radius of curvature R(z).

2.4. Incorporating TEM00 Deviations Into the Propagation Model

There exist deviations from theoretical Gaussian beams in the output of real-life lasers and fibers that are also considered in the proposed tool [15,16,17,18]. For example, low-power light beams from helium neon lasers can be a close approximation to TEM00, but it must be noted that as the power of the light source becomes higher, the pumping mechanism becomes more complex and the order of the emitted mode becomes higher and thus, the more the beam deviates from the ideal TEM00. Indeed, although the beam from a well-controlled helium neon laser and other gas lasers is very close to that of a theoretical TEM00 mode, for most lasers, even specifying a TEM00 mode, the output contains some components of higher-order modes that do not propagate according to the theoretical TEM00. Hence, it is necessary to introduce in the proposed tool these deviations from the theoretical Gaussian beam.
Particularly, to accommodate this variance, the M2 quality factor (known as the M-squared factor) is generally used in the Gaussian light field to describe the deviation of the laser beam from a theoretical Gaussian beam or theoretical TEM00, and it has also been considered in the presented proposal. This factor can be deduced from the beam parameter product (BPP), which is the result of the TEM00 beam divergence θ 0 and waist size W(z). In particular, the M2 can be obtained measuring the BPP of a real laser beam (product of the beam’s minimum diameter and far-field divergence) and dividing this value by the BPP of an ideal TEM00 at the same frequency. Hence, the M-squared factor or M2 factor is a long-recognized figure of merit used to determine the deviation of the laser beam from a theoretical Gaussian TEM00 that can be defined as:
M 2 = W 0 R θ 0 R W 0 θ 0
where W 0 R and θ 0 R are the beam waist and far-field divergence of the real beam, respectively. The M2 factor for a theoretical Gaussian beam is equal to the unity, M2 = 1 and the beam-waist and beam-divergence product is given by:
W 0 θ 0 = λ π
It follows then that for a real laser beam or deviated TEM00:
W 0 R θ 0 R = M 2 λ π
All real laser beams have M2 values greater than one (M2 > 1), although very high-quality beams present values very close to the unity. In all cases, the M2 factor affects the characteristics of a laser beam and cannot be neglected in optical designs. Specific values for M2 can be found in many common TEM00 light sources. For instance, for helium neon lasers M2 is generally lower than 1.1 (M2 < 1.1). For ion lasers, the M2 typical range goes from 1.1 to 1.3 (1.1 < M2 < 1.3). In addition, collimated TEM00 diode laser beams usually have an M2 ranging from 1.1 to 1.7 (1.1 < M2 < 1.7) and in high-energy multimode lasers, M2 can be as high as 25 or 30 (25 < M2 < 30).
Once the deviations from the theoretical TEM00 have been analyzed, the propagation model must be redefined accordingly. Specifically, the propagation equations for a real laser beam must now consider the following expressions for the beam radius:
W R ( z ) = W 0 R [ 1 + ( z z 0 R ) 2 ] 1 / 2 = W 0 R [ 1 + ( z λ M 2 π W 2 0 R ) 2 ] 1 / 2
and the radius of curvature:
R R ( z ) = z [ 1 + ( z 0 R z ) 2 ] 1 / 2 = z [ 1 + ( π W 2 0 R z λ M 2 ) 2 ]
where the Rayleigh distance is transformed to:
z 0 R = π W 2 0 R λ M 2

2.5. Laguerre–Gaussian Modes

The TEM00 mode can be regarded as a particular case of Laguerre–Gaussian modes [19]. When a laser light beam has a circularly symmetric profile or its associated cavity is cylindrically symmetric (as in the case of optical fibers too) the propagation model is efficiently described by the Laguerre–Gaussian modal decomposition. The related expressions are expressed in cylindrical coordinates considering generalized Laguerre polynomials. Particularly, each TEM mode is denoted using two integers p (positive or zero) and l, representing the radial index and the azimuthal index, respectively. The general case corresponds to TEMpl modes and to TEM00 when p and l are zero. The specific form of the complex amplitude for the general case of Laguerre–Gaussian TEMpl modes corresponds to the following expression:
U ( r ) = C L G l p W 0 W ( z ) exp ( r 2 W ( z ) ) | l | exp ( r 2 W 2 ( z ) ) L p | l | exp ( 2 r 2 W 2 ( z ) ) × exp ( j k r 2 2 R ( z ) ) exp ( j l ϕ ) exp ( j ψ ( z ) )
where L p | l | indicates the generalized Laguerre polynomials, C l p L G corresponds to a normalization constant, W(z) and R(z) were described in Equations (6) and (8), respectively, for the TEM00, and the Guoy phase along z is:
ψ ( z ) = ( N + 1 ) tan 1 ( z z 0 )
with N = |l| + 2p. As it can be inferred, this expression contains the TEM00 complex amplitude discussed above.
The physical effect of the rotational mode number l described by Equation (26), beyond the contribution to the Laguerre polynomials L p | l | , is essentially given by the phase factor exp ( j l ϕ ) , which determines that the light beam is advanced (or retarded) by a 2πl phase shift in a single rotation around the beam (in ϕ ). This can be regarded as an optical vortex of topological charge l that is related to the orbital angular momentum (OAM) of light in that TEMpl mode. The light in an optical vortex is twisted around the propagation axis z and due to this twisting, the light waves at the optical axis are cancelled by the combined effect. As a result, an optical vortex typically presents a ring form with a black hole in the central space when projected on a flat surface. The topological charge l specifies the integer number of twists of light for the same frequency, which can be a positive or negative number depending on the sense of rotation. It must be noted that the greater value of l, the greater the number of twists and thus, the faster the rotation or spinning around the propagation axis z. The twisting of light around the optical axis brings OAM with the optical mode, inducing torque on electric dipoles. It can be observed that OAM differs from the well-known spin angular momentum that generates a circular polarization, and it can be appreciated in the orbiting motion of trapped particles. When an optical vortex and a light plane wave interfere the phase is shown as concentric spirals with topological charge equal to number of arms in the spiral. Hence, optical vortices and, thus, OAM, are directly related to the generation of Laguerre–Gaussian TEM modes.
As introduced above, TEM00 and higher order modes, that is, TEMpl modes, are given by the Laguerre–Gaussian modal decomposition of order (p,l). Hence, TEM00 is a particular case of these modes when the order corresponds to (0,0). When l = 0, the OAM of the Laguerre–Gaussian mode is zero, exp ( j l ϕ ) = 1 , and thus, the light is not twisted around the optical axis as it propagates in the optical axis. In contrast, when l > 0, the light is rotated around the optical axis as it propagates, given the dependence of the mode with ϕ ruled by l. It must be pointed out that the consideration of OAM of modes has a large set of emerging applications in areas such as quantum computing (e.g., encoding and storage of information) and optical communications (especially for multiplexing in high-capacity fiber networks) [20,21,22].

2.6. TEM00 Photon Statistics

Photon statistics describe the probability distribution p ( n , T ) , to detect n photons in the time interval T [14,23]. The photon statistics of a TEM00 laser beam quantum mechanically correspond to that of a monochromatic field of a single-mode laser described by the coherent states | α :
| α = e 1 2 | α | 2 n α n n ! | n
Hereby, the intensity verifies:
I = n ^ = | α 2 | = n ¯
The probability to detect a specific number of photons is given by:
p ( n ) = | n | α | 2 = | α | 2 n n ! e | α | 2 = n ¯ n n ! e n ¯
which corresponds to a Poisson-distribution with mean n ¯ and variance:
Δ n 2 = n ¯
Finally, it must be pointed out that, contrarily to the classical description, the laser light reveals fluctuations. These fluctuations are due to the particle’s nature of the light that generates shot noise. For a very high number of photons, the relative fluctuations Δ n 2 / n ¯ 2 tend to zero. Thereby, the classical theory can be considered as a limit case.

3. Proposal

The presented computer-assisted TEM00 mode profiling system, CATEM00, proposes a design based on two of the basic characterization dimensions of light beams. On the one hand, the light beams are designed in the spatial dimension and, on the other hand, they are designed in the angular dimension. As has been studied, the TEM00 beams are characterized by spatial and angular confinement, i.e., paraxial waves. This confinement allows the generation of bounded beams both in area and in angle of impact, so relevant in most of the applications in which they are used. Hence, the suggested assisted design system offers first the possibility of beam profiling in spatial terms, considering intensity, power, and divergence and, secondly, the analysis of angular confinement based on wave fronts. In addition, to facilitate the design, the beam design software is integrated into an intuitive graphical user interface.

3.1. Spatial and Angular Beam Computer-Aided Designer: CATEM00

Both intensity and power of TEM00 beams are analyzed from the point of view of the evolution of the optical axis distance, the Rayleigh distance and the frequency of the light source. In addition, for the intensity configuration, the possibility to study the variation in the optical axis of the maximum intensity in consideration not only of the exact expression, but also of the far-field version, is included. Moreover, in the case of the design of the optical power, it is possible to analyze the evolution of the power fraction against the width of the Gaussian beam and obtain said power fraction in a fixed section whose radius coincides with the selected width. In relation to the divergence, the computer profiling system allows modeling the divergence cone of the TEM00 mode considering both its general version and its useful far-field version. The exact divergence and far-field cones are represented by performing the revolution symmetry around the optical axis of the Gaussian beam width and the divergence parameter, respectively.
On the other hand, for the design of the angular confinement, the wave fronts of the TEM00 mode are profiled from the point of view of Rayleigh distance to the optical source and the frequency of light emission. It is important to note that the design of the wave fronts must respect the paraxial approach of the Gaussian beam to ensure the validity of the modeling. Finally, as a way of introducing the modeling of the TEM00 and comparative mode in terms of spatial and angular confinement, a section is proposed within the flat and spherical wave characterization design system. Figure 3 shows the scheme of the design areas of the computer-assisted profiler.
Furthermore, Figure 4 shows the flow chart of the operation of the presented tool CATEM00. First, to carry out the design, it is necessary to specify the parameter to be optimized, i.e., the spatial characterization parameter: intensity, power. or divergence, and angular characterization parameter, i.e., the wave fronts or the study of fundamental waves. Once the parameter/study to be analyzed has been selected, the input information must be entered into the design tool. This information includes the frequency, the amplitude, the Rayleigh distance, the distance in the optical axis, the radius of the fixed section and the distance around the x-y axes where it is desired to carry out the paraxial approximation for angular confinement. Once the input information has been entered correctly in format, the beam designer can be executed and the obtained TEM00 beam profile can be observed. If no further analysis of the beam is desired, the tool allows to re-design based on either the selection of the parameter to be designed, the desired analysis or the introduction of parameters for the current design. In addition, it is possible to save the obtained design including the graphic results and the analyzed parameter, the selected analysis and the specific input data. Additionally, the tool provides a more detailed characterization of the obtained results, allowing a rotation of the 3D profile in distance, azimuth and elevation and extracting qualitative data from any point of the graphs, inserting legends or intensity bars to allow a better visual compression of the strength or intensity of the parameter in every point. In addition, after each detailed design analysis, the tool allows one to perform a new analysis, to save the data, or to increase the current results’ analysis.

3.2. Graphical User Interface

The beam designer is integrated into a graphical user interface to allow greater ease of use of the tool and interpretation of results. As Figure 5 shows, this interface consists of two fundamental panels. On the one hand, the main panel integrates the parameters’ selection, analysis, input, execution, main visualization and, on the other hand, the auxiliary panel allows the plan and elevation visualization of the main visualization. The visualization in plan and elevation is the default configuration that is given to the 3D representations, although it is also possible to configure the desired visualization dynamically by means of the in-depth analysis options of the tool. In addition, as can be seen in the upper left corner of the graphic interface, there is a tool menu for the detailed analysis of the results.
Figure 6 illustrates the parameter selection and analysis panels with all available options. As it can be observed, each parameter study is associated with various options in the analysis panel including all the choices specified in Figure 3 regarding intensity, power, divergence, and wave fronts. Thus, in total, the spatial and angular beam designer allows working with 16 profiling options. Additionally, Figure 7 shows the data entry panel that, as introduced in the previous section, is composed of the information: frequency, amplitude, Rayleigh distance, distance on the optical axis, the radius of the fixed section, the distance around the x–y axes, and the M2 factor. It should be noted that, depending on the type of parameter and beam analysis to be designed, the tool automatically enables those parameters and, furthermore, the factor that must be adjusted successively to achieve optimization is highlighted. The launch panel for the execution of CATEM00 is also shown in Figure 7.

4. Results and Discussion

4.1. Intensity Design

As explained, the TEM00 profiler allows the profiling of Gaussian beam intensity in four different aspects. The first characterization is the relation of the intensity with the distance to the light source or optical axis location z. In Figure 8, it can be appreciated the output of the profiler for the Gaussian beam intensity in two different positions of the optical axis (1 m and 10 m) considering the Rayleigh range and the frequency fixed. As observed, as the light advances in the optical z-axis, the Gaussian beam decreases its focality, i.e., light is spread, corresponding to an increment in the beam with W(z). As a function of the distance to the optical axis ρ = (x2 + y2)1/2 it can also be visualized that for each value of z, the intensity depends on the radial distance ρ.
In ρ = 0 (z axis), the intensity reaches its maximum and it decreases as ρ increases. Furthermore, it can be analyzed that this spreading effect is wider as the beam separates from the light source. It is important to realize that W(z) increases when advancing in the optical axis z and so does the intensity spreading along the transverse radio (as pointed out above, the definition of W(z), Equation (6), is related to the intensity spreading).
Additionally, the proposed Gaussian beam profiler allows the simulation and study of the intensity in relation to the depth of the focus (2z0) or Rayleigh distance (z0) of the light source. Figure 9 presents three outputs of the proposed profiler, representing the Gaussian beam intensity in a fixed position of the beam axis and a fixed frequency varying the Rayleigh range (0.15 m, 5 m, and 80 m). It can be studied that if z0 decreases or increases from the considered distance in the optical z-axis, the focality of the light source is decreased, which means that the resulting Gaussian beam intensity profile is the narrowest when z0 = z. Hence, in z0 the highest focality is achieved for the TEM00 mode.
Moreover, the Gaussian beam intensity in a fixed position of the beam axis and for a given Rayleigh range varying the frequency can be analyzed with the TEM00 beam profiler. As observed in Figure 10, if the frequency increases, the intensity profile gets narrower, which means that if the frequency is higher, the intensity is more confined around the optical axis. From the definition of the beam waist or spot size of the light source (Equation (7)), the focality (related to depth of focus, 2z0) becomes greater when W0 (half of the spot size) increases or when wavelength decreases (increment in frequency value).
Furthermore, the Gaussian beam intensity in the beam axis z considering a given frequency and Rayleigh range against a wide range of distance values z (2D representation) can be obtained with CATEM00. The representation contains a marker point indicating the Rayleigh range and include the approximate intensity or trend that the Gaussian beam intensity follows for far-field conditions (typically 10 times the Rayleigh range).
On the one hand, the exact intensity in the axis corresponds to Equation (11) when ρ = 0:
I ( 0 , z ) = I 0 [ W 0 W ( z ) ] 2 = I 0 1 + ( z / z 0 ) 2
As studied, the Gaussian beam intensity is a function of the distance to the optical axis ρ = (x2 + y2)1/2. Hence, from Equation (32), it can be deduced that the intensity in the z axis ρ = 0 reaches its maximum when z = 0:
max [ I ( ρ , z ) ] = I ( 0 , 0 ) = I 0
and it decays as z increases, reaching half of the maximum intensity value at the Rayleigh range z = z0:
I ( 0 , z 0 ) = I 0 2
Figure 11 shows that the beam intensity at points in the beam axis (ρ = 0) decreases gradually as the distance to the focus increases, reaching half its peak value at z = ±z0. It can also be appreciated that the overall peak intensity I(0,0) = I0 occurs at the beam center (z = 0 and ρ = 0). On the other hand, the approximate expression for intensity in far-field corresponds to:
I ( 0 , z ) = I 0 [ W 0 W ( z ) ] 2 = I 0 1 + ( z / z 0 ) 2 I 0 z 0 2 / z 2 : z > > z 0
Thus, it can be visualized that the intensity decreases with the distance following an inverse-square law in far-field conditions.

4.2. Power Design

It can be deduced from Equation (13) that the Gaussian beam total power is independent of z (in the absence of attenuation factors) and it corresponds to half of the maximum intensity multiplied by the beam area in the waist position. Nevertheless, a key aspect in the Gaussian beam power analysis resides in its spatial confinement as it propagates along the optical axis. The spatial confinement of the power of the Gaussian beam can be studied with the CATEM00 regarding its relation to the optical distance in z, the Rayleigh distance, and the frequency, as also considered previously for the intensity.
Figure 12 shows the power fraction in the laser position for all x and y directions for different positions in the z-axis and fixed values for the Rayleigh range and the frequency. As observed, the profiler indicates the power circular sections (described by its radius from the optical axis) that contain the diverse power fractions (e.g., power fraction equivalent to 1), and it can be appreciated that the associated radius increases as the Gaussian beam is studied in farther distances from the light source for all the power fractions.
A similar conclusion can be derived from the study of the power fraction in the light source position for all x and y directions for different values for the Rayleigh range and fixed frequency and position in the optical axis, given by the computer-aided tool. As illustrated in Figure 13, the power fraction confined in a fixed section area increases or decreases depending on the fact that the Rayleigh distance is higher or lower than the considered distance in z, achieving the highest spatial confinement when z0 = z.
Additionally, Figure 14 shows the performance of CATEM00 in the analysis of the power fraction in relation to the frequency. Specifically, the power fraction in the laser position for all x and y coordinates (wide range) for different values for the frequency and fixed Rayleigh distance and position in the optical axis is studied. It can be inferred that, analogously, to the previous intensity versus frequency analysis, the higher the frequency is, the higher the power fraction contained in a predefined section.
Moreover, from the computer-aided profiler it can be deduced that when the considered circular transverse section radius is ρ0 = W(z) and ρ0 = 1.5W(z), the contained power is approximately 86.47% and 99%, respectively. Figure 15 shows the power fraction for different values of the beam width W(z). A graphical indication of the power fraction that corresponds to a beam width equal to a selected cross-section is also presented.

4.3. Divergence Design

Regarding the divergence analysis, the beam profiler offers the possibility to study the exact divergence cone and the divergence cone in far-field conditions. Firstly, Figure 16 shows the exact divergence cone for different frequencies values. As known, the exact divergence cone is obtained by the revolution symmetry of the beam width W(z), given by Equation (6), and it contains 86.47% of the emitted power by the light source, approximately. Considering this expression, it can be deduced that the exact divergence cone gets narrower as the frequency is increased.
Secondly, as it can be appreciated in Figure 17, the profiler can also feature the approximate divergence cone for far-field conditions. In this approximation, given by Equation (16), the divergence increases linearly with z, defining a cone with half angle θ 0 that approximately contains 86.74% of the power. As observed, the divergence is directly proportional to the relation λ 2 W 0 , from which it can be deduced that in order to obtain very directional beams, short wavelengths, or large spot sizes are necessary. Thus, as it can be appreciated in Figure 17 that the higher the frequency (i.e., shorter wavelengths), the lower divergence. Analogous conclusions can be obtained from the analysis of the depth of focus with the proposed tool. The depth of the focus is proportional to the beam area at the light source location and it is inversely proportional to the wavelength, as introduced in Equation (17). Hence, a large depth of focus, or high spatial confinement of light along the z-axis, requires large spot sizes or short wavelengths.

4.4. Wave Fronts Design

Finally, the TEM00 profiler allows the characterization of the angular confinement. In this sense, it has to be remembered here that the Gaussian beam model is founded on the paraxial approximation and so, it is only valid when the wave fronts are tilted less than 30° approximately from the axis beam. This is translated to the fact that the TEM00 wave fronts given by the Gaussian beam model is precise when the beam waist W0 is larger than 2λ/π, or equivalently, for the radius to the optical source much greater than z, specified by the paraxial approximation of waves. Hence, the presented tool offers the possibility to change the distance to the optical axis, to set up the visualization of the wave front model respecting the paraxial approximation. Equation (19) determines the surfaces of the TEM00 wave fronts, where the radius of curvature R(z) of the wave fronts of a Gaussian beam helps in the analysis of the angular confinement. On the one hand, in z = 0, R(z) = ∞, which corresponds to the radius of curvature of plane wave fronts. On the other hand, in z = ∞, R(z) = z, which corresponds to the radius of curvature of spherical wave fronts that, considering the paraxial approximation, is associated to plane waves in far-field. Finally, in z = z0, the minimum radius of curvature, or equivalently, the greater angular confinement, is achieved, and this radius corresponds to the spot size, i.e., R(z) = Rmin = 2z0. Figure 18 shows the output of the computer-aided profiler for wave fronts considering different values for the frequency, confirming the previous analysis. It can be appreciated that the greater the frequency is, the lower the curvature radius of the wave fronts and, so, the greater the tilt of the wave fronts.
Moreover, CATEM00 allows the analysis of the dependence of the TEM00 wave fronts and the Rayleigh range z0. Figure 19 depicts the output for the wave fronts for two different Rayleigh distances, where it can be inferred that the higher the Rayleigh distance, the greater the achieved angular confinement. That is, the normal vectors to the wave fronts form smaller angles with the optical axis-z along z.

5. Conclusions

The proposed computer-aided TEM00 beam tool CATEM00 provides a complete and freely-available solution for the visual design in 3D of multiple lasers and single-mode fiber output beam. Based on the fundamentals of Gaussian beam propagation in wave theory, the design is possible not only in the characterization of the associated light sources in relation to their spatial confinement, considering intensity, power, and divergence, but also in relation to its angular confinement, in the form of limited wave fronts according to the paraxial approximation. It should be noted that both the design and the knowledge of the dependence of the beams, both in space and at an angle, make it possible to propose and detect the quality of the beam with precision, which is very critical in various applications, especially in biomedical related to surgery [1]. The tool has been developed in MATLAB code and the support platform for its execution is provided free of charge, in a way that the distribution of the work can be extensive. In addition, the graphical user interface that integrates the designer makes the propagation details transparent and allows simple modeling based on the characteristic parameters of Rayleigh range, amplitude, distance on the optical axis z, radial distance to the optical axis z, and radius of the fixed section to be considered. Finally, all the necessary quantitative values can be obtained from the configurable graphic displays in azimuth, elevation, and distance. In future works, the integration of further transversal modes of light sources and lenses will be considered. Moreover, the tool will be also offered in other programming languages to make it as accessible as possible to the scientific community.

Supplementary Materials

The proposed designer software is available online in [12] as well as in the supplementary files attached to this paper (see https://www.mdpi.com/2073-8994/12/1/83/s1). Please, refer to these sources to access the tool.

Author Contributions

Conceptualization: R.P.d.P., S.G.-G., J.E.M.-E., and A.M.; software: R.P.d.P., S.G.-G., J.E.M.-E., and A.M.; validation and formal analysis: R.P.d.P., S.G.-G., J.E.M.-E., and A.M.; writing—original draft preparation: R.P.d.P., S.G.-G., and J.E.M.-E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by University of Jaén, with Innovation Project grant number PID10_201819. The APC was funded by Innovation Project grant number PID10_201819 of the University of Jaén.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alsultanny, Y.A. Laser beam analysis using image processing. J. Comput. Sci. 2006, 2, 109–113. [Google Scholar] [CrossRef] [Green Version]
  2. González-Gómez, C. Laser Beam Shaping. Master’s Thesis, Brno University of Technology, Brno, Czech Republic, 2012. [Google Scholar]
  3. Laskin, A.; Kaiser, P.; Laskin, V.; Ostrun, A. Laser beam shaping for biomedical microscopy techniques. In Biophotonics Solutions for Better Health Care, V. Proceeding of the SPIE Photonics Europe Conference; SPIE: Brussels, Belgium, 2016. [Google Scholar]
  4. Webb, C.E.; Jones, J. Handbook of Laser Technology and Applications; Institute of Physics Publishing: Bristol, UK; Philadelphia, PA, USA, 2004. [Google Scholar]
  5. Karnakis, D.M.; Fieret, J.; Rumsby, P.T.; Gower, M.C. Microhole drilling using reshaped pulsed Gaussian laser beams. In Proceedings of the SPIE—The International Society for Optical Engineering, International Symposium on Optical Science and Technology: Laser Beam Shaping II, San Diego, CA, USA, 29 July–3 August 2001. [Google Scholar]
  6. Gaussian Beam Calculator. Available online: http://www.originalcode.com/downloads/GBC8p2.zip (accessed on 12 November 2019).
  7. FreeBMP. Available online: http://www.freebpm.com (accessed on 12 November 2019).
  8. Gaussian Beam Intensity/Irradiance Calculator. Available online: http://www.calctool.org/CALC/phys/optics/gauss_power_dist (accessed on 17 December 2019).
  9. Smith, S.T. MATLAB: Advanced GUI Development; Dog Ear Publishing: Indianapolis, Indiana, 2006. [Google Scholar]
  10. Prado, R.P.; García-Cárdenas, M.; García-Galán, S.; Muñoz-Expósito, J.E. Interactive tool for learning propagation in single-mode optical fibers in telecommunication engineering. Comput. Appl. Eng. Educ. 2019, 27, 789–813. [Google Scholar] [CrossRef]
  11. Jimenez, J.; Belmonte, A.; Garrido, J.; Ruz, M.L.; Vazquez, F. Software tool for acausal physical modelling and simulation. Symmetry 2019, 11, 1199. [Google Scholar] [CrossRef] [Green Version]
  12. Computer-Aided Laser-Fiber Output Beam 3D Spatial and Angular Designer (CATEM00). Available online: https://drive.google.com/file/d/17hhZS2NGDx4EWxHvdi2wWc7NiXAtBLO-/view?usp=sharing (accessed on 20 December 2019).
  13. Ophir SP907 Beam Profiling Camera. Available online: https://www.ophiropt.com/laser--measurement/beam-profilers/products/Beam-Profiling/Camera-Based-Beam-Profiling-with-BeamGage/SP907-Beam-Profiling-Camera (accessed on 20 December 2019).
  14. Saleh, B.E.A.; Teich, M.C. Fundamentals of Photonics; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  15. Siegman, A. Lasers; University Science Books: Sausalito, CA, USA, 1986. [Google Scholar]
  16. Sun, H. Thin lens equation for a real laser beam with weak lens aperture truncation. Opt. Eng. 1998, 37, 2906–2913. [Google Scholar] [CrossRef]
  17. Freiberg, R.J.; Halsted, A.S. Properties of low order transverse modes in Argon ion lasers. Appl. Opt. 1969, 2, 355–362. [Google Scholar] [CrossRef] [PubMed]
  18. Born, M.; Wolf, E. Principles of Optics, 7th ed.; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  19. Allen, L.; Beijersbergen, M.W.; Spreeuw, R.J.; Woerdman, J.P. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Phys. Rev. 1992, 45, 8185–8189. [Google Scholar] [CrossRef] [PubMed]
  20. Yan, Y. High-capacity millimeter-wave communications with orbital angular momentum multiplexing. Nature Commun. 2014, 5, 4876. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Yan, Y. Terabit-scale orbital angular momentum mode division multiplexing in fibers. Science 2013, 340, 1545–1548. [Google Scholar]
  22. Krenn, M.; Handsteiner, J.; Fink, M.; Fickler, R.; Ursin, R.; Malik, M.; Zeilinger, A. Twisted light transmission over 143 km. Proc. Natl. Acad. Sci. USA 2016, 113, 13648–13653. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. University of Erlangen. Photon Statistics—Advanced Laboratory Course in Physics; University of Erlangen: Bavaria, Germany, 2017. [Google Scholar]
Figure 1. Examples of detected TEM00 intensity with the beam camera in Spiron SP907 at the Telecommunication Engineering Department, University of Jaén.
Figure 1. Examples of detected TEM00 intensity with the beam camera in Spiron SP907 at the Telecommunication Engineering Department, University of Jaén.
Symmetry 12 00083 g001
Figure 2. Symmetry in spatial and angular dimensions of TEM00 mode and main parameters.
Figure 2. Symmetry in spatial and angular dimensions of TEM00 mode and main parameters.
Symmetry 12 00083 g002
Figure 3. Design areas of the computer-aided profiler.
Figure 3. Design areas of the computer-aided profiler.
Symmetry 12 00083 g003
Figure 4. Computer-aided design tool flow chart.
Figure 4. Computer-aided design tool flow chart.
Symmetry 12 00083 g004
Figure 5. Graphical user interface of CATEM00.
Figure 5. Graphical user interface of CATEM00.
Symmetry 12 00083 g005
Figure 6. Parameter and analysis panels with all available design options.
Figure 6. Parameter and analysis panels with all available design options.
Symmetry 12 00083 g006
Figure 7. Input data panel.
Figure 7. Input data panel.
Symmetry 12 00083 g007
Figure 8. Intensity vs. position in z design. (a) z = 1 m; (b) z = 10 m.
Figure 8. Intensity vs. position in z design. (a) z = 1 m; (b) z = 10 m.
Symmetry 12 00083 g008
Figure 9. Intensity vs. Rayleigh Distance z0 design. (a) z0 = 0.15 m; (b) z0 = 15 m; (c) z0 = 80 m.
Figure 9. Intensity vs. Rayleigh Distance z0 design. (a) z0 = 0.15 m; (b) z0 = 15 m; (c) z0 = 80 m.
Symmetry 12 00083 g009aSymmetry 12 00083 g009b
Figure 10. Intensity dependence with frequency design. (a) f = 569.947 THz; (b) f = 193.414 THz.
Figure 10. Intensity dependence with frequency design. (a) f = 569.947 THz; (b) f = 193.414 THz.
Symmetry 12 00083 g010aSymmetry 12 00083 g010b
Figure 11. Exact/far-field intensity design.
Figure 11. Exact/far-field intensity design.
Symmetry 12 00083 g011
Figure 12. Power vs. position in optical axis z design. (a) z = 1 m; (b) z = 5 m.
Figure 12. Power vs. position in optical axis z design. (a) z = 1 m; (b) z = 5 m.
Symmetry 12 00083 g012
Figure 13. Power vs. Rayleigh Distance z0 design. (a) z0 = 30 m; (b) z0 = 2 m; (c) z0 = 0.1 m.
Figure 13. Power vs. Rayleigh Distance z0 design. (a) z0 = 30 m; (b) z0 = 2 m; (c) z0 = 0.1 m.
Symmetry 12 00083 g013aSymmetry 12 00083 g013b
Figure 14. Power dependence with frequency design. (a) f = 100.414 THz; (b) f = 900.414 THz.
Figure 14. Power dependence with frequency design. (a) f = 100.414 THz; (b) f = 900.414 THz.
Symmetry 12 00083 g014aSymmetry 12 00083 g014b
Figure 15. Fixed section power vs. width design.
Figure 15. Fixed section power vs. width design.
Symmetry 12 00083 g015
Figure 16. Exact divergence cone design. (a) f = 193.414 THz; (b) f = 4093.414 THz.
Figure 16. Exact divergence cone design. (a) f = 193.414 THz; (b) f = 4093.414 THz.
Symmetry 12 00083 g016
Figure 17. Far-field divergence cone design. (a) f = 193.414 THz; (b) f = 4093.414 THz.
Figure 17. Far-field divergence cone design. (a) f = 193.414 THz; (b) f = 4093.414 THz.
Symmetry 12 00083 g017
Figure 18. Wave fronts vs. frequency design. (a) f = 193.414 THz; (b) f = 493.414 THz.
Figure 18. Wave fronts vs. frequency design. (a) f = 193.414 THz; (b) f = 493.414 THz.
Symmetry 12 00083 g018aSymmetry 12 00083 g018b
Figure 19. Wave fronts vs. Rayleigh z0 distance design. (a) z = 0.9 m; (b) z0 = 2 m.
Figure 19. Wave fronts vs. Rayleigh z0 distance design. (a) z = 0.9 m; (b) z0 = 2 m.
Symmetry 12 00083 g019aSymmetry 12 00083 g019b

Share and Cite

MDPI and ACS Style

Prado, R.P.d.; García-Galán, S.; Muñoz-Expósito, J.E.; Marchewka, A. Computer-Aided Laser-Fiber Output Beam 3D Spatial and Angular Design. Symmetry 2020, 12, 83. https://doi.org/10.3390/sym12010083

AMA Style

Prado RPd, García-Galán S, Muñoz-Expósito JE, Marchewka A. Computer-Aided Laser-Fiber Output Beam 3D Spatial and Angular Design. Symmetry. 2020; 12(1):83. https://doi.org/10.3390/sym12010083

Chicago/Turabian Style

Prado, Rocío Pérez de, Sebastián García-Galán, José Enrique Muñoz-Expósito, and Adam Marchewka. 2020. "Computer-Aided Laser-Fiber Output Beam 3D Spatial and Angular Design" Symmetry 12, no. 1: 83. https://doi.org/10.3390/sym12010083

APA Style

Prado, R. P. d., García-Galán, S., Muñoz-Expósito, J. E., & Marchewka, A. (2020). Computer-Aided Laser-Fiber Output Beam 3D Spatial and Angular Design. Symmetry, 12(1), 83. https://doi.org/10.3390/sym12010083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop