Source author record

M. James Jee

M. James Jee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2026arXiv

MrMARTIAN: A Multi-resolution Mass Reconstruction Algorithm Combining Free-form and Analytic Components

We present ${\tt MrMARTIAN}$ (Multi-resolution MAximum-entropy Reconstruction Technique Integrating Analytic Node), a new hybrid strong lensing (SL) modeling algorithm. By incorporating physically motivated analytic nodes into the free-form method ${\tt MARS}$, ${\tt MrMARTIAN}$ enables stable and flexible mass reconstructions while mitigating oversmoothing in the inner mass profile. Its multi-resolution framework increases the degrees of freedom in regions with denser strong lensing constraints, thereby enhancing computational efficiency for a fixed number of free parameters. We evaluate the performance of ${\tt MrMARTIAN}$ using publicly available simulated SL data and find that it consistently outperforms ${\tt MARS}$ in recovering both mass and magnification. In particular, it delivers significantly more stable reconstructions when multiple images are sparsely distributed. Finally, we apply ${\tt MrMARTIAN}$ to the galaxy cluster MACS J0416.1-2403, incorporating two analytic nodes centered on the northeastern and southwestern BCGs. Our mass model, constrained by 412 multiple images, achieves an image-plane rms scatter of ~0".11, the smallest to date for this dataset.

preprint2022arXiv

Comparison of spatial distributions of Intracluster light and Dark Matter

In a galaxy cluster, the relative spatial distributions of dark matter, member galaxies, gas, and intracluster light (ICL) may connote their mutual interactions over the cluster evolution. However, it is a challenging problem to provide a quantitative measure for the shape matching between two multi-dimensional scalar distributions. We present a novel methodology, named the {\em Weighted Overlap Coefficient (WOC)}, to quantify the similarity of 2-dimensional spatial distributions. We compare the WOC with a standard method known as the Modified Hausdorff Distance (MHD). We find that our method is robust, and performs well even with the existence of multiple sub-structures. We apply our methodology to search for a visible component whose spatial distribution resembled with that of dark matter. If such a component could be found to trace the dark matter distribution with high fidelity for more relaxed galaxy clusters, then the similarity of the distributions could also be used as a dynamical stage estimator of the cluster. We apply the method to six galaxy clusters at different dynamical stages simulated within the GRT simulation, which is an N-body simulation using the galaxy replacement technique. Among the various components (stellar particles, galaxies, ICL), the ICL+ brightest cluster galaxy (BCG) component most faithfully trace the dark matter distribution. Among the sample galaxy clusters, the relaxed clusters show stronger similarity in the spatial distribution of the dark matter and ICL+BCG than the dynamically young clusters. While the MHD results show weaker trend with the dynamical stages.

preprint2022arXiv

GREAT3 results I: systematic errors in shear estimation and the impact of real galaxy morphology

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the Sérsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

preprint2022arXiv

MARS: A New Maximum Entropy-Regularized Strong Lensing Mass Reconstruction Method

Free-form strong-lensing (SL) mass reconstructions typically suffer from overfitting, which manifest itself as false-positive small-scale fluctuations. We present a new free-form MAximum-entropy ReconStruction (${\tt MARS}$) method without the assumption that light traces mass (LTM). The ${\tt MARS}$ algorithm enables us to achieve excellent convergence in source positions ~0.001$"$, minimize spurious small-scale fluctuations, and provide a quasi-unique solution independently of initial conditions. Our method is tested with the publicly available synthetic SL data ${\tt FF-SIMS}$ and the comparison with the truth shows that the reconstruction quality is on a par with those of the best-performing LTM methods published in the literature, which have been demonstrated to outperform the existing free-form methods. In terms of the radial mass profile reconstruction, we achieve $<1$% agreements with the truth for the regions constrained by the multiple images. Finally, we apply ${\tt MARS}$ to A1689 and find that the cluster mass in the SL regime is dominated by the primary halo centered on the brightest cluster galaxy and the weaker secondary halo also coincident with the bright cluster member ~160 kpc northeast. Within the SL field, the A1689 radial profile is well-described by a Navarro-Frenk-White (NFW) profile with $c_{200}=5.53\pm0.77$ and $r_s=538^{+90}_{-100}$ kpc and we find no evidence that A1689 is over-concentrated.

preprint2022arXiv

Snowmass2021: Vera C. Rubin Observatory as a Flagship Dark Matter Experiment

Establishing that Vera C. Rubin Observatory is a flagship dark matter experiment is an essential pathway toward understanding the physical nature of dark matter. In the past two decades, wide-field astronomical surveys and terrestrial laboratories have jointly created a phase transition in the ecosystem of dark matter models and probes. Going forward, any robust understanding of dark matter requires astronomical observations, which still provide the only empirical evidence for dark matter to date. We have a unique opportunity right now to create a dark matter experiment with Rubin Observatory Legacy Survey of Space and Time (LSST). This experiment will be a coordinated effort to perform dark matter research, and provide a large collaborative team of scientists with the necessary organizational and funding supports. This approach leverages existing investments in Rubin. Studies of dark matter with Rubin LSST will also guide the design of, and confirm the results from, other dark matter experiments. Supporting a collaborative team to carry out a dark matter experiment with Rubin LSST is the key to achieving the dark matter science goals that have already been identified as high priority by the high-energy physics and astronomy communities.

