Shortcut to Multipartite Entanglement Generation: A Graph Approach to Boson Subtractions

We propose a graph method for systematically searching for schemes that can generate multipartite entanglement in linear bosonic systems with heralding. While heralded entanglement generation offers more tolerable schemes for quantum tasks than postselected ones, it is generally more challenging to find appropriate circuits for multipartite systems. We show that our graph mapping from boson subtractions provides handy tactics to overcome the limitations in circuit designs. We present a practical strategy to mitigate the limitation through the implementation of our graph technique. Our physical setup is based on the sculpting protocol, which utilizes an $ N$ spatially overlapped subtractions of single bosons to convert Fock states of evenly distributed bosons into entanglement. We have identified general schemes for qubit N-partite GHZ and W states, which are significantly more efficient than previous schemes. In addition, our scheme for generating the superposition of $N=3$ GHZ and W entangled states illustrates that our approach can be extended to derive more generalized forms of entangled states. Furthermore, we have found an N-partite GHZ state generation scheme for qudits, which requires substantially fewer particles than previous proposals. These results demonstrate the power of our approach in discovering optimized solutions for the generation of intricate heralded entangled states. As a proof of concept, we propose a linear optical scheme for the generation of the Bell state by heralding detections. We expect our method to serve as a promising tool in generating diverse entanglement.


