Nonlocal realism tests and quantum state tomography in Sagnac-based type-II polarization-entanglement SPDC-source.

We have experimentally created a robust, ultrabright and phase-stable polarization-entangled state close to maximally entangled Bell-state with %98-fidelity using the type-II spontaneous parametric down-conversion (SPDC) process in periodically-poled KTiOPO4 (PPKTP) collinear crystal inside a Sagnac interferometer (SI). Bell inequality measurement, Freedman's test, as the different versions of CHSH inequality, and also visibility test which all can be seen as the nonlocal realism tests, imply that our created entangled state shows a strong violation from the classical physics or any hidden-variable theory. We have obtained very reliable and very strong Bell violation as S = 2.78 ± 0.01 with high brightness V HV = % ( 99.969 ± 0.003 ) and V DA = % ( 96.751 ± 0.002 ) and very strong violation due to Freedman test as δ F = 0.01715 ± 0.00001 . Furthermore, using the tomographic reconstruction of quantum states together a maximum-likelihood-technique (MLT) as the numerical optimization, we obtain the physical non-negative definite density operator which shows the nonseparability and entanglement of our prepared state. By having the maximum likelihood density operator, we calculate some important entanglement-measures and entanglement entropies. The Sagnac configuration provides bidirectional crystal pumping yields to high-rate entanglement source which is very applicable in quantum communication, sensing and metrology as well as quantum information protocols, and has potential to be used in quantum illumination-based LIDAR and free-space quantum key distribution (QKD).

Introduction

Generation and characterization of polarization-entangled Bell-states are very important from fundamental and applied point of view in the different contexts of the photonic-based quantum technologies such as quantum metrology, realizing protocols of quantum computing, and also quantum communication (for a review, see [1], [2], [3], [4]). To our knowledge, the cleanest, easiest, and most inexpensive and accessible entanglement source, notably the polarization entanglement source, can be realized via the process of spontaneous parametric down-conversion (SPDC) in bulk nonlinear crystals (NLCs) such as BBO, PPKTP or BaBO. However, there are different methods to generate the polarization entangled photon pairs, such as the non-post selective two-photon resonant excitation scheme on a single semiconductor quantum dot [5] and the high-harmonic generation method plus a mixed co-propagating elliptically polarized waves in an isotropic media [6]. Moreover, very recently, it has been proposed [7] an interesting method to generate the long-distant polarization-spatial-mode hyper-entanglement for quantum communication using the PPKTP crystal over 11 km in multicore fiber. The authors have shown that the total entanglement purification efficiency can be estimated about 3-order more than the experiment using two pairs of entangled states with SPDC. It should be reminded that the probabilistic nature of their generation process in SPDC or four wave mixing leads to creation of zero or multiple photon pairs yielding a Poissonian distribution which limits their applications such as in complex algorithms where many qubits and gate operations are required.

SPDC is the well-known nonlinear optical process, in which a classical pump laser beam at higher frequency is incident onto an optical nonlinear material which under the so-called phase-matching (PM) conditions, i.e., energy-momentum conservation, a twin-photons beam, i.e., a time-energy entangled [8] pair of signal-idler photons at lower frequencies, is generated out of quantum vacuum [1], [2]. One of the advantages of nonclassical sources based on the SPDC in bulk NLCs is that they are more inexpensive, user-friendly, and work at room-temperature, which means that need no cooling, and have capability to be scalable in order to be commercialized. During these three decades, the SPDC in NLCs has been a variety range of applications in quantum metrology, quantum communication, quantum computing, quantum information, and quantum thermodynamics (for a review, see Ref. [4]). Among these, one can remark, for example, EPR realization [9], the entanglement generation [10], [11], [12], [13], [14], [15], quantum state teleportation [16], quantum ellipsometry [17], [18], quantum illumination [19], [20], [21], [22], [23], [24], [25], [26], quantum spectroscopy [27], [28], [29], [30], squeezing generation [31], quantum imaging [32], [33], [34], [35], [36], quantum communication [37], [38], [39], [40], [41], [42], [43], [44], [45], [46], [47], [48], [49], [50], and nonlocal realism tests [9], [51], [52], [53], [54], [55].

