
==== Front
Nature
Nature
Nature
0028-0836
1476-4687
Nature Publishing Group UK London

38720079
7357
10.1038/s41586-024-07357-5
Article
Fusion of deterministically generated photonic graph states
http://orcid.org/0000-0002-3378-7976
Thomas Philip
http://orcid.org/0000-0001-5353-0863
Ruscio Leonardo
http://orcid.org/0000-0002-0327-5646
Morin Olivier olivier.morin@mpq.mpg.de

Rempe Gerhard
https://ror.org/01vekys64 grid.450272.6 0000 0001 1011 8465 Max-Planck-Institut für Quantenoptik, Garching, Germany
8 5 2024
8 5 2024
2024
629 8012 567572
12 9 2023
26 3 2024
© The Author(s) 2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Entanglement has evolved from an enigmatic concept of quantum physics to a key ingredient of quantum technology. It explains correlations between measurement outcomes that contradict classical physics and has been widely explored with small sets of individual qubits. Multi-partite entangled states build up in gate-based quantum-computing protocols and—from a broader perspective—were proposed as the main resource for measurement-based quantum-information processing1,2. The latter requires the ex-ante generation of a multi-qubit entangled state described by a graph3–6. Small graph states such as Bell or linear cluster states have been produced with photons7–16, but the proposed quantum-computing and quantum-networking applications require fusion of such states into larger and more powerful states in a programmable fashion17–21. Here we achieve this goal by using an optical resonator22 containing two individually addressable atoms23,24. Ring25 and tree26 graph states with up to eight qubits, with the names reflecting the entanglement topology, are efficiently fused from the photonic states emitted by the individual atoms. The fusion process itself uses a cavity-assisted gate between the two atoms. Our technique is, in principle, scalable to even larger numbers of qubits and is the decisive step towards, for instance, a memory-less quantum repeater in a future quantum internet27–29.

Using an optical resonator containing two individually addressable atoms in a single cavity, fusion of deterministically generated photonic graph states to create ring and tree graph states with up to eight qubits is demonstrated.

Subject terms

Quantum information
Single photons and quantum effects
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcMain

The characteristics and capabilities of highly entangled graph states1,4 have been widely explored in theoretical quantum information science. These states form a useful subclass of multi-partite entangled states and have the common feature in that they can be represented by a graph comprising vertices and edges (Fig. 1a). A variety of quantum-information protocols have already been implemented in proof-of-principle experiments with graph states made of entangled photons from spontaneous parametric down-conversion (SPDC) sources8,30–32. However, the intrinsically low efficiency of the probabilistic SPDC process remains a substantial obstacle for scalability to large qubit numbers. An alternative and, in principle, deterministic approach using a sequence of single photons emitted from a single memory spin was recognized early on5,7,33 but could not be realized owing to technological shortcomings. The strategy was implemented only recently but with remarkable progress9–15 that finally led to an outperformance of SPDC systems in the achievable number of entangled photons16.Fig. 1 Toolbox for generating photonic graph states.

a, Common examples of graph states. Qubits are represented by vertices, whereas edges connecting them reflect their entanglement topology. b, Two independent graphs can be merged by means of a cavity-assisted fusion gate. c, Correlation measurements on the Bell state |Ψ+⟩=(|01⟩S+|10⟩S)2. Correlations in the measurement bases XX, YY and ZZ certify entanglement and good agreement with the ideal |Ψ+⟩ state. Error bars represent the 1σ standard error. d, Experimental apparatus. Two highly reflective mirrors form an asymmetric high-finesse cavity in which we optically trap two 87Rb atoms at a distance of d = (9 ± 6) μm. The vSTIRAP control laser can be applied either globally or atom-selectively using a high-NA objective. In the process, photons are generated leaving the cavity predominantly through the right mirror. e, Atomic-level scheme illustrating different steps in the protocol. (i) Photon emission from |2, ±2⟩. (ii) Before every photon emission, the qubit is mapped from |1, ±1⟩ to |2, ±2⟩ using Raman lasers at 790 nm. (iii) The same Raman laser system is used to perform single qubit gates on the qubit states |1, ±1⟩ (Raman lasers are not shown in panel d).

Source Data

Although these experiments were limited to elementary photonic graph states such as Greenberger–Horne–Zeilinger (GHZ) and linear cluster states (Fig. 1a), several emitter qubits can, in principle, be combined using quantum logic operations to fully exploit their capabilities17–19,27,28. Once implemented, this would enable architectures that can generate more complex types of graph state for which many powerful quantum-information protocols such as measurement-based quantum computers and quantum repeaters have been proposed2,3,6,27,28. Although recent proposals have successfully identified resource-efficient protocols for such architectures20,21,34, the emitters still need to satisfy some demanding conditions: a suitable energy-level structure for spin–photon entanglement, efficient emission of indistinguishable photons, coherent control of the emitter qubit and high-fidelity entangling gates between emitters. Despite the individual demonstration of these components, no experiment has yet achieved the successful integration of all of them into the same physical system.

Here we demonstrate fusion of photonic graph states produced from two individually addressable atoms in an optical cavity. First, we implement an atom–atom entangling gate based on two-photon interference in the cavity mode23,24. Extending previous work16, we then show that two graph states separately generated from both emitters can be fused into a larger graph (Fig. 1b). In particular, we demonstrate the generation of two important multi-qubit states, namely ring and tree graph states (see Fig. 1a). Both types of state have been identified as valuable resources for protection against qubit loss and/or computational errors in the framework of measurement-based quantum computation and communication25–28,35.

Our experimental setup is schematically shown in Fig. 1d and consists of two 87Rb atoms trapped in a high-finesse optical cavity. Both atoms are positioned at anti-nodes of the cavity mode to ensure strong light–matter coupling with a cooperativity of C = 1.8 and hence to enable efficient generation of single photons by means of a vacuum-stimulated Raman adiabatic passage (vSTIRAP)36. The cooperativity C = g2/(2κγ) is defined in terms of the cavity quantum electrodynamics parameters (g, κ, γ)/2π = (5.4, 2.7, 3.0) MHz. Here g denotes the coupling rate of a single atom to the cavity mode for the D2 line transition ∣F=1,mF=±1⟩↔∣F′=2,mF′=±2⟩, κ is the total cavity-field decay rate and γ is the atomic-polarization decay rate. We use the atomic-state notation |F, mF⟩ (F′,mF′), in which F (F′) denotes the total angular momentum of the ground (excited) state and mF (mF′) its projection along the quantization axis. The latter is given by a magnetic field oriented along the y axis (cavity axis), giving rise to a Zeeman splitting with Larmor frequency ωL/2π = 100 kHz. To minimize crosstalk between the two emitters, the cavity resonance is detuned by Δ/2π = −150 MHz with respect to the F=1,mF=±1↔F′=1,mF′=±1 transition37. The photons are outcoupled from the cavity and directed towards a polarization-resolving detection setup consisting of a polarizing beam splitter and a pair of single-photon detectors. The vSTIRAP control laser can be applied either globally along the x direction acting on both atoms simultaneously or atom-selectively using an acousto-optic deflector (AOD) combined with a high-NA objective on the z axis. Furthermore, atomic-state manipulation such as optical pumping or coherent driving of Raman transitions can be carried out on both atoms simultaneously using global laser beams.

