Generating nonclassical states of motion using spontaneous emission

Nonclassical motional states of matter are of interest both from a fundamental perspective but also for their potential technological applications as resources in various quantum processing tasks such as quantum teleportation, sensing, communication, and computation. In this work we explore the motional effects of a harmonically trapped, excited two-level emitter coupled to a one-dimensional (1D) photonic system. As the emitter decays it experiences a momentum recoil that entangles its motion with the emitted photon pulse. In the long-time limit the emitter relaxes to its electronic ground state, while its reduced motional state remains entangled with the outgoing photon. We find photonic systems where the long-time reduced motional state of the emitter, though mixed, is highly nonclassical and in some cases approaches a pure motional Fock state. Motional recoil engineering can be simpler to experimentally implement than complex measurement and feedback based methods to engineer novel quantum mechanical states of motion.


Introduction
Engineering nonclassical quantum states of motion has long attracted much attention particularly if the motional properties of the object in question admit the potential to generate macroscopic superpositions. Early work, both theoretical and experimental, on the creation of nonclassical motional states focused on the motional states of individual trapped ions in the late 90s [1,2,3,4,5,6,7]. Much of this work was performed in the Lamb-Dicke (LD) regime where the photon recoil energy is much smaller than a motional quanta, and thus the ion suffers little motional excitation. Beyond this regime (so called large LD regime), schemes for cooling via electromagnetically induced transparency [8] and the existence of bistable motional states [9] have been proposed. Studies involving trapped ions are typically done in a free-space setting. Alternatively, strong coupling of quantum systems to 1D electromagnetic environments, such as optical waveguides, provides a new avenue to study light-matter interactions. Theoretical studies have shed light on the scattering of single-photon pulses [10], waveguide-mediated entanglement [11], ultra-strong coupling [12], and spatial light rectification [13,14,15,16]. However, little theoretical research has been done on the optomechanical effects of mobile emitters ina 1D photonic system, even while much experimental progress has been made towards strong coupling of optically trapped emitters [17,18,19]. A notable exception is the recent work of Li et al. [20], which studied the modification of single-photon reflection spectra due to motional recoil of a trapped two-level emitter near a 1D waveguide. Complementary studies [21,22,23] use a standard phenomenological optomechanical coupling to include motional effects, including the creation of nonclassical motional states [24,25]. In the following we focus on microscopic optomechanical systems where an emitter, e.g. an atom or nanoparticle containing a emitter, is held in a harmonic motional trap near to a 1D photonic system. The motional coupling is fundamentally mediated by the recoil suffered by the emitter upon emission/absorption of a photon, and no phenomenology need be assumed for the optomechanical interaction.
In this article we explore the motional effects of spontaneous emission for an excited, harmonically trapped quantum emitter both within and beyond the LD regime. Optical decay into the modes of a bidirectional waveguide couples the emitter's motional state to the outgoing photon. In the long-time limit all the light escapes away along the waveguide, and the emitter relaxes to its electronic ground state. However, during the emission process the motional state is disturbed indefinitely. In the long-time limit entanglement between the motion and the outgoing photonic field produces reduced motional states that are often highly mixed and lack nonclassical features. We examine this pedagogical case in Sec. 2. We then explore two alternate physical architectures where the reduced motional state can become highly nonclassical. In Sec. 3 we consider an emitter next to a perfect mirror embedded in the waveguide. Including the motional coupling significantly modifies the well-known inhibition of decay due to interference from the field reflected off the mirror. We find that this effect can be used to generate pure single-excitation Fock states of motion. The second physical architecture involves coupling the harmonically trapped emitter to the circulating modes of a toroidal optical cavity, which is symmetrically coupled to left-and right-propagating modes in the waveguide. This architecture has been experimentally realized with tapered optical nanofibers by Rauschenbeutel and co-workers who have demonstrated a variety of effects including strong coupling [26], optical switching [27], and optical isolation [28]. To date the effects of recoil on the trapped emitter have not been experimentally investigated. We explore this architecture in Sec. 4, and identify a parameter regime in which the reduced motional state is both highly entangled with the emitted photonic field and also individually nonclassical.
For each physical scenario a quantum system couples to a bidirectional waveguide that supports a continuum of modes for each propagation direction. These modes behave as an electromagnetic reservoir into which system excitations can decay. This open-systems description of light-matter interactions is conveniently treated within the context of input-output theory [29]. We consider the case where there are no initial optical fields present and there are no intrinsic losses, i.e. the combined motional/optical system evolves unitarily. Then, tracing over the two waveguide channels gives the dynamics of the reduced quantum system in terms of a Markov master equation, where the subscripts { , r} refer to the left-and right-propagating waveguide modes, respectively. The time evolution proceeds according to a system HamiltonianĤ sys and two jump operatorsL i , one for each propagating mode. The jump operators describe the coupling of the quantum system to the left-and right-propagating modes. For each of the physical scenarios below we solve the master equation numerically for the system stateρ(t) and then trace out additional system degrees of freedom to find the reduced state of the motion,ρ m (t) = Tr H 0 [ρ(t)], where H 0 is the Hilbert space of all other system degrees of freedom.

