We have studied the low-lying excitations of a chain of coupled circuit-QED systems in the ultrastrong coupling regime, and report several intriguing properties of its two nearly degenerate ground states. The ground states are Schrödinger cat states at a truly large scale, involving maximal entanglement between the resonators and the qubits, and are mathematically equivalent to Majorana bound states. With a suitable design of physical qubits, they are protected against local fluctuations and constitute a non-local qubit. Further, they can be probed and manipulated coherently by attaching an empty resonator to one end of the circuit-QED chain.Confronted with formidable difficulties in solving strongly interacting many-body systems, it has been desired to find good quantum simulators. It may seem natural to simulate a many-body system with another tunable system of massive particles such as ultracold atomic gases [1]. In fact, any controllable quantum system, notably quantum computer if ever practical, can simulate efficiently many-body systems [2]. Indeed it has been recognized that photons confined in coupledcavities simulate closely the quantum behaviors of stronglycorrelated many-body systems [3,4]. Subsequent studies have revealed that Bose-Hubbard model [5], interacting spin models [4,6], and other exotic quantum phases [7] can be simulated efficiently using the coupled-cavities. Further, recent advances in solid-state devices such as circuit-QED systems [8,9] and micro-cavities [10,11] and ongoing efforts to fabricate large-scale cavity arrays [12] make the array of coupled cavities a promising candidate for an efficient quantum simulator.Meanwhile, the ultrastrong coupling regime of the cavity-QED system, where the light-matter coupling energy is comparable to or even higher than the energy of the cavity field, has been envisioned [13] and experimentally demonstrated [14]. The ultrastrong coupling brings about fundamentally different physics deeply connected to the high degree of entanglement between the "matter" and the photon [15][16][17][18][19]. However, the effect of ultrastrong coupling on the lowenergy excitations of an array of coupled cavity-QED systems remains unclear, and is our main concern in this work.In this paper, we investigate the low-lying excitations of a one-dimensional (1D) array of circuit-QED systems (cQEDs), with each cQED being in the ultrastrong coupling regime; see Fig. 1. It turns out that the array permits two nearly degenerate ground states separated by a finite energy gap from the continuum of higher-energy states. We find several intriguing properties of the two ground states: (i) They are Schrödinger cat states at a truly large scale, and involve maximal entanglement between the resonators and the qubits. (ii) With a suitable design of physical qubits, the two ground states are protected against local fluctuations and constitute a non-local qubit [20]. (iii) They are mathematically equivalent to the long-searched Majorana bound states [21]. (iv) They can be probed and mani...