Abstracts of Interest

Selected by:


Abstract: 2203.13143
Full Text: [ PostScript, PDF]

Title:New determination of the production cross section for secondary positrons and electrons in the Galaxy

Download PDF
Abstract: The cosmic-ray fluxes of electrons and positrons ($e^{\pm}$) are measured with high precision by the space-borne particle spectrometer AMS-02. To infer a precise interpretation of the production processes for $e^{\pm}$ in our Galaxy, it is necessary to have an accurate description of the secondary component, produced by the interaction of cosmic-ray proton and helium with the interstellar medium atoms. We determine new analytical functions of the Lorentz invariant cross section for the production of $\pi^\pm$ and $K^\pm$ by fitting data from collider experiments. We also evaluate the invariant cross sections for several other channels, involving for example hyperon decays, contributing at the few \% level on the total cross section. For all these particles, the relevant 2 and 3 body decay channels are implemented, with the polarized $\mu^\pm$ decay computed with next-to-leading order corrections. The cross section for scattering of nuclei heavier than protons is modeled by fitting data on $p+C$ collisions. The total differential cross section $d\sigma/dT_{e^\pm}(p+p\rightarrow e^\pm+X)$ is predicted from 10 MeV up to 10 TeV of $e^\pm$ energy with an uncertainty of about 5-7\% in the energies relevant for AMS-02 positron flux, thus dramatically reducing the precision of the theoretical model with respect to the state of the art. Finally, we provide a prediction for the secondary Galactic $e^\pm$ source spectrum with an uncertainty of the same level. As a service for the scientific community, we provide numerical tables and a script to calculate energy-differential cross sections.

Comments: 25 pages, 19 figures. Submitted to PRD. The tables of the energy differential cross sections of electrons and positrons with a script to read them can be found here: this https URL. Results are under peer-review process


Abstract: 2203.12866
Full Text: [ PostScript, PDF]

Title:COSE$ν$: A Collective Oscillation Simulation Engine for Neutrinos

Download PDF
Abstract: We introduce the implementation details of the simulation code \cosenu, which numerically solves a set of non-linear partial differential equations that govern the dynamics of neutrino collective flavor conversions. We systematically provide the details of both the finite difference method supported by Kreiss-Oliger dissipation and the finite volume method with seventh order weighted essentially non-oscillatory scheme. To ensure the reliability of the code, we perform the comparison of the simulation results with theoretically obtainable solutions. In order to understand and characterize the error accumulation behavior of the implementations when neutrino self-interactions are switched on, we also analyze the evolution of the deviation of the conserved quantities for different values of simulation parameters. We report the performance of our code with both CPUs and GPUs. The public version of the \cosenu~package is available at \url{this https URL}.

Comments: 17 pages, 13 figures


Abstract: 2203.12696
Full Text: [ PostScript, PDF]

Title:Rocks, Water and Noble Liquids: Unfolding the Flavor Contents of Supernova Neutrinos

Download PDF
Abstract: Measuring core-collapse supernova neutrinos, both from individual supernovae within the Milky Way and from past core collapses throughout the Universe (the diffuse supernova neutrino background, or DSNB), is one of the main goals of current and next generation neutrino experiments. Detecting the heavy-lepton flavor (muon and tau types, collectively $\nu_x$) component of the flux is particularly challenging due to small statistics and large backgrounds. While the next galactic neutrino burst will be observed in a plethora of neutrino channels, allowing to measure a small number of $\nu_x$ events, only upper limits are anticipated for the diffuse $\nu_x$ flux even after decades of data taking with conventional detectors. However, paleo-detectors could measure the time-integrated flux of neutrinos from galactic core-collapse supernovae via flavor-blind neutral current interactions. In this work, we show how combining a measurement of the average galactic core-collapse supernova flux with paleo detectors and measurements of the DSNB electron-type neutrino fluxes with the next-generation water Cherenkov detector Hyper-Kamiokande and the liquid noble gas detector DUNE will allow to determine the mean supernova $\nu_x$ flux parameters with precision of order ten percent.

