Quasi-molecular bosonic complexes -- a pathway to atomic analog of SQUID with controlled sensitivity

Recent experimental advances in realizing degenerate quantum dipolar gases in optical lattices and the flexibility of experimental setups in attaining various geometries offer the opportunity to explore exotic quantum many-body phases stabilized by anisotropic, long-range dipolar interaction. Moreover, the unprecedented control over the various physical properties of these systems, ranging from the quantum statistics of the particles, to the inter-particle interactions, allow one to engineer novel devices. In this paper, we consider dipolar bosons trapped in a stack of one-dimensional optical lattice layers, previously studied in [1]. Building on our prior results, we provide a description of the quantum phases stabilized in this system which include composite superfluids, solids, and supercounterfluids, most of which are found to be threshold- less with respect to the dipolar interaction strength. We also demonstrate the effect of enhanced sensitivity to rotations of a SQUID-type device made of two composite superfluids trapped in a ring-shaped optical lattice layer with weak links.


Introduction
The advances in experimental methods and techniques developed in atomic, molecular and optical physics experiments, such as atomic spectroscopy and interferometry, have led to the development of novel technologies, which have in turn impacted a diverse set of fields, including precision measurement and sensor technologies. Prominent examples include atomic clocks, atomic magnetometers, and atom interferometers [2]. The recent theoretical and experimental development in atomtronics, allows one to apply techniques from experimental atomic, molecular and optical physics, condensed matter theory, and quantum information to realize ideal analogues to various electronic systems with the atoms playing the role of the electrons. As such, these systems are set to play a significant role in the development of novel quantum technologies and quantum computation [3].
In the past decade, a variety of atomtronics devices such as atomtronic batteries, diodes, transistors, atom circuits have been proposed theoretically [3,4,5,6,7,8,9,10,11,12,13] and in certain cases realized experimentally [14,15,16]. More recently, a considerable effort has been devoted to the realization of atomic SQUID (superconducting quantum interference device), which are experimentally created by a Bose-Einstein condensate trapped in a toroidal trap with weak links provided by a variety of techniques such as rotating potential barriers, painted potentials, or Laguerre-Gaussian modes of a laser beam. These systems are analogous to superconducting rings with Josephson junctions, where the response of a superconductor to an external magnetic field is replaced by the response of a superfluid to rotation. The neutral versions of the SQUID were first realized in superfluids of 4 He and 3 He [17,18]. The equivalent of the rf-SQUID has been created in [19,20,21], while preliminary experiments for a dc-SQUID have been carried out in [22,23]. These experiments have paved the way to the atomtronic rotation sensors analogous to magnetic field sensors formed by superconducting devices. While the atom SQUID experiments listed above were performed in the continuum, recent theoretical proposals call for the realization of atomic rf-SQUIDS in ring-shaped optical lattices [24,25].
Ultracold atoms in optical lattices provide an ideal platform for engineering systems and devices in a highly controllable manner. The unprecedented level of control and flexibility of these experimental setups allow one to construct a variety of system geometries and manipulate inter-particle interactions in an almost ideal, decoherencefree setup. Moreover the experimental realization of atomic and molecular systems featuring long-range and anisotropic dipolar interactions [26,27,28,29,30,31,32,33,34,35,36] may lead to the experimental observation of the predicted novel phenomena such as p-wave superfluidity, superfluidity of multimers, solids and supersolids. In what follows we discuss quantum phases of hard-core dipolar bosons trapped in a multi-tube geometry, where each tube is a one-dimensional optical lattice and all tubes are parallel to each other and belong to the same plane. An external field is used to align the dipole moments perpendicular to the tubes. This system features a variety of solid phases, composite superfluids, and supercounterfluids, most of which are predicted to be stabilized at infinitesimal values of dipolar interactions [1]. In particular we use bosonization and renormalization group techniques, in association with large scale quantum Monte Carlo simulations, to study and confirm the threshold-less nature of these phases.
It is important to note that such phases can potentially be used for creating nonlinear elements (and their networks) which can be viewed as generalized Josephson junctions with interesting unexpected features. Unusual properties of a contact between paired and single atomic superfluids were first noticed in Ref. [37]. In general, a contact can be created between arbitrary phases. However, in order to have Josephson effect through such a contact, certain conditions following from the requirement of particle conservation should be met.
This paper is organized as follows. In Section 2 we describe the microscopic Hamiltonian of the system. Section 3 discusses the quantum phases stabilized by this model. Section 4 focuses the theoretical treatment of the bilayer case, presenting bosonization formalism along with numerical results from quantum Monte Carlo simulations. In Section 5 we summarize results from bosonization and renormalization group analysis for the multilayer system presented in [1]. Section 6 is devoted to discussing the nature of Josephson effect between different composite superfluids made of the same type of bosons. Specifically, we derive a dependence of the rotational "flux" quantization on the number of components forming the composites, and point out that this dependence may be utilized for achieving higher sensitivity in rotational sensors. Finally, in Section 7 we conclude.

