Literature DB >> 34767450

Programmable and sequential Gaussian gates in a loop-based single-mode photonic quantum processor.

Yutaro Enomoto1, Kazuma Yonezu1, Yosuke Mitsuhashi1, Kan Takase1, Shuntaro Takeda1,2.   

Abstract

A quantum processor to import, process, and export optical quantum states is a common core technology enabling various photonic quantum information processing. However, there has been no photonic processor that is simultaneously universal, scalable, and programmable. Here, we report on an original loop-based single-mode versatile photonic quantum processor that is designed to be universal, scalable, and programmable. Our processor can perform arbitrarily many steps of programmable quantum operations on a given single-mode optical quantum state by time-domain processing in a dynamically controlled loop-based optical circuit. We use this processor to demonstrate programmable single-mode Gaussian gates and multistep squeezing gates. In addition, we prove that the processor can perform universal quantum operations by injecting appropriate ancillary states and also be straightforwardly extended to a multimode processor. These results show that our processor is programmable, scalable, and potentially universal, leading to be suitable for general-purpose applications.

Entities:  

Year:  2021        PMID: 34767450      PMCID: PMC8589304          DOI: 10.1126/sciadv.abj6624

Source DB:  PubMed          Journal:  Sci Adv        ISSN: 2375-2548            Impact factor:   14.136


INTRODUCTION

A quantum processor capable of performing universal quantum operations on a given arbitrary optical quantum state is a common core technology enabling various photonic quantum information processing. Requirements for such a general-purpose quantum processor are (i) the basic functionality to import, process, and export optical quantum states; (ii) ability to perform a universal set of operations; (iii) scalability to efficiently perform arbitrary steps of operations; and (iv) programmability to perform various tasks. Such a processor is at the heart of continuous-variable quantum computing (, ) as well as pre-/postprocessing in quantum communication and sensing protocols (–). It is also an essential building block to realize Hamiltonian simulation for various quantum systems (–) and quantum neural networks for machine learning applications (–). Thus far, most of the elementary building blocks for photonic quantum processors have been demonstrated separately. In general, universal quantum processors are required to perform arbitrary Gaussian operations (such as squeezing and phase shift gates) for linear transformation and at least one non-Gaussian operation (such as a cubic phase gate) for nonlinear transformation. Sequential application of Gaussian and non-Gaussian operations is proven to constitute arbitrary transformations (). Gaussian operations alone are already useful for many quantum communication and sensing protocols (, , ), while non-Gaussian operations are proven to be necessary for quantum computational speedup over classical computers (). On the experimental side, optical Gaussian and non-Gaussian operations are efficiently implementable with “measurement and feedforward” schemes (–). Proof-of-principle demonstrations of such single-step Gaussian gates have been reported already (–). Recent research interests are moving on to integrating these building blocks to realize general-purpose quantum processors. However, to date, no photonic processor has been demonstrated that simultaneously satisfies all the requirements (i) to (iv) for the general-purpose processor. In the context of continuous-variable quantum computing, multistep quantum operations have been demonstrated with the measurement and feedforward schemes (, ), but these early demonstrations relied on spatial-mode encoding, making them difficult to scale up. These schemes were later extended to temporal-mode encoding to achieve high scalability (–), and programmable multistep quantum computations have been recently reported with temporal mode cluster states (, ). However, these implementations were designed specifically for the purpose of one-way quantum computing, which internally prepares input quantum states for computation and returns only classical calculation results instead of output quantum states. Thus, these implementations lack essential ingredients for general-purpose quantum processors such as an input port to take in external states or optical feedforward to finalize quantum operations. Here, we report on an original loop-based single-mode versatile photonic quantum processor that is designed to satisfy all the requirements for the general-purpose processor. Our loop-based processor can store an externally injected single-mode optical state (, ) and perform arbitrarily many steps of quantum operations by repeating a measurement and feedforward process in the loop. This time-domain processing is achieved by developing dynamical and synchronized control of the circuit parameters at a nanosecond time scale, such as transmissivity of beam splitters (BSs), measurement bases, and availability of feedforward. The dynamical controllability also allows us to fully program gate sequences without changing the hardware configurations. In this work, we use this processor to demonstrate programmable single-mode Gaussian gates and multistep squeezing gates. Moreover, we prove that the processor can perform a non-Gaussian gate by injecting appropriate ancillary states. The non-Gaussian gate is not implemented here, but its in-principle applicability and the multistep operationability demonstrated in this work show that our processor can perform universal single-mode quantum operations by sequential Gaussian and non-Gaussian gates. These results together show that our processor is scalable, programmable, and potentially universal, leading to be suitable for general-purpose applications. This work can be straightforwardly extended to a multimode processor () and will accelerate the realization of various quantum optical applications as a universal processing unit.

