Search

Hot Jupiters should initially form at considerable distances from host stars and subsequently migrate towards inner regions, supported directly by transit timing variation (TTV). We report the TTV of K2-237b, using reproduced timings fitted from \textit{Kepler} K2 and \textit{TESS} data. The timings span from 2016 to 2021, leading to an observational baseline of 5 years. The timing evolution presents a significant bias to a constant period scenario. The model evidence is evaluated utilizing the Bayesian Information Criterion (BIC), which favours the scenario of period decay with a $\Delta$BIC of 14.1. The detected TTV induces a period decay rate ($\dot{P}$) of -1.14$\pm$0.28$\times$10$^{-8}$ days per day ($-$0.36 s/year). Fitting the spectral energy distribution, we find infrared excess at the significance level of 1.5 $\sigma$ for WISE W1 and W2 bands, and 2 $\sigma$ level for W3 and W4 bands. This potentially reveals the existence of a stellar disk, consisting of hot dust at 800$\pm$300 K, showing a $L_{dust}/L_{\ast}$ of 5$\pm$3$\times$10$^{-3}$. We obtain a stellar age of 1.0$^{+1.4}_{-0.7}$$\times$10$^{9}$ yr from isochrone fitting. The properties of K2-237b potentially serve as a direct observational support to the planet disk migration though more observation are needed., Comment: 5 pages, Accepted for publication in MNRAS Letters

Conventional galaxy mass estimation methods suffer from model assumptions and degeneracies. Machine learning, which reduces the reliance on such assumptions, can be used to determine how well present-day observations can yield predictions for the distributions of stellar and dark matter. In this work, we use a general sample of galaxies from the TNG100 simulation to investigate the ability of multi-branch convolutional neural network (CNN) based machine learning methods to predict the central (i.e., within $1-2$ effective radii) stellar and total masses, and the stellar mass-to-light ratio $M_*/L$. These models take galaxy images and spatially-resolved mean velocity and velocity dispersion maps as inputs. Such CNN-based models can in general break the degeneracy between baryonic and dark matter in the sense that the model can make reliable predictions on the individual contributions of each component. For example, with $r$-band images and two galaxy kinematic maps as inputs, our model predicting $M_*/L$ has a prediction uncertainty of 0.04 dex. Moreover, to investigate which (global) features significantly contribute to the correct predictions of the properties above, we utilize a gradient boosting machine. We find that galaxy luminosity dominates the prediction of all masses in the central regions, with stellar velocity dispersion coming next. We also investigate the main contributing features when predicting stellar and dark matter mass fractions ($f_*$, $f_{\rm DM}$) and the dark matter mass $M_{DM}$, and discuss the underlying astrophysics., Comment: 17 pages, 12 figures, published on MNRAS

We propose a new method for measuring the spatial density distribution of the stellar halo of the Milky Way. Our method is based on a pairwise statistic of the distribution of stars on the sky, the angular two-point correlation function (ATPCF). The ATPCF utilizes two dimensional data of stars only and is therefore immune to the large uncertainties in the determination of distances to stars. We test our method using mock stellar data coming from various models including the single power-law (SPL) and the broken power-law (BPL) density profiles. We also test the influence of axisymmetric flattening factors using both constant and varying values. We find that the ATPCF is a powerful tool for recovering the spatial distributions of the stellar halos in our models. We apply our method to observational data from the type ab RR Lyrae catalog in the Catalina Survey Data Release 1. In the 3-parameter BPL model, we find that $s_{1}=2.46_{-0.20}^{+0.18}, s_{2}=3.99_{-1.33}^{+0.75}$ and $r_{0}=31.11_{-5.88}^{+7.61}$, which are in good agreement with previous results. We also find that introducing an extra parameter, the radially varying flattening factor, greatly improves our ability to model accurately the observed data distribution. This implies perhaps that the stellar halo of the Milky Way should be regarded as oblate., Comment: Accepted for publication in MNRAS

Knowing the conserved quantities that a galaxy's stellar orbits conform to is important in helping us understand the stellar distribution and structures within the galaxy. Isolating integrals of motion and resonances are particularly important, non-isolating integrals less so. We compare the behavior and results of two methods for counting the number of conserved quantities, one based on the correlation integral approach and the other a more recent method using machine learning. Both methods use stellar orbit trajectories in phase space as their only input, and we create such trajectories from theoretical spherical, axisymmetric and triaxial model galaxies. The orbits have known isolating integrals and resonances. We find that neither method is fully effective in recovering the numbers of these quantities, nor in determining the number of non-isolating integrals. From a computer performance perspective, we find the correlation integral approach to be the faster. Determining the algebraic formulae of (multiple) conserved quantities from the trajectories has not been possible due to the lack of an appropriate symbolic regression capability. Notwithstanding the shortcomings we have noted, it may be that the methods are usable as part of a trajectory analysis tool kit., Comment: Accepted for publication by RAA, 18 pages, 9 figures, 4 tables

