Next Article in Journal
Antioxidant Dimethyl Fumarate Temporarily but Not Chronically Improves Intracortical Microelectrode Performance
Previous Article in Journal
A Planar-Type Micro-Biopsy Tool for a Capsule-Type Endoscope Using a One-Step Nickel Electroplating Process
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design of Single-Mode Single-Polarization Large-Mode-Area Multicore Fibers

1
Department of Electrical and Computer Engineering, Tarbiat Modares University (TMU), Tehran 1411713116, Iran
2
Department of Engineering Physics, Polytechnique Montréal, Montréal, QC H3C 3A7, Canada
3
School of Electrical and Electronic Engineering, University of Adelaide, Adelaide 5001, Australia
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(10), 1901; https://doi.org/10.3390/mi14101901
Submission received: 1 September 2023 / Revised: 28 September 2023 / Accepted: 29 September 2023 / Published: 4 October 2023

Abstract

:
In laser science and industry, considerable effort is directed toward designing fibers for fiber laser and fiber amplifier applications, each of which offers a particular advantage over the others. Evanescently coupled multicore fibers, however, have been studied less extensively due to the relatively small mode area in the single-mode regime. Here, by proposing a new structure with stress-applying parts in a 37-core fiber and optimizing this structure through a comprehensive framework, we present 21 solutions characterized by large-mode-area and high beam quality in the single-mode, single-polarization regime. Different fiber designs are optimal for different output parameters. In one design, the mode area can significantly increase to above 880 μ m 2 , which is comparable with that of photonic-crystal fibers. Moreover, besides the single-mode operation, the beam quality factor ( M 2 factor) of the fundamental mode is considered an output parameter in the bent state and is improved up to 1.05 in another design. A comprehensive tolerance analysis is then performed to assess the performance of the designs under deviations from normal conditions. Moreover, in spite of the shifts in the loss of modes, the proposed high beam quality LMA fibers maintain single-polarization, single-mode operation across a wide range of core pitches, bending orientation angles, and bending radius deviations. Our results highlight the potential of multicore fibers for the efficient operation of fiber lasers and amplifiers.

1. Introduction

High-power fiber lasers have been used as alternatives for semiconductor lasers and other related technologies due to their high output power, beam quality, efficient heat dissipation, and compactness [1,2,3,4,5,6]. These types of fibers are also used in sensors. One of the prominent advantages of fiber sensors based on photonic-crystal is their easy manufacturing compared to other sensors [7,8,9,10]. The most important characteristic of fibers designed for fiber lasers and fiber amplifiers is LMA (Table 1 indicates the acronyms.), which is essential for increasing the delivered power, reducing nonlinear effects, and ensuring single-mode operation to enhance beam quality and mode stability. In addition, to achieve compact and low-footprint devices, the fiber is coiled, necessitating that the bending loss of the FM be minimized as much as possible [1,2,11,12,13]. Furthermore, to maintain a simple laser setup and to avoid a reduction in the degree of polarization, the use of single-polarization fibers, or polarization-preserving fibers, is necessary [14,15,16,17,18].
Increasing the EMA in conventional fibers by increasing the core diameter can severely degrade the beam quality and mode stability and lead to the excitation of HOMs. Moreover, bending the fiber highly increases FM loss in conventional fibers. Regarding these issues, there is an inherent trade-off between the conflicting FOMs of the fiber [19]. Several different structures, such as LMA fibers with asymmetric bend compensation, segmented cladding fibers, PCFs, LCFs, and MCFs, have been proposed to reduce the trade-off between the mode area, single-mode operation, and FM loss of bent fibers [20,21,22,23,24,25,26,27,28,29,30]. Additionally, in polarization-preserving fibers, the birefringence induced by stress using SAPs not only removes one of the polarizations of all modes but also decreases the FM-guided polarization leakage [15]. Another difficulty that arises as a result of both bending the fiber and applying birefringence is the deformation of the Gaussian profile of the FM. Although single-mode operation—achieved by increasing HOM loss and decreasing FM loss—results in high output beam quality and mode stability, the deformation of the FM as a result of fiber bending makes beam focusing more challenging [1,2,16,17,31,32]. Among all the proposed structures, multicore fibers present another suitable option, offering diffraction-limited output beam quality, a less complicated cladding structure compared to PCFs and LCFs, and a higher numerical aperture than step-index fibers [20,33]. The air holes in the cladding of a PCF must be precisely spaced and aligned. Even small imperfections can lead to additional losses [34]. However, the EMA of all-solid multicore fibers is low in the SM-operating regime [20,33,35]. It would be rational to think that by combining the advantages of multicore fibers with single-polarization fibers, one can improve the characteristics of the fiber. However, the conflict between output parameters increases due to the complicated nature of the fiber and the interplay between different physical parameters. To strike a balance between EMA, single-mode and single-polarization operation, and high beam quality, artificial intelligence techniques can be utilized [36,37,38,39,40,41], among which, the framework proposed in [42] is more compatible with our work. Furthermore, all solid dual-core fibers and dual-core photonic-crystal fibers have been designed and optimized to enhance the large-mode-area characteristic in single-mode operation [43].
In this paper, we design single-mode, single-polarization, large-mode-area multicore fibers (SM-SP-LMA-MCFs) that yield diffraction-limited output beam quality. In addition to utilizing 37 coupled cores to ensure a large mode area, SAPs are embedded on both sides of the core region to induce birefringence, thereby eliminating one of the polarizations. For the first time, to the best of our knowledge, apart from improving SM operation to improve beam quality, the M 2 factor of the FM in the bent state is introduced as an output parameter to the simulations to prevent the deformation of the FM in the coiled fiber. Then, considering all figures of merit and using a comprehensive multi-objective optimization framework, the fiber is optimized, and the results are analyzed with regard to the effects of FOMs. This paper is organized as follows. Section 2 describes the proposed design, the method of analysis, and the optimization framework. Section 3 presents the results, and the sensitivity of the optimal structures to variations in different structural parameters is meticulously investigated. Finally, the paper is summarized in Section 4.

