Synthetic spin-orbit coupling and topological polaritons in Janeys-Cummings lattices

The interaction between a photon and a qubit in the Janeys-Cummings (JC) model generates a kind of quasiparticle called polariton. While they are widely used in quantum optics, difficulties in engineering controllable coupling of them severely limit their applications to simulate spinful quantum systems. Here we show that, in the superconducting quantum circuit context, polariton states in the single-excitation manifold of a JC lattice can be used to simulate a spin-1/2 system, based on which tunable synthetic spin-orbit coupling and novel topological polaritons can be generated and explored. The lattice is formed by a sequence of coupled transmission line resonators, each of which is connected to a transmon qubit. Synthetic spin-orbit coupling and effective Zeeman field of the polariton can both be tuned by modulating the coupling strength between neighbouring resonators, allowing for the realization of a large variety of polaritonic topological semimetal bands. Methods for detecting the polaritonic topological edge states and topological invariants are also proposed. Therefore, our work suggests that the JC lattice is a versatile platform for exploring spinful topological states of matter, which may inspire developments of topologically protected quantum optical and information processing devices.


INTRODUCTION
The Janeys-Cummings (JC) model proposed in 1963 [1] is a seminal theoretical model treating light-matter interaction with full quantum theory, i.e., the interaction of a quantized electromagnetic field with a two-level atom. This model has been widely applied to many quantum platforms for studying the interaction of a quantized bosonic field with a qubit, which now has become the cornerstone in quantum optics and quantum computation [2][3][4][5][6][7]. Furthermore, an interconnected array of multiple JC systems can form a JC lattice [8][9][10][11], which provides an innovative quantum optical platform for studying condensed matter physics. This is highlighted by previous works which show that coupled JC systems can be used to realize the Bose-Hubbard model and investigate superfluidto-Mott-insulator phase transition [12][13][14]. However, spinful lattices have not been simulated in this platform due to the difficulty in engineering a tunable coupling between different cavities.
On the other hand, the search of topological states of matter in artificial systems recently has become a rapidly growing field of research [15][16][17][18][19][20][21][22][23][24]. Topological states are characterized by topological invariants which are robust to the smooth changes in system parameters and disorders, where topological edge states can be employed for robust quantum transport [25,26]. Therefore, they hold tremendous promise for fundamental new states of matter as well as for dissipationless quantum transport devices and topological quantum computation [27]. One of the key ingredients for generating such states is to realize tunable spin-orbit coupling (SOC). Significant theoretical and experimental progress on realizing synthetic SOC recently have been achieved in ultracold atom systems [28][29][30]. This progress stimulate great research interests to explore topological states with ultracold atoms trapped in optical lattices [31][32][33][34]. However, the limited trapping time and the site addressing difficulty increase the experimental complexity, i.e., it is generally difficult to have modulable coupling between two neighbouring sites in an optical lattice.
Here, we find that the JC lattice system can be used to realize various topological spin lattice models, where synthetic polaritonic SOC and Zeeman field can be induced with in situ tunability, which provides a flexible platform to explore topological states of matter with great controllability. Specifically, we consider realizing the JC lattice in the context of superconducting quantum circuits, where each JC lattice site is constructed by a transmission line resonator (TLR) coupled to a two-level transmon qubit. We find that the dressed polariton states in the single-excitation manifold in each JC lattice site can simulate a spin-1/2 system. Particularly, synthetic SOC and Zeeman field for polaritons can be induced and manipulated by only engineering the coupling strength between neighbouring resonators. Meanwhile, we show that, based on tunable synthetic SOC and Zeeman field, nodal-loop semimetal bands [35][36][37] and topological polaritons can be realized and explored in the simulated JC lattice. Moreover, through calculating the topological winding number, we find that this tunable system has a rich topological phase diagram.
Our proposal to explore the topological states in the JC lattice system is different from previous ones based on the optical lattices [28][29][30]. In particular, our proposal has a number of advantages. (i) Unlike ultracold atoms, polaritons are quasiparticles which are hybrids of photons and qubit excitations. Topological polaritons emerge from the topological structure of light-matter interaction, where photons and qubit excitations are topologically trivial by themselves, but combining together, they become hybrid topological states. Therefore using polariton for quantum simulation enriches our controlling methods -both photonic and atomic means take effects. (ii) The systematic parameters in JC lattice systems can be tuned at a single-site level, which allows us to generate a wide variety of SOC forms. (iii) The geometry of JC lattice can be artificially designed and lattice boundaries are easy to be created for observing topological edge states, thus various topological lattice models and topological effects can be constructed and probed. (iv) The particle number putting in a JC lattice can be deterministically controlled. With such an advantage, we present a method using single-particle quantum dynamics to probe topological winding numbers and topological polariton edge states. (v) In the quantum optics platform, JC lattice systems previously have generated multiple important applications, including masers, lasers, photon transistors, and quantum information processors. Meanwhile, there are indeed several disadvantages in our proposed JC lattice system, such as the limited system size, parameter fluctuations, and decoherence. However, the essential physics of the simulated topological polariton states, such as the topological invariants and edge states, can still be detected under these realistic circumstances. Therefore, the topological JC lattice system in superconducting quantum circuits offers the possibility to develop functional topological spin quantum devices.