We combine the kinematics of 159 globular clusters (GCs) provided by the Gaia Early Data Release 3 (EDR3) with other observational data to classify the GCs, and to estimate the mass of the Milky Way (MW). We use the age-metallicity relation, integrals of motion, action space and the GC orbits to identify the GCs as either formed in-situ (Bulge and Disk) or ex situ (via accretion). We find that $45.3\%$ have formed in situ, $38.4\%$ may be related to known merger events: Gaia-Sausage-Enceladus, the Sagittarius dwarf galaxy, the Helmi streams, the Sequoia galaxy, and the Kraken galaxy. We also further identify three new sub-structures associated with the Gaia-Sausage-Enceladus. The remaining $16.3\%$ of GCs are unrelated to the known mergers and thought to be from small accretion events. We select 46 GCs which have radii $8.0

Single-line spectroscopic binaries recently contribute to the stellar-mass black hole discovery, independently of the X-ray transient method. We report the identification of a single-line binary system LTD064402+245919, with an orbital period of 14.50 days. The observed component is a subgiant with a mass of 2.77$\pm$0.68M$_{\odot}$, radius 15.5$\pm$2.5R$_{\odot}$, effective temperature $T_{\rm eff}$ 4500$\pm$200K, and surface gravity log\emph{g} 2.5$\pm$0.25dex. The discovery makes use of the LAMOST time-domain (LAMOST-TD) and ZTF survey. Our general-purpose software pipeline applies the Lomb-Scargle periodogram to determine the orbital period and uses machine-learning to classify the variable type from the folded light curves. We apply a combined model to estimate the orbital parameters from both the light and radial velocity curves, taking constraints on the primary star mass, mass function, and detection limit of secondary luminosity into consideration. We obtain a radial velocity semi-amplitude of 44.6$\pm$1.5 km s$^{-1}$, mass ratio of 0.73$\pm$0.07, and an undetected component mass of 2.02$\pm$0.49M$_{\odot}$ when the type of the undetected component is not set. We conclude that the inclination is not well constrained, and that the secondary mass is larger than 1M$_{\odot}$ when the undetected component is modelled as a compact object. According to our investigations using an MCMC simulation, increasing the spectra SNR by a factor of 3 would enable the secondary light to be distinguished (if present). The algorithm and software in this work are able to serve as general-purpose tools for the identification of compact objects quiescent in X-rays., Comment: Accepted by ApJ

The stellar mass-to-light ratio gradient in SDSS $r-$band $\nabla (M_*/L_r)$ of a galaxy depends on its mass assembly history, which is imprinted in its morphology and gradients of age, metallicity, and stellar initial mass function (IMF). Taking a MaNGA sample of 2051 galaxies with stellar masses ranging from $10^9$ to $10^{12}M_\odot$ released in SDSS DR15, we focus on face-on galaxies, without merger and bar signatures, and investigate the dependence of the 2D $\nabla (M_*/L_r)$ on other galaxy properties, including $M_*/L_r$-colour relationships by assuming a fixed Salpeter IMF as the mass normalization reference. The median gradient is $\nabla M_*/L_r\sim -0.1$ (i.e., the $M_*/L_r$ is larger at the centre) for massive galaxies, becomes flat around $M_*\sim 10^{10} M_{\odot}$ and change sign to $\nabla M_*/L_r\sim 0.1$ at the lowest masses. The $M_*/L_r$ inside a half light radius increases with increasing galaxy stellar mass; in each mass bin, early-type galaxies have the highest value, while pure-disk late-type galaxies have the smallest. Correlation analyses suggest that the mass-weighted stellar age is the dominant parameter influencing the $M_*/L_r$ profile, since a luminosity-weighted age is easily affected by star formation when the specific star formation rate (sSFR) inside the half light radius is higher than $10^{-3} {\rm Gyr}^{-1}$. With increased sSFR gradient, one can obtain a steeper negative $\nabla (M_*/L_r)$. The scatter in the slopes of $M_*/L$-colour relations increases with increasing sSFR, for example, the slope for post-starburst galaxies can be flattened to $0.45$ from the global value $0.87$ in the $M_*/L$ vs. $g-r$ diagram. Hence converting galaxy colours to $M_*/L$ should be done carefully, especially for those galaxies with young luminosity-weighted stellar ages, which can have quite different star formation histories., Comment: 12 pages, 10 figures. Accepted for publication in MNRAS

We present a novel, iterative method using an empirical Bayesian approach for modeling the limb darkened WASP-121b transit from the TESS light curve. Our method is motivated by the need to improve $R_{p}/R_{\ast}$ estimates for exoplanet atmosphere modeling, and is particularly effective with the limb darkening (LD) quadratic law requiring no prior central value from stellar atmospheric models. With the non-linear LD law, the method has all the advantages of not needing atmospheric models but does not converge. The iterative method gives a different $R_{p}/R_{\ast}$ for WASP-121b at a significance level of 1$\sigma$ when compared with existing non-iterative methods. To assess the origins and implications of this difference, we generate and analyze light curves with known values of the limb darkening coefficients (LDCs). We find that non-iterative modeling with LDC priors from stellar atmospheric models results in an inconsistent $R_{p}/R_{\ast}$ at 1.5$\sigma$ level when the known LDC values are as those previously found when modeling real data by the iterative method. In contrast, the LDC values from the iterative modeling yields the correct value of $R_{p}/R_{\ast}$ to within 0.25$\sigma$. For more general cases with different known inputs, Monte Carlo simulations show that the iterative method obtains unbiased LDCs and correct $R_{p}/R_{\ast}$ to within a significance level of 0.3$\sigma$. Biased LDC priors can cause biased LDC posteriors and lead to bias in the $R_{p}/R_{\ast}$ of up to 0.82$\%$, 2.5$\sigma$ for the quadratic law and 0.32$\%$, 1.0$\sigma$ for the non-linear law. Our improvement in $R_{p}/R_{\ast}$ estimation is important when analyzing exoplanet atmospheres., Comment: Accepted by AJ

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., Comment: 13 pages, 11 figures, accepted by MNRAS