Trapped atom directly coupled to a waveguide
We first consider a trapped quantum emitter strongly coupled to a bidirectional waveguide as shown in Fig. 1a. While this scenario reveals no remarkable nonclassical features in the unconditional long-time motional state of the emitter, it serves as an instructive launching point for the investigations below. Further, by considering the conditional state of motion given detection of the emitted photon, as shown in Fig. 1b, we gain insight into physical settings that can indeed produce nonclassical signatures in the reduced motional state. Based on this we later study the case when the emitter is placed in front of a mirror as depicted in Fig. 1c. The emitter is taken to be a two-level system with internal electronic eigenstates |g and |e separated by transition frequency ω 0 . The emitter is held in a harmonic trap with a zero-point motion amplitude x zpm , and bosonic operatorsv, andv † , that respectively annihilate and create motional quanta at frequency ω m and satisfy [v,v † ] = 1. We work in the interaction picture with respect to the electronic degree of freedom, so the Hamiltonian is that of harmonic motion in the trap,Ĥ sys = ω mv †v . The emitter decays symmetrically into the left-and right-propagating waveguide modes with total optical decay rate Γ. Conservation of momentum requires that the emitter experience a momentum recoil opposite to the direction of photon emission. The jump operators that describe the process of decay and momentum recoil within the context of the master equation, Eq. (1), areL whereσ − = |g e| is the electronic lower operator. The momentum recoil that accompanies optical emission is generated by the dimensionless position operator x =v +v † with a strength determined by the Lamb-Dicke parameter η ≡ k 0 x zpm , where k 0 = 2π/λ 0 is the magnitude of the emitted photon's wavevector. An emitter initially prepared in the excited electronic state and motional ground state, |ψ(t 0 ) = |e ⊗ |0 m , will decay into the waveguide modes. For times long compared to the optical decay rate, t 1/Γ, the initial excitation decays entirely into the waveguide. The emitter's long-time state after the light has been traced out, ρ(t) = |g g| ⊗ρ m (t), is given by the master equation. Typically, after the emitted light propagates away the reduced motional stateρ m (t) is mixed due to entanglement between the outgoing photon and the motional degree of freedom. Nevertheless, mixed motional states can in principle exhibit nonclassicality; we seek such a physical situation.
Decay into the left-and right-propagating waveguide modes results in both positive and negative displacements in momentum. These displacements are correlated with the direction of the outgoing photon, resulting in a mixture of positive and negative momentum recoils when the photonic degrees of freedom are traced out. For an initially motionless trapped emitter we find two qualitatively different regimes. If the coupling to the waveguide is large compared with the motional frequency Γ/ω m > 1, then the resulting displacement is an impulsive kick, which leads to long time states whose motional Wigner function exhibits two cleanly separated peaks, see Fig 2(a). As the coupling to the waveguide decreases, the decay and resulting motional displacement occur over a duration which is either commensurate with, or longer than, the motional period. The resulting motional Wigner function becomes smeared, see Figs 2(b) and 2(c). However, no parameter values (η, Γ/ω m ), were found that result in a long-time motional Wigner function with any negativity. That is, the atomic motion is left in a highly mixed "classical" motional state. The reduced motional state purity, Tr[ρ 2 m (t)], is shown in 3. To further investigate the level of correlation between the motion and emitted optical fields we compute the motional entropy S m ≡ −Tr[ρ m lnρ m ], as a function of (η, Γ/ω m ). Since the joint state of the motion-excitation-optical field is pure at all times, and in the long-time limit the emitter is in the electronic ground state, the motional entropy is an entanglement measure for t → ∞. Large motional entropy indicates a large correlation between the motion and the outgoing photon field. From Fig. 3, we observe that the motion is heavily entangled with the outgoing photon pulses, preventing the reduced motional state from displaying negativity in its Wigner function.
Information about the direction of momentum recoil contained in the outgoing photon field can be retrieved by measuring it. A "click" at a photo-detector placed at each end of the waveguide determines whether the photon was emitted to the left or right. The corresponding recoil is simply a single coherent displacement of the momentum to the right or the left, respectively. In order to see nonclassical behavior, the "whichway" information in the outgoing photon field must be coherently erased. This can be achieved by placing a beamsplitter before the detectors as depicted in Fig. 1b exiting either port of the beamsplitter is no longer correlated with a single momentum kick but rather a superposition of recoils to the right and left. A click at either detector thus heralds a motional state in a superposition of momentum recoils. Additionally, detection of the photon disentangles it from emitter, and resulting conditional motional state is pure.
To investigate the conditional motional state, we simulated the evolution using a stochastic master equation for photon counting [29], which gives the emitter's state while both waveguide modes are continuously monitored. As long as only vacuum is measured at both detectors, the joint electronic-motional state remains in the initial state, |e ⊗ |0 m , according to the quantum Zeno effect. When the photon is detected, the electronic state is projected into |g , and the motional state receives a superposition of momentum kicks with a relative phase that depends on which of the two detectors clicks. Examples of the resulting conditional Wigner functions are shown in Fig. 2(d)-(f), revealing near perfect motional cat states.
It is well known that one can harness the nonunitarity of quantum measurements to generate nonclassical states. However, in practice the quality of the measurement   has a large effect on the nonclassicality. Meanwhile measurements at the quantum limit have become common but remain difficult to engineer. Thus, informed by the simple conditional model above where eliminating "which-way" information is the key to nonclassicality, we now focus our attention on more elaborate settings where nonclassicality can arise even for unconditional dynamics.

