Engineering the spectra of photon triplets generated from micro/nanofiber

Chuan Qu(瞿川), Dongqin Guo(郭東琴), Xiaoxiao Li(李笑笑), Zhenqi Liu(劉振旗), Yi Zhao(趙義),Shenghai Zhang(張勝海), and Zhengtong Wei(衛(wèi)正統(tǒng))

The College of Basic Department,Information Engineering University,Zhengzhou 450000,China

Keywords: photon triplets,micro/nanofiber,spectrum engineering

1.Introduction

Photons are well suited for the implementation of major quantum information processing (QIP) tasks, such as quantum computation,[1,2]quantum teleportation,[3]quantum key distribution,[4]and quantum metrology.[5]These photonicsbased QIP creates demands for sources of single photons[6,7]and of multiple photons[5,8]in quantum-entangled states.Indeed, over the past few decades, spontaneous parametric down-conversion (SPDC) in second-order nonlinear crystals,as well as spontaneous four-wave mixing (SFWM) in thirdorder nonlinear fibers,have emerged as the primary choices for entangled photon pair sources in many QIP experiments.[9,10]Importantly, photon triplet states have inherent advantages in generating Greenberger-Horne-Zeilinger (GHZ) states[11]and the heralded generation of photon pairs.[12]Several methods have been reported for generating photon triplets, including cascaded photon pair processes[13]and quantumdot molecule schemes.[14]Nevertheless, these methods are plagued by extremely low collection efficiencies.Third-order spontaneous parametric down-conversion (TOSPDC), the inverse process to third-harmonic generation (THG), where a pump photon is annihilated to simultaneously give birth to a photon triplet governed by energy and momentum conservation, may lead to the generation of GHZ states without post-selection.[15-17]Further,TOSPDC introduces three-mode squeezing operators to directly facilitate the realization of non-Gaussian states,[18,19]while two-mode squeezing operators lead to Gaussian squeezed states.Despite TOSPDC’s promising prospects for QIP, it confronts some technological challenges due to the weak third-order nonlinearity of typical optical materials and difficulties in achieving phase matching.[20]

A micro/nanofiber is an optical fiber with a core diameter approaching the submicron scale, a size comparable to the wavelength of the transmitted light.Here, the effective nonlinear-optical coefficients are significantly enhanced due to the reduced mode area and field enhancement that results from tight confinement.[21,22]The typical method for fabricating micro/nanofiber involves heating and stretching standard optical fibers until they reach a predetermined diameter.[23]In addition, the use of direct mode cutoff feedback can significantly enhance the accuracy and precision of real-time diameter control during the fiber-pulling process.[23,24]Note that the cladding of the original fiber acts as the core confining light,while the surrounding air serves as the new cladding.It is such a large step in the refractive index that tightly confines the mode inside the core.Furthermore,phase matching is actually dependent on the chromatic dispersion of the fiber.The ability of a micro/nanofiber to support multiple transmission modes enables phase matching of TOSPDC through so-called intermodal phase matching, where the pump operates in a highorder mode while the photon triplets are in the fundamental mode.Also,waveguide dispersion contributes greatly to chromatic dispersion, suggesting that photon-triplet wavelengths can be widely tailored by changing the micro/nanofiber diameter.The pigtail of the micro/nanofiber is retained as a standard optical fiber, which is beneficial for accessing the fiber quantum network with minimal coupling loss.

In this work, we study the spectrum engineering of photon triplets generated from micro/nanofibers and longperiod micro/nanofiber gratings.A multitude of theoretical and experimental studies have focused on realizing phase matching for degenerate signal frequencies to improve the efficiency of photon triplet generation, using methods such as fiber dispersion tuning,[25,26]nonlinear phase modulation enhancing[27]and quasi phase matching (QPM).[28-30]Some studies have validated these design schemes experimentally through THG.[17,20]Nonetheless, spectrum engineering, especially for non-degenerate photon triplets, has rarely been reported.In fact, phase mismatching at one-third pump frequency gives rise to non-degenerate photon triplets without a decrease in efficiency.An increased phase mismatching corresponds to a broader signal bandwidth,indicating that the collection efficiency of photon triplets decreases and the noise increases.Further,the photon triplets need to be separated into three channels for practical applications in QIP.To the best of our knowledge,we provide a frequency-division scheme with high heralding efficiency for the first time.Moreover,we propose a QPM scheme to generate tunable-wavelength photon triplets in a long-period micro/nanofiber grating.The results presented can also be extended to many optical materials and waveguide geometries where TOSPDC occurs.

