An Analytical Solution for Fiber Optic Links with Photonic-Assisted Millimeter Wave Upconversion Due to MZM Nonlinearities

Photonic-assisted technique employed in the generation of stable and high-purity signals in the millimeter wave band has been investigated using dual-parallel Mach-Zehnder modulator. Indeed, this type of generation performs a high order RF multiplier, whose operation principle steams out from the nonlinear characteristic of the external modulators. Its application in directdetection fiber optic link was investigated based on an exact formulation that properly takes into account both the dispersion of a single mode fiber and the modulator nonlinearity. Analytic expression for the link transfer function dependence with respect to both modulator bias point and modulation index were obtained. It is worthwhile to mention that RF tuning is enabled since no optical filtering is required. As an illustration of the modeling usefulness, it is analytically shown that a 60 GHz millimeter signal, generated by photonic quadruplicating a 15 GHz signal using dual-parallel Mach-Zehnder modulator, can be transmitted through a 80 km link of standard single mode fiber with less than 0.1 dB deviation with respect to ideal flat response. This performance is much better than previously published results.


I. INTRODUCTION
End-user modern applications are requiring faster transmission rates and better signal stability, which could be mapped as improvements on network systems with wide bandwidth and low signal distortion.Increasing convergence in wireless services suggests a migration to millimeter wave spectrum, resulting on smaller antennas and greater beam directivity.These are essential requirements in on-board radar systems and satellite communication links, where size and weight are precious resources.Pure electronic signal generation in microwave and millimeter wave (MMW) bands is a complicated and relatively expensive process due to factors such as the complexity of the electronic modules, the extremely small size of some components, the high losses by propagation in wired lines or wave-guides, and the difficulty in isolating circuits from electromagnetic interference [1].Low distortion, high stability and low phase noise are desirable characteristics in a carrier signal used in An Analytical Solution for Fiber Optic Links with Photonic-Assisted Millimeter Wave Upconversion Due to MZM Nonlinearities Romildo Henrique de Souza, Olympio Lucchini Coutinho, José Edimar Barbosa Oliveira Instituto Tecnológico de Aeronáutica -ITA, São José dos Campos -SP, Brazil, romildo@ita.br,olympio@ita.br,edimar@ita.brAntônio Alves Ferreira Júnior, José Antônio Justino Ribeiro Instituto Nacional de Telecomunicações -INATEL, Santa Rita do Sapucaí -MG, Brazil, antonioa@inatel.br,justino@inatel.brwireless networks and radar systems for electronic warfare [2], for instance.
Photonic approaches could be employed to generate high frequency carrier signal with such criteria by up-conversion of lower frequency reference signals available from electronic sources.
Furthermore, photonic techniques offer tunability and relatively lower cost compared to purely electronic systems [3].Additionally, a fiber optic link allows providing the signal to a remotely located site before reconversion into the electrical domain, resulting in integration between optical processing and transmission (Fig. 1).Fig. 1.Schematic representation of a general photonic-assisted frequency multiplication system used in remote applications (O/E -Optical/Electrical).
Among several proposed photonic-assisted frequency multiplication systems, those using Mach-Zehnder Modulators (MZM) reveal being a convenient and scalable solution, whose operation principle steams out from the nonlinear characteristic of the external modulators.Researches for obtaining sixteen-fold frequency signals, with 20 dB of Sideband Suppression Ratio (SSR), represent state of the art in this area [4].Examples of experimental demonstrations include an octupler at 42.4 GHz, with electric SSR greater than 30 dB [5], and a quadrupler at 72 GHz with 36 dB of optical SSR [6].
Dual Parallel MZM (DP-MZM) topologies do not require optical filtering, since parameters are properly settled in order to suppress optical carrier and certain harmonics.Commercial availability of devices with two or more modulators encapsulated in the same package opens the possibilities to explore different configurations.For example, consider two push-pull modulators in parallel and then connected to a third MZM as presented in Fig. 2a.Two electrical signals in quadrature are applied each one to a modulator, which are both biased at the Maximum Transmission Bias Point (MATB), while the third MZM operates at Minimum Transmission Bias Point (MITB) [6], [7].An alternative topology (Fig. 2b) employs two dual drive modulators in parallel, eliminating the use of a third cascade MZM [8].Both modulators are biased at MITB and the electric input signals on the same modulator are in quadrature, while π out-of-phase from the other modulator.However, analysis of these topologies have been performed by using approximate methods restricted to small signal regime, where modulation index were arbitrarily limited to low values [6], [7], [8].This paper is mainly concerned with the modeling of the effect of the MZM nonlinearity on the performance of optical up-conversion.Special attention is given to the MMW generation which relies on Dual-Parallel MZM (DP-MZM) and exact expressions for the harmonic amplitude dependence on both modulation index and bias point were deduced [9].Based on such formulation we were able to precisely predict the harmonic distortion suppression ratio within a range of optical modulation index far wider than that used in previous publications [6], [7], [8].Additionally, a fiber optic spool is inserted on the modulator output for remote transmission.Under such conditions, chromatic dispersion influences system multiplication factor and thus must be considered on the final project design [10].
At this scope, the fundamentals of optical up-conversion are presented in Section II, which addresses the optical and RF spectra at the output of a Dual-Drive Mach-Zehnder Modulator (DD-MZM) at distinct bias point conditions.Section III deduces the wave propagation equation in fiber optic links under chromatic dispersion and attenuation effects and relates them with a properly transfer function to account amplitude and phase response for each optical spectral component.A DD-MZM combined with a spool of fiber is the simplest frequency multiplication topology, which is rigorous analyzed in Section IV, aiming to establish the exact solution methodology in order to realize the effect of fiber chromatic dispersion on the multiplication factor.Section V applies the formulation on a dual parallel MZM scenario, resulting in higher multiplication factors and consequently reaching millimeter wave generation.Thus, this exact solution is applied to a back-to-back (BTB) system in (a)

