Chaos and two-level dynamics of the Atomtronic Quantum Interference Device

We study the Atomtronics Quantum Interference Device employing a semiclassical perspective. We consider an $M$ site ring that is described by the Bose-Hubbard Hamiltonian. Coherent Rabi oscillations in the flow of the current are feasible, with an enhanced frequency due to to chaos-assisted tunneling. We highlight the consequences of introducing a weak-link into the circuit. In the latter context we clarify the phase-space considerations that are involved in setting up an effective"systems plus bath"description in terms of Josephson-Caldeira-Leggett Hamiltonian.


Introduction
Atomtronics is a new quantum technology [1,2,3], with potential for novel quantum computing implementations [4,5,6,7]. Theory and experiments with Atomtronic superfluid circuits are in the focus of current research [8,9,10,11]. A major objective is to realize a Quantum Interference Device (AQUID) that possibly includes one or two weak-links [12]. This is analogous to a superconducting circuit, or to its low dimensional version (fluxon, Josephson vortex qubit) [13,14]. However the design considerations of such device are still somewhat vague.
We study an Atomtronic superfluid circuit that is described by the Bose-Hubbard Hamiltonian (BHH) [6]. Namely, we consider N bosons in an M site rotating ring such that the model parameters are (N, M, K, U, Φ), where K is the hopping frequency between the sites, U is the on-site interaction, and the rotation is formally equivalent to having an Aharonov-Bohm flux Φ. If a weak-link is introduced, there is an additional parameter α that characterizes the relative strength of the coupling.
For the purpose of qubit realization, the objective is to single out a two-level system (TLS) that is quasi-isolated from all the other microscopic degrees of freedom (DOFs). In the present context there are two flow-states that differ by their "winding number" m, meaning that they are characterized by a different value of the persistent current (I m ). The flow-states are required to be meta-stable, meaning that each of them will not decay in time. If they are quasi-degenerate, one would like to witness coherent Rabi oscillations. During a Rabi-based protocol the system evolves into a superposition of macroscopically distinct flow-states [15].

Model and Methods
We consider N Bosons in an M site rotating ring such that the model parameters are (N, M, K, U, Φ), where K is the hopping frequency between the sites, U is the onsite interaction, and the rotation is formally equivalent to having an Aharonov-Bohm flux Φ. If a weak-link is introduced, there is an additional parameter α ∼ K ′ /K that characterizes the relative strength of the coupling. Accordingly the ring is described by the Bose-Hubbard Hamiltonian (BHH): where a j (a † j ) are bosonic annihilation (creation) operators on the jth site and n j = a † j a j is the corresponding number operator. Periodic boundaries are imposed, meaning that a M ≡ a 0 . The parameter U takes into account the finite scattering length for the atomic two-body collisions on the same site. The hopping parameters are constant K j = K except in the weak-link where it is K ′ . The ring is pierced by an artificial (dimensionless) magnetic flux Φ, which can be experimentally induced for neutral atoms as a Coriolis flux by rotating the lattice at constant velocity [18,19], or as a synthetic gauge flux by imparting a geometric phase directly to the atoms via suitably designed laser fields [20,21,22]. The presence of the flux Φ in Eq.(1) has been taken into account through the Peierls substitution: In the quantum analysis, we diagonalize Eq.(1), and display the spectrum as in Fig.1a. For each eigenstate E α we calculate the fragmentation measure M as defined in Appendix B, while the average current is obtained using the following formula: In a classical context the average is taken over time for a very long trajectory.

Semiclassical perspective
For the purpose of semiclassical analysis it is convenient to write the BHH using actionangle variables: a → √ ne iϕ . Accordingly the Hamiltonian describes an M degrees of freedom (DOFs) system, namely, Since the total number of particles N = n j is a constant of the motion of the system, the Hamiltonian above describes d = M−1 coupled pendula: the dimer (M = 2) is the so-called bosonic Josephson junction; while the trimer (M = 3) is the minimal superfluid circuit. Our main focus is on BHH circuits with M = 3, 4, 5 sites, but we shall refer to M ≥ 6 rings too. The interaction is characterized by the dimensionless parameter The classical dynamics is governed bẏ where z ≡ (ϕ 1 , · · · , ϕ M , n 1 , · · · , n M ) are the canonical coordinates. The notation ∂ ν stands for derivative with respect to z ν , and J is the symplectic matrix. It is important to emphasize that upon re-scaling the only dimensionless parameters that affect the classical trajectories are (u, Φ) and K ′ /K. The effective Planck constant is = 1/N. The latter parameter, does not appear in the "classical" equations of motion Eq.(5), but only in the full quantum treatment of Eq.(1).