preprint2021arXiv

Discovery of a double radio relic in ZwCl1447.2+2619: A rare testbed for shock acceleration models with a peculiar surface brightness ratio

We report a discovery of a double radio relic in the cluster merger ZwCl1447.2+2619 ($z=0.376$) with uGMRT observations at $420\rm~MHz$ and $700\rm~MHz$. The linear sizes of the northern and southern relics are $\sim0.3~$Mpc and $\sim1.2~$Mpc, respectively, which is consistent with the theoretical expectation that a larger relic is produced in the less massive subcluster side. However, ZwCl1447.2+2619 is unlike other known double radio relic systems, where the larger relics are much more luminous by several factors. In this merger the higher surface brightness of the smaller northern relic makes its total radio luminosity comparable to that of the much larger southern relic. The surface brightness ratio $\sim0.1$ between the two radio relics differs significantly from the relation observed in other double radio relic systems. From our radio spectral analysis, we find that both relics signify similar weak shocks with Mach numbers of $2.9\pm0.8$ and $2.0\pm0.7$ for the northern and southern relics, respectively. Moreover, the northern relic is connected to a discrete radio source with an optical counterpart, which indicates the possible presence of cosmic ray injection and re-acceleration. Therefore, we propose that this atypical surface brightness ratio can be explained with the particle acceleration efficiency precipitously dropping in the weak shock regime and/or with re-acceleration of fossil cosmic rays. Our multi-wavelength analysis and numerical simulation suggest that ZwCl1447.2+2619 is a post-merger, which has experienced a near head-on collision $\sim0.7\rm~Gyr$ ago.

preprint2021arXiv

Multiwavelength Analysis of A1240, the Double Radio Relic Merging Galaxy Cluster Embedded in a ~80 Mpc-long Cosmic Filament

We present a multiwavelength study of the double radio relic cluster A1240 at z=0.195. Our Subaru-based weak lensing analysis detects three mass clumps forming a ~4 Mpc filamentary structure elongated in the north-south orientation. The northern ($M_{200}=2.61_{-0.60}^{+0.51}\times10^{14} M_{\odot}$) and middle ($M_{200}=1.09_{-0.43}^{+0.34}\times10^{14} M_{\odot}$) mass clumps separated by ~1.3 Mpc are associated with A1240 and co-located with the X-ray peaks and cluster galaxy overdensities revealed by Chandra and MMT/Hectospec observations, respectively. The southern mass clump ($M_{200}=1.78_{-0.55}^{+0.44}\times10^{14} M_{\odot}$), ~1.5 Mpc to the south of the middle clump, coincides with the galaxy overdensity in A1237, the A1240 companion cluster at z=0.194. Considering the positions, orientations, and polarization fractions of the double radio relics measured by the LOFAR study, we suggest that A1240 is a post-merger binary system in the returning phase with the time-since-collision ~1.7 Gyr. With the SDSS DR16 data analysis, we also find that A1240 is embedded in the much larger-scale (~80 Mpc) filamentary structure whose orientation is in remarkable agreement with the hypothesized merger axis of A1240.

preprint2020arXiv

Baryonic feedback measurement from KV450 cosmic shear analysis

While baryonic feedback is one of the most important astrophysical systematics that we need to address in order to achieve precision cosmology, few weak lensing studies have directly measured its impact on the matter power spectrum. We report measurement of the baryonic feedback parameter with the constraints on its lower and upper limits from cosmic shear. We use the public data from the Kilo-Degree Survey and the VISTA Kilo-Degree Infrared Galaxy Survey spanning 450 deg$^2$. Estimating both cosmological and feedback parameters simultaneously, we obtain $A_{\rm b}=1.01_{-0.85}^{+0.80}$, which shows a consistency with the dark matter-only (DMO) case at the $\sim1.2~σ$ level and a tendency toward positive feedback; the $A_{\rm b}=0$ ($0.81$) value corresponds to the DMO (OWLS AGN) case. Despite this full constraint of the feedback parameter, our $S_8~(\equiv σ_8 \sqrt{Ω_m / 0.3})$ measurement ($0.739^{+0.036}_{- 0.035}$) shifts by only $\sim6$% of the statistical error, compared to the previous measurement. When we assume the flat $Λ$CDM cosmology favored by the Nine-Year Wilkinson Microwave Anisotropy Probe (Planck) result, the feedback parameter is constrained to be $A_{\rm b}=1.21_{-0.54}^{+0.61}$ ($1.60_{-0.52}^{+0.53}$), which excludes the DMO case at the $\sim2.2~σ$ ($\sim3.1~σ$) level.

preprint2020arXiv

Constraining the Mass of the Emerging Galaxy Cluster SpARCS1049+56 at z=1.71 with Infrared Weak Lensing

