Showing 1–50 of 1,474 results for author: Chen, Y

NOEMA formIng Cluster survEy (NICE): Characterizing eight massive galaxy groups at $1.5 < z < 4$ in the COSMOS field

Authors: Nikolaj B. Sillassen, Shuowen Jin, Georgios E. Magdis, Emanuele Daddi, Tao Wang, Shiying Lu, Hanwen Sun, Vinod Arumugam, Daizhong Liu, Malte Brinch, Chiara D'Eugenio, Raphael Gobat, Carlos Gómez-Guijarro, Michael Rich, Eva Schinnerer, Veronica Strazzullo, Qinghua Tan, Francesco Valentino, Yijun Wang, Mengyuan Xiao, Luwenjia Zhou, David Blánquez-Sesé, Zheng Cai, Yanmei Chen, Laure Ciesla , et al. (19 additional authors not shown)

Abstract: The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are c… ▽ More The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are confirmed with ALMA, and one is confirmed by H$α$ from Subaru/FMOS. We constructed the integrated FIR SEDs for the eight groups, obtaining total IR SFR $=260-1300~{\rm M_\odot}$~yr$^{-1}$. We adopted six methods to estimate the dark matter masses, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass density. We found the radial stellar mass density are consistent with a NFW profile, supporting that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are $\log(M_{\rm h}/{\rm M_\odot})=12.8-13.7$ with uncertainty of 0.3 dex. From halo mass estimates, we derive baryonic accretion rate ${\rm BAR}=(1-8)\times10^{3}\,{\rm M_{\odot}/yr}$ for this sample. We find a quasi-linear correlation between the integrated SFR/BAR and the theoretical halo mass limit for cold streams, $M_{\rm stream}/M_{\rm h}$, with ${\rm SFR/BAR}=10^{-0.46\pm0.22}\left({M_{\rm stream}/M_{\rm h}}\right)^{0.71\pm0.16}$ with a scatter of $0.40\,{\rm dex}$. Further, we compare halo masses and stellar masses with simulations, and find all structures are consistent with being progenitors of $M_{\rm h}(z=0)>10^{14}\,{\rm M_{\odot}}$ galaxy clusters, and the most massive central galaxies have stellar masses consistent with brightest cluster galaxies (BCGs) progenitors in the TNG300 simulation. The results strongly suggest these structures are forming massive galaxy clusters via baryonic and dark matter accretion. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 44 pages (27pp appendix), 32 figures, 18 tables, accepted for publication in A&A

Mirage Sources and Large TeV Halo-Pulsar Offsets: Exploring the Parameter Space

Authors: Yiwei Bao, Ruo-Yu Liu, Gwenael Giacinti, Hai-Ming Zhang, Yang Chen

Abstract: We investigate the asymmetric propagation of 100 TeV electrons (whose radiation mainly concentrates on 20--30 TeV) in turbulent magnetic fields around pulsars, using GPU-accelerated simulations to explore their trajectories and interactions within pulsar wind nebulae and the interstellar medium. Key results include the identification of ``mirage'' sources indicating significant offsets in high-ene… ▽ More We investigate the asymmetric propagation of 100 TeV electrons (whose radiation mainly concentrates on 20--30 TeV) in turbulent magnetic fields around pulsars, using GPU-accelerated simulations to explore their trajectories and interactions within pulsar wind nebulae and the interstellar medium. Key results include the identification of ``mirage'' sources indicating significant offsets in high-energy emissions from their originating pulsars, challenging the results of traditional symmetric diffusion models. By varying parameters like source distance, magnetic field strength, and electron injection spectral index, the study delineates their effects on observable phenomena such as the probability that a source has at least one mirage around it, as well as the source separation. Our results offer insights into some puzzling sources observed recently by the Large High Altitude Air Shower Observatory (LHAASO), and shed light on the cosmic-ray transport mechanism in the interstellar medium. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Mirages and Large TeV Halo-Pulsar Offsets from Cosmic Ray Propagation

Authors: Yiwei Bao, Gwenael Giacinti, Ruo-Yu Liu, Hai-Ming Zhang, Yang Chen

Abstract: The study of extended $γ$-ray sources usually assumes symmetric diffusion of cosmic rays. However, recent observations of multiple sources near single pulsars and significant offsets between TeV halo centroids and their parent pulsars suggest that this assumption is overly simplistic. In this Letter, we demonstrate that asymmetric propagation of cosmic rays near their accelerators may create multi… ▽ More The study of extended $γ$-ray sources usually assumes symmetric diffusion of cosmic rays. However, recent observations of multiple sources near single pulsars and significant offsets between TeV halo centroids and their parent pulsars suggest that this assumption is overly simplistic. In this Letter, we demonstrate that asymmetric propagation of cosmic rays near their accelerators may create multiple TeV sources instead of a single symmetric source. This mechanism also explains the large offsets between TeV halo centroids and their pulsars. We demonstrate that several perplexing detected sources can be naturally explained without invoking additional invisible accelerators. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Baryon Acoustic Oscillations analyses with Density-Split Statistics

Authors: Tengpeng Xu, Yan-Chuan Cai, Yun Chen, Mark Neyrinck, Liang Gao, Qiao Wang

Abstract: Accurate modeling for the evolution of the Baryon Acoustic Oscillations (BAO) is essential for using it as a standard ruler to probe cosmology. We explore the non-linearity of the BAO in different environments using the density-split statistics and compare them to the case of the conventional two-point correlation function (2PCF). We detect density-dependent shifts for the position of the BAO with… ▽ More Accurate modeling for the evolution of the Baryon Acoustic Oscillations (BAO) is essential for using it as a standard ruler to probe cosmology. We explore the non-linearity of the BAO in different environments using the density-split statistics and compare them to the case of the conventional two-point correlation function (2PCF). We detect density-dependent shifts for the position of the BAO with respect to its linear version using halos from N-body simulations. Around low/high-densities, the scale of the BAO expands/contracts due to non-linear peculiar velocities. As the simulation evolves from redshift 1 to 0, the difference in the magnitude of the shifts between high- and low-density regions increases from the sub-percent to the percent level. In contrast, the scale of the BAO does not evolve in the total 2PCF in the same redshift range. The width of the BAO around high density regions increases as the universe evolves, similar to the known broadening of the BAO in the 2PCF due to non-linear evolution. In contrast, the width is smaller and stable for low density regions. We discuss possible implications for the reconstructions of the BAO in light of our results. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 16 pages, 10 figures

A Machine Learning Method for Monte Carlo Calculations of Radiative Processes

Authors: William Charles, Alexander Y. Chen

Abstract: Radiative processes such as synchrotron radiation and Compton scattering play an important role in astrophysics. Radiative processes are fundamentally stochastic in nature, and the best tools currently used for resolving these processes computationally are Monte Carlo (MC) methods. These methods typically draw a large number of samples from a complex distribution such as the differential cross sec… ▽ More Radiative processes such as synchrotron radiation and Compton scattering play an important role in astrophysics. Radiative processes are fundamentally stochastic in nature, and the best tools currently used for resolving these processes computationally are Monte Carlo (MC) methods. These methods typically draw a large number of samples from a complex distribution such as the differential cross section for electron-photon scattering, and then use these samples to compute the radiation properties such as angular distribution, spectrum, and polarization. In this work we propose a machine learning (ML) technique for efficient sampling from arbitrary known probability distributions that can be used to accelerate Monte Carlo calculation of radiative processes in astrophysical scenarios. In particular, we apply our technique to inverse Compton radiation and find that our ML method can be up to an order of magnitude faster than traditional methods currently in use. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 15 pages, 6 Figures

Concurrent Accretion and Migration of Giant Planets in their Natal Disks with Consistent Accretion Torque

Authors: Ya-Ping Li, Yi-Xian Chen, Douglas N. C. Lin

Abstract: Migration commonly occurs during the epoch of planet formation. For emerging gas giant planets, it proceeds concurrently with their growth through the accretion of gas from their natal protoplanetary disks. Similar migration process should also be applied to the stellar-mass black holes embedded in active galactic nucleus disks. In this work, we perform high resolution 3D and 2D numerical hydrodyn… ▽ More Migration commonly occurs during the epoch of planet formation. For emerging gas giant planets, it proceeds concurrently with their growth through the accretion of gas from their natal protoplanetary disks. Similar migration process should also be applied to the stellar-mass black holes embedded in active galactic nucleus disks. In this work, we perform high resolution 3D and 2D numerical hydrodynamical simulations to study the migration dynamics for accreting embedded objects over the disk viscous timescales in a self-consistent manner. We find that an accreting planet embedded in a predominantly viscous disk has a tendency to migrate outward, in contrast to the inward orbital decay of non-accreting planets. 3D and 2D simulations find the consistent outward migration results for the accreting planets. Under this circumstance, the accreting planet's outward migration is mainly due to the asymmetric spiral arms feeding from the global disk into the Hill radius. This is analogous to the unsaturated corotation torque although the imbalance is due to material accretion within the libration timescale rather than diffusion onto the inner disk. In a disk with a relatively small viscosity, the accreting planets clear deep gaps near their orbits. The tendency of inward migration is recovered, albeit with suppressed rates. By performing a parameter survey with a range of disks' viscosity, we find that the transition from outward to inward migration occurs with the effective viscous efficiency factor $α\sim 0.003$ for Jupiter-mass planets. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures, Accepted by ApJ

Spatially resolved analysis of Stellar Populations in NGC 2992: Impact of AGN feedback

Authors: Xiaoyu Xu, Junfeng Wang, Zhiyuan Li, Yanmei Chen

