By using extended bosonic coherent states, the solution to the Jaynes-Cummings model without the rotating-wave approximation can be mapped to that of a polynomial equation with a single variable. The solutions to this polynomial equation can give all eigenvalues and eigenfunctions of this model with all values of the coupling strength and the detuning exactly, which can be readily applied to recent circuit quantum electrodynamic systems operating in the ultra-strong coupling regime.PACS numbers: 42.50. Pq, 03.65.Ge, 85.25.Cp, 03.67.Lx The Jaynes-Cummings (JC) model[1] describes the interaction of a two-level atom (qubit) with a single bosonic mode. It is a fundamental one in quantum optics. Based on the assumption of nearresonance and relatively weak atom-cavity coupling, the rotating-wave approximation (RWA) is usually employed, and analytically exact solution can be trivially obtained.Recently, the JC model is closely related to condensed matter physics. It can be realized in some solid-state systems recently, such as one Josephson charge qubit coupling to an electromagnetic resonator [2], the superconducting quantum interference device coupled with a nanomechanical resonator [3,4], and the most recently LC resonator magnetically coupled to a superconducting qubit [5][6][7]. In traditional quantum optics where the coupling between the two-level "natural" atom and the single bosonic mode is quite weak, RWA is the most useful approximation. However, in the circuit quantum electrodynamic (QED), the artificial atoms may interact very strongly with on-chip resonant circuits [6][7][8][9][10], the RWA can not describe well the strong coupling regime [6]. Therefore, the JC model without the RWA is the focus of current interests [11][12][13][14][15][16][17][18][19][20][21].However, due to the inclusion of the counter-rotating terms, the Bosonic number is not conserved, the Bosonic Fock space has infinite dimensions, so any solution without the RWA is highly nontrivial. In the recent years, several non-RWA approaches have been proposed in the Dicke model [12], the quantum Zeno effect [22], and the spin-boson model [23]. Especially, by using extended bosonic coherent states, the present authors have solved the Dicke model without the RWA exactly in the numerical sense [12]. The most simple N = 1 Dicke model is just the JC model. This numerically exact solutions to the JC model are also described in detail in Refs. [18,19]. To the best of our knowledge, an analytical exact solutions are still lacking in the literature to date.In this paper, we propose a new method to solve exactly the JC model without the RWA by means of extended bosonic coherent states. The correlations among bosons are added step by step until further corrections will not change the results. Different from our previous work where the pure coherent state are fixed, the eigenvalue α of the pure coherent state in this paper is tunable. By solving Schr Without the RWA, the Hamiltonian of a qubit interacting with a single bosonic mode reads ( = 1)where a ...