Generating Multipartite Entanglement via Structured-Bath Engineering

Abstract

We present how to dissipatively prepare genuinely multipartite-entangled steady states with mutually independent local reservoirs by harnessing an Abelian strong symmetry and temporal bath correlations. In interacting Ising chains, we demonstrate that, for a broad class of system-environment couplings, flat baths exponentially erase Ising-axis correlations. Structured baths overcome this obstruction by converting temporal bath correlations into classical spatial correlations. Strong symmetry then promotes these correlations to the quantum domain, enabling global multipartite entanglement. Within weak-coupling theory, we obtain exact steady states despite the generator's nonquadratic nature: strong symmetry reduces their construction to a classical problem for arbitrary Ising couplings. This reduction reveals an exact connection between classical distinguishability and entanglement. We then show how to engineer the bath spectrum to increase the correlation length and enable the preparation of GHZ-like states from separable initial conditions. We test the mechanism beyond the analytical regime using exact Markovian embeddings and tensor-network simulations. Our approach is particularly suited to superconducting and trapped-ion platforms, and we propose a concrete implementation for the latter using sympathetically cooled motional modes.

Can spatially independent, single-site reservoirs prepare genuinely multipartite-entangled steady states in interacting Ising systems? This paper studies how frequency-selective dissipation and a conserved global $\mathbb Z_2$ parity combine to do so, including preparation from separable states in a fixed parity sector. For specified symmetry-respecting local couplings and arbitrary Ising interactions, flat-spectrum Lindblad baths exponentially erase connected Ising-axis correlations. Structured baths instead bias configuration-dependent spin flips. In the secular weak-coupling description of a uniform nearest-neighbor chain, stationary populations are classical Gibbs weights with $\beta_{\rm eff}=(4J)^{-1}\ln[S(4J)/S(-4J)]$, even for nonthermal reservoirs. For parity-conserving couplings $L_j=\sigma_j^z$, population dynamics and coherences between globally reversed configurations obey the same classical master equation, without secularization at the Redfield level. This yields $\rho_{\rm ss}=\rho_{\rm cl}[1+\langle U_z(0)\rangle U_z]$ for arbitrary Ising couplings. Within a fixed parity sector, bipartite negativity equals half the total-variation distance between the stationary distribution and its partial spin reversal. Every proper reduced state is separable. Genuine multipartite entanglement occurs exactly when $\max_{\mathbf s}P_{\mathbf s}>1/4$; for nearest-neighbor periodic chains, the required correlation length approaches $0.38N$. Lorentzian-reservoir Markovian embeddings, tensor-network calculations, and small-system simulations test these predictions beyond weak coupling. A trapped-ion implementation is proposed, not demonstrated. Finite reservoir occupation limits spectral asymmetry to $(\bar n+1)/\bar n$; certification assumes known parity. Redfield steady-state completeness is generic rather than unconditional, and preparation-time and dissipative-gap scaling remain open.