arxiv: v1 [math-ph] 11 Dec 2016

Transcription

1 The tructure of uadratic Gau um in Talbot effect Carlo R. Fernández-Poua Deartment of Communication Engineering, Univeridad Miguel Hernández, Av. Univeridad /n, E03203 Elche, Sain Dated: November 29, 2016 arxiv: v1 [math-h] 11 Dec 2016 The field diffracted from a one-dimenional, coherently illuminated eriodic tructure at fractional Talbot ditance can be decribed a a coherent um of hifted unit cell weighted by a et of hae given by uadratic Gau um. We reort on the comutation of thee um by ue of the roertie of a recently introduced integer, which i contructed here directly from the two corime number that that define the fractional Talbot lane. Uing integer, the comutation i reduced, u to a global hae, to the trivial comletion of the exonential of the uare of a um. In addition, it i hown that the Gau um can be reduced to two cae, deending only on the arity of integer. Exlicit and imler exreion for the two form of integer are alo rovided. The Gau um are reented a a Dicrete Fourier Tranform air between eriodic euence of length howing erfect eriodic autocorrelation. The relationhi with one-dimenional multilevel hae tructure i exemlified by the tudy of Talbot array illuminator. Thee reult rereent a imle mean for the deign and analyi of ytem emloying the fractional Talbot effect. Keyword: Talbot effect, uadratic Gau um, erfect eriodic autocorrelation, Talbot array illuminator. I. INTRODUCTION When a diffraction grating i illuminated by a coherent lane wavefront, a erie of elf-image can be oberved at certain ditance beyond the object, at integer multile of a fundamental Talbot ditance [1], [2]. A related henomenon occur at fractional value / of thi ditance, where the field i a coherent ueroition of hifted and weighted coie of the grating unit cell [3]. Thee henomena are coneuence of the interference between diffraction order, which acuire a uadratic hae under Frenel roagation, and have alo been oberved in other domain, uch a the angular ectrum [4] and, through the ace-time duality [5], [6], in time [7] and in the otical ectrum [8], [9]. A imle diffractive effect, they have found wideread alication [10], [11]. Much of thi work wa triggered by the ublication, twenty year ago, of a eminal aer by Berry and Klein [12], who howed the relationhi of the weighting factor at fractional Talbot ditance with the claical uadratic Gau um of number theory. The olution to thee Gau um wa reented in [12] by conidering three earated cae, relying on the reult of a reviou invetigation [13]. One of thee three olution contained an error which wa corrected later by Matutami and Ônihi [14]. More recently [15], it ha been hown that a roer rewording of the original threefold form of the Gau um lead to a comact decrition in term of the roduct of two hae, one that i contant and only deend on integer and and a econd that i uadratic and roortional to a new integer. The obervation in [15] can be tated a the fact that the Dicrete Fourier Tranform DFT of a uadratic hae euence i itelf a uadratic hae euence. In addition to it intrinic c.oua@umh.e beauty, thi correondence find immediate alication a a imle rule to imlement the hae in hyical ytem emloying fractional Talbot effect. However, the characterization of the integer that arie in the Gau um wa obtained in [15] by direct comarion, relying on the original comutation in [13] and [14]. In thi aer, we reviit the comutation of thee Gau um from a new erective. Uing tandard techniue from number theory, integer i contructed in an abtract way from the two corime integer and that define the fractional Talbot effect. The roertie hown by thi integer allow for an almot trivial derivation of the main roerty of the Gau um, namely it uadratic deendence. It i alo hown that both and the reult of the Gau um can be exreed only in two cae aociated to the arity of integer and, in addition, that both member of the DFT air that decribe the Gau um how erfect eriodic autocorrelation [16], [17]. Thee reult not only imlify the decrition of elf-imaging henomena, but alo ermit a comact analyi of ytem emloying the fractional Talbot effect. A an examle, the relationhi with multilevel hae tructure i exemlified by the detailed tudy of one-dimenional Talbot array illuminator TAI [18], [19], [20], [21], [22], [23]. In our reentation we will emloy the decrition of Talbot effect in the time domain [7], where temorally modulated eriodic wave ubtitute for diffraction grating and dierion for diffraction. Although our reult can be extended to different domain a ointed out before, the ue of the temoral formalim i referred here becaue it conform to the tandard definition of dicretetime tranform of ignal theory, which i the natural language uon which the main reult are reented. The lan of the aer i the following. In Section II we review the baic reult of [12] about the tructure of Talbot elf-image of eriodic object. In Section III, we contruct the integer and derive everal roertie, and