RESULTS

Working principle of the loop-shaped processor

The basic action of our processor is to let the input wave packet circulate in a loop and repeatedly perform measurement-induced Gaussian/non-Gaussian gates with the help of ancillary states (). Figure 1A shows the conceptual diagram of the loop-shaped optical circuit. The input state to the circuit is in the form of a single wave packet. The wave packet of each ancillary state is sequentially injected into the circuit via an optical switch (Switch-1) that transmits either the input or the ancillary states. These states are then coupled to the loop by a variable BS (VBS), the transmissivity of which is dynamically controlled. In the loop part, displacement and phase shift operations are performed. The amount of displacement and phase shift operation is also dynamically controlled, which is represented by time-dependent parameters of g(t) and θ(t) (Fig. 1A). Downstream of the loop part, there is another optical switch (Switch-2). One output connection of Switch-2 leads to a homodyne detector (HD) with a dynamically tunable measurement basis; is measured, where and are the quadrature operators of the light field, and ϕ is called the homodyne angle. The measurement outcome of the HD is fed forward to the displacement operation inside the loop when necessary. The other connection of Switch-2 is coupled to free space, which serves as the output port of the whole processor.
Fig. 1.

Conceptual schematic of the loop-based circuit.

(A) Physical configuration of the processor. The parameter of each optical and electrical component is dynamically changed as a single wave packet passes through it. PS, phase shift operation; Disp., displacement operation. (B) Equivalent circuit model of the processor. Single time-dependent components are mapped to multiple time-independent components. Temporal modes of the ancillary states are mapped to spatial modes. The amount of displacement operation can actually be dependent on all the past measurement results, although this figure implies that it is determined only by the last measurement result. The measurement basis can also be dependent on the past results.

Conceptual schematic of the loop-based circuit.

(A) Physical configuration of the processor. The parameter of each optical and electrical component is dynamically changed as a single wave packet passes through it. PS, phase shift operation; Disp., displacement operation. (B) Equivalent circuit model of the processor. Single time-dependent components are mapped to multiple time-independent components. Temporal modes of the ancillary states are mapped to spatial modes. The amount of displacement operation can actually be dependent on all the past measurement results, although this figure implies that it is determined only by the last measurement result. The measurement basis can also be dependent on the past results. Thanks to the dynamical feature of each component, the loop-shaped circuit can be virtually transformed into a sequential optical circuit, which is depicted in Fig. 1B. This figure essentially shows that single time-dependent components can be viewed as multiple time-independent components. In the case of the VBS, for example, the single VBS is interpreted as the multiple BSs having independent constant values of reflectivity, which are named BS-0 to BS-(n + 1). As can be seen in Fig. 1B, the processor is capable of applying measurement-induced operations arbitrarily many times to the input state. Moreover, such sequential operations are sufficient for universal single-mode operations. Except for already built-in displacement operations, the processor in Fig. 1A can implement any single-mode Gaussian gate by sequential operations of a phase shift, measurement-based squeezing with an ancillary squeezed state, and another phase shift (). A cubic phase gate, as an example of non-Gaussian gates, can also be performed by using proper ancillae based on a scheme that we propose in this work (refer to Materials and Methods). This scheme is an alternative to the one in () and has the advantage that a cubic phase gate is performed with the single-loop optical path instead of the dual loop. Therefore, the processor can perform universal single-mode operations by sequential Gaussian/non-Gaussian gates () with appropriate ancillary states. In addition to the above property, the operations by the processor are programmable as well since the parameters of the circuit are dynamically controlled from outside.

