Search

We study the stability and characteristics of two-dimensional (2D) quasi-isotropic quantum droplets (QDs) of fundamental and vortex types, formed by binary Bose-Einstein condensate with magnetic quadrupole-quadrupole interactions (MQQIs). The magnetic quadrupoles are built as pairs of dipoles and antidipoles polarized along the x-axis. The MQQIs are induced by applying an external magnetic field that varies along the x-axis. The system is modeled by the Gross-Pitaevskii equations including the MQQIs and Lee-Huang-Yang correction to the mean-field approximation. Stable 2D fundamental QDs and quasi-isotropic vortex QDs with topological charges S<4 are produced by means of the imaginary-time-integration method for configurations with the quadrupoles polarized parallel to the systems two-dimensional plane. Effects of the norm and MQQI strength on the QDs are studied in detail. Some results, including an accurate prediction of the effective area, chemical potential, and peak density of QDs, are obtained in an analytical form by means of the Thomas-Fermi approximation. Collisions between moving QDs are studied by means of systematic simulations., Comment: 11pages, 10 figures,to be published in New Journal of Physics

We investigate vortex states in Bose-Einstein condensates under the combined action of the spin-orbit coupling (SOC), gradient magnetic field, and harmonic-oscillator trapping potential. The linear version of the system is solved exactly. Through the linear-spectrum analysis, we find that, varying the SOC strength and magnetic-field gradient, one can perform energy-level inversion. With suitable parameters, initial higher-order vortex states can be made the ground state (GS). The nonlinear system is solved numerically, revealing that the results are consistent with the linear predictions in the case of repulsive inter-component interactions. On the other hand, inter-component attraction creates the GS in the form of mixed-mode states in a vicinity of the GS phase-transition points. The spin texture of both vortex- and mixed-mode GSs reveals that they feature the structure of 2D (baby) skyrmions., Comment: to be published in Physical Review A

We explore the defect formation and universally critical dynamics in two-dimensional (2D) two-component Bose-Einstein condensates(BECs) subjected to two types of potential traps: a homogeneous trap and a harmonic trap.We focus on the non-equilibrium universal dynamics of the miscible-immiscible phase transition with both linear and nonlinear quenching types.Although there exists spatial independence of the critical point, we find that the inhomogeneity of trap doesn't affect the phase transition of system and the critical exponents can still be explained by the homogeneous Kibble-Zurek mechanism. By analyzing the Bogoliubov excitations, we establish a power-law relationship between the domain correlation length, the phase transition delay, and the quench time.Furthermore, through real-time simulations of phase transition dynamics, the formation of domain defects and the delay of phase transition in non-equilibrium dynamics are demonstrated, along with the corresponding universal scaling of correlation length and phase transition delay for various quench time and quench coefficients, which align well with our analytical predictions.Our study confirms that the universality class of two-component BECs remains unaffected by dimensionality, while the larger nonlinear coefficients effectively suppress non-adiabatic excitations, offering a novel perspective for addressing adiabatic evolution., Comment: 9 pages, 8 figures

We construct strongly anisotropic quantum droplets with embedded vorticity in the 3D space, with mutually perpendicular vortex axis and polarization of atomic magnetic moments. Stability of these anisotropic vortex quantum droplets (AVQDs) is verified by means of systematic simulations. Their stability area is identified in the parametric plane of the total atom number and scattering length of the contact interactions. We also construct vortex-antivortex-vortex bound states and find their stability region in the parameter space. The application of a torque perpendicular to the vorticity axis gives rise to robust intrinsic oscillations or rotation of the AVQDs. The effect of three-body losses on the AVQD stability is considered too. The results show that the AVQDs can retain the topological structure (vorticity) for a sufficiently long time if the scattering length exceeds a critical value., Comment: 7 pages, 6 figures, and 80 References, Physical Review Letters, in press