The most important way to characterize the quantum state of the system is quantum state tomography (QST) which leads to reconstruction of the density matrices operator of the system (DMOS). Density matrice includes all information such as probabilities and coherences. Also, the evolution of the density operator enables us to evaluate the environmental effects on the quantum state of the system, and decoherence effects. It has been shown theoretically and experimentally [56], [57], [58] that using the theory of the measurement of qubits, which includes a tomographic reconstruction of the DMOS due to the linear set of 16 polarization-measurements together with the numerical optimization method, the so-called maximum-likelihood-technique (MLT), one can obtain and reconstruct the non-negative definite density operator. By having the density operator, all quantum coherence properties and entanglement entropies can be measured and calculated. Moreover, very recently, it has been theoretically developed and experimentally demonstrated a new method, resource- and computationally efficient, for quantum state tomography in Fock basis via Wigner-function reconstruction and semidefinite programming [59]. It has been shown that obtained density operator from this method is robust against the noise of measurement and relies on no approximate state displacements, and requires all physical properties. Furthermore, a new novel method has been developed and demonstrated for scalable on-chip QST [60], which is based on expanding a multi-photon state to larger dimensionality. It leads to scale linearly with the number of qubits and provides a tomographically complete set of data with no reconfigurability.

Optical interferometry is the most popular technique in optical metrology and physical optics which has been applied in a variety range of applications [61], [62], [63], [64], [65], [66], [67], [68]. Among the different configurations of interferometers, in recent decades, the Sagnac-based sources because of their stability and high rate have been used in quantum optical experiments in a variety range of application such as quantum key distribution (QKD) [43], [69] and nonclassical radiation sources [70], [71], [72], [73]. Here, it worths to remind that in a Sagnac interferometer (SI), a coherent light beam is split on two branches, and thus, the two split beams pass the common-path but in opposite directions. When the two beams return to the entrance port of SI, if their wave-functions coincide, thus, coherence condition is satisfied and they interfere.

In this paper, we are motivated by the above-mentioned research to investigate the nonlocal realism tests such as Bell's inequality test and Freedman's test as well as visibility test in an ultrabright, robust and high-rate polarization-entangled photons source via the type-II SPDC in a bulk PPKTP collinear crystal inside the common-path SI. Sagnac configuration enables us to Bidirectional pumping the PPKTP crystal inside which yields to robust and ultrabright polarization-entanglement source with high pair photon flux. After preparing and proving the entangled Bell-state, via the quantum state tomography assisted maximum-likelihood numerical optimization technique, we obtain the physical density matrix of the system, which results in entanglement entropies and entanglement measures.

Experimental results

Experimental setup

Fig. 1(a) shows the schematic of experimental setup for type-II Bell-polarization-entanglement generation via the SPDC in PPKTP crystal inside the common-path SI. A 48 mW-beam of a diode-laser passes through an optical isolator, and then, is reflected by mirror M1. In order to generate polarization entanglement in Sagnac configuration, one should bidirectionally pump the collinear PPKTP crystal with both vertical (V) and horizontal (H) polarization. So, to this, one should polarize the pump laser beam at via an HWP. Then, -polarized laser beam is reflected by M2 toward PBS, i.e., the entrance port of the Sagnac. In the reflected/transmitted port of PBS, the V/H-polarization part of the laser beam can pass. Both V and H polarized beam are reflected by M4 and M3 in Sagnac, and then are focused on the center of PPKTP crystal by concave lens with focal length 30 cm. The crystal is placed on a heater which is fixed on a high-precision translation stage with fine-tilt capability. The crystal temperature is fixed at 30 °C which can be controlled with 0.1 °C precision. Note that to have a maximum rate and efficiency, it is necessary to focus the beam on center of the crystal which is achievable by tuning the lens and the crystal position. Each polarized beam of pump in Sagnac can generate a signal-idler pair photons at wavelength 810 nm via the type-II SPDC process in the PPKTP crystal [see Fig. 1(a)].

Figure 1

