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

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

Near Mean-motion Resonances in the Systems Observed by Kepler: Affected by Mass Accretion and Type I Migration

Authors: Su Wang, Jianghui Ji

Abstract: The Kepler mission has released over 4496 planetary candidates, among which 3483 planets have been confirmed as of April 2017. The statistical results of the planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. The observations indicate that a plenty of planet pairs could have firstly been captured into mean motion resonances (MMRs) in planetary fo… ▽ More The Kepler mission has released over 4496 planetary candidates, among which 3483 planets have been confirmed as of April 2017. The statistical results of the planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. The observations indicate that a plenty of planet pairs could have firstly been captured into mean motion resonances (MMRs) in planetary formation. Subsequently, these planets depart from exact resonant locations to be near MMRs configurations. Through type I migration, two low-mass planets have a tendency to be trapped into first-order MMRs (2:1 or 3:2 MMRs), however two scenarios of mass accretion of planets and potential outward migration play an important role in reshaping their final orbital configurations. Under the scenario of mass accretion, the planet pairs can cross 2:1 MMRs and then enter into 3:2 MMRs, especially for the inner pairs. With such formation scenario, the possibility that two planets are locked into 3:2 MMRs can increase if they are formed in a flat disk. Moreover, the outward migration can make planets have a high likelihood to be trapped into 3:2 MMRs. We perform additional runs to investigate the mass relationship for those planets in three-planet systems, and we show that two peaks near 1.5 and 2.0 for the period ratios of two planets can be easily reproduced through our formation scenario. We further show that the systems in chain resonances (e.g., 4:2:1, 3:2:1, 6:3:2 and 9:6:4 MMRs), have been observed in our simulations. This mechanism can be applicable to understand the formation of systems of Kepler-48, Kepler-53, Kepler-100, Kepler-192, Kepler-297, Kepler-399, and Kepler-450. △ Less

Submitted 7 October, 2017; originally announced October 2017.

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

The scattering outcomes of Kepler circumbinary planets: planet mass ratio

Authors: Yan-Xiang Gong, Jianghui Ji

Abstract: Recent studies reveal that the free eccentricities of Kepler-34b and Kepler-413b are much larger than their forced eccentricities, implying that the scattering events may take place in their formation. The observed orbital configuration of Kepler-34b cannot be well reproduced in disk-driven migration models, whereas a two-planet scattering scenario can play a significant role of shaping the planet… ▽ More Recent studies reveal that the free eccentricities of Kepler-34b and Kepler-413b are much larger than their forced eccentricities, implying that the scattering events may take place in their formation. The observed orbital configuration of Kepler-34b cannot be well reproduced in disk-driven migration models, whereas a two-planet scattering scenario can play a significant role of shaping the planetary configuration. These studies indicate that circumbinary planets discovered by Kepler may have experienced scattering process. In this work, we extensively investigate the scattering outcomes of circumbinary planets focusing on the effects of planet mass ratio. We find that the planetary mass ratio and the the initial relative locations of planets act as two important parameters which affect the eccentricity distribution of the surviving planets. As an application of our model, we discuss the observed orbital configurations of Kepler-34b and Kepler-413b. We first adopt the results from the disk-driven models as the initial conditions, then simulate the scattering process occurred in the late evolution stage of circumbinary planets. We show that the present orbital configurations of Kepler-34b and Kepler-413b can be well reproduced when considering two unequal-mass planet ejection model. Our work further suggests that some of the currently discovered circumbinary single-planet systems may be the survivals of original multiple-planet systems. The disk-driven migration and the scattering events occurring in the late stage both play an irreplaceable role in sculpting the final systems. △ Less

Submitted 13 September, 2017; originally announced September 2017.

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

Journal ref: AJ, 154:179 (2017)

Using Chebyshev polynomials interpolation to improve the computation efficiency of gravity near an irregular-shaped asteroid

Authors: Shoucun Hu, Jianghui Ji

Abstract: In asteroid rendezvous missions, the dynamical environment near the asteroid's surface should be made clear prior to the mission launch. However, most of the asteroids have irregular shapes, which lower the efficiency of calculating their gravitational field by adopting the traditional polyhedral method. In this work, we propose a method to partition the space near the asteroid adaptively along th… ▽ More In asteroid rendezvous missions, the dynamical environment near the asteroid's surface should be made clear prior to the mission launch. However, most of the asteroids have irregular shapes, which lower the efficiency of calculating their gravitational field by adopting the traditional polyhedral method. In this work, we propose a method to partition the space near the asteroid adaptively along three spherical coordinates and use Chebyshev polynomials interpolation to represent the gravitational acceleration in each cell. Moreover, we compare four different interpolation schemes to obtain the best precision with the identical initial parameters. An error-adaptive octree division is combined to improve the interpolation precision near the surface. As an example, we take the typical irregular-shaped near-Earth asteroid 4179 Toutatis to show the advantage of this method, as a result, we show that the efficiency can be increased by hundreds to thousands times with our method. In a word, this method can be applicable to other irregular-shaped asteroids and can greatly improve the evaluation efficiency. △ Less