From Rybicki's analysis using the Fourier slice theorem, mathematically it is possible to reproduce uniquely an edge-on axisymmetric galaxy's 3D light distribution from its 2D surface brightness. Utilizing galaxies from a cosmological simulation, we examine the ability of Syer and Tremaine's made-to-measure method and Schwarzschild's method for stellar dynamical modeling to do so for edge-on oblate axisymmetric galaxies. Overall, we find that the methods do not accurately recover the 3D distributions, with the made-to-measure method producing more accurate estimates than Schwarzschild's method. Our results have implications broader than just luminosity density, and affect other luminosity-weighted distributions within galaxies, for example, age and metallicity., Comment: Accepted for publication by Research in Astronomy and Astrophysics. 26 pages, 15 figures, 4 tables

Binaries play key roles in determining stellar parameters and exploring stellar evolution models. We build a catalog of 88 eclipsing binaries with spectroscopic information, taking advantage of observations from both the Large Sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) and the Palomar Transient Factory (PTF) surveys. A software pipeline is constructed to identify binary candidates by examining their light curves. The orbital periods of binaries are derived from the Lomb-Scargle method. The key distinguishing features of eclipsing binaries are recognized by a new filter \textit{Flat Test}. We classify the eclipsing binaries by applying Fourier analysis on the light curves. Among all the binary stars, 13 binaries are identified as eclipsing binaries for the first time. The catalog contains information: position, primary eclipsing magnitude and time, eclipsing depth, the number of photometry and radial velocity observations, largest radial velocity difference, binary type, the effective temperature of observable star $T_{\rm eff}$, and surface gravity of observable star log \emph{g}. The false-positive probability is calculated by using both a Monte Carlo simulation and real data from the SDSS Stripe 82 Standard Catalog. The binaries in the catalog are mostly with a period of less than one day. The period distribution shows a 0.22-day cut-off which is consistent with the low probability of an eclipsing binary rotating with such a period., Comment: 31 pages, 10 figures, In Press (accepted by ApJS)

We investigate the relationship between the optical variability of galaxies and their distances from the centre of the Virgo Cluster using Palomar Transient Factory data. We define the ratio between the standard deviation of the galaxy brightness and the mean value of the standard deviation as a measure of a galaxy's optical variability. A sample of 814 Virgo galaxies with 230263 observations shows a monotonically decreasing trend of optical variability with increasing clustercentric distance. The variability level inside the cluster is 3.2$\sigma$ higher than the level outside. We fit the variability with a linear function and find that the data reject a distance-independent model. We examine 217 background galaxies for comparison and find no significant trend in galaxy variability. We assess the relation with Monte Carlo simulation by rebuilding the brightness of each galaxy. The simulation shows a monotonically decreasing relation for member galaxy variability and a distance-independent relation for background galaxies. Our result is consistent with the theory that the cold gas flowing inwards the cluster centre fuels AGN activity. This work is a new implementation of the method using optical variability to investigate the relation between galaxies evolution and their environment., Comment: 5 pages, 4 figures, in Press (accepted by MNRAS Letters)

We enhance the Syer & Tremaine made-to-measure (M2M) particle method of stellar dynamical modelling to model simultaneously both kinematic data and absorption line strength data thus creating a `chemo-M2M' modelling scheme. We apply the enhanced method to four galaxies (NGC 1248, NGC 3838, NGC 4452, NGC 4551) observed using the SAURON integral-field spectrograph as part of the ATLAS3D programme. We are able to reproduce successfully the 2D line strength data achieving mean chi^2 per bin values of ~1 with >95\% of particles having converged weights. Because M2M uses a 3D particle system, we are also able to examine the underlying 3D line strength distributions. The extent to which these distributions are plausible representations of real galaxies requires further consideration. Overall we consider the modelling exercise to be a promising first step in developing a `chemo-M2M' modelling system and in understanding some of the issues to be addressed. Whilst the made-to-measure techniques developed have been applied to absorption line strength data, they are in fact general and may be of value in modelling other aspects of galaxies., Comment: 18 pages

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 $\Omega_{\rm p}=40\ \rm{km\ s^{-1}\ kpc^{-1}}$, and bar angle $\theta_{\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., Comment: 10 pages, 7 figures, 2 tables. Accepted for publication in MNRAS