Thermal-Condensate Collisional Effects on Atomic Josephson Junction Dynamics

1. Introduction

A Josephson junction consists of a thin barrier or insulator separating two weakly-coupled superfluids [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18] or superconductors [19,20,21]. Firstly studied with superconductors [19,20,21], it has been shown that a current flows across the junction due to tunneling of Cooper pairs from one side to the other of the junction. This current could be direct or alternating in the absence or presence of applied external potential, respectively. Josephson effects have been studied extensively, both theoretically and experimentally, with superfluid helium (e.g. [1,2,3,4]), trapped ultracold atomic gases [5,6,7,8,9,10,11,12,13,14,15,16,17,18,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45] (see also [46,47,48,49,50,51,52,53]), and exciton-polariton condensates [54,55]. In particular, one way to induce the Josephson dynamics is by imposing an initial population imbalance between the two superfluids, which induces an initial acceleration of the superfluid across the junction. Previous studies at $T=0$ (e.g. [7,15,23,25,37,40]) have revealed the existence of different main dynamical regimes in an atomic Josephson junction: the Josephson plasma, the dissipative and the self-trapping regime. While individual works had studied transitions between two of those regimes, the 3 dynamical regimes were first clearly unified into a single Josephson phase diagram in our previous work [38]. The first and the latter regime are characterized by coherent dynamics while in the second one, phase slips, well-known from condensed matter systems (see, e.g. [1,56]), lead to the generation of topological defects (the nature of which is set by the underlying system geometry [36,37,38,57]) and sound waves, with the net effect of inducing dissipation of the Josephson oscillations.

The aim of the present work is to investigate the critical effect that collisions have on these 3 dynamical regimes when a condensate is coupled to a dynamical thermal cloud. Our study focuses on an elongated three-dimensional geometry, inspired by the experiment of [14,15], where the phase-slip was induced by the generation of vortex rings. In the present work we will refer to such dissipative regime as the vortex-induced dissipative regime, in order to distinguish such behaviour from the dissipation of the otherwise coherent dynamics arising purely from thermal effects.