(b)
Section VI, in order to compare with approximated solution found in the literature.Finally, conclusion and brief discussion are presented in Section VII.

II. FUNDAMENTALS OF DD-MZM UP-CONVERTER
Although the field of electro-optic modulation is well established, its application for accomplishing optical MMW signal generation, in which the MZM works as a frequency multiplier, had begun rather recently.Therefore, in order to be able to properly model such device, we consider worthwhile to present a brief review of the main characteristics of DD-MZM's operation.To this end, the MZM electronic drives should be assembled in such a way to enable its operation at any desired bias point on its transference curve.Then, the DD-MZM will be part of a fiber optic link, as depicted in Fig. 3.
The operation of such type of modulator relies on the interference between two optical guided modes whose propagation constant are controlled by a RF electric field pattern conveniently distributed inside an electro-optic substrate.The modulating electric field pattern, which is spatially formatted by thin film electrodes, behaves as quasi-TEM wave.It depends on RF reflection coefficients at the ends of the electrodes and on the optical and RF velocities mismatching (see Appendix for further detail).
The field pattern is driven by a combination of DC and RF sources, which provides the bias and RF modulating field, respectively.In this publication, it is assumed that the impedance and velocity matching are simultaneously attained, hence the RF phase shift induced by the voltage waveforms seen by the optical wave is given by (A-6).Under such condition, the RF inputs of the DD-MZM shown in Fig. 3  With respect to the optical carrier shown in Fig. 3, it is assumed that it steams out from a DFB laser with null linewidth and no intensity or phase noise, which yields an optical power  0 and corresponding electrical field strength  0 and frequency  0 .
Assuming a 50/50 splitting ratio and no insertion loss on the modulator, following rather standard procedures, the electric optical field at the DD-MZM output is obtained as: where An examination of (3), ( 4) and ( 5) reveals that, in general, the output signal has a combination of amplitude and phase modulation.However, when both RF phase shifts ( 1 and  2 ), and DC bias ( 1 and  2 ) are properly chosen, it is possible to achieve different modulation formats which ranges from pure intensity to pure phase modulation.Moreover these generalized expressions turns out to be very helpful when we undertake the modeling of a DP-MZM, which is made up of two adequately electronically driven DD-MZM's that are optically associated in a parallel configuration.Beyond that, bearing in mind that the DD-MZM's performance evaluation as an up-converter requires the knowledge of its spectra characteristics, i.e., both the optical and RF spectrum, the setup illustrated in Fig. 3 includes both an optical spectrum analyzer (OSA) and an electronic spectrum analyzer (ESA).
Using Jacobi-Anger expansion in (3), the electrical field can be represented as: with where   () is the n-th Bessel function of first kind.
As previously stated, this section aims at the modeling of an electro-optical up-conversion based on DD-MZM, and validating the obtained results through their comparison to well accepted results available in the literature, whenever it is possible.To this end, we have to particularize the schematic depicted in Fig. 3, in a way that it reduces to well established configurations [11], [12].Such references present valuable results regarding to both experimental and theoretical on the DD-MZM modulator performance when it operates either in single sideband (SSB) or double sideband (DSB) regime.
In order to be able to compare our results to those given in [12], we particularize the configuration shown in Fig. 3 by imposing  2 =  2 = 0.Under these conditions, (4) and ( 5) become: Furthermore, the above given equations provide the following expression for the optical power at the output of the DD-MZM, i.e., in the so-called back-to-back (BTB) situation: Equation ( 11) besides revealing the nonlinear dependence of the DD-MZM response with respect to the driving electronic signal, also points out that it can be strongly affected by the bias point, which is specified by a combination of phase shift and DC bias, namely the pair ( 1 ,  1 ).In general, the literature emphasizes the conditions which yield maximum transmission (MATB), quadrature (QB) and minimum transmission (MITB).In this paper we are interested in those situations in which ( 1 ,  1 ) meet the following requirements: (π, 0) for MATB; (π, π/2) for QB; and (π, π) for MITB.The localizations of such points over the modulator characteristic curve are shown in Fig. 4. Once again using the Jacobi-Anger expansion in (11), we obtain the optical power spectral component at frequency  0 +   :

