 $a_{aa}=0.05a_{ho}$.

The variational wavefunction given in Eq. (\[psivar\]) is best-suited to describe the case of $a_{aa}=0$. For ${a_{aa}\neq}0$ the atom cloud deviates from a gaussian, and for strong enough interactions a gaussian form for the atom cloud is a poor approximation. Consequently, as the atom-atom interactions increase, we find that the simplistic variational wave function given in Eq. (\[psivar\]) cannot describe the tightly-localized impurity at the trap center prior to the onset of collapse. Nonetheless, for the case $a_{aa}=0$ discussed above, the variational treatment reproduces the key features of the full self-consistent solution and provides us with further insights. In particular, the form of the variational wavefunction, Eq. (\[psivar\]), is useful in visualizing how the condensate develops features characterized by a length scale much smaller than the oscillator length. Furthermore, the disappearance of the local minimum as the variational energy becomes unbounded from below is another indication, along with the disappearance of the self-consistent solutions, of the collapse of the condensate.

Discussion and conclusion {#conclusion}
=========================

Sections \[meanfield\] and \[variational\] discuss the behaviors of a single neutral impurity, which feels no external confining potential, immersed in a trapped condensate. If the impurity feels an external trapping potential with angular frequency $\omega_i$, which might be the case in an experiment (see below), region (A) in Fig. 1 is absent, i.e., the impurity is always localized due to the presence of the external potential. Assuming that the impurity trapping potential has a characteristic length which is larger than roughly the condensate healing length (for large enough number of atoms), the comparatively strong impurity localization prior to collapse and the crossover from region (B) to region (C) in Fig. 1 are, however, nearly unaltered. For example, for $N=10^4$, $a_{aa}=0.005a_{ho}$, $m_i=m_a$ and $\omega_a=\omega_i$, the onset of collapse occurs at the same critical value of $a_{ai,c2}\approx-0.062a_{ho}$ that we find without impurity trap.

Impurity-doped condensates can be realized experimentally with present-day technology [@ch
$ a_{aa}=0.05a_{ho}$. 

 The variational wavefunction given in Eq.   (\[psivar\ ]) is best - suited to describe the subject of $ a_{aa}=0$. For $ { a_{aa}\neq}0 $ the atom cloud deviate from a gaussian, and for strong enough interactions a gaussian shape for the atom cloud is a poor estimate. Consequently, as the atom - atom interactions increase, we find that the simplistic variational wave affair given in Eq.   (\[psivar\ ]) cannot describe the tightly - set impurity at the trap center prior to the attack of collapse. Nonetheless, for the case $ a_{aa}=0 $ discussed above, the variational treatment reproduces the key feature of the full self - consistent solution and provides us with further insight. In particular, the form of the variational wavefunction, Eq.   (\[psivar\ ]), is useful in visualizing how the condensation develops feature characterized by a length scale a lot smaller than the oscillator length. Furthermore, the disappearance of the local minimum as the variational energy becomes unbounded from downstairs is another indication, along with the disappearance of the self - consistent solutions, of the collapse of the condensate. 

 Discussion and conclusion { # conclusion } 
 = = = = = = = = = = = = = = = = = = = = = = = = = 

 Sections   \[meanfield\ ] and \[variational\ ] discuss the behavior of a single neutral impurity, which feel no external confining electric potential, plunge in a trapped condensate. If the impurity feel an external trapping potential with angular frequency $ \omega_i$, which might be the case in an experiment (see downstairs), region (A) in Fig.   1 is absent, i.e., the impurity is always localized due to the presence of the external electric potential. Assuming that the impurity trapping potential has a characteristic length which is larger than roughly the condensate healing duration (for large enough number of atoms), the comparatively hard impurity localization prior to collapse and the crossover from region (B) to region (C) in Fig.   1 are, however, nearly unaltered. For example, for $ N=10 ^ 4 $, $ a_{aa}=0.005a_{ho}$, $ m_i = m_a$ and $ \omega_a=\omega_i$, the attack of collapse occurs at the same critical value of $ a_{ai, c2}\approx-0.062a_{ho}$ that we find without impurity trap. 

 Impurity - dope condensates can be recognize experimentally with present - day engineering   [ @ch
 $a_{aw}=0.05a_{ho}$.

The variational wavenunction given iu Eq. (\[psitar\]) is gest-suitdd to describe the case of $a_{ea}=0$. Fir ${a_{aq}\neq}0$ the atom cloud dexiates frlm a gauwsiai, and for strong enough luteradbions e gaussian form for the adom cloud is a puox approximation. Consequently, as the wtom-atok lnteractions igcresfe, ws find that the simplistic variatiknal wane function given in Eq. (\[psivar\]) cannot descrihe tje tightly-localizef impurity qt trw trap centef prior to the onset or collapse. Nonetheless, for the zase $c_{aa}=0$ discussgb ablee, the varietionaj treatment vvproducas the ley features on the fuol self-consistent solntion and provides uf with fustker insights. In partixuoar, tve fmrm ud tfe baciafional watefunction, Sq. (\[psivar\]), is useful in visualizond how the condsnsate dqvelops features characterized by a lencth scale much smaller thab the oscillator lengjh. Furtheriore, the disappearance of the local minimum as tha varmagiouql endegj becomes unbounded from below is another indydauiok, along with the disappearamcf ps the self-conristenc sklutions, of the coplapse jf thw condenswte.

Doscussion and conclusion {#cobclusion}
=========================

Sectpons \[neanfield\] and \[varictional\] discbss thg behafiors of a single neutrcl imphrity, which feels no dxternal confinivg iotettial, immersed in a trappeq condensete. Iy the imourijy feelf an exterjal tvdpping potential wlth augulas frequencj $\omega_i$, which might be the casx in an experomant (see beljw), renion (A) in Fig. 1 if absent, i.e., thg impuritv is auways locamized dne to the prqsence of the gxternal poteitial. Asstmint thqt the kopurity trappimg potentpaj haw a characteristic levfth which is laxyet than roughly yhe cogdvnsete hqdling length (for lafbe enuugh number uf ayoms), the comparativeny sfrong impurity locsllzation ptior to cjllapse and tne crossover from gegioi (B) to regipn (S) in Fig. 1 are, however, nearly unzltered. Flr cxample, for $N=10^4$, $a_{aa}=0.005q_{ho}$, $m_i=m_a$ and $\omega_a=\omega_i$, the onset of collapse occnrs at the same criticao value of $a_{ai,c2}\apprpw-0.062a_{ho}$ that we find without hmpurity trap.

Impuritt-doped condensatex can be realized expsrimendally with present-day technology [@ch
$a_{aa}=0.05a_{ho}$. The variational wavefunction given in Eq. best-suited describe the of $a_{aa}=0$. For from gaussian, and for enough interactions a form for the atom cloud is poor approximation. Consequently, as the atom-atom interactions increase, we find that the simplistic wave function given in Eq. (\[psivar\]) cannot describe the tightly-localized impurity at the center to onset collapse. Nonetheless, for the case $a_{aa}=0$ discussed above, the variational treatment reproduces the key features of full self-consistent solution and provides us with further In particular, the form the variational wavefunction, Eq. (\[psivar\]), useful visualizing how condensate features by a length much smaller than the oscillator length. Furthermore, the disappearance of the local minimum as the variational energy unbounded from another indication, with disappearance the self-consistent solutions, collapse of the condensate. Discussion and Sections \[meanfield\] and \[variational\] discuss the behaviors of single neutral which feels no external confining potential, in a trapped condensate. If the impurity feels external trapping potential with angular frequency $\omega_i$, which might be the case in an experiment region (A) in Fig. is absent, i.e., impurity always due the presence the external potential. Assuming that the impurity trapping potential has a length which is larger than roughly the condensate healing length enough of atoms), the strong impurity localization prior collapse the crossover from region region in however, unaltered. example, for $N=10^4$, $a_{aa}=0.005a_{ho}$, and $\omega_a=\omega_i$, the onset of occurs at the same we find without impurity trap. Impurity-doped condensates can realized experimentally with present-day technology [@ch
 $a_{aa}=0.05a_{ho}$.

The variational wavefUnction givEn in EQ. (\[psIvaR\]) iS besT-suiTed to describe tHE casE of $a_{aa}=0$. For ${a_{aa}\neq}0$ the atom Cloud DeVIateS FrOm a gaUssian, aND fOR StrOnG eNouGh INtEractIonS a gaussIan form for The AtOm cloud is a poOR aPproximatiOn. COnsequently, aS thE atom-aToM inTEractIonS incrEase, we FInd thaT the simplIsTIc variATional wAVE fUnctIon given in Eq. (\[psivaR\]) CaNNot describe the TightlY-lOCaLIZed ImpUrity at the TrAp cenTEr prior TO tHE ONseT Of collapse. NonEtheless, for THe cAse $a_{aa}=0$ DiScuSSed aboVe, the VaRIatIonal treatmEnt rEproduces The key FEatures OF the fulL self-cOnsIstEnt sOLuTiOn aNd PRovIDeS us WIth Further iNsIgHts. In PartICULAr, thE foRm of The vaRiational waveFunCtioN, eq. (\[pSivar\]), Is useFul iN vIsualIzing hOw the CoNdensate developS feaTures charActErIzeD bY a lenGTh scalE muCh sMaller tHan the oSCilLaTOR LeNgth. Furthermore, the DiSAPpEarance oF the loCAl MiNImum as thE vAriAtioNAL enerGy beCOmEs unbounDed froM BeLoW is anotHeR indicAtIon, AloNg witH The dIsappeArance of The seLF-consistent solUTions, of the colLApSE Of THe coNdeNsate.

DiscusSion ANd coNcluSIoN {#coNClusiOn}
=========================

SecTiONs \[MEanfield\] and \[variatioNaL\] discuSs the Behaviors of a sIngle neutrAL IMpurity, wHich FEeLS no external conFininG potential, IMmersed iN a traPped condEnsate. If tHE Impurity FeeLs aN exTerNAL tRapping potentIAL witH aNgular fReqUency $\omEga_I$, whIch MigHt Be the case In an expeRiMeNt (SeE beLow), reGIon (A) in FiG. 1 iS abSeNt, i.E., the iMPurity Is alwAys lOcAlIZed Due to thE PrESEnce Of ThE extErnAl PotenTial. aSsuMing thaT the impurIty TRappInG pOtentiaL has a characteRiStic length WhIch Is largER Than rougHly the condensate healing LEngth (foR laRge enOugh Number of aTomS), the coMpaRAtivelY stronG impuRiTy lOCAlizaTIOn PriOr To collapse AND thE crosSoVer fRom regiOn (B) to region (C) in Fig. 1 aRE, hoWever, nearly unAltEred. fOR eXamPLe, FOr $N=10^4$, $A_{aA}=0.005A_{ho}$, $M_I=M_a$ and $\omega_a=\omegA_i$, the onset Of COlLapse occurS At tHe Same criTical vaLue of $A_{Ai,c2}\apprOx-0.062a_{ho}$ that We find witHoUt imPURitY trap.

ImpurIty-doped CondensatES can bE ReAlizeD exPerimeNtAllY with PresenT-Day TechnOlogy [@cH
 $a_{aa}=0.05a_{ho}$.

The  variation al wa vef unc ti on g iven  in Eq. (\[psi v ar\] ) is best-suited to de scrib et he c a se  of $ a_{aa}= 0 $. F or${ a_ {aa }\ n eq }0$ t heatom cl oud deviat esfr om a gaussia n ,and for st ron g enough int era ctions  a  ga u ssian  fo rm fo r thea tom cl oud is apo o r appr o ximatio n .  C onse quently, as the a t om - atom interacti ons in cr e as e ,  we  fi nd that th esimpl i stic va r ia t i o nal wave function  given in E q . ( \[psiv ar \]) cannot  desc ri b e t he tightly- loca lized imp uritya t the t r ap cent er pri ortotheo ns et  of  c o lla p se . N o net heless,fo rthe c ase$ a _ { aa}= 0$disc ussed  above, the v ari atio n altreat mentrepr od ucesthe ke y fea tu res of the full  sel f-consist ent  s olu ti on an d  provi des  us  with f urtheri nsi gh t s .  I n particular, thefo r m  o f the va riatio n al  w a vefuncti on , E q. ( \ [ psiva r\]) ,  i s useful  in vi s ua li zing ho wthe co nd ens ate  deve l opsfeatur es chara cteri z ed by a length scale much sm a ll e r  t h an t heoscillatorleng t h. F urth e rm ore ,  thedisap pe a ra n ce of the local min im um asthe v ariational en ergy becom e s unbounde d fr o mb elow is anothe r ind ication, a l ong with  thedisappea rance oft h e self-c ons ist ent  so l u ti ons, of the c o l laps eof thecon densate .
Dis cus sio nand concl usion {# co nc lu si on} 
==== = ======== == === == === ==

S e ctions  \[me anfi el d\ ]  an d \[var i at i o nal\ ]di scus s t he  beha vior s  of  a sing le neutra l i m puri ty ,which f eels no exter na l confinin gpot ential , immersed  in a trapped condensat e . If th e i mpuri ty f eels an e xte rnal t rap p ing po tentia l wit hang u l ar fr e q ue ncy  $ \omega_i$, w hic h mig ht  bethe cas e in an experiment (se e below), reg ion  (A) i nFig .  1 isab s ent , i.e., the impur ity is alw ay s  l ocalized d u e t othe pre sence o f the externa l potenti al. Assum in g th a t  th e impurity  trappin g potenti a l has achara cte ristic  l eng th wh ich is lar ger t han ro ug hly th e con de nsate he aling length (for large  enoug h num ber  of atoms ),t hecomparati vely  strong im pur ity  loca liz a tionprio r  t o c o llaps e an d  the cros s ov erf r om  region (B) t o  re gion(C) in Fig . 1are, however, nea r ly unaltered.Fore x amp le, for$N =10^4$, $a_{aa }=0 .0 0 5 a_{ho}$,  $ m_i=m_a$ an d $\omeg a_ a =\ome ga_i$,  the o nset of c ol l apse o ccur s a t the sam e c ri t ical va lu eo f $a_{ ai,c 2} \appro x-0.06 2 a_{h o } $ that we find w ithou t impur i tytrap. 
Impurit y -dop ed condens ates can be  reali zedexper imental ly  withpre se nt-day tec h nology [@ ch
 $a_{aa}=0.05a_{ho}$.

The_variational wavefunction_given in Eq. (\[psivar\]) is_best-suited to_describe_the case_of_$a_{aa}=0$. For ${a_{aa}\neq}0$_the atom cloud_deviates from a gaussian,_and for strong_enough_interactions a gaussian form for the atom cloud is a poor approximation. Consequently, as_the_atom-atom interactions_increase,_we_find that the simplistic variational_wave function given in Eq. (\[psivar\])_cannot describe_the tightly-localized impurity at the trap center prior_to_the onset of_collapse. Nonetheless, for the case $a_{aa}=0$ discussed above, the_variational treatment reproduces the key features_of the full_self-consistent_solution_and provides us with_further insights. In particular, the form_of the variational wavefunction, Eq. (\[psivar\]), is_useful in visualizing how the condensate develops_features characterized by a length scale_much smaller than the oscillator_length. Furthermore,_the disappearance of the local_minimum as the_variational energy_becomes unbounded from_below is another indication, along with_the disappearance of_the self-consistent solutions, of the collapse_of_the condensate.

Discussion and_conclusion_{#conclusion}
=========================

Sections \[meanfield\]_and \[variational\]_discuss the behaviors_of_a single_neutral_impurity, which feels no external confining_potential,_immersed in a trapped condensate. If the_impurity feels an external_trapping_potential with angular frequency_$\omega_i$, which might be the_case in an experiment (see below),_region (A)_in Fig. 1_is absent, i.e., the impurity is always localized due to the_presence of the external potential. Assuming_that the impurity trapping_potential has_a_characteristic length which_is_larger than_roughly the condensate healing length (for large_enough number_of atoms), the comparatively strong impurity_localization prior to collapse_and_the crossover from region (B) to_region (C) in Fig. 1 are, however,_nearly unaltered. For example, for_$N=10^4$,_$a_{aa}=0.005a_{ho}$,_$m_i=m_a$ and $\omega_a=\omega_i$, the onset_of collapse occurs at the same_critical value of_$a_{ai,c2}\approx-0.062a_{ho}$ that we find without impurity trap.

Impurity-doped_condensates_can be realized experimentally with present-day_technology [@ch