Researcher profile

Shude Mao

Shude Mao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
39works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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}$&#39;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.

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&#39;&#39; 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.

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.

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.

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.

preprint2022arXiv

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

In order to exhume the buried signatures of &#34;missing planetary caustics&#34; 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}$.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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 &#34;dip&#34; 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 &#34;$\pm$&#34; 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.

preprint2020arXiv

A discrete chemo-dynamical model of M87&#39;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.

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.

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)

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&#39;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&#39;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&#39; 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.

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.

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

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.

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.

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.

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 &#34;elbow&#34; 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.

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&#39; 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&#39;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&#39;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.

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 &#34;skewness&#34; ($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.

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&#39;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 &#34;mass ratio desert&#34; at -3.7 \lesssim \log q \lesssim -3.0 for the 31 published KMTNet planets.

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.

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.

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&#39;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.

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.

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.