Janeys-Cummings lattice
The method for implementing a one-dimensional (1D) JC lattice in superconducting quantum circuits [4] is as follow. As shown in Fig. 1A, every unit cell consists of a TLR resonantly coupled with a transmon, forming a JC model [38]. The neighbouring TLRs are connected by a combination of a SQUID and a small inductor L in series, which can actually be regarded as the counterpart of a semitransparent mirror in the cavity QED system, allowing photons to hop across neighbouring cavities (see Methods). As a result, setting = 1 hereafter, the system Hamiltonian of this JC lattice is where N is the number of the unit cells; l is the JC type interacting Hamiltonian in lth unit cell. The condition g l ω l has to be met for justifying the JC coupling. σ + l = |e g| and σ − l = |g e| are the raising and lowering operators of the lth transmon qubits.â l andâ † l are the annihilation and creation operators of the photon in lth TLR. And J l (t) is the inter-TLR hopping strengths between lth and (l + 1)th unit cells. Different from optical cavities, the time-dependence of J l (t) here can be induced by adding a time-varying external magnetic flux threading through the SQUIDs (see Methods).
is the counter-rotating term that can be neglected by using the rotating wave approximation (RWA). The lowest three eigenstates of the JC Hamiltonian h l for the lth unit are |0g l , |↑ l = (|0e l + |1g l ) / √ 2 and |↓ l = (|0e l − |1g l ) / √ 2, where |ng l and |ne l (n = 0, 1, 2, · · · ) are the states containing n photons while the transmon is at the ground and excited state, respectively. And their eigenenergies are E l,0g = 0 and E l,↑(↓) = ω l + (−)g l . Here, we exploit the two single-excitation eigenstates | ↑ l and | ↓ l to simulate the effective electronic spin-up and spin-down state. As each of the two states consists of half "photon" and half "atom", they are regarded as a whole and was termed as "polariton".

