Search

Wave (fuzzy) dark matter ($\psi$DM) consists of ultralight bosons, featuring a solitonic core within a granular halo. Here we extend $\psi$DM to two components, with distinct particle masses $m$ and coupled only through gravity, and investigate the resulting soliton-halo structure via cosmological simulations. Specifically, we assume $\psi$DM contains $75$ per cent major component and $25$ per cent minor component, fix the major-component particle mass to $m_{\rm major}=1\times10^{-22}\,{\rm eV}$, and explore two different minor-component particle masses with $m_{\rm major}:m_{\rm minor}=3:1$ and $1:3$, respectively. For $m_{\rm major}:m_{\rm minor}=3:1$, we find that (i) the major- and minor-component solitons coexist, have comparable masses, and are roughly concentric. (ii) The soliton peak density is significantly lower than the single-component counterpart, leading to a smoother soliton-to-halo transition and rotation curve. (iii) The combined soliton mass of both components follows the same single-component core-halo mass relation. In dramatic contrast, for $m_{\rm major}:m_{\rm minor}=1:3$, a minor-component soliton cannot form with the presence of a stable major-component soliton; the total density profile, for both halo and soliton, is thus dominated by the major component and closely follows the single-component case. To support this finding, we propose a toy model illustrating that it is difficult to form a soliton in a hot environment associated with a deep gravitational potential. The work demonstrates the extra flexibility added to the multi-component $\psi$DM model can resolve observational tensions over the single-component model while retaining its key features., Comment: 19 pages, 24 figures, 1 table, accepted for publication in MNRAS

Low mass galaxies in the Local Group are dominated by dark matter and comprise the well studied “dwarf Spheroidal" (dSph) class, with typical masses of 109−10M⊙ and also the equally numerous “ultra faint dwarfs" (UFD), discovered recently, that are distinctly smaller and denser with masses of only 107−8M⊙. This bimodality amongst low mass galaxies contrasts with the scale free continuity expected for galaxies formed under gravity, as in the standard Cold Dark Matter (CDM) model for heavy particles. Within each dwarf class we find the core radius Rc is inversely related to velocity dispersion σ, quite the opposite of standard expectations, but indicative of dark matter in a Bose-Einstein state, where the Uncertainty Principle requires Rc × σ is fixed by Planks constant, h. The corresponding boson mass, mb = h/Rcσ, differs by one order of magnitude between the UDF and dSph classes, with 10−21.4eV and 10−20.3eV respectively. Two boson species is reinforced by parallel relations seen between the central density and radius of UDF and dSph dwarfs respectively, which both match the steep prediction, ρc ∝ Rc−4, for soliton cores in the ground state. Furthermore, soliton cores accurately fit the stellar profiles of UDF and dSph dwarfs where prominent, dense cores appear surrounded by low density halos, as predicted by our simulations. Multiple bosons may point to a String Theory interpretation for dark matter, where a discrete mass spectrum of axions is generically predicted to span many decades in mass, offering a unifying "Axiverse" interpretation for the physical "diversity" of dark matter dominated dwarf galaxies.

