Abstract visibility
Abstract text size

Papers for Thursday, Aug 05 2021

Papers with local authors

Oliver H. E. Philcox, Jiamin Hou, Zachary Slepian

26 pages, 9 figures, submitted to Phys. Rev. D. Code available at this https URL

0 votes
Paper 1 — arXiv:2108.01670
0 votes
Paper 1 — arXiv:2108.01670

We present an $8.1\sigma$ detection of the non-Gaussian 4-Point Correlation Function (4PCF) using a sample of $N_{\rm g} \approx 8\times 10^5$ galaxies from the BOSS CMASS dataset. Our measurement uses the $\mathcal{O}(N_{\rm g}^2)$ NPCF estimator of Philcox et al. (2021), including a new modification to subtract the disconnected 4PCF contribution (arising from the product of two 2PCFs) at the estimator level. This approach is unlike previous work and ensures that our signal is a robust detection of gravitationally-induced non-Gaussianity. The estimator is validated with a suite of lognormal simulations, and the analytic form of the disconnected contribution is discussed. Due to the high dimensionality of the 4PCF, data compression is required; we use a signal-to-noise-based scheme calibrated from theoretical covariance matrices to restrict to $\sim$ $100$ basis vectors. The compression has minimal impact on the detection significance and facilitates traditional $\chi^2$-like analyses using a suite of mock catalogs. The significance is stable with respect to different treatments of noise in the sample covariance (arising from the limited number of mocks), but decreases to $4.7\sigma$ when a minimum galaxy separation of $14 h^{-1}\mathrm{Mpc}$ is enforced on the 4PCF tetrahedra (such that the statistic can be modelled more easily). The detectability of the 4PCF in the quasi-linear regime implies that it will become a useful tool in constraining cosmological and galaxy formation parameters from upcoming spectroscopic surveys.

The hot Jupiter WASP-12b has been found to be on a decaying 1.09-day orbit. The mean density of the planet inferred from transit and radial-velocity data is near its Roche limit; just how near depends on the planet's uncertain internal structure. There is also spectroscopic evidence of mass loss. We accurately calculate the Roche density on the assumption of a synchronously rotating $n=1$ polytrope, and find this to be only $15-20\%$ below the observational estimates for the mean density. We discuss the implied shape of the planet, its lifetime before complete disruption, and its current rate of mass loss based on our improved tidal model. The present mass-loss rate is at least as sensitive to the opacity and temperature profiles of the planet's atmosphere as to its internal structure, however.

Federica B. Bianco, Željko Ivezić, R. Lynne Jones, Melissa L. Graham, Phil Marshall, Abhijit Saha, Michael A. Strauss, Peter Yoachim, Tiago Ribeiro, Timo Anguita, Franz E. Bauer, Eric C. Bellm, Robert D. Blum, William N. Brandt, Sarah Brough, Màrcio Catelan, William I. Clarkson, Andrew J. Connolly, Eric Gawiser, John Gizis, Renee Hlozek, Sugata Kaviraj, Charles T. Liu, Michelle Lochner, Ashish A. Mahabal, Rachel Mandelbaum, Peregrine McGehee, Eric H. Neilsen Jr., Knut A. G. Olsen, Hiranya Peiris, Jason Rhodes, Gordon T. Richards, Stephen Ridgway, Megan E. Schwamb, Dan Scolnic, Ohad Shemmer, Colin T. Slater, Anže Slosar, Stephen J. Smartt, Jay Strader, Rachel Street, David E. Trilling, Aprajita Verma, A. K. Vivas, Risa H. Wechsler, Beth Willman

Submitted as the opening paper of the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

0 votes
Paper 6 — arXiv:2108.01683
0 votes
Paper 6 — arXiv:2108.01683

Vera C. Rubin Observatory is a ground-based astronomical facility under construction, a joint project of the National Science Foundation and the U.S. Department of Energy, designed to conduct a multi-purpose 10-year optical survey of the southern hemisphere sky: the Legacy Survey of Space and Time. Significant flexibility in survey strategy remains within the constraints imposed by the core science goals of probing dark energy and dark matter, cataloging the Solar System, exploring the transient optical sky, and mapping the Milky Way. The survey's massive data throughput will be transformational for many other astrophysics domains and Rubin's data access policy sets the stage for a huge potential users' community. To ensure that the survey science potential is maximized while serving as broad a community as possible, Rubin Observatory has involved the scientific community at large in the process of setting and refining the details of the observing strategy. The motivation, history, and decision-making process of this strategy optimization are detailed in this paper, giving context to the science-driven proposals and recommendations for the survey strategy included in this Focus Issue.

Jiamin Hou, Robert N. Cahn, Oliver H. E. Philcox, Zachary Slepian

31 pages, 17 figures

0 votes
Paper 10 — arXiv:2108.01714
0 votes
Paper 10 — arXiv:2108.01714

We derive analytic covariance matrices for the $N$-Point Correlation Functions (NPCFs) of galaxies in the Gaussian limit. Our results are given for arbitrary $N$ and projected onto the isotropic basis functions of Cahn & Slepian (2020), recently shown to facilitate efficient NPCF estimation. A numerical implementation of the 4PCF covariance is compared to the sample covariance obtained from a set of lognormal simulations, Quijote dark matter halo catalogues, and MultiDark-Patchy galaxy mocks, with the latter including realistic survey geometry. The analytic formalism gives reasonable predictions for the covariances estimated from mock simulations with a periodic-box geometry. Furthermore, fitting for an effective volume and number density by maximizing a likelihood based on Kullback-Leibler divergence is shown to partially compensate for the effects of a non-uniform window function.

All other papers

Yutaka Fujita, Nozomu Kawakatu, Hiroshi Nagai

16 pages, 12 figures

Massive molecular gas has been discovered in giant elliptical galaxies at the centers of galaxy clusters. To reveal its role in AGN feedback in those galaxies, we construct a semi-analytic model of gas circulation. This model especially focuses on the massive molecular gas (interstellar cold gas on a scale of ~10 kpc) and the circumnuclear disk (~< 0.5 kpc). We consider the destruction of the interstellar cold gas by star formation and the gravitational instability for the circumnuclear disk. Our model can reproduce the basic properties of the interstellar cold gas and the circumnuclear disk such as their masses. We also find that the circumnuclear disk tends to stay at the boundary between stable and unstable states. This works as an 'adjusting valve' that regulates mass accretion toward the supermassive black hole. On the other hand, the interstellar cold gas serves as a 'fuel tank' in the AGN feedback. Even if the cooling of the galactic hot gas is prevented, the interstellar cold gas can sustain the AGN activities for ~> 0.5 Gyr. We also confirm that the small entropy of the hot gas (~< 30 keV cm^2) or the short cooling time (~< 1 Gyr) is a critical condition for the existence of the massive molecular gas in the galaxy. The dissipation time of the interstellar cold gas may be related to the critical cooling time.

Yu-Ching Chen, Hsiang-Chih Hwang, Yue Shen, Xin Liu, Nadia L. Zakamska, Qian Yang, Jennifer I. Li

20 pages, 10 figures, 3 tables, submitted to ApJ, comments are welcome

Dual supermassive black holes (SMBHs) at $\sim$kpc scales are the progenitor population of SMBH mergers and play an important role in understanding the pairing and dynamical evolution of massive black holes in galaxy mergers. Because of the stringent resolution requirement and the apparent rareness of these small-separation pairs, there are scarce observational constraints on this population, with few confirmed dual SMBHs at $<10$kpc separations at $z>1$. Here we present results from a pilot search for kpc-scale dual quasars selected with Gaia Data release 2 (DR2) astrometry and followed up with Hubble Space Telescope (HST) Wide Field Camera 3 dual-band (F475W and F814W) snapshot imaging. Our targets are quasars primarily selected with the varstrometry technique, i.e., light centroid jitter caused by asynchronous variability from both members in an unresolved quasar pair, supplemented by sub-arcsec pairs already resolved by Gaia DR2. We find an overall high fraction of HST-resolved pairs among the varstrometry-selected quasars (unresolved in Gaia DR2), $\sim 30-50\%$, increasing toward high redshift ($\sim 60-80\%$ at $z>1.5$). We discuss the nature of the 43 resolved sub-arcsec pairs based on HST and supplementary data. A substantial fraction ($\sim 40\%$) of these pairs are likely physical quasar pairs or gravitationally lensed quasars. We also discover a triple quasar candidate and a quadruply lensed quasar, which is among the smallest-separation quadruple lenses. These results provide important guidelines to improve varstrometry selection and follow-up confirmation of $\sim$kpc-scale dual SMBHs at high redshift.

We here introduce indicator functions, which identify regions of a given density in order to characterize the density dependence of clustering. After a general introduction to this tool, we show that indicator-function power spectra are biased versions of the linear spectrum on large scales. We provide a calculation from first principles for this bias, we show that it reproduces simulation results, and we provide a simple functional form for the translinear portion of the indicator-function spectra. We also outline two applications: first, these spectra facilitate surgical excision of non-linearity and thus significantly increase the reach of linear theory. Second, indicator-function spectra permit calculation of theoretical covariance matrices for counts-in-cells, facilitating parameter estimation with complementary CIC methods.

Alfredo Mejia-Narvaez, Gustavo Bruzual, Sebastian F. Sanchez, Leticia Carigi, Jorge Barrera-Ballesteros, Mabel Valerdi, Renbin Yan, Niv Drory

27 pages, 17 figures, Submitted to ApJ