2 2 in Section IV we comute the Gau um of Talbot effect and derive their autocorrelation roertie. Section V i devoted to the analyi of TAI and Section VI reent our concluion. The aer i comleted with three aendice which reent ome reult and definition from number theory and two lengthy comutation. II. THE TEMPORAL TALBOT EFFECT Let u conider a modulated otical wave, whoe electric field i given by Et Et exjω 0 t, with ω 0 the freuency of the otical carrier and et the otical comlex enveloe or baeband rereentation of the wave. The enveloe i aumed eriodic with fundamental eriod T, and rereented a the reetition of function wt: Et + n δ t nt wt, 1 where denote convolution. Here we aume that wt i contained in a eriod, 0 wt < T, and o it decribe the unit cell of the eriodic train. The ectrum of the otical enveloe Et i contained in a et of euallyaced angular freuencie, ω n 2πn/T with n integer. Thi field feed a linear delay line with lowet-order dierion. Without lo of generality, the tranfer function that decribe the linear roagation of the enveloe Et can be aumed uadratic, Hω ex jφω 2 /2, where φ i the dierion coefficient, eual to the derivative of the grou delay with angular freuency at the carrier value, ω 0. The temoral Talbot effect manifet itelf a the recovery of the initial train, or of a eriodic tructure directly related to the original train, at certain value of dierion φ characterized by the eriod T and two oitive and mutually corime integer and. The defining euation i: 2π φ T 2. 2 At thee value, Hω n ex jπσ φ n2, where σ φ i the ign of φ. The outut enveloe E t can be comuted in the ectral domain by ue of the Poion formula in 1 and the value of Hω n, E t 1 T + n e jπσ φ n2 e j2πn t T wt. 3 Note that the value of the um for σ φ 1 i the comlex conjugated to the um for σ φ 1, ince the change in the ign of the exonential can be comenated with an inverion n n. Thi firt art of the euation can be decomoed by changing the um over n to a um over the air of integer l 0,..., 1 and m,..., +, with n l + m. The um over m can be erformed uing again the Poion formula, and we get [12]: E T t δ t e e jσ φξ n δ n0 + k t n T δ t kt wt, 4 where e x rereent the arity of integer x, o that e x 0 when x i even and e x 1 when x i odd. In 4, e thu rereent the arity of the roduct of integer. The hae in 4 are given by the following Gau um: e jσ φξ n e πjσ φ m2 1 m e 2πj nm e πjσ φ 1+em2 e 2πj nm. 5 In the final art of the euation we have ued that 1 m 1 em2 σ φ. Thi econd form of the Gau um will be the tarting oint of our comutation in Section IV. According to 4, the outut i comoed of a reetitive tructure of the ame eriod T a the original train, exlicitly recognized in the econd term of the multile convolution. Thi train i hifted by half a eriod if i odd, a hown by the firt delta term. The reetitive tructure that define the outut train i decribed by the two lat term in 4, and i comoed of the coherent um of relica of the inut unit cell, wt, mutually hifted by T/ and weighted by factor exjσ φ ξ n /. We oint out that, exect in the cae that wt i contained in an interval of length T/, the lat line in 4 doe not decribe the unit cell of the outut, ince thi exreion may extend beyond an interval of length T. When 1, the outut train i imilar to the original. The reulting wave i referred to a an integer Talbot image of order or index, or imly a elf-image of order of the inut train. If, in addition, i odd, it i uually referred to a a hifted integer Talbot image. In general, the reulting train for > 1 are called fractional Talbot image of order /, and hifted when the roduct i odd. To roceed with our analyi we devote the following ection to the characterization of hae exjξ n in term of an integer contructed from and. III. A PARITY-DEPENDENT MODULAR INVERSE The main reult concerning the roertie of integer i tated in the following theorem: Theorem 1. Given two corime and oitive integer, and, there exit a uniue integer uch that it verifie the following three roertie:

3 3 a lie in the range 1 2 1, b i a olution of the modular euation: 1 + e mod 2, 6 c ha ooite arity to. In addition, uch a integer alo verifie: d i corime with. The theorem will be roved in four te. Firt, we how the exitence of integer verifying a and b by the exlicit contruction of olution to 6. Second, we will how that thee olution alo verify c. Third, we how that any other olution of 6 i either out of the range a or doe not verify c. Finally, we deduce that and mut be corime. For the definition of the ymbol ued in the roof we refer to Aendix A. Exitence. For even, E. 6 reduce to 1 mod 2, and ince integer and are mutually rime, o are and 2. Then the euation can olved a: even, 7 2 where [1/a] n i the invere of a mod n. In articular, the invere of mod 2 lie in the range When i odd, integer 2 and are corime, and a olution to 6 i given by: 2 odd. 8 2 Indeed, in thi cae the roduct i: 2 1 mod 1 + a, 9 2 for a certain integer a. But thi a mut be odd, ince in 8 i even and i odd, and therefore 9 can be written a: 1 + mod 2, 10 thu olving euation 6. Now, ince the invere of 2 mod i contained in the interval 1 [1/2] 1, integer lie in the range tated in the theorem. Parity. A for roerty c, the tatement i trivial for the olution found for odd, ince the exlicit contruction 8 how that in thi cae i even. For even we firt note that mut be odd. Now, the olution of 6 i of the form: 1 + 2b 11 for certain integer b. The right hand ide of thi euation i odd, and therefore mut be odd. Note that the exlicit olution 7 and 8 together with roerty c oint out that i an modular invere of whoe concrete form deend on the arity of. TABLE I. Value of integer in the range a a function of and Uniuene. Integer 7 and 8 are the uniue olution to 6 in the range with ooite arity to. In fact, given two olution of 6 for given and, namely given by 7 or 8, and a different olution, their difference verifie 0 mod 2. Here we need to analyze earately the cae. For even, and ince and 2 are corime, thi imlie that 0 mod 2 and o the uniuene of in the range follow. The ame argument alie for odd if we further aume that i alo odd, ince again and 2 are corime. The only remaining cae i odd and even. The euation above, 0 mod 2, only imlie that 0 mod, and o both 8 and + are olution to 6, a it canbe eaily checked. In articular, may lie in the range For examle, for 2 and 3, both 2 and 5 are olution of 6 in the range 1, 5. However, integer a given by 8 i elected by the fact that it ha the ooite arity to, contrary to +. Thi comlete the roof. Corime. Finally, to how that integer i indeed corime with, let u uoe that and hare a common factor, ay c. We will how that c ±1. Firt, note that c i alo a factor or the roduct. Then, dividing 6 by thi factor we get: c 1 c + e c mod c Now, ince the left hand ide of thi euation i an integer and c i a factor of, then 1/c hould be integer, and thu c ±1. The roof of Theorem 1 i comlete. We oint out that the exreion for derived in [14] and ued in [15] are different from the one given by 7 and 8. The euivalence i hown in Aendix B. The comutation of the value of, a a function of and, can roceed from the ytematic olution of 7 and 8 or from an euivalent exreion uch a thoe in Aendix B. A table for the lowet value of and wa reented in [15] and reroduced in Table I ee, however, note [24]. The calculu of