(a) Schematic of SI to create the ultrabright and robust type-II polarization-entangled Bell-state at wavelength 810 nm using the SPDC process in bulk PPKTP NLC. LD: 405 nm diode-laser, I: optical isolator for 405 nm, M1(M2): high-reflective mirror, HWP: half-wave plate [note that HWP in Sagnac arm is a dual wavelengths HWP at 405/810 nm], QWP: quarter-wave plate, L: concave lens with focal length 30 cm, M3(M4): high-reflective broadband mirror at 405 nm and 810 nm, PBS: polarizing beam-splitter [Note that in detectors arm PBS is set for 810 nm while in entrance port of Sagnac a dual wavelengths PBS is used], F: 780 nm-bandpass and 810 nm-narrowband filters, a SPCM: Excelitas single-photon counting module, MMF/SMF multi-/single-mode optical fiber, C: coupler lens module. Here, we used the quTools coincidence time-tagger with 81 ps resolution. PPKTP: AR coted 1 × 2 × 25 mm periodically-polled-KTP collinear crystal with Λ0 = 10.025 μm which is kept at temperature 30 °C by a temperature controller with 0.1 °C resolution. The coincidence time window and acquisition integration time are, respectively, 5 ns and 400 ms. PPKTP is cut for collinear type-II SPDC QPM condition. (b) Side-view of the experimental setup.

To coincide two wavefunctions and in order to generate the maximally polarization entangled Bell-state one should tune and rotate the QWP, acts as a compensator, at proper angle φ around the z-axis. In our experiment, to achieve the best entangled state we have set the QWP at angle . As shown in Fig. 1(a), after generation of pair photons in Sagnac, pair photons and , respectively, are reflected by the back-reflect port of PBS and dichroic mirror (DM), and then, pass through a HWP plus PBS, which act as a polarization-measurement box, and after passing through a 780 nm bandpass filter and a 810 nm narrowband filter with 12 nm FWHM are fed via multi-mode fibers (MMF) or single-mode fibers (SMF) to Excelitas single-photon counting modules (SPCMs) with 200 cps and 400 cps dark noises. The coincidence count (CC) rates are measured by quTools coincidence time-tagger with 81 ps resolution. The experimental setup is shown in Fig. 1(b). By checking the coincidence count (CC) rate at different angles of HWP in front of the SPCMs, one can tune the compensator in order to have maximum entanglement and generate an entangled state near to maximally type-II Bell-state.

Quasi-phase matching (QPM) condition

Periodic polling technique [74] is a technique based on micro-lithography for fabrication of a metallic thin film nano-layers with specified alternate wavevector orientation in a birefringent optical material in bulk or integrated optics [1] [see Fig. 2]. The period Λ which spaces regularly the domains is multiple of the desired operation wavelength. The periodically poled crystals, for example potassium titanyl phosphate (KTP), lithium niobate (LiNbO3), and lithium tantalate (LiTaO3), are frequently used as an efficient nonlinear optical materials at second-harmonic generation (SHG) and SPDC than the other types of crystals without periodic polling. Thermal pulsing, pulsed electric field and electron bombardment, or other techniques are usually used to reposition the atoms in the lattice, creating oriented domains which can be achieved either during the growth of the crystal, or subsequently. This structure is usually designed to achieve the so-called quasi-phase-matching (QPM) condition, a generalized energy-momentum conservation, in the material to extend the wavelength regime to μm-wavelengths. In QPM condition, one should simultaneously satisfy relations where , and , and is the refractive index which depends on the wavelength and temperature, [75]. As is seen, the last term in Eq. (2) is the wavevector due to the periodic polling which its temperature dependency provides the long range of wavelengths for pair generation via the SPDC process. Note that here the crystal length depends on the temperature (for more details, see Ref. [75]).

Figure 2

Schematic of a periodically-poled KTiOPO4 crystal (PPKTP) with length L. In the z-direction, the wave-vector is periodically modulated via exerting a high voltage on the coated metal on crystal with periodic length Λ at temperature 25 °C.

Let us remind that in a type-II () periodically poled PPKTP source, as a biphoton source, an extraordinary (e) polarized pump photon is down-converted to two ordinary (o) and extraordinary photons, the signal-idler biphoton, under the QPM condition [1]. This type of NLC has potential to generate the type-II polarization-entangled Bell-state (for more details about the state of SPDC, see Refs. [2], [12], [13]).

Nonlocality tests

In this subsection by verifying some well-known tests, we will show that the generated entangled state violates the classical physics prediction. In the proceeding, we will present our experimental results of visibility test, Bell parameter measurement and Freedman test as the nonlocal realism tests for the generated state via the SPDC in PPKTP crystal in Sagnac configuration.

Bell inequality measurement and visibility

We are interested in measuring CHSH (Clauser, Horne, Shimony, and Holt) inequality [76] which is slightly different from the original Bell's inequality. In 1964, Bell showed that all hidden-variable theories (HVTs) obey his inequality, while quantum mechanics, for example entangled state, violates [12], [77], [78].