Optical Transmittance
The result given by ( 12) has been widely used in both theoretical and experimental investigations of high speed electro-optic modulators performance [11].
Before undertaking any further analysis, let us point out a few features of the dependence of the optical power spectrum with respect to a chosen bias point, as predicted by (12).First of all, one realizes that both the amplitude of a given optical spectral component and the spectrum symmetry with respect to the carrier can be controlled by properly choosing the bias point, defined by ( 1 ,  1 ), as shown in Fig. 3.As an illustration, we choose ( 1 ,  1 ) = (π/2, π/2), which implies that the +1 st order sideband is null, whereas the -1 st order sideband is given by  0  1 2 ().Thereby, such bias point yields the so-called single sideband regime (SSB).A rather different feature is obtained when we choose ( 1 ,  1 ) = (π, 0), because the ±1 st order sidebands, are suppressed while the remaining components are all separated by even multiple of the MMW frequency   .This configuration is called double sideband (DSB) and its spectra, at both optical and RF domains, are depicted in Fig. 5. Now, we return to (12) aiming to predict the electrical current at the output of a photodetector (PD), when the setup shown in Fig. 3 is rearranged to perform BTB measurement.Obviously, under such assumption the fiber optic chromatic dispersion is not taken into account, therefore the MMW waveform and its corresponding spectral composition can be readily determined based on (12).One concludes that the electric current waveform is given by the following expression: where ℜ is the PD responsitivity.
In previous publication we have used the Jacobi-Anger expansion in (13) and then carried out a general analysis of the dependence of electrical signal spectral content, at the output of a DD-MZM, with respect to the bias point, i.e., ( 1 ,  1 ) [13].Such undertaking, beyond enabling a quick retrieve of a wide range of results available in the literature regarding to the single tone operation of DD-MZM, also provided valuable insights on the subject of bias control techniques for lithium niobate modulators [13].Part of the obtained results which are concerned to the optical and electronic power spectra with the DD-MZM operating in the DSB regime are given in Table I.This type of configuration played a major role during the earliest years of optical research on microwave optical generation by electro-optic up-conversion [14].It is worthwhile to point out that (13) can be rewritten in the rather familiar expression shown below: which provides a manner for quick double-checking of all predictions so far presented in this publication.Such features will modify the resulting electrical signal at the photodetector output and decrease system dynamic range.Once physical properties of dispersive materials depend on the transmitted signal frequency, phase factor β(ω) can be expanded in Taylor series around optical carrier angular frequency ω 0 [15]: Relying on slowly varying envelope approximation, the solution of the electric field wave equation assuming z-propagation and x-linear polarization is [15]: where (, ,  0 ) is the electrical field spatial distribution at frequency ω 0 and (, ) is the slowly varying wave envelope.
Plugging (16) in the wave equation and using the expanded phase factor β(ω), given by ( 15), first order perturbation theory leads to [15]: where are the phase factor coefficients,   is the attenuation factor,  is the nonlinear parameter due to Kerr effect and R(t) is the nonlinear response function due to Raman scattering, considering it as the most intense nonlinear phenomena in that case.
Let us assume intensity modulation of an optical carrier with angular frequency ω 0 by an electric signal ω RF , being transmitted through a fiber optic link.The electric field spectral components of the optical signal will be at frequencies  0 ±   , for any integer n, in a DSB modulation.Without losing accuracy and following general approach, the power series in (15) could be limited until the third term to obtain a reasonable phase factor for each spectrum component: where  0 = ( 0 ) gives rise to a fixed phase shift; is the inverse of the group velocity, which makes the second term in (18)  Neglecting fiber nonlinearities, once it is out of the scope of this paper, we set  = 0 in (17): As β 1 relates to a constant propagation delay, it turns to be convenient to apply the following transformation of variables [15],  =  −  1  and  = , in order to simplify (19): The above equation describes electromagnetic wave propagation in fiber optic media and could be denoted by a transfer function H(ω) for a spool of fiber with length  = .Each spectral component has different propagation delays, i.e., phase factor depends on the component frequency, which leads to distinct phases for each component.Thus, transfer function H(ω)is defined by an amplitude response U(ω) and by a phase response ϕ(ω): Although silica attenuation varies with signal frequency, here it will be assumed as a constant inside the 1550 nm optical transmission window.Doing so, amplitude response U(ω) becomes: −   20 (22) where   is the fiber attenuation factor in dB/km.
On the other hand, phase response for each optical spectral component of the modulated signal after propagation through a L-length fiber is: It could be noted that there is a quadratic relation between phase and frequency   for each component,   =     2 , where the proportionality factor   is given by: Rewriting (24) in terms of the linear dispersion parameter D, we have: where  0 is the optical signal wavelength and c is the speed of light.
On a direct detection receiver, quadratic photodetector output has a RF fundamental component which electrical power  1 is proportional to  2 ( 1   2 ).This is null when argument is a π⁄2 odd multiple [16].Therefore, when a signal with frequency   is launched into a fiber optic link, the fundamental electric power component will ideally vanish at a first distance  1 =  2 0 2   2 and then at each distance interval Δ =   0 2   2 .The electric power flows from fundamental component to higherorder harmonics and this phenomenon is known as signal fading [13], [17].Fig. 6 presents a numerical example of the variation of the fundamental electrical power on the photodetector output when 30 GHz and 60 GHz signals are modulated by a 1550 nm optical carrier and transmitted through a standard single mode dispersive fiber (D = 17 ps/nm-km).
Fig. 6.Example of the fundamental electrical power dependence on fiber length for 30 GHz and 60 GHz signal due to chromatic dispersion effect on an intensity modulation/direct detection fiber optic link [17].

