Single plane minimal tomography of double slit qubits

TL;DR

This paper explores implementing minimal quantum tomography for double slit qubits using a single measurement plane, demonstrating that a tetrahedron measurement can be realized through specific geometric and propagation conditions.

Contribution

It shows how to perform tetrahedron measurements of double slit qubits on a single plane, simplifying quantum state tomography procedures.

Findings

Tetrahedron measurement can be implemented with Gaussian slits and free propagation.

A specific detection plane and slit geometry enable minimal tomography.

Practical measurement sets are proposed for experimental realization.

Abstract

The determination of the density matrix of an ensemble of identically prepared quantum systems by performing a series of measurements, known as quantum tomography, is minimal when the number of outcomes is minimal. The most accurate minimal quantum tomography of qubits, sometimes called a tetrahedron measurement, corresponds to projections over four states which can be represented on the Bloch sphere as the vertices of a regular tetrahedron. We investigate whether it is possible to implement the tetrahedron measurement of double slit qubits of light, using measurements performed on a single plane. Assuming Gaussian slits and free propagation, we demonstrate that a judicious choice of the detection plane and the double slit geometry allows the implementation of a tetrahedron measurement. Finally, we consider possible sets of values which could be used in actual experiments.

Tables

Table 1. Table 1: Tetrahedron tomography: some numerical solutions

$\zeta$	$w_1$	$w_2$	$w_3$	$w_3$
3.4678	-1.0287	-0.268044	0.268044	1.0287
6.08028	-0.943335	-0.367661	0.367661	0.943335
7.70501	-1.89747	0.0782496	0.332645	0.628523
8.5243	-1.30348	-0.153589	0.363175	0.805324
8.55362	-1.124	-0.111174	0.111174	1.124
10.6561	-0.738203	-0.487839	0.272693	1.140046

Table 2. Table 2: Optimal geometry for nearest SIC tetrahedron tomography of double slit qubits

$\lambda$	$a$	$2d_{opt}$	$z_0$	$x_3$	$x_4$
650 nm	100 $\mu$m	553 $\mu$m	33.53 cm	63 $\mu$m	242 $\mu$m
780 nm	62.5 $\mu$m	346 $\mu$m	10.9 cm	39 $\mu$m	151 $\mu$m
826 nm	60 $\mu$m	332 $\mu$m	9.5 cm	38 $\mu$m	145 $\mu$m
810 nm	40 $\mu$m	221 $\mu$m	3.5 cm	25 $\mu$m	97 $\mu$m

Equations

$$\frac{{\partial }^2\psi (x,z)}{\partial x^2}+2ik\frac{\partial \psi (x,z)}{\partial z}=0.$$

$$\psi (\xi ,\zeta )$$

$$=\int d{\xi }^{\prime }\frac{e^{i(\xi -{\xi }^{\prime }{)}^2/(2\zeta )}}{\sqrt{2\pi i\zeta }}\psi ({\xi }^{\prime },0).$$

$$\rho (\zeta =0)=\sum _{i,j=1,2}{\rho }_{ij}|{\psi }_i(\zeta =0)⟩⟨{\psi }_j(\zeta =0)|.$$

$$p(\xi ,\zeta )$$

$$=\sum _{i,j=1,2}{\psi }_i(\xi ,\zeta ){\psi }_j^{\ast }(\xi ,\zeta ){\rho }_{ij}.$$

$$p(\xi ,\zeta )=I(\xi ,\zeta )\mathrm{Tr}\left(\hat{\Pi }(\xi ,\zeta )\rho (\zeta =0)\right),$$

$$\hat{\Pi }(\xi ,\zeta )=\frac{1}{I(\xi ,\zeta )}\left(\begin{array}{cc}{\psi }_1(\xi ,\zeta ){\psi }_1^{\ast }(\xi ,\zeta )& {\psi }_1(\xi ,\zeta ){\psi }_2^{\ast }(\xi ,\zeta )\\ {\psi }_2(\xi ,\zeta ){\psi }_1^{\ast }(\xi ,\zeta )& {\psi }_2(\xi ,\zeta ){\psi }_2^{\ast }(\xi ,\zeta )\end{array}\right).$$