Polaritonic spin-orbit coupling
We proceed to show that a spin-1/2 chain model with tunable Zeeman field and SOC can be simulated with the JC lattice by only adjusting the pulse shape of the coupling strengths J l (t) between neighbouring TLRs. Firstly, since each cell contains two pseudo-spin states, there are totally four intercell neighbouring hopping. In order to control each hopping separately, selective frequency addressing is employed (see Methods), i.e., we assign each of the four hopping with its unique hopping frequency. To achieve this, we adopt two sets of unit cells, A-type and B-type, which are different in the sense that they have different eigen-frequencies and coupling strengths of the JC model. Then we arrange them in an alternate way, as shown in Fig. 1A. Setting started with an A-type one, when l is odd (even), ω l = ω A (ω B ) and g l = g A (g B ).
Then, based on the current experimental reaches [4], we set ω A /2π = 6 GHz, ω B /2π = 5.65 GHz, g A /2π = 300 MHz, and g B /2π = 270 MHz. With this, the energy intervals of the four hopping are {|E l,α − E l+1,α |/2π} = {220, 320, 380, 920} MHz with α, α ∈ {↑, ↓}. The differences between every two of them are no less than 20 times of the effective hopping strength t 0 /2π = 3 MHz, thus they can be selective addressed in frequency. Correspondingly, the J l (t) contains four tunes and can be written as where is the sign of the hopping phase, 4t 0,lαα , ω d lαα and s lαα ϕ lαα are the amplitudes, frequencies and phases, corresponding to the hopping |α l → |α l+1 , respectively. Note that due to the alternate arrangement, based on the definition above, the J l (t)s take only two different forms -when l is odd (even), J l = J AB (J BA ) where the subscript "AB" ("BA") refers to that the A-type (B-type) unit cell is on the left. Experimentally, this time-dependent coupling strength J l (t) can be realized by adding external magnetic fluxes with both dc and ac components threading through the SQUIDs (see Methods). In this way, selective hopping can be induced, i.e., only when a frequency of the driving flux matches a particular hopping energy interval, that hopping can be triggered, otherwise it will not take into effect. Meanwhile, both the strengths and phases of the hopping can be controlled by the amplitudes and phases of the ac magnetic fluxes. However, only controlling these two is not enough for realizing the topological states, we still need to induce and adjust the spin splitting. Then we theoretically find out, and numerically prove that this spin splitting can be induced by just adding a detuning to the spinflipped transition tunes. Concretely, while the spin-preserved transition frequencies are set as ω d lαα = |E l,α − E l+1,α |, we set the spin-flipped transition frequencies with a detuning 2m Fig. 1B. Thus, in the rotating frame (explained later), there will be a spin splitting m for each cell, as shown in Fig. 1C. We now show how the time-dependent coupling strength in Eq. (2) can induce a designable spin transition process in a certain rotating frame. First, we map the Hamiltonian in Eq. (1) into the single excitation direct product subspace span{|0g, · · · , 0g, α lth , 0g, · · · , 0g } l=1,··· ,N ; α=↑,↓ . Hereafter, when there is no ambiguity, we use |α l to denote |0g, · · · , 0g, α lth , 0g, · · · , 0g , and |G to denote |0g, · · · , 0g .
Then, we define a rotating frame by a unitary operator Neglecting the fast rotating terms (see Methods), one obtains where S z l = |↑ l ↑| − |↓ l ↓|, t 0,lαα is the effective coupling strength andĉ † l,α = |α l G| is the creation operator for polariton with "spin" α in lth unit cell. As a result, the Hamiltonian (4) represents a general 1D spin-1/2 tight-binding lattice model with m, t 0,lαα and ϕ lαα being the equivalent Zeeman energy, hopping strength and hopping phase, respectively. These three variables can all be experimentally tuned in wide ranges by the frequencies, amplitudes and phases of the external ac magnetic fluxes. Notably, the effective on-site potential m can be tuned as either positive or negative depending on the detuning direction.
Meanwhile, it is worth noticing that although we have introduced two kinds of unit cells (A-type and B-type) in the laboratory frame, by adjusting the coupling parameters J AB and J BA , the translation symmetry in the rotating frame is still preserved. In other words, the smallest repeating unit in the rotating frame contains only one unit cell as shown in Fig. 1C. By now, the adjustable spin-preserved tunneling (α = α ) and SOC terms (α = α ) can both be induced, hence our superconducting quantum circuit setup can naturally be used to simulate a tunable SOC topological polariton insulator.
In order to justify the individual frequency addressing of the inter-cell transitions, we numerically simulate the dynamics of a system containing only two unit cells, with an A-type one on the left and a B-type one on the right. We test every hopping of the four transitions {|α A ↔ |α B } one by one, adding only one corresponding frequency ω d 1αα in J AB (t). As expected, when we pick out a resonant frequency ω d 1αα with m = 0, there will be a Rabi oscillation between the two corresponding target states. One example of these Rabi oscillations was shown in Fig. 1D, which is of the least fidelity among the four. Even in this worst case, and in the third Rabi cycle, the fidelity still reaches a high value of 0.9979, which justify the RWA. In addition, our numerical simulation also shows that, in the present of the unmatched driving, all the initial non-target states remain almost unchanged, thus justify our individual frequency addressing method.