This paper is organized as follows.In Section 2, we introduce the quantum theory analysis of TOSPDC and give the expression of joint spectral amplitude.In Section 3,we investigate the ellipse locus of joint spectra, and based on this, in Section 4,we propose a frequency-division scheme to separate non-degenerate photon triplets into three channels.Also, we study the tunable-wavelength photon triplets based on QPM in Section 5.Conclusions and some perspectives are drawn in Section 6.

2.Photon triplet states

The TOSPDC process is a third-order optical nonlinear process, originating from the third-order susceptibilityχ(3).The annihilation of individual photons from the pump modes gives birth to photon triplets, as shown in Fig.1.The three emitted signal modes are referred to as signal-1 (r),signal-2 (s), and idler (i) with angular frequenciesωr,ωsandωi, respectively.Pump angular frequency is denoted asωp.This TOSPDC process occurs by satisfying the energy conservationωr+ωs+ωi=ωpand phase-matching condition Δβ=βp-βr-βs-βi-βNL=0,whereβj(j=p,r,s,i)is the mode propagation constant for the four participating fields andβNLis the nonlinear contribution resulting from cross and selfphase modulation.[22]Δβis known as the phase mismatch.The light-matter interaction Hamiltonian for the TOSPDC is given by[20,31]

Fig.1.The TOSPDC process in micro/nanofiber.A pump photon decays into photon triplets.Additionally, the intensity distributions of the two mode fields are shown below.

whereχ(3)is the cubic susceptibility, ?0is the vacuum permittivity and the integral is evaluated over the cubic interaction volumeVint.The subscript i denotes the idler, and the other i signifies the imaginary unit.We describe classically the strong pump fields, in terms of monochromatic pump, its positive-frequency components can be written as Therefore, we obtain the photon triplet states in terms of the fiber lengthL

and joint spectral intensity (JSI)|?(ωr,ωs,ωi)|2is related to the probability of photon triplets emitted at frequencies ofωr,ωs, andωi.In Eq.(6), the product termγ2Iωrωsωi/ω2phas a slowly-varying dependence on frequency within the spectral range of interest,[15]such that we neglect this dependence and characterize the spectral properties of photon triplets with the joint spectral amplitude(JSA)?(ωr,ωs,ωi),given by

In this case, the JSI|?(ωr,ωs,ωi)|2=L2sinc2(ΔβL/2) and thus phase matching Δβ=0 corresponds to the most efficient photon-triplet emission.

Generally,the fulfillment of phase matching resorts to the so-called intra-modal phase matching,in which the pump is in a higher-order mode while the triplet photons are all in the fundamental mode(HE11mode).Here,the micro/nanofiber used to generate photon triplets has a submicron diameter,such that changing diameter drastically influences the contribution of waveguide dispersion to fiber dispersion.Figure 2(a) shows the effective refractive index of various modes versus the micro/nanofiber diameter.Here, the black curve denotes the fundamental mode with a wavelength of 1551 nm(angle frequencyω1), while the colorful curves represent higher-order modes with a wavelength of 517 nm (angle frequencyω3,ω3=3ω1).In the case of continuous-wave(CW)pump with low peak power,the nonlinear phase mismatchβNLis negligible and the phase matching is rewritten asneff(ω3)=neff(ω1)corresponding to the intersections between black curve and colorful curves.Even if multiple higher-order modes enable one to fulfill the phase matching of the TOSPDC process,most of them are inaccessible due to their poor overlap integralsfprsiand the difficulty in coupling pump into higher-order modes.The favorable regime occurs when the visible pump light is guided in the HE12mode,and the micro/nanofiber diameter is 767 nm,corresponding to the black dot in Fig.2(a).The pump wavelength ofλp0= 517 nm and the diameter ofd0=767 nm are referred to as reference values.Moreover,the intensity distributions of the HE11(ω1)mode and the HE12(ω3)mode are shown in Figs.2(b)and 2(c),respectively.It turns out that the overlap integral for the combination of HE11(ω1) and HE12(ω3) is the most efficient for generating photon triplets.[20]