System plus bath perspective
The conventional approach for analyzing a SQUID/AQUID is based on a "system plus bath" perspective. This perspective becomes meaningful once a weak-link is introduced, which is like having a "slow DOF". In order to motivate the conventional phenomenology one can regard the BHH Eq.(3) as describing masses that are connected by nonlinear springs. If one spring is very "weak", then at low energies the equal-partition theorem justifies an harmonic approximation for the small vibrations of the other springs. Accordingly we can regard the system has having one non-linear DOF ("pendulum") coupled to phonons ("harmonic bath"). The canonical coordinates that describe the weak-link are the phase difference ϕ = (ϕ M − ϕ 1 ), and its conjugate n = (n M − n 1 )/2. Hence we obtain the the Josephson Circuit Hamiltonian (JCH) with E C = U, and E L = [(N/M)/(M − 1)]K, and E J = (N/M)K ′ . The bath Hamiltonian has the standard Caldeira-Leggett form For small M the "bath" is merely a set of several oscillators, and possibly can be neglected, because the ω m are typically large compared with the natural frequency of the junction. For large M one can characterize the bath oscillators by an Ohmic spectral function The detailed derivation and the explicit expressions for the bath parameters in terms of the BHH parameters are presented in Appendix C, and will be further discussed in a later section. We note that in [13,6] the finite-temperature partition-function of the BHH ring has been introduced, and the reduced "system plus bath" action has been deduced. From the reduced action one could figure-out what is the effective JCH. In the present approach to the same system, we do not assume finite temperature, but merely re-arrange the Hamiltonian in a way that allows a "system plus bath" description. This is a valid procedure even if the ring is prepared (say) in a micro-canonical state with some arbitrary energy E. One may say that in our treatment E/M plays the role of the temperature.
Within the framework of the JCH treatment, the possibility of having metastable flow-states is controlled by the parameter For Φ = π the condition for having at least two local minima in the potential floor of the JCH, is α > α c , where α c = 1. Disregarding small quantum fluctuations, the two local minima can support a quasi-degenerate pair of flow-states. If the bath is ignored, then from the WKB approximation it follows that the tunnel splitting is given by some variation of the following expression [13] where C is a numerical prefactor. We would like to emphasize that there are several variations of this formula, depending on the relative size of (E C , E L , E J ), but they are all based on the assumption that Eq.(6) is a valid description.

Two-level system perspective
The objective is obviously to realize a two-level system (TLS) that is quasi-isolated from all the other microscopic DOFs [23,24,25,26,27,28,7]. In the present context there are two flow-states that differ by their "winding number" m, meaning that they are characterized by a different value of the persistent current (I m ). We label these states as and , and write the TLS hamiltonian as We refer to ∆ s as the splitting: if we draw the eigen-energies versus Φ we get an avoided crossing. The flow-states are required to be meta-stable, meaning that each of them will not decay in time. If they are quasi-degenerate, one would like to witness coherent Rabi oscillations. The quasi-degeneracy is controlled by Φ, and happens for Φ = 0 (say m = ±1) or for Φ = π (say m = 0, 1). During the Rabi oscillation the system evolves into a superposition of these macroscopically distinct flow-states. Such superposition is commonly termed "cat state". The conventional procedure to engineer a TLS is as follows: (i) To introduce a ring with a weak-link; (ii) To ensure that the weak-link DOF is only weakly-coupled to all the other ring DOFs; (iii) To analyze the operation of the device using the "system plus bath" paradigm of Caldeira and Leggett. The introduction of a weak-link allows the reduction of the many-body BHH Eq.(1) into the simpler JCH Eq.(6). The JCH consists of a single pendulum-like DOF that interacts weakly with harmonic-oscillators ("phonons"). The relative strength of the weak-link is characterized by a parameter α of Eq.(9). If α < α c the weak-link destroys the meta-stability (it is effectively like having a disconnected ring), hence we require α > α c . Furthermore, the introduction of a weak-link allows control over the coupling ∆ s between the flow-states. Without a weak-link this coupling might be too small for operational purpose, meaning that the time period (2π/∆ s ) of coherent Rabi oscillations might become too large for practical implementations.
If the bath is taken into account then there are two effects. One is "dressing" of the bare parameters, and the other is "noise". It is well known from the work of Caldeira and Leggett that coherent Rabi oscillations can be observed provided η < η c , where η c is of order unity (η c = π for the spin-boson model). We shall come back to this issue when we discuss the large M limit.