Thanks to the single-atom addressing beam, both atoms serve as independent emitters, each capable of generating individual spin–photon entanglement in parallel. An atom residing in a coherent superposition of the states |2, ±2⟩ can undergo a two-photon transition (vSTIRAP) to |1, ±1⟩, emitting a photon into the cavity mode (see Fig. 1e(i)). We choose the states |0⟩S ≡ |1, +1⟩ and |1⟩S ≡ |1, −1⟩ as the atomic qubit basis (‘S’ for ‘spin’). In the emission process, conservation of angular momentum gives rise to entanglement between the polarization of the photon and the atomic spin state. |0⟩ ≡ |R⟩ and |1⟩ ≡ |L⟩ define the photonic qubit, with R/L corresponding to right/left circular polarization, respectively. This process can be repeated after a Raman transfer from |1, ±1⟩ back to |2, ±2⟩. Using a specifically designed alternating sequence of photon emissions and atomic qubit rotations, elementary graph states such as GHZ or linear cluster states can be obtained16.

The global beam provides the possibility to entangle the two emitters, thus merging the graphs to which they are connected. The underlying mechanism involves two-photon interference in the cavity mode and resembles the type II fusion gate38. Although not strictly identical to fusion in its original form, we here refer to our implementation as a ‘cavity-assisted fusion gate’. The quality of this process depends crucially on the indistinguishability of the photons, which is ensured here by both atoms coupling to the same cavity mode and vSTIRAP control laser. To quantitatively characterize this process, we use the cavity-assisted fusion gate to entangle the two atoms and analyse the correlations in the obtained two-qubit state. Hence we initialize both atoms by optical pumping to the state |2, 0⟩. Next, we carry out the fusion by applying a global vSTIRAP control pulse generating two entangled spin–photon pairs. As the photons interfere in the cavity mode, the which-atom information is erased. Therefore, subsequent measurement in the Z(R/L) basis, with one photon being detected in each detector, projects the atoms onto the Bell state |Ψ+⟩=(|01⟩S+|10⟩S)/2 and indicates the success of the entangling operation. A detection of both photons in the same detector projects the atoms onto a product state (|00⟩S or |11⟩S), which means failure.

For the successful preparation of |Ψ+⟩, we observe strong correlations when measuring the atoms in the bases XX, YY and ZZ (Methods). The probability of each measurement outcome in the different bases is plotted in Fig. 1c. From this, we obtain a state fidelity F=0.915±0.005 with respect to the ideal state. This number varies between 0.851 ± 0.006 and 0.963 ± 0.008, depending on the choice of post-selection criteria for the photon arrival times (Methods). The scenario described above is the simplest case of the fusion mechanism shown in Fig. 1b, in which the emitter qubits do not share a bond with any other qubit before the fusion. The resulting state |Ψ+⟩ can be interpreted as a logical qubit redundantly encoded19,38 in the basis {|0⟩L ≡ |10⟩S, |1⟩L ≡ |01⟩S} (‘L’ for ‘logic’). In the graph-state picture, we express this as a vertex containing two circles. As we will see below, the same principle applies when the two atoms are part of a graph state and do share bonds with other qubits. In this case, the two emitter vertices are merged, preserving the bonds attached to them, as shown in Fig. 1b. If the fusion fails in case both photons end up in the same detector (RR or LL), the emitter vertices are removed from the graph. Although this implies a failure of the protocol, the portion of the graph generated up to this point can still be recovered.

We now use the two-atom Bell pair as a starting point for preparing various photonic graph states. As a first example, we demonstrate the generation of ring graph states consisting of up to eight qubits (Fig. 2e). Basically, we first grow a linear cluster state with the emitters at the ends of the chain, which we then fuse to obtain a ring. The generation steps are shown in Fig. 2a using the graph state representation. We start with a two-atom Bell state |Ψ+⟩, which we obtain from the cavity-assisted fusion. Next, we apply N photon-generation cycles interleaved with π/2 pulses, resulting in a linear cluster state with the atoms at its ends. In each cycle, we first perform the atomic qubit transfer to |2, ±2⟩ (see Fig. 1e(ii)). We then generate one photon from each atom (Fig. 1e(i)), applying the control pulse using the addressing system. The two photons are temporally separated by T ≈ 20 μs to allow enough time for the AOD to direct the addressing beam to the second atom. Afterwards, we perform a π/2 rotation on the atomic qubits by means of the intermediate state |2, 0⟩ using a sequence of Raman pulses (similar to ref. 16). Each cycle lasts 225 μs and has the effect of adding two photonic qubits to the linear cluster state. In the final step, we perform the qubit transfer followed by a photon-production pulse by means of the global beam. This realizes the fusion on the emitter qubits, effectively merging both ends of the chain. For N = 2 and N = 3 photon-generation cycles, this produces either a box-shaped or a hexagon-shaped graph as shown in Fig. 2b,e, respectively. Here, again, the two atoms carry a logical qubit redundantly encoded in |10⟩S and |01⟩S. This specific protocol only produces ring graph states of even parity, that is, an even number of vertices. However, as we show in Methods, ring graph states of odd parity can be obtained simply by adding a global π/4 rotation after the initial fusion gate. In the following, we will focus on the protocol as presented above and demonstrate the generation of the box and hexagon graphs, consisting of four and six vertices, respectively.Fig. 2 Ring graph states.

a, Experimental protocol showing the individual steps such as initialization (purple), vSTIRAP control pulse (blue), photons (red), qubit transfer (green) and qubit rotation (orange). The colour coding follows Fig. 1e. Each of the three lines represents operations carried out on atom 1, atom 2 or globally, that is, by means of a global laser beam acting on both atoms simultaneously. The multi-qubit quantum state is shown below in its graph representation at different stages of the sequence. The full protocol including atom readout takes 1 ms for the box and 1.2 ms for the hexagon graph states. b,e, Graph representation of the generated ‘box’ and ‘hexagon’ graph states. c,f, Measured expectation values for the stabilizers corresponding to the box-shaped/hexagon-shaped graphs in b and e. The colour of the bars (purple and blue) identifies the two sets of stabilizers a and b corresponding to the measurement settings Ma and Mb. d,g, Entanglement witness measurement obtained from coincidences in the measurement settings Ma and Mb. Product correlators Ga (purple) and Gb (blue) and fidelity interval defined by the lower (upper) bound F− (F+) in black. All error bars represent the 1σ standard error.

Source Data

To characterize the experimentally generated state and compare it with the ideal graph state, we measure its corresponding stabilizer operators. The stabilizers to a given graph are defined as Si=Xi∏j∈NiZj, in which Ni is the neighbourhood of vertex i. As the cavity-assisted fusion gate produces vertices that are encoded by two physical qubits, we use the concept of ‘redundantly encoded graph states’19. These are equivalent to regular graph states up to a local unitary transformation on the redundant physical qubits. The stabilizers to the graphs in Fig. 2b,e are shown on the x axis of Fig. 2c,f. To obtain the expectation value of a given stabilizer, we measure coincidences of the corresponding subset of qubits, for which each qubit is detected in either the Z or the X basis. For photonic qubits, the detection basis is set dynamically by means of an electro-optic polarization modulator. The readout of the atomic qubit is realized by an appropriate Raman rotation to set the basis, followed by up to three photon-generation attempts with the detection basis set to R/L (Methods).

