Correlation and Power Distribution of Intercore Crosstalk Field Components of Polarization-Coupled Weakly Coupled Single-Mode Multicore Fibres

: The correlation and power distribution of intercore crosstalk (ICXT) ﬁeld components of weakly coupled multicore ﬁbers (WC-MCFs) are important properties that determine the statistics of the ICXT and ultimately impact the performance of WC-MCF optical communication systems. Using intensive numerical simulation of the coupled mode equations describing ICXT of a single-mode WC-MCF with intracore birefringence and linear propagation, we assess the mean, correlation, and power distribution of the four ICXT ﬁeld components of unmodulated polarization-coupled homogeneous and quasi-homogeneous WC-MCFs with a single interfering core in a wide range of birefringence conditions and power distribution among the ﬁeld components at the interfering core input. It is shown that, for homogeneous and quasi-homogeneous WC-MCFs, zero mean uncorrelated ICXT ﬁeld components with similar power levels are observed for birefringence correlation length and birefringence beat length in the ranges of [ 0.5m,10m ] and [ 0.1m,10m ] , respectively, regardless of the distribution of power between the four ﬁeld components at the interfering core input. These results suggest that, for the range of birefringence correlation length and birefringence beat length typically found in singlecore telecommunication ﬁbers, we may model the ICXT power ﬂuctuations by a χ 2 distribution with four degrees of freedom in homogeneous and quasi-homogeneous polarization-coupled single-mode WC-MCFs, regardless of the distribution of power between the four ﬁeld components at the interfering core input.


Introduction
Intercore crosstalk (ICXT) is one of the most important physical impairments of weakly coupled multicore fibers (WC-MCFs). The ICXT can limit the use of WC-MCFs in modern optical communication systems (OCSs) by imposing restrictions to the transmission distance [1], core count [2], and use of higher order modulation formats [3]. Nevertheless, WC-MCFs have already been proven as a potential solution to improving the capacity of next-generation OCSs using space division multiplexing [4,5] by increasing their capacities by one order of magnitude comparatively to OCSs using singlecore fibers [5]. In fact, very recently, experiments replicating real-use environments [6,7] and a field-deployed testbed [8] were used to prove the effective potential of WC-MCFs.
In order to devise mitigating strategies for the effects of ICXT on the performance of WC-MCF systems, several works have proposed ICXT models based on coupled-mode theory and coupled-power theory to evaluate the mean ICXT power and other stochastic parameters related to ICXT [9][10][11][12][13][14][15][16][17][18][19]. Relying on those models, the impact of the ICXT on the performance of different types of OCSs were assessed [20][21][22]. However, although those models have considered effects such as fiber bending, twist, and structure fluctuations in the accumulation of ICXT along the WC-MCF, the impact of birefringence on ICXT accumulation has been much less considered.
WC-MCFs can be classified as follows [23]: (i) polarization-coupled fibers, for which random polarization effects dominate intercore coupling; (ii) scalar-coupled fibers, for which random structural fiber fluctuations effects dominate intercore coupling; and (iii) mixed-coupled fibers, where both effects, random polarization and random structure fiber fluctuations, are important. Up to now, the majority of studies have considered scalar-coupled WC-MCFs and only a few studies have considered birefringence effects in WC-MCFs [23][24][25][26].
In fact, only very recently, it was shown that random polarization-mode coupling due to intracore birefringence can play a crucial role in explaining the ICXT accumulation in polarization-coupled and mixed-coupled WC-MCF [23] and is closely related to the polarization-mode dispersion [26].
The stochastic properties of ICXT power fluctuations have been broadly studied for scalar-coupled WC-MCFs, both theoretically and experimentally. Several studies have considered and others have shown that, for unmodulated homogeneous and quasihomogeneous WC-MCFs with a single interfering core and multiple interfering cores, under quite general conditions, the four ICXT field components at the WC-MCF output are uncorrelated gaussian components with zero mean and balanced power [11,19,27]. Hence, the ICXT power fluctuations can be modelled by a chi-square χ 2 distribution with four degrees of freedom [11,19,27].
For polarization-coupled WC-MCFs, a preliminary study of the probability density function (PDF) of the ICXT power fluctuations has been presented [23] where a χ 2 distribution with four degrees of freedom has been shown using numerical simulation for a homogeneous WC-MCF and two pairs of values of birefringence correlation length and birefringence beat length. This shows that, for polarization-coupled WC-MCFs, at least for some specific conditions, the four ICXT field components are zero mean uncorrelated gaussian-distributed. In this context, it is relevant to assess when the mean, the correlation, and the power distribution of the ICXT field components hold those properties under quite general practical conditions of birefringence in polarization-coupled homogeneous and quasi-homogeneous WC-MCFs. This allows us to conclude the validity of modelling the ICXT power fluctuations by a χ 2 distribution with four degrees of freedom in polarization-coupled WC-MCFs.
In this article, the mean, correlation, and power distribution of the ICXT field components in unmodulated polarization-coupled single-mode WC-MCFs are studied by intensive numerical simulation, considering linear propagation and stochastic intracore birefringence effects along a quite large range of birefringence conditions in homogeneous and quasi-homogeneous WC-MCFs. We show that, for both homogeneous and quasihomogeneous WC-MCFs, the four ICXT field components are zero mean and uncorrelated and have similar power levels for birefringence correlation length and birefringence beat length in the ranges of [0.5 m, 10 m] and [0.1 m, 10 m], respectively.