Abstract: In NGC 2992, a galaxy-scale ionized gas outflow driven by AGN has long been recognized, yet its impact on the host galaxy has remained elusive. In this paper, we utilize data from the archival Very Large Telescope (VLT)/MUSE to present a spatially resolved analysis of stellar populations in this galaxy. Two different stellar population templates are employed to fit the stellar continuum, allowing… ▽ More In NGC 2992, a galaxy-scale ionized gas outflow driven by AGN has long been recognized, yet its impact on the host galaxy has remained elusive. In this paper, we utilize data from the archival Very Large Telescope (VLT)/MUSE to present a spatially resolved analysis of stellar populations in this galaxy. Two different stellar population templates are employed to fit the stellar continuum, allowing us to determine the light-weighted stellar age, metallicity, the fraction of the young stellar population (age $<100$ Myr, $P_{\rm Y}$), and the average age and metallicity of $P_{\rm Y}$. Our results reveal the presence of a very young stellar population ($\leq40$ Myr) within the dust lane and nearly along the galaxy's major axis. The light-weighted stellar age and the fraction of $P_{\rm Y}$ show negative trends along the major and minor axes. The average age and metallicity of $P_{\rm Y}$ present positive trends with increasing distance, except along the northern direction of the major axis. Within the circumnuclear region ($<1$ kpc), the distribution of the young stellar population is spatially anti-correlated with the AGN outflow cone. The highest fraction of $P_{\rm Y}$ is observed at the outskirts of the nuclear radio bubble in the northern region near the nucleus. Considering the coupling efficiency and timescales, we propose that the AGN outflow in this galaxy may exert both negative and positive feedback on its host. Additionally, the star formation and the AGN activities could be attributed to the interaction between NGC 2992 and NGC 2993. △ Less

Submitted 13 June, 2024; originally announced June 2024.

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

Constraints on Ultra Heavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes… ▽ More In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $γ$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Simulation Models for Exploring Magnetic Reconnection

Authors: Michael Shay, Subash Adhikari, Naoki Beesho, Joachim Birn, Jorg Buechner, Paul Cassak, Li-Jen Chen, Yuxi Chen, Giulia Cozzani, Jim Drake, Fan Guo, Michael Hesse, Neeraj Jain, Yann Pfau-Kempf, Yu Lin, Yi-Hsin Liu, Mitsuo Oka, Yuri A. Omelchenko, Minna Palmroth, Oreste Pezzi, Patricia H. Reiff, Marc Swisdak, Frank Toffoletto, Gabor Toth, Richard A. Wolf

Abstract: Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these… ▽ More Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM). △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: Chapter 5.2 of ISSI Book on Magnetic Reconnection, submitted to Space Science Reviews

Prevalence of non-standard collapsing of strong Langmuir turbulence in solar corona plasmas

Authors: Yaokun Li, Haomin Sun, Hao Ning, Sulan Ni, Xiangliang Kong, Jiansen He, Yao Chen

Abstract: We present a fully-kinetic simulation of the full life cycle of strong Langmuir turbulence (SLT) excited by electron beams that are accelerated under the solar corona conditions. We find that (1) most packets ($\sim$80%) are affected by their neighbors during their collapse, as a result, their spatial scale variations present non-standard evolutionary features, i.e., deviating away from what was p… ▽ More We present a fully-kinetic simulation of the full life cycle of strong Langmuir turbulence (SLT) excited by electron beams that are accelerated under the solar corona conditions. We find that (1) most packets ($\sim$80%) are affected by their neighbors during their collapse, as a result, their spatial scale variations present non-standard evolutionary features, i.e., deviating away from what was predicted by the Zakharov model; (2) the collapsing cavity is too shallow to trap the wave packet due to the growth of the Coulomb force, as a result a majority ($\sim$70%) of the packet energy runs away and a secondary localization may occur. The study indicates that the non-standard Langmuir collapse may play an important role in coronal plasmas interacting with an intense electron beam, that may be eventually confirmed by humanity's first mission to fly through the corona. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Modeling the refractive index profile n(z) of polar ice for ultra-high energy neutrino experiments

Authors: S. Ali, P. Allison, S. Archambault, J. J. Beatty, D. Z. Besson, A. Bishop, P. Chen, Y. C. Chen, B. A. Clark, W. Clay, A. Connolly, K. Couberly, L. Cremonesi, A. Cummings, P. Dasgupta, R. Debolt, S. de Kockere, K. D. de Vries, C. Deaconu, M. A. DuVernois, J. Flaherty, E. Friedman, R. Gaior, P. Giri, J. Hanson , et al. (45 additional authors not shown)

Abstract: We develop an in-situ index of refraction profile using the transit time of radio signals broadcast from an englacial transmitter to 2-5 km distant radio-frequency receivers, deployed at depths up to 200 m. Maxwell's equations generally admit two ray propagation solutions from a given transmitter, corresponding to a direct path (D) and a refracted path (R); the measured D vs. R (dt(D,R)) timing di… ▽ More We develop an in-situ index of refraction profile using the transit time of radio signals broadcast from an englacial transmitter to 2-5 km distant radio-frequency receivers, deployed at depths up to 200 m. Maxwell's equations generally admit two ray propagation solutions from a given transmitter, corresponding to a direct path (D) and a refracted path (R); the measured D vs. R (dt(D,R)) timing differences provide constraints on the index of refraction profile near South Pole, where the Askaryan Radio Array (ARA) neutrino observatory is located. We constrain the refractive index profile by simulating D and R ray paths via ray tracing and comparing those to measured dt(D,R) signals. Using previous ice density data as a proxy for n(z), we demonstrate that our data strongly favors a glaciologically-motivated three-phase densification model rather than a single exponential scale height model. Simulations show that the single exponential model overestimates ARA neutrino sensitivity compared to the three-phase model. △ Less

Submitted 11 June, 2024; v1 submitted 2 June, 2024; originally announced June 2024.

Filament eruption by multiple reconnections

Authors: Y. Liu, G. P. Ruan, B. Schmieder, J. H. Guo, Y. Chen, R. S. Zheng, J. T. Su, B. Wang

Abstract: Filament eruption is a common phenomenon in solar activity, but the triggering mechanism is not well understood. We focus our study on a filament eruption located in a complex nest of three active regions close to a coronal hole. The filament eruption is observed at multiple wavelengths: by the GONG, the STEREO, the SUTRI, and the AIA and Helioseismic and Magnetic Imager (HMI) on board the SDO. Th… ▽ More Filament eruption is a common phenomenon in solar activity, but the triggering mechanism is not well understood. We focus our study on a filament eruption located in a complex nest of three active regions close to a coronal hole. The filament eruption is observed at multiple wavelengths: by the GONG, the STEREO, the SUTRI, and the AIA and Helioseismic and Magnetic Imager (HMI) on board the SDO. Thanks to high temporal-resolution observations, we were able to analyze the evolution of the fine structure of the filament in detail. The filament changes direction during the eruption, which is followed by a halo coronal mass ejection detected by the LASCO on board the SOHO. A Type III radio burst was also registered at the time of the eruption. To investigate the process of the eruption, we analyzed the magnetic topology of the filament region adopting a nonlinear force-free-field (NLFFF) extrapolation method and the polytropic global magnetohydrodynamic (MHD) modeling. We modeled the filament by embeddingatwisted fluxropewiththe regularized Biot-Savart Laws (RBSL) method in the ambient magnetic f ield. The extrapolation results show that magnetic reconnection occurs in a fan-spine configuration resulting in a circular flare ribbon. The global modeling of the corona demonstrates that there was an interaction between the filament and open field lines, causing a deflection of the filament in the direction of the observed CME eruption and dimming area. The modeling supports the following scenario: magnetic reconnection not only occurs with the filament itself (the flux rope) but also with the background magnetic field lines and open field lines of the coronal hole located to the east of the flux rope. This multiwavelength analysis indicates that the filament undergoes multiple magnetic reconnections on small and large scales with a drifting of the flux rope. △ Less

Submitted 2 June, 2024; originally announced June 2024.

Electron acceleration and transport in the 2023-03-06 solar flare

Authors: Alexey Kuznetsov, Zhao Wu, Sergey Anfinogentov, Yang Su, Yao Chen

Abstract: We investigated in detail the M5.8 class solar flare that occurred on 2023-03-06. This flare was one of the first strong flares observed by the Siberian Radioheliograph in the microwave range and the Advanced Space-based Solar Observatory in the X-ray range. The flare consisted of two separate flaring events (a "thermal" and a "cooler" ones), and was associated with (and probably triggered by) a f… ▽ More We investigated in detail the M5.8 class solar flare that occurred on 2023-03-06. This flare was one of the first strong flares observed by the Siberian Radioheliograph in the microwave range and the Advanced Space-based Solar Observatory in the X-ray range. The flare consisted of two separate flaring events (a "thermal" and a "cooler" ones), and was associated with (and probably triggered by) a filament eruption. During the first part of the flare, the microwave emission was produced in an arcade of relatively short and low flaring loops. During the second part of the flare, the microwave emission was produced by energetic electrons trapped near the top of a large-scale flaring loop; the evolution of the trapped electrons was mostly affected by the Coulomb collisions. Using the available observations and the GX Simulator tool, we created a 3D model of the flare, and estimated the parameters of the energetic electrons in it. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Frontiers in Astronomy and Space Sciences - Stellar and Solar Physics

A two-phase model of galaxy formation: III. The formation of globular clusters

Authors: Yangyao Chen, Houjun Mo, Huiyuan Wang

Abstract: We develop a model of globular cluster (GC) formation within the cosmological hierarchy of structure formation. The model is rooted in the `two-phase' scenario of galaxy formation developed in Paper-I, where the fast accretion of dark matter halos at high redshift leads to the formation of self-gravitating, turbulent gas clouds that subsequently fragment into dynamically hot systems of dense sub-c… ▽ More We develop a model of globular cluster (GC) formation within the cosmological hierarchy of structure formation. The model is rooted in the `two-phase' scenario of galaxy formation developed in Paper-I, where the fast accretion of dark matter halos at high redshift leads to the formation of self-gravitating, turbulent gas clouds that subsequently fragment into dynamically hot systems of dense sub-clouds with masses $\sim 10^6$-$10^7 M_\odot$. Here we elaborate on the formation, evolution, and fate of these sub-clouds, and show that some of the sub-clouds can be compactified via two distinctive channels into a 'supernova-free' regime to form two distinct populations of GCs. The model is simple, characterized by a small number of free parameters underpinned by physical considerations, and can be efficiently implemented into cosmological N-body simulations to generate a coherent sample of halos, galaxies, and GCs. Our model can reproduce a range of observations on GCs, including the mass function, the size-mass relation, the frequency per unit host galaxy/halo mass, the bimodal metallicity distribution, and the spatial profile. Predictions for GCs are made for both the local Universe and for redshift up to $z \approx 10$, and can be tested by upcoming observations. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 35 pages, 19 figures, 1 table. Submitted to MNRAS