I. INTRODUCTION
Entanglement is an essential aspect of quantum information science, investigation of which has resulted in new fundamental understandings about the nature [1].Practically, entanglement is recognized as a valuable resource in the field of quantum information processing, with potential applications in areas such as cryptography [2,3] and computing [4,5].To study and utilize entanglement, it is a prerequisite to find reliable procedures to construct entangled quantum systems.One of the promising approaches to this task is to exploit the indistinguishability of quantum particles [6].Various works suggested theoretical and experimental entanglement generation schemes based on the identicality of particles and postselection.Along these lines, Refs.[7][8][9] showed that two spatially overlapped indistinguishable particles can carry bipartite entanglement.The quantitative relation of particle indistinguishability and spatial overlap to the bipartite entanglement was rigorously analyzed in Refs.[10][11][12].For the case of multipartite entanglement, schemes for GHZ and W states with identical particles have been theoretically suggested [13][14][15][16][17][18][19][20][21][22] and experimented [23,24].Ref. [22] presented a comprehensive graph-theoretic approach to embrace the schemes to generate the entanglement of identical particles in linear quantum networks (LQNs) with postselection.
On the other hand, considering that the postselected schemes are highly sensitive to particle loss in circuits [25,26] and the multipartite correlations can be created by the postselection bias [27,28], there have been several attempts to generate the entanglement of identical particles without postselection.Specifically, heralded generation of entangled states of photons was studied for bipartite [29][30][31][32] and multipartite systems [33][34][35][36][37][38].While postse-lected schemes involve identifying desired outcomes after conducting operations, heralded schemes use real-time signals or measurements in ancillary spatial modes as "heralds" to sort out successful runs during the experiment (see, e.g., Ref. [37] for a more detailed explanation).This property of heralded operations renders more tolerable schemes from photon loss [25], however it is usually more challenging to find proper circuits to obtain the heralded entanglement of an arbitrary N -partite systems than the postselected ones [37].
In this work, we introduce a systematic method to overcome the difficulty of obtaining heralded entanglement generation schemes for an arbitrary Npartite system.Our method employs the sculpting protocol introduced in Ref. [39], which generates an N -partite entangled state by applying N single-boson subtraction operators (which we name the "sculpting operator") to a 2N boson initial state.By setting the initial state to have the even distribution of the bosons in different 2N states, the spatially overlapped sculpting operation (i..e, a singleboson subtraction operator Â(l) is a summation of subtractions on different spatial modes) generates the N -partite entanglement.And various sculpting operators result in different entangled states.Since linear bosonic systems with heralding detectors can realize boson subtraction operators [40][41][42][43], we can design N -partite genuine entanglement generation schemes with this theoretical process.
In the sculpting protocol, the difficulty of designing a circuit for an N -partite entangled state is translated to the difficulty of finding a suitable sculpting operator to "chisel" the state.However, former research on the sculpting protocol [39,44] provides only proof-of-concept schemes demonstrating the method and lacks any systematic way of linking the features of sculpting operators to the expected final states.Our work shows that a graph picture of the sculpting protocol provides an straightforward strategy of finding appropriate sculpting operators for entanglement, hence appropriate heralding schemes themselves.We map multi-boson systems with sculpting operators into bipartite graphs (bigraphs), for which we develop techniques to understand key properties of the entanglement generation process.Our list of correspondence relations between sculpting protocols and graphs is a variation of that given in Ref. [10], which provided a systematic method to analyze and design LQNs for obtaining entanglement with postselection.In this graph picture, we have found a special type of bigraphs, which we name effective perfect matching (EPM) bigraphs.These bigraphs are highly useful because they can directly correspond to sculpting operators that generate entanglement.states has primarily relied on trial and error methods via the direct route depicted by the blue to red boxes, or by taking a detour through the yellow box.However, we present a more systematic approach to designing circuits for various entangled states by mapping the elements of many-boson systems onto graphs, as illustrated by the routes of the green arrows.
scribes our graph strategy to search for genuinely entangled states.
With our graph-theoretic approach, we present sculpting operators that generate qubit N -partite GHZ and W states, and an N = 3 Type 5 entangled state (the superposition of GHZ and W states [45]).The GHZ and W schemes are significantly more efficient than those given in Ref. [39].And contrary to the schemes in Refs.[25,37], they work for arbitrary number of parties.In addition, our N = 3 Type 5 entangled state generation scheme illustrates that our approach can be extended to find more generalized forms of entangled states.
To top it off, by generalizing the bigraph used to obtain qudit GHZ states, we also present a qudit N -partite GHZ state generation scheme with dN bosons.To our knowledge, our scheme requires much less bosons than any known heralded schemes as in Ref. [46].Our theoretical schemes can be realized in any many-boson system, e.g., linear optical systems with polarization qubit encoding and heralded detections.These outcomes showcase the effectiveness of our method in finding simple solutions for the generation of intricate entangled states.
Our work is organized as follows: Sec.II reviews the sculpting protocol introduced in Ref. [39].Sec.III explains our dictionary of mapping the sculpting protocol to bigraphs.We also show that the perfect matchings (PMs) of bigraphs determine the final state after the sculpting operation.Sec.IV gives sculpting operators that generate qubit Npartite GHZ and W states.Using the qubit GHZ generation graph, Sec.V presents a qudit GHZ state generation scheme.Sec.VI explains how linear optical systems with polarization qubit encoding and heralding detectors can build our sculpting schemes.To showcase our method, we present a simple Bell state generation example.Sec.VII summarizes the significance of our results and discusses possible follow-up researches.
For the video summary of this paper, go to https: //youtu.be/iQ4aWQJuRZI.