Hamiltonian
The system under consideration is described by the model where J > 0 is the intra-tube tunneling amplitude, a † xz (a xz ) is the creation (annihilation) operator for a hard-core boson at site (x, z), where z = 0, 1, 2, ..., N − 1 labels the tubes, which are parallel to each other and belong to the same plane, and x = 0, 1, 2, ..., L is the coordinate along a tube. The summation over the nearest neighboring sites is denoted by x,x ,z .... We consider zero inter-tube tunneling, reflected in the absence of terms a † xz a x z with z = z, so that the chemical potential µ z in layer z, controlling its filling factor ν z , can be chosen to be different in different tubes. Thus, the only coupling between tubes is provided by the dipolar poetntial. Here n xz = a † xz a xz is particle number operator. Figure 1a) shows a sketch of the system considered. For the geometry considered here where the polarization axis of the dipoles is along z-axis, the interaction V (x, z) takes the form where V d > 0 sets the energy scale determined by the inter-tube distance d z as V d ∝ 1/d 3 z , and α denotes the ratio of squares of distances between neighboring sites along x and z, respectively. In what follows we use α = 1. In the case of polarization axis perpendicular to the xz-plane, instead, the dipolar interaction is purely repulsive.
Here we consider hard-core bosons (that is, n xz = 0, 1 in Eq.(1)), which guarantees the stability of the system. For polar molecules trapped within this geometry, the hardcore condition can be satisfied by using a deep enough optical lattice in the z−direction, which prevents system collapse due to the attractive forces between aligned dipoles [38]. For magnetic atoms, one can reach the hard-core limit by exploiting Feshbach resonances which can be used to tune the attractive interactions. Moreover the use of a deep optical lattice along z together with varying d z can effectively suppress inter-tube tunneling.
In the following we show that in this geometry, where the interaction along the zaxis is attractive, any arbitrarily small V d can stabilize superfluidity of quasi-molecular complexes, provided that the bosons satisfy the hard-core condition.

Quantum phases in multi-layer geometries
In the absence of inter-tube hopping between the layers and the inter-tube interactions, hard-core bosons form N independent superfluids (N -SF), which are characterized by N quasi-condensates exhibiting algebraic quantum order. Within the low energy description the second quantized operators ψ z (x) built on a xz -operators can be replaced by phases φ z (x) as ψ z (x) → exp(iφ z (x)), where φ z is a variable belonging to the tube z. In the N -SF phase the quantum correlators exp(−iφ z (x)) exp(iφ z (0)) ∝ δ z,z |x| −b , where b > 0 is determined by the Luttinger parameter.
As shown in [1], upon introducing an infinitesimal inter-tube interaction, the N -SF phase may be destabilized in favor of composite superfluids, supercounterfluids and insulating phases. In a composite superfluid phase (CSF) the one-particle algebraic order vanishes, that is, the power law is replaced by the exponential decay exp(−iφ z (x)) exp(iφ z (0)) ∝ δ z,z exp(−|x|/λ), where λ is some length determined by interaction V d . In the situation of zero threshold λ → ∞ as V d → 0 [1]. This, in fact, determines the practical experimental limitation -a system size L should be larger than λ in order to observe CSF.
It is important to note that in the CSF the algebraic order is preserved only in the corresponding products of the one-particle operators. Here we provide a brief description of the composite quantum many-body phases.  Figure 1. Sketches of 1D quantum phases in the multi-tube geometry as described in Sec.3: a) N -SF phase existing for fillings of bosons (solid circles) different in each tube so that there are independent one-particle algebraic orders in each tube. The solid arrow shows the direction of tunneling along x; b) The CSF phase realized between tubes with equal filling. Dotted lines indicate binding of bosons from different layers and the fuzzy cloud depicts (three-particle) algebraic order for all layers; c) The SCF phase realized between layers with complementary fillings (upper two) so that the algebraic (two-particle) SCF order (shown by the fuzzy cloud) is formed only on these two layers. Dotted lines indicate binding between particles and holes (open circles). The bottom layer has filling incommensurate with both the upper two layers and, thus, features a one-particle algebraic order; d) The CB phase. Dotted lines indicate binding between bosons from different layers at ν = 1/2.

