Showing 1–50 of 84 results for author: Ji, J

The Formation of Ultra-short-period Planets under the Influence of the Nearby Planetary Companions

Authors: Jia Jun Zhu, Su Wang, Jianghui Ji, Yao Dong

Abstract: Ultra-short-period (USP) planets, defined as those with orbital periods shorter than 1 day, provide valuable insights into planetary evolution under strong stellar tidal interactions. In this work, we investigate the formation of USP planets in two-planet systems consisting of an inner terrestrial planet accompanied by an outer hot Jupiter (HJ). Our simulation results show USP planets can form thr… ▽ More Ultra-short-period (USP) planets, defined as those with orbital periods shorter than 1 day, provide valuable insights into planetary evolution under strong stellar tidal interactions. In this work, we investigate the formation of USP planets in two-planet systems consisting of an inner terrestrial planet accompanied by an outer hot Jupiter (HJ). Our simulation results show USP planets can form through a process driven by secular perturbations from the outer companion, which induce eccentricity excitation, tidal dissipation, and subsequent orbital decay of the inner planet. The probability of USP formation is governed by key factors, including the mass ratio between two planets, their orbital eccentricities, and the tidal dissipation process. 6.7\% of our simulations form USP planets, and USP planets form most efficiently when the mass ratio is around 4 $M_{\oplus}{\rm /}M_{\rm J}$, with the inner planet less than 8 $M_{\oplus}$. Furthermore, the eccentricity of the outer HJ plays a crucial role-moderate eccentricities ($e_{\rm outer}<0.1$) favor USP formation, whereas higher eccentricities ($e_{\rm outer}>0.1$) enhance the likelihood of orbital instability, often resulting in a lonely HJ. USP planets form more efficiently when the tidal dissipation function of the inner planet is comparable to the values estimated for terrestrial planets in the solar system. Comparison with observed planetary systems reveals that systems with large mass ratios or nearly circular outer planets tend to produce short-period (SP) planets instead of USP planets. Our findings offer a potential explanation for the most commonly observed system architectures, which predominantly feature either an HJ with an inner SP planet or a lonely HJ. △ Less

Submitted 28 September, 2025; originally announced September 2025.

Comments: 15 pages, 8 figures, accepted for publication in ApJ

Direct Imaging for the Debris Disk around $ε$ Eridani with the Cool-Planet Imaging Coronagraph

Authors: Chunhui Bao, Jianghui Ji, Gang Zhao, Yiming Zhu, Jiangpei Dou, Su Wang, Yao Dong

Abstract: We analyze the inner debris disk around $ε$ Eridani using simulated observations with the Cool-Planet Imaging Coronagraph (CPI-C). Using the radiative transfer code MCFOST, we generate synthetic scattered-light images and spectral energy distributions for three disk models that differ in inclination and radial extent, and compare these results with the anticipated performance of CPI-C. CPI-C can r… ▽ More We analyze the inner debris disk around $ε$ Eridani using simulated observations with the Cool-Planet Imaging Coronagraph (CPI-C). Using the radiative transfer code MCFOST, we generate synthetic scattered-light images and spectral energy distributions for three disk models that differ in inclination and radial extent, and compare these results with the anticipated performance of CPI-C. CPI-C can resolve disk structures down to $\sim$3 au, offering substantially finer spatial resolution than existing HST/STIS and Spitzer/IRS observations. Recovered inclinations and radial extents closely match the input models, constraining the disk geometry and informing potential planet-disk interactions in the $ε$ Eri system. Although the cold Jupiter-like planet $ε$ Eri b is not detected in our simulations, polarimetric methods may enable detection of its reflected light. These results highlight the capability of next-generation coronagraphs to probe cold dust in nearby planetary systems. △ Less

Submitted 20 September, 2025; originally announced September 2025.

Comments: 17 pages, 9 figures, accepted for publication in Research in Astronomy and Astrophysics

Discovery and dynamics of a Sedna-like object with a perihelion of 66 au

Authors: Ying-Tung Chen, Patryk Sofia Lykawka, Yukun Huang, JJ Kavelaars, Wesley C. Fraser, Michele T. Bannister, Shiang-Yu Wang, Chan-Kao Chang, Matthew J. Lehner, Fumi Yoshida, Brett Gladman, Mike Alexandersen, Edward Ashton, Young-Jun Choi, A. Paula Granados Contreras, Takashi Ito, Youngmin JeongAhn, Jianghui Ji, Myung-Jin Kim, Samantha M. Lawler, Jian Li, Zhong-Yi Lin, Hong-Kyu Moon, Surhud More, Marco Muñoz-Gutiérrez , et al. (8 additional authors not shown)

Abstract: Trans-Neptunian objects (TNOs) with large perihelion distances ($q > 60$ au) and semi-major axes ($a > 200$ au) provide insights into the early evolution of the solar system and the existence of a hypothetical distant planet. These objects are still rare and their detection is challenging, yet they play a crucial role in constraining models of solar system formation. Here we report the discovery o… ▽ More Trans-Neptunian objects (TNOs) with large perihelion distances ($q > 60$ au) and semi-major axes ($a > 200$ au) provide insights into the early evolution of the solar system and the existence of a hypothetical distant planet. These objects are still rare and their detection is challenging, yet they play a crucial role in constraining models of solar system formation. Here we report the discovery of a Sedna-like TNO, 2023\,KQ$_{14}$, nicknamed `Ammonite', with $q = 66$ au, $a = 252$ au, and inclination $i=11^\circ$. Ammonite's orbit does not align with those of the other Sedna-like objects and fills the previously unexplained `$q$-gap' in the observed distribution of distant solar system objects. Simulations demonstrate that Ammonite is dynamically stable over 4.5 billion years. % with less than 1\% variation in its semi-major axis. Our analysis suggests that Ammonite and the other Sedna-like objects may have shared a primordial orbital clustering around 4.2 billion years ago. Furthermore, Ammonite's stable orbit favors larger orbits ($\sim$ 500 au) rather than closer ones for a large hypothetical planet in present-day trans-Neptunian space. △ Less

Submitted 4 August, 2025; originally announced August 2025.

Comments: Accepted manuscript of an article published open access in Nature Astronomy under a Creative Commons Attribution-NonCommercial-NoDerivatives (CC BY-NC-ND 4.0) license. The final published version is available at https://doi.org/10.1038/s41550-025-02595-7

Journal ref: Nature Astronomy, 2025

Introduction to the Chinese Space Station Survey Telescope (CSST)

Authors: CSST Collaboration, Yan Gong, Haitao Miao, Hu Zhan, Zhao-Yu Li, Jinyi Shangguan, Haining Li, Chao Liu, Xuefei Chen, Haibo Yuan, Jilin Zhou, Hui-Gen Liu, Cong Yu, Jianghui Ji, Zhaoxiang Qi, Jiacheng Liu, Zigao Dai, Xiaofeng Wang, Zhenya Zheng, Lei Hao, Jiangpei Dou, Yiping Ao, Zhenhui Lin, Kun Zhang, Wei Wang , et al. (97 additional authors not shown)

Abstract: 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 inst… ▽ More 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. △ Less

Submitted 19 September, 2025; v1 submitted 6 July, 2025; originally announced July 2025.

Comments: 48 pages, 12 figures, 1 table. Accepted for publication in SCIENCE CHINA Physics, Mechanics & Astronomy

A Be star-black hole binary with a wide orbit from LAMOST time-domain survey

Authors: Qian-Yu An, Yang Huang, Wei-Min Gu, Yong Shao, Zhi-Xiang Zhang, Tuan Yi, B. D. Lailey, T. A. A. Sigut, Kyle Akira Rocha, Meng Sun, Seth Gossage, Shi-Jie Gao, Shan-Shan Weng, Song Wang, Bowen Zhang, Xinlin Zhao, Senyu Qi, Shilong Liao, Jianghui Ji, Junfeng Wang, Jianfeng Wu, Mouyuan Sun, Xiang-Dong Li, Jifeng Liu

Abstract: Binary systems consisting of an early type star and a black hole (BH) are crucial for understanding various astrophysical phenomena, particularly the origins of detected gravitational wave sources. Be binary systems are expected to represent a key evolutionary stage in hosting BHs. However, while hundreds of Be X-ray binaries are known, the only confirmed BH candidate in a Be binary remains highly… ▽ More Binary systems consisting of an early type star and a black hole (BH) are crucial for understanding various astrophysical phenomena, particularly the origins of detected gravitational wave sources. Be binary systems are expected to represent a key evolutionary stage in hosting BHs. However, while hundreds of Be X-ray binaries are known, the only confirmed BH candidate in a Be binary remains highly controversial. We report the discovery of ALS 8814, a Be star-BH binary with a moderately eccentric ($e = 0.23$) and wide orbit ($P = 176.6$ days), revealed by the radial velocity (RV) measurement of the visible Be star. Our analysis, combining flux-calibrated spectra in the Balmer discontinuity region and spectral template matching, yields a mass of $11.2^{+1.4}_{-1.2}$ $M_\odot$ for the Be star. The minimum mass of the unseen companion, assuming an edge-on inclination ($i = 90^{\circ}$), is $9.8\pm 0.7\,M_\odot$. We rule out the presence of non-degenerate companions in ALS 8814, indicating that it can only be a BH. This discovery represents a robust case of a Be-BH binary, identified purely through precise RV measurements from a single set of lines. The extremely low peculiar velocity of ALS 8814 suggests that the BH is formed via a direct core-collapse with a negligible natal kick, implying an almost perfect alignment between the Be star's spin and the orbital plane. In this context, the binary's inclination angle is estimated to be 22$^{\circ}$-49$^{\circ}$ by analyzing the shallow double-peaked profile of the H$α$ emission line. This inclination range corresponds to a BH mass estimate between $15\,M_\odot$ and $58\,M_\odot$. As the only unambiguous Be-BH binary system known to date, ALS 8814 provides valuable constraints on the BH formation in a binary system with a high-mass companion. △ Less

Submitted 29 May, 2025; originally announced May 2025.

Comments: 76 pages, 29 figures, to be submitted

Closeby Habitable Exoplanet Survey (CHES). IV. Synergy between astrometry and direct imaging missions of the Habitable World Observatory for detecting Earth-like planets

Authors: Chunhui Bao, Jianghui Ji, Dongjie Tan, Guo Chen, Xiumin Huang, Su Wang, Yao Dong

Abstract: The detection and characterization of habitable planets around nearby stars persist as one of the foremost objectives in contemporary astrophysics. This work investigates the synergistic integration of astrometric and direct imaging techniques by capitalizing on the complementary capabilities of the Closeby Habitable Exoplanet Survey (CHES) and Habitable Worlds Observatory (HWO). Planetary brightn… ▽ More The detection and characterization of habitable planets around nearby stars persist as one of the foremost objectives in contemporary astrophysics. This work investigates the synergistic integration of astrometric and direct imaging techniques by capitalizing on the complementary capabilities of the Closeby Habitable Exoplanet Survey (CHES) and Habitable Worlds Observatory (HWO). Planetary brightness and position vary over time due to phase effects and orbital architectures, information that can be precisely provided by CHES's astrometric measurements. By combining the precise orbital constraints from CHES with the imaging capabilities of HWO, we evaluate the improvements in detection efficiency, signal-to-noise ratio and overall planet yield. Completeness is quantified as the fraction of injected planets that are successfully detected, while yields are estimated for various scenarios using terrestrial planet occurrence rates derived from the Kepler dataset. Our results indicate that prior astrometric data significantly enhance detection efficiency. Under the adopted detection limit, our analysis indicates that prior CHES observations can increase completeness by approximately 10% and improve detection efficiency by factors ranging from two to thirty. The findings underscore the importance of interdisciplinary approaches in the search for and characterization of habitable worlds. △ Less