Comments: 18 pages, 5 figures


Abstract: 2203.11936
Full Text: [ PostScript, PDF]

Title:High-Energy Extragalactic Neutrino Astrophysics

Download PDF
Abstract: The detection of an astrophysical flux of neutrinos in the TeV-PeV energy range by the IceCube observatory has opened new possibilities for the study of extreme cosmic accelerators. The apparent isotropy of the neutrino arrival directions favors an extragalactic origin for this flux, potentially created by a large population of distant sources. Recent evidence for the detection of neutrino emission from extragalactic sources include the active galaxies TXS 0506+056 and NGC 1068. We here review the current status of the search for the sources of the high-energy neutrino flux, concentrating on its extragalactic contribution. We discuss the implications of these observations for multimessenger studies of cosmic sources and present an outlook for how additional observations by current and future instruments will help address fundamental questions in the emerging field of high-energy neutrino astronomy.

Comments: 24 pages, 5 figures, 1 table, invited review article accepted for publication in ARNPS


Abstract: 2203.13111
Full Text: [ PostScript, PDF]

Title:Unusually high HCO+/CO ratios in and outside supernova remnant W49B

Download PDF
Abstract: Galactic supernova remnants (SNRs) and their environments provide the nearest laboratories to study SN feedback. We performed molecular observations toward SNR W49B -- the most luminous Galactic SNR in the X-ray band, aiming to explore signs of multiple feedback channels of SNRs on nearby molecular clouds (MCs). We found very broad HCO+ lines with widths of dv = 48--75 km/s in the SNR southwest, providing strong evidence that W49B is perturbing MCs at a systemic velocity of $V_{LSR}=61$--65 km/s, and placing W49B at a distance of $7.9\pm 0.6$ kpc. We observed unusually high intensity ratios of HCO+ J=1-0/CO J=1-0 not only at shocked regions ($1.1\pm 0.4$ and $0.70\pm 0.16$), but also in quiescent clouds over 1 pc away from the SNR's eastern boundary (> 0.2). By comparing with the magnetohydrodynamics (MHD) shock models, we interpret that the high ratio in the broad-line regions can result from a cosmic-ray (CR) induced chemistry in shocked MCs, where the CR ionization rate is enhanced to around 10--100 times of the Galactic level. The high HCO+/CO ratio outside the SNR is probably caused by the radiation precursor, while the luminous X-ray emission of W49B can explain a few properties in this region. The above results provide observational evidence that SNR can strongly influence the molecular chemistry in and outside the shock boundary via their shocks, CRs, and radiation. We propose that the HCO+/CO ratio is a potentially useful tool to probe SNR's multi-channel influence on MCs.

Comments: 19 pages, 11 figures. Submitted to ApJ. Comments welcome


Abstract: 2203.12747
Full Text: [ PostScript, PDF]

Title:Seven Years of SN 2014C: a Multi-Wavelength Synthesis of an Extraordinary Supernova

Download PDF
Abstract: SN 2014C was originally classified as a Type Ib supernova, but at phase {\phi} = 127 d post-explosion strong H{\alpha} emission was observed. SN 2014C has since been observed in radio, infrared, optical and X-ray bands. Here we present new optical spectroscopic and photometric data spanning {\phi} = 947 - 2494 d post-explosion. We address the evolution of the broadened H{\alpha} emission line, as well as broad [O III] emission and other lines. We also conduct a parallel analysis of all publicly available multi-wavelength data. From our spectra, we find a nearly constant H{\alpha} FWHM velocity width of {\sim}2000 km/s that is significantly lower than that of other broadened atomic transitions ({\sim}3000 - 7000 km/s) present in our spectra ([O I] {\lambda}6300; [O III] {\lambda}{\lambda}4959,5007; He I {\lambda}7065; [Ca II] {\lambda}{\lambda}7291,7324). The late radio data demand a fast forward shock ({\sim}10,000 km/s at {\phi} = 1700 d) in rarified matter that contrasts with the modest velocity of the H{\alpha}. We propose that the infrared flux originates from a toroidal-like structure of hydrogen surrounding the progenitor system, while later emission at other wavelengths (radio, X-ray) likely originates predominantly from the reverse shock in the ejecta and the forward shock in the quasi-spherical progenitor He wind. We propose that the H{\alpha} emission arises in the boundary layer between the ejecta and torus. We also consider the possible roles of a pulsar and a binary companion.

