Miscibility in coupled dipolar and non-dipolar Bose-Einstein condensates

We perform a full three-dimensional study on miscible-immiscible conditions for coupled dipolar and non-dipolar Bose-Einstein condensates (BEC), confined within anisotropic traps. Without loosing general miscibility aspects that can occur for two-component mixtures, our main focus was on the atomic erbium-dysprosium ($^{168}$Er-$^{164}$Dy) and dysprosium-dysprosium ($^{164}$Dy-$^{162}$Dy) mixtures. Our analysis for pure-dipolar BEC was limited to coupled systems confined in pancake-type traps, after considering a study on the stability regime of such systems. In case of non-dipolar systems with repulsive contact intneeractions we are able to extend the miscibility analysis to coupled systems with cigar-type symmetries. For a coupled condensate with repulsive inter- and intra-species two-body interactions, confined by an external harmonic trap, the transition from a miscible to an immiscible phase is verified to be much softer than in the case the system is confined by a symmetric hard-wall potential. Our results, presented by density plots, are pointing out the main role of the trap symmetry and inter-species interaction for the miscibility. A relevant parameter to measure the overlap between the two densities was defined and found appropriate to quantify the miscibility of a coupled system.


Introduction
Following the first experimental realization of Bose-Einstein condensation (BEC) in 1995 [1], substantial progress has been verified in experimental and theoretical investigations with ultracold quantum gases, which can be traced by several paper reviews on the subject (as references in [2]), as well as by a recent book written by experts in the field [3]. One can also follow the advances in quantum simulations and control of BECs in reviews such as [4,5]. From studies with ultra-cold atoms, we can improve our understanding on quantum properties of a large variety of bosonic and fermionic systems, as well as molecular configurations with different atomic species. By controlling atomic properties potential technological applications exist from ultra-precise clocks till quantum computation. In this respect, by considering BEC mixing with different atomic and molecular configurations, also with degenerate complex atoms (alkaline-earth, lanthanides) and Fermi gases, one can investigate the crossover of BEC properties with Bardeen-Cooper-Schrieffer superconductivity, superfluid to Mott insulator transition in bosonic and fermionic systems, quantum phases of matter in optical lattices, ground-state fermionic molecules [6].
BEC with two-components was first produced with different hyperfine states of 87 Rb [7]. This is a simple example of a multicomponent system made with a mixture of two-species of bosons. Following that, one should notice several investigations using binary mixtures. This can be exemplified with works considering the dynamics of binary mixtures with bosons and fermions [8]; studies on the dynamics of phase separation and on how to control it [9,10]; also considering mixtures with different isotopes of the same atomic species [11], or with different atomic species [12]. As a relevant characteristics of multicomponent ultracold gases, we have their miscibility behavior which will depend on the nature of the interatomic interactions between different species. Miscible or immiscible two-component BEC systems can be distinguished by the spatial overlap or separation of the respective wave-functions of each component. Their phase separations were observed in spinor BECs of sodium in all hyperfine states of F=1 [13]. The advances in the experimental investigations with multi-component BECs have activated a large amount of theoretical descriptions applied to condensed mixtures having spatially segregated phases, by studying their properties related to static and dynamical stability [14][15][16][17][18][19][20][21][22][23][24][25][26].
Dipolar atomic and molecular systems, as well as mixtures with different dipolar atoms, have been explored in theoretical works connected with BEC since 2000 [27,28], followed by several other investigations motivated by the increasing experimental possibilities in cold-atom laboratories to probe the theoretical analysis and eventual proposals. The theoretical effort in this direction can be traced back by the following sample works, [29][30][31][32][33][34][35][36][37][38][39][40][41]. A more complete bibliography on the subject, covering experimental and theoretical approaches, can be found in recent reviews and dissertations, as given in [41][42][43][44][45].
When considering the investigations with dipolar systems in cold-atom laboratories, the pioneer work is the experiment with chromium 52 Cr, reported in [46]. Among the works following that, we have the investigations on the stability of dipolar gases [34], as well as on collisions of dipolar molecules [47]. With ultra-cold atoms having non-negligible magnetic moments, we have experiments with dysprosium described in [48], considering the dysprosium isotopes 162,164 Dy; and the investigation with erbium 168 Er in [42]. More recently, quantum droplets have being observed in a strongly dipolar condensed gas of 164 Dy [49], with new features being verified for dipolar BECs, due to the competition between isotropic short-range contact interaction and anisotropic long-range dipole-dipole interaction (DDI).
By considering the quite interesting recent investigations in ultracold laboratories with two-component dipolar BECs, studies on stability and miscibility properties are of interest due to the number of control parameters that can be explored in new experimental setups. The parameters are given by the strength of dipoles, number of atoms in each component, inter-and intra-species scattering lengths, as well as confining trap geometries or optical lattices. Among the theoretical studies cited above (most concentrated on stability of dipolar condensates), we have some of them are particularly related to miscibility of coupled BECs and structure formation, as [19,31,37,40].
In the present paper, our main proposal is to discuss miscibility conditions for general three-dimensional (3D) atomic BEC systems, which are constituted by two-coupled dipolar or non-dipolar species confined by asymmetric cylindrical harmonic traps. Due to stability requirements, the dipolar systems that we are considering will be confined in pancake-type traps. In case of non-dipolar coupled systems, we also discuss the miscibility by considering cigar-type symmetries. For our study on the miscibility, we start with a brief discussion by considering the homogeneous case. In order to simplify the formalism and a possible experimental realization, both species are assumed to be confined by a cylindrical trap with the same aspect ratio. Without losing the general conclusions related to miscibility of two-atomic dipolar BEC species, most of our study will focus on the particular 168 Er-164 Dy and 164 Dy-162 Dy mixtures, motivated by the actual experimental possibilities [42,48].
Our numerical results for the coupled dipolar Gross-Pitaevskii (GP) equation are presented by using different parameter configurations for the trapping properties, as well as for the inter-and intra-species twobody contact and dipolar interactions. The parameter region of stability for the dipolar system is discussed for different trap-aspect ratio and number of atoms in each species. As we are going to evidence, for a given mixture of two condensates confined by harmonic traps, the main parameters of the system that are possible to be manageable in an experimental realization with focus on the miscibility are the trap-aspect ratio and the twobody scattering lengths (these ones, controlled via Feshbach resonance techniques [5,50]).
Within our full-3D model for the coupled densities, when considering pure dipolar systems trapped in pancake-shaped harmonic potentials, we are also discussing some unusual local minimum structures and fluctuations in the densities, which are verified when the system is near the instability border (considering the critical aspect ratio and atom numbers). These structures are verified for coupled systems that are partially immiscible ( 168 Er-164 Dy, in the present case), as well as when it is completely miscible, such as 164 Dy-162 Dy. Such structures, verified for well defined trap-aspect ratio and number of atoms in stable configurations, suggest possible experimental studies with two-component dipolar BECs, considering miscible and immiscible systems.
The next sections are organized as follows. In section 2, we present the general 3D mean-field formalism (in full-dimension and dimensionless) for trapped two-component dipolar BECs, together with the definition of relevant parameters, as well as the numerical approach we are considering. In section 3, we first write down the miscibility conditions for homogeneous coupled systems, followed by the definition of an appropriate miscibility parameter, which is found appropriate to measure the overlap between densities of a general coupled system. Our numerical results are organized in two sections, in order to characterize the main relevant conditions for the observation of miscibility in coupled BEC systems. The role of the trap symmetry for the miscibility, considering different dipolar and non-dipolar mixed systems, is analyzed in section 4. In view of the particular relevance of the inter-species two-body interactions on the miscibility, the corresponding results are presented and discussed in section 5. Finally, in section 6, we present a summar with our principal conclusions and perspectives.

