Multifaceted phase ordering kinetics of an antiferromagnetic spin-1 condensate

TL;DR

This paper investigates the long-term phase domain coarsening in a quasi-one-dimensional spin-1 antiferromagnetic condensate after a magnetic field quench, revealing scaling laws consistent with established models and novel intermediate regimes.

Contribution

It demonstrates the applicability of binary mixture and vector field models to phase ordering kinetics in spin-1 condensates and identifies new intermediate scaling regimes.

Findings

Correlation length growth follows predicted scaling laws.

Distinct regimes for binary mixture and vector field models are identified.

Atypical scaling laws emerge in intermediate regions.

Abstract

We study phase domain coarsening in the long time limit after a quench of magnetic field in a quasi one-dimensional spin-1 antiferromagnetic condensate. We observe that the growth of correlation length obeys scaling laws predicted by the two different models of phase ordering kinetics, namely the binary mixture and vector field. We derive regimes of clear realization for both of them. We demonstrate appearance of atypical scaling laws, which emerge in intermediate regions.

Equations

$$g^{(1)}(x,t)=\Big{\langle}\,{\int{\rm d}x^{\prime}\ \phi(x^{\prime}+x,t)^{*}\phi(x^{\prime},t)\,}\Big{\rangle}.$$

$$H=\int \mathrm{d}x\left[{\vec{\psi }}^{†}\left(-\frac{{\hslash }^2}{2m}{\nabla }^2+q{c}_2\rho f_z^2\right)\vec{\psi }+\frac{c_0}{2}{n}^2+\frac{c_2}{2}{\mathbf{F}}^2\right],$$

$$i\hslash {\dot{\psi }}_{m_F}=$$

$$+q{c}_2\rho |m_F|{\psi }_{m_F}+c_2|m_F|{\psi }_{-m_F}^{\ast }{\psi }_0^2+2c_2{\delta }_{m_F,0}{\psi }_0^{\ast }{\psi }_1{\psi }_{-1},$$

$${\dot{n}}_{m_F}+\frac{\hslash }{m}{v}_{m_F}\nabla n_{m_F}=\frac{\hslash }{m}{n}_{m_F}\nabla v_{m_F}+2\frac{c_2}{\hslash }(|m_F|-2)n_0\sqrt{n_1{n}_{-1}}\sin (\theta )$$

$${\dot{v}}_{m_F}+4\frac{m}{\hslash }{v}_{m_F}\nabla v_{m_F}=2\frac{m}{\hslash }\nabla \left(\frac{{\nabla }^2\sqrt{n_{m_F}}}{\sqrt{n_{m_F}}}\right)+\nabla U_{m_F}-c_2(2-|m_F|)\nabla \left(n_0\sqrt{\frac{n_1{n}_{-1}}{n_{m_F}}}\cos (\theta )\right),$$

$$\frac{dL}{dt}\sim \frac{c_0}{L^2}⟹L(t)\sim t^{1/3}.$$

$${\left(\frac{dL}{dt}\right)}^2\sim \frac{c_0}{L}⟹L(t)\sim t^{2/3}.$$

Full text

Introduction

The theory of phase ordering kinetics (PhOK) states that the growth of order occurs through the coarsening of phase domains, when the system is quenched from disordered to ordered phase. The typical length scale of the phase domain increases then with time. According to the dynamic scaling hypothesis it is due to a global change of scale. In a homogeneous system, the typical size of a phase domain is characterized by the correlation length $L(t)$, which can be defined as the half-width of the equal-time correlation function for the local order parameter $\phi(x,t)$,

[TABLE]

The angle brackets in Eq.(1) indicate the average over initial conditions representing a disordered state. These initial conditions defined in the momentum representation are such $\langle\phi(k,0)^{*}\phi(k^{\prime},0)\rangle=\Delta\,\delta(k-k^{\prime})$, where $\Delta$ sets the size of initial fluctuations in the field $\phi$. The scaling hypothesis implies that due to the existence of the unique characteristic length scale $L(t)$, the correlation function is just one-parameter dependent, i.e. $g^{(1)}(x,t)=f\left(x/L(t)\right)$. Therefore, the central interest of the theory is in the time evolution of the correlation length, and whether this function follows universal scaling laws $L(t)\sim t^{1/z_{d}}$. Further, if yes, what is the value of the dynamical exponent $z_{d}$. The importance attributed to $z_{d}$ stems from its universal character. Since the exponent does not depend on microscopic properties, it reveals general features of an entire class of systems [1].

