Multiple bosons undergoing coherent evolution in a coupled network of sites constitute a so-called quantum walk system. The simplest example of such a two-particle interference is the celebrated Hong-Ou-Mandel interference. When scaling to larger boson numbers, simulating the exact distribution of bosons has been shown, under reasonable assumptions, to be exponentially hard. We analyze the feasibility and expected performance of a globally connected superconducting resonator based quantum walk system, using the known characteristics of state-of-the-art components. We simulate the sensitivity of such a system to decay processes and to perturbations and compare with coherent input states.Superconducting Josephson devices are a remarkable quantum information processing platform, with single and two qubit coherences close to and even surpassing fault tolerant thresholds [1, 2]. In harmonic superconducting resonators, various non-classical states have been formed on-demand and complex entangled states between such resonators have also been demonstrated, including NOON states of high order [3][4][5]. It is very important to benchmark the quality of entanglement achieved and to steadily expand the size of the systems under study. As such entanglement grows larger, it becomes a resource for potential sensing applications [6,7] and, more fundamentally, challenges various models of spontaneous collapse [8,9] or correlated error [10].A relatively simple (experimentally) platform for such explorations is the recently reformulated problem of quantum walks and associated BosonSampling [11]. Indistinguishable bosons are placed in a coupled array of resonators and allowed to interfere via the noninteracting coherent quantum walk of the particles among the resonators. Assuming a closed system (without gauge fields), the boson network then evolves in time with the Hamiltonianwhere the summation is over all nodes in the graph, J ij is the coupling strength between node i and j, andâ i is the ladder operator for resonator i. Terms for which i = j can be included to account for varying resonator oscillation energies [12].The resulting distribution of occupation probabilities in the different resonators is (thought to be) both exponentially hard to compute and verify classically [11]. This means the problem belongs to the complexity class #P, which is considered larger and more difficult than the notorious NP class. Note, however, that BosonSampling is expected to be non-universal in the computational sense. Universality with a multi-particle quantum walk hardware was only recently established in the presence of interactions between the bosons [13].Although initially this form of quantum simulation was thought to have no practical applications, it has recently been shown to be capable of simulating the vibronic spectra of molecules, if a non-trivial initial state can be created [14], in contrast, however, with the simple single photon inputs of BosonSampling.Experimental implementations of such multi-boson interference experiments ha...