Formalism for coupled BEC with dipolar interactions
For a mixed system with two atomic species identified by i=1, 2, having their masses, number of particles and local time-dependent wave-functions given by m i , N i , and y y º ( ) t r, i i , respectively, the general form of the mean-field GP equation, for the trapped system with dipolar interactions, can be described by [29], defining the magnetic-type dipolar interaction between particles i and j. The nonlinear contact interactions between the particles are given by , where a 11 , a 22 and = a a 12 21 are the two-body scattering lengths for intra (a ii ) and inter (a 12 ) species, with m ij being the reduced mass In the above, both wave-function components are normalized as For the confining trap potentials we assume harmonic cylindrical shapes, with frequencies w i and aspect ratios l i , such that where it will be assumed that each species i is confined by an angular frequency w i along the x-y plane, r º + ˆr e r e ; 1 1 2 2 and with l w i i along the z-directionr e 3 3 . The trap will be spherically symmetric for l = 1; i will have a cigar shape for l < 1; i and a pancake shape when l > 1 i . For the magnetic-type dipolar interaction between particles i and j, with respective dipole momentum strength given by m m m º D ij i j 0 (m 0 being the permeability in free space and m i the dipole moment of the species i), we have where -¢ r r determines the relative position of dipoles and θ is the angle between -¢ r r and the direction of polarization.
Let us rewrite (1) in dimensionless quantities, with the first component defining the scales for length, with  w º ( ) l m 1 1 and energy, w 1 . Within these units, we introduce new dimensionless variables and redefine the parameters such that  With the above, and also by redefining the wave-functions for the atomic species i, , the expression (1) can be rewritten as the following dimensionless coupled expressions:

Numerical approach
For the numerical approach, used to obtain our results when solving the full-3D coupled equations (6) and (7), we have employed the split-step Crank-Nicolson method, which is detailed in similar nonlinear studies, as in [17,51], where one can find more extended analysis and details on computer techniques convenient for nonlinear coupled equations, facing stability and accuracy of the results. In view of the particular integrodifferential structure of the coupled nonlinear differential equations when having dipolar interactions, we had to combine our approach in solving coupled differential equations with a standard method for evaluating dipolar integrals in momentum space [29,52,53].
By looking for stable solutions, the 3D numerical simulations were carried out in imaginary time with a grid size having 128 points for each dimension, where we have D = D = D = x y z 0.2 for the space-steps and D = t 0.004 for the time-step. The results were quite stable, verified by taking half of the mentioned grid sizes. As a preliminary calculation, which also help us to check the numerical code, we reproduce the stability diagram obtained in [29] for a spherically symmetric trap (l = 1), where only one atomic species was used, 52 Cr, for the coupled system, with the two species ( = m m 1 2 ) having opposite polarizations along the z direction. Besides that, in our analysis, we have also verified the stability of the numerical results by studying the effect of varying λ. In this regard, once verified that pancake-type configurations are required for stable dipolar configurations, the miscibility in cigar-type traps is being analyzed only for non-dipolar systems.
3. Miscibility of coupled systems 3.1. Homogeneous case with hard-wall barriers In order to characterize the transition between miscible and immiscible states, let us consider a simpler case, for the homogeneous 3D system with and hard-wall barriers, following the simplified energetic approach presented in [15]. Within this approximation, the miscible-immiscible transition (MIT) can be characterized by a threshold parameter, which is defined by the relation between the two-body repulsive interactions. The criterium for miscibility, also quoted in some recent works on binary BEC mixtures (see, for instance, [19,21]), was previously obtained from stability analysis of the excitation spectrum in [14]. It can be easily generalized to include all the two-body repulsive interactions.
As in the present case we are interested in mixed configuration with contact and dipolar interactions, we should first obtain a simple relation for the dipolar interactions appearing in the formalism. For that, we follow the approach given in [29] to deal with the integro-differential formalism where we have a divergence of the integrand in the zero limit for the inter-particle distances. The convolution theorem and Fourier transforms are applied for the magnetic dipolar potential and for the density components y | | i 2 . The Fourier transform of the dipolar potential (4) is given in terms of a cut-off parameter, which is of the order of the atomic radius. As this parameter is much smaller than a significant length scale of the system, one can safely consider the limit where it is zero, such that the dipolar potential between the atomic species i and j is given by where q k is the angle between the wave-vector k and the dipole moment.
Therefore, by including together the contact and dipolar interactions, the condition for MIT can be written as where D = 0 defines the critical value for the transition from miscible (D > 0) to immiscible (D < 0) systems. The mixed coupled state will have a lower total energy when the mutual repulsion between atoms is large enough such that we have D < 0, which is characterizing the system in an immiscible phase. The system will be in a miscible phase when D > 0, with a critical border for the transition given by D = 0. The MIT occurs when the inter-species and intra-species interactions are balanced. As verified, the two-body interactions are just being scaled by the dipolar interactions, in this simple case that the system is confined by hard-wall barriers. The observation that the dipolar interactions are not playing a significant role results from the fact that the DDIs averaged out for homogeneous gases. However, for numerical comparisons with more general miscibility conditions, it is still quite useful, as it does not depend on the number of atoms and condensate size.

Miscibility in general coupled systems
The miscible to immiscible transition for a coupled system, as defined by the critical limit Δ in (9), is changed from a first-order transition that occurs for homogeneous system to a second-order one when the kinetic energy is taken into account, as discussed in [20]. Therefore, in order to evaluate the miscibility of a coupled system for a more general case, we define a parameter which gives an approximate measure of how much overlapping we have between the densities of both components. By the definition of a parameter related to the densities we are modifying the one suggested in [20], where the overlap between the wave-functions was considered. Both definitions are equivalent in particular cases. Our parameter to define the miscibility of a coupled system is given by As f 1 and f 2 are both normalized to one, also having the same center, this expression implies that h = 1 for the complete overlap between the two densities, decreasing as the overlap diminishes. Therefore, we can define the system as almost completely immiscible when h  1 (close to zero); and, almost completely miscible when η is close to one. This parameter is extending to general non-homogeneous mixtures the MIT criterium (9) discussed for homogeneous systems. With η, intermediate cases can be determined when the system is partially immiscible or partially miscible. From our observation, which will follow from the analysis of results obtained for the densities in the next sections, when  h 0.5 the system shows already a clear space separation, with the components having their maxima in well separated points in the space, such that we can already define the coupled system as immiscible. For η between 0.5 and 0.8, the two densities start having an increasing overlap with their maxima approaching each other, such that we can define the system within this interval as partially miscible. The maxima of the densities are close together for  h 0.8, when we assume the system is miscible.
As a side remark, we should observe that we are treating two condensates with different atomic species symmetrically distributed around the center. As a general expected behavior, the density of the more massive species should be closer to the center, with the other density being pushed out. This can be explained even before activating particular interactions between the atoms, as the kinetic energy of the massive species is smaller than the corresponding kinetic energy of the less massive species.
Next two sections are dedicated to present our main results on the miscibility of coupled dipolar and nondipolar BECs. When choosing non-zero dipolar parameters, the results are exemplified by two mixtures for a better characterization of the miscibility properties. Within the actual experimental possibilities in BEC laboratories, we consider the erbium-dysprosium ( 168 Er-164 Dy) and dysprosium-dysprosium ( 164 Dy-162 Dy) mixtures.