II. SCULPTING PROTOCOL FOR QUBIT ENTANGLEMENT
In this section, we formalize the sculpting protocol [39] that converts the boson identity into entanglement.While N -partite entangled state was constructed in Ref. [39] based on 2N modes with the dual-rail qubit encoding, we re-explain it based on N spatial modes and consider the qubit state as a two-dimensional internal degree of freedom of bosons.This way of expression not only embraces the dual-rail encoding, but also provides a more intuitive description of qubit states in the system.
Since in our setup each boson in jth spatial mode (j ∈ {1, 2, • • • , N }) has a two-dimensional internal degree of freedom s (∈ {0, 1}), boson creation and annihilation operators are denoted as â † j,s and âj,s respectively.Then, as an input state, we distribute 2N bosons into N spatial modes so that each spatial mode has two bosons with orthogonal internal states (0 and 1, see Fig. 2).Therefore, the initial state is given by Following the former works [39,44], we call it the maximally symmetric state of 2N bosons.Rewriting |Sym N ⟩ in the mode occupation representation as and the particle number distribution in the 2N states as a vector, we see that |Sym N ⟩ is majorized by all the other Fock states of 2N bosons (i.e., since each mode has one particle, the particle number distribution vector of |Sym N ⟩ is written as ⃗ n Sym = (1, 1, • • • , 1).This vector is majorized by any particle number distribution vector of 2N dimension.See Ref. [47], Sec.II for the rigorous definitions and analyses).Several research papers showed that this kind of state is very resourceful in many quantum computation protocols [47][48][49].
What we need to obtain in the sculpting protocol is a state in which each spatial mode has one boson whose internal state corresponds to the qubit state.For such a final state, we need to annihilate N single-bosons from the initial state |Sym N ⟩.The N single-boson subtraction operators, which we name the sculpting operator, is expressed in the most general form as We see that the one-boson subtraction operator Â(l) can superpose among different spatial modes.Such an operation has been implemented probabilistically in several bosonic experimental setups [40][41][42][43].
For a later convenience, we rewrite Â(l) as where α Applying the sculpting operator ÂN to the initial state |Sym N ⟩, we obtain the final state, The correlation of |Sym N ⟩ from the indistinguishability of bosons and the spatial overlap among bosons that comes from ÂN establish the entanglement property of |Ψ⟩ f in .Therefore, we can state that the indistinguishability of identical particles and spatial overlap are two essential elements for the entanglement generation with sculpting protocol, as in postselected schemes [9,12,24,50].
There is an essential restriction that we need to impose on a sculpting operator ÂN : it must be a sum of operators that annihilates one particle per spatial mode so that the final total state |Ψ⟩ f in must consist of states with one particle per spatial mode.Then the internal state of one particle encodes the qubit information.This also means that any term of ÂN that subtracts both particles from a given spatial mode must vanish.Otherwise, as the total number of subtracted particles is fixed, there remains a term with two particles in the same spatial mode in |Ψ⟩.We call this restriction the no-bunching restriction.
We can find such an ÂN that satisfies the restriction by controlling the probability amplitudes of it.
For the simplest N = 2 example, Â2 is written as and the final state is given by In the above equation, α 1 β 1 â1ψ1 â1ϕ1 |Sym 2 ⟩ must vanish, otherwise the term has two particles in the second spatial mode which violates the no-bunching condition.
In the dual rail encoding setup [39,44], the nobunching condition appears as a seemingly different form.In that setup, repetitive annihilations on the same spatial mode naturally vanish.However, since two spatial modes combine to constitute one subsystem for the case, valid final states are only restricted to those with one boson per two spatial modes.This exactly corresponds to the no-bunching restriction in our setup.
All things considered, we summarize the sculpting protocol as follows: Sculpting protocol 1.Initial state: We prepare the maximally symmetric state |Sym N ⟩ of 2N bosons, i.e., each boson has different states (either spatial or internal) with each other as Eq. ( 36) (see Fig. 2).
2. Operation: We apply the sculpting operator ÂN of the form (4) to the initial state |Sym N ⟩.The sculpting process must satisfy the no-bunching condition.

Final state:
The final state can be fully separable, partially separable, or genuinely entangled.
It is worth mentioning that the degree of entanglement for an N -partite pure state can be categorized into three classes: fully separable, partially separable, and genuinely entangled [51] (see Ref. [22], 3.1 for a quick summary of these concepts).A state |ψ⟩ It is genuinely entangled if it cannot be separable under any bipartition of H.It is partially separable when it is neither fully separable nor genuinely entangled.The states we target in the work are genuinely entangled states such as GHZ and W states.
Another crucial remark is that a sculpting operator can generate an N -partite entangled state with K ancillary spatial modes.Then the above sculpting protocol can be slightly varied to an (N +K)-partite intial state with the sculpting operator ÂN+K .We can see such cases in Sec.IV B and IV C.
Most of the technical difficulty to find ÂN for a specific entanglement state comes from Step 2, for it is critical to control the probability amplitudes so that the sculpting operator satisfies the no-bunching restriction.There have been no systematic technique to find a suitable ÂN that simultaneously satisfies the no-bunching condition and generates nontrivial entanglement state [39].As we will explain in the following sections, our graph techniques facilitate a powerful tool to overcome this limitation.