Emitter in front of a mirror
As we have seen above, nonclassical motional states can arise when "which-way" information in the emitted photon is erased coherently. We first consider a straightforward way to achieve this by placing a perfect mirror near the trapped emitter to ensure that all emitted photons eventually exit in one direction down the waveguide, without relying on a conditional detection scheme. In the Lamb-Dicke regime, η 1, motional coherence after spontaneous emission near a mirror has been observed [30].
To simulate the master equation for this situation, we use the SLH-formalism for quantum networks [31,32], to find the appropriate Hamiltonian and jump operators. This formalism was designed to study quantum networks in the context of input-output theory and gives methods for reducing a complex set of interconnected network elements to a simpler description. Following the procedure detailed in Appendix A we obtain the which are used in the master equation, Eq. (1). Here φ = 2πd/λ 0 is the phase acquired by the field as it propagates the distance d between the emitter and mirror. The quantized position operator of the emitter,x =v +v † , appears in an oscillatory term in both the jump operator, Eq. (4), and the effective Hamiltonian. The latter describes an energy shift due to interaction with the reflected field, Eq. (6). We again consider the decay of an initially excited emitter in the ground motional state, |ψ(t 0 ) = |e ⊗ |0 m . The time evolution of the unconditional motional state of the emitter is found from the master equation in Eq. (1) with the jump operators and Hamiltonian in Eqs. (4) and (6). We plot the Wigner functions for the long-time motional state in Fig. 4, for the emitter located at three positions relative to the mirror given by φ = {0, π/4, π/2}. We observe a drastic difference from the previous case of emission into the bidirectional waveguide in Sec. 2. In all situations we observe that the motional state, though highly mixed, is also highly nonclassical due to significant negativity in the Wigner function. We thus have confirmed that if one arranges an optical setup where one forces the emitted light to re-interact with the emitter then the motional state can become nonclassical.
In the absence of motional recoil (η = 0), the optical decay can be modulated by placing the emitter at different positions with respect to the mirror. When the emitter is placed at a node, φ = (2n + 1)π/2 for integer n, the reflected and emitted fields destructively interfere and optical decay is entirely suppressed. This effect has been studied experimentally in free space by Blatt et al. [33,34,35,36], who examined the position-dependent modulation of the fluorescence of a single ion positioned in front of a mirror.
For nonzero Lamb-Dicke parameter, η > 0, the electronic and motional degrees of freedom are coupled, evident in the jump operator Eq. (4). The extent of the emitter's spatial wave function beyond the node, even when in the motional ground state, breaks the perfect destructive interference and allows decay in to the waveguide. This effect can be used to engineer pure, single-excitation Fock states of motion |1 m . The emitter is placed at a node, φ = π/2, and the Lamb-Dicke parameter η is varied, see Fig. (5). As η → 0, the long-time motional state becomes a pure state, as seen in Fig. 5(d). To verify that the pure state is indeed a single-excitation Fock state of motion, |1 m , in Fig.  5(e) we plot the trace distance, 1 2 Tr[ (ρ m −σ) 2 ], whereρ m is the long-time motional state andσ = |1 1| m . Again, as η → 0, the state asymptotically approaches |1 m .