Robust Nitrogen and Oxygen Abundances of Haro 3 from Optical and Infrared Emission

Authors: Yuguang Chen, Tucker Jones, Ryan L. Sanders, Dario Fadda, Jessica Sutter, Robert Minchin, Nikolaus Z. Prusinski, Sunny Rhoades, Keerthi Vasan GC, Charles C. Steidel, Erin Huntzinger, Paige Kelly, Danielle A. Berg, Fabio Bresolin, Rodrigo Herrera-Camus, Ryan J. Rickards Vaught, Guido Roberts-Borsani, Peter Senchyna, Justin S. Spilker, Daniel P. Stark, Benjamin Weiner, D. Christopher Martin, Mateusz Matuszewski, Rosalie C. McGurk, James D. Neill

Abstract: Accurate chemical compositions of star-forming regions are a critical diagnostic tool to characterize the star formation history and gas flows which regulate galaxy formation. However, the abundance discrepancy factor (ADF) between measurements from the "direct" optical electron temperature ($T_e$) method and from the recombination lines (RL) represents $\sim0.2$ dex systematic uncertainty in oxyg… ▽ More Accurate chemical compositions of star-forming regions are a critical diagnostic tool to characterize the star formation history and gas flows which regulate galaxy formation. However, the abundance discrepancy factor (ADF) between measurements from the "direct" optical electron temperature ($T_e$) method and from the recombination lines (RL) represents $\sim0.2$ dex systematic uncertainty in oxygen abundance. The degree of uncertainty for other elements is unknown. We conduct a comprehensive analysis of O$^{++}$ and N$^+$ ion abundances using optical and far-infrared spectra of a star-forming region within the nearby dwarf galaxy Haro 3, which exhibits a typical ADF. Assuming homogeneous conditions, the far-IR emission indicates an O abundance which is higher than the $T_e$ method and consistent with the RL value, as would be expected from temperature fluctuations, whereas the N abundance is too large to be explained by temperature fluctuations. Instead a component of highly obscured gas is likely required to explain the high far-IR to optical flux ratios. Accounting for this obscured component reduces both the IR-based metallicities and the inferred magnitude of temperature fluctuations, such that they cannot fully explain the ADF in Haro 3. Additionally, we find potential issues when predicting the RL fluxes from current atomic data. Our findings underscore the critical importance of resolving the cause of abundance discrepancies and understanding the biases between different metallicity methods. This work represents a promising methodology, and we identify further approaches to address the current dominant uncertainties. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 18 pages, 7 figures, and 2 tables. Submitted to ApJ

Forecasting the Population of Globular Cluster Streams in Milky Way-type Galaxies

Authors: Sarah Pearson, Ana Bonaca, Yingtian Chen, Oleg Y. Gnedin

Abstract: Thin stellar streams originating from globular clusters are among the most sensitive tracers of low-mass dark-matter subhalos. Joint analysis of the entire population of stellar streams will place the most robust constraints on the dark-matter subhalo mass function, and therefore the nature of dark matter. Here we use a hierarchical model of globular cluster formation to forecast the total number,… ▽ More Thin stellar streams originating from globular clusters are among the most sensitive tracers of low-mass dark-matter subhalos. Joint analysis of the entire population of stellar streams will place the most robust constraints on the dark-matter subhalo mass function, and therefore the nature of dark matter. Here we use a hierarchical model of globular cluster formation to forecast the total number, masses and radial distribution of dissolved globular cluster in Milky Way-like galaxies. Furthermore, we generate mock stellar streams from these progenitors' orbital histories taking into account the clusters' formation and accretion time, mass, and metallicity. Out of $\sim$10,000 clusters more massive than $10^4$ M$_{\odot}$, $\sim$9000 dissolved in the central bulge and are fully phase-mixed at the present, while the remaining $\sim$1000 survive as coherent stellar streams. This suggests that the current census of $\sim$80 globular cluster streams in the Milky Way is severely incomplete. Beyond 15 kpc from the Galactic center we are missing hundreds of streams, of which the vast majority are from accreted GCs. Deep Rubin photometry $(g\lesssim27.5)$ would be able to detect these streams, even the most distant ones beyond $> 75$ kpc. We also find that M31 will have an abundance of streams at galactocentric radii of 30-100 kpc. We conclude that future surveys will find a multitude of stellar streams from globular clusters which can be used for dark matter subhalo searches. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: 21 pages, 13 figures, submitted to AAS journals, comments welcome

Ly$α$ Halo Properties and Dust in the Circumgalactic Medium of $z \sim 2$ Star-forming Galaxies

Authors: Zhiyuan Song, Naveen A. Reddy, Yuguang Chen, Alice E. Shapley, Saeed Rezaee, Andrew Weldon, Tara Fetherolf, Alison L. Coil, Bahram Mobasher, Charles C. Steidel

Abstract: We present Keck Cosmic Web Imager IFU observations around extended Ly$α$ halos of 27 typical star-forming galaxies with redshifts $2.0 < z < 3.2$ drawn from the MOSFIRE Deep Evolution Field survey. We examine the average Ly$α$ surface-brightness profiles in bins of star-formation rate (SFR), stellar mass ($M_*$), age, stellar continuum reddening, SFR surface density ($\rm Σ_{SFR}$), and… ▽ More We present Keck Cosmic Web Imager IFU observations around extended Ly$α$ halos of 27 typical star-forming galaxies with redshifts $2.0 < z < 3.2$ drawn from the MOSFIRE Deep Evolution Field survey. We examine the average Ly$α$ surface-brightness profiles in bins of star-formation rate (SFR), stellar mass ($M_*$), age, stellar continuum reddening, SFR surface density ($\rm Σ_{SFR}$), and $\rm Σ_{SFR}$ normalized by stellar mass ($\rm Σ_{sSFR}$). The scale lengths of the halos correlate with stellar mass, age, and stellar continuum reddening; and anti-correlate with star-formation rate, $\rm Σ_{SFR}$, and $\rm Σ_{sSFR}$. These results are consistent with a scenario in which the down-the-barrel fraction of Ly$α$ emission is modulated by the low-column-density channels in the ISM, and that the neutral gas covering fraction is related to the physical properties of the galaxies. Specifically, we find that this covering fraction increases with stellar mass, age, and $E(B-V)$; and decreases with SFR, $\rm Σ_{SFR}$ and $\rm Σ_{sSFR}$. We also find that the resonantly scattered Ly$α$ emission suffers greater attenuation than the (non-resonant) stellar continuum emission, and that the difference in attenuation increases with stellar mass, age, and stellar continuum reddening, and decreases with $\rm Σ_{sSFR}$. These results imply that more reddened galaxies have more dust in their CGM. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 12 pages, 7 figures, 2 tables, submitted to ApJ

Solar Imaging Data Analytics: A Selective Overview of Challenges and Opportunities

Authors: Yang Chen, Ward Manchester, Meng Jin, Alexei Pevtsov

Abstract: We give a gentle introduction to solar imaging data, focusing on the challenges and opportunities of data-driven approaches for solar eruptions. The various solar phenomenon prediction problems that might benefit from statistical methods are presented. Available data products and software are described. State-of-the-art solar eruption forecasting models with data-driven approaches are summarized a… ▽ More We give a gentle introduction to solar imaging data, focusing on the challenges and opportunities of data-driven approaches for solar eruptions. The various solar phenomenon prediction problems that might benefit from statistical methods are presented. Available data products and software are described. State-of-the-art solar eruption forecasting models with data-driven approaches are summarized and discussed. Based on the characteristics of the datasets and state-of-the-art approaches, we point out several promising directions to explore from statistical modeling and computational perspectives. △ Less

Submitted 2 July, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Data quality control system and long-term performance monitor of the LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, H. X. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (263 additional authors not shown)

Abstract: The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To… ▽ More The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Characterization of carbon dioxide on Ganymede and Europa supported by experiments: Effects of temperature, porosity, and mixing with water

Authors: L. Schiltz, B. Escribano, G. M. Muñoz Caro, S. Cazaux, C. del Burgo Olivares, H. Carrascosa, I. Boszhuizen, C. González Díaz, Y. -J. Chen, B. M. Giuliano, P. Caselli

Abstract: The surfaces of icy moons are primarily composed of water ice that can be mixed with other compounds, such as carbon dioxide. The carbon dioxide (CO$_2)$ stretching fundamental band observed on Europa and Ganymede appears to be a combination of several bands that are shifting location from one moon to another. We investigate the cause of the observed shift in the CO$_2$ stretching absorption band… ▽ More The surfaces of icy moons are primarily composed of water ice that can be mixed with other compounds, such as carbon dioxide. The carbon dioxide (CO$_2)$ stretching fundamental band observed on Europa and Ganymede appears to be a combination of several bands that are shifting location from one moon to another. We investigate the cause of the observed shift in the CO$_2$ stretching absorption band experimentally. We also explore the spectral behaviour of CO$_2$ ice by varying the temperature and concentration.} %H$_2$O:CO$_2$ deposition ratios. We analyzed pure CO$_2$ ice and ice mixtures deposited at 10 K under ultra-high vacuum conditions using Fourier-transform infrared (FTIR) spectroscopy and temperature programmed desorption (TPD) experiments. Laboratory ice spectra were compared to JWST observation of Europa's and Ganymede's leading hemispheres. The simulated IR spectra were calculated using density functional theory (DFT) methods, exploring the effect of porosity in CO$_2$ ice. Pure CO$_2$ and CO$_2$-water ice show distinct spectral changes and desorption behaviours at different temperatures, revealing intricate CO$_2$ and H$_2$O interactions. The number of discernible peaks increases from two in pure CO$_2$ to three in CO$_2$-water mixtures. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Regions of suppressed diffusion around supernova remnants?

