Spin modulation instabilities and phase separation dynamics in trapped two-component Bose condensates

. In the study of trapped two-component Bose gases, a widely used dynamical protocol is to start from the ground state of a one-component condensate and then switch half the atoms into another hyperﬁne state. The slightly different intra-component and inter-component interactions can then lead to highly non-trivial dynamics, especially in the density mismatch between the two components, commonly referred to as ‘spin’ density. We study and classify the possible subsequent dynamics, over a wide variety of parameters spanned by the trap strength and by the inter-to intra-component interaction ratio. A stability analysis suited to the trapped situation provides us with a framework to explain the various types of dynamics in different regimes.

In a number of two-component BEC experiments reported over more than a decade, a standard technique has been to start from the equilibrium state of a single-component BEC, e.g., populating a single hyperfine state of 87 Rb, and then using a π/2 pulse to switch half the atoms to a different hyperfine state [1][2][3][4][5][6][7][8][9].This results in a binary condensate where the two intra-species interactions (g 11 and g 22 ) and one inter-species interaction (g 12 ) are all slightly different from each other, but the starting state is the ground state determined by g 11 alone.
Since it has been realized several times in several different laboratory setups, this is a paradigm non-equilibrium initial state for binary condensate dynamics.A thorough and general analysis of the dynamics subsequent to such a π/2 pulse is thus clearly important.In this article we present such an analysis, clarifying the combined role of the inter-species interaction (g 12 ) and the strength λ of the trapping potential.We provide a stability analysis mapping out regions of the λ-g 12 parameter space hosting different types of dynamics.Since it is now routine to monitor real-time dynamics in such experiments (e.g.[6]), we also directly analyze the real-time evolution after a π/2 pulse.
It is widely known that the ground state of a uniform twospecies BEC is phase separated or miscible depending on whether or not the inter-species repulsion dominates over the self-repulsions of the two species, i.e., if g 11 g 22 < (g 12 ) 2  (1) then the ground state is phase separated [15].This criterion is also a key ingredient in understanding dynamical features such as pattern dynamics in the density difference between the two species -such "spin patterns" emerge when the phase separation condition is satisfied.This can be understood as the onset of a modulation instability [16][17][18], identified by the appearance of an unstable mode in the excitation spectrum around a reference stationary state.For a homogeneous situation, linear stability analysis shows that modulation instability sets in when the condition of Eq. ( 1) is satisfied [16][17][18] .
The situation is different in the presence of a trapping potential.Phase separation in the ground state, as well as the appearance of modulation instability when starting from a mixed state, now requires larger inter-species repulsion [14,19].This suggests that the region of parameter space where pattern dynamics occurs also depends on the trap.A trap is almost always present in cold-atom experiments, and it is easy to imagine experiments where the trapping potential is not extremely shallow but varies between tight and shallow limits.It is thus necessary to examine the relevance of Eq. (1) for trapped binary BECs.To this end, we explore different trap strengths spanning several orders of magnitude, and identify the appropriate extensions of Eq. ( 1) for the type of spin dynamics resulting from the π/2 protocol described above.
We focus on the effects of two parameters.First, we study effects of changing cross-species interaction g 12 , thus generalizing Eq. ( 1) for trapped situations.Second, we explore the role of the relative strength of the trap with respect to the interactions.Our analysis, performed for a one-dimensional (1D) geometry, sheds light on the situation where g 11 and g 22 are close but unequal: (a) the stability analysis is performed for g 11 = g 22 and their difference serves only to select appropriate instability modes; (b) the simulations are performed with g 22 /g 11 = 1.01.
In Section II, we introduce the formalism and geometry.In Section III, we show results from a linear stability analysis for a sequence of trap strengths, and identify and analyze relevant modulation instabilities.Through an analysis of unstable modes, we present a classification of the parameter space into dynamically distinct regions, in relation to the prototypical initial state explained above.This may be regarded as a dynamical "phase diagram".A remarkable aspect is that the "phase transition" line most relevant to spin pattern dynamics does not arise from the first modulation instability (studied in Ref. [14]).This first instability mode is antisymmetric in space, and as a result is not naturally excited in a symmetric trap with symmetric initial conditions.Complex dynamics (not due to collective modes but rather due to modulation instability) is generated only when the first spatially symmetric mode becomes unstable, which occurs at a higher value of g 12 .
In Section IV we provide a relatively detailed account of the time evolution.For each trap strength λ, for values of g 12 not much larger than g 11 , we observe simple collective modes.Above a threshold value of g 12 , the oscillation amplitude becomes sharply stronger, and at the same time the motion becomes notably aperiodic, signaling that the dynamics is more complex than a combination of a few modes.Dynamical spin patterns start appearing at this stage and become more pronounced as g 12 is increased further.The threshold value at which the dynamics changes sharply corresponds to the second modulation instability line rather than the first, as we demonstrate through careful choice of parameters in each region of the phase diagram derived from stability analysis.Some further connections between the stability analysis and dynamical features, relating to the length scale of generated patterns, appear in Section V.In the concluding Section VI we place our results in context and point out open questions.