We study the stability and characteristics of two-dimensional (2D) circular quantum droplets (QDs) with embedded hidden vorticity (HV), i.e., opposite angular momenta in two components, formed by binary Bose-Einstein condensates (BECs) trapped in a radially-periodic potential. The system is modeled by the Gross-Pitaevskii (GP) equations with the Lee-Huang-Yang (LHY) terms, which represent the higher-order self-repulsion induced by quantum fluctuations around the meanfield state, and a potential which is a periodic function of the radial coordinate. Ring-shaped QDs with high winding numbers (WNs) of the HV type, which are trapped in particular circular troughs of the radial potential, are produced by means of the imaginary-time-integration method. Effects of the depth and period of the potential on these QD states are studied. The trapping capacity of individual circular troughs is identified. Stable compound states in the form of nested multiring patterns are constructed too, including ones with WNs of opposite signs. The stably coexisting ringshaped QDs with different WNs can be used for the design of BEC-based data-storage schemes., Comment: 14 pages, 6 figures,to be published in Physical Review E

A new scheme for producing semidiscrete self-trapped vortices (\textquotedblleft swirling photon droplets\textquotedblright ) in photonic crystals with competing quadratic ($\chi ^{(2)}$) and self-defocusing cubic ($\chi ^{(3)}$) nonlinearities is proposed. The photonic crystal is designed with a striped structure, in the form of spatially periodic modulation of the $\chi ^{(2)}$ susceptibility, which is imposed by the quasi-phase-matching technique. Unlike previous realizations of semidiscrete optical modes in composite media, built as combinations of continuous and arrayed discrete waveguides, the semidiscrete vortex droplets are produced here in the fully continuous medium. This work reveals that the system supports two types of semidiscrete vortex droplets, \textit{viz}., onsite- and intersite-centered ones, which feature, respectively, odd and even numbers of stripes, $\mathcal{N}$. Stability areas for the states with different values of $\mathcal{N}$ are identified in the system's parameter space. Some stability areas overlap with each others, giving rise to multistability of states with different $\mathcal{N}$. The coexisting states are mutually degenerate, featuring equal values of the Hamiltonian and propagation constant. An experimental scheme to realize the droplets is outlined, suggesting new possibilities for the long-distance transmission of structured light carrying orbital angular momentum in nonlinear media., Comment: 10 pages, 7 figures, and 76 references

In this study, we investigate the formation of domain defects and the universal critical real-time dynamics in a two-component Bose-Einstein condensate with nonlinear quenching across the miscible-immiscible phase transition. By analyzing the Bogoliubov excitations, we obtain the power-law relations among the defect density, the phase transition delay and the quench time near the phase transition. Moreover, by simulating the real-time dynamics across the miscible-immiscible phase transition, we clearly show the formation of domain defects and the delay of the phase transition. Furthermore, we find that the domain defects are suppressed by large nonlinear coefficients and long quench times. To accurately characterize the domain defects, we quantify the defect excitations using the correlation length and the domain number. In addition, by combining the power-law relations between the phase transition delay and the quench time, we extract the critical exponents for different nonlinear coefficients. Our study not only confirms that the critical exponents do not sensitively depend on the nonlinear quenches but also provides a dynamic path toward the suppression of nonadiabatic excitation., Comment: 7 pages, 6 figures

Creation of stable intrinsically anisotropic self-bound states with embedded vorticity is a challenging issue. Previously, no such states in Bose-Einstein condensates (BECs) or other physical settings were known. Dipolar BEC suggests a unique possibility to predict stable anisotropic vortex quantum droplets (AVQDs). We demonstrate that they can be created with the vortex' axis oriented \emph{perpendicular} to the polarization of dipoles. The stability area and characteristics of the AVQDs in the parameter space are revealed by means of analytical and numerical methods. Further, the rotation of the polarizing magnetic field is considered, and the largest angular velocities, up to which spinning AVQDs can follow the rotation in clockwise and anti-clockwise directions, are found. Collisions between moving AVQDs are studied too, demonstrating formation of bound states with a vortex-antivortex-vortex structure. A stability domain for such stationary bound states is identified. Unstable dipolar states, that can be readily implemented by means of phase imprinting, quickly transform into robust AVQDs, which suggests a straightforward possibility for the creation of these states in the experiment., Comment: 6 pages, 6 figures, and 48 references

