BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shude Mao

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.EPastro-ph.SRastro-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

99works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

99 published item(s)

preprint2026arXiv

Searching for Galaxy Cluster-Scale Strong lenses from the DESI Legacy Imaging Surveys

Galaxy cluster-scale strong gravitational lensing systems are rare yet valuable tools for investigating the properties of dark matter and dark energy, as well as providing the opportunity to study the distant universe at flux levels and spatial resolutions that would otherwise be unavailable. Large-scale imaging surveys present unprecedented opportunities to expand the sample of cluster lenses. In this study, we adopt a deep learning-based approach to identify cluster lenses from the DESI Legacy Imaging Surveys, utilizing the catalog of galaxy cluster candidates identified by Zou et al. (2021). Our lens-finder employs a ResNet-18 architecture, trained with mock images of cluster lenses as positives and observational images of cluster scale non-lenses as negatives. We do an iterative operation to increase the completeness of our work, namely adding the found true positive samples back to the training set and training again for several times. Human inspection is conducted to further refine the candidates, categorizing them into grades (A, B, C) according to the significance of the strongly lensed arcs. Reviewing all 540,432 objects in Zou's catalog, we discover 485 high-confidence cluster lens candidates with a cluster $M_{500}$ range of $10^{13.67\sim14.97}M_{\odot}$ and a Brightest Central Galaxy (BCG) redshift range of $0.04\sim0.89$. After excluding the lens candidates listed in previous studies, we identify 247 newly discovered cluster lens candidates, including 16 grade A, 90 grade B, and 141 grade C. This catalog of cluster lens candidates is publicly available online, and follow-up observations are encouraged to confirm and conduct thorough investigations of these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Introduction to the Chinese Space Station Survey Telescope (CSST)

The Chinese Space Station Survey Telescope (CSST) is an upcoming Stage-IV sky survey telescope, distinguished by its large field of view (FoV), high image quality, and multi-band observation capabilities. It can simultaneously conduct precise measurements of the Universe by performing multi-color photometric imaging and slitless spectroscopic surveys. The CSST is equipped with five scientific instruments, i.e. Multi-band Imaging and Slitless Spectroscopy Survey Camera (SC), Multi-Channel Imager (MCI), Integral Field Spectrograph (IFS), Cool Planet Imaging Coronagraph (CPI-C), and THz Spectrometer (TS). Using these instruments, CSST is expected to make significant contributions and discoveries across various astronomical fields, including cosmology, galaxies and active galactic nuclei (AGN), the Milky Way and nearby galaxies, stars, exoplanets, Solar System objects, astrometry, and transients and variable sources. This review aims to provide a comprehensive overview of the CSST instruments, observational capabilities, data products, and scientific potential.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Revisiting the atmosphere of HAT-P-70b with CARMENES high-resolution transmission spectroscopy

Owing to hot and inflated envelopes that facilitate atmospheric studies, ultra-hot Jupiters (UHJs) have attracted much attention. Significant progress has been achieved, from enlarging the sample size to broadening the studies to encompass diverse stellar types and ages. Here, we present a transmission spectroscopy study of HAT-P-70b, an UHJ orbiting a young A-type star, through high-resolution observations with CARMENES at the 3.5m Calar Alto telescope. By using the line-by-line technique, we confirm the previous detections of Ha, Na I, and Ca II, report a new tentative detection of K I, and impose an upper limit on the He triplet absorption. Through cross-correlation analysis, we identify the Ca II and Fe I absorptions, both blue-shifted by approximately 5 km/s, indicating a day-to-night side atmospheric wind. Additionally, we find a new tentative detection of K I. We do not see any significant atmospheric molecular signal in the near-infrared data. Putting HAT-P-70b in the context of UHJs from the literature, it turns out that (1) Ha absorption is more common on gas giants orbiting stars younger than 1 Gyr, with a relative detection probability of $P_{\rm Age<1\,Gyr}({\rm Ha})/P_{\rm Age\geq1\,Gyr}({\rm Ha})\sim 3$; (2) any UHJ is likely to exhibit Fe I absorption if it has Ca II.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Early-type galaxy density profiles from IllustrisTNG: III. Effects on outer kinematic structure

Early-type galaxies (ETGs) possess total density profiles close to isothermal, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at $1.5 R_{\mathrm{eff}}$ (effective radius) that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on these kinematic features through mock ETGs in IllustrisTNG. We create mock Integral-Field-Unit observations to extract projected stellar kinematic features for 207 ETGs with stellar mass $M_{\ast}\geqslant 10^{11} \mathrm{M_{\odot}}$ in TNG100-1. The mock observations reproduce the key outer ($1.5R_{\mathrm{eff}}$) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope $γ_{\mathrm{outer}}$ and the kurtosis $h_{4}$'s gradient. We find that $h_{4}$ is uncorrelated with stellar orbital anisotropy beyond $R_{\mathrm{eff}}$; instead we find that the variations in $γ_{\mathrm{outer}}$ and outer $h_{4}$ (a good proxy for $h_{4}$ gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the origin of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates the two quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

KMT-2021-BLG-0171Lb and KMT-2021-BLG-1689Lb: Two Microlensing Planets in the KMTNet High-cadence Fields with Followup Observations

Follow-up observations of high-magnification gravitational microlensing events can fully exploit their intrinsic sensitivity to detect extrasolar planets, especially those with small mass ratios. To make followup more uniform and efficient, we develop a system, HighMagFinder, based on the real-time data from the Korean Microlensing Telescope Network (KMTNet) to automatically alert possible ongoing high-magnification events. We started a new phase of follow-up observations with the help of HighMagFinder in 2021. Here we report the discovery of two planets in high-magnification microlensing events, KMT-2021-BLG-0171 and KMT-2021-BLG-1689, which were identified by the HighMagFinder. We find that both events suffer the ``central-resonant'' caustic degeneracy. The planet-host mass-ratio is $q\sim4.7\times10^{-5}$ or $q\sim 2.2\times10^{-5}$ for KMT-2021-BLG-0171, and $q\sim2.5\times10^{-4}$ or $q\sim 1.8\times10^{-4}$ for KMT-2021-BLG-1689. Together with two events reported by Ryu et al. (2022), four cases that suffer such degeneracy have been discovered in the 2021 season alone, indicating that the degenerate solutions may have been missed in some previous studies. We also propose a new factor for weighting the probability of each solution from the phase-space. The resonant interpretations for the two events are disfavored under this consideration. This factor can be included in future statistical studies to weight degenerate solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

OGLE-2018-BLG-0799Lb: a $q \sim 2.7 \times 10^{-3}$ Planet with Spitzer Parallax

We report the discovery and analysis of a planet in the microlensing event OGLE-2018-BLG-0799. The planetary signal was observed by several ground-based telescopes, and the planet-host mass ratio is $q = (2.65 \pm 0.16) \times 10^{-3}$. The ground-based observations yield a constraint on the angular Einstein radius $θ_{\rm E}$, and the microlensing parallax vector $\vecπ_{\rm E}$, is strongly constrained by the Spitzer data. However, the 2019 Spitzer baseline data reveal systematics in the Spitzer photometry, so there is ambiguity in the magnitude of the parallax. In our preferred interpretation, a full Bayesian analysis using a Galactic model indicates that the planetary system is composed of an $M_{\rm planet} = 0.26_{-0.11}^{+0.22}~M_{J}$ planet orbiting an $M_{\rm host} = 0.093_{-0.038}^{+0.082}~M_{\odot}$, at a distance of $D_{\rm L} = 3.71_{-1.70}^{+3.24}$ kpc. An alternate interpretation of the data shifts the localization of the minima along the arc-shaped microlens parallax constraints. This, in turn, yields a more massive host with median mass of $0.13 {M_{\odot}}$ at a distance of 6.3 kpc. This analysis demonstrates the robustness of the osculating circles formalism, but shows that further investigation is needed to assess how systematics affect the specific localization of the microlens parallax vector and, consequently, the inferred physical parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

OGLE-2019-BLG-1470LABc: Another Microlensing Giant Planet in a Binary System?

We report the discovery and analysis of a candidate triple-lens single-source (3L1S) microlensing event, OGLE-2019-BLG-1470. This event was first classified as a normal binary-lens single-source (2L1S) event, but a careful 2L1S modelling showed that it needs an additional lens or source to fit the observed data. It is found that the 3L1S model provides the best fit, but the binary-lens binary-source (2L2S) model is only disfavoured by $Δχ^2 \simeq 18$. All of the feasible models include a planet with planet-to-host mass-ratios $10^{-3} \lesssim q \lesssim 10^{-2}$. A Bayesian analysis based on a Galactic model indicates that the planet is super-Jovian, and the projected host-planet separation is about 3 $\mathrm{au}$. Specifically, for the best-fit 3L1S model, the two stars have masses of $M_1=0.57^{+0.43}_{-0.32}M_{\odot}$, and $M_2=0.18^{+0.15}_{-0.10}M_{\odot}$, with projected separation of $1.3^{+0.5}_{-0.5}$ $\mathrm{au}$, and the planetary mass is $M_3=2.2^{+1.8}_{-1.3}M_{\rm{Jupiter}}$. For the 2L2S model, the masses of the host star and the planet are $0.55^{+0.44}_{-0.31}M_{\odot}$ and $4.6^{+3.7}_{-2.6}M_{\rm{Jupiter}}$, respectively. By investigating the properties of all known microlensing planets in binary systems, we find that all planets in binary systems published by the KMTNet survey are located inside the resonant caustics range with $q \gtrsim 2 \times 10^{-3}$, indicating the incompleteness of the KMTNet sample for planets in binary systems. Thus, planets in binary systems cannot be included in the current study of the KMTNet mass-ratio function, and a systematic search for planetary anomalies in KMTNet microlensing light curves of binary systems is needed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Sub-percentage measure of distances to redshift of 0.1 by a new cosmic ruler

Distance-redshift diagrams probe expansion history of the Universe. We show that the stellar mass-binding energy (massE) relation of galaxies proposed in our previous study offers a new distance ruler at cosmic scales. By using elliptical galaxies in the main galaxy sample of the Sloan Digital Sky Survey Data Release 7, we construct a distance-redshift diagram over the redshift range from 0.05 to 0.2 with the massE ruler. The best-fit dark energy density is 0.675+-0.079 for flat Lambda-CDM, consistent with those by other probes. At the median redshift of 0.11, the median distance is estimated to have a fractional error of 0.34%, much lower than those by supernova (SN) Ia and baryonic acoustic oscillation (BAO) and even exceeding their future capability at this redshift. The above low-z measurement is useful for probing dark energy that dominates at the late Universe. For a flat dark energy equation of state model (flat wCDM), the massE alone constrains w to an error that is only a factor of 2.2, 1.7 and 1.3 times larger than those by BAO, SN Ia, and cosmic microwave background (CMB), respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Systematic KMTNet Planetary Anomaly Search, Paper I: OGLE-2019-BLG-1053Lb, A Buried Terrestrial Planet

In order to exhume the buried signatures of "missing planetary caustics" in the KMTNet data, we conducted a systematic anomaly search to the residuals from point-source point-lens fits, based on a modified version of the KMTNet EventFinder algorithm. This search reveals the lowest mass-ratio planetary caustic to date in the microlensing event OGLE-2019-BLG-1053, for which the planetary signal had not been noticed before. The planetary system has a planet-host mass ratio of $q = (1.25 \pm 0.13) \times 10^{-5}$. A Bayesian analysis yields estimates of the mass of the host star, $M_{\rm host} = 0.61_{-0.24}^{+0.29}~M_\odot$, the mass of its planet, $M_{\rm planet} = 2.48_{-0.98}^{+1.19}~M_{\oplus}$, the projected planet-host separation, $a_\perp = 3.4_{-0.5}^{+0.5}$ au, and the lens distance of $D_{\rm L} = 6.8_{-0.9}^{+0.6}$ kpc. The discovery of this very low mass-ratio planet illustrates the utility of our method and opens a new window for a large and homogeneous sample to study the microlensing planet-host mass-ratio function down to $q \sim 10^{-5}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Systematic KMTNet Planetary Anomaly Search. IV. Complete Sample of 2019 Prime-Field

We report the complete statistical planetary sample from the prime fields ($Γ\geq 2~{\rm hr}^{-1}$) of the 2019 Korea Microlensing Telescope Network (KMTNet) microlensing survey. We develop the optimized KMTNet AnomalyFinder algorithm and apply it to the 2019 KMTNet prime fields. We find a total of 14 homogeneously selected planets and report the analysis of three planetary events, KMT-2019-BLG-(1042,1552,2974). The planet-host mass ratios, $q$, for the three planetary events are $6.34 \times 10^{-4}, 4.89 \times 10^{-3}$ and $6.18 \times 10^{-4}$, respectively. A Bayesian analysis indicates the three planets are all cold giant planets beyond the snow line of their host stars. The 13 planets are basically uniform in $\log q$ over the range $-5.0 < \log q < -1.5$. This result suggests that the planets below $q_{\rm break} = 1.7 \times 10^{-4}$ proposed by the MOA-II survey may be more common than previously believed. This work is an early component of a large project to determine the KMTNet mass-ratio function, and the whole sample of 2016--2019 KMTNet events should contain about 120 planets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

TESS discovery of a sub-Neptune orbiting a mid-M dwarf TOI-2136

