Date: Thu, 9 May 2024 | Total: 4
We compare two different numerical methods to integrate in time spatially delocalized initial densities using the Schr\"odinger-Poisson equation system as the evolution law. The basic equation is a nonlinear Schr\"odinger equation with an auto-gravitating potential created by the wave function density itself. The latter is determined as a solution of Poisson's equation modelling, e.g., non-relativistic gravity. For reasons of complexity, we treat a one-dimensional version of the problem whose numerical integration is still challenging because of the extreme long-range forces (being constant in the asymptotic limit). Both of our methods, a Strang splitting scheme and a basis function approach using B-splines, are compared in numerical convergence and effectivity. Overall, our Strang-splitting evolution compares favourably with the B-spline method. In particular, by using an adaptive time-stepper rather large one-dimensional boxes can be treated. These results give hope for extensions to two spatial dimensions for not too small boxes and large evolution times necessary for describing, for instance, dark matter formation over cosmologically relevant scales.
We study the late-time dynamics of two particles confined in one spatial dimension and subject to two-body losses. The dynamics is exactly described by a non-Hermitian Hamiltonian that can be analytically studied both in the continuum and on a lattice. The asymptotic decay rate and the universal power-law form of the decay of the number of particles are exactly computed in the whole parameter space of the problem. When in the initial state the two particles are far apart, the average number of particles in the setup decays with time $t$ as $t^{-1/2}$; a different power law, $t^{-3/2}$, is found when the two particles overlap in the initial state. These results are valid both in the continuum and on a lattice, but in the latter case a logarithmic correction appears.
We consider a mixture of bosons and spin-polarized fermions in two dimensions at zero temperature with a tunable Bose-Fermi attraction. By adopting a diagrammatic T-matrix approach, we analyze the behavior of several thermodynamic quantities for the two species as a function of the density ratio and coupling strength, including the chemical potentials, the momentum distribution functions, the boson condensate density, and the Tan's contact parameter. By increasing the Bose-Fermi attraction, we find that the condensate is progressively depleted and Bose-Fermi pairs form, with a small fraction of condensed bosons surviving even for strong Bose-Fermi attraction. This small condensate proves sufficient to hybridize molecular and atomic states, producing quasi-particles with unusual Fermi liquid features. A nearly universal behavior of the condensate fraction, the bosonic momentum distribution, and Tan's contact parameter with respect to the density ratio is also found.
We analyze the pinning of vortices for a stationary rotating dipolar supersolid along the low-density paths between droplets as a function of the rotation frequency. We restrict ourselves to the stationary configurations of vortices with the same symmetry as that of the array of droplets. Our approach exploits the fact that the wave function of each droplet acquires a linear phase on the coordinates, and hence the relative phases between neighboring droplets allows us to predict the position of the vortices. For a confined system, the estimate accurately reproduces the Gross-Pitaevskii results in the spatial regions where the neighboring droplets are well defined.