Integrator for general spin-s Gross-Pitaevskii systems

We provide an algorithm, i-SPin 2, for evolving general spin-s Gross-Pitaevskii or nonlinear Schrödinger systems carrying a variety of interactions, where the 2s+1 components of the "spinor" field represent the different spin-multiplicity states. We consider many nonrelativistic interactions up to quartic order in the Schrödinger field (both short and long range, and spin-dependent and spin-independent interactions), including explicit spin-orbit couplings. The algorithm allows for spatially varying external and/or self-generated vector potentials that couple to the spin density of the field. Our work can be used for scenarios ranging from laboratory systems such as spinor Bose-Einstein condensates (BECs), to cosmological or astrophysical systems such as self-interacting bosonic dark matter. As examples, we provide results for two different setups of spin-1 BECs that employ a varying magnetic field and spin-orbit coupling, respectively, and also collisions of spin-1 solitons in dark matter. Our symplectic algorithm is second-order accurate in time, and is extensible to the known higher-order-accurate methods.

Dark Matters on the Scale of Galaxies

The cold dark-matter model successfully explains both the emergence and evolution of cosmic structures on large scales and, when we include a cosmological constant, the properties of the homogeneous and isotropic Universe. However, the cold dark-matter model faces persistent challenges on the scales of galaxies. Indeed, N-body simulations predict some galaxy properties that are at odds with the observations. These discrepancies are primarily related to the dark-matter distribution in the innermost regions of the halos of galaxies and to the dynamical properties of dwarf galaxies. They may have three different origins: (1) the baryonic physics affecting galaxy formation is still poorly understood and it is thus not properly included in the model; (2) the actual properties of dark matter differs from those of the conventional cold dark matter; (3) the theory of gravity departs from General Relativity. Solving these discrepancies is a rapidly evolving research field. We illustrate some of the solutions proposed within the cold dark-matter model, and solutions when including warm dark matter, self-interacting dark matter, axion-like particles, or fuzzy dark matter. We also illustrate some modifications of the theory of gravity: Modified Newtonian Dynamics (MOND), MOdified Gravity (MOG), and f(R) gravity.

Axion Stars

The dark matter particle can be a QCD axion or axion-like particle. A locally over-densed distribution of axions can condense into a bound Bose–Einstein condensate called an axion star, which can be bound by self-gravity or bound by self-interactions. It is possible that a significant fraction of the dark matter axion is in the form of axion stars. This would make some efforts searching for the axion as the dark matter particle more challenging, but at the same time it would also open up new possibilities. Some of the properties of axion stars, including their emission rates and their interactions with other astrophysical objects, are not yet completely understood.

Galaxy formation with BECDM - I. Turbulence and relaxation of idealized haloes

We present a theoretical analysis of some unexplored aspects of relaxed Bose-Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schrödinger-Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic r-3 NFW-like profile. We find a fundamental relation of the core-to-halo mass with the dimensionless invariant Ξ ≡ |E|/M3/(Gm/ħ)2 or Mc/M ≃ 2.6Ξ1/3, linking the soliton to global halo properties. For r ≥ 3.5 rc core radii, we find equipartition between potential, classical kinetic and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behaviour driven by the continuous reconnection of vortex lines due to wave interference. We analyse the turbulence 1D velocity power spectrum and find a k-1.1 power law. This suggests that the vorticity in BECDM haloes is homogeneous, similar to thermally-driven counterflow BEC systems from condensed matter physics, in contrast to a k-5/3 Kolmogorov power law seen in mechanically-driven quantum systems. The mode where the power spectrum peaks is approximately the soliton width, implying that the soliton-sized granules carry most of the turbulent energy in BECDM haloes.

Relativistic Axions from Collapsing Bose Stars

The substructures of light bosonic (axionlike) dark matter may condense into compact Bose stars. We study the collapse of critical-mass stars caused by attractive self-interaction of the axionlike particles and find that these processes proceed in an unexpected universal way. First, nonlinear self-similar evolution (called "wave collapse" in condensed matter physics) forces the particles to fall into the star center. Second, interactions in the dense center create an outgoing stream of mildly relativistic particles which carries away an essential part of the star mass. The collapse stops when the star remnant is no longer able to support the self-similar infall feeding the collisions. We shortly discuss possible astrophysical and cosmological implications of these phenomena.