N atomic superfluid (N -SF)
The N -SF phase, with N independent atomic superfluids, is stable in the presence of finite V d if the interacting tubes have fillings ν z that are both different and noncomplementary. In this case, there exist N independent one-particle algebraic orders (at T = 0) in the spatial correlators for each tube as discussed above. This scenario is depicted in Fig. 1 a).

Composite superfluid (CSF)
If n ≤ N tubes in the system have the same incommensurate filling factor, any intertube attraction destabilizes the independent atomic superfluids in favor of a superfluid of n-boson composites. This phase is characterized by an exponential decay of correlators exp(iφ z (x)) exp(−iφ z (0)) as discussed above, coupled with the survival of the algebraic order in the n-body density matrix Ψ † (x)Ψ(0) where Ψ(x) = ψ 1 (x)ψ 2 (x)...ψ n (x), with ψ i , i = 1, 2, ..., n denoting the fields characterized by the same value of the filling factor ν. This defines the CSF -a superfluid of quasi-molecular complexes which are formed by n bosons on n tubes. A sketch of this phase is shown in Fig. 1b). It is useful to note that, due to the non-local nature of the interaction, such a phase can form not necessarily between adjacent tubes. The only requirement is that the filling factors are the same in a set of n tubes.

Supercounterfluid (SCF)
If n tubes in the system have complementary fillings, the N -SF phase can be destabilized in favor of the SCF, previously introduced for a two-component bosonic system in Ref. [39]. SCF can exist in a lattice when the filling factors ν 1 and ν 2 for two components complement each other to an integer filling. In the case of hard-core bosons it is ν 1 + ν 2 = 1. Then, the repulsive interaction can induce binding of atoms of the component "1" to holes of the component "2". This corresponds to the algebraic order in the product ψ * 1 (x)ψ 2 (x), while the individual fields ψ 1 and ψ 2 exhibit no algebraic order.
This property can be extended to a general case of an N -component bosonic mixture where the super-flow of M < N components is (partially) compensated by the counter-flow of the remaining components, provided that the filling factors of the two sets are complementary. The SCF phase is sketched in Fig. 1c) for the case of two (upper) tubes being complementary to each other and the third one (at the bottom) exhibiting independent one-particle algebraic order. In other words, the SCF algebraic order exists for the fields ψ * 3 ψ 2 and the field ψ 1 exhibits one-particle algebraic order, where z = 1, 2, 3 from the bottom upward. The corresponding filling factors obey the relations ν 1 = ν 2 = ν 3 , ν 2 + ν 3 = 1.

Composite insulators
The system described by the model (1) can support insulating phases at various commensurate fillings ν z = 1/2, 1/3, 2/3, 1/5, .... In a single tube, that is, N = 1, such phases can occur only if V d exceeds corresponding thresholds. The situation is different for the cases ν z = 1/2, N > 2. As shown in Ref. [1], the checker-board (CB) insulator is formed even in the limit V d → 0. It is interesting to note that such an insulator is formed even if the inter-layer interaction is purely attractive. A sketch of the CB solid is shown in Fig. 1d).