Submitted 22 August, 2017; originally announced August 2017.

Comments: 17 pages, 10 figures, 3 tables, accepted to Research in Astronomy and Astrophysics

Thermophysical characteristics of the large main-belt asteroid (349) Dembowska

Authors: LiangLiang Yu, Bin Yang, Jianghui Ji, Wing-Huen Ip

Abstract: (349) Dembowska, a large, bright main-belt asteroid, has a fast rotation and oblique spin axis. It may have experienced partial melting and differentiation. We constrain Dembowska's thermophysical properties, e.g., thermal inertia, roughness fraction, geometric albedo and effective diameter within 3$σ$ uncertainty of $Γ=20^{+12}_{-7}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.25^{+0.60}_{-0.25}$,… ▽ More (349) Dembowska, a large, bright main-belt asteroid, has a fast rotation and oblique spin axis. It may have experienced partial melting and differentiation. We constrain Dembowska's thermophysical properties, e.g., thermal inertia, roughness fraction, geometric albedo and effective diameter within 3$σ$ uncertainty of $Γ=20^{+12}_{-7}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.25^{+0.60}_{-0.25}$, $p_{\rm v}=0.309^{+0.026}_{-0.038}$, and $D_{\rm eff}=155.8^{+7.5}_{-6.2}\rm~km$, by utilizing the Advanced Thermophysical Model (ATPM) to analyse four sets of thermal infrared data obtained by IRAS, AKARI, WISE and Subaru/COMICS at different epochs. In addition, by modeling the thermal lightcurve observed by WISE, we obtain the rotational phases of each dataset. These rotationally resolved data do not reveal significant variations of thermal inertia and roughness across the surface, indicating the surface of Dembowska should be covered by a dusty regolith layer with few rocks or boulders. Besides, the low thermal inertia of Dembowska show no significant difference with other asteroids larger than 100 km, indicating the dynamical lives of these large asteroids are long enough to make the surface to have sufficiently low thermal inertia. Furthermore, based on the derived surface thermophysical properties, as well as the known orbital and rotational parameters, we can simulate Dembowska's surface and subsurface temperature throughout its orbital period. The surface temperature varies from $\sim40$ K to $\sim220$ K, showing significant seasonal variation, whereas the subsurface temperature achieves equilibrium temperature about $120\sim160$ K below $30\sim50$ cm depth. △ Less

Submitted 10 August, 2017; originally announced August 2017.

Comments: 10 pages, 14 figures, accepted to MNRAS

Journal ref: MNRAS, 472, 2388-2397 (2017)

Surface thermophysical properties investigation of the potentially hazardous asteroid (99942) Apophis

Authors: LiangLiang Yu, Jianghui Ji, Wing-Huen Ip

Abstract: In this work, we investigate the surface thermophysical properties (thermal emissivity, thermal inertia, roughness fraction and geometric albedo) of asteroid (99942) Apophis, using the currently available thermal infrared observations of CanariCam on Gran Telescopio CANARIAS and far-infrared data by PACS of Herschel, on the basis of the Advanced thermophysical model. We show that the thermal emiss… ▽ More In this work, we investigate the surface thermophysical properties (thermal emissivity, thermal inertia, roughness fraction and geometric albedo) of asteroid (99942) Apophis, using the currently available thermal infrared observations of CanariCam on Gran Telescopio CANARIAS and far-infrared data by PACS of Herschel, on the basis of the Advanced thermophysical model. We show that the thermal emissivity of Apophis should be wavelength dependent from $8.70~μm$ to $160~μm$, and the maximum emissivity may arise around $20~μm$ similar to that of Vesta. Moreover, we further derive the thermal inertia, roughness fraction, geometric albedo and effective diameter of Apophis within a possible 1$σ$ scale of $Γ=100^{+100}_{-52}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.78\sim1.0$, $p_{\rm v}=0.286^{+0.030}_{-0.026}$, $D_{\rm eff}=378^{+19}_{-25}\rm~m$, and 3$σ$ scale of $Γ=100^{+240}_{-100}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.2\sim1.0$, $p_{\rm v}=0.286^{+0.039}_{-0.029}$, $D_{\rm eff}=378^{+27}_{-29}\rm~m$. The derived low thermal inertia but high roughness fraction may imply that Apophis could have regolith on its surface, and less regolith migration process has happened in comparison with asteroid Itokawa. Our results show that small-size asteroids could also have fine regolith on the surface, and further infer that Apophis may be delivered from the Main Belt by Yarkovsky effect. △ Less

