Showing 1–50 of 432 results for author: Pohl, M

SpaRG: Sparsely Reconstructed Graphs for Generalizable fMRI Analysis

Authors: Camila González, Yanis Miraoui, Yiran Fan, Ehsan Adeli, Kilian M. Pohl

Abstract: Deep learning can help uncover patterns in resting-state functional Magnetic Resonance Imaging (rs-fMRI) associated with psychiatric disorders and personal traits. Yet the problem of interpreting deep learning findings is rarely more evident than in fMRI analyses, as the data is sensitive to scanning effects and inherently difficult to visualize. We propose a simple approach to mitigate these chal… ▽ More Deep learning can help uncover patterns in resting-state functional Magnetic Resonance Imaging (rs-fMRI) associated with psychiatric disorders and personal traits. Yet the problem of interpreting deep learning findings is rarely more evident than in fMRI analyses, as the data is sensitive to scanning effects and inherently difficult to visualize. We propose a simple approach to mitigate these challenges grounded on sparsification and self-supervision. Instead of extracting post-hoc feature attributions to uncover functional connections that are important to the target task, we identify a small subset of highly informative connections during training and occlude the rest. To this end, we jointly train a (1) sparse input mask, (2) variational autoencoder (VAE), and (3) downstream classifier in an end-to-end fashion. While we need a portion of labeled samples to train the classifier, we optimize the sparse mask and VAE with unlabeled data from additional acquisition sites, retaining only the input features that generalize well. We evaluate our method - Sparsely Reconstructed Graphs (SpaRG) - on the public ABIDE dataset for the task of sex classification, training with labeled cases from 18 sites and adapting the model to two additional out-of-distribution sites with a portion of unlabeled samples. For a relatively coarse parcellation (64 regions), SpaRG utilizes only 1% of the original connections while improving the classification accuracy across domains. Our code can be found at github.com/yanismiraoui/SpaRG. △ Less

Submitted 24 September, 2024; originally announced October 2024.

Future frame prediction in chest cine MR imaging using the PCA respiratory motion model and dynamically trained recurrent neural networks

Authors: Michel Pohl, Mitsuru Uesaka, Hiroyuki Takahashi, Kazuyuki Demachi, Ritu Bhusal Chhatkuli

Abstract: Lung radiotherapy treatment systems are subject to a latency that leads to uncertainty in the estimated tumor location and high irradiation of healthy tissue. This work addresses future frame prediction in chest dynamic MRI sequences to compensate for that delay using RNNs trained with online learning algorithms. The latter enable networks to mitigate irregular movements, as they update synaptic w… ▽ More Lung radiotherapy treatment systems are subject to a latency that leads to uncertainty in the estimated tumor location and high irradiation of healthy tissue. This work addresses future frame prediction in chest dynamic MRI sequences to compensate for that delay using RNNs trained with online learning algorithms. The latter enable networks to mitigate irregular movements, as they update synaptic weights with each new training example. Experiments were conducted using four publicly available 2D thoracic cine-MRI sequences. PCA decomposes the time-varying deformation vector field (DVF), computed with the Lucas-Kanade optical flow algorithm, into static deformation fields and low-dimensional time-dependent weights. We compare various algorithms to forecast the latter: linear regression, least mean squares (LMS), and RNNs trained with real-time recurrent learning (RTRL), unbiased online recurrent optimization, decoupled neural interfaces and sparse 1-step approximation (SnAp-1). That enables estimating the future DVFs and, in turn, the next frames by warping the initial image. Linear regression led to the lowest mean DVF error at a horizon h = 0.32s (the time interval in advance for which the prediction is made), equal to 1.30mm, followed by SnAp-1 and RTRL, whose error increased from 1.37mm to 1.44mm as h increased from 0.62s to 2.20s. Similarly, the structural similarity index measure (SSIM) of LMS decreased from 0.904 to 0.898 as h increased from 0.31s to 1.57s and was the highest among the algorithms compared for the latter horizons. SnAp-1 attained the highest SSIM for h $\geq$ 1.88s, with values of less than 0.898. The predicted images look similar to the original ones, and the highest errors occurred at challenging areas such as the diaphragm boundary at the end-of-inhale phase, where motion variability is more prominent, and regions where out-of-plane motion was more prevalent. △ Less

Submitted 8 October, 2024; originally announced October 2024.

Comments: 28 pages, 16 figures

Spectral Graph Sample Weighting for Interpretable Sub-cohort Analysis in Predictive Models for Neuroimaging

Authors: Magdalini Paschali, Yu Hang Jiang, Spencer Siegel, Camila Gonzalez, Kilian M. Pohl, Akshay Chaudhari, Qingyu Zhao

Abstract: Recent advancements in medicine have confirmed that brain disorders often comprise multiple subtypes of mechanisms, developmental trajectories, or severity levels. Such heterogeneity is often associated with demographic aspects (e.g., sex) or disease-related contributors (e.g., genetics). Thus, the predictive power of machine learning models used for symptom prediction varies across subjects based… ▽ More Recent advancements in medicine have confirmed that brain disorders often comprise multiple subtypes of mechanisms, developmental trajectories, or severity levels. Such heterogeneity is often associated with demographic aspects (e.g., sex) or disease-related contributors (e.g., genetics). Thus, the predictive power of machine learning models used for symptom prediction varies across subjects based on such factors. To model this heterogeneity, one can assign each training sample a factor-dependent weight, which modulates the subject's contribution to the overall objective loss function. To this end, we propose to model the subject weights as a linear combination of the eigenbases of a spectral population graph that captures the similarity of factors across subjects. In doing so, the learned weights smoothly vary across the graph, highlighting sub-cohorts with high and low predictability. Our proposed sample weighting scheme is evaluated on two tasks. First, we predict initiation of heavy alcohol drinking in young adulthood from imaging and neuropsychological measures from the National Consortium on Alcohol and NeuroDevelopment in Adolescence (NCANDA). Next, we detect Dementia vs. Mild Cognitive Impairment (MCI) using imaging and demographic measurements in subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Compared to existing sample weighting schemes, our sample weights improve interpretability and highlight sub-cohorts with distinct characteristics and varying model accuracy. △ Less

Submitted 5 October, 2024; v1 submitted 1 October, 2024; originally announced October 2024.

Brain-Cognition Fingerprinting via Graph-GCCA with Contrastive Learning

Authors: Yixin Wang, Wei Peng, Yu Zhang, Ehsan Adeli, Qingyu Zhao, Kilian M. Pohl