4 4 i facilitated by the obervation that, regardle the arity of a given, 6 imlie that integer for given corime number and i the ame that for + 2 and. Thi roerty reflect the exact eriodicity, without half-interval hift, of the temoral Talbot effect. In the cae of odd, it alo follow from 6 that the integer for the air, and for the air + and, coincide. Unfortunately, the reult doe not hold for even. Secific olution can be found for everal erie of integer, a hown by the following reult. Prooition 2. a for any : 1 mod e. b for any : ± 1 mod 2 ± 1. c for any n, : ±1 + 2n ±2n mod 2. Proof. Firt we note that in thee three erie integer and are indeed corime, ince in all cae they are defined by a relation of the form a + b ±1 for certain integer a and b, and thu any common factor between and mut be eual to one following a imilar reaoning to that in 12. Second, it i a imle obervation that the arity of i the ooite to that of in all cae. Moreover, all rooed value of lie in the range 1 2 1, excet in cae c for which the reidue mod 2 ha been taken. Then, it uffice to check that the value of verify 6 by direct ubtitution. The reult i trivial for a, wherea for b and c it can be deduced from the obervation that 6 i euivalent to any of the following euation: ± 1 2 mod In articular, 13 define the elf-dual integer ± 1 for which the value of and coincide. Thi comlete the roof. The calculu of i alo imlified by the following reult, that relate it value for comlementary Talbot line of order / and /. Eentially it i a um rule that alie to air of integer in each row of Table I: Prooition 3. If i the integer correonding to and, with <, then the integer correonding to and, verifie; a 2 for odd, b for even and 1 1, and c 3 for even and Proof. Firt notice that, if i corime with, then i alo corime with. Moreover, integer a tated above ha the ame arity a, o condition c in Theorem 1 i fulfilled. A for condition a in that Theorem, for odd + 2, and therefore if i in the range o i. For even, both b and c above give in the range Note that cannot coincide with ince they mut be corime. To comlete the roof it uffice to check 6. Direct comutation give: + [n ], 14 with n 1, 2, 3 according to the tatement of the rooition. When even,, and n are odd, and thu mod 2. For odd, and n are even, and again mod 2. Therefore, ince by hyothei 6 i verified by, it i alo verified by. The roof i comlete. IV. THE GAUSS SUMS Integer contructed in the reviou ection rovide a imle route to evaluate the Gau um 4. Our comutation i baed on the obervation that 5 involve the exreion: ex πjσ φ 1 + e m 2 15 which can be conidered a a -euence ince, for even and odd, i eriodic in the running index m with eriod. We recall that, given a euence of comlex number, x n with n 0,..., 1, the Dicrete Fourier Tranform DFT, denoted by F, i defined a the -euence X m given by: 1 X m F{x n } m x n e 2πjnm/. 16 n0 The IDFT i the invere tranform to 16, defined a: x n F 1 {X m } n 1 1 We tate the Gau um 5 a a DFT air: X m e 2πjnm/. 17 Prooition 4. For any, corime and the integer contructed in Theorem 1, the following DFT air hold: x n e jσ φξ n e jσ φξ 0 ex πjσ φ n2 DF T X m [ ] ex jπσ φ 1 + e m 2, 18 with: e jξ0 e j π 4 1 e jξ0 e j π 4 e j π 4 1 odd, e j π 4 even, 19 and a b i the Jacobi ymbol of an arbitrary integer a and an odd and oitive integer b.