We present the discovery of TOI-2136b, a sub-Neptune planet transiting every 7.85 days a nearby M4.5V-type star, identified through photometric measurements from the TESS mission. The host star is located $33$ pc away with a radius of $R_{\ast} = 0.34\pm0.02\ R_{\odot}$, a mass of $0.34\pm0.02\ M_{\odot}$ and an effective temperature of $\rm 3342\pm100\ K$. We estimate its stellar rotation period to be $75\pm5$ days based on archival long-term photometry. We confirm and characterize the planet based on a series of ground-based multi-wavelength photometry, high-angular-resolution imaging observations, and precise radial velocities from CFHT/SPIRou. Our joint analysis reveals that the planet has a radius of $2.19\pm0.17\ R_{\oplus}$, and a mass measurement of $6.4\pm2.4\ M_{\oplus}$. The mass and radius of TOI2136b is consistent with a broad range of compositions, from water-ice to gas-dominated worlds. TOI-2136b falls close to the radius valley for low-mass stars predicted by the thermally driven atmospheric mass loss models, making it an interesting target for future studies of its interior structure and atmospheric properties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The effects of plasma on the magnification and time delay of strongly lensed fast radio bursts

The number of identified Fast Radio Bursts (FRBs) is increasing rapidly with current and future facilities. Strongly lensed FRBs are expected to be found as well, which can provide precise time delays and thus have rich applications in cosmology and fundamental physics. However, the radio signal of lensed FRBs will be deflected by plasma in lens galaxies in addition to gravity. Such deflections by both gravity and plasma will cause frequency dependent time delays, which are different from the dispersion delay and the geometric delay caused by gravitational lensing. Depending on the lensing and plasma models, the frequency-time delay relation of the lensed images can show distinguishing behaviours either between the multiple images, or from the dispersion relation. Such phenomena cannot be neglected in future studies, especially at low radio frequency, as plasma exists in lens galaxies in general. More importantly, such information provides not only a potential way to search for lensed FRBs, but also constraints on the mass and plasma distributions in lens galaxy. In particular, plasma may make the missing central images observable at low radio frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Fornax3D project: discovery of ancient massive merger events in the Fornax cluster galaxies NGC 1380 and NGC 1427

We report the discovery of ancient massive merger events in the early-type galaxies NGC 1380 and NGC 1427, members of the Fornax galaxy cluster. Both galaxies have been observed by the MUSE IFU instrument on the VLT, as part of the Fornax3D project. By fitting recently-developed population-orbital superposition models to the observed surface brightness as well as stellar kinematic, age, and metallicity maps, we obtain the stellar orbits, age and metallicity distributions of each galaxy. We then decompose each galaxy into multiple orbital-based components, including a dynamically hot inner stellar halo component which is identified as the relic of past massive mergers. By comparing to analogues from cosmological galaxy simulations, chiefly TNG50, we find that the formation of such a hot inner stellar halo requires the merger with a now-destroyed massive satellite galaxy of $3.7_{-1.5}^{+2.7} \times 10^{10}$ Msun (about $1/5$ of its current stellar mass) in the case of NGC 1380 and of $1.5_{-0.7}^{+1.6} \times10^{10}$ Msun (about $1/4$ of its current stellar mass) in the case of NGC 1427. Moreover, we infer that the last massive merger in NGC 1380 happened $\sim10$ Gyr ago based on the stellar age distribution of the re-grown dynamically cold disk, whereas the merger in NGC 1427 ended $t\lesssim 8$ Gyr ago based on the stellar populations in its hot inner stellar halo. The major merger event in NGC 1380 is the first one with both merger mass and merger time quantitatively inferred in a galaxy beyond the Local Volume. Moreover, it is the oldest and most massive merger uncovered in nearby galaxies so far.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The North/South Asymmetry of the Galaxy: Possible Connection to the Vertical Phase Space Snail

The Galaxy is found to be in disequilibrium based on recent findings of the North/South (N/S) asymmetry and the phase mixing signatures, such as a phase spiral (snail) structure in the vertical phase space ($z-V_{z}$). We show that the N/S asymmetry in a tracer population of dwarfs may be quantitatively modeled with a simple phase snail model superimposed on a smooth equilibrium background. As the phase snail intersects with the $z$ axis, the number density is enhanced, and the velocity dispersion ($σ_{z}$) is decreased relative to the other side of the Galactic plane. Fitting only to the observed asymmetric N/S $σ_{z}$ profiles, we obtain reasonable parameters for the phase space snail and the potential utilized in modeling the background, despite the complex dependence of the model on the potential parameters and the significant selection effects of the data. Both the snail shape and the N/S number density difference given by our best-fit model are consistent with previous observations. The equilibrium background implies a local dark matter density of $0.0151^{+0.0050}_{-0.0051}$ ${\rm M}_{\odot}\,{\rm pc}^{-3}$. The vertical bulk motion of our model is similar to the observation, but with a $\sim$1.2 $\rm km\,s^{-1}$ shift. Our work demonstrates the strong correlation between the phase space snail and the N/S asymmetry. Future observational constraints will facilitate more comprehensive snail models to unravel the Milky Way potential and the perturbation history encoded in the snail feature.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Validation of 13 Hot and Potentially Terrestrial TESS Planets

The James Webb Space Telescope (JWST) will be able to probe the atmospheres and surface properties of hot, terrestrial planets via emission spectroscopy. We identify 18 potentially terrestrial planet candidates detected by the Transiting Exoplanet Survey Satellite (TESS) that would make ideal targets for these observations. These planet candidates cover a broad range of planet radii ($R_{\rm p} \sim 0.6 - 2.0 R_\oplus$) and orbit stars of various magnitudes ($K_s = 5.78 - 10.78$, $V = 8.4 - 15.69$) and effective temperatures ($T_{\rm eff }\sim 3000 - 6000$ K). We use ground-based observations collected through the TESS Follow-up Observing Program (TFOP) and two vetting tools -- DAVE and TRICERATOPS -- to assess the reliabilities of these candidates as planets. We validate 13 planets: TOI-206 b, TOI-500 b, TOI-544 b, TOI-833 b, TOI-1075 b, TOI-1411 b, TOI-1442 b, TOI-1693 b, TOI-1860 b, TOI-2260 b, TOI-2411 b, TOI-2427 b, and TOI-2445 b. Seven of these planets (TOI-206 b, TOI-500 b, TOI-1075 b, TOI-1442 b, TOI-2260 b, TOI-2411 b, and TOI-2445 b) are ultra-short-period planets. TOI-1860 is the youngest ($133 \pm 26$ Myr) solar twin with a known planet to date. TOI-2260 is a young ($321 \pm 96$ Myr) G dwarf that is among the most metal-rich ([Fe/H] = $0.22 \pm 0.06$ dex) stars to host an ultra-short-period planet. With an estimated equilibrium temperature of $\sim 2600$ K, TOI-2260 b is also the fourth hottest known planet with $R_{\rm p} < 2 \, R_\oplus$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hot and counter-rotating star-forming disk galaxies in IllustrisTNG and their real-world counterparts

A key feature of a large population of low-mass, late-type disk galaxies are star-forming disks with exponential light distributions. They are typically also associated with thin and flat morphologies, blue colors, and dynamically cold stars moving along circular orbits within co-planar thin gas disks. However, the latter features do not necessarily always imply the former, in fact, a variety of different kinematic configurations do exist. In this work, we use the cosmological hydrodynamical IllustrisTNG Simulation to study the nature and origin of dynamically hot, sometimes even counter-rotating, star-forming disk galaxies in the lower stellar mass range (between $5\times 10^9\,\mathrm{M_{\odot}}$ and $2\times 10^{10}\,\mathrm{M_{\odot}}$). We find that being dynamically hot arises in most cases as an induced transient state, for example due to galaxy interactions and merger activities, rather than as an age-dependent evolutionary phase of star-forming disk galaxies. The dynamically hot but still actively star-forming disks show a common feature of hosting kinematically misaligned gas and stellar disks, and centrally concentrated on-going star formation. The former is often accompanied by disturbed gas morphologies, while the latter is reflected in low gas and stellar spins in comparison to their dynamically cold, normal disk counterparts. Interestingly, observed galaxies from MaNGA with kinematic misalignment between gas and stars show remarkably similar general properties as the IllustrisTNG galaxies, and therefore are plausible real-world counterparts. In turn, this allows us to make predictions for the stellar orbits and gas properties of these misaligned galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

KMT-2018-BLG-1988Lb: microlensing super-Earth orbiting a low-mass disk dwarf

We reexamine high-magnification microlensing events in the previous data collected by the KMTNet survey with the aim of finding planetary signals that were not noticed before. In this work, we report the planetary system KMT-2018-BLG-1988L that was found from this investigation. The planetary signal appears as a deviation with $\lesssim 0.2$~mag from a single-lens light curve and lasted for about 6 hours. The deviation exhibits a pattern of a dip surrounded by weak bumps on both sides of the dip. The analysis of the lensing light curve indicates that the signal is produced by a low mass-ratio ($q\sim 4\times 10^{-5}$) planetary companion located near the Einstein ring of the host star. The mass of the planet, $M_{\rm planet}=6.8^{+4.7}_{-3.5}~M_\oplus$ and $5.6^{+3.8}_{-2.8}~M_\oplus$ for the two possible solutions, estimated from the Bayesian analysis indicates that the planet is in the regime of a super-Earth. The host of the planet is a disk star with a mass of $M_{\rm host} = 0.47^{+0.33}_{-0.25}~M_\odot$ and a distance of $D_{\rm L}= 4.2^{+1.8}_{-.14}$~kpc. KMT-2018-BLG-1988Lb is the seventeenth microlensing planet with a mass below the upper limit of a super-Earth. The fact that 14 out of 17 microlensing planets with masses $\lesssim 10~M_\oplus$ were detected during the last 5 years since the full operation of the KMTNet survey indicates that the KMTNet database is an important reservoir of very low-mass planets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

KMT-2021-BLG-0912Lb: A microlensing super Earth around a K-type star

The light curve of the microlensing event KMT-2021-BLG-0912 exhibits a very short anomaly relative to a single-lens single-source form. We investigate the light curve for the purpose of identifying the origin of the anomaly. We model the light curve under various interpretations. From this, we find four solutions, in which three solutions are found under the assumption that the lens is composed of two masses (2L1S models), and the other solution is found under the assumption that the source is comprised of a binary-star system (1L2S model). The 1L2S model is ruled out based on the contradiction that the faint source companion is bigger than its primary, and one of the 2L1S solutions is excluded from the combination of the relatively worse fit, blending constraint, and lower overall probability, leaving two surviving solutions with the planet/host mass ratios of $q\sim 2.8\times 10^{-5}$ and $\sim 1.1\times 10^{-5}$. A subtle central deviation supports the possibility of a tertiary lens component, either a binary companion to the host with a very large or small separation or a second planet lying near the Einstein ring, but it is difficult to claim a secure detection due to the marginal fit improvement, lack of consistency among different data sets, and difficulty in uniquely specifying the nature of the tertiary component. With the observables of the event, it is estimated that the masses of the planet and host are $\sim (6.9~M_\oplus, 0.75~M_\odot)$ according to one solution and $\sim (2.8~M_\oplus, 0.80~M_\odot)$ according to the other solution, indicating that the planet is a super Earth around a K-type star, regardless of the solution.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Light Curve Calculations for Triple Microlensing Systems

We present a method to compute the magnification of a finite source star lensed by a triple lens system based on the image boundary (contour integration) method. We describe a new procedure to obtain continuous image boundaries from solutions of the tenth-order polynomial obtained from the lens equation. Contour integration is then applied to calculate the image areas within the image boundaries, which yields the magnification of a source with uniform brightness. We extend the magnification calculation to limb-darkened stars approximated with a linear profile. In principle, this method works for all multiple lens systems, not just triple lenses. We also include an adaptive sampling and interpolation method for calculating densely covered light curves. The C++ source code and a corresponding Python interface are publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Microlensing Predictions: Impact of Galactic Disc Dynamical Models

Galactic model plays an important role in the microlensing field, not only for analyses of individual events but also for statistics of the ensemble of events. However, the Galactic models used in the field varies, and some are unrealistically simplified. Here we tested three Galactic disc dynamic models, the first is a simple standard model that was widely used in this field, whereas the other two consider the radial dependence of the velocity dispersion, and in the last model, the asymmetric drift. We found that for a typical lens mass $M_{\rm L}=0.5M_{\odot}$, the two new dynamical models predict $\sim16\%$ or $\sim5\%$ less long-timescale events (e.g., microlensing timescale $t_{\rm E}>300$ days) and $\sim 5\%$ and $\sim 3.5\%$ more short-timescale events ($t_{\rm E}<3$ days) than the standard model. Moreover, the microlensing event rate as a function of Einstein radius $θ_{\rm E}$ or microlensing parallax $π_{\rm E}$ also shows some model dependence (a few percent). The two new models also have an impact on the total microlensing event rate. This result will also to some degree affect the Bayesian analysis of individual events, but overall, the impact is small. However, we still recommend that modelers should be more careful when choosing the Galactic model, especially in statistical works involving Bayesian analyses of a large number of events. Additionally, we find the asymptotic power-law behaviors in both $θ_{\rm E}$ and $π_{\rm E}$ distributions, and we provide a simple model to understand them.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Systematic KMTNet Planetary Anomaly Search, Paper II: Six New $q<2\times 10^{-4}$ Mass-ratio Planets