Nodal-loop topological polaritons
Our protocol provides a tunable platform using polaritons to study topological matters. Here, we take the nodal-loop semimetal as an application sample to demonstrate how to simulate a specific condensed matter model in our proposed setup. To fit our simulation setup, we reform the Hamiltonian of the original 3D nodal-loop model in Ref. [36], without losing any topological properties. Firstly, we relabel the coordinates to set the hopping terms with SOC to be along the x axis. Secondly, we consider the Fourier transformations along y and z directions with quasi-momenta k y and k z and treat them as system parameters. Then, according to Eq. (4), we can simulate this 3D nodal-loop model in our 1D system with the other two dimensions being the parametric dimensions. In this direct simulation, to engineer the four transitions between different paloriton states, we need four different tunes in the inter-cell coupling strength J nod l (t). The above implementation can further be simplified as following. We first make a unitary transformation to the original Hamiltonian so that The transformed 1D lattice Hamiltonian from Eq. (4), without losing any physical properties, is where m (k y , k z ) = M + 2d(cos k y + cos k z ) with M being the effective Zeeman energy and d being the effective hopping energy along y and z directions. Then, we set the parameters of the JC model to be ω A = ω B = 2π × 6 GHz, g A /2π = 200 MHz and g B /2π = 100 MHz. Correspondingly, the coupling strengths can be set to contain only two tunes as where ω d 1 /2π = 100 MHz and ω d 2 /2π = (300 − 2m) MHz. In this way, when transforming into the rotating frame of U and applying the RWA, the Hamiltonian (1) of our JC-lattice system will takes the form of Eq. (5), which accomplishes the quantum simulation of the topological nodal-loop semimetals (see Methods).
Experimentally, one can choose t 0 /2π = 3 MHz and set the detuning within the range of m ∈ 2π × [−20, 20] MHz, such that any value of the variable tunes ω d lαα (α = α ) still maintain a frequency difference no less than 20t 0 away from other tunes, and thus the crosstalk caused by unwanted tunes will be negligible. Last, for testing the validity of our theoretical protocol, numerical simulations will be given in the Experimental detection section after a brief introduction of the characteristics of the topological nodal-loop polaritons.

Characteristics of the topological polaritons
To investigate the bulk characteristic of the topological nodal-loop polaritons, we first consider the periodic boundary condition to obtain the Hamiltonian in Eq. (5) in the momentum space as where b y = −2t 0 cos k x , b z = 2t 0 sin k x + m (k y , k z ) and S y = i(| ↓ ↑ | − | ↑ ↓ |). This system has two energy bands which will touch when E = 0. The touching points form closed lines, the so-called nodal-loops, in momentum space as shown in Fig. 2A. These two loops appear in the k x = π/2 and/or k x = −π/2 planes. By fixing k x = ±π/2, we plot these two energy bands in the k y -k z space as shown in Fig. 2B and 2C, where one can see that the touching is right along the nodal-loops. The topological index characterizing each nodal loop is a winding number defined as where This shows that the quantized winding number is either 1 or 0, corresponding to the cases whether a nodal loop is enclosing the straight line, which is along k x direction of the fixed (k x , k y ) point, or not [36,37]. Hence, the nodal-loops divide the k y -k z space into regions with different winding numbers, as shown in Fig. 2D. For all straight lines along k x inside the nodal loop, each of them can be regarded as being corresponding to a topological 1D gapped subsystem with winding number 1.
Two striking topological characters of nodal-loop semimetal are the zero-energy modes inside the energy gap and their corresponding edge states. We take a slice of k y = 0, indicated by the green line in Fig. 2D, as an example to plot the energy spectrum with various k z for a finite chain with N = 20, as shown in Fig. 2E. The numerical result shows that there are two mid-gap degenerated zero-energy modes appear in the range of ν = 1 (the red area in Fig. 2D). The quantum states corresponding to the two mid-gap energies are edge states localized in the left and right end of the lattice, respectively. When N is large enough for ignoring the finite-size effects, their wave-functions can be expressed as where A = (1 − q 2 )/2(1 − q 2N ) and λ = ln(1/q) with q = tan(m π/8t 0 ). The phase factor i l−1 stems from the unitary transformation V to get the nodal-loop Hamiltonian in Eq. (5) that simplified our simulation, from the original Hamiltonian in Eq. (4). These analytical results of wavefunctions are in very good agreement with the numerical results for N = 20, as shown in Fig. 2F, thus justifies that we can use the JC lattice of experimentally capable size to simulate the topological features.
There are several phases in our simulated Hamiltonian where the phase transition is indicated by the emerging or vanishing of the nodal loops. Inferring from Eq. (8), the critical conditions are obtained as k x = −π/2 and k x = π/2, which are corresponding to one nodal loop in each of the two regions of −2t 0 −4d < M < −2t 0 +4d and 2t 0 −4d < M < 2t 0 +4d with d > 0 in the M -d plane. Therefore, in the area where the two regions overlap, there are two nodal loops. But in the area outside these two regions, there is no nodal loop so that the whole Brillouin zone will be in a purely trivial or nontrivial phase. Consequently, there are totally five different phases of different winding number configurations ν(k y , k z ) in the M -d plane, as shown in Fig. 3. According to the chosen parameters, t 0 and m, the area {M, d| − 6t 0 + 4d < M < 6t 0 − 4d} will include all the the five phases in the phase diagram of Fig. 3, i.e., all the possible phases can all be simulated in our proposed system.