5 5 For the definition and roertie of Jacobi ymbol we refer to Aendix A. We alo oint out that the reent reult reduce the number of cae to two, deending only on the arity of, intead to the conventional three cae tudied in Ref. [12], [13], [14] and [15]. Proof. We et σ φ 1, the general cae i conidered at the end of the roof. We have to how that: e jξn e jξ0 e πj n2 1 1 e πj 1+em2 e 2πj nm. 20 Firt we jutify the change of variable m k, with k 0,..., 1. Since and are corime, the value of the roduct k an the reidue of, i. e., the roduct k, with k 0,..., 1, an the et {0,..., 1} mod. The change m k can thu be erformed in 20 ince the contribution to the um of a certain m i invariant under hift of the form m m +, a it i immediate to how. With thi change 20 read: 1 1 k0 e πj 1+ek2 e 2πj nk. 21 We ue now 6 in the firt exonential of 21: e πj 1+ek2 e πj 1+e2 k 2 e πj k2 1 k e πj k2, 22 where we have ued that the roduct i even. With thee imlification, the exonent in 21 can be comleted to the uare of a um, 1 1 e πj n2 k0 and the roof i comlete if we how that: e jξ0 1 1 e πj k n2, 23 e πj m2. 24 The fact that the um in the right hand ide of 24 i indeed a hae can be hown by changing variable in the double um defining it modulu uared, following a imilar argument to that of an euivalent comutation in [12]. The exlicit evaluation of the um leading to 19 can be comleted by tandard techniue of the ame tye a thoe ued in [13], a i hown in Aendix C. We oberve that the uadratic hae exjξ n are eriodic and ymmetric, exjξ n exjξ n+ exjξ n exjξ n. 25 Now, uing the general relationhi for the DFT of the comlex-conjugated euence, F{x n} m F{x n } m, 26 it follow from the ymmetry roerty that the DFT of the comlex conjugated euence, ex jξ n i the comlex conjugated DFT, a tated in 18. The roof i comlete. We now comute the eriodic autocorrelation of the -euence x k exjξ k, 1 R ξ n ex jξ k exjξ k+n. 27 k0 Uing 18, thi give a um of the form: 1 R ξ n e jπ n2 e j2π kn, 28 k0 which can be evaluated in term of Kronecker delta uing that, for any integer a: 1 e j2πak/ δ a,0 mod. 29 k0 We oberve from 28 that the um i nonzero for integer n which are olution of n 0 mod. The olution are n 0 mod, ince and are corime, o: R ξ n δ n,0 mod. 30 The hae euence x k how erfect eriodic autocorrelation [16], and therefore it ectrum, a given by 18, ha magnitude. In fact, it can be hown that both -euence, x n and X m in 18, have erfect eriodic autocorrelation, ince both are examle of Chu contruction [17]. In that aer it wa hown that, given N and M corime integer, the M-euence ex±jπnm 2 /M M even ex±jπnmm + 1/M M odd, 31 have erfect eriodic autocorrelation. For even, the fact that x k exjξ k belong to thi et i immediate. For odd it follow after the hift in the running index k k + + 1/2, ince the reulting euence i in the form 31 u to a global factor. A for the Fourier-tranformed euence X m, it i already in form 31 for even. For odd, the euence read; X m e jπ+1m2 / e jπm2 / 1 m. 32 Uing again the hift m m + + 1/2 we are led, u a global hae and the factor, to: [ ex jπ m m ] 1 m ex [ jπ ] m m + 1, 33 which i the exreion in 31.