Emitter coupled to a waveguide via a toroidal cavity
We now consider a physical architecture that has been experimentally realized with tapered optical nanofibers by Rauschenbeutel and co-workers [28,37]. A trapped atom is strongly coupled to circulating optical modes in a "bottle resonator" fabricated in an optical fiber, which is evanescently in/out coupled to a tapered nanofiber. We model this architecture as depicted in Fig. 6. The trapped emitter couples with strength g to clockwise(counterclockwise) circulating photons in a toroidal resonator with frequency ω R and annihilation operatorsâ 1 (â 2 ). The emitter is detuned from the resonance frequency of the toroidal resonator modes by ∆ ≡ ω R − ω 0 , and the respective resonator modes are coupled to the left-and right-propagating waveguide modes with rate κ. Rather than breaking the directional symmetry as in Section 3, here the leftand right-propagating waveguide modes are treated on equal footing. The toroidal resonator facilitates a photon reinteracting with the emitter, which can then re-emit into either resonator mode, effectively scrambling the 'which way' information. In the interaction picture the Hamiltonian governing the coupling between the emitter and toroidal resonator modes iŝ a 1 a 2 g κ κ Figure 6. Schematic setup for a harmonically trapped two-level emitter coupled symmetrically with strength g to degenerate clockwise (â 1 ) and anticlockwise (â 2 ) circulating optical modes of a toroidal resonator. The resonator modes are outcoupled with strength κ to the left-and right-propagating modes of a 1D waveguide. The emitter is detuned by ∆ from the toroidal cavity modes. and the jump operators describing coupling from the toroidal resonator to the waveguide modes are,L We initialize the system with the emitter in the excited electronic state and ground motional state, and the toroidal cavities in vacuum. The open-system dynamics of the joint emitter-resonator system are given by the master equation in Eq (1) with the Hamiltonian and jump operators above. In the long-time limit we examine the reduced motional state,ρ m , by tracing out the electronic and resonator degrees of freedom after the emitter has decayed to the electronic ground state and all the light has exited into the waveguide. Example Wigner functions forρ m are show in Fig. (7) for fixed Lamb-Dicke parameter η = 2 and cavity decay κ = 2. Comparing the columns reveals that decreasing the emitter-resonantor coupling with respect to the motional frequency, g/ω m , eliminates the phase information. For small g/ω m the motional state performs many oscillations during a Jaynes-Cummings-type transfer of excitation from emitter to resonator.
Comparing the rightmost rows of Fig. 7 (g/ω m = 0.25) shows that increasing the emitter-resonator detuning produces rings of negativity in the Wigner function similar to those for Fock states of increasing n. However, the motional states cannot be a Fock state since they are highly mixed. Indeed, the motional states are in fact almost perfectly diagonal in the Fock basis, indicating a mixture over perfect Fock states. The number distributions for the diagonal elements ofρ m are shown in Fig. (8). As the detuning is increased, the relative amount of vacuum, n = 0, diminishes, and the resulting motional states are mixtures of nonclassical Fock states.  To compare nonclassicality of the long-time motional states, we compute the integrated negativity of the motional Wigner function, where A − are the areas in phase space in which the motional Wigner function is negative. We seek parameter values that leave W − as negative as possible indicating a highly nonclassical motional state. In Fig 9(a), we observe that nonclassical states require a small coupling between the emitter and the cavities (so that the light circulating in the toroidal cavities can interact with the emitter for longer), and with a periodic dependence on the detuning. Finally, nonclassicality of the long-time motional state has some robustness to situations where the trapped emitter is not initially in the motional ground state. This is common in realistic physical settings where perfect ground-state cooling is not achieved. As shown in Fig. 10, integrated negativity persists when the emitter's motional degree of freedom is initially in a thermal state with mean occupation numbern ≡ Tr[ρ m (0)v †v ].