Comparison of the bosonization results for N = 2 tubes with ab initio QMC
In this Section we review the bosonization approach to the CSF with N = 2 [40]so called paired superfluid (PSF). The bosonization framework [41] is the low energy description of a microscopic model. We provide detailed comparison with the QMC simulations [1] using the two-worm algorithm [42]. The excellent agreement between the two sets of results provides a strong ground for the extension of the bosonization treatment to arbitrary N , as previously presented in [1]. We will explicitly demonstrate that the PSF phase is stabilized by arbitrarily small V d .

Renormalization Group (RG) equations for N =2
The PSF in 1D is characterized by algebraic order in the product Ψ(x) = ψ 1 ψ 2 , while the individual fields ψ 1 and ψ 2 are disordered. We start by reviewing the bosonization approach for a single tube. The fields ψ z are represented in terms of the density ρ z (x) and the conjugated phase with θ z (x, t) being the density angle [41]. In the low energy limit, 1D superfluid can be described by the hydrodynamics of Luttinger liquid. For a single tube (that is, N = 1) its classical action in imaginary time τ can be written as where V s stands for speed of sound and K is the Luttinger parameter. It should be mentioned that these parameters are related to details of microscopic theory and finding these relations is not a trivial task. However, for the case of hard-core bosons these relations are known. In particular, in the absence of interaction (except for the hardcore condition), K = 1 and V s = J (in the atomic units) (see in [43]). As we will discuss below the relation K = 1 turns out to be crucial for the threshold-less nature of some of the composite phases. This important aspect has already been emphasized in Ref. [40] for the case of PSF in N = 2 system. In presence of a periodic lattice and interactions between bosons (in addition to the hard-core one), the higher harmonics in Eq.(3) must be taken into account. Then, the density-density interaction V (on different sites) introduces a variety of non-linear terms ∼ V cos(2m (θ 1 + πν 1 x)) cos(2m(θ 1 + πν 1 x)) inducing the so called backscattering (see in Ref. [43]), with V being the backscattering amplitudes determined by interaction (as V ∝ V d in our case). These terms may become relevant if ν 1 is commensurate with lattice period and K is below some critical value K c . For ν 1 = 1/2 the lowest relevant harmonic is the m = ±1 with the corresponding critical value K c = 1/2. Thus, a weak interaction V between hard-core bosons does not destroy 1D superfluidity in a single tube because K = 1 > 1/2. However, as the interaction becomes stronger, K decreases and eventually the insulating CB phase appears. As discussed in Ref. [1], the situation is dramatically different for the case N > 2 since the critical value of K is unity, K c = 1. Hence for N > 2 an arbitrarily small dipole interaction between tubes induces the CB phase.
Here, we consider the case N = 2 and ν 1 = ν 2 = 1/2, previously studied in detail by bosonization in Ref. [40]. We test these predictions by means of QMC simulations and find perfect agreement. For the sake of consistency we repeat the derivations so that the "language" of the comparison with QMC becomes clear.
The low energy action S for N = 2 tubes with equal fillings ν 1,2 = ν corresponding to the Hamiltonian (1) and which determines the partition function ∼ exp(−S) can be expressed in terms of the variables θ 1 and θ 2 as where the inter-tube interaction ∼ V ρ 1 ρ 2 is due to the dipolar interaction. The relationship between V and the interaction V d given by Eq. (2) is, in general, nontrivial. For weak interactions, the low energy limit is sensitive to the lowest harmonic of interaction for the so called forward scattering interaction V f s (to be used for the m = 0 harmonic in Eq. (3)) as well as to the backscattering amplitude V bs [1]. In the limit V d → 0 it is natural to expect that both V f s and V bs are ∝ V d . In particular, the inter-tube forward scattering interaction term takes the form where "-" is for the case of attractive interaction between the tubes (that is, the dipoles are along z-axis, see Eq. 2) and "+" for the repulsive interaction realized when the dipoles are perpendicular to the xz-plane. This terms should be added to Eq.(4) for both tubes. In what follows we consider the case of attractive interaction, that is, V f s < 0. This additional term S 12 changes the Luttinger parameters for channels θ + = θ 1 +θ 2 and θ − = θ 1 − θ 2 as so that K + > 2 and K − < 2. Ignoring weak change of V s , one finds where V + ∼ V − ∼ V d are the backscattering terms taken for m = ±1. For ν different from 1/2, the term V + can be dropped so that the first equation describes PSF which is a Luttinger liquid of pairs. In other words, there is algebraic order in the field Ψ = ψ 1 ψ 2 .
The situation with the channel S − is completely different. Since the critical value of K − is 2, the backscattering term becomes relevant so that the fluctuations of the angle θ − are suppressed. This implies that the superfluidity in the counter-flow channel is destroyed because density and phase are conjugate variables. Thus, suppression of fluctuations in θ − immediately implies strong fluctuations in the relative phase φ − = φ 1 − φ 2 so that the individual fields ψ 1 ∼ exp(iφ 1 ) and ψ 2 ∼ exp(iφ 2 ) lose their coherence.
The corresponding renormalization group (RG) equations (ignoring weak renormalization of the speed of sound [40]) are: where g = K − /2 and u ∝ V − with the proportionality coefficient determined by the small-distance cutoff. These equations are the standard Kosterlitz-Thoulless [44] RG equations, where l sets a typical space-time scale (in units V s = 1). At small scales, the flow of K − (l) starts at the initial value set by the forward scattering, Eq. (5). Since K − (0) < 2, Eq. (9) indicates that V − becomes more and more relevant as l increases. Simultaneously, the effective Luttinger parameter g(l) flows to zero as l → ∞. The comparison of the prediction of Eqs. (8) and (9) with ab initio simulations can be done for g(l). As discussed in the next section, the worm algorithm has a direct access to g(l) where l is determined by the ratio of the system size L to some microscopic scale L 0 , that is, l ∝ L/L 0 , provided that the extension along imaginary time /T = L/V s establishes the space-time symmetry. We will conclude this section by presenting the solution of the RG equations.
A general solution of the system of equations (8) and (9) can be expressed in terms of two constants of integration, η, l 0 > 0, determined by the initial values of u and g, which in their turn are set by the microscopic model (1). If η is real, the solution has a form On the other hand, if η = i|η|, the solution becomes This solution is valid if ξ > −1.
As will be shown in the section 4.2 the class of solutions pertaining to our simulations data belongs to the separatrix, that is, η → 0 and u = √ 2|ξ| where This equation will be used for the analysis of the QMC simulations where K − will be determined directly in a progression of system sizes for different values of the interaction.