Fig.2.(a) The dependence of the effective refractive index neff on the micro/nanofiber diameter for various higher-order modes at a wavelength of 517 nm and for the fundamental mode at a wavelength of 1551 nm.(b) and (c) The intensity distributions in terms of HE11(ω1) mode and HE12(ω3)mode,respectively.

3.Joint spectra with ellipse locus

In the process of TOSPDC, the pump wavelength and micro/nanofiber diameter may deviate slightly from their reference values, whereas the rules of energy conservation and phase matching can still be satisfied resulting from nondegenerate photon-triplet frequencies.In addition, nondegenerate TOSPDC with phase matching also maintains efficient photon-triplet emission, unlike the process of THG.Figure 3 shows the results of simulated photon-triplet JSI in the space of{ωs,ωi,ωr}with brighter one representing higher probabilities of emission.As for the monochromatic pump,the JSI is pasted upon the plane ofωs+ωi+ωr=ωpdue to energy conservation.Thus, we project JSI onto the three coordinate planes to obtain each marginal distribution.Here,the pump wavelength is shifted by-0.2 nm withλp=516.8 nm.Obviously, the non-degenerate JSI has a form of closed-loop belt and is absent from an emission maximum at the loop center with frequency ofωp/3.

Fig.3.JSI in the three-dimensional space of{ωs,ωi,ωr}in terms of the frequency non-degenerate configuration.The violet plane stands for the plane of ωs+ωi+ωr=ωp.Three marginal distributions are projected on the corresponding coordinate planes.

In the case of non-degenerate frequency,R ＞0 gives rise to an ellipse-belt JSI in shape, while in the case of degenerate frequency,R= 0 gives rise to an ellipse-cake JSI in shape.The highest emission probabilities occur at frequencies satisfying the ellipse equation in Eq.(10) and they get identical probabilities for both of degenerate and nondegenerate cases above.Furthermore, degenerate photon triplets get a Gaussian-approximation output spectrum, while non-degenerate photon triplets get a concave output spectrum,as shown in Figs.4(b)and 4(d).The two highest values of the concave spectrum correspond to the left/right extreme points of ellipse;using analytic geometry methods,the bandwidth of the concave spectrum readsNote that there are multiple weak ellipse belts,attributing to the assistant peaks in sinc function plots,thus they rapidly oscillate and decay.

Fig.4.JSIs in the {ωs,ωi} plane, corresponding to the marginal distribution, for the non-degenerate regime (a) and degenerate regime (c).Both JSIs are normalized by the total conversion probability.Additionally, their normalized output signal spectra are shown in panels (b) and(d), respectively.In panel (a), the red dashed curve, originating from Eq.(10),outlines an ellipse locus representing the highest emission probability.Two white arrows denote the major axis δ- and minor axis δ+,and they have 135° and 45° in intersection angle with the ωs axis, respectively.The location of the ellipse center is（2ωp/3,0）in {δ+,δ-}coordinate plane.

4.Frequency-division scheme to separate photon triplets into three channels

A critical procedure for a TOSPDC photon-triplet source to be practically applied to quantum information technology is to separate photon triplets into three frequency channels.Hence,detecting one of the photon triplets heralds the remaining photon pairs.Fortunately,in terms of non-degenerate photon triplets,the fixed eccentricityin the ellipse locus of JSI provides an option of the frequency-division scheme with high-enough heralding efficiency.The frequency-division scheme we give is shown in Fig.5, wherein the simulation parameters of JSI are identical to that in Fig.4(a).Here, the diagonal line in equation ofωs=ωiintersects the ellipse of perfect phase matching at two points,accordingly,we can divide the whole spectrum into three channels.

We assume an ideal rectangular filter and that channels 1,2 and 3 are of identical bandwidth.Therefore, we can obtain visually the channels in which s-signal and i-signal photons are located from the marginal JSI in the{ωs,ωi}plane.Identifying the channel where r-signal photons are located is crucial for increasing efficiency.Ifωs+ωi=2ωp/3+νwithνrepresenting the frequency deviation,we haveωr=ωp/3-ν.For example,in terms of point A in Fig.5,νindicates the vertical distance from the white dashed lineωs+ωi=2ωp/3.Then we move the point（ωp/3,ωp/3）with-νto obtain the channel of r-signal photons.Particularly, the heralding efficiency in terms of Fig.5 is up to 94.4%.