IV. MZM UP-CONVERSION IN FIBER OPTIC LINKS
In order to develop the link model which takes into account the fiber chromatic dispersion, we should emphasize some assumptions we are taken for granted with respect to link setup shown in Fig. 3. First, as far as the optical signal is concerned, all components behave linearly, except the PD which displays quadratic response; secondly, the source of RF nonlinearity is fully concentrated into the DD-MZM; and thirdly, the chromatic dispersion steams out from the optical fiber [11], [12].Additionally, due to our interest in electro-optic up-conversion, it is convenient to DD-MZM to operate at MATB bias point.
Under the above mentioned assumptions, we used ( 7) and ( 8) to determine the optical electrical field in the time domain, right at the DD-MZM's output.Using this result in combination with Maxwell's equations, we undertook the time domain analysis of the propagation of the optical field along a dispersive and linear spool of SMF.In order to take into account the chromatic dispersion's effect, we will consider ( 21), ( 22) e ( 23) to find the electrical field at the optical link output: Assuming that the optical power is evenly distributed over the PD cross section, the time dependent RF current turns out to be: where   is a constant (in ohms per square meter) which depends on both the fiber effective crosssection and the optical wave impedance, and mathematical symbol * denotes the complex conjugate.
Substituting (26) in (27), photodetector current becomes: where each current harmonic component is: for   and  + being, respectively, the phase of components   and ( + )  .
Clearly, harmonic electric current components depend on Γ parameter and modulation index m, which are function of optical fiber parameters (length and dispersion coefficient), modulator halfwave voltage and bias voltage, and RF signal characteristics (amplitude, frequency and phases).Laser power, photodetector responsitivity and fiber attenuation are common terms for all the harmonics.Dominant harmonic and the relative intensity to the other harmonics are defined by this set of parameters.Adequate scenarios give rise to an output signal which frequency is a multiple of the reference signal frequency.
For DSB modulation (∆ =  1 = , ∆ =  1 ,  2 =  2 = 0), (32) simplify to [9]: In this case, there are only even harmonics and the lowest order multiplier is a doubling.Fig. 7a e 7b show the relative RF power for second and fourth harmonics as a function of both Γ parameter and modulation index m, for DSB modulation and MATB operation.Given a modulation index, it is possible to analyze system performance in terms of Γ parameter by cutting the 3D-graphs, as done in Fig. 8a and 8b, for  = 1 and  = 2.5, respectively.On the first cut-plot, second harmonic is more intense than the fourth in a large Γ interval and the difference between them reaches 20 dB.However,  the 60 GHz millimeter wave band (57-64 GHz) [18].In such frequency range it is widely recognized that the generation of high purity and low phase noise MMW signals using conventional electronic technology constitutes an expensive and challenging undertaking.As a consequence, in the past decade a great deal of research has been guided towards optical millimeter wave generation based on MZM's up conversion [6], [7], [8].
Signal generation in millimeter wave band using up-conversion techniques requires high multiplication factors.A photonic approach system based solely in DD-MZM is not adequate to reach that spectrum with easily available reference signals.Multiplication factors could be enhanced by associating DD-MZM in cascade and/or in parallel.Several topologies have being described in literature.While the simple ones control the basic parameters of the modulators in the system, complex configurations could employ Fiber Bragg Gratings (FBG), Erbium Doped Fiber Amplification (EDFA), Semiconductor Optical Amplifier (SOA) or other electrical or optical devices [3].The main goal is trying to suppress certain optical spectral components, resulting more clearance between RF harmonics at the output.As mentioned before, this natural spacing reduces the need of optical filters and improves system performance.Among such experiments, steams out those dealing with frequency quadrupling based on an integrated dual parallel MZM (DP-MZM) [6], [7], [8].The device's main configuration comprises three sub-MZM assembled as illustrated in Electrical Optical EMZM(t) (up)   EMZM(t) (down)   Fig. 9.The DD-MZM embedded in each arm of the main MZM exhibits features which were fully discussed in Section II.To the best of our knowledge, the design of frequency quadrupler based on DP-MZM has been developed relying on models which does not fully takes into account the MZM's RF nonlinearity.For example, one of our paper main motivation is to examine the assumption which claims the existence of a certain modulation index depth beyond which no closed mathematical expression can be obtained for the RF output power [6], [7], [8].To this aim, we represent the electric field at the DP-MZM output as the sum of the individual electric field from each DD-MZM (up and down): where After some algebraic operations, we get [9]: As expected, the RF nonlinearity effects on the optical field at DP-MZM's output, as given by (37), seems to be far more interesting than that previously investigated in Section II, which deals with DD-MZM.In order to properly understand this statement one should apply the Jacobi-Anger expansion to rewrite (37) as a complex Fourier series: where Then, based in (39), it can be concluded that at the DP-MZM's output the optical carrier is fully suppressed.Beyond that, all the remaining   coefficients are null, except those whose order obeys the relationship  = 4 − 2, for  = 0, ±1, ±2, ±3 ….In another words, the optical spectrum only contains components at frequencies such as  0 ± 2  ,  0 ± 6  ,  0 ± 10  , and so on (Fig. 10a).By comparing these characteristic to those illustrated in Fig. 5a, which are attached to DD-MZM, steams out one of the main reason why DP-MZM is better suited for frequency quadrupling.DC 4ω RF 8ω RF 12ω RF 16ω RF ω 0 +2ω RF ω 0 -2ω RF ω 0 +6ω RF ω 0 -6ω RF ω 0 -10ω RF ω 0 +10ω RF Fig. 10.Frequency response (the amplitude are not to scale) of a DP-MZM: (a) optical spectrum and (b) electronic spectrum.
Following the methodology described on previous section by considering a spool of fiber optic with length L, we obtain an expression for the exact current harmonic component at the photodetector output, in a similar way than (30): Plugging (39) in (40) and applying once again Graf's Addition Theorem, we have [10]: The term    2 is non-zero only when N is a multiple of 4. Therefore, the output electrical current will have components in the frequencies 4  , 8  , 12  , and so on (Fig. 10b).Once again, harmonic electric current components depend on Γ parameter and modulation index.The 3D plots of the fourth and eighth harmonics are presented in Fig. 11a and 11b.Cuttings on fixed m allow identify the dominant harmonic and thus, the system multiplication factor.At  = 1, Fig. 12a reveals that fourth harmonic is at least 74 dB higher than the eight harmonic, resulting in a high-purity signal with a large range insensibility to variations of link parameters.It is also noted a performance improvement when this DP-MZM quadrupling is compared to the DD-MZM doubling, by means of Fig. 12a and 8a: fading effect was very pronounced and then it is practically eliminated, for same modulation index.As an illustrative example, a 60 GHz MMW signal could be generated from a 15 GHz signal using this DP-MZM scheme, with less than 0.1 dB deviation with respect to ideal flat response of a standard single mode 80 km fiber optic link.This is a great length improvement when compared to (a) (b) the 1 km fading noticed in Fig. 6 [17].Increases in modulation index will raise intensity of the other harmonics, as could be seen in Fig. 12b, for  = 5.