2. Structure and Theory

2.1. Designof the Single-Mode Single-Polarization Large-Mode-Area Multicore Fiber

Figure 1 shows the two-dimensional schematic of the proposed SM-SP-LMA-MF with diffraction-limited output beam quality. To remove one polarization of the fundamental mode in a high beam quality fiber, boron-doped silica, as SAPs, have been used on both sides in a multicore fiber [15]. Instead of using a 19-core fiber proposed in [20], 37 cores are arranged in a hexagonal lattice to ensure a larger mode area. As illustrated in Figure 1, parameters r 1 to r 20 represent the radii of cores, Λ indicates the core pitch, d L and d R denote the distance of the left and right SAPs, and parameters h L 1 , h L 2 , l L 1 , L L 2 , h R 1 , h R 2 , l R 1 , and l R 2 represent the dimensions of the left and right SAPs. Since the exerted bending is considered to be only along the x-direction, the structure is designed to be symmetrical along the y-direction. The background material is silica with a refractive index of 1.4496, and the NA of each core is 1.08 at Λ = 1.064 μ m [20,42]. The boron-doped silica used in SAPs has a refractive index [15] of 1.4416. Since the core region can be co-doped with materials, such as fluorine and boron oxide (to have a refractive index equal to the refractive index of pure silica), the results are valid for both active and passive fibers [13]. More information about the thermal expansion coefficient, Young’s modulus, Poisson’s ratio, and photoelastic constant can be found in reference [44].

2.2. Theoretical Description

The mode calculation is implemented in two steps. In the first step, the refractive index change induced by SAPs is calculated, whereas, during the second step, the calculated refractive index tensor is used to solve the Helmholtz equation via FEM in association with the PML boundary condition, thereby calculating the eigenvalues and the mode profiles.
In the first step, the new refractive index profile induced by thermal expansion is evaluated. Using the following equilibrium equation (Equation (1)), the stresses induced by thermal expansion can be calculated using [44]
. σ = . ε x ε y γ x y α α 0 ( 1 + υ ) ( T T ref ) ,
σ = D ε ,
where T r e f denotes the reference temperature, T is the annealing temperature, υ is Poisson’s ratio, α represents the expansion coefficient, ε x ( y ) and γ x y are the strain components, σ is the stress tensor, and D is the elasticity matrix. The relationship between the applied stress and the change in refractive index ( Δ n) can be calculated according to the elasto-optical effect equation [44]:
Δ n = C σ .
Using Equation (2) and with the values of C 1 and C 2 given in references [15,44], the birefringence (B) is calculated using the following equation:
B = | ( C 2 C 1 ) ( σ x σ y ) | .
In the second step, using the refractive index tensor calculated in the previous step, we solve the vector Helmholtz equation [45]:
Δ × ( Δ × E ) k 0 2 n 2 E = 0 ,
where E represents the electric field vector, k 0 is the wavenumber in the free space, and n is the modified refractive index. The output parameters that are required to evaluate the performance of the designed fiber are the effective mode area, bending loss, and M 2 factor. The bending loss of each mode (dB/m) is calculated using the imaginary part of its effective refractive index [27]:
L = 8.686 k 0 I m ( n e f f ) ,
where n e f f is the effective index of the mode, and k 0 is the free space wavenumber (1/m). The EMA of the propagating mode can be calculated using [27]:
A e f f = | E | 2 d x d y 2 | E | 4 d x d y .
The M 2 factor of a propagating mode is calculated according to M x 2 and M y 2 using the following equations proposed in [46],
M 2 = M x 2 M y 2 ,
where
M j 2 = 4 π σ 0 j 2 σ j 2 For j = x or y .
Let j = x or y, and terms σ 0 j and σ j are the second-order moments in the spatial and frequency domains, respectively. Term σ 0 j can be calculated using [46]
σ 0 j 2 = ( j j ¯ ) 2 | E ( x , y ) | 2 d x d y | E ( x , y ) | 2 d x d y ,
where the first-order moment in the spatial domain of a given mode is given by [46]
j ¯ = j | E ( x , y ) | 2 d x d y | E ( x , y ) | 2 d x d y .
Based on the 2D-Fourier transform theory, the field distribution in the frequency domain E( f x , f y ) is obtained by [46]
E ( f x , f y ) = e x p [ i 2 π ( f x x + f y y ) ] d x d y ,
where f x and f y are the transverse spatial frequencies. According to the frequency domain distribution of the field, the second-order moments in the frequency domain can be given by [46]
σ 2 = ( f j f ¯ j ) | E ( f x , f y ) | 2 d f x d f y | E ( x , y ) | 2 d f x d f y ,
where the first-order moment in the frequency domain of a given mode is given by [46]
f ¯ j = f j | E ( f x , f y ) | 2 d f x d f y | E ( x , y ) | 2 d f x d f y .

2.3. Optimization Framework

Using MOGWO in combination with an image processing technique—to classify all modes according to their order—we can optimize the proposed structure and offer designs that can provide worthy characteristics [42]. In addition to calculating the area and loss of modes using the FEM simulator, the effect of stress induced by thermal expansion on the refractive index profile and the M 2 factor have been included in the framework.
All the structural parameters, including the radii of the cores, the core pitch, SAP dimensions, and SAP distances from the core, are defined in a PM to be optimized.
P M = [ r 1 / Λ , , r 20 / Λ , h L 1 / Λ , h L 2 / Λ , l L 1 / Λ , l L 2 / Λ , h R 1 / Λ , h R 2 / Λ , l R 1 / Λ , l R 2 / Λ , d L / Λ , d R / Λ , Λ ]
Moreover, all fabrication and bending limitations, along with the desired range for the objectives, are incorporated in a CM, see Table 2. The beam quality factor in the x and y directions can be improved by decreasing the M 2 -factor and the term | M x 2 M y 2 | .
Figure 2 shows the calculation process of the output parameters. After calculating the refractive index change induced by SAPs for the candidate design, the FEM simulator calculates the mode profiles and their effective refractive indices. Mode order recognition is performed using image-processing techniques. Using the information from the image-processing part and FEM simulator, the EMA, L F X , L F Y / L F X , L H , M 2 , and | M x 2 M y 2 | are calculated. As the optimizer attempts to increase all the objectives, the objectives that should be decreased are multiplied by 1 , such as L F X , M 2 , and | M x 2 M y 2 | . If the merit factors are not in the acceptable range defined by the CM, a very low value ( 1 × 1 20 ) will be assigned to them. Hence, the optimizer can find that such candidate designs are not feasible and, consequently, they are neglected during the optimization process.