Kinetic models for the derivation of particular scaling functions were extensively studied and established in 1990s. The models were devoted to classical systems, and the general framework associated with the Hohenberg and Halperin A-J classification of dynamic critical phenomena [2] made the theory more universal. Several systems belonging to the same dynamical class of models exhibit the same universal scaling laws determined by the physical mechanism embedded into them. The phase ordering kinetics, for example, is governed by the H model for binary liquids, in which hydrodynamic processes dominates. On the other hand, the B model associated with vector fields is controlled by diffusion processes mainly. Even though both models are conservative, the scaling exponents corresponding to each of them are different. Nowadays, the subject of the PhOK has been revived in quantum systems, a good example of which are ultra-cold atomic gases [3, 4, 5, 6, 7, 8]. In particular advanced studies concern nonthermal fixed points [9, 10, 11, 12, 13]. The vast majority of works is devoted to two- and three-dimensional spin-1 ferromagnetic condensates [14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27], but also for one-dimensional spinor [28, 29, 30] or binary Bose condensates [31, 32, 33], including driven-dissipative systems [34, 35, 36].

In this paper the PhOK is investigated in a quasi one-dimensional antiferromagnetic spin-1 condensate. The sudden quench of a weak magnetic field leads to the transition from the antiferromagnetic state to a state where domains of atoms with different spin projections separate. We are particularly interested in the superfluid order for atoms in the $m_{F}=0$ Zeeman state. We have performed numerical calculations within the truncated Wigner approximation [37] and made the following observations. The scaling hypothesis holds for any regime of parameters, however, on the longest time scale, the whole variety of scaling exponents is observed. Their values are $3$ or $3/2$ both in the sense of hydrodynamic (H) models [14, 15, 17, 21, 22], and even $4$ (with a logarithmic correction) in the sense of the vector’s field (B) model — another class of the model of the Hohenberg and Halperin classification. All this occurs depending on the system’s parameters. The scaling exponent $z_{d}=4$ has not yet been reported in the PhOK studies with quantum systems.

Using parameters for which we observe a change of scaling exponent passing from the short to long time scales, the characteristic length $L(t)$ can be even a subject to multi-scaling behaviour [38, 39, 30]. The reason for such a multifaceted PhOK is that the spin-1 antiferromagnetic condensate effectively behaves as a binary mixture or vector fields model which is assigned to the H or B model, respectively. Therefore, the PhOK of the system can exhibit features characteristic for the H or B models independently, or both of them simultaneously. A study of the interplay between the models can be made. This is interesting because to date studies of the PhOK attributed system’s behaviour mostly to a particular scaling law characterised by a particular model. A natural questions arises whether distinct models, and hence their physical mechanisms, are mutually compatible, or under what conditions they co-occur. Here, we characterize and classify the appearance of various scaling exponents. We calculate borders for the limit cases in which the B and H models can be realized in their forms. In the region around borders, both models compete, leading to various scaling laws in which scaling exponents smoothly change among the two limiting cases.

I Model and methods

A spinor Bose-Einstein condensate of $N$ sodium atoms is considered [40, 41]. The system is represented by the vector $\vec{\psi}=(\psi_{1},\psi_{0},\psi_{-1})^{T}$, whose components describe atoms in the corresponding Zeeman levels numerated by the magnetic number $m_{F}=0,\pm 1$. We assume a ring-shaped quasi-one-dimensional geometry with periodic boundary conditions [42, 43, 44], where transverse degrees of freedom are confined in a strong potential with frequency $\omega_{\perp}$. The Hamiltonian of the system is

[TABLE]