II. GEOMETRY AND FORMALISM
The relevant time-resolved experiments have been performed in both quasi-1D geometries (highly elongated traps with strong radial trapping) [6] and in a 3D BEC of cylindrical symmetry with the radial variable playing analogous role as the 1D coordinate [2,5].Since the basic phenomena are very similar, we expect the same theoretical framework to describe the essential features of each case.For definiteness, in this work we show results for 1D geometry.We expect the general picture emerging from this work to be qualitatively true also for other geometries exhibiting the same type of spin dynamics.
We describe the dynamics in the mean field framework at zero temperature, i.e., by two coupled Gross-Pitaevskii equations [30][31][32]: Condensate wave functions ψ 1 (x, t) and ψ 2 (x, t) are normalized to unity, and λ is the strength of the harmonic trap.Factors of particle number and radial trapping frequency are absorbed as appropriate into the effective 1D interaction parameters g ij [6,32,33].We consider purely non-dissipative dynamics, i.e., we do not attempt to model experimental loss rates with a phenomenological dissipative term as done in, e.g., Refs.[5][6][7].
The equations above are in dimensionless form because we measure lengths in units of trap oscillator length and time in units of inverse trapping frequency, for a hypothetical trap of unit strength (λ = 1).The scale for trap strengths is itself fixed by imposing g 11 = 1.With this convention, small values of λ correspond to a BEC in the Thomas-Fermi limit.For comparison, we note that the parameters of the experiment of Ref. [6] corresponds to λ of order 10 −5 in these units.Of course one can switch between different units via the transformation: x → x/l, t → t/l 2 , λ → λl 2 , g → gl, and ψ → ψ √ l, where l is an appropriately chosen scale.
The initial state after a π/2 pulse involves both components occupying the ground state of a single-component system of interaction 2g 11 , because the atoms were all in the first hyperfine state before the pulse.We model this initial situation as a two-component BEC with g 11 = g 22 = g 12 .The π/2 pulse may then be regarded as a sudden change (a quantum quench [29]) of the interaction parameters g 22 and g 12 .
We use g 11 = g 22 for the stability analysis of Section III.For the explicit time evolution reported in Section IV, we use g 11 and g 22 values close but unequal: g 11 = 1, g 22 = 1.01.This choice of close values is convenient for illustrating the structure of the phase diagram, especially for shallow traps.In rubidium experiments the difference between g 11 and g 22 is somewhat larger (in the common case using 87 Rb hyperfine states Numerical simulations presented in Section IV were performed using a semi-implicit Crank-Nicolson method [34,35].

III. STABILITY ANALYSIS AND DYNAMICAL "PHASE DIAGRAM"
We provide in this section a stability analysis for g 11 = g 22 that maps out the regions of λ-g 12 parameter space which support pattern formation instabilities.
Ideally, one might like to perform a stability analysis around the initial state.However, in contrast to the homogeneous case [16], we are faced with the situation that the initial state is not a stationary state of the final Hamiltonian.The choice of reference state is therefore a somewhat subtle aspect of the present analysis.
We use as reference state ψ 0 (x) the lowest-energy spatially symmetric stationary state of the case g 11 = g 22 , with parameter g 12 set to its final value.(For large g 12 , this is not the ground state for these parameters, which is phaseseparated.)This reference state has the advantage of looking relatively similar to our actual initial state (two components totally overlapping in space), and of being a stationary state of the Hamiltonian for which we analyze linear stability.Our reference state can be regarded as placing both components in the single-component ground state for interaction g 11 + g 12 .We are not aware of a suitable stationary state even more similar to the actual initial state.We will see that our stability analysis around this reference state will predict remarkably well the main observed time-evolution features described in Section IV.
Note that it is not natural to use g 11 = g 22 , because stationary states for such a case typically do not overlap completely in space.Instead, in our approach the difference between g 11 and g 22 will play the important role of selecting certain instability modes over others.For this reason, inferences from the present analysis apply only to small relative differences between g 11 and g 22 .
We linearize Eqs. ( 2) and (3) around the reference stationary state ψ 0 (x): where µ is the chemical potential corresponding to the reference state.By keeping only terms of the first order in δψ 1 (x, t) and δψ 2 (x, t), we obtain a system of linear equations which can be cast in the form: Here M is a matrix differential operator which, upon discretization or upon expansion in a set of orthogonal functions, becomes the so-called stability matrix (e.g., [22,24]).We analyze below the eigenmodes of the stability matrix, which we have obtained by numerically calculating the reference stationary state ψ 0 (x) and expanding in the basis of harmonic trap (non-interacting) eigenstates.Since we use g 11 = g 22 for the stability analysis, eigenmodes will have well-defined "species parity", i.e. will all be either even [δψ 1 (x, t) = δψ 2 (x, t)] or odd [δψ 1 (x, t) = −δψ 2 (x, t)] with respect to the interchange of species.Even  x x , while in-out-in modes are persistent only above the higher g s 12 line.The spatially symmetric instability (g s 12 line) is the one relevant for experimental situations with symmetric traps.Squares mark values used in the dynamical simulations of Figs. 3 and 4 (Table 1).modes describe in-phase motion of the two components and simply correspond to the excitation spectrum of a singlecomponent BEC with interaction constant g 11 + g 12 .Odd modes are more interesting -they describe out-of-phase motion of two components and are therefore reflected in the spin dynamics.Additionally, due to the spatial inversion symmetry x → −x, the solutions will also have well-defined spatial parity, and we can distinguish spatially symmetric and antisymmetric modes.
Typical eigenspectra are presented in Fig. 1.In the case of a tight trap λ = 0.2, we notice two modes whose frequencies are nearly constant.These are even modes encoding single-component or in-phase physics.The lower one is the dipole (Kohn) mode with frequency equal to the trap frequency λ.The second nearly-constant mode is the breathing mode, which for elongated traps takes value close to ω 2 = 3λ 2 .The breathing mode (oscillations of cloud size) is visible in the plots of Fig. 3 (Section IV) as a fast oscillation of the total condensate widths.
The two lowest-lying eigenmodes are odd modes encoding out-of-phase physics.For g 12 1, their frequencies are significantly below the breathing mode, and therefore lead to relatively slow oscillations in the spin density.This will also be visible in the real-time dynamics presented in Section IV (first two columns of Figs. 3 and 4).The forms of the corresponding eigenvectors are shown in the lower right of Fig. 1.The nature of the eigenvectors shows that the motion related to the lowest mode corresponds to the left-right oscillations of the two species, while the next odd mode corresponds to spatially symmetric spin motion.The frequencies of these two modes become imaginary at certain values of g 12 , thus leading to the onset of modulation instabilities.The antisymmetric mode becomes unstable at smaller value of g 12 (g a 12 ≈ 1.6 for λ = 0.2) in comparison to the symmetric mode (g s 12 ≈ 2.4 for λ = 0.2).In a spatially symmetric trap, there is no natural mechanism for exciting the spatially antisymmetric mode.On the other hand, any difference between g 11 and g 22 naturally excites the second (spatially symmetric) mode.Thus, the second mode, occurring at larger g 12 , is the relevant instability for understanding the dynamics observed in experiments and explored numerically in Section IV.
We find similar excitation spectra for trap strengths λ spanning several orders of magnitude.The spatially antisymmetric mode becomes unstable before the spatially symmetric mode, and both instabilities get closer to 1 as the trap gets shallower.For example, for λ = 10 −3 (also shown in Fig. 1) the lowest instability sets in for g a 12 ≈ 1.02, while the next one appears at g s 12 ≈ 1.05.The distinction between two instabilities becomes ever smaller as we go toward a uniform system λ → 0, where the phase-separation condition Eq. ( 1) becomes exact.Nevertheless, even for shallow traps, the issue is not purely academic as the precision in experimental measurement and control of scattering lengths continues to improve [6,36].
In Fig. 2 (main panel), the results of the stability analyses are combined to present a dynamical "phase diagram".The two lines show the two instabilities (g a 12 and g s 12 ) as a function of trap strength λ.For very shallow traps, the two transition lines merge as g s 12 ≈ g a 12 ≈ 1.The lower transition line (g a 12 ) was previously introduced in Ref. [14].However, for a trap and initial state with left-right spatial symmetry, this is not the relevant dynamical transition line, because the first even mode only becomes unstable at some higher g 12 value, given by the g s 12 line.
In the next Section, we will see that spin pattern dynamics is indeed only generated when the inter-component repulsion g 12 exceeds the second instability line (g 12 > g s 12 ), and that crossing the first instability (g a 12 < g 12 < g s 12 ) is not enough for pattern formation in a spatially symmetric trap.

IV. DYNAMICAL FEATURES ACROSS THE PARAMETER SPACE
In this Section we present and analyze the dynamics obtained from direct numerical simulation of the Gross-Pitaevskii equations ( 2) and (3), after the system is initially prepared in the ground state of the situation g 11 = g 22 = g 12 = 1.The subsequent dynamics is performed with g 11 = 1, g 22 = 1.01, and several different values of g 12 for each trap strength λ.
It is difficult to show the full richness of pattern dynamics through plots of a few quantities.We choose to show the dynamics through two types of plots (Figs. 3 and 4).Fig. 3 shows the time dependence of the root mean square widths of the two components while Fig. 4  corresponds to a different trap strength (λ), and we approach the shallow trap (Thomas-Fermi) limit going from top to bottom.
For each λ the four values of g 12 from Table I are used for Figs. 3 and 4. We have chosen g 12 values such that the first panel in each row is in the parameter region where there are no instabilities, the second one is in the region where the only instability is the antisymmetric one, and the third on each row is at g 12 values just above the second, relevant, instability.The fourth panel on each row is at higher g 12 values.The choice of g 12 values with respect to instability lines is clear in the tighter traps of the top three rows, as also shown by squares in Fig. 2. For shallow traps (lower rows), the instability lines are too close together and too close to g 12 = 1, so making such choices is not meaningful.In the following, as we compare features of the different columns, we implicitly exclude the lowest row (smallest λ).This is also indicated by the fact that the schematic instability lines in Figs. 3 and 4 are not extended to the lowest row.
Broadly speaking, we note that there is only regular (collective-mode) dynamics in the second-column figures (g a 12 < g 12 < g s 12 ) even though an instability is present.There is generally a sharp difference between the second and third figure in each row, indicating that the second instability (g s 12 ) is the relevant one.The fourth panel on each row is at higher g 12 values, showing more rich dynamics.
In Fig. 3, we show time-dependence of the individual widths (w 1 , w 2 ) and also of the total root mean square width, w(t) = (w 2 1 (t) + w 2 2 (t))/2.Consistent with our observation that spatially symmetric modes (and not the antisymmetric ones) are naturally excited in the current setup, the dynamics shows signatures of the two most prominent spatially symmetric modes noted in Fig. 1.The breathing mode is the easiest to notice and most ubiquitous -it shows up in almost every parameter choice as oscillations in the total density (inphase in the two components), with a typical period given by 2π/ √ 3λ ≈ 3.63/λ.This follows from the frequency of this mode being almost constant near √ 3λ.We also see out-of-phase motion of the two components, associated with the lower spatially-symmetric mode in Fig. 1, which has odd species parity.In the first two columns of Fig. 3, corresponding to smaller values of g 12 such that this mode has small real frequencies, this is excited as a regular   2. While the first two columns look qualitatively the same and show regular oscillatory dynamics, in the third column we observe aperiodic motion of stronger amplitude that we relate to the onset of spin pattern dynamics.The spin dynamics is even more pronounced in the fourth column.
'spin' mode.For example, at λ = 10 −3 and g 12 = 1.04, we observe an out-of-phase oscillation with the period of approximately ≈ 30, much slower than the breathing mode.In addition, the oscillation period of the out-of-phase motion is slower in the second than in the first column of each row, corresponding to the decreasing frequency of the mode, as seen in stability analysis (Fig. 1).Once g 12 becomes large enough that the instability threshold for this mode is crossed, the oscillation amplitudes increase sharply and the width dynamics becomes strongly aperiodic and irregular (third column of Fig. 3).This signifies the onset of pattern dynamics, as opposed to the excitation of a regular collective mode around a stable state.Irregularity of the width dynamics at stronger g 12 is even more apparent in the fourth column of Fig. 3.
It is noteworthy that the spatially antisymmetric modes play no role and do not show up in these dynamical simulations.We see no signature of the Kohn mode.Nor do we see any sharp change associated with the instability of the antisymmetric mode, i.e., there is no sharp difference between the first two columns of Fig. 3.
In Fig. 4 we show the dynamics of the "spin density" The case of very shallow traps (last row), resembles the data in Refs.[17,18].As in Fig. 3, the first two columns show regular oscillations, corresponding to collective modes without instability.A sharp change occurs, not across the first instability line (between 1st and 2nd column), but instead across the second instability line (2nd and 3rd columns), especially for tighter traps (top three rows) where comparison with instability lines is meaningful.The sharp change can be noted through the color scales, which is dramatically different between second and third columns in the upper rows.

V. LENGTH SCALES OF PATTERNS
In homogeneous stability analysis, the length scale of patterns is inferred from the wavevector (momentum) at which an instability first occurs.Since we perform our stability analysis specifically for trapped systems, we do not have a momentum quantum number.Nevertheless, the eigenvectors of the unstable modes contain information about the form of patterns generated in the dynamics of the trapped system.This is illustrated in Fig. 5, where the eigenvectors of the lowest unstable even-parity modes are shown for several values of g 12 , together with the spin patterns generated in the non-equilibrium dynamics.There is a close match between the distance between nodes of the eigenvectors (rough analog of 'wavevector') and the length scales involved in the patterns.
In Fig. 4 we see that the patterns contain more spatial structure in shallow traps.The top two rows (tight traps) only show in-out-in type of patterns.This can be understood from the idea that the interactions induce length scales ('healing lengths') in the problem, which are smaller for larger interactions, and which set the length scale of spatial structures.For tight traps, the healing length set by the interactions is large or comparable to the cloud size, so that only global dynamical patterns are generated.In such traps, generation of complex patterns with many spatial oscillations would require much higher values of (g 12 /g 11 − 1).For shallow traps, the healing length becomes much smaller than the cloud size; as a result one can have a multitude of dynamical spin structures in the system, of the type seen in experiments and prior simulations [5,6,18].This heuristic explanation can be made quantitative by counting the number of nodes appearing in the eigenmodes (as in Fig. 5).

VI. CONCLUSIONS AND OPEN PROBLEMS
In this article we have analyzed a widely used dynamical protocol for two-component BECs, which involves starting from the ground state of one component and switching half the atoms to a different component through a π/2 pulse.We have presented a stability analysis suitable to the trapped situation, and also presented results from extensive dynamical simulations.Through an analysis of unstable modes, we have presented a classification of the parameter space into a number of dynamically distinct regions, in relation to the prototypical initial state.This may be regarded as a dynamical "phase diagram".
In the 'stable' regime of parameter space (no modulation instabilities), our stability analysis explains the observed slow spin oscillations compared to the fast breathing mode oscillations of the total density.We demonstrate that the important "phase transition" line for spatially symmetric situations relevant to most experiments is not the first instability (studied in Ref. [14]), but a second transition line.The first instability is antisymmetric in space, and as a result is not naturally excited in a symmetric trap.
Our stability analysis is performed relative to a stationary state of the situation g 11 = g 22 .The π/2 pulse of the experiments (in the cases where g 11 =g 22 ) can be considered as turning on a nonzero (g 11 − g 22 ), i.e., turning on 'buoyancy' such that one component gains more energy by being in the interior of the trap compared to the other.This helps to select instability modes which are symmetric in space.
Since we have used a stability analysis with g 11 = g 22 to analyze dynamics with g 11 =g 22 , an obvious question is how the ratio g 22 /g 11 affects the regime of applicability of this scheme.We expect that features of this (g 11 = g 22 ) stability analysis are useful for dynamical predictions as long as g 12 /g 11 −1 is roughly more than g 22 /g 11 −1.For example, for shallow traps (small λ), the instabilities occur at g 12 /g 11 − 1 values comparable to 0.01, which is why the placement of parameters in the three dynamical regions of the 'phase diagram' (Fig. 2) is not meaningful for the smallest λ values (lowest rows of Figs. 3 and 4).
For the stability analysis we used a reference stationary state which is of course not the initial state: the initial state is the ground state for g 11 = g 22 = g 12 , while the reference state is the lowest-energy spatially symmetric stationary state corresponding to the final value of g 12 .The instability lines found in this stability analysis would describe even better a situation where the dynamics is triggered by a small quench of g 12 , as opposed to the changes of g 12 that we consider here, which can be relatively large.We have looked at some examples of this type of dynamics and indeed find instabilities matching the stability analysis extremely well.However, although the initial state in the π/2 dynamics is somewhat different from the reference state of our stability analysis, our results show that this stability analysis does provide an excellent overall picture of the dynamics generated by the π/2 protocol.
Our work opens up a number of questions deserving further study.First of all, we have thoroughly explored the λ-g 12 parameter space, while assuming that the intra-component interactions g 11 and g 22 are unequal but close in value.The regime of large difference (g 11 − g 22 ) clearly might have other interesting dynamical features which are yet to be explored.
Second, in this work we have restricted ourselves to the mean field regime.While the mean field description captures well the richness of pattern formation phenomena (c.f.Refs.[5,14,17,18,20] in addition to present work), it may be worth asking whether quantum effects beyond mean field might have interesting consequences for the patter dynamics generated by a π/2 pulse.For bosons in elongated traps, regimes other than mean field (such as Lieb-Liniger or Tonks regimes) may occur naturally in experiments [37][38][39][40][41]. Dynamics subsequent to a π/2 pulse in strongly interacting 1D gases outside the mean field regime is an open area of investigation.
Third, we have assumed a spatially symmetric trap and an initial condition with spatial symmetry, and this plays a crucial role in the selection of instability channels.In a reallife experiment, the trap will have some left-right asymmetry.Also, thermal and quantum fluctuations can initiate spatially antisymmetric excitations.The extent to which a small spatial asymmetry affects spin dynamics remains unexplored; in such a case we would have some type of competition between two types of instabilities.Ref. [14] has studied dynamical effects of fluctuations (noise), but the effects of thermal and quantum fluctuations is yet to be studied in the context of a π/2 protocol.
Finally, one could consider time evolution and spatiotemporal patterns generated by a π/2 pulse in the presence of an optical lattice, described by the dynamics of a two-component Bose-Hubbard model.This is a situation easy to imagine realizing experimentally.One could speculate complex interplay between spin dynamics and the spatial arrangement of Mott and superfluid regions.

FIG. 1 .
FIG. 1. (Color online.)Results from stability analysis.Squared eigenvalues ω 2 of the stability matrix M are plotted against g12, for a tight trap (top) and for a shallow trap (bottom left).The arrows show the values of g12 for onset of the two instabilities, namely g a 12 (onset of spatially antisymmetric modulation instability) and g s 12 (onset of spatially symmetric instability).Typical eigenvectors corresponding to these two modes are shown in the panels on lower right.
1 ); however our insights should be relevant to a broad regime of possible experiments.A full exploration of the regime of arbitrary differences (g 11 − g 22 ) remains an open task, beyond the scope of the present manuscript.