Single-step Gaussian gate

To demonstrate the programmability of the processor, we perform several different single-step Gaussian gates by the measurement-based scheme in (), which first couples an input state to an ancillary squeezed state, then measures one part of the state, and lastly performs feedforward to the other part. As representatives among all the Gaussian gates, we choose squeezing gates and quadratic phase gates (QPGs). Input-output relations of these gates can be expressed in a matrix form as or withwhere and are the quadrature operators (refer to Materials and Methods for notations and conventions) of the input state, while and are those of the output state. Let us call r and κ as a squeezing parameter and a QPG parameter, respectively. We choose a squeezing parameter as r = 0.44 and 0.69 and a QPG parameter as κ = 0.46 and 0.75 for the demonstration. The operation condition of the system can be described by using Fig. 1B; the roundtrip number n is one for the single-step gate. We inject a vacuum field into the system as an input resource, and then, by using the displacement and the phase operations in the loop, we prepare three different types of input coherent states ∣α⟩ with α = 0, xin, and ipin (xin, pin ∈ ℝ). Let us call these input states as the vacuum state, X-coherent state, and P-coherent state for later convenience. Ancilla-1 is the squeezed vacuum state. The amount of phase shifts θ1 and θ2 and the transmissivity of the VBS are dynamically set to appropriate values according to the desired squeezing parameter r or QPG parameter κ (Table 1). In all cases, the homodyne angle ϕ1 is 0°. After the input state is processed by the squeezer and the phase shifters, the reflectivity of the BS is changed to zero [BS-(n + 1) in Fig. 1B] so that the wave packet comes out of the loop.
Table 1.

Parameter list of the optical components for the demonstration of Gaussian gates.

The first column identifies the target gate. R(BS-1) is the reflectivity of BS-1 in the virtual diagram shown in Fig. 1B. θ1 and θ2 are the amount of phase shifts, which can also be found in Fig. 1B. g1 is the ratio of the amount of displacement with respect to the value of the measured quadrature. The homodyne angle ϕ1 is not listed here since it is 0° in all cases.

r or κ R(BS-1) (%) θ1 (°) θ2 (°) g (dB)
r = 0.4441001.6
r = 0.6925004.8
κ = 0.4641−32.757.31.6
κ = 0.7525−26.663.44.8

Parameter list of the optical components for the demonstration of Gaussian gates.

The first column identifies the target gate. R(BS-1) is the reflectivity of BS-1 in the virtual diagram shown in Fig. 1B. θ1 and θ2 are the amount of phase shifts, which can also be found in Fig. 1B. g1 is the ratio of the amount of displacement with respect to the value of the measured quadrature. The homodyne angle ϕ1 is not listed here since it is 0° in all cases. Figure 2A describes the experimental setup for the demonstration. Each component is dynamically controlled, and the timing of the control is synchronized at a nanosecond time scale by using a timing controller. Figure 2B shows how such a dynamical and synchronized control is typically achieved in the demonstration. What makes it possible for the optical circuit to perform quantum logic gates in the time domain is to incorporate dynamical and synchronized control of the input switch (Switch-1) and feedforward system into our previous work (), which is the major technical advance in this work. In our setup shown in Fig. 2A, with Switch-2 omitted, the output of the processor is measured by the HD, where the output state is characterized by a tomography process with the local oscillator (LO) phase varied (). The detail of the setup is described in Materials and Methods.
Fig. 2.

Schematic of the experiment.

(A) Optical and electrical setup. Refer to Materials and Methods for the details. EOM, electro-optic modulator; MZI, Mach-Zehnder interferometer; QWP, quarter-wave plate. (B) Typical timing chart of each component. This plot represents the case of a QPG with κ = 0.75 on the X-coherent state (Table 1). The gate sequence starts at around 0 ns by taking the input pulse in the loop part with the reflectivity of the VBS set to zero. The time axis is presented in the laboratory frame, which leads to apparent time deviations in particular for (v) feedforward.

Schematic of the experiment.