It is known that stable 2D solitons of the semi-vortex (SV) and mixed-mode (MM) types are maintained by the interplay of the cubic attractive nonlinearity and spin-orbit coupling (SOC) in binary Bose-Einstein condensates. We introduce a double-layer system, in which two binary condensates, stabilized by the SOC, are linearly coupled by tunneling. By means of the numerical methods, it is found that symmetric two-layer solitons undergo the spontaneous-symmetry-breaking (SSB) bifurcation of the subcritical type. The bifurcation produces families of asymmetric 2D solitons, which exist up to the value of the total norm equal to the norm of the Townes solitons, above which the collapse occurs. This situation terminates at a critical value of the inter-layer coupling, beyond which the SSB bifurcation is absent, as the collapse sets in earlier. Symmetric 2D solitons that are destabilized by the SSB demonstrate dynamical symmetry breaking, in combination with intrinsic oscillations of the solitons, or transition to the collapse, if the soliton's norm is sufficiently large. Asymmetric MMs produced by the SSB instability start spontaneous drift, in addition to the intrinsic vibrations. Consideration of moving 2D solitons is a nontrivial problem because SOC breaks the Galilean invariance. It is found that the system supports moving MMs up to a critical value of the velocity, beyond which they undergo delocalization., Comment: 17 pages, 16 figures, to be published in CNSNS

We report solutions for stable compound solitons in a three-dimensional quasi-phase-matched photonic crystal with the quadratic ($\chi ^{(2)}$) nonlinearity. The photonic crystal is introduced with a checkerboard structure, which can be realized by means of the available technology. The solitons are built as four-peak vortex modes of two types, rhombuses and squares (intersite- and onsite-centered self-trapped states, respectively). Their stability areas are identified in the system's parametric space (rhombuses occupy an essentially broader stability domain), while all bright vortex solitons are subject to strong azimuthal instability in uniform $\chi^{(2)}$ media. Possibilities for experimental realization of the solitons are outlined., Comment: 6 pages, 6 figures, 55 references, Physical Review Letters, in press

The formation and dynamics of full vortex gap solitons (FVGSs) is investigated in two-component Bose-Einstein condensates with spin-orbit coupling (SOC), Zeeman splitting (ZS), and competing cubic and quintic nonlinear terms, while the usual kinetic energy is neglected, assuming that it is much smaller than the SOC and ZS terms. Unlike previous SOC system with the cubic-only attractive nonlinearity, in which solely semi-vortices may be stable, with the vorticity carried by a single component, the present system supports stable FVGS states, with the vorticity present in both components (such states are called here full vortex solitons, to stress the difference from the half-vortices). They populate the bandgap in the system's linear spectrum. In the case of the cubic self-attraction and quintic repulsion, stable FVGSs with a positive effective mass exist near the top of the bandgap. On the contrary, the system with cubic self-repulsion and quintic attraction produces stable FVGSs with a negative mass near the bottom of the bandgap. Mobility and collisions of FVGSs with different topological charges are investigated too., Comment: 14 pages,9 figures, 77 references. Communication in Nonlinear Science and Numerical Simulation, in press

We elaborate a scheme of trapping-expulsion management (TEM), in the form of the quadratic potential periodically switching between confinement and expulsion, as a means of stabilization of two-dimensional dynamical states against the backdrop of the critical collapse driven by the cubic self-attraction with strength g. The TEM scheme may be implemented, as spatially or temporally periodic modulations, in optics or BEC, respectively. The consideration is carried out by dint of numerical simulations and variational approximation (VA). In terms of the VA, the dynamics amounts to a nonlinear Ermakov equation, which, in turn, is tantamount to a linear Mathieu equation. Stability boundaries are found as functions of g and parameters of the periodic modulation of the trapping potential. Below the usual collapse threshold, which is known, in the numerical form, as g < 5.85 (in the standard notation), the stability is limited by the onset of the parametric resonance. This stability limit, including the setup with the self-repulsive sign of the cubic term (g < 0), is accurately predicted by the VA. At g > 5.85, the collapse threshold is found with the help of full numerical simulations. The relative increase of the critical value of g above 5.85 is ~ 1.5%, which is a meaningful result, even if its size is small, because the collapse threshold is a universal constant, which is difficult to change., Comment: to be published in Physical Review E