Comments: 33 pages, 12 figures


Abstract: 2203.11502
Full Text: [ PostScript, PDF]

Title:Phenomenological modelling of the Crab Nebula's broad-band energy spectrum and its apparent extension

Authors:Ludmilla Dirson (1 and 2), Dieter Horns (1) ((1) Institut für Experimentalphysik, Universität Hamburg (2) Université de Strasbourg, CNRS, Observatoire astronomique de Strasbourg)
Download PDF
Abstract: The Crab Nebula emits bright non-thermal radiation from radio to the most energetic photons. The underlying physical model of a relativistic wind from the pulsar terminating in a hydrodynamic standing shock remains unchanged since the early 1970s. In this model, an increase of the toroidal magnetic field downstream from the shock is expected. We introduce a detailed radiative model to calculate non-thermal synchrotron and inverse Compton as well as thermal dust emission self-consistently to compare quantitatively with observational data. Special care is given to the radial dependence of electron and seed field density. The radiative model is used to estimate the parameters of electrons and dust in the nebula. A combined fit based upon a $\chi^2$ minimisation reproduces successfully the complete data set used. For the best-fitting model, the energy density of the magnetic field dominates over the particle energy density up to a distance of $\approx 1.3~r_s$ ($r_s$: distance of the termination shock from the pulsar). The very high energy (VHE: $E>100$ GeV) gamma-ray spectra set the strongest constraints on the radial dependence of the magnetic field: $B(r)=(264\pm9)~\mu\mathrm{G} (r/r_s)^{-0.51\pm0.03}$. The reconstructed magnetic field and its radial dependence indicates a ratio of Poynting to kinetic energy flux $\sigma\approx 0.1$ at the termination shock, $\approx 30$ times larger than estimated up to now. Consequently, the confinement of the nebula would require additional mechanisms to slow down the flow through e.g. excitation of small-scale turbulence with possible dissipation of magnetic field.

Comments: submitted to A&A, abbr. abstract


Abstract: 2203.12908
Full Text: [ PostScript, PDF]

Title:The Optical Two and Three-Dimensional Fundamental Plane Correlations for More than 180 Gamma-Ray Burst Afterglows with Swift/UVOT, RATIR, and the SUBARU Telescope

Download PDF
Abstract: Gamma-ray bursts (GRBs) are fascinating events due to their panchromatic nature. We study optical plateaus in GRB afterglows via an extended search into archival data. We comprehensively analyze all published GRBs with known redshifts and optical plateaus observed by many ground-based telescopes (e.g., Subaru Telescope, RATIR) around the world and several space-based observatories such as the Neil Gehrels Swift Observatory. We fit 502 optical light curves (LCs), showing the existence of the plateau in 181 cases. This sample is 77\% larger than the previous one \citep{dainotti2021b}, and it is the largest compilation so far of optical plateaus. We discover the 3D fundamental plane relation at optical wavelengths using this sample. This correlation is between the rest-frame time at the end of the plateau emission, $T^{*}_{\rm opt}$, its optical luminosity, $L_{\rm opt}$, and the peak in the optical prompt emission, $L_{\rm peak, opt}$, thus resembling the three-dimensional (3D) X-ray fundamental plane relation \citep{dainotti2016}. We correct our sample for redshift evolution and selection effects, discovering that this correlation is indeed intrinsic to GRB physics. We investigate the rest-frame end time distributions in X-rays and optical ($T^{*}_{\rm opt}$, $T^{*}_{\rm X}$), and conclude that the plateau is achromatic only when selection biases are not considered. We also investigate if the 3D optical correlation may be a new discriminant between optical GRB classes and {\bf find that there is no significant separation between the classes compared to the Gold sample plane after correcting for evolution.}