(A) Optical and electrical setup. Refer to Materials and Methods for the details. EOM, electro-optic modulator; MZI, Mach-Zehnder interferometer; QWP, quarter-wave plate. (B) Typical timing chart of each component. This plot represents the case of a QPG with κ = 0.75 on the X-coherent state (Table 1). The gate sequence starts at around 0 ns by taking the input pulse in the loop part with the reflectivity of the VBS set to zero. The time axis is presented in the laboratory frame, which leads to apparent time deviations in particular for (v) feedforward. Figure 3 summarizes the results of operation of the squeezing gates with r = 0.44 and 0.69. Figure 3A shows the output Wigner functions of the squeezing gates performed on the X-coherent input state, as an example out of three different input states. The Wigner function of the input state is also shown. According to Eq. 1, the gate operation changes the standard deviations along x and p axes by factors of e and e−, respectively. In the same way, the mean values along x and p axes are also changed by the same factors of e and e−, respectively. Thus, the shapes of the output Wigner functions are squeezed in the direction of the p axis, while the positions are shifted along the x axis, following the matrix transformation of Eq. 1. It can be seen from these plots that the positions and shapes of the ellipses of the experimental outputs agree well with those of the theoretical ones (Fig. 3B), which indicates that the processor performs the squeezing gates as expected. The experimental and theoretical outputs also agree well for the vacuum and P-coherent input states (Supplementary Materials). Figure 3 (C and D) shows the comprehensive evaluation of the gates by quantitatively comparing the experimental and theoretical output states. Figure 3C shows the matrix elements of the squeezing gates S(r) inferred from transformation of mean values of the X- and P-coherent input states. The (1,1) and (2,1) elements are estimated by and , respectively, with the X-coherent input state injected into the processor. Here, ⟨⋯⟩ denotes the mean value. Similarly, the (1,2) and (2,2) elements are estimated by and , respectively, with the P-coherent input state. The matrix elements derived from the experiment coincide well with the theoretical expectations. Figure 3D shows variances along major and minor axes of squeezing ellipses of the output states and tilt angles of the ellipses. The vacuum input state is used for the estimation of the variances. These results also agree well with the theoretical expectations. The fidelity F between the experimental and ideal output states is F = 0.86(2) and 0.75(2) for r = 0.44 and 0.69, respectively, while the realistic model of our setup including optical loss and finite squeezing predicts F = 0.89 and 0.80, respectively. As a whole, the results shown in Fig. 3 together demonstrate that the squeezing gates are properly performed by the processor.
Fig. 3.

Output state characterization of the single-step squeezing gates.

(A) Representative Wigner functions of input and output states of the squeezing gates of the experiment. An x-shaped symbol denotes the center of the distribution, while a solid line is a contour of the relative height of with respect to the center height representing the standard deviation. The convention ℏ = 2 is adopted here. (B) Theoretical predictions for the corresponding Wigner functions. A set of an x-shaped symbol and a solid line is derived from a realistic model with optical loss and finite squeezing, while another set of a circular dot and a dashed line is from an ideal model with no optical loss and infinite squeezing. (C) Matrix elements of a squeezing gate as a function of a target squeezing parameter r. (D) Parameters of a squeezing ellipse as a function of a target squeezing parameter. The variances of squeezing and antisqueezing are shown in the top figure, while the squeezing angles are in the bottom one.

Output state characterization of the single-step squeezing gates.

(A) Representative Wigner functions of input and output states of the squeezing gates of the experiment. An x-shaped symbol denotes the center of the distribution, while a solid line is a contour of the relative height of with respect to the center height representing the standard deviation. The convention ℏ = 2 is adopted here. (B) Theoretical predictions for the corresponding Wigner functions. A set of an x-shaped symbol and a solid line is derived from a realistic model with optical loss and finite squeezing, while another set of a circular dot and a dashed line is from an ideal model with no optical loss and infinite squeezing. (C) Matrix elements of a squeezing gate as a function of a target squeezing parameter r. (D) Parameters of a squeezing ellipse as a function of a target squeezing parameter. The variances of squeezing and antisqueezing are shown in the top figure, while the squeezing angles are in the bottom one. In the same manner, Fig. 4 summarizes the results of the operation of the QPGs with κ = 0.46 and 0.75, showing that the QPGs are also properly performed. Comparing Fig. 4 (A and B) tells us that the positions and shapes of the output Wigner functions agree well with what the theoretical model predicts (refer to the Supplementary Materials for the vacuum and P-coherent input states). The agreement with the model can also be confirmed by Fig. 4 (C and D). The fidelity takes the value of 0.87(3) and 0.77(1) for κ = 0.46 and 0.75, respectively, while the realistic model predicts F = 0.89 and 0.80, respectively.
Fig. 4.