Submitted 5 May, 2025; originally announced May 2025.

Comments: 22 pages, 9 figures, accepted for publication in AJ

Asymmetry and Dynamical Constraints in 2-Limbs Retrieval of WASP-39 b Inferring from JWST Data

Authors: Zixin Chen, Jianghui Ji, Guo Chen, Fei Yan, Xianyu Tan

Abstract: Transmission spectroscopy has provided unprecedented insight into the makeup of exoplanet atmospheres. A transmission spectrum contains contributions from a planet's morning and evening limbs, which can differ in temperature, composition and aerosol properties due to atmospheric circulation. While high-resolution ground-based observations have identified limb asymmetry in several ultra-hot/hot exo… ▽ More Transmission spectroscopy has provided unprecedented insight into the makeup of exoplanet atmospheres. A transmission spectrum contains contributions from a planet's morning and evening limbs, which can differ in temperature, composition and aerosol properties due to atmospheric circulation. While high-resolution ground-based observations have identified limb asymmetry in several ultra-hot/hot exoplanets, space-based studies of limb asymmetry are still in their early stages. The prevalence of limb asymmetry across a broad range of exoplanets remains largely unexplored. We conduct a comparative analysis of retrievals on transmission spectra, including traditional 1D approaches and four 2D models that account for limb asymmetry. Two of these 2D models include our newly proposed dynamical constraints derived from shallow-water simulations to provide physically-motivated temperature differences between limbs. Our analysis of WASP-39 b using JWST observations and previous combined datasets (HST, VLT, and Spitzer) strongly favors 2D retrievals over traditional 1D approaches, confirming significant limb asymmetry in this hot Jupiter. Within our 2D framework, unconstrained models recover larger temperature contrasts than dynamically-constrained models, with improved fits to specific spectral features, although Bayesian evidence cannot definitively distinguish between these 2D approaches. Our results support the presence of homogeneous C/O in both the morning and evening atmospheres, but with temperature differences leading to variations in clouds and hazes. Using this treatment, we can study a larger sample of hot Jupiters to gain insights into atmospheric limb asymmetries on these planets. △ Less

Submitted 1 April, 2025; originally announced April 2025.

Comments: 16 pages, 6 figures, accepted for publication in AJ

Closeby Habitable Exoplanet Survey (CHES). III. Retrieval of Planetary Masses in Binaries Using the N-body Model with RV and Astrometry Synergy

Authors: Xiumin Huang, Jianghui Ji, Chunhui Bao, Dongjie Tan, Su Wang, Yao Dong, Guo Chen

Abstract: Given that secular perturbations in a binary system not only excite high orbital eccentricities but also alter the planetary orbital inclination, the classical Keplerian orbital model is no longer applicable for orbital retrieval. The combination of a dynamical model and observational data is essential for characterizing the configuration and planetary mass in close binaries. We calculate the theo… ▽ More Given that secular perturbations in a binary system not only excite high orbital eccentricities but also alter the planetary orbital inclination, the classical Keplerian orbital model is no longer applicable for orbital retrieval. The combination of a dynamical model and observational data is essential for characterizing the configuration and planetary mass in close binaries. We calculate the theoretical radial velocity (RV) signal in the N-body framework and observe a drift in the RV semi-amplitude, which leads to a reduction in the $m$sin$i$ detection threshold by 20 $M_{\oplus}$, with $\sim$ 100% detection probability in the $m_1$sin$i_1$-$a_1$ parameter space. High-precision RV data with an accuracy of 1 m/s can detect such dynamical effects. For four close-in binaries-GJ 86, GJ 3021, HD 196885, and HD 41004, the deviation between the minimum mass derived from the Keplerian and N-body models is found to be $> 0.2 ~ M_{\mathrm{Jup}}$. High-precision astrometric data are also necessary to resolve the 3D orbits and true masses exoplanets. We generate astrometric simulation data with accuracies corresponding to Gaia (57.8 $μ$as) and the Closeby Habitable Exoplanet Survey (CHES) (1 $μ$as), respectively. Joint orbit fitting is performed for RV + Gaia and RV + CHES synergy methods. Compared with the fitting results from the astrometry-only model, the synergy models more effectively constrain the range of orbital inclinations. Using simulation data, we derive precise uncertainties for the true planetary mass, which are critical for determining the evolution of planets around binary stars and multi-planet systems. △ Less

Submitted 21 March, 2025; originally announced March 2025.

Comments: 18 pages, 8 figures, accepted for publication in ApJ

Formation of Ultra-short-period Planet in Hot Jupiter Systems: Application to WASP-47

Authors: Su Wang, Mengrui Pan, Yao Dong, Gang Zhao, Jianghui Ji

Abstract: The WASP-47 system is notable as the first known system hosting both inner and outer low-mass planetary companions around a hot Jupiter, with an ultra-short-period (USP) planet as the innermost planetary companion. The formation of such an unique configuration poses challenges to the lonely hot Jupiter formation model. Hot Jupiters in multiple planetary systems may have a similar formation process… ▽ More The WASP-47 system is notable as the first known system hosting both inner and outer low-mass planetary companions around a hot Jupiter, with an ultra-short-period (USP) planet as the innermost planetary companion. The formation of such an unique configuration poses challenges to the lonely hot Jupiter formation model. Hot Jupiters in multiple planetary systems may have a similar formation process with warm Jupiter systems, which are more commonly found with companions. This implies that the WASP-47 system could bridge our understanding of both hot and warm Jupiter formation. In this work, we propose a possible formation scenario for the WASP-47 system based on its orbital configuration. The mean motion resonance trapping, giant planet perturbations, and tidal effects caused by the central star are key factors in the formation of USP planets in multiple planetary systems with hot Jupiters. Whether a planet can become an USP planet or a short period super-earth (SPSE) planet depends on the competition between eccentricity excitation by nearby giant planet perturbations and the eccentricity damping due to tidal effects. The $Q_p'$ value of the innermost planet is essential for the final planetary configuration. Our results suggest that a $Q_p'$ in the range of [1, 10] is favorable for the formation of the WASP-47 system. Based on the formation scenario, we estimate an occurrence rate of 8.4$\pm$2.4\% for USP planets in systems similar to WASP-47. △ Less

Submitted 2 March, 2025; originally announced March 2025.

Comments: 14 pages, 7 figures, accepted for publication in ApJL

Closeby Habitable Exoplanet Survey (CHES). II. An Observation Strategy for the Target Stars

Authors: Dongjie Tan, Jianghui Ji, Chunhui Bao, Xiumin Huang, Guo Chen, Su Wang, Yao Dong, Haitao Li, Junbo Zhang, Liang Fang, Dong Li, Lei Deng, Jiacheng Liu, Zi Zhu

Abstract: The Closeby Habitable Exoplanet Survey (CHES) constitutes a mission intricately designed to systematically survey approximately 100 solar-type stars located within the immediate proximity of the solar system, specifically within a range of 10 parsecs. The core objective of this mission is the detection and characterization of potentially habitable Earth-like planets or super-Earths within the habi… ▽ More The Closeby Habitable Exoplanet Survey (CHES) constitutes a mission intricately designed to systematically survey approximately 100 solar-type stars located within the immediate proximity of the solar system, specifically within a range of 10 parsecs. The core objective of this mission is the detection and characterization of potentially habitable Earth-like planets or super-Earths within the habitable zone of these stars. The CHES mission obtains high-precision astrometric measurements of planets orbiting the target stars by observing angular distance variations between the target star and reference stars. As a result, we surveyed the relevant parameters of both target and reference stars in detail, conducting a thorough analysis and calculation of the required observation accuracy, the number of observations, and the priority assigned to each target star. Observational emphasis will be concentrated on targets considered of higher priority, ensuring the effectiveness of their observation capabilities. Through this approach, we formulate a five-year observation strategy that will cover all the target stars within a six-month timeframe. The strategy not only fulfills the required observing capability but also exhibit high efficiency simultaneously, providing an executable program for future mission. Over the span of the mission's five-year duration, a cumulative observation time of 29,220 hours will be available. Approximately 86 percent of this, totaling 25,120 hours, is allocated for the observation of target stars. This allocation leaves approximately 4,100 hours for extended scientific observation programs. We have also performed simulated observations based on this strategy and verified its observational capability for exoplanets. △ Less

Submitted 12 August, 2024; originally announced August 2024.

Comments: 20 pages, 12 figures, accepted for publication in AJ

Evidence for Nightside Water Emission Found in Transit of Ultrahot Jupiter WASP-33b

Authors: Yuanheng Yang, Guo Chen, Fei Yan, Xianyu Tan, Jianghui Ji

Abstract: To date, the dayside thermal structure of ultrahot Jupiters (UHJs) is generally considered to be inverted, but their nightside thermal structure has been less explored. Here we explore the impact of nightside thermal emission on high-resolution infrared transmission spectroscopy, which should not be neglected, especially for UHJs. We present a general equation for the high-resolution transmission… ▽ More To date, the dayside thermal structure of ultrahot Jupiters (UHJs) is generally considered to be inverted, but their nightside thermal structure has been less explored. Here we explore the impact of nightside thermal emission on high-resolution infrared transmission spectroscopy, which should not be neglected, especially for UHJs. We present a general equation for the high-resolution transmission spectrum that includes planetary nightside thermal emission. This provides a new way to infer the thermal structure of the planetary nightside with high-resolution transmission spectroscopy. Using the cross-correlation technique, we find evidence for the presence of an H$_2$O emission signature on the UHJ WASP-33b during the transit, indicating an inverted temperature structure on its nightside. Such a result suggests a stronger heat transport through the circulation than currently expected. An alternative explanation is that the rotating visible hemisphere during transit leads to the potential contribution of the limb and dayside atmospheres to the detected emission signature. In the future, the combination of high-resolution full-phase curve spectroscopic observations and general circulation models will hopefully solve this puzzle and provide a complete picture of the three-dimensional nature of the chemistry, circulation, and thermal structure of UHJs. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Comments: 17 pages, 9 figures, 1 table. Accepted for publication in ApJL

Closeby Habitable Exoplanet Survey (CHES). I. Astrometric Noise and Planetary Detection Efficiency due to Stellar Spots and Faculae

Authors: Chunhui Bao, Jianghui Ji, Dongjie Tan, Guo Chen, Xiumin Huang, Su Wang, Yao Dong

Abstract: The Closeby Habitable Exoplanet Survey (CHES) is dedicated to the astrometric exploration for habitable-zone Earth-like planets orbiting solar-type stars in close proximity, achieving unprecedented micro-arcsecond precision. Given the elevated precision, thorough consideration of photocenter jitters induced by stellar activity becomes imperative. This study endeavors to model the stellar activity… ▽ More The Closeby Habitable Exoplanet Survey (CHES) is dedicated to the astrometric exploration for habitable-zone Earth-like planets orbiting solar-type stars in close proximity, achieving unprecedented micro-arcsecond precision. Given the elevated precision, thorough consideration of photocenter jitters induced by stellar activity becomes imperative. This study endeavors to model the stellar activity of solar-type stars, compute astrometric noise, and delineate the detection limits of habitable planets within the astrometric domain. Simulations were conducted for identified primary targets of CHES, involving the generation of simulated observed data for astrometry and photometry, accounting for the impact of stellar activity. Estimation of activity levels in our samples was achieved through chromospheric activity indices, revealing that over 90% of stars exhibited photocenter jitters below 1 $μ\mathrm{as}$. Notably, certain proximate stars, such as $α$ Cen A and B, displayed more discernible noise arising from stellar activity. Subsequent tests were performed to evaluate detection performance, unveiling that stellar activity tends to have a less pronounced impact on planetary detectability for the majority of stars. Approximately 95% of targets demonstrated a detection efficiency exceeding 80%. However, for several cold stars, e.g., HD 32450 and HD 21531, with the habitable zones close to the stars, a reduction in detection efficiency was observed. These findings offer invaluable insights into the intricate interplay between stellar activity and astrometric precision, significantly advancing our understanding in the search for habitable planets. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 18 pages, 10 figures, accepted for publication in AJ