FIG. 2 .
FIG. 2. (Color online.)The dynamical "phase diagram" showing the critical values of g12 for the onset of the two types of modulation instability versus the trap strength λ.The instability lines are shown as straight lines joining numerically determined values.The oscillation schematics on the right (and corresponding arrows) indicate that left-right-left oscillation modes are persistent everywhere above the g a12 line, while in-out-in modes are persistent only above the higher g s 12 line.The spatially symmetric instability (g s 12 line) is the one relevant for experimental situations with symmetric traps.Squares mark values used in the dynamical simulations of Figs.3 and 4(Table1).

5 FIG. 4
FIG. 4. (Color online.)Spin dynamics subsequent to the π/2 protocol, shown via the density difference λ −1/2 |ψ1(x, t)| 2 − |ψ2(x, t)| 2 .Traps and g12 values are the same as in Fig. 3 and Table I: λ decreases from 10 −1 to 10 −5 from top to bottom and g12 values are indicated near top of each panel.As in Fig. 3, the black solid line and the red dashed line indicate the instability lines from the "phase diagram" of Fig. 2. Note the sharp change of color-scale ranges between second and third columns in the upper rows, indicating that the dynamics changes dramatically only across the second instability line.
is the black solid line intermediate between the other two.From top to bottom: tight to shallow traps.For each trap strength, four values of g12 (indicated near top of each panel) from Table I are used.The two lines separating first and second column (red dashed) and second and third column (black solid) indicate the 'positions' of instability lines, from Figure FIG.3.(Color online.)Time evolution of root-mean-square widths after π/2 pulse (interaction quench).First component width w1(t) is shown as blue dashed line, second component width w2(t) is shown as red solid line (gray solid without color), the total width w(t) = (w 2 1 (t) + w 2 2 (t))/2