Quantum Monte Carlo Results
For simplicity, we restrict the interactions to nearest-neighbor inter-layer attraction. This corresponds to small filling factor, ν << 1, when the dipole interactions along the tubes can be ignored.
Our QMC simulations are focused on finding statistics of windings of particle worldlines. The system satisfies periodic boundary condition along space and along imaginary time. Then, these windings are well defined integers W x (z), z = 1, 2, along space and W τ (z), z = 1, 2, along time. As shown in Ref. [45], the mean square fluctuations of W x and W τ determine the superfluid stiffness and superfluid compressibility, respectively.
The Luttinger parameter K − can be expressed in terms of W x and W τ as follows (see [1]): Within the simulations, the RG scale l is set by the system size L. Practically, we have kept L ∝ β so that (W x (0) − W x (1)) 2 = (W τ (0) − W τ (1)) 2 , in order to ensure space-time symmetry. This guarantees that the system remains in its ground state for a given size. Fig. 2 presents the raw data for K − /2 obtained from simulations performed at various system sizes. These results allow us to identify the class of solutions to the RG flow best describing the data, and to determine the integration constant l 0 . We find that within the statistical errors of the simulations, the curves of K − /2 vs L/L 0 for various 0 < V d /J < 1 collapse on to one master curve, shown in figure 3 which is well described by the separatrix solution given by Eq. (12). The master curve has just one fit parameter l 0 . Thus, for a given value of the interaction V d /J, we have found a unique L 0 = L 0 (V d /J) so that the set of the data K − (L) (for a given V d /J) from Fig. 2 can be best superimposed on the separatrix solution where l/l 0 = L/L 0 . The result of this procedure for 11 values of V d /J (with the total of ∼100 data points) is presented in Fig. 3 where the solid line is the separatrix (12).
The inset of figure 3 shows a plot of L 0 as a function of V d /J. While L 0 is a rescaling parameter, it can also be interpreted as the size of a bound dimer. As can be seen from the inset of Fig. 3, L 0 diverges with decreasing inter-tube interaction, where a is a constant ( a = 3.82). This divergence proves that the critical value for the formation of the dimer superfluid is V d = 0. The excellent agreement between the numerical data and the RG solution (12) over almost 50 orders of magnitude of the (effective) distances as well as the dependence (14) proves that the PSF phase is stabilized for infinitesimal attraction between the layers and that the bosonization can describe very accurately the system properties. At this point a comment about the meaning of the master curve is in order. Using small V d in a system of finite size L may not result in an experimentally "visible" formation of PSF because K − is too close to K − = 2. However, for larger size L, K − will "flow" along the master curve (solid line), i.e. the separatrix, toward smaller values. Such an approach -matching the numerical solution with the RG flow in order to find the critical point of the Berezinskii-Kosterlitz-Thouless transition has been pioneered in Ref. [46].