VI. DUAL PARALLEL MZM FREQUENCY QUADRUPLING
It is opportune to point out that a great deal of previously published results on quadrupling are based on time domain analysis which leads to expression for the optical field in the form of infinite series [6], [7], [8], [14].This fact makes the calculations very cumbersome, even when one attempts to fully take into consideration the RF nonlinearity at moderated level of optical modulation [12], [13].Due to that, calculations are performed under the assumption that the series has a finite number of terms.Although a few approximations have been widely used, still there is the need for a more rigorous criterion that leads to an optimum value for the number of terms of the series.
In this publication we tackled the matter above mentioned using the theoretical background described in the previous section.Setting  = 0 in (41), we exclude fiber optic from the system and (41) resumes to an exact and very concise expression for the current components at the PD's output [9]: First of all, we emphasize that (41) and (42) do not relay in any assumptions about the optical modulation index.Furthermore, they point out that the RF current only has spectral components at certain given frequencies, namely: 4  , 8  , 12  , and so on.Such predictions are in full agreement with previously published values [6], [7], [8].Moreover, the rather wide frequency range between two consecutives spectral components, which is twice as larger than that previously presented for the DD-MZM setup, constitutes one of the reasons why the DP-MZM configuration has become attractive as frequency up-converter.Regarding the magnitude of the already mentioned spectral components, their approximated expressions, presented by different authors, and the exact one, (42), are shown in Table II.Our numerical calculations regarding to the DP-MZM up-converter were constrained to BTB configuration, and were undertaken based on approximated and exact analytical expressions given in Table II.Since the exact expressions do not rely on any mathematical assumption with respect to the RF power drive, we developed the calculations for modulation index as high as 13, as it is shown in Fig. 13.Actually, the various approximations introduced in [6], [7] and [8] still need a clear cut criterion to determine an upper bound limit for the modulation index.According to Fig. 13a, the two formulations are in very good agreement for modulation index as high as 5.0.For example, for  = 2.4, the RF Sideband Suppression Ratio (RFSSR) using the exact and approximated relationships are found to be 36.18dB and 36.10 dB, respectively.The discrepancy between the two predictions only becomes relevant when the index exceeds 7.0.These results reveal that the limit of validation of the approximation lies far beyond 2.4, as have been previously stated [6], [7], [8].Similarly, the comparison between higher order harmonics amplitudes predicted by the two formulations also exhibits very good agreement for modulation index even higher than 5.0, as shown in Fig. 6b