Output state characterization of the single-step QPGs.

(A) Representative Wigner functions of input and output states of the experiment. (B) Theoretical predictions for the corresponding Wigner functions. (C) Matrix elements of a QPG. (D) Parameters of a squeezing ellipse. The results are shown with the same rule as Fig. 3 except that the latter two plots are shown as functions of κ instead of r.

Output state characterization of the single-step QPGs.

(A) Representative Wigner functions of input and output states of the experiment. (B) Theoretical predictions for the corresponding Wigner functions. (C) Matrix elements of a QPG. (D) Parameters of a squeezing ellipse. The results are shown with the same rule as Fig. 3 except that the latter two plots are shown as functions of κ instead of r. The deviations between the experimental results and theoretical expectations found in Figs. 3 and 4 can be attributed to several reasons: suboptimal feedforward gain, systematic errors in the phase shift operations, and optical losses that are not taken into account in the theoretical model. It is worth emphasizing here that the same circuit implements all the above gate operations without any changes of the hardware configuration. Therefore, these results demonstrate the programmability of the present processor and capability of applying arbitrary Gaussian gates to a single-mode input.

Multistep squeezing gate

To show the capability of applying gates multiple times, we perform multistep squeezing gates. For the demonstration, a squeezing gate with r = 0.44 is applied n times (n = 1,2, and 3) to three different input states of the vacuum, X-coherent, and P-coherent states. The operation condition and the experimental setup are the same as that of the previous single-step squeezing gate except that settings for the squeezing process are held over multiple roundtrips of the input wave packet. We choose to apply the same squeezing gate over the multiple steps for the demonstration, rather than to apply different gates for each step, due to the limitation of our current electro-optic modulator (EOM) driver for EOM-2 (refer to the Supplementary Materials for more information). Figure 5 shows the results of the demonstration of the multistep squeezing gates. For the X-coherent input state, the Wigner functions of the experimental output of the n-step squeezer are shown in Fig. 5A, while the model expectations are in Fig. 5B. The Wigner functions for the vacuum and P-coherent input states are shown in the Supplementary Materials. They graphically show that the multistep gates are properly performed. The transformation of the mean value and variance is characterized as shown in Fig. 5 (C and D), which also shows good agreement with the theoretical expectations. The fidelity takes the value of 0.86(2),0.67(3), and 0.48(3) for n = 1,2, and 3, respectively, while the realistic model predicts F = 0.89,0.73, and 0.54, respectively. The deviations in the (1,1) element in Fig. 5C can be attributed to the spatial mode mismatch during the multiple roundtrips in the loop part. The results demonstrate that the loop-shaped processor is able to perform multistep gates on a given single-mode input in a programmable manner. It should be emphasized here that this setup can operate arbitrarily many steps of Gaussian gates in principle, although a squeezing gate is performed up to three times in the demonstration.
Fig. 5.

Output state characterization of the multistep squeezing gates.

(A) Representative Wigner functions of input and output states of the experiment. (B) Theoretical predictions for the corresponding Wigner functions. (C) Matrix elements of a multistep gate. (D) Parameters of a squeezing ellipse. The results are shown with the same rule as Fig. 3 except that the latter two plots are shown as functions of n, the number of repeats.

Output state characterization of the multistep squeezing gates.

(A) Representative Wigner functions of input and output states of the experiment. (B) Theoretical predictions for the corresponding Wigner functions. (C) Matrix elements of a multistep gate. (D) Parameters of a squeezing ellipse. The results are shown with the same rule as Fig. 3 except that the latter two plots are shown as functions of n, the number of repeats.

