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.

We offer a standing wave explanation for the rising proper motions of stars at the center of the globular cluster 47-Tucanae, amounting to $\simeq 0.44\%$ of the total mass. We show this can be explained as a solitonic core of dark matter composed of light bosons, $ m \geq 10^{-18} eV $, corresponding to $ \leq 0.27 pc$, as an alternative to a single black hole (BH) or a concentration of stellar BH remnants proposed recently. This is particularly important as having a concentrated stellar BH remnant with the above radii is very challenging without the heavy core since the three body encounters would prevent the BHs to be that concentrated. We propose this core develops from dark matter captured in the deep gravitational potential of this globular cluster as it orbits the dark halo of our galaxy. This boson may be evidence for a second light axion, additional to a lighter boson of $10^{-22} eV$, favored for the dominant dark matter implied by the large dark cores of dwarf spheroidal galaxies. The identification of two such light bosonic mass scales favors the generic string theory prediction of a wide, discrete mass spectrum of axionic scalar fields.

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.

A wavelike solution for the non-relativistic universal dark matter (wave-DM) is rapidly gaining interest, following distinctive predictions of pioneering simulations of cosmic structure as an interference pattern of coherently oscillating bosons. A prominent solitonic standing wave is predicted at the center of every galaxy, representing the ground state, that has been identified with the wide, kpc scale dark cores of common dwarf-spheroidal galaxies, providing a boson mass of, $\simeq 10^{-22}$ eV. A denser soliton is predicted for Milky Way sized galaxies where momentum is higher, so the de Broglie scale of the soliton is smaller, $\simeq 100$ pc, of mass $\simeq 10^9 M_\odot$. Here we show the central motion of bulge stars in the Milky Way implies the presence of such a dark core, where the velocity dispersion rises inversely with radius to a maximum of $\simeq 130$ km/s, corresponding to an excess central mass of $\simeq 1.5\times 10^9 M_\odot$ within $\simeq 100$ pc, favouring a boson mass of $\simeq 10^{-22}$ eV. This quantitative agreement with such a unique and distinctive prediction is therefore strong evidence for a light bosonic solution to the long standing Dark Matter puzzle, such as the axions generic in String Theory.

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$.

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.

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.

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.

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.

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$.

We make use of a catalog of 1600 Pan-STARRS1 groups produced by the probability friends-of-friends algorithm to explore how the galaxy properties, i.e. the specific star formation rate (SSFR) and quiescent fraction, depend on stellar mass and group-centric radius. The work is the extension of Lin et al. (2014). In this work, powered by a stacking technique plus a background subtraction for contamination removal, a finer correction and more precise results are obtained than in our previous work. We find that while the quiescent fraction increases with decreasing group-centric radius the median SSFRs of star-forming galaxies in groups at fixed stellar mass drop slightly from the field toward the group center. This suggests that the major quenching process in groups is likely a fast mechanism. On the other hand, a reduction in SSFRs by ~0.2 dex is seen inside clusters as opposed to the field galaxies. If the reduction is attributed to the slow quenching effect, the slow quenching process acts dominantly in clusters. In addition, we also examine the density-color relation, where the density is defined by using a sixth-nearest neighbor approach. Comparing the quiescent fractions contributed from the density and radial effect, we find that the density effect dominates over the massive group or cluster galaxies, and the radial effect becomes more effective in less massive galaxies. The results support mergers and/or starvation as the main quenching mechanisms in the group environment, while harassment and/or starvation dominate in clusters.

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.

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.

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.

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.

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.

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.

A variety of methods have been proposed to define and to quantify galaxy environments. While these techniques work well in general with spectroscopic redshift samples, their application to photometric redshift surveys remains uncertain. To investigate whether galaxy environments can be robustly measured with photo-z samples, we quantify how the density measured with the nearest neighbor approach is affected by photo-z uncertainties by using the Durham mock galaxy catalogs in which the 3D real-space environments and the properties of galaxies are exactly known. Furthermore, we present an optimization scheme in the choice of parameters used in the 2D projected measurements which yield the tightest correlation with respect to the 3D real-space environments. By adopting the optimized parameters in the density measurements, we show that the correlation between the 2D projected optimized density and real-space density can still be revealed, and the color-density relation is also visible out to $z \sim 0.8$ even for a photo-z uncertainty ($\sigma_{\Delta_{z}/(1+z)}$) up to 0.06. We find that at the redshift $0.3 < z < 0.5$ a deep ($i \sim 25$) photometric redshift survey with $\sigma_{\Delta_{z}/(1+z)} = 0.02$ yields a comparable performance of small-scale density measurement to a shallower $i \sim$ 22.5 spectroscopic sample with $\sim$ 10% sampling rate. Finally, we discuss the application of the local density measurements to the Pan-STARRS1 Medium Deep survey, one of the largest deep optical imaging surveys. Using data from $\sim5$ square degrees of survey area, our results show that it is possible to measure local density and to probe the color-density relation with 3$\sigma$ confidence level out to $z \sim 0.8$ in the PS-MDS. The color-density relation, however, quickly degrades for data covering smaller areas.