3. Results and Discussion

The parameters of 21 optimal designs, characterized by a large mode area, low bending loss for the x-polarized FM ( L F X ), high bending loss for the y-polarized FM ( L F Y ), and HOMs ( L H ), along with high beam quality (indicated by a low M 2 factor and | M x 2 M y 2 | ) are reported in Table A1, as detailed in Appendix A. Simulations are performed for a bending radius of 20 cm and a wavelength of 1064 nm. Given the many conflicting output parameters, the optimization problem has several solutions. These solutions are a set of solutions that trade-off these output parameters. All solutions have a standard output range for all parameters, and each design has at least one advantage over all other designs, and some designs may have multiple advantages. For example, design no. 9 has a higher EMA when compared to nos. 2, 3, 4, 5, 6, 7, 8, 10, 11, 12, 14, 16, 17, 18, 19, 20, and 21, as well as higher LH in comparison with nos. 1, 2, 5, 13, 14, 15, 16, and 21. While keeping other objectives in the desirable range by putting constraints on them, the EMA, L F Y / L F X , and L H can be increased separately to the range of 880 ( μ m 2 ), 4.95, and 485 dB/m for design nos. 15 and 17.
Likewise, without degrading the values of other FOMs to the undesirable regime, LFY, M 2 , and | M x 2 M y 2 | are decreased to the range of 0.12 dB/m, 1.05, and 0 for design n o s . 19 and 5. Figure 3 illustrates the distributed birefringence induced by SAPs at the fiber cross-section. The refractive index in the y-direction ( n y ) is so decreased that the effective refractive index of the y-polarized fundamental mode is in the same range as the cladding index. As a result, index-matched coupling occurs between y-polarized FM and cladding, which increases the loss of the y-polarized fundamental mode. As can be seen in Table 2, the height of SAPs is increased to the maximum range, around 2 Λ , to increase birefringence. Nevertheless, a further increase in the height of SAPs reduces the loss of all modes, including L F Y and L H . In fact, since the refractive index of the boron-doped silica is below the refractive index of cores and fused silica, they act as a barrier to the modes. Similarly, d L and d R —which imply the distance of SAPs from the center of the fiber—are in the range of 3.2 Λ –3.82 Λ , which means that they are placed near the core region to increase birefringence. Nonetheless, decreasing the distance gives rise to a reduction in the core area and, hence, there is a decrease in EMA. The average value of the left, middle, and right core radii are 0.1573 Λ , 0.1490 Λ , and 0.1462 Λ , respectively. Indeed, to compensate for the effect of bending, the refractive index in the bending direction is decreased, which is in agreement with the results reported in the literature [13,21,26].
The designs representing the best structures with the highest EMAs, L F Y / L F X , L H , and L F Y , as well as the lowest, M 2 , L F X , and | M x 2 M y 2 | , are highlighted in bold in Table A1 in Appendix A. Design no. 19—with the highest beam quality—has an M 2 factor of 1.05 in the bent state. Designing a fiber with an M 2 factor lower than 1.1 in the 20 cm bending radius is a problem. Despite this, a multicore fiber with a proper core arrangement and SAP dimensions can achieve a low deformation of the x-polarized fundamental mode. As can be seen in Figure 4, design no. 19 has a low bending loss and exhibits a Gaussian beam profile for the x-polarization of the fundamental mode. According to Table A1 in Appendix A, design no. 15—with the highest EMA—has the lowest acceptable L H . This structure has the highest core pitch (5.8 μ m) among all designs, so the fundamental mode occupies a larger area. Given the fact that there is a trade-off between the EMA and single-mode operation, this design has the lowest acceptable L H . Figure 5 provides clear evidence about L F X , L F Y , and L H in design nos. 15 and 17. Regarding the values of L F X , L F Y , and L H , we can see that the leakier the mode, the higher its loss. Due to the low leakage of the HOM power in the cladding, the L H of design no. 15 is the lowest among all structures. Therefore, this structure has the least effective single-mode operation. On the other hand, as can be seen in Figure 5f, the optical field strongly penetrates the cladding region and, therefore, the lower confinement level of the HOM in design no. 17 is a testament to its higher L H and effective single-mode operation.
Designing high-quality practical fibers requires considering all fabrication variations. Additionally, the versatility of the design enables it to operate in different settings, which is of great interest to end-users. In the following, the sensitivity of the fiber to the core pitch, the bending orientation, and the bending radius is investigated. Taking into account fiber output values in addition to sensitivities, users can select their desired design among the proposed designs.
Because unwanted changes in the bending orientation are inevitable during the bending process, evaluating the fiber characteristics under different bending orientation angles is of great importance [31]. Having a large Gaussian-like mode profile of the x-polarized fundamental mode confined to the core region is crucial for sustaining efficiency across various bending orientation angles, in addition to preserving single-mode and single-polarized operations. Figure 6 illustrates the effect of the bending orientation angle ( θ ) on the FOM performances of 21 optimal designs. We note that the EMAs for all designs are kept at the same level when the bending orientation angles change. However, both L F X and L F Y increase with the increase in the bending orientation angle. Because the fiber is designed to work optimally at θ = 0 , the optimal refractive index profile of the fiber’s cross-section is for this condition. Thus, changes in the value of θ result in deviations in the refractive index profile. The larger the changes in the bending orientation angle, the more pronounced the modifications are in the refractive index of the fiber cross-section, leading to more leakage of both polarizations of the fundamental mode. In addition, since the SAPs—both of which have lower refractive indices than background material—are positioned in the bending direction, they act as strong barriers to the modes’ expansion into the cladding region.
Nevertheless, by rotating the fiber along its axis, the barrier-like behaviors of SAPs in some areas of the fiber cross-section decrease and, hence, increase mode losses. In the same way, the L H increases as the bending orientation angle increases. However, in the case of design no. 1, it is not an upward curve, rather, it reaches a peak at θ = 9 . While x-polarized L P 11 has the lowest loss among higher-order modes for θ = 0 in all designs, the loss of x-polarized L P 11 exceeds that of y-polarized L P 11 , even for θ > 9 in design no. 1. Therefore, since we have introduced the L H as the minimum higher-order mode loss, L H corresponds to the loss of the y-polarized LP11, even for θ > 9 . Figure 7 provides valuable information about the mode profile and bending loss of the x- and y-polarized LP11, even in the bending orientation angles of 0 , 9 , and 20 for design no. 1. As the bending orientation angle increases, n y increases, thus decreasing the loss of the y-polarization of L P 11 , which is associated with L H for θ > 9 . Lastly, for all designs, because the mode profile of x-polarized FM slightly changes by changing the fiber orientation from a vertical position, the M 2 factor increases. For most of the designs, this parameter is in a standard range. Considering the sensitivities of all these parameters to bending orientation angle, we note that the large-mode-area fiber maintains both single-mode and single-polarized operation under varying bending orientations, at the cost of an increase in L F X .
Figure 8 shows the dependence of fiber characteristics on the bending radius. As can be seen from Figure 8a, EMA is minimally affected by the bending radius changes for the majority of designs. Among all designs, design nos. 15 and 21 with high core pitches are more sensitive to the bending radius. Despite this, the EMAs of these designs remain large even in small bending radii. The profiles of the x-polarized FMs in both designs with bending radius values of 10, 20, and 30 cm are shown in Figure 9.
The confinement increases with the bending radius in both designs, yet the localization increases in design no. 21 and decreases in design no. 15. This contrast leads to a reduction of EMA in the former and an improvement in the latter. Moreover, the smaller the fiber bending radius, the higher the L F X , L F Y , and L H for all designs. Thus, in the bending radius range of 10 to 30 cm, all designs remain in the single-mode single-polarization operation, but in a different loss range. For small bending radii, L F X is weakly guided but its difference with the y-polarization of the fundamental mode loss and HOMs is so high that the single-mode single-polarization operation is assured. On the other hand, as can be seen in Figure 8b,c, for large bending radii, where L F X and L F Y are low, the loss of the x-polarized fundamental mode is low enough to obtain the single-mode single-polarization operation.
Figure 10 indicates the dependence of output parameters on the core pitch ( Λ ). In this case, since the value of all structural parameters is defined as a ratio to the pitch, increasing the core pitch leads to an increase in the SM SP LMA MCF scale as a whole. As a result, the larger the core pitch, the larger the area the modes occupy, and the greater the EMA (Figure 10a). However, due to the weak coupling between the modes in two adjacent cores and larger core sizes in the larger pitch, the bending loss for all modes decreases as well, as seen in Figure 10b–d. This result is in agreement with the step-index fibers, where an increase in the core size results in the decreasing mode loss and appearance of HOMs. In addition to changing the refractive index, doping the fiber with rare-earth dopants can also change the spectroscopic properties of the core medium, such as the absorption and fluorescence spectra [47]. Furthermore, the quantum defect caused by the rare-earth dopants in the core medium results in significant heating of the core during high-power operations [48].

