-
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.
-
Filamentary Hierarchies and Superbubbles: Galactic Multiscale MHD Simulations of GMC to Star Cluster Formation
Authors:
Bo Zhao,
Ralph E. Pudritz,
Rachel Pillsworth,
Hector Robinson,
James Wadsley
Abstract:
There is now abundant observational evidence that star formation is a highly dynamical process that connects filament hierarchies and supernova feedback from galaxy scale kpc filaments and superbubbles, to giant molecular clouds (GMCs) on 100 pc scales and star clusters (1 pc). Here we present galactic multi-scale MHD simulations that track the formation of structure from galactic down to sub pc s…
▽ More
There is now abundant observational evidence that star formation is a highly dynamical process that connects filament hierarchies and supernova feedback from galaxy scale kpc filaments and superbubbles, to giant molecular clouds (GMCs) on 100 pc scales and star clusters (1 pc). Here we present galactic multi-scale MHD simulations that track the formation of structure from galactic down to sub pc scales in a magnetized, Milky Way like galaxy undergoing supernova driven feedback processes. We do this by adopting a novel zoom-in technique that follows the evolution of typical 3-kpc sub regions without cutting out the surrounding galactic environment, allowing us to reach 0.28 pc resolution in the individual zoom-in regions. We find a wide range of morphologies and hierarchical structure including superbubbles, turbulence, kpc atomic gas filaments hosting multiple GMC condensations that are often associated with superbubble compression; down to smaller scale filamentary GMCs and star cluster regions within them. Gas accretion and compression ultimately drive filaments over a critical, scale - dependent, line mass leading to gravitational instabilities that produce GMCs and clusters. In quieter regions, galactic shear can produce filamentary GMCs within flattened, rotating disk-like structures on 100 pc scales. Strikingly, our simulations demonstrate the formation of helical magnetic fields associated with the formation of these disk like structures.
△ Less
Submitted 28 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.
-
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.
-
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.
-
Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A
Authors:
The LHAASO Collaboration,
Zhen Cao,
F. Aharonian,
Q. An,
A. 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
, et al. (256 additional authors not shown)
Abstract:
We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at…
▽ More
We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components.
△ Less
Submitted 26 March, 2024; v1 submitted 15 March, 2024;
originally announced March 2024.
-
Exploring the dust grain size and polarization mechanism in the hot and massive Class 0 disk IRAS 16293-2422 B
Authors:
Joaquin Zamponi,
María José Maureira,
Hauyu Baobab Liu,
Bo Zhao,
Dominique Segura-Cox,
Chia-Lin Ko,
Paola Caselli
Abstract:
Multiwavelength dust continuum and polarization observations arising from self-scattering have been used to investigate grain sizes in young disks. However, the polarization by self-scattering is low in face-on optically thick disks and puts some of the size constraints from polarization on hold, particularly for the younger and more massive disks. The 1.3 mm emission detected toward the hot (…
▽ More
Multiwavelength dust continuum and polarization observations arising from self-scattering have been used to investigate grain sizes in young disks. However, the polarization by self-scattering is low in face-on optically thick disks and puts some of the size constraints from polarization on hold, particularly for the younger and more massive disks. The 1.3 mm emission detected toward the hot ($\gtrsim$400 K) Class 0 disk IRAS 16293-2422 B has been attributed to self-scattering, predicting grain sizes between 200-2000 $μ$m. We investigate the effects of grain size in the resultant flux and polarization fractions from self-scattering using a hot and massive Class 0 disk model and compare with observations. We compared new and archival high-resolution observations between 1.3 and 18 mm to a set of synthetic models. We have developed a new public tool to automate this process called Synthesizer. This is an easy-to-use program to generate synthetic observations from numerical simulations. Optical depths are in the range of 130 to 2 from 1.3 to 18 mm, respectively. Predictions from significant grain growth populations, including millimetric grains are comparable to the observations at all wavelengths. The polarization fraction produced by self-scattering reaches a maximum of $\sim$0.1% at 1.3 mm for a maximum grain size of 100 $μ$m, being an order of magnitude lower than that observed with ALMA. From the comparison of Stokes I fluxes, we conclude that significant grain growth could be present in the young Class 0 disk IRAS 16293 B, particularly in the inner hot region ($<10$ au, $T>$ 300 K) where refractory organics evaporate. The polarization produced by self-scattering in our model is not high enough to explain the observations at 1.3 and 7 mm, and effects like dichroic extinction or polarization reversal of elongated aligned grains remain other possible but untested scenarios.
△ Less
Submitted 4 November, 2023;
originally announced November 2023.
-
Does or did the supernova remnant Cassiopeia A operate as a PeVatron?
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:
For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE;…
▽ More
For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_γ\geq 100$~TeV) $γ$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
Very high energy gamma-ray emission beyond 10 TeV from GRB 221009A
Authors:
Zhen Cao,
F. Aharonian,
Q. An,
A. 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 highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the t…
▽ More
The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the trigger. The intrinsic energy spectrum of gamma-rays can be described by a power-law after correcting for extragalactic background light (EBL) absorption. Such a hard spectrum challenges the synchrotron self-Compton (SSC) scenario of relativistic electrons for the afterglow emission above several TeV. Observations of gamma-rays up to 13 TeV from a source with a measured redshift of z=0.151 hints more transparency in intergalactic space than previously expected. Alternatively, one may invoke new physics such as Lorentz Invariance Violation (LIV) or an axion origin of very high energy (VHE) signals.
△ Less
Submitted 22 November, 2023; v1 submitted 13 October, 2023;
originally announced October 2023.
-
Single View Refractive Index Tomography with Neural Fields
Authors:
Brandon Zhao,
Aviad Levis,
Liam Connor,
Pratul P. Srinivasan,
Katherine L. Bouman
Abstract:
Refractive Index Tomography is the inverse problem of reconstructing the continuously-varying 3D refractive index in a scene using 2D projected image measurements. Although a purely refractive field is not directly visible, it bends light rays as they travel through space, thus providing a signal for reconstruction. The effects of such fields appear in many scientific computer vision settings, ran…
▽ More
Refractive Index Tomography is the inverse problem of reconstructing the continuously-varying 3D refractive index in a scene using 2D projected image measurements. Although a purely refractive field is not directly visible, it bends light rays as they travel through space, thus providing a signal for reconstruction. The effects of such fields appear in many scientific computer vision settings, ranging from refraction due to transparent cells in microscopy to the lensing of distant galaxies caused by dark matter in astrophysics. Reconstructing these fields is particularly difficult due to the complex nonlinear effects of the refractive field on observed images. Furthermore, while standard 3D reconstruction and tomography settings typically have access to observations of the scene from many viewpoints, many refractive index tomography problem settings only have access to images observed from a single viewpoint. We introduce a method that leverages prior knowledge of light sources scattered throughout the refractive medium to help disambiguate the single-view refractive index tomography problem. We differentiably trace curved rays through a neural field representation of the refractive field, and optimize its parameters to best reproduce the observed image. We demonstrate the efficacy of our approach by reconstructing simulated refractive fields, analyze the effects of light source distribution on the recovered field, and test our method on a simulated dark matter mapping problem where we successfully recover the 3D refractive field caused by a realistic dark matter distribution.
△ Less
Submitted 1 December, 2023; v1 submitted 8 September, 2023;
originally announced September 2023.
-
On the importance of disc chemistry in the formation of protoplanetary disc rings
Authors:
C. A. Nolan,
B. Zhao,
P. Caselli,
Z. Y. Li
Abstract:
Radial substructures have now been observed in a wide range of protoplanetary discs (PPDs), from young to old systems, however their formation is still an area of vigorous debate. Recent magnetohydrodynamic (MHD) simulations have shown that rings and gaps can form naturally in PPDs when non-ideal MHD effects are included. However these simulations employ ad-hoc approximations to the magnitudes of…
▽ More
Radial substructures have now been observed in a wide range of protoplanetary discs (PPDs), from young to old systems, however their formation is still an area of vigorous debate. Recent magnetohydrodynamic (MHD) simulations have shown that rings and gaps can form naturally in PPDs when non-ideal MHD effects are included. However these simulations employ ad-hoc approximations to the magnitudes of the magnetic diffusivities in order to facilitate ring growth. We replace the parametrisation of these terms with a simple chemical network and grain distribution model to calculate the non-ideal effects in a more self-consistent way. We use a range of grain distributions to simulate grain formation for different disc conditions. Including ambipolar diffusion, we find that large grain populations (> 1μm), and those including a population of very small polyaromatic hydrocarbons (PAHs) facilitate the growth of periodic, stable rings, while intermediate sized grains suppress ring formation. Including Ohmic diffusion removes the positive influence of PAHs, with only large grain populations still producing periodic ring and gap structures. These results relate closely to the degree of coupling between the magnetic field and the neutral disc material, quantified by the non-dimensional Elsasser number Λ (the ratio of magnetic forces to Coriolis force). For both the ambipolar-only and ambipolar-ohmic cases, if the total Elsasser number is initially of order unity along the disc mid-plane, ring and gap structures may develop.
△ Less
Submitted 24 August, 2023;
originally announced August 2023.
-
Long-term multiwavelength monitoring and reverberation mapping of NGC 2617 during a changing-look event
Authors:
V. L. Oknyansky,
M. S. Brotherton,
S. S. Tsygankov,
A. V. Dodin,
A. M. Tatarnikov,
P. Du,
D. -W. Bao,
M. A. Burlak,
N. P. Ikonnikova,
V. M. Lipunov,
E. S. Gorbovskoy,
V. G. Metlov,
A. A. Belinski,
N. I. Shatsky,
S. G. Zheltouhov,
N. A. Maslennikova,
J. -M. Wang,
S. Zhai,
F. -N. Fang,
Y. -X. Fu,
H. -R. Bai,
D. Kasper,
N. A. Huseynov,
J. N. McLane,
J. Maithil
, et al. (10 additional authors not shown)
Abstract:
We present the results of photometric and spectroscopic monitoring campaigns of the changing look AGN NGC~2617 carried out from 2016 until 2022 and covering the wavelength range from the X-ray to the near-IR. The facilities included the telescopes of the SAI MSU, MASTER Global Robotic Net, the 2.3-m WIRO telescope, Swift, and others. We found significant variability at all wavelengths and, specifi…
▽ More
We present the results of photometric and spectroscopic monitoring campaigns of the changing look AGN NGC~2617 carried out from 2016 until 2022 and covering the wavelength range from the X-ray to the near-IR. The facilities included the telescopes of the SAI MSU, MASTER Global Robotic Net, the 2.3-m WIRO telescope, Swift, and others. We found significant variability at all wavelengths and, specifically, in the intensities and profiles of the broad Balmer lines. We measured time delays of ~ 6 days (~ 8 days) in the responses of the H-beta (H-alpha) line to continuum variations. We found the X-ray variations to correlate well with the UV and optical (with a small time delay of a few days for longer wavelengths). The K-band lagged the B band by 14 +- 4 days during the last 3 seasons, which is significantly shorter than the delays reported previously by the 2016 and 2017--2019 campaigns. Near-IR variability arises from two different emission regions: the outer part of the accretion disc and a more distant dust component. The HK-band variability is governed primarily by dust. The Balmer decrement of the broad-line components is inversely correlated with the UV flux. The change of the object's type, from Sy1 to Sy1.8, was recorded over a period of ~ 8 years. We interpret these changes as a combination of two factors: changes in the accretion rate and dust recovery along the line of sight.
△ Less
Submitted 23 August, 2023; v1 submitted 9 August, 2023;
originally announced August 2023.
-
The First LHAASO Catalog of Gamma-Ray Sources
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:
We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022.…
▽ More
We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with an extended size smaller than $2^\circ$ and a significance of detection at $> 5σ$. Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4σ$ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog.
△ Less
Submitted 27 November, 2023; v1 submitted 26 May, 2023;
originally announced May 2023.
-
Measurement of ultra-high-energy diffuse gamma-ray emission of the Galactic plane from 10 TeV to 1 PeV with LHAASO-KM2A
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 diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer ar…
▽ More
The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner ($15^{\circ}<l<125^{\circ}$, $|b|<5^{\circ}$) and outer ($125^{\circ}<l<235^{\circ}$, $|b|<5^{\circ}$) Galactic plane are detected with $29.1σ$ and $12.7σ$ significance, respectively. The outer Galactic plane diffuse emission is detected for the first time in the very- to ultra-high-energy domain ($E>10$~TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of $-2.99\pm0.04$, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of $\sim3$ than the prediction. A similar spectrum with an index of $-2.99\pm0.07$ is found in the outer Galaxy region, and the absolute flux for $10\lesssim E\lesssim60$ TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations.
△ Less
Submitted 19 August, 2023; v1 submitted 9 May, 2023;
originally announced May 2023.
-
Models of Rotating Infall for the B335 Protostar
Authors:
Neal J. Evans II,
Yao-Lun Yang,
Joel D. Green,
Bo Zhao,
James Di Francesco,
Jeong-Eun Lee,
Jes K. Jørgensen,
Minho Choi,
Philip C. Myers,
Diego Mardones
Abstract:
Models of the protostellar source, B335, are developed using axisymmetric three-dimensional models to resolve conflicts found in one-dimensional models. The models are constrained by a large number of observations, including ALMA, Herschel, and Spitzer data. Observations of the protostellar source B335 with ALMA show red-shifted absorption against a central continuum source indicative of infall in…
▽ More
Models of the protostellar source, B335, are developed using axisymmetric three-dimensional models to resolve conflicts found in one-dimensional models. The models are constrained by a large number of observations, including ALMA, Herschel, and Spitzer data. Observations of the protostellar source B335 with ALMA show red-shifted absorption against a central continuum source indicative of infall in the HCO$^+$ and HCN $J = 4\rightarrow 3$ transitions. The data are combined with a new estimate of the distance to provide strong constraints to three-dimensional radiative transfer models including a rotating, infalling envelope, outflow cavities, and a very small disk. The models favor ages since the initiation of collapse between $3 \times 10^4$ and $4 \times 10^4$ yr for both the continuum and the lines, resolving a conflict found in one-dimensional models. The models under-predict the continuum emission seen by ALMA, suggesting an additional component such as a pseudo-disk. The best-fitting model is used to convert variations in the 4.5 $μm$ flux in recent years into a model for a variation of a factor of 5-7 in luminosity over the last 8 years.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Increasing mass-to-flux ratio from the dense core to the protostellar envelope around the Class 0 protostar HH 211
Authors:
Hsi-Wei Yen,
Patrick Koch,
Chin-Fei Lee,
Naomi Hirano,
Nagayoshi Ohashi,
Jinshi Sai,
Shigehisa Takakuwa,
Ya-Wen Tang,
Ken'ichi Tatematsu,
Bo Zhao
Abstract:
To study transportation of magnetic flux from large to small scales in protostellar sources, we analyzed the Nobeyama 45-m N2H+ (1-0), JCMT 850 um polarization, and ALMA C18O (2-1) and 1.3 mm and 0.8 mm (polarized) continuum data of the Class 0 protostar HH 211. The magnetic field strength in the dense core on a 0.1 pc scale was estimated with the single-dish line and polarization data using the D…
▽ More
To study transportation of magnetic flux from large to small scales in protostellar sources, we analyzed the Nobeyama 45-m N2H+ (1-0), JCMT 850 um polarization, and ALMA C18O (2-1) and 1.3 mm and 0.8 mm (polarized) continuum data of the Class 0 protostar HH 211. The magnetic field strength in the dense core on a 0.1 pc scale was estimated with the single-dish line and polarization data using the Davis-Chandrasekhar-Fermi method, and that in the protostellar envelope on a 600 au scale was estimated from the force balance between the gravity and magnetic field tension by analyzing the gas kinematics and magnetic field structures with the ALMA data. Our analysis suggests that from 0.1 pc to 600 au scales, the magnetic field strength increases from 40-107 uG to 0.3-1.2 mG with a scaling relation between the magnetic field strength and density of $B \propto ρ^{0.36\pm0.08}$, and the mass-to-flux ratio increases from 1.2-3.7 to 9.1-32.3. The increase in the mass-to-flux ratio could suggest that the magnetic field is partially decoupled from the neutral matter between 0.1 pc and 600 au scales, and hint at efficient ambipolar diffusion in the infalling protostellar envelope in HH 211, which is the dominant non-ideal magnetohydrodynamic effect considering the density on these scales. Thus, our results could support the scenario of efficient ambipolar diffusion enabling the formation of the 20 au Keplerian disk in HH 211.
△ Less
Submitted 20 November, 2022;
originally announced November 2022.
-
PRODIGE -- Envelope to Disk with NOEMA II. Small-scale temperature structure and a streamer feeding the SVS13A protobinary using CH3CN and DCN
Authors:
T. -H. Hsieh,
D. M. Segura-Cox,
J. E. Pineda,
P. Caselli,
L. Bouscasse,
R. Neri,
A. Lopez-Sepulcre,
M. T. Valdivia-Mena,
M. J. Maureira,
Th. Henning,
G. V. Smirnov-Pinchukov,
D. Semenov,
Th. Möller,
N. Cunningham,
A. Fuente,
S. Marino,
A. Dutrey,
M. Tafalla,
E. Chapillon,
C. Ceccarelli,
B. Zhao
Abstract:
Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accreti…
▽ More
Aims. We present high sensitivity and high-spectral resolution NOEMA observations of the Class 0/I binary system SVS13A, composed of the low-mass protostars VLA4A and VLA4B with a separation of ~90 au. VLA4A is undergoing an accretion burst that enriches the chemistry of the surrounding gas. This gives us an excellent opportunity to probe the chemical and physical conditions as well as the accretion process. Methods. We observe the (12K-11K) lines of CH3CN and CH313CN, the DCN (3-2) line, and the C18O (2-1) line toward SVS13A using NOEMA. Results. We find complex line profiles at disk scales which cannot be explained by a single component or pure Keplerian motion. By adopting two velocity components to model the complex line profiles, we find that the temperatures and densities are significantly different between these two components. This suggests that the physical conditions of the emitting gas traced via CH3CN can change dramatically within the circumbinary disk. In addition, combining our observations of DCN (3-2) with previous ALMA high-angular-resolution observations, we find that the binary system (or VLA4A) might be fed by an infalling streamer from envelope scales (~700 au). If this is the case, this streamer contributes to the accretion of material onto the system with a rate of at least 1.4x10-6 Msun yr-1. Conclusions. We conclude that the CH3CN emission in SVS13A traces hot gas from a complex structure. This complexity might be affected by a streamer that is possibly infalling and funneling material into the central region.
△ Less
Submitted 25 December, 2022; v1 submitted 9 November, 2022;
originally announced November 2022.
-
The role of magnetic fields in the formation of protostars, disks, and outflows
Authors:
Yusuke Tsukamoto,
Anaëlle Maury,
Benoît Commerçon,
Felipe O. Alves,
Erin G. Cox,
Nami Sakai,
Tom Ray,
Bo Zhao,
Masahiro N. Machida
Abstract:
We present our current understanding of the formation and early evolution of protostars, protoplanetary disks, and the driving of outflows as dictated by the interplay of magnetic fields and partially ionized gas in molecular cloud cores. In recent years, the field has witnessed enormous development through sub-millimeter observations which in turn have constrained models of protostar formation. A…
▽ More
We present our current understanding of the formation and early evolution of protostars, protoplanetary disks, and the driving of outflows as dictated by the interplay of magnetic fields and partially ionized gas in molecular cloud cores. In recent years, the field has witnessed enormous development through sub-millimeter observations which in turn have constrained models of protostar formation. As a result of these observations % that the observations provided, the state-of-the-art theoretical understanding of the formation and evolution of young stellar objects is described. In particular, we emphasize the importance of the coupling, decoupling, and re-coupling between weakly ionized gas and the magnetic field on appropriate scales. This highlights the complex and intimate relationship between gravitational collapse and magnetic fields in young protostars.
△ Less
Submitted 27 September, 2022;
originally announced September 2022.
-
PRODIGE -- Envelope to disk with NOEMA I. A 3000 au streamer feeding a Class I protostar
Authors:
M. T. Valdivia-Mena,
J. E. Pineda,
D. M. Segura-Cox,
P. Caselli,
R. Neri,
A. López-Sepulcre,
N. Cunningham,
L. Bouscasse,
D. Semenov,
Th. Henning,
V. Piétu,
E. Chapillon,
A. Dutrey,
A. Fuente,
S. Guilloteau,
T. H. Hsieh,
I. Jiménez-Serra,
S. Marino,
M. J. Maureira,
G. V. Smirnov-Pinchukov,
M. Tafalla,
B. Zhao
Abstract:
Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-f…
▽ More
Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-forming region NGC 1333. Our goal is to study how the mass infall proceeds from envelope to disk scales in this source. Results. We discover a streamer delivering material toward Per-emb-50 in H$_2$CO and C$^{18}$O emission. The streamer's emission can be well described by the analytic solutions for an infalling parcel of gas along a streamline with conserved angular momentum, both in the image plane and along the line of sight velocities. The streamer has a mean infall rate of $1.3 \times 10^{ -6}$ M$_{ \odot}$ yr$^{ -1}$, $5 -10$ times higher than the current accretion rate of the protostar. SO and SO$_2$ emission reveal asymmetric infall motions in the inner envelope, additional to the streamer around Per-emb-50. Furthermore, the presence of SO$_2$ could mark the impact zone of the infalling material. Conclusions. The streamer delivers sufficient mass to sustain the protostellar accretion rate and might produce an accretion burst, which would explain the protostar's high luminosity with respect to other Class I sources. Our results highlight the importance of late infall for protostellar evolution: streamers might provide a significant amount of mass for stellar accretion after the Class 0 phase.
△ Less
Submitted 1 August, 2022;
originally announced August 2022.
-
Flux Variations of Cosmic Ray Air Showers Detected by LHAASO-KM2A During a Thunderstorm on 10 June 2021
Authors:
LHAASO Collaboration,
F. Aharonian,
Q. An,
Axikegu,
L. X. Bai,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
J. T. Cai,
Zhe Cao,
Zhen Cao,
J. Chang,
J. F. Chang,
E. S. Chen,
Liang Chen,
Liang Chen,
Long Chen,
M. J. Chen,
M. L. Chen,
S. H. Chen,
S. Z. Chen,
T. L. Chen,
X. J. Chen
, et al. (248 additional authors not shown)
Abstract:
The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations…
▽ More
The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations of trigger rates (increases or decreases) are found to be strongly dependent on the primary zenith angle. The flux of secondary particles increases significantly, following a similar trend with that of the shower events. To better understand the observed behavior, Monte Carlo simulations are performed with CORSIKA and G4KM2A (a code based on GEANT4). We find that the experimental data (in saturated negative fields) are in good agreement with simulations, assuming the presence of a uniform upward electric field of 700 V/cm with a thickness of 1500 m in the atmosphere above the observation level. Due to the acceleration/deceleration and deflection by the atmospheric electric field, the number of secondary particles with energy above the detector threshold is modified, resulting in the changes in shower detection rate.
△ Less
Submitted 6 December, 2022; v1 submitted 25 July, 2022;
originally announced July 2022.
-
Monitoring AGNs with H$β$ Asymmetry. III. Long-term Reverberation Mapping Results of 15 Palomar-Green Quasars
Authors:
Dong-Wei Bao,
Michael S. Brotherton,
Pu Du,
Jacob N. McLane,
T. E. Zastrocky,
Kianna A. Olson,
Feng-Na Fang,
Shuo Zhai,
Zheng-Peng Huang,
Kai Wang,
Bi-Xuan Zhao,
Sha-Sha Li,
Sen Yang,
Yong-Jie Chen,
Jun-Rong Liu,
Zhu-Heng Yao,
Yue-Chang Peng,
Wei-Jian Guo,
Yu-Yang Songsheng,
Yan-Rong Li,
Bo-Wei Jiang,
David H. Kasper,
William T. Chick,
My L. Nguyen,
Jaya Maithil
, et al. (20 additional authors not shown)
Abstract:
In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire l…
▽ More
In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire light curves and the measurements from individual observing seasons are provided. Reverberation mapping of 9 of our targets has been attempted for the first time, while the results for 6 others can be compared with previous campaigns. We measure the H$β$ time lags over periods of years and estimate their black hole masses. The long duration of the campaign enables us to investigate their broad line region (BLR) geometry and kinematics for different years by using velocity-resolved lags, which demonstrate signatures of diverse BLR geometry and kinematics. The BLR geometry and kinematics of individual objects are discussed. In this sample, the BLR kinematics of Keplerian/virialized motion and inflow is more common than outflow.
△ Less
Submitted 1 July, 2022;
originally announced July 2022.
-
Hot methanol in the [BHB2007] 11 protobinary system: hot corino versus shock origin? : FAUST V
Authors:
C. Vastel,
F. Alves,
C. Ceccarelli,
M. Bouvier,
I. Jimenez-Serra,
T. Sakai,
P. Caselli,
L. Evans,
F. Fontani,
R. Le Gal,
C. J. Chandler,
B. Svoboda,
L. Maud,
C. Codella,
N. Sakai,
A. Lopez-Sepulcre,
G. Moellenbrock,
Y. Aikawa,
N. Balucani,
E. Bianchi,
G. Busquet,
E. Caux,
S. Charnley,
N. Cuello,
M. De Simone
, et al. (41 additional authors not shown)
Abstract:
Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program F…
▽ More
Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program FAUST (Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars) to study the methanol line emission towards the [BHB2007] 11 protobinary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has been previously detected. Twelve methanol lines have been detected with upper energies in the range [45-537] K along with one 13CH3OH transition. The methanol emission is compact and encompasses both protostars, separated by only 28 au and presents three velocity components, not spatially resolved by our observations, associated with three different spatial regions, with two of them close to 11B and the third one associated with 11A. A non-LTE radiative transfer analysis of the methanol lines concludes that the gas is hot and dense and highly enriched in methanol with an abundance as high as 1e-5. Using previous continuum data, we show that dust opacity can potentially completely absorb the methanol line emission from the two binary objects. Although we cannot firmly exclude other possibilities, we suggest that the detected hot methanol is resulting from the shocked gas from the incoming filaments streaming towards [BHB2007] 11 A and B, respectively. Higher spatial resolution observations are necessary to confirm this hypothesis.
△ Less
Submitted 21 June, 2022;
originally announced June 2022.
-
The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out
Authors:
Paola Caselli,
Jaime E. Pineda,
Olli Sipilä,
Bo Zhao,
Elena Redaelli,
Silvia Spezzano,
Maria José Maureira,
Felipe Alves,
Luca Bizzocchi,
Tyler L. Bourke,
Ana Chacón-Tanarro,
Rachel Friesen,
Daniele Galli,
Jorma Harju,
Izaskun Jiménez-Serra,
Eric Keto,
Zhi-Yun Li,
Marco Padovani,
Anika Schmiedeke,
Mario Tafalla,
Charlotte Vastel
Abstract:
Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particula…
▽ More
Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved.
△ Less
Submitted 27 February, 2022;
originally announced February 2022.
-
FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417
Authors:
Satoshi Ohashi,
Claudio Codella,
Nami Sakai,
Claire J. Chandler,
Cecilia Ceccarelli,
Felipe Alves,
Davide Fedele,
Tomoyuki Hanawa,
Aurora Durán,
Cécile Favre,
Ana López-Sepulcre,
Laurent Loinard,
Seyma Mercimek,
Nadia M. Murillo,
Linda Podio,
Yichen Zhang,
Yuri Aikawa,
Nadia Balucani,
Eleonora Bianchi,
Mathilde Bouvier,
Gemma Busquet,
Paola Caselli,
Emmanuel Caux,
Steven Charnley,
Spandan Choudhury
, et al. (47 additional authors not shown)
Abstract:
We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the…
▽ More
We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed.
△ Less
Submitted 18 January, 2022;
originally announced January 2022.
-
Peta-electron volt gamma-ray emission from the Crab Nebula
Authors:
The LHAASO Collaboration,
Zhen Cao,
F. Aharonian,
Q. An,
Axikegu,
L. X. Bai,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
H. Cai,
J. T. Cai,
Zhe Cao,
J. Chang,
J. F. Chang,
B. M. Chen,
E. S. Chen,
J. Chen,
Liang Chen,
Liang Chen,
Long Chen,
M. J. Chen,
M. L. Chen,
Q. H. Chen
, et al. (250 additional authors not shown)
Abstract:
The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ pet…
▽ More
The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ petaelectronvolt (PeV). The ultra-high-energy photons exhibit the presence of a PeV electron accelerator (a pevatron) with an acceleration rate exceeding 15% of the absolute theoretical limit. Assuming that unpulsed $γ$-rays are produced at the termination of the pulsar's wind, we constrain the pevatron's size, between $0.025$ and $0.1$ pc, and the magnetic field $\approx 110 μ$G. The production rate of PeV electrons, $2.5 \times 10^{36}$ erg $\rm s^{-1}$, constitutes 0.5% of the pulsar's spin-down luminosity, although we do not exclude a non-negligible contribution of PeV protons to the production of the highest energy $γ$-rays.
△ Less
Submitted 11 November, 2021;
originally announced November 2021.
-
Accretion Flows or Outflow Cavities? Uncovering the Gas Dynamics around Lupus 3-MMS
Authors:
Travis J. Thieme,
Shih-Ping Lai,
Sheng-Jun Lin,
Pou-Ieng Cheong,
Chin-Fei Lee,
Hsi-Wei Yen,
Zhi-Yun Li,
Ka Ho Lam,
Bo Zhao
Abstract:
Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flo…
▽ More
Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flows in these early stages of protostellar formation still remain poorly characterized due to their low intensity and possibly some extended structures are disregarded as being part of the outflow cavity. We use ALMA archival data of a young Class 0 protostar, Lupus 3-MMS, to uncover four extended accretion flow-like structures in C$^{18}$O that follow the edges of the outflows. We make various types of position-velocity cuts to compare with the outflows and find the extended structures are not consistent with the outflow emission, but rather more consistent with a simple infall model. We then use a dendrogram algorithm to isolate five sub-structures in position-position-velocity space. Four out of the five sub-structures fit well ($>$95%) with our simple infall model, with specific angular momenta between $2.7-6.9\times10^{-4}\,$km$\,$s$^{-1}\,$pc and mass-infall rates of $0.5-1.1\times10^{-6}\,M_{\odot}\,$yr$^{-1}$. Better characterization of the physical structure in the supposed "outflow-cavities" is important to disentangle the true outflow cavities and accretion flows.
△ Less
Submitted 7 November, 2021;
originally announced November 2021.
-
The young protostellar disk in IRAS16293-2422 B is hot and shows signatures of gravitational instability
Authors:
Joaquin Zamponi,
María José Maureira,
Bo Zhao,
Hauyu Baobab Liu,
John D. Ilee,
Duncan Forgan,
Paola Caselli
Abstract:
Deeply embedded protostars are actively fed from their surrounding envelopes through their protostellar disk. The physical structure of such early disks might be different from that of more evolved sources due to the active accretion. We present 1.3 and 3\,mm ALMA continuum observations at resolutions of 6.5\,au and 12\,au respectively, towards the Class 0 source IRAS 16293-2422 B. The resolved br…
▽ More
Deeply embedded protostars are actively fed from their surrounding envelopes through their protostellar disk. The physical structure of such early disks might be different from that of more evolved sources due to the active accretion. We present 1.3 and 3\,mm ALMA continuum observations at resolutions of 6.5\,au and 12\,au respectively, towards the Class 0 source IRAS 16293-2422 B. The resolved brightness temperatures appear remarkably high, with $T_{\rm b} >$ 100\,K within $\sim$30\,au and $T_{\rm b}$ peak over 400\,K at 3\,mm. Both wavelengths show a lopsided emission with a spectral index reaching values less than 2 in the central $\sim$ 20\,au region. We compare these observations with a series of radiative transfer calculations and synthetic observations of magnetohydrodynamic and radiation hydrodynamic protostellar disk models formed after the collapse of a dense core. Based on our results, we argue that the gas kinematics within the disk may play a more significant role in heating the disk than the protostellar radiation. In particular, our radiation hydrodynamic simulation of disk formation, including heating sources associated with gravitational instabilities, is able to generate the temperatures necessary to explain the high fluxes observed in IRAS 16293B. Besides, the low spectral index values are naturally reproduced by the high optical depth and high inner temperatures of the protostellar disk models. The high temperatures in IRAS 16293B imply that volatile species are mostly in the gas phase, suggesting that a self-gravitating disk could be at the origin of a hot corino.
△ Less
Submitted 14 September, 2021;
originally announced September 2021.
-
No impact of core-scale magnetic field, turbulence, or velocity gradient on sizes of protostellar disks in Orion A
Authors:
Hsi-Wei Yen,
Bo Zhao,
Patrick M. Koch,
Aashish Gupta
Abstract:
We compared the sizes and fluxes of a sample of protostellar disks in Orion A measured with the ALMA 0.87 mm continuum data from the VANDAM survey with the physical properties of their ambient environments on the core scale of 0.6 pc estimated with the GBT GAS NH3 and JCMT SCUPOL polarimetric data. We did not find any significant dependence of the disk radii and continuum fluxes on a single parame…
▽ More
We compared the sizes and fluxes of a sample of protostellar disks in Orion A measured with the ALMA 0.87 mm continuum data from the VANDAM survey with the physical properties of their ambient environments on the core scale of 0.6 pc estimated with the GBT GAS NH3 and JCMT SCUPOL polarimetric data. We did not find any significant dependence of the disk radii and continuum fluxes on a single parameter on the core scale, such as the non-thermal line width, magnetic field orientation and strength, or magnitude and orientation of the velocity gradient. Among these parameters, we only found a positive correlation between the magnitude of the velocity gradient and the non-thermal line width. Thus, the observed velocity gradients are more likely related to turbulent motion but not large-scale rotation. Our results of no clear dependence of the disk radii on these parameters are more consistent with the expectation from non-ideal MHD simulations of disk formation in collapsing cores, where the disk size is self-regulated by magnetic braking and diffusion, compared to other simulations which only include turbulence and/or a magnetic field misaligned with the rotational axis. Therefore, our results could hint that the non-ideal MHD effects play a more important role in the disk formation. Nevertheless, we cannot exclude the influences on the observed disk size distribution by dynamical interaction in a stellar cluster or amounts of angular momentum on the core scale, which cannot be probed with the current data.
△ Less
Submitted 31 May, 2021;
originally announced June 2021.
-
Multi-Wavelength Monitoring and Reverberation Mapping of a Changing Look Event in the Seyfert Galaxy NGC 3516
Authors:
V. L. Oknyansky,
M. S. Brotherton,
S. S. Tsygankov,
A. V. Dodin,
D. -W. Bao,
B. -X. Zhao,
P. Du,
M. A. Burlak,
N. P. Ikonnikova,
A. M. Tatarnikov,
A. A. Belinski,
A. A. Fedoteva,
N. I. Shatsky,
E. O. Mishin,
S. G. Zheltouhov,
S. A. Potanin,
J. -M. Wang,
J. N. McLane,
H. A. Kobulnicky,
D. A. Dale,
T. E. Zastrocky,
J. Maithil,
K. A. Olson,
C. Adelman,
Z. Carter
, et al. (4 additional authors not shown)
Abstract:
We present the results of photometric and spectroscopic monitoring campaigns of the changing look AGN NGC 3516 carried out in 2018 to 2020 covering the wavelength range from the X-ray to the optical. The facilities included the telescopes of the CMO SAI MSU, the 2.3-m WIRO telescope, and the XRT and UVOT of Swift. We found that NGC 3516 brightened to a high state and could be classified as Sy1.5 d…
▽ More
We present the results of photometric and spectroscopic monitoring campaigns of the changing look AGN NGC 3516 carried out in 2018 to 2020 covering the wavelength range from the X-ray to the optical. The facilities included the telescopes of the CMO SAI MSU, the 2.3-m WIRO telescope, and the XRT and UVOT of Swift. We found that NGC 3516 brightened to a high state and could be classified as Sy1.5 during the late spring of 2020. We have measured time delays in the responses of the Balmer and He II 4686 lines to continuum variations. In the case of the best-characterized broad H-beta line, the delay to continuum variability is about 17 days in the blue wing and is clearly shorter, 9 days, in the red, which is suggestive of inflow. As the broad lines strengthened, the blue side came to dominate the Balmer lines, resulting in very asymmetric profiles with blueshifted peaks during this high state. During the outburst the X-ray flux reached its maximum on 1 April 2020 and it was the highest value ever observed for NGC 3516 by the Swift observatory. The X-ray hard photon index became softer, about 1.8 in the maximum on 21 Apr 2020 compared to the mean about 0.7 during earlier epochs before 2020. We have found that the UV and optical variations correlated well (with a small time delay of 1-2 days) with the X-ray until the beginning of April 2020, but later, until the end of Jun. 2020, these variations were not correlated. We suggest that this fact may be a consequence of partial obscuration by Compton-thick clouds crossing the line of sight.
△ Less
Submitted 11 May, 2021; v1 submitted 22 April, 2021;
originally announced April 2021.
-
Geminga SNR: Possible candidate of local cosmic-ray factory
Authors:
Bing Zhao,
Wei Liu,
Qiang Yuan,
Hong-Bo Hu,
Xiao-Jun Bi,
Han-Rong Wu,
Xun-Xiu Zhou,
Yi-Qing Guo
Abstract:
The precise measurements of energy spectra and anisotropy could help us uncover the local cosmic-ray accelerators. Our recent works have shown that spectral hardening above $200$ GeV in the energy spectra and transition of large-scale anisotropy at $\sim 100$ TeV are of local source origin. Less than $100$ TeV, both spectral hardening and anisotropy explicitly indicate the dominant contribution fr…
▽ More
The precise measurements of energy spectra and anisotropy could help us uncover the local cosmic-ray accelerators. Our recent works have shown that spectral hardening above $200$ GeV in the energy spectra and transition of large-scale anisotropy at $\sim 100$ TeV are of local source origin. Less than $100$ TeV, both spectral hardening and anisotropy explicitly indicate the dominant contribution from nearby sources. In this work, we further investigate the parameter space of sources allowed by the observational energy spectra and anisotropy amplitude. To obtain the best-fit source parameters, a numerical package to compute the parameter posterior distributions based on Bayesian inference, which is applied to perform an elaborate scan of parameter space. We find that by combining the energy spectra and anisotropy data, the permissible range of location and age of local source is considerably reduced. When comparing with the current local SNR catalog, only Geminga SNR could be the proper candidate of the local cosmic-ray source.
△ Less
Submitted 15 April, 2021;
originally announced April 2021.
-
East Asian VLBI Network Observations of Active Galactic Nuclei Jets: Imaging with KaVA+Tianma+Nanshan
Authors:
Yuzhu Cui,
Kazuhiro Hada,
Motoki Kino,
Bong Won Sohn,
Jongho Park,
Hyun Wook Ro,
Satoko Sawada-Satoh,
Wu Jiang,
Lang Cui,
Mareki Honma,
Zhi Qiang Shen,
Fumie Tazaki,
Tao An,
Ilje Cho,
Guang Yao Zhao,
Xiao Peng Cheng,
Kotaro Niinuma,
Kiyoaki Wajima,
Ying Kang Zhang,
Noriyuki Kawaguchi,
Juan Carlos Algaba,
Shoko Koyama,
Tomoya Hirota,
Yoshinori Yonekura,
Nobuyuki Sakai
, et al. (52 additional authors not shown)
Abstract:
The East Asian very-long-baseline interferometry (VLBI) Network (EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China, Japan, and Korea. EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions. After the combination of the Korean VLBI Network (KVN) and the VLBI Explorati…
▽ More
The East Asian very-long-baseline interferometry (VLBI) Network (EAVN) is a rapidly evolving international VLBI array that is currently promoted under joint efforts among China, Japan, and Korea. EAVN aims at forming a joint VLBI Network by combining a large number of radio telescopes distributed over East Asian regions. After the combination of the Korean VLBI Network (KVN) and the VLBI Exploration of Radio Astrometry (VERA) into KaVA, further expansion with the joint array in East Asia is actively promoted. Here we report the first imaging results (at 22 and 43 GHz) of bright radio sources obtained with KaVA connected to Tianma 65-m and Nanshan 26-m Radio Telescopes in China. To test the EAVN imaging performance for different sources, we observed four active galactic nuclei (AGN) having different brightness and morphology. As a result, we confirmed that Tianma 65-m Radio Telescope (TMRT) significantly enhances the overall array sensitivity, a factor of 4 improvement in baseline sensitivity and 2 in image dynamic range compared to the case of KaVA only. The addition of Nanshan 26-m Radio Telescope (NSRT) further doubled the east-west angular resolution. With the resulting high-dynamic-range, high-resolution images with EAVN (KaVA+TMRT+NSRT), various fine-scale structures in our targets, such as the counter-jet in M87, a kink-like morphology of the 3C273 jet and the weak emission in other sources, are successfully detected. This demonstrates the powerful capability of EAVN to study AGN jets and to achieve other science goals in general. Ongoing expansion of EAVN will further enhance the angular resolution, detection sensitivity and frequency coverage of the network.
△ Less
Submitted 14 April, 2021; v1 submitted 12 April, 2021;
originally announced April 2021.
-
Calibration of the Air Shower Energy Scale of the Water and Air Cherenkov Techniques in the LHAASO experiment
Authors:
F. Aharonian,
Q. An,
Axikegu,
L. X. Bai,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
H. Cai,
J. T. Cai,
Z. Cao Z. Cao,
J. Chang,
J. F. Chang,
X. C. Chang,
B. M. Chen,
J. Chen,
L. Chen,
L. Chen,
L. Chen,
M. J. Chen,
M. L. Chen,
Q. H. Chen,
S. H. Chen,
S. Z. Chen
, et al. (233 additional authors not shown)
Abstract:
The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galact…
▽ More
The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galactic cosmic rays due to the geomagnetic field. This deflection angle $Δ$ is inversely proportional to the energy of the cosmic rays. The precise measurements of the shifts by WCDA allows us to calibrate its energy scale for energies as high as 35 TeV. The energy scale measured by WCDA can be used to cross calibrate the energy reconstructed by WFCTA, which spans the whole energy range up to 10 PeV. In this work, we will demonstrate the feasibility of the method using the data collected from April 2019 to January 2020 by the WFCTA array and WCDA-1 detector, the first of the three water Cherenkov ponds, already commissioned at LHAASO site.
△ Less
Submitted 13 April, 2021; v1 submitted 11 April, 2021;
originally announced April 2021.
-
FAUST II. Discovery of a Secondary Outflow in IRAS 15398-3359: Variability in Outflow Direction during the Earliest Stage of Star Formation?
Authors:
Yuki Okoda,
Yoko Oya,
Logan Francis,
Doug Johnstone,
Shu-ichiro Inutsuka,
Cecilia Ceccarelli,
Claudio Codella,
Claire Chandler,
Nami Sakai,
Yuri Aikawa,
Felipe Alves,
Nadia Balucani,
Eleonora Bianchi,
Mathilde Bouvier,
Paola Caselli,
Emmanuel Caux,
Steven Charnley,
Spandan Choudhury,
Marta De Simone,
Francois Dulieu,
Aurora Durán,
Lucy Evans,
Cécile Favre,
Davide Fedele,
Siyi Feng
, et al. (44 additional authors not shown)
Abstract:
We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-…
▽ More
We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398-3359, by 1200 au. The arc-like structure is blue-shifted with respect to the systemic velocity. A velocity gradient of 1.2 km/s over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398-3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution.
△ Less
Submitted 18 January, 2021;
originally announced January 2021.
-
Construction and On-site Performance of the LHAASO WFCTA Camera
Authors:
F. Aharonian,
Q. An,
Axikegu,
L. X. Bai,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
H. Cai,
J. T. Cai,
Z. Cao,
Z. Cao,
J. Chang,
J. F. Chang,
X. C. Chang,
B. M. Chen,
J. Chen,
L. Chen,
L. Chen,
L. Chen,
M. J. Chen,
M. L. Chen,
Q. H. Chen,
S. H. Chen
, et al. (234 additional authors not shown)
Abstract:
The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this…
▽ More
The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this application. Eighteen SiPM-based cameras with square light funnels have been built for WFCTA. The telescopes have collected more than 100 million cosmic ray events and preliminary results indicate that these cameras are capable of working under moonlight. The characteristics of the light funnels and SiPMs pose challenges (e.g. dynamic range, dark count rate, assembly techniques). In this paper, we present the design features, manufacturing techniques and performances of these cameras. Finally, the test facilities, the test methods and results of SiPMs in the cameras are reported here.
△ Less
Submitted 4 July, 2021; v1 submitted 29 December, 2020;
originally announced December 2020.
-
Ground Experiments and Performance Evaluation of the Low-Frequency Radio Spectrometer Onboard the Lander of Chang'e-4 Mission
Authors:
Xinying Zhu,
Yan Su,
Yicai Ji,
Hongbo Zhang,
Bo Zhao,
Xiping Xue,
Shun Dai,
Junduo Li,
Chunlai Li
Abstract:
The Low-Frequency Radio Spectrometer (LFRS) is a scientific payload onboard the Chang'e-4 lunar lander launched in December 2018. The LFRS provides in-situ measurements of the low-frequency radio phenomena on the far-side of the Moon for the first time in human history. To evaluate the performance of the LFRS, a series of ground experiments are conducted using a prototype model of the LFRS. It is…
▽ More
The Low-Frequency Radio Spectrometer (LFRS) is a scientific payload onboard the Chang'e-4 lunar lander launched in December 2018. The LFRS provides in-situ measurements of the low-frequency radio phenomena on the far-side of the Moon for the first time in human history. To evaluate the performance of the LFRS, a series of ground experiments are conducted using a prototype model of the LFRS. It is not easy to perform the experiments because the EMI \footnote{ Abbreviation: EMI, Electro Magnetic Interference} from the Chang'e-4 lunar lander itself and the environment is very intense. The results after EMI mitigation show that the sensitivity of the LFRS may be 1e-18Wm^2Hz^-1.
△ Less
Submitted 8 December, 2020;
originally announced December 2020.
-
Monitoring AGNs with Hβ Asymmetry. II. Reverberation Mapping of Three Seyfert Galaxies Historically Displaying Hβ Profiles with Changing Asymmetry: Mrk 79, NGC 3227, and Mrk 841
Authors:
Michael S. Brotherton,
Pu Du,
Ming Xiao,
Dong-Wei Bao,
Bixuan Zhao,
Jacob N. McLane,
Kianna A. Olson,
Kai Wang,
Zheng-Peng Huang,
Chen Hu,
David H. Kasper,
William T. Chick,
My L. Nguyen,
Jaya Maithil,
Derek Hand,
Yan-Rong Li,
Luis C. Ho,
Jin-Ming Bai,
Wei-Hao Bian,
Jian-Min Wang
Abstract:
We report the results of reverberation mapping three bright Seyfert galaxies, Mrk 79, NGC 3227, and Mrk 841, from a campaign conducted from December 2016 to May 2017 with the Wyoming Infrared Observatory (WIRO) 2.3-meter telescope. All three of these targets have shown asymmetric broad H$β$ emission lines in the past, although their emission lines were relatively symmetric during our observations.…
▽ More
We report the results of reverberation mapping three bright Seyfert galaxies, Mrk 79, NGC 3227, and Mrk 841, from a campaign conducted from December 2016 to May 2017 with the Wyoming Infrared Observatory (WIRO) 2.3-meter telescope. All three of these targets have shown asymmetric broad H$β$ emission lines in the past, although their emission lines were relatively symmetric during our observations. We measured Hβ time lags for all three targets and estimated masses of their black holes -- for the first time in the case of Mrk 841. For Mrk 79 and NGC 3227, the data are of sufficient quality to resolve distinct time lags as a function of velocity and to compute two-dimensional velocity-delay maps. Mrk 79 shows smaller time lags for high-velocity gas but the distribution is not symmetric, and its complex velocity-delay map could result from the combination of both inflowing and outflowing Hβ emitting disks that may be part of a single larger structure. NGC 3227 shows the largest time lags for blueshifted gas and the two-dimensional velocity-delay map suggests a disk with some inflow. We compare our results with previous work and find evidence for different time lags despite similar luminosities, as well as evolving broad line region structures.
△ Less
Submitted 11 November, 2020;
originally announced November 2020.
-
The observation of the Crab Nebula with LHAASO-KM2A for the performance study
Authors:
F. Aharonian,
Q. An,
Axikegu,
L. X. Bai,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
H. Cai,
J. T. Cai,
Z. Cao,
Z. Cao,
J. Chang,
J. F. Chang,
X. C. Chang,
B. M. Chen,
J. Chen,
L. Chen,
L. Chen,
L. Chen,
M. J. Chen,
M. L. Chen,
Q. H. Chen,
S. H. Chen
, et al. (234 additional authors not shown)
Abstract:
As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the…
▽ More
As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power.
The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered.
△ Less
Submitted 13 October, 2020;
originally announced October 2020.
-
The Interplay between Ambipolar Diffusion and Hall Effect on Magnetic Field Decoupling and Protostellar Disc Formation
Authors:
Bo Zhao,
Paola Caselli,
Zhi-Yun Li,
Ruben Krasnopolsky,
Hsien Shang,
Ka Ho Lam
Abstract:
Non-ideal MHD effects have been shown recently as a robust mechanism of averting the magnetic braking "catastrophe" and promoting protostellar disc formation. However, the magnetic diffusivities that determine the efficiency of non-ideal MHD effects are highly sensitive to microphysics. We carry out non-ideal MHD simulations to explore the role of microphysics on disc formation and the interplay b…
▽ More
Non-ideal MHD effects have been shown recently as a robust mechanism of averting the magnetic braking "catastrophe" and promoting protostellar disc formation. However, the magnetic diffusivities that determine the efficiency of non-ideal MHD effects are highly sensitive to microphysics. We carry out non-ideal MHD simulations to explore the role of microphysics on disc formation and the interplay between ambipolar diffusion (AD) and Hall effect during the protostellar collapse. We find that removing the smallest grain population ($\lesssim$10 nm) from the standard MRN size distribution is sufficient for enabling disc formation. Further varying the grain sizes can result in either a Hall-dominated or an AD-dominated collapse; both form discs of tens of AU in size regardless of the magnetic field polarity. The direction of disc rotation is bimodal in the Hall dominated collapse but unimodal in the AD-dominated collapse. We also find that AD and Hall effect can operate either with or against each other in both radial and azimuthal directions, yet the combined effect of AD and Hall is to move the magnetic field radially outward relative to the infalling envelope matter. In addition, microphysics and magnetic field polarity can leave profound imprints both on observables (e.g., outflow morphology, disc to stellar mass ratio) and on the magnetic field characteristics of protoplanetary discs. Including Hall effect relaxes the requirements on microphysics for disc formation, so that prestellar cores with cosmic-ray ionization rate of $\lesssim$2--3$\times10^{-16}$ s$^{-1}$ can still form small discs of $\lesssim$10 AU radius. We conclude that disc formation should be relatively common for typical prestellar core conditions, and that microphysics in the protostellar envelope is essential to not only disc formation, but also protoplanetary disc evolution.
△ Less
Submitted 16 September, 2020;
originally announced September 2020.
-
Hall Effect in Protostellar Disc Formation and Evolution
Authors:
Bo Zhao,
Paola Caselli,
Zhi-Yun Li,
Ruben Krasnopolsky,
Hsien Shang,
Ka Ho Lam
Abstract:
The Hall effect is recently shown to be efficient in magnetized dense molecular cores, and could lead to a bimodal formation of rotationally supported discs (RSDs) in the first core phase. However, how such Hall dominated systems evolve in the protostellar accretion phase remains unclear. We carry out 2D axisymmetric simulations including Hall effect and Ohmic dissipation, with realistic magnetic…
▽ More
The Hall effect is recently shown to be efficient in magnetized dense molecular cores, and could lead to a bimodal formation of rotationally supported discs (RSDs) in the first core phase. However, how such Hall dominated systems evolve in the protostellar accretion phase remains unclear. We carry out 2D axisymmetric simulations including Hall effect and Ohmic dissipation, with realistic magnetic diffusivities computed from our equilibrium chemical network. We find that Hall effect only becomes efficient when the large population of very small grains (VSGs: $\lesssim$10 nm) is removed from the standard MRN size distribution. With such an enhanced Hall effect, however, the bimodality of disc formation does not continue into the main accretion phase. The outer part of the initial $\sim$40 AU disc formed in the anti-aligned configuration (${\bf Ω\cdot B}<0$) flattens into a thin rotationally supported Hall current sheet as Hall effect moves the poloidal magnetic field radially inward relative to matter, leaving only the inner $\lesssim$10--20 AU RSD. In the aligned configuration (${\bf Ω\cdot B}>0$), disc formation is suppressed initially but a counter-rotating disc forms subsequently due to efficient azimuthal Hall drift. The counter-rotating disc first grows to $\sim$30 AU as Hall effect moves the magnetic field radially outward, but only the inner $\lesssim$10 AU RSD is long-lived like in the anti-aligned case. Besides removing VSGs, cosmic ray ionization rate should be below a few 10$^{-16}$ s$^{-1}$ for Hall effect to be efficient in disc formation. We conclude that Hall effect produces small $\lesssim$10--20 AU discs regardless of the polarity of the magnetic field, and that radially outward diffusion of magnetic fields remains crucial for disc formation and growth.
△ Less
Submitted 16 September, 2020;
originally announced September 2020.
-
The supermassive black hole masses of reverberation-mapped active galactic nuclei
Authors:
Li-Ming Yu,
Wei-Hao Bian,
Xue-Guang Zhang,
Bi-Xuan Zhao,
Chan Wang,
Xue Ge,
Bing-Qian Zhu,
Yu-Qin Chen
Abstract:
Using different kinds of velocity tracers derived from the broad H$β$ profile (in the mean or rms spectrum) and the corresponding virial factors $f$, the central supermassive black hole (SMBH) masses ($M_{\rm BH}$) are calculated for a compiled sample of 120 reverberation-mapped (RM) AGNs. For its subsample of RM AGNs with measured stellar velocity dispersion ($σ_{\rm \ast}$), the multivariate lin…
▽ More
Using different kinds of velocity tracers derived from the broad H$β$ profile (in the mean or rms spectrum) and the corresponding virial factors $f$, the central supermassive black hole (SMBH) masses ($M_{\rm BH}$) are calculated for a compiled sample of 120 reverberation-mapped (RM) AGNs. For its subsample of RM AGNs with measured stellar velocity dispersion ($σ_{\rm \ast}$), the multivariate linear regression technique is used to calibrate the mean value $f$, as well as the variable FWHM-based $f$. It is found that, whether excluding the pseudo-bulges or not, $M_{\rm BH}$ from the H$β$ line dispersion in the mean spectrum ($σ_{\rm Hβ,mean}$) has the smallest offset rms with respect to the $M_{\rm BH}-σ_{\ast}$ relation. For the total sample excluding SDSS-RM AGNs, with respect to $M_{\rm BH}$ from $σ_{\rm \ast}$ or that from the H$β$ line dispersion in the rms spectrum ($σ_{\rm Hβ,rms}$), it is found that we can obtain $M_{\rm BH}$ from the $σ_{\rm Hβ,mean}$ with the smallest offset rms of 0.38 dex or 0.23 dex, respectively. It implies that, with respect to the H$β$ FWHM, we prefer $σ_{\rm Hβ,mean}$ to calculate $M_{\rm BH}$ from the single-epoch spectrum. Using the FWHM-based $f$, we can improve $M_{\rm BH}$ calculation from FWHM(H$β$) and the mean $f$, with a decreased offset rms from 0.52 dex to 0.39 dex with respect to $M_{\rm BH}$ from $σ_{\rm \ast}$ for the subsample of 36 AGNs with $σ_{\rm \ast}$. The value of 0.39 dex is almost the same as that from $σ_{\rm Hβ,mean}$ and the mean $f$.
△ Less
Submitted 14 August, 2020;
originally announced August 2020.
-
A protostellar system fed by a streamer of 10,500 au length
Authors:
Jaime E. Pineda,
Dominique Segura-Cox,
Paola Caselli,
Nichol Cunningham,
Bo Zhao,
Anika Schmiedeke,
Maria José Maureira,
Roberto Neri
Abstract:
Binary formation is an important aspect of star formation. One possible route for close-in binary formation is disk fragmentation$^{[1,2,3]}$. Recent observations show small scale asymmetries (<300 au) around young protostars$^{[2,4]}$, although not always resolving the circumbinary disk, are linked to disk phenomena$^{[5,6]}$. In later stages, resolved circumbinary disk observations$^{[7]}$ (<200…
▽ More
Binary formation is an important aspect of star formation. One possible route for close-in binary formation is disk fragmentation$^{[1,2,3]}$. Recent observations show small scale asymmetries (<300 au) around young protostars$^{[2,4]}$, although not always resolving the circumbinary disk, are linked to disk phenomena$^{[5,6]}$. In later stages, resolved circumbinary disk observations$^{[7]}$ (<200 au) show similar asymmetries, suggesting the origin of the asymmetries arises from binary-disk interactions$^{[8,9,10]}$. We observed one of the youngest systems to study the connection between disk and dense core. We find for the first time a bright and clear streamer in chemically fresh material (Carbon-chain species) that originates from outside the dense core (>10,500 au). This material connects the outer dense core with the region where asymmetries arise near disk scales. This new structure type, 10x larger than those seen near disk scales, suggests a different interpretation of previous observations: large-scale accretion flows funnel material down to disk scales. These results reveal the under-appreciated importance of the local environment on the formation and evolution of disks in early systems$^{[13,14]}$ and a possible initial condition for the formation of annular features in young disks$^{[15,16]}$.
△ Less
Submitted 27 July, 2020;
originally announced July 2020.
-
FAUST I. The hot corino at the heart of the prototypical Class I protostar L1551 IRS5
Authors:
E. Bianchi,
C. J. Chandler,
C. Ceccarelli,
C. Codella,
N. Sakai,
A. López-Sepulcre,
L. T. Maud,
G. Moellenbrock,
B. Svoboda,
Y. Watanabe,
T. Sakai,
F. Ménard,
Y. Aikawa,
F. Alves,
N. Balucani,
M. Bouvier,
P. Caselli,
E. Caux,
S. Charnley,
S. Choudhury,
M. De Simone,
F. Dulieu,
A. Durán,
L. Evans,
C. Favre
, et al. (41 additional authors not shown)
Abstract:
The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I p…
▽ More
The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I protostars has become of paramount importance. Here we report the discovery of a hot corino towards the prototypical Class I protostar L1551 IRS5, obtained within the ALMA Large Program FAUST. We detected several lines from methanol and its isopotologues ($^{13}$CH$_{\rm 3}$OH and CH$_{\rm 2}$DOH), methyl formate and ethanol. Lines are bright toward the north component of the IRS5 binary system, and a possible second hot corino may be associated with the south component. The methanol lines non-LTE analysis constrains the gas temperature ($\sim$100 K), density ($\geq$1.5$\times$10$^{8}$ cm$^{-3}$), and emitting size ($\sim$10 au in radius). All CH$_{\rm 3}$OH and $^{13}$CH$_{\rm 3}$OH lines are optically thick, preventing a reliable measure of the deuteration. The methyl formate and ethanol relative abundances are compatible with those measured in Class 0 hot corinos. Thus, based on the present work, little chemical evolution from Class 0 to I hot corinos occurs.
△ Less
Submitted 20 July, 2020;
originally announced July 2020.
-
Rapid Elimination of Small Dust Grains in Molecular Clouds
Authors:
Kedron Silsbee,
Alexei Ivlev,
Paola Caselli,
Olli Sipila,
Bo Zhao
Abstract:
We argue that impact velocities between dust grains with sizes less than $\sim 0.1$ $μm$ in molecular cloud cores are dominated by drift arising from ambipolar diffusion. This effect is due to the size dependence of the dust coupling to the magnetic field and the neutral gas. Assuming perfect sticking in collisions up to $\approx 50$ m/s, we show that this effect causes rapid depletion of small gr…
▽ More
We argue that impact velocities between dust grains with sizes less than $\sim 0.1$ $μm$ in molecular cloud cores are dominated by drift arising from ambipolar diffusion. This effect is due to the size dependence of the dust coupling to the magnetic field and the neutral gas. Assuming perfect sticking in collisions up to $\approx 50$ m/s, we show that this effect causes rapid depletion of small grains - consistent with starlight extinction and IR/microwave emission measurements, both in the core center ($n \sim 10^{6}$ cm$^{-3}$) and envelope ($n \sim 10^{4}$ cm$^{-3}$). The upper end of the size distribution does not change significantly if only velocities arising from this effect are considered. We consider the impact of an evolved dust size distribution on the gas temperature, and argue that if the depletion of small dust grains occurs as would be expected from our model, then the cosmic ray ionization rate must be well below $10^{-16}$ s$^{-1}$ at a number density of $10^{5}$ cm$^{-3}$.
△ Less
Submitted 29 September, 2020; v1 submitted 23 June, 2020;
originally announced June 2020.
-
Effect of grain size distribution and size-dependent grain heating on molecular abundances in starless and pre-stellar cores
Authors:
O. Sipilä,
B. Zhao,
P. Caselli
Abstract:
We present a new gas-grain chemical model to constrain the effect of grain size distribution on molecular abundances in starless and pre-stellar cores. We introduce grain-size dependence simultaneously for cosmic-ray (CR)-induced desorption efficiency and for grain equilibrium temperatures. We keep explicit track of ice abundances on a set of grain populations. We find that the size-dependent CR d…
▽ More
We present a new gas-grain chemical model to constrain the effect of grain size distribution on molecular abundances in starless and pre-stellar cores. We introduce grain-size dependence simultaneously for cosmic-ray (CR)-induced desorption efficiency and for grain equilibrium temperatures. We keep explicit track of ice abundances on a set of grain populations. We find that the size-dependent CR desorption efficiency affects ice abundances in a highly non-trivial way that depends on the molecule. Species that originate in the gas phase follow a simple pattern where the ice abundance is highest on the smallest grains (the most abundant in the distribution). Some molecules, such as HCN, are instead concentrated on large grains throughout the time evolution, while others (like $\rm N_2$) are initially concentrated on large grains, but at late times on small grains, due to grain-size-dependent competition between desorption and hydrogenation. Most of the water ice is on small grains at high medium density ($n({\rm H_2}) \gtrsim 10^6 \, \rm cm^{-3}$), where the water ice fraction, with respect to total water ice reservoir, can be as low as $\sim 10^{-3}$ on large (> 0.1 $μ$m) grains. Allowing the grain equilibrium temperature to vary with grain size induces strong variations in relative ice abundances in low-density conditions where the interstellar radiation field and in particular its ultraviolet component are not attenuated. Our study implies consequences not only for the initial formation of ices preceding the starless core stage, but also for the relative ice abundances on the grain populations going into the protostellar stage. In particular, if the smallest grains can lose their mantles due to grain-grain collisions as the core is collapsing, the ice composition in the beginning of the protostellar stage could be very different to that in the pre-collapse phase.
△ Less
Submitted 17 June, 2020;
originally announced June 2020.
-
Ionization: a possible explanation for the difference of mean disk sizes in star-forming regions
Authors:
M. Kuffmeier,
B. Zhao,
P. Caselli
Abstract:
Surveys of protoplanetary disks in star-forming regions of similar age revealed significant variations in average disk mass between some regions. For instance, disks in the Orion Nebular Cluster (ONC) and Corona Australis (CrA) are on average smaller than disks observed in Lupus, Taurus, Chamaeleon I or Ophiuchus. In contrast to previous models that study truncation of disks at a late stage of the…
▽ More
Surveys of protoplanetary disks in star-forming regions of similar age revealed significant variations in average disk mass between some regions. For instance, disks in the Orion Nebular Cluster (ONC) and Corona Australis (CrA) are on average smaller than disks observed in Lupus, Taurus, Chamaeleon I or Ophiuchus. In contrast to previous models that study truncation of disks at a late stage of their evolution, we investigate whether disks may already be born with systematically smaller disk sizes in more massive star-forming regions as a consequence of enhanced ionization rates. Assuming various cosmic-ray ionization rates, we compute the resistivities for ambipolar diffusion and Ohmic dissipation with a chemical network, and perform 2D non-ideal magnetohydrodynamical protostellar collapse simulations. A higher ionization rate leads to stronger magnetic braking, and hence to the formation of smaller disks. Accounting for recent findings that protostars act as forges of cosmic rays and considering only mild attenuation during the collapse phase, we show that a high average cosmic-ray ionization rate in star-forming regions like the ONC or CrA can explain the detection of smaller disks in these regions. Our results show that on average a higher ionization rate leads to the formation of smaller disks. Therefore, smaller disks in regions of similar age can be the consequence of different levels of ionization, and may not exclusively be caused by disk truncation via external photoevaporation. We strongly encourage observations that allow measuring the cosmic-ray ionization degrees in different star-forming regions to test our hypothesis.
△ Less
Submitted 5 June, 2020; v1 submitted 29 May, 2020;
originally announced June 2020.
-
Formation and Evolution of Disks around Young Stellar Objects
Authors:
Bo Zhao,
Kengo Tomida,
Patrick Hennebelle,
John J. Tobin,
Anaelle Maury,
Tomoya Hirota,
Álvaro Sánchez-Monge,
Rolf Kuiper,
Anna Rosen,
Asmita Bhandare,
Marco Padovani,
Yueh-Ning Lee
Abstract:
Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking t…
▽ More
Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
Transition from ordered pinched to warped magnetic field on a 100 au scale in the Class 0 protostar B335
Authors:
Hsi-Wei Yen,
Bo Zhao,
Patrick Koch,
Ruben Krasnopolsky,
Zhi-Yun Li,
Nagayoshi Ohashi,
Hsien Shang,
Shigehisa Takakuwa,
Ya-Wen Tang
Abstract:
We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarizat…
▽ More
We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarization percentages are systematically higher at 1.3 mm than 0.87 mm by a factor of ~1.7, suggesting that the polarized emission originates from magnetically aligned dust grains. We inferred the magnetic field orientations from the observed polarization orientations. We found that the magnetic field changes from ordered and highly pinched to more complicated and asymmetric structures within the inner 100 au scale of B335, and the magnetic field connects to the center along the equatorial plane as well as along the directions which are ~40-60 degrees from the equatorial plane. We performed non-ideal MHD simulations of collapsing dense cores. We found that similar magnetic field structures appear in our simulations of dense cores with the magnetic field and rotational axis slightly misaligned by 15 degrees but not in those with the aligned magnetic field and rotational axis. Our results suggest that the midplane of the inner envelope within the inner 100 au scale of B335 could be warped because of the misaligned magnetic field and rotational axis, and the magnetic field could be dragged by the warped accretion flows.
△ Less
Submitted 9 March, 2020;
originally announced March 2020.
-
Impact of magneto-rotational instability on grain growth in protoplanetary disks: I. Relevant turbulence properties
Authors:
Munan Gong,
Alexei V. Ivlev,
Bo Zhao,
Paola Caselli
Abstract:
Turbulence in the protoplanetary disks induces collisions between dust grains, and thus facilitates grain growth. We investigate the two fundamental assumptions of the turbulence in obtaining grain collisional velocities -- the kinetic energy spectrum and the turbulence autocorrelation time -- in the context of the turbulence generated by the magneto-rotational instability (MRI). We carry out nume…
▽ More
Turbulence in the protoplanetary disks induces collisions between dust grains, and thus facilitates grain growth. We investigate the two fundamental assumptions of the turbulence in obtaining grain collisional velocities -- the kinetic energy spectrum and the turbulence autocorrelation time -- in the context of the turbulence generated by the magneto-rotational instability (MRI). We carry out numerical simulations of the MRI as well as driven turbulence, for a range of physical and numerical parameters. We find that the convergence of the turbulence $α$-parameter does not necessarily imply the convergence of the energy spectrum. The MRI turbulence is largely solenoidal, for which we observe a persistent kinetic energy spectrum of $k^{-4/3}$. The same is obtained for solenoidal driven turbulence with and without magnetic field, over more than 1 dex near the dissipation scale. This power-law slope appears to be converged in terms of numerical resolution, and to be due to the bottleneck effect. The kinetic energy in the MRI turbulence peaks at the fastest growing mode of the MRI. In contrast, the magnetic energy peaks at the dissipation scale. The magnetic energy spectrum in the MRI turbulence does not show a clear power-law range, and is almost constant over approximately 1 dex near the dissipation scale. The turbulence autocorrelation time is nearly constant at large scales, limited by the shearing timescale, and shows a power-law drop close to $k^{-1}$ at small scales, with a slope steeper than that of the eddy crossing time. The deviation from the standard picture of the Kolmogorov turbulence with the injection scale at the disk scale height can potentially have a significant impact on the grain collisional velocities.
△ Less
Submitted 12 February, 2020;
originally announced February 2020.
-
An extended size-luminosity relation for the reverberation-mapped AGNs: the role of the accretion rate
Authors:
Li-Ming Yu,
Bi-Xuan Zhao,
Wei-Hao Bian,
Chan Wang,
Xue Ge
Abstract:
For a compiled sample of 120 reverberation-mapped AGNs, the bivariate correlations of the broad-line regions (BLRs) size ($R_{\rm BLR}$) with the continuum luminosity at 5100 Å ($L_{5100}$) and the dimensionless accretion rates ($\dot{\mathscr{M}}$) are investigated. Using our recently calibrated virial factor $f$, and the velocity tracer from the H$β$ Full-width at half-maximum (FWHM(H$β$)) or th…
▽ More
For a compiled sample of 120 reverberation-mapped AGNs, the bivariate correlations of the broad-line regions (BLRs) size ($R_{\rm BLR}$) with the continuum luminosity at 5100 Å ($L_{5100}$) and the dimensionless accretion rates ($\dot{\mathscr{M}}$) are investigated. Using our recently calibrated virial factor $f$, and the velocity tracer from the H$β$ Full-width at half-maximum (FWHM(H$β$)) or the line dispersion ($σ_{\rm Hβ}$) measured in the mean spectra, three kinds of SMBH masses and $\dot{\mathscr{M}}$ are calculated. An extended \RL relation including $\dot{\mathscr{M}}$ is found to be stronger than the canonical $R_{\rm BLR}({\rm Hβ}) - L_{\rm 5100}$ relation, showing smaller scatters. The observational parameters, $R_{\rm Fe}$ (the ratio of optical Fe II to H$β$ line flux) and the line profile parameter $D_{\rm Hβ}$ ($D_{\rm Hβ}=\rm FWHM(Hβ)/σ_{\rm Hβ}$), have relations with three kinds of $\dot{\mathscr{M}}$. Using $R_{\rm Fe}$ and $D_{\rm Hβ}$ to substitute $\dot{\mathscr{M}}$, extended empirical $R_{\rm BLR}({\rm Hβ}) - L_{\rm 5100}$ relations are presented. $R_{\rm Fe}$ is a better "fix" for the $R_{\rm BLR}({\rm Hβ}) - L_{\rm 5100}$ offset than the H$β$ shape $D_{\rm Hβ}$. The extended empirical $R_{\rm BLR}({\rm Hβ}) - L_{\rm 5100}$ relation including $R_{\rm Fe}$ can be used to calculate $R_{\rm BLR}$, and thus the single-epoch SMBH mass $M_{\rm BH}$. Our measured accretion rate dependence is not consistent with the simple model of the accretion disk instability leading the BLRs formation. The BLR may instead form from the inner edge of the torus, or from some other means in which BLR size is positively correlated with accretion rate and the SMBH mass.
△ Less
Submitted 25 November, 2019;
originally announced November 2019.
-
Gas flow and accretion via spiral streamers and circumstellar disks in a young binary protostar
Authors:
Felipe O. Alves,
Paola Caselli,
Josep M. Girart,
Dominique Segura-Cox,
Gabriel A. P. Franco,
Anika Schmiedeke,
Bo Zhao
Abstract:
The majority of stars are part of gravitationally bound stellar systems, such as binaries. Observations of protobinary systems constrain the conditions that lead to stellar multiplicity and subsequent orbital evolution. We report high-angular resolution observations of the circumbinary disk around [BHB2007] 11, a young binary protostar system. The two protostars are embedded in circumstellar disks…
▽ More
The majority of stars are part of gravitationally bound stellar systems, such as binaries. Observations of protobinary systems constrain the conditions that lead to stellar multiplicity and subsequent orbital evolution. We report high-angular resolution observations of the circumbinary disk around [BHB2007] 11, a young binary protostar system. The two protostars are embedded in circumstellar disks that have radii of 2 to 3 astronomical units and probably contain a few Jupiter masses. These systems are surrounded by a complex structure of filaments connecting to the larger circumbinary disk. We also observe accretion and radio jets associated with the protobinary system. The accretion is preferentially onto the lower-mass protostar, consistent with theoretical predictions.
△ Less
Submitted 2 October, 2019;
originally announced October 2019.