ICXT Model and Computation of Mean ICXT Power in Single-Mode WC-MCFs
This section presents the model used to obtain by numerical simulation the ICXT field of a single-mode WC-MCF considering intracore birefringence and linear propagation. In order to focus the analysis on the statistical properties of the ICXT field generated by each interfering core of a WC-MCF, we consider the interference between two cores only, designated within the paper, generically, by m (interfering core) and n (interfered core). Figure 1 shows a schematic illustration of the two-core WC-MCF cross section representing the interfering core m (blue) with radius a m , the interfered core n (red) with radius a n , the corresponding core pitch Λ mn , and the cladding of the MCF with diameter d c . The usual range of values for the WC-MCF size parameters are the following: (i) core radius [3.8 µm, 4.7 µm] [17,20,28,29]; (ii) core pith [30 µm, 45 µm] [8,17,20,[28][29][30]; and (iii) core cladding diameter [125 µm, 260 µm] [5] (p. 22), [8,20,25,[29][30][31][32]].  Figure 1. Schematic illustration of a two-core MCF cross section and the cartesian coordinates: core m (interfering core) with radius a m and core n (interfered core) with radius a n , the corresponding core pitch Λ mn , and the cladding of the MCF with diameter d c .

Evolution of Interfering and Interfered Fields in a Single-Mode WC-MCF
The coupled-mode equations (CME) used to model the crosstalk between cores m and n of a single-mode WC-MCF assuming intracore birefringence and linear propagation along the longitudinal coordinate of the MCF, z, with the propagation constants perturbed by fiber bending and twisting is given by [23] where β c,m = k 0 · n (int) eff,m and β c,n = k 0 · n (int) eff,n , with n (int) eff,m and n (int) eff,n being the intrinsic effective refractive index of cores m and n, respectively, i.e., the effective refractive index of the cores not perturbed by bend, twist, or any structure fluctuation along the MCF, and k 0 = 2π/λ, where λ is the wavelength. β b,m (z) and β b,n (z) are the perturbations of the propagation constants induced by bending and twist at the coordinate z of cores m and n, respectively. It should be emphasized that n (int) eff,m and n (int) eff,n are the effective refractive indexes for the fundamental (single) modes propagating through cores m and n, respectively.
Defining θ m (z) and θ n (z) as the accumulated phase due to propagation from the MCF input up to coordinate z in cores m and n, respectively, given by θ m (z) = z 0 β c,m (z )dz and θ n (z) = z 0 β c,n (z )dz , defining A n,p (z) = e θ n (z) · E n,p (z) ; A m,p (z) = e θ m (z) · E m,p (z) , and replacing (7) in (1) results in the following CMEs: In the CMEs (8), the terms f nm,p (z) and f mn,p (z) are four random processes introduced to model the phase fluctuations induced locally by the random structure fluctuations along the MCF on the contributions to the field of cores n and m, respectively. These random processes are usually characterized by their correlation lengths. We assume that the processes f nm,p (z) and f mn,p (z) are independent. The influence of this assumption on the ICXT statistics is expected to be reduced. This statement is supported by the fact that, when the birefringence effects can be neglected and the depletion of the interfering field is very low, only the random process that models the phase fluctuation induced on the interfered core affects the ICXT field [18]. Additionally, very small differences in the mean ICXT power estimates obtained by numerical simulation results were observed when independent or equal random processes are considered.
Koshiba et al. analyzed the effects of correlation length of the random process that models the phase fluctuation on the mean ICXT power when that random process includes only the random effect of structure fluctuation and does not include that of macrobend and twist, which are considered deterministic [12,13]. They investigated different types of autocorrelation function of the random process that models the phase fluctuation and found that the exponential autocorrelation function given by is adequate to estimate the mean ICXT power, where l c is the correlation length of the random process introduced to model the phase fluctuations induced by the random structural fluctuations of the MCF [12,13,33]. In the following, we consider that the random processes f nm,p (z) and f mn,p (z) follow the same statistics, with their autocorrelation functions holding R f nm,p (z) = R f mn,p (z) = R f (z). The power spectral density corresponding to the exponential autocorrelation function has a Lorentzian shape, and it is given by which is obtained by applying the Fourier transform operator defined as (6) to expression (9).

Computation of the Mean ICXT Power in a WC-MCF
We define the normalized ICXT power, XT, as the ratio between the power at the output of the interfered core n and the power at the input of the interfering core m (with no power launched in core n), i.e., where L is the MCF length. The mean of XT is denoted by XT , where · is the statistical mean operator. In the following, we refer to XT as the mean ICXT power. The field in the two directions, x and y, at the interfering core input is described as are the corresponding in-phase (I) and in-quadrature (Q) components. In the following, the index q denotes the quadrature components of the field in direction p, i.e., q ∈ {I, Q}. The power distribution at the input of the interfering core m between the two directions x and y, and the two quadrature components, I and Q, is controlled by two parameters, ζ x,y and ζ I,Q : (i) ζ x,y ∈ [0, 1] controls the power distribution between the two directions x and y, and (ii) ζ I,Q ∈ [0, 1] controls the power distribution between the I and Q components (assumed the same in both polarization directions). Therefore, the distribution of the field at the input of core m, A m (0), between the two directions x and y can be written as The ICXT field in the two directions at the interfered core output is described as A n,x (L) = A n,x,I (L) + A n,x,Q (L) and A n,y (L) = A n,y,I (L) + A n,y,Q (L). The mean ICXT power can be written as where XT p,q is the ICXT power component that describes the normalized portion of the crosstalk power at the output of the interfered core n in the q component of the p direction of the field, A 2 n,p,q (L), with the mean of the ICXT power component given by WC-MCFs can be classified in terms of the dominant crosstalk mechanism as follows [23]: (i) polarization-coupled fibers, for which Re R 2 β (z) is much narrower than R f (z), where Re{ · } denotes the real part operator; (ii) scalar-coupled fibers, for which The expressions were derived to assess the mean ICXT power at a generic propagation distance, z, for the following regimes [23]: (i) birefringence correlation length much lower than the birefringence beat length (L C L B ); (ii) birefringence correlation length comparable with the birefringence beat length (L C ≈ L B ); and (iii) birefringence correlation length much larger than the birefringence beat length (L C L B ). These are general expressions that can be particularized for polarization-coupled, scalar-coupled, and mixed-coupled fibers. In this article, we consider polarization-coupled WC-MCFs, imposing l c L B , L C .

Description of the Numerical Simulator and Assessment of the Number of Samples Required to Obtain Stabilized Estimates of the Mean ICXT Power
The objectives of this section are the following: (i) to describe the implementation of the numerical simulator and to discriminate the range of values considered for its parameters, and (ii) to assess the number of samples necessary to obtain stabilized estimates of the mean ICXT power and its four components using the simulator.

Brief Description about the Simulator Implementation
Four independent random processes β n,(1,2) (z) and β m,(1,2) (z) with exponential autocorrelation functions given by expressions (4) are generated by solving the stochastic differential Equation (3) using the numerical method proposed in Reference [39], which is a third-order Runge-Kutta method adapted to solve this type of stochastic differential equations. The statistical accuracy of the simulated random processes β n,(1,2) (z) and β m,(1,2) (z) was confirmed by comparing their autocorrelation and power spectral density functions obtained by simulation with the theoretical expressions (4) and (5), respectively.
Four independent random processes f nm,p (z) and f mn,p (z) with exponential autocorrelation functions given by expression (9) are also considered. As those random processes model the phase fluctuations induced by the random structure fluctuations [12,33], the random processes f nm,p (z) and f mn,p (z) can be written as f nm,p (z) = exp  · w m,p (z) and f mn,p (z) = exp  · w n,p (z) , where w m,p (z) and w n,p (z) are four independent Wiener processes that model the phase fluctuations induced by the random structure fluctuations. In order to obtain the exponential autocorrelation function with correlation length l c , the four Wiener processes are generated by the numerical integration of zero mean white noise with power spectral density of 2/l c [40] (pp. 71-74). The statistical accuracy of the simulated random processes f nm,p (z) and f mn,p (z) was confirmed by comparing their autocorrelation and power spectral density functions obtained by simulation with the theoretical expressions (9) and (10), respectively.
In the numerical simulation, a fifth-order Runge-Kutta method is used to solve the CMEs (8). Two highlights should be mentioned. (i) The step size ∆z used to solve numerically Equation (8) should be judiciously selected in order to simulate both the structure fluctuations and birefringence effects rigorously. On one hand, it should hold ∆z 1/|∆β c,mn |, with ∆β c,mn = β c,m − β c,n , since for correlation lengths of the random structural fluctuations of the order of 1 mm or higher as usually found in WC-MCFs [13,15,29,33] and non-homogeneous WC-MCFs, the faster variation of the right-hand sides of CMEs (1) results from the complex exponentials of the difference of phases θ m and θ n , which is imposed by ∆β c,mn . On the other hand, ∆z L C and ∆z √ 2L B /(6π) should also hold, simultaneously, to rigorously simulate the birefringence effects. (ii) In order to be used in the numerical solving of CMEs (1), samples of four independent random processes β n,(1,2) (z) and β m,(1,2) (z) and four independent Wiener processes w m,p (z) and w n,p (z) with adequate step sizes in z have to be generated.
The mean ICXT power estimated by numerical simulation was obtained considering two situations: polarization-coupled homogeneous ∆β c,mn = 0 rad · m −1 and quasihomogeneous ∆β c,mn = −701 rad · m −1 WC-MCF with length of L = 100 m. ∆β c,mn is adjusted by setting adequately the intrinsic effective refractive index of the core n, n (int) eff,n . Pairs for the birefringence correlation length and the birefringence beat length (L C , L B ) have been tested in the ranges of 0.5 m L C 10 m and 0.1 m L B 10 m, respectively. These ranges have been settled in order to include the set of values found typical for L C and L B in singlecore telecommunication fibers [37,41], [42] (p. 6), [35] (pp. 251, 261, 267, 269). To the best of authors' knowledge, measured values of the birefringence correlation length and birefringence beat length have not been reported yet for WC-MCFs.
According to the ranges considered for ∆β c,mn , L C and L B , ∆z in the range of 10 −5 − 10 −4 m was used in all simulations, averaging over 200 MCF realizations, with each realization corresponding to different samples of the eight noise random processes β n,(1,2) (z) and β m,(1,2) (z), and w n,p (z) and w m,p (z). Table 1 shows the parameters kept constant in the numerical results. We stress that we considered l c = 10 6 m and that the bending radius and twisting period were adjusted so that the bending and twisting effects on the ICXT field can be neglected. The parameters κ mn and κ nm were adjusted to κ mn = κ nm = 4.5 × 10 −5 m −1 and κ mn = κ nm = 4.5 × 10 −2 m −1 , according to the homogeneous or quasi-homogeneous regimes, respectively, in order to ensure a mean ICXT power range of about [−90 dB, −50 dB] in all situations tested.
In order to ensure acceptable computation times to estimate the statistical parameters of the ICXT field components through numerical simulation, the MCF length considered in the simulation study is L = 100 m. In fact, to generate β n,(1,2) (z) and β m,(1,2) (z) by solving the stochastic differential Equation (3), to propagate the random signal through the MCF solving the CMEs (8), and to obtain the ICXT power given by (11), a computational time of about 15 minutes per MCF realization is necessary. As 200 MCF realizations are required to obtain a confident estimate of the mean ICXT, the previous process must be repeated 200 hundreds of times for each (L B , L C ) pair, leading to a total time of about 50 hours per pair. These computational times are reached, running the simulations in a PC having a CPU Intel Core I9 10900KF with a clock speed of 5.3 GHz, a memory cache of 20 MB, and a RAM of 64 GB. These computational times increase proportionally to the MCF length. Therefore, in practical terms, acceptable computational times are hardly compatible to MCF lengths longer than 100 m.

Dependence of Estimates of Mean ICXT Power on the Number of Samples
Extended simulation tests were performed in order to assess the number of averaging samples necessary to stabilize the mean ICXT power estimator given by (13). A wide range of situations, considering 13 different pairs of birefringence correlation length and birefringence beat length, (L C , L B ), in the ranges indicated in Section 3.1, both for homogeneous ∆β c,mn = 0 rad · m −1 and quasi-homogeneous ∆β c,mn = −701 rad · m −1 single-mode WC-MCFs have been tested. For each pair (L C , L B ), different power distributions between the two directions x and y, and I and Q components (assumed the same in both polarization directions) of the interfering core (core m) have been tested by assuming different pairs ζ x,y , ζ I,Q , namely, (0, 0), (0, 0.5), (0, 1), (0.5, 0), (0.5, 0.5), (0.5, 1), (1, 0), (1, 0.5), and (1, 1). Figure 2 shows two paradigmatic examples among the universe of situations tested, illustrating the mean ICXT power estimator given by (13) in dB, i.e., 10 · log 10 XT , as a function of the number of samples, as well as each of the four components of Equation (13), XT p,q also in dB, i.e. 10 · log 10 XT p,q : (i) Figure 2a considers a homogeneous WC-MCF ∆β c,mn = 0 rad · m −1 with L C = 10 m, L B = 1 m, ζ x,y = 0.5 and ζ I,Q = 0.5, and (ii) Figure 2b considers a quasi-homogeneous WC-MCF ∆β c,mn = −701 rad · m −1 with L C = 1 m, L B = 1 m, ζ x,y = 0.5, and ζ I,Q = 1.  Figure 2 shows that the fluctuation of the mean ICXT power estimates occurring for 200 averaging samples is very low and, consequently, the stabilization achieved can be considered adequate. Figure 2 shows also a quasi-uniform distribution between the four components of the mean ICXT power at the interfered core output. Moreover, differences of about 6 dB between the mean ICXT power of each component and the total mean ICXT power are observed, which are in accordance with Equation (13). The fluctuations observed for 200 averaging samples are attributed to the inherent randomness of the estimate of the mean from a finite number of samples.
In the universe of situations tested, the level of fluctuations of the estimate of the mean ICXT power observed with an increase in the number of averaging samples is similar to the one shown in Figure 2. Therefore, all situations tested revealed that 200 averaging samples are adequate to ensure the stabilization of the estimators of mean ICXT power obtained by simulation.
To validate the numerical simulator, several situations were tested to assess the mean ICXT power as a function of L C for different L B , namely, 1 m and 10 m, both for homogeneous ∆β c,mn = 0 rad · m −1 and quasi-homogeneous ∆β c,mn = −701 rad · m −1 polarization-coupled WC-MCF. L C in the range indicated in Section 3.1 was considered. Excellent agreement between simulation estimates and analytical results obtained using the expressions proposed in [23] was achieved in all situations tested.

Parameters Used to Characterize the Statistical Properties of the ICXT Field Components
We use the following three parameters to quantify the mean deviation from 0; the level of correlation; and the power unbalancing of the ICXT field components, A n,p,q (L): (i) The normalized mean of the fields components at the WC-MCF output defined by and σ An given by with −∞ ≤ A n,p,q (L) ≤ +∞; (ii) The correlation coefficients of the ICXT field components defined as [43], pp. 153-154.
(iii) The unbalancing parameter defined in percentage as with XT p,q defined as (14), XT defined as (11), and max · defines the maximum of the four values between braces. This parameter assesses the unbalance in the mean power of the four ICXT field components at the WC-MCF output (in the two directions x and y, and in the two quadrature components I and Q). From Equations (13) and (18), we can easily infer the range for U as 0 % ≤ U ≤ 300 %: (i) when the mean ICXT power is uniformly distributed among the two directions x and y and the two quadrature components I and Q, XT p,q = XT /4, resulting in U = 0 %; (ii) on the other hand, when the mean ICXT power is entirely in one of the two directions x or y and one of the two quadrature components I or Q, the maximum of the four XT p,q terms between braces is equal to XT , resulting in U = 300 %.

Numerical Results
The objective of this section is, using the simulator developed in Section 3.1, to assess some of the most important statistical properties of the four ICXT field components A n,p,q (L), namely the means, correlations, and unbalancing power between them. Those properties are assessed for a wide range of situations, namely different (L C , L B ) pairs in the ranges indicated in Section 3.1 and different pairs of ζ x,y , ζ I,Q in the range of 0 ≤ ζ x,y ≤ 1 and 0 ≤ ζ I,Q ≤ 1, both for homogeneous ∆β c,mn = 0 rad · m −1 and quasihomogeneous ∆β c,mn = −701 rad · m −1 single-mode WC-MCF. More than 200 different situations are tested.
In the following two subsections, we present some of the most representative simulation results achieved considering L C = 10 m and L B = 10 m (L C ≈ L B ), L C = 1 m and L B = 10 m (L C L B ), and L C = 10 m and L B = 0.1 m (L C L B ) for homogeneous ∆β c,mn = 0 rad · m −1 and quasi-homogeneous ∆β c,mn = −701 rad · m −1 single-mode WC-MCF. Figure 3 shows the simulation results of the normalized mean of the A n,x,I (L) field component given by (15) (Figure 3a), the correlation coefficient between the field components A n,x,I (L) and A n,y,Q (L) given by (17) (Figure 3b), and the mean ICXT power of the term A n,x,I (L) given by (14) (in dB) (Figure 3c) for a homogeneous WC-MCF ∆β c,mn = 0 rad · m −1 with L C = 10 m and L B = 10 m (regime L C ≈ L B ).  Figure 3a shows that the normalized mean of A n,x,I (L) is considerably small in the entire range of power distribution at the interfering core input. Figure 3b shows that a very low correlation between A n,x,I (L) and A n,y,Q (L) is observed regardless of the power distribution at the interfering core input. Figure 3c shows that the mean ICXT power of the term A n,x,I (L) has small variations (lower than 1.5 dB) with the power distribution at the interfering core input. Similar observations hold for the other ICXT field components. Figure 4 shows the simulation results for the normalized mean of the A n,y,I (L) field component (Figure 4a), the correlation coefficient between the field components A n,y,I (L) and A n,y,Q (L) (Figure 4b), and the mean ICXT power of the term A n,y,I (L) in dB (Figure 4c) with L C = 1 m and L B = 10 m (regime L C L B ).  Figure 4a shows that the normalized mean of A n,y,I (L) is considerably small in the entire range of power distribution at the interfering core input. Figure 4b shows also that a very low correlation between A n,y,I (L) and A n,y,Q (L) is observed regardless of the power distribution at the interfering core input. Figure 4c shows that the mean ICXT power of the term A n,y,I (L) has small variations (lower than 1 dB) with the power distribution at the interfering core input. Similar observations hold for the other ICXT field components.

Results for Single-Mode Homogeneous WC-MCF
In order to analyze the impact of L B reduction on the variation in the mean, correlation, and mean power of ICXT field components, we maintain L C = 10 m and reduce the beat length to L B = 0.1 m. Figure 5 shows the simulation results of the normalized mean of the A n,x,I (L) field component (Figure 5a), the correlation coefficient between the field components A n,x,I (L) and A n,y,I (L) (Figure 5b), and the mean ICXT power of A n,y,Q (L) in dB (Figure 5c) with L C = 10 m and L B = 0.1 m (regime L C L B ). Figure 5a shows that the normalized mean of A n,x,I (L) is considerably small in the whole range of power distribution at the interfering core input. Figure 5b shows a very low correlation between A n,x,I (L) and A n,y,I (L) regardless of the power distribution at the interfering core input. Figure 5c shows that the mean ICXT power of the term A n,y,Q (L) has small variations (lower than 0.7 dB) with the power distribution at the interfering core input. Similar observations hold for the other ICXT field components. In order to analyze the dependence of the power unbalance of the four ICXT field components, A n,p,q (L), for homogeneous WC-MCF ∆β c,mn = 0 rad · m −1 , on the power distribution at the interfering core input assigned to one of the two directions x or y and one of the two quadrature components I or Q, Figure 6 shows the unbalancing parameter U for the three previous pairs of values (L B , L C ) and ∆β c,mn = 0 rad · m −1 , as a function of ζ x,y , ζ I,Q in the range 0 ≤ ζ x,y ≤ 1 and 0 ≤ ζ I,Q ≤ 1. Figure 6 shows that values of the U parameter lower than 20 % are achieved for all situations tested. Such a low unbalancing level indicates that a quasi-uniform power distribution between the four components of the ICXT field occurs for homogeneous WC-MCF regardless of the power distribution among the four field components at the interfering core input.

Results for Quasi-Homogeneous Single-Mode WC-MCF
In this subsection, we carry out an analysis similar to the one performed in Section 5.1 but for quasi-homogeneous single-mode WC-MCF with ∆β c,mn = −701 rad · m −1 . We consider the three regimes: (i) L C ≈ L B , (ii) L C L B , and (iii) L C L B . Particularly, for L C ≈ L B , Figure 7 shows the simulation results for the normalized mean of the A n,y,I (L) field component (Figure 7a), the correlation coefficient between the field components A n,y,I (L) and A n,y,Q (L) (Figure 7b), and the mean ICXT power of the term A n,y,I (L) (in dB) (Figure 7c), assuming L C = 10 m and L B = 10 m. Similar observations to the ones indicated in Figure 3-5 for homogeneous WC-MCF can be settled from Figure 7. In fact, (i) Figure 7a shows that the normalized mean of A n,y,I (L) is considerably small in the entire range of power distribution at the interfering core input, (ii) Figure 7b shows that a very low correlation between A n,y,I (L) and A n,y,Q (L) is observed regardless of the power distribution at the interfering core input, and (iii) Figure 7c shows that the mean ICXT power of the term A n,y,I (L) has small variations (lower than 0.9 dB) with the power distribution at the interfering core input. Similar conclusions holds for the other field components. In order to analyze the dependence of the power unbalance of the four ICXT field components, A n,p,q (L), for quasi-homogeneous WC-MCF ∆β c,mn = −701 rad · m −1 , on the power distribution at the interfering core input assigned to one of the two directions x or y and one of the two quadrature components I or Q, Figure 8 shows the unbalancing parameter U for the three previous pairs of values (L B , L C ), and ∆β c,mn = −701 rad · m −1 , as a function of ζ x,y , ζ I,Q in the range 0 ≤ ζ x,y ≤ 1 and 0 ≤ ζ I,Q ≤ 1. Figure 8 shows that values of the U parameter lower than 25 % are achieved for all situations tested. Such a low unbalancing level indicates that a quasi-uniform power distribution between the four components of the ICXT field occurs for quasi-homogeneous WC-MCF regardless of the power distribution among the four field components at the interfering core input.
We emphasize that the small variations of the normalized mean results around 0 shown in Figures 3a-5a and Figure 7a Figure 7b would be reduced by increasing the number of simulated samples (due to the previous indicated reason) as well as by increasing the MCF length to lengths higher than 100 m due to the decorrelation field effect caused by the increase in the birefringence effects. These behaviors have been observed in a few other ICXT estimates obtained from numerical simulations with up to 500 samples and for WC-MCF lengths of up to 1000 m.

Conclusions
We assessed by intensive numerical simulation the mean, correlation, and power distribution of the ICXT field components of unmodulated polarization-coupled homogeneous and quasi-homogeneous single-mode WC-MCFs with a single interfering core considering three regimes: (i) L C ≈ L B , (ii) L C L B , and (iii) L C L B . For that, the CME model described in Reference [23] to assess the ICXT in a single-mode WC-MCF with intracore birefringence and linear propagation along the longitudinal coordinate of the MCF was used.
We have found that, as it has been reported elsewhere for scalar-coupled single-mode WC-MCFs, the normalized mean of the four ICXT field components and their correlations are very small, with a nearly uniform distribution of the mean power between the four ICXT field components and a very weak dependence on the power distribution of the four field components at the interfering core input. These conclusions were drawn for homogeneous and quasi-homogeneous single-mode WC-MCFs, with birefringence correlation length and birefringence beat length in the ranges of [0.5 m, 10 m] and [0.1 m, 10 m], respectively.
These results suggest that, for the range of birefringence correlation length and birefringence beat length typically found in singlecore telecommunication fibers, we may model the ICXT power fluctuations by a χ 2 distribution with four degrees of freedom in homogeneous and quasi-homogeneous polarization-coupled single-mode WC-MCFs, regardless of the distribution of power between the four field components at the interfering core input. Funding: This research was partially funded by FCT/MCTES through national funds and, when applicable, co-funded EU funds under the project UIDB/EEA/50008/2020 and project DigCore/UIDP/ 50008/2020.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflicts of interest.

Abbreviations
The following abbreviations are used in this manuscript: