Shift-symmetric Horndeski theories admit an interesting class of Schwarzschild black hole solutions exhibiting time-dependent scalar hair. By making use of Lema\^{i}tre coordinates, we analyze perturbations around these types of black holes, and demonstrate that scalar perturbations around black hole backgrounds inevitably have gradient instabilities. Taken together with previously established results, this newly-discovered instability rules out black holes with time-dependent scalar hair in Horndeski theories.

Extremely light bosonic wave dark matter ($\psi$DM) is an emerging dark matter candidate contesting the conventional cold dark matter paradigm and a model subject to intense scrutiny of late. This work for the first time reports testable salient features pertinent to gravitational lenses of $\psi$DM halos. $\psi$DM halos are distinctly filled with large-amplitude, small-scale density fluctuations with $\delta\rho/\rho_{\rm halo}\sim 1$ in form of density granules. This halo characteristically yields fuzzy gravitational lenses that can yield ubiquitous flux ratio anomalies of few tens percent levels often found in lensed quasar multiple images. These fuzzy lenses may also produce rare hexad images and octad images, formed when the sources are located in some well-defined zones of highly distorted caustics. We have found new critical features appearing in the highly de-magnified lens center when the halo has sufficiently high surface density near a very compact massive core.

Quasinormal modes describe the return to equilibrium of a perturbed system, in particular the ringdown phase of a black hole merger. But as globally-defined quantities, the quasinormal spectrum can be highly sensitive to global structure, including distant small perturbations to the potential. In what sense are quasinormal modes a property of the resulting black hole? We explore this question for the linearized perturbation equation with two potentials having disjoint bounded support. We give a composition law for the Wronskian that determines the quasinormal frequencies of the combined system. We show that over short time scales the evolution is governed by the quasinormal frequencies of the individual potentials, while the sensitivity to global structure can be understood in terms of echoes. We introduce an echo expansion of the Green's function and show that, as expected on general grounds, at any finite time causality limits the number of echoes that can contribute. We illustrate our results with the soluble example of a pair of $\delta$-function potentials. We explicate the causal structure of the Green's function, demonstrating under what conditions two very different quasinormal spectra give rise to very similar ringdown waveforms.

We show that a black hole surrounded by scalar dark matter develops scalar hair. This is the generalization of a phenomenon pointed out by Jacobson, that a minimally coupled scalar with a non-trivial time dependence far away from the black hole would endow the black hole with hair. In our case, the time dependence arises from the oscillation of a scalar field with a non-zero mass. We systematically explore the scalar profile around the black hole for different scalar masses. In the small mass limit, the scalar field has a $1/r$ component at large radius $r$, consistent with Jacobson's result. In the large mass limit (with the Compton wavelength of order of the horizon or smaller), the scalar field has a $1/r^{3/4}$ profile yielding a pile-up close to the horizon, while distinctive nodes occur for intermediate masses. Thus, the dark matter profile around a black hole, while challenging to measure, contains information about the dark matter particle mass. As an application, we consider the case of the supermassive black hole at the center of M87, recently imaged by the Event Horizon Telescope. Its horizon size is roughly the Compton wavelength of a scalar particle of mass $10^{-20}$ eV. We consider the implications of the expected scalar pile-up close to the horizon, for fuzzy dark matter at a mass of $10^{-20}$ eV or below.

We reconsider the derivation of soft theorems associated with nonlinearly-realized symmetries in cosmology. Utilizing the path integral, we derive a generalized consistency relation that relates a squeezed $(N+1)$-point correlation function to an $N$-point function, where the relevant soft mode is at early rather than late time. This generalized (early-late-time) version has wider applicability than the standard consistency relation where all modes are evaluated at late times. We elucidate the conditions under which the latter follows from the former. A key ingredient is the physical mode condition: that the nonlinear part of the symmetry transformation must match the time dependence of the dominant, long wavelength physical mode. This is closely related to, but distinct from, the adiabatic mode condition. Our derivation sheds light on a number of otherwise puzzling features of the standard consistency relation: (1) the underlying nonlinearly-realized symmetries (such as dilation and special conformal transformation SCT) originate as residual gauge redundancies, yet the consistency relation has physical content---for instance, it can be violated; (2) the standard consistency relation is known to fail in ultra-slow-roll inflation, but since dilation and SCT remain good symmetries, there should be a replacement for the standard relation; (3) in large scale structure applications, it is known that the standard consistency relation breaks down if the long wavelength power spectrum is too blue. The early-late-time consistency relation helps address these puzzles. We introduce a toy model where explicit checks of this generalized consistency relation are simple to carry out. Our methodology can be adapted to cases where violations of the standard consistency relation involve additional light degrees of freedom beyond the inflaton.

Wave dark matter ($\psi$DM), which satisfies the Schr\"odinger-Poisson equation, has recently attracted substantial attention as a possible dark matter candidate. Numerical simulations have in the past provided a powerful tool to explore this new territory of possibility. Despite their successes to reveal several key features of $\psi$DM, further progress in simulations is limited, in that cosmological simulations so far can only address formation of halos below $\sim 2\times 10^{11} M_\odot$ and substantially more massive halos have become computationally very challenging to obtain. For this reason, the present work adopts a different approach in assessing massive halos by constructing wave-halo solutions directly from the wave distribution function. This approach bears certain similarity with the analytical construction of particle-halo (cold dark matter model). Instead of many collisionless particles, one deals with one single wave that has many non-interacting eigenstates. The key ingredient in the wave-halo construction is the distribution function of the wave power, and we use several halos produced by structure formation simulations as templates to determine the wave distribution function. Among different models, we find the fermionic King model presents the best fits and we use it for our wave-halo construction. We have devised an iteration method for constructing the nonlinear halo, and demonstrate its stability by three-dimensional simulations. A Milky-Way-sized halo has also been constructed, and the inner halo is found flatter than the NFW profile. These wave-halos have small-scale interferences both in space and time producing time-dependent granules. While the spatial scale of granules varies little, the correlation time is found to increase with radius by one order of magnitude across the halo.

We use holographic techniques to compute inflationary non-Gaussianities for general single-field inflation, including models with a non-trivial sound speed. In this holographic approach, the inflationary dynamics is captured by a relevant deformation of the dual conformal field theory (CFT) in the UV, while the inflationary correlators are computed by conformal perturbation theory. In this paper, we discuss the effects of higher derivative operators, such as $(\partial_\mu\phi\partial^\mu\phi)^{m}$, which are known to induce a non-trivial sound speed and source potentially large non-Gaussianities. We compute the full inflationary bispectra from the deformed CFT correlators. We also discuss the squeezed limit of the bispectra from the viewpoint of operator product expansions. As is generic in the holographic description of inflation, our power spectrum is blue tilted in the UV region. We extend our bispectrum computation to the IR region by resumming the conformal perturbations to all orders. We provide a self-consistent setup which reproduces a red tilted power spectrum, as well as all possible bispectrum shapes in the slow-roll regime.

We examine the nonlinear structure of gravitationally collapsed objects that form in our simulations of wavelike cold dark matter ($\psi$DM), described by the Schr\"{o}dinger-Poisson (SP) equation with a particle mass $\sim 10^{-22} {\rm eV}$. A distinct gravitationally self-bound solitonic core is found at the center of every halo, with a profile quite different from cores modeled in the warm or self-interacting dark matter scenarios. Furthermore, we show that each solitonic core is surrounded by an extended halo composed of large fluctuating dark matter granules which modulate the halo density on a scale comparable to the diameter of the solitonic core. The scaling symmetry of the SP equation and the uncertainty principle tightly relate the core mass to the halo specific energy, which, in the context of cosmological structure formation, leads to a simple scaling between core mass ($M_c$) and halo mass ($M_h$), $M_c \propto a^{-1/2} M_h^{1/3}$, where $a$ is the cosmic scale factor. We verify this scaling relation by (i) examining the internal structure of a statistical sample of virialized halos that form in our 3D cosmological simulations, and by (ii) merging multiple solitons to create individual virialized objects. Sufficient simulation resolution is achieved by adaptive mesh refinement and graphic processing units acceleration. From this scaling relation, present dwarf satellite galaxies are predicted to have kpc sized cores and a minimum mass of $\sim 10^8 {M_\odot}$, capable of solving the small-scale controversies in the cold dark matter model. Moreover, galaxies of $2\times10^{12} {M_\odot}$ at $z=8$ should have massive solitonic cores of $\sim 2\times10^9 {M_\odot}$ within $\sim 60 {\rm pc}$. Such cores can provide a favorable local environment for funneling the gas that leads to the prompt formation of early stellar spheroids and quasars.

