Quantum Tomography of the Photon-Plasmon Conversion Process in a Metal Hole Array
Lei Tang, Kaimin Zheng, Jiale Guo, Yi Ouyang, Yang Wu, Chuanqing Xia, Long Li, Fang Liu, Yong Zhang, Lijian Zhang, Min Xiao
QQuantum Tomography of the Photon -Plasmon Conversion Process in a Metal Hole Array L EI T ANG , K AIMIN Z HENG , J IALE G UO , Y I O UYANG , Y ANG W U , C HUANQING X IA , L ONG L I , F ANG L IU , Y ONG Z HANG , L IJIAN Z HANG , AND M IN X IAO National Laboratory of Solid State Microstructures, College of Engineering and Applied Sciences, and School of Physics, Nanjing University, Nanjing 210093, China Department of Physics, University of Arkansas, Fayetteville, Arkansas 72701, USA β These authors contribute equally to this work *[email protected], *[email protected], *[email protected]
Abstract:
In the past decades, quantum plasmonics has become an active area due to its potential applications in on-chip plasmonic devices for quantum information processing. However, the fundamental physical process, i.e., how a quantum state of light evolves in the photon-plasmon conversion process, has not been clearly understood. Here, we report a complete characterization of the plasmon-assisted extraordinary optical transmission process through quantum process tomography. By inputting various coherent states to interact with the plasmonic structure and detecting the output states with a homodyne detector, we reconstruct the process tensor of the photon-plasmon conversion process. Both the amplitude and phase information of the process are extracted, which explains the evolution of the quantum-optical state after the coupling with plasmons. Our experimental demonstration constitutes a fundamental block for future on-chip applications of quantum plasmonic circuits.
1. Introduction
Plasmonics provides the capabilities to localize and manipulate electromagnetic excitations within sub-wavelength scales, and therefore has great potentials in the miniaturization and scalability of photonic devices[1-5]. In recent years, the newly emerging field of quantum plasmonics[6,7] has become an active research area, which is mostly motivated by its potential applications in integrated plasmonic circuits for quantum information processing. Quantum plasmons excited by various quantum sources including single photon[8-10], two-photon entanglement[11,12] and squeezing[13] have been experimentally demonstrated. The fundamental physical characteristics of quantum plasmons have been investigated in intense light-matter interactions[14-16], quantum confinement effects[17], survival of entanglement and squeezing[11,18-20], decoherence and loss[8,11,18], perfect absorption of entangled photons[21] etc. Quantum plasmonic devices such as detectors[22,23], interferometers[9,24,25], and controlled-NOT gate[26] have been developed. The applications have been further extended to quantum plasmonic sensing and quantum plasmonic networks beyond the classical limit[27,28]. However, large-scale usage of quantum plasmons in quantum information processing still faces many technical challenges to be solved, such as reducing loss and enhancing efficiency of plasmon-based quantum optical logic gates, controlling the phase of quantum interference devices and understanding at what scale the quantum theory has to be considered in a plasmonic device. More importantly, the fundamental physical problem, how an arbitrary quantum state of light evolves through the coupling with plasmons, has not been investigated systematically. Physical understanding of quantum plasmons can be acquired in two steps through the evolution of various degrees of freedom of the optical field during the photon-plasmon conversion process. The first one is about the degree of freedom of single photons, e.g. polarization, spatial, temporal and spectral, which have been studied extensively[11,18,29]. However, the quantum fluctuations of photon-number degree of reedom and, moreover, the quasi-probability distribution in the phase space mark the difference between classical and non-classical optical fields, and are crucial for various quantum information applications[27,30,31]. Therefore, a complete understanding of the evolution of quantum fluctuations and coherence inside the plasmonic structure is indispensable. In this work we address this problem by applying the coherent-state quantum process tomography[32] (CSQPT) to a typical plasmon-assisted process, i.e. extraordinary optical transmission[33,34](EOT), in a metal-hole array. This transmission is largely attributed to the excitation of surface plasmon polaritons (SPPs) in the literature. Although subsequent theoretical and experimental work has claimed the presence of other surface waves[35] in this process, SPPs still play an important role in EOT process at visible and near-infrared frequencies. And it has been shown that such process can maintain certain quantum properties, for example, two-photons entanglement[11] and squeezing[13]. Here, we focus the character of far-fields in EOT process, by inputting a set of coherent states through a metal hole array and performing a tomographic reconstruction of the output states, we provide a complete characterization of the plasmon-assisted EOT process that spans the effect of the process in the photon-number Hilbert space. In particular, the reconstruction allows to predict the output state of any input states including non-classical ones.
2. Coherent state quantum process tomography
A quantum process π can be described by a positive, trace-preserving linear map that transfers the input states to the output states over Hilbert space Ξ . Complete characterization of a quantum process therefore means to know the effect of the process on arbitrary quantum states. Assuming {πΜ π } is the complete set spanning the single-mode Hilbert space, any input state can be decomposed as πΜ ππ = β π ππ πΜ π . (1) Therefore complete characterization of a quantum process is equivalent to determine the output state π(πΜ π ) for each πΜ π . Once this information is acquired, one can predict the output state πΜ ππ’π‘ = π(πΜ ππ ) through πΜ ππ’π‘ = β π π π(πΜ π ) π . (2) The challenge associated with this approach is the construction of the appropriate complete set. For the optical field a natural candidate is the photon-number, i . e . Fock basis {|πβ©β¨π|} . Under this basis, the process can be expressed as a rank-4 tensor π ππππ that relates the density matrix of the output state πΜ ππ’π‘ and that of the input state πΜ ππ : π ππππ’π‘ = β π ππππ π ππππππ , (3) where π ππππ = β¨π|πΜ ππ |πβ©, π ππππ’π‘ = β©π|πΜ ππ’π‘ |πβͺ and π ππππ = β©π|π(|πβͺβ©π|)|πβͺ . Characterization of a quantum process can be achieved through the reconstruction of the tensor, which is known as the quantum process tomography (QPT). However, direct reconstruction of π ππππ requires superpositions of different Fock states as the probe states, which may be infeasible with current techniques. Luckily, there is another set of complete basis, the coherent states {|πΌβͺ}, which is readily generated with the output of a laser. Any quantum state of light πΜ can be written as a linear combination of density matrices of coherent states |πΌβͺ πΜ ππ = 2 β« π ππ (πΌ)|πΌβͺβ©πΌ| π πΌ , (4) where π ππ (πΌ) is the stateβs Glauber-Sudarshan P function and the integration is performed over the entire complex plane. The output state can be expressed as: πΜ ππ’π‘ = 2 β« π ππ (πΌ) π(|πΌβͺβ©πΌ|) π πΌ . (5) Thus, it is sufficient to know the output states of every coherent state |πΌβͺ . In particular the process tensor in Fock basis can be reconstructed using[32] π ππππ = 2 β« π ππ (πΌ) β©π| π(|πΌβͺβ©πΌ|)|πβͺ π πΌ , (6) where π ππ (πΌ) is the P function of |πβ©β¨π| . In practical setups, experimental imperfections and statistical fluctuations are unavoidable. Therefore, instead of directly apply Eq. (6), maximum likelihood estimation (MLE)[36-38] can be used to mitigate of effect of noise. . Sample and experimental setup In the experiment, we utilize a typical plasmon-assisted EOT system to investigate the fundamental mechanism inside the photon-plasmon-photon conversion process. The sample is a hexagonal metal-hole array in a 100-nm-thick Au film sputtered on a glass substrate, which is fabricated by using a focused ion beam system (FEI Helios 600i). Each circular hole has a diameter of 460 nm and the array period is 759 nm. The shape of the transmission spectrum of such hexagonal hole array does not depend on the input polarization [39]. Because of the plasmon resonance, the transmissivity of our sample is greatly enhanced to be 62.0% at 1080 nm as shown in Fig. 1(a). The metal holes occupy ~33.4% of the area in our sample. The transmission efficiency π π΅ normalized to the aperture area is 1.86, which indicates that the EOT effect happens [40]. The discrepancy between the experimental data and theoretical simulations is due to the imperfections in the fabrication. We prepare a set of probe states and transmit them through the sample. The quadrature distribution of the output states are measured with a homodyne detector, which allows to reconstruct the Wigner function and density matrix of the output state[41,42] as well as the process tensor. Fig. 1. Transmission spectrum of gold plasmon sample and Experimental setup. (a) The blue slope is EOT transmissivity slope and the red dot is theory transmissivity slope (FDTD simulation). Their characteristic peaks are both at 1080nm. Inset: The electron microscope photo of our metal-hole arrays fragment. Its full size is 65ΞΌmΓ65ΞΌm. Period is 759nm. Hole diameter is 460nm. (b) A coherent infrared light at 1080nm is amplitude modulated with an electrooptical modulator and passes through a calibrated neutral density filter to prepare the probe states that incident on the plasmonic sample. The output state is measured with a homodyne detector. The relative phases between the probe states and local oscillator of the homodyne detector are set with a piezoelectric transducer. EOM is amplitude electrooptical modulator. ND is neutral attenuation piece. BS is beam splitter. PZT is piezoelectric transducer. OSC is high frequency signal generator. LPF is low frequency filter.
The schematic experimental setup is shown in Fig. 1(b). The signal field is a coherent infrared light at 1080nm, which is amplitude-modulated by an electro-optical modulator (EOM). The modulation frequency is set at 2 MHz and the applied voltage on the EOM is tunable between 0V and 10V. The modulated signal field is attenuated by a neutral density filter to prepare the probe states |πΌβ© with average photon number |πΌ| less than 10 and then is focused on the metal hole array sample. A two-lens system is used to reshape the incident beam to match the sample area. For the experiment we use 9 different probe states with the modulation voltage f the EOM from 0V to 8V with 1V increment, which give πΌ =
0, 0.1375, 0.2750, 0.4125, 0.5500, 0.6875, 0.8250, 0.9625, and 1.100, respectively. The transmitted light after the photon-plasmon-photon conversion process is measured with a homodyne detection system. The signal field passes through the metal hole sample and then interferes with a strong local field. The relative phase between the two fields is scanned by a piezoelectric transducer (PZT) from 0 to 2Ο and tracked with an ancillary beam. A pair of photodiodes are used to detect the interference intensity and a subtraction is performed on the AC signals. The subtracted AC signal is sent to a lock-in system. After frequency mixing and low-frequency filtering, the signal is collected with an oscilloscope. Such homodyne detection method can efficiently increase the signal-to-noise ratio in our measurement. We can obtain the quadrature information at different angles by slowly scanning the PZT with a frequency of 2Hz.
4. Quantum state tomography
First, we remove the sample and perform the quantum state tomography of the input coherent states at various modulation voltages to calibrate the input states and test the system. By scanning the phase of the local field between 0 and 2Ο, we collect 5Γ10 data points, which are equally divided into 20 phase sections. The collected data are used to reconstruct the density matrices πΜ ππ in Fock basis of the input states with MLE[43]. The phase-space quasiprobablity distribution, known as the Wigner function can be calculated from πΜ ππ π(π, π) = β« β¨π + π β² |πΜ ππ | ββ π β π β² β©exp (βππβ²π)ππβ² , (7) where π and π are two quadratures. The Wigner functions of the input states at 2V, 6V, and 8V modulation voltages are shown in Fig. 2(a), (c) and (e), respectively. Clearly, the center of the Wigner function moves away from the origin of the coordinates as increasing the modulation voltage since the amplitude of the input state increases with the modulation voltage. The Wigner functions possess Gaussian shapes and the distances between the centers and the origin are 0.2750, 0.8250, and 1.100, which corresponds to coherent stats |πΌβ© with πΌ = Fig. 2. Wigner function of input and output light. (a/c/e) The Wigner functions of input fields W(X,Y) input with 2V, 6V and 8V modulation voltage of the EOM, respectively. (b/d/f) The Wigner functions of output fields W(X,Y) output with 2V, 6V and 8V modulation voltage of the EOM, respectively. πΌ ππ corresponds to the input state and πΌ ππ’π‘ corresponds to the output state.
5. Quantum process tomography of metal-hole arrays
The changes in the Wigner functions of specific input states do not reveal all the information about the conversion process. We further apply CSQPT on the EOT process with 9 coherent states. We collect the outcomes of the homodyne detectors for each output state, and apply MLE algorithm[44] to reconstruct process tensor π ππππ . Here, π ππππ is a four-dimensional tensor. Since the EOT process has influences on both the amplitude and phase characteristics of the input states, we investigate them separately. We extract the diagonal elements of the tensor π ππππ to investigate the amplitude information (Fig. 3a). We choose m = n in the input state and k = l in the output state, which corresponds to he diagonal elements of the density matrices and describes the effect of the process on the photon-number distribution of the state. By analyzing the diagonal elements π ππππ of the process tensor, one can see the evolution of photon numbers from input field ( m ) to output field ( k ). For a given input photon number m , the output photon number k has a Binomial-like distribution, which hints the Bernoulli transformation of the process. By assuming a linear loss process with transmissivity T = 62.0%, we perform numerical simulations as shown in Fig. 3(b). By comparing the experimental results in Fig. 3(a) and the theoretical simulations in Fig. 3(b), we can get the fidelity F = 99.68% between the two matrices. Clearly, the EOT effect on the amplitude of the input state can be understood as a linear loss. Fig. 3. The results of CSQPT. (a)The diagonal elements of the process tensor π ππππ with input field index m and output field index k for an EOT process. (b)Numerical results of a linear loss process with transmissivity of 62.0%. (c/d/e) The off-diagonal elements of the process tensor Im{In[π ]} (c), Im{In[π ]} (d), Im{In[π ]} (e), where m and n represent the input field index in the Fock basis. In addition to changes in the photon-number distributions, the reconstructed process tensor also reveals the effect on the off-diagonal elements of the density matrix, i.e. the coherence between different photon-number components. In particular, we can acquire the effect of the process on the phases of the input state. Such information can be extracted from the tensor elements mapping the input state density matrix to certain off-diagonal element of the output state. For example, the phase value the π element of the output state is determined by the phases Im{ln[π ]} of π process tensor elements[45]. To elaborate this relation, we decompose the process π into a phase shift superoperator and a phase-symmetric process π β² πΜ ππ’π‘ = π(πΜ ππ ) = πΜ(π)π β² (πΜ ππ )πΜ β (π) = π πππΜ β πΜ π β² (πΜ ππ )π βπππΜ β πΜ , (8) where π is a constant phase shift. From Eq. (3), we can get the process tensor π ππππ =|πβ² ππππ |π ππ ππππ β π π(πβπ)π and Im{ln[π ππππ ]} = π ππππ + (π β π)π , (9) where π ππππ is the phase of πβ² ππππ , which accounts for photon-number-dependent, i.e. nonlinear, phase shifts. Fig. 3(c), (d) and (e) show several phase values Im{ln[π ππππ ]} of the EOT process, which are extracted from π , π , and π , respectively. Note for a phase-symmetric process, only the tensor elements π ππππ with π β π = π β π are non-zero[32]. The figures show that for given k and l , the phase changes are the same for various m and n , which well matches the prediction rom the Wigner functions (Fig. 2). These results indicate that there is only a constant linear phase shift involved in the process. In addition, the average values of the phases of π , π , and π are 0.9200, 1.8182, and 2.7849, respectively, which are linearly dependent on π β π . In short, the metal-hole arrays apply a constant phase shift with a value about 0.92 radian to the input state.
6. Discussion and conclusion
Although to fully understand the photon-plasmon-photon conversion process requires a detailed microscopic model, the reconstructed process tensor suggests that the conversion process can be effectively described with a beam-splitter model plus a constant phase. With this model, the input field operator πΜ ππ is transferred to the output field operator πΜ ππ’π‘ through the equation πΜ ππ’π‘ = β1 β ππΜ πππ‘β + π ππ βππΜ ππ , (10) where πΜ πππ‘β is the field operator of the lumped Markovian environment in vacuum, π is the overall transmissivity and π is the constant phase added by the process. The linear-loss model has been proposed in Refs. 8, 19 and tested with squeezed light[19] and heralded single photons[46]. Our tomography results go beyond the previous works by confirming validity of the model for arbitrary input states. Moreover, the new results also highlight the effect of the coupling process on the phase of the input state, which is closely related to quantum coherence. The model given in Eq. (10) can be applied to both classical and quantum optical field, and explains the extraordinary high transmissivity of light and entanglement survival[11,18-20]. To summarize, we have performed complete quantum tomography to reconstruct a photon-plasmon conversion process, i.e., EOT in a metal-hole array. Such reconstruction procedure discovers the fundamental characteristics including both the amplitude and phase information of a typical plasmonic process, which allows us to precisely estimate the interaction of arbitrary classical or quantum optical fields with such plasmonic structure, as well as to develop the microscopic model to interpret the process. Our experimental demonstration provides a fundamental understanding of a plasmon-assisted EOT process, which paves a way for the proper design of a quantum plasmonic components for future applications in on-chip quantum information processing. It should be noted that our experimental observations based on single input mode can be readily extended to the case with multiple optical modes. Since the applications of quantum plasmonics have been further extended to quantum plasmonic sensing and quantum plasmonic networks beyond the classical limit, the complete quantum tomography of various quantum plasmonic process will benefit the optimization of on-chip plasmonic devices. Funding
This work was supported by National Key R&D Program of China (2017YFA0303700, 2016YFA0302500), National Natural Science Foundation of China (NSFC) (91636106, 11874213, 11621091, 11474159, 61490711, 11874213).
References
1. J. A. Schuller, E. S. Barnard, W. Cai, Y. C. Jun, J. S. White, and M. L. Brongersma, "Plasmonics for extreme light concentration and manipulation," Nat. Mater. , 193β204 (2010). 2. E. Ozbay, "Plasmonics: Merging Photonics and Electronics at Nanoscale Dimensions," Science , 189β193 (2006). 3. Y. Zhu, D. Wei, Z. Kuang, Q. Wang, Y. Wang, X. Huang, Y. Zhang, and M. Xiao, "Broadband Variable Meta-Axicons Based on Nano-Aperture Arrays in a Metallic Film," Sci. Rep. , 11591 (2018). 4. Y. Wang, X. Fang, Z. Kuang, H. Wang, D. Wei, Y. Liang, Q. Wang, T. Xu, Y. Zhang, and M. Xiao, "On-chip generation of broadband high-order Laguerre-Gaussian modes in a metasurface," Opt. Lett. , 2463β2466 (2017). . D. Wei, Y. Wang, D. Liu, Y. Zhu, W. Zhong, X. Fang, Y. Zhang, and M. Xiao, "Simple and Nondestructive On-Chip Detection of Optical Orbital Angular Momentum through a Single Plasmonic Nanohole," ACS Photonics , 996β1002 (2017). 6. M. S. Tame, K. R. McEnery, Ε. K. Γzdemir, J. Lee, S. A. Maier, and M. S. Kim, "Quantum plasmonics," Nat. Phys. , 329β340 (2013). 7. F. Marquier, C. Sauvan, and J.-J. Greffet, "Revisiting Quantum Optics with Surface Plasmons and Plasmonic Resonators," ACS Photonics , 2091β2101 (2017). 8. M. S. Tame, C. Lee, J. Lee, D. Ballester, M. Paternostro, A. V. Zayats, and M. S. Kim, "Single-Photon Excitation of Surface Plasmon Polaritons," Phys. Rev. Lett. , 190504 (2008). 9. Y.-J. Cai, M. Li, X.-F. Ren, C.-L. Zou, X. Xiong, H.-L. Lei, B.-H. Liu, G.-P. Guo, and G.-C. Guo, "High-Visibility On-Chip Quantum Interference of Single Surface Plasmons," Phys. Rev. Applied , 014004 (2014). 10. M.-C. Dheur, E. Devaux, T. W. Ebbesen, A. Baron, J.-C. Rodier, J.-P. Hugonin, P. Lalanne, J.-J. Greffet, G. Messin, and F. Marquier, "Single-plasmon interferences," Sci. Adv. , e1501574 (2016). 11. E. Altewischer, M. P. van Exter, and J. P. Woerdman, "Plasmon-assisted transmission of entangled photons," Nature , 304β306 (2002). 12. X. F. Ren, G. P. Guo, Y. F. Huang, C. F. Li, and G. C. Guo, "Plasmon-assisted transmission of high-dimensional orbital angular-momentum entangled state," EPL , 753 (2006). 13. D. Wang, C. Xia, Q. Wang, Y. Wu, F. Liu, Y. Zhang, and M. Xiao, "Feedback-optimized extraordinary optical transmission of continuous-variable entangled states," Phys. Rev. B , 121406 (2015). 14. A. V. Akimov, A. Mukherjee, C. L. Yu, D. E. Chang, A. S. Zibrov, P. R. Hemmer, H. Park, and M. D. Lukin, "Generation of single optical plasmons in metallic nanowires coupled to quantum dots," Nature , 402β406 (2007). 15. A. Huck, S. Kumar, A. Shakoor, and U. L. Andersen, "Controlled Coupling of a Single Nitrogen-Vacancy Center to a Silver Nanowire," Phys. Rev. Lett. , 096801 (2011). 16. A. I. FernΓ‘ndez-DomΓnguez, S. I. Bozhevolnyi, and N. A. Mortensen, "Plasmon-Enhanced Generation of Nonclassical Light," ACS Photonics , 3447β3451 (2018). 17. P. Tassin, T. Koschny, M. Kafesaki, and C. M. Soukoulis, "A comparison of graphene, superconductors and metals as conductors for metamaterials and plasmonics," Nat. Photonics , 259β264 (2012). 18. S. Fasel, F. Robin, E. Moreno, D. Erni, N. Gisin, and H. Zbinden, "Energy-Time Entanglement Preservation in Plasmon-Assisted Light Transmission," Phys. Rev. Lett. , 110501 (2005). 19. A. Huck, S. Smolka, P. Lodahl, A. S. SΓΈrensen, A. Boltasseva, J. Janousek, and U. L. Andersen, "Demonstration of Quadrature-Squeezed Surface Plasmons in a Gold Waveguide," Phys. Rev. Lett. , 246802 (2009). 20. B. J. Lawrie, P. G. Evans, and R. C. Pooser, "Extraordinary Optical Transmission of Multimode Quantum Correlations via Localized Surface Plasmons," Phys. Rev. Lett. , 156802 (2013). 21. C. Altuzarra, S. Vezzoli, J. Valente, W. Gao, C. Soci, D. Faccio, and C. Couteau, "Coherent Perfect Absorption in Metamaterials with Entangled Photons," ACS Photonics , 2124β2128 (2017). 22. R. W. Heeres, S. N. Dorenbos, B. Koene, G. S. Solomon, L. P. Kouwenhoven, and V. Zwiller, "On-Chip Single Plasmon Detection," Nano Lett. , 661β664 (2010). 23. A. L. Falk, F. H. L. Koppens, C. L. Yu, K. Kang, N. de Leon Snapp, A. V. Akimov, M.-H. Jo, M. D. Lukin, and H. Park, "Near-field electrical detection of optical plasmons and single-plasmon sources," Nat. Phys. , 475β479 (2009). 24. R. W. Heeres, L. P. Kouwenhoven, and V. Zwiller, "Quantum interference in plasmonic circuits," Nat. Nanotechnol. , 719β722 (2013). 25. J. S. Fakonas, H. Lee, Y. A. Kelaita, and H. A. Atwater, "Two-plasmon quantum interference," Nat. Photonics , 317β320 (2014). 26. S. M. Wang, Q. Q. Cheng, Y. X. Gong, P. Xu, C. Sun, L. Li, T. Li, and S. N. Zhu, "A 14 Γ 14 ΞΌm2 footprint polarization-encoded quantum controlled-NOT gate based on hybrid waveguide," Nat. Commun. , 11490 (2016). 27. C. Lee, F. Dieleman, J. Lee, C. Rockstuhl, S. A. Maier, and M. Tame, "Quantum Plasmonic Sensing: Beyond the Shot-Noise and Diffraction Limit," ACS Photonics , 992β999 (2016). 28. M. W. Holtfrerich, M. Dowran, R. Davidson, B. J. Lawrie, R. C. Pooser, and A. M. Marino, "Toward quantum plasmonic networks," Optica, , 985β988 (2016). 29. S. Fasel, M. Halder, N. Gisin, and H. Zbinden, "Quantum superposition and entanglement of mesoscopic plasmons," New J. Phys. , 13β13 (2006). 30. W. Fan, B. J. Lawrie, and R. C. Pooser, "Quantum plasmonic sensing," Phys. Rev. A , 053812 (2015). 31. M. Dowran, A. Kumar, B. J. Lawrie, R. C. Pooser, and A. M. Marino, "Quantum-enhanced plasmonic sensing," Optica, , 628β633 (2018). 32. M. Lobino, D. Korystov, C. Kupchak, E. Figueroa, B. C. Sanders, and A. I. Lvovsky, "Complete Characterization of Quantum-Optical Processes," Science , 563β566 (2008). 33. T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and P. A. Wolff, "Extraordinary optical transmission through sub-wavelength hole arrays," Nature , 667β669 (1998). 34. F. J. GarcΓa de Abajo, "Colloquium: Light scattering by particle and hole arrays," Rev. Mod. Phys. , 1267β1290 (2007). 35. F. van Beijnum, C. RΓ©tif, C. B. Smiet, H. Liu, P. Lalanne, and M. P. van Exter, "Quasi-cylindrical wave contribution in experiments on extraordinary optical transmission," Nature , 411β414 (2012). 6. A. I. Lvovsky, "Iterative maximum-likelihood reconstruction in quantum homodyne tomography," J. Opt. B: Quantum Semiclass. Opt. , S556βS559 (2004). 37. J. FiurΓ‘Ε‘ek and Z. Hradil, "Maximum-likelihood estimation of quantum processes," Phys. Rev. A , 020101 (2001). 38. M. JeΕΎek, J. FiurΓ‘Ε‘ek, and Z. Hradil, "Quantum inference of states and processes," Phys. Rev. A , 012305 (2003). 39. T. Thio, H. F. Ghaemi, H. J. Lezec, P. A. Wolff, and T. W. Ebbesen, "Surface-plasmon-enhanced transmission through hole arrays in Cr films," J. Opt. Soc. Am. B, JOSAB , 1743β1748 (1999). 40. C. Genet and T. W. Ebbesen, "Light in tiny holes," Nature , 39β46 (2007). 41. G. M. DβAriano, C. Macchiavello, and M. G. A. Paris, "Detection of the density matrix through optical homodyne tomography without filtered back projection," Phys. Rev. A , 4298β4302 (1994). 42. U. Leonhardt, M. Munroe, T. Kiss, T. Richter, and M. G. Raymer, "Sampling of photon statistics and density matrix using homodyne detection," Opt. Commun. , 144β160 (1996). 43. A. I. Lvovsky and M. G. Raymer, "Continuous-variable optical quantum-state tomography," Rev. Mod. Phys. , 299β332 (2009). 44. S. Rahimi-Keshari, A. Scherer, A. Mann, A. T. Rezakhani, A. I. Lvovsky, and B. C. Sanders, "Quantum process tomography with coherent states," New J. Phys. , 013006 (2011). 45. C. Kupchak, S. Rind, B. Jordaan, and E. Figueroa, "Quantum Process Tomography of an Optically-Controlled Kerr Non-linearity," Sci. Rep. , 16581 (2015). 46. G. Di Martino, Y. Sonnefraud, S. KΓ©na-Cohen, M. Tame, Ε. K. Γzdemir, M. S. Kim, and S. A. Maier, "Quantum Statistics of Surface Plasmon Polaritons in Metallic Stripe Waveguides," Nano Lett.12