Submitted 20 March, 2017; originally announced March 2017.

Comments: 15 pages, 6 figures, 6 tables, accepted to Research in Astronomy and Astrophysics

Journal ref: Research in Astronomy and Astrophysics, 17, 070 (2017)

Terrestrial Planets Formation under Migration: the Systems near 4:2:1 Mean Motion Resonance

Authors: Zhao Sun, Jianghui Ji, Su Wang, Sheng Jin

Abstract: In this work, we extensively investigate the formation of near 4:2:1 mean motion resonances (MMRs) configuration by performing two sets of N-body simulations. We model the eccentricity damping, gas drag, type I and type II planetary migration of planetesimals, planetary embryos and giant planets in the first sets. For the simulations of giant planets with type II migration, the massive terrestrial… ▽ More In this work, we extensively investigate the formation of near 4:2:1 mean motion resonances (MMRs) configuration by performing two sets of N-body simulations. We model the eccentricity damping, gas drag, type I and type II planetary migration of planetesimals, planetary embryos and giant planets in the first sets. For the simulations of giant planets with type II migration, the massive terrestrial planets, with a mass up to several Earth masses, are likely produced in the systems. We further show that by shepherding and/or scattering mechanisms through Jovian planet's type II migration, the terrestrial planets and giant planets in the systems can be evolved into a near 4:2:1 MMRs. Moreover, the models are applicable to the formation of Kepler-238 and 302 systems. In the second set, we study the 4:2:1 MMRs formation in the terrestrial planetary systems, where the planets undergo type I migration and eccentricity damping. By considering type I migration, $\sim$ 17.1\% of the simulations indicate that terrestrial planets are evolved into 4:2:1 MMRs. However, this probability should depend on the initial conditions of planets. Hence, we conclude that both type I and type II migration can play a crucial role in close-in terrestrial planet formation. △ Less

Submitted 11 January, 2017; originally announced January 2017.

Comments: 18 pages, 9 figures, 3 tables, accepted to MNRAS

Journal ref: MNRAS, 467, 619-632 (2017)

Modeling Dust Emission of HL Tau Disk Based on Planet-Disk Interactions

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

Abstract: We use extensive global two-dimensional hydrodynamic disk gas+dust simulations with embedded planets, coupled with three dimensional radiative transfer calculations, to model the dust ring and gap structures in the HL Tau protoplanetary disk observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We include the self-gravity of disk gas and dust components and make reasonable choices… ▽ More We use extensive global two-dimensional hydrodynamic disk gas+dust simulations with embedded planets, coupled with three dimensional radiative transfer calculations, to model the dust ring and gap structures in the HL Tau protoplanetary disk observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We include the self-gravity of disk gas and dust components and make reasonable choices of disk parameters, assuming an already settled dust distribution and no planet migration. We can obtain quite adequate fits to the observed dust emission using three planets with masses 0.35, 0.17, and 0.26 $M_{Jup}$ at 13.1, 33.0, and 68.6 AU, respectively. Implications for the planet formation as well as the limitations of this scenario are discussed. △ Less

Submitted 3 January, 2016; originally announced January 2016.

Comments: Accepted for publication in ApJ

TianQin: a space-borne gravitational wave detector

Authors: Jun Luo, Li-Sheng Chen, Hui-Zong Duan, Yun-Gui Gong, Shoucun Hu, Jianghui Ji, Qi Liu, Jianwei Mei, Vadim Milyukov, Mikhail Sazhin, Cheng-Gang Shao, Viktor T. Toth, Hai-Bo Tu, Yamin Wang, Yan Wang, Hsien-Chi Yeh, Ming-Sheng Zhan, Yonghe Zhang, Vladimir Zharov, Ze-Bing Zhou

Abstract: TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gr… ▽ More TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gravitational waves within a few months of observing time. We describe the preliminary mission concept for TianQin, including the candidate source and experimental designs. We present estimates for the major constituents of the experiment's error budget and discuss the project's overall feasibility. Given the current level of technology readiness, we expect TianQin to be flown in the second half of the next decade. △ Less

Submitted 7 December, 2015; originally announced December 2015.

Comments: 32 pages, 4 figures, to be published in Class. Quantum Grav