Journal ref: The Astronomical Journal, Volume 167, Number 6, article id.286 (2024)

Forming Giant Planets Around Late-M Dwarfs: Pebble Accretion and Planet-Planet Collision

Authors: Mengrui Pan, Beibei Liu, Anders Johansen, Masahiro Ogihara, Su Wang, Jianghui Ji, Sharon X. Wang, Fabo Feng, Ignasi Riba

Abstract: We propose a pebble-driven core accretion scenario to explain the formation of giant planets around the late-M dwarfs of $M_{\star}{=}0.1{-}0.2 \ M_{\odot}$. In order to explore the optimal disk conditions for giant planet, we perform N-body simulations to investigate the growth and dynamical evolution of both single and multiple protoplanets in the disks with both inner viscously heated and outer… ▽ More We propose a pebble-driven core accretion scenario to explain the formation of giant planets around the late-M dwarfs of $M_{\star}{=}0.1{-}0.2 \ M_{\odot}$. In order to explore the optimal disk conditions for giant planet, we perform N-body simulations to investigate the growth and dynamical evolution of both single and multiple protoplanets in the disks with both inner viscously heated and outer stellar irradiated regions. The initial masses of the protoplanets are either assumed to be equal to $0.01 \ M_{\oplus}$ or calculated based on the formula derived from streaming instability simulations. Our findings indicate that massive planets are more likely to form in disks with longer lifetimes, higher solid masses, moderate to high levels of disk turbulence, and larger initial masses of protoplanets. In the single protoplanet growth cases, the highest planet core mass that can be reached is generally lower than the threshold necessary to trigger rapid gas accretion, which impedes the formation of giant planets. Nonetheless, in multi-protoplanet cases, the cores can exceed the pebble isolation mass barrier aided by frequent planet-planet collisions. This consequently speeds up their gas accretion and promotes giant planet formation, making the optimal parameter space to grow giant planets substantially wider. Taken together, our results suggest that even around very low-mass stellar hosts, the giant planets with orbital periods of ${\lesssim}100$ days are still likely to form when lunar-mass protoplanets first emerge from planetesimal accretion and then grow rapidly by a combination of pebble accretion and planet-planet collisions in disks with a high supply of pebble reservoir ${>}50 \ M_{\oplus}$ and turbulent level of $α_{\rm t} {\sim} 10^{-3}{-}10^{-2}$. △ Less

Submitted 16 November, 2023; originally announced November 2023.

Comments: 22 pages, 13 figures, 3 tables, Accepted for publication in A&A

PyMsOfa: A Python Package for the Standards of Fundamental Astronomy (SOFA) Service

Authors: Jianghui Ji, Dongjie Tan, Chunhui Bao, Xiumin Huang, Shoucun Hu, Yao Dong, Su Wang

Abstract: The Standards of Fundamental Astronomy (SOFA) is a service provided by the International Astronomical Union (IAU) that offers algorithms and software for astronomical calculations, which was released in two versions by FORTRAN 77 and ANSI C, respectively. In this work, we implement the python package PyMsOfa for SOFA service by three ways: (1) a python wrapper package based on a foreign function l… ▽ More The Standards of Fundamental Astronomy (SOFA) is a service provided by the International Astronomical Union (IAU) that offers algorithms and software for astronomical calculations, which was released in two versions by FORTRAN 77 and ANSI C, respectively. In this work, we implement the python package PyMsOfa for SOFA service by three ways: (1) a python wrapper package based on a foreign function library for Python (ctypes), (2) a python wrapper package with the foreign function interface for Python calling C code (cffi), and (3) a python package directly written in pure python codes from SOFA subroutines. The package PyMsOfa has fully implemented 247 functions of the original SOFA routines. In addition, PyMsOfa is also extensively examined, which is exactly consistent with those test examples given by the original SOFA. This python package can be suitable to not only the astrometric detection of habitable planets of the Closeby Habitable Exoplanet Survey (CHES) mission (Ji et al. 2022), but also for the frontiers themes of black holes and dark matter related to astrometric calculations and other fields. The source codes are available via https://github.com/CHES2023/PyMsOfa. △ Less

Submitted 17 October, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: 7 pages, 5 figures, accepted to Research in Astronomy and Astrophysics

Peculiar orbital characteristics of Earth quasi-satellite 469219 Kamo`oalewa: implications for the Yarkovsky detection and orbital uncertainty propagation

Authors: Shoucun Hu, Bin Li, Haoxuan Jiang, Gang Bao, Jianghui Ji

Abstract: 469219 Kamo`oalewa is selected as one of the primary targets of Tianwen-2 mission, which is currently believed to be the most stable quasi-satellite of Earth. Here we derive a weak detection of the Yarkovsky effect for Kamo`oalewa, giving $A_2 = -1.075\pm0.447\times 10^{-13} \rm{au/d}^2$, with the available ground-based optical observations from Minor Planet Center and a relatively conservative we… ▽ More 469219 Kamo`oalewa is selected as one of the primary targets of Tianwen-2 mission, which is currently believed to be the most stable quasi-satellite of Earth. Here we derive a weak detection of the Yarkovsky effect for Kamo`oalewa, giving $A_2 = -1.075\pm0.447\times 10^{-13} \rm{au/d}^2$, with the available ground-based optical observations from Minor Planet Center and a relatively conservative weighting scheme. Due to the quasi-satellite resonance with Earth, we show that the detection of Yarkovsky effect by orbital fitting with astrometric observations becomes difficult as its orbital drift shows a slow oscillatory growth resulting from the Yarkovsky effect. In addition, we extensively explore the characteristics of orbital uncertainty propagation and find that the positional uncertainty mainly arises from the geocentric radial direction in 2010-2020, and then concentrates in the heliocentric transverse direction in 2020-2030. Furthermore, the heliocentric transverse uncertainty is clearly monthly dependent, which can arrive at a minimum around January and a maximum around July as the orbit moves towards the leading and trailing edges, respectively, in 2025-2027. Finally, we investigate a long-term uncertainty propagation in the quasi-satellite regime, implying that the quasi-satellite resonance with Earth may play a crucial role in constraining the increase of uncertainty over time. Such interesting feature further implies that the orbital precision of Kamo`oalewa is relatively stable at its quasi-satellite phase, which may also be true for other quasi-satellites of Earth. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 15 pages, 11 figures, accepted for publication in AJ

Evolution of the Planetary Obliquity: The Eccentric Kozai-Lidov Mechanism Coupled with Tide

Authors: Xiumin Huang, Jianghui Ji, Shangfei Liu, Ruobing Dong, Su Wang

Abstract: The planetary obliquity plays a significant role in determining physical properties of planetary surfaces and climate. As direct detection is constrained due to the present observation accuracy, kinetic theories are helpful to predict the evolution of the planetary obliquity. Here the coupling effect between the eccentric Kozai-Lidov (EKL) effect and the equilibrium tide is extensively investigate… ▽ More The planetary obliquity plays a significant role in determining physical properties of planetary surfaces and climate. As direct detection is constrained due to the present observation accuracy, kinetic theories are helpful to predict the evolution of the planetary obliquity. Here the coupling effect between the eccentric Kozai-Lidov (EKL) effect and the equilibrium tide is extensively investigated, the planetary obliquity performs to follow two kinds of secular evolution paths, based on the conservation of total angular momentum. The equilibrium timescale of the planetary obliquity $t_{\mathrm{eq}}$ varies along with $r_{t}$, which is defined as the initial timescale ratio of the tidal dissipation and secular perturbation. We numerically derive the linear relationship between $t_{\mathrm{eq}}$ and $r_{t}$ with the maximum likelihood method. The spin-axis orientation of S-type terrestrials orbiting M-dwarfs reverses over $90^\circ$ when $r_{t} > 100$, then enter the quasi-equilibrium state between $40^\circ$ and $60^\circ$, while the maximum obliquity can reach $130^\circ$ when $r_{t} > 10^4 $. Numerical simulations show that the maximum obliquity increases with the semi-major axis ratio $a_1$/$a_2$, but is not so sensitive to the eccentricity $e_2$. The likelihood of obliquity flip for S-type terrestrials in general systems with $a_2 < 45$ AU is closely related to $m_1$. The observed potential oblique S-type planets HD 42936 b, GJ 86 Ab and $τ$ Boot Ab are explored to have a great possibility to be head-down over the secular evolution of spin. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 18 pages, 12 figures, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 956, Number 1, article id.45 (2023)

Mid-IR Observations of IRAS, AKARI, WISE/NEOWISE and Subaru for Large Icy Asteroid (704) Interamnia: a New Perspective of Regolith Properties and Water Ice Fraction

Authors: Haoxuan Jiang, Jianghui Ji, Liangliang Yu, Bin Yang, Shoucun Hu, Yuhui Zhao

Abstract: (704) Interamnia is one of the largest asteroids that locates in the outer main-belt region, which may contain a large amount of water ice underneath its surface. We observe this asteroid using 8.2 m Subaru telescope at mid-infrared wavebands, and utilize thermophysical model for realistic surface layers (RSTPM) to analyze mid-infrared data from Subaru along with those of IRAS, AKARI and WISE/NEOW… ▽ More (704) Interamnia is one of the largest asteroids that locates in the outer main-belt region, which may contain a large amount of water ice underneath its surface. We observe this asteroid using 8.2 m Subaru telescope at mid-infrared wavebands, and utilize thermophysical model for realistic surface layers (RSTPM) to analyze mid-infrared data from Subaru along with those of IRAS, AKARI and WISE/NEOWISE. We optimize the method to convert the WISE magnitude to thermal infrared flux with temperature dependent color corrections, which can provide significant references for main-belt asteroids at a large heliocentric distance with low surface temperature. We derive best-fitting thermal parameters of Interamnia - a mean regolith grain size of $190_{-180}^{+460}~\rm μm$, with a roughness of $0.30_{-0.17}^{+0.35}$ and RMS slope of $27_{-9}^{+13}$ degrees, thereby producing thermal inertia ranging from 9 to $92~\rm Jm^{-2}s^{-1/2}K^{-1}$ due to seasonal temperature variation. The geometric albedo and effective diameter are evaluated to be $0.0472_{-0.0031}^{+0.0033}$ and $339_{-11}^{+12}~\rm km$, respectively, being indicative of a bulk density of $1.86\pm0.63~\rm g/cm^3$. The low thermal inertia is consistent with typical B/C-type asteroids with $D\geq100$ km. The tiny regolith grain size suggests the presence of a fine regolith on the surface of Interamnia. Moreover, the seasonal and diurnal temperature distribution indicates that thermal features between southern and northern hemisphere appear to be very different. Finally, we present an estimation of volume fraction of water ice of $9\%\sim66\%$ from the published grain density and porosity of carbonaceous chondrites. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 17 pages, 10 figures, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 944, Number 2, article id.202 (2023)

The X-shaped radio galaxy J0725+5835 is associated with an AGN pair

Authors: Xiaolong Yang, Jialu Ji, Ravi Joshi, Jun Yang, Tao An, Ran Wang, Luis C. Ho, David H. Roberts, Lakshmi Saripalli