Authors: Yiwei Bao, Pasquale Blasi, Yang Chen

Abstract: The recent discovery of the so-called TeV halos has attracted much attention. The morphology of the emission requires that the region is characterized by severe suppression of the diffusion coefficient. This finding raises many questions as to its origin: 1) is the suppressed diffusion to be attributed to instabilities induced by the same radiating particles? 2) or does it actually show that the d… ▽ More The recent discovery of the so-called TeV halos has attracted much attention. The morphology of the emission requires that the region is characterized by severe suppression of the diffusion coefficient. This finding raises many questions as to its origin: 1) is the suppressed diffusion to be attributed to instabilities induced by the same radiating particles? 2) or does it actually show that the diffusion coefficient is small throughout the disc of the Galaxy? In both cases, one would expect that the surroundings of supernova remnants (SNRs) should also show evidence of reduced diffusion coefficient, since most remnants are located in the disc and are expected to be sites of effective particle acceleration. Should we expect the existence of regions of extended $γ$-ray emission from these regions as well? Here we investigate the transport of cosmic rays (CRs) escaped from SNRs in order to assess the viability of the idea of having a cocoon of suppressed diffusion around them. A comparison of our results with the $γ$-ray emission from the regions around HB9 and W28 does not provide solid evidence of reduced diffusivity. However, if indeed the phenomenon of reduced diffusivity occurs around SNRs surrounded by molecular clouds, our calculations show that the effects on the grammage of Galactic CRs can be significant. △ Less

Submitted 18 May, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: published in ApJ

Journal ref: 2024 ApJ 966 224

Discovery of Very-high-energy Gamma-ray Emissions from the Low Luminosity AGN NGC 4278 by LHAASO

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) i… ▽ More The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $γ$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$σ$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

A Generic Model for Persistent Radio Source around Fast Radio Bursts

Authors: Yushan Chen, Hao Tong

Abstract: The repeated fast radio burst FRB 121102A and FRB 190520B has been reported, along with a spatially coincident, compact, persistent radio emission. In this paper, we present a parameterized one-zone model, with a basic scenario that a relativistic magnetized wind from the pulsar sweeps up the surroundings, e.g. freely expanding supernova ejecta, giving rise to a power-law distribution of electron… ▽ More The repeated fast radio burst FRB 121102A and FRB 190520B has been reported, along with a spatially coincident, compact, persistent radio emission. In this paper, we present a parameterized one-zone model, with a basic scenario that a relativistic magnetized wind from the pulsar sweeps up the surroundings, e.g. freely expanding supernova ejecta, giving rise to a power-law distribution of electron filled between the forward shock and the termination shock. We show that via appropriate adjustment of the model parameters, we can obtain the synchrotron radio emission properties from the one-zone model bright enough to account for observation, simply and analytically fitting the observed spectra well. Through dynamical evolution of the model, we can also obtain time-varying of relevant properties. This parameterized model does not depend on concrete physical models such as central engine, instead we can constraint physical model via comparison between parameters and observation, indicating the information about the central engine and surroundings. We also discuss the synchrotron self-Compton emission in our scenario in the end, but find no clue on the counterparts at other waveband. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: 12 pages,5 figures, published by ApJ

Journal ref: ApJ, Volume 966, 2024, Issue 2, id.179, 9 pp

Galaxies with Biconical Ionized Structure in MaNGA - I. Sample Selection and Driven Mechanisms

Authors: Zhi-Jie Zhou, Yan-Mei Chen, Run-Quan Guan, Yong Shi, Qiu-Sheng Gu, Dmitry Bizyaev

Abstract: Based on the integral field unit (IFU) data from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we develop a new method to select galaxies with biconical ionized structures, building a sample of 142 edge-on biconical ionized galaxies. We classify these 142 galaxies into 81 star-forming galaxies, 31 composite galaxies, and 30 AGNs (consisting of 23 Seyferts and 7 LI(N)ERs) acco… ▽ More Based on the integral field unit (IFU) data from Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we develop a new method to select galaxies with biconical ionized structures, building a sample of 142 edge-on biconical ionized galaxies. We classify these 142 galaxies into 81 star-forming galaxies, 31 composite galaxies, and 30 AGNs (consisting of 23 Seyferts and 7 LI(N)ERs) according to the {\nii}-BPT diagram. The star-forming bicones have bar-like structures while AGN bicones display hourglass structures, and composite bicones exhibit transitional morphologies between them due to both black hole and star-formation activities. Star-forming bicones have intense star-formation activities in their central regions, and the primary driver of biconical structures is the central star formation rate surface density. The lack of difference in the strength of central black hole activities (traced by dust attenuation corrected {\oiii}$λ$5007 luminosity and Eddington ratio) between Seyfert bicones and their control samples can be naturally explained as that the accretion disk and the galactic disk are not necessarily coplanar. Additionally, the biconical galaxies with central LI(N)ER-like line ratios are edge-on disk galaxies that show strong central dust attenuation. The radial gradients of {\ha} surface brightness follow the $r^{-2.35}$ relation, roughly consistent with $r^{-2}$ profile, which is expected in the case of photoionization by a central point-like source. These observations indicate obscured AGNs or AGN echoes as the primary drivers of biconical structures in LI(N)ERs. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Comments: 12 pages, 9 figures, 1 table, Accepted for publication in MNRAS

A Multi-Peak Solar Flare with a High Turnover Frequency of The Gyrosynchrotron Spectra from the Loop-Top Source

Authors: Zhao Wu, Alexey Kuznetsov, Sergey Anfinogentov, Victor Melnikov, Robert Sych, Bing Wang, Ruisheng Zheng, Xiangliang Kong, Baolin Tan, Zongjun Ning, Yao Chen

Abstract: The origin of multiple peaks in lightcurves of various wavelengths remains illusive during flares. Here we discuss the flare of SOL2023-05-09T03:54M6.5 with six flux peaks as recorded by a tandem of new microwave and Hard X-ray instruments. According to its microwave spectra, the flare represents a high-turnover frequency (>15 GHz) event. The rather-complete microwave and HXR spectral coverage pro… ▽ More The origin of multiple peaks in lightcurves of various wavelengths remains illusive during flares. Here we discuss the flare of SOL2023-05-09T03:54M6.5 with six flux peaks as recorded by a tandem of new microwave and Hard X-ray instruments. According to its microwave spectra, the flare represents a high-turnover frequency (>15 GHz) event. The rather-complete microwave and HXR spectral coverage provides a rare opportunity to uncover the origin of such event together with simultaneous EUV images. We concluded that (1) the microwave sources originates around the top section of the flaring loops with a trend of source spatial dispersion with frequency;(2) the visible movement of the microwave source from peak to peak originates from the process of new flaring loops appearing sequentially along the magnetic neutral line; 3) the optically-thin microwave spectra are hard with the indices varying from -1.2 to -0.4, and the turnover frequency always exceeds 15 GHz; 4) higher turnover/peak frequency corresponds to stronger peak intensity and harder optically-thin spectra. Using the Fokker-Planck and GX simulator codes we obtained a good fit to the observed microwave spectra and spatial distribution of the sources at all peaks, if assuming the radiating energetic electrons have the same spatial distribution and single-power-law spectra but with the number density varying in a range of 30%. We conclude that the particle acceleration in this flare happens in a compact region nearing the looptop. These results provide new constraints on the acceleration of energetic electrons and the underlying flare intermittent reconnection process. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Comments: 23 pages, 11 figures

Four-hundred Very Metal-poor Stars Studied with LAMOST and Subaru. III. Dynamically Tagged Groups and Chemodynamical Properties

Authors: Ruizhi Zhang, Tadafumi Matsuno, Haining Li, Wako Aoki, Xiang-Xiang Xue, Takuma Suda, Gang Zhao, Yuqin Chen, Miho N. Ishigaki, Jianrong Shi, Qianfan Xing, Jingkun Zhao

Abstract: Very metal-poor (VMP) stars record the signatures of early accreted galaxies, making them essential tools for unraveling the early stages of Galaxy formation. Understanding the origin of VMP stars requires comprehensive studies of their chemical compositions and kinematics, which are currently lacking. Hence, we conduct a chemodynamical analysis of 352 VMP stars selected from one of the largest un… ▽ More Very metal-poor (VMP) stars record the signatures of early accreted galaxies, making them essential tools for unraveling the early stages of Galaxy formation. Understanding the origin of VMP stars requires comprehensive studies of their chemical compositions and kinematics, which are currently lacking. Hence, we conduct a chemodynamical analysis of 352 VMP stars selected from one of the largest uniform high-resolution VMP star samples, jointly obtained from LAMOST and Subaru. We apply a friends-of-friends clustering algorithm to the master catalog of this high-resolution sample, which consists of 5778 VMP stars. It results in 131 dynamically tagged groups with 89 associated with known substructures in the Milky Way, including Gaia-Sausage-Enceladus (GSE), Thamnos, Helmi streams, Sequoia, Wukong, Pontus, and the very metal-poor disk (VMPD). Our findings are: (i) the VMPD shows lower Zn abundances than the rest, which indicates that it could be a relic of small stellar systems; (ii) Sequoia shows moderately high r-process abundances; (iii) Helmi streams show deficiencies in carbon and light neutron-capture elements; (iv) the fraction of carbon-enhanced metal-poor stars with no enhancement in heavy elements (CEMP-no stars) seems low in the VMPD and the Helmi streams; and (v) a subgroup in GSE exhibits a very high fraction of r-process enhanced stars, with four out of five showing [Eu/Fe]> +1.0. The abundance patterns of other elements in VMP substructures largely match the whole VMP sample. We also study large-scale correlations between abundance ratios and kinematics without classifying stars into substructures, but it does not yield significant correlations once the overall chemical evolution is considered for most elements. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Comments: 30 pages, 24 figures, published in ApJ

Alfvén Wave Conversion to Low Frequency Fast Magnetosonic Waves in Magnetar Magnetospheres

Authors: Dominic Bernardi, Yajie Yuan, Alexander Y. Chen