where $m$ is the atomic mass, $n=\sum n_{m_{F}}=\sum\psi_{m_{F}}^{\dagger}\psi_{m_{F}}$ is the local atom density and ${\bf F}=(\psi^{\dagger}f_{x}\psi,\psi^{\dagger}f_{y}\psi,\psi^{\dagger}f_{z}\psi)$ is the spin density with the spin-1 matrices $f_{x,y,z}$. The spin-independent and spin-dependent interaction coefficients, $c_{0}$ and $c_{2}$, are both positive for sodium atoms. Namely, they are $c_{0}=2\hbar\omega_{\perp}\,(2a_{2}+a_{0})/3$ and $c_{2}=2\hbar\omega_{\perp}\,(a_{2}-a_{0})/3$, where $a_{S}$ is the $s$-wave scattering length for pairs of colliding atoms with total spin $S$ [40]. The term $q\,c_{2}\rho$ is the quadratic Zeeman energy, where the dimensionless parameter $q$ can be controlled using magnetic field or the microwave dressing [40], and $\rho$ is the mean density of the system. The Hamiltonian conserves the total atom number $N=\sum_{m_{F}}N_{m_{F}}$ and the magnetization $M=N_{1}-N_{-1}$. The ground state of the system in the thermodynamic limit [45] is presented in Fig. 1.

In the present study we consider the non-negative magnetization and choose the $2C$ state as an initial state. Next, the quadratic Zeeman shift is set to a fixed value $q>q_{c}$, and the evolution starts. The reason to consider quench from $2C$ towards $2C+\rho_{0}$, $\rho_{+}+\rho_{0}$ phases is that the symmetry breaking occurs in this direction. We describe dynamics of the system on the mean-field level by solving the time-dependent Gross-Pitaevskii (GP) equations

[TABLE]

where $n_{m_{F}}=|\psi_{m_{F}}|^{2}$ and $\delta_{m_{F},0}$ is the Kronecker delta function, see e.g. [47]. To obtain the initial state for an arbitrary chosen value of $M$, all atoms are prepared in the polar ground state $\vec{\psi}_{\rm pgs}=(0,\psi_{\rm pgs}(x),0)^{T}$. This state is then subject to double rotations: ($i$) a spin-1 rotation $e^{if_{y}\pi/2}$, which produces the intermediate state $\frac{\psi_{\rm pgs}}{\sqrt{2}}(1,0,-1)^{T}$, and ($ii$) a rotation $e^{-i\sigma_{y}\theta}$ through angle $\theta$ that is performed on the $m_{F}=\pm 1$ levels around the y-Pauli matrix $\sigma_{y}$. The above procedure leads to $\vec{\psi}_{M}=\frac{\psi_{\rm pgs}}{\sqrt{2}}(\sin\theta+\cos\theta,0,\sin\theta-\cos\theta)^{T}$, and the desired state for a given $M$ is constructed when $2\theta~{}=~{}\arcsin{\left(\frac{M}{N}\right)}$. The state for arbitrary $M$ can be also prepared experimentally by applying the two subsequent electromagnetic pulses [48]. In our calculations, stochastic white noise with variance $\Delta~{}=~{}\frac{1}{2}\frac{\text{particle}}{\text{momentum mode}}$ is added to all Zeeman components of the initial $\vec{\psi}_{M}$ to seed the formation of symmetry-breaking phase domains. The calculations are made for an ensemble of $N_{r}=300$ realizations, the number of grid points $\mathcal{N}=2^{10},2^{11}$ on the box size $L$ large enough to avoid finite size effects, i.e., $L=7\times 10^{2}$, $10^{3}$, $5\times 10^{3}$, $10^{4}$ $\mu$m, which corresponds to $\rho=14.3,10,2,1$ $\mu$m*-1*, respectively. We set the number of atoms to $N=10^{4}$.

As illustrated in Fig. 2 in the example for the $m_{F}=0$ component, the growth of phase domains and their coarsening is observed. The evolution of the density and phase can be divided into three stages: (i) creation of domains seed followed by spin domain formation around $t=10\tau$, (ii) early dynamics characterized by fast reduction of the number of domains, and (iii) further dynamics leading to domains merging at the longest time scale.