Comments: 24 pages, 7 figures and 3 table. Submitted to ApJ. The part in bold is the addition as a response to the referee queries


Abstract: 2203.12551
Full Text: [ PostScript, PDF]

Title:A Possible Subclassification of Fast Radio Bursts

Download PDF
Abstract: Although fast radio bursts (FRBs) have been an active field in astronomy and cosmology, their origin is still unknown to date. One of the interesting topics is the classification of FRBs, which is closely related to the origin of FRBs. Different physical mechanisms are required by different classes of FRBs. In the literature, they usually could be classified into non-repeating and repeating FRBs. Well motivated by the observations, here we are interested in the possible subclassification of FRBs. By using the first CHIME/FRB catalog, we propose to subclassify non-repeating (type I) FRBs into type Ia and Ib FRBs. The distribution of type Ia FRBs is delayed with respect to the cosmic star formation history (SFH), and hence they are probably associated with old stellar populations, while the distribution of type Ib FRBs tracks SFH, and hence they are probably associated with young stellar populations. Accordingly, the physical criteria for this subclassification of type I FRBs have been clearly determined. We find that there are some tight empirical correlations for type Ia FRBs but not for type Ib FRBs, and vice versa. These make them different in physical properties. Similarly, we suggest that repeating (type II) FRBs could also be subclassified into type IIa and IIb FRBs. A universal subclassification scheme is given at the end. This subclassification of FRBs might help us to reveal quite different physical mechanisms behind them, and improve their applications in astronomy and cosmology.

Comments: 17 pages, 4 tables, 11 figures, revtex4


Abstract: 2203.11787
Full Text: [ PostScript, PDF]

Title:In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

Download PDF
Abstract: The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.



Abstract: 2203.11464
Full Text: [ PostScript, PDF]

Title:Luminosity Selection for Gamma Ray Burst

Download PDF
Abstract: There exists an inevitable scatter in intrinsic luminosity of Gamma Ray Bursts(GRBs). If there is relativistic beaming in the source, viewing angle variation necessarily introduces variation in the intrinsic luminosity function (ILF). Scatter in the ILF can cause a selection bias where distant sources that are detected have a larger median luminosity than those detected close by. Median luminosity divides any given population into equal halves. When the functional form of a distribution is unknown, it can be a more robust diagnostic than any that use trial functional forms. In this work we employ a statistical test based on median luminosity and apply it to test a class of models for GRBs. We assume that the GRB jet has a finite opening angle and that the orientation of the GRB jet is random relative to the observer.We calculate $L_{median}$ as a function of redshift by simulating GRBs empirically, theoretically and use the luminosity vs redshift {\it Swift} data in order to compare the theoretical results with the observed ones.
The method accounts for the fact that at some redshifts there may be some GRBs that go undetected. We find that $L_{median}$ is extremely insensitive to the on-axis (i.e. maximal) luminosity of the jet.

Comments: 9 pages, 9 figures, 3 tables, accepted for publication in A&A. arXiv admin note: substantial text overlap with arXiv:2010.04810


Abstract: 2203.10669
Full Text: [ PostScript, PDF]

Title:MeerKAT uncovers the physics of an Odd Radio Circle

Download PDF
Abstract: Odd Radio Circles (ORCs) are recently-discovered faint diffuse circles of radio emission, of unknown cause, surrounding galaxies at moderate redshift ($z ~ 0.2-0.6). Here we present detailed new MeerKAT radio images at 1284 MHz of the first ORC, originally discovered with the Australian Square Kilometre Array Pathfinder, with higher resolution (6 arcsec) and sensitivity (~ 2.4 uJy/bm).
In addition to the new images, which reveal a complex internal structure consisting of multiple arcs, we also present polarisation and spectral index maps. Based on these new data, we consider potential mechanisms that may generate the ORCs.