Experimental detection methods
Polaritonic topological edge state detection. According to Eq. (10), or as shown in Fig. 2F, the polariton in the left or right edge state is maximally distributed in the leftmost or rightmost JC lattice site. Their internal spins are in the superposition states (| ↑ l + | ↓ l ) / √ 2 and (| ↑ l − | ↓ l ) / √ 2, respectively. Therefore one can find that, in the beginning, the left and right polaritonic edge states only have qubit excitation and photon components, respectively. Taking the detection of left edge state as an example, initially, the polariton in the leftmost JC lattice site is prepared into |0e 1 , i.e., the leftmost qubit (resonator) has been prepared in the excited (vacuum) state. The qubits and resonators in the other sites are prepared into the ground and vacuum states which means the initial systematic state is |ψ(t = 0) = |0e 1 |0g 2 · · · |0g N . After that, we let the above initial state evolve for a time about 0.5 µs. If the JC lattice is in the topological nontrivial phase supporting the left edge states, the final density distribution of the polaritons will maximally populate the leftmost site. The reason is that the initial state |ψ(t = 0) has a large overlap with the left edge state. It will evolve mainly via the edge state wave packet and maximally localized in the leftmost site. While if the system is in the topological trivial phase and has no edge states, the initial state will be a superposition of different bulk sates. The final density distribution will not have maximal dis- tribution in the leftmost site. Similarly, one also can prepare the JC lattice into |ψ(t = 0) = |0g 1 · · · |0g N −1 |1g N and detect [39] the right polaritonic topological edge state based on observing its time evolution.
In Fig. 4A and 4B, we have numerically calculated the time evolution of the polaritonic density when the JC lattice is in the topological trivial and nontrivial nodal-loop semimetal phases, respectively. For the trivial case, the wave packet has a ballistic spread versus time, which is a typical feature of bulk Bloch state. It shows that there is no edge state localization and the system is in topological trivial phase. For the nontrivial case, the density of the polaritons will always maximally localize in the leftmost JC lattice site, which indicates the existence of left topological edge states demonstrating that the system is in topological nontrivial phase. The time evolution of the qubit excitation and the photon population for the topological nontrivial case are also numerically calculated in Fig. 4C and 4D, which shows that the localized qubit excitation and the photon in the leftmost site have a Rabi-like oscillation feature inherited from the JC model.
Polaritonic topological invariant detection. Another important hallmark for topological nontrivial nodal-loop polaritonic semimetal phase is the nontrivial polaritonic topological winding number. Here we show that such polaritonic topological invariant also can be dynamically detected. Our method is based on a previous work which shows that the topological winding number rotted in the momentum space can be detected through measuring the dynamical chiral center in the real space [40]. The chiral operator for our topological polaritonic model is S x l = |↑ l ↓| + |↓ l ↑|. Then the chiral center operator for the JC lattice is defined asP d = polaritonic topological winding number can be related with the time-averaged dynamical chiral center associated with the single-polariton dynamics, i.e., where T is the evolution time, |ψ c (t) = exp(−iH nod t)|ψ c (0) is the time evolution of the initial single-polariton state |ψ c (0) = |0g 1 · · · | ↑ N/2 · · · |0g N , where one of the middle JC lattice site has been put one polariton in, with its spin prepared in the state | ↑ .
In Fig. 5A and 5B, we have numerically calculated the dynamical chiral centerP d (t) = ψ c (t)|P d |ψ c (t) for topological trivial and nontrivial cases, respectively. According to Eq. (11), one can find that the topological winding number is equal to twice the oscillation center ofP d (t). As shown in Fig. 5A,P d (t) oscillates around the average value 0, which gives the polaritonic topological winding number ν = 0. In contrast, the result for the topological nontrivial case in Fig. 5B shows thatP d (t) oscillates around 0.5, which yields the polaritonic topological winding number ν = 1. Experimentally, the states of qubits and resonators can be measured with fidelity higher than 0.99 in superconducting circuits. In our case, one only need to measure the qubit excitation and photon populations for getting their imbalance and deriving the chiral center, without requiring full quantum state tomography. In this way, the topological winding number can be easily and unambiguously detected based on monitoring singlepolariton quantum dynamics in a JC lattice.