In the hierarchical structure formation model of the universe, galaxy clusters are assembled through a series of mergers. Accordingly, it is expected that galaxy clusters in the early universe are actively forming and dynamically young. Located at a high redshift of z=1.71, SpARCS1049+56 offers a unique look into the galaxy cluster formation process. This cluster has been shown to be rich in cluster galaxies and to have intense star formation. Its high redshift pushes a weak-lensing analysis beyond the regime of the optical spectrum into that of the infrared. Equipped with deep Hubble Space Telescope Wide Field Camera 3 UVIS and IR observations, we present a weak-lensing characterization of SpARCS1049+56. As few IR weak-lensing studies have been performed, we discuss the details of PSF modeling and galaxy shape measurement for an IR weak-lensing procedure and the systematics that come with the territory. It will be critical to understand these systematics in future weak-lensing studies in the IR with the next generation space telescopes such as JWST, Euclid, and WFIRST. Through a careful analysis, the mass distribution of this young galaxy cluster is mapped and the convergence peak is detected at a 3.3 sigma level. The weak-lensing mass of the cluster is estimated to be $3.5\pm1.2\times10^{14}\ \text{M}_\odot$ and is consistent with the mass derived from a mass-richness scaling relation. This mass is extreme for a cluster at such a high redshift and suggests that SpARCS1049+56 is rare in the standard $Λ$CDM universe.

preprint2020arXiv

Discovery of a Radio Relic in the Massive Merging Cluster SPT-CL 2023-5535 from the ASKAP-EMU PILOT SURVEY

The ASKAP-EMU survey is a deep wide-field radio continuum survey designed to cover the entire southern sky and a significant fraction of the northern sky up to $+30^{\circ}$. Here, we report a discovery of a radio relic in the merging cluster SPT-CL 2023-5535 at z=0.23 from the ASKAP-EMU pilot 300 sq. deg survey (800-1088 MHz). The deep high-resolution data reveal a $\sim2$ Mpc-scale radio halo elongated in the east-west direction, coincident with the intracluster gas. The radio relic is located at the western edge of this radio halo stretched $\sim0.5$ Mpc in the north-south orientation. The integrated spectral index of the radio relic within the narrow bandwidth is $α^{\scriptstyle \rm 1088~MHz}_{\scriptstyle \rm 800~MHz}=-0.76 \pm 0.06$. Our weak-lensing analysis shows that the system is massive ($M_{200}=1.04\pm0.36\times 10^{15} M_{\odot}$) and composed of at least three subclusters. We suggest a scenario, wherein the radio features arise from the collision between the eastern and middle subclusters. Our discovery illustrates the effectiveness of the ASKAP-EMU survey in detecting diffuse emissions in galaxy clusters and when completed, the survey will greatly increase the number of merging cluster detections with diffuse radio emissions.

preprint2020arXiv

Reconstruction of Radio Relics and X-ray Tails in an Off-axis Cluster Merger: Hydrodynamical Simulations of A115

Although a merging galaxy cluster is a useful laboratory to study many important astrophysical processes and fundamental physics, only limited interpretations are possible without careful analysis of the merger scenario. However, the study is demanding because a thorough comparison of multi-wavelength observations with detailed numerical simulations is required. In this paper, we present such a study for the off-axis binary merger A115. The system possesses a number of remarkable observational features, but no convincing merger scenario, explaining the shape and location of the radio relic in harmony with the orientation of the cometary X-ray tails, has been presented. Our hydrodynamical simulation, with adaptive mesh refinement, suggests that the cometary X-ray tail of A115 might be a slingshot tail and can arise $\sim0.3$ Gyrs after the impact and before the two subclusters reach their apocenters. This scenario can predict the location and orientation of the giant radio relic, which is parallel to the northern X-ray tail. In addition, our study indicates that diffusive shock acceleration alone cannot generate the observed radio power unless aided by stronger magnetic fields and/or more significant presence of fossil electrons.

preprint2019arXiv

Candidate Brown-dwarf Microlensing Events with Very Short Timescales and Small Angular Einstein Radii

Short-timescale microlensing events are likely to be produced by substellar brown dwarfs (BDs), but it is difficult to securely identify BD lenses based on only event timescales $t_{\rm E}$ because short-timescale events can also be produced by stellar lenses with high relative lens-source proper motions. In this paper, we report three strong candidate BD-lens events found from the search for lensing events not only with short timescales ($t_{\rm E} \lesssim 6~{\rm days}$) but also with very small angular Einstein radii ($θ_{\rm E}\lesssim 0.05~{\rm mas}$) among the events that have been found in the 2016--2019 observing seasons. These events include MOA-2017-BLG-147, MOA-2017-BLG-241, and MOA-2019-BLG-256, in which the first two events are produced by single lenses and the last event is produced by a binary lens. From the Bayesian analysis conducted with the combined $t_{\rm E}$ and $θ_{\rm E}$ constraint, it is estimated that the lens masses of the individual events are $0.051^{+0.100}_{-0.027}~M_\odot$, $0.044^{+0.090}_{-0.023}~M_\odot$, and $0.046^{+0.067}_{-0.023}~M_\odot/0.038^{+0.056}_{-0.019}~M_\odot$ and the probability of the lens mass smaller than the lower limit of stars is $\sim 80\%$ for all events. We point out that routine lens mass measurements of short time-scale lensing events require survey-mode space-based observations.

preprint2019arXiv

Chandra Observations of the Spectacular A3411-12 Merger Event