III. GRAPH PICTURE OF BOSON SYSTEMS WITH SCULPTING OPERATORS
In this section, we present a list of correspondence relations between the fundamental elements of the sculpting protocol and those of graphs.With the mapping, we can replace key physical properties and restrictions on the sculpting operators with those on graphs, which renders a handy guideline to the operator-finding process for genuinely entangled states.
Ref. [22] proposed a list of correspondence relations between linear quantum networks (LQNs) and graphs for providing a systematic method to analyze and design networks for obtaining entanglement without postselection.Since our sculpting protocol also consists of linear transformations of boson subtraction operators, a similar graph mapping dictionary can be imposed to find a suitable ÂN that generates genuine entanglement.Indeed, we can map spatially overlapped subtraction operators into graph elements with a variation of the correspondence relations in Ref. [22], which leads to a practical graph-theoretic method to analyze our system.
The correspondence relations of elements between bosonic systems with sculpting operators and bigraphs can be enumerated as follows:

Boson systems
Bipartite Graph with sculpting operator j and ψ (l) j are defined as in Eq. ( 4).A brief glossary in graph theory can be found in Ref. [22], Appendix A.
In our graph picture, a subtraction operator Â(l) = N j=1 α (l) j âj,ψ (l) j is denoted as an unlabelled vertex in V .Dynamical variables specifying the operator, such as spatial distributions and internal states, are encoded as weighted edges connecting V to labelled vertices in U .Below, unlabelled and labelled vertices are drawn as dots (•) and circles ( j ) respectively.The array of subtraction operators (dots) are on the right hand side of the array of spatial modes (cir-cles).We can consider a more comprehensive mapping including creation operators, which is given in Appendix A. However, Table I suffices to analyze the crucial properties of sculpting operators for generating entanglement.
When Â2 is expanded as each term corresponds to a possible collective path (a possible connection of dots to circles in which each dot is uniquely connected to one circle) of the annihilation operators, e.g., the bigraph (10) has four possibilities for two annihilation operators to be applied to the spatial modes.Therefore, the expansion of Â2 ( 12) is expressed with collective paths as i.e., Â2 is a superposition of the above four collective paths.
For the sculpting operator Â2 to obey the nobunching condition, we must set the amplitudes so that the first two collective paths in (13) vanish when they are applied to We can achieve such a sculpting operator by setting Then it is direct to check that i.e., by fixing amplitudes as Eq. ( 14), we obtain a Bell state as the final state.
From the above N = 2 example, we can understand the role of the no-bunching condition in the graph picture.Since the bigraph expression of ÂN such as (10) is expanded with a summation of all the possible collective paths of annihilation operators as Eq. ( 13), we have to control the complex weights of the edges so that any collective path with more than two edges in the same circle does not contribute to the final state.This property can be understood with the concept of perfect matchings (PMs), which are independent sets of edges in which every vertex of U is connected to exactly one vertex of V (see Ref. [22], Appendix A), as follows: Property 1.For a specific sculpting operator ÂN , the final state |Ψ⟩ f in = ÂN |Sym N ⟩ must be fully determined by the addition of the perfect matchings (PMs) of the bigraph corresponding to ÂN .Indeed, we see that the two collective paths in Eq. ( 15) are the two perfect matchings of the bigraph (10).The above property is useful for understanding given sculpting operators in several aspects, which we explain in Appendix B with a general sculpting-operator-finding strategy based on this property.From now on, a bigraph that corresponds to a sculpting operator is called a sculpting bigraph.