Comments: Accepted by MNRAS


Abstract: 2203.11626
Full Text: [ PostScript, PDF]

Title:Assessing the Impact of Hydrogen Absorption on the Characteristics of the Galactic Center Excess

Download PDF
Abstract: We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation-transport modelling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic Center. The main benefits of the new gas model are, a), the ability to reproduce the negative line signals seen with the H$I$4PI survey and, b), the accounting for gas that primarily manifests itself through absorption.
We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter. We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H$I$ model is used to estimate the cosmic-ray induced diffuse gamma-ray emission. The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence for a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favoured by the data. This finding is robust under the variation of various parameters, for example the excitation temperature of atomic hydrogen, and a number of tests for systematic issues.

Comments: ApJ in press, 19 pages


Abstract: 2203.11223
Full Text: [ PostScript, PDF]

Title:Star Formation in the Central Molecular Zone of the Milky Way

Download PDF
Abstract: The Central Molecular Zone (CMZ) is a ring-like accumulation of molecular gas in the innermost few hundred parsecs of the Milky Way, generated by the inward transport of matter driven by the Galactic bar. The CMZ is the most extreme star-forming environment in the Galaxy. The unique combination of large-scale dynamics and extreme interstellar medium conditions, characterised by high densities, temperatures, pressures, turbulent motions, and strong magnetic fields, make the CMZ an ideal region for testing current star and planet formation theories. We review the recent observational and theoretical advances in the field, and combine these to draw a comprehensive, multi-scale and multi-physics picture of the cycle of matter and energy in the context of star formation in the closest galactic nucleus.

Comments: To appear in Protostars and Planets VII; Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura; Comments welcome


Abstract: 2203.13405
Full Text: [ PostScript, PDF]

Title:Centrally Concentrated H I Distribution Enhances Star Formation in Galaxies

Download PDF
Abstract: We use a sample of 13,511 nearby galaxies from the ALFALFA and SDSS spectroscopic surveys to study the relation between the spatial distribution of H I 21 cm emission and star formation rate (SFR). We introduce a new non-parametric quantity $K$, measured from the curve-of-growth of the line, to describe the shape of the integrated H I profile. The value of $K$ increases from double-horned to single-peaked profiles, depending on projection effects and the spatial and velocity distribution of the gas. Using carefully chosen samples to control for the competing factors that influence the integrated line profile, we argue that useful inferences can be made on the spatial distribution of the gas. We find that galaxies with a high value of $K$ tend to have more centrally concentrated H I distribution within the optical disk of the galaxy at fixed conditions, and that larger values of $K$ are associated with higher levels of total and central SFR. The results suggest that the global concentration of H I plays an important role in facilitating the conversion of neutral atomic hydrogen to molecular hydrogen gas, which, in turn, affects the star formation activity throughout the optical disk. Our sample is biased against quiescent galaxies, and thus the conclusions may not hold for galaxies with low SFR or low H I content.

Comments: Accepted by ApJ


Abstract: 2203.13404
Full Text: [ PostScript, PDF]

Title:Statistical Analysis of H I Profile Asymmetry and Shape for Nearby Galaxies

Download PDF
Abstract: We present a uniform analysis of the integrated profile of the H I emission line of 29,958 nearby ($z < 0.06$) galaxies extracted from the ALFALFA 21 cm survey. We apply the curve-of-growth technique to derive a database of spectral parameters and robust estimates of their associated uncertainties. Besides the central velocity and total flux, the main catalog provides new measures of line width, profile asymmetry, and profile shape. For a subsample of 13,511 galaxies with optical properties available from the Sloan Digital Sky Survey, we compute inclination angle-corrected line widths, rotation velocities empirically calibrated from spatially resolved observations, and dynamical masses based on H I sizes estimated from the H I mass. To facilitate subsequent scientific applications of the database, we also compile a number of ancillary physical properties of the galaxies, including their optical morphology, stellar mass, and various diagnostics of star formation activity. We use the homogeneous catalog of H I parameters to examine the statistical properties of profile asymmetry and shape. Across the full sample, which covers a wide range of stellar masses and environments, statistically significant H I profile asymmetry is detected in $\sim 20\%$ of the galaxy population. The global H I profiles are $35.2 \pm 0.3\%$ single-peaked, $26.9 \pm 0.3\%$ flat-topped, and $37.9 \pm 0.3\%$ double-horned. At a given inclination angle, double-horned profiles are preferentially associated with galaxies of higher stellar mass or optical concentration, while galaxies of lower mass or concentration tend to have single-peaked profiles.