We present deep Chandra observations of A3411-12, a remarkable merging cluster that hosts the most compelling evidence for electron re-acceleration at cluster shocks to date. Using the $Y_X-M$ scaling relation, we find $r_{500} \sim 1.3$ Mpc, $M_{500} = (7.1 \pm 0.7) \times 10^{14} \ M_{\rm{\odot}}$, $kT=6.5\pm 0.1$ keV, and a gas mass of $M_{\rm g,500} = (9.7 \pm 0.1) \times 10^{13} M_\odot$. The gas mass fraction within $r_{500}$ is $f_{\rm g} = 0.14 \pm 0.01$. We compute the shock strength using density jumps to conclude that the Mach number of the merging subcluster is small ($M \leq 1.15_{-0.09}^{+0.14}$). We also present pseudo-density, projected temperature, pseudo-pressure, and pseudo-entropy maps. Based on the pseudo-entropy map we conclude that the cluster is undergoing a mild merger, consistent with the small Mach number. On the other hand, radio relics extend over Mpc scale in the A3411-12 system, which strongly suggests that a population of energetic electrons already existed over extended regions of the cluster.

preprint2019arXiv

OGLE-2016-BLG-1227L: A Wide-separation Planet from a Very Short-timescale Microlensing Event

We present the analysis of the microlensing event OGLE-2016-BLG-1227. The light curve of this short-duration event appears to be a single-lens event affected by severe finite-source effects. Analysis of the light curve based on single-lens single-source (1L1S) modeling yields very small values of the event timescale, $t_{\rm E}\sim 3.5$ days, and the angular Einstein radius, $θ_{\rm E}\sim 0.009$ mas, making the lens a candidate of a free-floating planet. Close inspection reveals that the 1L1S solution leaves small residuals with amplitude $ΔI\lesssim 0.03$ mag. We find that the residuals are explained by the existence of an additional widely-separated heavier lens component, indicating that the lens is a wide-separation planetary system rather than a free-floating planet. From Bayesian analysis, it is estimated that the planet has a mass of $M_{\rm p} = 0.79^{+1.30}_{-0.39} M_{\rm J}$ and it is orbiting a low-mass host star with a mass of $M_{\rm host}=0.10^{+0.17}_{-0.05} M_\odot$ located with a projected separation of $a_\perp=3.4^{+2.1}_{-1.0}$ au. The planetary system is located in the Galactic bulge with a line-of-sight separation from the source star of $D_{\rm LS}=1.21^{+0.96}_{-0.63}$ kpc. The event shows that there are a range of deviations in the signatures of host stars for apparently isolated planetary lensing events and that it is possible to identify a host even when a deviation is subtle.

preprint2019arXiv

OGLE-2018-BLG-1700L: Microlensing Planet in Binary Stellar System

We report the discovery of a planet in a binary that was discovered from the analysis of the microlensing event OGLE-2018-BLG-1700. We identify the triple nature of the lens from the fact that the complex anomaly pattern can be decomposed into two parts produced by two binary-lens events, in which one binary pair has a very low mass ratio of $\sim 0.01$ between the lens components and the other pair has a mass ratio of $\sim 0.3$. We find two sets of degenerate solutions, in which one solution has a projected separation between the primary and its stellar companion less than the angular Einstein radius $\thetae$ (close solution), while the other solution has a separation greater than $\thetae$ (wide solution). From the Bayesian analysis with the constraints of the event time scale and angular Einstein radius together with the location of the source lying in the far disk behind the bulge, we find that the planet is a super-Jupiter with a mass of $4.4^{+3.0}_{-2.0}~M_{\rm J}$ and the stellar binary components are early and late M-type dwarfs with masses $0.42^{+0.29}_{-0.19}~M_\odot$ and $0.12^{+0.08}_{-0.05}~M_\odot$, respectively, and the planetary system is located at a distance of $D_{\rm L}=7.6^{+1.2}_{-0.9}~{\rm kpc}$. The planet is a circumstellar planet according to the wide solution, while it is a circumbinary planet according to the close solution. The projected primary-planet separation is $2.8^{+3.2}_{-2.5}~{\rm au}$ commonly for the close and wide solutions, but the primary-secondary binary separation of the close solution, $0.75^{+0.87}_{-0.66}~{\rm au}$, is widely different from the separation, $10.5^{+12.1}_{-9.2}~{\rm au}$, of the wide solution.

preprint2019arXiv

Precise Mass Determination of SPT-CL J2106-5844, the Most Massive Cluster at z>1

We present a detailed high-resolution weak-lensing (WL) study of SPT-CL J2106-5844 at z=1.132, claimed to be the most massive system discovered at z > 1 in the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey. Based on the deep imaging data from the Advanced Camera for Surveys and Wide Field Camera 3 on-board the Hubble Space Telescope, we find that the cluster mass distribution is asymmetric, composed of a main clump and a subclump ~640 kpc west thereof. The central clump is further resolved into two smaller northwestern and southeastern substructures separated by ~150 kpc. We show that this rather complex mass distribution is more consistent with the cluster galaxy distribution than a unimodal distribution as previously presented. The northwestern substructure coincides with the BCG and X-ray peak while the southeastern one agrees with the location of the number density peak. These morphological features and the comparison with the X-ray emission suggest that the cluster might be a merging system. We estimate the virial mass of the cluster to be $M_{200c} = (10.4^{+3.3}_{-3.0}\pm1.0)~\times~10^{14}~M_{\odot}$, where the second error bar is the systematic uncertainty. Our result confirms that the cluster SPT-CL J2106-5844 is indeed the most massive cluster at z>1 known to date. We demonstrate the robustness of this mass estimate by performing a number of tests with different assumptions on the centroids, mass-concentration relations, and sample variance.

preprint2016arXiv