Coherent dynamics in the absence of a weak-link
The stationary orbitals of a single particle in a clean ring are the momentum states with wavenumber k = (2π/M)m, where m is an integer modulo M. Coherent flow-states have N particles condensed into the same momentum orbital: Implying a macroscipically large current In the absence of interaction (U = 0) these coherent flow-states are the eigenstates of the BHH. For Φ = π the m = 0 and m = 1 flow-states are degenerate in energy. If we add not-too-strong interaction they become coupled and may form a doublet whose dynamics is generated by the TLS hamiltonian Eq.(11). The energy-difference δE ≡ E − E is determined by the deviation δΦ ≡ (Φ − π), and the coupling ∆ s is determined by the strength of the interaction. An example for such doublet if provided in Fig.1.
Assuming that we have a TLS doublet of flow-states with energy splitting ∆ s , one would expect to witness pure Rabi oscillations. If the system has been prepared (say) in a flow-state with clockwise current, the subsequent evolution would be implying alternating current with frequency ∆ s , namely, If we add weak-link or weak-disorder, the flow-states remain stable, provided the perturbation is not strong compared with the interaction. This is the essence of superfluidity. The stability is due to the non-zero interaction U. The interaction stabilizes the flow-states: instead of being located on a flat potential floor, the flow-states are located in local minima of the potential floor. Local minima are structurally-stable with respect to the added disorder, i.e. the local minima do not diminished by a weak perturbation. The common conception is that the the two minima are separated by a "forbidden region". This is the same reasoning that leads to Eq.(10), but here we refer to the multi-dimensional phase-space of the BHH Eq.(3) and not to the reduced single DOF description of Eq.(6). Nevertheless, both perspectives connect smoothly. Namely, a rough way to write Eq. (10), that illuminates the semiclassics is where the prefactor C M has some dependence on M. This version highlights the distinction between the "classical" parameters (u, α) and the "quantum" parameter = 1/N. In the absence of a weak-link one formally makes the sunstitution α → (M−1) as implied by Eq.(9). The energy scale ∆ 0 ≡ (E L E C ) 1/2 is like the "attempt frequency" of the Gamow-formula. In a later section we identify ∆ 0 as the frequency spacing between the phononic modes.
In the JCH based picture, the splitting ∆ s is exponentially small in N due to the existence of a classically "forbidden region" between the two local minima, which necessitates tunneling. This very small ∆ s creates difficulties in witnessing coherent two-level dynamics in such configuration. In order to have a bigger ∆ s a smaller α is required. But is should not be smaller than α c = 1 else the meta-stability is diminished. Note also that there is a trade-off between the weakness of the link and the quality of the superposition state [27].
The question arises whether one can manage without introducing a weak-link. In fact there is a loophole. In order to realize this loophole, one should be aware, following [17], that there are novel flow-states that are not supported by local minima of the potential, but by a "stability island" or by a "chaotic pond", or by an "Arnold web" region. We summarize all these possibilities in Appendix A -the exact details are not important. The important point is that the phase-space locations, where the flowstates reside, are not separated by a "forbidden region". Instead they are separated by a "chaotic-sea". A visualization of this possibility is provided by the quantum spectrum in Fig.2, which should be contrasted with that of Fig.1. The way we plot the quantum spectrum (following [17]) is in one-to-one correspondence with a section of the classical phase-space: In Fig.1 the two flow-states at the bottom are separated by a "forbidden region" where no states can reside; In contrast to that, in Fig.2, between the two metastable states there are many other states with roughly the same energy that reside in the "chaotic-sea".
If the coupling between the quasi-degenerate eigenstates is mediated by a chaotic sea, then ∆ s is much larger. This is known as chaos-assisted tunneling [29,30,31,32]. Possibly the term tunneling is not the best description for the mathematics that is involved. The rough idea is that the quantum-coupling between the two metastable states is mediated by some intermediate state in the chaotic sea. The coupling is roughly estimated using second-order perturbation theory as ∆ s ∼ U 2 /∆, where ∆ is the detuning from exact resonance. This expression does not contain a WKB suppression and are separated by a forbidden-region. The tunnel-coupling allows coherent Rabi oscillations with extremely slow frequency ∆ s . If we did not slightly perturbed Φ, the diagonalization would give zero current cat-states (symmetric and anti-symmetric superposition of the pertinent flow-states). In (b) the initial state is an m = 1 coherent state, and the system has exactly Φ=π rotation. This initial state has large overlap with the pair of quasi-degenerate cat eigenstates. Consequently we observe Rabi oscillations of the current with frequency ∆ s that is determined by the tunnel coupling.  exponent, so it is not small, but nevertheless it is very sensitive to the model parameters, as in the theory of universal conductance fluctuations.
In Fig.2 we provide a numerical demonstration of chaos-assisted Rabi oscillations. In this example the device is non-rotating (Φ=0), and the Rabi oscillations are between the metastable m = ±1 flow-states. The dependence of ∆ s on the number of particles for "chaos assisted tunneling" is contrasted with "under the barrier tunneling" in Fig.3.
Summarizing this section, we observe that the coupling between metastable flowstates can be via chaos-assisted tunneling, implying a relatively large ∆ s when compared with the conventional expectation. A weak-link in a few-site ring is not essential for getting large ∆ s . In fact its introduction is likely to be harmful for the device operation (see next section).