Fano effect arising from interference between two dissipation channels to the radiation continuum enables to tune the photon statistics. Understanding the role of Fano effect and exploiting to achieve strong photon correlations are of both fundamental and applied significance. We present an analytical description of Fano-enhanced photon correlations based on cavity quantum electrodynamics to show that, the Fano effect in atom-cavity systems can improve the degree of antibunching by over four orders of magnitude. The enhancement factors and the optimal conditions are explicitly given, and found related to the Fano parameter $q$. Remarkably, the Fano enhancement manifests robustness against the decoherence processes, and can survive in the weak coupling regime. We expect our work to provide insight in tuning the photon statistics through Fano effect, which offers a new route to enhance the photon correlations, as well as the possibility of generating nonclassical light in a wider diversity of systems without the need of strong light-matter interaction., Comment: 6 pages, 3 figures

We investigate the single-photon blockade (1PB) in the quasichiral regime of atom-photon interaction that mediates via dissipative environment, where the effective atom-photon interaction is asymmetrical but achiral. The synthetic magnetic current in the closed-loop coupling breaks down the reciprocity of atom-photon interaction, resulting in an asymmetrical and even completely unidirectional effective coupling between two selected quantum states. As an example, we couple the single-atom cavity-QED (cQED) system to a strongly dissipative auxiliary cavity. We find that in the quasichiral regime, the unconventional photon blockade (UPB) always incorporates with the conventional photon blockade (CPB) in the condition of maximum chirality. Furthermore, we show that 1PB in the quasichiral regime combines the advantages of UPB and CPB, demonstrating the perfect single-photon purity, higher efficiency, smooth time dynamics as well as lower requirement of modes coupling to achieve UPB. Our work paves the way for 1PB towards practical applications and reveals the intriguing quantum-optics phenomena in the quasichiral light-matter interaction., Comment: 9 pages, 6 figures

When the quasi-phase matching (QPM) parameters of the $\chi^{(2)}$ nonlinear crystal rotate along a closed path, geometric phase will be generated in the signal and idler waves that participate in the nonlinear frequency conversion. In this paper, we study two rotation schemes, full-wedge rotation, and half-wedge rotation, of the QPM parameters in the process of fully nonlinear three-wave mixing. These two schemes can effectively suppress the uncertainty in creating the geometric phase in the nonlinear frequency conversion process when the intensity of the pump is depleted. The finding of this paper provides an avenue toward constant control of the geometric phase in nonlinear optics applications and quantum information processing., Comment: 12 Pages, 37 References, and 7 Figures, published on Optics Express

Quantum compiling aims to construct a quantum circuit V by quantum gates drawn from a native gate alphabet, which is functionally equivalent to the target unitary U. It is a crucial stage for the running of quantum algorithms on noisy intermediate-scale quantum (NISQ) devices. However, the space for structure exploration of quantum circuit is enormous, resulting in the requirement of human expertise, hundreds of experimentations or modifications from existing quantum circuits. In this paper, we propose a variational quantum compiling (VQC) algorithm based on reinforcement learning (RL), in order to automatically design the structure of quantum circuit for VQC with no human intervention. An agent is trained to sequentially select quantum gates from the native gate alphabet and the qubits they act on by double Q-learning with \epsilon-greedy exploration strategy and experience replay. At first, the agent randomly explores a number of quantum circuits with different structures, and then iteratively discovers structures with higher performance on the learning task. Simulation results show that the proposed method can make exact compilations with less quantum gates compared to previous VQC algorithms. It can reduce the errors of quantum algorithms due to decoherence process and gate noise in NISQ devices, and enable quantum algorithms especially for complex algorithms to be executed within coherence time., Comment: 21 pages, 10 figures

We demonstrate a possibility of the creation of stable optical solitons combining one continuous and one discrete coordinate, with embedded vorticity, in an array of planar waveguides with intrinsic cubic-quintic nonlinearity. The same system may be realized in terms of the spatiotemporal light propagation in an array of tunnel-coupled optical fibers with the cubic-quintic nonlinearity. In contrast with zero-vorticity states, semidiscrete vortex solitons do not exist without the quintic term in the nonlinearity. Two types of the solitons, \emph{viz.}, intersite- and onsite-centered ones (IC and OC, respectively), with even and odd numbers $N$ of actually excited sites in the discrete direction, are identified. We consider the modes carrying the embedded vorticity $S=1$ and $2$. In accordance with their symmetry, the vortex solitons of the OC type exhibit an intrinsic core, while the IC solitons with a small $N$ may have a coreless structure. Facilitating their creation in the experiment, the modes reported in the present work may be much more compact states than their counterparts considered in other systems, and they feature strong anisotropy. They can be set in motion in the discrete direction, provided that the coupling constant exceeds a certain minimum value. Collisions between moving vortex solitons are considered too., Comment: 14 pages, 12 figures, 74 References,Published on Advanced Photonics Research