The experimentally measured expectation values of the stabilizers are shown in Fig. 2c,f for the box and hexagon graphs, respectively. Furthermore, in the case of an even number of vertices, it is possible to divide the stabilizers into two sets a and b, which can be measured with two local measurement settings Ma and Mb. Similar to in ref. 39, we introduce the operators Ga and Gb as the product ∏i∈a/b(1 + Si)/2 obtained from the measurement setting Ma/b. Their expectation values can be used to compute the fidelity lower bound F−=Ga+Gb−1 and upper bound F+=GaGb (Methods). Both bounds define a constraint for the fidelity given by the inequality F−≤F≤F+.

The results are shown in Fig. 2d,g. For the box-shaped graph, we find the fidelity to fall within the interval given by 0.59±0.03≤F≤0.80±0.02. Because the lower bound of this interval exceeds the threshold of 0.5, we have genuine multi-partite entanglement. In the case of the hexagon graph, the data do not prove genuine entanglement, as the lower bound falls below 0.5. Here we have 0.34−0.07+0.06≤F≤0.67±0.03. We nonetheless observe an overlap with the ideal graph state in terms of the stabilizer expectation values ⟨Si⟩. Also, we emphasize that the true state fidelity is probably higher than the obtained lower bound F−. Alternative characterization methods have been developed for a more precise estimation of the fidelity40,41 but are unsuitable for our current detection setup.

As a second example, we demonstrate the generation of a tree graph state26 consisting of eight qubits (Fig. 3c). In this scenario, we fuse two independent graphs into a larger one. To do this, we first generate two GHZ states, each represented by a star-shaped graph. These will eventually form two branches of the tree graph after merging them using the cavity-assisted fusion gate. The experimental protocol is shown in Fig. 3a. After initialization to the |2, 0⟩ state, each atom emits a photon on successively sending the vSTIRAP control pulse onto the atoms using the addressing system, generating two atom–photon entangled pairs (Fig. 3a, (1)). Two further photon-production cycles are carried out, each cycle consisting of a global Raman transfer to |2, ±2⟩ and a photon-generation pulse (Fig. 3a, (2,3)). Next, after a free evolution time of t0, a π/2 rotation is applied to both atomic qubits simultaneously. At this stage, two GHZ states of the form1 (∣+⟩S∣0++⟩+e−iϕ1,2(t0)∣−⟩S∣1−−⟩)/2

have been generated, each consisting of one atom and three photons. Here ±=0±1/2. Note that two of the photons in the above expression have experienced a Hadamard rotation, which—in the experiment—is absorbed into the setting of the measurement basis. This as well as a π/2 pulse on the atomic qubit has the effect that the respective qubit is ‘pushed out’ from the graph, thus forming a so-called ‘leaf’ qubit (see, for example, ref. 19).Fig. 3 Tree graph states.

a, Experimental protocol similar to Fig. 2a. The length of the individual pulses is the same as for the ring-state protocol. The full protocol including atom readout takes 0.8 ms. b, Quantum phase synchronization between the two branches of the tree. The data show the parity of each four-qubit GHZ state as a function of free evolution time t0 before the π/2 rotation. The parity of both branches is tuned to give a relative π phase difference. The dashed vertical line indicates the working point used for the experiment. Data in light red show the parity when slightly varying the photon separation T (see main text). c, Graph representation of the generated multi-qubit state with the chosen labelling convention. d, Measured stabilizer expectation values similar to Fig. 2. e, Product correlators Ga (purple) and Gb (blue) and fidelity interval (black). All error bars represent the 1σ standard error.

Source Data

The second term in equation (1) carries a phase factor ϕ1,2(t0), which arises from the free evolution of atoms 1 and 2, as denoted by the subscript. We write ϕ1(t0) = 2ωLt0 and ϕ2(t0) = 2ωL(t0 − T) as functions of the π/2 pulse timing given by t0, in which ϕ2(t0) contains the photon separation T as an extra parameter. Figure 3b shows the parity of each of the GHZ states and the oscillating behaviour as a function of t0. We show that, by adjusting t0 and T, we can actively tune this phase to be 0 and π for the two branches, respectively.

In the next step, the two branches are fused into a tree graph. To do this, we apply a global vSTIRAP control pulse, leading to the simultaneous emission of two photons that are detected in the R/L basis. As before, the protocol succeeds if one photon is detected in each detector, that is, in the R and L polarization states, respectively. This step can be thought of as a projection of the atomic qubits on the subspace {|01⟩S, |10⟩S}, given by the operator |01⟩S⟨01|S + |10⟩S⟨10|S. We then obtain the state2 |ψtree⟩=122[|10⟩S(|0++⟩+|1−−⟩)⊗2+|01⟩S(|0++⟩−|1−−⟩)⊗2],

which is an eigenstate to a set of stabilizers corresponding to a tree graph state of depth two, in which the root qubit can be seen as redundantly encoded (Fig. 3c). Owing to the redundant encoding, we again use the modified stabilizers S1 = −Z1Z2 and S2 = X1X2Z3Z6 for the physical qubits of the root vertex. If necessary, the atoms can be disentangled from the photonic state by performing an atom-to-photon state transfer16. In certain cases, however, the protocol may require the emitters to be part of the graph. An example is the one-way quantum repeater28, in which an emitter forms the root qubit of a tree graph.

The measured expectation values for the stabilizers are shown in Fig. 3d. We find all stabilizers to be above 0.7 and some of them close to 1, in good agreement with the ideal state for which ⟨Si⟩ = 1. Furthermore, we are able to prove genuine multi-partite entanglement by collecting eight-qubit coincidences. We find that the entanglement fidelity is constrained by the upper and lower bounds with 0.69−0.05+0.04≤F≤0.85−0.03+0.02, thus exceeding the classical threshold of 0.5.

The fidelities of the generated entangled states are limited by various sources of error. For single-emitter protocols, we have identified spontaneous scattering in the photon-emission process and imperfect Raman rotations as the main error mechanisms16. In this work, we attribute most of the infidelity to the cavity-assisted fusion gate, which is affected by spontaneous scattering as well as imperfect photon indistinguishability. A more detailed discussion can be found in Methods.

The generation of the presented graph states relies on a high overall source-to-detector efficiency, which—in this work—is close to 0.5 for a single photon emission. Hence, with a coincidence rate on the order of one per minute, we can collect hundreds of events in a few hours of measurement (Methods).