$${\psi }_k(\xi ,0)=\frac{1}{\sqrt[4]{\pi }}{e}^{-(\xi -(-1{)}^k\delta {)}^2/2},k=1,2.$$

$${\psi }_k(\xi ,\zeta )=\int d{\xi }^{\prime }\frac{e^{i(\xi -{\xi }^{\prime }{)}^2/(2\zeta )}}{\sqrt{2\pi i\zeta }}{\psi }_k({\xi }^{\prime },0)=\frac{e^{-\frac{1-i\zeta }{2(1+{\zeta }^2)}(\xi -(-1{)}^k\delta {)}^2}}{\sqrt{\sqrt{\pi }(1+i\zeta )}}.$$

$$I(\xi ,\zeta ;\delta )=\sum _{k=1,2}|{\psi }_k(\xi ,\zeta ){|}^2=\frac{2e^{-\frac{{\xi }^2+{\delta }^2}{1+{\zeta }^2}}}{\sqrt{\pi (1+{\zeta }^2)}}\cosh \left(\frac{2\delta \xi }{1+{\zeta }^2}\right).$$

$$\hat{\Pi }(\xi ,\zeta )=\frac{1}{\cosh w}\left(\begin{array}{cc}{e}^{-w}& e^{-iw\zeta }\\ e^{iw\zeta }& e^w\end{array}\right)=\frac{1}{2}\left(\begin{array}{cc}1-\tanh (w)& \text{sech}(w)e^{-iw\zeta }\\ \text{sech}(w)e^{iw\zeta }& 1+\tanh (w)\end{array}\right).$$

$$\mathbf{s}(\xi ,\zeta )=\text{sech}(w)\left(\cos (w\zeta ),\sin (w\zeta ),-\sinh (w)\right).$$

$${\mathbf{s}}_1$$

$${\mathbf{s}}_2$$

$${\mathbf{s}}_3$$

$${\mathbf{s}}_4$$

$$\delta =\frac{1}{2}\sqrt{\frac{(1+{\zeta }^2)(w_1^2-w_2^2)}{\ln \cosh w_1-\ln \cosh w_2}},$$

Full text

††journal: Optics Communications

1 Introduction

Quantum effects can be exploited to process information in different and sometimes more efficient ways [1] than those allowed by classical physics. For instance, quantum physics promises faster computers [2] and more secure communications [3]. Different systems [4] (photons, atoms, spins, and superconducting and nanomechanical structures) are being used to build quantum devices. Among the optical implementations of quantum information technologies, the most popular one uses photon polarization as a natural two-level system (qubit) [5]. Implementations based on the spatial and temporal [6] degrees of freedom of light have also been used. For example, orbital angular momentum eigenstates have been employed to define (spatial) quantum d-level systems (qudits) [7]. Spatial qudits can also be produced when photons are made to pass through an aperture with $d$ pixels [8] or with $d$ slits [9, 10]. The photon transverse position has been used to prepare, measure and control spatial qudit states [9, 11]. In particular, several methods to estimate the quantum state of slit qubits (and qudits) have been reported [10, 12, 13, 14, 15, 16].

Due to their common mathematical structure, classical waves and quantum mechanics share many physical effects, such as the Gouy geometrical phase [17], to mention just one. Recently, ideas originated in quantum estimation have been applied to classical contexts. For example, the quantum Fisher information has been used to show that suitable measurements allow the resolution of incoherent sources separated by distances which violate the Rayleigh criterion [18, 19]. The physical system that we consider in this work, a double-slit qubit, corresponds to the classical Young’s interference experiment. We present a proposal for minimal tomographic reconstruction of double-slit qubits employing only free propagation. Although our results have been presented in a quantum mechanics language, analogous results hold for the classical two-slit interference setup when we change photon detection by intensity measurement.