Results from N-tube bosonization
The formalism of the previous section can be extended to the multilayer case. For a thorough discussion of the N -tube bosonization we refer readers to [1]. The bosonization analysis for the general case is more involved since it is not possible to rewrite the action in terms of a sum of actions corresponding to decoupled sectors as it was done for N = 2. Likewise the N = 2 case, though, the goal is to study the flow of the amplitudes associated to the various harmonics (as defined in Eq. 3). Here, we present a short qualitative discussion.

CSF.
Similarly to the N = 2 case, the CSF phase occurs for arbitrarily small inter-tube attraction V d as long as the filling factors in the tubes are the same. Due to the spatially non-local nature of dipolar interactions, the CSF phase 3.2 can form between tubes with the same filling factors regardless of their positions. For example, in a system of N = 6 tubes where ν 1 = ν 2 = ν 5 = ν (here we consider ν = 1/2) with all other values ν 3 = ν 4 = ν 6 = ν, a CSF will form between tubes z = 1, 2, 5 while the remaining tubes will form 3 independent SF. Consequently, this system features algebraic orders of individual fields ψ 3 (x), ψ 4 (x), ψ 6 (x), as well as of the CSF field Ψ 125 (x) = ψ 1 (x)ψ 2 (x)ψ 5 (x), with the individual fields ψ 1 (x), ψ 2 (x), ψ 5 (x) being disordered.
The important aspect of the solutions found in Ref. [1] is that, in the presence of intra-tube repulsion V 0 > 0 provided by dipole forces, matching of the filling factors rather than the sign of the inter-tube interaction can determine the nature of possible phases. The CSF phase can also occur due to the inter-tube repulsion (for dipoles polarized perpendicular to the xz-plane) as long as the intra-tube repulsion is dominating, and provided the filling factors of the tubes forming CSF are the same.
Such a peculiarity -attraction due to repulsive interaction -is the property of hard-core bosons in 1D [1], as it can be seen from the renormalization of the forward scattering amplitude discussed in Eq. (5). If both the intra-and inter-tube repulsions exist, this equation becomes Thus, the Luttinger parameter K − for the channel θ − becomes K − < 2 if V 0 > V f s , where V f s is due to inter-tube repulsion. This, as discussed above, implies the suppression of the superflow in this channel (for equal filling factors). [It should be noted that, despite K + also being below 2, the superfluidity in this channel is protected by the filling factor being different from 1/2, as it was discussed for N = 2 layers].

SCF.
The supercounterfluid phases can occur between tubes with filling factors complementary to unity. These phases are thresholdless in V d as well [1]. If V 0 > |V f s |, the SCF occurs for attractive inter-tube interaction as can be seen from Eq.(15), provided the filling factors are complementary and different from 1/2.
An interesting situation occurs when there are groups of tubes with complementary fillings. Say, there are two sets of tubes, M 1 > 1 and M 2 > 1 such that in the first one the filling factor in each tube is ν = 1/2 and in the second one it is 1 − ν. Then the RG equations show the formation of a pair of composite superfluids-one per each group of tubes with equal filling factors-and that these composite superfluids further bind in the counterflow regime.
In the absence of inter-tube interactions, hard-core bosons can form a CB insulator at filling factor ν = 1/2 only if the repulsive interaction is strong enough. The situation changes dramatically in the presence of inter-tube interaction. RG equations show that regardless of the sign of the inter-tube interaction even in the limit V d → 0, a CB solid can be stabilized. Accordingly, there is no algebraic order in any channel.

