Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the general goal of quantum information processing is to manipulate the information coded in a particular quantum system , while simultaneously trying to isolate it from undesired perturbations due to the external environment . for instance , in the fields of quantum computation  @xcite and quantum simulation  @xcite , the processing stage usually relies on a designed unitary evolution which , preserving the coherence of the system , provides a quantum state that contains the outcome of the computation , or the properties of the quantum phase of matter under study . therefore , the experimental progress in these fields has been typically associated to technological developments that minimize experimental imperfections and maximize environmental isolation  @xcite . however , the coupling of a system to its surrounding environment is not necessarily a disadvantage . in fact , the dissipation of a quantum system , if judiciously exploited , can act as a resource for quantum information processing . this interesting change of paradigm started with the recognition that dissipation can assist the generation of entanglement between distant atoms , either in free space  @xcite , or trapped inside cavities  @xcite . by measuring the presence of a photon spontaneously emitted by the atoms  @xcite , or its absence  @xcite , it is possible to project an initially uncorrelated atomic state into a maximally entangled one , as has been recently experimentally demonstrated  @xcite . these schemes can be improved further by tailoring the atom - cavity interaction  @xcite or by modifying the measurement techniques  @xcite , and can also be generalized to provide a route towards universal quantum computation  @xcite . since the dissipation does not always render the desired result in these approaches , the measurement is required to select only the successful outcomes . therefore , the combination of dissipation and measurement can be considered as a probabilistic resource for quantum information . another promising approach is to assess whether , by designing the system - environment coupling in a particular way , the dissipation could provide the desired quantum state with certainty . this _ quantum reservoir engineering _ finds its roots in the theory of laser cooling of atoms  @xcite and , more related to this work , of trapped ions  @xcite . laser - cooling schemes try to design the system - environment coupling such that the dissipation yields a stationary state with reduced kinetic energy . these ideas can be taken a step further in order to non - probabilistically produce stationary states that display non - classical aspects , or a certain amount of entanglement  @xcite . for instance , by mimicking superradiance phenomena  @xcite , certain engineered decay channels yield partial entanglement in the stationary mixed state  @xcite . however , it would be highly desirable to devise dissipative protocols that provide _ maximally - entangled pure states _ asymptotically , the so - called dark states , showing fidelities comparable to those obtained by more standard unitary schemes  @xcite . this approach has been recently pursued by different groups which have shown that , provided that one can engineer a dissipation that acts quasi - locally on different parts of the system , it is possible to design dissipative protocols to produce a number of paradigmatic dark states  @xcite , or to perform universal dissipative quantum computation  @xcite . for instance , by engineering a particular dissipation that acts on pairs of adjacent atoms in an optical lattice , the system can be driven dissipatively into a superfluid phase  @xcite . in addition , as initially proposed for rydberg atoms  @xcite and first realized in experiments of trapped ions  @xcite , by concatenating multi - qubit gates with a controlled dissipation of ancillary qubits , a variety of multipartite entangled states ( i.e. bell , greenberger - horne - zeilinger , and dicke states ) have been dissipatively generated with reasonably high fidelities . let us note that this stroboscopic time - evolution corresponds to a markovian master equation in the limit of many gates and dissipation steps , where the dissipative jump operators correspond to multi - qubit jump operators . from a fundamental point of view , it is of interest to address whether similar dissipative protocols can also work in an _ ( i ) analog _ and _ ( ii ) global _ fashion ( i.e. using always - on couplings without individual addressing ) . additionally , from a more pragmatic point of view , such global analog schemes would not be limited by the accumulation of the errors in each step of the stroboscopic protocols . however , finding particular schemes to provide multi - qubit jump operators in an analog manner seems to be a tremendous task from both a theoretical and experimental point of view . therefore , we impose a further constraint , the jump operators should be composed of _ ( iii ) single - qubit operators_.    in this article , we develop an instance of a global analog dissipative protocol that generates multipartite entangled states corresponding to ground states of a quantum spin chain . the underlying idea is to engineer jump operators that are a particular sum of single - spin operators , in analogy to the models of collective spontaneous emission  @xcite . we show that , if the dissipation of the spins is mediated by a common harmonic mode , the spin chain sees a structured markovian environment . by structured environment we refer to a spectral density exhibiting a sharp peak at a certain frequency ( so called breit - wigner resonance ) , corresponding to a weakly damped harmonic oscillator . this structure will allow us to design jump operators in such a way that the stationary state of the system corresponds to the ground state of the quantum spin chain . in particular , we consider an isotropic xy spin chain , which describes a critical phase of matter . we explain in detail how to implement this protocol with trapped ions , relying on tools that have already been achieved experimentally . these tools are the so - called state - dependent forces  @xcite , and sympathetic resolved - sideband cooling developed for quantum computation  @xcite . we test the scheme for realistic parameters , and show that the fidelities that can be achieved are comparable to the unitary protocols that produce multipartite entangled states  @xcite . our procedure performs well for small chains of trapped ions , which is actually the regime in most of the current experiments . this article is organized as follows . in sec . [ dissipative_protocol ] , we describe the dissipative protocol . we start by introducing the model under consideration in sec . [ model ] , and then move onto an analytic discussion , supported by numerical results , of how to engineer a structured environment that allows us to prepare multipartite entangled states dissipatively ( sec .  [ cooling_spin_wave ] ) . in sec . [ spin_gs_cooling ] we explain how the method can be used to cool the system to the many - body ground state of the quantum spin chain . in sec . [ sec : ions ] we show how crystals of trapped atomic ions are ideally suited for an experimental implementation of our ideas . section  [ ion_crystal ] describes the trapped - ion setup , a two - species coulomb crystal . in secs . [ ion_xy_model ] ,  [ ion_spin_boson ] and  [ ion_damping ] we introduce the trapped - ion toolbox required to implement the dissipative protocol . section  [ numerical ] contains a numerical analysis of the trapped - ion dissipative protocol . finally , we present our conclusions and an outlook in sec . [ conclusions ] . and @xmath0 going respectively up and down the spectrum of spin - chain excitations with quasi - momentum @xmath1 and energy @xmath2 ( the model parameters are defined in the text ) . ] let us start by introducing the model under consideration : an interacting spin-@xmath3 chain coupled to a damped bosonic mode ( see fig . [ fig_scheme]*(a ) * ) . this model , which shall be referred to as the _ damped spin - boson chain _ ( dsbc ) in the rest of the manuscript , is a many - body generalization of a system considered recently for the dissipative generation of two - qubit entanglement  @xcite . the dsbc is described by the following master equation ( @xmath4 ) @xmath5+\mathcal{d}_{\rm b}(\rho),\ ] ] where @xmath6 is the total density matrix , and @xmath7 is the dissipative liouvillian . we consider a finite chain of @xmath8 spins evolving under the isotropic xy hamiltonian  @xcite @xmath9 where @xmath10 is an antiferromagnetic interaction strength , and we introduced the raising and lowering operators @xmath11 . the bosonic hamiltonian consists of a single mode @xmath12 where @xmath13 and @xmath14 are respectively the creation and annihilation operators ( the oscillator frequency @xmath15 can actually correspond to a detuning @xmath16 , as will become clear from the experimental implementation ) . the spins and the boson are coupled through @xmath17 in the above expression , we have introduced the site - dependent coupling strength @xmath18 , and the pauli matrices @xmath19 . finally , we also consider weak damping of the bosonic mode via a lindblad - type dissipator @xmath20 where @xmath21 is the damping rate . we note that all these individual ingredients can be realized with the state - of - the - art technology of trapped - ion crystals , as described in sec . [ sec : ions ] . the objective of this work is to understand the interplay of these dynamics , such that the liouvillian @xmath7 generates stationary entanglement between the initially uncorrelated spins . so far , we can already appreciate two of the properties of the protocol outlined in the introduction : _ ( i ) _ it is _ analog _ , since the couplings  - are considered to be switched on during the whole protocol . _ ( ii ) _ it is _ global _ , since the couplings  - address all the spins of the chain . let us note that the spin - chain hamiltonian   alone could already generate entangled states unitarily at given instants of time . nonetheless , these correlations have a transient nature , whereas the interest of this work is the onset of stationary correlations . therefore , we explore the interplay of this hamiltonian part with the irreversible dynamics introduced by the damped boson through the spin - boson coupling  .    in fig . [ fig_scheme]*(a ) * , we represent schematically the dsbc . the spin chain is considered to be our system @xmath22 , whereas the bosonic mode , together with the markovian reservoir where it dissipates , provide an effective structured environment @xmath23 for the spin chain . at first sight , the dephasing - like spin - boson coupling   seems to introduce a source of decoherence in the spin chain , thus hindering rather than assisting the generation of stationary entanglement . in the following section , we will show that this nave intuition is not always valid , and explain the subtle mechanism that allows for generation of stationary entanglement in a particular regime of the system parameters . we consider an initial state of the system @xmath24 , where @xmath25 is an arbitrary state of the harmonic mode , and the uncorrelated spin state @xmath26 has a single spin excitation at the left edge of the chain . in this section , we discuss how the system liouvillian @xmath7 acts on this state , allowing for the onset of stationary entanglement .      as a consequence of the xy hamiltonian , the initially localized spin excitation is exchanged between different neighbors . in fact , the spin chain corresponds exactly to a tight - binding model where the excitation hops along the sites of the underlying chain according to @xmath27 in this expression , @xmath28 results from the projection of the xy hamiltonian onto the single - excitation subspace @xmath29 , where @xmath30 is the projector onto an @xmath8-dimensional subspace spanned by the states @xmath31 , such that @xmath32 labels the position of the spin excitation in the chain . this tight - binding model gives rise to the following band structure ( see fig . [ fig_scheme]*(b ) * ) @xmath33 where each energy is associated to a `` spin - wave excitation '' @xmath34 the goal of this section is to find the conditions to generate dissipatively one of these multipartite entangled states , @xmath35 . we consider a spatial modulation of the spin - boson coupling strength   in the form @xmath36 , where @xmath37 . in the single - excitation subspace , the spin - boson coupling becomes @xmath38 where the ladder operators @xmath39 are responsible for ascending or descending along the ladder of spin - wave excitations ( see fig . [ fig_scheme]*(b ) * ) . accordingly , the spin - boson coupling   connects different spin - wave excitations , while simultaneously creating or annihilating bosonic quanta . moreover , since the harmonic mode dissipates irreversibly into a markovian environment  , the ladder operators @xmath40 introduce irreversible dynamics in the spin chain . in terms of these ladder operators , the single - excitation spin - chain hamiltonian   reads @xmath41 where for simplicity we reset the zero of energy at the bottom of the band . we show below how to generate mesoscopic entangled states by controlling the relative strengths of the dissipative processes that take the system up and down the _ spin - wave ladder_.      since we are interested in the stationary properties of the dsbc , we consider long times @xmath42 , such that the bosonic mode has enough time to relax under the action of the damping . besides , for @xmath43 , this relaxation is much faster than the process of energy exchange between the boson and the spins . in this limit , the boson thermalizes individually , and the spin chain evolves on a much slower timescale under such a bosonic background .    to obtain the effect of the boson background on the slower spin dynamics  @xcite , we must `` integrate out '' the bosonic degrees of freedom from the spin - boson coupling  . we start with a state of the form @xmath44 , where @xmath45 is the vacuum of the harmonic mode , and fulfills @xmath46 . we then expand the state at a later time @xmath47 , with @xmath48 , in powers of the coupling constant @xmath49 between the spin system and the oscillator , keeping up to second order in @xmath49 . tracing over the mode , we obtain : @xmath50 where the `` hat '' indicates that we work in interaction picture with respect to @xmath51 , and where we have introduced @xmath52 $ ] . using the explicit form of @xmath53 and @xmath54 , after some algebra and one integral in @xmath55 we find : @xmath56.\end{gathered}\ ] ] here , we introduced the _ collective jump operators _ @xmath57 with @xmath58 $ ] , being @xmath59 the energy difference between two neighboring spin - wave excitations in the spin - wave ladder ( see fig . [ fig_scheme]*(b ) * ) . the integrand in eq . ( [ first_integral ] ) contains some terms that oscillate rapidly in time . we now perform the remaining integral , assuming that the time interval @xmath60 is short compared to the time scale given by @xmath61 , but long enough such that @xmath62 . this condition restricts the values of @xmath49 for which this treatment is valid . to perform the integral , we group the frequencies @xmath63 in such a way that within each group the frequencies are equal , and the difference between the frequencies in different groups is finite ( we note that this is not possible in the limit of infinite sizes , where the energies become a continuum ) . it is worth noticing that the presence of degenerate frequencies is typical in this model and makes the grouping necessary . if one keeps only the dominant terms , which are the terms in the integrand that are constant in time , back in schrdinger picture one finds the following master equation governing the coarse - grained evolution over the time scales of interest : @xmath64+\mathcal{d}_{\rm s } ( \rho_{\rm s}),\ ] ] with @xmath65\,,\end{gathered}\ ] ] and where @xmath66\,.\end{gathered}\ ] ] in the expressions above , the sum over @xmath67 runs over the different transition frequencies in the system , @xmath68 denotes an anticommutator , and the lindblad operators are defined as : [ jump_operator ] j_^+ = _ n/_n= l_n^+ = j_^-^where the sum is over all the values of @xmath69 such that @xmath70 . the action of the different lindblad operators is weighted by the spectral density of the effective environment , @xmath71}.\ ] ] according to these expressions , the spins are subjected to a lorentzian reservoir centered at the boson detuning @xmath15 with a width given by the boson damping rate @xmath72 . therefore , the dissipation on the spins is not equal at all frequencies , but stronger at frequencies matching that of the bosonic mode ( i.e. they do not see a totally flat environment , but a structured one  @xcite ) .    as announced previously , the ladder operators @xmath40 are responsible for introducing the irreversible dynamics in the spin chain . in particular , they determine the collective jump operators  , where the adjective collective emphasizes that they act over all the spins in the chain . however , we remark that these jump operators are a sum of single - spin operators , as opposed to the multi - spin nature of some other engineered dissipation protocols considered recently  @xcite . with this discussion , we show the third property of the protocol outlined in the introduction : jump operators are composed of sums of _ ( iii ) single - qubit operators_. this draws an analogy to the models of collective spontaneous emission  @xcite , but we will show that the special dissipation mediated by the boson mode allows us to use them to prepare dissipatively the ground state of the spin chain .    , the jump operators that descend the ladder govern the dissipative dynamics ( i.e. @xmath73 ) . this can be understood in terms of the effective lorentzian spectral density @xmath74 of the environment seen by the spin chain , which overlaps minimally with the processes that climb up the ladder . * ( b , d ) * in the regime of negative oscillator detuning @xmath75 , the jump operators that climb up the ladder dominate the dissipative dynamics ( i.e. @xmath76 ) , since in this case the spectral density of the environment seen by the spin chain overlaps minimally with the processes that descend the ladder . ] we are searching for the steady state @xmath77 of eq .  , @xmath78 , with three properties : being pure , displaying multipartite entanglement , and being unique . a pure steady state @xmath79 necessarily belongs to a hamiltonian eigenspace  @xcite , which in our case is given by the spin - wave excitations @xmath80 . additionally , the steady state should also belong to the kernel of the dissipator , @xmath81 . by inspection of the jump operators , we identify two regimes where this can happen , which depend on the boson detuning with respect to the spin - wave energy difference : @xmath82 @xmath83 for the regime of positive detunings , we can then approximate the evolution using that @xmath84 is negligibly small for all @xmath69 , while all the @xmath85 are non - negligible . in that case the non - unitary part of the evolution is _ s _ s 2_j_a ( ) ( j_^- _ s j_^+ - \ { j_^+ j_^- , _ s } ) , taking the system to the lowest state in the manifold of spin waves ( the hamiltonian part does not induce transitions between spin waves ) . thus , the unique pure steady state corresponds to the lowest spin wave @xmath86 . conversely , for negative detunings , the steady state corresponds to the highest spin - wave @xmath87 . following  @xcite , it can be shown that in these limits the particular form of the lindblad operators does not allow for mixed steady states . we thus get the desired result : a unique pure state displaying stationary multipartite entanglement .    in fig . [ fig_2 ] , we represent schematically the two regimes of interest . when conditions   are fulfilled , for @xmath88 , the dissipative processes that go down the spin - wave ladder dominate since @xmath89 . effectively , the dissipation induced by the damped bosonic mode `` cools '' the spin state to the lowest - energy spin wave , @xmath90 ( see fig .  [ fig_2]*(a ) * ) . the bath spectral density peaks at the frequencies corresponding to the dissipative processes descending the spin - wave ladder ( see fig . [ fig_2]*(c ) * ) and therefore , the bath absorbs very efficiently the energy dissipated by the spin chain during its relaxation to the spin wave @xmath91 . conversely , when the system parameters fulfill conditions  , for negative detunings @xmath92 , it is the dissipative processes that climb up the ladder which dominate , @xmath93 . in this limit , the bosonic mode drives the spin state to the highest - energy spin wave , @xmath94 ( see fig . [ fig_2]*(b ) * ) . the bath spectral density is maximal for the ascending dissipative processes ( see fig .  [ fig_2]*(d ) * ) , and the bath provides the required energy to climb up the spin - wave ladder . we may thus conclude that in both regimes , the structured reservoir singles out a unique spin wave as the steady state of the chain , assisting the generation of stationary multipartite entanglement .    in the light of these results , we can also understand the regime @xmath95 . in this limit , the reservoir spectral density   becomes essentially flat , and both up / down processes contribute equally @xmath96 . it is straightforward to see that if @xmath97 is constant , the dissipator in eq . ( [ complete_dissipator ] ) satisfies @xmath98 , and the totally mixed state @xmath99 becomes a steady state . this regime corresponds to the nave argument of section  [ model ] , following which one expects that the dephasing - like term   can only decohere the spin chain . it is now clear that there are other regimes ,  - , where we can profit from a structured environment to assist the generation of stationary entanglement . we now comment on the time required to reach the aforementioned steady states . in the regime where the effective liouvillian   was derived , it is given by @xmath100 , where @xmath101 . for a particular experimental setup , where the spin couplings can not reach arbitrarily large values , this preparation time may turn out to be too long for practical purposes . we emphasize , however , that the same dissipative preparation of entangled states can be obtained in a non - markovian regime where @xmath102 , or @xmath103 . since simple master equations can not be derived analytically for this regime , we shall explore it numerically . let us note that although the spin - wave excitations   are genuinely multipartite entangled , the weights for the excitation of each of the spins in the system are different . a small modification of the scheme also allows for the dissipative generation of @xmath8-partite w - like states , where the w - state is defined as |w=(|+| + | ) . let us leave for a moment the constraint to consider purely @xmath104 couplings , and assume that it is possible to add external transverse fields that act locally at the two edges of the chain @xmath105 in the subspace with one spin excitation , the ground state of this hamiltonian is of the form |_n = _ i=1^n ( -1)^i |i , which is locally equivalent to the above @xmath8-partite w - state ( we note that the actual w - state can be obtained if the spin - spin coupling is ferromagnetic instead of antiferromagnetic ) . as opposed to the ground state of the spin system considered so far , the w - state has the property that it is fully symmetric under particle interchanges . in the state @xmath106 , the excitation is equally distributed over all the sites , and each of the particles is equally entangled with the rest this state can be prepared with the same method described before , at the price of introducing individual addressability in the trapped - ion proposal of sec . [ sec : ions ] . * approach to the asymptotic state : * time evolution of the error @xmath107 , where @xmath108 is the fidelity with the inhomogeneous spin - wave excitation @xmath109 for a chain of three spins . the curves correspond to the resonant case @xmath110 , while the other parameters are @xmath111 , @xmath112 ( red ) , @xmath113 , @xmath114 ( blue ) , and @xmath115 ( green ) . for these numerical calculations , the oscillator was truncated to three levels , with trace preserved to order @xmath116 . ] +    in this section , we analyze numerically the validity of our previous analytical derivations by integrating directly the master equation of the damped spin - boson chain  . we obtain the fidelity with the desired stationary entangled state @xmath117 where @xmath26 , we consider that the boson mode is initially in the vacuum state @xmath118 , and @xmath119 is the particular entangled state that the protocol targets . we test the robustness of this fidelity against a non - optimal choice of the system parameters  - , a limited protocol time @xmath120 , and an increasing particle number @xmath8 .    in figure [ fig : fidelity_dynamics ] , we consider the time evolution of the fidelity for the preparation of the lowest - energy spin - wave excitation of an @xmath121 spin chain @xmath122 for these calculations , we set the oscillator detuning @xmath15 to the value @xmath123 which is resonant with the spin - wave transitions , and vary the ratio between the spin - boson coupling strength @xmath49 and the damping rate @xmath72 . the plot indicates that , for similar values of the final fidelity , the convergence is much faster when @xmath49 and @xmath72 are of the same order . this implies that for practical purposes it is most convenient to work in the deeply non - markovian regime . the optimal ratio @xmath124 , however , depends in general on the number of spins and the value of @xmath15 . it is important to note that the resonance condition can be relaxed to some extent . in figure [ fig : fidelity_n3 ] we analyze the fidelity as a function of the detuning @xmath15 and of the coupling @xmath49 , setting @xmath125 . from our previous analysis , we expect the fidelity to be maximal when @xmath126 , and when @xmath127 approach zero to fulfill the requirement @xmath128 . this behavior is totally captured by the numerical results in fig . [ fig : fidelity_n3]*(a ) * , displaying the asymptotic fidelity . the maximum fidelity , highlighted with a white dot , is indeed found in this limit and is equal to 1 within the numerical errors . we remark that , even when the above conditions are only partially fulfilled , the dynamics can provide the desired state with fidelities well above @xmath129 . therefore , we can claim that our dissipative protocol is considerably robust against parameter imperfections . nevertheless , the time required to approach the desired target state might be prohibitively long . therefore , we also characterize the fidelity for a fixed protocol duration . in fig . [ fig : fidelity_n3]*(b ) * , we show the value of the fidelity at a finite fixed time @xmath130 . as can be seen in this figure , the optimal choice for the oscillator detuning is still @xmath123 , while the optimal value for @xmath131 has moved away from zero . in fact , the best choice now results from the interplay between the condition @xmath132 , and the required convergence speed . we emphasize that it is still possible to find parameters such that the achieved fidelities remain very high , @xmath133 . we now address the dissipative generation of w - like states by adding boundary transverse fields  . for the @xmath121 spin chain under consideration , the target state is @xmath134 in fig . [ fig : fidelity_n3]*(c)-(d ) * , we show the asymptotic and finite - time fidelities for the dissipative generation of such a tripartite entangled state with homogeneously - distributed excitations . we observe high fidelities comparable to the inhomogeneous spin - wave state  , and moreover , an analogous robustness with respect to a non - optimal choice of the system parameters . let us note that , for @xmath121 , there are two different transitions in the spectrum of these w - type spin - waves : one is resonant at a frequency of @xmath135 , and the remaining one occurs at @xmath136 . therefore , it is impossible to set an oscillator detuning at resonance with both transitions simultaneously , and thus the conditions   are not completely fulfilled . this explains the fact that the w - type fidelities achieved for finite times are lower than those corresponding to the inhomogeneous spin waves in fig . [ fig : fidelity_n3]*(b ) * , where the two transitions have the same frequency @xmath137 .    in this section , we have analyzed numerically the validity of the scheme for the generation of the simplest case of multipartite entanglement , namely tripartite entangled states . a question that should be carefully addressed is whether the same scheme can be used to generate @xmath8-partite entangled states , and how large can the attained fidelities be as @xmath8 is increased . this is the topic of the following section . we start by assessing the fulfillment of the necessary conditions  - , which rely on the energetic difference between the dissipative processes that climb up / down the spin - wave ladder , as the number of spins @xmath8 is increased . for very large spin chains @xmath138 , the energy differences between neighboring spin waves @xmath139 . this implies that the energetic argument selecting only processes going up or down the ladder can no longer hold . indeed , @xmath140 , so that the protocol ceases to be operative in the thermodynamic limit of infinitely long chains . however , for mesoscopic spin chains , the ratio @xmath141 can be controlled to an acceptable degree . in fig . [ fig : nsites]*(a ) * , we show that for positive detunings and @xmath142 , @xmath143 for the most of the spin - wave excitations , whereas slightly higher values are attained for the extremal spin excitations ( where @xmath144 ) . for ratios @xmath141 on this order , we expect that the dissipative preparation of entangled states still works with acceptable fidelities .    to be more specific , we note that the presence of different transition frequencies @xmath63 is generic for @xmath145 , which represents an obstacle for the efficiency of the procedure . as the number of spins is increased , the convergence also becomes slower because of the larger number of steps down / up the ladder towards the target state , and the lower transition frequencies which require a decrease in @xmath49 and @xmath72 . in fig . [ fig : nsites]*(b ) * , we show numerical results for the dependence of the protocol error @xmath146 achieved for a finite time @xmath147 as a function of the number of sites and optimizing the values of @xmath49 , @xmath72 and @xmath15 . we observe that the errors obtained for the dissipative state preparation of the inhomogeneous spin wave @xmath148 are below 10% for chains of up to ten sites . in comparison , the creation of the w - like states is worse , and such high fidelities can only be achieved for short chains of up to five sites .     * the dissipative protocol in mesoscopic chains : * * ( a ) * ratio @xmath149 of the relative strength of the dissipative processes climbing up and going down the ladder as a function of quasimomentum index @xmath69 for different number of spins @xmath8 . the values @xmath150 correspond to negative oscillator detunings @xmath151 which select the processes up the ladder , whereas @xmath152 corresponds to positive oscillator detunings @xmath153 selecting the cooling processes down the ladder . * ( b ) * error @xmath154 for the dissipative generation of the target state in a finite time @xmath147 as a function of the number of spins @xmath8 in the chain . the blue circles represent the errors corresponding to the inhomogeneous spin waves @xmath155 . the red diamonds represent the error for the w - like states @xmath156 . for the numerical calculation , the hilbert space of the oscillator was truncated to four levels . ] so far , our analysis has been restricted to an @xmath8-dimensional subspace of the spin chain , since the initial spin state @xmath26 contains a single excitation , and the number of spin excitations is preserved by the complete dsbc liouvillian  . in the following we explore the full @xmath157-dimensional hilbert space of the spin chain .      in order to treat the full hilbert space of the spins , we fermionize the spin chain via the so - called jordan - wigner transformation  @xcite , namely @xmath158 where @xmath159 are fermionic creation - annihilation operators . the spin - chain   and spin - boson   hamiltonians can be expressed in terms of the jordan - wigner fermions as follows @xmath160 where we recall that the spin - boson couplings are @xmath36 with @xmath37 . the next step is to express these hamiltonian terms in the spin - wave basis introduced in eq .  , which leads to the following expressions @xmath161 where the fermionic operators in momentum space are @xmath162 . in combination with the dissipative part given by eq .  , the spin - boson system is analogous to a _ damped single - mode holstein model _ , a dissipative version of the familiar holstein model describing electron - phonon interactions  @xcite . we can obtain the same formal expressions as in the single - excitation problem by rewriting the ladder operators in second - quantized form @xmath163 accordingly , in the regime where the boson degrees of freedom can be integrated out , we obtain a purely fermionic master equation which coincides with eqs .  - , but with the collective jump operators now expressed in terms of the jordan - wigner ladder operators @xmath164 as a result , the dissipative dynamics restricted to the single - excitation sector can be generalized directly to the full hilbert space with arbitrary numbers of spin excitations . let us now consider an initial state with an arbitrary number of spin excitations @xmath165 distributed along the chain @xmath166 we would like to determine the steady - state of the spin - boson system if the conditions   are fulfilled . in this limit , the lindblad operators are only of the form @xmath167 such a jump operator has the effect of lowering the energy of the fermionic quasiparticles . due to the pauli exclusion principle , the initial @xmath168 excitations can not all occupy the lowest energy level , with quasimomentum @xmath169 . instead , the stationary state must be of the following form @xmath170 this is precisely the ground state of the original isotropic xy model  , if supplemented with a homogeneous transverse field @xmath171 . in this hamiltonian , the transverse field plays the role of an effective chemical potential @xmath172 , which is determined by the initial number of spin excitations . due to the to the dissipative process , the excitations are distributed in the lowest available single - particle states ( see fig . [ fig_6]*(a ) * ) . this means , for each number of spins up in the initial state prepared , there is a value of the transverse field @xmath173 such that the asymptotic state corresponds to the ground state of @xmath174 . conversely , a given choice of @xmath175 determines the number of spins that should be up in the initial state so that the dissipative dynamics take the system into the desired ground state . spin excitations to the lowest non - occupied state . the stationary state corresponds to the ground state of the isotropic xy hamiltonian , with the particular filling determined by an additional transverse field acting as a chemical potential . * ( b ) * the same as in * ( a ) * but selecting the processes that climb up the ladder . the stationary state would correspond to the zero - temperature state of the same xy model in a transverse field , but with reversed coupling strengths . ] we could also explore the steady state when the conditions   are fulfilled . in this limit , the only lindblad operators are of the form : @xmath176 and pump all the excitations to the highest - energy available single - particle states @xmath177 corresponding to the ground - state of the xy hamiltonian @xmath178 with a ferromagnetic spin - spin coupling , and an inverted chemical potential ( see fig . [ fig_6]*(b ) * ) .       * scaling of the ground - state - cooling protocol : * error @xmath179 in producing the desired ground - state @xmath180 after a fixed time @xmath147 and for a number of spins ranging from @xmath181 to @xmath182 . for each number of sites , we plot in blue circles the average over the different possible numbers of excitations ( @xmath183 $ ] ) , and with yellow and red squares the best and worst cases , respectively . for the numerical calculation , the hilbert space of the oscillator was truncated to three levels . ] let us now confirm the validity of the above argument by numerical analysis for small spin chains with different numbers of initial excitations . in fig . [ fig : xx ] , we display the error in producing the ground - state   taking a fixed protocol time of @xmath147 . the highest errors occur when the number of excitations is 1 or @xmath184 , since the transition frequencies are lowest at these points ( this is the worst case , with the errors plotted in fig . [ fig : xx ] by means of red squares ) . with yellow squares , we represent the lowest attained errors . note that the average errors of the protocol ( blue circles , averaged over all the possible excitation numbers ) are as low as @xmath185 for chains up to @xmath182 spins , which would thus allow us to prepare dissipatively the ground state of the spin model with fidelities of @xmath186 . let us close this section by noting that , in the general case with an arbitrary number of spin excitations , a different choice for the site - dependence of the coupling coefficients @xmath18 may allow for jump operators that go down / up the spin ladder in larger steps , and therefore favor a faster convergence to the asymptotic state ( this appears to be specially well - suited for chains around the half - filling condition ) . however , for up to @xmath182 spins , we have not found any sizable advantage , and the highest fidelities were always achieved for transitions between neighboring spin waves @xmath187 . let us describe a specific platform where the dsbc model   can be realized . we consider a system of @xmath188 ions confined in a linear radio - frequency trap  @xcite , and arranged forming a one - dimensional chain ( see fig . [ fig_ion_scheme]*(a ) * ) . the vibrational degrees of freedom around the equilibrium positions can be expressed in terms of the so - called normal modes  @xcite , namely @xmath189 where @xmath190 refers to the trap axis , @xmath191 labels the different modes with frequencies @xmath192 , and @xmath193 are bosonic operators that annihilate / create phonons in a particular mode . in a linear chain of ions with equal mass , the phonon branches have the structure shown in fig .  [ fig_ion_scheme]*(b ) * , such that the different modes spread around the trap frequencies lying in the range @xmath1941 - 10mhz , @xmath1951 - 10mhz , and @xmath1960.1 - 1mhz . a property of the longitudinal branch that will be used in this work is the existence of an energy gap between the lowest - energy mode ( i.e. the center - of - mass mode ) and the following one .    as shown in fig . [ fig_ion_scheme]*(a ) * , @xmath8 ions of the chain correspond to a particular atomic species with hyperfine structure , whereas the remaining @xmath197 ions , that will be used for cooling , do not necessarily have hyperfine splitting ( we note that though we focus on hyperfine qubits , our proposal could be generalized to optical or zeeman qubits ) . for simplicity , in the following we describe the case @xmath198 . for the first @xmath8 ions , we select two hyperfine levels of the ground - state manifold , which are referred to as spin states @xmath199 , and have a transition frequency @xmath200 in the microwave range with a negligible linewidth @xmath201 ( fig .  [ fig_ion_scheme]*(c ) * ) . for the remaining @xmath202-th ion , we select two internal states @xmath203 from a certain transition with frequency @xmath204 , and linewidth @xmath20510 - 100@xmath206khz ( fig . [ fig_ion_scheme]*(d ) * ) . this is the so - called resolved sideband limit @xmath207 , which will allow for an efficient laser cooling close to the vibrational ground state for the longitudinal center - of - mass mode @xcite . let us note that these small linewidths can be obtained by a variety of methods depending on the particular ion species ( e.g. for weakly - allowed dipole transitions , for quadrupole - allowed transitions with an additional laser that admixes the excited state with that of a dipole - allowed transition , or by a raman configuration )  @xcite . the complete degrees of freedom are thus described by the following master equation @xmath208+\tilde{\mathcal{d}}(\rho),\ ] ] where we have introduced the hamiltonians @xmath209 such that @xmath210 , and @xmath211 ( we use @xmath212 instead of @xmath213 to emphasize the pseudo - spin character of these internal states ) . additionally , we have to take into account the dissipation from the dipole - allowed transition considering recoil effects  @xcite . for trapped ions , this can be described by the following super - operator @xmath214 , where @xmath215 where @xmath216 . this is the usual dissipator of a two - level atom in an electromagnetic environment . the spontaneous emission also has effects on the vibrational degrees of freedom due to the photon recoil . to leading order in the lamb - dicke parameter @xmath217 for each normal mode , the jump operators create and annihilate phonons as given by @xcite : @xmath218\tilde{\tau}_{n+1}^++\text{h.c . } , \end{split}\ ] ] where we have introduced @xmath219 , with @xmath220 the coefficients describing the displacement of the @xmath221 ion in mode @xmath69 along direction @xmath222 . once the vibrational and atomic degrees of freedom have been introduced , let us summarize the role that each of them will play in the realization of the damped spin - boson chain  : _ ( i ) _ the hyperfine levels of the @xmath8 ions will simulate the spins in the dsbc model . _ ( ii ) _ the @xmath223 vibrational modes along the @xmath224 and @xmath225 axes will be used to mediate the spin - spin interactions leading to the isotropic xy model in eq .  . _ ( iii ) _ the longitudinal center - of - mass mode will play the role of the harmonic oscillator in the dsbc model . _ ( iv ) _ the atomic levels of the @xmath202-th ion will be used to laser - cool the longitudinal center - of - mass mode . we note that the position of this auxiliary ion is irrelevant for the purpose of cooling , but it modifies the spin couplings in the chain . moreover , it is also possible to use more than one auxiliary ion as an additional knob to control the damping rate . in the following sections , we will detail the different elements required for the ion - trap implementation of the dsbc . we emphasize that all of them can be realized with state - of - the - art technology . we start by addressing how to implement the isotropic xy model   in the ion chain by means of a variation of the mlmer - srensen gate @xcite . the idea is to use the so - called spin - dependent forces , a tool that has already been implemented in several laboratories for different purposes  @xcite . by combining a pair of raman laser beams ( see fig . [ fig_ion_scheme]*(c ) * ) , each of which is tuned to the red / blue first vibrational sideband of the hyperfine transition  @xcite , it is possible to obtain the following state - dependent force @xmath226 where the `` hat '' indicates the interaction picture with respect to @xmath227 . in the above expression , we have assumed that the effective wavevector of the raman beams @xmath228 points along the @xmath229 axis of the trap ( see fig . [ fig_ion_scheme]*(a ) * ) . additionally , we assume that the strength of both raman beams is equal , such that they share a common rabi frequency @xmath230 , and have opposite detunings , such that we can define a common @xmath231 . we have introduced the laser lamb - dicke parameter @xmath232 . the constraints on these parameters are @xmath233 to neglect other terms in addition to such a state - dependent force . finally , we have defined the sum and difference of the raman beam phases @xmath234 and the pauli spin operator @xmath235 .    in order to obtain the desired isotropic xy model , we apply two orthogonal state - dependent forces with the following phases and directions @xmath236 altogether , the hamiltonian of the ion chain in the interaction picture becomes @xmath237    performing a magnus expansion  @xcite , the evolution turns out to be given by the unitary operator @xmath238 , where @xmath239 displays the following form to second order @xmath240.\ ] ] we want to derive an effective hamiltonian from this expression @xmath241 ( where the superscript `` ti '' stands for the trapped - ions realization ) . the first - order contribution leads to a couple of orthogonal state - dependent displacements , which can only be neglected in the limit @xmath242  @xcite . in addition , from the non - commutativity of the @xmath243 and @xmath244 forces , we get a residual spin - phonon coupling in the second - order term of the magnus expansion . in order to neglect it , we have to impose a further constraint , namely @xmath245 , which can be fulfilled if the trap frequencies along the @xmath246 axes are sufficiently different @xmath247 . under these constraints , an interacting quantum spin chain is obtained @xmath248 therefore , by adjusting the strengths of the raman beams , and their detunings , it is possible to find a regime where @xmath249 , and the above spin hamiltonian corresponds to the desired isotropic xy model since @xmath250 . let us remark that the resulting xy model has the peculiarity of displaying long - range couplings . in fact , when @xmath251 , and the raman lasers are far - detuned from the whole vibrational branch , it can be shown that these couplings decay with a dipolar law . the trapped - ion hamiltonian   becomes a realization of the xy interaction   in the dsbc model   after the identifications @xmath252 , and @xmath253 . for typical nearest - neighbor distances of @xmath2541 - 10@xmath173 m , the spin - spin couplings attain strengths in the @xmath2551 - 10 khz . let us note that modifications of the scheme of state - dependent forces have been proposed to yield other quantum spin models  @xcite , some of which have also been realized experimentally  @xcite . [ ion_spin_boson ]    let us now address how to implement the spin - boson coupling   in the ion chain . the idea is again to use a @xmath256-beam configuration ( see fig . [ fig_ion_scheme]*(c ) * ) , but in a different regime . rather than tuning the two - photon frequencies @xmath257 to the vibrational sidebands of the hyperfine transition , we impose that @xmath258 . moreover , we consider that the laser beams form a standing wave along the trap @xmath259-axis ( see fig . [ fig_ion_scheme]*(a ) * ) . under these constraints , the laser - ion interaction leads to a crossed - beam differential ac - stark shift , which can be interpreted as another state - dependent force in the @xmath260 basis @xmath261 where the `` hat '' refers to the interaction picture with respect to @xmath262 . here , the parameters are defined in analogy to those in eq .  , with three important differences : _ ( i ) _ @xmath263 is not the two - photon rabi frequency of the hyperfine transition , but rather a differential ac - stark shift coming from processes where a photon is exchanged between the pair of laser beams in the @xmath256 configuration . _ ( ii ) _ the detunings are changed to @xmath264 . _ ( iii ) _ since the ion chain lies along the @xmath259 axis , there is an additional site - dependent phase when shining the lasers such that @xmath265 . when adjusting the difference of the laser phases to be @xmath266 , then we obtain @xmath267 . let us highlight that the standing - wave nature of this spin - dependent force is not essential for the dissipative protocol . however , it makes the connection with the dsbc studied in sec . [ dissipative_protocol ] more transparent . we will show in sec . [ numerical ] that essentially the same fidelities can be achieved for a traveling - wave configuration , which has also been experimentally demonstrated  @xcite . we now exploit the particular properties of the longitudinal phonon branch ( see fig . [ fig_ion_scheme]*(b ) * ) . more precisely , we use the presence of a gap between the axial center - of - mass mode and the rest of the modes along the same direction . for the case of equal ions , @xmath268 , while for different ion species or isotopes , the exact value of the gap will depend on the mass ratio and the position of the cooling ions . we assume that @xmath269 is red - detuned with respect to the axial center - of - mass mode @xmath270 , such that @xmath271 . since the next vibrational modes are separated by a large energy gap , the laser coupling to the remaining phonon branch is highly off - resonant and can be thus neglected . accordingly , this state - dependent force gives _ l_z=_i_i1^z_i^zb_1ze^-i_1zt+ h.c .. we can finally move to a different picture where the above hamiltonian becomes time - independent @xmath272 , where [ spin_boson_ion ]    h_b^ti&=_i_1zb_1z^b_1z , + h_sb^ti&=_i_i1^z_i^z(b_1z+b_1z^ ) . therefore , it is clear that the center - of - mass mode plays the role of the harmonic oscillator @xmath273 in the dsbc , with the laser detuning corresponding to the oscillator detuning @xmath274 in the dsbc model  . in addition , the strength of the state - dependent force determines the spin - boson coupling @xmath275 in the dsbc model  . let us now comment on realistic values for the trapped - ion parameters . since the laser detuning should only fulfill @xmath2761 - 10@xmath206mhz , it will be easy to reach the required condition of @xmath2771 - 10@xmath206khz . on the other hand , we know from the previous section that @xmath278 is necessary to have an accurate pumping to the desired entangled state . nevertheless , @xmath49 should not be too small that the total preparation time becomes prohibitively long . a suitable choice could be @xmath279 0.1 - 0.5 khz . since @xmath280 , it will suffice to set the rabi frequency in the @xmath2811 - 5@xmath282khz . finally , we should adjust the laser wavevector such that @xmath283 . once the implementation of the xy model   and the spin - boson coupling   has been described , let us turn into the last required ingredient : an effective damping of the bosonic mode  . as discussed previously , the idea is to exploit the atomic levels of the @xmath284-th ion to laser - cool the longitudinal center - of - mass mode . since the resonance frequencies are very different , @xmath285 , the cooling lasers do not affect the spin dynamics of the other @xmath8 ions . however , since the vibrational modes are collective , it is possible to sympathetically cool the vibrations of the crystal by only acting on the @xmath284-th ion . this sympathetic laser cooling  @xcite has been already realized experimentally in small crystals for quantum computation  @xcite . our starting point is the master equation   with the dissipation superoperators  - . to control the effective damping , we introduce a laser beam red detuned from the atomic transition @xmath286 , such that @xmath287 ( see fig . [ fig_ion_scheme]*(d ) * ) . let us note that the particular laser - beam configuration will depend on the particular ion species , and the scheme to attain the resolved - sideband limit  @xcite . we consider that the laser is in a traveling - wave configuration , such that its wavevector is aligned parallel to the trap axis @xmath288 ( see fig . [ fig_ion_scheme]*(a ) * ) . after expanding in series of the lamb - dicke parameter @xmath289 , the laser - ion interaction leads to two different terms , the so - called carrier term @xmath290 and the red and blue sideband terms @xmath291 where we introduced @xmath292 , and the corresponding rabi frequency @xmath293 . by controlling appropriately these sideband terms , and their interplay with the atomic spontaneous emission , it is possible to tailor the damping of the longitudinal modes . let us rearrange the full liouvillian   as a sum of two terms @xmath294 , where @xmath295+\tilde{\mathcal{d}}^0(\bullet),\\ \mathcal{l}_1(\bullet)&=-{{\rm i}}[\tilde{h}_{\tilde{\rm l}_z}^1,\bullet]+\tilde{\mathcal{d}}^1(\bullet ) , \end{split}\ ] ] which is justified for small lamb - dicke parameters . to obtain the effective damping of the longitudinal center - of - mass mode , we use a similar formalism as in @xcite . in this case , the fastest timescale in the problem is given by the decay rate @xmath296 of the atomic states of the @xmath202-th ion . therefore , we must `` integrate out '' these atomic degrees of freedom , which can be accomplished by projecting the density matrix of the full system into @xmath297 , where @xmath298 fulfills @xmath299 , and @xmath300 is the reduced density matrix for the @xmath8 spins and the chain vibrational modes . for this we use the expression @xmath301 where the `` hats '' refer to the interaction picture with respect to @xmath262 . after some algebra , one gets the effective damping of the longitudinal modes @xmath302+{\mathcal{d}}_{\rm b}(\rho_{n}),\ ] ] where we have introduced the bosonic dissipator @xmath303 together with the effective rates for laser cooling and heating @xmath304 here , the so - called diffusion coefficient accounting for the recoil heating is @xmath305 , and the spectral functions of the sideband terms are @xmath306 which can be obtained by means of the quantum regression theorem . the formalism and cooling rates coincide with those of a single trapped ion  @xcite , with the difference that in the diffusion coefficient @xmath307 and forces @xmath308 one must consider the normal vibrational modes at the position of the cooled ion . the idea is to work in the resolved sideband regime @xmath207 , and set the laser - cooling parameters to optimize the cooling of the longitudinal center - of - mass mode @xmath309 . in this regime , the cooling dominates in  , and we are left with the desired damping of the bosonic mode @xmath310 once more , after the identifications @xmath273 , and @xmath311 , we obtain an analogous dsbc damping  . let us finally comment on the required cooling rates . ideally , we should achieve the regime @xmath3120.1 - 0.5 khz . since such cooling is not particularly fast , it is possible to find the right detuning @xmath287 such that @xmath3130.1 - 0.5 khz . besides , since we are in the resolved sideband limit , the heating can be made much smaller , @xmath314 . the stationary state of a harmonic oscillator under the action of @xmath315 in eq .   is a thermal state with mean number of quanta @xmath316 . in the limit @xmath317 , such that the cooling is the dominant effect and the heating only presents a small correction , the equilibrium state of the center - of - mass mode is close to the vacuum state .      in this last section , we will explore numerically how the dissipative protocol described in the part  [ dissipative_protocol ] of this manuscript can be implemented with the trapped - ion dscb model in eqs .  ,   and  , namely @xmath318+\mathcal{d}_{\rm b}^{\rm ti}(\rho).\ ] ] therefore , the results discussed below correspond to the actual spin - spin couplings @xmath319 in an ion trap , which are not nearest - neighbors couplings , but display a dipolar decay with the cube of the distance . furthermore , we have taken into account that in harmonic ion traps , the interparticle distance is not constant over the chain . more importantly , we have considered the effects of the always - present heating term in the trapped - ion effective damping  .     * heating in the trapped - ion dissipative protocol : * * ( a ) * fidelity with the desired target state for a fixed time @xmath147 and as a function of the ratio @xmath320 between the sympathetic heating and cooling rates in a chain of six ions in which the central four represent sites of the isotropic xy spin chain . the blue dots correspond to the case of @xmath321 initial spin excitations , and the orange squares to the case of @xmath322 excitations . the oscillator was truncated to five levels . * ( b ) * same as in * ( a ) * , but considering a shorter time @xmath323 , and a traveling - wave configuration of the spin - boson coupling in eq .  . the achieved fidelities for this experimentally simpler configuration are similar to the ones in * ( a)*. ]    in fig . [ fig:6 ion chain]*(a ) * , we show the optimal fidelity for the dissipative generation of the ideal ground state at different fillings  , by integrating numerically the dipolar and inhomogeneous trapped - ion dsbc model   as a function of the ratio between heating and cooling rate @xmath324 . we study a chain of @xmath325 ions , in which @xmath326 have a hyperfine structure and play the role of the spins in the xy chain , while the two peripheral ions are assumed to be used for sympathetic cooling ( this choice of auxiliary cooling ions at the ends of the chain makes the interparticle spacing slightly more homogeneous for the `` system '' ions ) . as can be seen in the figure , when @xmath327 , the fidelities with the target ground state approach 100@xmath328 , which allows us to conclude that the differences due to the inhomogeneity of the chain and the dipolar range of interactions with respect to the ideal homogeneous and nearest - neighbor dsbc   do not compromise the fidelities severely . the situation is different as finite heating rates @xmath329 are considered . in fig . [ fig:6 ion chain]*(a ) * , we observe that the fidelity with the target ground state decays essentially linear with @xmath324 . in the worst case considered , where @xmath330 leads to a mean phonon number of @xmath331 , the fidelities obtained are reduced to roughly @xmath332 . we note that resolved sideband cooling of single vibrational modes has been experimentally achieved reaching @xmath333 below 0.1 for single trapped ions  @xcite . sympathetic cooling of crystals up to @xmath326 ions has also been achieved  @xcite , reaching values as low as @xmath334 with pulsed techniques . we note that our scheme does not demand such an accurate ground - state cooling ( see fig . [ fig:6 ion chain]*(a ) * ) , and we expect that the required mean phonon numbers @xmath335 can also be achieved in a continuous cooling scheme . another source of experimental imperfections is given by slow drifts of the trap frequencies which will modify the detuning @xmath336 in different experimental runs , such that the conditions   will not be perfectly fulfilled . nonetheless , in fig . [ fig : fidelity_n3 ] we have shown that the dissipative character of the protocol endows it with a natural robustness with respect to non - optimal choices of the parameters . therefore , one can expect large fidelities even in the presence of slow drifts of the trap frequency within the khz - range . in any case , to alleviate imperfections caused by both magnetic - field dephasing and trap frequency drifts , we have considered a faster protocol where the required time is @xmath3371 - 10@xmath206ms . in fig . [ fig:6 ion chain]*(b ) * , we show that the achieved fidelities in this case are still above 0.8 ( for heating / cooling ratios below 0.2 ) . moreover , in this figure we have also considered substituting the standing - wave force in eq .   by a traveling - wave configuration , which is experimentally less demanding . the need of the @xmath260 state - dependent force   forbids the use of magnetic - field insensitive states ( i.e. clock states ) . therefore , the hyperfine spins will be subject to fluctuations induced by non - shielded external magnetic fields , typically leading to dephasing in a timescale of @xmath3381 - 10@xmath206ms . these timescales are comparable to the protocol times @xmath3391 - 10@xmath206ms , where we have taken @xmath3401 - 10@xmath206khz . nevertheless , these magnetic - field fluctuations act globally for standard radio - frequency traps @xmath341 , where @xmath342 is a fluctuation of the resonance frequency due to the zeeman shift ( note that this might not be the case for micro - fabricated surface traps , where fluctuating magnetic - field gradients can also arise ) . since the dynamics of the dsbc   conserves the number of spin excitations , this magnetic - field noise only introduces a global fluctuating phase , and thus does not decohere the state of the system . therefore , as far as the induced spin dynamics occurs within any subspace with a conserved number of excitations , the dissipative protocol is robust to global magnetic - field noise . * anisotropy in the trapped - ion dissipative protocol : * * ( a ) * fidelity with the desired target state for a fixed time @xmath343 and as a function of the anisotropy ratio @xmath344 for a chain of six ions in which the central four ions represent sites of the slightly - anisotropic xy spin chain . the blue dots correspond to the case of @xmath321 initial spin excitations , and the orange squares to the case of @xmath322 excitations . the oscillator was truncated to five levels . ] it is , however , important to stabilize the intensity of the state - dependent forces to the sweet spot where @xmath249 , and make sure that @xmath345 to neglect any term that does not conserve the number of excitations . the effect of off - resonant terms in the scheme leading to the xy interactions , that could give rise to non - conservation of the number of spin excitations , should be much smaller than in stroboscopic proposals @xcite since our protocol does not require fast gates @xcite . in fig . [ fig:7 ion chain ] , we show the achieved fidelities for anisotropies in @xmath346 . the results have been obtained numerically for the same regime as in fig . [ fig:6 ion chain]*(b ) * , and indicate that if @xmath347 differ by 1 part in @xmath348 , the fidelities are not severely reduced . the results in fig . [ fig:7 ion chain ] also show a marked difference in the robustness depending on the number of spin excitations in the initial state . this behaviour , at first surprising , can be understood by looking at the small undesired hamiltonian terms responsible for the anisotropy in an interaction picture with respect to the ideal isotropic hamiltonian . the anisotropy term creates or destroys fermionic quasiparticles in pairs , and rotates with a frequency that is the sum of the energies of the two quasiparticles . for the case with only one spin up , it is possible to create from the ideal target state pairs of fermions with energy adding up to zero , so that these error terms do not rotate . in the case with two spin excitations , on the contrary , the lower half of the fermionic spectrum is filled , so that any hamiltonian term creating or destroying a pair of fermions from the ideal target state oscillates in time , and therefore its effect is strongly suppressed . we have presented a method to dissipatively generate multipartite - entangled states corresponding to the ground states of small spin chains with isotropic xy hamiltonian in a transverse field . the protocol for dissipative state preparation is interesting from a fundamental point of view , as it illustrates how local and even purely markovian noise can assist entanglement production . indeed , the jump operators required are sums of terms acting on individual sites , in contrast with the few - body ( quasi - local ) nature of the operators in @xcite . when the noise deviates from exact markovianity , the required time for achieving the steady state is reduced , which illustrates the value of non - markovian effects for practical purposes . we would like to emphasize that the general idea underlying the procedure presented in this work is not specific to the model considered , and might be generalized to other spin hamiltonians . as an example , we introduced a variation that can be used to generate states locally equivalent to w - states . it would also be interesting to modify our dissipative protocol to prepare the ground state of a gapped spin model . for longer ion chains , the combination of more state - dependent forces   with different wavevectors could improve the scalability of the protocol . we have also explained in detail how to implement this method in small chains of trapped ions . in our proposal , two internal levels of the ions embody the spin system and a collective motional mode represents the damped oscillator . the preparation procedure requires the implementation of spin - spin interactions using the so - called mlmer - srensen scheme , the action of a state - dependent force , and sympathetic cooling , all ingredients within the capabilities of present ion - trap technology . numerical simulations including different sources of errors indicate that the protocol can produce the target states with fidelities comparable to the more standard coherent protocols  @xcite . the method presented is indeed robust against a number of experimental imperfections , and the implementation is simplified since only global addressing of the ions is required .    our results , however , may find a realization in different experimental setups . for instance , the application of our ideas to arrays of superconducting qubits in stripline resonators seems feasible . morever , since this model can also be understood as a quadratic fermonic model with an additional chemical potential fixing the number of particles , the results may also be interesting for fermionic atoms in optical lattices . the realization of this kind of dissipative system paves the way for the implementation of more complex scenarios to study the interplay of coherent and incoherent dynamics giving rise to noise - induced criticality @xcite . indeed , recent work in the area of non - equilibrium quantum phase transitions @xcite has shown a number of fascinating results that could be demonstrated in systems of trapped ions using tools similar to the ones in our protocol . d. e. browne , m. b. plenio , and s. f. huelga , http://prl.aps.org/abstract/prl/v91/i6/e067901[phys . * 91 , * 067901 ( 2003 ) ] ; l .- m . duan and h. j. kimble , http://prl.aps.org/abstract/prl/v90/i25/e253601[phys . rev . lett . * 90 , * 253601 ( 2003 ) . ] a. beige , d. braun , b. tregenna , and p. l. knight , http://prl.aps.org/abstract/prl/v85/i8/p1762_1[phys . * 85 , * 1762 ( 2000 ) ] ; j. pachos and h. walther , http://prl.aps.org/abstract/prl/v89/i18/e187903[phys . rev . lett . * 89 , * 187903 ( 2002 ) ] ; a. beige , d. braun , s. bose , s.f . huelga , p.l . knight , m.b . plenio and v. vedral , http://www.tandfonline.com/doi/abs/10.1080/09500340008232183[j . mod . opt . * 47 * , 2583- 2598 ( 2000 ) ] . j. f. poyatos , j. i. cirac , and p. zoller , http://prl.aps.org/abstract/prl/v77/i23/p4728_1[phys . * 77 , * 4728 ( 1996 ) ] ; m. b. plenio and s. f. huelga , http://prl.aps.org/abstract/prl/v88/i19/e197901[phys . rev . lett . * 88 , * 197901 ( 2002 ) ] . a. gonzalez - tudela , d. martin - cano , e. moreno , l. martin - moreno , c. tejedor , and f.j . garcia - vidal , http://prl.aps.org/abstract/prl/v106/i2/e020501[phys . * 106 , * 020501 ( 2011 ) ] ; a. bermudez , t. schaetz , and m. b. plenio , http://prl.aps.org/abstract/prl/v110/i11/e110502[phys . lett . * 110 * , 110502 ( 2013 ) ] ; m. j. kastoryano , f. reiter , and a. s. srensen , http://prl.aps.org/abstract/prl/v106/i9/e090502[phys . rev . * 106 * , 090502 ( 2011 ) ] . j. t. barreiro , m. mller , p. schindler , d. nigg , t. monz , m. chwalla , m. hennrich , c. f. roos , p. zoller , and r. blatt , http://www.nature.com/nature/journal/v470/n7335/full/nature09801.html[nature * 470 , * 486 ( 2011 ) ] ; p. schindler , m. mller , d. nigg , j. t. barreiro , e. a. martinez , m. hennrich , t. monz , s. diehl , p. zoller , and r. blatt , http://www.nature.com/nphys/journal/v9/n6/abs/nphys2630.html[nature phys . * 9 * , 361 ( 2013 ) ] .    c. monroe , d.m . meekhof , b.e . king , and d.j . wineland , http://www.sciencemag.org/content/272/5265/1131.abstract?sid=cc4c0d27-8270-4f9d-bd28-4e7af11fb993[science * 272 , * 1131 ( 1996 ) ] ; p. c. haljan , k .- a . brickman , l. deslauriers , p. j. lee , and c. monroe , http://prl.aps.org/abstract/prl/v94/i15/e153602[phys . * 94 , * 153602 ( 2005 ) ] ; m.j . mcdonnell , j.p . home , d.m . lucas , g. imreh , b.c . keitch , d.j . szwer , n.r . thomas , s.c . webster , d.n . stacey , a.m. steane , http://prl.aps.org/abstract/prl/v98/i6/e063603[phys . lett . , * 98 , * 063603 ( 2007 ) ] ; u. poschinger , a. walther , k. singer , and f. schmidt - kaler , http://prl.aps.org/abstract/prl/v105/i26/e263602[phys . rev . lett . * 105 , * 263602 ( 2010 ) . ] sackett , d. kielpinski , b.e . king , c. langer , v. meyer , c.j . myatt , m. rowe , q.a . turchette , w.m . itano , d.j . wineland , and c. monroe , http://www.nature.com/nature/journal/v404/n6775/abs/404256a0.html[nature * 404 , * 256 ( 2000 ) ] ; d. leibfried , b. demarco , v. meyer , d. lucas , m. barrett , j. britton , w. m. itano , b. jelenkovic , c. langer , t. rosenband , and d. j. wineland , http://www.nature.com/nature/journal/v422/n6930/abs/nature01492.html [ nature * 422 , * 412 ( 2003 ) ] ; j. benhelm , g. kirchmair , c. roos , and r. blatt , http://www.nature.com/nphys/journal/v4/n6/full/nphys961.html [ nat . phys . * 4 , * 463 ( 2008 ) ] ; a. noguchi , s. haze , k. toyoda , and s. urabe , http://prl.aps.org/abstract/prl/v108/i6/e060503[phys . * 108 , * 060503 ( 2011 ) ] .    c. ospelkaus , u. warring , y. colombe , k. r. brown , j. m. amini , d. leibfried , and d. j. wineland , http://www.nature.com/nature/journal/v476/n7359/full/nature10290.html[nature * 476 , * 181 ( 2011 ) ] ; a. khromova , ch . piltz , b. scharfenberger , t. f. gloger , m. johanning , a. f. varon , ch . wunderlich , http://prl.aps.org/abstract/prl/v108/i22/e220502[phys . * 108 , * 220502 ( 2012 ) . ] a. friedenauer , h. schmitz , j. t. glueckert , d. porras , and t. schaetz , http://www.nature.com/nphys/journal/v4/n10/full/nphys1032.html[nat . phys . * 4 , * 757 ( 2008 ) ] ; k. kim , m .- s . chang , s. korenblit , r. islam , e. e. edwards , j. k. freericks , g .- d . lin , l .- m . duan , and c. monroe , http://www.nature.com/nature/journal/v465/n7298/full/nature09071.html[nature * 465 , * 590 ( 2010 ) ] ; r. islam , e. e. edwards , k. kim , s. korenblit , c. noh , h. carmichael , g .- d.lin , l .- m . duan , c .- c . joseph wang , j. k. freericks , and c. monroe , http://www.nature.com/ncomms/journal/v2/n7/full/ncomms1374.html[nat . * 2 , * 377 ( 2011 ) ] ; b. p. lanyon , c. hempel , d. nigg , m. mller , r. gerritsma , f. zhringer , p. schindler , j. t. barreiro , m. rambach , g. kirchmair , m. hennrich , p. zoller , r. blatt , and c. f. roos , http://www.sciencemag.org/content/334/6052/57[science * 334 , * 57 ( 2011 ) ] ; j. w. britton , b. c. sawyer , a. c. keith , c .- c . j. wang , j. k. freericks , h. uys , m. j. biercuk and j. j. bollinger , http://www.nature.com/nature/journal/v484/n7395/full/nature10981.html[nature * 484 , * 489 ( 2012 ) . ]    h. schmitz , r. matjeschk , ch . schneider , j. glueckert , m. enderlein , t. huber , and t. schaetz , http://prl.aps.org/abstract/prl/v103/i9/e090504[phys . * 103 , * 090504 ( 2009 ) ] ; f. zhringer , g. kirchmair , r. gerritsma , e. solano , r. blatt , and c. f. roos , http://prl.aps.org/abstract/prl/v104/i10/e100503[phys . rev . lett . * 104 , * 100503 ( 2010 ) . ] h. rohde , s. t. gulde , c. f. roos , p. a. barton , d. leibfried , j. eschner , f. schmidt - kaler , r. blatt , http://iopscience.iop.org/1464-4266/3/1/357/[j . b * 3 , * s3 ( 2001 ) ] ; m. d. barrett , b. demarco , t. schaetz , v. meyer , d. leibfried , j. britton , j. chiaverini , w. m. itano , b. jelenkovic , j. d. jost , c. langer , t. rosenband , and d. j. wineland , http://pra.aps.org/abstract/pra/v68/i4/e042302[phys . a * 68 , * 042302 ( 2003 ) ] . j. p. home , m. j. mcdonnell , d. j. szwer , b. c. keitch , d. m. lucas , d. n. stacey , and a. m. steane , http://pra.aps.org/abstract/pra/v79/i5/e050305[phys . a * 79 , * 050305(r ) ( 2009 ) ] ; j. p. home , d. hanneke , j. d. jost , j. m. amini , d. leibfried , and d. j. wineland , http://www.sciencemag.org/content/325/5945/1227[science * 325 , * 1227 ( 2009 ) ] . a. imamoglu , http://pra.aps.org/abstract/pra/v50/i5/p3650_1[phys . a * 50 , * 3650 ( 1994 ) ] ; b. m. garraway , http://pra.aps.org/abstract/pra/v55/i3/p2290_1[phys . a * 55 , * 2290 ( 1997 ) . ] chin , a. rivas , s.f . huelga and m.b . * 51 * , 092109 ( 2010 ) ] . t. holstein , http://www.sciencedirect.com/science?_ob=articlelisturl&_method=list&_articlelistid=-67535988&_sort=r&_st=13&view=c&_acc t = c000057787&_version=1&_urlversion=0&_userid=2532525&md5=e2f1d3979b0d8f1c0b747b3c2c93f84e&searchtype = a[ann . * 8 , * 325 ( 1959 ) ] ; ( ibid ) http://www.sciencedirect.com/science/article/pii/000349165990003x[*8 , * 343 ( 1959 ) ]    d. j. wineland , c. monroe , w. m. itano , d. leibfried , b. e. king , and d. m. meekhof , http://nvlpubs.nist.gov/nistpubs/jres/103/3/cnt103-3.htm[j . res . * 103 , * 259 ( 1998 ) ] , _ and references therein_.        j. dalibard and c. cohen - tannoudji , http://iopscience.iop.org/0022-3700/18/8/019[j . b : at . mol . * 18 , * 166 ( 1985 ) ] ; y. castin , h. wallis , and j. dalibard , http://www.opticsinfobase.org/josab/abstract.cfm?uri=josab-6-11-2046[j . b. * 6 , * 2046 ( 1989 ) . ]          w. magnus , http://onlinelibrary.wiley.com/doi/10.1002/cpa.3160070404/abstract[commun . pure appl . * 7 , * 649 ( 1954 ) ] ; s. blanes , f. casas , j. a. oteo and j. ros , http://iopscience.iop.org/0143-0807/31/4/020[eur . * 31 , * 907 ( 2010 ) . ] m. pons , v , ahufinger , ch . wunderlich , a. sanpera , s. braungardt , a. sen , u. sen , and m. lewenstein , http://prl.aps.org/abstract/prl/v98/i2/e023003[phys . * 98 , * 023003 ( 2007 ) ] ; a. bermudez , d. porras , and m. a. martin - delgado , http://pra.aps.org/abstract/pra/v79/i6/e060303[phys . a * 79 , * 060303 ( r ) ( 2009 ) ] ; p. hauke , f. m. cucchietti , a. mller - hermes , m.c . bauls , j. i. cirac , and m. lewenstein , http://iopscience.iop.org/1367-2630/12/11/113037/[new j. phys . * 12 , * 113037 ( 2010 ) ] ; p. a. ivanov and f. schmidt - kaler , http://iopscience.iop.org/1367-2630/13/12/125008[new j. phys . * 13 , * 125008 ( 2011 ) ] ; a. bermudez , j. almeida , f. schmidt - kaler , a. retzker , and m. b. plenio , http://prl.aps.org/abstract/prl/v107/i20/e207209[phys . * 107 , * 207209 ( 2011 ) ] ; r. schmied , j. h. wesenberg , and d. leibfried , http://iopscience.iop.org/1367-2630/13/11/115011/[new j. phys . * 13 , * 115011 ( 2011 ) ] ; g .- d . lin , c. monroe , and l .- m . duan , http://prl.aps.org/abstract/prl/v106/i23/e230402[phys . * 106 , * 230402 ( 2011 ) ] ; a. bermudez , j. almeida , k. ott , h. kaufmann , s. ulm , u. poschinger , f. schmidt - kaler , a. retzker , and m. b. plenio , http://iopscience.iop.org/1367-2630/14/9/093042/[new j. phys . * 14 , * 093042 ( 2012 ) . ] d. kielpinski , b. e. king , c. j. myatt , c. a. sackett , q. a. turchette , w. m. itano , c. monroe , d. j. wineland , and w. h. zurek , http://pra.aps.org/abstract/pra/v61/i3/e032310[phys . rev . a * 61 , * 032310 ( 2000 ) ] ; g. morigi and h. walther , http://epjd.epj.org/index.php?option=com_article&access=standard&itemid=129&url=/articles/epjd/abs/2001/02/d0131/d0131.h tml[eur . j. d * 13 , * 261 ( 2001 ) ] . g. schn and a.d . zaikin , http://www.sciencedirect.com/science/article/pii/037015739090156v[phys . rep . * 198 * , 237 ( 1990 ) ] . l. capriotti , a. cuccoli , a. fubini , v. tognetti , and r. vaia , http://prl.aps.org/abstract/prl/v94/i15/e157001[phys . * 94 * , 157001 ( 2005 ) ] . m. a. cazalilla , f. sols , and f. guinea , http://prl.aps.org/abstract/prl/v97/i7/e076401[phys . * 97 * , 076401 ( 2006 ) ] . f. verstraete , m. m. wolf , and j. i. cirac , http://www.nature.com/nphys/journal/v5/n9/abs/nphys1342.html[nature phys . * 5 * , 633 ( 2009 ) ] . j. eisert and t. prosen , e - print http://arxiv.org/abs/1012.5013[arxiv:1012.5013 ( 2012 ) ] . s. morrison , a. s. parkins , http://prl.aps.org/abstract/prl/v100/i4/e040403[phys . lett . * 100 * , 040403 ( 2008 ) ] . t. prosen and e. ilievski , http://prl.aps.org/abstract/prl/v107/i6/e060403[phys . * 107 * , 060403 ( 2011 ) ] ; m. j. bhaseen , j. mayoh , b. d. simons , and j. keeling , http://pra.aps.org/abstract/pra/v85/i1/e013817[phys . a * 85 * , 013817 ( 2012 ) ] ; v. m. bastidas , c. emary , g. schaller , and t. brandes , http://pra.aps.org/abstract/pra/v86/i6/e063627[phys . a * 86 * , 063627 ( 2012 ) . ]
Abstract:

Answer: we propose a dissipative method to prepare the ground state of the isotropic xy spin hamiltonian in a transverse field . our model consists of a spin chain with nearest - neighbour interactions and an additional collective coupling of the spins to a damped harmonic oscillator . the latter provides an effective environment with a lorentzian spectral density and can be used to drive the chain asymptotically towards its multipartite - entangled ground state at a rate that depends on the degree of non - markovianity of the evolution . we also present a detailed proposal for the experimental implementation with a chain of trapped ions . the protocol does not require individual addressing , concatenated pulses , or multi - particle jump operators , and is capable of generating the desired target state in small ion chains with very high fidelities .