A brief account of quantum tomography, section 2, provides context to formulate the problem of state estimation for qubits defined by a two-slit setup (section 3). In section 4, we demonstrate how the minimal quantum tomography can be performed using measurements on a single plane. Finally, conclusions are drawn, and additional remarks, concerning possible experimental implementations, are made.

2 Quantum Tomography

Quantum tomography is an a posteriori process that allows a thorough description of the quantum state of an assembly of identically prepared systems, based on data obtained with measurement apparatuses [20]. The origin of quantum tomography of systems of continuous variables can be traced back to Pauli [21], who considered the problem of the reconstruction of the wavefunction of a spinless quantum particle, given its coordinate and momentum probability densities [22]. In general, the probability density and the probability current (not coordinate and momentum probability densities) allow the reconstruction of pure states [23]. In the case of mixed states, it is possible to reconstruct the Wigner quasiprobability function from the probability distributions along straight lines in phase space [24, 25]. Experimental state reconstruction in a diversity of quantum systems (including molecular vibrational modes, one-mode and two-mode states of light, trapped ions, and helium atoms) have been reported [20, 26].

Stokes [27] arguably is the father of tomographic methods for the reconstruction of systems with finite-dimension Hilbert spaces. However, the first systematic approach to state estimation is due to Fano [28], who introduced the notion of a quorum, a set of observables sufficient to determine the quantum state of a system. Any quantum tomography —spin tomography [29, 30, 31], for example— can be performed using different quora. The elements of a quorum are not necessarily associated with observables, but with positive semidefinite operators $\mathbb{P}_{m}$ which resolve the identity, $\mathbb{I}=\sum_{m=1}^{k}\mathbb{P}_{m}$. This set of operators, collectively known as a positive operator valued measure (POVM), describe generalized measurements.

When the statistics of a POVM can completely determine the quantum state of a system, it is said to be informationally complete (IC). An IC-POVM must contain at least $d^{2}$ operators $P_{m},$ to be able to estimate the $d^{2}-1$ reals parameters that uniquely determine the density matrix of a $d$-level system. Two prominent examples of IC measurements are mutually unbiased (MU) measurements and symmetric informationally complete (SIC) POVMs. Bases such that the angles between arbitrary pairs of elements of different bases are all equal are known as MU bases (MUBs) [32]; a set of $d+1$ MUBs of a system of $d$ levels is informationally complete [33]. When the Hilbert-Schmidt inner products between every pair of different operators of an IC-POVM are all equal, this POVM is a SIC-POVM [34]. Despite the lack of a formal proof, it is widely believed that SIC-POVMs exist in any Hilbert space of finite dimension [35]. The SIC-POVM for qubits is also known as a tetrahedron measurement. It comprises four subnormalized projectors over pure states. In the Bloch sphere, each pure state is associated with a unit-length Bloch vector. The tips of the four Bloch vectors which characterize the SIC-POVM are the vertices of a regular tetrahedron [36].

The determination of the quality of a given tomographic method is an important but complex problem. The notion of optimality of a tomographic method depends on the assumptions made about the experimental setup (individual, collective, fixed or adaptative measurements), the reconstruction method employed (linear inversion, maximum likelihood estimation, etc), the particular quantifier of accuracy or efficiency (e.g., minimum squared error, maximum fidelity), and on how the average over the whole state of states is performed. Under particular assumptions, it has been shown that SIC-POVMs are optimal among all minimal (those with the minimal number of outcomes) IC measurements [36, 37], while MU measurements are optimal among all choices of IC projective measurements [38]. Moreover, several figures of merit and assumptions show that (optimal) MU measurements are more accurate than the corresponding SIC-POVMs [39, 40, 38, 41]. Perhaps the most fundamental of these figures of merit is the quantum tomographic transfer function (the trace of the inverse of the Fisher matrix, averaged over pure states using the Haar measure), which gives the average optimal tomographic accuracy per sampling event for all unbiased state estimators, in the limit of a large number of sampling events [41].