In conclusion, we have generated ring graph states of up to 6 (8) logical (physical) qubits and a tree graph state made up of 7 (8) logical (physical) qubits by coupling two emitters by means of a cavity-assisted fusion gate. The latter constitutes the, in our view, decisive step towards scalable architectures of coupled single-photon sources for creating arbitrary photonic graph states. These could be realized with several atom–cavity systems that are embedded in a distributed architecture and connected by optical fibre links3. Alternatively, one could increase the number of emitters within the same cavity, for instance making use of arrays of optical tweezers. Both approaches are conceptually similar, whereas the latter takes advantage of hosting several emitters in the same hardware device. A larger number of emitters would enable tree states of higher depths or repeater graph states, which are proposed as useful tools to overcome photon loss in long-distance transmission lines6,27,28. Similar schemes can be used to generate two-dimensional cluster states to enable fault-tolerant quantum-computing protocols such as one-way or fusion-based quantum computation2,25,35. Finally, yet notably, the photons of the graph state could be individually steered to and stored in a distributed set of heralded quantum memories42, thereby bringing the flying entanglement to a standstill in a material system. In the context of multi-partite quantum networks43, this approach would offer many fascinating possibilities29 beyond those of a two-party quantum-communication link.

Methods

Experimental setup

The apparatus used in our work consists of a single-sided high-finesse cavity in which we optically trap two rubidium atoms. Most experimental details about the setup including the cavity quantum electrodynamics parameters have already been described elsewhere16. In the following, we provide further information that is important for the current work.

The atoms are trapped in a two-dimensional optical standing-wave potential formed by two pairs of counter-propagating laser beams. The first is a retro-reflected laser at a wavelength of λ = 1,064 nm along the x axis. The second propagates inside the cavity mode along the y axis with λ = 772 nm. The atoms are loaded from a magneto optical trap to the cavity centre using a second 1,064-nm running-wave laser. The light scattered by the atom during laser cooling is imaged by means of the objective onto an electron-multiplying charge-coupled device camera to spatially resolve the position of the atoms. After each loading attempt, we find a random number of atoms n at random positions. The experimental control software identifies atom pairs with a suitable relative distance d. If no such atom pair is present, a new loading attempt starts immediately. Otherwise, a tightly focused resonant laser beam, propagating through the objective and steered by the AOD, removes the n − 2 unwanted atoms. The x component of the centre-of-mass position of the atom pair (x2 + x1)/2 is then actively stabilized to the centre of the cavity mode by acting on the relative phase of the 1,064-nm counter-propagating laser beams. The y components y1 and y2 are controlled by modulating the optical power of the 772-nm intra-cavity trap until the atoms are found in a desired position.

Fusion gate and post-selection

For a fusion gate to be successful, two photons have to be detected, as described in the main text. Mathematically, this can be understood by considering two atom–photon entangled states of the form3 ψAP=12F=1,mF=1L−F=1,mF=−1R.

The relative minus sign in the above equation arises from the Clebsch–Gordan coefficients in the two emission paths. Applying the projector ⟨R|⟨L| to the product state |ψAP⟩ ⊗ |ψAP⟩ corresponds to the detection of an R and an L photon, signalling a successful fusion. This leaves us with the |Ψ+⟩ Bell state. Here we implicitly assumed that the two photons occupy the same spatiotemporal mode function. In the experiment, however, their temporal wave packet may not be perfectly indistinguishable, leading to an incomplete erasure of which-path information. Such imperfection can arise from spontaneous scattering by means of the excited state or from unbalanced atom–cavity or atom–laser coupling. This effect becomes visible when post-selecting on the arrival time of the photons. The influence of the arrival time on the fidelity of the atom–atom Bell state is summarized in Extended Data Fig. 1. Panel a shows the intensity profile of the photon temporal wave function as a function of tR,L, with tR and tL being the arrival times of the R-polarized and L-polarized photons produced in the fusion process, respectively. Events in which a photon arrives outside the time interval marked by the dashed lines are discarded. This interval contains about 98% of all single-photon counts. Panel b is a two-dimensional density plot of the number of two-photon events versus arrival times tR and tL. We can see that most events lie in the vicinity of the point tR = tL = 200 ns. The dashed line encloses the region defining the post-selection criteria, which we specify in more detail below. Panel c is a density plot similar to b showing the fidelity as a function of tR and tL. We find that the fidelity is highest near the diagonal of the plot, that is tR ≈ tL. This motivates our choice of the post-selection region enclosed by the dashed line. Pixels for which we did not acquire enough data to compute the fidelity are shown in white. The fidelity is computed using the formula4 F=141+⟨XX⟩+⟨YY⟩−⟨ZZ⟩.

Here XX, YY and ZZ are two-qubit operators consisting of the respective Pauli operators. In panels d–f, we analyse their expectation values ⟨XX⟩, ⟨YY⟩ and ⟨ZZ⟩ as a function of arrival time difference |tR − tL|. We plot the expectation value both for |tR − tL| = τ (orange) and |tR − tL| ≤ τ (purple), that is, the cumulative expectation value. We find all correlators to be in good agreement with the ideal case, for which we expect ⟨XX⟩ = ⟨YY⟩ = 1 and ⟨ZZ⟩ = −1. The high fidelity of the two-atom Bell state is also an indicator of a high photon indistinguishability. The dashed lines indicate the maximum value of τ, that is, |tR − tL| ≤ τ, chosen for the data presented in Fig. 1c.

Post-selection criteria

For the data in Extended Data Fig. 1c, as well as the data presented in the main text, we apply two post-selection steps. The first step consists of restricting the absolute detection time of the photons to a predefined interval of 1 μs width (see dashed lines in Extended Data Fig. 1a). This step applies to both single-photon and two-photon events. The second post-selection condition involves the relative arrival time difference |tR − tL| in the case of two-photon events and therefore only applies to photons generated in the fusion process. The diagonal dashed lines in Extended Data Fig. 1b,c mark the condition |tR − tL| ≤ τmax = 250 ns. Events in which the photons are detected with a relative delay larger than τmax are discarded. In about 80% of experimental runs, the two photons fall within the interval of τmax.

As stated in the main text, the atom–atom Bell-state fidelity ranges between 0.851(6) and 0.963(8). The first number refers to the scenario in which no post-selection on the photon arrival time is applied. The second number is obtained when restricting the photon arrival times to tR,L ≤ 500 ns and |tR − tL| ≤ 20 ns. In this case, the post-selection ratio is about 15%.

The above numbers refer to the scenario in which the atom is initialized to |F = 2, mF = 0⟩ before photon generation. However, in the ring and tree states protocol, the last fusion step consists of a two-photon emission from |F = 2, mF = ±2⟩. In this case, the photon wave packet is slightly longer, as the mF = ±2 Zeeman sublevels couple to different excited states in the emission process. Here we apply the same 1-μs time interval as for the mF = 0 case, as at least 95% of the photon wave packet is enclosed by this window. However, for the two-photon events in the fusion process, we choose a maximum time difference of τmax = 400 ns to accommodate for a post-selection fraction of about 80%, similar to the mF = 0 case.

Atom readout

At the end of the generation sequence for tree and ring graph states, the atomic qubits are still entangled with the photons previously generated. One way to measure the atomic qubits is to perform an atom-to-photon state transfer, as done in ref. 16. Here the qubit is mapped from |F = 1, mF = ±1⟩ to |F = 2, mF = ±1⟩ before photon production. In this way, the qubit is fully transferred to the photon, which can then be measured optically. In this work, however, we chose another technique to measure the atomic qubit. For a Z measurement, we transfer the qubit to |F = 2, mF = ±2⟩ and generate a photon measuring it in the R/L basis. Detecting an R (L) photon projects the atomic qubit onto the state |0⟩S (|1⟩S). When measuring the qubit in X or Y, we set the basis directly on the atomic qubit with a π/2 pulse whose phase is tuned according to the basis. The advantage of this scheme is that it can be repeated until success in the case of photon loss, thus increasing the overall efficiency of the state readout. However, as errors are more likely to occur after many repetitions, we limit the number of attempts to three.