The initial exponential growth followed by the phase domains formation is well understood [46] so far relying on unstable modes. Here, in turn, PhOK in the long time limit driven by phase domains merging is described. The description relies on the evolution of the correlation length determined from the first-order correlation function $g^{(1)}_{N}(x,t)=\langle\int{\rm d}x^{\prime}\ \psi_{0}(x^{\prime}+x,t)^{*}\psi_{0}(x^{\prime}+y,t)\rangle$, where $\psi_{0}(x,t)$ is a solution of GPEs for the $m_{F}=0$ Zeeman component. The computed correlation length $l_{h}$ has such a property that $g^{(1)}_{N}(l_{h},t)=h$ (in this paper $h=\frac{1}{2}$). Our initial condition, although trivialized, determines the following. No atoms in the $m_{F}=0$ state implies that its $l_{1/2}=0$ at $t=0$. Since atoms entirely occupied $m_{F}=\pm 1$ states, their correlation length is of the order of the system size. The situation reverses dramatically, when the system is abruptly quenched. Therefore, the analysis of the superfluid order from the correlation length of the $m_{F}=0$ spin component seems suitable for the PhOK investigation in our system. The units chosen for the characteristic length and time in the spin-1 system are $\xi_{s}$ and $\tau=\hbar/(c_{2}\rho)$, respectively.

II Scaling hypothesis and various scaling laws: numerical results

The scaling hypothesis states that during the PhOK process a single length scale $L(t)$ is expected which increases in time as $L(t)\sim t^{1/z_{d}}$, where $1/z_{d}$ is the scaling exponent. The whole variation of the correlation function (1) that occurs during PhOK becomes independent of time, when it is scaled by $L(t)$. We verify the hypothesis in our system and confirm that the correlation function $g_{N}^{(1)}(x,t)$ is one parameter-dependent after the proper re-scaling, i.e., $g_{N}^{(1)}(x,t)=f\left(x/L(t)\right)$. However, the value of the scaling exponent $1/z_{d}$ strongly depends on the system parameters.

Let us refer to the representative examples in Fig. 3, i.e., $g_{N}^{(1)}(x,t)$ (upper panels) and $f\left(x/L(t)\right)$ (bottom panels) with the properly chosen the $L(t)$ scaling. Our numerical calculations show that the scaling exponent is $1/3$ and $3/2$ (or equivalently the dynamical exponent $z_{d}=3$ and $z_{d}=3/2$) for $M=0$ and for a nonzero magnetization value when $q$ is relatively large, respectively, see an example in Fig. 3(c),(d) and Fig. 3(g),(h). The two scaling exponents have already been reported for spin-1 Bose-Einstein condensates with ferromagnetic interactions characterized by negative $c_{2}$ [49, 17]. These results can be explained on the basis of the hydrodynamic (binary mixture) model, as discussed by us in Appendix A. In the case of macroscopic magnetization and a relatively low value of $q$ the best matching to numerical data is observed when using scaling exponents predicted by the vector model with $L(t)\sim(t/\rm{ln}(t))^{1/4}$, see an example in Fig. 3(b),(f). The temporal scaling $L(t)\sim t^{1/2}$ can also be observed, see panels (a) and (e).

In the next two subsections, we discuss in details the particular scaling laws characteristic for the H and B models, and the extent of their occurrences. Our numerical findings are summarized in the last subsection in the form of the phase diagram for the spin-1 antiferromagnetic BECs.

II.1 In the limit of zero magnetization

Let us start the analysis with the zero magnetization case for which the scaling laws associated with the binary mixture or H model in the Hohenberg and Halperin classification of dynamic critical phenomena can be observed.

In the model the fluid flow contributes to the transport of the order parameter [1], in general. This is why hydrodynamical processes like inertial or viscous growth along with diffusion mechanism are included. Each mechanism dominates at a different stage of the domains formation, and eventually one wins at the longest time scale. The H model predicts the following scaling laws: $\sim t^{1/3}$ when the diffusive transport of the order parameter is dominant, $\sim t^{2/3}$ when the inertia of fluids are important, and $\sim t^{1}$ if the viscous process wins [1, 50].