We present a new special relativistic hydrodynamics (SRHD) code capable of handling coexisting ultra-relativistically hot and non-relativistically cold gases. We achieve this by designing a new algorithm for conversion between primitive and conserved variables in the SRHD solver, which incorporates a realistic ideal-gas equation of state covering both the relativistic and non-relativistic regimes. The code can handle problems involving a Lorentz factor as high as 106 and optimally avoid the catastrophic cancellation. In addition, we have integrated this new SRHD solver into the code gamer (https://github.com/gamer-project/gamer) to support adaptive mesh refinement and hybrid OpenMP/MPI/GPU parallelization. It achieves a peak performance of 7 × 107 cell updates per second on a single Tesla P100 GPU and scales well to 2048 GPUs. We apply this code to two interesting astrophysical applications: (a) an asymmetric explosion source on the relativistic blast wave and (b) the flow acceleration and limb brightening of relativistic jets.

Fuzzy Dark Matter (FDM), consisting of ultralight bosons ($m_{\rm b} \sim 10^{-22}\ \rm eV$), is an intriguing alternative to Cold Dark Matter. Numerical simulations that solve the Schr\"odinger-Poisson (SP) equation show that FDM halos consist of a central solitonic core, which is the ground state of the SP equation, surrounded by an envelope of interfering excited states. These excited states also interfere with the soliton, causing it to oscillate and execute a confined random walk with respect to the halo center of mass. Using high-resolution numerical simulations of a $6.6 \times 10^9 M_{\odot}$ FDM halo with $m_{\rm b} = 8 \times 10^{-23}\ \rm eV$ in isolation, we demonstrate that the wobbling, oscillating soliton gravitationally perturbs nuclear objects, such as supermassive black holes or dense star clusters, causing them to diffuse outwards. In particular, we show that, on average, objects with mass $\lesssim 0.3 \%$ of the soliton mass ($M_{\rm sol}$) are expelled from the soliton in $\sim 3\ \rm Gyr$, after which they continue their outward diffusion due to gravitational interactions with the soliton and the halo granules. More massive objects ($\gtrsim 1 \% M_{\rm sol}$), while executing a random walk, remain largely confined to the soliton due to dynamical friction. We also present an effective treatment of the diffusion, based on kinetic theory, that accurately reproduces the outward motion of low mass objects and briefly discuss how the observed displacements of star clusters and active galactic nuclei from the centers of their host galaxies can be used to constrain FDM., Comment: 20 pages, 9 figures, minor revisions, accepted for publication in ApJ

Fuzzy dark matter (FDM) has been a promising alternative to standard cold dark matter. The model consists of ultralight bosons with mass $m_b \sim 10^{-22}$ eV and features a quantum-pressure-supported solitonic core that oscillates. In this work, we show that the soliton density oscillations persist even after significant tidal stripping of the outer halo. We report two intrinsic yet distinct timescales associated, respectively, with the ground-state soliton wavefunction $\tau_{00}$ and the soliton density oscillations $\tau_\text{soliton}$, obeying $\tau_\text{soliton} /\tau_{00} \simeq 2.3$. The central star cluster (SC) in Eridanus II has a characteristic timescale $\tau_\text{soliton} / \tau_\text{SC} \sim 2$ to $3$ that deviates substantially from unity. As a result, we demonstrate, both analytically and numerically with three-dimensional self-consistent FDM simulations, that the gravitational heating of the SC owing to soliton density oscillations is negligible irrespective of $m_b$. We also show that the subhalo mass function to form Eridanus II does not place a strong constraint on $m_b$. These results are contrary to the previous findings by Marsh & Niemeyer (2019)., Comment: 19 pages, 19 figures; Accepted for publication in PRD

The evolution of cosmic perturbations in a two-component ultralight axionic universe is investigated. We present the first spectral computation of perturbations in multi-component universes. A particular case composed of light extreme axions and free massive particles offers a possibility for the formation of very high-redshift massive galaxies, which are typically required to host massive early quasars. Our computation retains the information of perturbed velocities for individual axion components, opening a new avenue for setting up initial conditions for future axion dark matter simulations., Comment: 10 pages, 8 figures

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 yields ubiquitous flux ratio anomalies of a few tens of percent, as is typically found for lensed quasars, and may also produce rare hexad and octad images, for sources located in well-defined caustic zones. 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., Comment: 6 pages, 4 figures

Wave dark matter ($\psi$DM) predicts a compact soliton core and a granular halo in every galaxy. This work presents the first simulation study of an elliptical galaxy by including both stars and $\psi$DM, focusing on the systematic changes of the central soliton and halo granules. With the addition of stars in the inner halo, we find the soliton core consistently becomes more prominent by absorbing mass from the host halo than that without stars, and the halo granules become "non-isothermal", "hotter" in the inner halo and "cooler" in the outer halo, as opposed to the isothermal halo in pure $\psi$DM cosmological simulations. Moreover, the composite (star+$\psi$DM) mass density is found to follow a $r^{-2}$ isothermal profile near the half-light radius in most cases. Most striking is the velocity dispersion of halo stars that increases rapidly toward the galactic center by a factor of at least 2 inside the half-light radius caused by the deepened soliton gravitational potential, a result that compares favorably with observations of elliptical galaxies and bulges in spiral galaxies. However in some rare situations we find a phase segregation turning a compact distribution of stars into two distinct populations with high and very low velocity dispersions; while the high-velocity component mostly resides in the halo, the very low-velocity component is bound to the interior of the soliton core, resembling stars in faint dwarf spheroidal galaxies., Comment: 15 pages, 11 figures

Wave dark matter ($\psi \rm{DM}$) composed of extremely light bosons ($m_{\psi} \sim 10^{-22}\,\rm eV$), with quantum pressure suppressing structures below a kpc-scale de Broglie wavelength, has become a viable dark matter candidate. Compared to the conventional free-particle $\psi {\rm DM}$ (${\rm FP} \psi {\rm DM}$), the extreme-axion $\psi \rm{DM}$ model (${\rm EA} \psi {\rm DM}$) proposed by Zhang & Chiueh (2017) features a larger cut-off wavenumber and a broad spectral bump in the matter transfer function. Here we conduct cosmological simulations to compare the halo abundances and assembly histories at $z=4-11$ between three different scenarios: ${\rm FP} \psi {\rm DM}$, ${\rm EA} \psi {\rm DM}$, and cold dark matter (CDM). We show that ${\rm EA} \psi {\rm DM}$ produces significantly more abundant low-mass haloes than ${\rm FP} \psi {\rm DM}$ with the same $m_{\psi}$, and therefore could alleviate the tension in $m_{\psi}$ required by the Ly$\alpha$ forest data and by the kpc-scale dwarf galaxy cores. We also find that, compared to the CDM counterparts, massive ${\rm EA} \psi {\rm DM}$ haloes are on average $3-4$ times more massive at $z=10-11$ due to their earlier formation, undergo a slower mass accretion at $7 \lesssim z \lesssim 11$, and then show a rapidly rising major merger rate exceeding CDM by $\sim 50\%$ at $4 \lesssim z \lesssim 7$. This fact suggests that ${\rm EA} \psi {\rm DM}$ haloes may exhibit more prominent starbursts at $z \lesssim 7$., Comment: 5 pages, 4 figures, only minor updates (accepted for publication in MNRAS letters)

Although still under debate, observations generally suggest that dwarf spheroidal (dSph) galaxies exhibit large constant-density cores in the centers, which can hardly be explained by dissipationless cold dark matter simulations without baryonic feedback. Wave dark matter (${\psi {\rm DM}}$), characterized by a single parameter, the dark matter particle mass $m_{\psi}$, predicts a central soliton core in every galaxy arising from quantum pressure against gravity. Here we apply Jeans analysis assuming a soliton core profile to the kinematic data of eight classical dSphs so as to constrain $m_{\psi}$, and obtain $m_{\psi}=1.18_{-0.24}^{+0.28}\times10^{-22}{\,\rm eV}$ and $m_{\psi}=1.79_{-0.33}^{+0.35}\times10^{-22}{\,\rm eV}~(2\sigma)$ using the observational data sets of Walker et al. (2007) and Walker et al. (2009b), respectively. We show that the estimate of $m_{\psi}$ is sensitive to the dSphs kinematic data sets and is robust to various models of stellar density profile. We also consider multiple stellar subpopulations in dSphs and find consistent results. This mass range of $m_{\psi}$ is in good agreement with other independent estimates, such as the high-redshift luminosity functions, the reionization history, and the Thomson optical depth to the cosmic microwave background., Comment: 12 pages, 11 figures, 2 tables, add appendix B about the joint analysis, submitted to MNRAS

Simulations of ultralight, $\sim 10^{-22}\,\rm eV$, bosonic dark matter exhibit rich wave-like structure, including a soliton core within a surrounding halo that continuously self-interferes on the de Broglie scale. We show here that as an inherent consequence, the soliton undergoes a confined random walk at the base of the halo potential. This is significant for the fate of the ancient central star cluster in Eridanus II, as the agitated soliton gravitationally shakes the star cluster in and out of the soliton on a time scale of $\sim 100\,\rm Myr$, so complete tidal disruption of the star cluster can occur within $\sim 1\,\rm Gyr$. This destructive effect can be mitigated by tidal stripping of the halo of Eridanus II, thereby reducing the agitation, depending on its orbit around the Milky Way. Our simulations show the Milky Way tide affects the halo much more than the soliton, so the star cluster in Eridanus II can survive for over $5\,\rm Gyr$ within the soliton if it formed after significant halo stripping., Comment: Accepted for publication in Phys. Rev. Lett. 6 pages, 4 figures

Accurate and precise measurements of the Hubble constant are critical for testing our current standard cosmological model and revealing possibly new physics. With Hubble Space Telescope (HST) imaging, each strong gravitational lens system with measured time delays can allow one to determine the Hubble constant with an uncertainty of $\sim 7\%$. Since HST will not last forever, we explore adaptive-optics (AO) imaging as an alternative that can provide higher angular resolution than HST imaging but has a less stable point spread function (PSF) due to atmospheric distortion. To make AO imaging useful for time-delay-lens cosmography, we develop a method to extract the unknown PSF directly from the imaging of strongly lensed quasars. In a blind test with two mock data sets created with different PSFs, we are able to recover the important cosmological parameters (time-delay distance, external shear, lens mass profile slope, and total Einstein radius). Our analysis of the Keck AO image of the strong lens system RXJ1131-1231 shows that the important parameters for cosmography agree with those based on HST imaging and modeling within 1-$\sigma$ uncertainties. Most importantly, the constraint on the model time-delay distance by using AO imaging with $0.045"$resolution is tighter by $\sim 50\%$ than the constraint of time-delay distance by using HST imaging with $0.09"$when a power-law mass distribution for the lens system is adopted. Our PSF reconstruction technique is generic and applicable to data sets that have multiple nearby point sources, enabling scientific studies that require high-precision models of the PSF., Comment: 20 pages, 15 figures

Using cosmological particle hydrodynamical simulations and uniform ultraviolet backgrounds, we compare Lyman-$\alpha$ forest flux spectra predicted by the conventional cold dark matter (CDM) model, the free-particle wave dark matter (FP$\psi$DM) model and extreme-axion wave dark matter (EA$\psi$DM) models of different initial axion field angles against the BOSS Lyman-$\alpha$ forest absorption spectra with a fixed boson mass $m_b\sim 10^{-22}$eV. We recover results reported previously (Ir\v{s}i\v{c} et al. 2017b; Armengaud et al. 2017) that the CDM model agrees better with the BOSS data than the FP$\psi$DM model by a large margin, and we find the difference of total $\chi^2$'s is $120$ for $420$ data bins. These previous results demand a larger boson mass by a factor $>10$ to be consistent with the date and are in tension with the favoured value determined from local satellite galaxies. We however find that such tension is removed as some EA$\psi$DM models predict Lyman-$\alpha$ flux spectra agreeing better with the BOSS data than the CDM model, and the difference of total $\chi^2$'s can be as large as $24$ for the same bin number. This finding arises with no surprise since EA$\psi$DM models have unique spectral shapes with spectral bumps in excess of the CDM power near the small-scale cutoff typical of $\psi$DM linear matter power spectra as well as more extended cutoffs than FP$\psi$DM (Zhang & Chiueh 2017a,b).

Acting as powerful gravitational lenses, the strong lensing galaxy clusters of the deep Hubble Frontier Fields (HFF) program permit access to lower-luminosity galaxies lying at higher redshifts than hitherto possible. We analyzed the HFF to measure the volume density of Lyman-break galaxies at $z > 4.75$ by identifying a complete and reliable sample up to $z \simeq 10$. A marked deficit of such galaxies was uncovered in the highly magnified regions of the clusters relative to their outskirts, implying that the magnification of the sky area dominates over additional faint galaxies magnified above the flux limit. This negative magnification bias is consistent with a slow rollover at the faint end of the UV luminosity function, and indicates a preference for Bose-Einstein condensate dark matter with a light boson mass of $m_\mathrm{B} \simeq 10^{-22} \, \mathrm{eV}$ over standard cold dark matter. We emphasize that measuring the magnification bias requires no correction for multiply lensed images (with typically three or more images per source), whereas directly reconstructing the luminosity function will lead to an overestimate unless such images can be exhaustively matched up, especially at the faint end that is accessible only in the strongly lensed regions. In addition, we detected a distinctive downward transition in galaxy number density at $z \gtrsim 8$, which may be linked to the relatively late reionization reported by Planck. Our results suggests that JWST will likely peer into an "abyss" with essentially no galaxies detected in deep NIR imaging at $z > 10$., 46 pages, 23 figures, accepted for publication in ApJ

The cold dark matter (CDM) paradigm successfully explains the cosmic structure over an enormous span of redshifts. However, it fails when probing the innermost regions of dark matter halos and the properties of the Milky Way’s dwarf galaxy satellites. Moreover, the lack of experimental detection of Weakly Interacting Massive Particle (WIMP) favors alternative candidates such as light axionic dark matter that naturally arise in string theory. Cosmological N-body simulations have shown that axionic dark matter forms a solitonic core of size of ≃150 pc in the innermost region of the galactic halos. The oscillating scalar field associated to the axionic dark matter halo produces an oscillating gravitational potential that induces a time dilation of the pulse arrival time of ≃400 ns/(m B /10 −22 eV) for pulsar within such a solitonic core. Over the whole galaxy, the averaged predicted signal may be detectable with current and forthcoming pulsar timing array telescopes. IDM acknowledges financial support from University of the Basque Country UPV/EHU under the program "Convocatoria de contratacion para la especializacion de personal investigador doctor en la UPV/EHU 2015", from the Spanish Ministerio de Economia y Competitividad through research project FIS2010-15492, and from the Basque Government through research project IT-956-16. RL is supported by the Spanish Ministry of Economy and Competitiveness through research projects FIS2010-15492 and Consolider EPI CSD2010-00064, and by the University of the Basque Country UPV/EHU under program UFI 11/55. IDM and RL also acknowledge support from the COST Action CA1511 Cosmology and Astrophysics Network for Theoretical Advances and Training Actions (CANTATA). TJB acknowledges generous hospitality from the Institute for Advanced Studies in Hong Kong and helpful conversations with Nick Kaiser and Kfir Blum. SHHT is supported by CRF Grant HKUST4/CRF/13G and GRF 16305414 issued by the Research Grants Council (RGC) of the Government of the Hong Kong SAR. Chipbond Technology Corporation is acknowledged for donating the GPU cluster with which this work was conducted. This work was supported in part by the National Science Council of Taiwan under grants NSC100-2112-M-002-018-MY3 and NSC99-2112-M-002-009-MY3

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., Comment: 14 pages, 7 figures, accepted by PRD

A 17.4-GHz 1-b sampling spectrum analyzer is built to demonstrate the working principle for a much faster spectrometer, where noise dithering and sine transformation are two key features of the principle. While the noise dithering can suppress the harmonics of signals, the sine transformation is shown to restore the noise. Measurements of this instrument reveal that the newly proposed novel scheme works as it should, with no technical showstopper for a much faster device.

Linear perturbations of the wave dark matter, or $\psi$ dark matter ($\psi$DM), of particle mass $\sim 10^{-22}$eV in the radiation-dominant era are analyzed, and the matter power spectrum at the photon-matter equality is obtained. We identify four phases of evolution for $\psi$DM perturbations, where the dynamics can be vastly different from the counterparts of cold dark matter (CDM). While in late stages after mass oscillation long-wave $\psi$DM perturbations are almost identical to CDM perturbations, some subtle differences remain, let alone intermediate-to-short waves that bear no resemblance with those of CDM throughout the whole evolutionary history. The dissimilarity is due to quantum mechanical effects which lead to severe mode suppression. We also discuss the axion model with a cosine field potential. The power spectrum of axion models are generally almost identical to those of $\psi$DM, but in the extreme case when the initial axion angle is near the field potential top, this axion model predict a power excess over a range of wave number and a higher spectral cutoff than $\psi$DM as if $\psi$DM had a higher particle mass., Comment: 20 pages, 7 figures

Light Axionic Dark Matter, motivated by string theory, is increasingly favored for the "no-WIMP era". Galaxy formation is suppressed below a Jeans scale, of $\simeq 10^8 M_\odot$ by setting the axion mass to, $m_B \sim 10^{-22}$eV, and the large dark cores of dwarf galaxies are explained as solitons on the de-Broglie scale. This is persuasive, but detection of the inherent scalar field oscillation at the Compton frequency, $\omega_B= (2.5{\rm \, months})^{-1}(m_B/10^{-22}eV)$, would be definitive. By evolving the coupled Schr\"odinger-Poisson equation for a Bose-Einstein condensate, we predict the dark matter is fully modulated by de-Broglie interference, with a dense soliton core of size $\simeq 150pc$, at the Galactic center. The oscillating field pressure induces General Relativistic time dilation in proportion to the local dark matter density and pulsars within this dense core have detectably large timing residuals, of $\simeq 400nsec/(m_B/10^{-22}eV)$. This is encouraging as many new pulsars should be discovered near the Galactic center with planned radio surveys. More generally, over the whole Galaxy, differences in dark matter density between pairs of pulsars imprints a pairwise Galactocentric signature that can be distinguished from an isotropic gravitational wave background., Comment: 6 pages, 3 figures. Accepted for publication in Physics Review Letter

Sub-horizon perturbations under the extreme initial condition of the axion model are investigated, where initial axion angles start near the potential maximum. This work focuses on a few new features found in the extreme axion model but absent in the free-particle model. A particularly novel new feature is the spectral excess relative to the CDM model in some wave number range, where the excess may be so large that landscapes of high-redshift universe beyond $z=10$ can be significantly altered. For axions of particle mass $10^{-22}$ eV, this range of wave number corresponds to first galaxies of few times $10^9-10^{10} M_\odot$. We demonstrate that sub-horizon perturbations are accurately described by Mathieu's equation and subject to parametric instability, which explains this novel feature. Actually the axion model is not a special one; perturbations in a wide range of scalar field models can share the similar characteristic., 11 pages, 3 figures

Radio observations suggest that 3C 75, located in the dumbbell shaped galaxy NGC 1128 at the center of Abell 400, hosts two colliding jets. Motivated by this source, we perform three-dimensional hydrodynamical simulations using a modified version of the GPU-accelerated Adaptive-MEsh-Refinement hydrodynamical parallel code ($\mathit{GAMER}$) to study colliding extragalactic jets. We find that colliding jets can be cast into two categories: 1) bouncing jets, in which case the jets bounce off each other keeping their identities, and 2) merging jets, when only one jet emerges from the collision. Under some conditions the interaction causes the jets to break up into oscillating filaments of opposite helicity, with consequences for their downstream stability. When one jet is significantly faster than the other and the impact parameter is small, the jets merge; the faster jet takes over the slower one. In the case of merging jets, the oscillations of the filaments, in projection, may show a feature which resembles a double helix, similar to the radio image of 3C 75. Thus we interpret the morphology of 3C 75 as a consequence of the collision of two jets with distinctly different speeds at a small impact parameter, with the faster jet breaking up into two oscillating filaments.

A new approximation method for inverting the Poisson's equation is presented for a continuously distributed and finite-sized source in an unbound domain. The advantage of this image multipole method arises from its ability to place the computational error close to the computational domain boundary, making the source region almost error free. It is contrasted to the modified Green's function method that has small but finite errors in the source region. Moreover, this approximation method also has a systematic way to greatly reduce the errors at the expense of somewhat greater computational efforts. Numerical examples of three-dimensional and two-dimensional cases are given to illustrate the advantage of the new method., 7 pages, 6 figures

We present a list of galaxy-scale lens candidates including a highly probable interacting galaxy-scale lens in the Hyper Suprime-Cam (HSC) imaging survey. We combine HSC imaging with the blended-spectra catalog from the Galaxy And Mass Assembly (GAMA) survey to identify lens candidates, and use lens mass modeling to confirm the candidates. We find 46 matches between the HSC S14A_0b imaging data release and the GAMA catalog. Ten of them are probable lens systems according to their morphology and redshifts. There is one system with an interacting galaxy pair, HSC J084928+000949, that has a valid mass model. We predict the total mass enclosed by the Einstein radius of $\sim0.72$" ($\sim1.65$kpc) for this new expected lens system to be $\sim10^{10.59}M_{\odot}$. Using the photometry in the {\it grizy} bands of the HSC survey and stellar population synthesis modeling with a Salpeter stellar initial mass function, we estimate the stellar mass within the Einstein radius to be $\sim10^{10.46}\,M_{\odot}$. We thus find a dark matter mass fraction within the Einstein radius of $\sim25\%$. Further spectroscopy or high-resolution imaging would allow confirmation of the nature of these lens candidates. The particular system with the interacting galaxy pair, if confirmed, would provide an opportunity to study the interplay between dark matter and stars as galaxies build up through hierarchical mergers., 10 pages, 5 figures

We present the implementation and performance of a class of directionally unsplit Riemann-solver-based hydrodynamic schemes on Graphic Processing Units (GPU). These schemes, including the MUSCL-Hancock method, a variant of the MUSCL-Hancock method, and the corner-transport-upwind method, are embedded into the adaptive-mesh-refinement (AMR) code GAMER. Furthermore, a hybrid MPI/OpenMP model is investigated, which enables the full exploitation of the computing power in a heterogeneous CPU/GPU cluster and significantly improves the overall performance. Performance benchmarks are conducted on the Dirac GPU cluster at NERSC/LBNL using up to 32 Tesla C2050 GPUs. A single GPU achieves speed-ups of 101(25) and 84(22) for uniform-mesh and AMR simulations, respectively, as compared with the performance using one(four) CPU core(s), and the excellent performance persists in multi-GPU tests. In addition, we make a direct comparison between GAMER and the widely-adopted CPU code Athena (Stone et al. 2008) in adiabatic hydrodynamic tests and demonstrate that, with the same accuracy, GAMER is able to achieve two orders of magnitude performance speed-up., Comment: 16 pages, 5 figures, accepted for publication in International Journal of High Performance Computing Applications (IJHPCA)

We study the environment of wet, dry, and mixed galaxy mergers at 0.75 < z < 1.2 using close pairs in the DEEP2 Galaxy Redshift Survey. We find that the typical environment of mixed and dry merger candidates is denser than that of wet mergers, mostly due to the color-density relation. While the galaxy companion rate (Nc) is observed to increase with overdensity, using N-body simulations we find that the fraction of pairs that will eventually merge decreases with the local density, predominantly because interlopers are more common in dense environments. After taking into account the merger probability of pairs as a function of local density, we find only marginal environment dependence of the fractional merger rate for wet mergers over the redshift range we have probed. On the other hand, the fractional dry merger rate increases rapidly with local density due to the increased population of red galaxies in dense environments. We also find that the environment distribution of K+A galaxies is similar to that of wet mergers alone and of wet+mixed mergers, suggesting a possible connection between K+A galaxies and wet and/or wet+mixed mergers. We conclude that, as early as z ~ 1, high-density regions are the preferred environment in which dry mergers occur, and that present-day red-sequence galaxies in overdense environments have, on average, undergone 1.2+-0.3 dry mergers since this time, accounting for (38+-10)% of their mass accretion in the last 8 billion years. Our findings suggest that dry mergers are crucial in the mass-assembly of massive red galaxies in dense environments, such as Brightest Cluster Galaxies (BCGs) in galaxy groups and clusters., 14 pages, 7 figures, 2 tables, submitted to ApJ, revised according to the referee's comments

GAMER, a parallel Graphic-processing-unit-accelerated Adaptive-MEsh-Refinement hydrodynamic code, has been extended to support magnetohydrodynamics (MHD) with both the corner-transport-upwind (CTU) and MUSCL-Hancock schemes and the constraint transport (CT) technique. The divergent preserving operator for adaptive mesh refinement (AMR) has been applied to reinforce the divergence-free constraint on the magnetic field. GAMER-MHD has fully exploited the concurrent executions between the GPU MHD solver and other CPU computation pertinent to AMR. We perform various standard tests to demonstrate that GAMER-MHD is both second-order accurate and robust, producing results as accurate as those given by high-resolution uniform-grid runs. We also explore a new 3D MHD test, where the magnetic field assumes the Arnold-Beltrami-Childress (ABC) configuration, temporarily becomes turbulent with current sheets and finally settles to a lowest-energy equilibrium state. This 3D problem is adopted for the performance test of GAMER-MHD. The single-GPU performance reaches $1.2\times 10^8$ and $5.5\times 10^7$ cell-updates/sec for the single- and double-precision calculations, respectively, on Tesla P100. We also demonstrate a parallel efficiency of $\sim 70\%$ for both weak and strong scaling using $1,024$ XK nodes on the Blue Waters supercomputers., Comment: 29 pages, 22 figures, 1 table. GAMER-MHD will soon become available at https://github.com/gamer-project/gamer

We present the newly developed code, GAMER (GPU-accelerated Adaptive MEsh Refinement code), which has adopted a novel approach to improve the performance of adaptive mesh refinement (AMR) astrophysical simulations by a large factor with the use of the graphic processing unit (GPU). The AMR implementation is based on a hierarchy of grid patches with an oct-tree data structure. We adopt a three-dimensional relaxing TVD scheme for the hydrodynamic solver, and a multi-level relaxation scheme for the Poisson solver. Both solvers have been implemented in GPU, by which hundreds of patches can be advanced in parallel. The computational overhead associated with the data transfer between CPU and GPU is carefully reduced by utilizing the capability of asynchronous memory copies in GPU, and the computing time of the ghost-zone values for each patch is made to diminish by overlapping it with the GPU computations. We demonstrate the accuracy of the code by performing several standard test problems in astrophysics. GAMER is a parallel code that can be run in a multi-GPU cluster system. We measure the performance of the code by performing purely-baryonic cosmological simulations in different hardware implementations, in which detailed timing analyses provide comparison between the computations with and without GPU(s) acceleration. Maximum speed-up factors of 12.19 and 10.47 are demonstrated using 1 GPU with 4096^3 effective resolution and 16 GPUs with 8192^3 effective resolution, respectively.

An alternative bosonic dark matter model is examined in detail via high-resolution simulations. These bosons have particle mass of order $10^{-22}eV$ and are non-interacting. If they do exist and can account for structure formation, these bosons must be condensed into the Bose-Einstein state and described by a coherent wave function. This matter, also known as Fuzzy Dark Matter (Hu, Barkana & Gruzinov 2000),, is speculated to be able, first, to eliminate the sub-galactic halos to solve the problem of over-abundance of dwarf galaxies, and, second, to produce flat halo cores in galaxies suggested by some observations. Our simulation results show that although this extremely light bosonic dark matter indeed suppresses low-mass halos, it can, to the contrary of expectation, yield singular halo cores. The density profile of the singular halo is almost identical to the halo profile of Navarro, Frenk & White (1997). Such a profile seems to be universal, in that it can be produced via either accretion or merger., 21 pages, 10 figures

The Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) started scientific operation in early 2007. This work describes the optimization of the system performance for the measurements of the Sunyaev-Zel'dovich effect for six massive galaxy clusters at redshifts $0.09 - 0.32$. We achieved a point source sensitivity of $63\pm 7$ mJy with the seven 0.6m dishes in 1 hour of on-source integration in 2-patch differencing observations. We measured and compensated for the delays between the antennas of our platform-mounted interferometer. Beam switching was used to cancel instrumental instabilities and ground pick up. Total power and phase stability were good on time scales of hours, and the system was shown to integrate down on equivalent timescales of 300 hours per baseline/correlation, or about 10 hours for the entire array. While the broadband correlator leads to good sensitivity, the small number of lags in the correlator resulted in poorly measured bandpass response. We corrected for this by using external calibrators (Jupiter and Saturn). Using Jupiter as the flux standard, we measured the disk brightness temperature of Saturn to be $149^{+5}_{-12}$ K., 9 pages, 7 figures, 1 table, accepted for publication in ApJ

We present observations, analysis, and results for the first-year operation of Array for Microwave Background Anisotropy (AMiBA), an interferometric experiment designed to study cosmology via the measurement of cosmic microwave background (CMB). AMiBA is the first CMB interferometer operating at 3 mm to have reported successful results, currently with seven close-packed antennas of 60 cm diameter giving a synthesized resolution of around 6'. During 2007, AMiBA detected the Sunyaev-Zel'dovich effects (SZEs) of six galaxy clusters at redshift 0.091 ≤ z ≤ 0.322. An observing strategy with on-off-source switching is used to minimize the effects from electronic offset and ground pickup. Planets were used to test the observational capability of AMiBA and to calibrate the conversion from correlator time-lag data to visibilities. The detailed formalism for data analysis is given. We summarize our early tests including observations of planets and quasars, and present images, visibility profiles, the estimated central coordinates, sizes, and SZE amplitudes of the galaxy clusters. Scientific implications are summarized. We also discuss possible systematic effects in the results.

AMiBA is the largest hexapod astronomical telescope in current operation. We present a description of this novel hexapod mount with its main mechanical components -- the support cone, universal joints, jack screws, and platform -- and outline the control system with the pointing model and the operating modes that are supported. The AMiBA hexapod mount performance is verified based on optical pointing tests and platform photogrammetry measurements. The photogrammetry results show that the deformations in the inner part of the platform are less than 120 micron rms. This is negligible for optical pointing corrections, radio alignment and radio phase errors for the currently operational 7-element compact configuration. The optical pointing error in azimuth and elevation is successively reduced by a series of corrections to about 0.4 arcmin rms which meets our goal for the 7-element target specifications., Comment: Accepted for ApJ, 33 pages, 15 figures

We discuss the observation, analysis, and results of the first-year science operation of AMiBA, an interferometric experiment designed to study cosmology via the measurement of Cosmic Microwave Background (CMB). In 2007, we successfully observed 6 galaxy clusters (z < 0.33) through the Sunyaev-Zel'dovich effect. AMiBA is the first CMB interferometer operating at 86–102 GHz, currently with 7 close-packed antennas of 60 cm in diameter giving a synthesized resolution of around 6 arcminutes. An observing strategy with on-off-source modulation is used to remove the effects from electronic offset and ground pickup. Formalism of the analysis is given and preliminary science results are summarized. Tests for systematic effects are also addressed. We also discuss the expansion plan.

The Array for Microwave Background Anisotropy is a 7-element interferometer sited on Mauna Loa, Hawaii. The seven 60cm telescopes are mounted on a 6-meter platform, and operates at 3mm wavelength. In October 2006, the telescope was officially dedicated and renamed as the Y. T. Lee AMiBA. During 2007, scientific operations have begun, after a long process of calibration and testing. At the time of this meeting, six clusters of galaxies have been detected and mapped via the inverse Compton scattering of the Cosmic Microwave Background radiation, also known as the Sunyaev-Zel'dovich Effect.

Inspired by Barnard 68, a Bok globule, that undergoes stable oscillations, we perform multi-phase hydrodynamic simulations to analyze the stability of Bok globules. We show that a high-density soft molecular core, with an adiabatic index $\gamma$ = 0.7 embedded in a warm isothermal diffuse gas, must have a small density gradient to retain the stability. Despite being stable, the molecular core can still collapse spontaneously as it will relax to develop a sufficiently large density gradient after tens of oscillations, or a few $10^7$ years. However, during its relaxation, the core may abruptly collapse triggered by the impingement of small-amplitude, long-wavelength ($\sim$ 6 $-$ 36 pc) sound waves in the warm gas. This triggered collapse mechanism is similar to a sonoluminescence phenomenon, where underwater ultrasounds can drive air bubble coalescence. The collapse configuration is found to be different from both inside-out and outside-in models of low-mass star formation; nonetheless the mass flux is close to the prediction of the inside-out model. The condition and the efficiency for this core collapse mechanism are identified. Generally speaking, a broad-band resonance condition must be met, where the core oscillation frequency and the wave frequency should match each other within a factor of several. A consequence of our findings predicts the possibility of propagating low-mass star formation, for which collapse of cores, within a mass range short of one order of magnitude, takes place sequentially tracing the wave front across a region of few tens of pc over $10^7$ years., Comment: 9 pages, 8 figures