Detailed protocol description

In the following, we will describe the generation protocol for the ring and tree graph states with explicit expressions for each step. In the derivation, we do not explicitly include the free evolution of the atomic qubit. In the experiment, the phases that arise from the qubit oscillation are tracked by measuring the stabilizer operators as a function of certain timing parameters related to, for instance, Raman transfers and photon emissions. Notably, these phases may be tuned for each atom independently by varying the respective time of the photon-production pulse.

Ring states

We first describe the protocol of the ring graph states and choose the pentagon ring as a specific example. The box-shaped and hexagon-shaped graphs are obtained from a similar protocol, only omitting a single π/4 rotation. A sketch of the experimental sequence is given in Extended Data Fig. 2a.

The first step of the protocol is to entangle the two atoms and prepare them in the Bell state |Ψ+⟩=(|01⟩S+|10⟩S)/2. To obtain the pentagon graph, which has an odd number of vertices, we need to apply a global −π/4 pulse. This ‘pushes’ the two qubits apart, forming two separate vertices with an edge between them (Extended Data Fig. 2a, (2)). The corresponding state (omitting normalization constants) reads5 ⟶−π/4∣00⟩S+∣01⟩S+∣10⟩S−∣11⟩S=∣0+⟩S+∣1−⟩S.

Here we have substituted the transformations6 0S⟶−π/4cosθ20S+sinθ21S,1S⟶−π/4−sinθ20S+cosθ21S

and used θ = −π/4, as well as cosπ8=2+22 and sinπ8=2−22. Subsequently, each atom emits a photon, giving7 ⟶PP∣0⟩S∣0⟩(∣0⟩∣0⟩S+∣1⟩∣1⟩S)+∣1⟩S∣1⟩(∣0⟩∣0⟩S−∣1⟩∣1⟩S)=(∣0⟩S∣0⟩+∣1⟩S∣1⟩)∣0⟩∣0⟩S+(∣0⟩S∣0⟩−∣1⟩S∣1⟩)∣1⟩∣1⟩S,

followed by a π/2 rotation on the atomic qubits:8 ⟶π/2(∣+⟩S∣0⟩+∣−⟩S∣1⟩)∣0⟩∣+⟩S+(∣+⟩S∣0⟩−∣−⟩S∣1⟩)∣1⟩∣−⟩S=(∣0⟩S∣+⟩+∣1⟩S∣−⟩)∣0⟩∣+⟩S+(∣0⟩S∣−⟩+∣1⟩S∣+⟩)∣1⟩∣−⟩S,

which is equal to a four-qubit linear cluster state with the atoms at both ends of the chain. Note that the (global) π/2 pulse affects only the spin component of the multi-qubit state. We perform another photon production on both spins and obtain9 ⟶PP(∣0⟩S∣0+⟩+∣1⟩S∣1−⟩)∣0⟩(∣0⟩∣0⟩S+∣1⟩∣1⟩S)+(∣0⟩S∣0−⟩+∣1⟩S∣1+⟩)∣1⟩(∣0⟩∣0⟩S−∣1⟩∣1⟩S).

We apply a Z gate to qubit 6 and a Hadamard to qubits 1 and 6 (indices run from left to right).10 ⟶H1⊗Z6⊗H6(∣+⟩S∣0+⟩+∣−⟩S∣1−⟩)∣0⟩(∣−⟩∣0⟩S+∣+⟩∣1⟩S)+(∣+⟩S∣0−⟩+∣−⟩S∣1+⟩)∣1⟩(∣+⟩∣0⟩S+∣−⟩∣1⟩S).

Now we perform the fusion operation on qubits 1 and 6 and obtain11 ⟶Fusion122(∣10⟩S(∣0+0+⟩+∣1−0+⟩+∣0−1−⟩+∣1+1−⟩)+∣01⟩S(∣0+0−⟩−∣1−0−⟩+∣0−1+⟩−∣1+1+⟩))=122(∣0⟩L(∣0+0+⟩+∣1−0+⟩+∣0−1−⟩+∣1+1−⟩)+∣1⟩L(∣0+0−⟩−∣1−0−⟩+∣0−1+⟩−∣1+1+⟩)).

Here we have moved the second spin qubit to the first position and reintroduced the logical qubit encoding using |0⟩L and |1⟩L. Furthermore, we have added a normalization constant. The above expression represents the state that corresponds to the graph shown in Extended Data Fig. 2c. The measured stabilizer expectation values are shown in Extended Data Fig. 2b.

Tree states

We now describe the experimental protocol for generating the target state of the form12 |ψtree⟩=122[|0⟩(|0++⟩+|1−−⟩)⊗2+|1⟩(|0++⟩−|1−−⟩)⊗2].

We start by preparing both atoms in the |F = 2, mF = 0⟩ state, followed by three sequential photon-production cycles on each atom in parallel. From this, we obtain the tensor product of two GHZ states, each consisting of one atom and three photons (see also ref. 16). Omitting normalization constants, we can write the state as13 ∣ψ⟩=(∣0⟩S∣000⟩+∣1⟩S∣111⟩)⊗(∣0⟩S∣000⟩−∣1⟩S∣111⟩).

Note that the second term carries a relative minus sign with respect to the first term. This is reflected in the parity measurement shown in Fig. 3b. We now perform a Hadamard gate on all qubits except qubits 2 and 6 (indices run from left to right) and obtain14 ⟶(∣+⟩S∣0++⟩+∣−⟩S∣1−−⟩)⊗(∣+⟩S∣0++⟩−∣−⟩S∣1−−⟩).

For the atoms, the Hadamard is carried out with a Raman laser (see Fig. 1e), whereas for the photons, it is absorbed into the setting of the measurement basis.

We now merge both branches into one larger graph state by applying the fusion gate. Hence we generate two photons from the atoms with the global STIRAP control laser. Detecting one photon in R and one in L effectively projects the atoms onto the subspace {|01⟩S, |10⟩S}.15 ⟶∣01⟩S⟨01∣S+∣10⟩S⟨10∣S(∣10⟩S+∣01⟩S)∣0++⟩⊗2+(∣10⟩S−∣01⟩S)∣0++⟩∣1−−⟩+(∣10⟩S−∣01⟩S)∣1−−⟩∣0++⟩+(∣10⟩S+∣01⟩S)∣1−−⟩⊗2=∣10⟩S(∣0++⟩+∣1−−⟩)⊗2+∣01⟩S(∣0++⟩−∣1−−⟩)⊗2.

For convenience, we have moved the second spin qubit to the first position in the above expression, which allows us to express the two atoms as a logical qubit encoded in the basis {|0⟩L ≡ |10⟩S, |1⟩L ≡ |01⟩S}. Adding a normalization constant, we can then write the final state as16 |ψtree⟩=122[|0⟩L(|0++⟩+|1−−⟩)⊗2+|1⟩L(|0++⟩−|1−−⟩)⊗2].