We study the stability of zero-vorticity and vortex lattice quantum droplets (LQDs), which are described by a two-dimensional (2D) Gross-Pitaevskii (GP) equation with a periodic potential and Lee-Huang-Yang (LHY) term. The LQDs are divided in two types: onsite-centered and offsitecentered LQDs, the centers of which are located at the minimum and the maximum of the potential, respectively. The stability areas of these two types of LQDs with different number of sites for zerovorticity and vorticity with S = 1 are given. We found that the u-N relationship of the stable LQDs with a fixed number of sites can violate the Vakhitov-Kolokolov (VK) criterion, which is a necessary stability condition for nonlinear modes with an attractive interaction. Moreover, the u-N relationship shows that two types of vortex LQDs with the same number of sites are degenerated, while the zero-vorticity LQDs are not degenerated. It is worth mentioning that the offsite-centered LQDs with zero-vorticity and vortex LQDs with S = 1 are heterogeneous., Comment: 12 pages, 8 figures, Frontiers of Physics, in press

This brief review summarizes recent theoretical and experimental results which predict and establish the existence of quantum droplets (QDs), i.e., robust two- and three-dimensional (2D and 3D) self-trapped states in Bose-Einstein condensates (BECs), which are stabilized by effective selffirepulsion induced by quantum fluctuations around the mean-field (MF) states [alias the Lee-Huang--Yang (LHY) effect]. The basic models are presented, taking special care of the dimension crossover, 2D -> 3D. Recently reported experimental results, which exhibit stable 3D and quasi-2D QDs in binary BECs, with the inter-component attraction slightly exceeding the MF self-repulsion in each component, and in single-component condensates of atoms carrying permanent magnetic moments, are presented in some detail. The summary of theoretical results is focused, chiefly, on 3D and quasi-2D QDs with embedded vorticity, as the possibility to stabilize such states is a remarkable prediction. Stable vortex states are presented both for QDs in free space, and for singular but physically relevant 2D modes pulled to the center by the inverse-square potential, with the quantum collapse suppressed by the LHY effect., Comment: A brief review article to be published in Frontiers in Physics

We investigate dynamics of two-dimensional chiral solitons of semi-vortex (SV) and mixed-mode (MM) types in spin-orbit-coupled Bose-Einstein condensates with the Manakov nonlinearity, loaded in a dual-core (double-layer) trap. The system supports two novel manifestations of Josephson phenomenology: one in the form of persistent oscillations between SVs or MMs with opposite chiralities in the two cores, and another one demonstrating robust periodic switching (identity oscillations) between SV in one core and MM in the other, provided that the strength of the inter-core coupling exceeds a threshold value. Below the threshold, the system creates composite states, which are asymmetric with respect to the two cores, or suffer the collapse. Robustness of the chirality and identity oscillations against deviations from the Manakov nonlinearity is investigated too. These dynamical regimes are possible only in the nonlinear system. In the linear one, exact stationary and dynamical solutions for SVs and MMs of the Bessel type are found. They sustain Josephson self-oscillations in different modes, with no interconversion between them., Comment: to be published in Physical Review Research

We study discrete solitons in zigzag discrete waveguide arrays with different types of linear mixing between nearest-neighbor and next-nearest-neighbor couplings. The waveguide array is constructed from two layers of one-dimensional (1D) waveguide arrays arranged in zigzag form. If we alternately label the number of waveguides between the two layers, the cross-layer couplings (which couple one waveguide in one layer with two adjacent waveguides in the other layer) construct the nearestneighbor couplings, while the couplings that couple this waveguide with the two nearest-neighbor waveguides in the same layer, i.e., self-layer couplings, contribute the next-nearest-neighbor couplings. Two families of discrete solitons are found when these couplings feature different types of linear mixing. As the total power is increased, a phase transition of the second kind occurs for discrete solitons in one type of setting, which is formed when the nearest-neighbor coupling and next-nearest-neighbor coupling feature positive and negative linear mixing, respectively. The mobilities and collisions of these two families of solitons are discussed systematically throughout the paper, revealing that the width of the soliton plays an important role in its, Comment: 8 pages, 8 figures,to be published in Physics Letters A