4. Conclusions

In summary, we proposed a novel large-mode-area fiber, implementing stress-applying parts in the cladding region of a multicore fiber. This proposed structure was optimized using a multi-objective optimization technique, considering a bending radius of 20 cm and a wavelength of 1.064 μ m. Accordingly, 21 optimal designs were obtained, each demonstrating single-mode, single-polarization operation, large mode area, and high beam quality. Different designs offer trade-offs between various output parameters. We demonstrated that one of the optimal designs represents the highest recorded effective mode area of 880 μ m 2 . An x-polarized fundamental mode loss of 0.12 dB/m, a y-polarized FM loss of 2.05, and a higher-order mode loss of 485 dB/m are among the best-recorded merit factors of our optimal designs. Apart from single-mode operation, a beam quality factor of 1.05 guarantees the fiber’s beam quality in the bent state. Additionally, the sensitivity of all designs to the core pitch size, bending radius, and bending orientation angle were studied in detail. We find that the large effective mode area of the fiber can be guaranteed, even under different bending radii and bending orientation angles. Moreover, despite the changes in the loss of the x-polarized fundamental mode, it was observed that the fiber maintains single-mode and single-polarization operations under different core pitch sizes, bending orientation angles, and bending radii.

Author Contributions

K.R. did the main simulations and wrote the draft and prepared the figures. D.F. and H.T. defined the project and prepared the final manuscript. J.M. prepared the Latex version. S.M.M. did the optimization. D.F. supervised the project and the presented results. D.A. edited the final version. All authors have read and agreed to the published version of the manuscript.

Funding

The work reported in this work has received no external funding.

Data Availability Statement

The datasets used and/or analyzed during the current study are available upon reasonable request from the corresponding author.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A