Our numerical GPEs results demonstrate that the diffusive $\sim t^{1/3}$ and inertial hydrodynamic $\sim t^{2/3}$ scaling laws are revealed in the long time limit, while the scaling law  $\sim t^{1}$ from a pure viscous effect is absent. The exact derivation of scaling exponents from GPEs is difficult, but the estimated analysis can be performed — a similar one was already done on the hydrodynamic equations for the binary mixture [1]. We discuss this in Appendix A. In general, the estimation of scaling laws for the H model requires introduction of a non-zero surface tension. In one-dimensional system described by GPEs, it equals zero. The scaling analysis we performed show that the interaction coefficient $c_{0}$ plays the role of surface tension and in a result, he scaling analysis of the hydrodynamic form of GPEs gives the same scaling laws as the typical H model, see Appendix A. Consequently, the appearance of the two scaling exponents $1/3$ and $2/3$ in the systems described by the GPEs can be justified.

The transition from the diffusive to the intertial hydrodynamic scaling laws is clearly visible in Fig. 4 (b). We observe that the average domain wall width is comparable with the width of the phase domain itself at the transition point where the scaling exponent changes. If the size of the phase domain is larger than the size of the average domain wall, then the diffusion transport defines the physics of the system and the scaling exponent. We use that reasoning to estimate the transition point between $z_{d}^{-1}=1/3$ and $z_{d}^{-1}=2/3$ scaling laws. To proceed, let us assume that the fractional size of phase domain over the entire system is given by the fractional volume occupied by the ground state phase $\rho_{0}$, composed of spin domains, i.e., $x_{0}(q)=1-\frac{\sqrt{q_{1}}}{\sqrt{q}}$. This formula can be established from the analysis of equilibrium conditions for the coexistence of 2C and $\rho_{0}$ phases [46]. The width of domain wall between these phases turns out to be set by the $q$-dependent healing length $\xi_{2C}$, estimated by expanding the Bogoliubov dispersion relation of the $m_{F}=0$ Zeeman component in powers of small momentum $k$, $\epsilon_{k}^{(0)}=c_{2}\rho\sqrt{\left(\xi_{s}^{2}k^{2}+1-q\right)^{2}-\left(1-2q\right)}\approx c_{2}\rho\ q\,\left(1+\xi_{2C}^{2}k^{2}\right)+O(k^{4})$ with $\xi_{2C}^{2}=\xi_{s}^{2}\left(1-q\right)/q^{2}$ [46]. The relation $x_{0}(q_{th})\approx\xi_{2C}/\xi_{s}$ or equivalently $1-\frac{\sqrt{q_{1}}}{\sqrt{q_{th}}}\approx\sqrt{1-q_{th}}/q_{th}$, gives the desired condition for the transition point $q_{th}$ between the two different scaling laws. Whenever $q>q_{th}$, the diffusive transport governs domains coarsening. The resulting estimate for $q_{th}$ is shown in Fig. 4(b) by the vertical solid line, and in Fig. 7 by the dashed line.

At this juncture, it is useful to connect the appearance of domains of phase (PhOK) and also domains in the real space, to the occupations of atoms in particular Zeeman states. The system we consider is three component in general. When $M$ decreases at relatively low $q$ the $m_{F}=0$ component becomes macroscopically occupied at longer time scales, while the remaining two components $m_{F}=\pm 1$ are occupied marginally, both to almost the same extent. For the purposes of the dynamics one can expect that the properties of atoms remaining in the $m_{F}=\pm 1$ components become identical. The system starts then to behave like a binary mixture composed of two species of atoms: these in the $m_{F}=0$ state and those in $m_{F}=\pm 1$ treated together as the second species. Similarly, in the large $q$ limit occupation of the $m_{F}=-1$ component becomes marginal while the remaining two $m_{F}=0,1$ are of the only importance. In both regimes of parameters, the binary mixture description of PhOK could become suitable for our system, and therefore, the resulting scaling exponents could appear.