Recent BICEP2 detection of low-multipole B-mode polarization anisotropy in the cosmic microwave background radiation supports the inflationary universe scenario and suggests a large inflaton field range. The latter feature can be achieved with axion fields in the framework of string theory. We present such a helical model which naturally becomes a model with a single cosine potential, and which in turn reduces to the (quadratic) chaotic inflation model in the super-Planckian limit. The slightly smaller tensor/scalar ratio $r$ of models of this type provides a signature of the periodic nature of an axion potential. We present a simple way to quantify this distinctive feature. As axions are intimately related to strings/vortices and strings are ubiquitous in string theory, we explore the possibility that cosmic strings may be contributing to the B-mode polarization anisotropy observed.

We study a racetrack model in the presence of the leading alpha'-correction in flux compactification in Type IIB string theory, for the purpose of getting conceivable de-Sitter vacua in the large compactified volume approximation. Unlike the K\"ahler Uplift model studied previously, the alpha'-correction is more controllable for the meta-stable de-Sitter vacua in the racetrack case since the constraint on the compactified volume size is very much relaxed. We find that the vacuum energy density \Lambda for de-Sitter vacua approaches zero exponentially as the volume grows. We also analyze properties of the probability distribution of \Lambda in this class of models. As in other cases studied earlier, the probability distribution again peaks sharply at \Lambda=0. We also study the Racetrack K\"ahler Uplift model in the Swiss-Cheese type model.

In this paper, we describe the architecture and performance of the GraCCA system, a Graphic-Card Cluster for Astrophysics simulations. It consists of 16 nodes, with each node equipped with 2 modern graphic cards, the NVIDIA GeForce 8800 GTX. This computing cluster provides a theoretical performance of 16.2 TFLOPS. To demonstrate its performance in astrophysics computation, we have implemented a parallel direct N-body simulation program with shared time-step algorithm in this system. Our system achieves a measured performance of 7.1 TFLOPS and a parallel efficiency of 90% for simulating a globular cluster of 1024K particles. In comparing with the GRAPE-6A cluster at RIT (Rochester Institute of Technology), the GraCCA system achieves a more than twice higher measured speed and an even higher performance-per-dollar ratio. Moreover, our system can handle up to 320M particles and can serve as a general-purpose computing cluster for a wide range of astrophysics problems.

We investigate Kolmogorov wave turbulence in QCD or, in other words, we calculate the spectrum of gluons as a function of time, f_k(t), in the presence of a source which feeds in energy density in the infrared region at a constant rate. We find an early, an intermediate and a late time form for the gluon spectrum. Wave turbulence in QCD turns out to be somewhat different than the turbulence in the case of \phi^4-type theories studied by Zakharov, L'vov and Falkovich. The hope is that a good understanding of QCD wave turbulence might lead to a better understanding of the instability problem in the early stages of the evolution after a heavy ion collision.

We present a scenario for non-radiative accretion onto the supermassive black hole at the galactic center. Conducting MHD simulations with 1400^3 grid zones that break the axial and reflection symmetries of earlier investigations and extend inward from the Bondi radius, we find a quasi-hydrostatic radial density profile rho propto r^{-0.72} with superadiabatic gradient corresponding to an n~0.72 polytrope. Buoyancy generated by magnetic dissipation is resisted by the same fields so effectively that energy is advected inward: a state of magnetically-frustrated convection. This scenario is consistent with observational constraints on energetics andouter boundary conditions.

We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin. The magnetic field is updated separately in two-dimensional advection-constraint steps. The electromotive force (EMF) is computed in the advection step using the TVD scheme, and this same EMF is used immediately in the constraint step in order to preserve \grad.B=0 without the need to store intermediate fluxes. Operator splitting is used to extend the code to three dimensions, and Runge-Kutta is used to get second order accuracy in time. The advantages of this code are high resolution per grid cell, second order accuracy in space and time, enforcement of the \grad.B=0 constraint to machine precision, no memory overhead, speed, and simplicity. A 3-D Fortran implementation less than 400 lines long is made freely available. We also implemented a fully scalable message-passing parallel MPI version. We present tests of the code on MHD waves and shocks.