Arxiv:2009.13942V2 [Cond-Mat.Quant-Gas] 30 Sep 2020 Rpee[ Graphene I Eedsrbdnmrcly[ Numerically Described Were Dia Soengah,Booy N Brotc.Atrad,The Afterwards, Optics

arXiv:2009.13942v2 [cond-mat.quant-gas] 30 Sep 2020 rpee[ graphene i eedsrbdnmrcly[ numerically described were dia soengah,booy n brotc.Atrad,the Afterwards, optics. fiber and such fields biology, different oceanography, in confirmed as been have and ago, century a [ sys optics linear different [ in waves water observed including been have tems, they nonl and from interaction, effects ear nontrivial indicates wavepackets calized ato at of control reflection motion. the on quantum flexibility excel- large study with an incidence to is normal experiments lengths wave for Broglie platform degenerate de lent quantum large and with ultracold gases of atomic system the surface, the length, number, wave φ [ soli of potential the on [ appear re- surfaces atoms quantum ultra-cold where system of atomic flections cold to extended are ments [ faces [ mirror such atomic systems, various as in observed been has phenomenon Such t hp hnpoaaigoe ogdsacs n emerge and [ distances, unaltered long collisions over from propagating when shape its ir rwlswtotrahn lsia unn on [ point turning classical bar- a reaching potential without from wells or reflected riers are particles whereby nomenon unu eeto faBs-isencnest rmarapi a from condensate Bose-Einstein a of reflection Quantum ∗ † 13 ( [email protected] [email protected] k oiosaesl-enocn aeakt,wihmaintai which wavepackets, self-reinforcing are Solitons unu eeto sacasclycutrnutv phe- counterintuitive classically a is reflection Quantum .Snesgicn eeto cusudrtecondition, the under occurs reflection significent Since ]. (1 = ) 8 λ ,adS,B7gassrae [ surfaces glass BK7 Si, and ], and , 10 12 /k , ,smcnutrhtrsrcue n hnfilms thin and heterostructures semiconductor ], 2 11 20 ae ewe h retto nl ftesoliton the of o angle length orientation the to the up between drop rates potential of widths the to susceptible nl ewe h nietsltnadterflce soliton reflected the reflectio and quantum soliton is the incident It velocity, the between incident angle of barriers. regime potential certain and a drops potential from reflected ea ihicesn oeta its oee,frabarr a for However, widths. potential increasing with decay it esta h ea egho h atrwv and wave matter the of length decay the than less width ASnmes 37.m 37.k 05.45.Yv 03.75.Kk, 03.75.Lm, numbers: PACS be different the from arisen is variation. systems two in reflectances the ) k dk/dx x esuytednmcbhvo faBs-isencondensate Bose-Einstein a of behavior dynamic the study We swl so h aiu tutr,sc as such structure, various the on as well as ] , s 21 2 = 3 stenra itnefo h tmto atom the from distance normal the is – .INTRODUCTION I. .Tertcly oiosi olna me- nonlinear in solitons Theoretically, ]. 5 s ,nnpru ei [ media nanoporous ], π/λ ∼ 4 3 eed nted rgi wave Broglie de the on depends , SCST egHaw etrfrFnaetlTer,Xa 7 Xian Theory, Fundamental for Center Huanwu Peng NSFC-SPTP hax e aoaoyfrTertclPyisFrontiers, Physics Theoretical for Laboratory Key Shaanxi 14 dλ/dx 2 – colo cec,WhnUiest fTcnlg,Whn43 Wuhan Technology, of University Wuhan Science, of School 1 17 nttt fMdr hsc,NrhetUiest,Xi’an, University, Northwest Physics, Modern of Institute 16 ,pams[ plasmas ], 22 .Teeitneo uhlo- such of existence The ]. s n nltcly[ analytically and ] ≫ oge Wang, Dongmei 1 9 [ .Rcn develop- Recent ]. 9 ,weetelocal the where ], 18 6 , , 7 19 ,ruhsur- rough ], ,adnon- and ], Dtd coe ,2020) 1, October (Dated: 1 23 oXiong, Bo half ] 1 , mic in- 2 θ soliton ]. d n - 0 = δR eetdfo otwls[ walls soft ar they from when dissipation reflected some exist there that shown has work u eeto fsc rgtsltn r icse intensi discussed are [ solitons bright such of reflection tum gaussian the the [ or when surface tential structure solid the cloud inte with the interacts from condensate destroy arisen can vortices, interaction and atomic solitons as re- cloud such atomic condensate, for [ factor diffraction key a and be flection to interaction verified the is reflections, atoms tween single-atom and light the from [ condensate the through [ equations (GP) Gross-Pitaevskii interpreted well-known be the theoretically by can system soliton atomic such ultracold atten Especially, much in attracted have solitons condensate, Bose-Einstein , e.g., Recently arrays. platform uide various [ systems in optical investigated as such also were solitons atomic oio,epcal t tutr,sbett oeta w potential a to subject dark structure, 2D its a of especially influence soliton, the Moreover, presented. been not have potent sharp a from barrier[ t significantly condensate affect the can of soliton probability dark reflection 2D a of orientation the region [ associated addressed been have and dynamics chaotic collision its solitons, matter-wave fringes interference of minima local [ the from and shock, tum [ xeietly atrwv oio ihapaedifferen phase of a with soliton matter-wave a Experimentally, hntevlct ftecnest slw nanra re- abnormal an low, is [ experimentally condensate and observed are the process, flectivity reflection of the velocity on the effect when of significant a reflection has the quantum densate of the velocity incident with the involved Conventionally, condensates. been has work ical and 2, 15 33 28 s i.e., , h aksltnadtedffrneo h reflection the of difference the and soliton dark the f ∗ rbblt ipastecsn ftedeflection the of cosine the displays probability n hw htfrarpdyvrigptniladin and potential varying rapidly a for that shown rsnsa xoeta rn.Ti iceac of discrepancy This trend. exponential an presents nrcn er,agetda feprmna n theoret- and experimental of deal great a years, recent In π e potential, ier air fmte ae ntergo fpotential of region the in waves matter of haviors , , , θ n a Yang Tao and 29 34 a egnrtdraiyb prpit hs imprinting phase appropriate by readily generated be can 40 = R ,dniyegneig[ engineering density ], ]. .Frmte-aedr oios eettheoretic recent a solitons, dark matter-wave For ]. 42 36 ( BC otiigadr oio separately soliton dark a containing (BEC) π/ θ .Hwvr h seta atr o uhinfluence such for factors essential the However, ]. ) , 2 ∼ 37 , l ayn oeta:terl fdark of role the potential: varying dly δR o 2 cos .Det h upiigrbsns fbright of robustness surprising the to Due ]. R ia,702,China 710127, Xi’an, ( s θ ipasoclaigexponential oscillating displays , 117 China 710127, ,3 4, 3, 1, θ ) o oeta drop, potential a For . sisniiefrtepotential the for insensitive is 00 China 0070, 02,China 10127, † 24 35 31 41 ,mnei ls[ films mangetic ], .Tennierectto nthe in excitation nonlinear The ]. , .Alts okas hw that shows also work latest A ]. 32 30 ,a osqec fquan- of consequence a as ], ,sepn ioepotential dipole a sweeping ], 6 38 , R 11 , 39 ( .Mroe,apart Moreover, ]. θ ) .As h quan- the Also ]. 25 is ,adwaveg- and ], 26 l on ell , tion. con- well vely of s 27 po- be- ial he ce r- s, ]. e , 2 the quantum reflection has not yet investigated. Correspond- ingly, the underlying mechanism of a BEC with a dark soliton reflected separately from potential well and potential barrier is still not clear. In this paper, we investigate the dynamic behavior of a 2D disk-shaped Bose-Einstein condensate containing a dark soliton separately reflected from the potential drop and the potential barrier. We try to find the influence of the orientation angle of the dark soliton, the structure of potential and the width of the potential on the quantum reflection rate. Af- terwards, we further used dynamic images to reflect the physical mechanism through numerical simulation. The underlying mechanism of the quantum reflection is discussed intensively.

II. METHOD AND SIMULATION SCHEME

We consider a BEC containing N 23Na atoms, which are FIG. 1. A schematic view of a condensate containing a dark soliton, initially trapped by a harmonic trapping potential of the form whose center-of-mass is away from the potential step (a) and poten- V (x,y,z) = m(ω2x2 + ω2y2 + ω2z2)/2. Here we chose h x y z tial drop (b) with distance ∆x in the x-axis. The top view of the the trap frequencies to be ω = ωx = ωy = 2π 10Hz and × initial density distribution of the condensate under one set of specific ωz =2π 100Hz so that the trap frequenciesin the z direction parameters (θ = π/4) (c) and its correspond phase diagram (d). are much× larger than that in the x and y directions. Hence, the system can be treated as a disk-shaped 2D condensate. The dynamics of such system at zero temperature is well described by the 2D time-dependent Gross-Pitaevskii equation, variation of trap potential, the condensate cloud is driven to 2 move toward the region of lower potential where the potential ~ 2 i~∂tψ = ψ + V (x, y)ψ + g2DN ψ (1) barrier Fig.1(a) and potential drop Fig.1(b) occur, which allow −2m | | us to study the reflection process of the condensate containing where ψ(x,y,t) is the temporal condensate wave function and a dark soliton. Fig.1(c) and (d) show a typical density distri- 4 the number of atoms, N, ischosentobe 1.8 10 in our calcu- bution and phase diagram of the initial state of a condensate × lation. g2D =2√2π~ωzasaz is the 2D couplingconstantwith with θ = π/4 containing a dark soliton, respectively. When the s wave scattering length as = 2.9nm. The radial and ax- the solitons are initially imprinted in the condensate, we can ial oscillation lengths are a0 = ~/mω and az = ~/mωz freely adjust the angle of orientation θ between the dark soli- respectively. The ground state ofp dark soliton is obtainedp by tons and the positive x-axis. To derive the dimensionless GP the imaginary time evolution of GP equation, where a π phase ~ equations, we use x0 = , t0 = 1/ω, E0 = ~ω as unit step is compulsively added into the condensate wave function q mω of length, time, energy, respectively. Note that in the region so that two stable density halves are formed. To get a prevalent of potential drop, the velocity of the condensate wave is much insight aboutthe interaction betweendark soliton and potential larger than the incident velocity so a sufficiently small grid size surface, we consider general potentials with the form, is desired to avoid the distortion arisen from finite size scale. Vh(x, y, 0) x< 0 To explore the role of dark soliton in quantum reflection of V (x, y)=  h x [0, d] (2) BEC from the potentials, it is required that the lifetime of dark  0 x∈ > d soliton should be larger than the timescale of reflection pro-  cess. Since the dark soliton in 2D BEC is generally unstable, where h > 0 and h < 0 are the intensity of the potential step it will eventually decay into vortices due to snake instability. and potential drop, respectively. d is the width of the barrier This indicates that the complete quantum reflection should ap- or well. pear before the distortion of the dark soliton. My latest work We consider two different potential profiles: potential I shows that the lifetime of dark soliton either located in the side [Fig.1 (a)] is the potential barrier of width d; potential II [Fig.1 area of 2D BEC or with quantum noise is still larger than the (b)] is a potential well of width d. The sharp potential choosen timescale of reflection process. Also we note that the stability here have fulfill readily the quantum reflection condition and of 2D dark soliton in trapped BEC is still an interesting topics make the reflection of solitonsprominent. Thusit enablesus to on its own. identify the key physical processes that occur when BEC expe- riences reflection, and thus sharpen the understanding needed for experimental studies to explain BECs with solitons reflected from the more complex semiconductor surface poten- III. RESULTS AND DISCUSSION tial or nanostructure potential. In the middle of the harmonic trap, the center-of-mass of the condensate is placed at the po- As stated in [43, 44], the center-of-mass velocity and in- sition which away from the trap center (∆x, 0). Due to the teratomic interaction of condensates play an important role in 3 the formation of interference and dynamic excitation. These two factors determine the time scales of interference and distortion as well as the competition of two time scales. Here, the quantum reflection of the condensate comes with quantum interference and possible nonlinear excitation, which indicates the center-of-massvelocity and interaction are also critical factors. As the condensate moves toward the trap center, its aver- x age incident velocity at the edge of the potential well or barrier 2 2 is approximately, v¯ ω ∆x + ∆y. Previous research shows ≈ q that at low incident velocity and high density, the interference time of the BEC is extended while its distortion time is reduced significantly, as a result that dark solitons and votices are cre- ated in the process. Although the interaction has little effecton refection rate, it causes the fragmentation of the BEC and in the threshold of experimental measurement, the low reflection rate is obtained [10]. Here, we further study the quantum reflection of BECs containing dark solitons from rapidly varying potentials, and specifically concentrate on the role of structure pattern of the nonlinear excitations in such reflection process. To demonstrate explicitly our findings, we define the reflection rate as

L 0 2 1 2 R(t)= ψ(x,y,t) dxdy (3) Z L Z L N − 2 − 2 | | where L is the system size chosen for calculations. We de- note by Rs the asymptotic reflection probability after the condensate cloud has been completely reflected by the potential. When the height or depth of the potential, h , is enough large, | | the BEC undergoes the elastic reflection with Rs = 1. In principle, the quantum reflection from two different potentials along x direction is only dependent on the velocity in x direction. While, for interacting system , the reflection definitely FIG. 2. (Color online) Dynamics of a BEC containing a dark soliton depend on the velocity in y direction, i.e., ωy∆y. (a). The typical density plots show that the initial matter wave soliton Our latest work shows that when a BEC with a dark soli- with orientation angle θ (i) incidents separately on a potential drop ton reflects from a potential barrier, its asymptotic reflection (ii) and potential barrier (iii). The asymptotic reflection rate Rs(θ) probability, Rs, is highly sensitive to the orientation angle θ of the BEC, reflected separately from potential drop (b) and potential of the soliton [42]. Here we further consider the behavior of step (c). The red curve with squares and the dark curve with circles in such BEC, which is reflected from potential well. We find that (b) are from numerical solutions and correspond to the width of the potential drop d = 4.2x0 and d = 4.7x0, respectively, and identical although Rs is still sensitive to the orientation angle in the sys- parameters: h = −180E0, ∆x = 7.05x0. The line curves are the tem of potential well, the variation of Rs as the function of θ, analytic fitting from the ansatz Rs = R0 + Rse cos 2θ: R0 = 0.419 i.e., R(θ), is stronglydependenton the width of potentialdrop. and Rse = −0.044 (blue dashed line) and R0 = 0.365 and Rse = In Fig.2 (b), for d =4.2x0, when θ varies from 0 to π, Rs for 0.106 (green dotted line). Similarly, the curve with squares in (c) is potential well increases initially, reaches the maximum value from numerical results for d = 5.64x0, h = 40E0, ∆x = 7.05x0, at θ = π/2, and finally decrease to the initial value. While and the solid line is from the ansatz for potential barrier with R0 = for d =4.7x0, the variation of Rs in terms of θ is completely 0.935 and Rse = 0.045 opposite, which is similar to the case of potential barrier. This . is in contrast with the reflection behavior from potential barriers (see Fig.2 (c) and Fig.3 (b)), which is immune against the width of potential barrier. The essential mechanism for such discrepancy will be demonstrated in later characteristic flection probability is proportional to the normal velocity, e.g., dynamics. Rs(θ) sin θ for our system. This indicates that the inter- Through much simulation, we find that for both potential action between∼ reflected component and incident component wells and potential barriers and in a certain regime of inci- of the matter-wave dark soliton may be an essential factor for dent velocity, the asymptotic reflection probability always has the quantum reflection. We note that since the strong diffrac- the cosine form of the deflection angle between incident soli- tion emerges for quantum reflection from potential drops (see ton and reflected soliton, i.e., Rs cos2θ. It is distinguished Fig2.(a)(ii)), the deflection angle does not fulfill 2θ exactly ∼ from the reflection of oblique incident waves, where the re- and corresponding Rs(θ) does not fit as perfectly as potential 4 barriers (see Fig2.(a)(iii) and (c)). for two system display distinct behaviors: ∆Rs shows an oscillating exponential decay into zero for a potential drop while displays an exponential decay into a nonzero constant for a potential barrier. It can be also seen that the quantum reflection of matter-wavedarksoliton from the potentialdropis very sensitive to the small region of potential length and δRs can varies upto 0.25 (1 means a total reflection),while the reflectionfrom the potential barrier is insensitive to small d and δRs decreases monotonically to a constant with increasing d. This indicates that there exist two mechanism for the quantum reflection of the condensate containing a dark soliton [45]: for potential drops, the length of matter-wave dark soliton determines the sensitive width of potential drops; for potential barriers, the decay length of the matter wave in the region of the barrier qualifies the sensitive width of the barrier.

FIG. 3. (Color online) Top: The reflection probability, Rs, as the function of the widths of potential drop (a) and potential barrier (b). Line-and-marker curves show the case of the condensate containing dark solitons with three different orientations for θ = 0 (red circles), θ = π/4 (black squares) and θ = π/2 (blue triangles). Bottom: The variation of the reflection probability, δRs = Rs(θ = π/2) − Rs(θ = 0), as a function of the widths of potential drop (c) and potential barrier (d). The solid curve with squares in (c) is from numerical calculation for h = −180E0 and ∆x = 7.47x0, and dotted curve is the fitting from an oscillating exponential function, δRs = −0.74 exp (−0.32d) sin (3.42d − 9.38). Similarly, the solid curve with squares in (d) is from numerical results of the width of the barrier for h = 40E0, ∆x = 7.05x0, and the dotted line is from the ansatz for potential barrier with δRs = 0.09 exp(−1.67d) − 0.09.

FIG. 4. Top two rows: The temporal density distribution of the BEC Fig.3 shows the influence of the width of rapidly varying containing a θ = π/2 dark soliton and reflected from potential drop of different widths. The depth of the potential well is h = −180E0, potentials in the quantum reflection probability of the mat- and the widths of the potential well in (Ia)-(Ic) and (IIa)-(IIc) are ter wave including a dark soliton. We find that when the d = 4.2x0 and d = 4.7x0, respectively. Bottom two rows: The width of the potential drop is less than the length of the soli- temporal density distribution of the BEC containing a dark soliton ton (or the condensate), i.e., d < 11.34x0, the reflection with orientation angle θ = π/2 and θ = 0 interacting with the po- probabilities, Rs, for different orientation angle are sensitive tential barrier, respectively. Here, its width is d = 0.28x0 and its to d and displays oscillating exponential decay. While for height is h = 100E0. The regime between two solid white lines in d > 11.34x0, Rs are immune against d and approach iden- (Ib) and (IIb) as well as (IIIb) and (IVb) indicate the case of potential tically a nonzero constant. For a potential barrier, since the ki- well and the barrier, respectively. The points of time are t = 0.00t0 netic energyof the incidentwaveis less than the barrierheight, ((Ia), (IIa), (IIIa), (IVa)), t = 1.12t0 ((Ib), (IIb)), t = 1.84t0 ((Ic), 1 2 2 (IIc)), t = 2.40t0 ((IIIb),(IVb)) and t = 3.68t0 ((IIIc), (IVc)). i.e., 2 mω (∆x) x0. Therefore, we argue that tively to the transmitted wave withp velocity 2m(E V )/~. when d is less than the decay length of the matter wave in the Here V > 0 for potential barriers and Vp < 0 for potential− barrier, the orientation of the dark soliton has trivial effect on drops. For potential drops, V < 0 and E V E in our − ≫ Rs while becomes significant when d is larger than the decay simulation, indicating that the velocity of the transmitted wave length. Moreover, Fig.3 (c) and (d) show that the variation of is much larger than the velocity of the incident wave as well reflection probabilities δRs = Rs(θ = π/2) Rs(θ = 0) as the reflected wave. Due to the interatomic interaction and − 5 the absence of confinement potential in the area of drop, the transmitted wave can expand significantly in a short time. As a result, an interference pattern with very small fringe spacings occurs in a relatively large condensate wave packet. This is in contrast with the system of potential barriers, where the small condensate wave packet contains an interference pattern with large fringe spacings. Furthermore, the interatomic interaction causes the change of the reflected component of the condensate in shape, which is very different from the incident component. Also the transmitted component distinguishes from the reflected component, so the matter wave reflection is not a strict specular reflection. To explore the discrepancy of the reflection behavior in the potential drop and potential barrier, we study the time evolution of the condensate density distribution as shown in Fig.4. FIG. 5. (Color online) The reflection probability of the BEC dom- For d = 4.2x0, the transmitted wave and the reflected wave inated by expansion motion as a function of the orientation of the in the region of potential drop (see the bracketed area by two soliton. Different line-and-maker plots indicate the condensate in dif- white lines in Fig.4(Ib)) produce an interference pattern with ferent initial traps: ω = 4ω0 (squared line); ω = 6ω0 (circular line); ω = 8ω0 (diamond line) small fringe spacings. Due to the proper width of potential drop and interatomic interaction, the valley of the matter-wave packet, i.e., a low-density area, emerges in the drop region and increasesthe effective depth of the potentialdrop, thus enhanc- condensate containing a dark soliton with different orienta- ing the quantum reflection. By contrast, for d =4.7x0, the in- tion angle on its quantum reflection. To obtain the expansion- terferenceof the matter wave causes the peak of the wavepack- dominated BEC, we create BEC under different initial traps: age, i.e., a high-densityof area, to occur in theregion of poten- ω = 4ω0, ω = 6ω0 and ω = 8ω0. At time t = 0, we tial drop (see Fig.4 (IIb)), which reduce the effective depth of suddenly displace the harmonic trap by a distance 2.83x0 and the potential dropand therebysuppress the quantumreflection. change its trap frequency to ω0. Thus the condensate expands Therefore, it demonstrate that the interference pattern as well rapidly when it is accelerated toward the potential well. Fig.3 as interatomic interaction are critical factors for the quantum shows the reflection probability of the expansion-dominated reflection of the condensate containing a dark soliton from the BEC, R , with respect to different orientation angle θ. It can potential drop. Moreover, for the reflection of the BEC from s seen clearly that the reflection rate is still sensitive to the orien- potential barrier, the transmitted velocity in the region of bar- tation of the dark soliton; the variation of R is up to 0.05. The rier (see the bracketed areas in Fig.4(IIIb)and (IVb)) are much s behavior of expansion motion on R is opposite to the one of smaller than the incident velocity during the transmission pro- s c.m. motionon R . In the expansion-dominatedprocess, from cess, so there do not exist wave knots and complex interfer- s θ =0 to π/2, R here decreases gradually and approaches the ence pattern. As a result, the reflection behavior is immnune s minimum value at θ = π/2. After θ varies from π/2 to π, R against the small barrier width. But for different orientation s increase gradually. Moreover, the larger expanding speed, the angle of the dark soliton, the quantum reflection for both sys- smaller R for the same θ. tems is still sensitive. The bottom two rows show the reflection s process of the condensate with different orientation from the same potential barrier. When the soliton orientation is par- allel to the barrier surface (Fig.4(IIIa)) and during the reflec- IV. EXPERIMENTAL PARAMETERS AND SUMMARY tion process, the soliton structure tends to prevent the return of the reflected wave by barriers due to its self-reinforcement and With regard to the matter-wave dark soliton with adjustable thus reduces the reflection rate (Fig.4(IIIb)). Note that there orientation angles, it may be prior to produce it by combining exist two reflected waves from the potential barrier. While for the condensate in a double potential well with the phase im- θ =0 the influence of the dark soliton on the reflected waves printing [28, 29, 46]. The sensitive widths of potential wells aretrivialandcloserto thecaseofno soliton(seeFig.2in [42]) and barriers for the dark-soliton condensate are approximately and therebythe reflection rate will increase. Therefore, whena the length of the soliton and the decay length of the matter BEC containing a dark soliton is reflected froma potentialbar- wave in the region of potential barrier, respectively. Accord- rier, the reflection probability for θ = π/2 is smaller than the ing to our system, where a quasi-2D condensate is produced one for θ = 0 (see Fig.2(c)). From Fig.4 ((Ic), (IIc)), one can under the parameters ω = 2π 10Hz, ωz = 2π 100Hz, 4 × × also see that after the condensate contacts with the potential N = 1.8 10 , the length of the dark soliton is 11.34x0 well, the transmitted component of the matter wave contains 75µm and× the decay length is approximately the order of≈ a dark-soliton structure with the life time longer than the re- x0 6.6µm. Curret cold-atom experiments can achieve box- flected component, indicating the stability of the dark soliton like≈ potential with length varing from several micrometers to passing through a potential drop as topological excitation. hundreds micrometers [47, 48], so the potential drops can be Next we explore the effect of the expansion motion of the realized in current experiments. Also we can increase the size 6 of the condensate via the trap potentials and interatomic inter- a certain width up to the order of the width of the condensate, action to fulfill the experimental requirment. this suggests that a potential drop may be a desirable candidate for new devices used to detect the nonlinear excitation. Our re- In summary, we have investigated some crucial factors that sults may provide a general view about the micro-mechanism influence the reflection probability of a dark-soliton conden- of atoms on solid surface potential, for example, semiconduc- sate. In particular, the widths of the potential well and barrier tor surfaces, graphene surfaces or nanostructures. have a greatinfluence on the changeof reflection. Theyare de- termined, respectively, by the length of the dark soliton and the decay length of the matter wave in the region of the potential ACKNOWLEDGMENTS barrier. Furthermore, we find that for a rapidly changing potential, the quantum reflection probability is represented as the We thank Y.-J.Lin for the inspiriting discussion about ex- cosine of the deflection angle between the incident soliton and perimental setup for our system. This work is supported by the reflected soliton within a certain range of incident veloc- the NSFC under grants Nos. 11775178 and 11775177, the ity, namely R(θ) cos2θ. The sensitivity of reflection prob- Major Basic Research Program of Natural Science of shaanxi ability on the orientation∼ of soliton may permit experiments to Province under grant Nos. 2017KCT-12 and 2017ZDJC-32, detect the existence of a dark soliton and probe the possible and the Open Research Fund of Shaanxi key Laboratory for structure of soliton. Since the BEC with a dark-soliton profile Theoretical Physics Frontiers under grant No. SXKLTPF - is much susceptible to potential drops than potential barrier for K20190602.

[1] Harald Friedrich and Johannes Trost, “Working [12] Marcius Silvestre, Tarik P. Cysne, Daniela Szilard, Fe- with WKB waves far from the semiclassical limit,” lipe A. Pinheiro, and Carlos Farina, “Tuning quantum Physics Reports 397, 359 (2004). reflection in graphene with an external magnetic field,” [2] Lew D. Landau, Lew D. Landau, and Lew D. Landau, Phys. Rev. A 100, 033605 (2019). Lehrbuch der theoretischen Physik III. Quantenmechanik (Eu- [13] T. E. Judd, R. G. Scott, A. M. Martin, B. Kaczmarek, ropa Lehrmittel Verlag, 1992). and T. M. Fromhold, “Quantum reflection of ultracold atoms [3] A. Landragin, J.-Y. Courtois, G. Labeyrie, N. Vansteenkiste, from thin films, graphene and semiconductor heterostructures,” C. I. Westbrook, and A. Aspect, “Measurement New Journal of Physics 13, 083020 (2011). of the van der waals force in an atomic mirror,” [14] Thierry Dauxois and Michel Peyrard, “Physics of solitons,” Phys. Rev. Lett. 77, 1464–1467 (1996). (2010), 10.2277/0521854210. [4] B. Segev, R. Côté, and M. G. Raizen, “Quantum reflection from [15] S. L. Cornish, N. G. Parker, A. M. Martin, T. E. Judd, R. G. an atomic mirror,” Phys. Rev. A 56, R3350 (1997). Scott, T. M. Fromhold, and C. S. Adams, “Quantum reflection [5] R. Côté, B. Segev, and M. G. Raizen, “Retardation effects on of bright matter-wave solitons,” Phys.D 238, 1299 (2009). quantum reflection from an evanescent-wave atomic mirror,” [16] Dodd, K Roger, Eilbeck, J. Chris, D John, Morris, and C Hed- Phys. Rev. A 58, 3999 (1998). ley, Solitons and nonlinear wave equations (Academic Press, [6] T. A. Pasquini, M. Saba, G.-B. Jo, Y. Shin, W. Ketterle, 1982). D. E. Pritchard, T. A. Savas, and N. Mulders, “Low [17] Hirota and Ryogo, “Exact Nsoliton solutions of the wave equa- velocity quantum reflection of Bose-Einstein condensates,” tion of long waves in shallow-water and in nonlinear lattices,” Phys. Rev. Lett. 97, 093201 (2006). Journal of Mathematical Physics 14, 810 (1973). [7] G. Dufour, R. Guérout, A. Lambrecht, V. V. Nesvizhevsky, [18] E. A. Kuznetsov, A. M. Rubenchik, and V. E. Za- S. Reynaud, and A. Yu. Voronin, “Quantum re- kharov, “Soliton stability in plasmas and hydrodynamics,” flection of antihydrogen from nanoporous media,” Physics Reports 142, 103–165 (1987). Phys. Rev. A 87, 022506 (2013). [19] Lonngren and K. E., “Soliton experiments in plasmas,” [8] Bum Suk Zhao, H. Christian Schewe, Gerard Meijer, Plasma Physics 25, 943–982 (2000). and Wieland Schöllkopf, “Coherent reflection of He [20] D. Y. Tang, H. Zhang, L. M. Zhao, and X. Wu, “Observation of atom beams from rough surfaces at grazing incidence,” high-order polarization-locked vector solitons in a fiber laser,” Phys. Rev. Lett. 105, 133203 (2010). Phys. Rev. Lett. 101, 153904 (2008). [9] Fujio Shimizu, “Specular reflection of very slow [21] J. D. Gibbon, P. J. Caudrey, R. K. Bullough, and metastable neon atoms from a solid surface,” J. C. Eilbeck, “An N-soliton solution of a nonlinear Phys. Rev. Lett. 86, 987 (2001). optics equation derived by a general inverse method,” [10] R. G. Scott, A. M. Martin, T. M. Fromhold, and F. W. Sheard, Lettere Al Nuovo Cimento 8, 775–779 (1973). “Anomalous quantum reflection of Bose-Einstein condensates [22] Zabusky and J. Norman, “Solitons and bound states from a silicon surface: the role of dynamical excitations,” of the time-independent Schriödinger equation,” Phys. Rev. Lett. 95, 073201 (2005). Physical Review 168, 124–128 (1968). [11] T. A. Pasquini, Y. Shin, C. Sanner, M. Saba, A. Schi- [23] A. C. Newell, The Inverse Scattering Transform (Springer rotzek, D. E. Pritchard, and W. Ketterle, “Quantum Berlin Heidelberg, 1980). reflection from a solid surface at normal incidence,” [24] N. P. Proukakis, N. G. Parker, D. J. Frantzeskakis, Phys. Rev. Lett. 93, 223201 (2004). and C. S. Adams, “Analogies between dark solitons in atomic Bose-Einstein condensates and optical systems,” 7

J. Opt. B- Qun. S. O. 6, S380 (2004). sates due to interatomic interactions and quantum noise,” [25] Mingzhong Wu, Boris A. Kalinikos, and Carl E. Patton, Phys. Rev. A 74, 053605 (2006). “Generation of dark and bright spin wave envelope soliton [36] M. Liu, L. She, H. W. Xiong, and M. S. Zhan, “Quan- trains through self-modulational instability in magnetic films,” tum reflection of vortices in Bose-Einstein condensates,” Phys. Rev. Lett. 93, 157207 (2004). Phys. Rev. A 74, 043619 (2006). [26] S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Seng- [37] A. M. Martin, R. G. Scott, and T. M. Fromhold, “Transmis- stock, A. Sanpera, G. V. Shlyapnikov, and M. Lewen- sion and reflection of Bose-Einstein condensates incident on a stein, “Dark solitons in Bose-Einstein condensates,” gaussian tunnel barrier,” Phys. Rev. A 75, 065602 (2007). Phys. Rev. Lett. 83, 5198 (1999). [38] A. D. Martin, C. S. Adams, and S. A. Gardiner, “Bright matter- [27] J. Denschlag, J. E. Simsarian, D. L. Feder, Charles W. wave soliton collisions in a harmonic trap: regular and chaotic Clark, L. A. Collins, J. Cubizolles, L. Deng, E. W. Ha- dynamics,” Phys. Rev. Lett. 98, 020402 (2007). gley, K. Helmerson, and W. P. Reinhardt, “Generating soli- [39] A. D. Martin, C. S. Adams, and S. A. Gardiner, “Bright solitary- tons by phase engineering of a Bose-Einstein condensate,” matter-wave collisions in a harmonic trap: regimes of solitonlike Science 287, 97 (2000). behavior,” Phys. Rev. A 77, 013620 (2008). [28] J. Denschlag, J. E. Simsarian, D. L. Feder, C. W. Clark, [40] Chaohong Lee and Joachim Brand, “Enhanced quantum reflec- L. A. Collins, J. Cubizolles, L. Deng, E. W. Hagley, tion of matter-wave solitons,” Epl 73, 321 (2005). K. Helmerson, and W. P. Reinhardt, “Generating soli- [41] M. Sciacca, C. F. Barenghi, and N. G. Parker, “Matter-wave tons by phase engineering of a Bose-Einstein condensate,” in dark solitons in boxlike traps,” Phys. Rev. A 95, 013628 (2017). Laser Cooling: Bose-Einstein Condensation and Atom Laser–Proceedings[42] Q. of L. CCAST Cheng, (World W. Laboratory) K. Bai, Y. Workshop Z. Zhang, B. Xiong, and (1999). T. Yang, “Influence of a dark soliton on the reflec- [29] S. Stellmer, C. Becker, P. Soltan-Panahi, E.-M. Richter, tion of a Bose–Einstein condensate by a square barrier,” S. Dörscher, M. Baumert, J. Kronjäger, K. Bongs, and K. Sen- Laser Physics 29, 015501 (2018). gstock, “Collisions of dark solitons in elongated Bose-Einstein [43] Bo Xiong, Tao Yang, and Keith A. Benedict, “Dis- condensates,” Phys. Rev. Lett. 101, 120406 (2008). tortion of interference fringes and the resulting vor- [30] B. P. Anderson, P. C. Haljan, C. A. Regal, D. L. Feder, L. A. tex production of merging Bose-Einstein condensates,” Collins, C. W. Clark, and E. A. Cornell, “Watching dark soli- Phys. Rev. A 88, 043602 (2013). tons decay into vortex rings in a Bose-Einstein condensate,” [44] T. Yang, B. Xiong, and Keith A. Benedict, “Dynamical excita- Phys. Rev. Lett. 86, 2926–2929 (2001). tions in the collision of two-dimensional Bose-Einstein conden- [31] P. Engels and C. Atherton, “Stationary and nonstationary fluid sates,” Phys. Rev. A 87, 023603 (2013). flow of a Bose-Einstein condensate through a penetrable bar- [45] Alexander Jurisch and Jan-Michael Rost, “Anomalous rier,” Phys. Rev. Lett. 99, 160405 (2007). quantum reflection as a quasidynamical damping effect,” [32] J. J. Chang, P. Engels, and M. A. Hoefer, “Formation of disper- Phys. Rev. A 81, 043610 (2010). sive shock waves by merging and splitting Bose-Einstein con- [46] T. P. Meyrath, F. Schreck, J. L. Hanssen, C. S. Chuu, and M. G. densates,” Phys. Rev. Lett. 101, 170404 (2008). Raizen, “A high frequency optical trap for atoms using hermite- [33] G.-B. Jo, J.-H. Choi, C. A. Christensen, T. A. Pasquini, Y.-R. gaussian beams,” Optics Express 13, 2843–2851 (2005). Lee, W. Ketterle, and D. E. Pritchard, “Phase-sensitive recom- [47] Alexander L. Gaunt, Tobias F. Schmidutz, Igor Gotli- bination of two Bose-Einstein condensates on an atom chip,” bovych, Robert P. Smith, and Zoran Hadzibabic, “Bose- Phys. Rev. Lett. 98, 180401 (2007). einstein condensation of atoms in a uniform potential,” [34] A. Weller, J. P. Ronzheimer, C. Gross, J. Esteve, M. K. Phys. Rev. Lett. 110, 200406 (2013). Oberthaler, D. J. Frantzeskakis, G. Theocharis, and [48] T. P. Meyrath, F. Schreck, J. L. Hanssen, C.-S. Chuu, P. G. Kevrekidis, “Experimental observation of os- and M. G. Raizen, “Bose-einstein condensate in a box,” cillating and interacting matter wave dark solitons,” Phys. Rev. A 71, 041604 (2005). Phys. Rev. Lett. 101, 130401 (2008). [35] R. G. Scott, D. A. W. Hutchinson, and C. W. Gar- diner, “Disruption of reflecting Bose-Einstein conden-