Bell inequality constrains the degree of polarization correlation under measurements at different polarizer angles . Bell parameter S which is criterion for Bell inequality includes measure which incorporates all possible polarization CC-measurement outcomes, and varies from −1 (when the polarizations always agree) to +1 (when the polarizations always disagree). It is given by [76], [77] where (for more details and for a short review, see Refs. ([12], [77])). Thus, S parameter is given by [76], [77] where stand for the different polarizer angles. To obtain S, one needs to 16 CC-measurement. The most important feature of S-parameter is that it is theory-independent and has no clear physical meaning. It can be proved that for any HVTs and arbitrary angles, Bell parameter satisfies inequality which is equivalent to the visibility of coincidence probability in both H-V and D-A basis [54], [76], [77]. Surprisingly, it can be proved that quantum mechanics can violate this inequality by considering entangled-state. Also, it can be maximized for the maximally entangled Bell-states, such state given in Eq. (1). By choosing angles , , (V), and ; we obtain [12], [77] for any maximally entangled Bell-state, while other quantum states give lower values such that . Interestingly, for these angles, it can be proved that HVT gives the maximum S-value . The Bell inequality emphasizes that no realistic, local and complete theory, i.e., the classical physics, in the EPR context will ever agree with quantum mechanics. In the ideal case, for the Bell state (1), the coincidence rate can be theoretically calculates as . In experiment, if (), it shows the violation the Bell inequality and approves that the nature does not agree the HVTs or classical physics.

After our best tuning the state near to the maximally entangled Bell state, we start to measure the CC rate to obtain the visibility and S-parameters which will be given in the following.

In Fig. 3, we have plotted the behavior of the experimentally measured CC rate when the polarizer angle of channel A(1) [or half of the HWP angle] is fixed at angles , , and for the different values of the polarizer angles of channel B(2) which agrees with the theoretical quantum prediction. Table 1 shows the experimental results for visibility measurement. Using the experimentally measured data, we find and with 5-digit accuracy for the visibility in H/V and 45/135 (D/A) bases, respectively. Our measured visibility shows very strong and reliable violation from the classical physics prediction (). Our measured visibility implies on that our polarization-entanglement source is a high-brightness nonclassical source.

Figure 3

Experimentally measured CC rate of channel A (1) vs. polarizer angle (the half HWP-angle) of channel B (2). (a) A is fixed at angles 0∘ and 90∘. (b) A is fixed at angles 45∘ and 135∘. Blue and brown thick-dotted points; red and green thick-triangle are referred to polarizer angles A = 0∘, A = 45∘, A = 90∘, and A = 135∘, respectively. Here, HWP plus PBS in front of detectors act as a polarizer box such that θpol = 2θHWP.

Table 1

Results for experimental visibility measurement. We obtain and for the visibility in H/V and 45/135 bases, respectively.

θpolA	θpolB	Rc (Kcps)	ΔRc (Kcps)
0	0	0.5376	0.0361
0	90	14.7923	0.1848
90	0	14.7830	0.1834
90	90	0.2919	0.0120
45	45	0.2743	0.0287
45	135	15.3375	0.1781
135	45	0.2211	0.0292
135	135	14.6600	0.2496

To calculate the Bell parameter, one should measure 16 CC rates which are given in Table 2. The experimental measured rates for our prepared entangled state lead to Bell parameter , which shows very reliable and strong violation. This strong violation shows a high-degree of entanglement of the prepared state in our experiment which can be used as a strong nonclassical source to quantum sensing and metrology in order to enhance the signal-to-noise ratio and the sensitivity of measurement by exploiting the quantum entanglement.

Table 2

Experimental measured CC rates corresponding to the Bell inequality measurement. We have obtained S = 2.78 ± 0.01 which shows very reliable (with 78σstd) and strong violation from the HVTs or any classical prediction. R, R and R are, respectively, singles and coincidence count rates as a function of polarizer angles in channel A and B. Here, the integration time and coincidence time-window are, respectively, T = 0.4 s and τ = 5 ns. The accidental rates are calculated using τRR/T.