Fig. 5 .
Fig. 5. Frequency response (amplitudes not in scale) of a DD-MZM operating at MATB in the DSB regime: (a) optical spectrum; (b) electronic spectrum.

Fig. 7 .
Fig. 7. Amplitude of harmonics generated in a DD-MZM up-converter in a fiber optical link set for DSB modulation and MATB bias: (a) 2 nd and (b) 4 th harmonic dependence with respect to Γ parameter and modulation index m.

Fig. 8 .
Fig. 8. Cuts in the amplitude 3D-plots of 2 nd , 4 th and 6 th harmonics for: (a)  = 1 and (b) = 2.5 .V. EMPLOYMENT OF DP-MZM TO ENHANCE MULTIPLICATION FACTOR IN FIBER OPTIC LINKSNowadays, innumerous applications in the area of telecommunication are being developed inside

Fig. 11 .
Fig. 11.Amplitude of harmonics generated in a DP-MZM up-converter in a fiber optical link: (a) 4 th and (b) 8 th harmonic dependence with respect to Γ parameter and modulation index m.
and Fig.6c.VII.CONCLUSIONWireless networks and radar systems are some of the major technological applications that rely on microwave and millimeter wave carrier signals.Photonic-based frequency multiplication topologies show up as tunable and reliable techniques.This paper presented the effects of chromatic dispersion in fiber optic links associated with two Mach-Zehnder Modulators in parallel, the so-called Dual-Parallel MZM, to remotely provide microwave and millimeter wave signals from lower frequency reference signals.It is suitable to suppress certain optical spectral components through exploring the nonlinear characteristic of the external modulators.The result is sufficiently spaced RF components without using narrow optical or electrical filters.Analytical exact solution for the photodetector output current was derived with no restriction to small signal analysis.It allows a more accurate once harmonic intensities could be precisely described by means of link parameters and modulation index.It has been proved that an up-conversion remote system can be achieved without considerable fading penalties due to fiber dispersion, by properly choosing link parameters.Contrary to previously published formulations, which rely on time domain analysis, the analysis here presented was carried out in the frequency domain and yielded close expressions for the RF power at the DP-MZM's output.The obtained results reveal that previously published formulations do not properly take into account the RF nonlinearity.Thereby, one of the main advantages of the formulation is the fact that it enables straightforward calculations of parameters such as RFSSR, which otherwise has been only approximately calculated based on a wide range of methods.Phase noise and fiber nonlinearities effects are the next topic to be considered in the proposed model, along with practical experiments on future research.APPENDIX Nowadays, broadband travelling-wave electro-optic modulator based on lithium niobate (LiNbO 3 ), like MZM modulator, still play a major rule in high-speed and long-haul state-of-the art optical fiber communication networks.It is well known that RF frequency response of this type of modulator strongly depends on parameters such as the RF reflection coefficients, mismatching between optical and RF velocity, and electrodes losses, among others[19]-[24].Quite often one assumes that the mentioned constraints are simultaneously overcome, and reaches the rather awkward conclusion that the link performance does not depend on the modulator bandwidth.This appendix attempts to review models which accounts for effect of the previously identified mismatching issues on the modulator RF performance, based on unified nomenclature which enables a ready tracing back of results published by many authors.Before undertaking the analysis, we would like to point out that LiNbO 3 is a dispersive material which exhibits an inherent mismatch between optical and microwave velocities, whereas the microwave attenuation and impedance mismatching are strongly dependents on electrodes geometries.Bearing in mind that a MZM modulator is made up of phase modulators (PM) properly connected, our model relies on the standard schematic representation shown in Fig. A-a.In such illustration,   () stands for the RF drive voltage at the electrodes input,   is the internal impedance of signal generator,   is the load impedance at the electrodes right end, and SiO 2 is the buffer layer.It is well known that the travelling-wave structures under investigation can support a quasi TEM mode, which can be represented by an equivalent transmission line, as shown in Fig. A-b.
1 and  2 the DC bias voltage.

TABLE I
be linear in frequency and consequently results on a fixed propagation delay time;  2 =  2 ( 0 ), so third term leads to a linear propagation delay in frequency.It results on a distortion in the final waveform.