Abstract: Rapid shear motion of magnetar crust can launch Alfvén waves into the magnetosphere. The dissipation of the Alfvén waves has been theorized to power the X-ray bursts characteristic of magnetars. However, the process by which Alfvén waves convert their energy to X-rays is unclear. Recent work has suggested that energetic fast magnetosonic (fast) waves can be produced as a byproduct of Alfvén waves… ▽ More Rapid shear motion of magnetar crust can launch Alfvén waves into the magnetosphere. The dissipation of the Alfvén waves has been theorized to power the X-ray bursts characteristic of magnetars. However, the process by which Alfvén waves convert their energy to X-rays is unclear. Recent work has suggested that energetic fast magnetosonic (fast) waves can be produced as a byproduct of Alfvén waves propagating on curved magnetic field lines; their subsequent dissipation may power X-ray bursts. In this work, we investigate the production of fast waves by performing axisymmetric force-free simulations of Alfvén waves propagating in a dipolar magnetosphere. For Alfvén wave trains that do not completely fill the flux tube confining them, we find a fast wave dominated by a low frequency component with a wavelength defined by the bouncing time of the Alfvén waves. In contrast, when the wave train is long enough to completely fill the flux tube, and the Alfvén waves overlap significantly, the energy is quickly converted into a fast wave with a higher frequency that corresponds to twice the Alfvén wave frequency. We investigate how the energy, duration, and wavelength of the initial Alfvén wave train affect the conversion efficiency to fast waves. For modestly energetic star quakes, we see that the fast waves that are produced will become non-linear well within the magnetosphere, and we comment on the X-ray emission that one may expect from such events. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 12 pages, 13 figures

Comparison of Ion-Proton Differential Speed between ICMEs and Solar Wind near 1 au

Authors: Xuechao Zhang, Hongqiang Song, Chengxiao Zhang, Hui Fu, Leping Li, Jinrong Li, Xiaoqian Wang, Rui Wang, Yao Chen

Abstract: The elemental abundance of ICMEs and solar wind near 1 au is often adopted to represent the abundance in the corresponding coronal sources. However, the absolute abundance of heavy ions (relative to hydrogen) near 1 au might be different from the coronal abundance due to the ion-proton differential speed ($V_{ip}$). To illustrate the $V_{ip}$ characteristics and explore whether it influences the a… ▽ More The elemental abundance of ICMEs and solar wind near 1 au is often adopted to represent the abundance in the corresponding coronal sources. However, the absolute abundance of heavy ions (relative to hydrogen) near 1 au might be different from the coronal abundance due to the ion-proton differential speed ($V_{ip}$). To illustrate the $V_{ip}$ characteristics and explore whether it influences the absolute abundance analysis for ICMEs and solar wind, we perform a statistical study on the $V_{ip}$ for He$^{2+}$, C$^{5+}$, O$^{6+}$, and Fe$^{10+}$ in both ICMEs and solar wind based on measurements of Advanced Composition Explorer. The results show that the $V_{ip}$ is negligible within ICMEs and slow solar wind ($<$ 400 km s$^{-1}$), while obvious in the intermediate (400 -- 600 km s$^{-1}$) and fast wind ($>$ 600 km s$^{-1}$). Previous studies showed that the $V_{ip}$ in ICMEs keeps negligible during propagation from 0.3 to 5 au, but in solar wind it increases with the decreasing heliocentric distance. Therefore, it might be questionable to infer the absolute abundance of coronal sources through in-situ abundance near 1 au for solar wind. Fortunately, the ion-oxygen (O$^{6+}$) differential speed ($V_{io}$) is negligible for He$^{2+}$, C$^{5+}$, and Fe$^{10+}$ within both ICMEs and solar wind, and previous studies suggested that the $V_{io}$ does not vary significantly with the heliocentric distance. This indicates that various heavy ions always flow at the same bulk speed and their relative abundance (relative to oxygen) near 1 au can represent the coronal abundance for both ICMEs and solar wind. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 14 pages, 8 figures

Evidence for Plasmoid-mediated Magnetic Reconnection during a Small-scale Flare in the Partially Ionized Low Solar Atmosphere

Authors: Guanchong Cheng, Lei Ni, Zehao Tang, Yajie Chen, Yuhao Chen, Jialiang Hu, Jun Lin

Abstract: Magnetic reconnection plays a crucial role in the energy release process for different kinds of solar eruptions and activities. The rapid solar eruption requires a fast reconnection model. Plasmoid instability in the reconnecting current sheets is one of the most acceptable fast reconnection mechanisms for explaining the explosive events in the magnetohydrodynamics (MHD) scale, which is also a pot… ▽ More Magnetic reconnection plays a crucial role in the energy release process for different kinds of solar eruptions and activities. The rapid solar eruption requires a fast reconnection model. Plasmoid instability in the reconnecting current sheets is one of the most acceptable fast reconnection mechanisms for explaining the explosive events in the magnetohydrodynamics (MHD) scale, which is also a potential bridge between the macroscopic MHD reconnection process and microscale dissipations. Plenty of high resolution observations indicate that the plasmoid-like structures exist in the high temperature solar corona, but such evidences are very rare in the lower solar atmosphere with partially ionized plasmas. Utilizing joint observations from the Goode Solar Telescope (GST) and the Solar Dynamics Observatory (SDO), we discovered a small-scale eruptive phenomenon in NOAA AR 13085, characterized by clear reconnection cusp structures, supported by Nonlinear Force-Free Field (NLFFF) extrapolation results. The plasmoid-like structures with a size about 150 km were observed to be ejected downward from the current sheet at a maximum velocity of 24 km$\cdot$s$^{-1}$ in the H$α$ line wing images, followed by enhanced emissions at around the post flare loop region in multiple wave lengths. Our 2.5D high-resolution MHD simulations further reproduced such a phenomenon and revealed reconnection fine structures. These results provide comprehensive evidences for the plasmoid mediated reconnection in partially ionized plasmas, and suggest an unified reconnection model for solar flares with different length scales from the lower chromosphere to corona. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Illuminating Black Hole Shadow with Dark Matter Annihilation

Authors: Yifan Chen, Ran Ding, Yuxin Liu, Yosuke Mizuno, Jing Shu, Haiyue Yu, Yanjie Zeng

Abstract: The Event Horizon Telescope (EHT) has revolutionized our ability to study black holes by providing unprecedented spatial resolution and unveiling horizon-scale details. With advancements leading to the next-generation EHT, there is potential to probe even deeper into the black hole's dark region, especially the inner shadow characterized by low-intensity foreground emissions from the jet, thanks t… ▽ More The Event Horizon Telescope (EHT) has revolutionized our ability to study black holes by providing unprecedented spatial resolution and unveiling horizon-scale details. With advancements leading to the next-generation EHT, there is potential to probe even deeper into the black hole's dark region, especially the inner shadow characterized by low-intensity foreground emissions from the jet, thanks to a significant enhancement in dynamic range by two orders of magnitude. We demonstrate how such enhanced observations could transform supermassive black holes into powerful probes for detecting annihilating dark matter, which can form a dense profile in the vicinity of supermassive black holes, by examining the morphology of the black hole image. △ Less

Submitted 2 May, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: 19 pages, 6 figures, references updated

Soft X-ray prompt emission from a high-redshift gamma-ray burst EP240315a

Authors: Y. Liu, H. Sun, D. Xu, D. S. Svinkin, J. Delaunay, N. R. Tanvir, H. Gao, C. Zhang, Y. Chen, X. -F. Wu, B. Zhang, W. Yuan, J. An, G. Bruni, D. D. Frederiks, G. Ghirlanda, J. -W. Hu, A. Li, C. -K. Li, J. -D. Li, D. B. Malesani, L. Piro, G. Raman, R. Ricci, E. Troja , et al. (170 additional authors not shown)

Abstract: Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a,… ▽ More Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 41 pages, 8 figures, 7 tables

New Timing Results of MSPs from NICER Observations

Authors: Shijie Zheng, Dawei Han, Heng Xu, Kejia Lee, Jianping Yuan, Haoxi Wang, Mingyu Ge, Liang Zhang, Yongye Li, Yitao Yin, Xiang Ma, Yong Chen, Shuangnan Zhang

Abstract: Millisecond pulsars (MSPs) are known for their long-term stability. Using six years of observations from the Neutron Star Interior Composition Explorer (NICER), we have conducted an in-depth analysis of the X-ray timing results for six MSPs: PSRs B1937+21, B1821$-$24, J0437$-$4715, J0030+0451, J0218+4232, and J2124$-$3358. The timing stability parameter $σ_z$ has been calculated, revealing remarka… ▽ More Millisecond pulsars (MSPs) are known for their long-term stability. Using six years of observations from the Neutron Star Interior Composition Explorer (NICER), we have conducted an in-depth analysis of the X-ray timing results for six MSPs: PSRs B1937+21, B1821$-$24, J0437$-$4715, J0030+0451, J0218+4232, and J2124$-$3358. The timing stability parameter $σ_z$ has been calculated, revealing remarkable timing precision on the order of $10^{-14}$ for PSRs B1937+21 and J0437$-$4715, and $10^{-13}$ for PSRs B1821$-$24, J0218+4232, and J0030+0451 over a timescale of 1000 days. These findings underscore the feasibility of autonomous in-orbit timekeeping using X-ray observations of MSPs. In addition, the consistency of long-term spin-down noise in the X-ray and radio bands has been investigated by comparison with IPTA radio data. △ Less

Submitted 24 April, 2024; originally announced April 2024.

JWST detections of amorphous and crystalline HDO ice toward massive protostars

Authors: Katerina Slavicinska, Ewine F. van Dishoeck, Łukasz Tychoniec, Pooneh Nazari, Adam E. Rubinstein, Robert Gutermuth, Himanshu Tyagi, Yuan Chen, Nashanty G. C. Brunken, Will R. M. Rocha, P. Manoj, Mayank Narang, S. Thomas Megeath, Yao-Lun Yang, Leslie W. Looney, John J. Tobin, Henrik Beuther, Tyler L. Bourke, Harold Linnartz, Samuel Federman, Dan M. Watson, Hendrik Linz