Conclusion
We have investigated the production of nonclassical motional states in several physical architectures based on coupling a trapped two-level emitter to a 1D photonic waveguide. In each case the fundamental optomechanical coupling between the internal degree of freedom and the motional degree of freedom is provided by momentum recoil as a photon is emitted, as required by conservation of momentum. While nonclassical states can be readily produced by detecting the photon, we have shown that even unconditional motional states can exhibit significant nonclassical features in the longtime limit. Our method for the production of nonclassical motional states relies only on the fundamental coupling between an emitter and its electromagnetic environment making it relatively straightforward compared to procedures involving complicated measurement and feedback. An interesting extension to generate larger, potentially entangled motional states involves multiple emitters in a single trap next to a waveguide. Cooperative emission might yield a superradiant recoil [38] that can be harnessed to produce extremely large motional nonclassical motional states. to the emitter and is absorbed and reemitted in both directions. As the fields propagatesbetween the beamsplitters and the emitter they acquire propagation phases φ 1 and φ 2 , respectively. This FP setup can be modelled as an input-output network with interconnected linkages from each of the three components (two BS and the emitter). The SLH formalism, reviewed in Ref. [39], was designed to model such networks. Within the formalism each component is fully described by three objects: a unitary scattering matrix S, a vector of jump operators L, and a HamiltonianĤ. Collected together as G = (S, L,Ĥ), these are referred to as an SLH-triple. We leave operator hats off of S and L for notational convenience. When multiple components are interconnected, the formalism provides rules to collapse the internal connections such that the entire network can itself be described by a single SLH-triple [31,32,39,40]. Specifically, Rule 4 of Ref. [39] describes how to perform feedback reduction of an SLH-triple by eliminating an internal interconnect. Once all the interconnects have been eliminated, the unconditional dynamics of the network as a whole are given by a master equation, Eq. (1), using the resulting L andĤ. Below we will discuss a particular caveat associated with feedback reduction of optical quantum networks; at the moment we focus on a short review of the procedure within the context of our problem.
Before applying feedback reductions we recast the setup in Fig A1a, into the more abstract form depicted in Fig A1b. We label the input ports from n in = 1, 2, . . . 6, and output ports from n out = 1, 2, . . . 6. We will have to perform four reductions as the network at large has two physical inputs and two outputs (arrows in Fig A1b). Before reductions, the SLH-triple for this six-port input-output device operates on input operator vector x = (a R in ,ā L in , a L in ,ā R in , emitter L , emitter R ), and gives the output vector y = (a R out ,ā L out , a L out ,ā R out , emitter L , emitter R ), where emitter R/L denotes the emitter coupling to the right-and left-moving modes. Using this ordering we denote the scattering matrix for the setup as where S L{R} BS (α{β}) is the two-port scattering matrix for the left{right} beamsplitter. The full S matrix is fully represented as where Cβ ≡ cos β, Sβ ≡ sin β, etc. The last two diagonal entries correspond to the phase acquired by the left-and right-moving modes as they propagate once between  Figure A1. Optical network for a two-level emitter coupled to a 1D waveguide with an embedded Fabry-Perot cavity, modeled as two interconnected beamsplitters. a) Graphical schematic of the emitter in between two beamsplitters. The left(right) beamsplitter has reflectivity β(α). The straight horizontal optical path depicts the right-propagating waveguide mode coming in from the left while the curved optical path depicts the left-propagating mode coming in from the right. Solid arrows indicate external inputs and output, while outlined arrows indicate internal connections. b) Abstract network diagram of the same setup where we label the six input and output ports of the beamsplitters and emitter. Four internal interconnects link the internal inputs and outputs that are eliminated in the feedback reduction.
the beamsplitters. The vector of jump operatorsL describes the linear coupling of the components to the input/output fields. As the beamsplitters have no internal dynamics, the only nontrivialL is that describing the emission of photons. Within the master equation, it drives spontaneous emission of the emitter (and recoil), into the left-and right-moving modes (taking into account the propagation phases). The vector of jump operators is then Finally, we work in the interaction picture with respect to the electronic degrees of freedom, so the HamiltonianĤ is simply that of the harmonic motion in the trap H sys = ω mv †v . We must now interconnect the various internal network links in the optical network in Fig A1b. Some details how this is achieved is outlined in Appendix A.2. Once all the internal links are reduced one obtains a two-input, two-output SLH-triple. Here we focus on two specific cases of interest, the first being when the left BS is completely transparent, e.g. there is no left BS, β = 0, and the right BS is fully reflective, i.e. is a mirror, α = π/2. The second case attempts to model the emitter in a Fabry-Perot cavity with 50:50 mirrors. To model this we set each BS to be a 50:50 BS, e.g. β = α = π/4. First, however, we derive the general SLH-triple for arbitrary α, and β, but set the emitter at the mid-point between the BS e.g. (φ 1 = φ 2 ≡ φ), for simplicity. After all feedback and internal link reductions have been performed, we arrive at the generalized final two-input, two-output SLH-triple: H sys = ω mv †v + |e e| Γe 4iφ 2F 1 F 2 sin(β)(1 + sin 2 (α)) sin(2φ + 2kx) (A.6) + sin(α)(1 + sin 2 (β)) sin(2φ − 2kx) + 2 sin(β) sin(α) sin(4φ) where F 1 = 1 − e 4iφ sin(α) sin(β), and F 2 = e 4iφ − sin(α) sin(β). Note that the scattering matrix is unitary: S † S = SS † = I. By setting β = 0 and α = π/2, we immediately obtain the SLH-triple presented in Sec. 3, for the emitter in front of a mirror. In this case φ represents the phase accrued by the light traveling between the emitter and the mirror.  5   Table A1. Reduction of enumerating the input and output ports as one eliminates internal links in the network in Fig A1b). To derive this table we proceed as follow: the left most columns enumerate the original input(left) and output(right) ports. We first eliminate the internal link connecting 1 out → 6 in . We then relabel the nodes in the second to left most columns of the remaining ports which now range from 1 : 5. We next eliminate the returning link, which in the new enumeration links 5 out → 3 in . We proceed to eliminate all four internal links to arrive at the reducing link identifications: 1 out → 6 in : 5 out → 3 in : 3 out → 4 in : 3 out → 2 in . After these four reductions one is left with a 2 input, 2 output port device.
care to re-identify the enumeration of the ports. As an example we can track the enumeration of the reducing network as we eliminate (in order), the internal links 1 out → 6 in : 6 out → 3 in : 4 out → 5 in : 5 out → 2 in , where we have used the labelling in Fig A1b). The actual enumeration of input to output ports must collapse as the size of the device shrinks and we give such a reducing enumeration for the example elimination we described in the previous sentence in Table Appendix A.2.
We now indicate how to perform the first reduction, obtaining the reduced network once the internal link 1 out → 6 in , is eliminated. We refer the reader to the SLH composition rules in section 5.2 of [39], and in particular Rule 4 or Eqns (61) of [39]. This rules describes how to obtain the reduced G red = (S red , L red ,Ĥ red ) description from the original network G = (S, L,Ĥ), when the internal link connecting output port x to input port y, i.e. x → y, is eliminated. In our example, G describes a 6 × 6, input-output device while G red describes a 5 × 5 device, and we have where Sxȳ is the scattering matrix S, omitting the x th row and y th column. Sx y is the y th column of S with the x th row deleted, S xȳ is the x th row of S with the y th column deleted, and S xy is the (x, y) element of S. Also, Lx is L without the x th row, and L x is the x th row of L.
We begin with the full description of G = (S, L,Ĥ) for the 6 × 6 network: andĤ red =Ĥ sys , since L 1 = 0, and S 16 = 0. To perform the next reduction we must refer back to the Fig A1b). In the original enumeration of the 6 × 6 device we next eliminate the internal link 6 → 3, but in the new labeling of the 5 × 5 reduced device we refer to Table Appendix A.2, and using the second to left most column of the Out and In ports we next eliminate the internal 5 → 3 link. We follow the above procedure to finally arrive at a 2 × 2 SLH network.