Miscibility results-role of the trap symmetry
The results analyzed in this section for the miscibility properties of coupled atomic condensates are more concerned with the role of the symmetry of the harmonic trap, considering different configurations for the possible internal interactions (contact and dipolar) between the atoms of both species. The contact interactions are characterized by the atomic two-body scattering lengths, with dipolar interactions due to the magnetic dipole moment of each atomic species. In the following, we split the section in three subsections for clarity. In part A, we consider miscibility in the case we have no contact interactions, starting with a discussion on the stability of such dipolar systems in terms of the trap-aspect ratio and number of atoms. As shown, for realistic number of atoms, we need pancake-shaped traps. Next, in part B, our analysis is concentrated in non-dipolar coupled structures, where we can examine cigar-type and pancake-type BEC configurations. We conclude the section with part C, by analyzing our results for the case that we have both dipolar and contact interactions.

. Stability analysis
The homogeneous case of purely dipolar condensate is unstable due to anisotropy of the DDI. However, in a way similar as the case of homogeneous non-dipolar BEC that are unstable (for attractive two-body interactions, < a 0), the instability usually can be overcome by applying some external trap, which will help to stabilize the dipolar BEC by imprinting anisotropy to the density distribution [34]. This is applicable for single and multicomponent dipolar BECs. On the stabiliy of two condensates purely dipolar, separated by a distance, we can mention [38]. For the case of single component dipolar condensates, it was previously shown in [27,28] that, by increasing the aspect ratio λ one can obtain a more stable configuration due to the DDI becoming effectively more repulsive.
As shown by the stability diagrams given in the two panels of figure 1, pure dipolar condensates require pancake-type traps to be stable, with the coupled mixture becoming less stable when the dipolar strengths (inter and intra-species) are close to the same values (for some fixed ratio of number of atoms in both species). This effect is mainly due to the inter-species repulsion in comparison with the corresponding intra-species ones.
To become clear this effect, our results are given for the coupled equations (6) and (7)  . Also, there is a critical trap aspect ratio (l c ) for the stability, As one can verify from the upper panel, for the 168 Er-164 Dy mixture with N 1 = 3000, 8000, and 10 000 the, this critical aspect ratio starts from l » c 7, 11, and 13 respectively. On the other case, for the 164 Dy-162 Dy mixture with the same sets of N 1 (= 3000, 8000, and 10 000), the critical lower limit for stability starts with l » This stability analysis will be considered in the following study on the miscibility of two coupled mixtures. However, it can also be of interest for corresponding experimental investigations (for the specific mixtures we have considered or for other similar dipolar mixtures).