Numerous studies deal with the implementation of tomographic schemes in experiments involving photon polarization, including several schemes realizing this tetrahedron measurement [42, 43, 44]. In contrast, fewer papers have been devoted to the tomographic reconstruction of double slit qubits. For example, a MUB approach was employed to measure the state of two double slit qubits [12, 13]. In this case, detectors placed in the near and far-field of the slits, aided by double slit spatial filtering, allow the simultaneous measurement of the three Pauli operators of each qubit. It was subsequently discovered that spatial filtering is not essential, because Pauli operators can be measured using a lens and “point” detectors in the image and focal planes. It was also recognized that, for a fixed detection-plane to slit-plane distance, the measurement of the interference pattern corresponds to a continuous POVM of a single qubit; therefore, the elements of the density operator can be obtained from the interference pattern [15]. A spatial light modulator, which can be used to control amplitudes and phases, was the key device to implement the minimal SIC tomography proposed in [36] to reconstruct double slit qubits [16].

Though optimal MU measurements are more accurate, SIC-POVMs are simpler. Due to its simplicity, we consider the implementation of the tetrahedron measurement of double slit qubits of light in this paper. We show that a minimal SIC-POVM tomography of double slit qubits can be implemented by measures on a single plane, using only free propagation (without resorting to lenses, spatial light modulators, or other optical elements).

3 Two-slit diffraction

We consider the double slit setup sketched in Fig. 1, in which an electromagnetic plane wave propagates in the positive $z$ direction, towards a detection plane, after traversing a double slit screen. The distance between the centers of the two slits, of width $a,$ is $2d.$ The midpoint between the slits is chosen as the origin of the coordinates. We assume quasimonochromatic waves, of frequency $\omega$ and wavenumber $k,$ in the paraxial approximation with a given polarization and a single spatial dimension $x$ on the slit and the detection planes. A simple mathematical description in terms of wave functions, recently used in the investigation of optical superresolution [19, 18], gives the same result as a second quantized treatment in which only the single photon subspace of the full Fock space is taken into account.

The light coming from the i-th slit is described by the wave function $\psi_{i}(x,z=0)=\braket{x}{\psi_{i}(z=0)},$ $i=1,2.$ The paraxial wave equation

[TABLE]

connects the wave function at the slit plane with the corresponding wave function at the detection plane. It is convenient to switch to the dimensionless quantities $\xi={x}/{a},$ and $\zeta={z}/{z_{0}}.$ The paraxial wave equation can be simply written as $\partial_{\xi\xi}\psi+2i\partial_{\zeta}\psi=0,$ where $z_{0}$ was chosen as $z_{0}=ka^{2}=2\pi a^{2}/\lambda.$ Here, $\lambda$ is the wavelength of the light that illuminates the double slit. The general solution of the paraxial wave equation is a superposition of the elementary solutions $e^{i\kappa\xi-i\kappa^{2}(\zeta/2)}$; that is, $\psi(\xi,\zeta)=\int d\kappa\tilde{\psi}(\kappa)e^{i\kappa\xi-i\kappa^{2}(\zeta/2)}$. Noticing that $\tilde{\psi}(\kappa)$ is the Fourier transform of $\psi(\xi^{\prime},0)$, we have

[TABLE]

In the last step, the integral over $\kappa$ was performed.

Under the assumptions made in this section, the most general state of light just outside the slit screen is

[TABLE]

A photon is detected at the (transversal) position $\xi$ on the detection plane with probability (density)

[TABLE]

In the following section, we show that a minimal SIC-POVM tomographic reconstruction of $\rho$, using photon detectors in a single plane, $\zeta=$ constant, is possible.

4 Minimal double slit qubit tomography

The photon detection probability density, at the plane detection $\zeta$ and transversal position $\xi$, can be written in the suggestive form

[TABLE]