weak-link in a few site ring
In this section we discuss what happens if a weak-link is introduced into a ring that has a small number of sites (M = 3, 4). In particular we ask what remains of the JCH phenomenology. The first implication of the JCH phenomenology is the prediction of a critical α below which a quasi-degenerate doublet of flow states cannot exist. If we naively use Eq.(9) we deduce that the condition α > 1 for getting such doublet is K ′ /K > 1/2 for M = 3 and K ′ /K > 1/3 for M = 4. In order to inquire what is the actual threshold we plot quantum spectra for various values of u and K ′ /K. See   Quantum spectrum (upper panels) and phase-space landscape (lower panels). The quantum spectra are for an M =3 ring with N =45 particles, dimensionless interaction u = 2.5, and weak-link coupling ratio K ′ /K = 1, 0.8, 0.65, 0.4 (from left to right). Axes and and color code are the same as in Fig.1. In each case an n 3 −n 1 = 0 Poincare section is displayed. The section coordinates are Q = (n 1 − n 2 )/(2N ) and P = (ϕ 1 − ϕ 2 ). The energy is chosen to be slightly above the ground state. The solid black line marks the borders of the allowed phase-space regions. The outer regions are "forbidden" energetically. The color code represents the averaged current for each classical trajectory: red for larger clockwise current; blue for large anti-clockwise current; and yellow-to-green for very small current. Looking at Fig.4 we see that for rings with M = 3, 4 sites, the α border is slightly higher than expected. We have verified using Poincare sections (see below) that for large u the border is in agreement with α c = 1. For completeness we also show that for very large u (of order N 2 ) the value of M for the ground-state becomes of order M, reflecting the Mott transition [16].
To understand what determines the α c border we display in Fig.5 so-called Poincare sections of classical trajectories that are generated by the Hamiltonian Eq.(3). Namely, for display purpose a pair of canonical coordinates (Q, P ) is selected, and for each trajectory the sequence of points (Q(t j ), P (t j )) where it intersects a specified phasespace section is recorded. We see clearly that in the α < α c regime the two stability islands merge, reflecting that we no longer have the "double well" structure in phase space.
But this is not enough. The JCH should be trusted also when we analyze tunneling or phase-slips through the forbidden region. For this purpose it should describe correctly the dynamics up to some energy well above the barrier. This means the the threshold E u for chaotic motion should be above the threshold E b for barrier crossing. See illustration in Fig.6. We therefore plot in Fig.7b, a Poincare section for an energy that is slightly above E b . What we see is that trajectories that go across the barriers are chaotic rather than regular. This indicates that a JCH description of the dynamics is in fact not valid.
Let us try to understand the reason for the failure of the JCH description. In the vicinity of a single flow-state worst case scenario is that a phase difference π has to be supported by the ring. The harmonic approximation requires π/(M − 1) < π/2 on each bond. This is marginally satisfied for M = 3. But if we want the JCH description to be  The lower panels are for the weak-link energy E ϕ (t) measured relative to the top of the barrier. The system parameters are u=200, and K J /K=0.3, and Φ=π. In all the panels the initial condition is in the vicinity of the barrier, with equal populations n j = N/6. The actual starting point is with n 1 = (1/6 + δ)N and n 6 = (1/6 − δ)N . In (a) the junction energy is mostly below the barrier (δ = 0.01), and we see that the dynamics is in qualitative agreement with the JCH: we observe regular flow-motion with rare jumps to the opposing flow-motion due to an activation by the "bath" DOFs. In (b) the junction energy is above the barrier (δ = 0.08), and we still observe pendulumlike regular motion. In (c) the energy is above the chaos threshold (δ = 0.16), and we get irregular chaotic motion that is no longer described by the JCH. This should be contrasted with the M = 3 trajectories of Fig.7b where the chaos threshold E u coincides with E b , invalidating the JCH phenomenology.
valid over a 2π range of ϕ, then the requirement becomes 2π/(M − 1) < π/2, meaning we have to consider rings with M ≥ 5 sites. Similar claim has appeared in [13]. In Fig.8 we verify that for an M = 6 ring with weak-link the chaos border E u is indeed well above the barrier energy E b . Up to E u the dynamics looks like that of a pendulum that is slightly affected by the other "bath" DOFs. Above E u the motion becomes chaotic and the JCH description is no longer applicable.