IV. QUBIT ENTANGLEMENT: GHZ AND W STATES
In this section, we present sculpting operators that generate qubit N -partite GHZ state, N -partite W state and a superposition of N = 3 GHZ and W states, using Property 1.Our operator solutions are more efficient and more feasible to construct in many-boson systems than those given in Ref. [39], especially for the W state case.
To find the sculpting operators for those entangled states, we define a specially convenient type of bigraphs, which we dub effective PM bigraphs (EPM).We restrict our attention to sculpting bigraphs whose edge weights of internal states are only among {|0⟩, |1⟩, |+⟩, |−⟩} with |±⟩ ≡ 1 √ 2 (|0⟩ ± |1⟩).Among the creation and annihilation operators in the above basis, we can easily see the following identity holds, which directly results in the following identities: The above identities are translated into our bigraph language as j = 0, Here, the internal state edge weights {|0⟩, |1⟩, |+⟩, |−⟩} are denoted as edge colors {Black, Dotted, Red, Blue} respectively for the convenience.
The amplitude edge weights are omitted.The translation from Eq. ( 17) to ( 18) can be explained more clearly with directed bigraphs (see Appendix A).
Then, we define effective PM bigraphs (EPM bigraphs) as bigraphs whose edges always attach to the circles as one of the above forms.An example of EPM bigraphs is the N = 2 bigraph (10) with restrictions (14), i.e., 1 2 because all the edges are attached to the circles as the first form of (18).From the identities (18), we can see a crucial property of EPM bigraphs as follows: Property 2. If a sculpting bigraph is an EPM bigraph, then the final state is always fully determined by the PMs of the bigraph.
The Combination of Properties 1 and 2 provides a convenient strategy to find sculpting operators that generate a specific entangled state.Since we can express an entangled state with a addtion of PMs, if we can draw an effective PM bigraph which has the same PMs, the bigraph corresponds to a sculpting bigraph that generates the entangled state.We will show that various qubit N -partite genuinely entangled states can be generated with such bigraphs.

A. Qubit GHZ state
The sculpting bigraph that generates the Npartite GHZ state is given by where the edge weights represent the probability amplitudes and edge colors Red and Blue represent the internal states |+⟩ and |−⟩.This bigraph was also used in Ref. [22] to obtain the GHZ state in LQNs (see bigraph (30) of Ref. [22]).
The sculpting operator ÂN corresponding to ( 20) is where ⊕ N in the last line is defined as the addition mod N .It is simple to verify that the bigraph (20) corresponds to a sculpting operator that generates the GHZ state.First, the edges of (20) attach to circles as the first of the three graphs in (18).Therefore we see that only the PMs contribute to the final state.
Second, the bigraph has two PMs and (the above equality holds since the dots are identical), which constructs the GHZ state.
In the operator form, with the identity ( 16), we see that the final state is explicitly given by From the normalization factor, we directly see that the success probability becomes 1/2 N −1 .Note that we can find other sculpting bigraphs for the GHZ state based on (20).While the GHZ state is invariant under the permutation of spatial modes, the bigraph (20) is not.Therefore, any bigraph with the permuted vertex labels of (20) also generates the GHZ state, i.e., σ(1) σ( 2) under a permutation σ (∈ S N ).Since the GHZ state is also invariant under the qubit state flip, the exchange of blue and red edges also gives the GHZ state.However, such graph transformations are already included in the above permutation.
It is worth comparing our GHZ solution (21) with the solution given in Ref. [39], Most importantly, the consecutive annihilations in (26) do not remove particles from orthogonal modes.Hence they are very challenging to realize with experimental setups.In contrast, the procedure ( 21) is based on orthogonal modes, so that a single unitary change of basis is sufficient to prepare all the modes from which a single particle is to be removed.On top of that, each mode in ( 26) is a weighted superposition of all the initial ones.Understanding the operator from the graph picture, (26) corresponds to a bigraph with N N edges.On the other hand, (20) corresponds to a bigraph with only 2 N edges.Therefore, the scheme described by ( 20) is more effective in the sense that each annihilation operator used there is constructed by superposing just two modes with internal state basis changes.