Fig.5.Frequency-division scheme to separate non-degenerate photon triplets into three channels.The frequency-division channels 1, 2 and 3 are divided by the red lines.The coordinate equation for the white dashed line is ωs+ωi =2ωp/3 and the coordinate equations for the two solid lines are ωs+ωi=2ωp/3±R/.Thus,these shadow regions represent that photon triplets can not be absolutely separated, resulting in a slight decrease in heralding efficiency.

In Fig.5,these shadow regions correspond to that photon triplets may not be absolutely separated into three channels,while other regions correspond to complete separation.Obviously, a degenerate source of photon triplets will give rise to a low heralding efficiency.In fact, the decrease in heralding efficiency mainly arises from the thickness of the ellipse belt as well as the assistant peaks of the sinc function.Both of the above-mentioned flaws can be mitigated by increasing the fiber length (see Eq.(9)).In addition, the broad ellipse loci suffer from a lot of fluorescent background,such that the frequency-division bands collecting signal photons should be set narrow enough to increase the signal-to-noise ratio.

The frequency-division scheme is realised using a frequency divider.The fabrication of this frequency divider can utilize the same technologies as the wavelength-division multiplexer in fiber communication systems.Furthermore, in order to increase the heralding efficiency,the frequency-channel edge should be steep and the channel position in the spectrum should align with the preset value.

Due to a fixed eccentricity of the ellipse locus,the bandwidth of photon triplets depends on the minor semi-axis lengthR.Figure 6 showsR2versus pump wavelength shiftsλp-λp0in terms of various micro/nanofiber diameter deviationsd-d0.The black dashed line marksR=0,indicating degenerate photon triplets.A positiveR2corresponds to non-degenerate photon triplets, while a negativeR2represents an unaccessible generation of photon triplets in the absence of phase matching.Obviously, around the reference value, the simultaneous decrease in micro/nanofiber diameter and increase in pump wavelength will result in the disappearance of photon triplets.As shown in Fig.7,the ellipse width can be adjusted by changing the diameter and pump wavelength.Due to the square root function form ofR, the ellipse loci in Figs.7(a) and 7(b) become denser from the inside out.Furthermore, changing the pump wavelength will move the ellipse center,while changing the diameter will keep the ellipse center fixed.Importantly,even if the diameter deviation from the reference value,resulting from micro/nanofiber fabrication error, may broaden the bandwidths of the photon triplets or destroy phase matching,choosing an appropriate pump wavelength can compensate for this error and give rise to a narrowband,non-degenerate photon triplet source, as illustrated in Fig.7(c).In addition, as can be seen from Fig.7(a),a pump with a wider spectrum will result in broadening the thickness of the ellipse belt, thereby lowering the heralding efficiency.Therefore,it is necessary to choose pumps with narrow enough bands.

Fig.6.The R2 versus pump wavelength shifts Δλp=λp-λp0 in terms of various micro/nanofiber diameter deviations Δd=d-d0. λp0=517 nm and d0 =767 nm are the reference values for pump wavelength and micro/nanofiber diameter,respectively,to perform degenerate perfect phase matching.Point A corresponds to Δλp=0 and Δd=0.

Fig.7.(a) Ellipse loci corresponding to various Δλp of 0 nm, -1 nm,-2 nm, -3 nm and -4 nm in the order from the inside out.(b) Ellipse loci corresponding to various Δd of 0 nm,5 nm,10 nm,15 nm and 20 nm in the order from the inside out.(c)JSI with Δλp=-6.42 nm and Δd =-10 nm.Some plots marked with capitals A-F correspond to the dots in Fig.6.

5.Tunable-wavelength photon triplet generation

Many applications of quantum information technique require multiple wavelengths, thus it is demanding to control the phase matching for adjusting the photon triplet wavelength.Phase matching mainly depends on the dispersion characteristics of fibers, whereas after finishing the fabrication of micro/nanofiber,its waveguide and material dispersion both are fixed, indicating a poor tunability of photon-triplet wavelength, as illustrated in Fig.6.Here, long-period micro/nanofiber grating emerges as a viable candidate for realizing tunable-wavelength photon triplet source, in which the periodic variation of dispersion along fiber length provides a new controllable degree of freedom to tailor phase matching.In general, the fiber grating is modeled as adding cosine oscillation term to the original refractive index, i.e.,n(z)=n0+Δncos(2πz/Λ),wheren0is the original refractive index,Δndenotes modulation depth of refractive index andΛis the grating pitch.[32]Note that Δnperforms differently for different fiber modes.Therefore, we rewrite the phase mismatching as