DISCUSSION

In conclusion, our experiment has demonstrated that the loop-shaped processor is programmable, capable of performing an arbitrary Gaussian gate, and compatible with multistep gate operations. The processor can perform arbitrarily many steps of quantum gates in principle and thus work as a universal processor for a single-mode input, once a cubic phase state is provided as ancillary non-Gaussian resources. It must be emphasized here again that our implementation of quantum gates includes the dynamical feedforward system so that the processor imports and exports photonic quantum states; in contrast, the previous implementations (, ), which were designed specifically for the purpose of one-way quantum computing, did not have an input port or a feedforward system and thus cannot be a general-purpose photonic processor that receives, processes, and returns photonic quantum states. With Switch-2 as depicted in Fig. 1 to send the output state to the free space, the processor could be a core part of various applications. Moreover, the processor can be extended to a multimode input and output straightforwardly by using an additional loop optical path outside the current processor as shown in the original proposal (). In this way, our processor will be evolved into a general-purpose photonic processor, which implements arbitrary photonic quantum information processing.

MATERIALS AND METHODS

Experimental setup

Figure 2A illustrates the experimental setup of the demonstration. The basic working principle and the functionality of a part of the setup are explained in (). The ancillary squeezed vacuum is produced by an optical parametric oscillator (OPO) that is pumped with a continuous-wave laser beam having a wavelength of 430 nm. The other beams have a wavelength of 860 nm. A set of EOM-1 and polarizing BSs (PBSs) forms Switch-1. In a similar way, a set of EOM-2, PBSs, and a quarter-wave plate works as a VBS (refer to the Supplementary Materials for more information). Displacement operation is performed by coupling a phase-modulated coherent field (displacing beam) by a high-reflectivity mirror inside the loop. The carrier field of the displacing beam is undesired and thus removed by introducing a Mach-Zehnder interferometer. EOM-5 on the path of the displacing beam is used for the feedforward operation [(v) in Fig. 2] and the input state preparation [(iii) in Fig. 2], in which a vacuum field from outside the processor is displaced to produce a coherent state. The experimental details around the displacing beam are explained further in the Supplementary Materials. The phase shift operation is performed by EOM-3. The output of the processor is measured by the HD with the homodyne angle varied, where EOM-4 shifts the propagation phase of the LO beam. The 3-dB bandwidth of the HD is 200 MHz. All the dynamical components are synchronized by using a timing controller. Figure 2B represents how the parameter of each component is dynamically and synchronously changed. The rise/fall time of EOMs for switching is ~10 ns. This plot is produced by an auxiliary measurement where classical light fields are introduced into the optical circuit and detected at the output port, so that the light fields probe the dynamical changes of the components. This plot corresponds to the case of a QPG with κ = 0.75 on the X-coherent input state. Since the loop has a roundtrip length of 19.8 m, which corresponds to a roundtrip time of 66 ns, individual optical modes are defined every 66 ns. Thus, the parameter of each optical and electrical component is changed with unit time being 66 ns, except for (iii) state preparation (Supplementary Materials).

Data analysis

For quantum tomography of the output states, the outcome of the HD is acquired by an oscilloscope and then stored in a computer. Among a full stretch of a single time bin, 66 ns, the central 46 ns is used for defining a wave packet mode, while the remaining part is regarded as switching time and thus discarded. We obtain a single sample of a quadrature amplitude by applying a mode function f(t – t0) to the acquired time series. The mode function is defined bywhere γ = 6 × 107/s and t1 = 23 ns. t0 is set to the center of the time bin. The main benefit of using this mode function is its insensitivity to the overall offset (), which is difficult to remove in many experimental situations. The raw value of a quadrature is calibrated by measuring a vacuum field immediately after the main measurement, so that the raw value is converted into a physical value. For each value of the homodyne angle, we sample 1000 independent quadrature amplitudes by repeating the same operations of the processor. Twelve different values of the homodyne angle are chosen in increments of 15°. As a whole, the tomography of the output state is performed using 12 × 1000 samples. The Wigner function is estimated by the method of maximum likelihood with the assumption that the Wigner function is Gaussian (). This assumption is justified in our case because a Gaussian operation never transforms a Gaussian state into a non-Gaussian state. By splitting 1000 samples into 10 individual subsets and estimating the Wigner function for each subset, we obtain the mean values and statistical errors of the parameters of the Wigner functions.

