<<

arXiv:1905.03459v1 [quant-ph] 9 May 2019 bouevleo h etra u fmomenta the of if sum zero vectorial from the different be of can value group absolute a of the tons m htn.Frasnl htnwt h the with particles, massless single of a groups For for and . with theory”. ticles relativity of “compati- language mass standard a the of with pro- definition (1) ble reasonable Equation only [2], the Okun vides L.B. space. by Minkowski argued 4-dimensional as the Moreover, re- in of with rotations “length” is to which squared spect group, the the of to 4-momentum of proportional the line is second (1) the in Equation expression the as Lorentz-invariant, -oetm lal,dfie nti a,temass the way, this in defined Clearly, 4-momentum. sls than less is sapi fnnolna htn iheulfrequencies equal with angle some photons and noncollinear of pair a is Fgr ) o hc qain1 yields 1) Equation which for 1), (Figure nryadmmnu ftegop and group, the of momentum and energy where asmmnu”interrelation mass-momentum” m smxmlat maximal is n h aso h ru qaszr.Temass The zero. equals group the of mass the and At igephot single h ento fEuto 1 svldbt o par- for both valid is (1) Equation of definition The skon[,2,i h eaiitcpyistemass the physics relativistic the in 2], [1, known As fagopo atce sdtrie yte“energy- the by determined is particles of group a of θ m tecase (the 0 = 2 i c 4 ueae particles, numerates = ,a thst e u o ruso pho- of groups for But be. to has it as 0, = ~ ~ r on hnteLrnzivratms srltddirect related is mass parame Lorentz-invariant spontaneous the parameter of when process found the are in formed states biphoton rnvrewv etr feitdphotons. emitted of vectors wave transverse  P k m X θ h ocp fteLrnzivratms fagopo parti of group a of mass Lorentz-invariant the of concept The = max θ i i .INTRODUCTION 1. = ω m ewe ietoso hi aevectors wave their of directions between oet-nain asadetnlmn fbpoo state biphoton of entanglement and mass Lorentz-invariant ε ~ i π /c b i ( ω/c =  θ rpgto fpoosi collinear is photons of propagation ) ≡ tecase (the h ipeteapeo hskind this of example simplest The . = ) 2 m − K c qain()gvsimmediately gives (1) Equation , 2 ( P π  2 ≫ c = ) ~ tot X 2 ω i .V Vintskevich V. S. eemnn ihdge fetnlmn fabpoo sta biphoton a of entanglement of degree high a determining 1 · ε sin 3 c ~p 2 tot P ainlRsac nvriyHge colo Economics, of School Higher University Research National and ) c i ~  2 tot ω  2 and c θ 2 , 0ysisaaUis,Mso,110,Russia 101000, Moscow, Ulitsa, 20Myasnitskaya = 2 usa cdm fSine ocw Russia Moscow, Science, of Academy Russian 1  ocwIsiueo hsc n Technology, and Physics of Institute Moscow ≡ 2 4 ..PohrvGnrlPyisInstitute, Physics General Prokhorov A.M. . cλ ~p π ε tot ~ P tot 2 , ntttkiPr ,Dolgoprudny, 9, Per. Institutskii tot ocwRgo 470 Russia 141700, Region Moscow r h total the are − stetotal the is c 1 2 .A Grigoriev A. D. , ∗ ~p fedorovmv@gmail ~ ~ tot 2 ω P m m and i (1) (3) (2) ( ~ θ k is ω ) i ru fprilsi mle hntesedo , of speed the than smaller is particles of group if Evidently, aeo n ru fprilsisma eoiyo motion of mean by: its given particles is of group any of case n ntecsso olna rpgto ( propagation collinear of cases the in and qain()gives (3) Equation wn onnolnaiyo htnpoaainwsseen was velocity propagation [3]. mean photon experimentally of the noncollinearity of to decreasing owing states, two-photon For ( case general a in photons two of vectors Wave 1: Figure where ietdoinain( orientation directed - reso antd mle hnteeeto mass the than smaller magnitude ∼ of orders 5-6 mlms ftopoosi ahrsal .. at E.g., small. rather is photons two of mass imal h ytmudrcnieainmvsa hl iha with whole a as moves speed consideration under system the mass. electron the than larger even or becomes hr t envlct un eo ..ters frame rest the i.e. zero, turns velocity mean its where h aewv vectors wave same the o n betwith object any For xsec fanneoms en meitl that immediately means mass nonzero a of Existence 10 . − 1 com , 27 2 λ v .V Fedorov V. M. , ywt pootoa o h Schmidt the to) (proportional with ly .Hwvri ytmof system a in However g. stepoo aeegh ueial,temax- the Numerically, wavelength. photon the is mle hntesedo light of speed the than smaller N rcdw ovrin h conditions The conversion. down tric ie agradcnbcm oprbewith comparable become can and larger times lsi hw ob plcbeto applicable be to shown is cles v m = h enpoaainsedo a of speed propagation mean the 0 6= c 2 ε p m c tot ). tot max m ~ k = 2 1 , and 0 6= 3 ∼ and c ewt epc to respect with te ∗  4 1 × ~ k − 2 h aso h group the of mass the , 10 ε m c < v tot 2 − N 2 33 s  htnpiswith pairs photon ,wihi about is which g, 1 / hr saframe a is there c b 2 ,adcounter- and ), nageneral a In . . λ c < v 1 = (4) a µ ) . 2

(r.f.). Of course, for groups of particles in the rest frame 2. SPDC particles can move but only in such a way that the vec- torial sums of their and momenta compensate each other and give zeros in the sums, i.e. Let us consider a rather simple case of the collinear frequency-degenerate regime of SPDC with the type-I r.f. r.f. ~ptot =0, v =0. (5) phase matching. This means that the pump propagates in a nonlinear crystal as extraordinary and both emitted This clarifies the physical meaning of the Lorentz- photons as ordinary waves. If the pump is vertically po- invariant mass of a group of particles: multiplied by c2, larized, both emitted photons have the same horizontal the mass m coincides with the total energy of a group in polarization. Frequencies of the pump and of both emit- its rest frame ted photons are assumed to be equal, correspondingly, to ωp and ωp/2. Let 0z be the central pump-propagation r.f. εtot direction, and distribution of the pump in the trans- m = 2 . (6) c verse directions ( 0z) be Gaussian with the width wp: 2 ⊥ 2 Ep exp( ~r⊥p/2wp). The momentum-representation An example of the rest frame for a pair of photons is that wave∝ function− of emitted photons depends on the trans- shown in Figure 1c. verse momenta of emitted photon ~k⊥1 and ~k⊥2 and is As classical light fields consist of photons, the defini- known [6] to be given by tion of the Lorentz-invariant mass is applicable also to light pulses considered as relativistic objects [4, 5]. In this we will consider manifolds of photons produced 2 in nonlinear birefringent crystals under the action of a ~ ~ 2 k⊥1 + k⊥2 wp classical pump in the process of Spontaneous Parametric Ψ(~k⊥ ,~k⊥ )= N exp 1 2   Down-Conversion (SPDC). As known, such photons can " − 2 #× be entangled. The key question to be addressed below 2 Lλp ~ ~ is whether there is any connection between the Lorentz- sinc k⊥1 k⊥2 , (10) 8πno − invariant mass of SPDC photons and the degree of their     entanglement. As for the pump, its Lorentz-invariant mass has to be sin x found too, at least as the benchmark for comparison with where sinc = x , L is the length of a crystal (along the that of SPDC photons. For the pump we can use the z-axis), no is the refractive index of the ordinary wave in result of our previous work [5]: in Equation (3) of Ref. the crystal, N is the normalization factor. [5] the mean propagation speed of a diverging Gaussian It should be noted that, in principle, anisotropy of the pulse was found to be given by refractive index of the extraordinary pump wave, np, can 2 give rise to the additional term in the argument of the λp ~ ~ ′ v = c 1 , (7) sinc-function, fanisotr(k⊥1, k⊥1) = Lnp(k1 x + k2 x)/4π, − 8π2w2 ′ p ! where np is the derivative of the refractive index np over the angle ϑ between the pump wave vector kp and the where wp is the pump waist. By comparing this result optical axis of a crystal OA, with OA assumed to be with that of Equation (4) at m εtot we find the pump- located in the (x, z)-plane. If fanisotr is not small (com- pulse Lorentz-invariant mass ≪ pared to 1), it plays a very important role by provid- ing anomalously high degree of angular entanglement of λp SPDC photons [7],[8]. But in this work we assume that mp tot = εtot 2 . (8) 2πc wp the term fanisotr is small and can be ignored. Specif- ically, the condition justifying this assumption has the ~ ′ ′ By assuming that εtot = N ωpwhere N is the number of form L np wp, where typically np 0.1. In the photons in the pump pulse, we can find from Equation opposite| case| ≪ of strongly pronounced| anisotropy,| ∼ entan- (8) the Lorentz-invariant mass per one photon glement of SPDC photons was investigated in the works [7, 8], and the arising in this case peculiarities of deriva- ~ tion of the Lorentz-invariant mass will be considered else- mp = . (9) c wp where separately.

Qualitatively this result can be obtained from Equa- The wave function (10) can be used for finding ~ ~ tion (2) with the angle θ taken equal to the diffraction- mean values of the biphoton momentum ~ k1 + k2 h i divergence angle of the pump θ = λ /w 1. As and, finally, the Lorentz-invariant mass mbiph. Because p p ≪ wp λp, the Lorentz-invariant mass of the pump per of the axial symmetry of the expression (10), mean one≫ photon is much smaller that the Lorentz-invariant transversal components of wave vectors are equal zero, mass of two counter-propagating photons (3). ~k⊥ + ~k⊥ = 0, whereas for the mean sum of longitu- h 1 2i 3 dinal momenta in the paraxial approximation we get Altogether Equations (11), (15) and (17) give the follow- ing expression for the mean momentum of biphoton pairs k + k = h z 1 z 2i (multiplied by the c) ωp c ~ ~ 2 ~ ~ 2 (k⊥ 1 + k⊥ 2) + (k⊥ 1 k⊥ 2) = ~ c − 2ωp h i h − i c(pz 1 + pz 2) = ωp c δp, (18) h i h i − ωp c 2 2 ~q+ + ~q− , (11) where c − 2ωp h i h i   ~ c 1 4πno πL where ~q± = ~k⊥ ~k⊥ . 1 2 δp = 2 + ln . (19) ± 2ωp w Lλp 2noλp In terms of two 2D vectors ~q+ and ~q− averagings are  p   1 1 − 2 understood as ... = 4 d~q+d~q Ψ (...), with 4 be- h i | | As always λp wp,L , in any case cδp ~ωp and, ing the transition Jacobian from variables ~k⊥ , ~k⊥ to R 1 2 2 ~ ≪2 { } 2 ≪ 2 as εbiph ( ωp) , the difference εbiph c(pz 1 + pz 2) q+, q−, and with the squared wave function (10) taking ≡ ~ − h i the form equals approximately 2 ωpcδp. In this approximation the biphoton Lorentz-invariant mass per one pair of 2 2 2 2 2 Lλp 2 SPDC photons takes the form Ψ = N exp ~q+ wp sinc ~q− . (12) | | − 8πno   1/2  ~ 1 4πno πL In this expression terms depending on ~q and ~q− are m = + ln . (20) + biph c w2 Lλ 2n λ factorized, which is very convenient for averagings.  p p  o p  The first step of calculations consists in finding the normalization factor from the condition 3. ENTANGLEMENT 1 2 d~q+ d~q− Ψ =1. (13) 4 | | ~ ~ Z Z Transverse momenta of SPDC photons k⊥1 and k⊥2 are 2D continuous variables. As known [9], the degree Both integrals over ~q+ and ~q− in (12) are easily cal- 2 of entanglement in such variables can be evaluated by culated separately to give, correspondingly, π/wp and 2 ∞ the Schmidt entanglement parameter K defined as the 8π no dx sinc2(x) with the last integral over x equal Lλp 0 inverse trace of the squared reduced density matrix. Cal- to π/2. Combined together, these results give finally R culation of this parameter is rather simple and straight- forward for the so called double-Gaussian wave functions, wp Lλp having the form of a product of two Gaussian functions, N = 2 . (14) π r no one of which depends on the sum of variables and the other one - on their difference. For non-double-Gaussian The next steps are finding ~q 2 and ~q 2 in Equation + − wave functions functions calculation of the parameter K (11). The first of these two quantitiesh i h is determinedi by is a problem, and even more difficult problem is its direct the integral of the exponential function in Equation (12) experimental measurement. giving Another and much easier measurable entanglement pa- 1 rameter R was suggested for the first time in the work ~q 2 = . (15) + 2 [10]. This parameter is defined mathematically as the h i wp ratio of widths of the unconditional to conditional prob- 2 As for the second term, ~q− , it can be reduced in a ability densities depending on the variable of one of two similar way to the followingh integrali form particles (e.g., photons). Defined in this way the param- ∞ eter R was found to coincide exactly with the Schmidt 2 8πno 2 ~q− = x sinc (x) dx, (16) parameter K [11], and in other cases to be on the order of h i Lλp 0 K [12]. In experiment one has to split the original beam Z of biphotons for two channels and perform two kinds of where x is the integration variable equal to the argument measurements. At first, photons can by counted by a of the sinc-function in Equation (12). In principle, for- single scanning detector in only one of two channels to mally, the remaining integral in Equation (16) diverges plot the single-particle distribution and to find its single- logarithmically at x . But this divergence is re- particle width. In the second series of measurements one lated to the used above→ paraxial ∞ approximation. In fact, ωp has to use two detectors located in different channels, one this approximation is valid only as long as k⊥ , 1, 2 2c of them scanning and another one kept at a constant po- ωp Lλp ωp 2 | | ≪ q− , and x xmax = ( ) . For this rea- sition, and only coinciding signals have to be registered at | | ≪ c ≪ 8πno c son, the integral in Equation (16) can be estimated as the computer. In this way one gets the coincidence (or xmax 2 1 conditional) distribution and measures its width. The 1 sin x dx/x = 2 ln(xmax) to give ratio of found widths is the parameter R. R 2 4πno πL It’s known also that if the wave function has the form ~q− = ln . (17) h i Lλ 2n λ of a product of two terms, one of which depends on p  o p  4 the sum and the other one - on the difference of vari- to each other: the Lorentz-invariant mass is proportional ables and if the widths of these distributions are a and b, to the Schmidt entanglement parameter K. then the coincidence and single-particle widths are, cor- It’s true, however, that this connection is not universal. respondingly, min a, b and max a, and, consequently, { } { } In the case of very strong focusing of the pump L Ld = R = max a, b / min a, b . Applied to the wave function 2 ≫ { } { } wp/λp, opposite to that considered above, the degree of of the form (10) and variables k1,x, k2,x, these definitions entanglement is high too, K L/Ld 1, but the give a =1/wp, b =2π no/Lλp and Lorentz-invariant mass (20) appears∼ to be≫ independent p p of the degree of entanglement and determined only by (s) ∆k max 1/wp , 2π no/Lλp the inverse waist of the pump, mbiph = ~/c wp. K R = 1,x . (21) ∼ (c) ∼ n o ∆k1,x min 1/wp , 2π pno/Lλp n p o This derivation and the derived expression confirm the 4. CONCLUSION known result [7, 8] that the degree of entanglement is high either if wp Lλp or if wp Lλp. In the second of these two≪ limiting cases the≫ crystal is assumed to be The main result of the presented analysis concerns short comparedp to the diffractionp length of the pump, demonstration that there are conditions when biphotons 2 states are highly entangled and this high entanglement is L wp/λp Ld. Under this condition ≪ ≡ directly related to the relatively high Lorentz-invariant 2πwp√no mass, mbiph K 1 (23). We believe that this is K Rshort L 1. (22) a fundamentally∝ important≫ new knowledge. Though we ∼ ∼ Lλ ≫ p realize that it may be difficult to imagine any ways of With this expression for the entanglementp parameters K measuring the Lorentz-invariant mass of biphotons di- and R, the derived above formula (20) for the Lorentz- rectly and independently of measuring coincidence and invariant mass of biphoton pairs can be rewritten as single-particle widths of momentum-distributions, find- ing the entanglement parameter R K and then using ∼ ~ K2 πL 1/2 if for finding mbiph from Equation (23). m = 1+ ln Acknowledgement biph c w π 2n λ p   o p  ~K 1 πL ~ ln . (23) The work of S.V. Vintskevich was supported by the ≈ 2c wp π 2noλp ≫ 2c wp Foundation for the Advancement of Theoretical Physics s   BASIS (PhD Student Grant No. 17-15-603-1). D.A. 2 This result shows that in the case L wp/λp Ld both Grigoriev and M.V. Fedorov acknowledge support of the the Lorentz-invariant mass of the biphoton≪ state≡ (10) and Russian Foundation for Basic Researches, grant RFBR its degree of entanglement are high and they are related 18-02-00634.

[1] L. B. Okun. Energy and Mass in Relativity Theory. [8] M. V. Fedorov, M. A. Efremov, P. A. Volkov, World Scientific, Singapore, 2009. E. V. Moreva, S. S. Straupe, and S. P. Kulik. Phys. [2] L. B. Okun. Am. J. Phys., 77:430, 2009. Rev. A, 77:032336, 2008. [3] D. Giovannini, J. Romero, V. Poto´cek, G. Ferenczi, [9] M. V. Fedorov and N. I. Miklin. Contemporary Physics, F. Speirits, S. M. Barnett, D. Faccio, and M. J. Pad- 55:94, 2014. gett. Science, 347 (6224):857, 2017. [10] M. V. Fedorov, M. A. Efremov, A. E. Kazakov, [4] M. V. Fedorov andS. V. Vintskevich. Laser Phys., K. W. Chan, C. K. Law, and J. H. Eberly. Phys. Rev. A, 27:036202, 2017. 69:052117, 2004. [5] M. V. Fedorov, S. V. Vintskevich, and D. A. Grigoriev. [11] M. V. Fedorov, M. A. Efremov, P. A. Volkov, and Europhys. Lett., 117:64001, 2017. J. H. Eberly. J. Phys. B: At. Mol. Opt. Phys., 39:S467, [6] M. V. Fedorov. Phys.Rev. A, 93:033830, 2016. 2006. [7] M. V. Fedorov, M. A. Efremov, P. A. Volkov, [12] W. Mauerer and C. Silberhorn. AIP Conf. Proc., E. V. Moreva, S. S. Straupe, and S. P. Kulik. Phys. 1110:220, 2009. Rev. Lett. , 99:063901, 2007.