Showing 1–50 of 55 results for author: Cai, M

The Host Galaxy Fluxes of Active Galaxy Nuclei Are Generally Overestimated by the Flux Variation Gradient Method

Authors: Minxuan Cai, Zhen Wan, Zhenyi Cai, Lulu Fan, Junxian Wang

Abstract: In terms of the variable nature of normal active galaxy nuclei (AGN) and luminous quasars, a so-called flux variation gradient (FVG) method has been widely utilized to estimate the underlying non-variable host galaxy fluxes. The FVG method assumes an invariable AGN color, but this assumption has been questioned by the intrinsic color variation of quasars and local Seyfert galaxies. Here, using an… ▽ More In terms of the variable nature of normal active galaxy nuclei (AGN) and luminous quasars, a so-called flux variation gradient (FVG) method has been widely utilized to estimate the underlying non-variable host galaxy fluxes. The FVG method assumes an invariable AGN color, but this assumption has been questioned by the intrinsic color variation of quasars and local Seyfert galaxies. Here, using an up-to-date thermal fluctuation model to simulate multi-wavelength AGN variability, we theoretically demonstrate that the FVG method generally overestimates the host galaxy flux; that is, it is more significant for brighter AGN/quasars. Furthermore, we observationally confirm that the FVG method indeed overestimates the host galaxy flux by comparing it to that estimated through other independent methods. We thus caution that applying the FVG method should be performed carefully in the era of time-domain astronomy. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Journal ref: Universe 2024, 10, 282

AT2023lli: A Tidal Disruption Event with Prominent Optical Early Bump and Delayed Episodic X-ray Emission

Authors: Shifeng Huang, Ning Jiang, Jiazheng Zhu, Yibo Wang, Tinggui Wang, Shan-Qin Wang, Wen-Pei Gan, En-Wei Liang, Yu-Jing Qin, Zheyu Lin, Lin-Na Xu, Min-Xuan Cai, Ji-An Jiang, Xu Kong, Jiaxun Li, Long Li, Jian-Guo Wang, Ze-Lin Xu, Yongquan Xue, Ye-Fei Yuan, Jingquan Cheng, Lulu Fan, Jie Gao, Lei Hu, Weida Hu , et al. (20 additional authors not shown)

Abstract: High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The… ▽ More High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The bump represents the longest separation time from the main peak among known TDEs to date. The main UV/optical outburst declines as $t^{-4.10}$, making it one of the fastest decaying optically selected TDEs. Furthermore, we detected sporadic X-ray emission 30 days after the UV/optical peak, accompanied by a reduction in the period of inactivity. It is proposed that the UV/optical bump could be caused by the self-intersection of the stream debris, whereas the primary peak is generated by the reprocessed emission of the accretion process. In addition, our results suggest that episodic X-ray radiation during the initial phase of decline may be due to the patched obscurer surrounding the accretion disk, a phenomenon associated with the inhomogeneous reprocessing process. The double TDE scenario, in which two stars are disrupted in sequence, is also a possible explanation for producing the observed early bump and main peak. We anticipate that the multicolor light curves of TDEs, especially in the very early stages, and the underlying physics can be better understood in the near future with the assistance of dedicated surveys such as the deep high-cadence survey of the 2.5-meter Wide Field Survey Telescope (WFST). △ Less

Submitted 26 March, 2024; v1 submitted 3 March, 2024; originally announced March 2024.

Comments: 14 pages, 8 figures,accepted for publication by ApJL

Hot Jupiter Formation in Dense Star Clusters

Authors: Leonard Benkendorff, Francesco Flammini Dotti, Katja Stock, Maxwell Xu Cai, Rainer Spurzem

Abstract: Hot Jupiters (HJ) are defined as Jupiter-mass exoplanets orbiting around their host star with an orbital period < 10 days. It is assumed that HJ do not form in-situ but ex-situ. Recent discoveries show that star clusters contribute to the formation of HJ. We present direct $N$-body simulations of planetary systems in star clusters and analyze the formation of HJ in them. We combine two direct $N$-… ▽ More Hot Jupiters (HJ) are defined as Jupiter-mass exoplanets orbiting around their host star with an orbital period < 10 days. It is assumed that HJ do not form in-situ but ex-situ. Recent discoveries show that star clusters contribute to the formation of HJ. We present direct $N$-body simulations of planetary systems in star clusters and analyze the formation of HJ in them. We combine two direct $N$-body codes: NBODY6++GPU for the dynamics of dense star clusters with 32 000 and 64 000 stellar members and LonelyPlanets used to follow 200 identical planetary systems around solar mass stars in those star clusters. We use different sets with 3, 4, or 5 planets and with the innermost planet at a semi-major axis of 5 au or 1 au and follow them for 100 Myr in our simulations. The results indicate that HJs are generated with high efficiency in dense star clusters if the innermost planet is already close to the host star at a semi-major axis of 1 au. If the innermost planet is initially beyond a semi-major axis of 5 au, the probability of a potential HJ ranges between $1.5-4.5$ percent. Very dense stellar neighborhoods tend to eject planets rather than forming HJs. A correlation between HJ formation and angular momentum deficit (AMD) is not witnessed. Young Hot Jupiters ($t_{\rm age} < 100$ Myrs) have only been found, in our simulations, in planetary systems with the innermost planet at a semi-major axis of 1 au. △ Less

Submitted 21 January, 2024; originally announced January 2024.

A hybrid approach for solving the gravitational N-body problem with Artificial Neural Networks

Authors: Veronica Saz Ulibarrena, Philipp Horn, Simon Portegies Zwart, Elena Sellentin, Barry Koren, Maxwell X. Cai

Abstract: Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural Networks (ANNs) to replace expensive parts of the integration of planetary systems. Neural networks that include physical knowledge have grown in popularity in the l… ▽ More Simulating the evolution of the gravitational N-body problem becomes extremely computationally expensive as N increases since the problem complexity scales quadratically with the number of bodies. We study the use of Artificial Neural Networks (ANNs) to replace expensive parts of the integration of planetary systems. Neural networks that include physical knowledge have grown in popularity in the last few years, although few attempts have been made to use them to speed up the simulation of the motion of celestial bodies. We study the advantages and limitations of using Hamiltonian Neural Networks to replace computationally expensive parts of the numerical simulation. We compare the results of the numerical integration of a planetary system with asteroids with those obtained by a Hamiltonian Neural Network and a conventional Deep Neural Network, with special attention to understanding the challenges of this problem. Due to the non-linear nature of the gravitational equations of motion, errors in the integration propagate. To increase the robustness of a method that uses neural networks, we propose a hybrid integrator that evaluates the prediction of the network and replaces it with the numerical solution if considered inaccurate. Hamiltonian Neural Networks can make predictions that resemble the behavior of symplectic integrators but are challenging to train and in our case fail when the inputs differ ~7 orders of magnitude. In contrast, Deep Neural Networks are easy to train but fail to conserve energy, leading to fast divergence from the reference solution. The hybrid integrator designed to include the neural networks increases the reliability of the method and prevents large energy errors without increasing the computing cost significantly. For this problem, the use of neural networks results in faster simulations when the number of asteroids is >70. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: Accepted for publication in the Journal of Computational Physics

The Lobster Eye Imager for Astronomy Onboard the SATech-01 Satellite

Authors: Z. X. Ling, X. J. Sun, C. Zhang, S. L. Sun, G. Jin, S. N. Zhang, X. F. Zhang, J. B. Chang, F. S. Chen, Y. F. Chen, Z. W. Cheng, W. Fu, Y. X. Han, H. Li, J. F. Li, Y. Li, Z. D. Li, P. R. Liu, Y. H. Lv, X. H. Ma, Y. J. Tang, C. B. Wang, R. J. Xie, Y. L. Xue, A. L. Yan , et al. (101 additional authors not shown)