Josephson effect between composite superfluids
One can engineer coherent transport through a link between two different composite superfluid, where the transport is controlled by specific powers of the order parameters corresponding to the two systems across the link (cf. [37]). In the following we use examples to illustrate this point.
Consider a CSF characterized by the field Ψ k = ψ 1 ...ψ k formed out of k components. As discussed in Ref. [1] the transition from k independent condensates to a single condensed field Ψ k corresponds to the loss of coherence of each individual field ψ i , while coherence is preserved in Ψ k . To emphasize the number of the components we will call such a CSF as k-CSF. Here we describe the Josephson coupling between the k-CSF and another phase, for example n-CSF with n = k or a simple superfluid exhibiting one-particle algebraic order in a field ψ. This coupling should preserve the global U(1) invariance, which arises as a consequence of the particle number conservation in the system. This requirement excludes couplings of the type ∼ (Ψ * k Ψ n + c.c.) or ∼ (Ψ * k ψ + c.c.). We begin by considering k tubes where the region x < 0 is occupied by k-CSF while the region x > 0 contains k-SF, that is, k independent superfluids characterized by one-particle fields ψ i , i = 1, 2, ..., k each exhibiting algebraic order. If there is a weak link at x = 0, the corresponding (minimal) Josephson coupling is given by the product ∼ Ψ * k k i=1 ψ i + c.c.. In terms of the corresponding phases the Josephson energy can be written as If the fields ψ i are all identical, that is their phases are locked in as one and the same phase φ due to the tunneling between the tubes in the region x > 0, the above relation takes the form In other words, the coherence is preserved only in the transport of k bosons across the junction, and any transfer of a number of bosons different from k would cost finite energy and, therefore, cannot be coherent. It should be noted that higher-order couplings would allow coherent transport of any integer number of groups of k bosons. We, however, are focusing on the lowest order coupling. Next, consider a junction formed between two composite superfluids k-CSF and n-CSF with k > 1, n > 1 and n = k. In the following we consider the case where k and n have no non-trivial common divisor. The corresponding fields exhibiting algebraic orders can be represented as Ψ k = k i ψ i and Ψ n = n i ψ i . In this case, the minimal Josephson effect can exist between the products of the fields (Ψ n ) k and (Ψ k ) n with the corresponding Josephson energy being where Ψ k = |Ψ k | exp(iϕ k ), Ψ n = |Ψ k | exp(iϕ n ). This form preserves global U(1) invariance of the system. In other words, the coherent transport occurs by groups of nk atoms. It is worth mentioning that transferring n atoms from n-CSF to k-CSF would require a compensation for the difference k−n = 0 by either adding or subtracting atoms. Clearly, such compensation costs finite energy. Similarly, no transfer of k atoms from k-CSF to n-CSF is possible without spending energy. In contrast, a transfer of a group of k quasi-molecular complexes each made of n atoms from n-CSF to k-CSF costs no energy and, therefore, can be coherent. [Conversely, a group of n complexes each built out of k atoms can be coherently transferred from k-CSF to n-CSF].
If p > 1 is the greatest common divisor, the minimal coupling can be achieved between the products (Ψ n ) k/p and (Ψ k ) n/p with the Josephson energy In this case, the coherent transport proceeds by increments of kn/p bosons. More specifically, k/p groups of n bosons are transferred from n-CSF to k-CSF and n/p groups of k bosons are transferred from k-CSF to n-CSF.
The most direct implication of the above analysis is an unexpected quantization of the rotational flux. We consider it below for the case of k-CSF and n-CSF quantum liquids confined to a ring geometry and separated by two weak links.

