We develop an irregular lattice mass-spring model to simulate and study the deformation modes of a thin elastic ribbon as a function of applied end-to-end twist and tension. Our simulations reproduce all reported experimentally observed modes, including transitions from helicoids to longitudinal wrinkles, creased helicoids and loops with self-contact, and transverse wrinkles to accordion self-folds. Our simulations also show that the twist angles at which the primary longitudinal and transverse wrinkles appear are well described by various analyses of the Föppl–von Kármán equations, but the characteristic wavelength of the longitudinal wrinkles has a more complex relationship to applied tension than previously estimated. The clamped edges are shown to suppress longitudinal wrinkling over a distance set by the applied tension and the ribbon width, but otherwise have no apparent effect on measured wavelength. Further, by analyzing the stress profile, we find that longitudinal wrinkling does not completely alleviate compression, but caps the magnitude of the compression. Nonetheless, the width over which wrinkles form is observed to be wider than the near-threshold analysis predictions: the width is more consistent with the predictions of far-from-threshold analysis. However, the end-to-end contraction of the ribbon as a function of twist is found to more closely follow the corresponding near-threshold prediction as tension in the ribbon is increased, in contrast to the expectations of far-from-threshold analysis. These results point to the need for further theoretical analysis of this rich thin elastic system, guided by our physically robust and intuitive simulation model.

Ultralight scalar dark matter has been proposed to constitute a component of dark matter, though the minimal scenarios have increasingly become constrained. In this work, we analyze scenarios where the dark matter consists of more than one ultralight boson, each with different masses. This potentially leads to formation of gravitationally-bound Bose-Einstein condensates with structures that are very different from condensates composed of a single scalar field. By generalizing from the well-understood single-flavor case, we explore a large range of input parameters, subject to stability criteria, and determine the allowed parameter space for two-flavor condensates as a function of particle physics parameters, paying particular attention to cases where such condensates could compose galactic cores. We also analyze single-flavor condensates subject to external gravity from massive inner bodies and find that such systems may mimic the size of galactic cores as well.

Taking a comprehensive view, including a full range of boundary conditions, we reexamine QCD axion star solutions based on the relativistic Klein-Gordon equation (using the Ruffini-Bonazzola approach) and its nonrelativistic limit, the Gross-Pitaevskii equation. A single free parameter, conveniently chosen as the central value of the wave function of the axion star, or alternatively the chemical potential with range $-m<μ<0$ (where $m$ is the axion mass), uniquely determines a spherically symmetric ground state solution, the axion condensate. We clarify how the interplay of various terms of the Klein-Gordon equation determines the properties of solutions in three separate regions: the structurally stable (corresponding to a local energy minimum) dilute and dense regions, and the intermediate, structurally unstable transition region. From the Klein-Gordon equation, one can derive alternative equations of motion including the Gross-Pitaevskii and Sine-Gordon equations, which have been used previously to describe axion stars in the dense region. In this work, we clarify precisely how and why such methods break down as the binding energy increases, emphasizing the necessity of using the full relativistic Klein-Gordon approach. Finally, we point out that, even after including perturbative axion number violating corrections, solutions to the equations of motion, which assume approximate conservation of axion number, break down completely in the regime with strong binding energy, where the magnitude of the chemical potential approaches the axion mass.

We analyze the accuracy of the variational method in computing physical quantities relevant for gravitationally bound Bose-Einstein condensates. Using a variety of spherically symmetric variational ansätze found in existing literature, we determine physical quantities and compare them to numerical solutions. We conclude that a ``linear + exponential'' wave function proportional to $(1+\xi) exp(-\xi)$ (where $\xi$ is a dimensionless radial variable) is the best fit for attractive self-interactions along the stable branch of solutions, while for small particle number $N$ it is also the best fit for repulsive self-interactions. For attractive self-interactions along the unstable branch, a single exponential is the best fit for small $N$, while a sech wave function fits better for large $N$. The Gaussian wave function ansatz, which is used often in the literature, is exceedingly poor across most of the parameter space, with the exception of repulsive interactions for large $N$. We investigate a “double exponential” ansatz with a free constant parameter, which is computationally efficient and can be optimized to fit the numerical solutions in different limits. We show that the double exponential can be tuned to fit the sech ansatz, which is computationally slow. We also show how to generalize the addition of free parameters in order to create more computationally efficient ansätze using the double exponential. Determining the best ansatz, according to several comparison parameters, will be important for analytic descriptions of dynamical systems. Finally, we examine the underlying relativistic theory, and critically analyze the Thomas-Fermi approximation often used in the literature.

Stability properties of gravitationally bound condensates composed of ultralight axionic Fuzzy Dark Matter (FDM) are studied. Previous work has shown that astrophysical collisions could make self-gravitating condensates structurally unstable, making them prone to collapse and decay; in the context of FDM, we reexamine the relevant timescales using the time-dependent variational method. We show that FDM condensates can be made unstable through gravitational interactions with central black holes, for black hole masses in a phenomenologically relevant range. Instability could also be stimulated by galaxy collisions. The subsequent decay takes place over a period lasting as long as many thousands of years. We also discuss the possible relevance of FDM condensates to understanding the composition of Ultracompact Dwarf (UCD) Galaxies. Future observation of extremely massive black holes in the central regions of UCDs can constrain this interpretation.

In a previous work, we analyzed collapsing axion stars using the low-energy instanton potential, showing that the total energy is always bounded and that collapsing axion stars do not form black holes. In this paper, we provide a proof that the conclusions are unchanged when using instead the more general chiral potential for QCD axions.

If QCD axions form a large fraction of the total mass of dark matter, then axion stars could be very abundant in galaxies. As a result, collisions with each other, and with other astrophysical bodies, can occur. We calculate the rate and analyze the consequences of three classes of collisions, those occurring between a dilute axion star and: another dilute axion star, an ordinary star, or a neutron star. In all cases we attempt to quantify the most important astrophysical uncertainties; we also pay particular attention to scenarios in which collisions lead to collapse of otherwise stable axion stars, and possible subsequent decay through number changing interactions. Collisions between two axion stars can occur with a high total rate, but the low relative velocity required for collapse to occur leads to a very low total rate of collapses. On the other hand, collisions between an axion star and an ordinary star have a large rate, \(\Gamma_{\odot} \) ∼ 3000 collisions/year/galaxy, and for sufficiently heavy axion stars, it is plausible that most or all such collisions lead to collapse. We identify in this case a parameter space which has a stable region and a region in which collision triggers collapse, which depend on the axion number (N ) in the axion star, and a ratio of mass to radius cubed characterizing the ordinary star (\(\frac{M_S}{R_S^3}\)). Finally, we revisit the calculation of collision rates between axion stars and neutron stars, improving on previous estimates by taking cylindrical symmetry of the neutron star distribution into account. Collapse and subsequent decay through collision processes, if occurring with a significant rate, can affect dark matter phenomenology and the axion star mass distribution.

Axion stars, gravitationally bound states of low-energy axion particles, have a maximum mass allowed by gravitational stability. Weakly bound states obtaining this maximum mass have sufficiently large radii such that they are dilute, and as a result, they are well described by a leading-order expansion of the axion potential. Heavier states are susceptible to gravitational collapse. Inclusion of higher-order interactions, present in the full potential, can give qualitatively different results in the analysis of collapsing heavy states, as compared to the leading-order expansion. In this work, we find that collapsing axion stars are stabilized by repulsive interactions present in the full potential, providing evidence that such objects do not form black holes. In the last moments of collapse, the binding energy of the axion star grows rapidly, and we provide evidence that a large amount of its energy is lost through rapid emission of relativistic axions.