weak-link in a many site ring
Consider N bosons is ring of length L, such that the average density is ρ = N/L. The so called Lieb-Liniger parameter that controls the quantum aspect of the interaction is γ = mg/ρ. For γ ≫ 1 the hard-core bosons are like fermions, while for γ ≪ 1 we can use a "classical" description. In the latter case the "trajectories" obey the so-called Gross-Pitaevskii (GP) equation. In fact the parameter γ does not appear in the GP treatment of the model. The only dimensionless parameter of the GP description is We shall refer to it as the "classical" dimensionless parameter, while = 1/N can be regarded as the dimensionless Planck constant. Within the framework of the "classical" (GP) treatment the low excitations of the systems are phonons with sound velocity c = (gρ/m) 1/2 . For a finite length ring the spacing in the frequencies of the phononic modes is ∆ 0 = πc/L. If we add a periodic potential that divides the ring into M sites, we get a system that possibly can be described by the BHH Eq.(1). The analogue of the GP is the discrete nonlinear Schrodinger (DNLS) equation. The distance between the sites is a = L/M and the average number of particles per site isn = N/M. The effective parameters of the BHH are accordingly U = g/a and K = 1 where S 0 reflects the height of the barrier. The effective quantum parameter is This parameter controls the quantum Mott transition. Namely for γ * > 1 superfluidity is diminished ifn is close to integer. In addition we can define the "classical" dimensionless parameter which is analogous to u L of Eq.(17) as The u M parameter controls the DNLS equation, and determines the stability of the steady flow solutions, as well as the thresholds for self-trapping and soliton formation. Due to the discretization we have effectively M phononic modes, whose spectrum is charaterized by the cutoff frequency Where ∆ 0 is formally the same as for a continuous ring, but with m * . For a regular ring with a weak-link the reduction to an effective JCH provides the following expressions [33]: E C = g/L and E L = ρ/(mL) and E J = αE L . The parameter α is controlled by the tunnel-coupling, which is determined by the height of the barrier at the weak-link. Our derivation in Appendix C has provided similar expressions, but there are some differences. First of all the effective mass is of course m * and not m, and therefore the effective quantum parameter γ becomes γ * . A secondary difference is that E C = g/L is replaced by E C = g/L * , where the effective length over which the density varies might be as small as L * = a. The latter value reflects the extreme case of uniform distribution of the particles along the ring. Expression Eq.(9) for the parameter α can be written as where S 0 and S J reflect the heights of the barriers in regular bonds, and at the weak-link respectively. We turn our attention to the bath. The derivation in Appendix C shows that within the bilinear-coupling approximation the effective number of bath DOFs is In Ref. [33], regarding regular ring, it has been claimed that if (E J /N) ≪ ∆ 0 (called there "the small ring limit") then the bath can be ignored. In the context of the present Bose-Hubbard circuit this condition takes the form K ′ ≪ ω c , meaning that the bath should have high frequency cutoff compared with the hopping rate. But from the work on the spin-boson problem we know that the condition for witnessing coherent oscillation is η < π which implies that γ * should be large compered with unity. We identify that this is a problematic non-semiclassical regime where the Mott transition takes place. Namely, for γ * > 1 the superfluidity of the system depends sensitively on the filling ratio N/M. In a grand-canonical perspective the system has the tendency to become a Mott insulator.