Ring with k-CSF n-CSF junctions and velocity quantization
In a simple superfluid of bosons each of mass m, the quantization of the velocity circulation is characterized by the quantum of circulation Φ 0 = h/m, where m is the bare atomic mass. This relation follows from the irrotational nature of superfluidity, the phase φ of the superfluid being defined modulo 2π, and the requirement of the Galilean invariance. Accordingly, in the k-CSF phase (where the coherent transport occurs by groups of k bosons) the circulation quantum is given by Φ k = h/(mk) and Φ n = h/(mn) in the n-CSF phase. In the following we address how the quantization changes if there are two composite superfluids confined in the ring geometry and connected by weak links using two different approaches.
Consider a ring, Fig. 4, of radius r containing two composite superfluids k-CSF and n-CSF and rotating at angular velocity Ω around the axis through the ring center and perpendicular to the ring plane. The Josephson energy E J of the system as a function of the difference of the phases ϕ k and ϕ n is the sum of the energies written for two (identical) junctions Without loss of generality we can set ϕ k (a) = 0 and ϕ n (b) = 0. Then, integrating the velocity along the ring, we find I n = ϕ n (a)/(mn) along the n-CSF arc and I k = ϕ k (b)/(mk) along the k-CSF arc. If the arc-lengths L are large, the energy of the currents along the rim is ∝ 1/L and can be ignored and one should only consider the Josephson energy given by Eq. (20). This energy is periodic in the phases ϕ k (b) and ϕ n (a) with the corresponding periods 2pπ/n and 2pπ/k. Such a periodicity is the analog to the phase φ of a simple superfluid being defined modulo 2π. Thus, the quantum of the circulation along the ring can be written as Φ nk = I k + I n = hp/(mnk) = Φ 0 p/(nk). This result implies that the coherent current circulating along the ring can be induced at special values of the angular frequency with the period δΩ = p nk mr 2 .
In general, lowering the period leads to an increased sensitivity to rotations. A similar situation occurs in the case when n = k, when the maximum divisor is p = n = k and the factor p/(kn) becomes simply 1/n. Using two composite superfluids with, for example, k = n + 1 results in the reduction of the period by the factor 1/n(n + 1).

Summary and outlook
In summary, we considered a system of dipolar bosons trapped in a stack of onedimensional optical lattice layers (tubes) with all tubes lying on the same plane. The direction of polarization is perpendicular to the tubes and can lie in the plane or perpendicular to it. Within this geometry, in the first case the dipolar interactions between bosons are attractive in different tubes and repulsive within each tube. In the second case, the interactions are purely repulsive. Based on previous results [1], we provided a description of the various quantum phases stabilized in this system depending on interactions and the filling factor in each layer. These phases include independent superfluids, composite superfluids, composite solids and supercounterfluids.
Composite superfluids and supercounterfluids are found to be threshold-less with respect to the dipolar interaction strength as derived by a combination of bosonization and renormalization group techniques. We described in detail the bosonization formalism for the case of two tubes and verified its main prediction -the absence of the threshold for the formation of paired superfluid -by ab initio QMC simulations. We have also discussed how bosonization predicts that the result obtained for two tubes also holds for arbitrary number of tubes. Moreover, we proposed how the phases stabilized in this multilayer geometry can be used to create generalized Josephson junctions, with applications to higher sensitivity rotation sensors. In particular, utilizing junctions between composite superfluids made of k and n components trapped in ring-shaped optical lattice layers, one can engineer a sensor for rotations analogous to a SQUID. In this setup, the rotational flux quantum is significantly reduced, by a factor p/(kn), where p is the maximal common divisor. Hence, the accuracy of this rotation sensor can be significantly higher than in the cases of simple superfluids or identical composite superfluids.
In this work we set out to elucidate the multi-faceted potential of the composite phases of dipolar bosons with respect to fundamental questions and possible applications in various devices. One of the most interesting features is the non-trivial quantization of circulation of superfluid velocity. We have explored this question for the simplest geometry and we don't exclude that more interesting features can be uncovered for various networks and combinations of the composite phases. Moreover, questions associated with the possibility of practical realization and detection provide another interesting direction for future research.