Abstract: X-shaped radio galaxies (XRGs) are those that exhibit two pairs of unaligned radio lobes (main radio lobes and "wings"), one of the promising models for the peculiar morphology is jet re-orientation. To clarify it, we conducted the European VLBI Network (EVN) 5 GHz observation of an XRG J0725+5835, which resembles the archetypal binary AGNs 0402+379 in radio morphology but it is larger in angular… ▽ More X-shaped radio galaxies (XRGs) are those that exhibit two pairs of unaligned radio lobes (main radio lobes and "wings"), one of the promising models for the peculiar morphology is jet re-orientation. To clarify it, we conducted the European VLBI Network (EVN) 5 GHz observation of an XRG J0725+5835, which resembles the archetypal binary AGNs 0402+379 in radio morphology but it is larger in angular size. In our observation, two milliarcsec (mas) scale radio components with non-thermal radio emission are detected, each of them coincides with an optical counterpart with similar photometric redshift and (optical and infrared) magnitude, corresponding to dual active nuclei. Furthermore, with the improved VLA images, we find a bridge between the two radio cores and a jet bending in the region surrounding the companion galaxy, which further supports the interplay between the main and companion galaxies. In addition, we also report the discovery of an arcsec-scale jet in the companion. Given the projected separation of $\sim100$ kpc between the main and companion galaxies, XRG J0725+5835 is likely associated with a dual jetted-AGN system. In both EVN and VLA observations, we find signatures that the jet is changing its direction, which is likely responsible for the X-shaped morphology. On the origin of jet re-orientation, several scenarios are discussed. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: 12 pages, 8 figures, 2 table, published in ApJ

Report number: 2022ApJ...933...98Y

LSST Survey Strategies and Brown Dwarf Parallaxes

Authors: John E. Gizis, Peter Yoachim, R. Lynne Jones, Dylan Hilligoss, Jinbiao Ji

Abstract: The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has the potential to measure parallaxes for thousands of nearby ultracool dwarfs, enabling improved measurements of the brown dwarf luminosity function. We develop a simple model to estimate the number of L dwarfs and T dwarfs with parallaxes with signal-to-noise ratio greater than ten in the baseline LSST survey. High quality… ▽ More The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) has the potential to measure parallaxes for thousands of nearby ultracool dwarfs, enabling improved measurements of the brown dwarf luminosity function. We develop a simple model to estimate the number of L dwarfs and T dwarfs with parallaxes with signal-to-noise ratio greater than ten in the baseline LSST survey. High quality astrometry imposes scheduling constraints. We assess different possible observing strategies using quantitative metrics and make recommendations as part of the LSST community input process. We find that the new substellar parallax sample will represent a nearly order-of-magnitude increase on existing samples, with ~50-100 objects per spectral type bin for late-L to mid-T dwarfs. The sample size is robust (+/- 5%) against most survey strategy changes under consideration, although we do identify areas of tension with other uses of twilight time that could have larger impact. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: Accepted to the Astrophysical Journal Supplement, for special issue on Rubin LSST Science Optimization. https://iopscience.iop.org/journal/0067-0049/page/rubin_cadence

Extremely Inclined Orbit of S-type Planet $γ$ Cep Ab Induced by Eccentric Kozai-Lidov Mechanism

Authors: Xiumin Huang, Jianghui Ji

Abstract: $γ$ Cep Ab is a typical S-type planet, which occupies a nearly perpendicular planetary orbit relative to the binary. Here we use the Markov Chain Monte Carlo (MCMC) sampler to conduct full N-body fitting and derive self-consistent orbital solutions for this hierarchical system. Then we employ the Eccentric Kozai-Lidov (EKL) mechanism to explain the extremely inclined orbit of S-type planet $γ… ▽ More $γ$ Cep Ab is a typical S-type planet, which occupies a nearly perpendicular planetary orbit relative to the binary. Here we use the Markov Chain Monte Carlo (MCMC) sampler to conduct full N-body fitting and derive self-consistent orbital solutions for this hierarchical system. Then we employ the Eccentric Kozai-Lidov (EKL) mechanism to explain the extremely inclined orbit of S-type planet $γ$ Cep Ab. The EKL mechanism plays an essential role in exploring significant oscillations of the mutual inclination $i_{\mathrm{mut}}$ between the planet and the secondary star. We perform qualitative analysis and extensive numerical integrations to investigate the flip conditions and timescales of $γ$ Cep Ab's orbit. When the planetary mass is 15 $M_{\mathrm{Jup}}$, the planet can reach $i_{\mathrm{mut}} \sim$ 113$^{\circ}$ with the critical initial conditions of $i_{\mathrm{mut}} < 60^{\circ}$ and $e_1<0.7$. The timescale for the first orbital flip decreases with the increase of the perturbation Hamiltonian. Flipping orbits of $γ$ Cep Ab are confirmed to have a large possibility to retain stable based on surfaces of section and the secular stability criterion. Furthermore, we extend the application of EKL to general S-type planetary systems with $a_1/a_2\leq0.1$, where the most intense excitation of $i_{\mathrm{mut}}$ occurs when $a_1/a_2=0.1$ and $e_2 \sim 0.8$, and the variation of planetary mass mainly affect the flip possibility where $e_1\leq 0.3$. △ Less

Submitted 3 September, 2022; originally announced September 2022.

Comments: 21 pages, 14 figures, accepted for publication in AJ

Journal ref: The Astronomical Journal, Volume 164, Number 5, article id.177 (2022)

CHES: a space-borne astrometric mission for the detection of habitable planets of the nearby solar-type stars

Authors: Jianghui Ji, Haitao Li, Junbo Zhang, Liang Fang, Dong Li, Su Wang, Yang Cao, Lei Deng, Baoquan Li, Hao Xian, Xiaodong Gao, Ang Zhang, Fei Li, Jiacheng Liu, Zhaoxiang Qi, Sheng Jin, Yaning Liu, Guo Chen, Mingtao Li, Yao Dong, Zi Zhu, CHES consortium

Abstract: The Closeby Habitable Exoplanet Survey (CHES) mission is proposed to discover habitable-zone Earth-like planets of the nearby solar-type stars ($\sim 10~\mathrm{pc}$ away from our solar system) via micro-arcsecond relative astrometry. The major scientific objectives of CHES are: to search for Earth Twins or terrestrial planets in habitable zones orbiting 100 FGK nearby stars; further to conduct a… ▽ More The Closeby Habitable Exoplanet Survey (CHES) mission is proposed to discover habitable-zone Earth-like planets of the nearby solar-type stars ($\sim 10~\mathrm{pc}$ away from our solar system) via micro-arcsecond relative astrometry. The major scientific objectives of CHES are: to search for Earth Twins or terrestrial planets in habitable zones orbiting 100 FGK nearby stars; further to conduct a comprehensive survey and extensively characterize the nearby planetary systems. The primary payload is a high-quality, low-distortion, high-stability telescope. The optical subsystem is a coaxial three-mirror anastigmat (TMA) with a $1.2 \mathrm{~m}$-aperture, $0.44^{\circ} \times 0.44^{\circ}$ field of view and $500 \mathrm{~nm}-900 \mathrm{~nm}$ working waveband. The camera focal plane is composed of 81 MOSAIC scientific CMOS detectors each with $4 \mathrm{~K} \times 4 \mathrm{~K}$ pixels. The heterodyne laser interferometric calibration technology is employed to ensure micro-arcsecond level (1 $μ$as) relative astrometry precision to meet the requirements for detection of Earth-like planets. CHES satellite operates at the Sun-Earth L2 point and observes the entire target stars for 5 years. CHES will offer the first direct measurements of true masses and inclinations of Earth Twins and super-Earths orbiting our neighbor stars based on micro-arcsecond astrometry from space. This will definitely enhance our understanding of the formation of diverse nearby planetary systems and the emergence of other worlds for solar-type stars, and finally to reflect the evolution of our own solar system. △ Less

Submitted 9 June, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: 39 pages, 37 figures, Invited Review, accepted to Research in Astronomy and Astrophysics

Journal ref: Res. Astron. Astrophys., Volume 22, Number 7, article id. 072003 (2022)

The Terrestrial Planet Formation around M Dwarfs: In-situ, Inward Migration or Reversed Migration

Authors: Mengrui Pan, Su Wang, Jianghui Ji

Abstract: Terrestrial planets are commonly observed to orbit M dwarfs with close-in trajectories. In this work, we extensively perform N-body simulations of planetesimal accretion with three models of in-situ, inward migration and reversed migration to explore terrestrial formation in tightly compact systems of M dwarfs. In the simulations, the solid disks are assumed to be 0.01\% of the masses of host star… ▽ More Terrestrial planets are commonly observed to orbit M dwarfs with close-in trajectories. In this work, we extensively perform N-body simulations of planetesimal accretion with three models of in-situ, inward migration and reversed migration to explore terrestrial formation in tightly compact systems of M dwarfs. In the simulations, the solid disks are assumed to be 0.01\% of the masses of host stars and spread from 0.01 to 0.5 AU with the surface density profile scaling with $r^{-k}$ according to the observations. Our results show that in-situ scenario may produce $7.77^{+3.23}_{-3.77}$ terrestrial planets with an average mass of $1.23^{+4.01}_{-0.93} \ M_{\oplus}$ around M dwarfs. The number of planets tends to increase as the disk slope is steeper or with a larger stellar mass. Moreover, we show that $2.55^{+1.45}_{-1.55}$ planets with mass of $3.76^{+8.77}_{-3.46} \ M_{\oplus}$ are formed in the systems via inward migration, while $2.85^{+1.15}_{-0.85}$ planets with $3.01^{+13.77}_{-2.71} \ M_{\oplus}$ are yielded under reversed migration. Migration scenarios can also deliver plentiful water from the exterior of ice line to the interior due to more efficient accretion. The simulation outcomes of reversed migration model produce the best matching with observations, being suggestive of a likely mechanism for planetary formation around M dwarfs. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 510, 4134-4145 (2022)

Nii: a Bayesian orbit retrieval code applied to differential astrometry

Authors: Sheng Jin, Xiaojian Ding, Su Wang, Yao Dong, Jianghui Ji

Abstract: Here we present an open source Python-based Bayesian orbit retrieval code (Nii) that implements an automatic parallel tempering Markov chain Monte Carlo (APT-MCMC) strategy. Nii provides a module to simulate the observations of a space-based astrometry mission in the search for exoplanets, a signal extraction process for differential astrometric measurements using multiple reference stars, and an… ▽ More Here we present an open source Python-based Bayesian orbit retrieval code (Nii) that implements an automatic parallel tempering Markov chain Monte Carlo (APT-MCMC) strategy. Nii provides a module to simulate the observations of a space-based astrometry mission in the search for exoplanets, a signal extraction process for differential astrometric measurements using multiple reference stars, and an orbital parameter retrieval framework using APT-MCMC. We further verify the orbit retrieval ability of the code through two examples corresponding to a single-planet system and a dual-planet system. In both cases, efficient convergence on the posterior probability distribution can be achieved. Although this code specifically focuses on the orbital parameter retrieval problem of differential astrometry, Nii can also be widely used in other Bayesian analysis applications. △ Less

Submitted 12 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

FOSSIL: I. The Spin Rate Limit of Jupiter Trojans

Authors: Chan-Kao Chang, Ying-Tung Chen, Wesley C. Fraser, Fumi Yoshida, Matthew J. Lehner, Shiang-Yu Wang, JJ Kavelaars, Rosemary E. Pike, Mike Alexandersen, Takashi Ito, Young-Jun Choi, A. Paula Granados Contreras, Youngmin JeongAhn, Jianghui Ji, Myung-Jin Kim, Samantha M. Lawler, Jian Li, Zhong-Yi Lin, Patryk Sofia Lykawka, Hong-Kyu Moon, Surhud More, Marco Munoz-Gutierrez, Keiji Ohtsuki, Tsuyoshi Terai, Seitaro Urakawa , et al. (3 additional authors not shown)