where $I(\xi,\zeta)=\sum_{k=1,2}|\psi_{k}(\xi,\zeta)|^{2}$ is the intensity envelope, and $\hat{\Pi}(\xi,\zeta)$, written in the basis $\{\ket{\psi_{1}},\ket{\psi_{2}}\}$, is the projection operator

[TABLE]

Each projection operator projects over a pure state, $\varrho(\xi,\zeta)=\ket{\psi(\xi,\zeta)}\bra{\psi(\xi,\zeta)}$ $=\frac{1}{2}(\mathbb{I}+\boldsymbol{s}(\xi,\zeta)\cdot\boldsymbol{\sigma}),$ where $\mathbb{I}$ is the $2\times 2$ identity matrix, $\boldsymbol{\sigma}$ is a vector whose components are the three Pauli matrices, and $\boldsymbol{s}(\xi,\zeta)$ is the unit-length Bloch vector associated with the pure state $\varrho(\xi,\zeta).$ In a single plane, $\zeta=\zeta_{0}$ ($\zeta_{0}$ a constant), the family of Bloch vectors $\boldsymbol{s}(\xi,\zeta_{0})$ describes a curve over the Bloch sphere, as shown in Fig. 2, as $\xi$ takes values on the real axis.

The photon detection probability density $p(\xi,\zeta)$, equal to the product of the intensity envelope $I(\xi,\zeta)$ and the population $0\leq\braket{\psi(\xi,\zeta)}{\rho(0)}{\psi(\xi,\zeta)}\leq 1,$ is bounded by the intensity envelope. Therefore, $I(\xi,\zeta)$ is the maximum photon detection probability at a given point $\xi$ on the detection plane, as we measure over every possible state $\rho.$

The photon detection probability density $p(\xi,\zeta)$ depends on i) the four elements of the state $\rho(0)$ in the basis $\{\ket{\psi_{1}},\ket{\psi_{2}}\},$ ii) the two wave functions $\psi_{1}(\xi)$ and $\psi_{2}(\xi),$ and iii) the evolution operator from the slit plane to the detection plane. Besides the information on evolution, assumed by the usual state estimation schemes, we additionally assume complete knowledge of the wave functions $\psi_{1}(\xi)$ and $\psi_{2}(\xi).$ Photon detection on a fixed detector plane $\zeta_{0}$ is described by the continuous POVM $\mathbb{I}=\int d\xi\,I(\xi,\zeta_{0})\ket{\xi,\zeta_{0}}\bra{\xi,\zeta_{0}}.$ We want to use a minimal tomography, $\mathbb{M}=\sum_{i=1}^{4}I(\xi_{i},\zeta_{0})\Delta\xi\,\ket{\xi_{i},\zeta_{0}}\bra{\xi_{i},\zeta_{0}},$ which corresponds to four preselected positions of point detectors. The spatial width of the photon detectors, $\Delta\xi,$ is chosen to be small enough to satisfy $\ket{\xi\pm\Delta\xi,\zeta_{0}}\approx\ket{\xi,\zeta_{0}}.$ In order to go from the continuous measurement to the minimal tomography, most of the information contained in the former measurement must be discarded. Hence, the continuous POVM will generally give a better state estimation than the minimal tomography $\mathbb{M}.$ The minimal tomography, taking into account only the four point detectors, is modeled by $\mathbb{M}=\sum_{i=1}^{4}P(\xi_{i},\zeta_{0})\ket{\xi_{i},\zeta_{0}}\bra{\xi_{i},\zeta_{0}},$ where the weigths $P(\xi_{i},\zeta_{0})$ are given by $I(\xi_{i},\zeta_{0})\Delta\xi/\Sigma(\zeta_{0}),$ and $\Sigma(\zeta_{0})=\sum_{i=1}^{4}I(\xi_{i},\zeta_{0})\Delta\xi/2.$