6 6 V. TALBOT ARRAY ILLUMINATORS The reviou reult rovide a imle decrition of the weighted um of unit cell that define the fractional Talbot effect. Uing 20 in 4, the outut field in a Talbot lane with order / i decribed a the reetition of the ignal: w t ejσ φξ 0 1 n0 ex jπσ φ n2 w t n T, 34 comoed of the um of delayed inut unit cell, weighted by a uadratic hae factor and the amlitude factor 1/. In thi ection we ue thi reult to analyze a concrete examle. The Talbot array illuminator TAI are multilevel hae grating that concentrate coherent light in a eriodic array of har binary irradiance ditribution of rectangular form, at fractional Talbot ditance of order / after the array [18], [19]. Referring to Fig. 1, we denote the eriod of the hae grating a L and the width of the focuing ot a. The comreion ratio, /L, i a meaure of it concentration caacity, and eual the value of integer [19]. Phae Amlitude Amlitude Phae L X Z 1 Z 4 T FIG. 1. Scheme of a Talbot array illuminator for / 1/4. The TAI hae can be determined when the roagation i conidered in the revered way [21],[22], a hown in Fig. 2. Starting from a binary amlitude grating of oening ratio 1/, which i illuminated by a collimated and atially coherent wavefront, one earche for ditribution of uniform light intenity at the fractional Talbot lane after the amlitude grating. At thee lane the wavefront i comoed of a et of hae level with a eriodicity eual to the original eriod L of the amlitude grating. In thi invere TAI cheme, the deign roblem i imly to determine thee hae level at the fractional Talbot lane. Note that in the atial formalim revering the roagation direction z i euivalent to a comlex conjugation, ince the one-dimenional araxial wave euation i: 2ik 0 A z 2 A x 2, 35 where Ax, z i the comlex amlitude of the araxial wave and k 0 the wavenumber. Therefore, the TAI hae are the comlex conjugated to thoe acuired in roagation. Alternatively, one can reroduce in a multilevel hae grating the hae acuired by fractional / Talbot roagation, without comlex conjugation, and oerate the hae grating at the comlementary / fractional Talbot ditance [22]. The combined effect i the roagation over a / 1 Talbot ditance, and thu the original amlitude grating i reroduced with a halferiod hift. In our reentation, however, we ue the former oint of view a it doe not reuire the change of the Talbot order. Amlitude Phae Phae Amlitude X Z 1 Z 4 T FIG. 2. Scheme of the invere Talbot array illuminator for / 1/4. Within the time-domain formalim the general olution can be found a follow. A in the invere TAI roblem, let u conider the temoral analogue to a binary amlitude grating with oening ratio 1/, which i a ule of rectangular form and width T/. Uing 34, we oberve that the outut unit cell at a fractional Talbot lane of order / i comoed of relica of the baic rectangular ule, each with a uadratic hae factor roortional to the ratio /. According to the reviou develoment, the TAI hae hould be the comlex conjugate. Exlicitly, Φ / n ex jσ φ ξ n ex [ jσ φ ξ 0 + π ] n2. 36 To check thi reult, let u conider the direct TAI geometry. Cw laer light of unit amlitude, analogou to the lane wave that illuminate the grating in Fig. 1, i hae-modulated according to the level in 36, each in a time interval of duration T/. The inut unit cell can be reented a: 1 t w t ex jσ φ ξ n rect T n, 37 n0 where the rectangle function i defined a rectx 1 for x < 1/2 and zero otherwie. To comute the field after

7 7 fractional Talbot we ue 34 and 37 in 4, obtaining: E t m n0 e jσ φξ m ξ t n rect T n m. 38 We have omitted the δt e T term for imlicity, ince it i not relevant in the reent comutation. With the change n + m k thi read: E t k n0 e jσ φξ k n ξ t n rect T k. 39 The um in n i in the form of a convolution, but can be written a the autocorrelation R ξ k δ k,0 mod uing the ymmetry roerty 25. The autocorrelation i nonzero only for k 0 mo, o etting k n for n integer we finally get: E t + n t rect T n. 40 The outut enveloe i the exected rectangle ule train with eriod T and ule width T/. Thi reult can be decribed a follow. The multilevel hae modulation 37 define a erie of time bin of duration T/ within the fundamental eriod T, exlicitly hown by the ucceion of rect function. The hifting and weighting of unit cell roduced by Talbot effect, together with the eriodicity of the the train, tranlate into a hifting and weighting of bin, which lead to the convolution exreed in E. 38. Thi roce ha been exemlified in Fig. 3. If the multilevel hae are choen a the comlex conjugated to the Talbot roagation hae, the bin-wie multile interference become the erfect eriodic autocorrelation 30. The multile interference i therefore contructive for only only bin er eriod, and i detructive for the ret. In other word, the array illuminator relie on the autocorrelation roertie of the uadratic Gau um ariing in Talbot effect. The TAI oeration rincile can alo be analyzed in the ectral domain. The Fourier tranform of the unit cell 37 i: W ν T inc νt 1 e jσ φξ n e j2πνnt/. 41 n0 Since the unit cell i reeated with eriod T, the ectrum of the enveloe Et i comoed of ectral line ν m/t, with m integer: Eν 1 inc νt + m δ ν m T 1 e jσ φξ n e j2πmn/, 42 n0 where incxinπx/πx. The um over n become the DFT 18, and the reult can be reented a: Eν 1 [ ] ex jπσ φ νt 2 νt + inc δ ν m 1 m. 43 T m Thi i the ectrum of a train of uared ule of width T/, amlitude and eriod T, hifted by half a eriod when i odd, and diered by the uadratic hae factor in the firt line. The dierive line after the hae modulator imly comenate for thi hae factor, and render the ule train chir-free and thu of uare temoral rofile. The eak amlitude in 40 i, which amount to a gain factor of in eak ower with reect to the ower of the cw laer light. Notice alo that 37 can be interreted a the multilevel hae modulation of a erfect rectangular ule. In general, had we hae modulated an arbitrary ule contained in an interval of duration T/ and reeated alo in T/ hift, we would have obtained at the outut the ame ule with eriodicity T and a gain factor of in ower. Thi i the rincile of the noiele ule amlification by coherent addition demontrated in [25]. It i illutrative to analyze ome articular cae derived from 36. For intance, it contain the ix binary TAI decribed in [20], aociated to the fractional Talbot lane / 1/2 1, / 1/3 4, / 2/3 2, and the correonding hifted integer of the form + /. Moreover, for 1, and uing the firt reult of Prooition 2, hae 36 take the form: Φ 1/ n 1 ne exjσ φ πn 2 /, 44 which coincide, u to a global hae, with the erie of olution originally reented in [19]. The euivalence i immediate for even, and for odd, and uing their notation, it follow from the ubtitution n I + /2 with I half-integer ranging from N/2 to N/2 1. We alo mention that the et of TAI hae 36 i not the ame et of hae ued to induce the ectral Talbot effect [8] by ue of multilevel hae modulation [9]. The ectral Talbot effect conit in the generation of elf-imaging henomena in the ectrum of a freuency comb, i. e., in the broadband ectrum of a train of ultrahort ule of eriod T, where the et of ectral line earated by a freuency difference of 1/T lay the role of a diffraction grating in the uual configuration of the Talbot effect. In order to induce the ectral Talbot effect it i thu neceary to generate a uadratic hae tranformation in the recirocal domain, in thi cae in time. Such a tranformation i therefore a uadratic hae modulation of the form exjαt 2 /2 for certain value of the chir arameter α. For ultrahort ule, which are eentially located at definite temoral intant t n nt, the ectral Talbot effect i generated