The newly established luminosity functions of high-z galaxies at $4 \lesssim z \lesssim 10$ can provide a stringent check on dark matter models that aim to explain the core properties of dwarf galaxies. The cores of dwarf spheroidal galaxies are understood to be too large to be accounted for by free streaming of warm dark matter without overly suppressing the formation of such galaxies. Here we demonstrate with cosmological simulations that wave dark matter, $\psi$DM, appropriate for light bosons such as axions, does not suffer this problem, given a boson mass of $m_{\psi} \ge 1.2 \times 10^{-22}{\,\rm eV}$ ($2\sigma$). In this case, the halo mass function is suppressed below $\sim 10^{10}{\,M_\odot}$ at a level that is consistent with the high-z luminosity functions, while simultaneously generating the kpc-scale cores in dwarf galaxies arising from the solitonic ground state in $\psi$DM. We demonstrate that the reionization history in this scenario is consistent with the Thomson optical depth recently reported by Planck, assuming a reasonable ionizing photon production rate. We predict that the luminosity function should turn over slowly around an intrinsic UV luminosity of $M_{\rm UV} \gtrsim -16$ at $z \gtrsim 4$. We also show that for galaxies magnified $\mathord{>}10\times$ in the Hubble Frontier Fields, $\psi$DM predicts an order of magnitude fewer detections than cold dark matter at $z \gtrsim 10$ down to $M_{\rm UV} \sim -15$, allowing us to distinguish between these very different interpretations for the observed coldness of dark matter.

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. Current and upcoming imaging surveys will contain thousands of new lensed quasars, augmenting the existing sample by at least two orders of magnitude. To find such lens systems, we built a robot, CHITAH, that hunts for lensed quasars by modeling the configuration of the multiple quasar images. Specifically, given an image of an object that might be a lensed quasar, CHITAH first disentangles the light from the supposed lens galaxy and the light from the multiple quasar images based on color information. A simple rule is designed to categorize the given object as a potential four-image (quad) or two-image (double) lensed quasar system. The configuration of the identified quasar images is subsequently modeled to classify whether the object is a lensed quasar system. We test the performance of CHITAH using simulated lens systems based on the Canada-France-Hawaii Telescope Legacy Survey. For bright quads with large image separations (with Einstein radius $r_{\rm ein}>1.1"$) simulated using Gaussian point-spread functions, a high true-positive rate (TPR) of $\sim$90% and a low false-positive rate of $\sim$$3\%$ show that this is a promising approach to search for new lens systems. We obtain high TPR for lens systems with $r_{\rm ein}\gtrsim0.5"$, so the performance of CHITAH is set by the seeing. We further feed a known gravitational lens system, COSMOS 5921$+$0638, to CHITAH, and demonstrate that CHITAH is able to classify this real gravitational lens system successfully. Our newly built CHITAH is omnivorous and can hunt in any ground-based imaging surveys.

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.

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.

Scalar field dark matter likely is able to solve all small-scale cosmology problems facing the cold dark matter (CDM), and has become an emerging contender to challenge the CDM. It however requires a particle mass $\sim 1 - 2 \times10^{-22}$eV. We find such an extremely small particle mass can naturally arise from a non-QCD axion mechanism, under fairly general assumptions that a few species of self-interacting light particles of comparable masses and a massless gauge boson decouple from the bright sector since the photon temperature exceeds 200 GeV, and the axion is the dominant dark matter. These assumptions also set the axion decay constant scale to several $\times 10^{16}$ GeV. Given the above axion mass range, we further pin down the dark-sector particles to consist of only one species of fermion and anti-fermion, likely right-handed neutrinos. With a mass around $92-128$ eV, the dark-sector particles may constitute a minority population of dark matter. If the gauge boson lives on SU(2), a dilute instanton gas can contribute about $2.5\%$ of the total relativistic relics in the cosmic microwave background radiation.

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.

The conventional cold, particle interpretation of dark matter (CDM) still lacks laboratory support and struggles with the basic properties of common dwarf galaxies, which have surprisingly uniform central masses and shallow density profiles. In contrast, galaxies predicted by CDM extend to much lower masses, with steeper, singular profiles. This tension motivates cold, wavelike dark matter ($\psi$DM) composed of a non-relativistic Bose-Einstein condensate, so the uncertainty principle counters gravity below a Jeans scale. Here we achieve the first cosmological simulations of this quantum state at unprecedentedly high resolution capable of resolving dwarf galaxies, with only one free parameter, $\bf{m_B}$, the boson mass. We demonstrate the large scale structure of this $\psi$DM simulation is indistinguishable from CDM, as desired, but differs radically inside galaxies. Connected filaments and collapsed haloes form a large interference network, with gravitationally self-bound solitonic cores inside every galaxy surrounded by extended haloes of fluctuating density granules. These results allow us to determine $\bf{m_B=(8.1^{+1.6}_{-1.7})\times 10^{-23}~eV}$ using stellar phase-space distributions in dwarf spheroidal galaxies. Denser, more massive solitons are predicted for Milky Way sized galaxies, providing a substantial seed to help explain early spheroid formation. Suppression of small structures means the onset of galaxy formation for $\psi$DM is substantially delayed relative to CDM, appearing at $\bf{z\lesssim 13}$ in our simulations.

Aims: We performed a detailed photometric analysis of the lensed system UM673 (Q0142-100) and an analysis of the tentative lens models. Methods: High-resolution adaptive optics images of UM673 taken with the Subaru telescope in the H band were examined. We also analysed the J, H and K-band observational data of UM673 obtained with the 1.3m telescope at the CTIO observatory. Results: We present photometry of quasar components A and B of UM673, the lens, and the nearby bright galaxy using H-band observational data obtained with the Subaru telescope. Based on the CTIO observations of UM673, we also present J- and H-band photometry and estimates of the J, H and K-band flux ratios between the two UM673 components in recent epochs. The near-infrared fluxes of the A and B components of UM673 and their published optical fluxes are analysed to measure extinction properties of the lensing galaxy. We estimate the extinction-corrected flux ratio between components A and B to be about 2.14 mag. We discuss lens models for the UM673 system constrained with the positions of the UM673 components, their flux ratio, and the previously measured time delay

Measurements of the polarized radiation often reveal specific physical properties of emission sources, such as strengths and orientations of magnetic fields offered by synchrotron radiation and Zeeman line emission, and the electron density distribution by free-free emission. Polarization-capable, millimeter/sub-millimeter telescopes are normally equipped with either septum polarizers or ortho-mode transducers (OMT) for the detection of polarized radiation. While the septum polarizer is traditionally conceived to be limited to a significantly narrower bandwidth than the OMT, it does possess advantageous features for astronomical polarization measurements unparalleled by the OMT. Challenging the conventional bandwidth limit, we design an extremely wideband circular waveguide septum polarizer, covering $42\%$ bandwidth, from 77 GHz to 118 GHz, without any undesired resonance. Stokes parameters constructed from the measured data in between 77 GHz and 115 GHz show that the leakage from $I$ to $Q$ and $U$ is below $\pm 2\%$ and the $Q-U$ mutual leakage below $\pm 1\%$. Such a performance is comparable to other modern polarizers, but the bandwidth of this polarizer can be at least twice as wide. This extremely wide-band design removes the major weakness of the septum polarizer and opens up a new window for future astronomical polarization measurements.

(Abridged) In tandem with observational datasets, we utilize realistic mock catalogs, based on a semi-analytic galaxy formation model, constructed specifically for Pan-STARRS1 Medium Deep Surveys in order to assess the performance of the Probability Friends-of-Friends (PFOF, Liu et al.) group finder, and aim to develop a grouping optimization method applicable to surveys like Pan-STARRS1. Producing mock PFOF group catalogs under a variety of photometric redshift accuracies ({\sigma}{\Delta}z/(1+zs)), we find that catalog purities and completenesses from ``good' {\sigma}{\Delta}z/(1+zs)) ~ 0.01) to ``poor' {\sigma}{\Delta}z/(1+zs)) ~ 0.07) photo-zs gradually degrade respectively from 77% and 70% to 52% and 47%. To avoid model dependency of the mock for use on observational data we apply a ``subset optimization' approach, using spectroscopic-redshift group data from the target field to train the group finder for application to that field, as an alternative method for the grouping optimization. We demonstrate this approach using these spectroscopically identified groups as the training set, i.e. zCOSMOS groups for PFOF searches within PS1 Medium Deep Field04 (PS1MD04) and DEEP2 EGS groups for searches in PS1MD07. We ultimately apply PFOF to four datasets spanning the photo-z uncertainty range from 0.01 to 0.06 in order to quantify the dependence of group recovery performance on photo-z accuracy. We find purities and completenesses calculated from observational datasets broadly agree with their mock analogues. Further tests of the PFOF algorithm are performed via matches to X-ray clusters identified within the PS1MD04 and COSMOS footprints. Across over a decade in group mass, we find PFOF groups match ~85% of X-ray clusters in COSMOS and PS1MD04, but at a lower statistical significance in the latter.