In a nonlinear three-wave mixing process, the interacting waves can accumulate an adiabatic geometric phase (AGP) if the nonlinear coefficient of the crystal is modulated in a proper manner along the nonlinear crystal. This concept was studied so far only for the case in which the pump wave is much stronger than the two other waves, hence can be assumed to be constant. Here we extend this analysis for the fully nonlinear process, in which all three waves can be depleted and we show that the sign and magnitude of the AGP can be controlled by the period, phase and duty cycle of the nonlinear modulation pattern. In this fully nonlinear interaction, all the states of the system can be mapped onto a closed surface. Specifically, we study a process in which the eigenstate of the system follows a circular rotation on the surface. Our analysis reveals that the AGP equals to the difference between the total phase accumulated along the circular trajectory and that along its vertical projection, which is universal for the undepleted (linear) and depleted (nonlinear) cases. Moreover, the analysis reveals that the AGPs in the processes of sum-frequency generation and difference-frequency generation have opposite chirality. Finally, we utilize the AGP in the fully nonlinear case for splitting the beam into different diffraction orders in the far field., Comment: 9 pages, 36 references, 7 figures. Published on Physical Review A

We consider possibilities to grasp and drag one-dimensional solitons in two-component Bose- Einstein condensates (BECs), under the action of gravity, by tweezers induced by spatially confined spin-orbit (SO) coupling applied to the BEC, with the help of focused laser illumination. Solitons of two types are considered, semi-dipoles and mixed modes. We find critical values of the gravity force, up to which the solitons may be held or transferred by the tweezers. The dependence of the critical force on the magnitude and spatial extension of the localized SO interaction, as well as on the solitons norm and speed (in the transfer regime), are systematically studied by means of numerical methods, and analytically with the help of a quasi-particle approximation for the soliton. In particular, a noteworthy finding is that the critical gravity force increases with the increase of the transfer speed (i.e., moving solitons are more robust than quiescent ones). Nonstationary regimes are addressed too, by considering abrupt application of gravity to solitons created in the weightless setting. In that case, solitons feature damped shuttle motion, provided that the gravity force does not exceed a dynamical critical value, which is smaller than its static counterpart. The results may help to design gravimeters based on ultracold atoms., Comment: 13 pages, 9 figures, to be published in New Journal of Physics

It was recently found that, under the action of the spin-orbit coupling (SOC) and Zeeman splitting (ZS), binary BEC with intrinsic cubic nonlinearity supports families of gap solitons, provided that the kinetic energy is negligible in comparison with the SOC and ZS terms. We demonstrate that, also under the action of SOC and ZS, a similar setting may be introduced for BEC with two components representing different atomic states, resonantly coupled by microwave radiation, while the Poisson equation accounts for the feedback of the two-component atomic wave function onto the radiation. The microwave-mediated interaction induces an effective nonlinear trapping potential, which strongly affects the purport of the linear spectrum in this system. As a result, families of both gap and embedded solitons (those overlapping with the continuous spectrum) are found, being chiefly stable. The shape of the solitons features exact or broken skew symmetry. In addition to fundamental solitons (whose shape may or may not include a node), a family of dipole solitons is constructed too, which are even more stable than their fundamental counterparts. A nontrivial stability area is identified for moving solitons in the present system, which lacks Galilean invariance. Colliding solitons merge into a single one., Comment: to be published in Physical Review A

Optical diffraction limit has been a long-term scientific issue since Ernst Abbe first introduced the concept in 1873. It is a constraint on the smallest light spot that can be achieved. Substantial effort has been invested in the past decade to beat this limit by exploiting evanescent waves. But this method encounters serious near-field limitations. A more promising route to breaking the constraint is to explore optical superoscillation in the far field with engineered metamaterials. However, these particular structures involve with very complicated optimization-based design that requires precisely tailoring the interference of propagating waves with low spatial frequency. To overcome these limitations, here we explore a new approach based on the two-hundred-year-old discovery: Possion-Arago spots. We show for the first time that by using a single disc, constructive interference of propagating waves with high-spatial-frequency wavevectors can be realized, generating a diffraction-unlimited localized Possion-Arago spot with achievable size down to $\lambda/$20. Actually, such an element permits creation of an ultra-long nearly nondiffracting superoscillatory needle with appreciable field of view. This easy-to-fabrication element provides a promising route to overcome the diffraction limit, thus might open new avenues to exploit various applications in different fields., Comment: The paper includes 4 figures

