Multiplane light conversion (MPLC) offers an alternative to adaptive optics for coupling a turbulence-corrupted free-space optical beam into single-mode fiber or waveguide. Recent published test results suggest that such a conversion device has comparable or better performance than an adaptive optics system. To better understand the device characteristics, simulations were conducted to quantify the power loss for different turbulence strengths and the number of Hermite–Gaussian modes used in the conversion process. The specific case studied is a prototype free-space laser communication system developed by the U.S. Army Research Laboratory. Simulation and statistical results of the conversion performance are reported. Analysis of a beam power combiner after MPLC is also discussed. |
1.IntroductionMultiplane light conversion (MPLC) devices perform unitary spatial transformation1–3 of optical beams composed of a set of orthogonal modes, such as Hermite–Gaussian (HG), to a different set of orthogonal modes via specially designed phase plates. MPLC covers a range of applications, such as spatial mode division (de)multiplexing for high-capacity optical fiber communications,4 light-wave coherent detection,5 and optical phase-front precompensation for communications and directed energy.6 For free-space optical communications (FSOCs), MPLC decomposes turbulence-corrupted optical beam into a relatively small number (e.g., ) of low-order HG modes that are subsequently transformed into spatially separated fundamental Gaussian modes for coupling into single-mode fibers. Guided Gaussian beams in these fibers are then coherently combined to recover the transmitted optical signal.7 An example of a simplified FSOC receiver schematic with a four-output MPLC followed by a 4-to-1 power combiner is shown in Fig. 1. Beam pointing and tracking devices, such as fast-steering mirrors, for maintaining optical alignment to the MPLC are not shown, but their presence before the MPLC as a part of the FSOC receiver architecture is assumed. Combining the outputs into one single-mode fiber or waveguide eliminates the complexity of implementing multiple or photoreceivers, especially for large . Also, the likelihood of some of these photoreceivers becoming thermal-noise-limited is expected to be higher than the detection of the combined optical signal using a single photoreceiver. Note that the four-output MPLC and the 4-to-1 power combiner shown in Fig. 1 are for clarity and illustrative purpose only. In general, -output MPLC and -to-1 power combiner configurations with , 15, 21, or higher are typical. In this report, the cases of , 15, and 21 were investigated. Current commercially available MPLC devices are based on microoptics.8 Recently published FSOC test results suggest that the MPLC device has a comparable or better performance than adaptive optics.7 Adaptive optics systems mitigate turbulence-induced beam distortion and allow free-space coupling into single-mode fiber. These systems have high performance but typically require complex subsystems, such as a wave-front sensor and two-dimensional (2D) deformable reflector array. As a result of their complexity, relatively large size, and cost, they are not commonly employed in military FSOC applications that require mobility and low size, weight, and power (SWaP). MPLC, which is passive with no moving components and has a comparatively small footprint, is an alternative to adaptive optics that potentially fulfill military requirements for FSOC. In MPLC, the turbulence-corrupted optical beam can be expressed as a superposition or summation of many orthogonal HG spatial modes with corresponding complex coefficients , as described in Sec. 7.1 in the Appendix. The transformation of HG into Gaussian modes in MPLC was achieved via free-space propagation to a series of phase plates. These specially designed phase plates impart phase-front perturbations to the input beam. As the beam propagates in the free-space region between phase plates these phase-front perturbations are converted to spatial amplitude perturbations via an equivalent 2D Fourier transform, thereby shaping the beam intensity profile. Alternating between successive phase-front perturbation induced by the phase plate and free-space propagation eventually transforms the HG mode into the desired mode, such as Gaussian beam. The phase plates profiles, number of propagations, and distance between propagations require careful designs to achieve the desired output beam profiles and distribution,1,5 e.g., spatially separated Gaussian beams. The MPLC described in Ref. 1, for example, consists of a single extended reflective phase plate and a spherical mirror in a multipass cavity configuration that performs simultaneous spatial unitary transformation of three HG modes into three spatially separated Gaussian beams that couple into single-mode fibers array, or vice versa. Conversion between a larger number of spatially separated Gaussian and HG beams or arbitrary input beam shapes has been demonstrated.4,5,7 Typically, small intensity features or high-spatial frequencies in the turbulence-corrupted optical beam require more higher-order HG modes to reproduce or recover the energy of fine features of the original turbulence-corrupted optical beam. In practice, the MPLC device is limited to a relatively small due to implementation complexity, optics and fiber coupling losses, and cost. With only a limited number of HG modes for decomposition, smaller intensity features in the original beam are not well reproduced, yielding a net power loss, i.e., , where and are the optical beam powers of the original and the -mode HG reconstructed beams, respectively. As a result, conversion loss from free-space to optical fiber due to an insufficient number of modes can be significant. Due to the random nature of turbulence, such a conversion or power loss is expected to fluctuate in a similar way as the effect of scintillation. In addition, there are static losses from optics (e.g., Fresnel loss) in the MPLC device, but they are not affected by turbulence. To better understand the impact of (the refractive-index structure parameter) and (the number of HG modes used for transformation of the original beam) on power loss or the fidelity of the reconstructed beam, wave-optics simulations were performed to generate many different turbulence-corrupted beam profiles. The simulated turbulence-corrupted beam is then decomposed into HG modes, and the corresponding complex HG mode coefficients are computed, as described in Sec. 7.1 in the Appendix. Statistics of the power loss extracted from the values versus different and were estimated with many turbulence propagation simulation trials. The case studied here followed parameters of a prototype FSOC system designed and constructed at the U.S. Army Research Laboratory (ARL). Simulation results suggest that, for , the power loss or penalty below 1 dB is achievable with only 10 HG modes decomposition. The overall variation of the power loss or normalized variance is below 0.01 over 1200 simulation trials. A larger number of HG modes or lower yields a better power loss and stability as expected. The relatively low penalty and tolerance to high is in part due to the small receiver aperture diameter () relative to the diffraction-only beam diameter () used in ARL’s FSOC prototype system, i.e., a reduced ratio, where is the Fried parameter.9 At the MPLC outputs, the transformed Gaussian beams are coupled into single-mode fibers, which are then coherently combined into one single-mode fiber or waveguide, as depicted in Fig. 1, for . A silicon-on-insulator photonic integrated circuit (PIC) beam combiner, with its small footprint, is a candidate for achieving this function. For example, a compact low-loss -branch or -junction-based10,11 -to-1 tree combiner structure (: integer) with integrated thermo-optic phase shifters used in some optical phased array designs12 with an edge couplers array can be readily employed to accept single-mode fibers coupling. The efficiency of such a symmetric tree combiner is largest when all input powers are the same and all input optical phases are aligned, i.e., identical or modulo of . However, as shown later, the mode power or HG mode coefficients distribution can be highly nonuniform due to turbulence-induced beam distortion, yielding uneven distribution of the Gaussian beam powers. As a result, degraded beam power combining efficiency is expected. Analysis and simulations were conducted to quantify such a degradation or penalty versus and for 8-, 16-, and 32-to-1 -branch tree combiners. 2.HG Modes DecompositionIn this section, the definition of the power loss or ratio that quantifies the performance of the low-mode-number HG decomposition in MPLC is introduced. This is followed by a discussion on the spatial scale factor for HG mode required to maximize the power ratio or beam reconstruction fidelity. 2.1.Power Ratio DefinitionA power ratio is defined here as the ratio of the optical power of the reconstructed () to the original () beam and is given as where and are the reconstructed and original optical field spatial profiles after the receiver aperture, respectively. The reconstructed beam is a superposition of HG modes with complex coefficients, , computed from the original turbulence-corrupted beam at the receiver, . It can be shown that this power ratio definition is equivalent to the conventional mode overlap integral, as described in Sec. 7.2 in the Appendix. The power ratio and its turbulence statistics investigated here provide one of several performance metrics for free-space coupling into single-mode fiber via MPLC.2.2.Spatial Scale Factor of HG ModesHG modes are characterized by the mode number or index and a spatial scale factor or (HG mathematical form and plots are shown in Sec. 7.3 in the Appendix). Because the spatial extend of HG mode is determined by , careful selection of its value is critical for accurately reconstructing the original beam and achieving a maximum power ratio. The selection of is not entirely obvious, especially at strong turbulences. Intuitively, should not be too different from the original beam size. However, as the product ( is the free-space propagation distance) increases, the Fried parameter or atmospheric coherence diameter seen by the original propagated beam reduces, leading to more likelihood of breakup of the optical beam into many smaller beam spots and patches with random sizes, i.e., the beam energy shifting to higher spatial frequencies with more finer intensity features. The beam size then becomes an increasingly ill-defined quantity. As a result, more smaller scale higher-order HG modes are required to reconstruct these smaller features, leading to a decrease of the “optimal” scale factor . One approach to select an optimal value is by searching for a that maximizes , i.e., the reconstructed beam power or the power ratio. This value, however, is expected to vary due to the random nature of turbulence even though and remain constant. Many propagation simulations are needed to extract the statistics of , and the ensemble average value of or is used to inform the design of MPLC devices and experimental testing. 3.Simulation Method and ResultsWave-optics simulations13 were performed to generate a turbulence-corrupted 1550-nm optical beam profile at a receiver 2 km away. This was then used to compute complex HG mode coefficients , as described in Sec. 7.1 in the Appendix. These HG modes were combined to reconstruct the original beam from which the power ratio was computed. To extract relevant statistics, 1200 independent simulation runs were executed for values of , , and and for , 15, and 21. Using the Rytov variance14 (), the three considered here are typically categorized as weak fluctuations (), strong fluctuations (), and focusing/saturation regime (), respectively. The wave-optics simulation was based on split-step propagation with the angular spectrum form of the Fresnel diffraction integral. A source was specified at the input plane and propagated across alternating instances of free space and a turbulent phase screen. The angular spectrum form of the Fresnel diffraction integral via 2D fast Fourier transform was used for the propagation calculations. Propagation was bounded by a super-Gaussian absorbing boundary to suppress energy spreading beyond the simulation grid similar to the concept of data windowing and antialiasing filter.13 Five independent random phase screens separated by 400 m were placed along the 2-km propagation path. Accurate phase screen representation of the turbulent profile is critical for obtaining reliable scintillation data and statistics. These five random phase screens were generated using the 2D Fourier transform technique.13 These phase screens were enhanced with subharmonic components to improve accuracy and low spatial frequency representation.13 The modified von Kármán power spectral density (PSD) was used for the phase screen computation, taking into account both inner and outer scales and its ability to collapse to a Kolmogorov PSD. The inner and outer scales used in the simulations were 5 mm and 5 m, respectively. A Cartesian grid with mesh points was used to represent the optical beam profiles with a 1-mm mesh size and an interpolated mesh size after beam truncation through the 46-mm receiver aperture. The quadratic phase-front envelop of the turbulence-corrupted beam due to diffraction was removed by subtracting the quadratic phase-front of a zero-turbulence propagated beam from it. Further details on wave-optics simulation can be found in Ref. 13 and the references therein. The launched optical beam at a wavelength of has a Gaussian profile that was used in all simulations. The Gaussian beam has a diameter 22.5 mm with a full beam divergence of 0.6 mrad. The diffraction-only beam diameter at 2 km is due to beam spreading without turbulence. The receiver aperture diameter is 46 mm. As a result, only a small fraction of beam area () is collected by the receiver. These parameters are derived from a prototype FSOC terminal designed and constructed at ARL. Figure 2(a) shows a series of snapshots of simulated turbulence-corrupted beam (field amplitude) within a viewing window at 2 km for to . Fried parameters computed for spherical wave () along with values are listed on the plots: , 1.96, and 7.82 for , , and , respectively. is the receiver aperture diameter, and is the propagation distance; these two quantities are fixed. The beam after the 46-mm aperture is analyzed and hereafter is referred to as the original turbulence-corrupted beam to be used for HG modes decomposition and reconstruction. Optical powers () after the 46-mm receiver aperture collected for 1200 simulation trials are shown in Fig. 2(b) (left to right), for , , and , respectively. The mean () and normalized variance (, : standard deviation) of are shown in the legend. Note that increases by about a factor of 10 for from to but only by a factor of 2.4 from to , indicating saturation. Figure 3 shows a series of snapshots (leftmost column) of the simulated turbulent-corrupted beam (optical field amplitude) truncated by the 46-mm diameter receiver aperture at 2 km for values of , , and within a viewing window slightly larger than the receiver aperture. The corresponding HG-based reconstructed beam profiles for , 15, and 21 are shown in the second to fourth columns (left to right), respectively. As can be seen, the fidelity of the reconstructed beam decreases at higher but generally improves with larger . The values shown in the reconstructed beam plots are in the number of grid points, and they tend to decrease as increases. The original beam power defined earlier is referenced to the truncated beam power or total power over the 46-mm diameter receiver aperture. Note that, for shown in the top-left plot, the beam profile is relatively flat over the receiver aperture. Decomposition of a flat profile with the first few HG modes is challenging, as can be seen in the plot next to the top-left one in which a hole or dip near the center appears to “balance” the two adjacent peaks to achieve an effective flat profile. Figure 4 shows distribution of HG mode coefficient amplitudes, , corresponds to reconstructed beams for and values depicted in Fig. 3. As expected, for the fundamental HG or Gaussian mode is largest at low , but it spreads to other with higher turbulence strengths. The general trend of versus the mode number appears exponential-like, with more complex substructures at higher values. Note that HG mode coefficients generally decrease with increasing values. This is due to the reduction of optical power transmitted through the receiver aperture (46-mm diameter) as the beam distortion increases with the turbulence strength. The power ratio or conversion loss for 1200 independent simulation trials is shown in Fig. 5 for different and values. The power ratio shown was obtained using ensemble averaged generated from separate sets of 1200 simulation runs, where was determined via power ratio maximization in each run (see Sec. 7.3 in the Appendix). It is clear that the power ratio and its fluctuation improve with larger and lower . The mean () and normalized variance () values of the power ratio over 1200 samples are also shown in the plots. Figure 6 shows of versus the HG mode number. Note that for (HG mode 1) increases with , indicating that the contribution of becomes less reliable at high turbulence conditions, as expected. Figure 7 shows the convergence of for the power ratio versus the accumulated number of simulation runs up to 1200. As can be seen, the statistics of power ratio generally stabilized around 600 or more trials. This suggests that our 1200-run statistical results should be reliable for the current investigated case. Statistics of the power ratio, mean (), and normalized variance () for different and over 1200 samples are shown in Fig. 8, including the case of . Numbers highlighted in red are ensemble averaged power ratios in decibel units, which can be used as a guide in the link budget computation. The last column shows the ensemble average of the ratio , where is the receiver aperture radius. Bar graphs for in dB units and versus are shown on the right panel for a visual comparison. For the HG mode decomposition penalty to be , MPLC with a mode number as low as 10 can tolerate a up to about . The power ratio fluctuation or is below 0.01 in this case, indicating stability and robustness of the conversion to turbulence with , as can be seen in the table and bar graph. 4.Power CombinerThe Gaussian outputs of MPLC devices are coupled into single-mode fibers, which are then coherently combined into one single-mode fiber or waveguide. A silicon-on-insulator PIC -branch10,11 tree combiner with inputs and phase shifters for input phase alignment is considered and analyzed. Then, simulation results are discussed. Other PIC power splitting/combining structures, such as directional couplers, can also be considered. The -branch was selected here due its compactness, low excess loss, relatively broad and flat spectral response, and fabrication tolerance based on our experience of fabricated PIC -branches gained from a series of silicon photonics multiproject wafer runs over the years. The performance of the power combiner or its efficiency for different and values was computed using the complex HG coefficients values or obtained from simulations. For clarity and to elucidate the impact of turbulence, internal mode conversion losses of the MPLC and fiber coupling and other static losses are ignored. These losses are turbulence-independent and can be readily included in the combiner efficiency once they are available or measured. With these assumptions and the normalized HG modes, the input optical fields to the power combiner are equal to the HG mode coefficients or , as shown in Sec. 7.3 in the Appendix. Figure 9 shows an example for a schematic of a PIC 4-to-1 -branch symmetric tree combiner with main components, such as edge couplers, phase shifters, and -branches,10,11 to combine outputs of MPLC. Here again, was selected for clarity and illustrative purpose only. The PIC phase shifters are used to perform phase-alignment: equalizing or aligning optical phases of the input optical signals from the MPLC output fibers to maximize and stabilize the combined output optical power. These phase shifters can be based on thermo-optics using resistive heaters or plasma dispersion in the PN (P-type, N-type) junction or PIN (P-type, Intrinsic, N-type) waveguide via doping.15 Recent reported phase shifters16,17 based on optomechanical PIC electrostatic actuators are also candidates to meet low SWaP requirements. Electrical controls for the phase shifters are not shown for clarity. A feedback control loop (e.g., Ref. 18) that maximizes the combiner output can be used to drive the phase shifters array, which is not considered here. As described earlier and in Sec. 7.4 in the Appendix, the -branch coherent combiner output is maximum when all of its inputs have the same power and same phase or are phase aligned. As shown in Fig. 4, or mode power distribution can be highly nonuniform, yielding an uneven distribution of the Gaussian beam powers to the combiner inputs. Also, the (: integer) number of inputs of the -branch combiner does not always match , yielding suboptimal efficiency. The impact of these impairments, along with derivation of combiner efficiency, which is defined as the ratio of the combined output power () to the original receiver aperture truncated beam power (), is described in Sec. 7.4 in the Appendix. Figure 10 shows the simulated 8-to-1 -branch tree combiner efficiency without () and with () phase alignment computed for 1200 simulation trials for different (top to bottom rows) and (left to right columns). The legend shows the normalized standard deviation () and mean () of and over 1200 samples. The numbers in red are mean efficiencies in decibels. Efficiencies for the 16- and 32-to-1 combiners were also computed. Figure 11 shows stacked bar graphs for statistics of efficiencies, mean , and normalized variance of (left column) and (right column) for 8-, 16-, and 32-to-1 combiners for the investigated and values. Each stacked bar graph is partitioned in three parts from left to right for the 8-, 16-, and 32-to-1 -branch tree combiners, respectively. As can be seen, an 8-to-1 combiner yields best efficiency for 10 HG modes, whereas a 16-to-1 combiner should be considered for 15 and 21 HG modes. Overall mean combiner efficiencies vary between and without phase alignment and between to with phase alignment for the three and three values. The -branch combiner efficiency is almost always suboptimal, as shown in Sec. 7.4 in the Appendix, due to nonequal input powers from the MPLC outputs because of turbulence-induced beam distortion. Furthermore, scalability of the combiner is likely limited due to accumulated PIC component losses at high input counts. 5.DiscussionsThe investigations conducted here mainly focused on impact of turbulence on MPLC and power combiner performance. For clarify, turbulence-independent static losses due to MPLC optics, fiber coupling, and PIC component are not included. Other imperfections of MPLC, such as internal mode conversion loss and mode crosstalk in MPLC, are not considered to elucidate the turbulence-induced penalty. The combiner efficiency described earlier is defined and normalized to the original beam power. Thus, the power ratio of the reconstructed beam is already included in the combiner efficiency. As a result, the combiner efficiency can be interpreted as a coupling efficiency from the free-space optical beam after the receiver aperture to a single-mode fiber or waveguide. The power combiner analysis presented here assumed an ideal -branch with no excess loss and perfect 50:50 power split. In practice, both excess loss and splitting ratio deviation from ideal are nonzero and wavelength dependent. For the -branch combiner, the accumulated excess loss increases geometrically as due to the cascading of -branches. The feedback control of phase shifters to maximize the combiner output will require further studies, which is beyond the scope of this investigation. An efficient optimization algorithm, such as stochastic parallel gradient descent,18 can be considered. The fiber length difference of the connecting fiber array between MPLC and PIC generates signal distortion at the combiner output similar to intersymbol interference (ISI) in multimode fibers. For a fiber length accuracy within 5 mm or a maximum delay of about 25 ps, the maximum symbol rate would be 20 Gbaud for a maximum of a half-symbol-period delay. Digital signal processing, such as adaptive equalization,19 can be applied to compensate for fiber differential delay-induced ISI. The overall receiver performance, such as signal-to-noise ratio (SNR), using a -branch combiner were estimated and are described in Sec. 7.5 in the Appendix. To mitigate the impact of polarization variation in single-mode fibers, polarization-maintaining fibers or polarization diversity using a PIC polarization rotator-splitter20 can be employed. The output of the combiner shown in Fig. 9 can be connected directly to a PIC photodetector integrated on the same chip,21 e.g., leveraging mature passive and active PIC fabrication capabilities of photonics foundries, such as American Institute for Manufacturing Integrated Photonics.22 6.SummaryAnalysis and a series of simulations were conducted to quantify the performance of MPLC and HG decomposition/reconstruction of turbulence-corrupted free-space optical beams for different and values. As a case study, simulations were conducted for a specific launched beam size, divergence, and propagation distance for a prototype FSOC system that was designed and constructed at ARL. These results provide an initial assessment on potential performance of the MPLC device with the prototype system. For example, depending on the requirements of overall loss, a design parameter space of and that is subject to constraints of the power loss based on simulation results can be formulated. The method and approach described here can be readily applied to other FSOC laser transceivers. Performance analyses were also conducted for a combiner that serializes multiple parallel outputs of MPLC into a single optical signal. Efficiencies of a “traditional” -branch symmetric tree combiners (8-, 16-, and 32-to-1) were computed for different and values. The overall robust performance of the passive MPLC and PIC combiner makes it an attractive alternative to adaptive optics system for turbulence-tolerant FSOC applications. 7.Appendix: HG Mode Decomposition and Power Combiner Analysis7.1.Power Ratio DefinitionThe turbulence-corrupted beam optical field profile at the receiver, , is expressed in terms of an infinite sum of orthogonal HG modes, , as follows: where and are the optical field amplitude and phase profiles, respectively. Using the orthogonality of the HG mode complex coefficients are computed as follows:The reconstructed beam or using only HG modes is is satisfied as . Then, Taking the integral of the absolute square of the integrand on both sides yields The above inequality is due to triangular inequality for any two nonzero complex numbers and : . Recognizing that the double-integral above is equivalent to the total optical beam power, then where is the original beam power, is the reconstructed beam power, and is a residual power loss term. Because is a nonnegative number,The above analyzed equation shows that optical power of the reconstructed beam using a finite number of HG modes is always less than the optical power of the original beam. The power ratio quantity is defined as the ratio of the two optical powers as follows: This quantity can be readily computed via wave-optics simulation for various turbulence scenarios and values. It can be shown that the power ratio quantity above is identical to the conventional definition of the mode overlap integral commonly used to quantify coupling efficiency between two different optical modes, as described next. 7.2.Mode Overlap IntegralThe mode overlap integral is defined as follows:23 where and are optical powers of the original and reconstructed beams, respectively. The integral in the numerator using HG mode decomposition is given byUsing the orthogonality properties of , One can write The above is the same as because where is normalized as described below, yielding . This results in the following yielding a power ratio the same as .7.3.HG ModesThe HG mode24 is given as (ignoring the constant Guoy phase shift) where and are the spatial positions and is associated with the spatial extend of the mode. is the Hermite polynomial and is computed using the following recursive relations:24 where . The amplitude is given asThe expression above yields normalized HG mode such that where is the Kronecker delta. For fundamental or Gaussian mode, is simply the radius of the intensity profile because . Figure 12 shows for the first 36 modes. Note that the numbering of the mode in relation to the index shown in the figures is adopted here.The is same as the spatial scale factor described in the main text. Optimal value of or was obtained by maximizing the power ratio. Figure 13 shows the ratio , where is the receiver aperture radius obtained for each 1200 simulation run for different and . As can be seen, tends to reduce with higher and values. The legend shows the ratio of the ensemble averaged to or . With the normalized HG modes defined above, the MPLC outputs to the combiner are derived as follows. The HG modes are transformed to a series of spatially separated Gaussian modes at the MPLC output: where and is an -dependent spatially shifted Gaussian mode and is a transform (free-space propagation through a series of phase plates) that converts to . The complex coefficient is associated with the transform from to via ; it may also include any losses from the optics. Because is a static loss mostly from the MPLC implementation loss and is independent of turbulence, it was assumed here that all and is normalized, i.e., . Furthermore, assuming no coupling loss from the MPLC outputs to the fiber arrays or fiber arrays to the PIC chip, the ’th input signal optical field () to the power combiner is simply given by the HG modes coefficients: .7.4.Power Combiner AnalysisThe performance of the power combiner is analyzed, and a combiner efficiency is defined here for the -branch tree combiner. Consider a single ideal lossless -branch with two input optical fields, and , as shown in Fig. 14. The combined output optical field is given as15 The combined output power is given by where and are the input optical powers. Consider the case in which the input phases are aligned such that , i.e., constructive interference in the fundamental mode at the output waveguide. The combiner efficiency is defined as follows:Note that, if , the combined output and . Figure 14 shows a plot of versus the ratio . The efficiency reaches maximum at 1, where the output power is simply the sum of the two input powers only when both input powers are equal. When one of the input power is zero, the combined output power is half of the nonzero input power. The above analysis of a single -branch combiner can be extended to -to-1 combiner. For MPLC outputs to be combined using a -to-1 -branch combiner, the combined output optical field is given as where is the complex optical field to the ’th input of the -branch combiner. The expression means selecting whichever number is smaller in the parenthesis. For the phase-aligned case (i.e., all are identical) using the compensating phase shifters, the combined output field is simply the sum of the input field amplitudes given asDue to the factor in the -branch combiner, unity combiner efficiency can only be achieved when all input amplitudes are identical and all phases are aligned. 7.5.SNR AnalysisAnalysis of SNR for the -branch power combiner with direct detection is described. The SNR of a -branch combiner with direct-detection using a single photodetector is given as where is the detector responsivity in A/W, is the detector dark current in A, is the combined optical power at the -branch combiner output in watts, is the noise bandwidth in Hz, is the detector thermal noise PSD in , and is the electronic charge. The first, second, and third term in the denominator represents shot, dark current, and thermal noise power of the detector, respectively.is related to or the optical power collected at the receiver aperture via the efficiency of the MPLC and -branch combiner described in the main text. For the ideal case and no turbulence, , the SNR is given as Figure 15 shows the computed SNR versus with -branch for thermal noise current and . The other parameters are at 1550 nm,21 ,21 , and . ReferencesG. Labroille et al.,
“Efficient and mode selective spatial mode multiplexer based on multi-plane light conversion,”
Opt. Express, 22
(13), 15599 https://doi.org/10.1364/OE.22.015599 OPEXFF 1094-4087
(2014).
Google Scholar
J. Morizur et al.,
“Programmable unitary spatial mode manipulation,”
J. Opt. Soc. Am. A, 27
(11), 2524 https://doi.org/10.1364/JOSAA.27.002524 JOAOD6 0740-3232
(2010).
Google Scholar
D. Miller,
“All linear optical devices are mode converters,”
Opt. Express, 20
(21), 23985 https://doi.org/10.1364/OE.20.023985 OPEXFF 1094-4087
(2012).
Google Scholar
G. Rademacher et al.,
“Peta-bit-per-second optical communications system using a standard cladding diameter 15-mode fiber,”
Nat. Commun., 12 4238 https://doi.org/10.1038/s41467-021-24409-w NCAOBW 2041-1723
(2021).
Google Scholar
Y. Zhang et al.,
“An ultra-broadband polarization-insensitive optical hybrid using multiplane light conversion,”
J. Lightwave Technol., 38
(22), 6286 https://doi.org/10.1109/JLT.2020.3012108 JLTEDG 0733-8724
(2020).
Google Scholar
T. Weyrauch and M. Vorontsov,
“Atmospheric compensation with a speckle beacon in strong scintillation conditions: directed energy and laser communication applications,”
Appl. Opt., 44
(30), 6388 https://doi.org/10.1364/AO.44.006388 APOPAI 0003-6935
(2005).
Google Scholar
R. Calvo et al.,
“Alternative passive fiber coupling system based on multiplane light conversion for satellite-to-ground communications,”
Proc. SPIE, 11272 112720Q https://doi.org/10.1117/12.2546076 PSISDG 0277-786X
(2020).
Google Scholar
“Home - World leader for light shaping - Cailabs,”
https://www.cailabs.com/en/
().
Google Scholar
J. Schmidt, Numerical Simulation of Optical Wave Propagation, SPIE Press, Bellingham, WA
(2010). Google Scholar
Y. Zhang et al.,
“A compact and low loss Y-junction for submicron silicon waveguide,”
Opt. Express, 21
(1), 1310 https://doi.org/10.1364/OE.21.001310 OPEXFF 1094-4087
(2013).
Google Scholar
C. Lalau-Keraly et al.,
“Adjoint shape optimization applied to electromagnetic design,”
Opt. Express, 21
(18), 21693 https://doi.org/10.1364/OE.21.021693 OPEXFF 1094-4087
(2013).
Google Scholar
R. Fatemi, A. Khachaturian and A. Hajimiri,
“A non-uniform sparse 2-D large-FOV optical phased array with a low-power PWM drive,”
IEEE J. Solid-State Circuits, 54
(5), 1200 https://doi.org/10.1109/JSSC.2019.2896767 IJSCBC 0018-9200
(2019).
Google Scholar
J. Liu et al.,
“Mid and long-wave infrared free-space optical communication,”
Proc. SPIE, 11133 1113302 https://doi.org/10.1117/12.2530969 PSISDG 0277-786X
(2019).
Google Scholar
L. Andrews and R. Phillips, Laser Beam Propagation through Random Media, SPIE Press, Bellingham, WA
(1998). Google Scholar
L. Chrostowski and M. Hochberg, Silicon Photonics Design, Cambridge University Press, Cambridge, United Kingdom
(2015). Google Scholar
M. Pruessner et al.,
“Foundry-processed optomechanical photonic integrated circuits,”
OSA Contin., 4
(4), 1215 https://doi.org/10.1364/OSAC.419410
(2021).
Google Scholar
P. Edginer et al.,
“Silicon photonic microelectromechanical phase shifters for scalable programmable photonics,”
Opt. Lett., 46
(22), 5671 https://doi.org/10.1364/OL.436288 OPLEDP 0146-9592
(2021).
Google Scholar
Q. Hu et al.,
“Adaptive stochastic parallel gradient descent approach for efficient fiber coupling,”
Opt. Express, 28
(9), 13141 https://doi.org/10.1364/OE.390762 OPEXFF 1094-4087
(2020).
Google Scholar
T. Mori et al.,
“Modal dispersion compensation by using digital coherent receiver with adaptive equalization in multi-mode fiber transmission,”
Opt. Fiber Technol., 19 132 https://doi.org/10.1016/j.yofte.2012.12.001 1068-5200
(2013).
Google Scholar
W. Sacher et al.,
“Polarization rotator-splitters in standard active silicon photonics platforms,”
Opt. Express, 22
(4), 3777 https://doi.org/10.1364/OE.22.003777 OPEXFF 1094-4087
(2014).
Google Scholar
M. Byrd et al.,
“Mode-evolution-based coupler for high saturation power Ge-on-Si photodetectors,”
Opt. Lett., 42
(4), 851 https://doi.org/10.1364/OL.42.000851 OPLEDP 0146-9592
(2017).
Google Scholar
N. Fahrenkopf et al.,
“The AIM photonics MPW: a highly accessible cutting edge technology for rapid prototyping of photonic integrated circuits,”
IEEE J. Sel. Top. Quantum Electron., 25
(5), 8201406 https://doi.org/10.1109/JSTQE.2019.2935698 IJSQEN 1077-260X
(2019).
Google Scholar
R. Hunsperger, Integrated Optics: Theory and Technology, 2nd ed.Spring-Verlag, New York
(1984). Google Scholar
A. Siegman, Lasers, University Science Books, Sausalito, CA
(1986). Google Scholar
BiographyPaul M. Pellegrino received his BS degree and PhD in physics in 1989 and 1996, respectively. He has been employed with ARL as a physicist and manager for more than 25 years. Currently, he is the chief of the Integrated Photonics Branch within Army Research Directorate in the ARL. He has authored more than 90 technical papers in optics. He is an advocate for DoD’s continued support of integrated photonics manufacturing. He is a member of OPTICA and a SPIE fellow. Justin Bickford received his BS degree in electrical engineering from Rochester Institute of Technology in Rochester, New York, United States, in 2002. He received his MS degree and PhD in electrical engineering (applied physics focus) from the University of California, San Diego, California, United States, in 2004 and 2008, respectively. Since 2008, he has been an optoelectronics research engineer for the Combat Capabilities Development Command-Army Research Laboratory (CCDC-ARL) in Adelphi, Maryland, United States, to create advanced capabilities for our soldiers. Currently, he leads an applied research team in the branch of integrated photonics at ARL and advises PMs on DARPA projects and the DoD management of the AIM Photonics Institute. |