Comments: Accepted by ApJS


Abstract: 2203.10842
Full Text: [ PostScript, PDF]

Title:Radio detections of IR-selected runaway stellar bow shocks

Download PDF
Abstract: Massive stars moving at supersonic peculiar velocities through the interstellar medium (ISM) can create bow shocks, arc-like structures at the interface between the stellar wind and the ISM. Many such bow shocks have been detected and catalogued at IR wavelengths, but detections in other wavebands remain rare. Strikingly, while electrons are expected to be accelerated in the bow shock and their non-thermal emission may include synchrotron emission at low frequencies, only two massive runaway stellar bow shocks have to date been detected in the radio band. Here, we examine a sample of fifty IR-detected bow shocks from the E-BOSS catalogues in recently released radio images from the Rapid ASKAP Continuum Survey (RACS). We identify three confident and three likely counterparts, as well as three inconclusive candidates requiring confirmation via follow-up observations. These detections significantly increase the number of known radio massive stellar bow shocks and highlight the advantage of dedicated searches with current and next-generation radio telescopes. We investigate the underlying radio emission mechanism for these radio sources, finding a mix of free-free-dominated and synchrotron-dominated systems. We also discuss the non-detected targets by putting constraints on their emission properties and investigating their detectability with future observations. Finally, we propose several future avenues of research to advance the study and understanding of bow shocks at radio frequencies.

Comments: 17 pages, 12 figures, 3 tables. Resubmitted to MNRAS after final round of very minor referee comments


Abstract: 2203.11949
Full Text: [ PostScript, PDF]

Title:Cosmic-Ray Transport in Varying Galactic Environments

Download PDF
Abstract: We study the propagation of mildly-relativistic cosmic rays (CRs) in multiphase interstellar medium environments with conditions typical of nearby disk galaxies. We employ the techniques developed in Armillotta+21 to post-process three high-resolution TIGRESS magnetohydrodynamic simulations modeling local patches of star-forming galactic disks. Together, the three simulations cover a wide range of gas surface density, gravitational potential, and star formation rate (SFR). Our prescription for CR propagation includes the effects of advection by the background gas, streaming along the magnetic field at the local ion Alfvén speed, and diffusion relative to the Alfvén waves, with the diffusion coefficient set by the balance between streaming-driven Alfvén wave excitation and damping mediated by local gas properties. We find that the combined transport processes are more effective in environments with higher SFR. These environments are characterized by higher-velocity hot outflows (created by clustered supernovae) that rapidly advect CRs away from the galactic plane. As a consequence, the ratio of midplane CR pressure to midplane gas pressures decreases with increasing SFR. We also use the post-processed simulations to make predictions regarding potential dynamical impacts of CRs. The relatively flat CR pressure profiles near the midplane argue that they would not provide significant support against gravity for most of the ISM mass. However, the CR pressure gradients are larger than the other pressure gradients in the extra-planar region (|z|>0.5 kpc), suggesting that CRs may affect the dynamics of galactic fountains and/or winds. The degree of this impact is expected to increase in environments with lower SFR.

Comments: Accepted for publication in ApJ


Abstract: 2203.11757
Full Text: [ PostScript, PDF]

Title:Deep Learning-based Imaging in Radio Interferometry