Table A1. Structural parameters of the 21 SM-SP-LMA-MCF design.
Table A1. Structural parameters of the 21 SM-SP-LMA-MCF design.
#1#2#3#4#5#6#7#8#9#10#11#12#13#14#15#16#17#18#19#20#21
r 1 / Λ 0.130.130.110.180.160.150.140.150.150.150.130.110.150.110.140.170.10.130.150.130.1
r 2 / Λ 0.150.180.160.160.150.170.180.120.150.140.140.140.180.170.190.180.170.150.180.160.19
r 3 / Λ 0.160.170.160.210.180.180.170.180.190.170.190.170.190.180.150.190.170.180.210.140.14
r 4 / Λ 0.150.140.110.150.140.160.140.150.120.140.150.120.160.140.120.120.120.130.140.130.14
r 5 / Λ 0.150.150.10.190.140.180.130.150.140.140.140.120.150.160.170.130.140.110.160.120.11
r 6 / Λ 0.140.150.130.140.130.130.160.160.170.140.150.130.130.120.140.160.150.140.160.160.15
r 7 / Λ 0.170.150.140.150.130.130.160.140.140.150.160.170.140.180.10.120.130.190.160.150.16
r 8 / Λ 0.140.130.130.180.160.170.120.170.130.140.110.130.120.120.160.110.150.120.130.140.11
r 9 / Λ 0.140.120.120.170.120.140.120.110.140.140.130.160.110.160.110.150.130.140.160.130.11
r 10 / Λ 0.170.170.120.150.180.130.160.170.160.160.130.130.170.120.150.180.130.170.150.170.17
r 11 / Λ 0.180.180.20.20.180.20.170.210.170.170.180.150.190.140.180.150.180.170.180.150.15
r 12 / Λ 0.130.150.160.180.160.170.140.150.160.150.160.170.130.160.140.140.170.190.220.180.12
r 13 / Λ 0.120.150.130.190.130.180.140.150.140.140.150.130.170.170.140.160.130.130.160.150.12
r 14 / Λ 0.180.150.130.140.160.150.120.160.140.150.130.130.120.120.120.160.160.120.140.150.16
r 15 / Λ 0.160.180.180.190.170.170.190.160.180.190.190.180.180.180.140.170.190.20.170.180.2
r 16 / Λ 0.10.160.170.140.160.150.160.150.120.140.130.150.140.130.130.110.160.140.170.130.13
r 17 / Λ 0.170.110.120.140.140.160.120.160.150.130.110.130.130.110.150.140.120.130.140.120.14
r 18 / Λ 0.150.210.180.170.180.140.190.170.170.160.170.160.20.20.220.170.180.220.180.170.19
r 19 / Λ 0.170.160.180.180.160.170.170.130.140.140.160.140.180.190.140.120.170.160.150.170.17
r 20 / Λ 0.130.150.120.150.180.150.160.150.150.160.150.140.140.110.160.170.120.140.180.140.13
h l 1 / Λ 1.581.741.711.781.631.891.611.751.781.61.641.731.671.571.821.631.881.911.711.661.54
h l 2 / Λ 1.71.581.481.751.481.771.691.671.581.731.461.641.631.851.721.541.481.851.741.481.52
l l 1 / Λ 0.940.790.810.80.770.790.780.870.830.90.860.850.850.950.810.830.880.960.930.780.92
l l 2 / Λ 0.870.710.780.70.770.670.790.810.80.780.860.830.720.860.690.690.730.860.830.740.88
h r 1 / Λ 1.981.821.691.891.841.791.731.861.941.961.791.721.731.981.671.921.6321.571.81.84
h r 2 / Λ 1.741.811.651.981.651.831.781.811.721.841.641.881.691.811.761.771.731.881.851.61.8
l r 1 / Λ 0.880.840.950.810.90.980.950.870.860.870.790.90.940.810.890.910.960.870.850.880.84
l r 2 / Λ 0.90.9110.840.890.950.930.990.810.930.850.970.880.810.90.970.90.950.820.86
d l / Λ −3.7−3.4−3.5−3.7−3.4−3.7−3.4−3.4−3.5−3.4−3.2−3.5−3.4−3.4−3.6−3.3−3.6−3.5−3.5−3.5−3.6
d r / Λ 3.73.53.53.43.73.83.63.83.43.63.73.73.73.73.63.63.63.93.73.53.5
Λ ( μ m)5.15.24.845.14.25.14.65.15.15.155.15.25.85.24.54.54.15.15.7
E M A ( μ m 2 )779644502456642498617635705660650601710604880702493494428605656
L F X (dB/m)0.520.130.540.130.270.240.290.50.450.230.410.550.270.140.230.420.430.260.120.440.18
L F Y / L F X 3.723.852.384.292.392.053.373.523.393.913.283.743.163.714.953.873.354.294.473.244.68
L H (dB/m)59573641906724111911579851132415957163348539623318152
L F Y (dB/m)1.920.491.280.540.640.50.981.771.540.921.342.050.850.511.121.641.441.140.521.430.85
M 2 1.241.161.141.061.161.061.171.161.231.151.191.151.201.151.341.271.091.101.051.181.32
| M x 2 M y 2 |   ×   10 3 11026514633171116201814187334142115866