whereapresents the original phase mismatch without the effects of grating, whileb=Δn3ωp/c-Δn1(ωs+ωi+ωr)/c,reduced tob=ωp(Δn3-Δn1)/caccording to energy conservation.Here,Δn3and Δn1denote the modulation depth of refractive index for HE12and HE11modes,respectively.Further,the accumulated phase mismatch from 0 tozcan be written as

Figure 9 shows the JSIs ofq-order QPM by using longperiod micro/nanofiber gratings.The preset ones areq=2 andq=-2, respectively.JSI shows multiple concentric ellipse loci due to QPM, and each ellipse locus corresponds to aq.From the inside out,qincreases sequentially.Indeed,the minor semi-axis length of theq-order ellipse is given by

On the one hand,as shown in Fig.9(a),no photon triplets are generated in such a homogeneous micro/nanofiber configuration due to the absence of phase matching, while QPM enables the regeneration of photon triplets.On the other hand,as shown in Fig.9(c), QPM can reduce the width of ellipse and thus improve the signal-to-noise ratio.Here,q=0 corresponds to the homogeneous micro/nanofiber.

Fig.8.The plots of Bessel function Jq(x)for various q.Parity,Jq(-x)=(-1)qJq(x).

Fig.9.(a) and (c) JSIs for the q-order QPM.Panels (b) and (d) are the normalized output signal spectra of TOSPDC corresponding to panels (a) and (c), respectively.The diameter of micro/nanofiber is the reference value 767 nm, fiber length L=5 cm and grating pitch Λ =800 μm.Δn3-Δn1 =0.00198 and λp =519.75 nm for panels (a) and(b),Δn3-Δn1=0.00196 and λp=514.24 nm for panels(c)and(d).

To enhance the intensity of presetq-order ellipse locus and reduce the influences of other orders,the modulation depth of refractive index needs to be set at a proper value in whichbΛ/2πcorresponds to the first extreme point ofJq(x).Thus,it is necessary to change the modulation depth of the refractive index as needed.Here, we provide a feasible suggestion that mechanically induced long-period fiber gratings,as discussed in Refs.[33,34].Moreover,the thickness of outer ellipse belts is too thin to accumulate stronger intensity contributions.Figures 9(b) and 9(d) show that the outer ellipse loci contribute much to the overall output spectra,but these contributions are mainly concentrated outside the frequency-division band, at the same time,the central part of the spectra corresponding to the preset ellipse locus exhibits a very prominent intensity.Using the frequency-division scheme presented in Section 4,we can effectively eliminate the influences from undesired orders of QPM.

6.Conclusion

We theoretically investigate the engineering of photontriplet spectra generated from micro/nanofiber.Firstly, we provide the expression of JSA in the process of TOSPDC.Indeed,the JSI is correlated with the probability of photon-triplet emission.The most efficient emission occurs when phase matching is fulfilled,resorting to intra-modal phase matching

where the pump is in HE12mode while the photon triplets are all in HE11mode.Further, the two-dimensional JSI shows an ellipse locus with a fixed eccentricity ofAccordingly,we present a frequency-division scheme to separate photon triplets into three channels with high heralding efficiency.The decrease in heralding efficiency primarily results from the thickness of the ellipse belt as well as the assistant peaks of the sinc function.Both of them can be mitigated by increasing the fiber length.In particular, the width of the ellipse locus depends on the phase mismatching at one-third pump frequency,such that one can adjust the width of the ellipse locus via changing pump wavelength or micro/nanofiber diameter.Importantly, choosing an appropriate pump wavelength can compensate micro/nanofiber fabrication errors and give rise to a narrowband non-degenerate photon triplet source with a high signal-to-noise ratio.Notably, long-period micro/nanofiber gratings exhibiting periodic oscillation of dispersion along the fiber length provide a new controllable degree of freedom to tailor phase matching.Thus, QPM is introduced and plays a dominant role in the generation of tunable-wavelength photon triplets.We believe that this work provides a unique pathway towards tunable-wavelength photon triplet sources with high signal-to-noise ratios for applications in quantum information technologies.

Acknowledgments

Project supported by the National Natural Science Foundation of China (Grant No.61605249) and the Science and Technology Key Project of Henan Province of China (Grant Nos.182102210577 and 232102211086).