Structure of coupled pure-dipolar condensates
Once analyzed the stability for pure dipolar coupled systems, in this subsection we characterize the role of the trapping aspect ratio λ in the structure and miscibility of a pure-dipolar coupled condensate. Therefore, as in the preceding subsection, all the two-body scattering lengths ( = a 0 ij ) are fixed to zero and we consider both coupled systems with 168 Er-164 Dy and 164 Dy-162 Dy. Our numerical results for the structure of the coupled system can be visualized through density plots, with the immiscible and miscible regimes characterized by the parameter η defined by (10). In figure 2 we present 3D surface plots for the densities f | ( )| x y , , 0 2 in the (frames (a), (c), (e), (g)) and f | ( )| x z , 0, 2 (frames (b), (d), (f), (h)), given respectively in the (x, y) and (x, z) surfaces. These plots are providing a 3D visualization of the density overlapping and distribution of the two-components. More close to the center, we have the more massive component (in red) of the mixture, which is enhanced when the system is more immiscible. How close to the center is the other species (in green) will depend on how miscible is the mixture, being clearly identified for immiscible mixtures. This can be verified by comparing the set of panels (a)-(d) for 168 Er-164 Dy with the ones (e)-(h) for 164 Dy-162 Dy.
Corresponding to the panels of figure 2, we also have the one-dimensional (1D) plots for the densities in figure 3, which are given as functions of one of the dimensions, x or z, with the other two dimensions at the center. From these 1D densities, the amount of overlapping between the densities can better be observed, being more helpful when comparing different parameter configurations. Both systems, 168 Er-164 Dy and 164 Dy-162 Dy, represented in figures 2 and 3, are been shown for pancake-shaped traps with l = 7. The number of atoms that was considered was dictated by the stability of the mixtures (as one can follow from figure 1), such that we can study the density properties when the mixtures are stable but not far from the border where they can become unstable.
In  condensate is confined in a pancake-type cylindrical trap, with aspect ratio l = 7, such that the distribution along thez axis is more concentrated near the center (z = 0) than the other two directions. As also characterized by the miscibility parameter, h = 0.77, in the present case the coupled system is partially miscible. This is quite well represented in the left four panels given in figure 3. , in the frames (g), (h) of these figures. In this case, where we have two isotopes of the same atom, we notice that the system is completely miscible, having the same value close to one for the miscibility parameter h = 0.99. The complete overlap between the densities of the two components are clearly shown in panels (e)-(h) of figure 3.
In these figures 2 and 3, by using the aspect ratio l = 7, for each one of the mixtures we can identify a nontrivial structure emerging in the condensate near the boundary of stability, with a local minimum at the center (for both components) in the symmetrical x-y plane (see left panels of both figures). At the z-direction, we have a normal Gaussian shape, as seen in the right panels of the figures.
In both the mixtures shown in figures 2 and 3, the upper panels (a), (b) and (e), (f) are for parameters very close to the stability threshold, as indicated by the results given in figure 1. Therefore, the structure observed for the density of the component 1, which has the largest fraction of atoms in the coupled system, can be explained by fluctuations close to the instability regime. For each system, by going from smaller values, shown in panels (c), (d) and (g), (h), to larger values, shown in panels (a), (b) and (e), (f), of N 1 (keeping N 2 about the same), we are approaching the unstable regime, visualized by the occurrence of oscillation peaks in the density of the dominant component in the system. The observed number of four peaks around the center is related to stability requirements when the number of atoms is close to the maximum allowed limit for a given trap asymmetry. By small variations of the parameters near the stability, before the collapse of the system, it is possible to increase such number of peaks.
In both sets of systems shown in figures 2 and 3, while observing the peak oscillations occurring for the first component, one can observe only the biconcave-shaped condensate for the second component. In this regard, we should also note that biconcave-shaped structures with local minimum in dipolar BECs have already been reported in the cases of single component condensates. In [33], such structures are explained as due to roton instability for certain specific pancake-type trap aspect ratios λ (≈ 7, 11, 15, 19, ...,), being not observed for other values of λ. In our analysis, we confirm the values reported in [33]. However, when considering singlecomponent BECs, we have also verified these type of structures for other particular trap-aspect ratios, such as λ ≈ 8, 12, 16, and 20. However, for the cases of fully anisotropic traps, where two aspect ratios are considered, this kind of density fluctuation has also been observed in [39], for single component dipolar BECs. When considering two-component systems, we have verified this kind of structured states with no particular restriction on the values of λ. The biconcave structure occurs for particular aspect ratios of the trap, due to the repulsive interaction between the dipoles. In the present study for a coupled dipolar system, we observe that we have also the repulsive inter-species DDI in addition to the intra-species DDI, leading to the observed four peaks in the density oscillations. The biconcave structures with a local minimum and few peaks are occurring near unstable regimes, which makes the experimental reproduction non-trivial to be observed.
In order to illustrate the behavior of the coupled densities, as we go from the parameter region close to the border of stability to a more stable configuration, we show three illustrative density plots in figure 4 considering the particular case with l = 8 for the mixture 168 Er-164 Dy, with fixed number of dysprosium atoms = N 950 2 . From the panels (a) to (c) we are varying the number of erbium atoms from 3000 to 2000. We notice that, by going from a more stable region (shown in panel (c), where = N 1 2000) to a region near the instability (panel (a), where = N 1 3000) the shape of the density distribution changes from a Gaussian to a biconcave format. By approaching very closely the instability region (increasing N 1 ), one can observe the fluctuation in the density around the center.
Besides the small mass difference between the atoms, in these cases of pure-dipolar coupled systems, the main difference is given by the respective dipolar interactions, which is directly related to the magnetic moment dipoles of the components. While for the 168 Er-164 Dy we have quite different moment dipoles for the two atomic components (one is about half of the other), for the mixture 164 Dy-162 Dy we have about the same dipolar parameters, with the inter-and intra-species interactions balancing each other. In view of that, the first mixture (with erbium) is expected to be more immiscible. This is reflected in the corresponding values of η (= 0.77 for 168 Er-164 Dy, and = 0.99 for 164 Dy-162 Dy). The results for η are consistent with the approximate criterium that one could use for homogeneous mixtures, where we have D » -0.014 for 168 Er-164 Dy, whereas D » 0 for 164 Dy-162 Dy. However, with the parameter η we can have a more realistic indication of the partial overlap between the densities.