B. Qubit W state
A sculpting bigraph for N -partite W state can be conceived with one ancillary spatial mode as The initial state is prepared in a slightly varied way as one ancillary boson at the ancillary spatial mode.It is as manifest as for the GHZ state case to see the above sculpting bigraph (27) generates W state. First, the bigraph is an effective bigraph since edges attach to circles as the first and second graphs in (18).Second, the above bigraph has N PMs which correspond to the W state with an ancillary system.
In the operator form, with the first identity of Eq. ( 16) again, we have The success probability is |α N −1 β| 2 , whose maximal value becomes (N −1) N −1 We can also check that this bigraph can be used to generate W state in LQNs with postselection.Indeed, by drawing a bigraph that corresponds to the schemes suggested in Ref. [16,17], we obtain the same form of bigraph with (30).
Comparing with the W state generation scheme suggested in Ref. [39], we can easily see that our current scheme have accomplished an outstanding improvement.The scheme in Ref. [39] starts from 4N bosons in 2N modes and goes through two steps of sculpting to generate the final N -partite W state.On the other hand, using the graph mapping technique, we have obtained a much more efficient Npartite W-state generation scheme just with 2N + 1 bosons in N + 1 spatial modes and one simple step of sculpting.

C. N = 3 Type 5 states
The genuinely entangled states that we have discussed so far have some convenient symmetries, which admit relatively simple sculpting bigraphs for generating them.However, we can also conceive less symmetric entangled states with the support of ancillary modes.It is always achieved by any EPM bigraph that connects all the dots to A j with black or dotted edges.
As an example, we present a EPM bigraph for a tripartite system that generates a superposition of the N = 3 GHZ and W states, which is called the N = 3 Type 5 state in Ref. [52].The EPM bigraph used here includes three ancillae to construct such a sculpting operator: In the above bigraph, the amplitude weights are omitted under the assumption that they are nonzero and satisfy the normalization conditions.The corresponding sculpting operator of the above bigraph is given by Â3+A,B,C When the above sculpting operator is applied to the initial state See that the five states with nonzero amplitudes in the final line of the above equation constitute the set of bases that can transform to any tripartite state under local operations [45].The final state we just obtained is categorized as N = 3 Type 5 state in Ref. [53].This N = 3 example shows that the graph method has potential to design other general form of multipartite entangled states with heralding detectors.
V. QUDIT ENTANGLEMENT: GHZ STATE Our graph picture also provides a useful insight to find sculpting operators for the general qudit systems.We will present in this section a sculpting bigraph for the qudit N -partite GHZ state, which has a generalized form of the qubit GHZ bigraph (39).
The qudit state is represented by a d-dimensional internal degree of freedom s (∈ {0, 1, • • • , d}) of bosons.To construct N partite qudit genuinely entangled states, we initially distribute dN bosons into N spatial modes so that exactly d bosons with mutually orthogonal internal states belong to a spatial mode (see Fig. 3).Hence, the initial state is given by Here |Sym N,d ⟩ denotes the N spatial mode maximally symmetric state with a d-dimensional internal degree of freedom.The sculpting operator ÂN = must be set to extract (d−1) bosons per spatial mode so that one boson per spatial mode in the final state determines the qudit state of each subsystem.All in all, the sculpting protocol is modified for qudits as follows: Sculpting protocol of qudits  36).See Fig. 3.

Operation:
We apply the sculpting operator ÂN to the initial state |Sym N,d ⟩.
The sculpting operator must be set to extract (d − 1) bosons per spatial mode.

Final state:
The final state can be fully separable, partially separable, or genuinely entangled.
Now we provide a sculpting operator that generates the N -partite GHZ state of d-level systems, denoted as |GHZ N,d ⟩, by generalizing the qubit GHZ sculpting operator in Sec.IV A.
First, by generalizing the d = 2 basis set {|+⟩, |−⟩} for the internal states of the sculpting operators, we choose the arbitrary d-dimensional basis set {| k⟩} d−1 k=0 where (ω = e i 2πk d ) for the internal states of the sculpting operators.
Second, we use an overlap of (d − 1) copies of the graph (20) for the sculpting bigraph, i.e., the following bigraph corresponds to the sculpting operator Then, the final state is given by The second line is obtained by Eq. ( 44) and the third by Eq. ( 43).The success probability is 1/3 N −1 .
Remark.-To the best of our knowledge, our scheme needs much less bosons than any other heralded schemes for the qudit GHZ state.For example, the scheme in Ref. [46] that needs (2d + 1) particles for the d-level bipartite GHZ state (i.e., the Bell state) and 25 particles for the 3-level tripartite GHZ state in their optimized method.Our sculpting bigraphs just need 2d and 9 particles respectively.