We make use of four galaxy catalogs based on four different semi-analytical models (SAMs) implemented in the Millennium simulation to study the environmental effects and the model dependence of galaxy merger rate. We begin the analyses by finding that galaxy merger rate in the SAMs has mild redshift evolution, consistent with results of previous works. To study the environmental dependence of galaxy merger rate, we adopt two estimators, the local overdensity (1+{\delta}n) defined as the surface density from the nth-nearest-neighbor (n = 6 is chosen in this study) and the host halo mass Mh. We find that galaxy merger rate Fmg shows strong dependence on the local overdensity (1+{\delta}n) and the dependence is similar at all redshifts. For the overdensity estimator, the merger rate Fmg is found about twenty times larger in the densest regions than in under-dense ones in two of the four models while it is roughly four times higher in the other two. In other words, the discrepancies of the merger rate difference between two extremes can differ by a factor of ~ five depending on the SAMs adopted. On the other hand for the halo mass estimator, Fmg does not monotonically increase with the host halo mass Mh, but peaks in the $M_h$ range between 10^12 and 10^13 h-1 M{\Theta}, which corresponds to group environments. High merger rate in high local density regions corresponds primarily to the high merge rate in group environments......

We study brightness variations in the double lensed quasar UM673 (Q0142-100) with the aim of measuring the time delay between its two images. In the paper we combine our previously published observational data of UM673 obtained during the 2003 - 2005 seasons at the Maidanak Observatory with archival and recently observed Maidanak and CTIO UM673 data. We analyze the V, R and I-band light curves of the A and B images of UM673, which cover ten observational seasons from August 2001 to November 2010. We also analyze the time evolution of the difference in magnitudes between images A and B of UM673 over more than ten years. We find that the quasar exhibits both short-term (with amplitude of \sim 0.1 mag in the R band) and high-amplitude (\sim 0.3 mag) long-term variability on timescales of about several months and several years, respectively. These brightness variations are used to constrain the time delay between the images of UM673. From cross-correlation analysis of the A and B quasar light curves and error analysis we measure the mean time delay and its error of 89 \pm11 days. Given the input time delay of 88 days, the most probable value of the delay that can be recovered from light curves with the same statistical properties as the observed R-band light curves of UM673 is 95{+5/-16}{+14/-29} days (68 and 95 % confidence intervals). Analysis of the V - I color variations and V, R and I-band magnitude differences of the quasar images does not show clear evidence of the microlensing variations between 1998 and 2010.

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.