anglea [degree]		single count [Kcps]		coincidence rate [Kcps]
θpolA	θpolB	RA	RB	R_c	ΔRc
0	-22.5	207.6210	214.0125	1.7567	0.0724
0	22.5	203.7475	207.4312	2.1062	0.0534
0	67.5	200.1308	166.2038	10.8926	0.2150
0	112.5	200.9170	174.8037	10.4809	0.1840
-45	-22.5	178.0473	213.5483	2.3775	0.0671
-45	22.5	182.7763	207.4940	11.5540	0.1933
-45	67.5	179.1138	166.5318	10.3532	0.1431
-45	112.5	177.7340	174.5832	1.5611	0.0742
45	-22.5	175.5212	213.7574	11.4244	0.1974
45	22.5	179.6321	207.3254	2.2608	0.0876
45	67.5	177.1177	166.8515	1.9555	0.0765
45	112.5	175.7228	175.2128	11.4158	0.1852
90	-22.5	160.5896	212.8264	12.3135	0.2163
90	22.5	155.0147	206.7563	11.5542	0.1430
90	67.5	154.5420	167.3865	1.8846	0.0741
90	112.5	154.5125	174.6846	2.2487	0.0634

Here, angles are attributed to polarizer's angles that in our configuration are equivalent to .

Freedman's test

Beyond the standard CHSH Bell test which requires 16 measurements, the so-called Freedman's inequality [79], [80] requires only 3 coincidence measurements with very simpler mathematical calculations and experiment, more understandable from the philosophical aspects and simpler physical interpretation, notably for undergraduate students in quantum optics Laboratories.

Freedman's inequality can be quantified by Freedman parameter which is defined as (to see full details of derivation of Freedman's inequality, see Refs. [79], [80]) where is the coincidence count for polarizer angles a and b such that . The optimum values for angles are and [79], [80]. is the coincidence when both polarizers/analyzers are removed (note that in our experiment the combination of HWP and PBS in front of each detector plays the role of polarizer). It can be proved that for any realistic and local theory such as classical physics, the Freedman parameter is always nonpositive, i.e., . If in an experiment Freedman parameters reliably becomes positive, , thus, one says the HVTs or classical physics cannot interpret it. For example, it can be shown that an entangled state violates the Freedman's inequality and becomes positive which shows it cannot be explained by any HVTs or classical physics.

Let us consider Bell state in Eq. (1), it can be easily shown that the coincidence count can be obtained as [2], [12], [77], [79], [80] where () is the transmittance of each arm when its polarizer box (HWP plus PBS) is present. Fig. 4 shows the normalized experimental measured coincidence count rate vs. difference of the polarizers angle for two cases: one angle is fixed at () and the other one changes. By fitting the theoretical expression (7) to the experimental coincidence pints, we obtain the average transmittance as .

Figure 4

Normalized coincidence count rate versus Φ for (a) A-polarizer is fixed at 0∘, and (b) A-polarizer is fixed at 45∘, and B-polarizer changes such that Φ varies zero to 90∘. The red-dotted and green-triangle points obtained from our experiment while the blue- and red-solid lines show the theoretical fitted curve from Eq. (7) to the experimental points.

To calculate the Freedman's parameter , we measured the CC rates for different values of Φ [see Table 3]. Our experimental data provide Freedman parameter with 5-digit accuracy as when and and ; and when and and . As theory predicted for entangled state for and , Freedman parameter becomes positive which shows a strong violation. It should be mentioned that to calculate the error in Freedman parameter, we have used relation [80]

Table 3

Results for measuring the Freedman parameter δF. The rows in bold are corresponding to the maximum value of violation from the Freedman inequality. The optimized values corresponding to violation from Freedman's inequality are highlighted in bold.

angle [degree]			rates [cps]
θpolA	θpolB	Φab	RA	RB	Rc
0	0	0	202794.8	223713.2	15.50
0	11.25	11.25	203103.3	219832	741.09
0	22.5	22.5	202627.5	212547.6	2412.05
0	33.75	33.75	201121.3	201762.6	4984.87
0	45	45	200279	191063.5	7729.67
0	56.25	56.25	202510.3	179757.9	10767.86
0	67.5	67.5	202614.5	171022.5	13122.65
0	78.75	78.75	203349.2	166741.3	14627.07
0	90	90	204432.5	165988.7	15241.43
45	45	0	175596.8	189897.2	209.28
45	56.25	11.25	173757.4	179289.1	842.59
45	67.5	22.5	172162	171451.6	2645.83
45	78.75	33.75	172518.4	166549.4	4796.4
45	90	45	171925.9	167041	7764.02
45	101.25	56.25	173068.9	170308	10662.76
45	112.5	67.5	172572.9	178028.4	12819.16
45	123.75	78.75	172656.7	189516.9	14733.98
45	135	90	174888.5	199757.1	15405.61