Cosmic Shear Results from the Deep Lens Survey - II: Full Cosmological Parameter Constraints from Tomography

We present a tomographic cosmic shear study from the Deep Lens Survey (DLS), which, providing a limiting magnitude r_{lim}~27 (5 sigma), is designed as a pre-cursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on depth. Using five tomographic redshift bins, we study their auto- and cross-correlations to constrain cosmological parameters. We use a luminosity-dependent nonlinear model to account for the astrophysical systematics originating from intrinsic alignments of galaxy shapes. We find that the cosmological leverage of the DLS is among the highest among existing >10 sq. deg cosmic shear surveys. Combining the DLS tomography with the 9-year results of the Wilkinson Microwave Anisotropy Probe (WMAP9) gives Omega_m=0.293_{-0.014}^{+0.012}, sigma_8=0.833_{-0.018}^{+0.011}, H_0=68.6_{-1.2}^{+1.4} km/s/Mpc, and Omega_b=0.0475+-0.0012 for LCDM, reducing the uncertainties of the WMAP9-only constraints by ~50%. When we do not assume flatness for LCDM, we obtain the curvature constraint Omega_k=-0.010_{-0.015}^{+0.013} from the DLS+WMAP9 combination, which however is not well constrained when WMAP9 is used alone. The dark energy equation of state parameter w is tightly constrained when Baryonic Acoustic Oscillation (BAO) data are added, yielding w=-1.02_{-0.09}^{+0.10} with the DLS+WMAP9+BAO joint probe. The addition of supernova constraints further tightens the parameter to w=-1.03+-0.03. Our joint constraints are fully consistent with the final Planck results and also the predictions of a LCDM universe.

preprint2016arXiv

IDCS J1426.5+3508: Weak Lensing Analysis of a Massive Galaxy Cluster at $z=1.75$

We present a weak lensing study of the galaxy cluster IDCS J1426.5+3508 at $z=1.75$, which is the highest redshift strong lensing cluster known and the most distant cluster for which a weak lensing analysis has been undertaken. Using F160W, F814W, and F606W observations with the Hubble Space Telescope, we detect tangential shear at $2σ$ significance. Fitting a Navarro-Frenk-White mass profile to the shear with a theoretical median mass-concentration relation, we derive a mass $M_{200,\mathrm{crit}}=2.3^{+2.1}_{-1.4}\times10^{14}$ M$_{\odot}$. This mass is consistent with previous mass estimates from the Sunyaev-Zel'dovich (SZ) effect, X-ray, and strong lensing. The cluster lies on the local SZ-weak lensing mass scaling relation observed at low redshift, indicative of minimal evolution in this relation.

preprint2016arXiv

MC^2: Mapping the Dark Matter Distribution of the "Toothbrush" Cluster RX J0603.3+4214 with Hubble Space Telescope and Subaru Weak-lensing

The galaxy cluster RX J0603.3+4214 at z=0.225 is one of the rarest clusters boasting an extremely large (~2 Mpc) radio-relic. Because of the remarkable morphology of the relic, the cluster is nicknamed "Toothbrush Cluster". Although the cluster's underlying mass distribution is one of the critical pieces of information needed to reconstruct the merger scenario responsible for the puzzling radio-relic morphology, its proximity to the Galactic plane b~10 deg has imposed significant observational challenges. We present a high-resolution weak-lensing study of the cluster with Subaru/Suprime Cam and Hubble Space Telescope imaging data. Our mass reconstruction reveals that the cluster is comprised of complicated dark matter substructures closely tracing the galaxy distribution, however in contrast with the relatively simple binary X-ray morphology. Nevertheless, we find that the cluster mass is still dominated by the two most massive clumps aligned north-south with a ~3:1 mass ratio (M_{200}=6.29_{-1.62}^{+2.24} x 10^{14} Msun and 1.98_{-0.74}^{+1.24} x 10^{14} Msun for the northern and southern clumps, respectively). The southern mass peak is ~2' offset toward the south with respect to the corresponding X-ray peak, which has a "bullet"-like morphology pointing south. Comparison of the current weak-lensing result with the X-ray, galaxy, and radio-relic suggests that perhaps the dominant mechanism responsible for the observed relic may be a high-speed collision of the two most massive subclusters, although the peculiarity of the morphology necessitates involvement of additional sub-clusters. Careful numerical simulations should follow in order to obtain more complete understanding of the merger scenario utilizing all existing observations.

preprint2016arXiv

The Globular Cluster System of the Coma cD Galaxy NGC 4874 from Hubble Space Telescope ACS and WFC3/IR Imaging