Abstract: Rotation periods of 53 small (diameters $2 < D < 40$ km) Jupiter Trojans (JTs) were derived using the high-cadence light curves obtained by the FOSSIL phase I survey, a Subaru/Hyper Suprime-Cam intensive program. These are the first reported periods measured for JTs with $D < 10$ km. We found a lower limit of the rotation period near 4 hr, instead of the previously published result of 5 hr (Ryan e… ▽ More Rotation periods of 53 small (diameters $2 < D < 40$ km) Jupiter Trojans (JTs) were derived using the high-cadence light curves obtained by the FOSSIL phase I survey, a Subaru/Hyper Suprime-Cam intensive program. These are the first reported periods measured for JTs with $D < 10$ km. We found a lower limit of the rotation period near 4 hr, instead of the previously published result of 5 hr (Ryan et al. 2017; Szabo et al. 2017, 2020) found for larger JTs. Assuming a rubble-pile structure for JTs, a bulk density of 0.9 gcm$^{-3}$ is required to withstand this spin rate limit, consistent with the value $0.8-1.0$ gcm$^{-3}$ (Marchis et al. 2006; Mueller et al. 2010; Buie et al. 2015; Berthier et al. 2020) derived from the binary JT system, (617) Patroclus-Menoetius system. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: Accepted by PSJ on July 9th, 2021. 15 pages, 7 figures, and 3 tables

Thermophysical Modeling of 20 Themis Family Asteroids with WISE/NEOWISE Observations

Authors: Haoxuan Jiang, Jianghui Ji

Abstract: Themis family is one of the largest and oldest asteroid populations in the main-belt. Water-ice may widely exist on the parent body (24) Themis. In this work, we employ the Advanced Thermophysical Model as well as mid-infrared measurements from NASA's Wide-Field Infrared Survey Explorer to explore thermal parameters of 20 Themis family members. Here we show that the average thermal inertia and geo… ▽ More Themis family is one of the largest and oldest asteroid populations in the main-belt. Water-ice may widely exist on the parent body (24) Themis. In this work, we employ the Advanced Thermophysical Model as well as mid-infrared measurements from NASA's Wide-Field Infrared Survey Explorer to explore thermal parameters of 20 Themis family members. Here we show that the average thermal inertia and geometric albedo are ~$39.5\pm26.0 ~\rm J m^{-2} s^{-1/2} K^{-1}$ and $0.067\pm0.018$, respectively. The family members have a relatively moderate roughness fraction on their surfaces. We find that the relatively low albedos of Themis members are consistent with the typical values of B-type and C-type asteroids. As aforementioned, Themis family bears a very low thermal inertia, which indicates a fine and mature regolith on their surfaces. The resemblance of thermal inertia and geometric albedo of Themis members may reveal their close connection in origin and evolution. In addition, we present the compared results of thermal parameters for several prominent families. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 22 pages, 25 figures, accepted for publication in AJ

Journal ref: The Astronomical Journal, Volume 162, Number 2, article id.40 (2021)

A Gaussian Process Regression Reveals No Evidence for Planets Orbiting Kapteyn's Star

Authors: Anna Bortle, Hallie Fausey, Jinbiao Ji, Sarah Dodson-Robinson, Victor Ramirez Delgado, John Gizis

Abstract: Radial-velocity (RV) planet searches are often polluted by signals caused by gas motion at the star's surface. Stellar activity can mimic or mask changes in the RVs caused by orbiting planets, resulting in false positives or missed detections. Here we use Gaussian process (GP) regression to disentangle the contradictory reports of planets vs. rotation artifacts in Kapteyn's star (Anglada-Escude et… ▽ More Radial-velocity (RV) planet searches are often polluted by signals caused by gas motion at the star's surface. Stellar activity can mimic or mask changes in the RVs caused by orbiting planets, resulting in false positives or missed detections. Here we use Gaussian process (GP) regression to disentangle the contradictory reports of planets vs. rotation artifacts in Kapteyn's star (Anglada-Escude et al. 2014, Robertson et al. 2015, Anglada-Escude et al. 2016). To model rotation, we use joint quasi-periodic kernels for the RV and H-alpha signals, requiring that their periods and correlation timescales be the same. We find that the rotation period of Kapteyn's star is 125 days, while the characteristic active-region lifetime is 694 days. Adding a planet to the RV model produces a best-fit orbital period of 100~years, or 10 times the observing time baseline, indicating that the observed RVs are best explained by star rotation only. We also find no significant periodic signals in residual RV data sets constructed by subtracting off realizations of the best-fit rotation model and conclude that both previously reported "planets" are artifacts of the star's rotation and activity. Our results highlight the pitfalls of using sinusoids to model quasi-periodic rotation signals. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 15 pages, including 10 figures. Accepted for publication in Astronomical Journal

Non-Sinusoidal Transit Timing Variations for the Exoplanet HAT-P-12b

Authors: Devesh P. Sariya, Ing-Guey Jiang, Li-Hsin Su, Li-Chin Yeh, Tze-En Chang, V. V. Moskvin, A. A. Shlyapnikov, V. Ignatov, David Mkrtichian, Evgeny Griv, Vineet Kumar Mannaday, Parijat Thakur, D. K. Sahu, Swadesh Chand, D. Bisht, Zhao Sun, Jianghui Ji

Abstract: Considering the importance of investigating the transit timing variations (TTVs) of transiting exoplanets, we present a follow-up study of HAT-P-12b. We include six new light curves observed between 2011 and 2015 from three different observatories, in association with 25 light curves taken from the published literature. The sample of the data used, thus covers a time span of about 10.2 years with… ▽ More Considering the importance of investigating the transit timing variations (TTVs) of transiting exoplanets, we present a follow-up study of HAT-P-12b. We include six new light curves observed between 2011 and 2015 from three different observatories, in association with 25 light curves taken from the published literature. The sample of the data used, thus covers a time span of about 10.2 years with a large coverage of epochs (1160) for the transiting events of the exoplanet HAT-P-12b. The light curves are used to determine the orbital parameters and conduct an investigation of possible transit timing variations. The new linear ephemeris shows a large value of reduced chi-square = 7.93, and the sinusoidal fitting using the prominent frequency coming from a periodogram shows a reduced chi-square around 4. Based on these values and the corresponding O-C diagrams, we suspect the presence of a possible non-sinusoidal TTV in this planetary system. Finally, we find that a scenario with an additional non-transiting exoplanet could explain this TTV with an even smaller reduced chi-square value of around 2. △ Less

Submitted 16 December, 2020; originally announced December 2020.

Comments: 19 pages, 7 tables, 6 figures, accepted by RAA

Departure from the Exact Location of Mean Motion Resonances Induced by the Gas Disk in the Systems Observed by Kepler

Authors: Su Wang, D. N. C. Lin, Xiaochen Zheng, Jianghui Ji

Abstract: The statistical results of transiting planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. A large number of planet pairs are found near the exact location of mean motion resonances (MMRs). In this work, we find out that the depletion and structures of gas disk play crucial roles in driving planet pairs out of exact location of MMRs. Under such sce… ▽ More The statistical results of transiting planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. A large number of planet pairs are found near the exact location of mean motion resonances (MMRs). In this work, we find out that the depletion and structures of gas disk play crucial roles in driving planet pairs out of exact location of MMRs. Under such scenario, planet pairs are trapped into exact MMRs during orbital migration firstly and keep migrating in a same pace. The eccentricities can be excited. Due to the existence of gas disk, eccentricities can be damped leading to the change of orbital period. It will make planet pairs depart from the exact location of MMRs. With depletion timescales larger than 1 Myr, near MMRs configurations are formed easily. Planet pairs have higher possibilities to escape from MMRs with higher disk aspect ratio. Additionally, with weaker corotation torque, planet pairs can depart farther from exact location of MMRs. The final location of the innermost planets in systems are directly related to the transition radius from optically thick region to inner optically thin disk. While the transition radius is smaller than 0.2 AU at the late stage of star evolution process, the innermost planets can reach around 10 days. Our formation scenario is a possible mechanism to explain the formation of near MMRs configuration with the innermost planet farther than 0.1 AU. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 18 pages, 7 figures, accepted for publication in AJ

Critical spin periods of sub-km-sized cohesive rubble-pile asteroids: dependencies on material parameters

Authors: Shoucun Hu, Derek C. Richardson, Yun Zhang, Jianghui Ji

Abstract: In this work, we employ a soft-sphere discrete element method with a cohesion implementation to model the dynamical process of sub-km-sized cohesive rubble piles under continuous spinup. The dependencies of critical spin periods $T_c$ on several material parameters for oblate rubble piles with different bulk diameters $D$ are explored. Our numerical simulations show that both the increase of inter… ▽ More In this work, we employ a soft-sphere discrete element method with a cohesion implementation to model the dynamical process of sub-km-sized cohesive rubble piles under continuous spinup. The dependencies of critical spin periods $T_c$ on several material parameters for oblate rubble piles with different bulk diameters $D$ are explored. Our numerical simulations show that both the increase of interparticle cohesion and particle shape parameter in our model can strengthen the bodies, especially for the smaller ones. In addition, we find there exists some critical diameter $D_{cri,ρ}$ at which the variation trend of $T_c$ with the bulk density $ρ$ reverses. Though a greater static friction coefficient $μ_S$ can strengthen the body, this effect attains a minimum at a critical diameter $D_{cri,φ}$ close to $D_{cri,ρ}$. The continuum theory (analytical method) is used for comparison and two equivalent critical diameters are obtained. The numerical results were fitted with the analytical method and the ratio of the interparticle cohesion $c$ to the bulk cohesion $C$ is estimated to be roughly 88.3. We find this ratio keeps constant for different $c$ and $ρ$, while it strongly depends on the friction angle $φ$. Also, our numerical results further show that the dependency of $T_c$ on $φ$ is opposite from that predicted by the continuum theory when $D$ < $D_{cri,φ}$. Finally, we find that the two critical diameters happen to be close to the diameter when the mean normal stress of the body equals zero, which is the separation between the compressive regime and the tensile regime. △ Less

Submitted 26 September, 2020; originally announced September 2020.

A Tale of Planet Formation: From Dust to Planets

Authors: Beibei Liu, Jianghui Ji

Abstract: The characterization of exoplanets and their birth protoplanetary disks has enormously advanced in the last decade. Benefitting from that, our global understanding of the planet formation processes has been substantially improved. In this review, we first summarize the cutting-edge states of the exoplanet and disk observations. We further present a comprehensive panoptic view of modern core accret… ▽ More The characterization of exoplanets and their birth protoplanetary disks has enormously advanced in the last decade. Benefitting from that, our global understanding of the planet formation processes has been substantially improved. In this review, we first summarize the cutting-edge states of the exoplanet and disk observations. We further present a comprehensive panoptic view of modern core accretion planet formation scenarios, including dust growth and radial drift, planetesimal formation by the streaming instability, core growth by planetesimal accretion and pebble accretion. We discuss the key concepts and physical processes in each growth stage and elaborate on the connections between theoretical studies and observational revelations. Finally, we point out the critical questions and future directions of planet formation studies. △ Less

Submitted 1 November, 2020; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 48 pages, 8 figures, invited review for Research in Astronomy and Astrophysics

Near Mean Motion Resonance of Terrestrial Planet Pair Induced by Giant Planet: Application to Kepler-68 System

Authors: Mengrui Pan, Su Wang, Jianghui Ji

Abstract: In this work, we investigate configuration formation of two inner terrestrial planets near mean motion resonance (MMRs) induced by the perturbation of a distant gas-giant for the Kepler-68 system, by conducting thousands of numerical simulations. The results show that the formation of terrestrial planets is relevant to the speed of Type I migration, the mass of planets, and the existence of giant… ▽ More In this work, we investigate configuration formation of two inner terrestrial planets near mean motion resonance (MMRs) induced by the perturbation of a distant gas-giant for the Kepler-68 system, by conducting thousands of numerical simulations. The results show that the formation of terrestrial planets is relevant to the speed of Type I migration, the mass of planets, and the existence of giant planet. The mass and eccentricity of the giant planet may play a crucial role in shaping the final configuration of the system. The inner planet pair can be trapped in 5:3 or 7:4 MMRs if the giant planet revolves the central star with an eccentric orbit, which is similar to the observed configuration of Kepler-68. Moreover, we find that the eccentricity of the middle planet can be excited to roughly 0.2 if the giant planet is more massive than 5 $M_J$,otherwise the terrestrial planets are inclined to remain near-circular orbits. Our study may provide a likely formation scenario for the planetary systems that harbor several terrestrial planets near MMRs inside and one gas-giant exterior to them. △ Less

Submitted 25 June, 2020; originally announced June 2020.

Comments: 13 pages, 8 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 496, 4688-4699 (2020)

Determination of size, albedo and thermal inertia of 10 Vesta family asteroids with WISE/NEOWISE observations

Authors: Haoxuan Jiang, Jianghui Ji, Liangliang Yu

Abstract: In this work, we investigate the size, thermal inertia, surface roughness and geometric albedo of 10 Vesta family asteroids by using the Advanced Thermophysical Model (ATPM), based on the thermal infrared data acquired by mainly NASA's Wide-field Infrared Survey Explorer (WISE). Here we show that the average thermal inertia and geometric albedo of the investigated Vesta family members are 42… ▽ More In this work, we investigate the size, thermal inertia, surface roughness and geometric albedo of 10 Vesta family asteroids by using the Advanced Thermophysical Model (ATPM), based on the thermal infrared data acquired by mainly NASA's Wide-field Infrared Survey Explorer (WISE). Here we show that the average thermal inertia and geometric albedo of the investigated Vesta family members are 42 $\rm J m^{-2} s^{-1/2} K^{-1}$ and 0.314, respectively, where the derived effective diameters are less than 10 km. Moreover, the family members have a relatively low roughness fraction on their surfaces. The similarity in thermal inertia and geometric albedo among the V-type Vesta family member may reveal their close connection in the origin and evolution. As the fragments of the cratering event of Vesta, the family members may have undergone similar evolution process, thereby leading to very close thermal properties. Finally, we estimate their regolith grain sizes with different volume filling factors. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Comments: 29 pages, 40 figures, accepted for publication in AJ

Journal ref: The Astronomical Journal, Volume 159, Number 6, article id. 264 (2020)

Birth environment of circumbinary planets: are there CBPs on the inclined orbits?

Authors: Chuan-Tao Ma, Yan-Xiang Gong, Xiao-Mei Wu, Jianghui Ji

Abstract: The inclination distribution of circumbinary planets (CBPs) is an important scientific issue. It is of great significance in estimating the occurrence rate of CBPs and studying their formation and evolution. Although the CBPs currently discovered by the transit method are nearly coplanar, the true inclination distribution of CBPs is still unknown. Previous researches on CBPs mostly regarded them a… ▽ More The inclination distribution of circumbinary planets (CBPs) is an important scientific issue. It is of great significance in estimating the occurrence rate of CBPs and studying their formation and evolution. Although the CBPs currently discovered by the transit method are nearly coplanar, the true inclination distribution of CBPs is still unknown. Previous researches on CBPs mostly regarded them as an isolated binary-planet system, without considering the birth environment of their host binaries. It is generally believed that almost all stars are born in clusters. Therefore, it is necessary to consider the impact of the close encounters of stars on the CBP systems. This article discusses how the close encounters of fly-by stars affect the inclination of CBP. Based on extensive numerical simulations, we found that CBPs in close binary with a spacing of about 0.2 au are almost unaffected by star fly-bys. Their orbits remain coplanar. However, when the spacing of the binary stars is greater than 1 au, the 2-3 fly-bys of the intruding star can excite a considerable inclination even for the CBP near the unstable boundary of the binary. For the planets in the outer region, a single star fly-by can excite an inclination to more than 5 degrees. Especially, CBPs in near polar or retrograde orbits can naturally form through bianry-star encounters. If close binaries are born in open clusters, our simulations suggest that there may be high-inclination CBPs in binaries with spacing > 1 au. △ Less

Submitted 29 February, 2020; originally announced March 2020.

Comments: Published in MNRAS, vol. 493, p. 1907

Meso-scale Instability Triggered by Dust Feedback in Dusty Rings: Origin and Observational Implications

Authors: Pinghui Huang, Hui Li, Andrea Isella, Ryan Miranda, Shengtai Li, Jianghui Ji

Abstract: High spatial resolution observations of protoplanetary disks (PPDs) by ALMA have revealed many details that are providing interesting constraints on the disk physics as well as dust dynamics, both of which are essential for understanding planet formation. We carry out high-resolution, 2D global hydrodynamic simulations, including the effects of dust feedback, to study the stability of dusty rings.… ▽ More High spatial resolution observations of protoplanetary disks (PPDs) by ALMA have revealed many details that are providing interesting constraints on the disk physics as well as dust dynamics, both of which are essential for understanding planet formation. We carry out high-resolution, 2D global hydrodynamic simulations, including the effects of dust feedback, to study the stability of dusty rings. When the ring edges are relatively sharp and the dust surface density becomes comparable to the gas surface density, we find that dust feedback enhances the radial gradients of both the azimuthal velocity profile and the potential vorticity profile at the ring edges. This eventually leads to instabilities on meso-scales (spatial scales of several disk scale heights), causing dusty rings to be populated with many compact regions with highly concentrated dust densities on meso-scales. We also produce synthetic dust emission images using our simulation results and discuss the comparison between simulations and observations. △ Less

Submitted 20 January, 2020; originally announced January 2020.

Comments: 22 pages, 16 figures, Submitted to ApJ

The High Energy X-ray telescope (HE) onboard the Insight-HXMT astronomy satellite

Authors: C. Z. Liu, Y. F. Zhang, X. F. Li, X. F. Lu, Z. Chang, Z. W. Li, A. Z. Zhang, Y. J. Jin, H. M. Yu, Z. Zhang, M. X. Fu, Y. B. Chen, J. F. Ji, Y. P. Xu, J. K. Deng, R. C. Shang, G. Q. Liu, F. J. Lu, S. N. Zhang, Y. W. Dong, T. P. Li, M. Wu, Y. G. Li, H. Y. Wang, B. B. Wu , et al. (8 additional authors not shown)

Abstract: The Insight-Hard X-ray Modulation Telescope (Insight-HXMT) is a broad band X-ray and gamma-ray (1-3000 keV) astronomy satellite. The High Energy X-ray telescope (HE) is one of its three main telescopes. The main detector plane of HE is composed of 18 NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as primary detector to measure ~ 20-250 keV photons incident from the field of view (FOV) de… ▽ More The Insight-Hard X-ray Modulation Telescope (Insight-HXMT) is a broad band X-ray and gamma-ray (1-3000 keV) astronomy satellite. The High Energy X-ray telescope (HE) is one of its three main telescopes. The main detector plane of HE is composed of 18 NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as primary detector to measure ~ 20-250 keV photons incident from the field of view (FOV) defined by the collimators, and CsI(Na) is used as an active shield detector to NaI(Tl) by pulse shape discrimination. CsI(Na) is also used as an omnidirectional gamma-ray monitor. The HE collimators have a diverse FOV: 1.1°x 5.7° (15 units), 5.7°x 5.7° (2 units) and blocked (1 unit), thus the combined FOV of HE is about 5.7°x 5.7°. Each HE detector has a diameter of 190 mm, resulting in the total geometrical area of about 5100 cm_2. The energy resolution is ~15% at 60 keV. The timing accuracy is better than 10 μs and dead-time for each detector is less than 10 μs. HE is devoted to observe the spectra and temporal variability of X-ray sources in the 20-250 keV band either by pointing observations for known sources or scanning observations to unveil new sources, and to monitor the gamma-ray sky in 0.2-3 MeV. This paper presents the design and performance of the HE instruments. Results of the on-ground calibration experiments are also reported. △ Less

Submitted 10 October, 2019; originally announced October 2019.

Comments: Accepted by SCIENCE CHINA Physics, Mechanics & Astronomy

Revisiting Advanced Thermal Physical Model: New Perspectives on Thermophysical Characteristics of (341843) 2008 EV5 from Four-bands WISE Data with Sunlight-Reflection Model

Authors: Haoxuan Jiang, Liangliang Yu, Jianghui Ji

Abstract: In this paper, we investigate thermophysical characteristics of near-Earth asteroid (341843) 2008 EV5, based on our improved Advanced Thermal Physical Model (ATPM) by considering the contribution of sunlight-reflection for rough surface, along with four wavebands observations from Wide-field Infrared Survey Explorer (WISE) and the radar-derived shape model. Here we derive that 2008 EV5 has a relat… ▽ More In this paper, we investigate thermophysical characteristics of near-Earth asteroid (341843) 2008 EV5, based on our improved Advanced Thermal Physical Model (ATPM) by considering the contribution of sunlight-reflection for rough surface, along with four wavebands observations from Wide-field Infrared Survey Explorer (WISE) and the radar-derived shape model. Here we derive that 2008 EV5 has a relatively low thermal inertia of $Γ=110 _{-12}^{+40}\rm~J m^{-2} s^{-1/2} K^{-1}$ but a high roughness fraction. The geometric albedo and effective diameter are then constrained to be $p_v=0.095_{-0.003}^{+0.016}$ and $D_{\rm eff}=431_{-33}^{+6}\rm~m$, respectively. The low thermal inertia indicates that 2008 EV5 may have undergone sufficient space weathering over secular evolution. The high roughness may have resemblances to the appearances of Bennu and Ryugu recently observed by spacecrafts, where a great number of boulders are widely distributed on the asteroid's surface. Moreover, we numerically perform 1000 backward simulations of 2008 EV5's cloned orbits within 1 Myr to explore its origin, and present a probability of $\sim6.1\%$ that the asteroid originates from the main belt. Finally, we estimate that the mean grain size of the surface ranges from 0.58 to 1.3 mm, and infer that it is unlikely to find water ice on most area of 2008 EV5, but there may exist water ice on high-latitudes near polar region. △ Less

Submitted 21 September, 2019; originally announced September 2019.

Comments: 20 pages, 12 figures, accepted for publication in AJ

Journal ref: The Astronomical Journal, Volume 158, Number 5, article id. 205 (2019)

The Observability Of Vortex-Driven Spiral Arms In Protoplanetary Disk: Basic Spiral Properties

Authors: Pinghui Huang, Ruobing Dong, Hui Li, Shengtai Li, Jianghui Ji

Abstract: Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet-disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex… ▽ More Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet-disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex-driven spirals by the Rossby Wave Instability in 2D hydrodynamics simulations. We study how the surface density contrast, the number, and the shape of vortex-driven spirals depend on the properties of the vortex. We also compare vortex-driven spirals with those induced by planets. The surface density contrast of vortex-driven spirals in our simulations are comparable to those driven by a sub-thermal mass planet, typically a few to a few tens of Earth masses. In addition, different from the latter, the former is not sensitive to the mass of the vortex. Vortex-driven spiral arms are not expected to be detectable in current scattered light observations, and the prominent spirals observed in scattered light in a few protoplanetary disks, such as SAO 206462 (HD 135344B), MWC 758, and LkHa 330, are unlikely to be induced by the candidate vortices in them. △ Less

Submitted 2 September, 2019; originally announced September 2019.

Comments: 8 pages, 4 figures, Astrophysical Journal Letters accepted

In-orbit demonstration of X-ray pulsar navigation with the Insight-HXMT satellite

Authors: S. J. Zheng, S. N. Zhang, F. J. Lu, W. B. Wang, Y. Gao, T. P. Li, L. M. Song, M. Y. Ge, D. W. Han, Y. Chen, Y. P. Xu, X. L. Cao, C. Z. Liu, S. Zhang, J. L. Qu, Z. Chang, G. Chen, L. Chen, T. X. Chen, Y. B. Chen, Y. P. Chen, W. Cui, W. W. Cui, J. K. Deng, Y. W. Dong , et al. (91 additional authors not shown)

Abstract: In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and ob… ▽ More In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively. △ Less

Submitted 5 August, 2019; originally announced August 2019.

Comments: Accepted by the Astrophysical Journal Supplement

New constraints on the dust and gas distribution in the LkCa 15 disk from ALMA

Authors: Sheng Jin, Andrea Isella, Pinghui Huang, Shengtai Li, Hui Li, Jianghui Ji

Abstract: We search a large parameter space of the LkCa 15's disk density profile to fit its observed radial intensity profile of $^{12}$CO (J = 3-2) obtained from ALMA. The best-fit model within the parameter space has a disk mass of 0.1 $M_{\odot}$ (using an abundance ratio of $^{12}$CO/H$_2$ $=$ 1.4 $\times 10^{-4}$ in mass), an inner cavity of 45 AU in radius, an outer edge at $\sim$ 600 AU, and a disk… ▽ More We search a large parameter space of the LkCa 15's disk density profile to fit its observed radial intensity profile of $^{12}$CO (J = 3-2) obtained from ALMA. The best-fit model within the parameter space has a disk mass of 0.1 $M_{\odot}$ (using an abundance ratio of $^{12}$CO/H$_2$ $=$ 1.4 $\times 10^{-4}$ in mass), an inner cavity of 45 AU in radius, an outer edge at $\sim$ 600 AU, and a disk surface density profile follows a power-law of the form $ρ_r \propto r^{-4}$. For the disk density profiles that can lead to a small reduced $χ^2$ of goodness-of-fit, we find that there is a clear linear correlation between the disk mass and the power-law index $γ$ in the equation of disk density profile. This suggests that the $^{12}$CO disk of LkCa 15 is optically thick and we can fit its $^{12}$CO radial intensity profile using either a lower disk mass with a smaller $γ$ or a higher disk mass with a bigger $γ$. By comparing the $^{12}$CO channel maps of the best-fit model with disk models with higher or lower masses, we find that a disk mass of $\sim$ 0.1 $M_{\odot}$ can best reproduce the observed morphology of the $^{12}$CO channel maps. The dust continuum map at 0.87 mm of the LkCa 15 disk shows an inner cavity of the similar size of the best-fit gas model, but its out edge is at $\sim$ 200 AU, much smaller than the fitted gas disk. Such a discrepancy between the outer edges of the gas and dust disks is consistent with dust drifting and trapping models. △ Less

Submitted 1 July, 2019; originally announced July 2019.

Comments: 17 pages, 9 figures, accepted for publication in ApJ

3D analysis of spatial resolution of MIRO/Rosetta measurements at 67P/CG

Authors: L. Rezac, Y. Zhao, P. Hartogh, J. Ji, D. Marshall, X. Shi

Abstract: The MIRO instrument's remote sensing capability is integral in constraining water density, temperature and velocity fields in the coma of 67P/Churyumov-Gersimenko. Our aim is to quantify the contribution to the water density from all facets inside and outside the field-of-view (FOV) of MIRO, in both, nadir and limb geometries. This information is crucial for understanding the MIRO derived coma pro… ▽ More The MIRO instrument's remote sensing capability is integral in constraining water density, temperature and velocity fields in the coma of 67P/Churyumov-Gersimenko. Our aim is to quantify the contribution to the water density from all facets inside and outside the field-of-view (FOV) of MIRO, in both, nadir and limb geometries. This information is crucial for understanding the MIRO derived coma production rates and their relation to the nucleus characteristics, and inherent spatial resolution of the data. This study relies on a detailed 3D nucleus shape model, illumination conditions and the pointing information of the viewing geometry. With these parameters, we can evaluate the relative contribution of water density originating from facets directly inside the MIRO beam as well outside of the beam as a function of distance along the MIRO line-of-sight. We also calculate the ratio of in-beam versus out-of-beam number density. We demonstrate that despite the rather small MIRO field-of-view there is only a small fraction of molecules that originate from facets within the MIRO beam. This is true for nadir, but also translated into the limb observing geometry. The MIRO instrument cannot discriminate active from non-active regions directly from observations. This study also suggests that the beam averaged solar incidence angle, local time and mean normal vectors are not necessary related to molecules within the MIRO beam. These results also illustrate why the 1D spherical Haser model can be applied with relative success for analyzing the MIRO data (and generally any Rosetta measurements). The future possibilities of constraining gas activity distribution on the surface should use 3D codes extracting information from the MIRO spectral line shapes which contain additional information. The presented results are in fact applicable to all relevant instruments onboard Rosetta. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Comments: 16 pages, 9 figures

Constraining spatial pattern of early activity of comet 67P/C-G with 3D modeling of the MIRO observations

Authors: Y. Zhao, L. Rezac, P. Hartogh, J. Ji, R. Marschall, H. U. Keller

Abstract: Our aim is to investigate early activity (July 2014) of 67P/CG with 3D coma and radiative transfer modeling of MIRO measurements, accounting for nucleus shape, illumination, and orientation of the comet. We investigate MIRO line shape information for spatial distribution of water activity on the nucleus during the onset of activity. During this period we show that MIRO line shape have enough infor… ▽ More Our aim is to investigate early activity (July 2014) of 67P/CG with 3D coma and radiative transfer modeling of MIRO measurements, accounting for nucleus shape, illumination, and orientation of the comet. We investigate MIRO line shape information for spatial distribution of water activity on the nucleus during the onset of activity. During this period we show that MIRO line shape have enough information to clearly isolate contribution from Hapi and Inhotep independently, and compare it to the nominal case of activity from the entire illuminated surface. We also demonstrate that spectral line shapes differ from the 1D model for different viewing geometries and coma conditions relevant to this study. Specifically, line shapes are somewhat sensitive to the location of the terminator in the coma. At last, fitting the MIRO observations we show that the Imhotep region (possible extended source of H$_{2}$O due to CO$_{2}$ activity) contributes only a small fraction of the total number of water molecules into MIRO beam in the early activity. On the other hand, a strong enhancement of water activity from the Hapi region seems required to fit the MIRO line shapes. This is consistent with earlier Rosetta results. Nevertheless, within the assumption of our coma and surface boundary conditions, we cannot get a reasonable fit to all MIRO mapping observations in July 2014, which may illustrate that a more sophisticated coma model or more accurate temperature/velocity distribution is needed. △ Less

Submitted 27 February, 2019; originally announced February 2019.

Comments: 22 pages, 15 figures, submitted

Dynamical evolution and stability maps of the Proxima Centauri system

Authors: Tong Meng, Jianghui Ji, Yao Dong

Abstract: Proxima Centauri was recently discovered to host an Earth-mass planet of Proxima b, and a 215-day signal which is probably a potential planet c. In this work, we investigate the dynamical evolution of the Proxima Centauri system with the full equations of motion and semi-analytical models including relativistic and tidal effects. We adopt the modified Lagrange-Laplace secular equations to study th… ▽ More Proxima Centauri was recently discovered to host an Earth-mass planet of Proxima b, and a 215-day signal which is probably a potential planet c. In this work, we investigate the dynamical evolution of the Proxima Centauri system with the full equations of motion and semi-analytical models including relativistic and tidal effects. We adopt the modified Lagrange-Laplace secular equations to study the evolution of eccentricity of Proxima b, and find that the outcomes are consistent with those from the numerical simulations. The simulations show that relativistic effects have an influence on the evolution of eccentricities of planetary orbits, whereas tidal effects primarily affects the eccentricity of Proxima b over long timescale. Moreover, using the MEGNO (the Mean Exponential Growth factor of Nearby Orbits) technique, we place dynamical constraints on orbital parameters that result in stable or quasi-periodic motions for coplanar and non-coplanar configurations. In the coplanar case, we find that the orbit of Proxima b is stable for the semi-major axis ranging from 0.02 au to 0.1 au and the eccentricity being less than 0.4. This is where the best-fitting parameters for Proxima b exactly fall. Additional simulations show that the robust stability of this system would favor an eccentricity of Proxima b less than 0.45 and that of Proxima c below 0.65. In the non-coplanar case, we find that mutual inclinations of two planets must be lower than $50^{\circ}$ in order to provide stability. Finally, we estimate the mass of Proxima c to be $3.13~{M_ \oplus } \le {m_c} \le 70.7~{M_ \oplus }$ when $1.27~{M_ \oplus } \le {m_b} \le 1.6~{M_ \oplus }$, if ${i_{mutual}} \le {50^{\circ}}$ and $ΔΩ= 0^{\circ}$. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 482, 372-383 (2019)

Identifying Anticyclonic Vortex Features Produced by the Rossby Wave Instability in Protoplanetary Disks

Authors: Pinghui Huang, Andrea Isella, Hui Li, Shengtai Li, Jianghui Ji

Abstract: Several nearby protoplanetary disks have been observed to display large scale crescents in the (sub)millimeter dust continuum emission. One interpretation is that these structures correspond to anticyclonic vortices generated by the Rossby wave instability within the gaseous disk. Such vortices have local gas over-densities and are expected to concentrate dust particles with Stokes number around u… ▽ More Several nearby protoplanetary disks have been observed to display large scale crescents in the (sub)millimeter dust continuum emission. One interpretation is that these structures correspond to anticyclonic vortices generated by the Rossby wave instability within the gaseous disk. Such vortices have local gas over-densities and are expected to concentrate dust particles with Stokes number around unity. This process might catalyze the formation of planetesimals. Whereas recent observations showed that dust crescent are indeed regions where millimeter-size particles have abnormally high concentration relative to the gas and smaller grains, no observations have yet shown that the gas within the crescent region counter-rotates with respect to the protoplanetary disk. Here we investigate the detectability of anticyclonic features through measurement of the line-of-sight component of the gas velocity obtained with ALMA. We carry out 2D hydrodynamic simulations and 3D radiative transfer calculation of a protoplanetary disk characterized by a vortex created by the tidal interaction with a massive planet. As a case study, the disk parameters are chosen to mimic the IRS 48 system, which has the most prominent crescent observed to date. We generate synthetic ALMA observations of both the dust continuum and 12CO emission around the frequency of 345 GHz. We find that the anticyclonic features of vortex are weak but can be detected if both the source and the observational setup are properly chosen. We provide a recipe for maximizing the probability to detect such vortex features and present an analysis procedure to infer their kinematic properties. △ Less

Submitted 18 September, 2018; originally announced September 2018.

Comments: 14 pages, 8 figures, Accepted for publication in Astrophysical Journal

Insight-HXMT observations of the New Black Hole Candidate MAXI J1535-571: timing analysis

Authors: Y. Huang, J. L. Qu, S. N. Zhang, Q. C. Bu, Y. P. Chen, L. Tao, S. Zhang, F. J. Lu, T. P. Li, L. M. Song, Y. P. Xu, X. L. Cao, Y. Chen, C. Z. Liu, H. -K. Chang, W. f. Yu, S. S. Weng, X. Hou, A. K. H. Kong, F. G. Xie, G. B. Zhang, J. F. ZHOU, Z. Chang, G. Chen, L. Chen , et al. (98 additional authors not shown)

Abstract: We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventual… ▽ More We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system. △ Less

Submitted 15 August, 2018; originally announced August 2018.

Comments: 12 pages, 11 figures, Sumbitted to ApJ

Formation of S-type planets in close binaries: scattering induced tidal capture of circumbinary planets

Authors: Yan-Xiang Gong, Jianghui Ji

Abstract: Although several S-type and P-type planets in binary systems were discovered in past years, S-type planets have not yet been found in close binaries with an orbital separation not more than 5 au. Recent studies suggest that S-type planets in close binaries may be detected through high-accuracy observations. However, nowadays planet formation theories imply that it is difficult for S-type planets i… ▽ More Although several S-type and P-type planets in binary systems were discovered in past years, S-type planets have not yet been found in close binaries with an orbital separation not more than 5 au. Recent studies suggest that S-type planets in close binaries may be detected through high-accuracy observations. However, nowadays planet formation theories imply that it is difficult for S-type planets in close binaries systems to form in situ. In this work, we extensively perform numerical simulations to explore scenarios of planet-planet scattering among circumbinary planets and subsequent tidal capture in various binary configurations, to examine whether the mechanism can play a part in producing such kind of planets. Our results show that this mechanism is robust. The maximum capture probability is $\sim 10\%$, which can be comparable to the tidal capture probability of hot Jupiters in single star systems. The capture probability is related to binary configurations, where a smaller eccentricity or a low mass ratio of the binary will lead to a larger probability of capture, and vice versa. Furthermore, we find that S-type planets with retrograde orbits can be naturally produced via capture process. These planets on retrograde orbits can help us distinguish in situ formation and post-capture origin for S-type planet in close binaries systems. The forthcoming missions (PLATO) will provide the opportunity and feasibility to detect such planets. Our work provides several suggestions for selecting target binaries in search for S-type planets in the near future. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 478, 4565-4574 (2018)

Investigation of Thermal Inertia and Surface Properties for Near-Earth Asteroid (162173) 1999 JU3

Authors: Liangliang Yu, Jianghui Ji, Su Wang

Abstract: In order to obtain the substantial information about the surface physics and thermal property of the target asteroid (162173) 1999 JU3, which will be visited by Hayabusa 2 in a sample return mission, with the Advanced Thermal Physical Model (ATPM) we estimate the possible thermal inertia distribution over its surface, and infer the major material composition of its surface materials. In addition,… ▽ More In order to obtain the substantial information about the surface physics and thermal property of the target asteroid (162173) 1999 JU3, which will be visited by Hayabusa 2 in a sample return mission, with the Advanced Thermal Physical Model (ATPM) we estimate the possible thermal inertia distribution over its surface, and infer the major material composition of its surface materials. In addition, the effective diameter and geometric albedo are derived to be $D_{\rm eff}=1.13\pm0.03\rm~km$, $p_{\rm v}=0.042\pm0.003$, respectively, and the average thermal inertia is estimated to be about $(300\pm50)\rm~J\cdot m^{-2}\cdot s^{-0.5}\cdot K^{-1}$. According to the derived thermal inertia distribution, we infer that the major area on the surface of the target asteroid may be covered by loose materials, such as rock debris, sands, and so on, but few bare rocks may exist in a very small region. In this sense, the sample return mission of Hayabusa 2 is feasible, when it is performed successfully, it will certainly bring significant scientific information to the research of asteroids. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 15 pages, 7 figures, published in Chinese Astronomy and Astrophysics

Tidal Evolution of Kepler lower-mass Planets

Authors: Yao Dong, Jianghui Ji, Su Wang

Abstract: The planets with a radius $<$ 4 $R$$_\oplus$ observed by the Kepler mission exhibit a unique feature, and propose a challenge for current planetary formation models. The tidal effect between a planet and its host star plays an essential role in reconfiguring the final orbits of the short-period planets. In this work, based on various initial Rayleigh distributions of the orbital elements, the fina… ▽ More The planets with a radius $<$ 4 $R$$_\oplus$ observed by the Kepler mission exhibit a unique feature, and propose a challenge for current planetary formation models. The tidal effect between a planet and its host star plays an essential role in reconfiguring the final orbits of the short-period planets. In this work, based on various initial Rayleigh distributions of the orbital elements, the final semi-major axis distributions of the planets with a radius $<$ 4 $R_\oplus$ after suffering tidal evolutions are investigated. Our simulations have qualitatively revealed some statistical properties: the semi-major axis and its peak value all increase with the increase of the initial semi-major axis and eccentricity. For the case that the initial mean semi-major axis is less than 0.1 au and the mean eccentricity is larger than 0.25, the results of numerical simulation are approximately consistent with the observation. In addition, the effects of other parameters, such as the tidal dissipation coefficient, stellar mass and planetary mass, etc., on the final semi-major axis distribution after tidal evolution are all relatively small. Based on the simulation results, we have tried to find some clues for the formation mechanism of low-mass planets. We speculate that these low-mass planets probably form in the far place of protoplanetary disk with a moderate eccentricity via the type I migration, and it is also possible to form in situ. △ Less

Submitted 13 May, 2018; originally announced May 2018.

Comments: 14 pages, 10 figures, to appear in Chinese Astronomy and Astrophysics

Photometric Follow-up Transit (Primary Eclipse) Observations of WASP-43 b and TrES-3b and A Study on Their Transit Timing Variations

Authors: Zhao Sun, Jianghui Ji, Yao Dong

Abstract: Two photometric follow-up transit (primary eclipse) observations on WASP-43 b and four observations on TrES-3 b are performed using the Xuyi Near-Earth Object Survey Telescope. After differential photometry and light curve analysis, the physical parameters of the two systems are obtained and are in good match with the literature. Combining with transit data from a lot of literature, the residuals… ▽ More Two photometric follow-up transit (primary eclipse) observations on WASP-43 b and four observations on TrES-3 b are performed using the Xuyi Near-Earth Object Survey Telescope. After differential photometry and light curve analysis, the physical parameters of the two systems are obtained and are in good match with the literature. Combining with transit data from a lot of literature, the residuals ($O-C$) of transit observations of both systems are fitted with the linear and quadratic functions. With the linear fitting, the periods and transit timing variations (TTVs) of the planets are obtained, and no obvious periodic TTV signal is found in both systems after an analysis. The maximum mass of a perturbing planet located at the 1:2 mean motion resonance (MMR) for WASP-43 b and TrES-3 b is estimated to be 1.826 and 1.504 Earth mass, respectively. By quadratic fitting, it is confirmed that WASP-43 b may have a long-term TTV which means an orbital decay. The decay rate is shown to be $\dot{P} =(-0.005248 \pm 0.001714)$~s$\cdot$yr$^{-1}$, and compared with the previous results. Based on this, the lower limit of the stellar tidal quality parameter of WASP-43 is calculated to be $Q'_{\rm *} \ge 1.5\times10^5$, and the remaining lifetimes of the planets are presented for the different $Q'_{\rm *}$ values of the two systems, correspondingly. △ Less

Submitted 24 April, 2018; originally announced April 2018.

Comments: 28 pages, 6 figures, published in Chinese Astronomy and Astrophysics

Journal ref: Chinese Astronomy and Astrophysics 42 (2018) 101-128

The formation mechanism of 4179 Toutatis' elongated bi-lobed structure in a close Earth encounter scenario

Authors: Shoucun Hu, Jianghui Ji, Derek C. Richardson, Yuhui Zhao, Yun Zhang

Abstract: The optical images of near-Earth asteroid 4179 Toutatis acquired by Chang'e-2 spacecraft show that Toutatis has an elongated contact binary configuration, with the contact point located along the long axis. We speculate that such configuration may have resulted from a low-speed impact between two components. In this work, we performed a series of numerical simulations and compared the results with… ▽ More The optical images of near-Earth asteroid 4179 Toutatis acquired by Chang'e-2 spacecraft show that Toutatis has an elongated contact binary configuration, with the contact point located along the long axis. We speculate that such configuration may have resulted from a low-speed impact between two components. In this work, we performed a series of numerical simulations and compared the results with the optical images, to examine the mechanism and better understand the formation of Toutatis. Herein we propose an scenario that an assumed separated binary precursor could undergo a close encounter with Earth, leading to an impact between the primary and secondary, and the elongation is caused by Earth's tide. The precursor is assumed to be a doubly synchronous binary with a semi-major axis of 4 $R_p$ (radius of primary) and the two components are represented as spherical cohesionless self-gravitating granular aggregates. The mutual orbits are simulated in a Monte Carlo routine to provide appropriate parameters for our \emph{N}-body simulations of impact and tidal distortion. We employ the \emph{pkdgrav} package with a soft-sphere discrete element method (SSDEM) to explore the entire scenarios. The results show that contact binary configurations are natural outcomes under this scenario, whereas the shape of the primary is almost not affected by the impact of the secondary. However, our simulations further provide an elongated contact binary configuration best-matching to the shape of Toutatis at an approaching distance $r_p$ = 1.4 $\sim$ 1.5 $R_e$ (Earth radius), indicative of a likely formation scenario for configurations of Toutatis-like elongated contact binaries. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: 15 pages, 12 figures, accepted to MNRAS

Limits on radio emission from meteors using the MWA

Authors: Xiang Zhang, Paul Hancock, Hadrien A. R. Devillepoix, Randall B. Wayth, A. Beardsley, B. Crosse, D. Emrich, T. M. O. Franzen, B. M. Gaensler, L. Horsley, M. Johnston-Hollitt, D. L. Kaplan, D. Kenney, M. F. Morales, D. Pallot, K. Steele, S. J. Tingay, C. M. Trott, M. Walker, A. Williams, C. Wu, Jianghui Ji, Yuehua Ma

Abstract: Recently, low frequency, broadband radio emission has been observed accompanying bright meteors by the Long Wavelength Array (LWA). The broadband spectra between 20 and 60 MHz were captured for several events, while the spectral index (dependence of flux density on frequency, with $S_ν\propto ν^α$) was estimated to be $-4\pm1$ during the peak of meteor afterglows. Here we present a survey of meteo… ▽ More Recently, low frequency, broadband radio emission has been observed accompanying bright meteors by the Long Wavelength Array (LWA). The broadband spectra between 20 and 60 MHz were captured for several events, while the spectral index (dependence of flux density on frequency, with $S_ν\propto ν^α$) was estimated to be $-4\pm1$ during the peak of meteor afterglows. Here we present a survey of meteor emission and other transient events using the Murchison Widefield Array (MWA) at 72-103 MHz. In our 322-hour survey, down to a $5σ$ detection threshold of 3.5 Jy/beam, no transient candidates were identified as intrinsic emission from meteors. We derived an upper limit of -3.7 (95% confidence limit) on the spectral index in our frequency range. We also report detections of other transient events, like reflected FM broadcast signals from small satellites, conclusively demonstrating the ability of the MWA to detect and track space debris on scales as small as 0.1 m in low Earth orbits. △ Less

Submitted 19 April, 2018; originally announced April 2018.

Comments: Accepted for Publication in MNRAS on 11 April 2018, 11 pages, 11 figures, 5 tables