We consider a binary bosonic condensate with weak mean-field (MF) residual repulsion, loaded in an array of nearly one-dimensional traps coupled by transverse hopping. With the MF force balanced by the effectively one-dimensional attraction, induced in each trap by the Lee-Hung-Yang correction (produced by quantum fluctuations around the MF state), stable onsite-centered and intersite-centered semi-discrete quantum droplets (QDs) emerge in the array, as fundamental ones and self-trapped vortices, with winding numbers, at least, up to 5, in both tightly-bound and quasi-continuum forms. The application of a relatively strong trapping potential leads to squeezing transitions, which increase the number of sites in fundamental QDs, and eventually replace vortex modes by fundamental or dipole ones. The results provide the first realization of stable semi-discrete vortex QDs, including ones with multiple vorticity., Comment: 7 pages, 4 figures, and 75 references, Published on Physical Review Letters. See: https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.123.133901

We report detailed investigation of the existence and stability of mixed and demixed modes in binary atomic Bose-Einstein condensates with repulsive interactions in a ring-trap geometry. The stability of such states is examined through eigenvalue spectra for small perturbations, produced by the Bogoliubov-de Gennes equations, and directly verified by simulations based on the coupled Gross-Pitaevskii equations, varying inter- and intra-species scattering lengths so as to probe the entire range of miscibility-immiscibility transitions. In the limit of the one-dimensional (1D) ring, i.e., a very narrow one, stability of mixed states is studied analytically, including hidden-vorticity (HV) modes, i.e., those with opposite vorticities of the two components and zero total angular momentum. The consideration of demixed 1D states reveals, in addition to stable composite single-peak structures, double- and triple-peak ones, above a certain particle-number threshold. In the 2D annular geometry, stable demixed states exist both in radial and azimuthal configurations. We find that stable radially-demixed states can carry arbitrary vorticity and, counter-intuitively, the increase of the vorticity enhances stability of such states, while unstable ones evolve into randomly oscillating angular demixed modes. The consideration of HV states in the 2D geometry expands the stability range of radially-demixed states., Comment: to be published in New J. Phys

We introduce a two-component system which models a pseudospinor Bose-Einstein condensate (BEC), with a microwave field coupling its two components. The feedback of BEC of the field (the local-field effect) is taken into account by dint of the respective Poisson equation, which is solved using the Green's function. This gives rise to an effective long-range self-trapping interaction, which may act alone, or be combined with the contact cubic nonlinearity. The system is made discrete by loading the BEC into a deep optical-lattice potential. Numerical solutions demonstrate that onsite-centered fundamental solitons are stable in the cases of attractive or zero contact interactions, while offsite-centered solitons are unstable. In the case of the repulsive onsite nonlinearity, offsite solitons are stable, while their onsite-centered counterparts are stable only at sufficiently small values of the norm, where bistability between the off- and onsite-centered mode takes place. The shape of the onsite-centered solitons is very accurately predicted by a variational approximation (which includes essential technical novelties). Spatially-antisymmetric (\textquotedblleft twisted") solitons are stable at small values of the norm, being unstable at larger norms. In the strongly asymmetric version of the two-component system, which includes the Zeeman splitting, the system is reduced to a single discrete Gross-Pitaevskii equation, by eliminating the small higher-energy component., Comment: 12 Pages, 10 Figures, and 41 References

It was recently found that the spin-orbit (SO) coupling can help to create stable matter-wave solitons in spinor Bose-Einstein condensates in the two-dimensional (2D) free space. Being induced by external laser illumination, the effective SO coupling can be applied too in a spatially confined area. Using numerical methods and the variational approximation (VA), we build families of 2D solitons of the semi-vortex (SV) and mixed-mode (MM) types, and explore their stability, assuming that the SO-coupling strength is confined in the radial direction as a Gaussian. The most essential result is identification, by means of the VA and numerical methods, of the minimum size of the spatial confinement for which the 2D system maintains stable solitons of the SV and MM types., Comment: 9 pages, 6 figures, and 61 references, Commun Nonlinear Sci Numer Simulat 73 481 (2019)