We introduce CoSHA: a Code for Stellar properties Heuristic Assignment. In order to estimate the stellar properties, CoSHA implements a Gradient Tree Boosting algorithm to label each star across the parameter space ($T_{\text{eff}}$, $\log{g}$, $\left[\text{Fe}/\text{H}\right]$, and $\left[\alpha/\text{Fe}\right]$). We use CoSHA to estimate these stellar atmospheric parameters of $22\,$k unique stars in the MaNGA Stellar Library (MaStar). To quantify the reliability of our approach, we run both internal tests using the G\"ottingen Stellar Library (GSL, a theoretical library) and the first data release of MaStar, and external tests by comparing the resulting distributions in the parameter space with the APOGEE estimates of the same properties. In summary, our parameter estimates span in the ranges: $T_{\text{eff}}=[2900,12000]\,$K, $\log{g}=[-0.5,5.6]$, $\left[\text{Fe}/\text{H}\right]=[-3.74,0.81]$, $\left[\alpha/\text{Fe}\right]=[-0.22,1.17]$. We report internal (external) uncertainties of the properties of $\sigma_{T_{\text{eff}}}\sim48\,(325)\,$K, $\sigma_{\log{g}}\sim0.2\,(0.4)$, $\sigma_{\left[\text{Fe}/\text{H}\right]}\sim0.13\,(0.27)$, $\sigma_{\left[\alpha/\text{Fe}\right]}\sim0.09\,(0.14)$. These uncertainties are comparable to those of other methods with similar objectives. Despite the fact that CoSHA is not aware of the spatial distribution of these physical properties in the Milky Way, we are able to recover the main trends known in the literature with great statistical confidence. The catalogue of physical properties can be accessed in \url{this http URL}.

Zorawar Wadiasingh, Christian J. T. van der Merwe, Christo Venter, Alice K. Harding, Matthew G. Baring

Presented at the 37th International Cosmic Ray Conference (ICRC 2021)

Pulsar winds interacting with sources of external pressure are well-established as efficient and prolific TeV accelerators in our Galaxy. Yet, enabled by observations from Fermi-LAT, a growing class of non-accreting pulsars in binaries has emerged and these are likely to become apparent as TeV emitters in the CTA era. This class consists of the black widows and redbacks, binaries in which a millisecond pulsar interacts with its low-mass companion. In such systems, an intrabinary shock can form as a site of particle acceleration and associated nonthermal emission. We motivate why these sources are particularly interesting for understanding pulsar winds. We also describe our new multizone code which models the X-ray and gamma-ray synchrotron and inverse Compton spectral components for select spider binaries, including diffusion, convection, and radiative energy losses in an axially symmetric, steady-state approach. This new multizone code simultaneously yields energy-dependent light curves and orbital-phase-resolved spectra. It also better constrains the multiplicity of electron/positron pairs that have been accelerated up to TeV energies and are necessary to power orbitally-modulated synchrotron emission components between the X-rays and MeV/GeV bands potentially observed in some systems. This affords a more robust prediction of the expected high-energy and VHE gamma-ray flux. Nearby MSPs with hot or flaring companions may be promising targets for CTA, and it is possible that spider binaries could contribute to the observed AMS-02 energetic positron excess.

K. B. Schmidt, J. Kerutt, L. Wisotzki, T. Urrutia, A. Feltre, M. V. Maseda, T. Nanayakkara, R. Bacon, L. A. Boogaard, S. Conseil, T. Contini, E. C. Herenz, W. Kollatschny, M. Krumpe, F. Leclercq, G. Mahler, J. Matthee, V. Mauerhofer, J. Richard, J. Schaye

Accepted by A&A. Contains 54 pages, 21 figures, and appendices. Catalogs described in Tables B.1 and C.1 will be made available with published paper

[Abbreviated] Rest-frame UV emission lines probe physical parameters of the emitting star-forming galaxies and their environments. The strongest main UV line, Ly$\alpha$, has been instrumental in advancing the general knowledge of galaxy formation in the early universe. However, observing Ly$\alpha$ emission becomes increasingly challenging at $z \gtrsim 6$ when the neutral hydrogen fraction of the CGM and IGM increases. Secondary weaker UV emission lines provide important alternative methods for studying galaxy properties at high redshift. We present a large sample of rest-frame UV emission line sources at intermediate redshift for calibrating and exploring the connection between secondary UV lines and the emitting galaxies' physical properties and their Ly$\alpha$ emission. The sample of 2052 emission line sources with $1.5 < z < 6.4$ was selected through untargeted source detection in three-dimensional MUSE data cubes. We searched optimally extracted 1D spectra of the full sample for UV emission features via emission line template matching, resulting in a sample of more than 100 rest-frame UV emission line detections. We show that the detection efficiency of (non-Ly$\alpha$) UV emission lines increases with survey depth, and that the UV emission line strength often correlate with the strength of Ciii]. We measured the velocity offsets of resonant emission lines with respect to systemic tracers as well as the electron density and the gas-phase abundance. Lastly, using "PhotoIonization Model Probability Density Functions" we find that the UV line emitters generally have ionization parameter log10(U) $\approx$ -2.5 and metal mass fractions that scatter around Z $\approx$ 10$^{-2}$, that is Z $\approx$ 0.66Z$\odot$. Value-added catalogs of the full sample of MUSE objects studied in this work and a collection of UV line emitters from the literature are provided with this paper.

Loukas Vlahos, Heinz Isliker, Nikos Sioulas