GAMER is a GPU-accelerated Adaptive-MEsh-Refinement code for astrophysical simulations. In this work, two further extensions of the code are reported. First, we have implemented the MUSCL-Hancock method with the Roe's Riemann solver for the hydrodynamic evolution, by which the accuracy, overall performance and the GPU versus CPU speed-up factor are improved. Second, we have implemented the out-of-core computation, which utilizes the large storage space of multiple hard disks as the additional run-time virtual memory and permits an extremely large problem to be solved in a relatively small-size GPU cluster. The communication overhead associated with the data transfer between the parallel hard disks and the main memory is carefully reduced by overlapping it with the CPU/GPU computations.

A wideband analog correlator has been constructed for the Yuan-Tseh Lee Array for Microwave Background Anisotropy. Lag correlators using analog multipliers provide large bandwidth and moderate frequency resolution. Broadband IF distribution, backend signal processing and control are described. Operating conditions for optimum sensitivity and linearity are discussed. From observations, a large effective bandwidth of around 10 GHz has been shown to provide sufficient sensitivity for detecting cosmic microwave background variations.

The Sunyaev-Zel'dovich Effect (SZE) has been observed toward six massive galaxy clusters, at redshifts 0.091 \leq z \leq 0.322 in the 86-102 GHz band with the Y. T. Lee Array for Microwave Background Anisotropy (AMiBA). We modify an iterative method, based on the isothermal \beta-models, to derive the electron temperature T_e, total mass M_t, gas mass M_g, and integrated Compton Y within r_2500, from the AMiBA SZE data. Non-isothermal universal temperature profile (UTP) \beta models are also considered in this paper. These results are in good agreement with those deduced from other observations. We also investigate the embedded scaling relations, due to the assumptions that have been made in the method we adopted, between these purely SZE-deduced T_e, M_t, M_g and Y. Our results suggest that cluster properties may be measurable with SZE observations alone. However, the assumptions built into the pure-SZE method bias the results of scaling relation estimations and need further study.

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.

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.

We investigate the scaling relations between the X-ray and the thermal Sunyaev-Zel'dovich Effect (SZE) properties of clusters of galaxies, using data taken during 2007 by the Y.T. Lee Array for Microwave Background Anisotropy (AMiBA) at 94 GHz for the six clusters A1689, A1995, A2142, A2163, A2261, and A2390. The scaling relations relate the integrated Compton-y parameter Y_{2500} to the X-ray derived gas temperature T_{e}, total mass M_{2500}, and bolometric luminosity L_X within r_{2500}. Our results for the power-law index and normalization are both consistent with the self-similar model and other studies in the literature except for the Y_{2500}-L_X relation, for which a physical explanation is given though further investigation may be still needed. Our results not only provide confidence for the AMiBA project but also support our understanding of galaxy clusters.

