Jaromír Fiurášek

Generation of highly nonclassical quantum states of light is essential for optical quantum information processing and quantum metrology. Given the lack of sufficiently strong nonlinear interactions between optical fields, the commonly employed optical quantum-state preparation schemes are conditional, based on nonlinearity induced by heralding photon number measurement on a part of a multimode squeezed Gaussian state. The development and optimization of such probabilistic quantum-state engineering schemes represents one of the central challenges in current quantum optics. As technology advances and experiments progress toward the detection of higher numbers of photons, the maximization of the heralding probability becomes essential to ensure sufficiently high state-preparation rates. Here, we show that for conditional quantum state preparation schemes based on Gaussian states and photon number measurements, the maximization of the heralding probability can be achieved by finding solutions to a system of polynomial equations, which offers an efficient way to find the optimal configuration and allows us to apply techniques dedicated specifically to solving such systems of equations. Our approach can seamlessly incorporate bounds on available single-mode quadrature squeezing, which is highly experimentally relevant. We mainly consider the generation of finite superpositions of Fock states but show that the approach can be straightforwardly extended to the generation of squeezed superpositions of Fock states. We focus on Gaussian states with vanishing coherent displacements, hence, the conditionally generated states have well-defined photon number parity. We illustrate our general methodology through examples of generation of single-mode states with two and three heralding modes, and generation of two-mode states with two heralding modes.

We address the task of converting n uses of an unknown unitary transformation into a quantum state (i.e., storage) and later retrieval of the transformation. Specifically, we consider the case where the unknown unitary is selected with equal prior probability from two options. First, we prove that the optimal storage strategy involves the sequential application of the n uses of the unknown unitary, and it produces the optimal state for discrimination between the two possible unitaries. Next, we show that incoherent “measure-and-prepare” retrieval achieves the maximum fidelity between the retrieved operation and the original (qubit) unitary. We then identify the retrieval strategy that maximizes the probability of successfully and perfectly retrieving the unknown transformation. In the regime in which the fidelity between the two possible unitaries is large the probability of success scales as P succ = 1 − O ( n − 2 ) , which is a quadratic improvement with respect to the case in which the unitaries are drawn from the entire unitary group U ( d ) with uniform prior probability. Finally, we present an optical experiment for this approach and assess the storage and retrieval quality using quantum tomography of states and processes. The results are discussed in relation to nonoptimal measure-and-prepare strategy, highlighting the advantages of our protocol.

Entanglement of continuous-variable Gaussian states can be distilled by combination of de-Gaussifying operation such as single-photon subtraction and iterative heralded Gaussification. Here we present and analyze a simplified equivalent version of such entanglement distillation protocol, where the Gaussian measurements utilized in heralded Gaussification are eliminated and are absorbed into the preparation of suitable input Gaussian states of the simplified protocol. The simplified scheme contains less detectors and its overall success probability increases in comparison with the original scheme, while producing completely equivalent outputs. Our simplification of the entanglement distillation protocol closely parallels the recently proposed simplification of a scheme for breeding optical single-mode Gottesman–Kitaev–Preskill states (Aghaee Rad et al 2025 Nature 638 912 ). We investigate operation of the simplified entanglement distillation scheme for both pure and mixed input states and clarify how multicopy distillation of Gaussian entanglement emerges in a setup without any heralding Gaussian measurements.

Stellar rank of quantum state of light quantifies the amount of non-Gaussian resources required for its generation. One popular and practical approach to certification of stellar rank is based on measurement of click statistics with an array of binary detectors that can only distinguish the presence and absence of photons. Specifically, it was shown that measurements with an array of m + 1 detectors allow one to certify stellar rank m of approximate Fock state | m 〉 , even when the state is subjected to losses or certain noise. Here we address the question as to how many click detectors are in principle required to certify stellar ranks higher than one. We show that two click detectors arranged in a Hanbury Brown-Twiss setup suffice. Interestingly, detection of stellar ranks higher than one is greatly facilitated by making the total detection efficiency of the detectors sufficiently low but well calibrated. Losses affect the response of the detection scheme in a way that can be exploited to certify stellar rank of higher Fock states. We explicitly construct the corresponding stellar rank witnesses and discuss dependence of the stellar rank thresholds on parameters of the considered setup. Our results reveal that it is possible to certify stellar ranks higher than 1 even with a minimalistic scheme that provides only very coarse-grained information about the photon number statistics of the characterized state.

We investigate fundamental bounds on the ability to determine photon number distribution and other related quantities from tomographically incomplete measurements with an array of M detectors that can only distinguish the absence or presence of photons. We show that the lower and upper bounds on photon number probabilities can be determined by solving a linear program. We present and discuss numerical results for various input states, including thermal states, coherent states, squeezed states, and highly nonclassical single-photon subtracted squeezed vacuum states. Besides photon number probabilities, we also investigate bounds on the parity of photon number distribution that determines the value of Wigner function of the state at the origin of phase space. Moreover, we also discuss estimation of mean photon number as an example of a quantity described by an unbounded operator. Our approach and results can provide quantitative guidance on the number of detection channels required to determine the photon number distribution with a given precision.

Probabilistic heralded Gaussification of quantum states of light is an important ingredient of protocols for distillation of continuous variable entanglement and squeezing. An elementary step of heralded Gaussification protocol consists of interference of two copies of the state at a balanced beam splitter, followed by conditioning on outcome of a suitable Gaussian quantum measurement on one output mode. When iterated, the protocol either converges to a Gaussian state, or diverges. Here we report on experimental investigation of the convergence properties of iterative heralded Gaussification. We experimentally implement two iterations of the protocol, which requires simultaneous processing of four copies of the input state. We utilize the phase-randomized coherent states as the input states, which greatly facilitates the experiment, because these states can be generated deterministically and their mean photon number can easily be tuned. We comprehensively characterize the input and partially Gaussified states by balanced homodyne detection and quantum state tomography. Our experimental results are in good agreement with theoretical predictions and they provide new insights into the convergence properties of heralded quantum Gaussification.

Probabilistic heralded Gaussification of quantum states of light is an important ingredient of protocols for distillation of continuous variable entanglement and squeezing. An elementary step of the heralded Gaussification protocol consists of interference of two copies of the state at a balanced beam splitter, followed by conditioning on the outcome of a suitable Gaussian quantum measurement on one output mode. When iterated, the protocol either converges to a Gaussian state or diverges. Here, we report on experimental investigation of the convergence properties of iterative heralded Gaussification. We experimentally implement two iterations of the protocol, which requires simultaneous processing of four copies of the input state. We utilize the phase-randomized coherent states as the input states, which greatly facilitates the experiment because these states can be generated deterministically, and their mean photon number can easily be tuned. We comprehensively characterize the input and partially Gaussified states by balanced homodyne detection and quantum state tomography. Our experimental results are in good agreement with theoretical predictions and they provide new insights into the convergence properties of heralded quantum Gaussification.