Parker (1983) suggested a mechanism for the formation of current sheets (CSs) in the solar atmosphere. His main idea was that the tangling of coronal magnetic field lines by photospheric random flows facilitates the continuous formation of CSs in the solar atmosphere. This part of his idea represents one of the many ways by which the turbulent convection zone drives the formation of coherent structures and CSs in the solar atmosphere. Other mechanisms include emerging magnetic flux, interaction of current filaments, and explosive magnetic structures. However, there are two unproven assumptions in the initial idea of Parker for the coronal heating through nanoflares that must be re-examined. They are related to his suggestion that {ALL CSs formed are led to magnetic reconnection and that magnetic reconnection heats the plasma in the solar atmosphere. Let us discuss these two assumptions briefly in this short comment: (1) Are ALL coherent structures and CSs formed by the turbulent convection zone reconnecting? Does turbulence associated with non-reconnecting CSs play a role in the heating of the corona? (2) Does magnetic reconnection heat the plasma?

A. Sinha, V. Baghmanyan, G. Peron, Y. Gallant, S. Casanova, M. Holler, A. Mitchell (for the H.E.S.S. collaboration)

Contribution to the 37th ICRC, Berlin 2021

Cosmic Ray (CR) interactions with the dense gas inside Giant Molecular Clouds (GMCs) produce neutral pions, which in turn decay into gamma rays. Thus, the gamma ray emission from GMCs is a direct tracer of the cosmic ray density and the matter density inside the clouds. Detection of enhanced TeV emission from GMCs, i.e., an emission significantly larger than what is expected from the average Galactic cosmic rays illuminating the cloud, can imply a variation in the local cosmic ray density, due to, for example, the presence of a recent accelerator in proximity to the cloud. Such gamma-ray observations can be crucial in probing the cosmic ray distribution across our Galaxy, but are complicated to perform with present generation Imaging Atmospheric Cherenkov Telescopes (IACTs). These studies require differentiating between the strong cosmic-ray induced background, the large scale diffuse emission, and the emission from the clouds, which is difficult to the small field of view of present generation IACTs. In this contribution, we use H.E.S.S. data collected over 16 years to search for TeV emission from GMCs in the inner molecular galacto-centric ring of our Galaxy. We implement a 3D FoV likelihood technique, and simultaneously model the hadronic background, the galactic diffuse emission and the emission expected from known VHE sources to probe for excess TeV gamma ray emission from GMCs.

J.-M. Grießmeier, D. A. Smith, G. Theureau, T. J. Johnson, M. Kerr, L. Bondonneau, I. Cognard, M. Serylak

17 pages, 3 figures, 1 appendix. Accepted by A&A

The Fermi Large Area Telescope has detected over 260 gamma-ray pulsars. About one quarter of these are labeled as radio-quiet. In the population of nonrecycled gamma-ray pulsars, the fraction of radio-quiet pulsars is higher, about one half. Most radio observations of gamma-ray pulsars have been performed at frequencies between 300 MHz and 2 GHz. However, pulsar radio fluxes increase rapidly with decreasing frequency, and their radio beams often broaden at low frequencies. As a consequence, some of these pulsars might be detectable at low radio frequencies even when no radio flux is detected above 300 MHz. Our aim is to test this hypothesis with low-frequency radio observations. We have observed 27 Fermi-discovered gamma-ray pulsars with the international LOw Frequency ARray (LOFAR) station FR606 in single-station mode. We used the LOFAR high band antenna (HBA) band (110-190 MHz). On average, we use 9 h of observation per target after the removal of affected datasets, resulting in a sensitivity for pulse-averaged flux on the order of 1-10 mJy. We do not detect radio pulsations from any of the 27 sources, and we establish stringent upper limits on their low-frequency radio fluxes. These nondetections are compatible with the upper limits derived from radio observations at other frequencies. We also determine the pulsars' geometry from the gamma-ray profiles to see for which pulsars the low-frequency radio beam is expected to cross Earth. This set of observations provides the most constraining upper limits on the flux density at 150 MHz for 27 radio-quiet gamma-ray pulsars. In spite of the beam-widening expected at low radio frequencies, most of our nondetections can be explained by an unfavorable viewing geometry; for the remaining observations, especially those of pulsars detected at higher frequencies, the nondetection is compatible with insufficient sensitivity.

Yuichi Oyama

2 pages, 1 figure, Prepared for submission to Physical Review Letters (Comment)

The seven seconds' gap in the Kamiokande-II SN1987A neutrino data is reexamined.

Ricardo López-Valdivia, Kimberly R. Sokal, Gregory N. Mace, Benjamin T. Kidder, Maryam Hussaini, Larissa Nofi, L. Prato, Christopher M. Johns-Krull, Heeyoung Oh, Jae-Joon Lee, Chan Park, Jae Sok Oh, Adam Kraus, Kyle F. Kaplan, Joe Llama, Andrew W. Mann, Hwihyun Kim, Michael A. Gully-Santiago, Hye-In Lee, Soojong Pak, Narae Hwang, Daniel T. Jaffe

21 pages, 10 figures, 5 tables. Accepted for publication in ApJ

We present fundamental parameters for 110 canonical K- & M-type (1.3$-$0.13$M_\odot$) Taurus-Auriga young stellar objects (YSOs). The analysis produces a simultaneous determination of effective temperature ($T_{\rm eff}$), surface gravity ($\log$ g), magnetic field strength (B), and projected rotational velocity ($v \sin i$). Our method employed synthetic spectra and high-resolution (R$\sim$45,000) near-infrared spectra taken with the Immersion GRating INfrared Spectrometer (IGRINS) to fit specific K-band spectral regions most sensitive to those parameters. The use of these high-resolution spectra reduces the influence of distance uncertainties, reddening, and non-photospheric continuum emission on the parameter determinations. The median total (fit + systematic) uncertainties were 170 K, 0.28 dex, 0.60 kG, 2.5 km s$^{-1}$ for $T_{\rm eff}$, $\log$ g, B, and $v \sin i$, respectively. We determined B for 41 Taurus YSOs (upper limits for the remainder) and find systematic offsets (lower $T_{\rm eff}$, higher $\log$ g and $v \sin i$) in parameters when B is measurable but not considered in the fit. The average $\log$ g for the Class II and Class III objects differs by 0.23$\pm$0.05dex, which is consistent with Class III objects being the more evolved members of the star-forming region. However, the dispersion in $\log$ g is greater than the uncertainties, which highlights how the YSO classification correlates with age ($\log$ g), yet there are exceptionally young (lower $\log$ g) Class III YSOs and relatively old (higher $\log$ g) Class II YSOs with unexplained evolutionary histories. The spectra from this work are provided in an online repository along with TW Hydrae Association (TWA) comparison objects and the model grid used in our analysis.

Shang-Min Tsai, Matej Malik, Daniel Kitzmann, James R. Lyons, Alexander Fateev, Elspeth Lee, Kevin Heng

50 pages, 41 figures. Submitted to AAS journals

We present an update of the open-source photochemical kinetics code VULCAN (Tsai et al. 2017; https://github.com/exoclime/VULCAN) to include C-H-N-O-S networks and photochemistry. Additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross-sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an intercomparison of HD 189733b models between Moses et al. (2011), Venot et al. (2012), and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep troposphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and ammonia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP-33b, HD 189733b, GJ 436b, and 51 Eridani b. We look into the effects of the C/O ratio and chemistry of titanium/vanadium species for WASP-33b; we revisit HD 189733b for the effects of sulfur and carbon condensation; the effects of internal heating and vertical mixing ($K_{\textrm{zz}}$) are explored for GJ 436b; we test updated planetary properties for 51 Eridani b with S$_8$ condensates. We find sulfur can couple to carbon or nitrogen and impact other species such as hydrogen, methane, and ammonia. The observable features of the synthetic spectra and trends in the photochemical haze precursors are discussed for each case.

One of the clearest but unresolved questions for Europa is the thickness of its icy shell. Europa's surface is resplendent with geological features that bear on this question, and ultimately on its interior, geological history, and astrobiological potential. We characterize the size and topographic expression of circular and subcircular features created by endogenic thermal and tectonic disturbances on Europa: pits, uplifts, and small, subcircular chaos. We utilize the medium-resolution Galileo regional maps (RegMaps), as well as high-resolution regions, digital elevation models derived from albedo-controlled photoclinometry, and in some cases stereo-controlled photoclinometry. While limited in extent, the high-resolution images are extremely valuable for detecting smaller features and for overall geomorphological analysis. A peak in the size-distribution for all features is found at ~5-6 km in diameter and no pits smaller than 3.3 km in diameter were found in high resolution images. Additionally, there is a trend for larger pits to be deeper, and larger uplifts to be higher. Our data support a diapiric or intracrustal sill interpretation (as opposed to purely non-intrusive, melt-through models) and place a lower limit on ice shell thickness at the time of feature formation of 3-to-8 km, assuming isostasy and depending on the composition of the ice and underlying ocean.

Takashi J. Moriya, Robert M. Quimby, Brant E. Robertson

14 pages, 9 figures, 1 table, submitted to The Astrophysical Journal

Massive stars play critical roles for the reionization of the Universe. Individual massive stars at the reionization epoch (z > 6) are too faint to observe and quantify their contributions to reionization. Some massive stars, however, explode as superluminous supernovae (SLSNe) or pair-instability supernovae (PISNe) that are luminous enough to observe even at z > 6 and allow for the direct characterization of massive star properties at the reionization epoch. In addition, hypothetical long-sought-after PISNe are expected to be present preferentially at high redshifts, and their discovery will have a tremendous impact on our understanding of massive star evolution and the formation of stellar mass black holes. The near-infrared Wide Field Instrument on Nancy Grace Roman Space Telescope will excel at discovering such rare high-redshift supernovae. In this work, we investigate the best survey strategy to discover and identify SLSNe and PISNe at z > 6 with Roman. We show that the combination of the F158 and F213 filters can clearly separate both SLSNe and PISNe at z > 6 from nearby supernovae through their colors and magnitudes. The limiting magnitudes are required to be 27.0 mag and 26.5 mag in the F158 and F213 filters, respectively, to identify supernovae at z > 6. If we conduct a 10 deg2 transient survey with these limiting magnitudes for 5 years with a cadence of one year, we expect to discover 22.5 +- 2.8 PISNe and 3.1 +- 0.3 SLSNe at z > 6, depending on the cosmic star-formation history. The same survey is estimated to discover 76.1 +- 8.2 PISNe and 9.1 +- 0.9 SLSNe at 5 < z < 6. Such a supernova survey requires the total observational time of approximately 525 hours in 5 years. The legacy data acquired with the survey will also be beneficial for many different science cases including the study of high-redshift galaxies.

Gregory S. Vance (1), Kimberly L. Emig (2), Cecilia Lunardini (3), Rogier A. Windhorst (1) ((1) School of Earth and Space Exploration, Arizona State University, (2) National Radio Astronomy Observatory, (3) Department of Physics, Arizona State University)

26 pages, 4 figures, 7 tables; to be submitted to JCAP

Searches for spatial associations between high-energy neutrinos observed at the IceCube Neutrino Observatory and known astronomical objects may hold the key to establishing the neutrinos' origins and the origins of hadronic cosmic rays. While extragalactic sources like the blazar TXS 0506+056 merit significant attention, Galactic sources may also represent part of the puzzle. Here, we explore whether open clusters and supernova remnants in the Milky Way contribute measurably to the IceCube track-like neutrino events above 200 TeV. By searching for positional coincidences with catalogs of known astronomical objects, we can identify and investigate neutrino events whose origins are potentially Galactic. We use Monte Carlo randomization together with models of the Galactic plane in order to determine whether these coincidences are more likely to be causal associations or random chance. In all analyses presented, the number of coincidences detected was found to be consistent with the null hypothesis of chance coincidence. Our results imply that the combined contribution of Galactic open clusters and supernova remnants to the track-like neutrino events detected at IceCube is well under 30%. This upper limit is compatible with the results presented in other Galactic neutrino studies.

J. Shejeelammal, Aruna Goswami

Accepted for publication in the Astrophysical Journal (ApJ)

The observed abundance diversities among the CEMP stars can shed light on the formation and evolution of elements in the early Galaxy. In this work, we present results obtained from a detailed abundance analysis of a sample of seven extrinsic carbon stars. The analysis is based on high-resolution spectra obtained with HCT/HESP (R$\sim$60,000) and SUBARU/HDS (R$\sim$50,000). We present, for the first time, the elemental abundance results for the objects BD$-$19 132, BD$-$19 290, HE~1304$-$2111, HE~1354$-$2257, and BD+19 3109. Abundances of a few elements are available in literature for HE~1157$-$0518 and HD~202851, we present an update on the abundances of these elements along with new abundance estimates for several other elements. Our analysis confirms the object HD~202851 to be a CH star. While BD$-$19 132, HE~1354$-$2257, and BD+19 3109 are found to be CEMP-s stars, the objects BD$-$19 290, HE~1157$-$0518, and HE~1304$-$2111 are found to belong to CEMP-r/s group. The observed abundance patterns of the three CEMP-r/s stars are well reproduced with the i-process model predictions. While the objects BD+19 3109 and HD~202851 are confirmed binaries, the binary status of the remaining objects are not known. Analysis based on different elemental abundance ratios confirms low-mass former AGB companions for all the objects. Kinematic analysis shows that BD$-$19 290, HE~1157$-$0518, HE~1354$-$2257, and BD+19 3109 belong to the Galactic halo, whereas BD$-$19 132, HE~1304$-$2111, and HD~202851 are members of Galactic thin disk.

Recently the MAGIC telescope observed three TeV gamma-ray extended sources in the galactic plane in the neighborhood of radio SNR G24.7+0.6. Among them, the PWN HESS J1837-069 was detected earlier by the HESS observatory during its first galactic plane survey. The other two sources, MAGIC J1835-069 and MAGIC J1837-073 are detected for the first time at such high energies. Here we shall show that the observed gamma-rays from the SNR G24.7+0.6 and the HESS J1837-069 can be explained in terms of hadronic interactions of the PWN/SNR accelerated cosmic rays with the ambient matter. We shall further demonstrate that the observed gamma-rays from the MAGIC J1837$-$073 can be interpreted through hadronic interactions of runaway cosmic-rays from PWN HESS J1837-069 with the molecular cloud at the location of MAGIC J1837-073. No such association has been found between MAGIC J1835$-$069 and SNR G24.7+0.6 or PWN HESS J1837$-$069. We have examined the maximum energy attainable by cosmic-ray particles in the SNR G24.7+0.6/ PWN HESS J1837-069 and the possibility of their detection with future gamma-ray telescopes. The study of TeV neutrino emissions from the stated sources suggests that the HESS J1837$-$069 should be detected by IceCube Gen-2 neutrino telescope in a few years of observation.

Giovanni Morlino

Proceeding of the 37th International Cosmic Ray Conference (ICRC2021)

We investigate the process of particle acceleration at the termination shock that develops in the bubble excavated by winds of star clusters in the interstellar medium. We develop a theory of diffusive shock acceleration at such shock and we find that the maximum energy may reach the PeV region for very powerful clusters. We show how the maximum energy is limited by two different processes: the particle escape from the bubble boundary and the drop of energy gain for particles able to diffuse up to the center of the cluster. A crucial role in this problem is played by the dissipation of kinetic energy of the wind to magnetic perturbations which determines the diffusion regime of particles: in case the diffusion is close to Bohm than PeV energies can be reached.

Hyesung Kang

11 pages, 6 figures, submitted to JKAS

Radio relics could be generated by multiple shocks induced in the turbulent intracluster medium during galaxy mergers. Kang (2021) demonstrated that the re-acceleration of cosmic ray (CR) protons via diffusive shock acceleration (DSA) by multiple shocks could enhance the acceleration efficiency and flatten the CR spectrum, compared to a single episode of DSA. Here we examine the CR electron acceleration through multiple re-acceleration by considering energy losses and decompression of the particle distribution and magnetic fields in the postshock region between consecutive shock passages. We find that the accumulated effects of repeated re-acceleration are significant, if preceding shocks are stronger than the last shock and the shock passage interval is $\lesssim20$ Myr. In such cases, both the CR spectrum and the ensuing radiation spectrum behind the last shock are enhanced and become flatter than the canonical DSA power-law forms. As a result, the shock Mach number estimated from radio observations tends be higher than the actual Mach number of the last shock. Thus, multiple episodes of DSA may explain the enhanced acceleration efficiency for CR electrons and the discrepancy of shock Mach numbers, $M_{\rm X} \lesssim M_{\rm rad}$, inferred for some observed radio relics.

V. A. Allakhverdyan, A. D. Avrorin, A. V. Avrorin, V. M. Aynutdinov, R. Bannasch, Z. Bardačová, I. A. Belolaptikov, I. V. Borina, V. B. Brudanin, N. M. Budnev, V. Y. Dik, G. V. Domogatsky, A. A. Doroshenko, R. Dvornický, A. N. Dyachok, Zh.-A. M. Dzhilkibaev, E. Eckerová, T. V. Elzhov, L. Fajt, S. V. Fialkovski, A. R. Gafarov, K. V. Golubkov, N. S. Gorshkov, T. I. Gress, M. S. Katulin, K. G. Kebkal, O. G. Kebkal, E. V. Khramov, M. M. Kolbin, K. V. Konischev, K. A. Kopański, A. V. Korobchenko, A. P. Koshechkin, V. A. Kozhin, M. V. Kruglov, M. K. Kryukov, V. F. Kulepov, Pa. Malecki, Y. M. Malyshkin, M. B. Milenin, R. R. Mirgazov, D. V. Naumov, V. Nazari, W. Noga, D. P. Petukhov, E. N. Pliskovsky, M. I. Rozanov, V. D. Rushay, E. V. Ryabov, G. B. Safronov, B. A. Shaybonov, et al. (15 additional authors not shown)

Submitted to Proc. of the 37th International Cosmic Ray Conference (ICRC 2021), PoS-1144, July 12th -- 23rd, 2021, Online -- Berlin, Germany. 8 pages, 7 figures

Baikal-GVD is a neutrino telescope currently under construction in Lake Baikal. GVD is formed by multi-meganton subarrays (clusters). The design of Baikal-GVD allows one to search for astrophysical neutrinos already at early phases of the array construction. We present here preliminary results of a search for high-energy neutrinos with GVD in 2019-2020.

Chong Li, Keping Qiu, Bo Hu, Yue Cao

10 pages, 7 figures, Accepted for publication in ApJL

Using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we detect a giant HI filamentary structure in the sky region of 307$.\!\!^{\circ}$7 $<$ $\alpha$ $<$ 311$.\!\!^{\circ}$0 and 40$.\!\!^{\circ}$9 $<$ $\delta$ $<$ 43$.\!\!^{\circ}$4. The structure has a velocity range of $-$170 km s$^{-1}$ to $-$130 km s$^{-1}$, and a mean velocity of $-$150 km s$^{-1}$, putting it to a Galactocentric distance of 22 kpc. The HI structure has a length of 1.1 kpc, which appears to be so far the furthest and largest giant filament in the Galaxy and we name it Cattail. Its mass is calculated to be 6.5 $\times$ 10$^4$ M$_{\odot}$ and the linear mass density is 60 M$_{\odot}$ pc$^{-1}$. Its width is 207 pc, corresponding to an aspect ratio of 5:1. Cattail possesses a small velocity gradient (0.02 km s$^{-1}$ pc$^{-1}$) along its major axis. Together with the HI4PI data, we find that Cattail could have an even larger length, up to 5 kpc. We also identify another new elongated structure to be the extension into the Galactic first quadrant of the Outer Scutum-Centaurus (OSC) arm, and Cattail appears to be located far behind the OSC. The question about how such a huge filament is produced at the extreme Galactic location remains open. Alternatively, Cattail might be part of a new arm beyond the OSC, though it is puzzling that the structure does not fully follow the warp of the Galactic disk.

Hiroyuki Tashiro, Katsuya T. Abe, Teppei Minoda

6 pages, 3 figures

Baryonic gas falling onto a primordial black hole (PBH) emits photons via the free-free process. These photons can contribute the diffuse free-free background radiation in the frequency range of the cosmic microwave background radiation (CMB). We show that the intensity of the free-free background radiation from PBHs depends on the mass and abundance of PBHs. In particular, considering the growth of a dark matter (DM) halo around a PBH by non-PBH DM particles strongly enhances the free-free background radiation. This enhancement allows us to obtain the tight constraint on the PBH abundance with mass $M_{\rm PBH} > 10 M_\odot$. Comparison with the CMB free-free foreground component in high galactic latitude measured by Planck, our constraint on the PBH abundance fraction to the total DM, $f_{\rm{PBH}}$, is $f_{\rm{PBH}}<1.2\times10^{-5}(M_{\rm{PBH}}/100M_{\odot})^{-2.3}$. We also discuss the impact of the free-free radiation from PBHs on the measurement of 21-cm background radiation from high redshifts, including the EDGES anomaly.

Ekaterina Ilin, Katja Poppenhaeger, Sarah J. Schmidt, Silva P. Järvinen, Elisabeth R. Newton, Julián D. Alvarado-Gómez, J. Sebastian Pineda, James R. A. Davenport, Mahmoudreza Oshagh, Ilya Ilyin

23 pages, 18 figures, accepted to MNRAS, see conference poster this https URL for short summary

White-light flares are magnetically driven localized brightenings on the surfaces of stars. Their temporal, spectral, and statistical properties present a treasury of physical information about stellar magnetic fields. The spatial distributions of magnetic spots and associated flaring regions help constrain dynamo theories. Moreover, flares are thought to crucially affect the habitability of exoplanets that orbit these stars. Measuring the location of flares on stars other than the Sun is challenging due to the lack of spatial resolution. Here we present four fully convective stars observed with the Transiting Exoplanet Survey Satellite (TESS) that displayed large, long-duration flares in white-light which were modulated in brightness by the stars' fast rotation. This allowed us to determine the loci of these flares directly from the light curves. All four flares occurred at latitudes between 55 deg and 81 deg, far higher than typical solar flare latitudes. Our findings are evidence that strong magnetic fields tend to emerge close to the stellar rotational poles for fully convective stars, and suggest that the impact of flares on the habitability of exoplanets around small stars could be weaker than previously thought.

T. Saito, C. Delgado, O. Blanch, M. Artero, J. A. Barrio, F. Cassol, C. Diaz, D. Hadasch, D. Hoffmann, J. Houles, Y. Inome, M. Iori, L. Jouvin, D. Kerszberg, Y. Kobayashi, H. Kubo, G. Martinez, D. Mazin, E. Moretti, T. Nakamori, S. Nozaki, T. Oka, A. Okumura, M. Palatiello, M. Polo, J. Prast, S. Sakurai, Y. Sunada, M. Takahashi, R. Takeishi, L. A. Tejedor, T. Yamamoto, T. Yoshida (for the CTA LST project)

10 pages, 9 figures, Proceedings of the 37th International Cosmic Ray Conference (ICRC 2021), Berlin, Germany

The first Large Size Telescope (LST-1) of the Cherenkov Telescope Array has been operational since October 2018 at La Palma, Spain. We report on the results obtained during the camera commissioning. The noise level of the readout is determined as a 0.2 p.e. level. The gain of PMTs are well equalized within 2\% variation, using the calibration flash system. The effect of the night sky background on the signal readout noise as well as the PMT gain estimation are also well evaluated. Trigger thresholds are optimized for the lowest possible gamma-ray energy threshold and the trigger distribution synchronization has been achieved within 1~ns precision. Automatic rate control realizes the stable observation with 1.5\% rate variation over 3 hours. The performance of the novel DAQ system demonstrates a less than 10\% dead time for 15 kHz trigger rate even with sophisticated online data correction.

Galaxy redshift surveys are designed to map cosmic structures in three dimensions for large-scale structure studies. Nevertheless, limitations due to sampling and the survey window are unavoidable and degrade the cosmological constraints. We present an analysis of the VIMOS Public Extragalactic Redshift Survey (VIPERS) over the redshift range $0.6 < z < 1$ that is optimised to extract the cosmological parameters while fully accounting for the complex survey geometry. We employ the Gibbs sampling algorithm to iteratively draw samples of the galaxy density field in redshift space, the galaxy bias, the matter density, baryon fraction and growth-rate parameter $f\sigma_8$ based on a multivariate Gaussian likelihood and prior on the density field. Despite the high number of degrees of freedom, the samples converge to the joint posterior distribution and give self-consistent constraints on the model parameters. We validate the approach using VIPERS mock galaxy catalogues. Although the uncertainty is underestimated by the Gaussian likelihood on the scales that we consider by $50\%$, the dispersion of the results from the mock catalogues gives a robust error estimate. We find that the precision of the results matches those of the traditional analyses applied to the VIPERS data that use more constrained models. By relaxing the model assumptions, we confirm that the data deliver consistent constraints on the $\Lambda$CDM model. This work provides a case-study for the application of maximum-likelihood analyses for the next generation of galaxy redshift surveys.

Snehanshu Maiti, Kirit Makwana, Heshou Zhang, Huirong Yan

10 pages, 11 figures, submitted to APJ

This paper studies cosmic ray (CR) transport in magneto hydrodynamic (MHD) turbulence. CR transport is strongly dependent on the properties of the magnetic turbulence. We perform test particle simulations to study the interactions of CR with both total MHD turbulence and decomposed MHD modes. The spatial diffusion coefficients and the pitch angle scattering diffusion coefficients are calculated from the test particle trajectories in turbulence. Our results confirm that the fast modes dominate the CR propagation, whereas Alfv\'en and slow modes are much less efficient and have shown similar pitch angle scattering rates. We investigate the cross field transport on large and small scales. On large/global scales, normal diffusion is observed and the diffusion coefficient is suppressed by $M_A^\zeta$ compared to the parallel diffusion coefficients, with $\zeta$ closer to 4 in Alfv\'en modes than that in total turbulence as theoretically expected. For the CR transport on scales smaller than the turbulence injection scale, both the local and global magnetic reference frames are adopted. Super diffusion is observed on such small scales in all the cases. Particularly, CR transport in Alfv\'en modes show clear Richardson diffusion in the local reference frame. Our results have broad applications to CRs in various astrophysical environments.

W. Bednarek

19 pages, 5 figures, accepted to MNRAS

A part of early type stars is characterised by strong dipole magnetic field that is modified by the outflow of dense wind from the stellar surface. At some distance from the surface (above the Alfven radius), the wind drives the magnetic field into the reconnection in the equatorial region of the dipole magnetic field. We propose that electrons accelerated in these reconnection regions can be responsible for efficient comptonization of stellar radiation producing gamma-ray emission. We investigate the propagation of electrons in the equatorial region of the magnetosphere by including their advection with the equatorial wind. The synchrotron and IC spectra are calculated assuming that a significant part of the wind energy is transferred to relativistic electrons. As an example, the parameters of luminous, strongly magnetized star HD 37022 ($\Theta^1$ Ori C) are considered. The IC gamma-ray emission is predicted to be detected either in the GeV energy range by the Fermi-LAT telescope or in the sub-TeV energies by the Cherenkov Telescope Array. However, since the stellar winds are often time variable and the magnetic axis can be inclined to the rotational axis of the star, the gamma-ray emission is expected to show variability with the rotational period of the star and, on a longer time scale, with the stellar circle of the magnetic activity. Those features might serve as tests of the proposed scenario for gamma-ray emission from single, luminous stars.

S. Ustamujic, S. Orlando, M. Miceli, F. Bocchino, M. Limongi, A. Chieffi, C. Trigilio, G. Umana, F. Bufano, A. Ingallinera, G. Peres

18 pages, 12 figures; accepted for publication in Astronomy & Astrophysics

LBVs are massive evolved stars that suffer sporadic and violent mass-loss events. They have been proposed as the progenitors of some core-collapse SNe, but this idea is still debated due to the lack of direct evidence. Since SNRs can carry in their morphology the fingerprints of the progenitor stars as well as of the inhomogeneous CSM sculpted by the progenitors, the study of SNRs from LBVs could help to place core-collapse SNe in context with the evolution of massive stars. We investigate the physical, chemical and morphological properties of the remnants of SNe originating from LBVs, in order to search for signatures, revealing the nature of the progenitors, in the ejecta distribution and morphology of the remnants. As a template of LBVs, we considered the actual LBV candidate Gal 026.47+0.02. We selected a grid of models, which describe the evolution of a massive star with properties consistent with those of Gal 026.47+0.02 and its final fate as core-collapse SN. We developed a 3D HD model that follows the post-explosion evolution of the ejecta from the breakout of the shock wave at the stellar surface to the interaction of the SNR with a CSM characterized by two dense nested toroidal shells, parametrized in agreement with multi-wavelength observations of Gal 026.47+0.02. Our models show a strong interaction of the blast wave with the CSM which determines an important slowdown of the expansion of the ejecta in the equatorial plane where the two shells lay, determining a high degree of asymmetry in the remnant. After 10000 years of evolution the ejecta show an elongated shape forming a broad jet-like structure caused by the interaction with the shells and oriented along the axis of the toroidal shells.

H. Abdalla, F. Aharonian, F. Ait Benkhali, E.O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun, P. Brun, M. Bryan, M. Büchele, T. Bulik, T. Bylund, S. Caroff, A. Carosi, S. Casanova, T. Chand, S. Chandra, A. Chen, G. Cotter, M. Curyło, J. Damascene Mbarubucyeye, I.D. Davids, J. Davies, C. Deil, J. Devin, L. Dirson, A. Djannati-Atai, A. Dmytriiev, A. Donath, V. Doroshenko, L. Dreyer, C. Duffy, J. Dyks, K. Egberts, F. Eichhorn, S. Einecke, G. Emery, J.-P. Ernenwein, K. Feijen, S. Fegan, A. Fiasson, G. Fichet de Clairfontaine, G. Fontaine, et al. (176 additional authors not shown)

Accepted for publication in Astronomy & Astrophysics

We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for modeling the spectral energy distribution, an analysis of Fermi-LAT Pass 8 data was also included. We unambiguously detect N132D at very high energies (VHE) with a significance of 5.7 sigma. We report the results of a detailed analysis of its spectrum and localization based on the extended H.E.S.S. data set. The joint analysis of the extended H.E.S.S and Fermi-LAT data results in a spectral energy distribution in the energy range from 1.7 GeV to 14.8 TeV, which suggests a high luminosity of N132D at GeV and TeV energies. We set a lower limit on a gamma-ray cutoff energy of 8 TeV with a confidence level of 95%. The new gamma-ray spectrum as well as multiwavelength observations of N132D when compared to physical models suggests a hadronic origin of the VHE gamma-ray emission. SNR N132D is a VHE gamma-ray source that shows a spectrum extending to the VHE domain without a spectral cutoff at a few TeV, unlike the younger oxygen-rich SNR Cassiopeia A. The gamma-ray properties of N132D may be affected by an interaction with a nearby molecular cloud that partially lies inside the 95% confidence region of the source position. [Abridged]

R. F. Wimmer-Schweingruber, N. Janitzek, D. Pacheco, I. Cernuda, F. Espinosa Lara, R. Gómez-Herrero, G. M. Mason, R. C. Allen, Z. G. Xu, F. Carcaboso, A. Kollhoff, P. Kühl, J. L. Freiherr von Forstner, L. Berger, J. Rodriguez-Pacheco, G. C. Ho, G. B. Andrews, V. Angelini, A. Aran, S. Boden, S. I. Böttcher, A. Carrasco, N. Dresing, S. Eldrum, R. Elftmann, V. Evans, O. Gevin, J. Hayes, B. Heber, T. S. Horbury, S. R. Kulkarni, D. Lario, W. J. Lees, O. Limousin, O. Malandraki, C. Martín, H. O'Brien, M. Prieto Mateo, A. Ravanbakhsh, O. Rodriguez Polo, S. Sánchez Prieto, C. E. Schlemm, H. Seifert, J. C. Terasa, K. Tyagi, R. Vainio, A. Walsh, M. K. Yedla

Solar Orbiter strives to unveil how the Sun controls and shapes the heliosphere and fills it with energetic particle radiation. To this end, its Energetic Particle Detector (EPD) has now been in operation, providing excellent data, for just over a year. EPD measures suprathermal and energetic particles in the energy range from a few keV up to (near-) relativistic energies (few MeV for electrons and about 500 MeV/nuc for ions). We present an overview of the initial results from the first year of operations and we provide a first assessment of issues and limitations. During this first year of operations of the Solar Orbiter mission, EPD has recorded several particle events at distances between 0.5 and 1 au from the Sun. We present dynamic and time-averaged energy spectra for ions that were measured with a combination of all four EPD sensors, namely: the SupraThermal Electron and Proton sensor (STEP), the Electron Proton Telescope (EPT), the Suprathermal Ion Spectrograph (SIS), and the High-Energy Telescope (HET) as well as the associated energy spectra for electrons measured with STEP and EPT. We illustrate the capabilities of the EPD suite using the 10-11 December 2020 solar particle event. This event showed an enrichment of heavy ions as well as $^3$He, for which we also present dynamic spectra measured with SIS. The high anisotropy of electrons at the onset of the event and its temporal evolution is also shown using data from these sensors. We discuss the ongoing in-flight calibration and a few open instrumental issues using data from the 21 July and the 10-11 December 2020 events and give guidelines and examples for the usage of the EPD data. We explain how spacecraft operations may affect EPD data and we present a list of such time periods in the appendix. A list of the most significant particle enhancements as observed by EPT during this first year is also provided.

Alison Mitchell, Gavin Rowell, Silvia Celli, Sabrina Einecke

8 pages, 5 figures, Proceedings of the 37th International Cosmic Ray Conference (ICRC2021)

Interstellar clouds can act as target material for hadronic cosmic rays; gamma-rays produced through inelastic proton-proton collisions and spatially associated with the clouds provide a key indicator of efficient particle acceleration. However, even for PeVatron sources reaching PeV energies, the system of cloud and accelerator must fulfil several conditions in order to produce a detectable gamma-ray flux. In this contribution, we characterise the necessary properties of both cloud and accelerator. Using available Supernova Remnant (SNR) and interstellar cloud catalogues, and assuming particle acceleration to PeV energies in a nearby SNR, we produce a ranked shortlist of the most promising target systems; those for which a detectable gamma-ray flux is predicted. We discuss detection prospects for future facilities including CTA and SWGO; and compare our predictions with known gamma-ray sources, including the Ultra-High-Energy sources recently detected by LHAASO. A range of model scenarios are tested, including variation in the diffusion coefficient and particle spectrum, under which the best candidate clouds in our shortlist are consistently bright. On average, a detectable gamma-ray flux is more likely for more massive clouds; for systems with lower separation distance between the SNR and cloud; and for slightly older SNRs, due to the time required for particles to traverse the separation distance.

Manuel Meyer, James Davies, Julian Kuhlmann

Proceedings of the 37th International Cosmic Ray Conference (ICRC 2021). The online documentation is available at this https URL

Axions and axion-like particles (ALPs) are hypothetical particles that occur in extensions of the Standard Model and are candidates for cold dark matter. They could be detected through their oscillations into photons in the presence of external electromagnetic fields. gammaALPs is an open-source python framework that computes the oscillation probability between photons and axions/ALPs. In addition to solving the photon-ALP equations of motion, gammaALPs includes models for magnetic fields in different astrophysical environments such as jets of active galactic nuclei, intra-cluster and intergalactic media, and the Milky Way. Users are also able to easily incorporate their own custom magnetic-field models. We review the basic functionality and features of gammaALPs, which is heavily based on other open-source scientific packages such as numpy and scipy. Although focused on gamma-ray energies, gammaALPs can be easily extended to arbitrary photon energies.

Karen L. Masters, Coleman Krawczyk, Shoaib Shamsi, Alexander Todd, Daniel Finnegan, Matthew Bershady, Kevin Bundy, Brian Cherinka, Amelia Fraser-McKelvie, Dhanesh Krishnarao, Sandor Kruk, Richard R. Lane, David Law, Chris Lintott, Michael Merrifield, Brooke Simmons, Anne-Marie Weijmans, Renbin Yan

13 pages, 9 figures. MNRAS accepted

The challenge of consistent identification of internal structure in galaxies - in particular disc galaxy components like spiral arms, bars, and bulges - has hindered our ability to study the physical impact of such structure across large samples. In this paper we present Galaxy Zoo: 3D (GZ: 3D) a crowdsourcing project built on the Zooniverse platform which we used to create spatial pixel (spaxel) maps that identify galaxy centres, foreground stars, galactic bars and spiral arms for 29831 galaxies which were potential targets of the MaNGA survey (Mapping Nearby Galaxies at Apache Point Observatory, part of the fourth phase of the Sloan Digital Sky Surveys or SDSS-IV), including nearly all of the 10,010 galaxies ultimately observed. Our crowd-sourced visual identification of asymmetric, internal structures provides valuable insight on the evolutionary role of non-axisymmetric processes that is otherwise lost when MaNGA data cubes are azimuthally averaged. We present the publicly available GZ:3D catalog alongside validation tests and example use cases. These data may in the future provide a useful training set for automated identification of spiral arm features. As an illustration, we use the spiral masks in a sample of 825 galaxies to measure the enhancement of star formation spatially linked to spiral arms, which we measure to be a factor of three over the background disc, and how this enhancement increases with radius.

Matej Petković, Luke Lucas, Tomaž Stepišnik, Panče Panov, Nikola Simidjievski, Dragi Kocev

Presented at the SMC-IT 2021 conference

The Mars Express (MEX) spacecraft has been orbiting Mars since 2004. The operators need to constantly monitor its behavior and handle sporadic deviations (outliers) from the expected patterns of measurements of quantities that the satellite is sending to Earth. In this paper, we analyze the patterns of the electrical power consumption of MEX's thermal subsystem, that maintains the spacecraft's temperature at the desired level. The consumption is not constant, but should be roughly periodic in the short term, with the period that corresponds to one orbit around Mars. By using long short-term memory neural networks, we show that the consumption pattern is more irregular than expected, and successfully detect such irregularities, opening possibility for automatic outlier detection on MEX in the future.

Manuel Meyer, Tanja Petrushevska (for the Fermi-LAT collaboration)

Proceedings of the 37th International Cosmic Ray Conference (ICRC 2021)

During a core-collapse supernova (SN), axion-like particles (ALPs) could be produced through the Primakoff process and subsequently convert into gamma rays in the magnetic field of the Milky Way. Using a sample of well studied extragalactic SNe at optical wavelengths, we estimate the time of the core collapse and search for a coincident gamma-ray burst with the Fermi Large Area Telescope (LAT). Under the assumption that at least one SN was contained within the LAT field of view, we exclude photon-ALP couplings within a factor of $\sim$5 of previous limits from SN1987A. With the increasing number of SNe observed with optical surveys, our results demonstrate the potential to probe ALP dark matter with combined optical and gamma-ray observations. We also provide preliminary results for the estimation of explosion times of 15 close-by SNe observed recently with ZTF. Our findings show that the explosion time can be estimated within one day (statistical uncertainty only) making them promising targets for a follow-up LAT analysis.

A. Gokus, A. Kreikenbohm, K. Leiter, T. Bretz, T. Dauser, D. Dorner, D. Elsaesser, F. Eppel, J. Hessdoerfer, M. Kadler, A. Kraus, M. Kreter, I. Kreykenbohm, M. Langejahn, K. Mannheim, P. Thalhammer, J. Wilms, A. Arbet-Engels, D. Baack, M. Balbo, A. Biland, J. Buss, L. Eisenberger, D. Hildebrand, R. Iotov, A. Kalenski, A. Mitchell, D. Neise, M. Noethe, A. Paravac, W. Rhode, B. Schleicher, V. Sliusar, R. Walter

8 pages, 4 figures; Proceedings of the 37th International Cosmic Ray Conference (ICRC2021). 12-23 July, 2021. Berlin, Germany. Online at this https URL

The blazar Mrk 421 shows frequent, short flares in the TeV energy regime. Due to the fast nature of such episodes, we often fail to obtain sufficient simultaneous information about flux variations in several energy bands. To overcome this lack of multi-wavelength (MWL) coverage, especially for the pre- and post-flare periods, we have set up a monitoring program with the FACT telescope (TeV energies) and the Neil Gehrels Swift Observatory (X-rays). On 2019 June 9, Mrk 421 showed a TeV outburst reaching a flux level of more than two times the flux of the Crab Nebula at TeV energies. We acquired simultaneous data in the X-rays with additional observations by XMM-Newton and INTEGRAL. For the first time, we can study a TeV blazar in outburst taking advantage of highly sensitive X-ray data from XMM-Newton and INTEGRAL combined. Our dataset is complemented by pointed radio observations by Effelsberg at GHz frequencies. We present our first results, including the {\gamma}-ray and X-ray light curves, a timing analysis of the X-ray data obtained with XMM-Newton , as well as the radio spectra before, during and after the flare.

Tirthendu Sinha, Saurabh Sharma, Neelam Panwar, N. Matsunaga, K. Ogura, N. Kobayashi, R. K. Yadav, A. Ghosh, R. Pandey, P. S. Bisht

Accepted for publication in ApJ

We present the results from our time-series imaging data taken with the 1.3m Devasthal fast optical telescope and 0.81m Tenagara telescope in $V$, $R_{c}$, $I_{c}$ bands covering an area of $\sim18^\prime.4\times 18^\prime.4$ towards the star-forming region Sh 2-190. This photometric data helped us to explore the nature of the variability of pre-main sequence (PMS) stars. We have identified 85 PMS variables, i.e., 37 Class II and 48 Class III sources. Forty-five of the PMS variables are showing periodicity in their light curves. We show that the stars with thicker discs and envelopes rotate slower and exhibit larger photometric variations compared to their disc-less counterparts. This result suggests that rotation of the PMS stars is regulated by the presence of circumstellar discs. We also found that the period of the stars show a decreasing trend with increasing mass in the range of $\sim$0.5-2.5 M$_\odot$. Our result indicates that most of the variability in Class II sources is ascribed to the presence of thick disc, while the presence of cool spots on the stellar surface causes the brightness variation in Class III sources. X-ray activities in the PMS stars were found to be at the saturation level reported for the main sequence (MS) stars. The younger counterparts of the PMS variables are showing less X-ray activity hinting towards a less significant role of a stellar disc in X-ray generation.

Yingkang Zhang, Tao An, Sándor Frey, Xiaolong Yang, Máté Krezinger, Oleg Titov, Alexey Melnikov, Pablo de Vicente, Fengchun Shu, Ailing Wang

11 pages, 5 figures, accepted by MNRAS

Jets of high-redshift active galactic nuclei (AGNs) can be used to directly probe the activity of the black holes in the early Universe. Radio sources with jets misaligned with respect to the line of sight are expected to dominate the high-redshift AGN population. In this paper, we present the high-resolution imaging results of a z=4.57 AGN J2102+6015 by analyzing its multi-epoch dual-frequency very long baseline interferometry (VLBI) data. The 8.4-GHz VLBI images reveal two major features along the east-west direction separated by $\sim$ 10 milli-arcsec (mas). From the spectral index map, both features show flat/inverted spectra. The separation between the two features remains almost unchanged over an observation period of $\sim$ 13 years, placing an upper limit of the separation speed as about 0.04 mas year$^{-1}$. Previous studies have classified the source as a GHz-peaked spectrum quasar. Our results indicate that J2102+6015 is most likely a young, compact symmetric object rather than a blazar-type core-jet source.

Gy.M. Szabó, D. Gandolfi, A. Brandeker, Sz. Csizmadia, Z. Garai, N. Billot, C. Broeg, D. Ehrenreich, A. Fortier, L. Fossati, S. Hoyer, L. Kiss, A. Lecavelier des Etangs, P.F.L. Maxted, I. Ribas, Y. Alibert, R. Alonso, G. Anglada Escudé, T. Bárczy, S.C.C. Barros, D. Barrado, W. Baumjohann, M. Beck, T. Beck, A. Bekkelien, X. Bonfils, W. Benz, L. Borsato, M-D. Busch, J. Cabrera, S. Charnoz, A. Collier Cameron, C. Corral Van Damme, M.B. Davies, L. Delrez, M. Deleuil, O.D.S. Demangeon, B.-O. Demory, A. Erikson, M. Fridlund, D. Futyan, A. García Muñoz, M. Gillon, M. Guedel, P. Guterman, K. Heng, K.G. Isaak, G. Lacedelli, J. Laskar, M. Lend, C. Lovis, A. Luntzer, D. Magrin, V. Nascimbeni, G. Olofsson, H.P. Osborn, R. Ottensamer, I. Pagano, E. Pallé, G. Peter, D. Piazza, et al. (20 additional authors not shown)

16 pages, 11 figures, 7 tables; accepted by Astronomy and Astrophysics

AU Mic is a young planetary system with a resolved debris disc showing signs of planet formation and two transiting warm Neptunes near mean-motion resonances. Here we analyse three transits of AU Mic b observed with the CHaracterising ExOPlanet Satellite (CHEOPS), supplemented with sector 1 and 27 Transiting Exoplanet Survey Satellite (TESS) photometry, and the All-Sky Automated Survey (ASAS) from the ground. The refined orbital period of AU Mic b is 8.462995 \pm 0.000003 d, whereas the stellar rotational period is P_{rot}=4.8367 \pm 0.0006 d. The two periods indicate a 7:4 spin--orbit commensurability at a precision of 0.1%. Therefore, all transits are observed in front of one of the four possible stellar central longitudes. This is strongly supported by the observation that the same complex star-spot pattern is seen in the second and third CHEOPS visits that were separated by four orbits (and seven stellar rotations). Using a bootstrap analysis we find that flares and star spots reduce the accuracy of transit parameters by up to 10% in the planet-to-star radius ratio and the accuracy on transit time by 3-4 minutes. Nevertheless, occulted stellar spot features independently confirm the presence of transit timing variations (TTVs) with an amplitude of at least 4 minutes. We find that the outer companion, AU Mic c may cause the observed TTVs.

We show that a recent constraint on the cosmic birefringence effect due to dark energy can be related to the constraints on the coupling of axion dark matter to photon, by relying on a simple model of two-axion alignment mechanism with periodic potentials. Owing to the alignment of the potentials, one linear combination of two fields provides a nearly flat direction and acts as dark energy, whereas the other combination provides a steep direction and acts as dark matter. This scenario solves the known conceptual issues of one-field model for dark energy and predicts the connection between seemingly disparate constraints on the dark sectors of our universe.

Souvik Bose, Luc Rouppe van der Voort, Jayant Joshi, Vasco M.J. Henriques, Daniel Nóbrega-Siverio, Juan Martínez-Sykora, Bart De Pontieu

20 pages, 10 figures, accepted for publication in Astronomy and Astrophysics. Abstract abridged due to arXiv requirements

Spectroscopic observations of the emission lines formed in the solar transition region (TR) commonly show persistent downflows of the order of 10--15 km/s. The cause of such downflows, however, is still not fully clear and has remained a matter of debate. We aim to understand the cause of such downflows by studying the coronal and TR responses to the recently reported chromospheric downflowing rapid red shifted excursions (RREs), and their impact on heating the solar atmosphere. We have used two sets of coordinated data from SST, IRIS, and SDO for analyzing the response of the downflowing RREs in the TR and corona. To provide theoretical support, we use an already existing 2.5D MHD simulation of spicules performed with the Bifrost code. We find ample occurrences of downflowing RREs and show several examples of their spatio-temporal evolution, sampling multiple wavelength channels ranging from the cooler chromospheric to hotter coronal channels. These downflowing features are thought to be likely associated with the returning components of the previously heated spicular plasma. Furthermore, the TR Doppler shifts associated with them are close to the average red shifts observed in this region which further implies that these flows could (partly) be responsible for the persistent downflows observed in the TR. We also propose two mechanisms (a typical upflow followed by a downflow and downflows along a loop), from the perspective of numerical simulation, that could explain the ubiquitous occurrence of such downflows. A detailed comparison between the synthetic and observed spectral characteristics, reveals a distinctive match, and further suggests an impact on the heating of the solar atmosphere. We present evidence that suggests that at least some of the downflowing RREs are the chromospheric counterparts of the TR and lower coronal downflows.

Geetanjali Sarkar, Raghvendra Sahai

89 pages, 104 figures, submitted to The Astrophysical Journal (ApJ)

Common Envelope (CE) systems are the result of Roche lobe overflow in interacting binaries. The subsequent evolution of the CE, its ejection and the formation of dust in its ejecta while the primary is on the Red Giant Branch, gives rise to a recently-identified evolutionary class -- dusty post-RGB stars. Their spectral energy distributions (SEDs) suggest that their mass-ejecta are similar to dusty post-Asymptotic Giant Branch (post-AGB) stars. We have modeled the SEDs of a select sample of post-RGB and post-AGB stars in the Large Magellanic Cloud (LMC), quantified the total dust mass (and gas mass assuming gas-to-dust ratio) in the disks and shells and set constraints on the dust grain compositions and sizes. We find that the shell masses in the post-RGBs are generally less than those in post-AGBs, with the caveat that a substantial amount of mass in both types of objects may lie in cold, extended shells. Our models suggest that circumstellar disks, when present, are geometrically thick, flared structures with a substantial opening angle, consistent with numerical simulations of CE evolution (CEE). Comparison of our model dust mass values with the predictions of dust production during CEE on the RGB suggest that CEE occurred near or at the top of the RGB for the post-RGB sources in our sample. A surprising result is that the ejected dust in some post-RGB sources appears to be carbon-rich, providing independent support for the hypothesis of binary interaction leading to the formation of dusty post-RGB objects.

Saeid Foroughi-Abari, Felix Kling, Yu-Dai Tsai

9 pages, 3 figures. Para uma FORMOSA aventura!

We identify potentially the world's most sensitive location to search for millicharged particles in the 10 MeV to 100 GeV mass range: the forward region at the LHC. We propose constructing a scintillator-based experiment, FORward MicrOcharge SeArch (FORMOSA) in this location, and estimate the corresponding sensitivity projection. We show that FORMOSA can discover millicharged particles in a large and unexplored parameter space, and study strongly interacting dark matter that cannot be detected by ground-based direct-detection experiments. The newly proposed LHC Forward Physics Facility (FPF) provides an ideal structure to host the full FORMOSA experiment.

Marco Chianese, Damiano F. G. Fiorillo, Rasmi Hajjar, Gennaro Miele, Ninetta Saviano

8 pages, 5 figures

Among the several strategies for indirect searches of dark matter, one very promising one is to look for the gamma-rays from decaying dark matter. Here we use the most up-to-date upper bounds on the gamma-ray flux from $10^5$ to $10^{11}$ GeV, obtained from CASA-MIA, KASCADE, KASCADE-Grande, Pierre Auger Observatory, and Telescope Array. We obtain global limits on dark matter lifetime in the range of masses $m_\mathrm{DM}=[10^7-10^{15}]~\mathrm{GeV}$. We provide the bounds for a set of decay channels chosen as representatives. The constraints derived here are new and cover a region of the parameter space not yet explored. We compare our results with the projected constraints from future neutrino telescopes, in order to quantify the improvement that will be obtained by the complementary high-energy neutrino searches.

Juan Barranco, Javier Chagoya, Alberto Diez-Tejedor, Gustavo Niz, Armando A. Roque

49+1 pages, 17 figures, 2 tables

We establish the existence of time-dependent solitons in a modified gravity framework, which is defined by the low energy limit of theories with a weakly broken galileon symmetry and a mass term. These are regular vacuum configurations of finite energy characterized by a single continuous parameter representing the amplitude of the scalar degree of freedom at the origin. When the central field amplitude is small the objects are indistinguishable from boson stars. In contrast, increasing the central value of the amplitude triggers the effect of higher derivative operators in the effective theory, leading to departures from the previous solutions, until the theory becomes strongly coupled and model-dependent. The higher order operators are part of the (beyond) Horndeski theory, hence the name of the compact objects. Moreover, a remnant of the galileon non-renormalization theorem guarantees that the existence and properties of these solutions are not affected by quantum corrections. Finally, we discuss the linear stability under small radial perturbations, the mass-radius relation, the compactness, the appearance of innermost stable circular orbits and photon spheres, and some astrophysical signatures (accretion disks, gravitational radiation and lensing) that may be relevant to falsify the model.

Jim Madsen, Laurent Mulot, Christian Spiering (for the IceCube Collaboration)

Presented at the 37th International Cosmic Ray Conference (ICRC 2021). See arXiv:2107.06966 for all IceCube contributions

We present an overview of the art project Aganta Kairos (To Fish the Metaphysical Time). This project celebrates the neutrino, the ghost particle, which scientists consider a cosmic messenger and the artist regards as a link between people who care about their relationship to the cosmos and question their origins. The artwork is based on a performance of celebration and seeks to build a human community that encompasses different knowledge domains and interpretations of the universe. This intersection of knowledge is realized during the performance of placing a plaque, held with witnesses, and during subsequent exhibitions. Images, sounds, videos, and sculpture testify to the diversity of approaches to questioning our origins, ranging from traditional western science to ancient shamanism. The sites were selected for their global coverage and, for the South Pole, Mediterranean, and Lake Baikal, their connection to ongoing neutrino experiments. In December 2020, a plaque was installed at the South Pole IceCube Laboratory, the seventh and final site. We provide examples of images and links to additional images and videos.

What is the chance we start a conversation with another civilization like our own? Our technological society produced signals that could be received by other extraterrestrial civilizations, within a sphere around us with a radius of $\sim 10^2$ light years. Given that, the Copernican principle provides a lower limit on the response time that we should expect from transmitters on Earth-like planets around Sun-like stars. If our civilization lives longer, the expected number of responses could increase. We explore the chance of detecting a response in the future, and show that a response should only be expected to arrive after a few millennia.

A view on the history and current status of dark matter and dark energy, at a fairly introductory level.

Lorenzo Di Pietro, Victor Gorbenko, Shota Komatsu

93 pages, 22 figures

We study the fundamentals of quantum field theory on a rigid de Sitter space. We show that the perturbative expansion of late-time correlation functions to all orders can be equivalently generated by a non-unitary Lagrangian on a Euclidean AdS geometry. This finding simplifies dramatically perturbative computations, as well as allows us to establish basic properties of these correlators, which comprise a Euclidean CFT. We use this to infer the analytic structure of the spectral density that captures the conformal partial wave expansion of a late-time four-point function, to derive an OPE expansion, and to constrain the operator spectrum. Generically, dimensions and OPE coefficients do not obey the usual CFT notion of unitarity. Instead, unitarity of the de Sitter theory manifests itself as the positivity of the spectral density. This statement does not rely on the use of Euclidean AdS Lagrangians and holds non-perturbatively. We illustrate and check these properties by explicit calculations in a scalar theory by computing first tree-level, and then full one-loop-resummed exchange diagrams. An exchanged particle appears as a resonant feature in the spectral density which can be potentially useful in experimental searches.

A foundational model has been developed based on trends built from empirical data of space exploration and computing power through the first six plus decades of the Space Age which projects earliest possible launch dates for human-crewed missions from cis-lunar space to selected Solar System and interstellar destinations. The model uses computational power, expressed as transistors per microprocessor, as a key broadly limiting factor for deep space missions' reach and complexity. The goal of this analysis is to provide a projected timeframe for humanity to become a multi-world species through off-world colonization, and in so doing all but guarantees the long-term survival of the human race from natural and human-caused calamities that could befall life on Earth. Be-ginning with the development and deployment of the first nuclear weapons near the end of World War II, humanity entered a 'Window of Peril' which will not be safely closed until robust off-world colonies become a reality. Our findings suggest the first human-crewed missions to land on Mars, selected Asteroid Belt objects, and selected moons of Jupiter and Saturn can occur before the end of the 21st century. Launches of human-crewed interstellar missions to exoplanet destinations within roughly 40 lightyears of the Solar System are seen as possible during the 23rd century and launch of intragalactic missions by the end of the 24th century. An aggressive and sustained space exploration program, which includes colonization, is thus seen as critical to the long-term survival of the human race.

Black hole spectroscopy is the proposal to observe multiple quasinormal modes in the ringdown of a binary black hole merger. In addition to the fundamental quadrupolar mode, overtones and higher harmonics may be present and detectable in the gravitational wave signal, allowing for tests of the no-hair theorem. We analyze in detail the strengths and weaknesses of the standard Rayleigh criterion supplied with a Fisher matrix error estimation, and we find that the criterion is useful, but too restrictive. Therefore we motivate the use of a conservative high Bayes factor threshold to obtain the black hole spectroscopy horizons of current and future detectors, i.e., the distance (averaged in sky location and binary inclination) up to which one or more additional modes can be detected and confidently distinguished from each other. We set up all of our searches for additional modes starting at $t = 10(M_1+M_2)$ after the peak amplitude in simulated signals of circular nonspinning binaries. An agnostic multimode analysis allows us to rank the subdominant modes: for nearly equal mass binaries we find $(\ell, m, n) = (2,2,1)$ and $(3,3,0)$ and, for very asymmetric binaries, $(3,3,0)$ and $(4,4,0)$, for the secondary and tertiary modes, respectively. At the current estimated rates for heavy stellar mass binary black hole mergers, with primary masses between 45 and 100 solar masses, we expect an event rate of mergers within the $(2,2,1)$ spectroscopy horizon of $0.03 - 0.10\ {\rm yr}^{-1}$ for LIGO at design sensitivity and $(0.6 - 2.4) \times 10^3\ {\rm yr}^{-1}$ for the future third generation ground-based detector Cosmic Explorer.

Yuuki Wada, Takahiro Matsumoto, Teruaki Enoto, Kazuhiro Nakazawa, Takayuki Yuasa, Yoshihiro Furuta, Daisuke Yonetoku, Tatsuya Sawano, Go Okada, Hidehito Nanto, Shohei Hisadomi, Yuna Tsuji, Gabriel Sousa Diniz, Kazuo Makishima, Harufumi Tsuchiya

33 pages, 31 figures, submitted to Physical Review Research

In 2015 the Gamma-Ray Observation of Winter Thunderstorms (GROWTH) collaboration launched a mapping observation campaign for high-energy atmospheric phenomena related to thunderstorms and lightning discharges. This campaign has developed a detection network of gamma rays with up to 10 radiation monitors installed in Kanazawa and Komatsu cities, Ishikawa Prefecture, Japan, where low-charge-center winter thunderstorms frequently occur. During four winter seasons from October 2016 to April 2020, in total 70 gamma-ray glows, minute-lasting bursts of gamma rays originating from thunderclouds, were detected. Their average duration is 58.9 sec. Among the detected events, 77% were observed in nighttime. The gamma-ray glows can be classified into temporally-symmetric, temporally-asymmetric, and lightning-terminated types based on their count-rate histories. An averaged energy spectrum of the gamma-ray glows is well fitted with a power-law function with an exponential cutoff, whose photon index, cutoff energy, and flux are $0.613\pm0.009$, $4.68\pm0.04$ MeV, and $(1.013\pm0.003)\times10^{-5}$ erg cm$^{-2}$ s$^{-1}$ (0.2-20.0MeV), respectively. The present paper provides the first catalog of gamma-ray glows and their statistical analysis detected during winter thunderstorms in the Kanazawa and Komatsu areas.

A. Biekert, C. Chang, C. W. Fink, M. Garcia-Sciveres, E. C. Glazer, W. Guo, S. A. Hertel, S. Kravitz, J. Lin, M. Lisovenko, R. Mahapatra, D. N. McKinsey, J. S. Nguyen, V. Novosad, W. Page, P. Patel, B. Penning, H. D. Pinckney, M. Pyle, R. K. Romani, A. S. Seilnacht, A. Serafin, R. J. Smith, P. Sorensen, B. Suerfu, A. Suzuki, V. Velan, G. Wang, S. L. Watkins, V. G. Yefremenko, L. Yuan, J. Zhang

16 pages, 18 figures

Superfluid $^4$He is a promising target material for direct detection of light ($<$ 1 GeV) dark matter. Possible signal channels available for readout in this medium include prompt photons, triplet excimers, and roton and phonon quasiparticles. The relative yield of these signals has implications for the sensitivity and discrimination power of a superfluid $^4$He dark matter detector. Using a 16 cm$^3$ volume of 1.75 K superfluid $^4$He read out by six immersed photomultiplier tubes, we measured the scintillation from electronic recoils ranging between 36.3 and 185 keV$_\mathrm{ee}$, yielding a mean signal size of $1.12^{+0.02}_{-0.03}$ phe/keV$_\mathrm{ee}$, and nuclear recoils from 53.2 to 1090 keV$_\mathrm{nr}$. We compare the results of our relative scintillation yield measurements to an existing semi-empirical model based on helium-helium and electron-helium interaction cross sections. We also study the behavior of delayed scintillation components as a function of recoil type and energy, a further avenue for signal discrimination in superfluid $^4$He.