We present new HST optical and near-infrared (NIR) photometry of the rich globular cluster (GC) system of NGC 4874, the cD galaxy in the core of the Coma cluster (Abell 1656). NGC 4874 was observed with the HST Advanced Camera for Surveys in the F475W (g) and F814W (I) passbands and the Wide Field Camera 3 IR Channel in F160W (H). The GCs in this field exhibit a bimodal optical color distribution with more than half of the GCs falling on the red side at g-I > 1. Bimodality is also present, though less conspicuously, in the optical-NIR I-H color. Consistent with past work, we find evidence for nonlinearity in the g-I versus I-H color-color relation. Our results thus underscore the need for understanding the detailed form of the color-metallicity relations in interpreting observational data on GC bimodality. We also find a very strong color-magnitude trend, or "blue tilt," for the blue component of the optical color distribution of the NGC 4874 GC system. A similarly strong trend is present for the overall mean I-H color as a function of magnitude; for M_814 < -10 mag, these trends imply a steep mass-metallicity scaling with $Z\propto M_{\rm GC}^{1.4\pm0.4}$, but the scaling is not a simple power law and becomes much weaker at lower masses. As in other similar systems, the spatial distribution of the blue GCs is more extended than that of the red GCs, partly because of blue GCs associated with surrounding cluster galaxies. In addition, the center of the GC system is displaced by 4+/-1 kpc towards the southwest from the luminosity center of NGC 4874, in the direction of NGC 4872. Finally, we remark on a dwarf elliptical galaxy with a noticeably asymmetrical GC distribution. Interestingly, this dwarf has a velocity of nearly -3000 km/s with respect to NGC 4874; we suggest it is on its first infall into the cluster core and is undergoing stripping of its GC system by the cluster potential.

preprint2015arXiv

The return of the merging galaxy subclusters of El Gordo?

Merging galaxy clusters with radio relics provide rare insights to the merger dynamics as the relics are created by the violent merger process. We demonstrate one of the first uses of the properties of the radio relic to reduce the uncertainties of the dynamical variables and determine the 3D configuration of a cluster merger, ACT-CL J0102-4915, nicknamed El Gordo. From the double radio relic observation and the X-ray observation of a comet-like gas morphology induced by motion of the cool core, it is widely believed that El Gordo is observed shortly after the first core-passage of the subclusters. We employ a Monte Carlo simulation to investigate the three-dimensional (3D) configuration and dynamics of El Gordo. Using the polarization fraction of the radio relic, we constrain the estimate of the angle between the plane of the sky and the merger axis to be $α= 21~{\rm degree} \pm^9_{11}$. We find the relative 3D merger speed of El Gordo to be $2400\pm^{400}_{200}~{\rm km}~{\rm s}^{-1}$ at pericenter. The two possible estimates of the time-since-pericenter are $0.46\pm^{0.09}_{0.16}$ Gyr and $0.91\pm^{0.22}_{0.39}$ Gyr for the outgoing and returning scenario respectively. We put our estimates of the time-since-pericenter into context by showing that if the time-averaged shock velocity is approximately equal to or smaller than the pericenter velocity of the corresponding subcluster in the center of mass frame, the two subclusters are more likely to be moving towards, rather than away, from each other, post apocenter. We compare and contrast the merger scenario of El Gordo with that of the Bullet Cluster, and show that this late-stage merging scenario explains why the southeast dark matter lensing peak of El Gordo is closer to the merger center than the southeast cool core.

preprint2014arXiv

Hubble Space Telescope/Advanced Camera for Surveys Confirmation of the Dark Substructure in A520

We present the results from a weak gravitational lensing study of the merging cluster A520 based on the analysis of Hubble Space Telescope/Advanced Camera for Surveys (ACS) data. The excellent data quality allows us to reach a mean number density of source galaxies of ~109 per sq. arcmin, which improves both resolution and significance of the mass reconstruction compared to a previous study based on Wide Field Planetary Camera 2 (WFPC2) images. We take care in removing instrumental effects such as the trailing of charge due to radiation damage of the ACS detector and the position-dependent point spread function (PSF). This new ACS analysis confirms the previous claims that a substantial amount of dark mass is present between two luminous subclusters. We examine the distribution of cluster galaxies and observe very little light at this location. We find that the centroid of the dark peak in the current ACS analysis is offset to the southwest by ~1 arcmin with respect to the centroid from the WFPC2 analysis. Interestingly, this new centroid is in better spatial agreement with the location where the X-ray emission is strongest, and the mass-to-light ratio estimated with this centroid is much higher 813+-78 M_sun/L_Rsun than the previous value; the aperture mass based on the WFPC2 centroid provides a slightly lower, but consistent mass. Although we cannot provide a definite explanation for the presence of the dark peak, we discuss a revised scenario, wherein dark matter with a more conventional range sigma_DM/m_DM < 1 cm^2/g of self-interacting cross-section can lead to the detection of this dark substructure. If supported by detailed numerical simulations, this hypothesis opens up the possibility that the A520 system can be used to establish a lower limit of the self-interacting cross-section of dark matter.

preprint2014arXiv

MC^2: Constraining the Dark Matter Distribution of the Violent Merging Galaxy Cluster CIZA J2242.8+5301 by Piercing through the Milky Way

CIZA J2242.8+5301 is a merging system with a prominent (~2 Mpc long) radio relic, which together with the morphology of the X-ray emission provides strong evidence for a violent collision along the N-S axis. We present our constraints on the dark matter distribution of this unusual system using Subaru and CFHT imaging data. Measuring a high S/N lensing signal from this cluster is potentially a challenging task because of its proximity to the Milky Way plane (|b|~5 deg). We overcome this challenge with careful observation planning and systematics control, which enables us to successfully map the dark matter distribution of the cluster with high fidelity. The resulting mass map shows that the mass distribution is highly elongated along the N-S merger axis inferred from the orientation of the radio relics. Based on our mass reconstruction, we identify two sub-clusters, which coincide with the cluster galaxy distributions. We determine their masses using MCMC analysis by simultaneously fitting two NFW halos without fixing their centroids. The resulting masses of the northern and southern systems are $M_{200}=11.0_{-3.2}^{+3.7}\times10^{14} M_{\sun}$ and $9.8_{-2.5}^{+3.8}\times10^{14} M_{\sun}$, respectively, indicating that we are witnessing a post-collision of two giant systems of nearly equal mass. When the mass and galaxy centroids are compared in detail, we detect ~ 1' (~190 kpc) offsets in both northern and southern sub-clusters. We find that the galaxy luminosity-mass offset for the northern clump is statistically significant at the ~2 sigma level whereas the detection is only marginal for the southern sub-cluster in part because of a relatively large mass centroid error. We conclude that it is yet premature to uniquely attribute the galaxy-mass misalignment to SIDM and discuss caveats.