Miscibility in non-dipolar coupled systems 4.2.1. Role of the trapping aspect ratio
The effect of the geometry due to the external harmonic trap potential on the phase separation of the mixtures is investigated in the next, by considering different aspect ratios λ for the cases when the nonlinearity is given at least by repulsive two-body interactions. The miscibility parameter, given by the factor η, indicating the amount of mixing in the densities of the two-component, is presented as a function of λ in figure 5, for the case that twobody contact interactions are fixed such that = = a a a 40 11 22 0 and = a a 50 12 0 , with fixed number of atoms for both components. As discussed when this parameter was defined, a complete overlap between the densities implies in h = 1, being zero in the other limit of a complete immiscible system. For the given non-zero contact interactions, we observe that a complete miscible coupled dipolar system cannot be obained within a stable configuration, as the maximum value verified for the miscibility parameter is below 0.6 for both dipolar systems and number of atoms that we have considered. However, as we have already pointed out in the previous subsection, a value of η near 1 is possible to be obtained in case of complete dipolar systems (when the contact interactions are set to zero).
The purpose of this subsection is mainly to discuss full non-dipolar systems, with results given in figure 5 for non-zero dipolar interactions (for the 168 Er-164 Dy, as well as for 164 Dy-162 Dy). These two kind of mixtures are of interest to show how the dipolar interactions can modify the miscibility behavior of a coupled BEC system. As noticed by the red-solid line with empty circles given in this figure (with = d 0 ij ), for non-dipolar systems we can also obtain stable almost immiscible states in cigar-type configurations with l < 1. With non-zero dipolar parameters, and considering the given contact interactions, stable condensates are limited to pancake-type configurations with l > 6 for 168 Er-164 Dy mixture; and l > 15 for the 164 Dy-162 Dy mixture. As we can see from these results for a non-dipolar coupled system, the immiscibility is more evident when the trap is strongly deformed. It can happen for pancake-type traps, as well as for cigar-type traps. However, due to the stability of the condensates, it should be more favorable to build immiscible coupled condensates with pancake-type traps.

Miscibility in cigar-and pancake-type traps
For non-dipolar mixtures, we resume our analysis by the 1D density plots given in figure 6, followed by a corresponding 3D illustration in figure 7. For the parameters we have both components with the same number of atoms ( = = N N 6000 1 2 ), and with repulsive two-body interactions such that = = a a a 40 11 22 0 and = a a 50 12 0 . In figure 6, the panels (a)-(c) display the central densities ( = = y z 0) as functions of x; with the panels (d)-(f) showing the densities (for = = x y 0) as functions of z. From these density plots, we can also infer how well the miscibility parameter η, given by (10), can be quantitatively used to estimate the miscibility condition of a two-component condensate. In the upper panels,   extremes for the two components, in the case of a pancake-type trap with l = 20, shown in panels (c) and (f). On the other hand, for l = 1, shown in the panels (b) and (e), we can observe a strong superposition of the two components, with their maxima at the same position, with the miscibility parameter close to one (h = 0.99).
The two 3D illustrations shown in figure 7, for the strong deformed cases with l = 0.1 (cigar type) and l = 20 (pancake type) are also characterizing the corresponding density distributions.
By considering the above discussion on the usefulness of the miscibility parameter η, in the next density results that will be shown we rely on the value of this parameter as a relevant quantitative observable for the miscibility analysis and corresponding surface visualization of the densities. As observed, for that we can take h < 0.5 as being almost immiscible, with maxima for the two components in well distinct positions. In the other Figure 9. The miscibility η is shown as a function of a 12 in three frames, with fixed intra-species contact interactions, = = a a 11 22 40a 0 (a), 10a 0 (b) and 0 (c). The aspect ratio is l = 20, with number of atoms = = N N 6000 1 2 . The dipolar parameters used for the curves shown in each panel are given inside the frame (c). When the system turns out to be unstable, the lines are interrupted or no-results shown.

Miscibility in coupled BEC with dipolar and contact interactions
The miscibility behavior of the coupled condensate in terms of the aspect ratio λ is further investigated in figure 8 for dipolar systems in addition to the repulsive two-body contact interactions already used in figures 6, 7 ( = = a a a 40 11 22 0 , = a a 50 12 0 ). In this case, we consider the coupled system 168 Er-164 Dy, with the dipolar parameters given by = ). By considering the same atom number for the components = = N N 6000 1 2 , as verified before, stable results are possible only for pancake-type trapping potentials with l > 6. In general, pancake traps provide strong axial confinement and help to increase the repulsive interaction along radial directions which guides to phase separation between the mixtures. In figure 8, by increasing the values for the aspect ratio (such as l = 6, 10 and 20), we can verify how the structure of the densities for the two components varies in an immiscible condition. The above observations for dipolar and non-dipolar BECs explain their specific characteristics. In particular, for erbium-dysprosium coupled dipolar condensate in stable configurations (l > 6), we are clarifying that it presents immiscible density structures, which are more pronounced for larger values of the trap-aspect ratio.