II.2 Macroscopic magnetization

In the case of macroscopic magnetization $M$, we can talk about the regime of parameters for which occupations of all three Zeeman components are significant. Also then the variation of the scaled correlation length $l_{1/2}/\xi_{s}$ versus scaled time $t/\tau$ is expected, as shown in Fig. 5 for various system densities $\rho$. This time, an increase of the correlation length does not exhibit the scaling attached to the binary mixture models. Instead, we observe that $l_{1/2}$ reveals growth typical for vector field models [1]. Those models predict $L(t)\sim t^{1/2}$ (model A) and $L(t)\sim(t/{\rm ln}t)^{1/4}$ (model B) for non conserved and conserved order parameters, respectively.

Better understanding of the evolution of the order parameter in the case of considered system can be provided with definition of $n_{0}=N_{0}/N$ — the fractional number of atoms in the $m_{F}=0$ Zeeman component, where $N_{0}=\int{\rm d}x\,|\psi_{0}(x)|^{2}$.It is shown in Fig. 3(b). Notice that for early times $n_{0}$ is susceptible to large fluctuations. In the range $\sim 10\tau\dots 40\tau$, the order parameter is found to be evidently non conserved. The correlation length scaling exponent is $1/2$ there, see Fig. 3(c). However, for the longest time scale which we are interested in, the change of $n_{0}$ is less significant although its variance is far from being equal to zero. Despite that fact, the value of the scaling exponent changes to $1/4$ with logarithmic correction. The two scaling laws can also be observed by properly re-scaling the correlation function $g_{N}^{(1)}(x,t)$ as demonstrated in Fig. 3. Note, however, that the universal character can have only the scaling revealed in the long time limit, that is $z_{d}^{-1}=1/4$. The temporal $1/2$ scaling is mentioned by us only to demonstrate the multiscaling behaviour of PhOK characteristic for our system in early times. The similar observation and conclusion was also made for ferromagnetic spin-1 Bose-Einstein condensates [15].

The character of the system smoothly transforms to the binary mixture when the value of $q$ increases. The occupation of the $m_{F}=-1$ Zeeman component diminishes to zero and the remaining two, $m_{F}=0,1$ are of the only importance. The change of the scaling exponent to $1/3$ typical for the H model is observed by us independently of $M$.

II.3 Classification of scaling laws

To summarize the results at this stage: we show the variation of scaling laws in time for $q=0.5$, $0.75$, $1.0$, $1.25$ in Fig. 6 depending on the magnetization value. Our observations confirm the presence of mechanisms which are competitive with each other and they lead finally to different behaviour in the long times limit. Therefore, several values of the growth exponents can emerge during the later evolution, namely $z_{d}\to 3/2$, $z_{d}\to 3$ or even to $z_{d}\to 4$ with logarithmic correction. The change of scaling exponent while varying magnetization is typical for our system and appears in a wide range of the value of $q$.

Careful analysis of the numerical results indicates an important role of occupations of particular Zeeman states. We conclude that occupations of atoms in particular Zeeman components give an intuitive picture of which of the scaling law will be revealed in the long dynamics of our system. The vector (B) model fits the best when occupations of all the three Zeeman components are significant. It takes the place for macroscopic magnetization for values of $q$ a bit larger than the critical one, $q_{c}$. In turn, a high concentration of atoms can stimulate stronger interactions, the liquid character and thus the hydrodynamic description would be more accurate. It is when magnetization is small and in the large $q$ limit, $q\gg q_{c}$. In this regime, the hydrodynamic (H) model seems to match the best.

The border between both models can be deduced by matching when the number of atoms in the $m_{F}=-1$ component vanishes, i.e., $N_{-1}(M,q)/N\to 0$. This is illustrated in Fig. 7. The vector field model is realized below the gray solid border line, while the binary mixture model applies above it. The change between the two regions is smooth, and so $z_{d}$ does not perfectly reflect the particular model all around the border line.

III Discussion and conclusion