preprint2014arXiv

MC^2: Galaxy Imaging and Redshift Analysis of the Merging Cluster CIZA J2242.8+5301

X-ray and radio observations of CIZA J2242.8+5301 suggest that it is a major cluster merger. Despite being well studied in the X-ray, and radio, little has been presented on the cluster structure and dynamics inferred from its galaxy population. We carried out a deep (i<25) broad band imaging survey of the system with Subaru SuprimeCam (g & i bands) and the Canada France Hawaii Telescope (r band) as well as a comprehensive spectroscopic survey of the cluster area (505 redshifts) using Keck DEIMOS. We use this data to perform a comprehensive galaxy/redshift analysis of the system, which is the first step to a proper understanding the geometry and dynamics of the merger, as well as using the merger to constrain self-interacting dark matter. We find that the system is dominated by two subclusters of comparable richness with a projected separation of 6.9'^{+0.7}_{-0.5} (1.3^{+0.13}_{-0.10} Mpc). We find that the north and south subclusters have similar redshifts of z=0.188 with a relative line-of-sight velocity difference of 69+/-190 km/s. We also find that north and south subclusters have velocity dispersions of 1160^{+100}_{-90} km/s and 1080^{+100}_{-70} km/s, respectively. These correspond to masses of 16.1^{+4.6}_{-3.3}x10^14 M_sun and 13.0^{+4.0}_{-2.5}x10^14 M_sun, respectively. While velocity dispersion measurements of merging clusters can be biased we believe the bias in this system to be minor due to the large projected separation and nearly plane-of-sky merger configuration. CIZA J2242.8+5301 is a relatively clean dissociative cluster merger with near 1:1 mass ratio, which makes it an ideal merger for studying merger associated physical phenomena.

preprint2014arXiv

The rise and fall of star-formation in $\bf z\sim0.2$ merging galaxy clusters