This is equal to the expression in equation (12), with the only difference being that the root qubit is now redundantly encoded by the two atoms. Alternatively, it would be possible to remove one of the atoms from the state by an X basis measurement.

Coincidence rate

For each multi-qubit state, the typical generation and detection rate is between 0.4 and 2.3 coincidences per minute. The total number of events as well as the total measurement time are summarized in Extended Data Table 1 for each graph state generated in this work. These numbers include all post-selection steps as described above.

Entanglement witness and fidelity bounds

To quantify the agreement between the experimentally produced multi-photon state and the target state, we use an entanglement witness. This has the advantage that we can derive a lower bound of the fidelity without the need for full quantum-state tomography. The fidelity of a density matrix ρ with respect to the target state |ψ⟩ is defined as17 F=Tr{ρψψ}.

Using the stabilizers, we can express the projector to the target state as18 ψψ=∏i1+Si2=∏i∈a1+Si2∏j∈b1+Sj2=Ga⋅Gb.

Here we have written the projector as a product of two terms Ga and Gb associated with two sets of stabilizers a and b. Each set a/b can be measured with a single local measurement setting Ma/Mb. These only involve measurements in the X or Z basis for every qubit. We can then write the projector in terms of Ga and Gb, giving19 ψψ=Ga⋅Gb=Ga+Gb−1+1−Ga1−Gb

As the stabilizers Si take the values +1 or −1, the product terms Ga and Gb are either 1 or 0. We conclude that (1 − Ga)(1 − Gb) is non-negative. Omitting this term, we find the lower bound20 F−≡Ga+Gb−1≤F.

The above expression is applicable if the stabilizers can be divided into two sets a and b, each of which can be measured with a single measurement setting (Ma and Mb). In the context of our experiment, this applies to tree graph states as well as ring graph states of even parity, that is, an even number of vertices. To the best of our knowledge, there is no equivalent method for ring graph states of odd parity, such as the pentagon graph, and a fidelity lower bound cannot be derived.

We can further derive a fidelity upper bound based on the terms Ga and Gb. First, for any pure state |ψ⟩, we have21 ψGaGbψ≤ψGaGa†ψψGb†Gbψ,

by direct application of the Cauchy–Schwarz inequality. The terms (1 + Si)/2 are projectors, because Si2=1 and therefore22 1+Si22=1+2Si+Si24=1+Si2.

By construction, the stabilizers Si commute and therefore the projectors (1 + Si)/2 commute as well. Hence, because Ga/b are products of commuting projectors, Ga and Gb themselves are also projectors:23 Ga/b2=∏i∈a/b1+Si22=∏i∈a/b1+Si22=∏i∈a/b1+Si2=Ga/b.

Equation (21) can then be simplified as ⟨ψ∣GaGb∣ψ⟩≤⟨ψ∣Ga∣ψ⟩⟨ψ∣Gb∣ψ⟩.

Then, to generalize to mixed states, we write the mixed state ρ as a linear combination of pure states, that is, ρ=∑kpkψkψk, and apply the above inequality to each of them:24 ⟨GaGb⟩=∑kpk⟨ψk∣GaGb∣ψk⟩≤∑kpk⟨ψk∣Ga∣ψk⟩⟨ψk∣Gb∣ψk⟩.

We identify the right term as a scalar product of two vectors and again use the Cauchy–Schwarz inequality25 ∑kpkψkGaψkpkψkGbψk≤∑kpkψkGaψk∑k′pk′ψk′Gbψk′,

which shows the upper bound of the fidelity26 F=⟨GaGb⟩≤⟨Ga⟩⟨Gb⟩≡F+.

In the next section, we will use both fidelity bounds for a comparison between the experimental data and the expected fidelity.

Estimation of errors

In our previous work16, we identified some error mechanisms present in our system. For single-emitter protocols, the main error sources are spontaneous scattering in the photon-emission process (about 1% per photon) and imperfect Raman rotations (about 1% per π/2 pulse). In the following, we discuss several more mechanisms that could negatively affect the fidelity. In some cases, the effect of these mechanisms on the fidelity of multi-qubit entangled states is difficult to quantify because of the complexity of the entanglement topology and the protocols to generate it. Furthermore, measuring the fidelity of multi-qubit states is a non-trivial task and our measurement setup only allows us to extract a lower and an upper bound of the fidelity.

Fusion gate

For the two-emitter protocols developed in this work, the cavity-assisted fusion gate is probably the largest source of error. As shown in the main text, this mechanism can be used to prepare the |Ψ+⟩ Bell state with a fidelity ranging between 0.85 and 0.96, depending on how strictly we post-select on the arrival time of the photons. The fact that the fidelity decreases with a larger arrival time difference τ (see Extended Data Fig. 1) can be explained by an imperfect indistinguishability of the photons involved in the fusion process. For the standard value of τmax = 250 ns, the fidelity of the |Ψ+⟩ Bell state is 0.92. This number includes state readout of the two atoms, each of which is expected to introduce an error similar to a photon emission (roughly 1%). We conclude that the infidelity from the fusion process is on the order of 6%.

Decoherence

Another potential source of infidelity is atomic decoherence caused by magnetic-field noise or intensity fluctuations of the optical-trapping beams. We have measured the coherence time of the atomic qubit T2 to be approximately 1 ms. However, the atomic qubit is largely protected by a dynamical decoupling mechanism that is built into the protocol16, thereby extending the coherence time. The exact extent to which this mechanism takes effect depends on the specific timing parameters in the sequence and the frequency range in which the noise sources are most dominant (for example, magnetic-field fluctuations). Therefore, it is difficult to quantify how much the decoherence translates into infidelity of the final graph state. Furthermore, different types of graph state are more or less susceptible to noise44. It is therefore not straightforward to theoretically model the role of decoherence in the fidelity of the final multi-partite entangled state.

Qubit leakage

During the protocol, the emitter qubits are continuously transferred between different atomic states. These states are |1, ±1⟩, |2, ±2⟩ and |2, 0⟩, in which we again write the state as |F, mF⟩ with the quantum numbers F and mF. However, there seems to be a low probability that, during the emission process, the atom undergoes a transition to |1, 0⟩ (instead of |1, ±1⟩). This is readily explained by and consistent with the finding of spontaneous scattering during the vSTIRAP process, but may equally result from a contamination of σ+/σ− polarization components in the vSTIRAP control laser. The latter is in turn caused by either an imperfect polarization setting or longitudinal polarization components owing to the tight focus of the beam. The unwanted σ+/σ− components couple to the F′=1,mF′=±1 states and can thus drive a two-photon transition to |F = 1, mF = 0⟩. This process results in the atom leaving the qubit subspace but, unfortunately, such an event remains undetected. If the protocol resumes with a Raman π/2 pulse, the parasitic population in |1, 0⟩ is then partly transferred to |2, ±1⟩, as the corresponding transitions have the same resonance frequency. A subsequently emitted photon will then yield a random measurement outcome, which is detrimental to the fidelity of the state.