Miscibility results-role of the inter-species interaction
Besides the trap symmetry, another quite relevant parameter for the miscibility of coupled condensates is given by the inter-species interaction. In this case, when considering a system with fixed dipolarity, the appropriate parameter, which can be tuned via Feshbach resonance techniques [5], is given by the inter-species contact interaction. Therefore, our aim in the following analysis is to characterize the role of the inter-species scattering length a 12 for the miscibility. For that, we assume both species have the same intra-species scattering lengths, = = a a a 40 11 22 0 , choosing a particular large pancake-type trap with l = 20 and fixed number of atoms for . These parameters are dictated by the previous stability analysis, looking for stronger characterization of miscibility properties.
The general behavior for the miscibility can be verified by the parameter η, which is shown in figure 9 as a function of the ratio a a ). We should point out that the miscibility is going down slowly (increasing the immiscibility of the mixture) to an approximate constant value for  ¥ a 12 . Even in the non-dipolar case: with = a ii 40, 10 and 0, we have h~0.14, 0.10 and 0.04, respectively. This can be understood from the residual mixing of the two-component wave-functions near the unitary limit of the inter-species (when  ¥ a 12 ). In this limit for the inter-species, we should notice that we have the well-known Efimov effect [54], with increasing number of three-body bound and resonant states mixing the two component system. Already observed in cold-atom laboratories [55], this a pure quantum effect that occurs near zero two-body binding, when the effective potential goes as the inverse-square of the distance, and we have a long extension of the corresponding two-body wave functions. The role of this effect on coupled systems deserves further analysis in experiments, in particular when varying the mass-ratio of the mixture, as the resonant states have well-known theoretical predictions [56] also for mass-asymmetric mixtures.
Next, we can examine the role of the magnetic dipolar interactions in the miscibility of the two components. First, we should remind that these interactions between the magnetic dipoles are long-range ones, which go as the third power of the distance between the dipoles. They should act more effectively when the components are close together, but have residual effects due to long range behavior, which makes the miscibility be reduced to a non-zero constant value for  ¥ a 12 , in a way similar to the non-dipolar case in this limit, such that we have , and the intra-species scattering lengths are fixed such that = = a a a 40 11 22 0 . The components are being identified inside the upper panels. combined effects of two long-range interactions. In the other extreme, we can see the role of dipolar interactions in the miscibility by looking the region where = a 0 12 in the panel (c). When we have no inter-species dipolar interactions (with only one of the species having magnetic moment) there is no coupling between the systems, such that the overlap between the densities is partial. One of the densities is given by a linear trapped equation; with the other given by a nonlinear equation (confined by the same harmonic trap, but additionally having the repulsive dipolar interaction). As larger is the effect due to repulsive dipolar interactions, less miscible is the system. When both intra-species dipolar interaction are switched on, we have also the corresponding interspecies parameter coupling the system, with the net repulsive effect being averaged out.
Another aspect verified in the panels of figure 9 refers to the stability of the system, which can be observed for this specific case that we have a pancake-shaped trap with l = 20 and fixed number of atoms N i =6000. In general, from figure 9 one can observe that the immiscibility will significantly increase by increasing a 12 in all the verified situations. Also, by taking a fixed value for this inter-species contact interaction with > a a ii 12 , from the three panels one can verify that by decreasing the intra-species contact interactions = a a 11 22 , the immiscibility of the system increases (η decreases).
Density profiles of the coupled system are represented in the next figures 10-13, where we show results considering three values of the two-body inter-species scattering lengths, with the corresponding intra-species contact interactions a ii kept fixed at 40a 0 . In all these cases, in order to facilitate the comparison of the results, we assume the same three values for a 12 (= a 10 0 , a 43 0 and a 60 0 ) and maintain the trap-aspect ratio fixed to a pancake-type, l = 20, with the number of particles given by = = N N 6000 1 2 , as in the plotted results of figure 9. For a better quantitative comparison of the results, we present the densities in 1D plots, as functions of x , and the intra-species scattering lengths are fixed such that = = a a a 40 11 22 0 . The components are being identified inside the upper panels.
(left panels) and z (right panels). In figure 10, we have the case of non-dipolar systems; in figure 11 when considering just one of the components being dipolar. In figures 12 and 13, we consider two cases, when both components are dipolar, with parameters corresponding to the mixtures 168 Er-164 Dy and 164 Dy-162 Dy, respectively. First, by considering non-dipolar BECs, in figure 10, with three values for the inter-species scattering length, (a 12 varying between 10 and 60a 0 , with = = a a a 40 11 22 0 ). The transition is expected to occur close to » a a 12 11 ( » m m 1 2 ). This is consistent with the results shown in this figure, where we can verify that for > a a 12 11 there is already a clear separation between the two species, with h = 0.28 for = a a 60 12 0 . It is also shown that for = a a 10 12 0 the mixture is completely miscible (h = 1), having a sudden transition for = a a 43 12 0 (just above = a a 40 11 0 ) to a partial immiscible state. By further increasing a 12 , an almost complete immiscible state is reached (for = a a 60 12 0 , h = 0.28) . Next, in figure 11, we study the case where only the component 1 is dipolar. By comparing with the nondipolar case, we can observe the effect of the dipolar interaction in breaking the sudden transition verified when increasing the inter-species contact interaction, such that the transition is much softer than the one shown for the non-dipolar case.
In order to observe the effect in the miscibility when changing the inter-species scattering length for the cases that both components are dipolar, we present 1D plots for the densities in figures 12 and 13. For the kind of mixture considered in 12 , the 168 Er-164 Dy, we notice that the system is partially miscible even at < a a 12 11 , becoming almost immiscible for > a a 30 12 0 . However, in case of 164 Dy-162 Dy mixture, we observe that the system is almost miscible for small inter-species scattering lengths, becoming immiscible as we increase the inter-species scattering length. Therefore, we clearly noticed from these two figures, the immiscibility increasing when the inter-species scattering length become dominant in respect to the dipolar interactions.