DISCUSSION
We further investigate how our construction and detection methods are influenced by the quantum decoherence effects. The Lindblad master equation is adopted to take three main decoherence factors including the losses of the photon and the decay and dephasing of the transmon into account. The Lindblad master equation of our system can be written aṡ where ρ is the density operator of the whole system, γ is the decay rate or noise strength which are set to be the same here, Γ l,1 = a l , Γ l,2 = σ − l and Γ l,3 = σ z l are the photon-loss, transmon-loss and the transmon-dephasing operators in the lth lattice, respectively. In Fig. 6A, we plot the edge-site population P 1 (t) = tr ρ(t) a † 1 a 1 + σ + 1 σ − 1 after 0.5 µs and the oscillation center ν/2 of the trivial and nontrivial cases for different decay rates. It shows that the edge state population and the chiral center smoothly decrease when the decay rate increase. However, our detection method can tolerate the decay rate up to the order of 2π × 100 kHz, while the typical decay rates are only 2π × 5 kHz. We also investigate the minimum sites that are needed for the experimental detecting the oscillation-center. In the presence of the decay rates of 2π × 5 kHz, we plot the oscillation center for lattice with different number of sites in both the trivial and nontrivial cases, as shown in Fig. 6B. It shows that the four sites case corresponding to an oscillation center value about 0.40, is big enough for distinguishing the topologically trivial and nontrivial cases with current state-of-the-art technologies.
To conclude, we have introduced the concept of topological states into the JC lattice which is one of the most important building blocks in quantum optics and quantum information processing. We have studied the topological structure of light-matter interaction and shown that SOC physics and topological polaritons can be pursued in the JC lattice. Different from synthetic topological states in ultracold atomic, photonic and acoustic systems, topological polaritons are topological superposition states of photons and qubits. Tunable synthetic polaritonic SOC is induced by engineering the JC lattice couplings, which provides the basic ingredient for realizing spinful topological states of matter. We have also provided a method using single-particle quantum dynamics in real space to directly observe the polaritonic topological edge states and topological invariants.
Our work has a broad generalization and opens the door for exploring spinful topological states of matter using polaritons in JC lattice system. (i) In addition to mimicking spin-1/2, polariton states in multiple-excitation manifold have an extended spin degree of freedom and can also be used to mimic larger spin. It is challenging to realize SOC for the larger spin case in solid state materials and ultracold atoms. However, the method proposed in our work can be generalized to realize synthetic SOC for large-spin polaritons. This allows us to explore a large variety of topological states, including triple point topological states of matter [41]; (ii) Polariton states, the eigenstates of a JC model, can be referred as a synthetic dimension where each pair of states mimic a spin in a spinlattice model and the coupling between polariton states provides the hopping in the synthetic dimension. With such synthetic dimension, high-dimensional topological states of matter can be explored in a low-dimensional JC lattice, including topological states beyond three dimensions [42,43]; (iii) Polaritons have tunable strong nonlinear interaction, which allows us to study polaritonic fractional topological states of matter [44]; (iv) Besides superconducting circuit, JC lattices can also be realized in many quantum optical systems, including trapped ions [2], cavity quantum electrodynamics [3], nanoscopic lattice [6], optomechanical systems [7] and so on.

