achieving it. At the same time, the Bose-Hubbard model with attractive interactions (ABH) is theoretically known to have a phase transition from the superfluid phase to a highly entangled nonlocal superposition, but observation of this phase transition has remained out of experimental reach. In this theoretical work, we jointly address these two problems by (1) proposing an experimentally accessible quantum simulation of the ABH phase transition in an array of tunably coupled superconducting circuit microwave resonators and (2) incorporating the simulation into a highly scalable protocol that takes as input any microwave resonator state with negligible occupation of number states |0> and |1> and nonlocally superposes it across the whole array of resonators. The large-scale multipartite entanglement produced by the protocol is of the W-type, which is well-known for its robustness. The protocol utilizes the ABH phase transition to generate the multipartite entanglement of all of the resonators in parallel, and is therefore deterministic and permits an increase in resonator number without increase in protocol complexity; the number of resonators is limited instead by system characteristics such as resonator frequency disorder and inter-resonator coupling strength. Only one local and two global controls are required for the protocol. We numerically demonstrate the protocol with realistic system parameters, and estimate that current experimental capabilities can realize the protocol with high fidelity for greater than 40 resonators.

desirable, but not yet observed, are mechanical states that exhibit entanglement, wherein non-classical correlations exist between distinct modes; squeezing, wherein the quantum uncertainty of an observable quantity is reduced below the standard quantum limit; and Schr“odinger cats, wherein a single mode is cast in a quantum superposition of macroscopically distinct classical states. Also, while most investigations of electromechanical systems have focussed on single- or few-body scenarios, the many-body regime remains virtually unexplored. In such a regime quantum phase transitions naturally present themselves as a resource for quantum state generation, thereby providing a route toward entangling a large number of electromechanical systems in highly non-classical states. Here we show how to use existing superconducting circuit technology to implement a (quasi) quantum phase transition in an array of electromechanical systems such that entanglement, squeezing, and Schr“odinger cats become simultaneously observable across multiple microwave and micromechanical oscillators.