Abstract: This work aims to utilize the increased sensitivity and resolution of the JWST to quantify the HDO/H$_{2}$O ratio in ices toward young stellar objects (YSOs) and to determine if the HDO/H$_{2}$O ratios measured in the gas phase toward massive YSOs (MYSOs) are representative of the ratios in their ice envelopes. Two protostars observed in the Investigating Protostellar Accretion (IPA) program using… ▽ More This work aims to utilize the increased sensitivity and resolution of the JWST to quantify the HDO/H$_{2}$O ratio in ices toward young stellar objects (YSOs) and to determine if the HDO/H$_{2}$O ratios measured in the gas phase toward massive YSOs (MYSOs) are representative of the ratios in their ice envelopes. Two protostars observed in the Investigating Protostellar Accretion (IPA) program using JWST NIRSpec were analyzed: HOPS 370, an intermediate-mass YSO (IMYSO), and IRAS 20126+4104, a MYSO. The HDO ice toward these sources was detected above the 3$σ$ level and quantified via its 4.1 $μ$m band. The contributions from the CH$_{3}$OH combination modes to the observed optical depth in this spectral region were constrained via the CH$_{3}$OH 3.53 $μ$m band to ensure that the integrated optical depth of the HDO feature was not overestimated. H$_{2}$O ice was quantified via its 3 $μ$m band. From these fits, ice HDO/H$_{2}$O abundance ratios of 4.6$\pm$1.8$\times$10$^{-3}$ and 2.6$\pm$1.2$\times$10$^{-3}$ are obtained for HOPS 370 and IRAS 20126+4104, respectively. The simultaneous detections of both crystalline HDO and crystalline H$_{2}$O corroborate the assignment of the observed feature at 4.1 $μ$m to HDO ice. The ice HDO/H$_{2}$O ratios are similar to the highest reported gas HDO/H$_{2}$O ratios measured toward MYSOs as well as the hot inner regions of isolated low-mass protostars, suggesting that at least some of the gas HDO/H$_{2}$O ratios measured toward massive hot cores are representative of the HDO/H$_{2}$O ratios in ices. The need for an H$_{2}$O-rich CH$_{3}$OH component in the CH$_{3}$OH ice analysis supports recent experimental and observational results that indicate that some CH$_{3}$OH ice may form prior to the CO freeze-out stage in H$_{2}$O-rich ice layers. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A. 23 pages, 17 figures, 10 tables

LIGO operates with quantum noise below the Standard Quantum Limit

Authors: Wenxuan Jia, Victoria Xu, Kevin Kuns, Masayuki Nakano, Lisa Barsotti, Matthew Evans, Nergis Mavalvala, Rich Abbott, Ibrahim Abouelfettouh, Rana Adhikari, Alena Ananyeva, Stephen Appert, Koji Arai, Naoki Aritomi, Stuart Aston, Matthew Ball, Stefan Ballmer, David Barker, Beverly Berger, Joseph Betzwieser, Dripta Bhattacharjee, Garilynn Billingsley, Nina Bode, Edgard Bonilla, Vladimir Bossilkov , et al. (146 additional authors not shown)

Abstract: Precision measurements of space and time, like those made by the detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO), are often confronted with fundamental limitations imposed by quantum mechanics. The Heisenberg uncertainty principle dictates that the position and momentum of an object cannot both be precisely measured, giving rise to an apparent limitation called the Stan… ▽ More Precision measurements of space and time, like those made by the detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO), are often confronted with fundamental limitations imposed by quantum mechanics. The Heisenberg uncertainty principle dictates that the position and momentum of an object cannot both be precisely measured, giving rise to an apparent limitation called the Standard Quantum Limit (SQL). Reducing quantum noise below the SQL in gravitational-wave detectors, where photons are used to continuously measure the positions of freely falling mirrors, has been an active area of research for decades. Here we show how the LIGO A+ upgrade reduced the detectors' quantum noise below the SQL by up to 3 dB while achieving a broadband sensitivity improvement, more than two decades after this possibility was first presented. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Report number: LIGO-P2400059

APPLE: An Evolution Code for Modeling Giant Planets

Authors: Ankan Sur, Yubo Su, Roberto Tejada Arevalo, Yi-Xian Chen, Adam Burrows

Abstract: We introduce APPLE, a novel planetary evolution code designed specifically for the study of giant exoplanet and Jovian planet evolution in the era of Galileo, Juno, and Cassini. With APPLE, state-of-the-art equations of state for hydrogen, helium, ice, and rock are integrated with advanced features to treat ice/rock cores and metals in the gaseous envelope; models for helium rain and hydrogen/heli… ▽ More We introduce APPLE, a novel planetary evolution code designed specifically for the study of giant exoplanet and Jovian planet evolution in the era of Galileo, Juno, and Cassini. With APPLE, state-of-the-art equations of state for hydrogen, helium, ice, and rock are integrated with advanced features to treat ice/rock cores and metals in the gaseous envelope; models for helium rain and hydrogen/helium immiscibility; detailed atmosphere boundary tables that also provide self-consistent albedos and spectra; and options to address envelope metal gradients and stably-stratified regions. Our hope is that these purpose-built features of APPLE will help catalyze the development of the next generation of giant exoplanet and Jovian planet evolutionary models. △ Less

Submitted 11 June, 2024; v1 submitted 22 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ

JWST ERS Program Q3D: The pitfalls of virial BH mass constraints shown in a z = 3 quasar with an ultramassive host

Authors: Caroline Bertemes, Dominika Wylezalek, David S. N. Rupke, Nadia L. Zakamska, Sylvain Veilleux, Benjamin Beckmann, Andrey Vayner, Swetha Sankar, Yuzo Ishikawa, Nadiia Diachenko, Weizhe Liu, Yu-Ching Chen, Jerome Seebeck, Dieter Lutz, Guilin Liu

Abstract: We present JWST MIRI/NIRSpec observations of the extremely red quasar SDSS J165202.64+172852.3 at z~3, one of the most luminous quasars known to date, driving powerful outflows and hosting a clumpy starburst, amidst several interacting companions. We estimate the black hole (BH) mass of the system based on the broad H$α$ and H$β$ lines, as well as the Pa$β$ emission in the IR and MgII in the UV. W… ▽ More We present JWST MIRI/NIRSpec observations of the extremely red quasar SDSS J165202.64+172852.3 at z~3, one of the most luminous quasars known to date, driving powerful outflows and hosting a clumpy starburst, amidst several interacting companions. We estimate the black hole (BH) mass of the system based on the broad H$α$ and H$β$ lines, as well as the Pa$β$ emission in the IR and MgII in the UV. We recover a very broad range of mass estimates, with constraints ranging between log $M_{\rm BH}$=9 and 10.1, which is exacerbated if imposing a uniform BLR geometry at all wavelengths. Several factors may contribute to the large spread: measurement uncertainties (insufficient sensitivity to detect the broadest component of the faint Pa$β$ line, spectral blending, ambiguities in the broad/narrow component distinction), lack of virial equilibrium (in a system characterised by powerful outflows and rapid accretion), and uncertainties on the luminosity-inferred size of the broad line region, a.o. given central dust obscuration. We constrain the stellar mass via SED fitting, suggesting the host to be extremely massive at $10^{12.8\pm 0.5} M_\odot$ - ~2 dex above the characteristic mass of the Schechter fit to the z=3 stellar mass function. Notably, J1652's central BH might be interpreted as being either undermassive, overmassive, or in line with the BH mass-stellar mass relation, depending on the choice of assumptions. The recovered Eddington ratio varies accordingly, but exceeds 10% in any case. We put our results into context by providing an extensive overview and discussion of recent literature results and their associated assumptions. Our findings provide an important demonstration of the uncertainties inherent in virial BH mass estimates, which are of particular relevance in the JWST era given the growing number of studies on rapidly accreting quasars at high redshift. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 15 pages, 7 figures, submitted to A&A, comments welcome

MESSENGER observations of Mercury's planetary ion escape rates and their dependence on true anomaly angle

Authors: Weijie Sun, Ryan M. Dewey, Xianzhe Jia, Jim M. Raines, James A. Slavin, Yuxi Chen, Tai Phan, Gangkai Poh, Shaosui Xu, Anna Milillo, Robert Lillis, Yoshifumi Saito, Stefano Livi, Stefano Orsini

Abstract: This study investigates the escape of Mercury's sodium-group ions (Na+-group, including ions with m/q from 21 to 30 amu/e) and their dependence on true anomaly angle (TAA), i.e., Mercury's orbital phase around the Sun, using measurements from MESSENGER. The measurements are categorized into solar wind, magnetosheath, and magnetosphere, and further divided into four TAA intervals. Na+-group ions fo… ▽ More This study investigates the escape of Mercury's sodium-group ions (Na+-group, including ions with m/q from 21 to 30 amu/e) and their dependence on true anomaly angle (TAA), i.e., Mercury's orbital phase around the Sun, using measurements from MESSENGER. The measurements are categorized into solar wind, magnetosheath, and magnetosphere, and further divided into four TAA intervals. Na+-group ions form escape plumes in the solar wind and magnetosheath, with higher fluxes along the solar wind's motional electric field. The total escape rates vary from 0.2 to 1 times 10^{25} atoms/s with the magnetosheath being the main escaping region. These rates exhibit a TAA dependence, peaking near the perihelion and similar during Mercury's remaining orbit. Despite Mercury's tenuous exosphere, Na+-group ions escape rate is comparable to other inner planets. This can be attributed to several processes, including that Na+-group ions may include several ion species, efficient photoionization frequency for elements within Na+-group, etc. △ Less

Submitted 20 April, 2024; originally announced April 2024.

Comments: 21 pages, 4 figures

Carbon- and Oxygen-rich stars in MaStar: identification and classification

Authors: Lewis Hill, Claudia Maraston, Daniel Thomas, Renbin Yan, Yanping Chen, Guy S. Stringfellow, Richard R. Lane, José G. Fernández-Trincado