Galaxies can form in a sufficiently deep gravitational potential so that efficient gas cooling occurs. We estimate that such potential is provided by a halo of mass $M \gtsim M_{c} \approx 7.0 \times 10^{12} ~ (\Delta_{c}(z) (1+z)^{3})^{-1/2} \Msun$, where $\Delta_{c}(z)$ is the mean overdensity of spherically virialized objects formed at redshift $z$, and $M_{c} \approx 4.0 \times 10^{11} \Msun$ at $z = 0$. Based on this criterion, our galaxy samples are constructed from cosmology simulation data by using HiFOF to select subhalos in those FOF halos that are more massive than $M_{c}$. There are far more dark subhalos than galaxy-hosting subhalos. Several tests against observations have been performed to examine our galaxy samples, including the differential galaxy mass functions, the galaxy space density, the projected two point correlation functions (CF), the HODs, and the kinematic pair fractions. These tests show good agreements. Based on the consistency with observations, our galaxy sample is believed to correctly represent galaxies in real universe, and can be used to study other unexplored galaxy properties.

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.

The Y. T. Lee Array for Microwave Background (AMiBA) has reported the first science results on the detection of galaxy clusters via the Sunyaev Zel'dovich effect. The science objectives required small reflectors in order to sample large scale structures (20') while interferometry provided modest resolutions (2'). With these constraints, we designed for the best sensitivity by utilizing the maximum possible continuum bandwidth matched to the atmospheric window at 86-102GHz, with dual polarizations. A novel wide-band analog correlator was designed that is easily expandable for more interferometer elements. MMIC technology was used throughout as much as possible in order to miniaturize the components and to enhance mass production. These designs will find application in other upcoming astronomy projects. AMiBA is now in operations since 2006, and we are in the process to expand the array from 7 to 13 elements.

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.

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.

We present observations, analysis and results for the first-year operation of 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 (SZE) 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. Science results are summarized. We also discuss possible systematic effects in the results.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

We study the redshift evolution of galaxy pair fractions and merger rates for different types of galaxies using kinematic pairs selected from the DEEP2 Redshift Survey. By parameterizing the evolution of the pair fraction as (1+z)^{m}, we find that the companion rate increases mildly with redshift with m = 0.41+-0.20 for all galaxies with -21 < M_B^{e} < -19. Blue galaxies show slightly faster evolution in the blue companion rate with m = 1.27+-0.35 while red galaxies have had fewer red companions in the past as evidenced by the negative slope m = -0.92+-0.59. We find that at low redshift the pair fraction within the red sequence exceeds that of the blue cloud, indicating a higher merger probability among red galaxies compared to that among the blue galaxies. With further assumptions on the merger timescale and the fraction of pairs that will merge, the galaxy major merger rates for 0.1 < z <1.2 are estimated to be ~10^{-3}h^{3}Mpc^{-3}Gyr^{-1} with a factor of 2 uncertainty. At z ~ 1.1, 68% of mergers are wet, 8% of mergers are dry, and 24% of mergers are mixed, compared to 31% wet mergers, 25% dry mergers, and 44% mixed mergers at z ~ 0.1. The growth of dry merger rates with decreasing redshift is mainly due to the increase in the co-moving number density of red galaxies over time. About 22% to 54% of present-day L^{*} galaxies have experienced major mergers since z ~ 1.2, depending on the definition of major mergers. Moreover, 24% of the red galaxies at the present epoch have had dry mergers with luminosity ratios between 1:4 and 4:1 since z ~ 1. Our results also suggest that the wet mergers and/or mixed mergers may be partially responsible for producing red galaxies with intermediate masses while a significant portion of massive red galaxies is assembled through dry mergers at later times.

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.

Using data from the DEEP2 Galaxy Redshift Survey and HST/ACS imaging in the Extended Groth Strip, we select nearly 100 interacting galaxy systems including kinematic close pairs and morphologically identified merging galaxies. Spitzer MIPS 24 micron fluxes of these systems reflect the current dusty star formation activity, and at a fixed stellar mass (M_{*}) the median infrared luminosity (L_{IR}) among merging galaxies and close pairs of blue galaxies is twice (1.9 +/- 0.4) that of control pairs drawn from isolated blue galaxies. Enhancement declines with galaxy separation, being strongest in close pairs and mergers and weaker in wide pairs compared to the control sample. At z ~ 0.9, 7.1% +/- 4.3% of massive interacting galaxies (M_{*} > 2*10^{10} M_{solar}) are found to be ULIRGs, compared to 2.6% +/- 0.7% in the control sample. The large spread of IR luminosity to stellar mass ratio among interacting galaxies suggests that this enhancement may depend on the merger stage as well as other as yet unidentified factors (e.g., galaxy structure, mass ratio, orbital characteristics, presence of AGN or bar). The contribution of interacting systems to the total IR luminosity density is moderate (<= 36 %).

The technique of weak-lensing aperture mass densitometry, so called the zeta-statistic, has recently been popular in actual observations for measurement of individual cluster mass. It has however been anticipated that the line-of-sight projection by foreground and background matter can adversely affect the cluster mass determination with not only substantial error dispersion but also a sizable positive systematic bias. Additionally, the finite number of background galaxies even at a reasonable observing depth can also introduce Poisson noise to the mass estimate. In this paper, we quantitatively investigate the degree of errors separately contributed by the two sources to the mass determination of those galaxy clusters with M_{200}>10^{14}M_{\odot}. We find that the aperture mass of zeta-statistic turns out to be a mass estimator of much reduced systematic bias, due to the cancellation by the positively biased local background mass sheet. However, the error dispersion of M_{200} arising from both projection effect and Poisson noise is found to be still sizable (40%-90%), even for the shear-selected, clean sample where multiple clusters located within a suitable projected aperture are removed. We also investigate how to remedy this large-error problem in weak lensing measurements, and propose a plausible alternative mass estimator, M(<theta_{1000}), an aperture mass measured within about half the virial radius. The aperture mass M(<theta_{1000}) is free of bias and has a substantially reduced error dispersion, 39% for the worst case of high-z, low-mass clusters, that can be smaller than the error dispersion of M_{200} as much as a factor 3.

We present simulations of interferometric Sunyaev-Zel'dovich effect (SZE) and optical weak lenisng observations for the forthcoming AMiBA experiment, aiming at searching for high-redshift clusters of galaxies. On the basis of simulated sky maps, we have derived theoretical halo number counts and redshift distributions of selected halo samples for an AMiBA SZE survey and a weak lensing follow-up survey. By utilizing the conditional number counts of weak lensing halos with the faint SZE detection, we show that a combined SZE and weak lensing survey can gain an additional fainter halo sample at a given false positive rate, which cannot be obtained from either survey alone.

The cross-correlation of Sunyaev-Zel'dovich effect (SZ) and weak-lensing imaging surveys can be used to test how well hot baryons trace dark matter in clusters of galaxies. We examine this concept using mock SZ and weak-lensing surveys based on the forthcoming AMiBA experiment and generated from a pre-heated cosmological N-body/hydrodynamic simulation. A cross-correlation diagram between matched lensing convergence and Compton y maps exhibits butterfly-wing-like structures, corresponding to individual clusters, that encode rich information about the distributions of hot gas and dark matter. When the cluster redshift and temperature are available the slope of a wing reveals the cluster gas fraction and the width of the wing indicates how badly the hot gas traces dark matter. On the basis of simulated data we discuss systematic errors in the projected gas fraction estimates that would be obtained from such survey comparisons.

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.

We study bound object formation in a background universe accelerated by cosmological constant and by quintessence. If the acceleration lasts forever, due to the existence of event horizon, one would have naively expected the universe to approach a state of cold death. However, we find that many local regions in the universe can in fact be protected by their own gravity to form mini-universes, provided that their present matter densities exceed some critical value. In the case with cosmological constant ($\Lambda$CDM cosmology) the condition of forming mini-universe is that the ratio of present density parameters $\Omega_{i}^0/\Omega_\Lambda$ ought to be larger than a critical value 3.63. Such mini-universes typically weigh less than $2\times 10^{14}$ solar masses, with the lighter ones having tight and compact configurations. In the case with quintessence, the final ratio $\Omega_i/\Omega_q$ of a mini-universe is found to be always larger than $w_q^2-w_q$, where $w_q$ is the equation of state parameter.

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.

We derive the distribution of the dimensionless specific angular momentum of dark matter halos, $P(j)$, in the framework of the standard tidal torque theory, and explain the characteristic shape of $P(j)$ commonly observed in N-body simulations. A scalar quantity (shear scalar), $r$, is introduced for measuring the effective strength of the tidal torque force acting upon the halo from the surrounding matter. It is found that the ubiquitous and broad shape of $P(j)$ can mostly be attributed to the unique property of the shear scalar $r$, and also that $P(j)$ is insensitve to the underlying collapse dynamics. Our result demonstrates that although the shape of $P(j)$ is ubiquitous and close to log-normal, the distribution is not exactly log-normal but decays exponentially at the high angular momentum end, and drops slowly as a power-law with an index 2 at the low angular momentum end.

We examine the expected radio confusion on the thermal Sunyaev-Zel'dovich (SZ) effect of galaxy clusters at 28.5 GHz and 90 GHz based on the cluster radio luminosity function (CRLF) at 1.4 GHz. With the observationally determined spectral index (${\alpha}$, where $S_{\rm\nu}\propto\nu^{-\alpha}$) distribution, instead of a single average index $<\alpha>$, we convert the cluster radio luminosity function at 1.4 GHz to the high frequency ones and estimate the total radio flux in a cluster. At 28.5 GHz, radio confusion is up to 10~100% for small clusters ($M\sim 2\times 10^{14}M_{\odot}$) below redshift z=1, with more severe confusion for smaller and lower redshift clusters. By contrast, at 90 GHz the confusion is less than 10% for small clusters even at low redshifts.

Dynamics of the multi-component, multi-field quintessence and gravity is formulated as relativistic N-particle dynamics, embedded in a static viscus flat space and under the forces given by an interacting Lorentz scalar potential via exchange of field bosons. The Ratra-Peebles power-law potential of effective single-field quintessence can be derived from this microscopic perspective. In certain situations, the effective dynamics can be made identical to that of the single complex quintessence, except for that the overall U(1) symmetry is not manifestly broken. The present formulation provides a convenient gauge for analyzing the superhorizon perturbations and possibly for quantization of superhorizon fields and gravity.

We present a new collapse condition to describe the formation of dark halos via nonspherical gravitational clustering. This new nonspherical collapse condition is obtained by the logical generalization of the spherical model to the nonspherical one. By solving a diffusion-like random matrix equation with the help of the Monte Carlo method, we show that this nonspherical collapse condition yields the mass function derived by Sheth & Tormen (1999) which has been shown to be in excellent agreement with the recent N-body results of high resolution. We expect that this nonspherical collapse condition might provide us a deeper insight into the structure formation, and suggest that it should be widely applied to various cosmological issues such as the galaxy merging history, the galaxy bias, and so forth.

As part of a 4-year Cosmology and Particle Astrophysics (CosPA) Research Excellence Initiative in Taiwan, AMiBA $-$ a 19-element dual-channel 85-105 GHz interferometer array is being specifically built to search for high redshift clusters of galaxies via the Sunyaev-Zeldovich Effect (SZE). In addition, AMiBA will have full polarization capabilities, in order to probe the polarization properties of the Cosmic Microwave Background. AMiBA, to be sited on Mauna Kea in Hawaii or in Chile, will reach a sensitivity of $\sim 1$ mJy or 7$\mu$K in 1 hour. The project involves extensive international scientific and technical collaborations. The construction of AMiBA is scheduled to starting operating in early 2004.

The work examines the small-scale E and B patterns of the CMBR polarization anisotropy. Particularly we address the topological natures and spatial symmetries of the E and B modes, and how one may make use of the local measurements of Stokes parameters to separate these two modes. The analysis of a local map in the Fourier space for separating the E and B modes can be an ill-defined problem, due to the non-periodic boundary condition in the map for each mode. A strategy for the E-B separation in a local map can be through an appropriate projection of the polarization tensor field into a unique vector field, which naturally contains a curl-free E component and a divergence-free B component. An integral method, equivalent to the real-space, top-hat filtering, is proposed for extracting the smoothed "charge" and "current" of the E and B modes separately. The top-hat filter function exhibits, in its power spectrum, an oscillation of a well-defined frequency comparable to the filter size. It is pointed out that when the inverse filter size is chosen properly, matching the oscillation period in the E-mode power spectrum, the sensitivity for the E-mode detection can be made to be significantly improved.

We study Sunyaev-Zeldovich Effect (SZE) cluster counts in different cosmologies. It is found that even without the full knowledge of the redshift distribution of SZE clusters, one can still readily distinguish a flat universe with a cosmological constant from an open universe. We divide clusters into a low redshift group (with redshift z<0.5) and a high redshift group (with z>1), and compute the relation r=N(z<0.5)/N(z>1), where N is the number of flux-limited (S_lim) SZE clusters. With about the same total number of SZE clusters N(z>0), the r-value for a flat universe with a consmological constant and that for an open universe occupy different regions in the S_lim -- r plot for the most likely cosmological parameters 0.25<\Omega_0<0.35 and 0.2<\Gamma<0.3, where \Gamma is the shape parameter of the initial power spectrum of density fluctuations. Thus, with a deep SZE cluster survey, the ratio r can reveal, independent of the normalization of the power spectrum, whether we are living in a low-density flat universe or an open universe. Within the flat universe scenario, the SZE cluster-normalized \sigma_8 has also been investigated in this work, where sigma_8 is the r.m.s. density fluctuations within the top-hat scale of 8 h^{-1} Mpc. A functional relation sigma_8 \sim \Omega_0^{-0.13} is found. Combining it with the X-ray cluster-normalized \sigma_8\sim\Omega_0^{-0.52+0.13\Omega_0}, one can put tight constraints on both \Omega_0 and \sigma_8 simulataneously.

In this paper we have made an attempt to derive the radial profiles of density and temperature of intracluster gas based on the two well-established facts at present: the X-ray observed surface brightness of clusters described by the standard beta model and the (NFW) universal density profile as the underlying dark matter distribution. We have numerically solved the hydrostatic equation by demanding that the volumed-averaged baryon fraction of a cluster should asymptotically approach the universal value at its viral radius. We have shown that the radial temperature variation derived from these constraints differs significantly from the conventional polytropic equation of state: The gas temperature profile may show a dramatic increase or decrease with outward radius, depending sensitively on the beta parameter. A large beta value (typically greater than 0.8) is required in order to ensure that the X-ray temperature makes a drop at the virial radius. This indicates that either the NFW profile is inappropriate to apply for the overall dark matter distribution of clusters or other non-gravitational heating processes may play an important role in the dynamical evolution of clusters.

Increasing evidences provided primarily by the cluster lensing and numerical simulations of cluster formation indicate that galaxy clusters may contain compact cores that are substantially smaller than the cores revealed by the X-ray observations of hot intracluster gas. In this paper we present a model that describes how two distinct cores can grow simultaneously as a result of infall from the background dark matter. This model needs a pre-existing localized large fluctuation, which can be the non-Gaussian density peak of the primordial fluctuations, to seed the individual growing cluster. On the other hand, numerous recent observations also show that no strong evolution is detected for galaxy clusters within the redshift up to $z\sim 0.8$. We therefore present a comparison of observations with the saturated cluster configuration resulting from our infall model. In the saturated state, the predicted compact core mass is about few$\times 10^{11} M_\odot$ and the core size about 1 kpc for a cluster mass $\sim$ several $\times 10^{15} M_\odot$ within a radius of 3 Mpc. We have successfully reproduced the dark matter distribution revealed by the gravitational lensing, and the observed radial distributions of cluster galaxies, intracluster gas (i.e., the conventional $\beta$ model) and baryon fractions in a consistent way. This model, when combined with the observed X-ray surface brightness profiles of clusters, predicts that the overall temperature of intracluster gas has a tendency of radial decline with a mean polytropic index $\gamma\approx1.2$. Finally, the so-called $\beta$ discrepancy finds a natural explanation for such a polytropic gas in our model.

Rich and massive clusters of galaxies at intermediate redshift are capable of magnifying and distorting the images of background galaxies. A comparison of different mass estimators among these clusters can provide useful information about the distribution and composition of cluster matter and their dynamical evolution. Using a hitherto largest sample of lensing clusters drawn from literature, we compare the gravitating masses of clusters derived from the strong/weak gravitational lensing phenomena, from the X-ray measurements based on the assumption of hydrostatic equilibrium, and from the conventional isothermal sphere model for the dark matter profile characterized by the velocity dispersion and core radius of galaxy distributions in clusters. While there is an excellent agreement between the weak lensing, X-ray and isothermal sphere model determined cluster masses, these methods are likely to underestimate the gravitating masses enclosed within the central cores of clusters by a factor of 2--4 as compared with the strong lensing results. Such a mass discrepancy has probably arisen from the inappropriate applications of the weak lensing technique and the hydrostatic equilibrium hypothesis to the central regions of clusters as well as an unreasonably large core radius for both luminous and dark matter profiles. Nevertheless, it is pointed out that these cluster mass estimators may be safely applied on scales greater than the core sizes. Namely, the overall clusters of galaxies at intermediate redshift can still be regarded as the dynamically relaxed systems, in which the velocity dispersion of galaxies and the temperature of X-ray emitting gas are good indicators of the underlying gravitational potentials of clusters.

We present nine color CCD intermediate-band spectrophotometry of a two square degree field centered on the old open cluster M67, from 3890$\rm \AA$ to nearly 1$\mu$. These observations are taken as a part of the BATC (Beijing-Arizona-Taipei-Connecticut) Color Survey of the Sky, for both scientific and calibration reasons. With these data we show that the BATC survey can reach its goal of obtaining spectrophotometry to a zero point accuracy of 0.01 mag, and down to V = 21 with 0.3 mag random error. We fit the color-magnitude diagrams (CMDs) with Worthey's theoretical models. The net result is the excellent fit of the 4.0 Gyr, [Fe/H] = $-0.10$ model to our data, including a good fit to the main sequence (MS) turn-off. Our data are consistent with a toy model with 50\% of the stars in M67 being binaries and a random distribution of binary mass-ratios, although other models with different mass-ratio distributions cannot be ruled out. The spatial distribution and mass function (MF) of stars in M67 show marked effects of dynamical evolution and evaporation of stars from the cluster. Blue stragglers and binary stars are the most condensed within the cluster, with degree of condensation depending on mass.We find M67 to have an elongated shape, oriented at an angle of $15^{\circ}$ relative to the galactic plane. Within its tidal radius, the observed MF of M67 between 1.2 $\rm M_\odot$ and $\rm 0.8 M_\odot$ has a Salpeter slope $\rm \eta = -1.93 \pm 0.66$. For stars of mass below 0.8 $\rm M_\odot$, $\rm \eta \sim 0$. It is plausible that the leveling-off of the MF at lower masses is a result of evaporation of lower mass stars in this mass range at a rate of one every $\sim 10^7$ years. If so, it is plausible that the IMF of M67 has the canonical field value of $\rm \eta = -2.0$.