Revisit of rotational dynamics of Asteroid 4179 Toutatis from Chang'e-2's flyby

Authors: Yuhui Zhao, Shoucun Hu, Jianghui Ji

Abstract: This paper presents analysis of the rotational parameters of Toutatis based on the observational results from Chang'e-2's close flyby. The 3-D shape model derived from ground-based radar observation is used to calculate the 3-1-3 Euler angles at the flyby epoch, which are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$. The large amplitude of Toutatis' tu… ▽ More This paper presents analysis of the rotational parameters of Toutatis based on the observational results from Chang'e-2's close flyby. The 3-D shape model derived from ground-based radar observation is used to calculate the 3-1-3 Euler angles at the flyby epoch, which are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$. The large amplitude of Toutatis' tumbling attitude is demonstrated to be the result of the large deviation of the angular momentum axis and the rotational axis. Two rotational periods are evaluated to be $5.38\pm0.03$ days for rotation about the long axis and $7.40\pm0.03$ days for precession of the long axis about the angular momentum vector based on Fourier analysis. These results provide a further understanding of rotational state of Toutatis. △ Less

Submitted 13 November, 2015; originally announced November 2015.

Comments: 4 pages, 3 figures, accepted to IAU Symposium 318: Asteroids: New Observations, New Models

Journal ref: IAU Symposium 318, (2016)156-159

Thermophysical Characteristics of OSIRIS-REx Target Asteroid (101955) Bennu

Authors: LiangLiang Yu, Jianghui Ji

Abstract: In this work, we investigate the thermophysical properties, including thermal inertia, roughness fraction and surface grain size of OSIRIS-REx target asteroid (101955) Bennu by using a thermophysical model with the recently updated 3D radar-derived shape model (\cite[Nolan et al., 2013]{Nolan2013}) and mid-infrared observations (\cite[M$\ddot{u}$ller et al, 2012]{Muller2012}, \cite[Emery et al., 2… ▽ More In this work, we investigate the thermophysical properties, including thermal inertia, roughness fraction and surface grain size of OSIRIS-REx target asteroid (101955) Bennu by using a thermophysical model with the recently updated 3D radar-derived shape model (\cite[Nolan et al., 2013]{Nolan2013}) and mid-infrared observations (\cite[M$\ddot{u}$ller et al, 2012]{Muller2012}, \cite[Emery et al., 2014]{Emery2014}). We find that the asteroid bears an effective diameter of $510^{+6}_{-40}$ m, a geometric albedo of $0.047^{+0.0083}_{-0.0011}$, a roughness fraction of $0.04^{+0.26}_{-0.04}$, and thermal inertia of $240^{+440}_{-60}\rm~Jm^{-2}s^{-0.5}K^{-1}$ for our best-fit solution. The best-estimate thermal inertia suggests that fine-grained regolith may cover a large portion of Bennu's surface, where a grain size may vary from $1.3$ to $31$~mm. Our outcome suggests that Bennu is suitable for the OSIRIS-REx mission to return samples to Earth. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 3 pages, 1 figures, accepted to IAU Symposium 318: Asteroids: New Observations, New Models

Journal ref: IAU Symposium 318, (2016)218-220

Chang'e-2 spacecraft observations of asteroid 4179 Toutatis

Authors: Jianghui Ji, Yun Jiang, Yuhui Zhao, Su Wang, Liangliang Yu

Abstract: On 13 December 2012, Chang'e-2 completed a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 meters from the asteroid's surface. The observations show that Toutatis has an irregular surface and its shape resembles a ginger-root of a smaller lobe (head) and a larger lobe (body). Such bilobate shape is indicative of a contact binary origin for Toutatis. In additi… ▽ More On 13 December 2012, Chang'e-2 completed a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 meters from the asteroid's surface. The observations show that Toutatis has an irregular surface and its shape resembles a ginger-root of a smaller lobe (head) and a larger lobe (body). Such bilobate shape is indicative of a contact binary origin for Toutatis. In addition, the high-resolution images better than 3 meters provide a number of new discoveries about this asteroid, such as an 800-meter depression at the end of the large lobe, a sharply perpendicular silhouette near the neck region, boulders, indicating that Toutatis is probably a rubble-pile asteroid. Chang'e-2 observations have significantly revealed new insights into the geological features and the formation and evolution of this asteroid. In final, we brief the future Chinese asteroid mission concept. △ Less

Submitted 6 November, 2015; originally announced November 2015.

Comments: 9 pages, 5 figures, accepted to IAU Symposium 318: Asteroids: New Observations, New Models

Journal ref: IAU Symposium 318, (2016)144-152