We introduce a one-dimensional PT-symmetric system, which includes the cubic self-focusing, a double-well potential in the form of an infinitely deep potential box split in the middle by a delta-functional barrier of an effective height epsilon, and constant linear gain and loss, gamma, in each half-box. The setting may be realized in microwave photonics. Using numerical methods, we construct PT-symmetric and antisymmetric modes, which represent, respectively, the system's ground state and first excited state, and identify their stability. When they are unstable, the instability leads to blowup, except for the case of epsilon = 0, when an unstable symmetric mode transforms into a weakly oscillating breather, and an unstable antisymmetric mode relaxes into a stable symmetric one. At epsilon > 0, the stability area is much larger for the PT-antisymmetric state than for its symmetric counterpart. The stability areas shrink with with increase of the total power, P. In the linear limit, which corresponds to P --> 0, the stability boundary is found in an analytical form. The stability area of the antisymmetric state originally expands with the growth of gamma, and then disappears at a critical value of gamma., Comment: to be published in Phil. Trans. Roy. Soc. A (a theme issue on "Dissipative Structures in Matter out of Equilibrium: From Chemistry, Photonics and Biology", ed. by M. Tlidi, M. Clerc, and K. Panaiotov)

It was recently found that the Lee-Huang-Yang (LHY) correction to the mean-field Hamiltonian suppresses the collapse and creates stable localized modes (two-component "quantum droplets", QDs) in two and three dimensions. We construct two-dimensional\ self-trapped modes in the form of QDs with vorticity $S$ embedded into each component. The QDs feature a flat-top shape, which expands with the increase of $S$ and norm $N$. An essential finding, produced by a systematic numerical analysis and analytical estimates, is that the vortical QDs are \emph{stable} (which is a critical issue for vortex solitons in nonlinear models) up to $S=5$, for $N$ exceeding a certain threshold value. In the condensate of $^{39}$K atoms, in which QDs with $S=0$ and a quasi-2D shape were created recently, the vortical droplets may have radial size $\lesssim 30$ $\mathrm{\mu}$m, with the number of atoms in the range of $10^{4}-10^{5}$. It is worthy to note that \textit{hidden-vorticity} states in QDs with topological charges $% S_{+}=-S_{-}=1$ in its components, which are prone to strong instability in other settings, have their stability region too, although it may be located beyond applicability limits of the underlying model. Dynamics of elliptically deformed QDs, which form rotating elongated patterns or ones with strong oscillations of the eccentricity, as well as collisions of QDs, are also addressed., Comment: 13 pages, 7 figures, and 75 references, Physical Review A, in press

We study the spatiotemporal solitary modes that propagate in a hollow twisted cylinder waveguide pipe with a self-focusing Kerr nonlinearity. Three generic solitary modes, one belonging to the zero-harmonic (0H) and the other two belonging to the first-harmonic (1H), are found in the first rotational Brillouin zone. The 0H solitary modes can be termed as a quasi-1D (one-dimensional) temporal soliton. Their characteristics depend only on the energy flow. The 1H solitary mode can be termed a quasi-2D (two-dimensional) bullet, whose width is much narrower than the angular domain of the waveguide. In contrast to the 0H mode, the characteristics of the 1H solitary mode depend on both their energy flow and the rotating speed of the waveguide. We demonstrate numerically that the 1H solitary modes are stable when their energy flow is smaller than the threshold norm of the \emph{Townes soliton}. The boundaries of the bistable area for these two types of solitary modes are predicted by the analyses via two-mode approximation. This prediction is in accordance with the numerical findings. We also demonstrate analytically that the 1H solitary mode of this system can be applied to emulate the nonlinear dynamics of solitary modes with 1D Rashba spin-orbit (SO) coupling by optics. Two degenerated states of the 1H solitary mode, semi-dipole and mixed mode, are found from this setting via the mechanism of SO coupling. Collisions between the pair of these two types of solitary modes are also discussed in the paper. The pair of the 0H solitary mode features only the elastic collision, whereas the pair of 1H solitary modes can feature both elastic and inelastic collision when the total energy flow of the two modes are smaller or close to the threshold norm of the \emph{Townes soliton}., Comment: 6 figures, 52 references