Download PDF
Abstract: The sparse layouts of radio interferometers result in an incomplete sampling of the sky in Fourier space which leads to artifacts in the reconstructed images. Cleaning these systematic effects is essential for the scientific use of radiointerferometric images. Established reconstruction methods are often time-consuming, require expert-knowledge, and suffer from a lack of reproducibility. We have developed a prototype Deep Learning-based method that generates reproducible images in an expedient fashion. To this end, we take advantage of the efficiency of Convolutional Neural Networks to reconstruct image data from incomplete information in Fourier space. The Neural Network architecture is inspired by super-resolution models that utilize residual blocks. Using simulated data of radio galaxies that are composed of Gaussian components we train Deep Learning models whose reconstruction capability is quantified using various measures. The reconstruction performance is evaluated on clean and noisy input data by comparing the resulting predictions with the true source images. We find that source angles and sizes are well reproduced, while the recovered fluxes show substantial scatter, albeit not worse than existing methods without fine-tuning. Finally, we propose more advanced approaches using Deep Learning that include uncertainty estimates and a concept to analyze larger images.

Comments: Accepted for publication in Astronomy & Astrophysics


Abstract: 2203.10710
Full Text: [ PostScript, PDF]

Title:FAIR standards for astronomical data

Download PDF
Abstract: We present an overview of the "FAIR Guiding Principles for scientific data management and stewardship", first published in 2016, and how they relate to astronomical data management. In particular, we discuss the connection between the FAIR principles and IVOA standards, and how data management systems with these standards implemented are close to compliance. We then look at what extra components are required to make astronomical data FAIR. Finally, we give a case study of the All-Sky Virtual Observatory (Australia's node of the VO) and their implementation of the FAIR principles.

Comments: 10 pages, 1 figure, to be published in Proc. ADASS XXXI (2021)


Abstract: 2203.12610
Full Text: [ PostScript, PDF]

Title:AI Poincaré 2.0: Machine Learning Conservation Laws from Differential Equations

Authors:Ziming Liu (MIT), Varun Madhavan (IIT), Max Tegmark (MIT)
Download PDF
Abstract: We present a machine learning algorithm that discovers conservation laws from differential equations, both numerically (parametrized as neural networks) and symbolically, ensuring their functional independence (a non-linear generalization of linear independence). Our independence module can be viewed as a nonlinear generalization of singular value decomposition. Our method can readily handle inductive biases for conservation laws. We validate it with examples including the 3-body problem, the KdV equation and nonlinear Schrödinger equation.

Comments: 17 pages, 10 figures


Abstract: 2203.11945
Full Text: [ PostScript, PDF]

Title:sympy2c: from symbolic expressions to fast C/C++ functions and ODE solvers in Python

Download PDF
Abstract: Computer algebra systems play an important role in science as they facilitate the development of new theoretical models. The resulting symbolic equations are often implemented in a compiled programming language in order to provide fast and portable codes for practical applications. We describe sympy2c, a new Python package designed to bridge the gap between the symbolic development and the numerical implementation of a theoretical model. sympy2c translates symbolic equations implemented in the SymPy Python package to C/C++ code that is optimized using symbolic transformations. The resulting functions can be conveniently used as an extension module in Python. sympy2c is used within the PyCosmo Python package to solve the Einstein-Boltzmann equations, a large system of ODEs describing the evolution of linear perturbations in the Universe. After reviewing the functionalities and usage of sympy2c, we describe its implementation and optimization strategies. This includes, in particular, a novel approach to generate optimized ODE solvers making use of the sparsity of the symbolic Jacobian matrix. We demonstrate its performance using the Einstein-Boltzmann equations as a test case. sympy2c is widely applicable and may prove useful for various areas of computational physics. sympy2c is publicly available at this https URL

Comments: 28 pages, 5 figures, 5 tables, Link to package: this https URL, the described packaged sympy2c is used within arXiv:2112.08395


This page created: Mon Mar 28 13:27:15 ACDT 2022 by seinecke

For a printable title listing click here
For details on generating this page see the instructions. If there are problems with this page (and I expect there will be from time to time) contact Jose.

For previous lists of abstracts of interest click Previous abstracts of interest