Quantum state tomography

Density operator reconstruction

Now, we are going to calculate the density operator of the prepared entangled state to see how much our state is close to the maximally entangled Bell-state .

To reconstruct the physical nonnegative density operator, we follow the generalized method of qubit measurement in Ref. [56]. The introduced method in Ref. [56] combines the experimentally measured different polarization-projections together a numerical optimization, the-so-called maximum-likelihood-technique (MLT), to reconstruct the physical density matrix which requires constraints (i) Hermiticity ; (ii) non-negative semidefinite eigenvalue in the interval [0,1] (); (iii) ; and also the important condition (iv) . In experiment, to provide the different projections, one should use a QWP before the polarizers (HWP-plus-PBS) in front of each detector. Let us remind why only projection-measurement cannot provide the physical nonnegative density matrix. Usually in Lab due to the statistical errors and inaccuracies, experimental measured values are different from their expected theoretical ones, and thus, we cannot reconstruct the physical density operator only by experimental measuring and it needs to numerical optimization to satisfy all physical constraints. In the other words, the only tomographic measurement without optimization provides an unphysical density matrix. Using the Jones matrix and our explained experimental configuration, the projection state can be defined as [56] It is equivalent to the projection measurement operator with . Then, the projected state of each qubit, , can be written as [56] with where are, respectively, the angle of the fast-axis of HWP and QWP concerning the vertical axis-polarization.

Fig. 5 shows the real and imaginary parts of the maximum-likelihood physical density matrix obtained from the experimental tomographic data in Table 4. It shows that our prepared state is close to the maximally Bell state . After numerical optimization using the measured tomographic data, the maximum-likelihood density operator matrix is obtained as with eigenvalues () corresponding to .

Figure 5

Graphical representation of the real (top panel) and imaginary (bottom panel) parts of the maximum-likelihood physical density matrix obtained from the experimental tomographic data given in Table 4.

Table 4

Experimental measurement data for the tomographic analysis states to reconstruct the physical density matrix using the MLT. Here, we have used the notation , and .

projection modea		wave-plates angle		single count rate [cps]				coincidence count rate [cps]
ν	AB	(hA,qA)	(hB,qB)	RA	ΔRA	RB	ΔRB	R_c	ΔRA
1	HH	(45,0)	(45,0)	147948.9	1105.93	153163.1	480.21	393.90	15.28
2	HV	(45,0)	(0,0)	147634.8	709.19	166312.8	609.79	20311.20	154.15
3	VV	(0,0)	(0,0)	192185.9	1093.16	166866.8	935.68	431.10	40.73
4	VH	(0,0)	(45,0)	191927.2	791.79	167524.8	1061.66	20020.70	232.02
5	RH	(22.5,0)	(45,0)	167385.3	864.70	168598.9	678.13	11489.00	184.40
6	RV	(22.5,0)	(0,0)	167324.6	774.44	166527.4	797.42	9051.50	184.44
7	DV	(22.5,45)	(0,0)	168523.2	828.53	167088.7	871.87	11123.50	204.01
8	DH	(22.5,45)	(45,0)	168084.6	1171.37	168071.7	591.91	8941.10	193.73
9	DR	(22.5,45)	(22.5,0)	165216.2	633.79	164677.3	567.47	7956.20	228.51
10	DD	(22.5,45)	(22.5,45)	170360	624.23	165050.3	776.69	716.90	56.21
11	RD	(22.5,0)	(22.5,45)	168862.9	1573.20	165050.3	677.33	9472.30	151.59
12	HD	(45,0)	(22.5,45)	142171.6	890.40	165763	1349.81	10282.60	137.01
13	VD	(0,0)	(22.5,45)	190450.1	1090.98	166000.8	1074.42	8346.60	110.84
14	VL	(0,0)	(22.5,90)	192181.3	1090.04	168922.5	783.59	12727.60	175.61
15	HL	(45,0)	(22.5,90)	146866.6	803.78	169126.4	784.87	7103.00	147.25
16	RL	(22.5,0)	(22.5,90)	166698.1	1194.72	168640.5	1260.55	18817.20	238.57