References

  1. Chen, M.Y.; Li, Y.R.; Zhou, J.; Zhang, Y.K. Design of asymmetric large-mode area optical fiber with low-bending loss. J. Light. Technol. 2013, 31, 476–481. [Google Scholar] [CrossRef]
  2. Raja, G.T.; Varshney, S.K. Extremely large mode-area bent hybrid leakage channel fibers for lasing applications. IEEE J. Sel. Top. Quantum Electron. 2014, 20, 251–259. [Google Scholar] [CrossRef]
  3. Stutzki, F.; Jansen, F.; Otto, H.J.; Jauregui, C.; Limpert, J.; Tünnermann, A. Designing advanced very-large-mode-area fibers for power scaling of fiber-laser systems. Optica 2014, 1, 233–242. [Google Scholar] [CrossRef]
  4. Liu, Y.; Zhang, F.; Zhao, N.; Lin, X.; Liao, L.; Wang, Y.; Peng, J.; Li, H.; Yang, L.; Dai, N.; et al. Single transverse mode laser in a center-sunken and cladding-trenched Yb-doped fiber. Opt. Express 2018, 26, 3421–3426. [Google Scholar] [CrossRef]
  5. Jafari, S.; Fathi, D.; Taleb, H. Optimal pulse-shaping in actively q-switched ytterbium-doped fiber lasers. IEEE Access 2020, 8, 77716–77724. [Google Scholar] [CrossRef]
  6. Jafari, S.; Fathi, D.; Taleb, H.; Abbott, D. Suppression of the Stimulated Brillouin and Raman Scattering in Actively Q-switched Fiber Lasers through Temporal Pulse Shaping. Ann. Phys. 2021, 533, 2000541. [Google Scholar] [CrossRef]
  7. Mahmoodi, Y.; Fathi, D. New oil derivative refractive index sensor using ribbon of multi-walled carbon nanotubes. Optik 2023, 276, 170644. [Google Scholar] [CrossRef]
  8. Mousavi-Kiasari, S.M.G.; Rashidi, K.; Fathi, D.; Taleb, H.; Mirjalili, S.M.; Faramarzi, V. Computational design of highly-sensitive graphene-based multilayer SPR biosensor. Photonics 2022, 9, 688. [Google Scholar] [CrossRef]
  9. Mahmoodi, Y.; Fathi, D. High-performance refractive index sensor for oil derivatives based on MWCNT photonic crystal microcavity. Opt. Laser Technol. 2021, 138, 106865. [Google Scholar] [CrossRef]
  10. Majdi, M.; Fathi, D. Graphene Based Nano Bio-sensor: Sensitivity Improvement. Sci. Iran. 2017, 24, 3531–3535. [Google Scholar] [CrossRef]
  11. Limpert, J.; Schmidt, O.; Rothhardt, J.; Röser, F.; Schreiber, T.; Tünnermann, A.; Ermeneux, S.; Yvernault, P.; Salin, F. Extended single-mode photonic crystal fiber lasers. Opt. Express 2006, 14, 2715–2720. [Google Scholar] [CrossRef] [PubMed]
  12. Tsuchida, Y.; Saitoh, K.; Koshiba, M. Design of single-moded holey fibers with large-mode-area and low bending losses: The significance of the ring-core region. Opt. Express 2007, 15, 1794–1803. [Google Scholar] [CrossRef] [PubMed]
  13. Napierała, M.; Nasiłowski, T.; Bereś-Pawlik, E.; Berghmans, F.; Wójcik, J.; Thienpont, H. Extremely large-mode-area photonic crystal fibre with low bending loss. Opt. Express 2010, 18, 15408–15418. [Google Scholar] [CrossRef]
  14. Schmidt, O.; Rothhardt, J.; Eidam, T.; Röser, F.; Limpert, J.; Tünnermann, A.; Hansen, K.; Jakobsen, C.; Broeng, J. Single-polarization ultra-large-mode-area Yb-doped photonic crystal fiber. Opt. Express 2008, 16, 3918–3923. [Google Scholar] [CrossRef]
  15. Molardi, C.; Sun, B.; Yu, X.; Cucinotta, A.; Selleri, S. Polarization-Maintaining Large Mode Area Fiber Design for 2-μm Operation. IEEE Photonics Technol. Lett. 2016, 28, 2483–2486. [Google Scholar] [CrossRef]
  16. Ustimchik, V.; Rissanen, J.; Popov, S.; Chamorovskii, Y.; Nikitov, S.A. Anisotropic tapered polarization-maintaining large mode area optical fibers. Optics Express 2017, 25, 10693–10703. [Google Scholar] [CrossRef]
  17. Zhang, W.; Lou, S.; Wang, X.; Yan, S.; Tang, Z.; Xing, Z. A broadband single mode single polarization metal wires-embedded hollow core anti-resonant fiber for polarization filter. Opt. Fiber Technol. 2019, 53, 102011. [Google Scholar] [CrossRef]
  18. Liu, M.; Yuan, H.; Shum, P.; Shao, C.; Han, H.; Chu, L. Simultaneous achievement of highly birefringent and nonlinear photonic crystal fibers with an elliptical tellurite core. Appl. Opt. 2018, 57, 6383–6387. [Google Scholar] [CrossRef]
  19. Li, M.J.; Chen, X.; Liu, A.; Gray, S.; Wang, J.; Walton, D.T.; Zenteno, L.A. Limit of effective area for single-mode operation in step-index large mode area laser fibers. J. Light. Technol. 2009, 27, 3010–3016. [Google Scholar]
  20. Vogel, M.M.; Abdou-Ahmed, M.; Voss, A.; Graf, T. Very-large-mode-area, single-mode multicore fiber: Erratum. Opt. Lett. 2010, 35, 465. [Google Scholar] [CrossRef]
  21. Fini, J.M. Large mode area fibers with asymmetric bend compensation. Opt. Express 2011, 19, 21866–21873. [Google Scholar] [CrossRef]
  22. Xing, Z.; Wang, X.; Lou, S.; Zhang, W.; Tang, Z.; Yan, S. Bend resistant large mode area fiber with step-index core and single trench. Opt. Fiber Technol. 2019, 48, 15–21. [Google Scholar] [CrossRef]
  23. Coscelli, E.; Dauliat, R.; Poli, F.; Darwich, D.; Cucinotta, A.; Selleri, S.; Schuster, K.; Benoit, A.; Jamier, R.; Roy, P.; et al. Analysis of the modal content into large-mode-area photonic crystal fibers under heat load. IEEE J. Sel. Top. Quantum Electron. 2015, 22, 323–330. [Google Scholar] [CrossRef]
  24. Li, W.; Chen, D.; Qinling, Z.; Hu, L. Large-mode-area single-mode-output Neodymium-doped silicate glass all-solid photonic crystal fiber. Sci. Rep. 2015, 5, 1–5. [Google Scholar] [CrossRef] [PubMed]
  25. Ma, S.; Ning, T.; Li, J.; Pei, L.; Zhang, C.; Wen, X. Detailed study of bending effects in large mode area segmented cladding fibers. Appl. Opt. 2016, 55, 9954–9960. [Google Scholar] [CrossRef]
  26. Saini, T.S.; Kumar, A.; Sinha, R.K. Asymmetric large-mode-area photonic crystal fiber structure with effective single-mode operation: Design and analysis. Appl. Opt. 2016, 55, 2306–2311. [Google Scholar] [CrossRef]
  27. Kabir, S.; Razzak, S.A. An enhanced effective mode area fluorine doped octagonal photonic crystal fiber with extremely low loss. Photonics Nanostructures Fundam. Appl. 2018, 30, 1–6. [Google Scholar] [CrossRef]
  28. Ma, S.; Ning, T.; Lu, S.; Zheng, J.; Li, J.; Pei, L. Bending-resistant design of a large mode area segmented cladding fiber with resonant ring. J. Light. Technol. 2018, 36, 2844–2849. [Google Scholar] [CrossRef]
  29. Ma, S.; Ning, T.; Pei, L.; Li, J.; Zheng, J. Bend-resistant large mode area fiber with novel segmented cladding. Opt. Laser Technol. 2018, 98, 113–120. [Google Scholar] [CrossRef]
  30. Ma, S.; Ning, T.; Pei, L.; Li, J.; Zheng, J. Design and analysis of a modified segmented cladding fiber with parabolic-profile core. Laser Phys. Lett. 2018, 15, 035104. [Google Scholar] [CrossRef]
  31. Wang, X.; Lou, S.; Lu, W. Bending orientation insensitive large mode area photonic crystal fiber with triangular core. IEEE Photonics J. 2013, 5, 7100408. [Google Scholar] [CrossRef]
  32. Ji, J.; Lin, H.; Sidharthan, R.; Ho, D.; Zhou, Y.; Nilsson, J.; Yoo, S. Bendable large-mode-area fiber with a non-circular core. Appl. Opt. 2018, 57, 6388–6395. [Google Scholar] [CrossRef] [PubMed]
  33. Wang, P.; Cheng, G.; Yi, R.; Liu, X.; Shang, T.; Wang, Z.; Guo, L. Theoretical and experimental study of 37-core waveguides with large mode area. Appl. Opt. 2013, 52, 7981–7986. [Google Scholar] [CrossRef] [PubMed]
  34. Guan, C.; Yuan, L.; Shi, J. Supermode analysis of multicore photonic crystal fibers. Opt. Commun. 2010, 283, 2686–2689. [Google Scholar] [CrossRef]
  35. Fini, J.M. Large-mode-area multicore fibers in the single-moded regime. Opt. Express 2011, 19, 4042–4046. [Google Scholar] [CrossRef]
  36. Gonzalez, R.C.; Woods, R.E.; Eddins, S. Digital Image Processing Using MATLAB; Pearson Prentice Hall: Upper Saddle River, NJ, USA, 2004. [Google Scholar]
  37. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey wolf optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef]
  38. Mirjalili, S.; Mirjalili, S.M.; Hatamlou, A. Multi-verse optimizer: A nature-inspired algorithm for global optimization. Neural Comput. Appl. 2016, 27, 495–513. [Google Scholar] [CrossRef]
  39. Mirjalili, S.; Saremi, S.; Mirjalili, S.M.; Coelho, L.D.S. Multi-objective grey wolf optimizer: A novel algorithm for multi-criterion optimization. Expert Syst. Appl. 2016, 47, 106–119. [Google Scholar] [CrossRef]
  40. Mirjalili, S.M.; Merikhi, B.; Mirjalili, S.Z.; Zoghi, M.; Mirjalili, S. Multi-objective versus single-objective optimization frameworks for designing photonic crystal filters. Appl. Opt. 2017, 56, 9444–9451. [Google Scholar] [CrossRef]
  41. Safdari, M.J.; Mirjalili, S.M.; Bianucci, P.; Zhang, X. Multi-objective optimization framework for designing photonic crystal sensors. Appl. Opt. 2018, 57, 1950–1957. [Google Scholar] [CrossRef]
  42. Rashidi, K.; Mirjalili, S.M.; Taleb, H.; Fathi, D. Optimal design of large mode area photonic crystal fibers using a multiobjective gray wolf optimization technique. J. Light. Technol. 2018, 36, 5626–5632. [Google Scholar] [CrossRef]
  43. Huttunen, A.; Törmä, P. Optimization of dual-core and microstructure fiber geometries for dispersion compensation and large mode area. Opt. Express 2005, 13, 627–635. [Google Scholar] [CrossRef] [PubMed]
  44. Schreiber, T.; Schultz, H.; Schmidt, O.; Röser, F.; Limpert, J.; Tünnermann, A. Stress-induced birefringence in large-mode-area micro-structured optical fibers. Opt. Express 2005, 13, 3637–3646. [Google Scholar] [CrossRef]
  45. Adams, M.J. An Introduction to Optical Waveguides; John Wiley & Sons: Hobboken, NJ, USA, 1981. [Google Scholar]
  46. Guo-Quan, Z.; Yan, F. M2 factor of four-petal Gaussian beam. Chin. Phys. B 2008, 17, 3708. [Google Scholar] [CrossRef]
  47. Qiao, Y.; Wen, L.; Wu, B.; Ren, J.; Chen, D.; Qiu, J. Preparation and spectroscopic properties of Yb-doped and Yb–Al-codoped high silica glasses. Mater. Chem. Phys. 2008, 107, 488–491. [Google Scholar] [CrossRef]
  48. Hansen, K.R.; Alkeskjold, T.T.; Broeng, J.; Lægsgaard, J. Thermo-optical effects in high-power ytterbium-doped fiber amplifiers. Opt. Express 2011, 19, 23965–23980. [Google Scholar] [CrossRef]