The leakage mechanism described above is difficult to quantify, mainly because our experiment lacks an mF-selective state readout. We do however estimate that the longitudinal polarization components of the addressing beam have a relative amplitude on the order of about 1%, contributing to each single-photon emission. For the global beam, this effect is negligible owing to a larger focus.

Other sources of error

Other sources of error include drifts of the optical fibres, such as for the Raman beam, the global and addressing vSTIRAP beam or the optical traps, as well as the magnetic field. Furthermore, the position of the atoms is not fixed but varies from one loading attempt to another. In this work, we chose position criteria that are less strict than those in ref. 16, to increase the data rate. In combination with the drifts mentioned above, this leads to a variance in coupling between the atoms and the cavity, as well as the atoms and different laser beams. As a consequence, this may affect the fidelity of different processes, such as the fusion gate or Raman transfers. Furthermore, a drift of the magnetic field or the light shift induced by the optical trap can influence the phase of the atomic qubits at different stages of the protocol.

A way to reduce the overall infidelity would be to increase the cooperativity C. This would reduce the effect of spontaneous scattering, improve photon indistinguishability and thereby increase the fidelity of the fusion process and partly mitigate the qubit-leakage error. Photon emission through the D1 line of rubidium would have a similar effect, owing to a larger hyperfine splitting in the 52P1/2 excited state. Another strategy to improve the system would be a better control of the atom positions by using more advanced trapping techniques, such as optical tweezers. This would greatly reduce all errors associated with the variance of the atom positions. It would also allow longer trapping times and therefore higher data rates.

Error model

As an (oversimplified) ansatz to estimate the combined effect of the error mechanisms described above, we write the density matrix as a mixture of the ideal density matrix and white noise. This is a common approach to investigate, for instance, the robustness of entanglement witnesses against noise (see, for example, ref. 45). The density matrix then reads27 ρexp=1−pnoiseρideal+pnoise12n,

in which pnoise is the total error probability, ρideal is the ideal density matrix, 1 is the identity matrix and n is the number of qubits. We decompose pnoise into the different error contributions and write28 pnoise=1−1−pPNP1−pRNR1−pFNF.

Here pP denotes the probability of spontaneous scattering during photon emission, pR the error probability during a Raman rotation, pF the error probability for the fusion process and NP, NR and NF are the respective number of operations in the protocol. Note that we do not include mechanisms such as decoherence or qubit leakage in the above formula, as we are unable to assign a value to a specific step of the protocol.

In Extended Data Table 2, we compare the fidelity model to the measured lower and upper bounds as defined by equation (20) and equation (26), respectively. For the tree and box graph states, the predicted fidelities Fmodel are found to fall between the measured bounds, as expected. For the hexagon graph, Fmodel falls slightly above the upper bound but is still consistent with it when taking into account the statistical uncertainty (less than one standard deviation). As mentioned earlier, the model does not include the effect of qubit leakage, decoherence and drifts of, for instance, the magnetic field or optical fibres. Hence, it is likely that the predicted fidelities are slightly overestimated.

Online content

Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at 10.1038/s41586-024-07357-5.

Supplementary information

Peer Review File

Source data

Source Data Fig. 1

Source Data Fig. 2

Source Data Fig. 3

Extended data figures and tables

Extended Data Fig. 1 Atom–atom entanglement by means of the cavity-assisted fusion gate.

a, Histogram showing the total photon count rate as a function of tR and tL, in which tR/L is the arrival time of the right-hand-polarized/left-hand-polarized photon generated in the fusion process. Only events in which both photons were detected are shown. Dashed lines mark the acceptance window for post-selection. b, Density plot of the number of counts as a function of tR and tL. c, Density plot of the fidelity as a function of tR and tL. d–f, Expectation values of the correlators XX, YY and ZZ as a function of photon arrival time difference. The orange line shows the correlator for time difference |tR − tL| = τ, whereas the purple line is the cumulative correlator, meaning for events in which |tR − tL| ≤ τ. The dashed lines mark the maximum τ we choose for Fig. 1c. Error bars represent the 1σ standard error.

Extended Data Fig. 2 Protocol for the generation of the pentagon graph.

A two-atom graph state is obtained from the cavity-assisted fusion gate, followed by a −π/4 pulse. A chain is grown along one dimension using photon emissions and π/2 rotations on the atomic qubits. Both ends of the chain are merged to form a ring. Error bars represent the 1σ standard error.

Extended Data Table 1 Coincidence rate

Coincidence rate

Coincidence-rate statistics for the generated box, pentagon, hexagon and tree graph states.

Extended Data Table 2 Error estimation

Error estimation

We compare the predicted fidelities with the measured upper and lower bounds of the fidelity. For our model, we use the following error probabilities for the different steps: pP = 0.02, pR = 0.01 and pF = 0.06.

Extended data

is available for this paper at 10.1038/s41586-024-07357-5.

Supplementary information

The online version contains supplementary material available at 10.1038/s41586-024-07357-5.

Acknowledgements

P.T. thanks G. Styliaris for valuable discussions. This work was supported by the Bundesministerium für Bildung und Forschung through the Verbund QR.X (16KISQ019), by the Deutsche Forschungsgemeinschaft under Germany’s Excellence Strategy - EXC-2111 - 390814868 and by the European Union’s Horizon Europe research and innovation programme through the project QIA-Phase 1 under grant agreement no. 101102140. The opinions expressed in this document reflect only the author’s view and reflects in no way the European Commission’s opinions. The European Commission is not responsible for any use that may be made of the information it contains.

Author contributions

P.T. and L.R. performed the experiments. P.T. and O.M. analysed the data. G.R. and O.M. supervised the project. All authors were involved in the discussion of the results and writing of the manuscript.

Peer review

Peer review information

Nature thanks Barak Dayan, Peter Lodahl and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Funding

Open access funding provided by Max Planck Society.

Data availability