VI. HERALDED SCHEME OF SCULPTING PROTOCOLS IN LINEAR OPTICS: BELL STATE EXAMPLE
In Sec.IV and V, we have presented sculpting bigraphs that generate genuinely multipartite entangled states.Therefore, if we know how to build spatially overlapped subtraction operators by heralding, we can directly design heralded entanglement generation circuits by combining these operators as the structure of sculpting bigraphs.
There are general schemes in optics [40][41][42][43] to establish subtraction operators of bosons.Based on such methods, Ref. [39] proposed an optical scheme for constructing sculpting operators.More recently, Ref. [44] suggested an experimental scheme with arithmetic subtractions of trapped ions (which are near-deterministic operations established in Ref. [54] that work unitarily except when the system is in the vacuum) to generate the GHZ state with the sculpting operator (26) in Ref. [39].Ref. [39] and [44] used dual-rail encoding and binomial encoding respectively.
Here, we briefly explain an alternative linear optical circuit designed to implement spatially overlapped subtraction operators by heralding, which becomes building blocks of heralded schemes that generates entanglement.Then, as a proof of concept, we design a heralded Bell state generation scheme by utilizing the heralded subtraction operator and the Bell state bigraph (19).The optical elements used in the circuit can be applied to more general multipartite entagled states.For a more comprehensive explanation and solutions on the construction of heralded schemes from our sculpting schemes, see Ref. [55].
Since all the qubit sculpting schemes in our work correspond to EPM bigraphs whose final states are determined by the identities ( 16) and ( 17), we need to find heralded optical circuits that corresponds to those identities.The basic elements of our optical schemes are polarizing beam splitters (PBSs) and half-wave plates (HWPs), hence highly feasible.It also implies that our translation rule can be applied to any bosonic systems with operators that play the roles of PBSs and HWPs.
In our setup, we encode the internal boson states {|0⟩, |1⟩, |+⟩, |−⟩} as the polarization of photons {|D⟩, |A⟩, |H⟩, |V ⟩} where (D=diagonal, A=antidiagonal, H= horizotal, V =vertical).The polarized states have the following relations: Then the initial state is given by We first consider the heralded optical circuit for the identity (16), which can be rewritten as â± We can perform the above operation with the following optical circuit: (51) where the one-photon detection at 21 and 22 correspond to â+ and â− respectively.In the above circuit, PBSs transform photons as , and HWPs as {H, V } ↔ {D, A}.
Step-by-step explanation 1.The first HWP rotates the photon state basis: 2. The first PBS divides the photon paths according to the internal states: where 1 and 2 denote the upper and lower paths of the PBS.