CIZA J2242.8+5301 (`Sausage') and 1RXS J0603.3+4213 (`Toothbrush') are two low-redshift ($z\sim0.2$), massive ($\sim2\times10^{15}M_\odot$), post-core passage merging clusters, which host shock waves traced by diffuse radio emission. To study their star-formation properties, we uniformly survey the `Sausage' and `Toothbrush' clusters in broad and narrow band filters and select a sample of $201$ and $463$ line emitters, down to a rest-frame equivalent width ($13$Å). We robustly separate between H$α$ and higher redshift emitters using a combination of optical multi-band (B, g, V, r, i, z) and spectroscopic data. We build H$α$ luminosity functions for the entire cluster region, near the shock fronts, and away from the shock fronts and find striking differences between the two clusters. In the dynamically younger, $1$ Gyr old `Sausage' cluster we find numerous ($59$) H$α$ emitters above a star-formation rate (SFR) of $0.17$ M$_{\sun}$ yr$^{-1}$ surprisingly located in close proximity to the shock fronts, embedded in very hot intra-cluster medium plasma. The SFR density for the cluster population is at least at the level of typical galaxies at $z\sim2$. Down to the same star-formation rate, the possibly dynamically more evolved `Toothbrush' cluster has only $9$ H$α$ galaxies. The cluster H$α$ galaxies fall on the SFR-stellar mass relation $z\sim0.2$ for the field. However, the `Sausage' cluster has an H$α$ emitter density $>20$ times that of blank fields. If the shock passes through gas-rich cluster galaxies, the compressed gas could collapse into dense clouds and excite star-formation for a few $100$ Myr. This process ultimately leads to a rapid consumption of the molecular gas, accelerating the transformation of gas-rich field spirals into cluster S0s or ellipticals.

preprint2014arXiv

Weighing "El Gordo" with a Precision Scale: Hubble Space Telescope Weak-lensing Analysis of the Merging Galaxy Cluster ACT-CL J0102-4915 at z=0.87

(Abridged) We present a HST weak-lensing study of the merging galaxy cluster "El Gordo" (ACT-CL J0102-4915) at z=0.87 discovered by the Atacama Cosmology Telescope collaboration as the strongest SZ decrement in its ~1000 sq. deg survey. Our weak-lensing analysis confirms that ACT-CL J0102-4915 is indeed an extreme system consisting of two massive (~10^15 Msun each) subclusters with a projected separation of ~0.7 Mpc. This binary mass structure revealed by our lensing study is consistent with the cluster galaxy distribution and the dynamical study carried out with 89 spectroscopic members. We estimate the mass of ACT-CL J0102-4915 by simultaneously fitting two axisymmetric NFW profiles allowing their centers to vary. Our MCMC analysis shows that the masses of the northwestern (NW) and the southeastern (SE) components are M200c=(1.38+-0.22) x 10^15 Msun and (0.78+-0.20) x 10^15 Msun, respectively. The lensing-based velocity dispersions are consistent with their spectroscopic measurements. The centroids of both components are tightly constrained (~4") and close to the optical luminosity centers. The X-ray and mass peaks are spatially offset by ~8" (~62 kpc), which is significant at the ~2 sigma confidence level and confirms that the baryonic and dark matter in this cluster are disassociated. The dark matter peak, however, does not lead the gas peak in the direction expected if we are viewing the cluster soon after first core passage during a high speed merger. Under the assumption that the merger is happening in the plane of the sky, extrapolation of the two NFW halos to a radius r200a=2.4 Mpc yields a combined mass of M200a=(3.13+-0.56) x 10^15 Msun. This extrapolated total mass is consistent with our two-component-based dynamical analysis and previous X-ray measurements, projecting ACT-CL J0102-4915 to be the most massive cluster at z>0.6 known to date.

preprint2013arXiv

Cosmic shear results from the deep lens survey - I: Joint constraints on omega_m and sigma_8 with a two-dimensional analysis

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 4-meter telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the Omega_M-sigma_8 likelihood contour tightly, providing a joint constraint of Omega_M=0.262+-0.051 and sigma_8=0.868+-0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the Omega_M-sigma_8 degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain Omega_M=0.278+-0.018 and sigma_8=0.815+-0.020.

preprint2013arXiv

Morphology with Light Profile Fitting of Confirmed Cluster Galaxies at z=0.84

We perform a morphological study of 124 spectroscopically confirmed cluster galaxies in the z=0.84 galaxy cluster RX J0152.7-1357. Our classification scheme includes color information, visual morphology, and 1-component and 2-component light profile fitting derived from Hubble Space Telescope riz imaging. We adopt a modified version of a detailed classification scheme previously used in studies of field galaxies and found to be correlated with kinematic features of those galaxies. We compare our cluster galaxy morphologies to those of field galaxies at similar redshift. We also compare galaxy morphologies in regions of the cluster with different dark-matter density as determined by weak-lensing maps. We find an early-type fraction for the cluster population as a whole of 47%, about 2.8 times higher than the field, and similar to the dynamically young cluster MS 1054 at similar redshift. We find the most drastic change in morphology distribution between the low and intermediate dark matter density regions within the cluster, with the early type fraction doubling and the peculiar fraction dropping by nearly half. The peculiar fraction drops more drastically than the spiral fraction going from the outskirts to the intermediate-density regions. This suggests that many galaxies falling into clusters at z~0.8 may evolve directly from peculiar, merging, and compact systems into early-type galaxies, without having the chance to first evolve into a regular spiral galaxy.

preprint2012arXiv

Galaxy-Mass Correlations on 10 Mpc Scales in the Deep Lens Survey

We examine the projected correlation of galaxies with mass from small scales (<few hundred kpc) where individual dark matter halos dominate, out to 15 Mpc where correlated large-scale structure dominates. We investigate these profiles as a function of galaxy luminosity and redshift. Selecting 0.8 million galaxies in the Deep Lens Survey, we use photometric redshifts and stacked weak gravitational lensing shear tomography out to radial scales of 1 degree from the centers of foreground galaxies. We detect correlated mass density from multiple halos and large-scale structure at radii larger than the virial radius, and find the first observational evidence for growth in the galaxy-mass correlation on 10 Mpc scales with decreasing redshift and fixed range of luminosity. For a fixed range of redshift, we find a scaling of projected halo mass with rest-frame luminosity similar to previous studies at lower redshift. We control systematic errors in shape measurement and photometric redshift, enforce volume completeness through absolute magnitude cuts, and explore residual sample selection effects via simulations.

preprint2011arXiv

Towards Precision LSST Weak-Lensing Measurement - I: Impacts of Atmospheric Turbulence and Optical Aberration

The weak-lensing science of the LSST project drives the need to carefully model and separate the instrumental artifacts from the intrinsic lensing signal. The dominant source of the systematics for all ground based telescopes is the spatial correlation of the PSF modulated by both atmospheric turbulence and optical aberrations. In this paper, we present a full FOV simulation of the LSST images by modeling both the atmosphere and the telescope optics with the most current data for the telescope specifications and the environment. To simulate the effects of atmospheric turbulence, we generated six-layer phase screens with the parameters estimated from the on-site measurements. For the optics, we combined the ray-tracing tool ZEMAX and our simulated focal plane data to introduce realistic aberrations and focal plane height fluctuations. Although this expected flatness deviation for LSST is small compared with that of other existing cameras, the fast f-ratio of the LSST optics makes this focal plane flatness variation and the resulting PSF discontinuities across the CCD boundaries significant challenges in our removal of the systematics. We resolve this complication by performing PCA CCD-by-CCD, and interpolating the basis functions using conventional polynomials. We demonstrate that this PSF correction scheme reduces the residual PSF ellipticity correlation below 10^-7 over the cosmologically interesting scale. From a null test using HST/UDF galaxy images without input shear, we verify that the amplitude of the galaxy ellipticity correlation function, after the PSF correction, is consistent with the shot noise set by the finite number of objects. Therefore, we conclude that the current optical design and specification for the accuracy in the focal plane assembly are sufficient to enable the control of the PSF systematics required for weak-lensing science with the LSST.

preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.