It is well-known that the most accurate minimal POVM, $\sum_{i=1}^{4}(1/2)\left(\mathbb{I}+\boldsymbol{s}_{i}\cdot\boldsymbol{\sigma}\right)/2=\mathbb{I},$ projects over four pure states equally spaced on the Bloch sphere [36, 45]; the tips of the corresponding Bloch vectors, $\boldsymbol{s}_{i},\,i=1,2,3,4$, are the four vertices of a regular tetrahedron. These Bloch vectors satisfy the inner product conditions, $\boldsymbol{s}_{i}\cdot\boldsymbol{s}_{j}=\frac{4}{3}\delta_{i,j}-\frac{1}{3},$ $i,j=1,2,3,4.$ If the family of Bloch vectors $\boldsymbol{s}(\xi,\zeta_{0})$ contains four elements $\boldsymbol{s}(\xi_{i},\zeta_{0}),i=1,2,3,4,$ which satisfy these conditions, a minimal tomographic reconstruction is possible if photon detectors are placed on the plane $\zeta=\zeta_{0}$, with transversal positions $\xi_{i},i=1,2,3,4.$ Notice that the coefficients of the projectors are all equal. Therefore, the intensity envelope, $I(\xi,\zeta_{0}),$ must be exactly the same for the four values of $\xi_{i}.$

In this work, we focus on Gaussian slits, which can be thought of as an approximation to the rectangular slits used in actual experiments. Gaussian slits also describe the incidence of Gaussian beams on a biprism. Mathematically, Gaussian slits are modeled by the normalized wave functions

[TABLE]

where $2\delta$ is the distance between the centers of the two slits, and the standard deviation of the slits is 1 ($a$ in the original variables). For $\zeta\geq 0$, we do have (for k=1,2)

[TABLE]

Therefore, the intensity envelope is

[TABLE]

On the other hand, by setting $w=2\xi\delta/(1+\zeta^{2})$, we can write the projectors as

[TABLE]

Taking into account that $\varrho(\xi,\zeta)=\frac{1}{2}(I+\boldsymbol{s}(\xi,\zeta)\cdot\boldsymbol{\sigma}),$ we obtain

[TABLE]

To search four Bloch vectors $\boldsymbol{s}_{i}=\boldsymbol{s}(\xi_{i},\zeta_{0}),$ $i=1,2,3,4,$ whose tips are the vertices of a tetrahedron, we define the positive function $f(\boldsymbol{s}_{1},\boldsymbol{s}_{2},\boldsymbol{s}_{3},\boldsymbol{s}_{4})=\sum_{i>j}(\boldsymbol{s}_{i}\cdot\boldsymbol{s}_{j}+1/3)^{2},$ which vanishes only when the tetrahedron condition is met. Rough local minima of $f$, obtained with an algorithm that uses random values of $w$ and $\zeta$, are fed to a gradient algorithm, which provides a more accurate approximation. Some solutions, corresponding to the smallest values of $\zeta$ found, are listed in Table 1. Notice the existence of solutions which are not symmetric under reflection, $\xi\to-\xi.$

The four Bloch vectors, for the smallest $\zeta$ solution, are

[TABLE]

They are drawn in Fig. 2, along with the set of projectors obtained by varying the transversal position of the detector. It is easy to see that the $s_{x}=\pm 1/\sqrt{3}$ for symmetric solutions. Indeed, since $\boldsymbol{s}(\xi,\zeta)$ is a unit-length vector and $\boldsymbol{s}(-\xi,\zeta)=(s_{x}(\xi,\zeta),-s_{y}(\xi,\zeta),-s_{z}(\xi,\zeta)),$ we have $\boldsymbol{s}(-\xi,\zeta)\cdot\boldsymbol{s}(\xi,\zeta)=2s_{x}^{2}(\xi,\zeta)-1.$ Finally, when $\boldsymbol{s}(\xi,\zeta)$ and $\boldsymbol{s}(-\xi,\zeta)$ belong to a tetrahedron, its inner product is $-1/3=2s_{x}^{2}(\xi,\zeta)-1$. The desired result follows from this equality.

