Dark Matter could potentially manifest itself in the form of asymmetric dark stars. In this paper we entertain the possibility of probing such asymmetric bosonic dark matter stars by the use of atomic clocks. If the dark sector connects to the standard model sector via a Higgs or photon portal, the interior of boson stars that are in a Bose-Einstein condensate state can change the values of physical constants that control the timing of atomic clock devices. Dilute asymmetric dark matter boson stars passing through the Earth can induce frequency shifts that can be observed in separated Earth based atomic clocks. This gives the opportunity to probe a class of dark matter candidates that for the moment cannot be detected with any different conventional method.

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 non-relativistic limit, the Gross-Pitaevskii equation. A single free parameter, conveniently chosen as the central value of the wavefunction of the axion star, or alternatively the chemical potential with range $-m<\mu< 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.

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.

We analyze the accuracy of the variational method in computing physical quantities relevant for gravitationally bound Bose-Einstein condensates. Using a variety of variational ans\"atze found in existing literature, we determine physical quantities and compare them to exact numerical solutions. We conclude that a "linear+exponential" wavefunction 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 wavefunction fits better for large $N$. The Gaussian wavefunction 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 exact 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\"atze 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.

Coherent oscillation of axions or axion-like particles may give rise to long-lived clumps, called axion stars, because of the attractive gravitational force or its self-interaction. Such a kind of configuration has been extensively studied in the context of oscillons without the effect of gravity, and its stability can be understood by an approximate conservation of particle number in a non-relativistic effective field theory (EFT). We extend this analysis to the case with gravity to discuss the longevity of axion stars and clarify the EFT expansion scheme in terms of gradient energy and Newton's constant. Our EFT is useful to calculate the axion star configuration and its classical lifetime without any ad hoc assumption. In addition, we derive a simple stability condition against small perturbations. Finally, we discuss the consistency of other non-relativistic effective field theories proposed in the literature.

We develop a self-consistent, Gravitoelectromagnetic (GEM) formulation of a slowly rotating, self-gravitating and dilute Bose-Einstein condensate (BEC), intended for astrophysical applications in the context of dark matter halos. GEM self-consistently incorporates the effects of frame dragging to lowest order in $v/c$ via the Gravitomagnetic field. BEC dark matter has attracted attention as an alternative to Cold dark matter (CDM) and Warm dark matter (WDM) for some time now. The BEC is described by the Gross-Pitaevskii-Poisson (GPP) equation with an arbitrary potential allowing for either attractive or repulsive interactions. Owing to the difficulty in obtaining exact solutions to the GEM equations of motion without drastic approximations, we employ the variational method to examine the conditions under which rotating condensates, stable against gravitational collapse, may form in models with attractive and repulsive quartic interactions. We also describe the approximate dynamics of an imploding and rotating condensate by employing a collective coordinate description in terms of the condensate radius.

The method pioneered by Ruffini and Bonazzola (RB) to describe boson stars involves an expansion of the boson field which is linear in creation and annihilation operators. In the nonrelativistic limit, the equation of motion of RB is equivalent to the nonlinear Schr\"odinger equation. Further, the RB expansion constitutes an exact solution to a non-interacting field theory, and has been used as a reasonable ansatz for an interacting one. In this work, we show how one can go beyond the RB ansatz towards an exact solution of the interacting operator Klein-Gordon equation, which can be solved iteratively to ever higher precision. Our Generalized Ruffini-Bonazzola approach takes into account contributions from nontrivial harmonic dependence of the wavefunction, using a sum of terms with energy $k\,E_0$, where $k\geq1$ and $E_0$ is the chemical potential of a single bound axion. The method critically depends on an expansion in a parameter $\Delta \equiv \sqrt{1 - E_0{}^2/m^2} < 1$, where $m$ is the mass of the boson. In the case of the axion potential, we calculate corrections which are relevant for axion stars in the transition or dense branches of solutions. We find with high precision the local minimum of the mass, $M_{min}\approx 463\,f^2/m$, at $\Delta\approx0.27$, where $f$ is the axion decay constant. This point marks the crossover from the transition branch to the dense branch of solutions, and a corresponding crossover from structural instability to stability.

Axion particles can form macroscopic condensates, whose size can be galactic in scale for models with very small axion masses $m\sim10^{-22}$ eV, and which are sometimes referred to under the name of Fuzzy Dark Matter. Many analyses of these condensates are done in the non-interacting limit, due to the weakness of the self-interaction coupling of axions. We investigate here how certain results change upon inclusion of these interactions, finding a decreased maximum mass and a modified mass-radius relationship. Further, these condensates are, in general, unstable to decay through number-changing interactions. We analyze the stability of galaxy-sized condensates of axion-like particles, and sketch the parameter space of stable configurations as a function of a binding energy parameter. We find a strong lower bound on the size of Fuzzy Dark Matter condensates which are stable to decay, with lifetimes longer than the age of the universe.