Optical loss and OPO

Let us here describe the model parameters for our experimental setup, which are used to produce the theoretical predictions in Figs. 3 to 5. The quality of the squeezed state from the OPO is modeled in the following way. The amount of optical loss between the OPO and the VBS is measured to be 7%, including the internal loss inside the OPO. The roundtrip optical loss in the loop part is measured to be 6%. The readout loss is measured to be 19%, including the optical loss between the VBS and the HD as well as imperfection in the homodyne detection, such as the mode mismatching effect of the LO field and electrical noise. The squeezing level of the pure squeezed state from the OPO, before suffering from any optical loss, is estimated to be 8.0 dB by measuring the squeezed state from the OPO with the loop part bypassed and considering the losses of 7 and 19% in the path from the OPO to the homodyne detection. The realistic model including the optical loss is derived on the basis of these model parameters. In the ideal model, each value of the loss is set to zero, and the squeezing level of the ancillary state from the OPO is assumed to be infinity.

Continuous-variable gate conventions

In this subsection, let us summarize the conventions concerning the continuous-variable gates. and are the creation and annihilation operators of an optical mode, respectively, which satisfy . We use and for the corresponding quadrature amplitudes, which are defined by and . A squeezing gate is first characterized by a unitary operator . In the Heisenberg picture, the input-output relation of can be written in a matrix form aswhere input and output quadratures are labeled by subscripts “in” and “out,” respectively. Let us then define a QPG by a unitary operator . In a similar way, the input-output relation of can be written as Here, it is worth mentioning that the following matrix identity holdswhere κ = sinhr, cos2θ1 = tanhr, sin2θ1 = − 1/ cosh r, and θ2 = θ1 + π/2. This identity proves that a QPG can be realized by a successive operation of a phase shifter with θ1, a squeezing gate with r, and another phase shifter with θ2. Throughout the paper, we adopt ℏ = 2 so that the variance of quadrature amplitudes is equal to unity for the vacuum state.

Cubic phase gate

We explain the implementation scheme of a cubic phase gate with the loop-based processor here. Let us first define a cubic phase gate by a unitary operator . The quadrature operators of the output of the gate can be written in terms of those of the input as As originally proposed in (), a cubic phase gate can be implemented by a measurement-induced method using a cubic phase state and a squeezed state as ancillae. Figure 6 shows the overview of the method. The input wave packet has interaction with the cubic phase state at BS-1. One outgoing beam from BS-1 is measured by HD-1, where the quadrature operator is measured, while the other is transferred to BS-2. Then, the squeezed state is transferred to BS-2 as an ancilla. One outgoing beam from BS-2 is measured by HD-2, where the specific linear combination of the quadrature operators, , is measured. Let y denote the measurement outcome of HD-2. Here, the parameter μ is determined by the measurement outcome of HD-1 q aswhere R and T are the reflectivity and transmissivity of the BSs, respectively. The other beam from BS-2 becomes the output of the whole gate after the displacement operation is performed on it. The amount of the displacement is as followswhere , , and g3 = 3γ(T − R)/T. As a whole, the output of this circuit can be expressed with respect to the quadratures of input and ancillary states aswhere the subscripts “sq” and “CPG” denote the quadratures of the squeezed state and the cubic phase state, respectively. Here, approaches to zero as the squeezing level of the ancillary squeezed state approaches infinity. Under this condition, approaches to the ideal form (Eq. 7). Similarly, the second term of approaches to zero since it takes the form of the nullifier of a cubic phase state; the reflectivity and transmissivity of the BSs are chosen according to the given γ and the ancillary cubic phase state. The third term can be made arbitrarily small by carefully managing the quality of the two ancillary states, although the variance of is infinity for the ideal ancillary cubic phase state. In summary, this measurement-induced gate approaches the cubic phase gate described in Eqs. 7 and 8 in the limit of pure ancillary states. Apparently, the structure shown in Fig. 6 fits in the loop-shaped architecture (Fig. 1). Note that this method is slightly different from another proposal (), which does not fit in the circuit shown in Fig. 1 as it requires an additional outer loop ().
Fig. 6.