The observed far greater number of giant arcs than that predicted by the ΛCDM cosmology has been a long-standing puzzle. We adopt high-resolution cosmological simulations to assess this problem and quantitatively show that this issue can be resolved by inclusion of high-redshift source galaxies. Our ΛCDM results reveal a much higher giant-arc probability (>10-5) for source galaxies of z > 2 than that for z < 2. Together with the source redshift distribution derived from the NTT catalogs, we obtain the predicted giant-arc number counts, which turn out to be consistent with most observations. The fact that the high-redshift (z ≥ 1.5) galaxies are preferentially magnified by the foreground lenses with a surprisingly high efficiency explains why previous theoretical work substantially underestimated the giant-arc number density. In previous work, the source galaxies were assumed to be mostly around z = 1, an assumption commonly adopted in weak-lensing studies. In addition, although the lens mergers have been found to yield enhanced lensing cross sections, we find that merely 5% of giant arcs are associated with merging lenses. Our results support not only the ΛCDM cosmology but also the long-standing anticipation that galaxy clusters, mostly formed below z = 1, are powerful gravitational telescopes for probing the high-redshift protogalaxies. We illustrate how the patchy dust regions in interacting galaxies at z = 4 may appear as filamentary structures within a giant arc when observed at the submillimeter wave band.

We report an instantaneous heating scheme at z=2 that can produce bound objects obeying the observed mass-temperature and luminosity-temperature relations. Such heating reduces the Sunyaev–Zeldovich (SZ) flux by a factor 3–2 on scale of 4–6 arcminutes. It exacerbates the need for the matter-fluctuation normalization to σ8 to assume an exceedingly high value (~1.2), if the several-arcminute excess anisotropy detected by the CBI experiment is entirely caused by the SZ effect of inter-galactic hot gas. By contrast, we find that the several-arcminute excess detected by ACBAR can be consistent with the preheated SZ signals of σ8=1.

The Array for Microwave Background Anisotropy is a 7-element interferometer to be sited on Mauna Loa, Hawaii. The seven 1.2m telescopes are mounted on a 6-meter platform, and operates at 3mm wavelength. At the time of this meeting, the telescope is under construction at the Vertex factory in Germany. It is due to be delivered in the middle of 2004. A 2-element prototype instrument has already been deployed to Mauna Loa where initial tests are underway.

Due to the finite-size beam, interferometry observations of the CMB polarization tend to measure a mixture of E and B modes arising from uw-space smearing. We propose an interferometry observation strategy that can disentangle the polarization E and B modes in a clean manner. This strategy employs a similar principle as that of the single-dish observation proposed earlier, but requiring additional linear combinations of signals measured from different baselines. In this work, we consider a hexagonally packed interferometry array as an example to demonstrate the ultility and efficiency of the proposed method. The measurement from a unit cell of equilateral triangular baselines in interferometry measurement is found to be comparable to that from a single dish in single-dish measurements.

Mounting evidence from $x$-ray observations reveals that bound objects should receive some form of energy in the past injected from non-gravitaional sources. We report that an instantaneous heating scheme, for which gases in dense regions were subjected to a temperature jump of 1keV at $z=2$ whereas those in rarified regions remained intact, can produce bound objects obeying the observed mass-temperature and luminosity-temperature relations. Such preheating lowers the peak Sunyaev-Zeldovich (SZ) power by a factor of 2 and exacerbates the need for the normalization of matter fluctuations $\sigma_8$ to assume an extreme high value $(\sim 1.1)$ for the SZ signals to account for the excess anisotropy on 5-arcminute scale detected by the Cosmic Background Imager in the cosmic microwave background radiation.

Although it has previously been recognized that the CMB polarization measurement on a sky ring is able to separate the B and E-mode patterns, we show that the rotational symmetry of the CMB polarization measurement is indeed unique, in that it enables the B and E modes to be separately measured in a simple and clean manner. The separation of B and E modes can in principle be achieved even with a single-polarization detector of arbitrary beam pattern. Based on this premise, a specific observing strategy is suggested, where the telescope scans along a sky circle with the polarimeter axis always oriented to a fixed direction with respect to the radial direction of the sky circle. The observational strategy for B-mode measurements can further be refined by choosing appropriate beam sizes to suppress the accidental E-mode leakage. We present the expected variances of the measured E and B-mode surface-brightness fluctuations as functions of the sky annulus radius and beam size. The measured peak power is expected to be up to 10% of the total power contained in the CMB polarization fluctuations for both B and E modes. The temperature-polarization, or T-E, correlation can be measured by the present method as well. By adopting this strategy, the AMiBA experiment is expected to detect, in a single sky annulus, a 1\sigma$ B-mode signal in 340 hours.

Nonspherical collapse is modelled, under the Zeldovich approximation, by six-dimensional random walks of the initial deformation tensor field. The collapse boundary adopted here is a slightly-modified version of that proposed by Chiueh and Lee (2001). Not only the mass function agrees with the fitting formula of Sheth and Tormen (1999), but the bias function and conditional mass function constructed by this model are also found to agree reasonably well with the simulation results of Jing (1998) and Somerville et al. (2000), respectively. In particular, by introducing a small mass gap, we find a fitting formula for the conditional mass function, which works well even at small time intervals between parent and progenitor halos during the merging history., Comment: 28 pages, 11 figures, section 3.2 complemented, Fig.1 and Fig. 4 modified, version accepted for publication in ApJ