We summarize our work on constant roll inflationary models. It was understood recently that constant roll inflation, in a regime beyond the slow roll approximation, can give models that are in agreement with the observational constraints. We describe a new class of constant roll inflationary models and investigate the behavior of scalar perturbations in them. We also comment on other non-slow roll regimes of inflation.

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 \sim 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 ($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.

We investigate the decay of condensates of scalars in a field theory defined by $V({\cal A})=m{}^2\,f{}^2\,[1-\cos({\cal A}\,/\,f)]$, where $m$ and $f$ are the mass and decay constant of the scalar field. An example of such a theory is that of the axion, in which case the condensates are called axion stars. The axion field, $\cal A$, is self adjoint. As a result the axion number is not an absolutely conserved quantity. Therefore, axion stars are not stable and have finite lifetimes. Bound axions, localized on the volume of the star, have a coordinate uncertainty $\Delta x \sim R \sim 1/(m_a \Delta)$, where $R$ is the radius of the star and $\Delta = \sqrt{1-E_0{}^2\,/\,m_a^2}$. Here $m_a$ and $E_0$ are the mass, and the ground state energy of the bound axion. Then the momentum distribution of axions has a width of $\Delta p\sim m_a\,\Delta$. At strong binding, $\Delta={\cal O}(1)$, bound axions can easily transfer a sufficient amount of momentum to create and emit a free axion, leading to fast decay of the star with a transition rate $\Gamma \sim m_a$. However, when $\Delta\ll 1$, the momentum distribution is more restricted, and as shown in this paper, the transition rate for creating a free axion decreases as $\exp (-p \, \Delta x) \sim \exp (-\Delta^{-1})$. Then sufficiently large, weakly bound axion stars, produced after the big bang, survive until the present time. We plot the region of their stability, limited by decay through axion loss and by gravitational instability, as a function of the mass of the axion and the mass of the star.

We study the possibility that self-interacting bosonic dark matter forms star-like objects. We study both the case of attractive and repulsive self-interactions, and we focus particularly in the parameter phase space where self-interactions can solve well standing problems of the collisionless dark matter paradigm. We find the mass radius relations for these dark matter bosonic stars, their density profile as well as the maximum mass they can support.

We investigate some of the experimental, observational and theoretical consequences of hypothetical stable black holes in the mass range between the electro-weak scale and the Planck mass, 2.4$\times 10^{15}$ TeV. For the purpose of calculations we use Lovelock black holes in odd dimensions. If such black holes exist they contribute to dark matter. We show that the passage of the black holes through matter and the collision of black holes have a well defined experimental signature. Depending on their cross section and energy they also accumulate in stars and influence their development.

Black branes are studied in Einstein-Gauss-Bonnet (EGB) gravity. Evaporation drives black branes towards one of two singularities depending on the sign of $\alpha$, the Gauss-Bonnet coupling. For positive $\alpha$ and sufficiently large ratio $\sqrt{\alpha}/L$, where $L/2\pi$ is the radius of compactification, black branes avoid the Gregory-Laflamme (GL) instability before reaching a critical state. No black branes with the radius of horizon smaller than the critical value can exist. Approaching the critical state branes have a nonzero Hawking temperature. For negative $\alpha$ all black branes encounter the GL instability. No black branes may exist outside of the interval of the critical values, $0\leq\beta<3$, where $\beta=1- 8\alpha/r_h^2$ and $r_h$ is the radius of horizon of the black brane. The first order phase transition line of GL transitions ends in a second order phase transition point at $\beta=0$.

In a recent publication we developed a canonical quantization program describing the gravitational collapse of a spherical dust cloud in 2+1 dimensions with a negative cosmological constant $-\Lambda\equiv -l^{-2}<0$. In this paper we address the quantization of the Banados--Teitelboim--Zanelli (BTZ) black hole. We show that the mass function describing the black hole is made of two pieces, a constant non-vanishing boundary contribution and a discrete spectrum of the form $\mu_n = \frac{\hbar}{l}(n+ \frac 12)$. The discrete spectrum is obtained by applying the Wheeler--DeWitt equation with a particular choice of factor ordering and interpreted as giving the energy levels of the collapsed matter shells that form the black hole. Treating a black hole microstate as a particular distribution of shells among the levels, we determine the canonical entropy of the BTZ black hole. Comparison with the Bekenstein--Hawking entropy shows that the boundary energy is related to the central charge of the Virasoro algebra that generates the asymptotic symmetry group of the three-dimensional anti-de Sitter space AdS$_3$. This gives a connection between the Wheeler--DeWitt approach and the conformal field theory approach.

We investigate the thermodynamics of static black objects such as black holes, black strings and their generalizations to D dimensions (`black branes') in a gravitational theory containing the four dimensional Gauss-Bonnet term in the action, when D-4 of the dimensions are compactified on a torus. The entropies of black holes and black branes are compared to obtain information on the stability of these objects and to find their phase diagrams. We demonstrate the existence of a critical mass, which depends on the scale of the compactified dimensions, below which the black hole entropy dominates over the entropy of the black membrane.