Figure 1. The transverse cross-section of the proposed SM-SP-LMA-MCF.
Figure 1. The transverse cross-section of the proposed SM-SP-LMA-MCF.
Micromachines 14 01901 g001
Figure 2. The modified module for calculating output parameters.
Figure 2. The modified module for calculating output parameters.
Micromachines 14 01901 g002
Figure 3. The birefringence distribution induced by SAPs in SM-SP-LMA-MCFs.
Figure 3. The birefringence distribution induced by SAPs in SM-SP-LMA-MCFs.
Micromachines 14 01901 g003
Figure 4. The electric field distribution (V/m) of the x-polarized fundamental mode for design no. 19.
Figure 4. The electric field distribution (V/m) of the x-polarized fundamental mode for design no. 19.
Micromachines 14 01901 g004
Figure 5. Mode field distribution for (a) the x-polarized FM of design no. 15, (b) the x-polarized FM of design no. 17, (c) the y-polarized FM of design no. 15, (d) the y-polarized FM of design no. 17, (e) the HOM of design no. 15, and (f) the HOM of design no. 17.
Figure 5. Mode field distribution for (a) the x-polarized FM of design no. 15, (b) the x-polarized FM of design no. 17, (c) the y-polarized FM of design no. 15, (d) the y-polarized FM of design no. 17, (e) the HOM of design no. 15, and (f) the HOM of design no. 17.
Micromachines 14 01901 g005
Figure 6. Variations of the (a) EMA, (b) L F X , (c) L F Y , (d) L H , and (e) M 2 factor on the bending orientation angle.
Figure 6. Variations of the (a) EMA, (b) L F X , (c) L F Y , (d) L H , and (e) M 2 factor on the bending orientation angle.
Micromachines 14 01901 g006
Figure 7. The electric field distribution of (a) x-polarized LP11, observed with θ = 0 (b) y-polarized L P 11 , observed with θ = 0 (c) x-polarized L P 11 , observed with θ = 9 (d) y-polarized L P 11 , observed with θ = 9 (e) x-polarized L P 11 , observed with θ = 20 (f) y-polarized L P 11 , observed with θ = 20 .
Figure 7. The electric field distribution of (a) x-polarized LP11, observed with θ = 0 (b) y-polarized L P 11 , observed with θ = 0 (c) x-polarized L P 11 , observed with θ = 9 (d) y-polarized L P 11 , observed with θ = 9 (e) x-polarized L P 11 , observed with θ = 20 (f) y-polarized L P 11 , observed with θ = 20 .
Micromachines 14 01901 g007
Figure 8. (a) EMA (b) L F X , (c) L F Y , and (d) L H behavior versus the bending radius.
Figure 8. (a) EMA (b) L F X , (c) L F Y , and (d) L H behavior versus the bending radius.
Micromachines 14 01901 g008
Figure 9. Mode field distribution for the x-polarized FM of design no. 21 with (a) R = 10 cm, (b) R = 20 cm, and (c) R = 30 cm, and design no. 15 with (d) R = 10 cm, (e) R = 20 cm, and (f) R = 30 cm.
Figure 9. Mode field distribution for the x-polarized FM of design no. 21 with (a) R = 10 cm, (b) R = 20 cm, and (c) R = 30 cm, and design no. 15 with (d) R = 10 cm, (e) R = 20 cm, and (f) R = 30 cm.
Micromachines 14 01901 g009
Figure 10. The effect of the core pitch size on the (a) EMA, (b) L F X , (c) L F Y , (d) L H , and (e) L F Y / L F X .
Figure 10. The effect of the core pitch size on the (a) EMA, (b) L F X , (c) L F Y , (d) L H , and (e) L F Y / L F X .
Micromachines 14 01901 g010
Table 1. Acronyms.
Table 1. Acronyms.
AcronymsMeaningAcronymsMeaning
SAPsStress-applying partsLMALarge-mode-area
FMFundamental modeEMAEffective Mode Area
HOMsHigher-order modesFOMsFigure of merits
PCFsPhotonic-crystal fibersLCFsLeaky channel fibers
MCFsMulticore fibersSMSingle-mode
SPSingle-polarizationNANumerical aperture
FEMFinite element methodPMLPerfectly matched layer
EMAEffective mode areaCMConstraint module
MOGWOMulti-objective gray wolf optimizerPMParameter module
Table 2. Range of input and output Parameters.
Table 2. Range of input and output Parameters.
ParametersMinMax
Λ ( μ m)410
[ r 1 , r 2 , , r 19 , r 20 ] 0.1 Λ 0.4 Λ
[ h l 1 , h l 2 , h r 1 , h r 2 ] 02 Λ
[ l l 1 , l l 2 , l r 1 , l r 2 ] 01 Λ
d 1 −5 Λ −3.2 Λ
d 2 3.2 Λ 5 Λ
E M A ( μ m 2 )400-
L F X (dB/m)-0.6
L F Y / L F X 2-
L H (dB/m)15-
L F Y (dB/m)0.5-
M 2 11.4
| M x 2 M y 2 |   ×   10 3 00.2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rashidi, K.; Fathi, D.; Maleki, J.; Taleb, H.; Mirjalili, S.M.; Abbott, D. Design of Single-Mode Single-Polarization Large-Mode-Area Multicore Fibers. Micromachines 2023, 14, 1901. https://doi.org/10.3390/mi14101901

AMA Style

Rashidi K, Fathi D, Maleki J, Taleb H, Mirjalili SM, Abbott D. Design of Single-Mode Single-Polarization Large-Mode-Area Multicore Fibers. Micromachines. 2023; 14(10):1901. https://doi.org/10.3390/mi14101901

Chicago/Turabian Style

Rashidi, Kamyar, Davood Fathi, Javad Maleki, Hussein Taleb, Seyed Mohammad Mirjalili, and Derek Abbott. 2023. "Design of Single-Mode Single-Polarization Large-Mode-Area Multicore Fibers" Micromachines 14, no. 10: 1901. https://doi.org/10.3390/mi14101901

APA Style

Rashidi, K., Fathi, D., Maleki, J., Taleb, H., Mirjalili, S. M., & Abbott, D. (2023). Design of Single-Mode Single-Polarization Large-Mode-Area Multicore Fibers. Micromachines, 14(10), 1901. https://doi.org/10.3390/mi14101901

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