Photonic True-Time-Delay Beamformer For A Phased Array Antenna Receiver Based On Self-Heterodyne Detection
Photonic True-Time-Delay Beamformer For A Phased Array Antenna Receiver Based On Self-Heterodyne Detection
fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                1
             This paragraph of the first footnote will contain the date on which you                                V. C. Duarte is with Instituto de Telecomunicações and Physics Department,
          submitted your paper for review. This work is co-funded by FCT/MEC through                             University of Aveiro, 3810-193 Aveiro, Portugal (e-mail: vanessaduarte@
          national funds and when applicable co-funded by FEDER – PT2020 partnership                             av.it.pt).
          agreement under the project UID/EEA/50008/2013 and also supported by the                                  M. V. Drummond is with Instituto de Telecomunicações, University of
          European Commission through the project BEACON (FP7-SPACE-2013-1-                                      Aveiro, 3810-193, Aveiro, Portugal (e-mail: mvd@av.it.pt).
          607401).                                                                                                  R. N. Nogueira is with Instituto de Telecomunicações, University of Aveiro
                                                                                                                 and Watgrid Lda., 3810-193, Aveiro, Portugal (e-mail: rnogueira@av.it.pt).
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                2
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                3
                                   𝑣(𝑡) = 𝑎(𝑡) sin(𝜔RF 𝑡 + 𝜃(𝑡)),                                            (5)   where 𝐸MZDI (𝑓) is the delayed output optical signal, 𝑀(𝑓) is
                                                                                                                   the second output optical signal, which can be used for
          where 𝑎(𝑡) and 𝜃(𝑡) are the amplitude and phase variation of                                             operation monitoring, 𝜏(𝑓) is the transfer function of the optical
          the input RF signal, respectively. For a matter of simplicity, let                                       delay line and 𝐸(𝑓)𝑒𝑗𝛽 is the input optical signal. Without loss
          us consider infinite ER, i.e., 𝑘ER = 1, and an RF signal with                                            of generality, the amplitude and group delay responses of the
          reduced amplitude, such that the MZM operates in its linear                                              TODL can be derived by setting 𝛾 = 𝜋⁄2, as the variation of 𝛾
          domain. Equation (4) can be simplified to                                                                only results in a displacement of the frequency response along
                                                            1 𝜋                                                    the frequency. Such responses are given by
                                     𝐸(𝑡) = 𝐸in (𝑡) ⋅              𝑣(𝑡).                                    (6)
                                                           √2 𝑉𝜋                                                                                   1
                                                                                                                               ‖𝐻(𝑓)‖2 = (sin(𝜙) cos(2𝜋𝑓𝜏) + 1),                                              (9)
                                                                                                                                                   2
          B. Tunable Optical Delay Line                                                                                                           𝜏 sin 𝜙 cos(2𝜋𝑓𝜏)−cos 𝜙+1
                                                                                                                                  𝜏𝑔 (𝑓) = (                                         ).                     (10)
             After being phase shifted by 𝛽, the modulated optical signal                                                                         2      sin 𝜙 cos(2𝜋𝑓𝜏)+1
          is delayed by the TODL. The TODL output signals are given
          by                                                                                                       Fig. 2 shows the amplitude and group delay responses for
                                                                                                                   different values of 𝜙. Both responses have a periodical behavior
                     𝐸    (𝑓)   1 1 𝑗 𝑒𝑗𝜏(𝑓)               0 ] [1     𝑗 𝑒𝑗𝜙       0] [1 𝑗 ]
                    [ 𝑀𝑍𝐷𝐼 ] =    [   ][                               ][                                          with a period of 1/𝜏. The amplitude response is identical for
                       𝑀(𝑓)    2√2 𝑗 1   0                𝑒𝑗𝛾 𝑗       1 0         1 𝑗 1
                                                                                                            (7)    supplementary values of 𝜙, whereas the group delay response
                                                  𝑗𝛽
                                          ⋅ [𝐸(𝑓)𝑒 ] ,                                                             increases with 𝜙, proving that the delay imposed by the MZDI
                                                0
                                                                                                                   depends on the power coupling ratio between the arms of the
          resulting in                                                                                             MZDI.
                     𝐸    (𝑓)   1   𝑒𝑗𝜏(𝑓) (𝑒𝑗𝜙 − 1) − 𝑒𝑗𝛾 (𝑒 𝑗𝜙 + 1)                                                 As depicted in inset C of Fig. 1, the modulated optical signal
                    [ MZDI ] =    [                                       ] ⋅ 𝐸(𝑓)𝑒𝑗𝛽 ,                     (8)    has two sidebands, one of which is heterodynally detected.
                       𝑀(𝑓)    2√2 𝑗𝑒 𝑗𝜏(𝑓) (𝑒 𝑗𝜙 − 1) + 𝑗𝑒 𝑗𝛾 (𝑒 𝑗𝜙 + 1)
                                                                                                                   Without loss of generality, it is assumed that such sideband is
                                                                                                                   the one at the highest frequency. In order to have the MZDI
                                                                                                                   correctly operating as a TODL, one must center such sideband
                                                                                                                   at a maximum of the amplitude response, i.e., 𝑓RF ∙ 𝜏 = 𝑛,
                                                                                                                   where 𝑛 is an integer. If such is the case, the delayed optical
                                                                                                                   signal can be approximated by
                             -1                                                                                    where 𝑉𝜋,up , 𝑣up (𝑡) and 𝑉B,up are the half-wave voltage, RF
                                                                                                      () / 
                                                                                               0.6
                                                                                               0.4                 input signal and bias voltage of the upper MZM, respectively,
                                                                                                            g
                             -2
                                                                                                                   and 𝑉𝜋,low , 𝑣low (𝑡) and 𝑉B,low are the same variables of the
                                                                                               0.2
                                                                                                                   lower MZM. 𝑉𝜋,PS and 𝑉B,PS are the half-wave voltage and bias
                             -3                                                                0
                               0         50              100                    150                                voltage of the optical phase shifter. Considering VB,up =
                                                     [deg]                                                                        𝑉𝜋      𝜋     𝜋                 𝜋
                                                                                                                   𝑉B,low = −            ⋅ ,         𝑉           = , 𝑉𝜋,up = 𝑉𝜋,low = 𝑉𝜋,IQ , and
           Fig. 3. Attenuation and group delay response of the TODL as a                                                            𝜋      2    𝑉𝜋,PS B,PS         2
           function of 𝜙.                                                                                          assuming reduced a modulation index, (12) results in
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                4
                             1
                𝐼(𝑡) = − 𝑅𝐺sig 𝐺OLO 𝑃in
                                                         𝜋2
                                                                𝐴IQ 𝑎(𝑡) ⋅ (cos((𝜔RF −                           Equation (23) shows that besides the desired tone at 𝜔RF − 𝜔LO
                             4                       𝑉𝜋 𝑉𝜋,IQ                                        (18)        and the undesired tone at 𝜔RF + 𝜔LO , there is now another
               𝜔LO )𝑡 + 𝜃(𝑡) + 𝛽) − cos((𝜔RF + 𝜔LO )𝑡 + 𝜃(𝑡) − 𝛽)).                                              undesired tone at 𝜔LO . In order to avoid such that tone falls at
          The undesired tone at the highest frequency of 𝜔RF + 𝜔LO can                                           𝜔RF + 𝜔LO , one must have
          be filtered out through low-pass filtering, leading to
                                                                                                                                               𝜔RF − 𝜔LO ≠ 𝜔LO ,                                            (24)
                        1                  𝜋2                                                                                                                              𝜔RF
                𝐼(𝑡) = − 𝑅𝐺sig 𝐺OLO 𝑃in         𝐴 𝑎(𝑡)                                                           which, once again, rules out 𝜔LO =         . An identical analysis
                        4               𝑉𝜋 𝑉𝜋,IQ IQ                                                  (19)                                                 2
                                 ⋅ cos((𝜔RF − 𝜔LO )𝑡 + 𝜃(𝑡) + 𝛽).                                                can be performed for the IQM, in which a finite ER results in
                                                                                                                 the OLO having a residual tone also at 𝜔o . After coherent
          Equation (19) shows that the original RF signal is frequency                                           detection, besides the tones at 𝜔RF ± 𝜔LO and at 𝜔LO , there is
          converted from 𝜔RF down to 𝜔RF − 𝜔LO . It also shows that the                                          another very weak undesired tone at DC. Such tone can be
          optical phase shifting of 𝛽 is translated into an identical RF                                         easily filtered out using bandpass filtering.
          phase shifting due to the rejection of the received highest                                               Besides harmonic generation and finite ER, one must also
          frequency tone.                                                                                        consider the unavoidable drifting of the bias voltage of a
                                                                                                                 modulator, which cannot be neglected even when it is kept
          E. Impairments                                                                                         under control. Small but realistic drifts can be straightforwardly
             The model derived above assumes an ideal behavior of all                                            modeled in (6) as
          devices. In practice, the impairments caused by them need to be
                                                                                                                                                           1 𝜋
          taken into account in order to determine the robustness of the                                                          𝐸(𝑡) = 𝐸in (𝑡) ⋅               𝑣(𝑡) + Δ𝑉B ,                               (25)
                                                                                                                                                           √2 𝑉𝜋
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                5
          where Δ𝑉B is the deviation of the bias voltage of the MZM from                                         The first term introduces undesired tones at DC and 2𝜔RF . The
          the ideal value. This term, after coherent detection, produces an                                      second term only produces a tone at DC, whereas the third term
          undesired tone at 𝜔LO , similarly to the one produced due to the                                       represents the tones at 𝜔RF ± 𝜔LO . Consequently, in order to
          MZM having a finite ER. The same analysis is valid for drifts                                          suppress the first and second terms one must have 𝜔RF −
          in bias voltages of the IQM. Therefore, by carefully choosing                                          𝜔LO ≠ 0.
          𝜔IF in compliance with (21) and (24), third-order harmonic                                                Self-heterodyne detection has the advantage of enabling
          generation, finite ER of the modulators, and drifting of the bias                                      phase noise cancellation, as observed in (17). However, in the
          voltages of modulators do not impair the output signal, as the                                         derivation of (17) is assumed that the TODL does not delay the
          resulting spurious tones can be rejected by bandpass filtering.                                        modulated optical signal. However, the modulated optical
             In order to achieve ideal frequency shifting of the input laser                                     signal is delayed by 𝜏𝑔 (𝜙), 0 ≤ 𝜏𝑔 (𝜙) ≤ 𝜏, and therefore
          signal using a DP-MZM, it is required that all electrical and
                                                                                                                   ∗
          optical signals are perfectly balanced, which is unrealistic. In                                        𝐸in (𝑡 − 𝜏𝑔 (𝜙)) 𝐸in (𝑡) = 𝑃in exp (𝑗 (𝜔o 𝜏𝑔 (𝜙) + 𝜙PN (𝑡) −
          practice, even though most power is located at the desired tone                                                                                                                                   (31)
          𝜔o + 𝜔LO , there is a residual tone at 𝜔o − 𝜔LO that cannot be                                                                      𝜙PN (𝑡 − 𝜏𝑔 (𝜙)))).
          fully suppressed. The presence of the residual tone at 𝜔o − 𝜔LO
          can be modeled by generalizing (14) to                                                                 The first term results in a constant phase shift, and can be easily
                                     1      𝜋                                                                    accounted for by adjusting the optical phase shift 𝛽. The second
                𝐸OLO (𝑡) =                    𝐸in (𝑡)                                                            and third terms show that ideal phase noise cancellation is not
                                    2√2 𝑉𝜋,IQ
                                         ⋅ 𝐴IQ (𝑘SBR exp(𝑗𝜔IF 𝑡)
                                                                                                     (26)        achieved, as in practice 𝜙PN (𝑡) − 𝜙PN (𝑡 − 𝜏𝑔 (𝜙)) ≠ 0.
                                         + (1 − 𝑘SBR ) exp(−𝑗𝜔LO 𝑡)),                                            Consequently, the worst-case scenario corresponds to having
                                                                                                                 the highest possible delay, i.e., 𝜏(𝜙) = 𝜏. Instead of resorting
          where 𝑘SBR is related with the ratio between the amplitudes of
                                                                                                                 to a complex statistical assessment of 𝜙PN (𝑡) − 𝜙PN (𝑡 −
          the upper and lower sidebands, 0 ≤ 𝑘SBR ≤ 1, and ideally
          𝑘SBR = 1. After coherent detection and low-pass filtering, the                                         𝜏𝑔 (𝜙)), a simple insight can be obtained using the concept of
          output current is given by                                                                             laser coherence time. The coherence time is the time during
                                     1                         𝜋2
                                                                                                                 which the laser signal may be considered coherent, and is given
                     𝐼(𝑡) = − 𝑅𝐺sig 𝐺OLO 𝑃in                           𝐴IQ 𝑎(𝑡) ⋅                                by
                                     4                      𝑉𝜋 𝑉𝜋,IQ
                 (𝑘SBR cos((𝜔RF − 𝜔LO )𝑡 + 𝜃(𝑡) + 𝛽) − (1 −                                          (27)
                                                                                                                                                              1
                    𝑘SBR ) cos((𝜔RF − 𝜔LO )𝑡 + 𝜃(𝑡) − 𝛽)).                                                                                          𝜏CT =        ,                                          (32)
                                                                                                                                                              Δ𝜈
          Equation (27) shows that the residual tone of the OLO directly                                         where Δ𝜈 is the linewidth of the laser source. Therefore, it can
          impacts the output signal, as it promotes fading whose                                                 be assumed that 𝜙PN (𝑡) − 𝜙PN (𝑡 − 𝜏𝑔 (𝜙)) ≈ 0, as long as
          relevance depends on 𝑘SBR and 𝛽. The worst case occurs for                                             𝜏𝑔 (𝜙) < 𝜏CT . For realistic lasers, the linewidth can be safely
          𝑘SBR = 0.5, as
                                                                                                                 assumed to be lower than 10 MHz, leading to a coherence time
                     1                        𝜋2                                                                 of 100 ns. As shown in section III, such value is much higher
           𝐼(𝑡) = 𝑅𝐺sig 𝐺OLO 𝑃in                      𝐴IQ 𝑎(𝑡) ⋅ sin((𝜔RF − 𝜔LO )𝑡 +
                     4                     𝑉𝜋 𝑉𝜋,IQ                                                  (28)        than typical values of 𝜏.
                                            𝜃(𝑡)) sin(𝛽).
                                                                                                                 F. Power budget
          In this case, instead of translating into RF phase shifting, optical
                                                                                                                    In order to calculate the power budget of the proposed
          phase shifting only produces pure fading.
             Concerning coherent detection, the most relevant impairment                                         beamforming system, let us first consider Fig. 4, which presents
          of the BPD is the imbalance between its photodiodes. This can                                          a trivial extension of the basic architecture shown in Fig. 1 from
          be modeled in (16) as
                                    𝐼(𝑡) = 𝐼1 (𝑡) − 𝑏BPD 𝐼2 (𝑡)
                          𝑅                                   2
                         = (|𝐺sig 𝐸MZDI (𝑡) + 𝑗𝐺OLO 𝐸OLO (𝑡)|
                          2                                                                          (29)
                                           − 𝑏BPD |𝑗𝐺sig 𝐸MZDI (𝑡)
                                                                           2
                                                    + 𝐺OLO 𝐸OLO (𝑡)| ),
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                6
          one to 𝑁 AEs, again considering a single laser source. Using a                                         PAA receiver of a satellite receives multiple return user beams
          single laser source has the advantage that a single OLO and                                            within a frequency range of 27.5 GHz to 30 GHz and with a
          coherent detector are needed for multiple AEs. The                                                     maximum spectral bandwidth of 1 GHz of bandwidth. The
          disadvantage is that the power of the laser source must be                                             modulation format is typically chosen from the Digital Video
          shared among all AEs. Therefore, for a large number of AEs,                                            Broadcasting - Satellite - Second Generation (DVB-S2)
          one should consider more laser sources, for example by using                                           standard, which supports various amplitude and phase-shift
          multiple beamforming systems such as the one of Fig. 4, and                                            keying (PSK) formats.
          combining their outputs either by electrical or digital means.                                            An ideal CW laser with an average output power of 100 mW
          The laser signal is split by a 1:2 splitter, in which the upper                                        was generated. A 1 GBd quadrature PSK (QPSK) data signal
          output is in turn split by 𝑁 paths, once per AE, and the lower                                         pulse-shaped by a root-raised cosine (RRC) filter with a roll-off
          output is used for OLO generation. The output signal of the 𝑛th                                        factor of 𝛽RRC = 0.5 was generated with length of 29 symbols.
          MZM is given by the generalization of (6) to                                                           The data signal was modulated onto an ideal carrier with a
                                             𝜋                                                                   frequency of 𝑓RF = 30 GHz. Such signal is overdemanding for
                    𝐸𝑛 (𝑡) = 𝐸in (𝑡) ⋅ √𝑘1 𝑘2 𝑣𝑛 (𝑡),                (33)
                                             𝑉𝜋                                                                  the considered scenario, and thus provides a conservative
                                                                                                                 performance evaluation of the system. 𝑁 copies of such RF
          where 𝑘1 and 𝑘2 are the splitting ratios of the 1:2 and 1: 𝑁
          splitters. Concerning the generation of the OLO, (14) becomes                                          signal were produced, individually delayed, and fed to the 𝑁
                                                                                                                 MZMs. A modulation index of 25% relatively to 𝑉𝜋 was
                                               1 𝜋
                  𝐸OLO (𝑡) = √1 − 𝑘1                     𝐸in (𝑡) ⋅ 𝐴IQ exp(𝑗𝜔LO 𝑡).                  (34)        considered. The delay given to each copy was defined
                                               2 𝑉𝜋,IQ
                                                                                                                 according to the following PAA specifications. The steering
          Since the final purpose of the beamforming system is to                                                range of the PAA was limited to |Θ| ≤ 30o , and the distance
          constructively combine the contributions from all AEs, for a                                           between adjacent AEs was of 𝑑 = 0.85𝜆, where 𝜆 is the
          matter of simplicity let us assume ideal TODLs such that                                               wavelength of the input radio signal. The time interval between
          𝐸MZDI,𝑛 (𝑡) = 𝐸𝑛 (𝑡) = 𝐸(𝑡). The low-pass filtered output                                              the radio signals received by two adjacent AEs is given by 𝛿 =
          current amplitude is given by the generalization of (19) to                                            𝑑 sin(Θ)⁄𝑐 , where 𝑐 is the speed of light. Hence, 𝛿 varies
                          𝑅                                               𝜋2                                     between ±𝛿max, where 𝛿max = 14.17 ps, and thus the copy of
              𝐼 = −𝑁 𝐺sig 𝐺OLO 𝑃in √𝑘1 (1 − 𝑘1 )𝑘2 𝑘3                              𝐴IQ 𝐴sig,         (35)
                          4                                             𝑉𝜋 𝑉𝜋,IQ                                 the RF signal fed to the 𝑛th AE is delayed by (𝑛 − 1)𝛿, where
          where 𝑘3 is the coupling ratio of the 𝑁: 1 combiner. As 𝑘2 =                                           𝑛 = 1, … , 𝑁.
                       ⌈log2 𝑁⌉                                                                                     Regarding the value of 𝜏, even though the operation principle
          𝑘 3 = 2−         2      , for ⌈log 2 𝑁⌉ = log 2 𝑁 (35) becomes                                         of the TODL is analogous to the one presented in [13], in the
                        𝑅                               𝜋2                                                       present case it is not necessary to define 𝜏 according to equation
                   𝐼 = − 𝐺sig 𝐺OLO 𝑃in √𝑘1 (1 − 𝑘1 )         𝐴 𝐴 .                                   (36)
                        4                            𝑉𝜋 𝑉𝜋,IQ IQ sig                                             (7) of [13]. As only one of the sidebands is detected, there is
                                                                                                                 only the need to center one of them at a maximum of the
          The average output electrical power is given by                                                        TODL’s amplitude response. Therefore, the minimum required
                    1                                                                                            value of 𝜏 only depends on 𝛿max and on 𝑁, and is given by
             𝑃avg = 𝐼2 ⋅ 𝑅Ω
                    2                                                                                                                𝜏 = 𝛿max ∙ (𝑁 − 1).                      (39)
                                                          2                                          (37)
                    𝑅2              2    2    3
                                                   𝜋2          2
                  =    𝑘1 (1 − 𝑘1 )𝐺sig 𝐺OLO 𝑅Ω (        ) . 𝑃in 𝑃IQ 𝑃sig ,                                      The phase shifters 𝛽1 , … , 𝛽𝑁 , were adjusted such that all delayed
                     8                           𝑉𝜋 𝑉𝜋IQ
                                                                                                                 optical signals were constructively added, thereby maximizing
          where 𝑅Ω is the output load, 𝑃IQ is the electrical power fed to                                        the output electrical power . The phase shifters of the TODL
          the IQM, and 𝑃sig is the electrical power fed to each MZM.                                             𝛾1 , … , 𝛾𝑁 were set by default to an ideal value, enabling a
          Therefore, the power budget of the system is given by                                                  perfect alignment of the input optical signals with the frequency
                                                                                                                 response of the MZDIs.
                                                                                            2
                  𝑃out       𝑅2              2    2
                                                              𝜋2                                     (38)           Concerning the generation of the OLO, the IQM was fed with
                2          =    𝑘1 (1 − 𝑘1 )𝐺sig 𝐺OLO 𝑅Ω3 (         ) .
              𝑃in 𝑃IQ 𝑃sig   8                              𝑉𝜋 𝑉𝜋IQ
          This result shows that for a fixed set of input powers the output
          power is constant and independent of the number of AEs. It also
          shows that the best strategy to increase the output power is
          increasing the gain of the optical amplifiers associated with 𝐺sig
          and 𝐺OLO , and also increasing 𝑃in .
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                  7
                                                                                                                        OIP3 [dBm]
          detected signal was match-filtered by an RRC filter. The                                                                   50
          filtered signal was finally demodulated by a typical QPSK
                                                                                                                                     45
          demodulator, however without any carrier phase recovery
          technique nor algorithm.
                                                                                                                                     40
             Other than laser phase noise, no noise was considered at any                                                             10        20    30             40                                       50
          point of the system, meaning that signal degradation only                                                                         IQ modulation index[%]
          stemmed from phase noise, third-order intermodulation                                                    Fig. 6. OIP3 for different modulation indexes of the IQM.
          distortion (IMD3), and non-ideal frequency response of the
          TODLs. The error vector magnitude (EVM) [17] was the only                                                                  180
          figure of merit used, as for slightly distorted signals the bit error                                                                                                               4
          rate is consistently zero.                                                                                                 160                                                         6
                                                                                                                                                     2
                                                                                                                                     140                                                                   10
                                                                                                                                                                       4
          A. Laser Phase Noise                                                                                                                                                               8
             As observed in (31), in practice, laser phase noise is not                                                              120
                                                                                                                                                                                                      12
          ideally cancelled, and therefore its impact should be assessed.
                                                                                                                                                                              6
                                                                                                                     [degrees]
          A system with a single AE was considered. In order to consider                                                             100
                                                                                                                                                                                            10
          the worst-case scenario, the TODL was set to its maximum
                                                                                                                                                2
                                                                                                                                      80
          delay, i.e., 𝜏𝑔 (𝜙) = 𝜏, meaning that degradation was caused
                                                                                                                                                                                                         14
                                                                                                                                                                                       8
          only by non-ideal phase noise cancellation.                                                                                 60
                                                                                                                                                                                                         12
             As shown in Fig. 5 the EVM increases both with Δ𝜈 and 𝜏.                                                                                                              6
                                                                                                                                      40
          As expected, for 𝜏 = 0 the EVM is at its lowest, independently                                                                                                                                   10
          of the linewidth, as phase noise cancellation is ideal. Assuming                                                                               2                        4                        8
                                                                                                                                      20                                                                   6
          that the FSOLO is delayed on purpose by 𝜏/2, the maximum                                                                                                                                         4
                                                                                                                                                                                  2                        2
          absolute delay between the FSOLO and any delayed optical                                                                    0
                                                                                                                                       0    2            4          6              8             10          12
          signal is of 𝜏/2, instead of 𝜏. This means that, in practice, the                                                                                       /TRF
          degradation quantified in Fig. 5 can be achieved for a value of                                             Fig. 7. EVM (%) as a function of the maximum time delay of the
          𝜏 that is twice as indicated. As an example, let us consider a                                              TODL, 𝜏, and the phase shift 𝜙, responsible for tuning the time
          realistic linewidth of 1 MHz. For such a value, the EVM                                                     delay between 0 (0º) and 𝜏 (180º), for a system comprising a single
          estimated for an indicated time delay of 𝜏 = 6TRF is twice the                                              AE. The dots represent the angles 𝜙 which result in the highest
          minimum EVM. Given the considered maximum value of                                                          EVM for the corresponding values of τ.
          𝛿max = 14.17 ps, the maximum number of AEs supported by
          such value is of (𝑁 − 1)𝛿max = 2𝜏, which yields 𝑁 = 29 AEs.                                                                18
                                                                                                                                           RRC=0
          B. Intermodulation distortion                                                                                              16
                                                                                                                                           RRC=0.25
             In order to quantify the impact of modulation-imposed                                                                   14    RRC=0.5
          distortion, the IMD3 was characterized. Without loss of
          generality, a system with a single AE was once again                                                                       12    RRC=0.75
          considered. The MZM was modulated by two purely sinusoidal                                                                       RRC=1
                                                                                                                         EVM [%]
                                                                                                                                     10
          signals with frequencies of 30 and 30.05 GHz. Different
          modulation indexes were applied to the IQM. In order to avoid                                                               8
          overdrive, the maximum modulation index was of 50%.
             The input third-order intercept point (IIP3), was measured at                                                            6
          43.7 dBm, regardless of the modulation index of the IQM. As                                                                 4
          discussed in section II, the non-ideal generation of the OLO
          only produces undesired tones in the output electrical signal,                                                              2
          which can be filtered out after bandpass filtering. Hence, the
                                                                                                                                      0
          generation of the OLO does not impact the IMD3 of the system.                                                                0    2           4          6              8           10            12
          The measured output third-order intercept points (OIP3) are                                                                                            /TRF
          shown in Fig. 6. Increasing the modulation index of the IQM                                                 Fig. 8. EVM as a function of 𝜏 for a system comprising a single AE
          from 10% to 50% results in an increase of the OIP3 of 12 dB.                                                and 𝜙 = 90o. Different roll-off factors of the pulse-shaping filter
                                                                                                                      are considered. The insets correspond to the constellation diagrams
          C. Performance assessment for a single antenna element                                                      obtained for 𝜏 = 1.5 ∙ TRF and 𝜏 = 10 ∙ TRF , for 𝛽RRC = 0.5.
            The TODL has an ideal behavior only when configured for
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                8
6 1
                                                                                                                                                                                                  9
                                                                                                                                                   3
                                                                         14
22 20 18 16
                                                                                                                                                                                                         11
                                                                                                                                                                    5
                                         6
                                                                                                                                                                                          8
                                                                       12
                     5
                                                                                                                                                                                                  10
                                                                                                                                0.8
                                                                                                                                                                              7
                                                                                                                                                            4
                                                                                                                                      2
                                                                    10
                                                                                                                                                                        6
                                   4
                                                                                                                                                                                                          12
                     4                                                                                                                                                                        9
                                                                  8
                                                                                                                                0.6
                                                                                                                                                                                                          11
              [%]
max
                                                                                                                                                                                          8
                                                                                                12 10
                                                              6
                     3
                                                                                                                    /
                RF
                                                                                                                                                                                                         10
                     2
                                                                                                                                                                                     7
             /T
                                                                                                                                                                4
                                                                                                                                0.4                                              6
                                                                                                                                          2
                                                         4
                                                                                                   8                                                                                                        9
                     2
                                                                                                                                                                3                    5                     8
                                                                                                     6
                                                                                                                                0.2                                              4                      67
                     1                               2                                               4                                                                                                      5
                                                                                                                                                            2                        3                    4
                                                                                                                                                                                2                            3
                                                                          2                                                                                                                              2
                                                                                                                                 0
                     0                                                                                                                    10           15            20                  25            30
                      0                20  40           60          80
                                                                                                                                                                Number of AEs
                                        || [degrees]
           Fig. 9. EVM (%) as a function of the detuning between the frequency
                                                                                                                  Fig. 10. EVM (%) as a function of the number of AEs and of the
           response of the TODL and the input optical signal, ∆𝛾, and τ, for a
                                                                                                                  normalized time interval between the radio signals received by
           system comprising a single AE and a for 𝜙 = 90o.
                                                                                                                  two adjacent AEs , 𝛿/𝛿max .
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                                9
          gain without affecting the performance of the system.                                                  QPSK input signal centered at 30 GHz, and a unidimensional
            Fig. 10 shows that the EVM increases with 𝛿 and with the                                             PAA with a steering range of |Θ| ≤ 30o and with a distance
          number of AEs. For 𝛿 = 0 all MZDIs have 𝜙 = 0. Hence, the                                              between adjacent AEs of 0.85𝜆. It was confirmed that the
          delayed optical signals are not impaired by the TODLs, and                                             degradation imposed by non-ideal phase noise cancelation is
          consequently the output optical signal has the lowest EVM.                                             negligible for realistic laser linewidths. Concerning the
          However, as 𝛿 increases, so does the delay added by each                                               distortion introduced by the non-ideal frequency response of
          consecutive TODL, given by                                                                             each TODL, it was observed that it is caused by narrowband
                                                                                                                 filtering and not by group delay distortion. For multiple AEs the
                                         𝜏𝑛 = 𝛿 ∙ (𝑛 − 1),                                          (40)         distortion induced by narrowband filtering increases with the
                                                                                                                 number of AEs, as so does the value of 𝜏. The distortion
          where 𝑛 = 1, … , 𝑁. The bandwidth of the amplitude response                                            observed in the output signal is the average of the distortion
          of the 𝑛th MZDI depends on 𝜏𝑛 , reaching a minimum for 𝜏𝑛 =                                            imposed by all TODLs. This means that the worst-case
          𝜏/2. Consequently, the delayed optical signals are subject to                                          performance of the system can be conservatively estimated by
          narrowband filtering whose bandwidth depends on 𝜏𝑛 . The                                               considering a system with a single AE, the same value of 𝜏, and
                                                                                                                                              𝜋
          output optical signal is the addition of delayed optical signals                                       the worst case value of 𝜙 = . Reduced distortion was observed
                                                                                                                                              2
          with different levels of distortion. Hence, the distortion                                             for a system comprising up to 32 AEs, which means that for the
          observed in the output optical signal is the average of the                                            considered simulation scenario, the system is suitable for a
          distortions observed in the 𝑁 delayed optical signals. The                                             bidimensional PAA with 16 × 16 AEs.
          strongest distortion is observed when the average value of 𝜏𝑛 is                                          Within the scope of project BEACON [18], future work will
          closest to its worst-case value of 𝜏/2. This is generally the case                                     take two directions. On a theoretical side, as the presented work
          when 𝛿 = 𝛿max , as in this case 𝜏𝑛 varies from 0 up to 𝜏. Fig. 11                                      only addresses the impact of signal distortion, it is important to
          shows the variation of the EVM as a function of the number of                                          assess the relevance of noise and beam squinting, also in a
          AEs. The dashed line corresponds to the worst performance                                              scenario based on [16]. On a practical side, the beamformer will
          obtained in a system comprising a single AE with 𝜙 = 90o and                                           be demonstrated resorting to silicon photonic integrated
          with a value of 𝜏 given by (39). Such performance is                                                   circuits.
          consistently slightly worse than the one of a system comprising
          multiple AEs and having the same value of 𝜏. Consequently, in                                                                                 REFERENCES
          order to estimate the worst-case performance of a complete
                                                                                                                 [1]          R. C. Hansen, Phased Array Antennas. Wiley, 2009.
          system comprising multiple AEs, one can simply consider a                                              [2]          M. Sotom, B. Benazet, A. Le Kernec, and M. Maignan, “Microwave
          system with a single AE, with 𝜙 = 90o , and with the same                                                           photonic technologies for flexible satellite telecom payloads,”
          value of 𝜏. The insets of Fig. 11 show that the proposed                                                            Optical Communication, 2009. ECOC ’09. 35th European
                                                                                                                              Conference on. pp. 1–4, 2009.
          photonic beamforming system does not introduce severe                                                  [3]          W. Roh, J.-Y. Seol, J. Park, B. Lee, J. Lee, Y. Kim, J. Cho, K.
          distortion, even for as many as 32 AEs. Consequently, the                                                           Cheun, and F. Aryanfar, “Millimeter-wave beamforming as an
          proposed system is suitable for a unidimensional PAA with 32                                                        enabling technology for 5G cellular communications: theoretical
                                                                                                                              feasibility and prototype results,” Communications Magazine, IEEE,
          AEs, or alternatively, for a bidimensional PAA with 16 × 16                                                         vol. 52, no. 2. pp. 106–113, 2014.
          AEs.                                                                                                   [4]          J. Yao, “Microwave Photonics,” Lightwave Technology, Journal of,
                                                                                                                              vol. 27, no. 3. pp. 314–335, 2009.
                                                                                                                 [5]          J. Capmany and D. Novak, “Microwave photonics combines two
                                            IV. CONCLUSIONS                                                                   worlds,” Nat Phot., vol. 1, no. 6, pp. 319–330, Jun. 2007.
             In this paper a novel photonic TTD beamformer for PAA                                               [6]          W. Ng, A. A. Walston, G. L. Tangonan, J. J. Lee, I. L. Newberg,
                                                                                                                              and N. Bernstein, “The first demonstration of an optically steered
          receiver was presented. The proposed system is based on a set                                                       microwave phased array antenna using true-time-delay,” Lightwave
          of TODLs and on self-heterodyne coherent detection. Each                                                            Technology, Journal of, vol. 9, no. 9. pp. 1124–1131, 1991.
          TODL is based on a MZDI with a time delay of 𝜏, in which the                                           [7]          J. J. Lee, R. Y. Loo, S. Livingston, V. I. Jones, J. B. Lewis, H. Yen,
                                                                                                                              G. L. Tangonan, and M. Wechsberg, “Photonic wideband array
          delay can be varied between 0 and 𝜏 by defining the power ratio                                                     antennas,” Antennas and Propagation, IEEE Transactions on, vol.
          between the arms of the MZDI. Self-heterodyne detection                                                             43, no. 9. pp. 966–982, 1995.
          enables maximum sensitivity, laser phase noise cancellation,                                           [8]          R. Soref, “Optical dispersion technique for time-delay beam
                                                                                                                              steering,” Appl. Opt., vol. 31, no. 35, pp. 7395–7397, 1992.
          photonic RF phase shifting and photonic RF frequency                                                   [9]          H. Zmuda, R. A. Soref, P. Payson, S. Johns, and E. N. Toughlian,
          downconversion.                                                                                                     “Photonic beamformer for phased array antennas using a fiber
             A theoretical analysis of the system was derived to                                                              grating prism,” Photonics Technology Letters, IEEE, vol. 9, no. 2.
                                                                                                                              pp. 241–243, 1997.
          demonstrate its operation. The analytical results show that the
                                                                                                                 [10]         B. Ortega, J. L. Cruz, J. Capmany, M. V Andres, and D. Pastor,
          system is robust to realistic impairments such as third-order                                                       “Variable delay line for phased-array antenna based on a chirped
          harmonic generation, finite ER and bias deviations of the                                                           fiber grating,” Microwave Theory and Techniques, IEEE
          MZMs, imbalance between the photodiodes of the BPD and                                                              Transactions on, vol. 48, no. 8. pp. 1352–1360, 2000.
                                                                                                                 [11]         A. Meijerink, C. G. H. Roeloffzen, R. Meijerink, L. Zhuang, D. A.
          non-ideal phase noise cancellation. The power budget was                                                            I. Marpaung, M. J. Bentum, M. Burla, J. Verpoorte, P. Jorna, A.
          derived, showing that for a fixed set of input powers the output                                                    Hulzinga, and W. van Etten, “Novel Ring Resonator-Based
          power is independent of the number of AEs.                                                                          Integrated Photonic Beamformer for Broadband Phased Array
                                                                                                                              Receive Antennas—Part I: Design and Performance
             The system was numerically assessed considering a 1 GBd                                                          Analysis,” Lightwave Technology, Journal of, vol. 28, no. 1. pp. 3–
  0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.
This article has been accepted for publication in a future issue of this journal, but has not been fully edited. Content may change prior to final publication. Citation information: DOI 10.1109/JLT.2016.2620146, Journal of
                                                                                                       Lightwave Technology
          > REPLACE THIS LINE WITH YOUR PAPER IDENTIFICATION NUMBER (DOUBLE-CLICK HERE TO EDIT) <                                                                                                              10
                       18, 2010.                                                                                              With Separate Tuning of the Optical Carrier,” IEEE Photonics
          [12]         L. Zhuang, C. G. H. Roeloffzen, A. Meijerink, M. Burla, D. A. I.                                       Technology Letters, vol. 21, no. 22. pp. 1686–1688, 2009.
                       Marpaung, A. Leinse, M. Hoekman, R. G. Heideman, and W. van                               [16]         P. Inigo, O. Vidal, B. Roy, E. Alberty, N. Metzger, D. Galinier, J.
                       Etten, “Novel Ring Resonator-Based Integrated Photonic                                                 Anzalchi, G. Huggins, and S. Stirland, “Review of terabit/s satellite,
                       Beamformer for Broadband Phased Array Receive Antennas—Part                                            the next generation of HTS systems,” 2014 7th Advanced Satellite
                       II: Experimental Prototype,” J. Light. Technol., vol. 28, no. 1, pp.                                   Multimedia Systems Conference and the 13th Signal Processing for
                       19–31, 2010.                                                                                           Space Communications Workshop (ASMS/SPSC). pp. 318–322,
          [13]         M. V Drummond, P. P. Monteiro, and R. N. Nogueira, “Photonic                                           2014.
                       True-Time Delay Beamforming Based on Polarization-Domain                                  [17]         R. A. Shafik, S. Rahman, and R. Islam, “On the Extended
                       Interferometers,” Lightwave Technology, Journal of, vol. 28, no. 17.                                   Relationships Among EVM, BER and SNR as Performance
                       pp. 2492–2498, 2010.                                                                                   Metrics,” Electrical and Computer Engineering, 2006. ICECE ’06.
          [14]         M. V Drummond, R. N. Nogueira, P. P. Monteiro, M. A. Violas, C.                                        International Conference on. pp. 408–411, 2006.
                       Sterner, and P.-Y. Fonjallaz, “Tunable Optical Dispersion                                 [18]         “BEACON, photonics for telecom satellites.” [Online]. Available:
                       Compensator Based on Power Splitting Between Two Dispersive                                            http://www.space-beacon.eu/. [Accessed: 31-Aug-2016].
                       Media,” Lightwave Technology, Journal of, vol. 28, no. 8. pp. 1164–
                       1175, 2010.
          [15]         P. A. Morton and J. B. Khurgin, “Microwave Photonic Delay Line
0733-8724 (c) 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.