Discussion
We observe that a TLS modeling of quasi-degenerate flow-states in a few-site ring is feasible, meaning that coherent Rabi oscillations are not over-damped. This is true with or without a weak-link, and the frequency is possibly determined by chaos-assistance tunneling. In particular we have demonstrated numerically Rabi oscillations between metastable flow-states in a non-rotating (Φ = 0) circuit that consists of M = 4 sites.
We have determines what is the minimal value of α that does not endanger the meta-stability of the Φ = π flow-states. Clearly below this minimal value a weak-link is not useful. From a semi-classical perspective this value is the threshold for the merging of two stability islands. For large rings, assuming that the JCH phenomenology is valid, the minimal value is implied by the familiar condition α > α c with α c = 1 We note that in a super-conducting circuit, due to the Meisner effect, the effective inductance is larger, and α is typically large.
In the semiclassical perspective the flow-states are supported by a local minimum of the energy landscape (Landau stability), or by a region that is surrounded by KAM tori. In the latter case, for rings with M > 3 sites the stabilization is due to a manybody quantum localization effect, that suppresses the Arnold diffusion. Depending on the type of states involved, the coupling might be via a forbidden-region (as implied by the JCH phenomenology), or it might be mediated by a chaotic sea. In the latter case the chaos-assisted tunneling provides a weaker dependence on the number of particles involved.
The system plus bath perspective.-Formally the circuit has d = M−1 interacting DOFs, while in the approximated JCH version we have a single DOF (ϕ, n) that interacts with a "bath" that consists of a few DOFs. If the bath is ignored the motion in the single DOF phase-space is regular, and looks formally the same as that of a pendulum. If α > α c , a separatrix is formed, hence we have two stability-islands that can support the two quasi-degenerate flow-states. But if the bath is taken into account, the projected motion in the (ϕ, n) coordinates becomes "dressed" and "noisy", in the same sense as discussed by Caldeira, Leggett and followers. These effects endanger the coherent Rabi oscillations.
Large M ring.-For a regular ring with bosons one can define the Lieb-Liniger parameter γ. Having γ > 1 means that quantum effects become important (GP description becomes problematic), but nevertheless there is no quantum phasetransition. For the BHH ring (bose gas in an optical lattice), we have defined an effective γ * that corresponds to the effective mass in the lattice. As before γ * > 1 means that quantum effects are important. But here the consequences are much more dramatic. Namely, the quantum regime γ * > 1 is identified as the Mott-regime, where depending on the filling-ratio the ring can become a Mott-insulator. On the other hand the analysis shows that γ * > 1 is the condition for witnessing coherent Rabi oscillations. So there is clash here: on the one hand we want the ring to be in a superfluid phase (avoid Mott); on the other hand we want to have weak coupling to the bath in order to witness coherent oscillations. Thus for a many-site ring the requirement for observing coherent Rabi oscillation might be in clash with the quantum Mott transition.
Small M ring.-We wanted to understand how this standard JCH phenomenology is modified if the ring consists of a small number of sites. Then the "bath" consists of a small number of DOFs and the standard Caldeira-Leggett perspective becomes questionable. One direction [34] is to say that the interaction with chaotic DOFs is essentially like the interaction with infinitely many harmonic DOFs, hence coming back to Caldeira-Leggett phenomenology. This type of argument might work for rings with M ≥ 6 sites for which the effective number of bath DOFs is d bath ≥ 2. We did not take this route here. Rather we discussed the whole issue in a much more fundamental level, focusing on rings with M = 3, 4, 5 sites.
Arnold diffusion is in a sense the low dimensional version of having a "bath". The essence of Arnold diffusion is that a selected DOF does not perform an unperturbed integrable (pendulum like) motion. Rather the motion always "diffuses" due to the "noise" that is induced by the other DOFs. Hence we have here a formal equivalence with the "system plus bath" perspective. It follows rigorously that a necessary condition for the applicability of the "system plus bath" paradigm with regard to a circuit with a weak-link requires more than 3 sites. But this is not a sufficient condition. We have emphasized that a JCH modeling implies regular motion up to an energy that exceeds the barrier height. Such high threshold for chaos is apparently feasible only for rings that have more than 5 sites.