It must be mentioned that when the measurements are taken, only one wave-plate angle has to be changed between projection measurements. That is reason to choose this set of projection ν-arrangement.

To quantify and see how much the prepared state is close to the desired entangled state, one should calculate the fidelity in the proceeding.

Fidelity

Fidelity is a measure of the closeness of two quantum states which expresses the probability that one quantum state is similar and exactly close to the other quantum state, but, it is not a metric on the density operator space. Mathematically, the fidelity for two arbitrary density operators is defined as It can be easily shown that the fidelity is a symmetric quantity, and thus, . For two pure states (), the fidelity is simplified as , which is exactly the definition of the inner product in wavefunction space or vector space.

By considering and , we obtain between our prepared entangled state and the maximally entangled Bell state of Eq. (1) which shows a very good fidelity to prepare entangled state close to the desired entangled Bell state . In the following, we calculate some important entanglement entropies and measures using the reconstructed density matrix (13).

Entanglement entropy

In this subsection, we are going to calculate Concurrence, the entanglement of formation, tangle, logarithmic negativity, and different entanglement entropies such as linear entropy, Von-Neumann entropy, and Renyi 2-entropy which all can be derived from the density operator.

The Von Neumann entropy which quantifies the purity of a quantum state is given by [56] which for the pure state becomes zero. Our obtained Von Neumann entropy , shows our prepared state is the mixture-state. To quantify the degree of mixture of a quantum state, one can use the linear-entropy which is given by [56] that lies between zero and 1 (). We have obtained which shows our state is not a pure state.

To compute the quantum coherence properties of a mixed quantum state for a two-qubit system, one can calculate the concurrence, entanglement of formation, and tangle, which are equivalent measures of the entanglement of a mixed state. By considering the spin-flip matrix as we can define the non-Hermitian matrix with the left and right eigenstates, and eigenvalues with assumption , and thus, the concurrence is defined as [56] where and . We obtain the concurrence as . The tangle and the entanglement of formation are, respectively, defined as where is a monotonically increasing function. We find that and . Renyi 2-entropy, which quantifies the purity of the subsystem of a pure state, is another measure to quantify the entanglement, and is defined for the density operator of a subsystem which belongs to a bipartite pure state as [81] where the reducible density operators are given by The high value of Renyi 2-entropy shows a high degree of entanglement or low purity for the subsystem. In our case, we find , which implies on impurity of our state. Finally, the Logarithmic negativity, as upper bound to the distillable entanglement, which is obtained from the Peres-Horodecki criterion for separability [82], is defined as where is partial transpose with respect to subsystem A. It should be noted that, the negativity is independent of the transposed party because . Therefore, if the density matrix is separable, then, all eigenvalues of are non-negative. Otherwise, if the eigenvalues are negative; thus, is entangled. In our case, we find which shows the high degree of entanglement of our prepared state.

Summary, conclusion and outlooks

In this work, we experimentally prepared, measured and characterized high brightness, robust and stable type-II polarization-entangled state very close to the maximally entangled Bell-state with %98 fidelity via SPDC process in PPKTP NLC inside the SI. We have measured Bell and Freedman parameters as well as the visibility which are criteria for nonlocality test. They all showed the prepared state has high degree of entanglement. By reconstructing the physical density operator using the tomographic projection measurements together with the numerical optimization technique, we calculated some entanglement entropies which implied on the nonseparablity of the prepared state.

As outlooks, we are going to use this realized high-rate phase-stable entanglement source to implement the quantum LIDAR based on quantum illumination and to implement free-space QKD experiment.

Declarations

Author contribution statement

Ali Motazedifard: Proposed and developed the primary idea of the nonlocal tests in Sagnac configuration using PPKTP crystal; Conceived and designed the experiments; Contributed to build the experimental setup and performed the experiment; Calculated the theoretical analysis; Analyzed and interpreted the data; Contributed reagents, materials, analysis tools or data; Wrote the paper.

S.A. Madani: Contributed to build the experimental setup and performed the experiment; Did the numerical calculations; Analyzed the data; Contributed reagents, materials, analysis tools or data.

N.S. Vayaghan: Contributed to build the experimental setup.

J.J. Dashkasan: Contributed to build the experimental setup and drew Fig. 1(a).

Funding statement

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Data availability statement

Data included in article/supplementary material/referenced in article.

Declaration of interests statement

The authors declare no conflict of interest.

Additional information

No additional information is available for this paper.