8 8 Multilevel hae ex jξ n 0 T m 0 m 1 m 2 Talbot hae wt ex jξ 0 ex jξ 1 ex jξ 2 exjξ 0 exjξ 3 wt T/3 ex jξ 2 ex jξ 0 ex jξ 1 ex jξ 2 exjξ 1 exjξ 2 wt 2T/3 ex jξ 1 ex jξ 2 ex jξ 0 ex jξ 1 ex jξ 2 exjξ 2 exjξ time FIG. 3. Scheme of the multile interference in time bin in a Talbot array illuminator with / 1/3 and σ φ +1. The unit cell i contained in interval [0, T and divided in three time bin, hown with a gray background in the to row. Each bin contain the correonding multilevel hae. The hifting left and hae weighting right of bin induced by Talbot effect i hown in ucceive row, uing the cyclic tructure of bin inherited by the eriodicity of the unit cell. The reult of the bin-wie coherent um i rereented in the bottom. imly by etting α ±2π T 2. Therefore, the effect can be induced by the et of hae [9]: Φ / n ex ±jπ n2. 45 Thi euence of hae ha, in general, a eriod of 2, contrary to 36 where the eriod i alway. Thi i becaue in 45 ratio /, which i the order of the ectral Talbot effect, can contain two odd number. By contrat, ratio / in 36 i an irreducible fraction of integer with ooite arity, a hown by the third roerty of Theorem 1. Notice, however, that the two et of hae coincide at the elf-dual integer ± 1. In the ret of the ection we etablih the connection with a known characterization of the TAI hae u to an irrelevant global hae and an arbitrary hift of the hae grating by integer multile of T/ [23], which rovide an alternative form of Chu contruction [17]. Thi reult i thu ufficient when the illumination how thi ymmetry, a i the cae of the TAI roblem, but not in more comlex ytem where thi ymmetry i broken, uch a in Talbot line with tructured illumination. The TAI hae are written in [23] in term of the modular invere of mod, denoted a r: r. 46 We firt analyze the relationhi between thi r and our. We notice that, although in ome cae they coincide, thee integer may differ. For intance, for 5 and 8 we have r [1/5] 8 5 but [1/5] In general, integer r and are either eual mod 2 or differ by, ince the condition r 1 mod imlie either r mod 2 or r + mod 2. For even thi mean that, ince mut be odd, r i alo odd, indeendently of the concrete form of the relationhi, r or r +. For odd integer mut be even, and o r + i allowed only when r i odd, i.e., r + e r mod 2, ince thi aure that, a a function of r, i alway even. We are now in oition to exre the connection: Prooition 5. The TAI hae can be written a: ex jπ n2 ex jπ r m2 even ex jπ n2 ex jπ r mm 1 odd, 47 where tand for euality u to a global hae and a oible hift of index n. Proof. We mut analyze cae by cae. For even, the euality i obviou if r, and for r + we have: ex jπ n2 ex jπ r n2 1 n ex jπ r n2 1 nr ex jπ r nn ex jπ r n 2 ex jπ r,

9 9 where in the econd te we have ued that r i odd. The euivalence u to a global hae follow after the hift m n /2. For odd we ue the exlicit formula derived above, r + e r and comute for both cae. For r odd or r +, and uing the ame algebra a in 48, we have: ex jπ n2 ex jπ r nn ex jπ r mm 1 ex jπ after the hift m n [ 1/2], howing the euivalence with 47. Finally, for r even or r we have: ex jπ n2 ex jπ r n2 ex jπ r nn, 50 becaue the additional term ex jπrn 1 rn i unity ince r i even. In thi form, 50 i euivalent to 49 and lead again to 47. The roof i comlete. VI. CONCLUSIONS In thi aer we have reented a novel and imle comutation of the uadratic Gau um that define the Talbot effect. It i baed on the roertie of integer, a arity-deendent modular invere contructed from the two corime number and that define a general fractional Talbot lane. The comutation can be comactly decribed in two cae deending on the arity of, and how that the DFT of a uadratic hae euence i itelf a uadratic hae euence. In addition, both member of the DFT air have erfect eriodic autocorrelation. Thee reult were exemlified by the tudy of Talbot array illuminator, where the erfect autocorrelation arie a a multile interference roce in bin within the fundamental eriod of the grating. Alo, connection with reviou characterization of TAI have been rovided. The reult derived here rereent a imle mean for the deign and analyi of ytem baed on fractional Talbot effect. ACKNOWLEDGMENTS Thi aer wa written during a tay at INRS-EMT, Montreal, Canada, funded by Generalitat Valenciana, Sain through a BEST/2016/281 grant, and evolved from converation with Reza Maram, Lui Romero Corté and Joé Azaña, to whom I am indebted. Thi aer i dedicated to the memory of Prof. María Victoria Pérez, Carlo Gómez-Reino and Felie Mateo, with whom I learnt the beauty of Talbot effect. Aendix A: Some reult from number theory A baic reult in number theory [26], [27], [28], ometime referred to a Bézout lemma, tate that given two integer, a and b, with greatet common divior d, the modular euation ax + by n, i olvable for x and y if and only if integer n i multile of d. In articular, if a and b are corime, the reult guarantee that there exit olution to the euation: ax + by 1. A1 The air of olution x 0, y 0 of A1 are not uniue, and in fact determined mod b and a, reectively, ince x 0 + mb, y 0 ma for any integer m i alo a air of olution. It can be hown that thi erie exhaut the oible air of olution to A1. Let u define x 0 a the lowet oitive integer that belong to a air of olution to A1. x 0 i non-zero and bounded, x 0 < b. The modular multilicative invere, or imly the modular invere of a mod b, denoted a [1/a] b, i thi uniue integer x 0 that verifie: a 1 mod b A2 a b and lie in the range 1 [1/a] b b 1. The ame contruction can be alied to the modular invere [1/b] a, which in thi cae lie in the range 1 [1/b] a a 1. Bézout lemma alo imlie that, given a and b corime, the euation ax + by n admit olution for any integer n. In articular, thi mean that the modular euation ax n n mod b i olvable for any n [28,. 31]. A comlete et of olution i given by x n n[1/a] b mod b, including the cae n 0. If we et n in the range of reidue mod b, 0 n b 1 and reduce the value of x n to their reidue mod b, we have thu defined a bijective ma n x n within the reidue mod b. Thi i the ma that ha been ued to erform change of variable of the tye x y ax mod b for a, b corime in invariant exreion mod b in everal art of the aer. The Jacobi ymbol a b of an arbitrary integer a and an odd and oitive integer b, which are mutually rime, i defined a the roduct of the Legendre ymbol of the rime factor of b. To be ecific, if the decomoition of b in rime factor i b α1 1 α N N, then: a α1 αn a a A3 b 1 N where for rime j, the Legendre ymbol i: a +1 if there exit an integer x uch that x j 2 a mod j 1 otherwie. A4 Acoording to thi definition, the Legendre ymbol a/ j eual +1 if and only if a i a uadratic reidue mod j,

10 10 otherwie it i 1. In articular, 1/ j +1 for any j, and a/ j a / j if a a mod j. Proertie and value of the Legendre ymbol can be conulted in [26], [27], [28] and in the econd aendix of [13]. The Jacobi ymbol A3 verifie the ame roertie a the Legendre ymbol, and alo verifie that it value i +1 when a i a uadratic reidue mod b. The convere, however, i not necearily true, ince for a being a uadratic reidue of b it i neceary that a i uadratic reidue of all of it rime factor j, and o all factor in A3 mut be +1. Aendix B: Euivalence with alternative exreion of integer The comutation of exjξ n in [13], corrected in [14] and a reented in [15], i imilar to 19 by with a different exreion for. We how here the euivalence of the integer in 7 and 8 with the reidue mod 2 of the three tye of integer defined in [14], [15], here denoted a, by howing that they verify the condition tated in Theorem 1. In the firt cae we alo how the connection by exlicit maniulation on. a odd and odd. The integer i given by: [ ] [ ] mod B1 Uing that [1/2] + 1/2 [13] [27], we can imlify thi exreion to: [ ] mod 2. B2 2 Incidentally, we mention that the value in the original comutation [13] wa: [ ] 2 1 HB 4 mod 2. B3 4 Returning to B2, we recall the definition of the modular invere: 2 1 mod, B4 2 we can be further imlify B2 to: 2 [1 mod ] mod mod 2, 2 B5 which coincide with 8. Alternatively, we can how the euivalence by a direct check of the hyothee of Theorem 1. Firt we rehrae B4 a: a B6 2 for certain integer a, which hould be odd ince i odd and the left hand ide of the euation i even. Then the roduct i: a 2 1+ mod 2, B7 which i the defining euation 6 in Theorem 1. Thi econd form of the euivalence i comleted by oberving in B1 that the arity of i even, and i thu ooite to that of. b odd and even. The integer i defined a: [ ] 2 1 mod 2. B8 Uing the definition of modular invere, 1 mod 1 + b B9 for certain integer b, we have that the roduct i: 1 + b b + b 2 1 mod 2, B10 ince i even. Thi i the defining euation 6 for even. We finally notice from B9 that [1/] mut be odd, o i odd and thu of ooite arity to. c even and odd. The integer i: [ ] 2 1 mod 2. B11 The definition of modular invere B9 now imlie that integer b i odd, ince i odd and i even. Then, the roduct i: 1 + b mod 2, B12 which i again the defining euation 6 for odd. The roof i comleted by noticing that i now even ince i even. Aendix C: Comutation of hae exjξ 0 We comute the Gau um 24, which lead to the reult 19. A in [13], the comutation i reduced to the alication of the tandard form of the Gau um found in text of number theory [27,. 86]. Thi tandard form alie to um of uadratic hae of the tye exj2πan 2 /b where a i an arbitrary integer and b i odd. In the cae of odd, we are in thee condition ince i even. Then, /2 e jξ0 e jπ 12 /8 2 e jπ 12 /8 e jπ 1/4, C1

11 11 where in the econd te we have ued the numerator roduct rule of Jacobi ymbol and the value 2/ ex[ jπ 2 1/8]. Thi give the firt formula in 19. Now, we ue 6 and again the numerator roduct rule to reach an exreion involving and : 1 + e 1 1, C2 and o / / and the firt line of 19 i roved. For even, i odd and we have to invert the fraction in the exonential of 24. Uing the recirocity of the Gau um [13], [14], [26] we are led to: e jξ0 1 1 e jπ m2 e jπ/4 1 e jπ m2. C3 The ue of the tandard form give the firt exreion of 19 for even: e jξ0 e jπ/4. C4 Since the original Gau um in C3 i invariant under hift + 2, o it i the right hand ide of C4, and therefore: + 2n 1 n/2. C5 for any n 0. In articular, we can evaluate C4 uing any other rereentation of euivalent mod 2. With the rereentation [1/] 2 from Aendix B we can traightforwardly reduce thi reult to an exreion involving and. Firt, 2 [1/], C6 where we have ued the denominator roduct rule of Jacobi ymbol and the fact that they only take value ±1. And econd, we oberve that ince [1/] 2 i odd, [1/] mut alo be odd. Therefore, etting 2t+1 2 for certain integer t we get: e jπ/4 e jπ2t+12 /4 e jπ/4, C7 which together with C6 give the econd exreion in 19 for even. [1] H. F. Talbot, Philoohical Magazine 9, [2] L. Rayleigh, Philoohical Magazine 11, [3] J. T. Winthro and C. R. Worthington, Journal of the Otical Society of America 55, [4] J. Azaña and H. G. de Chatellu, Phyical Review Letter 112, [5] S. Akhmanov, A. Sukhorukov, and A. Chirkin, Soviet Journal of Exerimental and Theoretical Phyic 28, [6] B. H. Kolner, IEEE Journal of Quantum Electronic 30, [7] J. Azaña and M. A. Muriel, IEEE Journal of Selected Toic in Quantum Electronic 7, [8] J. Azaña, Otic Letter 30, [9] J. Carauitena, M. Beltrán, R. Llorente, J. Martí, and M. A. Muriel, Otic Letter 36, [10] K. Patorki, The elf-imaging henomenon and it alication, in Progre in Otic, Vol. 27, edited by E. Wolf Elevier, [11] J. Wen, Y. Zhang, and M. Xiao, Advance in Otic and Photonic 5, [12] M. V. Berry and S. Klein, Journal of Modern Otic 43, [13] J. H. Hannay and M. V. Berry, Phyica D: Nonlinear Phenomena 1, [14] S. Matutani and Y. Ônihi, Foundation of Phyic Letter 16, [15] L. Romero-Corté, H. G. de Chatellu, and J. Azaña, Otic Letter 41, [16] H. D. Luke, IEEE Tranaction on Aeroace and Electronic Sytem 24, [17] D. Chu, IEEE Tranaction on Information Theory 18, [18] A. W. Lohmann and J. A. Thoma, Alied Otic 29, [19] J. R. Leger and G. J. Swanon, Otic Letter 15, [20] V. Arrizón and J. Ojeda-Cataneda, Otic Letter 18, [21] P. Szwaykowki and V. Arrizón, Alied Otic 32, [22] V. Arrizón and J. Ojeda-Cataneda, Alied Otic 33, [23] C. Zhou, S. Stankovic, and T. Tchudi, Alied Otic 38, [24] The firt row 2 in the table of [15] i in error. It read but it hould read [25] R. Maram, J. Van Howe, M. Li, and J. Azaña, Nature Communication 5, [26] E. Hecke, Lecture on the Theory of Algebraic Number, 2nd ed., Graduate Text in Mathematic, Vol. 77 Sringer, [27] S. Lang, Algebraic Number Theory, 2nd ed., Graduate Text in Mathematic, Vol. 110 Sringer, [28] K. Ireland and M. Roen, A Claical Introduction to Modern Number Theory, 2nd ed., Graduate Text in Mathematic, Vol. 84 Sringer, 1990.