Abstract: The Lobster Eye Imager for Astronomy (LEIA), a pathfinder of the Wide-field X-ray Telescope of the Einstein Probe (EP) mission, was successfully launched onboard the SATech-01 satellite of the Chinese Academy of Sciences on 27 July 2022. In this paper, we introduce the design and on-ground test results of the LEIA instrument. Using state-of-the-art Micro-Pore Optics (MPO), a wide field-of-view (Fo… ▽ More The Lobster Eye Imager for Astronomy (LEIA), a pathfinder of the Wide-field X-ray Telescope of the Einstein Probe (EP) mission, was successfully launched onboard the SATech-01 satellite of the Chinese Academy of Sciences on 27 July 2022. In this paper, we introduce the design and on-ground test results of the LEIA instrument. Using state-of-the-art Micro-Pore Optics (MPO), a wide field-of-view (FoV) of 346 square degrees (18.6 degrees * 18.6 degrees) of the X-ray imager is realized. An optical assembly composed of 36 MPO chips is used to focus incident X-ray photons, and four large-format complementary metal-oxide semiconductor (CMOS) sensors, each of 6 cm * 6 cm, are used as the focal plane detectors. The instrument has an angular resolution of 4 - 8 arcmin (in FWHM) for the central focal spot of the point spread function, and an effective area of 2 - 3 cm2 at 1 keV in essentially all the directions within the field of view. The detection passband is 0.5 - 4 keV in the soft X-rays and the sensitivity is 2 - 3 * 10-11 erg s-1 cm-2 (about 1 mini-Crab) at 1,000 second observation. The total weight of LEIA is 56 kg and the power is 85 W. The satellite, with a design lifetime of 2 years, operates in a Sun-synchronous orbit of 500 km with an orbital period of 95 minutes. LEIA is paving the way for future missions by verifying in flight the technologies of both novel focusing imaging optics and CMOS sensors for X-ray observation, and by optimizing the working setups of the instrumental parameters. In addition, LEIA is able to carry out scientific observations to find new transients and to monitor known sources in the soft X-ray band, albeit limited useful observing time available. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted by RAA

First wide field-of-view X-ray observations by a lobster eye focusing telescope in orbit

Authors: C. Zhang, Z. X. Ling, X. J. Sun, S. L. Sun, Y. Liu, Z. D. Li, Y. L. Xue, Y. F. Chen, Y. F. Dai, Z. Q. Jia, H. Y. Liu, X. F. Zhang, Y. H. Zhang, S. N. Zhang, F. S. Chen, Z. W. Cheng, W. Fu, Y. X. Han, H. Li, J. F. Li, Y. Li, P. R. Liu, X. H. Ma, Y. J. Tang, C. B. Wang , et al. (53 additional authors not shown)

Abstract: As a novel X-ray focusing technology, lobster eye micro-pore optics (MPO) feature both a wide observing field of view and true imaging capability, promising sky monitoring with significantly improved sensitivity and spatial resolution in soft X-rays. Since first proposed by Angel (1979), the optics have been extensively studied, developed and trialed over the past decades. In this Letter, we repor… ▽ More As a novel X-ray focusing technology, lobster eye micro-pore optics (MPO) feature both a wide observing field of view and true imaging capability, promising sky monitoring with significantly improved sensitivity and spatial resolution in soft X-rays. Since first proposed by Angel (1979), the optics have been extensively studied, developed and trialed over the past decades. In this Letter, we report on the first-light results from a flight experiment of the Lobster Eye Imager for Astronomy ($LEIA$), a pathfinder of the wide-field X-ray telescope of the Einstein Probe mission. The piggyback imager, launched in July 2022, has a mostly un-vignetted field of view of $18.6^\circ \times 18.6^\circ $. Its spatial resolution is in the range of 4$-$7 arcmin in FWHM and the focal spot effective area is 2$-$3 cm$^2$, both showing only mild fluctuations across the field of view. We present images of the Galactic center region, Sco X-1 and the diffuse Cygnus Loop nebular taken in snapshot observations over 0.5$-$4 keV. These are truly wide-field X-ray images of celestial bodies observed, for the first time, by a focusing imaging telescope. Initial analyses of the in-flight data show excellent agreement between the observed images and the on-ground calibration and simulations. The instrument and its characterization are briefly described, as well as the flight experiment. The results provide a solid basis for the development of the present and proposed wide-field X-ray missions using lobster eye MPO. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 11 pages, 4 figures. Accepted for publication in Astrophysical Journal Letter

Search for relativistic fractionally charged particles in space

Authors: DAMPE Collaboration, F. Alemanno, C. Altomare, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, M. S. Cai, E. Casilli, E. Catanzani, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, M. Y. Cui, T. S. Cui, Y. X. Cui, H. T. Dai, A. De-Benedittis, I. De Mitri, F. de Palma, M. Deliyergiyev, A. Di Giovanni, M. Di Santo , et al. (126 additional authors not shown)

Abstract: More than a century after the performance of the oil drop experiment, the possible existence of fractionally charged particles FCP still remains unsettled. The search for FCPs is crucial for some extensions of the Standard Model in particle physics. Most of the previously conducted searches for FCPs in cosmic rays were based on experiments underground or at high altitudes. However, there have been… ▽ More More than a century after the performance of the oil drop experiment, the possible existence of fractionally charged particles FCP still remains unsettled. The search for FCPs is crucial for some extensions of the Standard Model in particle physics. Most of the previously conducted searches for FCPs in cosmic rays were based on experiments underground or at high altitudes. However, there have been few searches for FCPs in cosmic rays carried out in orbit other than AMS-01 flown by a space shuttle and BESS by a balloon at the top of the atmosphere. In this study, we conduct an FCP search in space based on on-orbit data obtained using the DArk Matter Particle Explorer (DAMPE) satellite over a period of five years. Unlike underground experiments, which require an FCP energy of the order of hundreds of GeV, our FCP search starts at only a few GeV. An upper limit of $6.2\times 10^{-10}~~\mathrm{cm^{-2}sr^{-1} s^{-1}}$ is obtained for the flux. Our results demonstrate that DAMPE exhibits higher sensitivity than experiments of similar types by three orders of magnitude that more stringently restricts the conditions for the existence of FCP in primary cosmic rays. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: 19 pages, 6 figures, accepted by PRD

Report number: 106, 063026

Journal ref: Physical Review D 106.6 (2022): 063026

Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

Authors: Katja Stock, Dimitri Veras, Maxwell X. Cai, Rainer Spurzem, Simon Portegies Zwart

Abstract: Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of th… ▽ More Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution simulations. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Comments: 16 pages, 18 figures, accepted for publication in MNRAS

Search for gamma-ray spectral lines with the DArk Matter Particle Explorer

Authors: Francesca Alemanno, Qi An, Philipp Azzarello, Felicia Carla Tiziana Barbato, Paolo Bernardini, Xiao-Jun Bi, Ming-Sheng Cai, Elisabetta Casilli, Enrico Catanzani, Jin Chang, Deng-Yi Chen, Jun-Ling Chen, Zhan-Fang Chen, Ming-Yang Cui, Tian-Shu Cui, Yu-Xing Cui, Hao-Ting Dai, Antonio De Benedittis, Ivan De Mitri, Francesco de Palma, Maksym Deliyergiyev, Margherita Di Santo, Qi Ding, Tie-Kuang Dong, Zhen-Xing Dong , et al. (121 additional authors not shown)

Abstract: The DArk Matter Particle Explorer (DAMPE) is well suitable for searching for monochromatic and sharp $γ$-ray structures in the GeV$-$TeV range thanks to its unprecedented high energy resolution. In this work, we search for $γ$-ray line structures using five years of DAMPE data. To improve the sensitivity, we develop two types of dedicated data sets (including the BgoOnly data which is the first ti… ▽ More The DArk Matter Particle Explorer (DAMPE) is well suitable for searching for monochromatic and sharp $γ$-ray structures in the GeV$-$TeV range thanks to its unprecedented high energy resolution. In this work, we search for $γ$-ray line structures using five years of DAMPE data. To improve the sensitivity, we develop two types of dedicated data sets (including the BgoOnly data which is the first time to be used in the data analysis for the calorimeter-based gamma-ray observatories) and adopt the signal-to-noise ratio optimized regions of interest (ROIs) for different DM density profiles. No line signals or candidates are found between 10 and 300 GeV in the Galaxy. The constraints on the velocity-averaged cross section for $χχ\to γγ$ and the decay lifetime for $χ\to γν$, both at 95% confidence level, have been calculated and the systematic uncertainties have been taken into account. Comparing to the previous Fermi-LAT results, though DAMPE has an acceptance smaller by a factor of $\sim 10$, similar constraints on the DM parameters are achieved and below 100 GeV the lower limits on the decay lifetime are even stronger by a factor of a few. Our results demonstrate the potential of high-energy-resolution observations on dark matter detection. △ Less

Submitted 6 December, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: 14 pages, 8 figures. Update the content to keep up with the published version

Journal ref: Science Bulletin, Volume 67, Issue 7, 15 April 2022, Pages 679-684

Neural Symplectic Integrator with Hamiltonian Inductive Bias for the Gravitational $N$-body Problem

Authors: Maxwell X. Cai, Simon Portegies Zwart, Damian Podareanu

Abstract: The gravitational $N$-body problem, which is fundamentally important in astrophysics to predict the motion of $N$ celestial bodies under the mutual gravity of each other, is usually solved numerically because there is no known general analytical solution for $N>2$. Can an $N$-body problem be solved accurately by a neural network (NN)? Can a NN observe long-term conservation of energy and orbital a… ▽ More The gravitational $N$-body problem, which is fundamentally important in astrophysics to predict the motion of $N$ celestial bodies under the mutual gravity of each other, is usually solved numerically because there is no known general analytical solution for $N>2$. Can an $N$-body problem be solved accurately by a neural network (NN)? Can a NN observe long-term conservation of energy and orbital angular momentum? Inspired by Wistom & Holman (1991)'s symplectic map, we present a neural $N$-body integrator for splitting the Hamiltonian into a two-body part, solvable analytically, and an interaction part that we approximate with a NN. Our neural symplectic $N$-body code integrates a general three-body system for $10^{5}$ steps without diverting from the ground truth dynamics obtained from a traditional $N$-body integrator. Moreover, it exhibits good inductive bias by successfully predicting the evolution of $N$-body systems that are no part of the training set. △ Less

Submitted 28 November, 2021; originally announced November 2021.

Comments: 7 pages, 2 figures, accepted for publication at the NeurIPS 2021 workshop "Machine Learning and the Physical Sciences"

Fast and Credible Likelihood-Free Cosmology with Truncated Marginal Neural Ratio Estimation

Authors: Alex Cole, Benjamin Kurt Miller, Samuel J. Witte, Maxwell X. Cai, Meiert W. Grootes, Francesco Nattino, Christoph Weniger

Abstract: Sampling-based inference techniques are central to modern cosmological data analysis; these methods, however, scale poorly with dimensionality and typically require approximate or intractable likelihoods. In this paper we describe how Truncated Marginal Neural Ratio Estimation (TMNRE) (a new approach in so-called simulation-based inference) naturally evades these issues, improving the $(i)$ effici… ▽ More Sampling-based inference techniques are central to modern cosmological data analysis; these methods, however, scale poorly with dimensionality and typically require approximate or intractable likelihoods. In this paper we describe how Truncated Marginal Neural Ratio Estimation (TMNRE) (a new approach in so-called simulation-based inference) naturally evades these issues, improving the $(i)$ efficiency, $(ii)$ scalability, and $(iii)$ trustworthiness of the inferred posteriors. Using measurements of the Cosmic Microwave Background (CMB), we show that TMNRE can achieve converged posteriors using orders of magnitude fewer simulator calls than conventional Markov Chain Monte Carlo (MCMC) methods. Remarkably, the required number of samples is effectively independent of the number of nuisance parameters. In addition, a property called \emph{local amortization} allows the performance of rigorous statistical consistency checks that are not accessible to sampling-based methods. TMNRE promises to become a powerful tool for cosmological data analysis, particularly in the context of extended cosmologies, where the timescale required for conventional sampling-based inference methods to converge can greatly exceed that of simple cosmological models such as $Λ$CDM. To perform these computations, we use an implementation of TMNRE via the open-source code \texttt{swyft}. △ Less

Submitted 8 November, 2022; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: v2: accepted journal version. v1: 37 pages, 13 figures. \texttt{swyft} is available at https://github.com/undark-lab/swyft, and demonstration code for cosmological examples is available at https://github.com/acole1221/swyft-CMB

Observations of Forbush Decreases of cosmic ray electrons and positrons with the Dark Matter Particle Explorer

Authors: Francesca Alemanno, Qi An, Philipp Azzarello, Felicia Carla Tiziana Barbato, Paolo Bernardini, XiaoJun Bi, MingSheng Cai, Elisabetta Casilli, Enrico Catanzani, Jin Chang, DengYi Chen, JunLing Chen, ZhanFang Chen, MingYang Cui, TianShu Cui, YuXing Cui, HaoTing Dai, Antonio De Benedittis, Ivan De Mitri, Francesco de Palma, Maksym Deliyergiyev, Margherita Di Santo, Qi Ding, TieKuang Dong, ZhenXing Dong , et al. (124 additional authors not shown)

Abstract: The Forbush Decrease (FD) represents the rapid decrease of the intensities of charged particles accompanied with the coronal mass ejections (CMEs) or high-speed streams from coronal holes. It has been mainly explored with ground-based neutron monitors network which indirectly measure the integrated intensities of all species of cosmic rays by counting secondary neutrons produced from interaction b… ▽ More The Forbush Decrease (FD) represents the rapid decrease of the intensities of charged particles accompanied with the coronal mass ejections (CMEs) or high-speed streams from coronal holes. It has been mainly explored with ground-based neutron monitors network which indirectly measure the integrated intensities of all species of cosmic rays by counting secondary neutrons produced from interaction between atmosphere atoms and cosmic rays. The space-based experiments can resolve the species of particles but the energy ranges are limited by the relative small acceptances except for the most abundant particles like protons and helium. Therefore, the FD of cosmic ray electrons and positrons have just been investigated by the PAMELA experiment in the low energy range ($<5$ GeV) with limited statistics. In this paper, we study the FD event occurred in September, 2017, with the electron and positron data recorded by the Dark Matter Particle Explorer. The evolution of the FDs from 2 GeV to 20 GeV with a time resolution of 6 hours are given. We observe two solar energetic particle events in the time profile of the intensity of cosmic rays, the earlier and weak one has not been shown in the neutron monitor data. Furthermore, both the amplitude and recovery time of fluxes of electrons and positrons show clear energy-dependence, which is important in probing the disturbances of the interplanetary environment by the coronal mass ejections. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: This article is dedicated to the 72nd anniversary of People's Republic of China

Oort cloud Ecology II: The chronology of the formation of the Oort cloud

Authors: Simon Portegies Zwart, Santiago Torres, Maxwell X. Cai, Anthony Brown

Abstract: We present a chronology of the formation and early evolution of the Oort cloud by simulations. These simulations start with the Solar System being born with planets and asteroids in a stellar cluster orbiting the Galactic center. Upon ejection from its birth environment, we continue to follow the evolution of the Solar System while it navigates the Galaxy as an isolated planetary system. We conclu… ▽ More We present a chronology of the formation and early evolution of the Oort cloud by simulations. These simulations start with the Solar System being born with planets and asteroids in a stellar cluster orbiting the Galactic center. Upon ejection from its birth environment, we continue to follow the evolution of the Solar System while it navigates the Galaxy as an isolated planetary system. We conclude that the range in semi-major axis between 100au and several 10$^3$\,au still bears the signatures of the Sun being born in a 1000MSun/pc$^3$ star cluster, and that most of the outer Oort cloud formed after the Solar System was ejected. The ejection of the Solar System, we argue, happened between 20Myr and 50Myr after its birth. Trailing and leading trails of asteroids and comets along the Sun's orbit in the Galactic potential are the by-product of the formation of the Oort cloud. These arms are composed of material that became unbound from the Solar System when the Oort cloud formed. Today, the bulk of the material in the Oort cloud ($\sim 70$\%) originates from the region in the circumstellar disk that was located between $\sim 15$\,au and $\sim 35$\,au, near the current location of the ice giants and the Centaur family of asteroids. According to our simulations, this population is eradicated if the ice-giant planets are born in orbital resonance. Planet migration or chaotic orbital reorganization occurring while the Solar System is still a cluster member is, according to our model, inconsistent with the presence of the Oort cloud. About half the inner Oort cloud, between 100 and $10^4$\,au, and a quarter of the material in the outer Oort cloud, $\apgt 10^4$\,au, could be non-native to the Solar System but was captured from free-floating debris in the cluster or from the circumstellar disk of other stars in the birth cluster. △ Less

Submitted 11 June, 2021; v1 submitted 26 May, 2021; originally announced May 2021.

Comments: Corrected proofs for publication in A&A

Journal ref: A&A 652, A144 (2021)

Measurement of the cosmic ray helium energy spectrum from 70 GeV to 80 TeV with the DAMPE space mission

Authors: F. Alemanno, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, M. S. Cai, E. Catanzani, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, M. Y. Cui, T. S. Cui, Y. X. Cui, H. T. Dai, A. D'Amone, A. De Benedittis, I. De Mitri, F. de Palma, M. Deliyergiyev, M. Di Santo, T. K. Dong, Z. X. Dong, G. Donvito , et al. (120 additional authors not shown)

Abstract: The measurement of the energy spectrum of cosmic ray helium nuclei from 70 GeV to 80 TeV using 4.5 years of data recorded by the DArk Matter Particle Explorer (DAMPE) is reported in this work. A hardening of the spectrum is observed at an energy of about 1.3 TeV, similar to previous observations. In addition, a spectral softening at about 34 TeV is revealed for the first time with large statistics… ▽ More The measurement of the energy spectrum of cosmic ray helium nuclei from 70 GeV to 80 TeV using 4.5 years of data recorded by the DArk Matter Particle Explorer (DAMPE) is reported in this work. A hardening of the spectrum is observed at an energy of about 1.3 TeV, similar to previous observations. In addition, a spectral softening at about 34 TeV is revealed for the first time with large statistics and well controlled systematic uncertainties, with an overall significance of $4.3σ$. The DAMPE spectral measurements of both cosmic protons and helium nuclei suggest a particle charge dependent softening energy, although with current uncertainties a dependence on the number of nucleons cannot be ruled out. △ Less

Submitted 21 May, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: 11 pages, 13 figures, published in Phys. Rev. Lett. Add one more digit for first three columns in Table S2

Journal ref: Phys. Rev. Lett. 126, 201102 (2021)

Inside-Out Planet Formation: VI. Oligarchic Coagulation of Planetesimals from a Pebble Ring?

Authors: Maxwell X. Cai, Jonathan C. Tan, Simon Portegies Zwart

Abstract: Inside-Out Planet Formation (IOPF) is a theory addressing the origin of Systems of Tightly-Packed Inner Planets (STIPs) via {\it in situ} formation and growth of the planets. It predicts that a pebble ring is established at the pressure maximum associated with the dead zone inner boundary (DZIB) with an inner disk magnetorotational instability (MRI)-active region. Using direct $N$-body simulations… ▽ More Inside-Out Planet Formation (IOPF) is a theory addressing the origin of Systems of Tightly-Packed Inner Planets (STIPs) via {\it in situ} formation and growth of the planets. It predicts that a pebble ring is established at the pressure maximum associated with the dead zone inner boundary (DZIB) with an inner disk magnetorotational instability (MRI)-active region. Using direct $N$-body simulations, we study the collisional evolution of planetesimals formed from such a pebble ring, in particular examining whether a single dominant planet emerges. We consider a variety of models, including some in which the planetesimals are continuing to grow via pebble accretion. We find that the planetesimal ring undergoes oligarchic evolution, and typically turns into 2 or 3 surviving oligarchs on nearly coplanar and circular orbits, independent of the explored initial conditions or form of pebble accretion. The most massive oligarchs typically consist of about $70\%$ of the total mass, with the building-up process typically finishing within $\sim 10^5$ years. However, a relatively massive secondary planet always remains with $\sim30-65\%$ of the mass of the primary. Such secondary planets have properties that are inconsistent with the observed properties of the innermost pairs of planets in STIPs. Thus, for IOPF to be a viable theory for STIP formation, it needs to be shown how oligarchic growth of a relatively massive secondary from the initial pebble ring can be avoided. We discuss some potential additional physical processes that should be included in the modeling and explored as next steps. △ Less

Submitted 13 December, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 15 pages, 10 figures, accepted for publication in MNRAS

Energy Partition in Four Confined Circular-Ribbon Flares

Authors: Z. M. Cai, Q. M. Zhang, Z. J. Ning, Y. N. Su, H. S. Ji

Abstract: In this study, we investigated the energy partition of four confined circular-ribbon flares (CRFs) near the solar disk center, which are observed simultaneously by SDO, GOES, and RHESSI. We calculated different energy components, including the radiative outputs in 1$-$8, 1$-$70, and 70$-$370 Å, total radiative loss, peak thermal energy derived from GOES and RHESSI, nonthermal energy in flare-accel… ▽ More In this study, we investigated the energy partition of four confined circular-ribbon flares (CRFs) near the solar disk center, which are observed simultaneously by SDO, GOES, and RHESSI. We calculated different energy components, including the radiative outputs in 1$-$8, 1$-$70, and 70$-$370 Å, total radiative loss, peak thermal energy derived from GOES and RHESSI, nonthermal energy in flare-accelerated electrons, and magnetic free energy before flares. It is found that the energy components increase systematically with the flare class, indicating that more energies are involved in larger flares. The magnetic free energies are larger than the nonthermal energies and radiative outputs of flares, which is consistent with the magnetic nature of flares. The ratio $\frac{E_{nth}}{E_{mag}}$ of the four flares, being 0.70$-$0.76, is considerably higher than that of eruptive flares. Hence, this ratio may serve as an important factor that discriminates confined and eruptive flares. The nonthermal energies are sufficient to provide the heating requirements including the peak thermal energy and radiative loss. Our findings impose constraint on theoretical models of confined CRFs and have potential implication for the space weather forecast. △ Less

Submitted 19 February, 2021; originally announced February 2021.

Comments: 20 pages, 14 figures, accepted for publication in Solar Physics, comments are welcome

Mapping Diffuse Emission in Lyman UV band

Authors: Li Ji, Zheng Lou, Jinlong Zhang, Keqiang Qiu, Shuangying Li, Wei Sun, Shuping Yan, Shuinai Zhang, Yuan Qian, Sen Wang, Klaus Werner, Taotao Fang, Tinggui Wang, Jürgen Barnstedt, Sebastian Buntrock, Mingsheng Cai, Wen Chen, Lauro Conti, Lei Deng, Sebastian Diebold, Shaojun Fu, Jianhua Guo, Lars Hanke, Yilin Hong, Christoph Kalkuhl , et al. (16 additional authors not shown)

Abstract: The CAFE (Census of warm-hot intergalactic medium, Accretion, and Feedback Explorer) and LyRIC (Lyman UV Radiation from Interstellar medium and Circum-galactic medium) have been proposed to the space agencies in China respectively. CAFE was first proposed in 2015 as a joint scientific CAS-ESA small space mission. LyRIC was proposed in 2019 as an independent external payload operating on the Chines… ▽ More The CAFE (Census of warm-hot intergalactic medium, Accretion, and Feedback Explorer) and LyRIC (Lyman UV Radiation from Interstellar medium and Circum-galactic medium) have been proposed to the space agencies in China respectively. CAFE was first proposed in 2015 as a joint scientific CAS-ESA small space mission. LyRIC was proposed in 2019 as an independent external payload operating on the Chinese Space Station. Both missions are dedicated to mapping the Lyman UV emission (ionized oxygen (O VI) resonance lines at 103.2 and 103.8 nm, and Lyman series) for the diffuse sources in our Galaxy and the circum-galactic mediums of the nearby galaxies. We present the primary science objectives, mission concepts, the enabling technologies, as well as the current status. △ Less

Submitted 15 December, 2020; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 9 pages, 6 figures, SPIE Astronomical Telescopes + Instrumentation 2020 conference Paper No. 11444-4

DeepGalaxy: Deducing the Properties of Galaxy Mergers from Images Using Deep Neural Networks

Authors: Maxwell X. Cai, Jeroen Bédorf, Vikram A. Saletore, Valeriu Codreanu, Damian Podareanu, Adel Chaibi, Penny X. Qian

Abstract: Galaxy mergers, the dynamical process during which two galaxies collide, are among the most spectacular phenomena in the Universe. During this process, the two colliding galaxies are tidally disrupted, producing significant visual features that evolve as a function of time. These visual features contain valuable clues for deducing the physical properties of the galaxy mergers. In this work, we pro… ▽ More Galaxy mergers, the dynamical process during which two galaxies collide, are among the most spectacular phenomena in the Universe. During this process, the two colliding galaxies are tidally disrupted, producing significant visual features that evolve as a function of time. These visual features contain valuable clues for deducing the physical properties of the galaxy mergers. In this work, we propose DeepGalaxy, a visual analysis framework trained to predict the physical properties of galaxy mergers based on their morphology. Based on an encoder-decoder architecture, DeepGalaxy encodes the input images to a compressed latent space $z$, and determines the similarity of images according to the latent-space distance. DeepGalaxy consists of a fully convolutional autoencoder (FCAE) which generates activation maps at its 3D latent-space, and a variational autoencoder (VAE) which compresses the activation maps into a 1D vector, and a classifier that generates labels from the activation maps. The backbone of the FCAE can be fully customized according to the complexity of the images. DeepGalaxy demonstrates excellent scaling performance on parallel machines. On the Endeavour supercomputer, the scaling efficiency exceeds 0.93 when trained on 128 workers, and it maintains above 0.73 when trained with 512 workers. Without having to carry out expensive numerical simulations, DeepGalaxy makes inferences of the physical properties of galaxy mergers directly from images, and thereby achieves a speedup factor of $\sim 10^5$. △ Less

Submitted 22 October, 2020; originally announced October 2020.

Comments: 7 pages, 7 figures. Accepted for publication at the 2020 IEEE/ACM Fifth Workshop on Deep Learning on Supercomputers (DLS)

On the survival of resonant and non-resonant planetary systems in star clusters

Authors: Katja Stock, Maxwell X. Cai, Rainer Spurzem, M. B. N. Kouwenhoven, Simon Portegies Zwart

Abstract: Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing… ▽ More Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing $N = 8k, 16k, 32k$, and $64k$ stars. The planetary systems, which consist of the four Solar system giant planets Jupiter, Saturn, Uranus, and Neptune, are initialized in different orbital configurations, to study the effect of the system architecture on the dynamical evolution of the entire planetary system, and on the escape rate of the individual planets. We find that the current orbital parameters of the Solar system giants (with initially circular orbits, as well as with present-day eccentricities) and a slightly more compact configuration, have a high resilience against stellar perturbations. A configuration with initial mean-motion resonances of 3:2, 3:2, and 5:4 between the planets, which is inspired by the Nice model, and for which the two outermost planets are usually ejected within the first $10^5$ yr, is in many cases stabilized due to the removal of the resonances by external stellar perturbation and by the rapid ejection of at least one planet. Assigning all planets the same mass of 1 Jovian mass almost equalizes the survival fractions. Our simulations reproduce the broad diversity amongst observed exoplanet systems. We find not only many very wide and/or eccentric orbits, but also a significant number of (stable) retrograde orbits. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 21 pages, 21 figures, accepted by MNRAS

Linking the formation and fate of exo-Kuiper belts within solar system analogues

Authors: Dimitri Veras, Katja Reichert, Francesco Flammini Dotti, Maxwell X. Cai, Alexander J. Mustill, Andrew Shannon, Catriona H. McDonald, Simon Portegies Zwart, M. B. N. Kouwenhoven, Rainer Spurzem

Abstract: Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially em… ▽ More Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS

Measurement of the cosmic-ray proton spectrum from 40 GeV to 100 TeV with the DAMPE satellite

Authors: Q. An, R. Asfandiyarov, P. Azzarello, P. Bernardini, X. J. Bi, M. S. Cai, J. Chang, D. Y. Chen, H. F. Chen, J. L. Chen, W. Chen, M. Y. Cui, T. S. Cui, H. T. Dai, A. D'Amone, A. De Benedittis, I. De Mitri, M. Di Santo, M. Ding, T. K. Dong, Y. F. Dong, Z. X. Dong, G. Donvito, D. Droz, J. L. Duan , et al. (129 additional authors not shown)

Abstract: The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This i… ▽ More The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays. △ Less

Submitted 30 September, 2019; v1 submitted 27 September, 2019; originally announced September 2019.

Comments: 37 pages, 5 figures, published in Science Advances

Journal ref: Science Advances, Vol. 5, no. 9, eaax3793 (2019)

Planetary systems in a star cluster I: the Solar system scenario

Authors: Francesco Flammini Dotti, M. B. N. Kouwenhoven, Maxwell Xu Cai, Rainer Spurzem

Abstract: Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planeta… ▽ More Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planetary systems escape from the star cluster within 50 Myr; those with low escape speeds often remain intact during and after the escape process. While most planetary systems inside the star cluster remain intact, a subset is strongly perturbed during the first 50 Myr. Over the course of time, 0.3 % - 5.3 % of the planets escape, sometimes up to tens of millions of years after a stellar encounter occurred. Survival rates are highest for Jupiter, while Uranus and Neptune have the highest escape rates. Unless directly affected by a stellar encounter itself, Jupiter frequently serves as a barrier that protects the terrestrial planets from perturbations in the outer planetary system. In low-density environments, Jupiter provides protection from perturbations in the outer planetary system, while in high-density environments, direct perturbations of Jupiter by neighbouring stars is disruptive to habitable-zone planets. The diversity amongst planetary systems that is present in the star clusters at 50 Myr, and amongst the escaping planetary systems, is high, which contributes to explaining the high diversity of observed exoplanet systems in star clusters and in the Galactic field △ Less

Submitted 21 August, 2019; originally announced August 2019.

Comments: 20 pages, 46 figures, accepted for pubblication on Monthly Notices of the Royal Astronomical Society on date 13 August 2019

The on-orbit calibration of DArk Matter Particle Explorer

Authors: G. Ambrosi, Q. An, R. Asfandiyarov, P. Azzarello, P. Bernardini, M. S. Cai, M. Caragiulo, J. Chang, D. Y. Chen, H. F. Chen, J. L. Chen, W. Chen, M. Y. Cui, T. S. Cui, H. T. Dai, A. D'Amone, A. De Benedittis, I. De Mitri, M. Ding, M. Di Santo, J. N. Dong, T. K. Dong, Y. F. Dong, Z. X. Dong, D. Droz , et al. (133 additional authors not shown)

Abstract: The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO… ▽ More The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation. △ Less

Submitted 3 July, 2019; originally announced July 2019.

Journal ref: Astroparticle Physics, Volume 106, p. 18-34 (2019)

Galactic Tide and Local Stellar Perturbations on the Oort Cloud: Creation of Interstellar Comets

Authors: Santiago Torres, Maxwell X. Cai, A. G. A. Brown, S. Portegies Zwart

Abstract: Comets in the Oort cloud evolve under the influence of internal and external perturbations, such as giant planets, stellar passages, and the galactic tidal field. We aim to study the dynamical evolution of the comets in the Oort cloud, accounting for external perturbations (passing stars and the galactic tide). We first construct an analytical model of stellar encounters. We find that individual p… ▽ More Comets in the Oort cloud evolve under the influence of internal and external perturbations, such as giant planets, stellar passages, and the galactic tidal field. We aim to study the dynamical evolution of the comets in the Oort cloud, accounting for external perturbations (passing stars and the galactic tide). We first construct an analytical model of stellar encounters. We find that individual perturbations do not modify the dynamics of the comets in the cloud unless very close (< 0.5pc) encounters occur. Using proper motions, parallaxes, and radial velocities from Gaia DR2, we construct an astrometric catalogue of 14,659 stars that are within 50pc from the Sun. For all these stars we calculate the time and the closest distance to the Sun. We find that the cumulative effect of relatively distant ($\leq1$ pc) passing stars can perturb the comets in the Oort cloud. Finally, we study the dynamical evolution of the comets in the Oort cloud under the influence of multiple stellar encounters within 2.5pc from the Sun and the galactic tidal field over $\pm10$Myr. We considered two models for the Oort cloud, compact (a $\leq$0.25 pc) and extended (a$ \leq0.5$ pc). We find that the cumulative effect of stellar encounters is the major perturber of the Oort cloud for a compact configuration while for the extended, the galactic tide is the major perturber. In both cases, the effect of passing stars and the galactic tide raises the semi-major axis of $\sim1.1$\% of the comets at the edge of the cloud up to interstellar regions ($a >0.5$pc). This leads to the creation of transitional interstellar comets, which might become interstellar objects due to external perturbations. This raises the question about the existence of a cloud of objects in the interstellar space which might overlap with our Oort cloud if we consider that other planetary systems face similar processes for the ejection of comets. △ Less

Submitted 15 August, 2019; v1 submitted 25 June, 2019; originally announced June 2019.

Comments: 13 pages, 13 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 629, A139 (2019)

Planets In Young Massive Clusters: On the survivability of planets in young massive clusters and its implication of planet orbital architectures in globular clusters

Authors: Maxwell Xu Cai, S. Portegies Zwart, M. B. N. Kouwenhoven, Rainer Spurzem

Abstract: As of August 2019, among the more than 4000 confirmed exoplanets, only one has been detected in a globular cluster (GC) M4. The scarce of exoplanet detections motivates us to employ direct $N$-body simulations to investigate the dynamical stability of planets in young massive clusters (YMCs), which are potentially the progenitors of GCs. In an $N=128{\rm k}$ cluster of virial radius 1.7 pc (compar… ▽ More As of August 2019, among the more than 4000 confirmed exoplanets, only one has been detected in a globular cluster (GC) M4. The scarce of exoplanet detections motivates us to employ direct $N$-body simulations to investigate the dynamical stability of planets in young massive clusters (YMCs), which are potentially the progenitors of GCs. In an $N=128{\rm k}$ cluster of virial radius 1.7 pc (comparable to Westerlund-1), our simulations show that most wide-orbit planets ($a\geq 20$~au) will be ejected within a timescale of 10 Myr. Interestingly, more than $70\%$ of planets with $a<5$~au survive in the 100 Myr simulations. Ignoring planet-planet scattering and tidal damping, the survivability at $t$ Myr as a function of initial semi-major axis $a_0$ in au in such a YMC can be described as $f_{\rm surv}(a_0, t)=-0.33 \log_{10}(a_0) \left(1 - e^{-0.0482t} \right) + 1$. Upon ejection, about $28.8\%$ of free-floating planets (FFPs) have sufficient speeds to escape from the host cluster at a crossing timescale. The other FFPs will remain bound to the cluster potential, but the subsequent dynamical evolution of the stellar system can result in the delayed ejection of FFPs from the host cluster. Although a full investigation of planet population in GCs requires extending the simulations to multi-Gyr, our results suggest that wide-orbit planets and free-floating planets are unlikely to be found in GCs. △ Less

Submitted 9 September, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: 12 pages, 11 figures, 1 table, accepted by MNRAS

Tropical and Extratropical General Circulation with a Meridional Reversed Temperature Gradient as Expected in a High Obliquity Planet

Authors: Wanying Kang, Ming Cai, Eli Tziperman

Abstract: Planets with high obliquity receive more radiation in the polar regions than at low latitudes, and thus, assuming an ocean-covered surface with sufficiently high heat capacity, their meridional temperature gradient was shown to be reversed for the entire year. The objective of this work is to investigate the drastically different general circulation of such planets, with an emphasis on the tropica… ▽ More Planets with high obliquity receive more radiation in the polar regions than at low latitudes, and thus, assuming an ocean-covered surface with sufficiently high heat capacity, their meridional temperature gradient was shown to be reversed for the entire year. The objective of this work is to investigate the drastically different general circulation of such planets, with an emphasis on the tropical Hadley circulation and the mid-latitude baroclinic eddy structure. We use a 3D dry dynamic core model, accompanied by an eddy-free configuration and a generalized 2D Eady model. When the meridional temperature gradient $T_y$ is reversed, the Hadley cell remains in the same direction, because the surface wind pattern and hence the associated meridional Ekman transport are not changed, as required by the baroclinic eddy momentum transport. The Hadley cell under reversed $T_y$ also becomes much shallower and weaker, even when the magnitude of the gradient is the same as in the normal case. The shallowness is due to the bottom-heavy structure of the baroclinic eddies in the reverse case, and the weakness is due to the weak wave activity. We propose a new mechanism to explain the mid-latitude eddy structure for both cases, and verify it using the generalized Eady model. With seasonal variations included, the annual mean circulation resembles that under perpetual annual mean setup. Approaching the solstices, a strong cross-equator Hadley cell forms in both cases, and about 2/3 of the Hadley circulation is driven by eddies, as shown by eddy-free simulations and using a decomposition of the Hadley cell. △ Less

Submitted 18 February, 2019; originally announced February 2019.

Survivability of planetary systems in young and dense star clusters

Authors: Arjen van Elteren, Simon Portegies Zwart, Inti Pelupessy, Maxwell Cai, Steve McMillan

Abstract: We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced… ▽ More We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced. Four, five, or size planets per star were introduced in orbit around the 500 solar-like stars with a maximum orbital separation of 400au. Our study focuses on the production of free-floating planets. A total of 357 become unbound from a total of 2522 planets in the initial conditions of the simulation. Of these, 281 leave the cluster within the crossing timescale of the star cluster; the others remain bound to the cluster as free-floating intra-cluster planets. Five of these free-floating intra-cluster planets are captured at a later time by another star. The two main mechanisms by which planets are lost from their host star, ejection upon a strong encounter with another star or internal planetary scattering, drive the evaporation independent of planet mass of orbital separation at birth. The effect of small perturbations due to slow changes in the cluster potential are important for the evolution of planetary systems. In addition, the probability of a star to lose a planet is independent of the planet mass and independent of its initial orbital separation. As a consequence, the mass distribution of free-floating planets indistinguishable from the mass distribution of planets bound to their host star. △ Less

Submitted 11 March, 2019; v1 submitted 12 February, 2019; originally announced February 2019.

Comments: A&A in press

Journal ref: A&A 624, A120 (2019)

Planetary Systems in Star Clusters: the dynamical evolution and survival

Authors: Francesco Flammini Dotti, Maxwell Xu Cai, Rainer Spurzem, M. B. N. Kouwenhoven

Abstract: Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary co… ▽ More Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary companions. To understand the origin and dynamical evolution of such exoplanet systems, it is necessary to carefully study the effect of their environments. Here, we combine theoretical estimates with state-of-the-art numerical simulations of evolving planetary systems similar to our own solar system in different star cluster environments. We combine the REBOUND planetary system evolution code, and the NBODY6++GPU star cluster evolution code, integrated in the AMUSE multi-physics environment. With our study we can constrain the effect of external perturbations of different environments on the planets and debris structures of a wide variety of planetary systems, which may play a key role for the habitability of exoplanets in the Universe. △ Less

Submitted 30 November, 2018; originally announced November 2018.

Comments: 2 pages, 2 figures, Proceedings for Symposia IAU 345 "Origins: from the Protosun to the First Steps of Life" at the XXX general assembly in Vienna

A MODEST review

Authors: Anna Lisa Varri, Maxwell Xu Cai, Francisca Concha-Ramírez, František Dinnbier, Nora Lützgendorf, Václav Pavlík, Sara Rastello, Antonio Sollima, Long Wang, Alice Zocchi

Abstract: We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the… ▽ More We present an account of the state of the art in the fields explored by the research community invested in 'Modeling and Observing DEnse STellar systems'. For this purpose, we take as a basis the activities of the MODEST-17 conference, which was held at Charles University, Prague, in September 2017. Reviewed topics include recent advances in fundamental stellar dynamics, numerical methods for the solution of the gravitational N-body problem, formation and evolution of young and old star clusters and galactic nuclei, their elusive stellar populations, planetary systems, and exotic compact objects, with timely attention to black holes of different classes of mass and their role as sources of gravitational waves. Such a breadth of topics reflects the growing role played by collisional stellar dynamics in numerous areas of modern astrophysics. Indeed, in the next decade, many revolutionary instruments will enable the derivation of positions and velocities of individual stars in the Milky Way and its satellites and will detect signals from a range of astrophysical sources in different portions of the electromagnetic and gravitational spectrum, with an unprecedented sensitivity. On the one hand, this wealth of data will allow us to address a number of long-standing open questions in star cluster studies; on the other hand, many unexpected properties of these systems will come to light, stimulating further progress of our understanding of their formation and evolution. △ Less

Submitted 17 October, 2018; originally announced October 2018.

Comments: 42 pages; accepted for publication in 'Computational Astrophysics and Cosmology'. We are much grateful to the organisers of the MODEST-17 conference (Charles University, Prague, September 2017). We acknowledge the input provided by all MODEST-17 participants, and, more generally, by the members of the MODEST community

The orbital eccentricity of small planet systems

Authors: Vincent Van Eylen, Simon Albrecht, Xu Huang, Mariah G. MacDonald, Rebekah I. Dawson, Maxwell X. Cai, Daniel Foreman-Mackey, Mia S. Lundkvist, Victor Silva Aguirre, Ignas Snellen, Joshua N. Winn

Abstract: We determine the orbital eccentricities of individual small Kepler planets, through a combination of asteroseismology and transit light-curve analysis. We are able to constrain the eccentricities of 51 systems with a single transiting planet, which supplement our previous measurements of 66 planets in multi-planet systems. Through a Bayesian hierarchical analysis, we find evidence that systems wit… ▽ More We determine the orbital eccentricities of individual small Kepler planets, through a combination of asteroseismology and transit light-curve analysis. We are able to constrain the eccentricities of 51 systems with a single transiting planet, which supplement our previous measurements of 66 planets in multi-planet systems. Through a Bayesian hierarchical analysis, we find evidence that systems with only one detected transiting planet have a different eccentricity distribution than systems with multiple detected transiting planets. The eccentricity distribution of the single-transiting systems is well described by the positive half of a zero-mean Gaussian distribution with a dispersion $σ_e = 0.32 \pm 0.06$, while the multiple-transit systems are consistent with $σ_e = 0.083^{+0.015}_{-0.020}$. A mixture model suggests a fraction of $0.76^{+0.21}_{-0.12}$ of single-transiting systems have a moderate eccentricity, represented by a Rayleigh distribution that peaks at $0.26^{+0.04}_{-0.06}$. This finding may reflect differences in the formation pathways of systems with different numbers of transiting planets. We investigate the possibility that eccentricities are "self-excited" in closely packed planetary systems, as well as the influence of long-period giant companion planets. We find that both mechanisms can qualitatively explain the observations. We do not find any evidence for a correlation between eccentricity and stellar metallicity, as has been seen for giant planets. Neither do we find any evidence that orbital eccentricity is linked to the detection of a companion star. Along with this paper we make available all of the parameters and uncertainties in the eccentricity distributions, as well as the properties of individual systems, for use in future studies. △ Less

Submitted 18 November, 2018; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: Accepted for publication in The Astronomical Journal

Stability of exomoons around the Kepler transiting circumbinary planets

Authors: Adrian S. Hamers, Maxwell X. Cai, Javier Roa, Nathan Leigh

Abstract: The Kepler mission has detected a number of transiting circumbinary planets (CBPs). Although currently not detected, exomoons could be orbiting some of these CBPs, and they might be suitable for harboring life. A necessary condition for the existence of such exomoons is their long-term dynamical stability. Here, we investigate the stability of exomoons around the Kepler CBPs using numerical $N$-bo… ▽ More The Kepler mission has detected a number of transiting circumbinary planets (CBPs). Although currently not detected, exomoons could be orbiting some of these CBPs, and they might be suitable for harboring life. A necessary condition for the existence of such exomoons is their long-term dynamical stability. Here, we investigate the stability of exomoons around the Kepler CBPs using numerical $N$-body integrations. We determine regions of stability and obtain stability maps in the (a_m,i_pm) plane, where a_m is the initial exolunar semimajor axis with respect to the CBP, and i_pm is the initial inclination of the orbit of the exomoon around the planet with respect to the orbit of the planet around the stellar binary. Ignoring any dependence on i_pm, for most Kepler CBPs the stability regions are well described by the location of the 1:1 mean motion commensurability of the binary orbit with the orbit of the moon around the CBP. This is related to a destabilizing effect of the binary compared to the case if the binary were replaced by a single body, and which is borne out by corresponding 3-body integrations. For high inclinations, the evolution is dominated by Lidov-Kozai oscillations, which can bring moons in dynamically stable orbits to close proximity within the CBP, triggering strong interactions such as tidal evolution, tidal disruption, or direct collisions. This suggests that there is a dearth of highly-inclined exomoons around the Kepler CBPs, whereas coplanar exomoons are dynamically allowed. △ Less

Submitted 2 August, 2018; v1 submitted 15 June, 2018; originally announced June 2018.

Comments: Accepted for publication by MNRAS. 12 pages, 7 figures

Direct detection of a break in the teraelectronvolt cosmic-ray spectrum of electrons and positrons

Authors: G. Ambrosi, Q. An, R. Asfandiyarov, P. Azzarello, P. Bernardini, B. Bertucci, M. S. Cai, J. Chang, D. Y. Chen, H. F. Chen, J. L. Chen, W. Chen, M. Y. Cui, T. S. Cui, A. D'Amone, A. De Benedittis, I. De Mitri, M. Di Santo, J. N. Dong, T. K. Dong, Y. F. Dong, Z. X. Dong, G. Donvito, D. Droz, K. K. Duan , et al. (133 additional authors not shown)

Abstract: High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to… ▽ More High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $γ$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs. △ Less

Submitted 29 November, 2017; originally announced November 2017.

Comments: 18 pages, 6 figures, Nature in press, doi:10.1038/nature24475

Journal ref: Nature, 552, 63-66 (2017)

The origin of interstellar asteroidal objects like 1I/2017 U1 'Oumuamua

Authors: Simon Portegies Zwart, Santiago Torres, Inti Pelupessy, Jeroen Bedorf, Maxwell Cai

Abstract: We study the origin of the interstellar object 1I/2017 U1 'Oumuamua by juxtaposing estimates based on the observations with simulations. We speculate that objects like 'Oumuamua are formed in the debris disc as left over from the star and planet formation process, and subsequently liberated. The liberation process is mediated either by interaction with other stars in the parental star-cluster, by… ▽ More We study the origin of the interstellar object 1I/2017 U1 'Oumuamua by juxtaposing estimates based on the observations with simulations. We speculate that objects like 'Oumuamua are formed in the debris disc as left over from the star and planet formation process, and subsequently liberated. The liberation process is mediated either by interaction with other stars in the parental star-cluster, by resonant interactions within the planetesimal disc or by the relatively sudden mass loss when the host star becomes a compact object. Integrating backward in time in the Galactic potential together with stars from the Gaia-TGAS catalogue we find that about 1.3Myr ago 'Oumuamua passed the nearby star HIP 17288 within a mean distance of $1.3$pc. By comparing nearby observed L-dwarfs with simulations of the Galaxy we conclude that the kinematics of 'Oumuamua is consistent with relatively young objects of $1.1$--$1.7$Gyr. We just met 'Oumuamua by chance, and with a derived mean Galactic density of $\sim 3\times 10^{5}$ similarly sized objects within 100\,au from the Sun or $\sim 10^{14}$ per cubic parsec we expect about 2 to 12 such visitors per year within 1au from the Sun. △ Less

Submitted 11 May, 2018; v1 submitted 9 November, 2017; originally announced November 2017.

Comments: MNRAS (in press)

The signatures of the parental cluster on field planetary systems

Authors: Maxwell Xu Cai, Simon Portegies Zwart, Arjen van Elteren

Abstract: Due to the high stellar densities in young clusters, planetary systems formed in these environments are likely to have experienced perturbations from encounters with other stars. We carry out direct $N$-body simulations of multi-planet systems in star clusters to study the combined effects of stellar encounters and internal planetary dynamics. These planetary systems eventually become part of the… ▽ More Due to the high stellar densities in young clusters, planetary systems formed in these environments are likely to have experienced perturbations from encounters with other stars. We carry out direct $N$-body simulations of multi-planet systems in star clusters to study the combined effects of stellar encounters and internal planetary dynamics. These planetary systems eventually become part of the Galactic field population the parental cluster dissolves, which is where most presently-known exoplanets are observed. We show that perturbations induced by stellar encounters lead to distinct signatures in the field planetary systems, most prominently, the excited orbital inclinations and eccentricities. Planetary systems that form within the cluster's half-mass radius are more prone to such perturbations. The orbital elements are most strongly excited in the outermost orbit, but the effect propagates to the entire planetary system through secular evolution. Planet ejections may occur long after a stellar encounter. The surviving planets in these reduced systems tend to have, on average, higher inclinations and larger eccentricities compared to systems that were perturbed less strongly. As soon as the parental star cluster dissolves, external perturbations stop affecting the escaped planetary systems, and further evolution proceeds on a relaxation time scale. The outer regions of these ejected planetary systems tend to relax so slowly that their state carries the memory of their last strong encounter in the star cluster. Regardless of the stellar density, we observe a robust anticorrelation between multiplicity and mean inclination/eccentricity. We speculate that the "Kepler dichotomy" observed in field planetary systems is a natural consequence of their early evolution in the parental cluster. △ Less

Submitted 23 November, 2017; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: 9 pages, 6 figures, 2 tables, accepted for publication in MNRAS

SiMon: Simulation Monitor for Computational Astrophysics

Authors: Penny Xuran Qian, Maxwell Xu Cai, Simon Portegies Zwart, Ming Zhu

Abstract: Scientific discovery via numerical simulations is important in modern astrophysics. This relatively new branch of astrophysics has become possible due to the development of reliable numerical algorithms and the high performance of modern computing technologies. These enable the analysis of large collections of observational data and the acquisition of new data via simulations at unprecedented accu… ▽ More Scientific discovery via numerical simulations is important in modern astrophysics. This relatively new branch of astrophysics has become possible due to the development of reliable numerical algorithms and the high performance of modern computing technologies. These enable the analysis of large collections of observational data and the acquisition of new data via simulations at unprecedented accuracy and resolution. Ideally, simulations run until they reach some pre-determined termination condition, but often other factors cause extensive numerical approaches to break down at an earlier stage. In those cases, processes tend to be interrupted due to unexpected events in the software or the hardware. In those cases, the scientist handles the interrupt manually, which is time-consuming and prone to errors. We present the Simulation Monitor (SiMon) to automatize the farming of large and extensive simulation processes. Our method is light-weight, it fully automates the entire workflow management, operates concurrently across multiple platforms and can be installed in user space. Inspired by the process of crop farming, we perceive each simulation as a crop in the field and running simulation becomes analogous to growing crops. With the development of SiMon we relax the technical aspects of simulation management. The initial package was developed for extensive parameter searchers in numerical simulations, but it turns out to work equally well for automating the computational processing and reduction of observational data reduction. △ Less

Submitted 29 June, 2017; originally announced June 2017.

Comments: 13 pages, 6 figures, Accepted for publication in Publications of the Astronomical Society of the Pacific (PASP)

The DArk Matter Particle Explorer mission

Authors: J. Chang, G. Ambrosi, Q. An, R. Asfandiyarov, P. Azzarello, P. Bernardini, B. Bertucci, M. S. Cai, M. Caragiulo, D. Y. Chen, H. F. Chen, J. L. Chen, W. Chen, M. Y. Cui, T. S. Cui, A. D'Amone, A. De Benedittis, I. De Mitri, M. Di Santo, J. N. Dong, T. K. Dong, Y. F. Dong, Z. X. Dong, G. Donvito, D. Droz , et al. (139 additional authors not shown)

Abstract: The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives… ▽ More The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals. △ Less

Submitted 14 September, 2017; v1 submitted 26 June, 2017; originally announced June 2017.

Comments: 45 pages, including 29 figures and 6 tables. Published in Astropart. Phys

Journal ref: Astropart. Phys. 95 (2017) 6-24

Stability of Multiplanetary Systems in Star Clusters

Authors: Maxwell Xu Cai, M. B. N. Kouwenhoven, Simon F. Portegies Zwart, Rainer Spurzem

Abstract: Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems ar… ▽ More Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems are assigned to solar-type stars with $\sim 1 M_{\odot}$ and evolved for 50~Myr. We found that following the depletion of protoplanetary disks, external perturbations and planet-planet interactions are two driving mechanisms responsible for the destabilization of planetary systems. The planet survival rate varies from $\sim 95\%$ in the N=2k cluster to $\sim 60\%$ in the N=32k cluster, which suggests that most planetary systems can indeed survive in low-mass clusters, except in the central regions. We also find that planet ejections through stellar encounters are cumulative processes, as only $\sim 3\%$ of encounters are strong enough to excite the eccentricity by $Δe \geq 0.5$. Short-period planets can be perturbed through orbit crossings with long-period planets. When taking into account planet-planet interactions, the planet ejection rate nearly doubles, and therefore multiplicity contributes to the vulnerability of planetary systems. In each ensemble, $\sim 0.2\%$ of planetary orbits become retrograde due to random directions of stellar encounters. Our results predict that young low-mass star clusters are promising sites for next-generation planet surveys, yet low planet detection rates are expected in dense globular clusters such as 47 Tuc. Nevertheless, planets in denser stellar environments are likely to have shorter orbital periods, which enhances their detectability. △ Less

Submitted 12 June, 2017; originally announced June 2017.

Comments: 19 pages, 13 figures, 4 tables, accepted for publication in MNRAS

Planetary systems in star clusters

Authors: M. B. N. Kouwenhoven, Qi Shu, Maxwell Xu Cai, Rainer Spurzem

Abstract: Thousands of confirmed and candidate exoplanets have been identified in recent years. Consequently, theoretical research on the formation and dynamical evolution of planetary systems has seen a boost, and the processes of planet-planet scattering, secular evolution, and interaction between planets and gas/debris disks have been well-studied. Almost all of this work has focused on the formation and… ▽ More Thousands of confirmed and candidate exoplanets have been identified in recent years. Consequently, theoretical research on the formation and dynamical evolution of planetary systems has seen a boost, and the processes of planet-planet scattering, secular evolution, and interaction between planets and gas/debris disks have been well-studied. Almost all of this work has focused on the formation and evolution of isolated planetary systems, and neglect the effect of external influences, such as the gravitational interaction with neighbouring stars. Most stars, however, form in clustered environments that either quickly disperse, or evolve into open clusters. Under these conditions, young planetary systems experience frequent close encounters with other stars, at least during the first 1-10 Myr, which affects planets orbiting at any period range, as well as their debris structures. △ Less

Submitted 4 September, 2016; originally announced September 2016.

Comments: Proceedings of Cosmic-Lab: Star Clusters as Cosmic Laboratories for Astrophysics, Dynamics, and Fundamental Physics - MODEST16

Evolution of star clusters on eccentric orbits

Authors: Maxwell Xu Cai, Mark Gieles, Douglas C. Heggie, Anna Lisa Varri

Abstract: We study the evolution of star clusters on circular and eccentric orbits using direct $N$-body simulations. We model clusters with initially $N=8{\rm k}$ and $N=16{\rm k}$ single stars of the same mass, orbiting around a point-mass galaxy. For each orbital eccentricity that we consider, we find the apogalactic radius at which the cluster has the same lifetime as the cluster with the same $N$ on a… ▽ More We study the evolution of star clusters on circular and eccentric orbits using direct $N$-body simulations. We model clusters with initially $N=8{\rm k}$ and $N=16{\rm k}$ single stars of the same mass, orbiting around a point-mass galaxy. For each orbital eccentricity that we consider, we find the apogalactic radius at which the cluster has the same lifetime as the cluster with the same $N$ on a circular orbit. We show that then, the evolution of bound particle number and half-mass radius is approximately independent of eccentricity. Secondly, when we scale our results to orbits with the same semi-major axis, we find that the lifetimes are, to first order, independent of eccentricity. When the results of Baumgardt and Makino for a singular isothermal halo are scaled in the same way, the lifetime is again independent of eccentricity to first order, suggesting that this result is independent of the Galactic mass profile. From both sets of simulations we empirically derive the higher order dependence of the lifetime on eccentricity. Our results serve as benchmark for theoretical studies of the escape rate from clusters on eccentric orbits. Finally, our results can be useful for generative models for cold streams and cluster evolution models that are confined to spherical symmetry and/or time-independent tides, such as Fokker-Planck models, Monte Carlo models, and (fast) semi-analytic models. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: 7 pages, 8 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society

Block Time Step Storage Scheme for Astrophysical N-body Simulations

Authors: Maxwell Xu Cai, Yohai Meiron, M. B. N. Kouwenhoven, Paulina Assmann, Rainer Spurzem

Abstract: Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to take the computing power of massively parallel clusters with more than $10^5$ GPU cores. While unprecedented accuracy and dynamical scales have been ac… ▽ More Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to take the computing power of massively parallel clusters with more than $10^5$ GPU cores. While unprecedented accuracy and dynamical scales have been achieved, the enormous amount of data being generated continuously poses great challenges for the subsequent procedures of data analysis and archiving. As an urgent response to these challenges, in this paper we propose an adaptive storage scheme for simulation data, inspired by the block time step integration scheme found in a number of direct $N$-body integrators available nowadays. The proposed scheme, namely the block time step storage scheme, works by minimizing the data redundancy with assignments of data with individual output frequencies as required by the researcher. As demonstrated by benchmarks, the proposed scheme is applicable to a wide variety of simulations. Despite the main focus of developing a solution for direct $N$-body simulation data, the methodology is transferable for grid-based or tree-based simulations where hierarchical time stepping is used. △ Less

Submitted 24 June, 2015; originally announced June 2015.

Comments: 13 pages, 10 figures, 5 tables, Accepted for publication in The Astrophysical Journal Supplement Series

Journal ref: ApJS, 219, 31, 2015

Planetary Systems in Star Clusters

Authors: Maxwell Xu Cai, Rainer Spurzem, M. B. N. Kouwenhoven

Abstract: In the solar neighborhood, where the typical relaxation timescale is larger than the cosmic age, at least 10\% to 15\% of Sun-like stars have planetary systems with Jupiter-mass planets. In contrast, dense star clusters, charactered by frequent close encounters, have been found to host very few planets. We carry out numerical simulations with different initial conditions to investigate the dynamic… ▽ More In the solar neighborhood, where the typical relaxation timescale is larger than the cosmic age, at least 10\% to 15\% of Sun-like stars have planetary systems with Jupiter-mass planets. In contrast, dense star clusters, charactered by frequent close encounters, have been found to host very few planets. We carry out numerical simulations with different initial conditions to investigate the dynamical stability of planetary systems in star cluster environments. △ Less

Submitted 27 January, 2015; v1 submitted 7 January, 2015; originally announced January 2015.

Comments: Paper submitted for proceedings of IAU Symp. No. 312, 2 pages, 3 figures

Magnetic Interactions in Pre-Main-Sequence Binaries

Authors: Fred C. Adams, Michael J. Cai, Daniele Galli, Susana Lizano, Frank H. Shu

Abstract: Young stars typically have strong magnetic fields, so that the magnetospheres of newly formed close binaries can interact, dissipate energy, and produce synchrotron radiation. The V773 Tau A binary system, a pair of T Tauri stars with a 51 day orbit, displays such a signature, with peak emission taking place near periastron. This paper proposes that the observed emission arises from the change in… ▽ More Young stars typically have strong magnetic fields, so that the magnetospheres of newly formed close binaries can interact, dissipate energy, and produce synchrotron radiation. The V773 Tau A binary system, a pair of T Tauri stars with a 51 day orbit, displays such a signature, with peak emission taking place near periastron. This paper proposes that the observed emission arises from the change in energy stored in the composite magnetic field of the system. We model the fields using the leading order (dipole) components and show that this picture is consistent with current observations. In this model, the observed radiation accounts for a fraction of the available energy of interaction between the magnetic fields from the two stars. Assuming antisymmetry, we compute the interaction energy $E_{\rm int}$ as a function of the stellar radii, the stellar magnetic field strengths, the binary semi-major axis, and orbital eccentricity, all of which can be measured independently of the synchrotron radiation. The variability in time and energetics of the synchrotron radiation depend on the details of the annihilation of magnetic fields through reconnection events, which generate electric fields that accelerate charged particles, and how those charged particles, especially fast electrons, are removed from the interaction region. However, the major qualitative features are well described by the background changes in the global magnetic configuration driven by the orbital motion. The theory can be tested by observing a collection of pre-main-sequence binary systems. △ Less

Submitted 20 October, 2011; originally announced October 2011.

Comments: 32 pages, 5 figures, accepted to ApJ

Stability of Magnetized Disks and Implications for Planet Formation

Authors: Susana Lizano, Daniele Galli, Mike J. Cai, Fred C. Adams

Abstract: This paper considers gravitational perturbations in geometrically thin disks with rotation curves dominated by a central object, but with substantial contributions from magnetic pressure and tension. The treatment is general, but the application is to the circumstellar disks that arise during the gravitational collapse phase of star formation. We find the dispersion relation for spiral density wav… ▽ More This paper considers gravitational perturbations in geometrically thin disks with rotation curves dominated by a central object, but with substantial contributions from magnetic pressure and tension. The treatment is general, but the application is to the circumstellar disks that arise during the gravitational collapse phase of star formation. We find the dispersion relation for spiral density waves in these generalized disks and derive the stability criterion for axisymmetric $(m=0)$ disturbances (the analog of the Toomre parameter $Q_T$) for any radial distribution of the mass-to-flux ratio $λ$. The magnetic effects work in two opposing directions: on one hand, magnetic tension and pressure stabilize the disk against gravitational collapse and fragmentation; on the other hand, they also lower the rotation rate making the disk more unstable. For disks around young stars the first effect generally dominates, so that magnetic fields allow disks to be stable for higher surface densities and larger total masses. These results indicate that magnetic fields act to suppress the formation of giant planets through gravitational instability. Finally, even if gravitational instability can form a secondary body, it must lose an enormous amount of magnetic flux in order to become a planet; this latter requirement represents an additional constraint for planet formation via gravitational instability and places a lower limit on the electrical resistivity. △ Less

Submitted 28 September, 2010; originally announced September 2010.

Comments: accepted in ApJ

The Equilibrium Structure of Prolate Magnetized Molecular Cores

Authors: Michael J. Cai, Ronald E. Taam

Abstract: The structure of molecular cloud cores supported by thermal pressure and a poloidal magnetic field is reinvestigated in the magnetohydrostatic and axisymmetric approximation. In addition to oblate configurations found in earlier work, solutions yielding prolate spheroidal shapes have also been obtained for a reference state described by a uniform sphere threaded by a uniform background magnetic… ▽ More The structure of molecular cloud cores supported by thermal pressure and a poloidal magnetic field is reinvestigated in the magnetohydrostatic and axisymmetric approximation. In addition to oblate configurations found in earlier work, solutions yielding prolate spheroidal shapes have also been obtained for a reference state described by a uniform sphere threaded by a uniform background magnetic field. The solutions for prolate configurations are found to be relevant for lower masses than for their oblate counterparts. Of particular importance is the result that the prolate cloud cores have radii less than a maximum given by $0.25 pc (a/0.2km/s)^2 (P_{\rm ext}/10^{-12} dyne/cm^2)^{-1/2}$, where $a$ is the sound speed and $P_{\rm ext}$ is the external pressure of the background medium. The existence of such solutions obviates the presence of toroidal fields in such modeled structures. △ Less

Submitted 15 December, 2009; originally announced December 2009.

Comments: 5 pages, 3 figures. Accepted for publication by ApJL

Migration of Extrasolar Planets: Effects from X-Wind Accretion Disks

Authors: Fred C. Adams, Mike J. Cai, Susana Lizano

Abstract: Magnetic fields are dragged in from the interstellar medium during the gravitational collapse that forms star/disk systems. Consideration of mean field magnetohydrodynamics (MHD) in these disks shows that magnetic effects produce subkeplerian rotation curves and truncate the inner disk. This letter explores the ramifications of these predicted disk properties for the migration of extrasolar plan… ▽ More Magnetic fields are dragged in from the interstellar medium during the gravitational collapse that forms star/disk systems. Consideration of mean field magnetohydrodynamics (MHD) in these disks shows that magnetic effects produce subkeplerian rotation curves and truncate the inner disk. This letter explores the ramifications of these predicted disk properties for the migration of extrasolar planets. Subkeplerian flow in gaseous disks drives a new migration mechanism for embedded planets and modifies the gap opening processes for larger planets. This subkeplerian migration mechanism dominates over Type I migration for sufficiently small planets (m_P < 1 M_\earth) and/or close orbits (r < 1 AU). Although the inclusion of subkeplerian torques shortens the total migration time by only a moderate amount, the mass accreted by migrating planetary cores is significantly reduced. Truncation of the inner disk edge (for typical system parameters) naturally explains final planetary orbits with periods P=4 days. Planets with shorter periods P=2 days can be explained by migration during FU-Ori outbursts, when the mass accretion rate is high and the disk edge moves inward. Finally, the midplane density is greatly increased at the inner truncation point of the disk (the X-point); this enhancement, in conjunction with continuing flow of gas and solids through the region, supports the in situ formation of giant planets. △ Less

Submitted 10 August, 2009; originally announced August 2009.

Comments: 15 pages, 2 figures, accepted to ApJ Letters

Polarized Emission of Sagittarius A*

Authors: Lei Huang, Siming Liu, Zhi-Qiang Shen, Ye-Fei Yuan, Mike J. Cai, Hui Li, Christopher L. Fryer

Abstract: We explore the parameter space of the two temperature pseudo-Newtonian Keplerian accretion flow model for the millimeter and shorter wavelength emission from Sagittarius A*. A general relativistic ray-tracing code is used to treat the radiative transfer of polarized synchrotron emission from the flow. The synchrotron self-Comptonization and bremsstrahlung emission components are also included. I… ▽ More We explore the parameter space of the two temperature pseudo-Newtonian Keplerian accretion flow model for the millimeter and shorter wavelength emission from Sagittarius A*. A general relativistic ray-tracing code is used to treat the radiative transfer of polarized synchrotron emission from the flow. The synchrotron self-Comptonization and bremsstrahlung emission components are also included. It is shown that the model can readily account for the millimeter to sub-millimeter emission characteristics with an accretion rate of ~6x10^17g.s^-1 and an inclination angle of ~40 deg. However, the corresponding model predicted near-infrared and X-ray fluxes are more than one order of magnitude lower than the observed 'quiescent' state values. While the extended quiescent-state X-ray emission has been attributed to thermal emission from the large-scale accretion flow, the NIR emission and flares are likely dominated by emission regions either within the last stable orbit of a Schwarzschild black hole or associated with outflows. With the viscous parameter derived from numerical simulations, there is still a degeneracy between the electron heating rate and the magnetic parameter. A fully general relativistic treatment with the black hole spin incorporated will resolve these issues. △ Less

Submitted 30 July, 2009; originally announced July 2009.

Comments: 27 pages, 6 figures, ApJ accepted

The Challenge of Sub-Keplerian Rotation for Disk Winds

Authors: Frank H. Shu, Susana Lizano, Daniele Galli, Mike J. Cai, Subhanjoy Mohanty

Abstract: Strong magnetization makes the disks surrounding young stellar objects rotate at rates that are too sub-Keplerian to enable the thermal launching of disk winds from their surfaces unless the rate of gas diffusion across field lines is dynamically fast. This under-appreciated implication of disk magnetization poses a considerable challenge for disk-wind theory. Strong magnetization makes the disks surrounding young stellar objects rotate at rates that are too sub-Keplerian to enable the thermal launching of disk winds from their surfaces unless the rate of gas diffusion across field lines is dynamically fast. This under-appreciated implication of disk magnetization poses a considerable challenge for disk-wind theory. △ Less

Submitted 12 June, 2008; originally announced June 2008.

Comments: 4 pages, accepted by ApJL

Linearly and Circularly Polarized Emission in Sagittarius A*

Authors: L. Huang, S. Liu, Z. -Q. Shen, M. J. Cai, H. Li, C. L. Fryer

Abstract: We perform general relativistic ray-tracing calculations of the transfer of polarized synchrotron radiation through the relativistic accretion flow in Sagittarius (Sgr) A*. Based on a two-temperature magneto-rotational-instability (MRI) induced accretion mode, the birefringence effects are treated self-consistently. By fitting the spectrum and polarization of Sgr A* from millimeter to near-infra… ▽ More We perform general relativistic ray-tracing calculations of the transfer of polarized synchrotron radiation through the relativistic accretion flow in Sagittarius (Sgr) A*. Based on a two-temperature magneto-rotational-instability (MRI) induced accretion mode, the birefringence effects are treated self-consistently. By fitting the spectrum and polarization of Sgr A* from millimeter to near-infrared bands, we are able to not only constrain the basic parameters related to the MRI and the electron heating rate, but also limit the orientation of the accretion torus. These constraints lead to unique polarimetric images, which may be compared with future millimeter and sub-millimeter VLBI observations. In combination with general relativistic MHD simulations, the model has the potential to test the MRI with observations of Sgr A*. △ Less

Submitted 25 February, 2008; originally announced February 2008.

Comments: 12 pages, 2 figures, ApJL accepted

X-winds in Action

Authors: Mike J. Cai, Hsien Shang, Hsiao-Hsuan Lin, Frank H. Shu

Abstract: The interaction of accretion disks with the magnetospheres of young stars can produce X-winds and funnel flows. With the assumption of axial symmetry and steady state flow, the problem can be formulated in terms of quantities that are conserved along streamlines, such as the Bernoulli integral (BI), plus a partial differential equation (PDE), called the Grad-Shafranov equation (GSE), that govern… ▽ More The interaction of accretion disks with the magnetospheres of young stars can produce X-winds and funnel flows. With the assumption of axial symmetry and steady state flow, the problem can be formulated in terms of quantities that are conserved along streamlines, such as the Bernoulli integral (BI), plus a partial differential equation (PDE), called the Grad-Shafranov equation (GSE), that governs the distribution of streamlines in the meridional plane. The GSE plus BI yields a PDE of mixed type, elliptic before critical surfaces where the flow speed equals certain characteristic wave speeds are crossed and hyperbolic afterward. The computational difficulties are exacerbated by the locations of the critical surfaces not being known in advance. To overcome these obstacles, we consider a variational principle by which the GSE can be attacked by extremizing an action integral, with all other conserved quantities of the problem explicitly included as part of the overall formulation. To simplify actual applications we adopt the cold limit of a negligibly small ratio of the sound speed to the speed of Keplerian rotation in the disk where the X-wind is launched. We also ignore the obstructing effects of any magnetic fields that might thread a disk approximated to be infinitesimally thin. We then introduce trial functions with adjustable coefficients to minimize the variations that give the GSE. We tabulate the resulting coefficients so that other workers can have analytic forms to reconstruct X-wind solutions for various astronomical, cosmochemical, and meteoritical applications. △ Less

Submitted 19 September, 2007; originally announced September 2007.

Comments: 31 pages, 8 figures. Accepted for publication in ApJ

Black Hole Shadow Image and Visibility Analysis of Sagittarius A*

Authors: Lei Huang, Mike Cai, Zhi-Qiang Shen, Feng Yuan

Abstract: The compact dark objects with very large masses residing at the centres of galaxies are believed to be black holes. Due to the gravitational lensing effect, they would cast a shadow larger than their horizon size over the background, whose shape and size can be calculated. For the supermassive black hole candidate Sgr A*, this shadow spans an angular size of about 50 micro arc second, which is u… ▽ More The compact dark objects with very large masses residing at the centres of galaxies are believed to be black holes. Due to the gravitational lensing effect, they would cast a shadow larger than their horizon size over the background, whose shape and size can be calculated. For the supermassive black hole candidate Sgr A*, this shadow spans an angular size of about 50 micro arc second, which is under the resolution attainable with the current astronomical instruments. Such a shadow image of Sgr A* will be observable at about 1 mm wavelength, considering the scatter broadening by the interstellar medium. By simulating the black hole shadow image of Sgr A* with the radiatively inefficient accretion flow model, we demonstrate that analyzing the properties of the visibility function can help us determine some parameters of the black hole configuration, which is instructive to the sub-millimeter VLBI observations of Sgr A* in the near future. △ Less

Submitted 13 March, 2007; v1 submitted 12 March, 2007; originally announced March 2007.

Comments: 11 pages, 18 figures (.eps), MNRAS accepted

Journal ref: Mon.Not.Roy.Astron.Soc.379:833-840,2007