SQUID induced time-dependent photon hopping
We here present how to induce the time-dependent photon hopping strength between two TLRs. As shown in Fig. 1A, neighbouring TLRs are connected by a common SQUID and an inductor L in series then to the ground. The SQUID actually serves as a single Josephson junction (JJ) but with effective Josephson inductance tunable by the external flux. Concretely, by applying external magnetic flux Φ ext = Φ dc + Φ ac threading through the SQUID, when Φ ac φ 0 with φ 0 being the reduced flux quanta, the effective inductance of the SQUID reads [45] where I c is the shared critical current of the two JJ in each SQUID. On the other hand, comparing to the inductance of the TLRs, both the SQUID and the inductor L have far smaller inductances, thus there is a voltage node but a current peak at both ends of each TLR. For these boundary conditions, after the conventional quantization of the TLRs [38,45], the flux density and the charge density wave-function of the lowestenergy mode in lth TLR can be expressed aŝ where ω l = π/ √ L l C l is the frequency of the photon with L l and C l being the inductance and capacitance, respectively. d l is the length and x l is the coordinate of the lth TLR. Meanwhile, because of the relatively low-impedances of the SQUID and the inductor, the currents from the ends of every TLR will flow directly through them to the ground, without crossing to their neighbouring TLRs. Hence, the interaction Hamiltonian between the lth and (l + 1)th TLR is just the summation energy of the SQUID and the inductor L where I ri(le) l =φ l (x l , t)d l /L l | x l =d l (0) is the current of right(left)-end of the lth TLR. Moreover, If we set where n is the number of the tunes in Φ ac , Φ dc = 4πn φ 0 where n is an arbitrary positive integer, and choose resonant or detuned frequencies of ω j , after the RWA, we obtain which right meets the form of Eqs. (1) and (2). This equation can be interpreted as describing the photons hopping between neighbouring unit cells, which means that the SQUID-L combination can actually serve as a counterpart of the semitransparent mirror in the optical cavity system.
whereĉ † l,α = |α l G|. Then in every J l (t), we add four tunes, corresponding to the four inter-cell hopping, cf. Eq. (17), each of which contains its independent tunable amplitude, frequency and phase as J l (t) = α,α 4t 0,lαα cos ω d lαα t + s lαα ϕ lαα , where s lαα is defined in Eq. (3) being the sign of each phase and the frequencies are set as where 2m is a detuning. By now, the form of J l (t) is determined leaving m, t 0,lαα and ϕ lαα to be chosen arbitrarily depending on the topological insulator model one simulates. Eventually, the target topological insulator model is in the rotating frame transformed by +e i(E l+1,↓ −E l,↑ +2m)tĉ l,↑ĉ † l+1,↓ + h.c., (22) where S z l = |↑ l ↑| − |↓ l ↓|. After doing the multiplication in this equation, the four resonant terms and their Hermitian conjugates will be absent of time t (frequency addressing). Meanwhile, if the conditions {t 0,lαα } α,α =↑,↓ {ω lαα , ω lαα − ω lββ } α,α ,β,β =↑,↓; ω lαα =ω lββ are satisfied, all the other terms are fast rotating term that can be dropped within RWA. As a result, we obtain the tight-binding model with tunable SOC in Eq. (4).

Simplified method of implementation
We here explain why the four hopping terms in the nodalloop semimetal modal can be induced by the coupling strength in Eq. (6) using only two tunes. Firstly, when we set parameters of the JC model to be ω A = ω B = 2π × 6 GHz, g A /2π = 200 MHz and g B /2π = 100 MHz, the energy intervals of the four hopping overlap and reduce into two set of intervals, i.e., a spin-conserved hopping interval |E l,↑ − E l+1,↑ | = |E l,↓ − E l+1,↓ | = 2π × 100MHz, (23) and a spin-flipped hopping interval |E l,↑ − E l+1,↓ | = |E l,↓ − E l+1,↑ | = 2π × 300MHz. (24) Therefore, in this way, one tune in J nod l (t) can induce two hopping while the two spin-conserved hopping still remain being controlled separated from the two spin-flipped hopping, thus they can be of different detuning.
Secondly, note that within this setting, according to the definition in Eq. (3), we have s l↑↑ = −s l↓↓ and s l↑↓ = −s l↓↑ , thus both the two spin-conserved hopping terms and spinflipped hopping terms can be induced by a same tune, but with opposite signs, i.e., i and −i. This is possible because the RWA selects different terms. These opposite signs are ideal for realizing the wanted SOC in our protocol. Anyway, this is only possible when we do the unitary transformation V to transform the hopping phase from the original 1 and −1 into the pure imaginary numbers i and −i. For the former case, one still has to use four tunes and induce them separately. Therefore, this transformation simplifies our simulation.