Abstract: Many longitudinal neuroimaging studies aim to improve the understanding of brain aging and diseases by studying the dynamic interactions between brain function and cognition. Doing so requires accurate encoding of their multidimensional relationship while accounting for individual variability over time. For this purpose, we propose an unsupervised learning model (called \underline{\textbf{Co}}ntra… ▽ More Many longitudinal neuroimaging studies aim to improve the understanding of brain aging and diseases by studying the dynamic interactions between brain function and cognition. Doing so requires accurate encoding of their multidimensional relationship while accounting for individual variability over time. For this purpose, we propose an unsupervised learning model (called \underline{\textbf{Co}}ntrastive Learning-based \underline{\textbf{Gra}}ph Generalized \underline{\textbf{Ca}}nonical Correlation Analysis (CoGraCa)) that encodes their relationship via Graph Attention Networks and generalized Canonical Correlational Analysis. To create brain-cognition fingerprints reflecting unique neural and cognitive phenotype of each person, the model also relies on individualized and multimodal contrastive learning. We apply CoGraCa to longitudinal dataset of healthy individuals consisting of resting-state functional MRI and cognitive measures acquired at multiple visits for each participant. The generated fingerprints effectively capture significant individual differences and outperform current single-modal and CCA-based multimodal models in identifying sex and age. More importantly, our encoding provides interpretable interactions between those two modalities. △ Less

Submitted 20 September, 2024; originally announced September 2024.

Evaluating the Quality of Brain MRI Generators

Authors: Jiaqi Wu, Wei Peng, Binxu Li, Yu Zhang, Kilian M. Pohl

Abstract: Deep learning models generating structural brain MRIs have the potential to significantly accelerate discovery of neuroscience studies. However, their use has been limited in part by the way their quality is evaluated. Most evaluations of generative models focus on metrics originally designed for natural images (such as structural similarity index and Frechet inception distance). As we show in a c… ▽ More Deep learning models generating structural brain MRIs have the potential to significantly accelerate discovery of neuroscience studies. However, their use has been limited in part by the way their quality is evaluated. Most evaluations of generative models focus on metrics originally designed for natural images (such as structural similarity index and Frechet inception distance). As we show in a comparison of 6 state-of-the-art generative models trained and tested on over 3000 MRIs, these metrics are sensitive to the experimental setup and inadequately assess how well brain MRIs capture macrostructural properties of brain regions (i.e., anatomical plausibility). This shortcoming of the metrics results in inconclusive findings even when qualitative differences between the outputs of models are evident. We therefore propose a framework for evaluating models generating brain MRIs, which requires uniform processing of the real MRIs, standardizing the implementation of the models, and automatically segmenting the MRIs generated by the models. The segmentations are used for quantifying the plausibility of anatomy displayed in the MRIs. To ensure meaningful quantification, it is crucial that the segmentations are highly reliable. Our framework rigorously checks this reliability, a step often overlooked by prior work. Only 3 of the 6 generative models produced MRIs, of which at least 95% had highly reliable segmentations. More importantly, the assessment of each model by our framework is in line with qualitative assessments, reinforcing the validity of our approach. △ Less

Submitted 12 September, 2024; originally announced September 2024.

Comments: Accepted MICCAI 20224

Latent 3D Brain MRI Counterfactual

Authors: Wei Peng, Tian Xia, Fabio De Sousa Ribeiro, Tomas Bosschieter, Ehsan Adeli, Qingyu Zhao, Ben Glocker, Kilian M. Pohl

Abstract: The number of samples in structural brain MRI studies is often too small to properly train deep learning models. Generative models show promise in addressing this issue by effectively learning the data distribution and generating high-fidelity MRI. However, they struggle to produce diverse, high-quality data outside the distribution defined by the training data. One way to address the issue is usi… ▽ More The number of samples in structural brain MRI studies is often too small to properly train deep learning models. Generative models show promise in addressing this issue by effectively learning the data distribution and generating high-fidelity MRI. However, they struggle to produce diverse, high-quality data outside the distribution defined by the training data. One way to address the issue is using causal models developed for 3D volume counterfactuals. However, accurately modeling causality in high-dimensional spaces is a challenge so that these models generally generate 3D brain MRIS of lower quality. To address these challenges, we propose a two-stage method that constructs a Structural Causal Model (SCM) within the latent space. In the first stage, we employ a VQ-VAE to learn a compact embedding of the MRI volume. Subsequently, we integrate our causal model into this latent space and execute a three-step counterfactual procedure using a closed-form Generalized Linear Model (GLM). Our experiments conducted on real-world high-resolution MRI data (1mm) demonstrate that our method can generate high-quality 3D MRI counterfactuals. △ Less

Submitted 9 September, 2024; originally announced September 2024.

Particle acceleration, escape and non-thermal emission from core-collapse supernovae inside non-identical wind-blown bubbles

Authors: Samata Das, Robert Brose, Martin Pohl, Dominique M. -A. Meyer, Iurii Sushch

Abstract: In the core-collapse scenario, the supernova remnants evolve inside the complex wind-blown bubbles, structured by massive progenitors during their lifetime. Therefore, particle acceleration and the emissions from these SNRs can carry the fingerprints of the evolutionary sequences of the progenitor stars. We time-dependently investigate the impact of the ambient environment of core-collapse SNRs… ▽ More In the core-collapse scenario, the supernova remnants evolve inside the complex wind-blown bubbles, structured by massive progenitors during their lifetime. Therefore, particle acceleration and the emissions from these SNRs can carry the fingerprints of the evolutionary sequences of the progenitor stars. We time-dependently investigate the impact of the ambient environment of core-collapse SNRs on particle spectra and the emissions. We use the RATPaC code to model the particle acceleration at the SNRs with progenitors having ZAMS masses of 20 Msol and 60 Msol. We have constructed the pre-supernova circumstellar medium by solving the hydrodynamic equations for the lifetime of the progenitor stars. Then, the transport equation for cosmic rays, and magnetic turbulence in test-particle approximation along with the induction equation for the evolution of large-scale magnetic field have been solved simultaneously with the hydrodynamic equations for the expansion of SNRs inside the pre-supernova CSM. The structure of the wind bubbles along with the magnetic field and the scattering turbulence regulate the spectra of accelerated particles for both SNRs. For the 60 Msol progenitor the spectral index reaches 2.4 even below 10 GeV during the propagation of the SNR shock inside the hot shocked wind. In contrast, we have not observed persistent soft spectra at earlier evolutionary stages of the SNR with 20 Msol progenitor, for which the spectral index becomes 2.2 only for a brief period. Later, the spectra become soft above ~10 GeV for both SNRs, as weak driving of turbulence permits the escape of high-energy particles from the remnants. The emission morphology of the SNRs strongly depends on the type of progenitors. For instance, the radio morphology of the SNR with 20 Msol progenitor is centre-filled at early stages whereas that for the more massive progenitor is shell-like. △ Less

Submitted 28 August, 2024; originally announced August 2024.

Journal ref: A&A, 689 (2024) A9

An indirect search for dark matter with a combined analysis of dwarf spheroidal galaxies from VERITAS

Authors: A. Acharyya, C. B. Adams, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, A. Falcone, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Holder, T. B. Humensky, W. Jin, M. N. Johnson, P. Kaaret, M. Kertzman , et al. (37 additional authors not shown)

Abstract: Understanding the nature and identity of dark matter is a key goal in the physics community. In the case that TeV-scale dark matter particles decay or annihilate into standard model particles, very-high-energy (VHE) gamma rays (greater than 100 GeV) will be present in the final state. The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an imaging atmospheric Cherenkov telescop… ▽ More Understanding the nature and identity of dark matter is a key goal in the physics community. In the case that TeV-scale dark matter particles decay or annihilate into standard model particles, very-high-energy (VHE) gamma rays (greater than 100 GeV) will be present in the final state. The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an imaging atmospheric Cherenkov telescope array that can indirectly detect VHE gamma rays in an energy range of 100 GeV to > 30 TeV. Dwarf spheroidal galaxies (dSphs) are ideal candidates in the search for dark matter due to their high dark matter content, high mass-to-light ratios, and their low gamma-ray fluxes from astrophysical processes. This study uses a legacy data set of 638 hours collected on 17 dSphs, built over 11 years with an observing strategy optimized according to the dark matter content of the targets. The study addresses a broad dark matter particle mass range, extending from 200 GeV to 30 PeV. In the absence of a detection, we set the upper limits on the dark matter velocity-weighted annihilation cross section. △ Less

Submitted 6 August, 2024; v1 submitted 23 July, 2024; originally announced July 2024.

Comments: 18 pages, 7 figures, 3 tables, accepted in PRD

A multi-wavelength study to decipher the 2017 flare of the blazar OJ 287

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. Brill, J. P. Caldwell, M. Carini, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, Q. Feng, J. P. Finley, J. Foote, L. Fortson, A. Furniss, G. Gallagher, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Hoang , et al. (49 additional authors not shown)

Abstract: In February 2017, the blazar OJ~287 underwent a period of intense multiwavelength activity. It reached a new historic peak in the soft X-ray (0.3-10 keV) band, as measured by Swift-XRT. This event coincides with a very-high-energy (VHE) $γ$-ray outburst that led VERITAS to detect emission above 100 GeV, with a detection significance of $10σ$ (from 2016 December 9 to 2017 March 31). The time-averag… ▽ More In February 2017, the blazar OJ~287 underwent a period of intense multiwavelength activity. It reached a new historic peak in the soft X-ray (0.3-10 keV) band, as measured by Swift-XRT. This event coincides with a very-high-energy (VHE) $γ$-ray outburst that led VERITAS to detect emission above 100 GeV, with a detection significance of $10σ$ (from 2016 December 9 to 2017 March 31). The time-averaged VHE $γ$-ray spectrum was consistent with a soft power law ($Γ= -3.81 \pm 0.26$) and an integral flux corresponding to $\sim2.4\%$ that of the Crab Nebula above the same energy. Contemporaneous data from multiple instruments across the electromagnetic spectrum reveal complex flaring behavior, primarily in the soft X-ray and VHE bands. To investigate the possible origin of such an event, our study focuses on three distinct activity states: before, during, and after the February 2017 peak. The spectral energy distributions during these periods suggest the presence of at least two non-thermal emission zones, with the more compact one responsible for the observed flare. Broadband modeling results and observations of a new radio knot in the jet of OJ~287 in 2017 are consistent with a flare originating from a strong recollimation shock outside the radio core. △ Less

Submitted 26 August, 2024; v1 submitted 16 July, 2024; originally announced July 2024.

Comments: 21 pages, 7 figures, accepted for publication in The Astrophysical Journal

Learning Human Detected Differences in Directed Acyclic Graphs

Authors: Kathrin Guckes, Alena Beyer, Margit Pohl, Tatiana von Landesberger

Abstract: Prior research has shown that human perception of similarity differs from mathematical measures in visual comparison tasks, including those involving directed acyclic graphs. This divergence can lead to missed differences and skepticism about algorithmic results. To address this, we aim to learn the structural differences humans detect in graphs visually. We want to visualize these human-detected… ▽ More Prior research has shown that human perception of similarity differs from mathematical measures in visual comparison tasks, including those involving directed acyclic graphs. This divergence can lead to missed differences and skepticism about algorithmic results. To address this, we aim to learn the structural differences humans detect in graphs visually. We want to visualize these human-detected differences alongside actual changes, enhancing credibility and aiding users in spotting overlooked differences. Our approach aligns with recent research in machine learning capturing human behavior. We provide a data augmentation algorithm, a dataset, and a machine learning model to support this task. This work fills a gap in learning differences in directed acyclic graphs and contributes to better comparative visualizations. △ Less

Submitted 8 June, 2024; originally announced June 2024.

A Shape Change Enhancing Hierarchical Layout for the Pairwise Comparison of Directed Acyclic Graphs

Authors: Kathrin Guckes, Marc Schäpers, Margit Pohl, Andreas Kerren, Tatiana von Landesberger

Abstract: Comparing directed acyclic graphs is essential in various fields such as healthcare, social media, finance, biology, and marketing. DAGs often result from contagion processes over networks, including information spreading, retweet activity, disease transmission, financial crisis propagation, malware spread, and gene mutations. For instance, in disease spreading, an infected patient can transmit th… ▽ More Comparing directed acyclic graphs is essential in various fields such as healthcare, social media, finance, biology, and marketing. DAGs often result from contagion processes over networks, including information spreading, retweet activity, disease transmission, financial crisis propagation, malware spread, and gene mutations. For instance, in disease spreading, an infected patient can transmit the disease to contacts, making it crucial to analyze and predict scenarios. Similarly, in finance, understanding the effects of saving or not saving specific banks during a crisis is vital. Experts often need to identify small differences between DAGs, such as changes in a few nodes or edges. Even the presence or absence of a single edge can be significant. Visualization plays a crucial role in facilitating these comparisons. However, standard hierarchical layout algorithms struggle to visualize subtle changes effectively. The typical hierarchical layout, with the root on top, is preferred due to its performance in comparison to other layouts. Nevertheless, these standard algorithms prioritize single-graph aesthetics over comparison suitability, making it challenging for users to spot changes. To address this issue, we propose a layout that enhances shape changes in DAGs while minimizing the impact on aesthetics. Our approach involves outwardly swapping changes, altering the DAG's shape. We introduce new drawing criteria. Our layout builds upon a Sugiyama-like hierarchical layout and implements these criteria through two extensions. We designed it this way to maintain interchangeability and accommodate future optimizations, such as pseudo-nodes for edge crossing minimization. In our evaluations, our layout achieves excellent results, with edge crossing aesthetics averaging around 0.8 (on a scale of 0 to 1). Additionally, our layout outperforms the base implementation by an average of 60-75\%. △ Less

Submitted 8 June, 2024; originally announced June 2024.

GuidelineExplorer -- Navigating through the Forrest of Actionable Guidelines on Node-Link Graph Visualization

Authors: Kathrin Guckes, Lisa Eisenhardt, Margit Pohl, Tatiana von Landesberger

Abstract: Creating graph visualizations involves many decisions, such as layout, node and edge appearance, and color choices. These decisions are challenging due to the multitude of options available. For instance, graph layout can be force-directed or orthogonal, and edges can be curved, tapered, partially drawn, or animated. Thus, research offers a multitude of guidelines to optimize graph visualizations… ▽ More Creating graph visualizations involves many decisions, such as layout, node and edge appearance, and color choices. These decisions are challenging due to the multitude of options available. For instance, graph layout can be force-directed or orthogonal, and edges can be curved, tapered, partially drawn, or animated. Thus, research offers a multitude of guidelines to optimize graph visualizations for human perception and usability. Guidelines can be actionable, providing direct instructions, or non-actionable, specifying what to avoid. This work focuses on actionable guidelines for node-link diagrams, aiding designers in making better decisions. Given the abundance of graph visualization research and the difficulty in navigating it, this work aims to collect and structure actionable guidelines for node-linkvisualizations. To demonstrate the general applicability of our approach to structuring actionable guidelines for node-link diagrams, we also included guidelines for visualizing graphs as matrices. It also proposes a visual interactive system, GuidelineExplorer, to apply guidelines directly to graphs, streamlining the design process and promoting collaboration within the research community. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Poynting flux transport channels formed in polar cap regions of neutron star magnetospheres

Authors: Jan Benáček, Andrey Timokhin, Patricio A. Muñoz, Axel Jessner, Tatiana Rievajová, Martin Pohl, Jörg Büchner

Abstract: Pair cascades in polar cap regions of neutron stars are considered to be an essential process in various models of coherent radio emissions of pulsars. The cascades produce pair plasma bunch discharges in quasi-periodic spark events. The cascade properties, and therefore also the coherent radiation, depend strongly on the magnetospheric plasma properties and vary significantly across and along the… ▽ More Pair cascades in polar cap regions of neutron stars are considered to be an essential process in various models of coherent radio emissions of pulsars. The cascades produce pair plasma bunch discharges in quasi-periodic spark events. The cascade properties, and therefore also the coherent radiation, depend strongly on the magnetospheric plasma properties and vary significantly across and along the polar cap. It is furthermore still uncertain from where the radio emission emanates in polar cap region. We investigate the generation of electromagnetic waves by pair cascades and their propagation in the polar cap for three representative inclination angles of a magnetic dipole, $0^\circ$, $45^\circ$, and $90^\circ$. 2D particle-in-cell simulations that include quantum-electrodynamic pair cascades are used in a charge limited flow from the star surface. We found that the discharge properties are strongly dependent on the magnetospheric current profile in the polar cap and that transport channels for high intensity Poynting flux are formed along magnetic field lines where the magnetospheric currents approach zero and where the plasma cannot carry the magnetospheric currents. There, the parallel Poynting flux component is efficiently transported away from the star and may eventually escape the magnetosphere as coherent radio waves. The Poynting flux decreases faster with the distance from the star in regions of high magnetospheric currents. Our model shows that no process of energy conversion from particles to waves is necessary for the coherent radio wave emission. Moreover, the pulsar radio beam does not have a cone structure, but rather the radiation generated by the oscillating electric gap fields directly escapes along open magnetic field lines in which no pair creation occurs. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: Submitted to Astronomy & Astrophysics

Supernova remnants of red supergiants: from barrels to Cygnus loops

Authors: D M A Meyer, P F Velazquez, M Pohl, K Egberts, M Petrov, M A Villagran, D F Torres, R Batzofin

Abstract: Core-collapse supernova remnants are the nebular leftover of defunct massive stars which have died during a supernova explosion, mostly while undergoing the red supergiant phase of their evolution. The morphology and emission properties of those remnants are a function of the distribution of circumstellar material at the moment of the supernova, the intrisic properties of the explosion, as well as… ▽ More Core-collapse supernova remnants are the nebular leftover of defunct massive stars which have died during a supernova explosion, mostly while undergoing the red supergiant phase of their evolution. The morphology and emission properties of those remnants are a function of the distribution of circumstellar material at the moment of the supernova, the intrisic properties of the explosion, as well as those of the ambient medium. By means of 2.5 dimensional numerical magnetohydrodynamics simulations, we model the long term evolution of supernova remnants generated by runaway rotating massive stars moving into a magnetised interstellar medium. Radiative transfer calculations reveal that the projected non-thermal emission of the supernova remnants decreases with time, i.e. older remnants are fainter than younger ones. Older (80 kyr) supernova remnants whose progenitors were moving with space velocity corresponding to a Mach number M = 1 (v_star = 20 km/s ) in the Galactic plane of the ISM (nISM = 1/cm3 ) are brighter in synchrotron than when moving with a Mach number M = 2 (v_star = 40 km/s ). We show that runaway red supergiant progenitors first induce an asymmetric non thermal 1.4 GHz barrel like synchrotron supernova remnants (at the age of about 8 kyr), before further evolving to adopt a Cygnus loop like shape (at about 80 kyr). It is conjectured that a significative fraction of supernova remnants are currently in this bilateral-to-Cygnus-loop evolutionary sequence, and that this should be taken into account in the data interpretation of the forthcoming Cherenkov Telescope Array (CTA) observatory. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: Accepted at Astronomie and Astrophysics

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 23 July, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 44 pages JCAP style (excluding author list and references), 19 figures; minor changes to match published version

Journal ref: JCAP 07 (2024) 047

Respiratory motion forecasting with online learning of recurrent neural networks for safety enhancement in externally guided radiotherapy

Authors: Michel Pohl, Mitsuru Uesaka, Hiroyuki Takahashi, Kazuyuki Demachi, Ritu Bhusal Chhatkuli

Abstract: In lung radiotherapy, infrared cameras can record the location of reflective objects on the chest to infer the position of the tumor moving due to breathing, but treatment system latencies hinder radiation beam precision. Real-time recurrent learning (RTRL), is a potential solution as it can learn patterns within non-stationary respiratory data but has high complexity. This study assesses the capa… ▽ More In lung radiotherapy, infrared cameras can record the location of reflective objects on the chest to infer the position of the tumor moving due to breathing, but treatment system latencies hinder radiation beam precision. Real-time recurrent learning (RTRL), is a potential solution as it can learn patterns within non-stationary respiratory data but has high complexity. This study assesses the capabilities of resource-efficient online RNN algorithms, namely unbiased online recurrent optimization (UORO), sparse-1 step approximation (SnAp-1), and decoupled neural interfaces (DNI) to forecast respiratory motion during radiotherapy treatment accurately. We use time series containing the 3D position of external markers on the chest of healthy subjects. We propose efficient implementations for SnAp-1 and DNI based on compression of the influence and immediate Jacobian matrices and an accurate update of the linear coefficients used in credit assignment estimation, respectively. The original sampling frequency was 10Hz; we performed resampling at 3.33Hz and 30Hz. We use UORO, SnAp-1, and DNI to forecast each marker's 3D position with horizons (the time interval in advance for which the prediction is made) h<=2.1s and compare them with RTRL, least mean squares, and linear regression. RNNs trained online achieved similar or better accuracy than most previous works using larger training databases and deep learning, even though we used only the first minute of each sequence to predict motion within that exact sequence. SnAp-1 had the lowest normalized root mean square errors (nRMSE) averaged over the horizon values considered, equal to 0.335 and 0.157, at 3.33Hz and 10.0Hz, respectively. Similarly, UORO had the highest accuracy at 30Hz, with an nRMSE of 0.0897. DNI's inference time, equal to 6.8ms per time step at 30Hz (Intel Core i7-13700 CPU), was the lowest among the RNN methods examined. △ Less

Submitted 3 March, 2024; originally announced March 2024.

Comments: 34 pages, 11 figures

Robust inference of the Galactic centre gamma-ray excess spatial properties

Authors: Deheng Song, Christopher Eckner, Chris Gordon, Francesca Calore, Oscar Macias, Kevork N. Abazajian, Shunsaku Horiuchi, Manoj Kaplinghat, Martin Pohl

Abstract: The gamma-ray Fermi-LAT Galactic centre excess (GCE) has puzzled scientists for over 15 years. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and poi… ▽ More The gamma-ray Fermi-LAT Galactic centre excess (GCE) has puzzled scientists for over 15 years. Despite ongoing debates about its properties, and especially its spatial distribution, its nature remains elusive. We scrutinize how the estimated spatial morphology of this excess depends on models for the Galactic diffuse emission, focusing particularly on the extent to which the Galactic plane and point sources are masked. Our main aim is to compare a spherically symmetric morphology - potentially arising from the annihilation of dark matter (DM) particles - with a boxy morphology - expected if faint unresolved sources in the Galactic bulge dominate the excess emission. Recent claims favouring a DM-motivated template for the GCE are shown to rely on a specific Galactic bulge template, which performs worse than other templates for the Galactic bulge. We find that a non-parametric model of the Galactic bulge derived from the VVV survey results in a significantly better fit for the GCE than DM-motivated templates. This result is independent of whether a GALPROP-based model or a more non-parametric ring-based model is used to describe the diffuse Galactic emission. This conclusion remains true even when additional freedom is added in the background models, allowing for non-parametric modulation of the model components and substantially improving the fit quality. When adopted, optimized background models provide robust results in terms of preference for a boxy bulge morphology of the GCE, regardless of the mask applied to the Galactic plane. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 21 pages, 11 figures

Report number: LAPTH-009/24

Journal ref: 2024MNRAS.530.4395S

Nonlinear feedback of the electrostatic instability on the blazar-induced pair beam and GeV cascade

Authors: Mahmoud Alawashra, Martin Pohl

Abstract: Relativistic pair beams produced in the cosmic voids by TeV gamma rays from blazars are expected to produce a detectable GeV-scale cascade that is missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields or, alternatively, an energy loss of the beam due to the beam-plasma instability. Here, we study how t… ▽ More Relativistic pair beams produced in the cosmic voids by TeV gamma rays from blazars are expected to produce a detectable GeV-scale cascade that is missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields or, alternatively, an energy loss of the beam due to the beam-plasma instability. Here, we study how the beam-plasma instability feeds back on the beam, using a realistic two-dimensional beam distribution. We find that the instability broadens the beam opening angles significantly without any significant energy loss, thus confirming a recent feedback study on a simplified one-dimensional beam distribution. However, narrowing diffusion feedback of the beam particles with Lorentz factors less than $10^6$ might become relevant even though initially it is negligible. Finally, when considering the continuous creation of TeV pairs, we find that the beam distribution and the wave spectrum reach a new quasi-steady state, in which the scattering of beam particles persists and the beam opening angle may increase by a factor of hundreds. Understanding the implications on the GeV cascade emission requires accounting for inverse-Compton cooling. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: Accepted for publication in the Astrophysical Journal. 16 pages and 12 figures

An Angular Diameter Measurement of $β$ UMa via Stellar Intensity Interferometry with the VERITAS Observatory

Authors: A. Acharyya, J. P. Aufdenberg, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. J. Chromey, J. D. Davis, Q. Feng, G. M. Foote, A. Furniss, W. Hanlon, C. E. Hinrichs, J. Holder, W. Jin, P. Kaaret, M. Kertzman, D. Kieda, T. K. Kleiner, N. Korzoun, T. LeBohec, M. A. Lisa, M. Lundy, N. Matthews, C. E McGrath , et al. (22 additional authors not shown)

Abstract: We use the VERITAS imaging air Cherenkov Telescope (IACT) array to obtain the first measured angular diameter of $β$ UMa at visual wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak ($β$ UMa), an A1-type subgiant… ▽ More We use the VERITAS imaging air Cherenkov Telescope (IACT) array to obtain the first measured angular diameter of $β$ UMa at visual wavelengths using stellar intensity interferometry (SII) and independently constrain the limb-darkened angular diameter. The age of the Ursa Major moving group has been assessed from the ages of its members, including nuclear member Merak ($β$ UMa), an A1-type subgiant, by comparing effective temperature and luminosity constraints to model stellar evolution tracks. Previous interferometric limb-darkened angular-diameter measurements of $β$ UMa in the near-infrared (CHARA Array, $1.149 \pm 0.014$ mas) and mid-infrared (Keck Nuller, $1.08 \pm 0.07$ mas), together with the measured parallax and bolometric flux, have constrained the effective temperature. This paper presents current VERITAS-SII observation and analysis procedures to derive squared visibilities from correlation functions. We fit the resulting squared visibilities to find a limb-darkened angular diameter of $1.07 \pm 0.04 {\rm (stat)} \pm 0.05$ (sys) mas, using synthetic visibilities from a stellar atmosphere model that provides a good match to the spectrum of $β$ UMa in the optical wave band. The VERITAS-SII limb-darkened angular diameter yields an effective temperature of $9700\pm200\pm 200$ K, consistent with ultraviolet spectrophotometry, and an age of $390\pm 29 \pm 32 $ Myr, using MESA Isochrones and Stellar Tracks (MIST). This age is consistent with $408 \pm 6$ Myr from the CHARA Array angular diameter. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Information acquisition in Adapt/Exchange decisions: When do people check alternative solution principles?

Authors: Romy Müller, Maria Pohl

Abstract: Many problems can be solved in two ways: either by adapting an existing solution, or by exchanging it for a new one. To investigate under what conditions people consider new solutions, we traced their information acquisition processes in a simulated mechanical engineering task. Within a multi-step optimisation procedure, participants could either adapt the properties of a currently used machine co… ▽ More Many problems can be solved in two ways: either by adapting an existing solution, or by exchanging it for a new one. To investigate under what conditions people consider new solutions, we traced their information acquisition processes in a simulated mechanical engineering task. Within a multi-step optimisation procedure, participants could either adapt the properties of a currently used machine component, or exchange this component for a new one. They had the opportunity to check whether the solutions met a set of requirements, which was varied systematically. We investigated whether participants would consistently check both solutions, or whether they would satisfice, ignoring the new solution as long as the current one was good enough. The results clearly refuted consistent checking, but only partly confirmed satisficing. On the one hand, participants indeed checked the new solution least often when the current one was applicable without problems. On the other hand, in this case the new solution still was not fully ignored. However, the latter finding could be traced back to a few participants who diverged from our anticipated strategy of first checking the current solution, and directly went for the new one. The results suggest that in Adapt/Exchange decisions, people do not usually check both solutions in an unbiased manner, but rely on existing solutions as long as they are good enough. △ Less

Submitted 30 December, 2023; originally announced January 2024.

Few Shot Part Segmentation Reveals Compositional Logic for Industrial Anomaly Detection

Authors: Soopil Kim, Sion An, Philip Chikontwe, Myeongkyun Kang, Ehsan Adeli, Kilian M. Pohl, Sang Hyun Park

Abstract: Logical anomalies (LA) refer to data violating underlying logical constraints e.g., the quantity, arrangement, or composition of components within an image. Detecting accurately such anomalies requires models to reason about various component types through segmentation. However, curation of pixel-level annotations for semantic segmentation is both time-consuming and expensive. Although there are s… ▽ More Logical anomalies (LA) refer to data violating underlying logical constraints e.g., the quantity, arrangement, or composition of components within an image. Detecting accurately such anomalies requires models to reason about various component types through segmentation. However, curation of pixel-level annotations for semantic segmentation is both time-consuming and expensive. Although there are some prior few-shot or unsupervised co-part segmentation algorithms, they often fail on images with industrial object. These images have components with similar textures and shapes, and a precise differentiation proves challenging. In this study, we introduce a novel component segmentation model for LA detection that leverages a few labeled samples and unlabeled images sharing logical constraints. To ensure consistent segmentation across unlabeled images, we employ a histogram matching loss in conjunction with an entropy loss. As segmentation predictions play a crucial role, we propose to enhance both local and global sample validity detection by capturing key aspects from visual semantics via three memory banks: class histograms, component composition embeddings and patch-level representations. For effective LA detection, we propose an adaptive scaling strategy to standardize anomaly scores from different memory banks in inference. Extensive experiments on the public benchmark MVTec LOCO AD reveal our method achieves 98.1% AUROC in LA detection vs. 89.6% from competing methods. △ Less

Submitted 15 April, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: Accepted in AAAI2024

AM$^3$: An Open-Source Tool for Time-Dependent Lepto-Hadronic Modeling of Astrophysical Sources

Authors: Marc Klinger, Annika Rudolph, Xavier Rodrigues, Chengchao Yuan, Gaëtan Fichet de Clairfontaine, Anatoli Fedynitch, Walter Winter, Martin Pohl, Shan Gao

Abstract: We present the AM$^3$ (``Astrophysical Multi-Messenger Modeling'') software, which has been successfully used in the past to simulate the multi-messenger emission, including neutrinos, from active galactic nuclei, including the blazar sub-class, gamma-ray bursts, and tidal disruption events. AM$^3$ is a documented state-of-the-art open source software that efficiently solves the coupled integro-di… ▽ More We present the AM$^3$ (``Astrophysical Multi-Messenger Modeling'') software, which has been successfully used in the past to simulate the multi-messenger emission, including neutrinos, from active galactic nuclei, including the blazar sub-class, gamma-ray bursts, and tidal disruption events. AM$^3$ is a documented state-of-the-art open source software that efficiently solves the coupled integro-differential equations for the spectral and temporal evolution of the relevant particle densities (photons, electrons, positrons, protons, neutrons, pions, muons, and neutrinos). AM$^3$ includes all relevant non-thermal processes (synchrotron, inverse Compton scattering, photon-photon annihilation, proton-proton and proton-photon pion production, and photo-pair production). The software self-consistently calculates the full cascade of primary and secondary particles, outperforming simple test-particle approaches, and allows for non-linear feedback and predictions in the time domain. It also allows to track separately the contributions of different radiative processes to the overall photon and neutrino spectra, including the different hadronic interaction channels. With its efficient hybrid solver combining analytical and numerical techniques, AM$^3$ combines efficiency and accuracy at a user-adjustable level. We describe the technical details of the numerical framework and present examples of applications to various astrophysical environments. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: The BSD-licensed source code can be found under https://gitlab.desy.de/am3/am3, for the user guide and documentation see https://am3.readthedocs.io/en/latest/

VERITAS contributions to the 38th International Cosmic Ray Conference

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. L. Christiansen, A. J. Chromey, A. Duerr, M. Errando, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, W. Hanlon, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, Z. Hughes, T. B. Humensky, W. Jin, M. N. Johnson, M. Kertzman , et al. (39 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. Compilation of papers presented by the VERITAS Collaboration at the 38th International Cosmic Ray Conference (ICRC), held July 26 through August 3, 2023 in Nagoya, Japan. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: html page. ICRC 2023, Nagoya, Japan

Kinetic simulations of non-relativistic high-Mach-number perpendicular shocks propagating in a turbulent medium

Authors: Karol Fulat, Artem Bohdan, Gabriel Torralba Paz, Martin Pohl

Abstract: Strong non-relativistic shocks are known to accelerate particles up to relativistic energies. However, for Diffusive Shock Acceleration electrons must have a highly suprathermal energy, implying a need for very efficient pre-acceleration. Most published studies consider shocks propagating through homogeneous plasma, which is an unrealistic assumption for astrophysical environments. Using 2D3V part… ▽ More Strong non-relativistic shocks are known to accelerate particles up to relativistic energies. However, for Diffusive Shock Acceleration electrons must have a highly suprathermal energy, implying a need for very efficient pre-acceleration. Most published studies consider shocks propagating through homogeneous plasma, which is an unrealistic assumption for astrophysical environments. Using 2D3V particle-in-cell simulations, we investigate electron acceleration and heating processes at non-relativistic high-Mach-number shocks in electron-ion plasma with a turbulent upstream medium. For this purpose slabs of plasma with compressive turbulence are separately simulated and then inserted into shock simulations, which requires matching of the plasma slabs at the interface. Using a novel procedure of matching electromagnetic fields and currents, we perform simulations of perpendicular shocks setting different intensities of density fluctuations ($\lesssim 10\%$) in the upstream. The new simulation technique provides a framework for studying shocks propagating in turbulent media. We explore the impact of the fluctuations on electron heating, the dynamics of upstream electrons, and the driving of plasma instabilities. Our results indicate that while the presence of the turbulence enhances variations in the upstream magnetic field, their levels remain too low to influence significantly the behavior of electrons at perpendicular shocks. △ Less

Submitted 22 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ, 17 pages, 10 figures

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

Metadata-Conditioned Generative Models to Synthesize Anatomically-Plausible 3D Brain MRIs

Authors: Wei Peng, Tomas Bosschieter, Jiahong Ouyang, Robert Paul, Ehsan Adeli, Qingyu Zhao, Kilian M. Pohl

Abstract: Generative AI models hold great potential in creating synthetic brain MRIs that advance neuroimaging studies by, for example, enriching data diversity. However, the mainstay of AI research only focuses on optimizing the visual quality (such as signal-to-noise ratio) of the synthetic MRIs while lacking insights into their relevance to neuroscience. To gain these insights with respect to T1-weighted… ▽ More Generative AI models hold great potential in creating synthetic brain MRIs that advance neuroimaging studies by, for example, enriching data diversity. However, the mainstay of AI research only focuses on optimizing the visual quality (such as signal-to-noise ratio) of the synthetic MRIs while lacking insights into their relevance to neuroscience. To gain these insights with respect to T1-weighted MRIs, we first propose a new generative model, BrainSynth, to synthesize metadata-conditioned (e.g., age- and sex-specific) MRIs that achieve state-of-the-art visual quality. We then extend our evaluation with a novel procedure to quantify anatomical plausibility, i.e., how well the synthetic MRIs capture macrostructural properties of brain regions, and how accurately they encode the effects of age and sex. Results indicate that more than half of the brain regions in our synthetic MRIs are anatomically accurate, i.e., with a small effect size between real and synthetic MRIs. Moreover, the anatomical plausibility varies across cortical regions according to their geometric complexity. As is, our synthetic MRIs can significantly improve the training of a Convolutional Neural Network to identify accelerated aging effects in an independent study. These results highlight the opportunities of using generative AI to aid neuroimaging research and point to areas for further improvement. △ Less

Submitted 6 October, 2023; originally announced October 2023.

A multi-wavelength investigation of PSR J2229+6114 and its pulsar wind nebula in the radio, X-ray, and gamma-ray bands

Authors: I. Pope, K. Mori, M. Abdelmaguid, J. D. Gelfand, S. P. Reynolds, S. Safi-Harb, C. J. Hailey, H. An, VERITAS Collaboration, :, P. Bangale, P. Batista, W. Benbow, J. H. Buckley, M. Capasso, J. L. Christiansen, A. J. Chromey, A. Falcone, Q. Feng, J. P. Finley, G. M Foote, G. Gallagher, W. F Hanlon, D. Hanna, O. Hervet , et al. (35 additional authors not shown)

Abstract: G106.3$+$2.7, commonly considered a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct ("head" & "tail") regions in the radio band. A discovery of very-high-energy (VHE) gamma-ray emission ($E_γ> 100$ GeV) followed by the recent detection of ultra-high-energy (UHE) gamma-ray emission ($E_γ> 100$ TeV) from the tail region suggests tha… ▽ More G106.3$+$2.7, commonly considered a composite supernova remnant (SNR), is characterized by a boomerang-shaped pulsar wind nebula (PWN) and two distinct ("head" & "tail") regions in the radio band. A discovery of very-high-energy (VHE) gamma-ray emission ($E_γ> 100$ GeV) followed by the recent detection of ultra-high-energy (UHE) gamma-ray emission ($E_γ> 100$ TeV) from the tail region suggests that G106.3$+$2.7 is a PeVatron candidate. We present a comprehensive multi-wavelength study of the Boomerang PWN (100" around PSR J2229+6114) using archival radio and Chandra data obtained from two decades ago, a new NuSTAR X-ray observation from 2020, and upper limits on gamma-ray fluxes obtained by Fermi and VERITAS observatories. The NuSTAR observation allowed us to detect a 51.67 ms spin period from the pulsar PSR J2229+6114 and the PWN emission characterized by a power-law model with $Γ= 1.52\pm0.06$ up to 20 keV. Contrary to the previous radio study by Kothes et al. 2006, we prefer a much lower PWN B-field ($B\sim3$ $μ$G) and larger distance ($d \sim 8$ kpc) based on (1) the non-varying X-ray flux over the last two decades, (2) the energy-dependent X-ray PWN size resulting from synchrotron burn-off and (3) the multi-wavelength spectral energy distribution (SED) data. Our SED model suggests that the PWN is currently re-expanding after being compressed by the SNR reverse shock $\sim 1000$ years ago. In this case, the head region should be formed by GeV--TeV electrons injected earlier by the pulsar propagating into the low density environment. △ Less

Submitted 6 October, 2023; originally announced October 2023.

LSOR: Longitudinally-Consistent Self-Organized Representation Learning

Authors: Jiahong Ouyang, Qingyu Zhao, Ehsan Adeli, Wei Peng, Greg Zaharchuk, Kilian M. Pohl

Abstract: Interpretability is a key issue when applying deep learning models to longitudinal brain MRIs. One way to address this issue is by visualizing the high-dimensional latent spaces generated by deep learning via self-organizing maps (SOM). SOM separates the latent space into clusters and then maps the cluster centers to a discrete (typically 2D) grid preserving the high-dimensional relationship betwe… ▽ More Interpretability is a key issue when applying deep learning models to longitudinal brain MRIs. One way to address this issue is by visualizing the high-dimensional latent spaces generated by deep learning via self-organizing maps (SOM). SOM separates the latent space into clusters and then maps the cluster centers to a discrete (typically 2D) grid preserving the high-dimensional relationship between clusters. However, learning SOM in a high-dimensional latent space tends to be unstable, especially in a self-supervision setting. Furthermore, the learned SOM grid does not necessarily capture clinically interesting information, such as brain age. To resolve these issues, we propose the first self-supervised SOM approach that derives a high-dimensional, interpretable representation stratified by brain age solely based on longitudinal brain MRIs (i.e., without demographic or cognitive information). Called Longitudinally-consistent Self-Organized Representation learning (LSOR), the method is stable during training as it relies on soft clustering (vs. the hard cluster assignments used by existing SOM). Furthermore, our approach generates a latent space stratified according to brain age by aligning trajectories inferred from longitudinal MRIs to the reference vector associated with the corresponding SOM cluster. When applied to longitudinal MRIs of the Alzheimer's Disease Neuroimaging Initiative (ADNI, N=632), LSOR generates an interpretable latent space and achieves comparable or higher accuracy than the state-of-the-art representations with respect to the downstream tasks of classification (static vs. progressive mild cognitive impairment) and regression (determining ADAS-Cog score of all subjects). The code is available at https://github.com/ouyangjiahong/longitudinal-som-single-modality. △ Less

Submitted 29 September, 2023; originally announced October 2023.

Journal ref: International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI) 2023

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Imputing Brain Measurements Across Data Sets via Graph Neural Networks

Authors: Yixin Wang, Wei Peng, Susan F. Tapert, Qingyu Zhao, Kilian M. Pohl

Abstract: Publicly available data sets of structural MRIs might not contain specific measurements of brain Regions of Interests (ROIs) that are important for training machine learning models. For example, the curvature scores computed by Freesurfer are not released by the Adolescent Brain Cognitive Development (ABCD) Study. One can address this issue by simply reapplying Freesurfer to the data set. However,… ▽ More Publicly available data sets of structural MRIs might not contain specific measurements of brain Regions of Interests (ROIs) that are important for training machine learning models. For example, the curvature scores computed by Freesurfer are not released by the Adolescent Brain Cognitive Development (ABCD) Study. One can address this issue by simply reapplying Freesurfer to the data set. However, this approach is generally computationally and labor intensive (e.g., requiring quality control). An alternative is to impute the missing measurements via a deep learning approach. However, the state-of-the-art is designed to estimate randomly missing values rather than entire measurements. We therefore propose to re-frame the imputation problem as a prediction task on another (public) data set that contains the missing measurements and shares some ROI measurements with the data sets of interest. A deep learning model is then trained to predict the missing measurements from the shared ones and afterwards is applied to the other data sets. Our proposed algorithm models the dependencies between ROI measurements via a graph neural network (GNN) and accounts for demographic differences in brain measurements (e.g. sex) by feeding the graph encoding into a parallel architecture. The architecture simultaneously optimizes a graph decoder to impute values and a classifier in predicting demographic factors. We test the approach, called Demographic Aware Graph-based Imputation (DAGI), on imputing those missing Freesurfer measurements of ABCD (N=3760) by training the predictor on those publicly released by the National Consortium on Alcohol and Neurodevelopment in Adolescence (NCANDA, N=540)... △ Less

Submitted 19 August, 2023; originally announced August 2023.

Comments: Accepted at the 6th workshop on PRedictive Intelligence in Medicine (PRIME 2023) - MICCAI 2023

Multiwavelength Observations of the Blazar PKS~0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

Authors: M. Pohl, K. Mori

Abstract: We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2 degrees away from the best-fit position of the IceCube neutrino event 211208A. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT… ▽ More We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2 degrees away from the best-fit position of the IceCube neutrino event 211208A. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution, and the gamma-ray data from Fermi-LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both measurements provide strong constraints on leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed spectral energy distribution (SED). The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models, and it would allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. A numerical lepto-hadronic model with external target photons reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: Proceedings paper presented at the 38th International Cosmic Ray Conference (ICRC2023), held 26 July - 3 August, 2023, in Nagoya, Japan

An Explainable Geometric-Weighted Graph Attention Network for Identifying Functional Networks Associated with Gait Impairment

Authors: Favour Nerrise, Qingyu Zhao, Kathleen L. Poston, Kilian M. Pohl, Ehsan Adeli

Abstract: One of the hallmark symptoms of Parkinson's Disease (PD) is the progressive loss of postural reflexes, which eventually leads to gait difficulties and balance problems. Identifying disruptions in brain function associated with gait impairment could be crucial in better understanding PD motor progression, thus advancing the development of more effective and personalized therapeutics. In this work,… ▽ More One of the hallmark symptoms of Parkinson's Disease (PD) is the progressive loss of postural reflexes, which eventually leads to gait difficulties and balance problems. Identifying disruptions in brain function associated with gait impairment could be crucial in better understanding PD motor progression, thus advancing the development of more effective and personalized therapeutics. In this work, we present an explainable, geometric, weighted-graph attention neural network (xGW-GAT) to identify functional networks predictive of the progression of gait difficulties in individuals with PD. xGW-GAT predicts the multi-class gait impairment on the MDS Unified PD Rating Scale (MDS-UPDRS). Our computational- and data-efficient model represents functional connectomes as symmetric positive definite (SPD) matrices on a Riemannian manifold to explicitly encode pairwise interactions of entire connectomes, based on which we learn an attention mask yielding individual- and group-level explainability. Applied to our resting-state functional MRI (rs-fMRI) dataset of individuals with PD, xGW-GAT identifies functional connectivity patterns associated with gait impairment in PD and offers interpretable explanations of functional subnetworks associated with motor impairment. Our model successfully outperforms several existing methods while simultaneously revealing clinically-relevant connectivity patterns. The source code is available at https://github.com/favour-nerrise/xGW-GAT . △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: Accepted by the 26th International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI 2023). MICCAI Student-Author Registration (STAR) Award. 11 pages, 2 figures, 1 table, appendix. Source Code: https://github.com/favour-nerrise/xGW-GAT

Multiwavelength Observations of the Blazar PKS 0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. Brill, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, T. B. Humensky , et al. (185 additional authors not shown)

Abstract: We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ra… ▽ More We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution (SED), and the gamma-ray data from Fermi -LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both X-ray and gamma-ray measurements provide strong constraints on the leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed SED. The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models and allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. We show a numerical lepto-hadronic model with external target photons that reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 21 pages, 3 figures, accepted by ApJ

A VERITAS/Breakthrough Listen Search for Optical Technosignatures

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Pedro Batista, Wystan Benbow, Aryeh Brill, M Capasso, Manel Errando, Abraham Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Sean Griffin, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, T. Humensky, Weidong Jin, Philip Kaaret , et al. (43 additional authors not shown)

Abstract: The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectabl… ▽ More The Breakthrough Listen Initiative is conducting a program using multiple telescopes around the world to search for "technosignatures": artificial transmitters of extraterrestrial origin from beyond our solar system. The VERITAS Collaboration joined this program in 2018, and provides the capability to search for one particular technosignature: optical pulses of a few nanoseconds duration detectable over interstellar distances. We report here on the analysis and results of dedicated VERITAS observations of Breakthrough Listen targets conducted in 2019 and 2020 and of archival VERITAS data collected since 2012. Thirty hours of dedicated observations of 136 targets and 249 archival observations of 140 targets were analyzed and did not reveal any signals consistent with a technosignature. The results are used to place limits on the fraction of stars hosting transmitting civilizations. We also discuss the minimum-pulse sensitivity of our observations and present VERITAS observations of CALIOP: a space-based pulsed laser onboard the CALIPSO satellite. The detection of these pulses with VERITAS, using the analysis techniques developed for our technosignature search, allows a test of our analysis efficiency and serves as an important proof-of-principle. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 15 pages, 7 figures

Evolution of Relativistic Pair Beams: Implications for Laboratory and TeV Astrophysics

Authors: Marvin Beck, Oindrila Ghosh, Florian Grüner, Martin Pohl, Carl B. Schroeder, Günter Sigl, Ryan D. Stark, Benno Zeitler

Abstract: Missing cascades from TeV blazar beams indicate that collective plasma effects may play a significant role in their energy loss. It is possible to mimic the evolution of such highly energetic pair beams in laboratory experiments using modern accelerators. The fate of the beam is governed by two different processes, energy loss through the unstable mode and energetic broadening of the pair beam thr… ▽ More Missing cascades from TeV blazar beams indicate that collective plasma effects may play a significant role in their energy loss. It is possible to mimic the evolution of such highly energetic pair beams in laboratory experiments using modern accelerators. The fate of the beam is governed by two different processes, energy loss through the unstable mode and energetic broadening of the pair beam through diffusion in momentum space. We chalk out this evolution using a Fokker-Planck approach in which the drift and the diffusion terms respectively describe these phenomena in a compact form. We present particle-in-cell simulations to trace the complete evolution of the unstable beam-plasma system for a generic narrow Gaussian pair beam for which the growth rate is reactive. We show that the instability leads to an energetic broadening of the pair beam, slowing down the instability growth in the linear phase, in line with the analytical and numerical solutions of the Fokker-Planck equation. Whereas in a laboratory experiment the change in the momentum distribution is an easily measured observable as a feedback of the instability, the consequence of diffusive broadening in an astrophysical scenario can be translated to an increase in the opening angle of the pair beam. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 20 pages, 10 figures

VERITAS discovery of very high energy gamma-ray emission from S3 1227+25 and multiwavelength observations

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Wystan Benbow, Aryeh Brill, Jodi Christiansen, Alisha Chromey, Manel Errando, Abe Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Greg Gallagher, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, Weidong Jin, Madalyn Johnson, Philip Kaaret , et al. (46 additional authors not shown)

Abstract: We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th… ▽ More We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th resulted in a strong 13$σ$ detection with a differential photon spectral index, $Γ$ = 3.8 $\pm$ 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target of opportunity observations with Swift, optical photometry, polarimetry and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of $τ_{obs}$ = 6.2 $\pm$ 0.9 hours, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cut-off. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 18 pages, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

Electron-Ion Temperature Ratio in Astrophysical Shocks

Authors: John C. Raymond, Parviz Ghavamian, Artem Bohdan, Dongsu Ryu, Jacek Niemiec, Lorenzo Sironi, Aaron Tran, Elena Amato, Masahiro Hoshino, Martin Pohl, Takanobu Amano, Federico Fiuza

Abstract: Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Her… ▽ More Collisionless shock waves in supernova remnants and the solar wind heat electrons less effectively than they heat ions, as is predicted by kinetic simulations. However, the values of T$_e$/T$_p$ inferred from the H alpha profiles of supernova remnant shocks behave differently as a function of Mach number or Alfvén Mach number than what is measured in the solar wind or predicted by simulations. Here we determine T$_e$/T$_p$ for supernova remnant shocks using H alpha profiles, shock speeds from proper motions, and electron temperatures from X-ray spectra. We also improve the estimates of sound speed and Alfvén speed used to determine Mach numbers. We find that the H alpha determinations are robust and that the discrepancies among supernova remnant shocks, solar wind shocks and computer-simulated shocks remain. We discuss some possible contributing factors, including shock precursors, turbulence and varying preshock conditions. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Search for Ultraheavy Dark Matter from Observations of Dwarf Spheroidal Galaxies with VERITAS

Authors: A. Acharyya, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, M. Baumgart, W. Benbow, J. H. Buckley, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. F. Hanlon, O. Hervet, J. Hoang, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (29 additional authors not shown)

Abstract: Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has bee… ▽ More Dark matter is a key piece of the current cosmological scenario, with weakly interacting massive particles (WIMPs) a leading dark matter candidate. WIMPs have not been detected in their conventional parameter space (100 GeV $\lesssim M_χ \lesssim$ 100 TeV), a mass range accessible with current Imaging Atmospheric Cherenkov Telescopes. As ultraheavy dark matter (UHDM; $M_χ \gtrsim$ 100 TeV) has been suggested as an under-explored alternative to the WIMP paradigm, we search for an indirect dark matter annihilation signal in a higher mass range (up to 30 PeV) with the VERITAS gamma-ray observatory. With 216 hours of observations of four dwarf spheroidal galaxies, we perform an unbinned likelihood analysis. We find no evidence of a $γ$-ray signal from UHDM annihilation above the background fluctuation for any individual dwarf galaxy nor for a joint-fit analysis, and consequently constrain the velocity-weighted annihilation cross section of UHDM for dark matter particle masses between 1 TeV and 30 PeV. We additionally set constraints on the allowed radius of a composite UHDM particle. △ Less

Submitted 24 February, 2023; v1 submitted 17 February, 2023; originally announced February 2023.

Comments: 10 pages, 7 figures

VERITAS and Fermi-LAT constraints on the Gamma-ray Emission from Superluminous Supernovae SN2015bn and SN2017egm

Authors: A. Acharyya, C. B. Adams, P. Bangale, W. Benbow, J. H. Buckley, M. Capasso, V. V. Dwarkadas, M. Errando, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, W. F Hanlon, O. Hervet, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (34 additional authors not shown)

Abstract: Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnet… ▽ More Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar's spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV - 100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 20 pages, 7 figures, 2 tables

VTSCat: The VERITAS Catalog of Gamma-Ray Observations

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, J. H. Buckley, A. Brill, M. Capasso, J. L. Christiansen, A. J. Chromey, M. K. Daniel, M. Errando, A. Falcone, K. A Farrell, Q. Feng, J. P. Finley, G. M Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, C. Giuri, O. Gueta , et al. (64 additional authors not shown)

Abstract: The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four… ▽ More The ground-based gamma-ray observatory VERITAS (Very Energetic Radiation Imaging Telescope Array System) is sensitive to photons of astrophysical origin with energies in the range between $\approx 85$ GeV to $\approx 30$ TeV. The instrument consists of four 12-m diameter imaging Cherenkov telescopes operating at the Fred Lawrence Whipple Observatory (FLWO) in southern Arizona. VERITAS started four-telescope operations in 2007 and collects about 1100 hours of good-weather data per year. The VERITAS collaboration has published over 100 journal articles since 2008 reporting on gamma-ray observations of a large variety of objects: Galactic sources like supernova remnants, pulsar wind nebulae, and binary systems; extragalactic sources like star forming galaxies, dwarf-spheroidal galaxies, and highly-variable active galactic nuclei. This note presents VTSCat: the catalog of high-level data products from all VERITAS publications. △ Less

Submitted 13 January, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: Version with corrected author list

Journal ref: Research Notes of the AAS, Volume 7, Number 1, 2023

Pre-acceleration in the Electron Foreshock II: Oblique Whistler Waves

Authors: Paul J. Morris, Artem Bohdan, Martin S. Weidl, Michelle Tsirou, Karol Fulat, Martin Pohl

Abstract: Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large scal… ▽ More Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large scale particle-in-cell simulation that generates a single shock with parameters motivated from supernova remnants. Within the foreshock, reflected electrons excite the oblique whistler instability and produce electromagnetic whistler waves, which co-move with the upstream flow and as non-linear structures eventually reach radii of up to 5 ion-gyroradii. We show that the inner electromagnetic configuration of the whistlers evolves into complex non-linear structures bound by a strong magnetic field around 4 times the upstream value. Although these non-linear structures do not in general interact with co-spatial upstream electrons, they resonate with electrons that have been reflected at the shock. We show that they can scatter, or even trap, reflected electrons, confining around $0.8\%$ of the total upstream electron population to the region close to the shock where they can undergo substantial pre-acceleration. This acceleration process is similar to, yet approximately 3 times more efficient than, stochastic shock drift acceleration. △ Less

Submitted 2 January, 2023; originally announced January 2023.

Comments: 16 pages, 11 figures, accepted for publication in ApJ

Generating Realistic Brain MRIs via a Conditional Diffusion Probabilistic Model

Authors: Wei Peng, Ehsan Adeli, Tomas Bosschieter, Sang Hyun Park, Qingyu Zhao, Kilian M. Pohl

Abstract: As acquiring MRIs is expensive, neuroscience studies struggle to attain a sufficient number of them for properly training deep learning models. This challenge could be reduced by MRI synthesis, for which Generative Adversarial Networks (GANs) are popular. GANs, however, are commonly unstable and struggle with creating diverse and high-quality data. A more stable alternative is Diffusion Probabilis… ▽ More As acquiring MRIs is expensive, neuroscience studies struggle to attain a sufficient number of them for properly training deep learning models. This challenge could be reduced by MRI synthesis, for which Generative Adversarial Networks (GANs) are popular. GANs, however, are commonly unstable and struggle with creating diverse and high-quality data. A more stable alternative is Diffusion Probabilistic Models (DPMs) with a fine-grained training strategy. To overcome their need for extensive computational resources, we propose a conditional DPM (cDPM) with a memory-efficient process that generates realistic-looking brain MRIs. To this end, we train a 2D cDPM to generate an MRI subvolume conditioned on another subset of slices from the same MRI. By generating slices using arbitrary combinations between condition and target slices, the model only requires limited computational resources to learn interdependencies between slices even if they are spatially far apart. After having learned these dependencies via an attention network, a new anatomy-consistent 3D brain MRI is generated by repeatedly applying the cDPM. Our experiments demonstrate that our method can generate high-quality 3D MRIs that share a similar distribution to real MRIs while still diversifying the training set. The code is available at https://github.com/xiaoiker/mask3DMRI_diffusion and also will be released as part of MONAI, at https://github.com/Project-MONAI/GenerativeModels. △ Less

Submitted 7 September, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Journal ref: MICCAI 2023

Joint Graph Convolution for Analyzing Brain Structural and Functional Connectome

Authors: Yueting Li, Qingyue Wei, Ehsan Adeli, Kilian M. Pohl, Qingyu Zhao

Abstract: The white-matter (micro-)structural architecture of the brain promotes synchrony among neuronal populations, giving rise to richly patterned functional connections. A fundamental problem for systems neuroscience is determining the best way to relate structural and functional networks quantified by diffusion tensor imaging and resting-state functional MRI. As one of the state-of-the-art approaches… ▽ More The white-matter (micro-)structural architecture of the brain promotes synchrony among neuronal populations, giving rise to richly patterned functional connections. A fundamental problem for systems neuroscience is determining the best way to relate structural and functional networks quantified by diffusion tensor imaging and resting-state functional MRI. As one of the state-of-the-art approaches for network analysis, graph convolutional networks (GCN) have been separately used to analyze functional and structural networks, but have not been applied to explore inter-network relationships. In this work, we propose to couple the two networks of an individual by adding inter-network edges between corresponding brain regions, so that the joint structure-function graph can be directly analyzed by a single GCN. The weights of inter-network edges are learnable, reflecting non-uniform structure-function coupling strength across the brain. We apply our Joint-GCN to predict age and sex of 662 participants from the public dataset of the National Consortium on Alcohol and Neurodevelopment in Adolescence (NCANDA) based on their functional and micro-structural white-matter networks. Our results support that the proposed Joint-GCN outperforms existing multi-modal graph learning approaches for analyzing structural and functional networks. △ Less

Submitted 27 October, 2022; originally announced November 2022.

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

Authors: H. Abe, S. Abe, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari , et al. (418 additional authors not shown)

Abstract: MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS expe… ▽ More MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA. △ Less

Submitted 6 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Nonthermal Electron Acceleration at Collisionless Quasi-perpendicular Shocks

Authors: T. Amano, Y. Matsumoto, A. Bohdan, O. Kobzar, S. Matsukiyo, M. Oka, J. Niemiec, M. Pohl, M. Hoshino

Abstract: Shock waves propagating in collisionless heliospheric and astrophysical plasmas have been studied extensively over the decades. One prime motivation is to understand the nonthermal particle acceleration at shocks. Although the theory of diffusive shock acceleration (DSA) has long been the standard for cosmic-ray acceleration at shocks, plasma physical understanding of particle acceleration remains… ▽ More Shock waves propagating in collisionless heliospheric and astrophysical plasmas have been studied extensively over the decades. One prime motivation is to understand the nonthermal particle acceleration at shocks. Although the theory of diffusive shock acceleration (DSA) has long been the standard for cosmic-ray acceleration at shocks, plasma physical understanding of particle acceleration remains elusive. In this review, we discuss nonthermal electron acceleration mechanisms at quasi-perpendicular shocks, for which substantial progress has been made in recent years. The discussion presented in this review is restricted to the following three specific topics. The first is stochastic shock drift acceleration (SSDA), which is a relatively new mechanism for electron injection into DSA. The basic mechanism, related in-situ observations and kinetic simulations results, and how it is connected with DSA will be discussed. Second, we discuss shock surfing acceleration (SSA) at very high Mach number shocks relevant to young supernova remnants (SNRs). While the original proposal under the one-dimensional assumption is unrealistic, SSA has now been proven efficient by a fully three-dimensional kinetic simulation. Finally, we discuss the current understanding of the magnetized Weibel-dominated shock. Spontaneous magnetic reconnection of self-generated current sheets within the shock structure is an interesting consequence of Weibel-generated strong magnetic turbulence. We argue that high Mach number shocks with both Alfven and sound Mach numbers exceeding 20-40 will likely behave as a Weibel-dominated shock. Despite a number of interesting recent findings, the relative roles of SSDA, SSA, and magnetic reconnection for electron acceleration at collisionless shocks and how the dominant particle acceleration mechanisms change depending on shock parameters remain to be answered. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Comments: To appear in Reviews of Modern Plasma Physics as an invited review

Multiple Instance Neuroimage Transformer

Authors: Ayush Singla, Qingyu Zhao, Daniel K. Do, Yuyin Zhou, Kilian M. Pohl, Ehsan Adeli

Abstract: For the first time, we propose using a multiple instance learning based convolution-free transformer model, called Multiple Instance Neuroimage Transformer (MINiT), for the classification of T1weighted (T1w) MRIs. We first present several variants of transformer models adopted for neuroimages. These models extract non-overlapping 3D blocks from the input volume and perform multi-headed self-attent… ▽ More For the first time, we propose using a multiple instance learning based convolution-free transformer model, called Multiple Instance Neuroimage Transformer (MINiT), for the classification of T1weighted (T1w) MRIs. We first present several variants of transformer models adopted for neuroimages. These models extract non-overlapping 3D blocks from the input volume and perform multi-headed self-attention on a sequence of their linear projections. MINiT, on the other hand, treats each of the non-overlapping 3D blocks of the input MRI as its own instance, splitting it further into non-overlapping 3D patches, on which multi-headed self-attention is computed. As a proof-of-concept, we evaluate the efficacy of our model by training it to identify sex from T1w-MRIs of two public datasets: Adolescent Brain Cognitive Development (ABCD) and the National Consortium on Alcohol and Neurodevelopment in Adolescence (NCANDA). The learned attention maps highlight voxels contributing to identifying sex differences in brain morphometry. The code is available at https://github.com/singlaayush/MINIT. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Modeling kilonova afterglows: Effects of the thermal electron population and interaction with GRB outflows

Authors: Vsevolod Nedora, Tim Dietrich, Masaru Shibata, Martin Pohl, Ludovica Crosato Menegazzi

Abstract: Given an increasing number of gamma-ray bursts accompanied by potential kilonovae there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distributions affect kilonova afterglow light curves. We employ semi-analytic afterglow model,… ▽ More Given an increasing number of gamma-ray bursts accompanied by potential kilonovae there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distributions affect kilonova afterglow light curves. We employ semi-analytic afterglow model, $\texttt{PyBlastAfterglow}$. We consider kilonova ejecta profiles from ab-initio numerical relativity binary neutron star merger simulations, targeted to GW170817. We do not perform model selection. We find that the emission from thermal electrons dominates at early times. If the interstellar medium density is high (${\simeq}0.1\,\ccm$) it adds an early time peak to the light curve. As ejecta decelerates the spectral and temporal indexes change in a characteristic way that, if observed, can be used to reconstruct the ejecta velocity distribution. For the low interstellar medium density, inferred for GRB 170817A, the emission from the non-thermal electron population generally dominates. We also assess how kilonova afterglow light curves change if the interstellar medium has been partially removed and pre-accelerated by laterally expanding gamma-ray burst ejecta. We find that the main effect is the emission suppression at early time ${\lesssim}10^{3}\,$days, and at its maximum it reaches ${\sim}40\%$ when the fast tail of the kilonova ejecta moves subsonically through the wake of laterally spreading gamma-ray burst ejecta. The subsequent rebrightening, when these ejecta break through and shocks form, is very mild (${\lesssim}10\%$), and may not be observable. △ Less

Submitted 13 January, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 20 pages, 14 figures

Bridging the Gap between Deep Learning and Hypothesis-Driven Analysis via Permutation Testing

Authors: Magdalini Paschali, Qingyu Zhao, Ehsan Adeli, Kilian M. Pohl

Abstract: A fundamental approach in neuroscience research is to test hypotheses based on neuropsychological and behavioral measures, i.e., whether certain factors (e.g., related to life events) are associated with an outcome (e.g., depression). In recent years, deep learning has become a potential alternative approach for conducting such analyses by predicting an outcome from a collection of factors and ide… ▽ More A fundamental approach in neuroscience research is to test hypotheses based on neuropsychological and behavioral measures, i.e., whether certain factors (e.g., related to life events) are associated with an outcome (e.g., depression). In recent years, deep learning has become a potential alternative approach for conducting such analyses by predicting an outcome from a collection of factors and identifying the most "informative" ones driving the prediction. However, this approach has had limited impact as its findings are not linked to statistical significance of factors supporting hypotheses. In this article, we proposed a flexible and scalable approach based on the concept of permutation testing that integrates hypothesis testing into the data-driven deep learning analysis. We apply our approach to the yearly self-reported assessments of 621 adolescent participants of the National Consortium of Alcohol and Neurodevelopment in Adolescence (NCANDA) to predict negative valence, a symptom of major depressive disorder according to the NIMH Research Domain Criteria (RDoC). Our method successfully identifies categories of risk factors that further explain the symptom. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: Accepted at the 5th workshop on PRedictive Intelligence in Medicine (PRIME 2022) - MICCAI 2022

Prediction of the motion of chest internal points using a recurrent neural network trained with real-time recurrent learning for latency compensation in lung cancer radiotherapy

Authors: Michel Pohl, Mitsuru Uesaka, Kazuyuki Demachi, Ritu Bhusal Chhatkuli

Abstract: During the radiotherapy treatment of patients with lung cancer, the radiation delivered to healthy tissue around the tumor needs to be minimized, which is difficult because of respiratory motion and the latency of linear accelerator systems. In the proposed study, we first use the Lucas-Kanade pyramidal optical flow algorithm to perform deformable image registration of chest computed tomography sc… ▽ More During the radiotherapy treatment of patients with lung cancer, the radiation delivered to healthy tissue around the tumor needs to be minimized, which is difficult because of respiratory motion and the latency of linear accelerator systems. In the proposed study, we first use the Lucas-Kanade pyramidal optical flow algorithm to perform deformable image registration of chest computed tomography scan images of four patients with lung cancer. We then track three internal points close to the lung tumor based on the previously computed deformation field and predict their position with a recurrent neural network (RNN) trained using real-time recurrent learning (RTRL) and gradient clipping. The breathing data is quite regular, sampled at approximately 2.5Hz, and includes artificial drift in the spine direction. The amplitude of the motion of the tracked points ranged from 12.0mm to 22.7mm. Finally, we propose a simple method for recovering and predicting 3D tumor images from the tracked points and the initial tumor image based on a linear correspondence model and Nadaraya-Watson non-linear regression. The root-mean-square error, maximum error, and jitter corresponding to the RNN prediction on the test set were smaller than the same performance measures obtained with linear prediction and least mean squares (LMS). In particular, the maximum prediction error associated with the RNN, equal to 1.51mm, is respectively 16.1% and 5.0% lower than the maximum error associated with linear prediction and LMS. The average prediction time per time step with RTRL is equal to 119ms, which is less than the 400ms marker position sampling time. The tumor position in the predicted images appears visually correct, which is confirmed by the high mean cross-correlation between the original and predicted images, equal to 0.955. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: 21 pages, 24 figures; accepted manuscript version

Journal ref: Computerized Medical Imaging and Graphics, Volume 91, 2021, p.101941