Conceptual schematic of the implementation method of a cubic phase gate.

The input state to this circuit is processed in a measurement-induced method using a cubic phase state and a squeezed state as ancillae. The homodyne angle of HD-2 is adaptively controlled according to the outcome of HD-1. The displacement operation is determined by both of the results from HD-1 and HD-2.

Conceptual schematic of the implementation method of a cubic phase gate.

The input state to this circuit is processed in a measurement-induced method using a cubic phase state and a squeezed state as ancillae. The homodyne angle of HD-2 is adaptively controlled according to the outcome of HD-1. The displacement operation is determined by both of the results from HD-1 and HD-2.
  8 in total

1.  Efficient classical simulation of continuous variable quantum information processes.

Authors:  Stephen D Bartlett; Barry C Sanders; Samuel L Braunstein; Kae Nemoto
Journal:  Phys Rev Lett       Date:  2002-02-14       Impact factor: 9.161

2.  Universal quantum computation with continuous-variable cluster states.

Authors:  Nicolas C Menicucci; Peter van Loock; Mile Gu; Christian Weedbrook; Timothy C Ralph; Michael A Nielsen
Journal:  Phys Rev Lett       Date:  2006-09-13       Impact factor: 9.161

3.  Demonstration of unconditional one-way quantum computations for continuous variables.

Authors:  Ryuji Ukai; Noriaki Iwata; Yuji Shimokawa; Seiji C Armstrong; Alberto Politi; Jun-ichi Yoshikawa; Peter van Loock; Akira Furusawa
Journal:  Phys Rev Lett       Date:  2011-06-17       Impact factor: 9.161

4.  Universal Quantum Computing with Measurement-Induced Continuous-Variable Gate Sequence in a Loop-Based Architecture.

Authors:  Shuntaro Takeda; Akira Furusawa
Journal:  Phys Rev Lett       Date:  2017-09-21       Impact factor: 9.161

5.  Deterministic generation of a two-dimensional cluster state.

Authors:  Mikkel V Larsen; Xueshi Guo; Casper R Breum; Jonas S Neergaard-Nielsen; Ulrik L Andersen
Journal:  Science       Date:  2019-10-18       Impact factor: 47.728

6.  Generation of time-domain-multiplexed two-dimensional cluster state.

Authors:  Warit Asavanant; Yu Shiozawa; Shota Yokoyama; Baramee Charoensombutamon; Hiroki Emura; Rafael N Alexander; Shuntaro Takeda; Jun-Ichi Yoshikawa; Nicolas C Menicucci; Hidehiro Yonezawa; Akira Furusawa
Journal:  Science       Date:  2019-10-18       Impact factor: 47.728

7.  Continuous variable quantum optical simulation for time evolution of quantum harmonic oscillators.

Authors:  Xiaowei Deng; Shuhong Hao; Hong Guo; Changde Xie; Xiaolong Su
Journal:  Sci Rep       Date:  2016-03-10       Impact factor: 4.379

8.  On-demand photonic entanglement synthesizer.

Authors:  Shuntaro Takeda; Kan Takase; Akira Furusawa
Journal:  Sci Adv       Date:  2019-05-17       Impact factor: 14.136

  8 in total
  1 in total

1.  Quantum computational advantage with a programmable photonic processor.

Authors:  Lars S Madsen; Fabian Laudenbach; Mohsen Falamarzi Askarani; Fabien Rortais; Trevor Vincent; Jacob F F Bulmer; Filippo M Miatto; Leonhard Neuhaus; Lukas G Helt; Matthew J Collins; Adriana E Lita; Thomas Gerrits; Sae Woo Nam; Varun D Vaidya; Matteo Menotti; Ish Dhand; Zachary Vernon; Nicolás Quesada; Jonathan Lavoie
Journal:  Nature       Date:  2022-06-01       Impact factor: 69.504

  1 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.