The dashed purple box subtracts â †
by heralding in the last PBS and send the remained one to the initial mode: In the last PBS, the detection of a photon with V in the upper mode 21 (H in the lower mode 22) heralds the final state â † H (â † V ).
In the above process, the dashed purple box plays the role of the heralded subtraction operator.We can deform the subtractor to design a heralded optical circuit for spatially overlapped subtraction operators such as (â 1,+ − â2,− ).Instead of attaching both wires of the subtractor to the same mode as in (51), we now attach two wires to different modes 1 and 2 as in (55) so that it plays the role of a spatially subtraction operator.A crucial difference of the subtraction operator in the above from that in (51) is the need of an HWP between two PBSs, which makes possible the subtraction of different internal states from difference input modes.
As a proof of concept that the above heralded operator plays the role of spatially overlapped subtraction operator, we provide a Bell-state generation scheme by employing the heralded subtraction operator based on the sculpting bigraph (19), which corresponds to in our optical setup.We implement two heralded subtraction operator of the form (55) following the structure of the sculpting bigraph ( 19) that corresponds to (57), which results in the following circuit: (58) By postselecting only the cases when each detector observes one photon, we generate Bell states as expected.Note that the wires of heralded subtraction operators are attached as the structure of the Bell sculpting bigraph (19) (see the wires in the dashed purple box).Appendix D explains the state evolution in detail.
In the case of W state and N = 3 Type 5 state, we require subtraction operators for the superposition of larger than 2 spatial modes in the ancillae.This can be achieved by generalizing (55).Ref. [55], Appendix A provides a more thorough analysis and optical circuits specific to these states.It is worth noting that any bosonic system with linear operators that transforms the spatial and internal states as PBSs and HWPs can execute the same entanglement generation schemes as those given here.
For the case of qudit entanglement generation with sculpting operators, we need a bosonic system that has a higher level of internal degree of freedom.For example, photons can have the orbital angular momentum (OAM) that can encode qudit information.Therefore, we can generate the GHZ qudit entanglement given in Sec.V with OAM beam splitters [56] and OAM-only Fourier transformation operators [57].OAM beam splitters change the outgoing paths of photons with respect to the internal states, hence a d-level generalization of the PBS.OAM-only Fourier transformation operators transform the computational basis of the internal states, hence a d-level generalization of the HWP.

VII. DISCUSSIONS
Our strategy for finding entanglement generation schemes based on linear bosonic systems with heralding involves a two-step process: first, we find a theoretical sculpting operator that generates an entangled state.Second, we construct a concrete experimental circuit for such a sculpting operator.For the first step of finding sculpting operators, we have exploited graph techniques by imposing the correspondence relations of bosonic systems to bigraphs.We have shown that the graph picture of bosonic systems facilitates a powerful tool to find proper sculpting operators.For the second step, we have explained that a spatially overlapped subtraction operator can be installed in linear optical networks with heralding.This operator allows us to design circuits for generating heralded entanglement by combining them based on sculpting bigraphs.As the simplest example, we have presented a Bell state circuit with the heralded subtraction operators, which can be extended to other sculpting schemes that we have proposed (see Ref. [55]).Our current results suggest several interesting future research directions.
First, our formalism and strategy can be extended to encompass more complex qudit systems.We have introduced EPM bigraphs for finding qubit solutions, which can be generalized to qudit cases.We can suggest a more complete demonstration on the construction of desigining qudit heralded schemes with more solutions for qudit entanglement.We can encode such qudit entangled states as the orbital angular momentum (OAM) of photons with OAM beam splitters [56] and OAM-only Fourier transformation operators [57].
Second, our sculpting protocol can identify other interesting multipartite entangled states.Since our graph approach provides a handy guideline to coming up with useful sculpting operators, we expect that it would be used to find heralded schemes for other crucial entanglements.For example, Ref [58] has recently found sculpting schemes for generating a special type of graph states, i.e., caterpillar graph states, with less photons than fusion gates [59].Therefore, one of our next goals will be to find, e.g., more general types of graph states, N -particle N -level singlet states [60], and cyclically symmetric states [61].
FIG.1: Previous research on generating entangled states has primarily relied on trial and error methods via the direct route depicted by the blue to red boxes, or by taking a detour through the yellow box.However, we present a more systematic approach to designing circuits for various entangled states by mapping the elements of many-boson systems onto graphs, as illustrated by the routes of the green arrows.

FIG. 2 :
FIG. 2: The initial state |Sym N ⟩ of 2N bosons in N spatial modes.Each spatial mode has two bosons, one in the internal state |0⟩ and the other in |1⟩.

TABLE I :
Correspondence relations of a sculpting operator to a sculpting bigraphIn the above table, α 1. Initial state: We prepare the maximally symmetric state |Sym N,d ⟩ of dN bosons, i.e., each boson has different states (either spatial or internal) with each other as Eq. ( AtStep 5, we sort out the states that click one of the two detectors in each PBS.When the two heralding detectors receive photons of the same polarization, we obtain a Bell state 1 √ 2 (|HH⟩ + |V V ⟩).When the two heralding detectors receive photons of the opposite polarization, we obtain a different Bell state1