We apply the automated AnomalyFinder algorithm of Paper I (Zang et al. 2021b) to 2018-2019 light curves from the $\simeq 13\,{\rm deg}^2$ covered by the six KMTNet prime fields, with cadences $Γ\geq 2\,{\rm hr}^{-1}$. We find a total of 11 planets with mass ratios $q<2\times 10^{-4}$, including six newly discovered planets, one planet that was reported in Paper I, and recovery of four previously discovered planets. One of the new planets, OGLE-2018-BLG-0977Lb, is in a planetary-caustic event, while the other five (OGLE-2018-BLG-0506Lb, OGLE-2018-BLG-0516Lb, OGLE-2019-BLG-1492Lb, KMT-2019-BLG-0253, and KMT-2019-BLG-0953) are revealed by a "dip" in the light curve as the source crosses the host-planet axis on the opposite side of the planet. These subtle signals were missed in previous by-eye searches. The planet-host separations (scaled to the Einstein radius), $s$, and planet-host mass ratios, $q$, are, respectively, $(s,q\times 10^5) = (0.88, 4.1)$, $(0.96\pm 0.10, 8.3)$, $(0.94\pm 0.07, 13)$, $(0.97\pm 0.07, 18)$, $(0.97\pm0.04,4.1)$, and $(0.74,18)$, where the "$\pm$" indicates a discrete degeneracy. The 11 planets are spread out over the range $-5<\log q < -3.7$. Together with the two planets previously reported with $q\sim 10^{-5}$ from the 2018-2019 non-prime KMT fields, this result suggests that planets toward the bottom of this mass-ratio range may be more common than previously believed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A discrete chemo-dynamical model of M87's globular clusters: Kinematics extending to ~ 400 kpc

We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/σ\sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/σ\sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An accurate analytic mass model for lensing galaxies

We develop an analytic mass model for lensing galaxies, based on a broken power-law (BPL) density profile, which is a power-law profile with a mass deficit or surplus in the central region. Under the assumption of an elliptically symmetric surface mass distribution, the deflection angle and magnification can be evaluated analytically for this new model. We compute the theoretical prediction for various quantities, including the volume and surface mass density profiles of the galaxies, and the aperture and luminosity-weighted line-of-sight velocity dispersions, and compare them to those measured from the Illustris simulation. We find an excellent agreement between our model prediction and the simulation, which validates our modeling. The high efficiency and accuracy of our model manifests itself as a promising tool for studying properties of galaxies with strong lensing.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Basis function expansions for galactic dynamics: Spherical versus cylindrical coordinates

Aims. The orbital structure of galaxies is strongly influenced by the accuracy of the force calculation during orbit integration. We explore the accuracy of force calculations for two expansion methods and determine which one is preferable for orbit integration. Methods. We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations.We specifically compare two methods, one was introduced by Hernquist \& Ostriker (HO), which uses a spherical coordinate system and was built specifically for the Hernquist model, and the other by Vasiliev \& Athanassoula (CylSP) has a cylindrical coordinate system. Our comparisons include the Dehnen profile, its triaxial extension (of which the Hernquist profile is a special case) and a multicomponent system including a bar and disk density distributions for both analytical models and $N$-body realizations. Results. For the generalized Dehnen density, the CylSP method is more accurate than the HO method for nearly all inner power-law indices and shapes at all radii. For $N$-body realizations of the Dehnen model, or snapshots of an $N$-body simulation, the CylSP method is more accurate than the HO method in the central region for the oblate, prolate, and triaxial Hernquist profiles if the particle number is more than $5\times 10^5$. For snapshots of the Hernquist models with spherical shape, the HO method is preferred. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Disentangling the formation history of galaxies via population-orbit superposition: method validation

We present population-orbit superposition models for external galaxies based on Schwarzschild's orbit-superposition method, by tagging the orbits with age and metallicity. The models fit the density distributions, as well as kinematic, age and metallicity maps from Integral Field Unit (IFU) spectroscopy observations. We validate the method and demonstrate its power by applying it to mock data, similar to those obtained by the Multi-Unit Spectroscopic Explorer (MUSE) IFU on the Very Large Telescope (VLT). These mock data are created from Auriga galaxy simulations, viewed at three different inclination angles ($\vartheta=40^o, 60^o, 80^o$). Constrained by MUSE-like mock data, our model can recover the galaxy's stellar orbit distribution projected in orbital circularity $λ_z$ vs. radius $r$, the intrinsic stellar population distribution in age $t$ vs. metallicity $Z$, and the correlation between orbits' circularity $λ_z$ and stellar age $t$. A physically motivated age-metallicity relation improves recovering the intrinsic stellar population distributions. We decompose galaxies into cold, warm and hot + counter-rotating components based on their orbit circularity distribution, and find that the surface density, mean velocity, velocity dispersion, age and metallicity maps of each component from our models well reproduce those from simulation, especially for projections close to edge-on. These galaxies exhibit strong global age vs. $σ_z$ relation, which is well recovered by our model. The method has the power to reveal the detailed build-up of stellar structures in galaxies, and offers a complement to local resolved, and high-redshift studies of galaxy evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Early-type galaxy density profiles from IllustrisTNG: I. Galaxy correlations and the impact of baryons