In this paper, we have explored the superfluid phase-ordering dynamics of an antiferromagnetic spin-1 condensate quenched from the antiferromagnetic state to a state where domains of atoms with different spin projections separate. We have found that the growth of domains is scale-invariant with various dynamic critical exponent $z_{d}$ that are typical for the B and H models. We classified the various scaling exponents due to the system parameters.

When the occupation of the $m_{F}=-1$ Zeeman state is significant, We conclude that the scaling exponents typical for vectors fields model appear ($z_{d}=4$ with logarithmic correction). At short times we observe also that temporally the correlation length $l_{1/2}$ exhibits the scaling $L(t)\sim t^{1/2}$. These all properties are characteristic for multiscaling of the PhOK.

On the other hand, when the occupation of the $m_{F}=-1$ component is marginal the scaling laws typical for the H model are realized, namely $z_{d}=3/2$ when $q$ and magnetization are small, and $z_{d}=3$ in the large $q$ limit for any $M$. While the appearance of scaling laws typical for the H model can be justified from the hydrodynamic form of GPEs, the occurrence of scaling typical for B model is not evident. The latter provides an interesting direction for future work.

To conclude: the antiferromagnetic spin-1 condensate captures the universal two-model feature of the PhOK in the long time limit. The system parameters ($q$, $M/N$) set the physics, and determine the dynamical scaling exponent corresponding to the model H or B. Switching between the models by changing the initial state or the $q$ parameter is allowed.

Acknowledgements

We wish to thank Piotr Deuar and Michał Matuszewski for useful discussions, and Oleh Hul for careful reading of manuscript. EW and JP acknowledge support from the National Science Centre (Poland) grants Nos. 2015/18/E/ST2/00760 and 2012/07/E/ST2/01389, respectively.

Author contributions statement

J.P. and E.W. equally contributed to the numerical calculations, analysis of the results and writing the manuscript. A.S. contributed to the numerical simulations of the $M/N=0.5$ case.

Appendix A Scaling analysis from spin-1 hydrodynamic equations

In the following, we present the hydrodynamic forms of the GPEs. The wave functions in Eq. (3) describing atoms in particular components are replaced by $\psi_{m_{F}}=\sqrt{n_{m_{F}}}\,e^{i\,\theta_{m_{F}}}$ – according to Madelung transformation. After some algebra, we obtain corresponding equations for the densities $n_{m_{F}}$, namely

[TABLE]

and the corresponding Navier-Stokes equations for velocities $v_{m_{F}}=\frac{\hbar}{m}\nabla\theta_{m_{F}}$

[TABLE]

where $\theta=\theta_{1}+\theta_{-1}-2\theta_{0}$, and $U_{1}=c_{0}n+qc_{2}\rho+c_{2}(n_{1}-n_{-1}+n_{0})$, $U_{0}=c_{0}n+c_{2}(n_{1}+n_{-1})$, $U_{-1}=c_{0}n+q\,c_{2}\rho+c_{2}(n_{-1}-n_{1}+n_{0})$.

We will perform here the scaling analysis corresponding to the one presented in [1] for the hydrodynamic (B) model. The chemical potential can be identified as $\tilde{\mu}=U_{0}$, see details in  [46], and therefore the derivative of chemical potential as $\phi\nabla\tilde{\mu}\to\nabla U_{0}$. Taking this as an assumption one can perform the scaling analysis for $m_{F}=0$ as in below.

\small1⃝ Scaling by diffusion mechanism:

In the diffusive regime it is assumed that the rate of change of the domain size $\frac{dL}{dt}$ is associated with the chemical potential gradient $\nabla U_{0}$. The scaling of $U_{0}$ is linked to the scaling of the local density $U_{0}\sim c_{0}n\sim c_{0}N/L$, therefore one has:

[TABLE]

Note, the role of the surface tension $\sigma$ typically introduced in analysis of binary mixture [1] is played here by the interaction coefficient $c_{0}$.

\small2⃝ Scaling by inertial growth:

If the inertial terms become important, the scaling law is given by the relation $(v_{0}\nabla)v_{0}\sim\nabla U_{0}$. Then, one recovers the same scaling as for the H model, namely

[TABLE]