To address this, we will use a fully self-consistent model, known as the Zaremba-Nikuni-Griffin, or `ZNG’, model [58,59,60]: this describes the atomic dynamics in terms of a dissipative Gross-Pitaevskii equation for the condensate, self-consistently coupled to a quantum Boltzmann equation governing the dynamics of the thermal cloud. Such description includes the full dynamical coupling between the condensate and thermal cloud, e.g. due to their coupled motion, but additionally includes the important effect of collisions which can either perturb the distribution of particles within the thermal cloud, or also transfer particles between the condensate and the thermal cloud. Beyond demonstrating the relevance of our studies [37] to experiments, and identifying the full Josephson junction dynamical regime phase diagram at $T=0$ [38], our previous work [39] has carefully characterised the effect of the thermal cloud dynamics on the damping and frequency shift of the main oscillation frequencies of the condensate, thermal and total population imbalance within the Josephson plasma and vortex-induced dissipation regimes: such prior study has been conducted within a simplified form of the ZNG model, in which only the self-consistent dynamical coupling of the condensate and thermal cloud was included, but all perturbing collisional effects were ignored: in other words, this corresponded to a Gross-Pitaevskii equation coupled to a collisionless Boltzmann equation. While this should in principle be a good approximation at very low temperatures, and for relatively short times, collisions are known to affect the overall system dynamics, with their effect increasing with increasing temperature [61,62].

In the present work we extend our earlier analysis in two ways: firstly, we revisit the results of [39], in order to investigate the role of such collisional processes on the system dynamics; secondly, we extend such analysis to the Macroscopic Quantum Self-Trapping regime (not previously studied within such formalism) presenting results for both collisionless and collisional implementations of the model. As such this paper presents a unified view of the role of thermal dissipation across all 3 distinct Josephson Junction dynamical regimes: Josephson Plasma (JP), Vortex-Induced Dissipative (VID) and Macroscopic Quantum Self-Trapping (MQST) regimes.

We note that the finite-temperature decay of the self-trapping regime (already experimentally observed in [11]) has also been previously studied in Ref. [34] in the context of an alternative model, known as the stochastic projected Gross-Pitaevskii equation [63]. Rather than making a clear distinction between condensate and thermal cloud as in ZNG, such model treats the low-lying modes of the system cumulatively as a classical field, while simultaneously ignoring the dynamics of high-lying modes, which are treated as a reservoir yielding a dissipative term and a stochastic noise source to the low-lying modes of the system. The relation between such models has been discussed, e.g. , in [60,64,65,66]. Importantly for our present work (and leaving aside the comparison of such models) we note that in carrying out their analysis, the work of Ref. [34] fixes the total particle number (up to the imposed energy cutoff) while varying temperature. Given that Josephson junction dynamics (e.g. frequencies) depend on the condensate number, whereas such approximation changes the condensate number with temperature, such work cannot allow for a systematic comparison of the frequency shifts at fixed condensate number.

Our present study focusses on the effects of the collisions on both the condensate and thermal imbalance dynamics, in all three dynamical regimes and for a temperature range $T\le 0.67T_c$, which is sufficiently far from the critical point. The paper is structured as follows. In Sect. Section 2 we describe the numerical methods employed for this study: in Sect. Section 3 we show the role of thermal cloud on different dynamical regimes, while in Sect. Section 4 we describe our main conclusions.

2. Methods

We model the system at finite temperature by numerically solving the full Zaremba-Nikuni-Griffin (ZNG) formalism [58,59,60,64], explicitly including self-consistently obtained collisional terms. This model has been successfully implemented for the study of topological defect dynamics [65,67,68,69,70], the study of first and second sound in highly elongated trap [71], the study of evaporative cooling [72], collective modes [61,73,74,75], and mixtures [62]. It consists of separating the system into two components: the condensate wavefunction $\psi (\mathbf{r},t)$ and the phase-space distribution function $f(\mathbf{r},\mathbf{p},t)$ representing the thermal cloud. As a result, two coupled equations are solved self-consistently; the dissipative Gross-Pitevskii equation (DGPE) and the Quantum Boltzmann equation (QBE) describing the condensate and thermal cloud dynamics, respectively. The time-dependent DGPE can be written as

$$i\hslash {\displaystyle \frac{\partial \psi (\mathbf{r},t)}{\partial t}}=\left[-{\displaystyle \frac{{\hslash }^2}{2M}}{\nabla }^2+V_{\mathrm{ext}}\left(\mathbf{r}\right)+g(n_{\mathrm{BEC}}(\mathbf{r},t)+2n_{\mathrm{th}}(\mathbf{r},t))-iR(\mathbf{r},t)\right]\psi (\mathbf{r},t)$$

(1)

where M is the particle mass (here ${}^6\mathrm{Li}$ molecule), $V_{\mathrm{ext}}$ the external trapping potential, $g=4\pi {\hslash }^{2}a/M$ the interaction strength with a the s-wave scattering length between molecules, $n^{\mathrm{th}}$ the thermal cloud density and $n^{\mathrm{BEC}}={\left|\psi \right|}^2$ the condensate density. The rightmost term within the bracket, $-iR(\mathbf{r},t)$, represents the dissipative or source (or ‘growth’) which facilitates the condensate number to change in time due to particle exhange with the thermal cloud.

The corresponding equation for the thermal cloud is the Quantum Boltzmann equation for the Wigner phase-space distribution function f in the form:

$${\displaystyle \frac{\partial f}{\partial t}}+\mathbf{v}·{\nabla }_{\mathbf{r}}f-\left({\nabla }_{\mathbf{r}}{V}_{\mathrm{eff}}^{\mathrm{th}}\right)·\left({\nabla }_{\mathbf{p}}f\right)={\mathrm{C}}_{12}[f,\psi ]+{\mathrm{C}}_{22}\left[f\right]$$

(2)

where $V_{\mathrm{eff}}^{\mathrm{th}}=V_{\mathrm{ext}}+2g(n^{\mathrm{BEC}}+n^{\mathrm{th}})$ is the generalized mean-field potential felt by the thermal particles and $n^{\mathrm{th}}$ the thermal cloud density obtained through the phase-space distribution f as $n^{\mathrm{th}}(\mathbf{r},t)=1/{(2\pi \hslash )}^3\int d\mathbf{p}f(\mathbf{p},\mathbf{r},t)$). Here, the first collisional term includes the collisions between the thermal and condensate particles which cause particle transfer between the two subsystems, while the second term includes collisions between the thermal particles which re-distribute the thermal atoms between the single-particle states. In the collisional limit considered here, even though the total number is numerically well conserved, the condensate number ($N_{BEC}\left(t\right)$) instead changes in time due to the particle exchange between the condensate and the thermal cloud. This is taken into account in the $C_{12}$ collisional integral in the Quantum Boltzmann equation; the latter is related to the source term $-iR$ in the dissipative GPE equation via $R(\mathbf{r},t)=(\hslash /2n_{\mathrm{BEC}})\int d\mathbf{p}/{(2\pi \hslash )}^3{C}_{12}(\mathbf{p},\mathbf{r},t)$. As a consequence, the thermal cloud particle number, $N^{th}\left(t\right)$, also changes in time.

The geometry considered in our study is an anisotropic elongated 3D harmonic trap , with a double-well potential of the form

$$\begin{array}{ccc}\hfill V_{\mathrm{ext}}(x,y,z)& =& {\displaystyle \frac{1}{2}}M\left({{\omega }_x}^2{x}^2+{{\omega }_y}^2{y}^2+{{\omega }_z}^2{z}^2\right)\hfill \\ & +& V_0{e}^{-2x^2/w^2}.\hfill \end{array}$$

(3)

The trap frequency values are ${\omega }_x=2\pi \times 15$ Hz, ${\omega }_y=2\pi \times 148$ Hz, and ${\omega }_z=2\pi \times 187.5$ Hz along the x, y and z directions, respectively, chosen to coincide with the experimental values of Ref. [14,15]. The parameters $V_0$ and w are the Gaussian barrier height and $1/e^2$ width, respectively. Our numerical studies are conducted in a grid of $\left[-48,48\right]l_x,\left[-8,8\right]l_x,\left[-8,8\right]l_x$ along the x, y and z directions respectively, where $l_x=\sqrt{\hslash /M{\omega }_x}$, based on $2048\times 128\times 128$ grid points for the thermal cloud. For the less spatially extended condensate, we use a corresponding half grid of size $\left[-24,24\right]l_x$, with 1024 grid points along the x axis. We initially find the equilibrium state by adding the tilted linear potential $-ϵx$ to the $V_{\mathrm{ext}}(x,y,z)$. A typical 2D integrated equilibrium profile is shown in Figure 1(a) for the condensate (i) and the thermal cloud (ii). As evident in this plot, the presence of the repulsive interaction between the condensate and the thermal cloud leads to the thermal particles populating regions of lower condensate density, i.e. at the edges of the condensate, or also in the barrier region.

To initiate the dynamics, we instantaneously set the parameter $ϵ$ to zero, thus creating an initial population imbalance, $z_0$, between the two wells. Our earlier work in the limit of a pure $T=0$ condensate revealed that different dynamical regimes emerge depending on the barrier parameters [38]: for relatively small barrier widths or heights the system transits from Josephson plasma (JP) to vortex-induced dissipative regime (VID) with increasing $z_0$, as soon as it exceeds a critical value $z^c$. In the opposite limit of relatively large barrier height or width the system transits from the Josephson plasma to the macroscopic quantum self-trapping regime (MQST) for $z_0$ beyond the self-trapping critical value $z^{st}$. For our studies here we fix the Gaussian barrier width at $w=3.8\xi$, where $\xi =\hslash /\sqrt{2\mu (T=0)M}=0.52\mu$m denotes the condensate healing length, and change the barrier height. To probe the effect of collisions on these dynamical regimes in the clearest manner, we consider two values for the main barrier height: (i) $V_0=104h{\nu }_x\simeq \mu (T=0)$ for the JP-VID transition, and (ii) $V_0=210h{\nu }_x\simeq 2\mu (T=0)$ for the JP-MQST transition. For each dynamical regime, when comparing the results at different temperatures the barrier height value, $V_0$ is kept fixed. In order to avoid the possibility of a transition to a different dynamical regime while increasing temperature [found to occur for fixed total number in [39], all finite temperature studies presented here are conducted at fixed condensate particle number, $N^{\mathrm{BEC}}\simeq 5\times {10}^4$. This means that at different temperatures, the total particle number takes different values, increasing with increasing temperature, in order to keep $N^{\mathrm{BEC}}$ fixed.

We define the condensate $z_{\mathrm{BEC}}\left(t\right)$, the thermal cloud $z_{\mathrm{th}}\left(t\right)$ and the total $z_{\mathrm{tot}}\left(t\right)$ population imbalances respectively as:

$$\begin{array}{ccc}\hfill z_{\mathrm{BEC}}\left(t\right)& =& {\displaystyle \frac{N_{\mathrm{BEC}}^R\left(t\right)-N_{\mathrm{BEC}}^L\left(t\right)}{N_{\mathrm{BEC}}\left(t\right)}}\hfill \\ \hfill \end{array}$$

(4)

$$\begin{array}{ccc}\hfill z_{\mathrm{th}}\left(t\right)& =& {\displaystyle \frac{N_{\mathrm{th}}^R\left(t\right)-N_{\mathrm{th}}^L\left(t\right)}{N_{\mathrm{th}}\left(t\right)}}\hfill \\ \hfill \end{array}$$

(5)

$$\begin{array}{ccc}\hfill z_{\mathrm{tot}}\left(t\right)& =& {\displaystyle \frac{N^R\left(t\right)-N^L\left(t\right)}{N_{\mathrm{tot}}}}\hfill \end{array}$$

(6)

where $N_{\mathrm{BEC}}^{\mathrm{R}/\mathrm{L}}$ and $N_{\mathrm{th}}^{bamR/L}$ are the number of the condensate and thermal particles on the right/left sides of the barrier, respectively, while $N_{\mathrm{BEC}}$ and $N_{\mathrm{th}}$ are the total thermal and condensate particle number. The presence of a nonzero population imbalance both for the condensate and the thermal particle means that we have a condensate and thermal (or normal) current in the system, with the current extracted as $I=-(N_{\mathrm{BEC}/th}/2)(d{z}_{\mathrm{BEC}/th}/dt)largeba$.

While a single numerical run in the absence of collisions generally accurately models the system dynamics, for a systematic monitoring of the results in the collisional case we generally need to consider a number of (typically five) simulations corresponding to the same physical conditions. Although in the JP and VID regimes the collisional results are found to exhibit negligible variation from run to run, in the MQST dynamics the self-trapping decay time is much more sensitive; in the latter case we thus also perform multiple MQST runs even for the collisionless regime (in addition to also doing so for the collisional regime).

The condensate imbalance data are analyzed by implementing a two-component fit of the form:

$$\begin{array}{cc}\hfill F\left(t\right)& =a_{J}cos(2\pi {\nu }_{J}t+{\varphi }_J)exp(-{\gamma }_{J}t)\hfill \\ \hfill & +a_{i}cos(2\pi {\nu }_{i}t+{\varphi }_i)exp(-{\gamma }_{i}t)\hfill \end{array}$$

(7)

with $a_{J/i}$ and ${\nu }_{J/i}$ being the amplitude and frequency values of each component respectively, while ${\gamma }_{J/i}$ are the damping rates, and ${\varphi }_{J/i}$ the phases of each component.

3. Results

In this section we will discuss the role of the thermal cloud dynamics and, in particular, the importance of collisions in each of the 3 dynamical Josephson junction regimes. This is shown in Figs. 2-4 for 3 different non-zero temperatures, with each figure showing the dynamics of the condensate (left column) and the thermal cloud (right column) without (red lines) and with (black lines) collisions, and for each of the 3 different regimes – from top to bottom showing Josephson Plasma, Vortex Induced Dissipative and Macroscopic Quantum Self-Trapping regimes. For reference, their corresponding $T=0$ dynamics have been previously discussed at length in Ref. [38].

Specifically, at $T=0$, the dynamics in the Josephson plasma regime exhibit coherent oscillations about a zero value with constant amplitude in both the condensate population imbalance (or its time derivative, which gives the current) and the Josephson phase (i.e. the phase difference between the two parts of the condensate). Such oscillations persist at relatively low temperatures $T\ll T_c$ (where $T_c$ is the critical temperature), with only very minor damping and collisions having a negligible effect. This can be seen for the specific case of $T=0.22T_c$ in Figure 2(a,i).

In the vortex-induced dissipative (VID) regime, when the dynamics following the initial population imbalance exceeds the local speed of sound, the system is known to generate (potentially successive) vortices, each accompanied by a $2\pi$ phase jump in the relative phase and some sound emission [37]. In the particular geometry, such vortices appear in the form of vortex rings in the transverse plane: such vortex rings can enter and propagate within the bulk condensate and be subsequently observed experimentally, as discussed in [37]. An image of such a vortex ring can be seen in Figure 1(b) showing both the condensate density containing the vortex ring (i) and the distribution of the thermal cloud (ii) which is maximised within the low-density region of the vortex ring core, and around the condensate. This feature, dominating the dynamics already at $T=0$ [38], significantly perturbs the subsequent condensate motion, followed by significantly reduced amplitude sinusoidal oscillations of the imbalance around a zero mean value. We observe similar behavior at $T=0.22T_c$ of $z_{\mathrm{BEC}}$ (t) in both the collisionless and collisional limits [Figure 2(b,i)], with collisions having no significant effect.

In contrast, the macroscopic quantum self-trapping (MQST) regime at $T=0$ is characterized by small-amplitude oscillations around a nonzero mean value and $2\pi$ phase jumps in the relative phase, but without the appearance of propagating vortices and sound waves [38]. In the collisionless limit, the condensate imbalance continues to oscillate around a nonzero mean value at low temperature, as shown in Figure 2(c,i) (red line), with only a very gradual increase in its oscillation amplitude with time. However, when collisions are included, a transition to a different dynamical regime at around $t=0.5$s for our specific parameters (see Figure 2(c,i), black line): such regime sees the imbalance oscillate around a zero time-average value, resembling the anharmonic oscillations regime characteristic of the transition to self-trapping [25]. Thus, inclusion of collisions in the model are in fact essential for accurately describing the self-trapping regime even at low temperatures.

The MQST regime is highly unstable in the presence of a thermal cloud, and its decay time can vary between runs, even in the collisionless limit. Nevertheless, at $T=0.22T_c$, the imbalance $z_{\mathrm{BEC}}$(t) in all five collisionless runs conducted did not decay within the maximum explored time interval of 0.72 s. In contrast, in the collisional limit, MQST decayed in all runs, even at this low temperature.

Collisions facilitate thermal equilibration between the two wells across all dynamical regimes, as evidenced in Figure 2(a)-(c)(ii). In both Josephson and dissipative regimes – particularly the latter – the thermal imbalance oscillates around a non-zero mean value under collisionless conditions. Collisions suppress this offset, shifting the equilibrium state toward zero thermal imbalance. This contrasts with previous collisionless studies [39], where low-temperature thermal particles relied solely on incoherent tunneling through the barrier. Our results demonstrate that collisional processes, specifically condensate-thermal particle exchange, enhance such equilibration. The MQST regime presents unique observational challenges due to its minimal initial thermal imbalance ($z_0^{\mathrm{th}}\simeq 0.01$). This small signal amplitude creates a high noise-to-signal ratio, rendering thermal imbalance oscillations nearly indistinguishable despite their physical presence. The suppression of measurable oscillations persists even when accounting for collisional effects in this regime.

Next we consider an intermediate temperature $T=70\mathrm{nK}\approx 0.5T_c$ , selected as the limiting value at which the barrier height practically becomes equal to the temperature, i.e. $V_0=104h{\nu }_x\simeq k_{B}T$. Figure 3 presents the time evolution of both the condensate imbalance $z_{\mathrm{BEC}}$ (t) and the thermal imbalance $z_{\mathrm{th}}$ (t) across all three regimes at this temperature. In the Josephson (Figure 3(a)) and dissipative regimes (Figure 3(b)), the $z_{\mathrm{BEC}}$ (t) profiles reveal two key features. First, increasing temperature leads to stronger damping of the condensate imbalance oscillations, consistent with previous findings [39]. In fact, the presence of normal current is expected to damp the Josephson oscillations. Second, collisions further enhance this damping effect, particularly in the dissipative regime, as will be discussed in detail later.

Similarly, collisions significantly reduce the amplitude of $z_{\mathrm{th}}$ (t) oscillations compared to the lower temperature case of $T=0.22T_c$, with the thermal imbalance becoming almost completely damped by $t=0.4$ s in the dissipative regime (see Figure 3(b)-(ii)). The inset in Figure 3(b)-(i) shows a zoomed-in view of the early-time evolution of $z_{\mathrm{BEC}}$(t), where two distinct ‘dips’ or ‘kinks’ appear even in the collisionless limit. These features are signatures of backflow caused by vortex ring generation, as described in [37]. Notably, collisions have little impact on the timing of these dips, which correspond to the vortex ring entering the local Thomas-Fermi surface, and they slightly influence the decay time (i.e. the time when $z_{\mathrm{BEC}}$ (t) assumes zero value). At this higher T in the MQST regime, the condensate imbalance decay time and temporal profile differs from run to run and these differences are even more pronounced in the collisional limit (as will be shown later). Therefore, the particular results shown in Figure 3(c) represent just one characteristic example among a range of possible behaviors.

Figure 4 instead shows the profiles of $z_{\mathrm{BEC}}$ (t) and $z_{\mathrm{th}}$(t) at $T=0.58T_c$, for which $T>V_0$. We have previously shown [39] that such a condition implies the thermal cloud can now pass over the barrier, thus able to execute its own independent oscillations. This new motion can in turn significantly affect the condensate dynamics. In the Josephson plasma regime (Figure 4(a)), collisions cause significantly increased damping of both (i) the condensate imbalance $z_{\mathrm{BEC}}$ (t) and (ii) the thermal imbalance $z_{\mathrm{th}}$(t) compared to the collisionless case. Interestingly, in the dissipative regime (Figure 4(b)), the condensate imbalance starts oscillating with a single dominant frequency, driven by the oscillating thermal component. Thus, in the collisionless limit and in the VID regime, the second component disappears from the condensate imbalance spectrum at temperature $T_1=V_0/k_B$. In this regime, collisions have only a minor effect on $z_{\mathrm{BEC}}$ (Figure 4(b)-(i)), while they strongly influence the thermal imbalance $z_{\mathrm{th}}$(t) (Figure 4(b)-(ii)). The MQST regime is included as an example of possible imbalance profiles, but it is important to note that these profiles can vary substantially from run to run.

3.1. Effect of Collisions on the Frequency and Damping of Dominant Modes in the Josephson Plasma and Vortex-Induced Dissipative Regimes

For each of the Josephson and dissipative regime, we fit the $z_{\mathrm{BEC}}$ (t) and $z_{\mathrm{tot}}$ (t) at each of the five runs with two-component frequency fit and take their mean value as best estimate of our results, with the error bars indicating the maximum deviation.

Consistent with the collisionless results of Ref. [39], the dominant frequencies of both JP and VID at $T\lesssim V_0/k_B\approx 0.5T_c$ correspond to the Josephson plasma frequency ${\nu }_J^{\mathrm{BEC}}$ and the frequency ${\nu }_1^{\mathrm{BEC}}\approx 2{\nu }_J^{\mathrm{BEC}}$ (attributed to second-order term in the tunneling hamiltonian). While the Josephson frequency completely dominates the BEC dynamics in the JP regime, the opposite is true in the VID regime, due to the significant sound energy emitted [38]. At the other extreme, for $T\approx 0.58T_c>V_0/k_B$, the thermal cloud starts exhibiting its own oscillations over the barrier: such motion emerges in a pronounced way in both dynamical regimes, with the sound waves becoming simultaneously significantly suppressed. In such a limit, we previously found [39] that the system dynamics can once again be described by two dominant frequencies, only now these are the Josephson frequency and a frequency associated with the thermal cloud motion, at a frequency approaching the trap frequency.

It is thus natural to ask how this picture is affected by the presence of collisions: The first observation to make here is that the addition of collisions appears to slightly shift the temperature at which the transition between the low and high temperature behaviours emerges to slightly lower values.

The results of the main component frequency value and damping rates as a function of temperature are shown in Figure 5. In particular, Figure 5(a) shows the dominant frequency value (upper subplots) and damping rates (lower subplots) of the condensate and total imbalance in the Josephson ‘plasma’ regime. In both cases the main frequency has a similar value and close to the expected Josephson plasma frequency; this is why it is called ${\nu }_J^{\mathrm{BEC}}$ and ${\nu }_J^{\mathrm{tot}}$. Figure 5(b) instead shows the corresponding plot for the dissipative regime where now the dominant frequency for $0<T\le 0.58T_c$ becomes ${\nu }_1^{\mathrm{BEC}}$ for the condensate imbalance and ${\nu }_1^{\mathrm{tot}}$ for the total imbalance. In the collisionless limit, the increased temperature shifts towards lower values both the condensate and total imbalance dominant frequency while enhancing their damping due to repulsive interaction with thermal component as in Ref. [39].

In the presence of collisions, two main effects are noted. Firstly, we observe that the inclusion of collisions affects differently the dominant frequency in the Josephson and dissipative regime both for the condensate and total imbalance. In fact, in the former case collisions soften the effect of the decrease of ${\nu }_J^{\mathrm{BEC}}$ and ${\nu }_J^{\mathrm{tot}}$ with temperature. For $T<0.5T_c$ the damping rate in the collisional limit has similar value to the collisionless limit but when the thermal cloud starts having more of an effect and exhibiting its own motion, i.e. for $T>0.5Tc=V_0/k_B$, the damping rate ${\gamma }_J^{\mathrm{BEC}}$ and ${\gamma }_J^{\mathrm{tot}}$ becomes larger in the collisional limit. In the dissipative regime instead, the collisions makes the ${\nu }_1^{\mathrm{BEC}}$(i) and ${\nu }_1^{\mathrm{tot}}$ (ii) achieving smaller values with respect to the collisionless limit (for the same $T/T_c$ ) but with a profile similar to the collisionless limit. Moreover, the damping rates ${\gamma }_1^{\mathrm{BEC}}$ and ${\gamma }_1^{\mathrm{tot}}$ are larger in the collisional limit almost the entire range of T.

3.2. Role of Collisions on the Macroscopic Quantum Self-Trapping

Next we consider the thermally-induced dissipation of Macroscopic Quantum Self-Trapping. The presence of the thermal component induces a normal Ohmic current, i.e. a normal conductance [34] in the system by causing dissipation and decay of the condensate imbalance or current oscillations in the self-trapping regime. This decay time can vary significantly from run to run, even in the collisionless limit, as illustrated in Figure 6 for four different temperatures: 28 nK ($0.22T_c$), 62 nK ($0.44T_c$), 70 nK ($0.5T_c$), and 88 nK ($0.58T_c$). At the lowest temperature (28 nK), all runs remain in the self-trapping regime within our probed temporal regime, as shown in Figure 6(a). However, as the temperature increases, the influence of the thermal cloud becomes more pronounced, leading to a decay of the self-trapping regime. This behaviour is evident at higher temperatures, where the imbalance begins to oscillate around zero, indicating a transition to a different dynamical regime.

To quantify this decay process within our simulated results, we extract the decay time from each run and calculate the mean value and associated error bars for both the collisional and collisionless limits, with the results presented in Figure 7. In the collisionless regime, we find that for temperatures below $0.4T_c$, the self-trapping state remains stable throughout the observed time interval (0.72 s). As the temperature increases, the self-trapping regime decays progressively earlier, with the decay time decreasing with temperature. In contrast, in the collisional regime, self-trapping is unstable across the entire temperature range, with a decay time which is shorter than the corresponding time in the collisionless regime. These findings highlight the significant impact of temperature and interactions on the stability of the self-trapping regime.

4. Conclusions

We investigated the impact of thermal-condensate collisions on the dynamical regimes of a bosonic Josephson junction by numerically solving the full Zaremba-Nikuni-Griffin (ZNG) model, which accounts self-consistently for all collisional processes at finite temperatures. We analyzed how temperature influences the frequency and damping of the dominant oscillation modes of both the condensate and total current in the Josephson plasma (JP) and vortex-induced dissipative (VID) regimes, as well as the decay time in the macroscopic quantum self-trapping (MQST) regime. Our results show that the decay time in the MQST regime is highly sensitive to collisions even at low temperatures, decreasing as temperature rises. In contrast, collisions in the JP and VID regimes primarily increase damping of the oscillation frequency. Notably, collisions cause a shift toward higher values of oscillation frequency in these regimes. Overall, our study demonstrates that collisions start playing a crucial role in describing the JP and VID regimes at intermediate to high temperatures, while they are already significant for the MQST regime at low temperatures. It would be interesting to compare our predictions against well-calibrated experimental data in the future.