Summary and conclusion
Motivated by recent experimental studies with dipolar systems, we focus our study in the two coupled mixtures given by 168 Er-164 Dy and 164 Dy-162 Dy. In the present work, our approach was to evidence the miscibility properties of a coupled condensate with two different species of atoms, having contact or dipolar pairwise interaction between them. First, we did an investigation related to the stability of dipolar coupled systems, as one varies the trap-aspect ratio and the number of atoms of both species. For realistic number of atoms in a mixed BEC system, we found necessary to consider pancake-type configurations (l > 1) for the coupled condensates, when the dominant nonlinear interaction is dipolar. In order to be more complete on the characterization of miscibility, we have extended our study to non-dipolar systems in cigar-type configurations, where it was possible to point out strong immiscibility for the coupled system.
For the study of miscibility, we extended an approach for the critical conditions of homogeneous coupled systems confined in hard-wall barriers. We observed that the critical MIT conditions remain unaffected by the dipolar interactions, once all the parameters of the previous definition are rescaled by incorporating the dipolar ones. In order to measure the miscibility of a more general confined coupled system, a relevant parameter η was defined in terms of the two-component densities. For a general system, this parameter is shown to be adequate to verify the miscibility of a coupled system than the usually simplified criterium obtained for homogeneous systems from energy consideration, where the kinetic energy is ignored.
By studying the miscibility for pure-dipolar coupled system (zero two-body contact interactions), we first compare the properties of the two mixtures given by 168 Er-164 Dy and 164 Dy-162 Dy. Besides the small mass difference between the atoms in both this two mixtures, one should notice that the main difference in their respective dipolar interactions ( ) a ij d is due to the differences in the magnetic moment dipole of erbium and dysprosium atoms. As verified, the two mixtures have quite different miscibility behavior, with 164 Dy-162 Dy being completely miscible (h = 0.99) and 168 Er-164 Dy partially miscible (h = 0.77), when we fix to the same values the other parameters (trap-aspect ratio and number of atoms). Such behavior is clearly due to the interspecies dipolar strength in comparison with the intra-species one. In this pure-dipolar case, we are also pointing out some non-trivial biconcave-shaped structures, with local minimum states for both components of the coupled system and with the emergence of density oscillations (manifested by a few peaks), when the system is near the stability border. This behavior is verified for both coupled mixtures that we have examined, with no direct relation with the miscibility of the species. Already reported in [39], this kind of nontrivial biconcave configuration has also being observed for single component dipolar systems, with the density fluctuations attributed to roton mode of the condensates. The role of the trap aspect ratio and inter-species contact interaction for the miscible-immiscible phase transition was studied for different configurations, from nondipolar to pure dipolar systems.
The present results on the miscibility are expected to be quite relevant in studies with dipolar and nondipolar coupled systems, in order to set parameters in experimental realizations. Within this perspective, when considering dipolar systems, we choose atomic mixtures with large repulsive dipolar strengths that are being investigated in BEC laboratories. In a possible straightforward extension of this work, we could alter the confining conditions of both condensates, with their center being separated by some distance, or by using different aspect ratios. As another perspective for future developments, we can mention studies of systems under rotations, where rich vortex structures will emerge, by following recent interest in the subject that can be traced from [57] and references therein.