The datasets generated and/or analysed during this study are available at 10.5281/zenodo.10717770 (ref. 46). Source data are provided with this paper.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Briegel HJ Raussendorf R Persistent entanglement in arrays of interacting particles Phys. Rev. Lett. 2001 86 910 10.1103/PhysRevLett.86.910 11177971
2. Raussendorf R Briegel HJ A one-way quantum computer Phys. Rev. Lett. 2001 86 5188 10.1103/PhysRevLett.86.5188 11384453
3. Barrett SD Kok P Efficient high-fidelity quantum computation using matter qubits and linear optics Phys. Rev. A 2005 71 060310(R) 10.1103/PhysRevA.71.060310
4. Hein, M. et al. in Quantum Computers, Algorithms and Chaos (eds Casati, G. Shepelyansky, D. L., Zoller, P. & Benenti, G.) 115–218 (IOS Press, 2006).
5. Lindner NH Rudolph T Proposal for pulsed on-demand sources of photonic cluster state strings Phys. Rev. Lett. 2009 103 113602 10.1103/PhysRevLett.103.113602 19792371
6. Azuma K Tamaki K Lo H-K All-photonic quantum repeaters Nat. Commun. 2015 6 6787 10.1038/ncomms7787 25873153
7. Wilk T Webster SC Kuhn A Rempe G Single-atom single-photon quantum interface Science 2007 317 488 490 10.1126/science.1143835 17588899
8. Walther P Experimental one-way quantum computing Nature 2005 434 169 176 10.1038/nature03347 15758991
9. Schwartz I Deterministic generation of a cluster state of entangled photons Science 2016 354 434 437 10.1126/science.aah4758 27608669
10. Istrati D Sequential generation of linear cluster states from a single photon emitter Nat. Commun. 2020 11 5501 10.1038/s41467-020-19341-4 33127924
11. Besse J-C Realizing a deterministic source of multipartite-entangled photonic qubits Nat. Commun. 2020 11 4877 10.1038/s41467-020-18635-x 32985501
12. Yang CW Yu Y Li J Bao X-H Pan J-W Sequential generation of multiphoton entanglement with a Rydberg superatom Nat. Photonics 2022 16 658 661 10.1038/s41566-022-01054-3
13. Ferreira, V. S., Kim, G., Butler, A., Pichler, H. & Painter, O., Deterministic generation of multidimensional photonic cluster states with a single quantum emitter. Nat. Phys. 10.1038/s41567-024-02408-0 (2024).
14. Cogan D Su Z-E Kenneth O Gershoni D Deterministic generation of indistinguishable photons in a cluster state Nat. Photonics 2023 17 324 329 10.1038/s41566-022-01152-2 37064524
15. Coste N High-rate entanglement between a semiconductor spin and indistinguishable photons Nat. Photonics 2023 17 582 587 10.1038/s41566-023-01186-0
16. Thomas P Ruscio L Morin O Rempe G Efficient generation of entangled multiphoton graph states from a single atom Nature 2022 608 677 681 10.1038/s41586-022-04987-5 36002484
17. Economou SE Lindner N Rudolph T Optically generated 2-dimensional photonic cluster state from coupled quantum dots Phys. Rev. Lett. 2010 105 093601 10.1103/PhysRevLett.105.093601 20868157
18. Russo A Barnes E Economou SE Generation of arbitrary all-photonic graph states from quantum emitters New J. Phys. 2019 21 055002 10.1088/1367-2630/ab193d
19. Hilaire P Vidro L Eisenberg HS Economou SE Near-deterministic hybrid generation of arbitrary photonic graph states using a single quantum emitter and linear optics Quantum 2023 7 992 10.22331/q-2023-04-27-992
20. Li B Economou SE Barnes E Photonic resource state generation from a minimal number of quantum emitters npj Quantum Inf. 2022 8 11 10.1038/s41534-022-00522-6
21. Löbl, M. C., Paesani, S. & Sørensen, A. S., Loss-tolerant architecture for quantum computing with quantum emitters. Quantum 8, 1302 (2024).
22. Reiserer A Rempe G Cavity-based quantum networks with single atoms and optical photons Rev. Mod. Phys. 2015 87 1379 10.1103/RevModPhys.87.1379
23. Casabone B Heralded entanglement of two ions in an optical cavity Phys. Rev. Lett. 2013 111 100505 10.1103/PhysRevLett.111.100505 25166644
24. Langenfeld S Thomas P Morin O Rempe G Quantum repeater node demonstrating unconditionally secure key distribution Phys. Rev. Lett. 2021 126 230506 10.1103/PhysRevLett.126.230506 34170169
25. Bartolucci S Fusion-based quantum computation Nat. Commun. 2023 14 912 10.1038/s41467-023-36493-1 36805650
26. Varnava M Browne DE Rudolph T Loss tolerance in one-way quantum computation via counterfactual error correction Phys. Rev. Lett. 2006 97 120501 10.1103/PhysRevLett.97.120501 17025947
27. Buterakos D Barnes E Economou SE Deterministic generation of all-photonic quantum repeaters from solid-state emitters Phys. Rev. X 2017 7 041023
28. Borregaard J One-way quantum repeater based on near-deterministic photon-emitter interfaces Phys. Rev. X 2020 10 021071
29. Epping M Kampermann H Bruß D Large-scale quantum networks based on graphs New J. Phys. 2016 18 053036 10.1088/1367-2630/18/5/053036
30. Yao X-C Experimental demonstration of topological error correction Nature 2012 482 489 494 10.1038/nature10770 22358838
31. Bell BA Experimental demonstration of a graph state quantum error-correction code Nat. Commun. 2014 5 3658 10.1038/ncomms4658 24752224
32. Zhong H-S 12-photon entanglement and scalable scattershot boson sampling with optimal entangled-photon pairs from parametric down-conversion Phys. Rev. Lett. 2018 121 250505 10.1103/PhysRevLett.121.250505 30608840
33. Schön C Solano E Verstraete F Cirac JI Wolf MM Sequential generation of entangled multiqubit states Phys. Rev. Lett. 2005 95 110503 10.1103/PhysRevLett.95.110503 16196992
34. Zilk, F. et al. in Proc. 2022 IEEE/ACM Third International Workshop on Quantum Computing Software (QCS) 57–67 (IEEE, 2023).
35. Bell TJ Pettersson LA Paesani S Optimizing graph codes for measurement-based loss tolerance PRX Quantum 2023 4 020328 10.1103/PRXQuantum.4.020328
36. Morin O Körber M Langenfeld S Rempe G Deterministic shaping and reshaping of single-photon temporal wave functions Phys. Rev. Lett. 2019 123 133602 10.1103/PhysRevLett.123.133602 31697544
37. Langenfeld S Morin O Körber M Rempe G A network-ready random-access qubits memory npj Quantum Inf. 2020 6 86 10.1038/s41534-020-00316-8
38. Browne DE Rudolph T Resource-efficient linear optical quantum computation Phys. Rev. Lett. 2005 95 010501 10.1103/PhysRevLett.95.010501 16090595
39. Tóth G Gühne O Detecting genuine multipartite entanglement with two local measurements Phys. Rev. Lett. 2005 94 060501 10.1103/PhysRevLett.94.060501 15783712
40. Flammia ST Liu Y-K Direct fidelity estimation from few Pauli measurements Phys. Rev. Lett. 2011 106 230501 10.1103/PhysRevLett.106.230501 21770489
41. Tiurev K Sørensen AS Fidelity measurement of a multiqubit cluster state with minimal effort Phys. Rev. Res. 2022 4 033162 10.1103/PhysRevResearch.4.033162
42. Brekenfeld M Niemietz D Christesen JD Rempe G A quantum network node with crossed optical fibre cavities Nat. Phys. 2020 16 647 651 10.1038/s41567-020-0855-3
43. Pompili M Realization of a multinode quantum network of remote solid-state qubits Science 2021 372 259 264 10.1126/science.abg1919 33859028
44. Dür W Briegel H-J Stability of macroscopic entanglement under decoherence Phys. Rev. Lett. 2004 92 180403 10.1103/PhysRevLett.92.180403 15169473
45. Tóth G Gühne O Entanglement detection in the stabilizer formalism Phys. Rev. A 2005 72 022340 10.1103/PhysRevA.72.022340
46. Thomas, P., Ruscio, L., Morin, O. & Rempe, G. Data of “Fusion of deterministically generated photonic graph states”. 10.5281/zenodo.10717770 (2024).