We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langleγ^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $σ_{γ^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $γ^{\prime}$ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction and in-situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, $γ^{\prime}$ is almost constant with redshift below $z = 2$. Through analyzing IllustrisTNG model variations we show that black hole kinetic winds are crucial to lowering $γ^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $γ^{\prime}$ are subdominant compared to AGN feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter halos are well-described by steeper-than-NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter only (DMO) run. Their inner density slopes anti-correlates (remain constant) with the halo mass in the FP (DMO) run, and anti-correlates with the halo concentration parameter $c_{200}$ in both types of runs. The dark matter halos of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

LHS 1815b: The First Thick-Disk Planet Detected By TESS

We report the first discovery of a thick-disk planet, LHS 1815b (TOI-704b, TIC 260004324), detected in the TESS survey. LHS 1815b transits a bright (V = 12.19 mag, K = 7.99 mag) and quiet M dwarf located $ 29.87\pm0.02 pc$ away with a mass of $0.502\pm0.015 M_{\odot}$ and a radius of $0.501\pm0.030 R_{\odot}$. We validate the planet by combining space and ground-based photometry, spectroscopy, and imaging. The planet has a radius of $1.088\pm 0.064 R_{\oplus}$ with a $3 σ$ mass upper-limit of $8.7 M_{\oplus}$. We analyze the galactic kinematics and orbit of the host star LHS1815 and find that it has a large probability ($P_{thick}/P_{thin} = 6482$) to be in the thick disk with a much higher expected maximal height ($Z_{max} = 1.8 kpc$) above the Galactic plane compared with other TESS planet host stars. Future studies of the interior structure and atmospheric properties of planets in such systems using for example the upcoming James Webb Space Telescope (JWST), can investigate the differences in formation efficiency and evolution for planetary systems between different Galactic components (thick and thin disks, and halo).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the local dark matter density with LAMOST DR5 and Gaia DR2

We apply the vertical Jeans equation to the kinematics of Milky Way stars in the solar neighbourhood to measure the local dark matter density. More than 90,000 G- and K-type dwarf stars are selected from the cross-matched sample of LAMOST DR5 and Gaia DR2 for our analyses. The mass models applied consist of a single exponential stellar disc, a razor thin gas disc and a constant dark matter density. We first consider the simplified vertical Jeans equation which ignores the tilt term and assumes a flat rotation curve. Under a Gaussian prior on the total stellar surface density, the local dark matter density inferred from Markov Chain Monte Carlo simulations is $0.0133_{-0.0022}^{+0.0024}\ {\rm M}_{\odot}\,{\rm pc}^{-3}$. The local dark matter densities for subsamples in an azimuthal angle range of $-10^{\circ} < ϕ< 5^{\circ}$ are consistent within their 1$σ$ errors. However, the northern and southern subsamples show a large discrepancy due to plateaux in the northern and southern vertical velocity dispersion profiles. These plateaux may be the cause of the different estimates of the dark matter density between the north and south. Taking the tilt term into account has little effect on the parameter estimations and does not explain the north and south asymmetry. Taking half of the difference of $σ_{z}$ profiles as unknown systematic errors, we then obtain consistent measurements for the northern and southern subsamples. We discuss the influence of the vertical data range, the scale height of the tracer population, the vertical distribution of stars and the sample size on the uncertainty of the determination of the local dark matter density.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Redshift Evolution of the Fundamental Plane Relation in the IllustrisTNG Simulation

We investigate the fundamental plane (FP) evolution of early-type galaxies in the IllustrisTNG-100 simulation (TNG100) from redshift $z=0$ to $z=2$. We find that a tight plane relation already exists as early as $z=2$. Its scatter stays as low as $\sim 0.08$ dex across this redshift range. Both slope parameters $b$ and $c$ (where $R \propto σ^b I^c$ with $R$, $σ$, and $I$ being the typical size, velocity dispersion, and surface brightness) of the plane evolve mildly since $z=2$, roughly consistent with observations. The FP residual $\rm Res$ ($\equiv\,a\,+\,b\log σ\,+\,c\log I\,-\,\log R$, where $a$ is the zero point of the FP) is found to strongly correlate with stellar age, indicating that stellar age can be used as a crucial fourth parameter of the FP. However, we find that $4c+b+2=δ$, where $δ\sim 0.8$ for FPs in TNG, rather than zero as is typically inferred from observations. This implies that a tight power-law relation between the dynamical mass-to-light ratio $M_{\rm dyn}/L$ and the dynamical mass $M_{\rm dyn}$ (where $M_{\rm dyn}\equiv 5σ^2R/G$, with $G$ being the gravitational constant) is not present in the TNG100 simulation. Recovering such a relation requires proper mixing between dark matter and baryons, as well as star formation occurring with correct efficiencies at the right mass scales. This represents a powerful constraint on the numerical models, which has to be satisfied in future hydrodynamical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SDSS-IV MaNGA: Stellar population correlates with stellar root-mean-square velocity $V_{\rm rms}$ gradients or total-density-profile slopes at fixed effective velocity dispersion $σ_{\rm e}$

Galaxy properties are known to correlate most tightly with the galaxy effective stellar velocity dispersion $σ_{\rm e}$. Here we look for {\em additional} trends at fixed $σ_{\rm e}$ using 1339 galaxies ($M_\ast \gtrsim 6\times10^9$ M$_\odot$) with different morphologies in the MaNGA (DR14) sample with integral-field spectroscopy data. We focus on the gradients ($γ_{\rm rms} \equiv σ(R_{\rm e}/4)/σ_{\rm e}$) of the stellar root-mean-square velocity ($V_{\rm rms} \equiv \sqrt{V^2 + σ^2}$), which we show traces the total mass density gradient $γ_{\rm tot}$ derived from dynamical models and, more weakly, the bulge fraction. We confirm that $γ_{\rm rms}$ increases with $σ_{\rm e}$, age and metallicity. We additionally find that these correlations still exist at fixed $σ_{\rm e}$, where galaxies with larger $γ_{\rm rms}$ are found to be older and more metal-rich. It means that mass density gradients contain information of the stellar population which is not fully accounted for by $σ_{\rm e}$. This result puts an extra constraint on our understanding of galaxy quenching. We compare our results with galaxies in the IllustrisTNG hydrodynamical simulations and find that, at fixed $σ_{\rm e}$, similar trends exist with age, the bulge fraction, and the total mass density slope but, unlike observations, no correlation with metallicity can be detected in the simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structural and Stellar Population Properties vs. Bulge Types in Sloan Digital Sky Survey Central Galaxies

This paper studies pseudo-bulges (P-bulges) and classical bulges (C-bulges) in Sloan Digital Sky Survey central galaxies using the new bulge indicator $ΔΣ_1$, which measures relative central stellar-mass surface density within 1 kpc. We compare $ΔΣ_1$ to the established bulge-type indicator $Δ\langleμ_e\rangle$ from Gadotti (2009) and show that classifying by $ΔΣ_1$ agrees well with $Δ\langleμ_e\rangle$. $ΔΣ_1$ requires no bulge-disk decomposition and can be measured on SDSS images out to $z = 0.07$. Bulge types using it are mapped onto twenty different structural and stellar-population properties for 12,000 SDSS central galaxies with masses 10.0 < log $M_*$/$M_{\odot}$ < 10.4. New trends emerge from this large sample. Structural parameters show fairly linear log-log relations vs. $ΔΣ_1$ and $Δ\langleμ_e\rangle$ with only moderate scatter, while stellar-population parameters show a highly non-linear "elbow" in which specific star-formation rate remains roughly flat with increasing central density and then falls rapidly at the elbow, where galaxies begin to quench. P-bulges occupy the low-density end of the horizontal arm of the elbow and are universally star-forming, while C-bulges occupy the elbow and the vertical branch and exhibit a wide range of star-formation rates at fixed density. The non-linear relation between central density and star-formation rate has been seen before, but this mapping onto bulge class is new. The wide range of star-formation rates in C-bulges helps to explain why bulge classifications using different parameters have sometimes disagreed in the past. The elbow-shaped relation between density and stellar indices suggests that central structure and stellar-populations evolve at different rates as galaxies begin to quench.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The influence of the secular perturbation of an intermediate-mass companion: I. Eccentricity excitation of disk stars at the Galactic center

There is a dense group of OB and Wolf-Rayet stars within a fraction of a parsec from the super-massive black hole (SMBH) at the Galactic Center. These stars appear to be coeval and relatively massive. A subgroup of these stars orbits on the same plane. If they emerged with low to modest eccentricity orbits from a common gaseous disk around the central super-massive black hole, their inferred lifespan would not be sufficiently long to account for the excitation of their high orbital eccentricity through dynamical relaxation. Here we analyze the secular perturbation on Galactic Center stars by an intermediate-mass companion (IMC) as a potential mechanism to account for these young disk stars' high eccentricity. This IMC may be either an intermediate-mass black hole (IMBH) or a compact cluster such as IRS-13E. If its orbital angular momentum vector is anti-parallel to that of the disk stars, this perturbation would be effective in exciting the eccentricity of stars with orbital precession rates which resonate with IMC's precession rate. If it orbits around the SMBH in the same direction as the disk stars, the eccentricity of the young stars can still be highly excited by the IMC during the depletion of their natal disk, possible associated with the launch of the Fermi bubble. In this scenario, IMC's precession rate decreases and its secular resonance sweeps through the proximity of the young stars. We carry out numerical simulations with various inclination angles between the orbits of IMC and the disk stars and show this secular interaction is a robust mechanism to excite the eccentricity and inclination of some disk stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Understanding the velocity distribution of the Galactic Bulge with APOGEE and Gaia

We revisit the stellar velocity distribution in the Galactic bulge/bar region with APOGEE DR16 and {\it Gaia} DR2, focusing in particular on the possible high-velocity (HV) peaks and their physical origin. We fit the velocity distributions with two different models, namely with Gauss-Hermite polynomial and Gaussian mixture model (GMM). The result of the fit using Gauss-Hermite polynomials reveals a positive correlation between the mean velocity ($\bar{V}$) and the "skewness" ($h_{3}$) of the velocity distribution, possibly caused by the Galactic bar. The $n=2$ GMM fitting reveals a symmetric longitudinal trend of $|μ_{2}|$ and $σ_{2}$ (the mean velocity and the standard deviation of the secondary component), which is inconsistent to the $x_{2}$ orbital family predictions. Cold secondary peaks could be seen at $|l|\sim6^\circ$. However, with the additional tangential information from {\it Gaia}, we find that the HV stars in the bulge show similar patterns in the radial-tangential velocity distribution ($V_{\rm R}-V_{\rm T}$), regardless of the existence of a distinct cold HV peak. The observed $V_{\rm R}-V_{\rm T}$ (or $V_{\rm GSR}-μ_{l}$) distributions are consistent with the predictions of a simple MW bar model. The chemical abundances and ages inferred from ASPCAP and CANNON suggest that the HV stars in the bulge/bar are generally as old as, if not older than, the other stars in the bulge/bar region.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

KMT-2016-BLG-1836Lb: A Super-Jovian Planet From A High-Cadence Microlensing Field

We report the discovery of a super-Jovian planet in the microlensing event KMT-2016-BLG-1836, which was found by the Korea Microlensing Telescope Network's high-cadence observations (Γ~ 4~{hr}^{-1}). The planet-host mass ratio q ~ 0.004. A Bayesian analysis indicates that the planetary system is composed of a super-Jovian M_{planet} = 2.2_{-1.1}^{+1.9} M_{J} planet orbiting an M or K dwarf M_{\rm host} = 0.49_{-0.25}^{+0.38} M_{Sun}, at a distance of D_{L} = 7.1_{-2.4}^{+0.8} kpc. The projected planet-host separation is 3.5^{+1.1}_{-0.9} AU, implying that the planet is located beyond the snowline of the host star. Future high-resolution images can potentially strongly constrain the lens brightness and thus the mass and distance of the planetary system. Without considering detailed detection efficiency, selection or publication biases, we find a potential "mass ratio desert" at -3.7 \lesssim \log q \lesssim -3.0 for the 31 published KMTNet planets.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

OGLE-2015-BLG-1649Lb: A gas giant planet around a low-mass dwarf

We report the discovery of an exoplanet in microlensing event OGLE-2015-BLG-1649. The planet/host-star mass ratio is $q =7.2 \times 10^{-3}$ and the projected separation normalized by the Einstein radius is $s = 0.9$. The upper limit of the lens flux is obtained from adaptive optics observations by IRCS/Subaru, which excludes the probability of a G-dwarf or more massive host star and helps to put a tighter constraint on the lens mass as well as commenting on the formation scenarios of giant planets orbiting low-mass stars. We conduct a Bayesian analysis including constraints on the lens flux to derive the probability distribution of the physical parameters of the lens system. We thereby find that the masses of the host star and planet are $M_{L} = 0.34 \pm 0.19 M_{\odot}$ and $M_{p} = 2.5^{+1.5}_{-1.4} M_{Jup}$, respectively. The distance to the system is $D_{L} = 4.23^{+1.51}_{-1.64}$kpc. The projected star-planet separation is $a_{\perp} = 2.07^{+0.65}_{-0.77}$AU. The lens-source relative proper motion of the event is quite high, at $\sim 7.1 \, {\rm mas/yr}$. Therefore, we may be able to determine the lens physical parameters uniquely or place much stronger constraints on them by measuring the color-dependent centroid shift and/or the image elongation with additional high resolution imaging already a few years from now.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

OGLE-2015-BLG-1771Lb: A Microlens Planet Orbiting an Ultracool Dwarf?

We report the discovery and the analysis of the short (tE < 5 days) planetary microlensing event, OGLE-2015-BLG-1771. The event was discovered by the Optical Gravitational Lensing Experiment (OGLE), and the planetary anomaly (at I ~ 19) was captured by The Korea Microlensing Telescope Network (KMTNet). The event has three surviving planetary models that explain the observed light curves, with planet-host mass ratio q \~ 5.4 * 10^{-3}, 4.5 * 10^{-3} and 4.5 * 10^{-2}, respectively. The first model is the best-fit model, while the second model is disfavored by Δχ^2 ~ 3. The last model is strongly disfavored by Δχ^2 ~ 15 but not ruled out. A Bayesian analysis using a Galactic model indicates that the first two models are probably composed of a Saturn-mass planet orbiting a late M dwarf, while the third one could consist of a super-Jovian planet and a mid-mass brown dwarf. The source-lens relative proper motion is mu_rel ~ 9 mas/yr, so the source and lens could be resolved by current adaptive-optics (AO) instruments in 2021 if the lens is luminous.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

SDSS-IV MaNGA: Internal mass distributions and orbital structures of early-type galaxies and their dependence on environment

In our earlier 2019 paper, we evaluated the reliability of Schwarzschild's orbit-superposition dynamical modelling method in estimating the internal mass distribution, intrinsic stellar shapes and orbit distributions of early-type galaxies (ETGs) taken from the Illustris cosmological simulation. We now apply the same techniques to galaxies taken from the integral-field survey Mapping Nearby Galaxies with APO (MaNGA), using a sample of 149 ETGs in the mass range of $10^{9.90}\sim10^{11.80} M_{\odot}$ and made up of 105 central and 44 satellite galaxies. We find that low-mass ETGs with $\log(M_*/M_{\odot})<11.1$ have an average dark matter fraction of $\sim0.2$ within one effective radius $R_{\rm e}$, tend to be oblate-like, and are dominated by rotation about their minor axis. High-mass ETGs with $\log(M_*/M_{\odot})>11.1$ have an average dark matter fraction of $\sim0.4$ within one effective radius $R_{\rm e}$, tend to be prolate-like, and are dominated by rotation about their major axis and by centrophilic orbits. The changes of internal structures within one $R_{\rm e}$ are dominated by the total stellar mass of the individual galaxies. We find no differences of internal structures between central and satellite ETGs for the same stellar masses. However, for similar stellar mass and colour distributions, we find that ETGs more prolate-like, or with more hot orbits, tend to have higher close neighbour counts at $r_p\sim40$ kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Spitzer + VLTI-GRAVITY Measure the Lens Mass of a Nearby Microlensing Event

We report the lens mass and distance measurements of the nearby microlensing event TCP J05074264+2447555. We measure the microlens parallax vector $π_{\rm E}$ using Spitzer and ground-based light curves with constraints on the direction of lens-source relative proper motion derived from Very Large Telescope Interferometer (VLTI) GRAVITY observations. Combining this $π_{\rm E}$ determination with the angular Einstein radius $θ_{\rm E}$ measured by VLTI GRAVITY observations, we find that the lens is a star with mass $M_{\rm L} = 0.495 \pm 0.063~M_{\odot}$ at a distance $D_{\rm L} = 429 \pm 21~{\rm pc}$. We find that the blended light basically all comes from the lens. The lens-source proper motion is $μ_{\rm rel,hel} = 26.55 \pm 0.36~{\rm mas\,yr^{-1}}$, so with currently available adaptive-optics (AO) instruments, the lens and source can be resolved in 2021. This is the first microlensing event whose lens mass is unambiguously measured by interferometry + satellite parallax observations, which opens a new window for mass measurements of isolated objects such as stellar-mass black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Spitzer Microlensing parallax reveals two isolated stars in the Galactic bulge

We report the mass and distance measurements of two single-lens events from the 2017 Spitzer microlensing campaign. The ground-based observations yield the detection of finite-source effects, and the microlens parallaxes are derived from the joint analysis of ground-based observations and Spitzer observations. We find that the lens of OGLE-2017-BLG-1254 is a $0.60 \pm 0.03 M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.11~\text{kpc}$, where $D_{\rm LS}$ is the distance between the lens and the source. The second event, OGLE-2017-BLG-1161, is subject to the known satellite parallax degeneracy, and thus is either a $0.51^{+0.12}_{-0.10} M_{\odot}$ star with $D_{\rm LS} = 0.40 \pm 0.12~\text{kpc}$ or a $0.38^{+0.13}_{-0.12} M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.19~\text{kpc}$. Both of the lenses are therefore isolated stars in the Galactic bulge. By comparing the mass and distance distributions of the eight published Spitzer finite-source events with the expectations from a Galactic model, we find that the Spitzer sample is in agreement with the probability of finite-source effects occurrence in single lens events.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Free-floating planets from core accretion theory: microlensing predictions

We calculate the microlensing event rate and typical time-scales for the free-floating planet (FFP) population that is predicted by the core accretion theory of planet formation. The event rate is found to be ~$1.8\times 10^{-3}$ of that for the stellar population. While the stellar microlensing event time-scale peaks at around 20 days, the median time-scale for FFP events (~0.1 day) is much shorter. Our values for the event rate and the median time-scale are significantly smaller than those required to explain the \cite{Sum+11} result, by factors of ~13 and ~16, respectively. The inclusion of planets at wide separations does not change the results significantly. This discrepancy may be too significant for standard versions of both the core accretion theory and the gravitational instability model to explain satisfactorily. Therefore, either a modification to the planet formation theory is required, or other explanations to the excess of short-time-scale microlensing events are needed. Our predictions can be tested by ongoing microlensing experiment such as KMTNet, and by future satellite missions such as WFIRST and Euclid.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Orbital classification in an N-body bar

The dynamics and evolution of any galactic structure are strongly influenced by the properties of the orbits that constitute it. In this paper, we compare two orbit classification schemes, one by Laskar (NAFF) , and the other by Carpintero and Aguilar (CA), by applying both of them to orbits obtained by following individual particles in a numerical simulation of a barred galaxy. We find that, at least for our case and some provisos, the main frequencies calculated by the two methods are in good agreement: for $80\%$ of the orbits the difference between the results of the two methods is less than $5\%$ for all three main frequencies. However, it is difficult to evaluate the amount of regular or chaotic bar orbits in a given system. The fraction of regular orbits obtained by the NAFF method strongly depends on the critical frequency drift parameter, while in the CA method the number of fundamental frequencies strongly depends on the frequency difference parameter $L_{\rm r}$ and the maximum integer used for searching the linear independence of the fundamental frequencies. We also find that, for a given particle, in general the projection of its motion along the bar minor axis is more regular than the other two projections, while the projection along the intermediate axis is the least regular.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SDSS-IV MaNGA IFS Galaxy Survey --- Survey Design, Execution, and Initial Data Quality

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SDSS-IV MaNGA: environmental dependence of stellar age and metallicity gradients in nearby galaxies

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The BOSS Emission-Line Lens Survey. III. : Strong Lensing of Ly$α$ Emitters by Individual Galaxies

We introduce the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$α$ EmitteR sYstems (BELLS GALLERY) Survey, which is a Hubble Space Telescope program to image a sample of galaxy-scale strong gravitational lens candidate systems with high-redshift Ly$α$ emitters (LAEs) as the background sources. The goal of the BELLS GALLERY Survey is to illuminate dark substructures in galaxy-scale halos by exploiting the small-scale clumpiness of rest-frame far-UV emission in lensed LAEs, and to thereby constrain the slope and normalization of the substructure-mass function. In this paper, we describe in detail the spectroscopic strong-lens selection technique, which is based on methods adopted in the previous Sloan Lens ACS (SLACS) Survey, BELLS, and SLACS for the Masses Survey. We present the BELLS GALLERY sample of the 21 highest-quality galaxy--LAE candidates selected from $\approx 1.4 \times 10^6$ galaxy spectra in the BOSS of the Sloan Digital Sky Survey III. These systems consist of massive galaxies at redshifts of approximately 0.5 strongly lensing LAEs at redshifts from 2--3. The compact nature of LAEs makes them an ideal probe of dark substructures, with a substructure-mass sensitivity that is unprecedented in other optical strong-lens samples. The magnification effect from lensing will also reveal the structure of LAEs below 100 pc scales, providing a detailed look at the sites of the most concentrated unobscured star formation in the universe. The source code used for candidate selection is available for download as a part of this release.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The BOSS Emission-Line Lens Survey. IV. : Smooth Lens Models for the BELLS GALLERY Sample

We present \textsl{Hubble Space Telescope} (\textsl{HST}) F606W-band imaging observations of 21 galaxy-Ly$α$ emitter lens candidates in the Baryon Oscillation Spectroscopic Survey (BOSS) Emission-Line Lens Survey (BELLS) for GALaxy-Ly$α$ EmitteR sYstems (BELLS GALLERY) survey. 17 systems are confirmed to be definite lenses with unambiguous evidence of multiple imaging. The lenses are primarily massive early-type galaxies (ETGs) at redshifts of approximately $0.55$, while the lensed sources are Ly$α$ emitters (LAEs) at redshifts from 2 to 3. Although the \textsl{HST} imaging data are well fit by smooth lens models consisting of singular isothermal ellipsoids in an external shear field, a thorough exploration of dark substructures in the lens galaxies is required. The Einstein radii of the BELLS GALLERY lenses are on average $60\%$ larger than those of the BELLS lenses because of the much higher source redshifts which will allow a detailed investigation of the radius evolution of the mass profile in ETGs. With the aid of the average $\sim 13 \times$ lensing magnification, the LAEs are resolved to comprise individual star-forming knots of a wide range of properties with characteristic sizes from less than 100 pc to several kpc, rest-frame far UV apparent AB magnitudes from 29.6 to 24.2, and typical projected separations of 500 pc to 2 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The InfraRed Imaging Spectrograph (IRIS) for TMT: latest science cases and simulations

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The role of major gas-rich mergers on the evolution of galaxies from the blue cloud to the red sequence

With the aim of exploring the fast evolutionary path from the blue cloud of star-forming galaxies to the red sequence of quiescent galaxies in the local universe, we select a local advanced merging infrared luminous and ultraluminous galaxy (adv-merger (U)LIRGs) sample and perform careful dust extinction corrections to investigate their positions in the SFR-$M_{\ast}$, u-r and NUV-r color-mass diagrams. The sample consists of 89 (U)LIRGs at the late merger stage, obtained from cross-correlating the IRAS Point Source Catalog Redshift Survey and 1 Jy ULIRGs samples with the Sloan Digital Sky Survey DR7 database. Our results show that $74\%\pm 5\%$ of adv-merger (U)LIRGs are localized above the $1\, σ$ line of the local star-forming galaxy main sequence. We also find that all adv-merger (U)LIRGs are more massive than and as blue as the blue cloud galaxies after corrections of Galactic and internal dust extinctions, with $95\%\pm 2\%$ and $81\%\pm 4\%$ of them outside the blue cloud on the u-r and NUV-r color-mass diagrams, respectively. These results, combined with the short timescale for exhausting the molecular gas reservoir in adv-merger (U)LIRGs ($3\times 10^{7}$ to $3\times 10^{8}$ years), imply that the adv-merger (U)LIRGs are likely at the starting point of the fast evolutionary track previously proposed by several groups. While the number density of adv-merger (U)LIRGs is only $\sim 0.1\%$ of the blue cloud star-forming galaxies in the local universe, this evolutionary track may play a more important role at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The UV-optical Color Gradients in Star-Forming Galaxies at 0.5<z<1.5: Origins and Link to Galaxy Assembly

The rest-frame UV-optical (i.e., NUV-B) color index is sensitive to the low-level recent star formation and dust extinction, but it is insensitive to the metallicity. In this Letter, we have measured the rest-frame NUV-B color gradients in ~1400 large ($\rm r_e>0.18^{\prime\prime}$), nearly face-on (b/a>0.5) main-sequence star-forming galaxies (SFGs) between redshift 0.5 and 1.5 in the CANDELS/GOODS-S and UDS fields. With this sample, we study the origin of UV-optical color gradients in the SFGs at z~1 and discuss their link with the buildup of stellar mass. We find that the more massive, centrally compact, and more dust extinguished SFGs tend to have statistically more negative raw color gradients (redder centers) than the less massive, centrally diffuse, and less dusty SFGs. After correcting for dust reddening based on optical-SED fitting, the color gradients in the low-mass ($M_{\ast} <10^{10}M_{\odot}$) SFGs generally become quite flat, while most of the high-mass ($M_{\ast} > 10^{10.5}M_{\odot}$) SFGs still retain shallow negative color gradients. These findings imply that dust reddening is likely the principal cause of negative color gradients in the low-mass SFGs, while both increased central dust reddening and buildup of compact old bulges are likely the origins of negative color gradients in the high-mass SFGs. These findings also imply that at these redshifts the low-mass SFGs buildup their stellar masses in a self-similar way, while the high-mass SFGs grow inside out.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Assessing the Jeans Anisotropic Multi-Gaussian Expansion method with the Illustris simulation

We assess the effectiveness of the Jeans-Anisotropic-MGE (JAM) technique with a state-of-the-art cosmological hydrodynamic simulation, the Illustris project. We perform JAM modelling on 1413 simulated galaxies with stellar mass M^* > 10^{10}M_{sun}, and construct an axisymmetric dynamical model for each galaxy. Combined with a Markov Chain Monte Carlo (MCMC) simulation, we recover the projected root-mean-square velocity (V_rms) field of the stellar component, and investigate constraints on the stellar mass-to-light ratio, M^*/L, and the fraction of dark matter f_{DM} within 2.5 effective radii (R_e). We find that the enclosed total mass within 2.5 R_e is well constrained to within 10%. However, there is a degeneracy between the dark matter and stellar components with correspondingly larger individual errors. The 1 sigma scatter in the recovered M^*/L is 30-40% of the true value. The accuracy of the recovery of M^*/L depends on the triaxial shape of a galaxy. There is no significant bias for oblate galaxies, while for prolate galaxies the JAM-recovered stellar mass is on average 18% higher than the input values. We also find that higher image resolutions alleviate the dark matter and stellar mass degeneracy and yield systematically better parameter recovery.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Combining time delays and image positions for quadruple lenses: a moment approach

Time delays in gravitational lenses can be used to determine the Hubble constant and the lens potential. In future surveys, many gravitational lenses can be discovered, and their time delays and image positions can in principle be measured. Using an elliptical power-law potential, we show that combinations of image positions and time delays for quadruple lenses yield simple analytical expressions that are connected with observable quantities. These relations can be used to obtain the approximate axis ratio q, the Einstein radius and the slope. We apply this method to RX J1131-1231, and show that our analytical results match the full numerical determinations approximately. Our approach can quickly determine rough values of lens parameters, which can then be used as initial guesses for further refinement through numerical modelling and may be useful for automated lens search in large surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Determining the local dark matter density with LAMOST data

Measurement of the local dark matter density plays an important role in both Galactic dynamics and dark matter direct detection experiments. However, the estimated values from previous works are far from agreeing with each other. In this work, we provide a well-defined observed sample with 1427 G \& K type main-sequence stars from the LAMOST spectroscopic survey, taking into account selection effects, volume completeness, and the stellar populations. We apply a vertical Jeans equation method containing a single exponential stellar disk, a razor thin gas disk, and a constant dark matter density distribution to the sample, and obtain a total surface mass density of $\rm {78.7 ^{+3.9}_{-4.7}\ M_{\odot}\ pc^{-2}}$ up to 1 kpc and a local dark matter density of $0.0159^{+0.0047}_{-0.0057}\,\rm M_{\odot}\,\rm pc^{-3}$. We find that the sampling density (i.e. number of stars per unit volume) of the spectroscopic data contributes to about two-thirds of the uncertainty in the estimated values. We discuss the effect of the tilt term in the Jeans equation and find it has little impact on our measurement. Other issues, such as a non-equilibrium component due to perturbations and contamination by the thick disk population, are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

How well can cold-dark-matter substructures account for the observed radio flux-ratio anomalies?

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Kinematics of the X-shaped Milky Way Bulge: Expectations from a Self-consistent N-body Model

We explore the kinematics (both the radial velocity and the proper motion) of the vertical X-shaped feature in the Milky Way with an N-body bar/bulge model. From the solar perspective, the distance distribution of particles is double-peaked in fields passing through the X-shape. The separation and amplitude ratio between the two peaks qualitatively match the observed trends towards the Galactic bulge. We confirm clear signatures of cylindrical rotation in the pattern of mean radial velocity across the bar/bulge region. We also find possible imprints of coherent orbital motion inside the bar structure in the radial velocity distribution along l=0 degree, where the near and far sides of the bar/bulge show excesses of approaching and receding particles. The coherent orbital motion is also reflected in the slight displacement of the zero-velocity-line in the mean radial velocity, and the displacement of the maximum/minimum in the mean longitudinal proper motion across the bulge region. We find some degree of anisotropy in the stellar velocity within the X-shape, but the underlying orbital family of the X-shape cannot be clearly distinguished. Two potential applications of the X-shape in previous literature are tested, i.e., bulge rotation and Galactic center measurements. We find that the proper motion difference between the two sides of the X-shape can be used to estimate the mean azimuthal streaming motion of the bulge, but not the pattern speed of the bar. We also demonstrate that the Galactic center can be located with the X-shape, but the accuracy depends on the fitting scheme, the number of fields, and their latitudinal coverage.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the diffuse soft X-ray emission from the nuclear region of M51

We present an analysis of the diffuse soft X-ray emission from the nuclear region of M51 combining both XMM-Newton RGS and Chandra data. Most of the RGS spectrum of M51 can be fitted with a thermal model with a temperature of $\sim0.5$ keV except for the OVII triplet, which is forbidden-line dominated. The Fe L-shell lines peak around the southern cloud, where the OVIII and NVII Lya lines also peak. In contrast, the peak of the OVII forbidden line is about 10$"$ offset from that of the other lines, indicating that it is from a spatially distinct component. The spatial distribution of the OVII triplet mapped by the Chandra data shows that most of the OVII triplet flux is located at faint regions near edges, instead of the southern cloud where other lines peak. This distribution of the OVII triplet is inconsistent with the photoionization model. Other mechanisms that could produce the anomalous OVII triplet, including a recombining plasma and charge exchange X-ray emission, are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Are high velocity peaks in the Milky Way bulge due to the bar?

Recently the commissioning APOGEE observations of the Galactic bulge reported that a significant fraction of stars ($\sim10%$) are in a cold ($σ_{\rm V} \approx 30$ km/s) high velocity peak (Galactocentric radial velocity $\approx 200$ km/s). These stars are speculated to reflect the stellar orbits in the Galactic bar. In this study, we use two $N$-body models of the Milky Way-like disk galaxy with different bar strengths to critically examine this possibility. The general trends of the Galactocentric radial velocity distribution in observations and simulations are similar, but neither our models nor the BRAVA data reveal a statistically significant cold high velocity peak. A Monte Carlo test further suggests that it is possible for a spurious high velocity peak to appear if there are only a limited number of stars observed. Thus, the reported cold high velocity peak, even if it is real, is unlikely due to stars on the bar-supporting orbits. Our models do predict an excess of stars with high radial velocity, but not in a distinct peak. In the distance--velocity diagram, the high velocity particles in different fields exist at a similar distance $\sim8.5 \pm 1$ kpc away from the Sun. This result may be explained with geometric intersections between the line-of-sight and the particle orbits; high velocity stars naturally exist approximately at the tangent point, without constituting a distinct peak. We further demonstrate that even without the presence of a bar structure, particle motions in an axisymmetric disk can also exhibit an excess of high velocity stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Empirical Study of Simulated Two-planet Microlensing Event

We undertake the first study of two-planet microlensing models recovered from simulations of microlensing events generated by realistic multi-planet systems in which 292 planetary events including 16 two-planet events were detected from 6690 simulated light curves. We find that when two planets are recovered, their parameters are usually close to those of the two planets in the system most responsible for the perturbations. However, in one of the 16 examples, the apparent mass of both detected planets was more than doubled by the unmodeled influence of a third, massive planet. This fraction is larger than, but statistically consistent with, the roughly 1.5% rate of serious mass errors due to unmodeled planetary companions for the 274 cases from the same simulation in which a single planet is recovered. We conjecture that an analogous effect due to unmodeled stellar companions may occur more frequently. For seven out of 23 cases in which two planets in the system would have been detected separately, only one planet was recovered because the perturbations due to the two planets had similar forms. This is a small fraction (7/274) of all recovered single-planet models, but almost a third of all events that might plausibly have led to two-planet models. Still, in these cases, the recovered planet tends to have parameters similar to one of the two real planets most responsible for the anomaly.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How well can cold-dark-matter substructures account for the observed lensing flux-ratio anomalies?

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Microlensing of Sub-parsec Massive Binary Black Holes in Lensed QSOs: Light Curves and Size-Wavelength Relation

Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circum-binary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared-optical-UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius vs. wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius-wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Large Synoptic Survey telescope (LSST) and Euclid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Predictions for Microlensing Planetary Events from Core Accretion Theory

We conduct the first microlensing simulation in the context of planet formation model. The planet population is taken from the Ida & Lin core accretion model for $0.3M_\odot$ stars. With $6690$ microlensing events, we find for a simplified Korea Microlensing Telescopes Network (KMTNet) the fraction of planetary events is $2.9\%$ , out of which $5.5\%$ show multiple-planet signatures. The number of super-Earths, super-Neptunes and super-Jupiters detected are expected to be almost equal. Our simulation shows that high-magnification events and massive planets are favored by planet detections, which is consistent with previous expectation. However, we notice that extremely high-magnification events are less sensitive to planets, which is possibly because the 10 min sampling of KMTNet is not intensive enough to capture the subtle anomalies that occur near the peak. This suggests that while KMTNet observations can be systematically analyzed without reference to any follow-up data, follow-up observations will be essential in extracting the full science potential of very high-magnification events. The uniformly high-cadence observations expected for KMTNet also result in $\sim 55\%$ of all detected planets being non-caustic-crossing, and more low-mass planets even down to Mars-mass being detected via planetary caustics. We also find that the distributions of orbital inclinations and planet mass ratios in multiple-planet events agree with the intrinsic distributions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The InfraRed Imaging Spectrograph (IRIS) for TMT: Overview of innovative science programs

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Next Generation Virgo cluster Survey. V. modelling the dynamics of M87 with the Made-to-Measure method

We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The velocity distribution in the solar neighbourhood from the LAMOST pilot survey

We use about 15,000 F/G nearby dwarf stars selected from the LAMOST pilot survey to map the U-V velocity distribution in the solar neighbourhood. An extreme deconvolution algorithm is applied to reconstruct an empirical multi-Gaussian model. In addition to the well known substructures, e.g., Sirius, Coma Berenices, Hyades-Pleiades over-densities, several new substructures are unveiled. A ripple-like structure from (U, V) = (-120, -5) to (103, -32)km/s is clearly seen in the U-V distribution. This structure seems associated with resonance induced by the Galactic bar, since it is extended in U while having a small dispersion in V at the same time. A ridge structure between (U, V) = (-60, 40) and (-15, 15) km/s is also found. Although similar substructures have been seen in the Hipparcos data, their origin is still unclear. Another compact over-density is seen at (U, V) = (-102, -24). With this large data sample, we find that the substructure located at V~70 km/s and the Arcturus group are essentially parallel in V, which may indicate that they originate from an unrelaxed disk component perturbed by the rotating bar.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A New Photometric Model of the Galactic Bar using Red Clump Giants

We present a study of the luminosity density distribution of the Galactic bar using number counts of red clump giants (RCGs) from the OGLE-III survey. The data were recently published by Nataf et al. (2013) for 9019 fields towards the bulge and have $2.94\times 10^6$ RC stars over a viewing area of $90.25 \,\textrm{deg}^2$. The data include the number counts, mean distance modulus ($μ$), dispersion in $μ$ and full error matrix, from which we fit the data with several tri-axial parametric models. We use the Markov Chain Monte Carlo (MCMC) method to explore the parameter space and find that the best-fit model is the $E_3$ model, with the distance to the GC is 8.13 kpc, the ratio of semi-major and semi-minor bar axis scale lengths in the Galactic plane $x_{0},y_{0}$, and vertical bar scale length $z_0$, is $x_0:y_0:z_0 \approx 1.00:0.43:0.40$ (close to being prolate). The scale length of the stellar density profile along the bar's major axis is $\sim$ 0.67 kpc and has an angle of $29.4^\circ$, slightly larger than the value obtained from a similar study based on OGLE-II data. The number of estimated RC stars within the field of view is $2.78 \times 10^6$, which is systematically lower than the observed value. We subtract the smooth parametric model from the observed counts and find that the residuals are consistent with the presence of an X-shaped structure in the Galactic centre, the excess to the estimated mass content is $\sim 5.8%$. We estimate the total mass of the bar is $\sim 1.8 \times 10^{10} M_\odot$. Our results can be used as a key ingredient to construct new density models of the Milky Way and will have implications on the predictions of the optical depth to gravitational microlensing and the patterns of hydrodynamical gas flow in the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Schwarzschild model of the Galactic bar with initial density from N-body simulations

Using the potential from N-body simulations, we construct the Galactic bar models with the Schwarzschild method. By varying the pattern speed and the position angle of the bar, we find that the best-fit bar model has pattern speed $Ω_{\rm p}=40\ \rm{km\ s^{-1}\ kpc^{-1}}$, and bar angle $θ_{\rm bar}=45^{\circ}$. $N$-body simulations show that the best-fit model is stable for more than 1.5 Gyrs. Combined with the results in Wang et al. (2012), we find that the bar angle and/or the pattern speed are not well constrained by BRAVA data in our Schwarzschild models. The proper motions predicted from our model are slightly larger than those observed in four fields. In the future, more kinematic data from the ground and space-based observations will enable us to refine our model of the Milky Way bar.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Biases in physical parameter estimates through differential lensing magnification

We study the lensing magnification effect on background galaxies. Differential magnification due to different magnifications of different source regions of a galaxy will change the lensed composite spectra. The derived properties of the background galaxies are therefore biased. For simplicity, we model galaxies as a superposition of an axis-symmetric bulge and a face-on disk in order to study the differential magnification effect on the composite spectra. We find that some properties derived from the spectra (e.g., velocity dispersion, star formation rate and metallicity) are modified. Depending on the relative positions of the source and the lens, the inferred results can be either over- or under-estimates of the true values. In general, for an extended source at strong lensing regions with high magnifications, the inferred physical parameters (e.g. metallicity) can be strongly biased. Therefore detailed lens modelling is necessary to obtain the true properties of the lensed galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Degeneracies in triple gravitational microlensing

We study microlensing light curves by a triple lens, in particular, by a primary star plus two planets. A four-fold degeneracy is confirmed in the light curves, similar to the close and wide degeneracy found in a double lens. Furthermore, we derive a set of equations for triple-lens in the external shear approximation. By using these external shear equations, we identify two kinds of continuous degeneracies which may confused double and triple lenses, i.e. the continuous external shear degeneracy among triple-lens systems and the double-triple lens degeneracy. These degeneracies are particularly important in high magnification events, and thus some caution needs to be applied when one infers the fraction of stars hosting multiple planets from microlensing. We study the dependence of the degeneracies on the lensing parameters (e.g., source trajectory) and give recipes about how the degeneracies should be explored with real data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of plasma on gravitational lensing

We study gravitational lensing when plasma surrounds the lens. An extra deflection angle is induced by the plasma in addition to the deflection generated by gravity. An inhomogeneous plasma distribution generates a greater effect than a homogeneous one, and may cause significant effects to be detected in low frequency radio observations (a few hundred MHz). In particular, the lensed image positions will be different between optical and radio observations. The change of position due to a plasma can reach a few tens of milli-arcsec, which is readily detectable. One can use the position difference in different frequencies to estimate the density of plasma in the lens. The magnification ratios between multiple images are mainly determined by other properties of the lens, and are only weakly affected by the plasma. More importantly, we find that the strong lensing time delay will be affected by the plasma. Estimation of the Hubble constant from the time delay in low radio frequency observation may be slightly biased due to plasma in the lens. Unfortunately, the ionosphere of the Earth strongly affects low frequency radio observations. Thus our ability to detect the effect depends on how well we are able to calibrate out the ionosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fe K lines in the nuclear region of M82

We study the spatial distribution of the Fe 6.4 and 6.7 keV lines in the nuclear region of M82 using the Chandra archival data with a total exposure time of 500 ks. The deep exposure provides a significant detection of the Fe 6.4 keV line. Both the Fe 6.4 and 6.7 keV lines are diffuse emissions with similar spatial extent, but their morphology do not exactly follow each other. Assuming a thermal collisional-ionization-equilibrium model, the fitted temperatures are around 5-6 keV and the Fe abundances are about 0.4-0.6 solar value. We also report the spectrum of a point source, which shows a strong Fe 6.7 keV line and is likely a supernova remnant or a superbubble. The fitted Fe abundance of the point source is 1.7 solar value. It implies that part of the iron may be depleted from the X-ray emitting gases as the predicted abundance is about 5 solar value assuming complete mixing. If this is a representative case of the Fe enrichment, a mild mass-loading of a factor of 3 will make the Fe abundance of the point source in agreement with that of the hot gas, which then implies that most of the hard X-ray continuum (2-8 keV) of M82 has a thermal origin. In addition, the Fe 6.4 keV line is consistent with the fluorescence emission irradiated by the hard photons from nuclear point sources.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational lensing effects on sub-millimetre galaxy counts

We study the effects on the number counts of sub-millimetre galaxies due to gravitational lensing. We explore the effects on the magnification cross section due to halo density profiles, ellipticity and cosmological parameter (the power-spectrum normalisation $σ_8$). We show that the ellipticity does not strongly affect the magnification cross section in gravitational lensing while the halo radial profiles do. Since the baryonic cooling effect is stronger in galaxies than clusters, galactic haloes are more concentrated. In light of this, a new scenario of two halo population model is explored where galaxies are modeled as a singular isothermal sphere profile and clusters as a Navarro, Frenk and White (NFW) profile. We find the transition mass between the two has modest effects on the lensing probability. The cosmological parameter $σ_8$ alters the abundance of haloes and therefore affects our results. Compared with other methods, our model is simpler and more realistic. The conclusions of previous works is confirm that gravitational lensing is a natural explanation for the number count excess at the bright end.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Serendipitous discovery of a massive cD galaxy at z=1.096: Implications for the early formation and late evolution of cD galaxies

We have made a serendipitous discovery of a massive cD galaxy at z=1.096 in a candidate rich cluster in the HUDF area of GOODS-South. This brightest cluster galaxy is the most distant cD galaxy confirmed to date. Ultra-deep HST/WFC3 images reveal an extended envelope starting from ~10 kpc and reaching ~70 kpc in radius along the semi-major axis. The spectral energy distributions indicate that both its inner component and outer envelope are composed of an old, passively-evolving stellar population. The cD galaxy lies on the same mass-size relation as the bulk of quiescent galaxies at similar redshifts. The cD galaxy has a higher stellar mass surface density but a similar velocity dispersion to those of more-massive, nearby cDs. If the cD galaxy is one of the progenitors of today's more massive cDs, its size and stellar mass have had to increase on average by factors of $3.4\pm1.1$ and $3.3\pm1.3$ over the past ~8 Gyrs, respectively. Such increases in size and stellar mass without being accompanied by significant increases in velocity dispersion are consistent with evolutionary scenarios driven by both major and minor dry mergers. If such cD envelopes originate from dry mergers, our discovery of even one example proves that some BCGs entered the dry merger phase at epochs earlier than z=1. Our data match theoretical models which predict that the continuance of dry mergers at z<1 can result in structures similar to those of massive cD galaxies seen today. Moreover, our discovery is a surprise given that the extreme depth of the HUDF is essential to reveal such an extended cD envelope at z>1 and, yet, the HUDF covers only a minuscule region of sky. Adding that cDs are rare, Our serendipitous discovery hints that such cDs may be more common than expected. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-dimensional Microlensing

We study three-dimensional microlensing where two lenses are located at different distances along the line of sight. We formulate the lens equation in complex notations and recover several previous results. There are in total either 4 or 6 images, with an equal number of images with positive and negative parities. We find that the sum of signed magnifications for six image configurations is unity. Furthermore, we show that the light curves can be qualitatively different from those for binary lensing in a single plane. In particular, the magnifications between a `U'-shaped caustic crossing can be close to unity, rather than having a minimum magnification of 3 as in the single plane binary lensing. There is only a small probability three-dimensional microlensing events will be seen in microlensing toward the Galactic centre. It is more likely they will be first seen in cosmological microlensing.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Model for the Milky Way Bar

We use Schwarzschild's orbit-superposition technique to construct self-consistent models of the Galactic bar. Using $χ^2$ minimisation, we find that the best-fit Galactic bar model has a pattern speed $Ω_{\rm p}=60 \rm{km s^{-1} kpc^{-1}}$, disk mass $\rm{M_{\rm d}=1.0\times10^{11}M_{\odot}}$ and bar angle $θ_{\rm bar}=20^{\circ}$ for an adopted bar mass $\rm{M_{\rm bar}=2\times10^{10}M_{\odot}}$. The model can reproduce not only the three-dimensional and projected density distributions but also velocity and velocity dispersion data from the BRAVA survey. We also predict the proper motions in the range $l=[-12^{\circ},12^{\circ}]$, $b=[-10^{\circ},10^{\circ}]$, which appear to be higher than observations in the longitudinal direction. The model is stable within a timescale of 0.5 Gyr, but appears to deviate from steady-state on longer timescales. Our model can be further tested by future observations such as those from GAIA.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Astrophysical Applications of Gravitational Microlensing

Since the first discovery of microlensing events nearly two decades ago, gravitational microlensing has accumulated tens of TBytes of data and developed into a powerful astrophysical technique with diverse applications. The review starts with a theoretical overview of the field and then proceeds to discuss the scientific highlights. (1) Microlensing observations toward the Magellanic Clouds rule out the Milky Way halo being dominated by MAssive Compact Halo Objects (MACHOs). This confirms most dark matter is non-baryonic, consistent with other observations. (2) Microlensing has discovered about 20 extrasolar planets (16 published), including the first two Jupiter-Saturn like systems and the only "cold Neptunes" yet detected. They probe a different part of the parameter space and will likely provide the most stringent test of core accretion theory of planet formation. (3) Microlensing provides a unique way to measure the mass of isolated stars, including brown dwarfs to normal stars. Half a dozen or so stellar mass black hole candidates have also been proposed. (4) High-resolution, target-of-opportunity spectra of highly-magnified dwarf stars provide intriguing "age" determinations which may either hint at enhanced helium enrichment or unusual bulge formation theories. (5) Microlensing also measured limb-darkening profiles for close to ten giant stars, which challenges stellar atmosphere models. (6) Data from surveys also provide strong constraints on the geometry and kinematics of the Milky Way bar (through proper motions); the latter indicates predictions from current models appear to be too anisotropic compared with observations. The future of microlensing is bright given the new capabilities of current surveys and forthcoming new telescope networks from the ground and from space. Some open issues in the field are identified and briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Determining neutron star masses with weak microlensing

The masses of stars including stellar remnants are almost exclusively known from binary systems. In this work, we study gravitational microlensing of faint background galaxies by isolated neutron stars (pulsars). We show that the resulting surface brightness distortions can be used to determine the masses of neutron star. Due to different evolutionary histories, isolated neutron stars may have different masses from those in binary systems, and thus provide unique insight into their equation of states under extreme physical conditions. We search for existing pulsar catalogs and find one promising pair of a nearby pulsar and a background galaxy. This method will become more practical for the next generation optical and radio surveys and telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Internal kinematics of groups of galaxies in the Sloan Digital Sky Survey data release 7

We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final data release of the Sloan Digital Sky Survey (SDSS). For groups of given mass we estimate the redshift-space cross-correlation function (CCF) with respect to a reference galaxy sample, xi(r_p, pi), the projected CCF, w_p(r_p), and the real-space CCF, xi(r). The VDP is then extracted from the redshift distortion in xi(r_p, pi), by comparing xi(r_p, pi) with xi(r). We find that the velocity dispersion (VD) within virial radius (R_200) shows a roughly flat profile, with a slight increase at radii below ~0.3 R_200 for high mass systems. The average VD within the virial radius, sigma_v, is a strongly increasing function of central galaxy mass. We apply the same methodology to N-body simulations with the concordance Lambda cold dark matter cosmology but different values of the density fluctuation parameter sigma_8, and we compare the results to the SDSS results. We show that the sigma_v-M_* relation from the data provides stringent constraints on both sigma_8 and sigma_ms, the dispersion in log M_* of central galaxies at fixed halo mass. Our best-fitting model suggests sigma_8 = 0.86 +/- 0.03 and sigma_ms = 0.16 +/- 0.03. The slightly higher value of sigma_8 compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our simulations. Our VD measurements also provide a direct measure of the dark matter halo mass for central galaxies of different luminosities and masses, in good agreement with the results obtained by Mandelbaum et al. (2006) from stacking the gravitational lensing signals of the SDSS galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Made-to-measure galaxy models - II Elliptical and Lenticular Galaxies

We take a sample of 24 elliptical and lenticular galaxies previously analysed by the SAURON project using three-integral dynamical models created with Schwarzschild's method, and re-analyse them using the made-to-measure (M2M) method of dynamical modelling. We obtain good agreement between the two methods in determining the dynamical mass-to-light (M/L) ratios for the galaxies with over 80% of ratios differing by < 10% and over 95% differing by < 20%. We show that (M/L)_M2M is approximately equal to (M/L)_Sch. For the global velocity dispersion anisotropy parameter delta, we find similar values but with fewer of the made-to-measure models tangentially anisotropic by comparison with their SAURON Schwarzschild counterparts. Our investigation is the largest comparative application of the made-to-measure method to date.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Made-to-measure galaxy models - III Modelling with Milky Way observations

We demonstrate how the Syer & Tremaine made-to-measure method of stellar dynamical modelling can be adapted to model a rotating galactic bar. We validate our made-to-measure changes using observations constructed from the existing Shen et al. (2010) N-body model of the Milky Way bar, together with kinematic observations of the Milky Way bulge and bar taken by the Bulge Radial Velocity Assay (BRAVA). Our results for a combined determination of the bar angle and bar pattern speed (~30 degrees and ~40 km/s/kpc) are consistent with those determined by the N-body model. Whilst the made-to-measure techniques we have developed are applied using a particular N-body model and observational data set, they are in fact general and could be applied to other Milky Way modelling scenarios utilising different N-body models and data sets. Additionally, we use the exercise as a vehicle for illustrating how N-body and made-to-measure methods might be combined into a more effective method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Effects of Line-of-Sight Structures on Lensing Flux-ratio Anomalies in a LCDM Universe

The flux-ratio anomalies observed in multiply-lensed quasar images are most plausibly explained as the result of perturbing structures superposed on the underlying smooth matter distribution of the primary lens. The cold dark matter cosmological model predicts that a large number of substructures should survive inside larger halos but, surprisingly, this population alone has been shown to be insufficient to explain the observed distribution of the flux ratios of quasar's multiple images. Other halos (and their own subhalos) projected along the line of sight to the primary lens have been considered as additional source of perturbation. In this work, we use ray tracing through the Millennium II simulation to investigate the importance of projection effects due to halos and subhalos of mass m>1E8 Msun/h and extend our analysis to lower masses, m>1E6 Msun/h, using Monte-Carlo halo distributions. We find that the magnitude of the violation depends strongly on the density profile and concentration of the intervening halos, but clustering plays only a minor role. For a typical lensing geometry (lens at redshift 0.6 and source at redshift 2), background haloes (behind the main lens) are more likely to cause a violation than foreground halos. We conclude that line-of-sight structures can be as important as intrinsic substructures in causing flux-ratio anomalies. The combined effect of perturbing structures within the lens and along the line of sight in the LCDM universe results in a cusp-violation probability of 20-30%. This alleviates the discrepancy between models and current data, but a larger observational sample is required for a stronger test of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the soft X-ray emission of M82

We present a spatial analysis of the soft X-ray and Hα emissions from the outflow of the starburst galaxy M82. We find that the two emissions are tightly correlated on various scales. The O VII triplet of M82, as resolved by X-ray grating observations of XMM-Newton, is dominated by the forbidden line, inconsistent with the thermal prediction. The O VII triplet also shows some spatial variations. We discuss three possible explanations for the observed O VII triplet, including the charge exchange at interfaces between the hot outflow and neutral cool gas, a collisional non-equilibrium-ionization recombining plasma, and resonance scattering.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Star formation activities in early-type brightest cluster galaxies

We identify a total of 120 early-type Brightest Cluster Galaxies (BCGs) at 0.1<z<0.4 in two recent large cluster catalogues selected from the Sloan Digital Sky Survey (SDSS). They are selected with strong emission lines in their optical spectra, with both Hα and [O II]λ3727 line emission, which indicates significant ongoing star formation. They constitute about ~ 0.5% of the largest, optically-selected, low-redshift BCG sample, and the fraction is a strong function of cluster richness. Their star formation history can be well described by a recent minor and short starburst superimposed on an old stellar component, with the recent episode of star formation contributing on average only less than 1 percent of the total stellar mass. We show that the more massive star-forming BCGs in richer clusters tend to have higher star formation rate (SFR) and specific SFR (SFR per unit galaxy stellar mass). We also compare their statistical properties with a control sample selected from X-ray luminous clusters, and show that the fraction of star-forming BCGs in X-ray luminous clusters is almost one order of magnitude larger than that in optically-selected clusters. BCGs with star formation in cooling flow clusters usually have very flat optical spectra and show the most active star formation, which may be connected with cooling flows.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A distortion of very--high--redshift galaxy number counts by gravitational lensing

The observed number counts of high-redshift galaxy candidates have been used to build up a statistical description of star-forming activity at redshift z >~ 7, when galaxies reionized the Universe. Standard models predict that a high incidence of gravitational lensing will probably distort measurements of flux and number of these earliest galaxies. The raw probability of this happening has been estimated to be ~ 0.5 percent, but can be larger owing to observational biases. Here we report that gravitational lensing is likely to dominate the observed properties of galaxies with redshifts of z >~ 12, when the instrumental limiting magnitude is expected to be brighter than the characteristic magnitude of the galaxy sample. The number counts could be modified by an order of magnitude, with most galaxies being part of multiply imaged systems, located less than 1 arcsec from brighter foreground galaxies at z ~ 2. This lens-induced association of high-redshift and foreground galaxies has perhaps already been observed among a sample of galaxy candidates identified at z ~ 10.6. Future surveys will need to be designed to account for a significant gravitational lensing bias in high-redshift galaxy samples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charge Exchange X-ray Emission of M82: Kα triplets of O VII, Ne IX, and Mg XI

Starburst galaxies are primary feedback sources of mechanical energy and metals, which are generally measured from associated X-ray emission lines assuming that they are from the thermal emission of the outflowing hot gas. Such line emission, however, can also arise from the charge exchange X-ray emission (CXE) between highly ionized ions and neutral species. To understand the feedback of energy and metals, it is crucial to determine the origin of the X-ray emission lines and to distinguish the contributions from the CXE and the thermal emission. The origin of the lines can be diagnosed by the Kα triplets of He-like ions, because the CXE favors the inter-combination and forbidden lines, while the thermal emission favors the resonance line. We analyze the triplets of O VII, Ne IX, and Mg XI observed in the XMM- Newton reflection grating spectra of the starburst galaxy M82. The flux contribution of the CXE is 90%, 50%, and 30% to the O VII, Ne IX, and Mg XI triplet, respectively. Averaged over all the three triplets, the contribution of the CXE is \sim 50% of the total observed triplet flux. To correctly understand the hot outflow of starburst galaxies, it is necessary to include the CXE. Based on the measured CXE contributions to the O VII, Ne IX, and Mg XI triplets, we estimate the relative abundances of O, Ne, and Mg of the outflow and find they are similar to the solar ratios.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Charge Exchange X-ray Emission of Nearby Star-forming Galaxies

Properties of hot gas outflows from galaxies are generally measured from associated X-ray line emission assuming that it represents atomic transitions in thermally excited hot gas. X-ray line emission, however, can also arise from the charge exchange between highly ionised ions and neutral species. The Kα triplet of He-like ions can be used as a powerful diagnostic, because the charge exchange X-ray emission (CXE) favours the inter-combination and forbidden lines, while the thermal emission favours the resonance line. We analyse the OVII triplet of a sample of nine nearby star-forming galaxies observed by the XMM-Newton reflection grating spectrometers. For most galaxies, the forbidden lines are comparable to or stronger than the resonance lines, which is in contrast to the thermal prediction. For NGC 253, M51, M83, M61, NGC 4631, and the Antennae galaxy, the observed line ratios are consistent with the ratio of the CXE; for M94 and NGC 2903, the observed ratios indicate multiple origins; for M82, different regions show different line ratios, also indicating multiple origins. We discuss other possible mechanisms that can produce a relatively strong forbidden line, such as a collisional non-equilibrium-ionization recombining/ionizing plasma, which are not favoured. These results suggest that the CXE may be a common process and contribute a significant fraction of the soft X-ray line emission for galaxies with massive star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Effects of supermassive binary black holes on gravitational lenses

Recent observations indicate that many if not all galaxies host massive central black holes (BHs). In this paper we explore the influence of supermassive binary black holes (SMBBHs) on their actions as gravitational lenses. When lenses are modelled as singular isothermal ellipsoids, binary black holes change the critical curves and caustics differently as a function of distance. Each black hole can in principle create at least one additional image, which, if observed, provides evidence of black holes. By studying how SMBBHs affect the cumulative distribution of magnification for images created by black holes, we find that the cross section for at least one such additional image to have a magnification larger than $10^{-5}$ is comparable to the cross section for producing multiple-images in singular isothermal lenses. Such additional images may be detectable with high-resolution and large dynamic range maps of multiply-imaged systems from future facilities, such as the SKA. The probability of detecting at least one image (two images) with magnification above $10^{-3}$ is $\sim 0.2 \fBH$ ($\sim 0.05 \fBH$) in a multiply-imaged lens system, where $\fBH$ is the fraction of galaxies housing binary black holes. We also study the effects of SMBBHs on the core images when galaxies have shallower central density profiles (modelled as non-singular isothermal ellipsoids). We find that the cross section of the usually faint core images is further suppressed by SMBBHs. Thus their presence should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Estimate of halo ellipticity as a function of radius with flexions

The cold dark matter theory predicts triaxial dark matter haloes. The radial distribution of halo ellipticity depends on baryonic processes and the nature of dark matter particles (collisionless or collisional). Here we show that we can use lensing flexion ratios to measure the halo ellipticity as a function of radius. We introduce a weight function and study the relationship between the first and second order statistics of flexion ratios, both of which can be used to reduce the bias in the estimate of ellipticity. we perform numerical tests for our method, and demonstrate that it can reduce the bias and determine the halo ellipticity as a function of radius. We also point out that the minimum mean flexion ratio can be used to trace the centres of galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lensing by a singular isothermal sphere and a black hole

Most galaxies host central supermassive black holes. As two galaxies merge, the black holes also merge. The final single black hole may suffer a kick due to asymmetric gravitational radiation and may not be at the centre of the galaxy; off-centre black holes may also be produced by other means such as sustained acceleration due to asymmetric jet power. We model the main galaxy as a singular isothermal sphere and the black hole as an off-centre point lens, and study the critical curves and caustics using complex notation. We identify the critical parameters that govern the transitions in the topology of critical curves, caustics and pseudo-caustics, and find the number of images can be two, three, four and five. We show examples of image configurations, including cases where three highly de-magnified images are found close to the centre. The perturbation on the image magnification due to the black hole scales linearly with its mass in the off-centre case, and quadratically when the black hole is at the centre. Such images are difficult to observe unless high-contrast and high-resolution imaging facilities (e.g., the Square Kilometer Array in the radio) become available.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measurement of halo properties with weak lensing shear and flexion

We constrain properties of cluster haloes by performing likelihood analysis using lensing shear and flexion data. We test our analysis using two mock cluster haloes: an isothermal ellipsoid (SIE) model and a more realistic elliptical Navarro-Frenk-White (eNFW) model. For both haloes, we find that flexion is more sensitive to the halo ellipticity than shear. The introduction of flexion information significantly improves the constraints on halo ellipticity, orientation and mass. We also point out that there is a degeneracy between the mass and the ellipticity of SIE models in the lensing signal.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rapidly Rotating Lenses: Repeating features in the lightcurves of short period binary microlenses

Microlensing is most sensitive to binary lenses with relatively large orbital separations, and as such, typical binary microlensing events show little or no orbital motion during the event. However, despite the strength of binary microlensing features falling off rapidly as the lens separation decreases, we show that it is possible to detect repeating features in the lightcurve of binary microlenses that complete several orbits during the microlensing event. We investigate the lightcurve features of such Rapidly Rotating Lens (RRL) events. We derive analytical limits on the range of parameters where these effects are detectable, and confirm these numerically. Using a population synthesis Galactic model we estimate the RRL event rate for a ground-based and space-based microlensing survey to be 0.32fb and 7.8fb events per year respectively, assuming year-round monitoring and where fb is the binary fraction. We detail how RRL event parameters can be quickly estimated from their lightcurves, and suggest a method to model RRL events using timing measurements of lightcurve features. Modelling RRL lightcurves will yield the lens orbital period and possibly measurements of all orbital elements including the inclination and eccentricity. Measurement of the period from the lightcurve allows a mass-distance relation to be defined, which when combined with a measurement of microlens parallax or finite source effects, can yield a mass measurement to a two-fold degeneracy. With sub-percent accuracy photometry it is possible to detect planetary companions, but the likelihood of this is very small.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Detectability of Orbital Motion in Stellar Binary and Planetary Microlenses

A standard binary microlensing event lightcurve allows just two parameters of the lensing system to be measured: the mass ratio of the companion to its host, and the projected separation of the components in units of the Einstein radius. However, other exotic effects can provide more information about the lensing system. Orbital motion in the lens is one such effect, which if detected, can be used to constrain the physical properties of the lens. To determine the fraction of binary lens lightcurves affected by orbital motion (the detection efficiency) we simulate lightcurves of orbiting binary star and star-planet (planetary) lenses and simulate the continuous, high-cadence photometric monitoring that will be conducted by the next generation of microlensing surveys that are beginning to enter operation. The effect of orbital motion is measured by fitting simulated lightcurve data with standard static binary microlensing models; lightcurves that are poorly fit by these models are considered to be detections of orbital motion. We correct for systematic false positive detections by also fitting the lightcurves of static binary lenses. For a continuous monitoring survey without intensive follow-up of high magnification events, we find the orbital motion detection efficiency for planetary events with caustic crossings to be 0.061+-0.010, consistent with observational results, and 0.0130+-0.0055 for events without caustic crossings (smooth events). Similarly for stellar binaries, the orbital motion detection efficiency is 0.098+-0.011 for events with caustic crossings and is 0.048+-0.006 for smooth events. These result in combined (caustic crossing and smooth) orbital motion detection efficiencies of 0.029+-0.005 for planetary lenses and 0.070+-0.006 for stellar binary lenses. We also investigate how various microlensing parameters affect the orbital motion detectability. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Major Dry-Mergers In Early-Type Brightest Cluster Galaxies

We search for ongoing major dry-mergers in a well selected sample of local Brightest Cluster Galaxies (BCGs) from the C4 cluster catalogue. 18 out of 515 early-type BCGs with redshift between 0.03 and 0.12 are found to be in major dry-mergers, which are selected as pairs (or triples) with $r$-band magnitude difference $\dm<1.5$ and projected separation $\rp<30$ kpc, and showing signatures of interaction in the form of significant asymmetry in residual images. We find that the fraction of BCGs in major dry-mergers increases with the richness of the clusters, consistent with the fact that richer clusters usually have more massive (or luminous) BCGs. We estimate that present-day early-type BCGs may have experienced on average $\sim 0.6 (\tmerge/0.3\Gyr)^{-1}$ major dry-mergers and through this process increases their luminosity (mass) by $15% (\tmerge/0.3\Gyr)^{-1} (\fmass/0.5)$ on average since $z=0.7$, where $\tmerge$ is the merging timescale and $\fmass$ is the mean mass fraction of companion galaxies added to the central ones. We also find that major dry-mergers do not seem to elevate radio activities in BCGs. Our study shows that major dry-mergers involving BCGs in clusters of galaxies are not rare in the local Universe, and they are an important channel for the formation and evolution of BCGs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The clustering of barred galaxies in the local Universe

We study the clustering properties of barred galaxies using data from the Sloan Digital Sky Survey (SDSS). We compute projected redshift-space two-point cross-correlation functions w_p(r_p) for a sample of nearly 1000 galaxies for which we have performed detailed structural decompositions using the methods described in Gadotti (2009). The sample includes 286 barred galaxies. The clustering of barred and unbarred galaxies of similar stellar mass is indistinguishable over all the scales probed (~20kpc - 30Mpc). This result also holds even if the sample is restricted to bars with bluer g-i colours (and hence younger ages). Our result also does not change if we split our sample of barred galaxies according to bar-to-total luminosity ratio, bar boxyness, effective surface brightness, length, or the shape of the surface density profile within the bar. There is a hint that red, elliptical bars are more strongly clustered than red and less elliptical bars, on scales >~ 1 Mpc, although the statistical significance is not high. We conclude that there is no significant evidence that bars are a product of mergers or interactions. We tentatively interpret the stronger clustering of the more elliptical bars as evidence that they are located in older galaxies, which reside in more massive haloes.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Rotation in Gravitational Lenses

Gravitational lensing deflects light. A single lens deflector can only shear images, but cannot induce rotations. Multiple lens planes can induce rotations. Such rotations can be observed in quadruply imaged sources, and can be used to distinguish between two proposed solutions of the flux anomaly problem: substructures in lensing galaxies vs large scale structure. We predict the expected amount of rotation due to large scale structure in strong lensing systems, and show how this effect can be measured using ~ mas VLBI astrometry of quadruple lenses with extended source structures. The magnitude of rotation is around one degree. The biggest theoretical uncertainty is the power spectrum of dark matter on very small scales. This procedure can potentially be turned around to measure the dark matter power spectrum on very small scales. We list the predicted RMS rotation angles for several quadruple lenses with known lens and source redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Galaxy formation in pre-processed dark halos

Recent N-body simulations show that the formation of a present-day, galaxy sized dark matter halo in the CDM cosmogony in general consists of an early fast collapse phase, during which the potential associated with a halo is established, followed by a slow accretion phase, during which mass is added rather gently in the outer region. Here we outline a scenario in which the fast collapse phase is accompanied with rapid formation of cold clouds and with starbursts that can eject a large amount of gas from the halo center.The loss of orbital energy of the cold clouds to the dark matter and the ejection of gas from halo center by starburst can significantly reduce the halo concentration. The outflow from the starburst can also heat the gas in the protogalaxy region. Subsequent formation of galaxies in the slow accretion regime is therefore in halos that have been pre-processed by these processes and may have properties different from that given by N-body simulations. This scenario can help to solve several outstanding problems in the standard Lambda CDM model of galaxy formation without compromising its success in allowing structure formation at high redshift. The predicted rotation curves are quite flat,alleviating the discrepancy of the Tully-Fisher relation predicted in the Lambda CDM model with observations. The flattened galaxy halos allow accreted minihalos to survive in their central regions longer, which may be helpful in producing the flux anomalies observed in some gravitational lensing systems. The preheating by the early starbursts effectively reduces the amount of gas that can be accreted into galaxy halos, which may explain why the baryon fraction in a spiral galaxy is in general much lower than the universal baryon fraction, f_B~0.16, in the standard Lambda CDM model.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

An analytic model for the galactic winds and mass outflows

Galactic winds and mass outflows are observed both in nearby starburst galaxies and in high-redshift star-forming galaxies. In this paper we develop a simple analytic model to understand the observed superwind phenomenon. Our model is built upon the model of McKee & Ostriker (1977) for the interstellar medium. It allows one to predict how properties of a superwind, such as wind velocity and mass outflow rate, are related to properties of its star-forming host galaxy, such as size, gas density and star formation rate. The model predicts a threshold of star formation rate density for the generation of observable galactic winds. Galaxies with more concentrated star formation produce superwinds with higher velocities. The predicted mass outflow rates are comparable to (or slightly larger than) the corresponding star formation rates. We apply our model to both local starburst galaxies and high-redshift Lyman break galaxies, and find its predictions to be in good agreement with current observations. Our model is simple, and so can be easily incorporated into numerical simulations and semi-analytical models of galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Tully-Fisher Relation and its Implications for Halo Density Profile and Self-interacting Dark Matter

We show that the Tully-Fisher relation observed for spiral galaxies can be explained in the current scenario of galaxy formation without invoking subtle assumptions, provided that galactic-sized dark haloes have shallow, core-like central profiles, with a core radius proportional to halo circular velocity. In such a system, both disk and halo contribute significantly to the maximum rotation of the disk, and the interaction between the disk and halo components acts to reduce the scatter in the Tully-Fisher relation. With model parameters chosen in plausible ranges, the model can well accommodate the zero-point, slope, scatter of the observed Tully-Fisher relation, as well as the large range of disk sizes. The model predicts that LSB disks obey a Tully-Fisher relation similar to that of normal disks, if disk mass-to-light ratio is properly taken into account. The halo profile required by the Tully-Fisher relation is as shallow as that required by the rotation curves of faint disks, but much shallower than that predicted by conventional CDM models. Our results cannot be explained by some of the recent proposals for resolving the conflict between conventional CDM models and the rotation-curve shapes of faint galaxies. If dark matter self-interaction (either scattering or annihilation) is responsible for the shallow profile, the observed Tully-Fisher relation requires the interaction cross section σ_X to satisfy <σ_{X}|v|>/m_{X}~10^{-16} cm^3/s/GeV.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Extended Source Effects in Astrometric Gravitational Microlensing

Extended source size effects have been detected in photometric monitoring of gravitational microlensing events. We study similar effects in the centroid motion of an extended source lensed by a point mass. We show that the centroid motion of a source with uniform surface brightness can be obtained analytically. For a source with circularly symmetric limb-darkening profile, the centroid motion can be expressed as a one-dimensional integral, which can be evaluated numerically. We found that when the impact parameter is comparable to the source radius, the centroid motion is significantly modified by the finite source size. In particular, when the impact parameter is smaller than the source radius, the trajectories become clover-leaf like. Such astrometric motions can be detected using space interferometers such as the Space Interferometry Mission. Such measurements offer exciting possibilities of determining stellar parameters such as stellar radius to excellent accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Probing the Structure of Lensing Galaxies with Quadruple Lenses: The Effect of the ``External'' Shear

We study a general elliptical potential of the form $ψ(x^2+y^2/q^2)~ (0<q\le 1)$ plus an additional shear (with an arbitrary direction) as models for the observed quadruple lenses. It is shown that a minimum additional shear is needed even just to reproduce the observed positions alone. We also obtain the dependence of the axial ratio, $q$, on the orientation of the major axis of potential. A general relation also exists between the shear, the position angle and axial ratio of the lensing galaxy. The relation shows a generic degeneracy in modelling quadruple lenses. In particular, it shows that only the ratio of the ellipticity, $ε\equiv (1-q^2)/(1+q^2)$, to the magnitude of shear, $γ$ can be determined. All these results are valid regardless of the radial profile of the potential. Our formalism applies when the galaxy position is observed, which is the case for seven of the eight known quadruple lenses. Application to these seven cases reveals two quadruple lenses CLASS 1608+656 and HST 12531-2914, requiring highly significant shear with magnitude $\approx 0.2$. For HST 12531-2914, there must be a misalignment between the major axis of light and the major axis of potential (mass). We conclude that detailed modelling of quadruple lenses can yield valuable quantitative information about the shape of lensing galaxies and their dark matter halos.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

On the Microlensing Optical Depth of the Galactic Bar

The microlensing probability (optical depth $τ$) toward the Galactic center carries information about the mass distribution of the Galactic bulge/bar, so can be used to constrain the very uncertain shape parameters of the bar. We find $tau$ depends on the bar mass, radial profile, angle, axis scale lengths and boxyness by a few simple analytical formulae, which shows: (1) $τ$ is proportional to the mass of the bar, $M$. (2) $τ$ falls along the minor axis with a strong gradient. (3) An oblate bulge can have more optical depth than a triaxial bar if the bar angle $α>45$ degress. (4) $τ$ is the largest if the angle $α$ and the axis ratio $y_0/x_0$ conspires so that $y_0/x_0=\tan α$. (5) At a fixed field on the minor axis but away from the center, boxy bars with a flat density profile tend to give a larger optical depth than ellipsoidal bars with a steep profile. (6) Main sequence sources should have a significantly lower (20-50\% lower) optical depth than red clump giants if main sequence stars are not observed as deep as the bright clump giants. An application to four COBE-constrained models (Dwek et al. 1994) shows most models produce optical depth $2σ$ lower than MACHO and OGLE observed values even with both a massive bar $2.8\times 10^{10} M_\odot$ and a full disk. The high $τ$ argues for a massive ($> 2\times 10^{10}M_\odot$) boxy bar with $y_0/x0\approx \tanα$ and $α<20$ deg and with a flat radial profile up to corotation.

0 citations0 reviewsNo code linkedOpen access