Abstract: Carbon- and Oxygen-rich stars populating the Thermally-Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ~1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (~3 Myr) TP-AGB lifetime. Here we exp… ▽ More Carbon- and Oxygen-rich stars populating the Thermally-Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ~1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (~3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (~60,000 spectra) to identify C,O-rich type stars. We define an optical colour selection with cuts of (g-r)>2 and (g-i)<1.55(g-r)-0.07, calibrated with known C- and O- rich spectra. This identifies C-,O-rich stars along clean, separated sequences. An analogue selection is found in V,R,I bands. Our equation identifies C- and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C and O-rich spectra. We find remarkably good fits for the majority of candidates and categorise them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 16 pages, 13 figures, MNRAS in press

The Population of Massive Stars in AGN Disks

Authors: Yi-Xian Chen, Douglas N. C. Lin

Abstract: Gravitational instability in the outskirts of Active Galactic Nuclei (AGN) disks lead to disk fragmentation and formation of super-massive (several 10^2Msun) stars with potentially long lifetimes. Alternatively, stars can be captured ex-situ and grow from gas accretion in the AGN disk. However, the number density distribution throughout the disk is limited by thermal feedback as their luminosities… ▽ More Gravitational instability in the outskirts of Active Galactic Nuclei (AGN) disks lead to disk fragmentation and formation of super-massive (several 10^2Msun) stars with potentially long lifetimes. Alternatively, stars can be captured ex-situ and grow from gas accretion in the AGN disk. However, the number density distribution throughout the disk is limited by thermal feedback as their luminosities provide the dominant heating source. We derive equilibrium stellar surface density profiles under two limiting contexts: in the case where the stellar lifetimes are prolonged due to recycling of hydrogen rich disk gas, only the fraction of gas converted into heat is removed from the disk accretion flow. Alternatively, if stellar composition recycling is inefficient and stars can evolve off the main sequence, the disk accretion rate is quenched towards smaller radii resembling a classical star-burst disk, albeit the effective removal rate depends not only on the stellar lifetime, but also the mass of stellar remnants. For AGNs with central Supermassive Black Hole (SMBH) masses of \sim 10^6 to 10^8Msun accreting at \sim 0.1 Eddington efficiency, we estimate a total number of 10^3 to 10^5 coexisting massive stars and the rate of stellar mergers to be 10^-3 to 1 per year. We motivate the detailed study of interaction between a swarm of massive stars through hydro and N body simulations to provide better prescriptions of dynamical processes in AGN disks, and to constrain more accurate estimates of the stellar population. △ Less

Submitted 30 May, 2024; v1 submitted 12 April, 2024; originally announced April 2024.

Comments: Published in ApJ

Monitoring AGNs with H$β$ Asymmetry. IV. First Reverberation Mapping Results of 14 AGNs

Authors: T. E. Zastrocky, Michael S. Brotherton, Pu Du, Jacob N. McLane, Kianna A. Olson, D. A. Dale, H. A. Kobulnicky, Jaya Maithil, My L. Nguyen, William T. Chick, David H. Kasper, Derek Hand, C. Adelman, Z. Carter, G. Murphree, M. Oeur, T. Roth, S. Schonsberg, M. J. Caradonna, J. Favro, A. J. Ferguson, I. M. Gonzalez, L. M. Hadding, H. D. Hagler, C. J. Rogers , et al. (19 additional authors not shown)

Abstract: We report first-time reverberation mapping results for 14 AGNs from the ongoing Monitoring AGNs with H$β$ Asymmetry campaign (MAHA). These results utilize optical spectra obtained with the Long Slit Spectrograph on the Wyoming Infrared 2.3m Telescope between 2017 November-2023 May. MAHA combines long-duration monitoring with high cadence. We report results from multiple observing seasons for 9 of… ▽ More We report first-time reverberation mapping results for 14 AGNs from the ongoing Monitoring AGNs with H$β$ Asymmetry campaign (MAHA). These results utilize optical spectra obtained with the Long Slit Spectrograph on the Wyoming Infrared 2.3m Telescope between 2017 November-2023 May. MAHA combines long-duration monitoring with high cadence. We report results from multiple observing seasons for 9 of the 14 objects. These results include H$β$ time lags, supermassive black hole masses, and velocity-resolved time lags. The velocity-resolved lags allow us to investigate the kinematics of the broad-line region. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 35 pages, 19 figures, accepted for publication in ApJ Supplement

Probing the shape of the brown dwarf desert around main-sequence A-F-G-type stars using post-common-envelope WD$-$BD binaries

Authors: Zhangliang Chen, Yizhi Chen, Chen Chen, Hongwei Ge, Bo Ma

Abstract: Brown dwarfs (BDs) possessing masses within the range $40-60 M_{\rm Jup}$ are rare around solar-type main-sequence (MS) stars, which gives rise to the brown dwarf desert (BDD). One caveat associated with previous studies of BDD is the relatively limited sample size of MS$-$BD binaries with accurately determined BD masses. We aim to produce a large sample of brown dwarf companions with precisely de… ▽ More Brown dwarfs (BDs) possessing masses within the range $40-60 M_{\rm Jup}$ are rare around solar-type main-sequence (MS) stars, which gives rise to the brown dwarf desert (BDD). One caveat associated with previous studies of BDD is the relatively limited sample size of MS$-$BD binaries with accurately determined BD masses. We aim to produce a large sample of brown dwarf companions with precisely determined mass around main-sequence A-F-G type stars using observations of post common-envelope white dwarf (WD)$-$BD binaries. We employ the rapid binary evolution code COMPAS to deduce the properties of MS$-$BD binary progenitors from post common-envelope WD$-$BD binaries. This method supplements the directly observed MS$-$BD binary sample, enriching the data available for analyzing BDD around main-sequence A-F-G type stars. Our study opens a new window for studying the shape of BDD around A-F-G type main-sequence stars in the short period regime. We find tentative evidence that the `driest' part of BDD around A-F-G type stars may extend into an orbital period of several hundred days, albeit with a small sample size. More post common-envelope WD$-$BD binaries detected in the future will advance our understanding of the BDD around A-F-G type stars. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 7 pages, 4 figures. Accepted for publication in A&A

The influence of thermonuclear bursts on polar caps of accreting X-ray millisecond pulsar MAXI J1816-195

Authors: Long Ji, Mingyu Ge, Yupeng Chen, Zhaosheng Li, Peng-Ju Wang, Shu Zhang, Shuang-Nan Zhang

Abstract: We report accretion-powered pulsations for the first time during thermonuclear bursts in hard X-rays, which were observed with Insight-HXMT in 2022 during the outburst of the accreting X-ray millisecond pulsar MAXI J1816-195. By stacking 73 bursts, we detected pulse profiles in 8-30 keV and 30-100 keV during bursts, which are identical to those obtained from the persistent (non-burst) emission. On… ▽ More We report accretion-powered pulsations for the first time during thermonuclear bursts in hard X-rays, which were observed with Insight-HXMT in 2022 during the outburst of the accreting X-ray millisecond pulsar MAXI J1816-195. By stacking 73 bursts, we detected pulse profiles in 8-30 keV and 30-100 keV during bursts, which are identical to those obtained from the persistent (non-burst) emission. On average, no significant phase lag was observed between burst and persistent pulse profiles. In addition, we suggest that the interaction with burst photons can be used as a direct diagnostic to distinguish contributions from the hot plasma near polar caps and the corona around the accretion disk, which are highly degenerate in their spectral shapes. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 3 figures, accepted for publication in ApJL

Alfven Wave Mode Conversion in Neutron Star Magnetospheres: A Semi-analytic Approach

Authors: Alexander Y. Chen, Yajie Yuan, Dominic Bernardi

Abstract: We write down the force-free electrodynamics (FFE) equations in dipole coordinates, and solve for normal modes corresponding to Alfvénic perturbations in the magnetosphere of a neutron star. We show that a single Alfvén wave propagating on dipole field lines spontaneously sources a fast magnetosonic (fms) wave at the next order in the perturbation expansion, without needing 3-wave interaction. The… ▽ More We write down the force-free electrodynamics (FFE) equations in dipole coordinates, and solve for normal modes corresponding to Alfvénic perturbations in the magnetosphere of a neutron star. We show that a single Alfvén wave propagating on dipole field lines spontaneously sources a fast magnetosonic (fms) wave at the next order in the perturbation expansion, without needing 3-wave interaction. The frequency of the sourced fms wave is twice the original Alfvén wave frequency, and the wave propagates spherically outwards. The properties of the outgoing fms wave can be computed exactly using the usual devices of classical electrodynamics. We extend the calculation to the closed zone of a rotating neutron star magnetosphere, and show that the Alfvén wave also sources a spherical fms wave but at the same frequency as the primary Alfvén wave. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: 15 pages, 5 figures, comments welcome!

A Monte-Carlo Simulation on Resonant Scattering of X-ray Line Emission in Supernova Remnants

Authors: Yiping Li, Gao-Yuan Zhang, Yang Chen, Lei Sun, Shuinai Zhang

Abstract: Resonant scattering (RS) of X-ray line emission in supernova remnants (SNRs) may modify the observed line profiles and fluxes and has potential impact on estimating the physical properties of the hot gas and hence on understanding the SNR physics, but has not been theoretically modeled ever. Here we present our Monte-Carlo simulation of RS effect on X-ray resonant-line emission, typified by O VII… ▽ More Resonant scattering (RS) of X-ray line emission in supernova remnants (SNRs) may modify the observed line profiles and fluxes and has potential impact on estimating the physical properties of the hot gas and hence on understanding the SNR physics, but has not been theoretically modeled ever. Here we present our Monte-Carlo simulation of RS effect on X-ray resonant-line emission, typified by O VII He$α$ r line, from SNRs. We employ the physical conditions characterized by the Sedov-Taylor solution and some basic parameters similar to those in Cygnus Loop. We show that the impact of RS effect is most significant near the edge of the remnant. The line profiles are predicted to be asymmetric because of different temperatures and photon production efficiencies of the expanding gas at different radii. We also predict the surface brightness of the line emission would decrease in the outer projected region but is slightly enhanced in the inner. The G-ratio of the OVII He$α$ triplet can be effectively elevated by RS in the outer region. We show that RS effect of the O VII He$α$ r line in the southwestern boundary region of Cygnus Loop is non-negligible. The observed OVII G-ratio $\sim$1.8 of the region could be achieved with RS taken into account for properly elevated O abundance from the previous estimates. Additional simulation performed for the SNRs in ejecta-dominated phase like Cas A shows that RS in the shocked ejecta may have some apparently effects on the observational properties of oxygen resonant lines. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Comments: 19 pages, 17 figures, accepted by ApJ

LHAASO-KM2A detector simulation using Geant4

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (254 additional authors not shown)

Abstract: KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with… ▽ More KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Thermal X-ray Emission in the Western Half of the LMC Superbubble 30 Dor C

Authors: Yi-Heng Chi, Han-Xiao Chen, Yang Chen, Yi-Fan Meng, Ping Zhou, Lei Sun, Wei Sun

Abstract: While 30 Dor C is a unique superbubble in the Large Magellanic Cloud for its luminous non-thermal X-ray emission, the thermal X-ray emission it emanates has not yet been thoroughly investigated and well constrained. Based on the separate ~1 Ms deep XMM-Newton and Chandra observations, we report the discovery of the thermally-emitting plasma in some portions of the western half of 30 Dor C. The the… ▽ More While 30 Dor C is a unique superbubble in the Large Magellanic Cloud for its luminous non-thermal X-ray emission, the thermal X-ray emission it emanates has not yet been thoroughly investigated and well constrained. Based on the separate ~1 Ms deep XMM-Newton and Chandra observations, we report the discovery of the thermally-emitting plasma in some portions of the western half of 30 Dor C. The thermal emission can be reproduced by a collisional-ionization-equilibrium plasma model with an average electron temperature of ~0.4 keV. We find a significant overabundance of the intermediate-mass elements such as O, Ne, Mg, and Si, which may be indicative of a recent supernova explosion in 30 Dor C. Dynamical properties in combination with the information of the OB association LH 90 suggest that the internal post-main-sequence stars dominate the power of the superbubble and blow it out in the past ~1 Myr. △ Less

Submitted 25 June, 2024; v1 submitted 6 April, 2024; originally announced April 2024.

Comments: 15 pages, 18 figures. Accepted by MNRAS

Are High $Σ_1$ Massive Blue Spiral Galaxies Rejuvenated Systems?

Authors: Cai-Na Hao, Xiaoyang Xia, Yong Shi, Rui Guo, Yanmei Chen, Shuai Feng, Junqiang Ge, Qiusheng Gu

Abstract: Quiescent galaxies generally possess denser cores than star-forming galaxies with similar mass. As a measurement of the core density, the central stellar mass surface density within a radius of 1 kpc ($Σ_1$) was thus suggested to be closely related to galaxy quenching. Massive star-forming galaxies with high $Σ_1$ do not fit into this picture. To understand the origin of such galaxies, we compare… ▽ More Quiescent galaxies generally possess denser cores than star-forming galaxies with similar mass. As a measurement of the core density, the central stellar mass surface density within a radius of 1 kpc ($Σ_1$) was thus suggested to be closely related to galaxy quenching. Massive star-forming galaxies with high $Σ_1$ do not fit into this picture. To understand the origin of such galaxies, we compare the spatially-resolved stellar population and star formation properties of massive ($ > 10^{10.5}{\rm M}_{\odot}$) blue spiral galaxies with high and low $Σ_1$, divided by $Σ_1 = 10^{9.4} M_\odot \, {\rm kpc}^{-2}$, based on the final release of MaNGA IFU data. We find that both high $Σ_1$ and low $Σ_1$ blue spirals show large diversities in stellar population and star formation properties. Despite the diversities, high $Σ_1$ blue spirals are statistically different from the low $Σ_1$ ones. Specifically, the radial profiles of the luminosity-weighted age and Mgb/${\rm \langle Fe \rangle}$ show that high $Σ_1$ blue spirals consist of a larger fraction of galaxies with younger and less $α$-element enhanced centers than their low $Σ_1$ counterparts, $\sim 55\%$ versus $\sim 30\%$. The galaxies with younger centers mostly have higher central specific star formation rates, which still follow the spaxel-based star formation main sequence relation though. Examinations of the H$α$ velocity field and the optical structures suggest that galactic bars or galaxy interactions should be responsible for the rejuvenation of these galaxies. The remaining $\sim 45\% $ of high $Σ_1$ blue spirals are consistent with the inside-out growth scenario. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 22 pages, 14 figures, accepted for publication in ApJ

CSST Strong Lensing Preparation: a Framework for Detecting Strong Lenses in the Multi-color Imaging Survey by the China Survey Space Telescope (CSST)

Authors: Xu Li, Ruiqi Sun, Jiameng Lv, Peng Jia, Nan Li, Chengliang Wei, Zou Hu, Xinzhong Er, Yun Chen, Zhang Ban, Yuedong Fang, Qi Guo, Dezi Liu, Guoliang Li, Lin Lin, Ming Li, Ran Li, Xiaobo Li, Yu Luo, Xianmin Meng, Jundan Nie, Zhaoxiang Qi, Yisheng Qiu, Li Shao, Hao Tian , et al. (7 additional authors not shown)

Abstract: Strong gravitational lensing is a powerful tool for investigating dark matter and dark energy properties. With the advent of large-scale sky surveys, we can discover strong lensing systems on an unprecedented scale, which requires efficient tools to extract them from billions of astronomical objects. The existing mainstream lens-finding tools are based on machine learning algorithms and applied to… ▽ More Strong gravitational lensing is a powerful tool for investigating dark matter and dark energy properties. With the advent of large-scale sky surveys, we can discover strong lensing systems on an unprecedented scale, which requires efficient tools to extract them from billions of astronomical objects. The existing mainstream lens-finding tools are based on machine learning algorithms and applied to cut-out-centered galaxies. However, according to the design and survey strategy of optical surveys by CSST, preparing cutouts with multiple bands requires considerable efforts. To overcome these challenges, we have developed a framework based on a hierarchical visual Transformer with a sliding window technique to search for strong lensing systems within entire images. Moreover, given that multi-color images of strong lensing systems can provide insights into their physical characteristics, our framework is specifically crafted to identify strong lensing systems in images with any number of channels. As evaluated using CSST mock data based on an Semi-Analytic Model named CosmoDC2, our framework achieves precision and recall rates of 0.98 and 0.90, respectively. To evaluate the effectiveness of our method in real observations, we have applied it to a subset of images from the DESI Legacy Imaging Surveys and media images from Euclid Early Release Observations. 61 new strong lensing system candidates are discovered by our method. However, we also identified false positives arising primarily from the simplified galaxy morphology assumptions within the simulation. This underscores the practical limitations of our approach while simultaneously highlighting potential avenues for future improvements. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: The paper is accepted by the AJ. The complete code could be downloaded with DOI of: 10.12149/101393. Comments are welcome

Transport of the magnetic flux away from a decaying sunspot via convective motions

Authors: Chenxi Zheng, Thierry Roudier, Brigitte Schmieder, Guiping Ruan, Jean-Marie Malherbe, Yang Liu, Yao Chen, Wenda Cao

Abstract: Aims. The aim of this paper is to consider relationship between the decay of sunspots and convection via the motion of the family of granules and how the diffusion mechanism of magnetic field operates in a decaying sunspot. Methods. We report the decay of a sunspot observed by the 1.6m Goode Solar Telescope (GST) with the TiO Broadband Filter Imager (BFI) and the Near-InfraRed Imaging Spectropolar… ▽ More Aims. The aim of this paper is to consider relationship between the decay of sunspots and convection via the motion of the family of granules and how the diffusion mechanism of magnetic field operates in a decaying sunspot. Methods. We report the decay of a sunspot observed by the 1.6m Goode Solar Telescope (GST) with the TiO Broadband Filter Imager (BFI) and the Near-InfraRed Imaging Spectropolarimeter (NIRIS). The analysis was aided by the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamic Observatory (SDO). In the first step, we followed the decay of the sunspot with HMI data over three days by constructing its evolving area and total magnetic flux. In the second step, the high spatial and temporal resolution of the GST instruments allowed us to analyze the causes of the decay of the sunspot. Afterward, we followed the emergence of granules in the moat region around the sunspot over six hours. The evolution of the trees of fragmenting granules (TFGs) was derived based on their relationship with the horizontal surface flows. Results. We find that the area and total magnetic flux display an exponential decrease over the course of the sunspot decay. We identified 22 moving magnetic features (MMFs) in the moats of pores, which is a signature of sunspot decay through diffusion. We note that the MMFs were constrained to follow the borders of TFGs during their journey away from the sunspot. Conclusions. The TFGs and their development contribute to the diffusion of the magnetic field outside the sunspot. The conclusion of our analysis shows the important role of the TFGs in sunspot decay. Finally, the the family of granules evacuates the magnetic field. △ Less

Submitted 31 March, 2024; originally announced April 2024.

Automated Identification and Segmentation of Hi Sources in CRAFTS Using Deep Learning Method

Authors: Zihao Song, Huaxi Chen, Donghui Quan, Di Li, Yinghui Zheng, Shulei Ni, Yunchuan Chen, Yun Zheng

Abstract: We introduce a machine learning-based method for extracting HI sources from 3D spectral data, and construct a dedicated dataset of HI sources from CRAFTS. Our custom dataset provides comprehensive resources for HI source detection. Utilizing the 3D-Unet segmentation architecture, our method reliably identifies and segments HI sources, achieving notable performance metrics with recall rates reachin… ▽ More We introduce a machine learning-based method for extracting HI sources from 3D spectral data, and construct a dedicated dataset of HI sources from CRAFTS. Our custom dataset provides comprehensive resources for HI source detection. Utilizing the 3D-Unet segmentation architecture, our method reliably identifies and segments HI sources, achieving notable performance metrics with recall rates reaching 91.6% and accuracy levels at 95.7%. These outcomes substantiate the value of our custom dataset and the efficacy of our proposed network in identifying HI source. Our code is publicly available at https://github.com/fishszh/HISF. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: 6 pages, 4 figures