The numerical solutions that we have found, described by $\mathbb{M}=\sum_{i=1}^{4}\mathbb{P}_{i}=\sum_{i=1}^{4}P(\xi_{i},\zeta_{0})\ket{\xi_{i},\zeta_{0}}\bra{\xi_{i},\zeta_{0}},$ are not true solutions, in the sense that the subnormalized projectors $\mathbb{P}_{i}$ do not add to the identity matrix. This problem originates on the fact that the weights $P(\xi_{i},\zeta_{0}),$ which are proportional to the intensity envelope $I(\xi_{i},\zeta_{0}),$ are not equal. To deal with this problem, one can artificially “balance” the actual counts, multiplying them by an appropriate constant, as have been done in some experimental realizations of minimal tomography. A better alternative would be to examine the behavior of the intensity envelope as a function of the transverse coordinate and the slit separation (see Fig. 3), with the aim to find solutions intrinsically balanced. When we write the intensity envelope $I(\xi,\zeta;\delta)$ as $2\exp({-\frac{(1+\zeta^{2})w^{2}}{4\delta^{2}}})\exp(-\frac{\delta^{2}}{1+\zeta^{2}})\cosh(w)/\sqrt{\pi(1+\zeta^{2})}$, we see that the symmetric solutions are balanced when the distance between the slits is chosen as

[TABLE]

where $w_{1}<w_{2}<0.$ Indeed, when $\delta$ is chosen according to (6), the value of the intensity envelope is the same at the four photon detector positions. In the case of the solution at $\zeta\approx 3.4678$, this tomography is a POVM for $\delta\approx 2.76444.$ In Fig. 3, we have plotted the intensity envelope as a function of the dimensionless transverse coordinate $\xi,$ for several values of the dimensionless distance between the centers of the slits, $\delta.$ The curves have been rescaled in such a way that, for $\zeta=0,$ its maximum value would be unity. Moreover, as a final check, the two slits can be seen to be non-overlapping, by plotting $\psi_{1}(\xi,\zeta=0,\delta)$ and $\psi_{2}(\xi,\zeta=0,\delta)$ (Fig. 4). In fact, $\braket{\psi_{1}(\xi,\zeta;\delta)}{\psi_{2}(\xi,\zeta;\delta)}=e^{-\delta^{2}}$ is of the order of $10^{-4}$ for $\delta\approx 2.76.$

5 Final remarks

In this paper, we consider the implementation of the SIC tetrahedron tomography for double slit qubits of light using measurements performed on a single plane. We have found that $\lambda$, the wavelength of the photons used in the experiment, and $a$, the width of the slits, can be used to determine the geometry of the experimental setup which implements the optimal state estimation, and does not require to artificially balance the measured probabilities. The particular details of the solution depend on the states which are propagated from the double slit screen. Assuming plane waves and Gaussian slits, the optimal geometry, for detectors placed nearest to the double slit plane, features a distance between the centers of the slits equal to 5.53 $a$. The distance between the slit-plane and the detection-plane is $3.47\times 2\pi a^{2}/\lambda$. The detectors must be placed at the transversal positions $x_{1}=-x_{4},$ $x_{2}=-x_{3},$ $x_{3}$ and $x_{4}$, measured from the midpoint between the slits (but on the detection plane). In Table 2, we give some typical values of the geometry, by using wavelengths and slit widths reported in the literature.

The results of this work can be applied not only to the classical version of this setup, but also to spatial qubits of matter. Moreover, the tomographic reconstruction of two double slit qubits can be carried out by joint measurements of two single qubit SIC tetrahedron projectors. An approach similar to the one used in this work makes it possible to investigate the implementation of the SIC tomography for d-slit qudits of light, using measurements performed on a minimum number of planes.

Acknowledgements

We thank Prof. M. C. Nemes for suggesting this problem, Dr. Breno Marques for stimulating discussions, and an anonymous referee for his/her constructive input. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.