Showing 1–50 of 437 results for author: Lang, M

Persistence of small polarons into the superconducting phase of Ba$_{1-x}$K$_x$BiO$_3$

Authors: Muntaser Naamneh, Eugenio Paris, Daniel McNally, Yi Tseng, Wojciech R. Pudelko, Dariusz J. Gawryluk, J. Shamblin, Eric OQuinn, Benjamin Cohen-Stead, Ming Shi, Milan Radovic, M. Lang, Thorsten Schmitt, Steven Johnston, Nicholas C. Plumb

Abstract: Bipolaronic superconductivity is an exotic pairing mechanism proposed for materials like Ba$_{1-x}$K$_x$BiO$_3$ (BKBO); however, conclusive experimental evidence for a (bi)polaron metallic state in this material remains elusive. Here, we combine resonant inelastic x-ray and neutron total scattering techniques with advanced modelling to study the local lattice distortions, electronic structure, and… ▽ More Bipolaronic superconductivity is an exotic pairing mechanism proposed for materials like Ba$_{1-x}$K$_x$BiO$_3$ (BKBO); however, conclusive experimental evidence for a (bi)polaron metallic state in this material remains elusive. Here, we combine resonant inelastic x-ray and neutron total scattering techniques with advanced modelling to study the local lattice distortions, electronic structure, and electron-phonon coupling ($e$-ph) in BKBO as a function of doping. Data for the parent compound ($x = 0$) indicates that the electronic gap opens in predominantly oxygen-derived states strongly coupled to a long-range ordered breathing distortion of the oxygen sublattice. Upon doping, short-range breathing distortions and sizable ($e$-ph) coupling persist into the superconducting regime ($x = 0.4$). Comparisons with exact diagonalization and determinant quantum Monte Carlo calculations further support this conclusion. Our results provide compelling evidence that BKBO's metallic phase hosts a liquid of small (bi)polarons derived from local breathing distortions of the lattice, with implications for understanding the low-temperature superconducting instability △ Less

Submitted 1 August, 2024; originally announced August 2024.

Who Let the Guards Out: Visual Support for Patrolling Games

Authors: Matěj Lang, Adam Štěpánek, Róbert Zvara, Vojtěch Řehák, Barbora Kozlíková

Abstract: Effective security patrol management is critical for ensuring safety in diverse environments such as art galleries, airports, and factories. The behavior of patrols in these situations can be modeled by patrolling games. They simulate the behavior of the patrol and adversary in the building, which is modeled as a graph of interconnected nodes representing rooms. The designers of algorithms solving… ▽ More Effective security patrol management is critical for ensuring safety in diverse environments such as art galleries, airports, and factories. The behavior of patrols in these situations can be modeled by patrolling games. They simulate the behavior of the patrol and adversary in the building, which is modeled as a graph of interconnected nodes representing rooms. The designers of algorithms solving the game face the problem of analyzing complex graph layouts with temporal dependencies. Therefore, appropriate visual support is crucial for them to work effectively. In this paper, we present a novel tool that helps the designers of patrolling games explore the outcomes of the proposed algorithms and approaches, evaluate their success rate, and propose modifications that can improve their solutions. Our tool offers an intuitive and interactive interface, featuring a detailed exploration of patrol routes and probabilities of taking them, simulation of patrols, and other requested features. In close collaboration with experts in designing patrolling games, we conducted three case studies demonstrating the usage and usefulness of our tool. The prototype of the tool, along with exemplary datasets, is available at https://gitlab.fi.muni.cz/formela/strategy-vizualizer. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Graph Neural Networks: A suitable Alternative to MLPs in Latent 3D Medical Image Classification?

Authors: Johannes Kiechle, Daniel M. Lang, Stefan M. Fischer, Lina Felsner, Jan C. Peeken, Julia A. Schnabel

Abstract: Recent studies have underscored the capabilities of natural imaging foundation models to serve as powerful feature extractors, even in a zero-shot setting for medical imaging data. Most commonly, a shallow multi-layer perceptron (MLP) is appended to the feature extractor to facilitate end-to-end learning and downstream prediction tasks such as classification, thus representing the de facto standar… ▽ More Recent studies have underscored the capabilities of natural imaging foundation models to serve as powerful feature extractors, even in a zero-shot setting for medical imaging data. Most commonly, a shallow multi-layer perceptron (MLP) is appended to the feature extractor to facilitate end-to-end learning and downstream prediction tasks such as classification, thus representing the de facto standard. However, as graph neural networks (GNNs) have become a practicable choice for various tasks in medical research in the recent past, we direct attention to the question of how effective GNNs are compared to MLP prediction heads for the task of 3D medical image classification, proposing them as a potential alternative. In our experiments, we devise a subject-level graph for each volumetric dataset instance. Therein latent representations of all slices in the volume, encoded through a DINOv2 pretrained vision transformer (ViT), constitute the nodes and their respective node features. We use public datasets to compare the classification heads numerically and evaluate various graph construction and graph convolution methods in our experiments. Our findings show enhancements of the GNN in classification performance and substantial improvements in runtime compared to an MLP prediction head. Additional robustness evaluations further validate the promising performance of the GNN, promoting them as a suitable alternative to traditional MLP classification heads. Our code is publicly available at: https://github.com/compai-lab/2024-miccai-grail-kiechle △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: Accepted at MICCAI 2024 - GRAIL Workshop

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 23 July, 2024; originally announced July 2024.

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

MedEdit: Counterfactual Diffusion-based Image Editing on Brain MRI

Authors: Malek Ben Alaya, Daniel M. Lang, Benedikt Wiestler, Julia A. Schnabel, Cosmin I. Bercea

Abstract: Denoising diffusion probabilistic models enable high-fidelity image synthesis and editing. In biomedicine, these models facilitate counterfactual image editing, producing pairs of images where one is edited to simulate hypothetical conditions. For example, they can model the progression of specific diseases, such as stroke lesions. However, current image editing techniques often fail to generate r… ▽ More Denoising diffusion probabilistic models enable high-fidelity image synthesis and editing. In biomedicine, these models facilitate counterfactual image editing, producing pairs of images where one is edited to simulate hypothetical conditions. For example, they can model the progression of specific diseases, such as stroke lesions. However, current image editing techniques often fail to generate realistic biomedical counterfactuals, either by inadequately modeling indirect pathological effects like brain atrophy or by excessively altering the scan, which disrupts correspondence to the original images. Here, we propose MedEdit, a conditional diffusion model for medical image editing. MedEdit induces pathology in specific areas while balancing the modeling of disease effects and preserving the integrity of the original scan. We evaluated MedEdit on the Atlas v2.0 stroke dataset using Frechet Inception Distance and Dice scores, outperforming state-of-the-art diffusion-based methods such as Palette (by 45%) and SDEdit (by 61%). Additionally, clinical evaluations by a board-certified neuroradiologist confirmed that MedEdit generated realistic stroke scans indistinguishable from real ones. We believe this work will enable counterfactual image editing research to further advance the development of realistic and clinically useful imaging tools. △ Less

Submitted 21 July, 2024; originally announced July 2024.

Comments: Accepted at MICCAI24 Simulation and Synthesis in Medical Imaging (SASHIMI) workshop

Enhancing the Utility of Privacy-Preserving Cancer Classification using Synthetic Data

Authors: Richard Osuala, Daniel M. Lang, Anneliese Riess, Georgios Kaissis, Zuzanna Szafranowska, Grzegorz Skorupko, Oliver Diaz, Julia A. Schnabel, Karim Lekadir

Abstract: Deep learning holds immense promise for aiding radiologists in breast cancer detection. However, achieving optimal model performance is hampered by limitations in availability and sharing of data commonly associated to patient privacy concerns. Such concerns are further exacerbated, as traditional deep learning models can inadvertently leak sensitive training information. This work addresses these… ▽ More Deep learning holds immense promise for aiding radiologists in breast cancer detection. However, achieving optimal model performance is hampered by limitations in availability and sharing of data commonly associated to patient privacy concerns. Such concerns are further exacerbated, as traditional deep learning models can inadvertently leak sensitive training information. This work addresses these challenges exploring and quantifying the utility of privacy-preserving deep learning techniques, concretely, (i) differentially private stochastic gradient descent (DP-SGD) and (ii) fully synthetic training data generated by our proposed malignancy-conditioned generative adversarial network. We assess these methods via downstream malignancy classification of mammography masses using a transformer model. Our experimental results depict that synthetic data augmentation can improve privacy-utility tradeoffs in differentially private model training. Further, model pretraining on synthetic data achieves remarkable performance, which can be further increased with DP-SGD fine-tuning across all privacy guarantees. With this first in-depth exploration of privacy-preserving deep learning in breast imaging, we address current and emerging clinical privacy requirements and pave the way towards the adoption of private high-utility deep diagnostic models. Our reproducible codebase is publicly available at https://github.com/RichardObi/mammo_dp. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Comments: Early Accept at MICCAI 2024 Deep-Breath Workshop

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, G. M. 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 17 July, 2024; v1 submitted 16 July, 2024; originally announced July 2024.

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

Progressive Growing of Patch Size: Resource-Efficient Curriculum Learning for Dense Prediction Tasks

Authors: Stefan M. Fischer, Lina Felsner, Richard Osuala, Johannes Kiechle, Daniel M. Lang, Jan C. Peeken, Julia A. Schnabel

Abstract: In this work, we introduce Progressive Growing of Patch Size, a resource-efficient implicit curriculum learning approach for dense prediction tasks. Our curriculum approach is defined by growing the patch size during model training, which gradually increases the task's difficulty. We integrated our curriculum into the nnU-Net framework and evaluated the methodology on all 10 tasks of the Medical S… ▽ More In this work, we introduce Progressive Growing of Patch Size, a resource-efficient implicit curriculum learning approach for dense prediction tasks. Our curriculum approach is defined by growing the patch size during model training, which gradually increases the task's difficulty. We integrated our curriculum into the nnU-Net framework and evaluated the methodology on all 10 tasks of the Medical Segmentation Decathlon. With our approach, we are able to substantially reduce runtime, computational costs, and CO2 emissions of network training compared to classical constant patch size training. In our experiments, the curriculum approach resulted in improved convergence. We are able to outperform standard nnU-Net training, which is trained with constant patch size, in terms of Dice Score on 7 out of 10 MSD tasks while only spending roughly 50% of the original training runtime. To the best of our knowledge, our Progressive Growing of Patch Size is the first successful employment of a sample-length curriculum in the form of patch size in the field of computer vision. Our code is publicly available at https://github.com/compai-lab/2024-miccai-fischer. △ Less

Submitted 11 July, 2024; v1 submitted 10 July, 2024; originally announced July 2024.

Comments: Accepted at MICCAI2024; Changes for Camera-Ready-Version for MICCAI2024 (missing in this arxiv submission): Replaced T-Test with Wilcoxon Signed Ranked Test, as DSC samples are not normally distributed => now only significant improvements and no significant decreases in performance for PGPS/PGPS+

Mask the Unknown: Assessing Different Strategies to Handle Weak Annotations in the MICCAI2023 Mediastinal Lymph Node Quantification Challenge

Authors: Stefan M. Fischer, Johannes Kiechle, Daniel M. Lang, Jan C. Peeken, Julia A. Schnabel

Abstract: Pathological lymph node delineation is crucial in cancer diagnosis, progression assessment, and treatment planning. The MICCAI 2023 Lymph Node Quantification Challenge published the first public dataset for pathological lymph node segmentation in the mediastinum. As lymph node annotations are expensive, the challenge was formed as a weakly supervised learning task, where only a subset of all lymph… ▽ More Pathological lymph node delineation is crucial in cancer diagnosis, progression assessment, and treatment planning. The MICCAI 2023 Lymph Node Quantification Challenge published the first public dataset for pathological lymph node segmentation in the mediastinum. As lymph node annotations are expensive, the challenge was formed as a weakly supervised learning task, where only a subset of all lymph nodes in the training set have been annotated. For the challenge submission, multiple methods for training on these weakly supervised data were explored, including noisy label training, loss masking of unlabeled data, and an approach that integrated the TotalSegmentator toolbox as a form of pseudo labeling in order to reduce the number of unknown voxels. Furthermore, multiple public TCIA datasets were incorporated into the training to improve the performance of the deep learning model. Our submitted model achieved a Dice score of 0.628 and an average symmetric surface distance of 5.8~mm on the challenge test set. With our submitted model, we accomplished third rank in the MICCAI2023 LNQ challenge. A finding of our analysis was that the integration of all visible, including non-pathological, lymph nodes improved the overall segmentation performance on pathological lymph nodes of the test set. Furthermore, segmentation models trained only on clinically enlarged lymph nodes, as given in the challenge scenario, could not generalize to smaller pathological lymph nodes. The code and model for the challenge submission are available at \url{https://gitlab.lrz.de/compai/MediastinalLymphNodeSegmentation}. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: Accepted for publication at the Journal of Machine Learning for Biomedical Imaging (MELBA) https://melba-journal.org/2024:008

Journal ref: Machine.Learning.for.Biomedical.Imaging. 2 (2024)

Multi-modal Transfer Learning between Biological Foundation Models

Authors: Juan Jose Garau-Luis, Patrick Bordes, Liam Gonzalez, Masa Roller, Bernardo P. de Almeida, Lorenz Hexemer, Christopher Blum, Stefan Laurent, Jan Grzegorzewski, Maren Lang, Thomas Pierrot, Guillaume Richard

Abstract: Biological sequences encode fundamental instructions for the building blocks of life, in the form of DNA, RNA, and proteins. Modeling these sequences is key to understand disease mechanisms and is an active research area in computational biology. Recently, Large Language Models have shown great promise in solving certain biological tasks but current approaches are limited to a single sequence moda… ▽ More Biological sequences encode fundamental instructions for the building blocks of life, in the form of DNA, RNA, and proteins. Modeling these sequences is key to understand disease mechanisms and is an active research area in computational biology. Recently, Large Language Models have shown great promise in solving certain biological tasks but current approaches are limited to a single sequence modality (DNA, RNA, or protein). Key problems in genomics intrinsically involve multiple modalities, but it remains unclear how to adapt general-purpose sequence models to those cases. In this work we propose a multi-modal model that connects DNA, RNA, and proteins by leveraging information from different pre-trained modality-specific encoders. We demonstrate its capabilities by applying it to the largely unsolved problem of predicting how multiple RNA transcript isoforms originate from the same gene (i.e. same DNA sequence) and map to different transcription expression levels across various human tissues. We show that our model, dubbed IsoFormer, is able to accurately predict differential transcript expression, outperforming existing methods and leveraging the use of multiple modalities. Our framework also achieves efficient transfer knowledge from the encoders pre-training as well as in between modalities. We open-source our model, paving the way for new multi-modal gene expression approaches. △ Less

Submitted 20 June, 2024; originally announced June 2024.

MSC Class: 68T07 (Primary)

TREE: Tree Regularization for Efficient Execution

Authors: Lena Schmid, Daniel Biebert, Christian Hakert, Kuan-Hsun Chen, Michel Lang, Markus Pauly, Jian-Jia Chen

Abstract: The rise of machine learning methods on heavily resource constrained devices requires not only the choice of a suitable model architecture for the target platform, but also the optimization of the chosen model with regard to execution time consumption for inference in order to optimally utilize the available resources. Random forests and decision trees are shown to be a suitable model for such a s… ▽ More The rise of machine learning methods on heavily resource constrained devices requires not only the choice of a suitable model architecture for the target platform, but also the optimization of the chosen model with regard to execution time consumption for inference in order to optimally utilize the available resources. Random forests and decision trees are shown to be a suitable model for such a scenario, since they are not only heavily tunable towards the total model size, but also offer a high potential for optimizing their executions according to the underlying memory architecture. In addition to the straightforward strategy of enforcing shorter paths through decision trees and hence reducing the execution time for inference, hardware-aware implementations can optimize the execution time in an orthogonal manner. One particular hardware-aware optimization is to layout the memory of decision trees in such a way, that higher probably paths are less likely to be evicted from system caches. This works particularly well when splits within tree nodes are uneven and have a high probability to visit one of the child nodes. In this paper, we present a method to reduce path lengths by rewarding uneven probability distributions during the training of decision trees at the cost of a minimal accuracy degradation. Specifically, we regularize the impurity computation of the CART algorithm in order to favor not only low impurity, but also highly asymmetric distributions for the evaluation of split criteria and hence offer a high optimization potential for a memory architecture-aware implementation. We show that especially for binary classification data sets and data sets with many samples, this form of regularization can lead to an reduction of up to approximately four times in the execution time with a minimal accuracy degradation. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Towards Learning Contrast Kinetics with Multi-Condition Latent Diffusion Models

Authors: Richard Osuala, Daniel M. Lang, Preeti Verma, Smriti Joshi, Apostolia Tsirikoglou, Grzegorz Skorupko, Kaisar Kushibar, Lidia Garrucho, Walter H. L. Pinaya, Oliver Diaz, Julia A. Schnabel, Karim Lekadir

Abstract: Contrast agents in dynamic contrast enhanced magnetic resonance imaging allow to localize tumors and observe their contrast kinetics, which is essential for cancer characterization and respective treatment decision-making. However, contrast agent administration is not only associated with adverse health risks, but also restricted for patients during pregnancy, and for those with kidney malfunction… ▽ More Contrast agents in dynamic contrast enhanced magnetic resonance imaging allow to localize tumors and observe their contrast kinetics, which is essential for cancer characterization and respective treatment decision-making. However, contrast agent administration is not only associated with adverse health risks, but also restricted for patients during pregnancy, and for those with kidney malfunction, or other adverse reactions. With contrast uptake as key biomarker for lesion malignancy, cancer recurrence risk, and treatment response, it becomes pivotal to reduce the dependency on intravenous contrast agent administration. To this end, we propose a multi-conditional latent diffusion model capable of acquisition time-conditioned image synthesis of DCE-MRI temporal sequences. To evaluate medical image synthesis, we additionally propose and validate the Fréchet radiomics distance as an image quality measure based on biomarker variability between synthetic and real imaging data. Our results demonstrate our method's ability to generate realistic multi-sequence fat-saturated breast DCE-MRI and uncover the emerging potential of deep learning based contrast kinetics simulation. We publicly share our accessible codebase at https://github.com/RichardObi/ccnet and provide a user-friendly library for Fréchet radiomics distance calculation at https://pypi.org/project/frd-score. △ Less

Submitted 17 July, 2024; v1 submitted 20 March, 2024; originally announced March 2024.

Comments: Early Accept at MICCAI2024

Two-step growth of high-quality single crystals of the Kitaev magnet $α$-RuCl$_{3}$

Authors: R. Namba, K. Imamura, R. Ishioka, K. Ishihara, T. Miyamoto, H. Okamoto, Y. Shimizu, Y. Saito, Y. Agarmani, M. Lang, H. Murayama, Y. Xing, S. Suetsugu, Y. Kasahara, Y. Matsuda, K. Hashimoto, T. Shibauchi

Abstract: The layered honeycomb magnet $α$-RuCl$_3$ is the most promising candidate for a Kitaev quantum spin liquid (KQSL) that can host charge-neutral Majorana fermions. Recent studies have shown significant sample dependence of thermal transport properties, which are a key probe of Majorana quasiparticles in the KQSL state, highlighting the importance of preparing high-quality single crystals of $α$-RuCl… ▽ More The layered honeycomb magnet $α$-RuCl$_3$ is the most promising candidate for a Kitaev quantum spin liquid (KQSL) that can host charge-neutral Majorana fermions. Recent studies have shown significant sample dependence of thermal transport properties, which are a key probe of Majorana quasiparticles in the KQSL state, highlighting the importance of preparing high-quality single crystals of $α$-RuCl$_3$. Here, we present a relatively simple and reliable method to grow high-quality single crystals of $α$-RuCl$_3$. We use a two-step crystal growth method consisting of a purification process by chemical vapor transport (CVT) and a main crystal growth process by sublimation. The obtained crystals exhibit a distinct first-order structural phase transition from the monoclinic ($C2/m$) to the rhombohedral ($R\bar{3}$) structure at $\sim150$ K, which is confirmed by the nuclear quadrupole resonance spectra with much sharper widths than previously reported. The Raman spectra show the absence of defect-induced modes, supporting the good crystallinity of our samples. The jumps in the thermal expansion coefficient and specific heat at the antiferromagnetic (AFM) transition at 7.6-7.7 K are larger and sharper than those of previous samples grown by the CVT and Bridgman methods and do not show any additional AFM transitions at 10-14 K due to stacking faults. The longitudinal thermal conductivity in the AFM phase is significantly larger than previously reported, indicating a very long mean free path of heat carriers. All the results indicate that our single crystals are of superior quality with good crystallinity and few stacking faults, which provides a suitable platform for studying the Kitaev physics. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 8 pages, 5 figures

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

The High Energy Light Isotope eXperiment program of direct cosmic-ray studies

Authors: HELIX Collaboration, S. Coutu, P. S. Allison, M. Baiocchi, J. J. Beatty, L. Beaufore, D. H. Calderon, A. G. Castano, Y. Chen, N. Green, D. Hanna, H. B. Jeon, S. B. Klein, B. Kunkler, M. Lang, R. Mbarek, K. McBride, S. I. Mognet, J. Musser, S. Nutter, S. OBrien, N. Park, K. M. Powledge, K. Sakai, M. Tabata , et al. (5 additional authors not shown)

Abstract: HELIX is a new NASA-sponsored instrument aimed at measuring the spectra and composition of light cosmic-ray isotopes from hydrogen to neon nuclei, in particular the clock isotopes 10Be (radioactive, with 1.4 Myr lifetime) and 9Be (stable). The latter are unique markers of the production and Galactic propagation of secondary cosmic-ray nuclei, and are needed to resolve such important mysteries as t… ▽ More HELIX is a new NASA-sponsored instrument aimed at measuring the spectra and composition of light cosmic-ray isotopes from hydrogen to neon nuclei, in particular the clock isotopes 10Be (radioactive, with 1.4 Myr lifetime) and 9Be (stable). The latter are unique markers of the production and Galactic propagation of secondary cosmic-ray nuclei, and are needed to resolve such important mysteries as the proportion of secondary positrons in the excess of antimatter observed by the AMS-02 experiment. By using a combination of a 1 T superconducting magnet spectrometer (with drift-chamber tracker) with a high-resolution time-of-flight detector system and ring-imaging Cherenkov detector, mass-resolved isotope measurements of light cosmic-ray nuclei will be possible up to 3 GeV/n in a first stratospheric balloon flight from Kiruna, Sweden to northern Canada, anticipated to take place in early summer 2024. An eventual longer Antarctic balloon flight of HELIX will yield measurements up to 10 GeV/n, sampling production from a larger volume of the Galaxy extending into the halo. We review the instrument design, testing, status and scientific prospects. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Presented at the 16th Topical Seminar on Innovative Particle and Radiation Detectors (IPRD23), Siena, Italy, to appear in JINST Proc

Towards Top-Down Reasoning: An Explainable Multi-Agent Approach for Visual Question Answering

Authors: Zeqing Wang, Wentao Wan, Qiqing Lao, Runmeng Chen, Minjie Lang, Keze Wang, Liang Lin

Abstract: Recently, several methods have been proposed to augment large Vision Language Models (VLMs) for Visual Question Answering (VQA) simplicity by incorporating external knowledge from knowledge bases or visual clues derived from question decomposition. Although having achieved promising results, these methods still suffer from the challenge that VLMs cannot inherently understand the incorporated knowl… ▽ More Recently, several methods have been proposed to augment large Vision Language Models (VLMs) for Visual Question Answering (VQA) simplicity by incorporating external knowledge from knowledge bases or visual clues derived from question decomposition. Although having achieved promising results, these methods still suffer from the challenge that VLMs cannot inherently understand the incorporated knowledge and might fail to generate the optimal answers. Contrarily, human cognition engages visual questions through a top-down reasoning process, systematically exploring relevant issues to derive a comprehensive answer. This not only facilitates an accurate answer but also provides a transparent rationale for the decision-making pathway. Motivated by this cognitive mechanism, we introduce a novel, explainable multi-agent collaboration framework designed to imitate human-like top-down reasoning by leveraging the expansive knowledge of Large Language Models (LLMs). Our framework comprises three agents, i.e., Responder, Seeker, and Integrator, each contributing uniquely to the top-down reasoning process. The VLM-based Responder generates the answer candidates for the question and gives responses to other issues. The Seeker, primarily based on LLM, identifies relevant issues related to the question to inform the Responder and constructs a Multi-View Knowledge Base (MVKB) for the given visual scene by leveraging the understanding capabilities of LLM. The Integrator agent combines information from the Seeker and the Responder to produce the final VQA answer. Through this collaboration mechanism, our framework explicitly constructs an MVKB for a specific visual scene and reasons answers in a top-down reasoning process. Extensive and comprehensive evaluations on diverse VQA datasets and VLMs demonstrate the superior applicability and interpretability of our framework over the existing compared methods. △ Less

Submitted 14 May, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: 16 pages, 5 figures

Ultrafast gap dynamics upon photodoping the Mott-insulating phase of a two-dimensional organic charge-transfer salt

Authors: Konstantin Warawa, Yassine Agarmani, Harald Schubert, Martin Dressel, Michael Lang, Hartmut G. Roskos, Mark D. Thomson

Abstract: We investigate experimentally the ultrafast changes in the spectral response of the Mott insulator $κ$-(BEDT-TTF)$_2$Cu[N(CN)$_2$]Cl ($κ$-Cl) upon photodoping with intense excitation at 1.6 eV and probing with continuum pulses simultaneously covering both the terahertz and infrared (IR) ranges (from 0 to 0.6 eV). A quantitative analysis of the differential reflectivity using a multi-band Lorentzia… ▽ More We investigate experimentally the ultrafast changes in the spectral response of the Mott insulator $κ$-(BEDT-TTF)$_2$Cu[N(CN)$_2$]Cl ($κ$-Cl) upon photodoping with intense excitation at 1.6 eV and probing with continuum pulses simultaneously covering both the terahertz and infrared (IR) ranges (from 0 to 0.6 eV). A quantitative analysis of the differential reflectivity using a multi-band Lorentzian model provide absolute changes in spectral weights and objective global time constants for the relaxation vs. temperature. The transient conductivity spectra deduced from the analysis suggest that the transient photoinduced spectral weight is dominated by a progressive closure of the Mott gap with increasing excitation density, i.e. due to changes in the inter-Hubbard-band absorption by the remaining singly occupied states. We critically examine this scenario compared to that proposed previously, whereby the low-energy spectral weight is attributed to a Drude-like response of photoexcited doublons/holons. We also consider the observed slowing down of the relaxation rate with increasing excitation density, and temperature dependence of the initial doublon/holon density in terms of the phonon-mediated gap recombination model. △ Less

Submitted 28 October, 2023; originally announced October 2023.

Slow and Non-Equilibrium Dynamics due to Electronic Ferroelectricity in a Strongly-Correlated Molecular Conductor

Authors: Tatjana Thomas, Yassine Agarmani, Steffi Hartmann, Mark Kartsovnik, Natalia Kushch, Stephen M. Winter, Sebastian Schmid, Peter Lunkenheimer, Michael Lang, Jens Mueller

Abstract: Using a combination of resistance fluctuation (noise) and dielectric spectroscopy we investigate the nature of relaxor-type electronic ferroelectricity in the organic conductor $κ$-(BETS)$_2$Mn[N(CN)$_2$]$_3$, a system representative for a wider class of materials, where strong correlations of electrons on a lattice of dimerized molecules results in an insulating ground state. The two complementar… ▽ More Using a combination of resistance fluctuation (noise) and dielectric spectroscopy we investigate the nature of relaxor-type electronic ferroelectricity in the organic conductor $κ$-(BETS)$_2$Mn[N(CN)$_2$]$_3$, a system representative for a wider class of materials, where strong correlations of electrons on a lattice of dimerized molecules results in an insulating ground state. The two complementary spectroscopies reveal a distinct low-frequency dynamics. By dielectric spectroscopy we detect an intrinsic relaxation that is typical for relaxor ferroelectrics below the metal-to-insulator transition at $T_{\rm{MI}}\sim 25\,$K. Resistance noise spectroscopy reveals fluctuating two-level processes above $T_{\rm MI}$ which strongly couple to the applied electric field, a signature of fluctuating polar nanoregions (PNR), i.e. clusters of quantum electric dipoles fluctuating collectively. The PNR preform above the metal insulator transition. Upon cooling through $T_{\rm MI}$, a drastic increase of the low-frequency $1/f$-type fluctuations and slowing down of the charge carrier dynamics is accompanied by the onset of strong non-equilibrium dynamics indicating a glassy transition of interacting dipolar clusters, the scaling properties of which are consistent with a droplet model. The freezing of nano-scale polar clusters and non-equilibrium dynamics is suggested to be a common feature of organic relaxor-type electronic ferroelectrics and needs to be considered in theoretical models describing these materials. △ Less

Submitted 3 November, 2023; v1 submitted 26 October, 2023; originally announced October 2023.

Combined experimental and theoretical studies on glasslike transitions in the frustrated molecular conductors $θ$-(BEDT-TTF)$_2MM'$(SCN)$_4$

Authors: Yohei Saito, Owen Ganter, Chao Shang, Kenichiro Hashimoto, Takahiko Sasaki, Stephen M. Winter, Jens Müller, Michael Lang

Abstract: We present results of the coefficient of thermal expansion for the frustrated quasi-two-dimensional molecular conductor $θ$-(BEDT-TTF)$_2$RbZn(SCN)$_4$ for temperatures 1.5 K $\leq T \leq$ 290 K. A pronounced first-order phase transition anomaly is observed at the combined charge-order/structural transition at 215 K. Furthermore, clear evidence is found for two separate glasslike transitions at… ▽ More We present results of the coefficient of thermal expansion for the frustrated quasi-two-dimensional molecular conductor $θ$-(BEDT-TTF)$_2$RbZn(SCN)$_4$ for temperatures 1.5 K $\leq T \leq$ 290 K. A pronounced first-order phase transition anomaly is observed at the combined charge-order/structural transition at 215 K. Furthermore, clear evidence is found for two separate glasslike transitions at $T_{\mathrm{g}}$ = 90-100 K and $T_{\mathrm{g}}^\dagger$ = 120-130 K, similar to previous findings for $θ$-(BEDT-TTF)$_2$CsZn(SCN)$_4$ and $θ$-(BEDT-TTF)$_2$CsCo(SCN)$_4$, reported in T. Thomas et al., Phys. Rev. B 105, L041114 (2022), both of which lack the charge-order/structural transition. Our findings indicate that these glasslike transitions are common features for the $θ$-(BEDT-TTF)$_2MM^\prime$(SCN)$_4$ family with $M$ = (Rb, Cs) and $M^\prime$ = (Co, Zn), irrespective of the presence or absence of charge order. These results are consistent with our model calculations on the glasslike dynamics associated with the flexible ethylene endgroups of the BEDT-TTF molecules for various $θ$-(BEDT-TTF)$_2MM^\prime$(SCN)$_4$ salts, predicting two different conformational glass transitions. Moreover, calculations of the hopping integrals show a substantial degree of dependence on the endgroups' conformation, suggesting a significant coupling to the electronic degrees of freedom. Our findings support the possibility that the glassy freezing of the ethylene endgroups could drive or enhance glassy charge dynamics. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 11 pages, 10 figures, 1 table

Journal ref: Phys. Rev. Research 6, 023003 (2024)

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.

Less Power for More Learning: Restricting OCaml Features for Effective Teaching

Authors: Max Lang, Nico Petzendorfer

Abstract: We present a framework for sandboxing and restricting features of the OCaml programming language to effectively automate the grading of programming exercises, scaling to hundreds of submissions. We describe how to disable language and library features that should not be used to solve a given exercise. We present an overview of an implementation of a mock IO system to allow testing of IO-related ex… ▽ More We present a framework for sandboxing and restricting features of the OCaml programming language to effectively automate the grading of programming exercises, scaling to hundreds of submissions. We describe how to disable language and library features that should not be used to solve a given exercise. We present an overview of an implementation of a mock IO system to allow testing of IO-related exercises in a controlled environment. Finally, we detail a number of security considerations to ensure submitted code remains sandboxed, allowing automatic grading to be trusted without manual verification. The source code of our implementation is publicly available. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: for associated git repository, see https://github.com/just-max/less-power

ACM Class: K.3.2

Supervised dimensionality reduction for multiple imputation by chained equations

Authors: Edoardo Costantini, Kyle M. Lang, Klaas Sijtsma

Abstract: Multivariate imputation by chained equations (MICE) is one of the most popular approaches to address missing values in a data set. This approach requires specifying a univariate imputation model for every variable under imputation. The specification of which predictors should be included in these univariate imputation models can be a daunting task. Principal component analysis (PCA) can simplify t… ▽ More Multivariate imputation by chained equations (MICE) is one of the most popular approaches to address missing values in a data set. This approach requires specifying a univariate imputation model for every variable under imputation. The specification of which predictors should be included in these univariate imputation models can be a daunting task. Principal component analysis (PCA) can simplify this process by replacing all of the potential imputation model predictors with a few components summarizing their variance. In this article, we extend the use of PCA with MICE to include a supervised aspect whereby information from the variables under imputation is incorporated into the principal component estimation. We conducted an extensive simulation study to assess the statistical properties of MICE with different versions of supervised dimensionality reduction and we compared them with the use of classical unsupervised PCA as a simpler dimensionality reduction technique. △ Less

Submitted 31 October, 2023; v1 submitted 4 September, 2023; originally announced September 2023.

Tailoring Stateless Model Checking for Event-Driven Multi-Threaded Programs

Authors: Parosh Aziz Abdulla, Mohamed Faouzi Atig, Frederik Meyer Bønneland, Sarbojit Das, Bengt Jonsson, Magnus Lång, Konstantinos Sagonas

Abstract: Event-driven multi-threaded programming is an important idiom for structuring concurrent computations. Stateless Model Checking (SMC) is an effective verification technique for multi-threaded programs, especially when coupled with Dynamic Partial Order Reduction (DPOR). Existing SMC techniques are often ineffective in handling event-driven programs, since they will typically explore all possible o… ▽ More Event-driven multi-threaded programming is an important idiom for structuring concurrent computations. Stateless Model Checking (SMC) is an effective verification technique for multi-threaded programs, especially when coupled with Dynamic Partial Order Reduction (DPOR). Existing SMC techniques are often ineffective in handling event-driven programs, since they will typically explore all possible orderings of event processing, even when events do not conflict. We present Event-DPOR , a DPOR algorithm tailored to event-driven multi-threaded programs. It is based on Optimal-DPOR, an optimal DPOR algorithm for multi-threaded programs; we show how it can be extended for event-driven programs. We prove correctness of Event-DPOR for all programs, and optimality for a large subclass. One complication is that an operation in Event-DPOR, which checks for redundancy of new executions, is NP-hard, as we show in this paper; we address this by a sequence of inexpensive (but incomplete) tests which check for redundancy efficiently. Our implementation and experimental evaluation show that, in comparison with other tools in which handler threads are simulated using locks, Event-DPOR can be exponentially faster than other state-of-the-art DPOR algorithms on a variety of programs and manages to completely avoid unnecessary exploration of executions. △ Less

Submitted 29 July, 2023; originally announced July 2023.

Comments: Published at ATVA23

ACM Class: D.3

Controlling stable Bloch points with electric currents

Authors: Martin Lang, Swapneel Amit Pathak, Samuel J. R. Holt, Marijan Beg, Hans Fangohr

Abstract: The Bloch point is a point singularity in the magnetisation configuration, where the magnetisation vanishes. It can exist as an equilibrium configuration and plays an important role in many magnetisation reversal processes. In the present work, we focus on manipulating Bloch points in a system that can host stable Bloch points - a two-layer FeGe nanostrip with opposite chirality of the two layers.… ▽ More The Bloch point is a point singularity in the magnetisation configuration, where the magnetisation vanishes. It can exist as an equilibrium configuration and plays an important role in many magnetisation reversal processes. In the present work, we focus on manipulating Bloch points in a system that can host stable Bloch points - a two-layer FeGe nanostrip with opposite chirality of the two layers. We drive Bloch points using spin-transfer torques and find that Bloch points can move collectively without any Hall effect and report that Bloch points are repelled from the sample boundaries and each other. We study pinning of Bloch points at wedge-shaped constrictions (notches) in the nanostrip and demonstrate that arrays of Bloch points can be moved past a series of notches in a controlled manner by applying consecutive current pulses of different strength. Finally, we simulate a T-shaped geometry and demonstrate that a Bloch point can be moved along different paths by applying current between suitable strip ends. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 12 pages, 7 figures

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

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)

Current Status of the Muon g-2 Interpretations within Two-Higgs-Doublet Models

Authors: Syuhei Iguro, Teppei Kitahara, Martin S. Lang, Michihisa Takeuchi

Abstract: In this article, we review and update implications of the muon anomalous magnetic moment (muon $g-2$) anomaly for two-Higgs-doublet models (2HDMs), which are classified according to imposed symmetries and their resulting Yukawa sector. In the minimal setup, the muon $g-2$ anomaly can be accommodated by the type-X (lepto-philic) 2HDM, flavor-aligned 2HDM (FA2HDM), muon-specific 2HDM ($μ$2HDM), and… ▽ More In this article, we review and update implications of the muon anomalous magnetic moment (muon $g-2$) anomaly for two-Higgs-doublet models (2HDMs), which are classified according to imposed symmetries and their resulting Yukawa sector. In the minimal setup, the muon $g-2$ anomaly can be accommodated by the type-X (lepto-philic) 2HDM, flavor-aligned 2HDM (FA2HDM), muon-specific 2HDM ($μ$2HDM), and $μτ$-flavor violating 2HDM. We summarize all relevant experimental constraints from high-energy collider experiments and flavor experiments, as well as the theoretical constraints from the perturbative unitarity and vacuum stability bounds, to these 2HDMs in light of the muon $g-2$ anomaly. We clarify the available parameter spaces of these 2HDMs and investigate how to probe the remaining parameter regions in future experiments. In particular, we find that, due to the updated $B_s\toμ^+ μ^-$ measurement, the remaining parameter region of the FA2HDM is almost equivalent to the one of the type-X 2HDM. Furthermore, based on collider simulations, we find that the type-X 2HDM is excluded and the $μ$2HDM scenario will be covered with the upcoming Run 3 data. △ Less

Submitted 14 November, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: 46 pages, 9 figures, 2 tables

Report number: P3H--23--025, TTP23--013

Journal ref: Phys.Rev.D 108 (2023) 11, 115012

Energetics and Dynamics of a stable Bloch point

Authors: Thomas Brian Winkler, Marijan Beg, Martin Lang, Mathias Kläui, Hans Fangohr

Abstract: Magnetic Bloch points (BPs) are highly confined magnetization configurations, that often occur in transient spin dynamics processes. However, opposing chiralities of adjacent layers for instance in a FeGe bilayer stack can stabilize such magnetic BPs at the layer interface. These BPs configurations are metastable and consist of two coupled vortices (one in each layer) with same circularity and opp… ▽ More Magnetic Bloch points (BPs) are highly confined magnetization configurations, that often occur in transient spin dynamics processes. However, opposing chiralities of adjacent layers for instance in a FeGe bilayer stack can stabilize such magnetic BPs at the layer interface. These BPs configurations are metastable and consist of two coupled vortices (one in each layer) with same circularity and opposite polarity. Each vortex is stabilized by opposite sign Dzyaloshinskii-Moriya interactions. This stabilization mechanism potentially opens the door towards BP-based spintronic applications. An open question, from a methodological point of view, is whether the Heisenberg (HB) model approach (atomistic model) as to be used to study such systems or if the -- computationally more efficient -- micromagnetic (MM) models can be used and still obtain robust results. We are modelling and comparing the energetics and dynamics of a stable BP obtained using both HB and MM approaches. We find that an MM description of a stable BP leads qualitatively to the same results as the HB description, and that an appropriate mesh discretization plays a more important role than the chosen model. Further, we study the dynamics by shifting the BP with an applied in-plane field and investigating the relaxation after switching the filed off abruptly. The precessional motion of coupled vortices in a BP state can be drastically reduced compared to a classical vortex, which may be also an interesting feature for fast and efficient devices. A recent study has shown that a bilayer stack hosting BPs can be used to retain information [1]. △ Less

Submitted 17 March, 2023; originally announced March 2023.

Comments: 19 pages, 6 figures

Pressure study on the interplay between magnetic order and valence-change crossover in EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$

Authors: Bernd Wolf, Theresa Lundbeck, Jan Zimmermann, Marius Peters, Kristin Kliemt, Cornelius Krellner, Michael Lang

Abstract: We present results of the magnetic susceptibility on high-quality single crystals of EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$ for Ge concentrations 0 $\leq x \leq$ 0.105 performed under varying hydrostatic (He-gas) pressure 0 $\leq p \leq$ 0.5 GPa. The work extends on recent studies at ambient pressure demonstrating the drastic change in the magnetic response from valence-change-crossover behavior for $x$ =… ▽ More We present results of the magnetic susceptibility on high-quality single crystals of EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$ for Ge concentrations 0 $\leq x \leq$ 0.105 performed under varying hydrostatic (He-gas) pressure 0 $\leq p \leq$ 0.5 GPa. The work extends on recent studies at ambient pressure demonstrating the drastic change in the magnetic response from valence-change-crossover behavior for $x$ = 0 and 0.058, to long-range antiferromagnetic (afm) order below $T_{\text{N}}$ = 47 K for $x$ = 0.105. The valence-change-crossover temperature $T'_{\text{V}}$ shows an extraordinarily strong pressure dependence of d$T'_{\text{V}}$/d$p$ = +(80 $\pm$ 10) K/GPa. In contrast, a very small pressure dependence of d$T_{\text{N}}$/d$p \leq$ +(1 $\pm$ 0.5) K/GPa is found for the afm order upon pressurizing the $x$ = 0.105 crystal from $p$ = 0 to 0.05 GPa. Remarkably, by further increasing the pressure to 0.1 GPa, a drastic change in the ground state from afm order to valence-change-crossover behavior is observed. Estimates of the electronic entropy, derived from analyzing susceptibility data at varying pressures, indicate that the boundary between afm order and valence-change crossover represents a first-order phase transition. Our results suggest a particular type of second-order critical endpoint of the first-order transition for $x$ = 0.105 at $p_{\text{cr}} \approx$ 0.06 GPa and $T_{\text{cr}} \approx$ 45 K where intriguing strong-coupling effects between fluctuating charge-, spin- and lattice degrees of freedom can be expected. △ Less

Submitted 2 June, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

From valence fluctuations to long-range magnetic order in EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$ single crystals

Authors: Marius Peters, Kristin Kliemt, Michelle Ocker, Bernd Wolf, Pascal Puphal, Matthieu Le Tacon, Michael Merz, Michael Lang, Cornelius Krellner

Abstract: EuPd$_2$Si$_2$ is a valence-fluctuating system undergoing a temperature-induced valence crossover at $T'_V\approx160\,$K. We present the successful single crystal growth using the Czochralski method for the substitution series EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$, with substitution levels $x\leq 0.15$. A careful determination of the germanium content revealed that only half of the nominal concentration… ▽ More EuPd$_2$Si$_2$ is a valence-fluctuating system undergoing a temperature-induced valence crossover at $T'_V\approx160\,$K. We present the successful single crystal growth using the Czochralski method for the substitution series EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$, with substitution levels $x\leq 0.15$. A careful determination of the germanium content revealed that only half of the nominal concentration is build into the crystal structure. From thermodynamic measurements it is established that $T'_V$ is strongly suppressed for small substitution levels and antiferromagnetic order from stable divalent europium emerges for $x\gtrsim 0.10$. The valence transition is accompanied by a pronounced change of the lattice parameter $a$ of order 1.8%. In the antiferromagnetically ordered state below $T_N = 47$ K, we find sizeable magnetic anisotropy with an easy plane perpendicular to the crystallographic c direction. An entropy analysis revealed that no valence fluctuations are present for the magnetically ordered materials. Combining the obtained thermodynamic and structural data, we construct a concentration-temperature phase diagram demonstrating a rather abrupt change from a valence-fluctuating to a magnetically-ordered state in EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$. △ Less

Submitted 24 May, 2023; v1 submitted 13 March, 2023; originally announced March 2023.

Journal ref: Phys. Rev. Materials 7, 064405 (2023)

3D Masked Autoencoders with Application to Anomaly Detection in Non-Contrast Enhanced Breast MRI

Authors: Daniel M. Lang, Eli Schwartz, Cosmin I. Bercea, Raja Giryes, Julia A. Schnabel

Abstract: Self-supervised models allow (pre-)training on unlabeled data and therefore have the potential to overcome the need for large annotated cohorts. One leading self-supervised model is the masked autoencoder (MAE) which was developed on natural imaging data. The MAE is masking out a high fraction of visual transformer (ViT) input patches, to then recover the uncorrupted images as a pretraining task.… ▽ More Self-supervised models allow (pre-)training on unlabeled data and therefore have the potential to overcome the need for large annotated cohorts. One leading self-supervised model is the masked autoencoder (MAE) which was developed on natural imaging data. The MAE is masking out a high fraction of visual transformer (ViT) input patches, to then recover the uncorrupted images as a pretraining task. In this work, we extend MAE to perform anomaly detection on breast magnetic resonance imaging (MRI). This new model, coined masked autoencoder for medical imaging (MAEMI) is trained on two non-contrast enhanced MRI sequences, aiming at lesion detection without the need for intravenous injection of contrast media and temporal image acquisition. During training, only non-cancerous images are presented to the model, with the purpose of localizing anomalous tumor regions during test time. We use a public dataset for model development. Performance of the architecture is evaluated in reference to subtraction images created from dynamic contrast enhanced (DCE)-MRI. △ Less

Submitted 10 March, 2023; originally announced March 2023.

Numerical simulation projects in micromagnetics with Jupyter

Authors: Martin Lonsky, Martin Lang, Samuel Holt, Swapneel Amit Pathak, Robin Klause, Tzu-Hsiang Lo, Marijan Beg, Axel Hoffmann, Hans Fangohr

Abstract: We report a case study where an existing materials science course was modified to include numerical simulation projects on the micromagnetic behavior of materials. The Ubermag micromagnetic simulation software package is used in order to solve problems computationally. The simulation software is controlled through Python code in Jupyter notebooks. Our experience is that the self-paced problem-solv… ▽ More We report a case study where an existing materials science course was modified to include numerical simulation projects on the micromagnetic behavior of materials. The Ubermag micromagnetic simulation software package is used in order to solve problems computationally. The simulation software is controlled through Python code in Jupyter notebooks. Our experience is that the self-paced problem-solving nature of the project work can facilitate a better in-depth exploration of the course contents. We discuss which aspects of the Ubermag and the project Jupyter ecosystem have been beneficial for the students' learning experience and which could be transferred to similar teaching activities in other subject areas. △ Less

Submitted 17 July, 2024; v1 submitted 3 March, 2023; originally announced March 2023.

Comments: Main article: 9 pages; supplementary material: 24 pages. Accepted for publication in American Journal of Physics

Field-induced effects in the spin liquid candidate PbCuTe$_{2}$O$_{6}$

Authors: Paul Eibisch, Christian Thurn, Arif Ata, Ulrich Tutsch, Yohei Saito, Steffi Hartmann, Bernd Wolf, Abanoub R. N. Hanna, A. T. M. Nazmul Islam, Shravani Chillal, Bella Lake, Michael Lang

Abstract: PbCuTe$_2$O$_6$ is considered as one of the rare candidate materials for a three-dimensional quantum spin liquid (QSL). This assessment was based on the results of various magnetic experiments, performed mainly on polycrystalline material. More recent measurements on single crystals revealed an even more exotic behavior, yielding ferroelectric order below $T_{\text{FE}}\approx 1\,\text{K}$, accomp… ▽ More PbCuTe$_2$O$_6$ is considered as one of the rare candidate materials for a three-dimensional quantum spin liquid (QSL). This assessment was based on the results of various magnetic experiments, performed mainly on polycrystalline material. More recent measurements on single crystals revealed an even more exotic behavior, yielding ferroelectric order below $T_{\text{FE}}\approx 1\,\text{K}$, accompanied by distinct lattice distortions, and a somewhat modified magnetic response which is still consistent with a QSL. Here we report on low-temperature measurements of various thermodynamic, magnetic and dielectric properties of single crystalline PbCuTe$_2$O$_6$ in magnetic fields $B\leq 14.5\,\text{T}$. The combination of these various probes allows us to construct a detailed $B$-$T$ phase diagram including a ferroelectric phase for $B \leq$ $8\,\text{T}$ and a $B$-induced magnetic phase at $B \geq$ $11\,\text{T}$. These phases are preceded by or coincide with a structural transition from a cubic high-temperature phase into a distorted non-cubic low-temperature state. The phase diagram discloses two quantum critical points (QCPs) in the accessible field range, a ferroelectric QCP at $B_{c1}$ = $7.9\,\text{T}$ and a magnetic QCP at $B_{c2}$ = $11\,\text{T}$. Field-induced lattice distortions, observed in the state at $T>$ $1\,\text{K}$ and which are assigned to the effect of spin-orbit interaction of the Cu$^{2+}$-ions, are considered as the key mechanism by which the magnetic field couples to the dielectric degrees of freedom in this material. △ Less

Submitted 9 May, 2023; v1 submitted 27 February, 2023; originally announced February 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

A simple and general approach for reversible condensation polymerization with cyclization

Authors: Michael Lang, Kiran Suresh Kumar

Abstract: We develop a simple recursive approach to treat reversible condensation polymerization with cyclization. Based upon a minimum set of balance equations, the law of mass action, Gaussian chain statistics, and the assumption of independent reactions, we derive exact analytical solutions for systems without cyclization, for systems containing only smallest loops, or systems that exclusively form loops… ▽ More We develop a simple recursive approach to treat reversible condensation polymerization with cyclization. Based upon a minimum set of balance equations, the law of mass action, Gaussian chain statistics, and the assumption of independent reactions, we derive exact analytical solutions for systems without cyclization, for systems containing only smallest loops, or systems that exclusively form loops. Exact numerical solutions are computed for the general case of a homopolymerization of flexible precursor polymers. All solutions were tested with Monte-Carlo simulations. A generalization for good solvent is discussed and it is shown that this generalization agrees with preceding work in the limit of low and high polymer volume fractions. The new aspect of our approach is its flexibility that allows for a rather simple generalization to more complex situations. These include different kinds of reversible linear polymerization, non-linear polymerization, first shell substitution effect, semiflexibility, or a low molecular weight cut-off for cyclization. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Journal ref: Macromolecules 54, 2021, 7021

Reversible stepwise condensation polymerization with cyclization: strictly alternating co-polymerization and homopolymerization based upon two orthogonal reactions

Authors: Michael Lang, Kiran Suresh Kumar

Abstract: In a preceding work [M. Lang, K. Kumar, A simple and general approach for reversible condensation polymerization with cyclization, Macromolecules 54 (2021) 7021], we have introduced a simple recursive scheme that allows to treat stepwise linear reversible polymerizations of any kind with cyclization. This approach is used to discuss the polymerization of linear Gaussian strands (LGS) with two diff… ▽ More In a preceding work [M. Lang, K. Kumar, A simple and general approach for reversible condensation polymerization with cyclization, Macromolecules 54 (2021) 7021], we have introduced a simple recursive scheme that allows to treat stepwise linear reversible polymerizations of any kind with cyclization. This approach is used to discuss the polymerization of linear Gaussian strands (LGS) with two different reactive groups $A$ and $B$ on either chain end that participate in two orthogonal reactions and the strictly alternating copolymerization of LGS that carry $A$ reactive groups with LGS equipped with type $B$ reactive groups. The former of these cases has not been discussed theoretically in literature, the latter only regarding some special cases. We provide either analytical expressions or exact numerical solutions for the general cases with and without cyclization. Weight distributions, averages, polydispersity, and the weight fractions of cyclic and linear species are computed. All numerical solutions were tested by Monte-Carlo simulations. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Journal ref: Macromolecules 54, 2021, 7036

Causal Analysis of Influence of the Solar Cycle and Latitudinal Solar-Wind Structure on Corotation Forecasts

Authors: Nachiketa Chakraborty, Harriet Turner, Mathew Owens, Matthew Lang

Abstract: Studying solar wind conditions is central to forecasting impact of space weather on Earth. Under the assumption that the structure of this wind is constant in time and corotates with the Sun, solar wind and thereby space weather forecasts have been made quite effectively. Such corotation forecasts are well studied with decades of observations from STEREO and near-Earth spacecrafts. Forecast accura… ▽ More Studying solar wind conditions is central to forecasting impact of space weather on Earth. Under the assumption that the structure of this wind is constant in time and corotates with the Sun, solar wind and thereby space weather forecasts have been made quite effectively. Such corotation forecasts are well studied with decades of observations from STEREO and near-Earth spacecrafts. Forecast accuracy depends upon the latitudinal separation (or offset $Δθ$) between source and spacecraft, forecast lead time ($Δt$) and the solar cycle via the sunspot number (SSN). The precise dependencies factoring in uncertain- ties however, are a mixture of influences from each of these factors. And for high precision forecasts, it is important to understand what drives the forecast accuracy and its uncertainty. Here we present a causal inference approach based on information theoretic measures to do this. Our framework can compute not only the direct (linear and non-linear) dependencies of the forecast mean absolute error (MAE) on SSN, $Δt$ and $Δt$, but also how these individual variables combine to enhance or diminish the MAE. We provide an initial assessment of this with potential of aiding data assimilation in the future. △ Less

Submitted 27 January, 2023; originally announced January 2023.

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

The new APD-Based Readout of the Crystal Barrel Calorimeter -- An Overview

Authors: CBELSA/TAPS Collaboration, :, C. Honisch, P. Klassen, J. Müllers, M. Urban, F. Afzal, J. Bieling, S. Ciupka, J. Hartmann, P. Hoffmeister, M. Lang, D. Schaab, C. Schmidt, M. Steinacher, D. Walther, R. Beck, K. -T. Brinkmann, V. Crede, H. Dutz, D. Elsner, W. Erni, E. Fix, F. Frommberger, M. Grüner , et al. (26 additional authors not shown)

Abstract: The Crystal Barrel is an electromagnetic calorimeter consisting of 1380 CsI(Tl) scintillators, and is currently installed at the CBELSA/TAPS experiment where it is used to detect decay products from photoproduction of mesons. The readout of the Crystal Barrel has been upgraded in order to integrate the detector into the first level of the trigger and to increase its sensitivity for neutral final s… ▽ More The Crystal Barrel is an electromagnetic calorimeter consisting of 1380 CsI(Tl) scintillators, and is currently installed at the CBELSA/TAPS experiment where it is used to detect decay products from photoproduction of mesons. The readout of the Crystal Barrel has been upgraded in order to integrate the detector into the first level of the trigger and to increase its sensitivity for neutral final states. The new readout uses avalanche photodiodes in the front-end and a dual back-end with branches optimized for energy and time measurement, respectively. An FPGA-based cluster finder processes the whole hit pattern within less than 100 ns. The important downside of APDs -- the temperature dependence of their gain -- is handled with a temperature stabilization and a compensating bias voltage supply. Additionally, a light pulser system allows the APDs' gains to be measured during beamtimes. △ Less

Submitted 16 January, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

$B_s \to μ^+ μ^-$ in a Two-Higgs-Doublet Model with flavour-changing up-type Yukawa couplings

Authors: Martin S. Lang, Ulrich Nierste

Abstract: We present a Two-Higgs-Doublet Model in which the structure of the quark Yukawa matrices is governed by three spurions breaking the flavour symmetries of the quark Yukawa sector. The model naturally suppresses flavour-changing neutral current (FCNC) amplitudes in the down-type sector, but permits sizable FCNC couplings in the up sector. We calculate the branching ratio of $B_s \to μ^+ μ^-$ to lead… ▽ More We present a Two-Higgs-Doublet Model in which the structure of the quark Yukawa matrices is governed by three spurions breaking the flavour symmetries of the quark Yukawa sector. The model naturally suppresses flavour-changing neutral current (FCNC) amplitudes in the down-type sector, but permits sizable FCNC couplings in the up sector. We calculate the branching ratio of $B_s \to μ^+ μ^-$ to leading and next-to-leading order of QCD for the case with FCNC couplings of the heavy neutral Higgs bosons to up-type quarks and verify that all counterterms follow the pattern dictated by the spurion expansion of the Yukawa matrices. We find correlations between $B_s \to μ^+ μ^-$, $b\to sγ$, and the Higgs masses. The $B_s - \bar B_s$ mixing amplitude is naturally suppressed in the model but can probe a portion of the parameter space with very heavy Higgs bosons. △ Less

Submitted 12 April, 2024; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: 35 pages, 11 figures

Report number: TTP22-074, P3H-22-128

LoTSS Jellyfish Galaxies IV: Enhanced Star Formation on the Leading Half of Cluster Galaxies and Gas Compression in IC3949

Authors: Ian D. Roberts, Maojin Lang, Daria Trotsenko, Ashley Bemis, Sara L. Ellison, Lihwai Lin, Hsi-An Pan, Alessandro Ignesti, Sarah Leslie, Reinout J. van Weeren

Abstract: With MaNGA integral field spectroscopy, we present a resolved analysis of star formation for 29 jellyfish galaxies in nearby clusters, identified from radio continuum imaging taken by the Low Frequency Array. Simulations predict enhanced star formation on the "leading half" of galaxies undergoing ram pressure stripping, and in this work we report observational evidence for this elevated star forma… ▽ More With MaNGA integral field spectroscopy, we present a resolved analysis of star formation for 29 jellyfish galaxies in nearby clusters, identified from radio continuum imaging taken by the Low Frequency Array. Simulations predict enhanced star formation on the "leading half" of galaxies undergoing ram pressure stripping, and in this work we report observational evidence for this elevated star formation. The dividing line (through the galaxy center) that maximizes this star formation enhancement is systematically tied to the observed direction of the ram pressure stripped tail, suggesting a physical connection between ram pressure and this star formation enhancement. We also present a case study on the distribution of molecular gas in one jellyfish galaxy from our sample, IC3949, using ALMA CO J=1-0, HCN J=1-0, and HCO$^+$ J=1-0 observations from the ALMaQUEST survey. The $\mathrm{H_2}$ depletion time (as traced by CO) in IC3949 ranges from $\sim\!1\,\mathrm{Gyr}$ in the outskirts of the molecular gas disk to $\sim\!11\,\mathrm{Gyr}$ near the galaxy center. IC3949 shows a clear region of enhanced star formation on the leading half of the galaxy where the average depletion time is $\sim\!2.7\,\mathrm{Gyr}$, in line with the median value for the galaxy on the whole. Dense gas tracers, HCN and HCO$^+$, are only detected at the galaxy center and on the leading half of IC3949. Our results favour a scenario in which ram pressure compresses the interstellar medium, promoting the formation of molecular gas that in turn fuels a localized increase of star formation. △ Less

Submitted 28 October, 2022; originally announced October 2022.

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

From magnetic order to valence-change crossover in EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$ using He-gas pressure

Authors: Bernd Wolf, Felix Spathelf, Jan Zimmermann, Theresa Lundbeck, Kristin Kliemt, Cornelius Krellner, Michael Lang

Abstract: We present results of magnetic susceptibility and thermal expansion measurements performed on high-quality single crystals of EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$ for 0 $\leq$ x $\leq$ 0.2 and temperatures 2 K $\leq T \leq$ 300 K. Data were taken at ambient pressure and finite He-gas pressure $p$ $\leq$ 0.5 GPa. For x = 0 and ambient pressure we observe a pronounced valence-change crossover centred arou… ▽ More We present results of magnetic susceptibility and thermal expansion measurements performed on high-quality single crystals of EuPd$_2$(Si$_{1-x}$Ge$_x$)$_2$ for 0 $\leq$ x $\leq$ 0.2 and temperatures 2 K $\leq T \leq$ 300 K. Data were taken at ambient pressure and finite He-gas pressure $p$ $\leq$ 0.5 GPa. For x = 0 and ambient pressure we observe a pronounced valence-change crossover centred around $T'_V$ $\approx$ 160 K with a non-magnetic ground state. This valence-change crossover is characterized by an extraordinarily strong pressure dependence of d$T'_V$ /d$p$ = (80 $\pm 10)$ K/GPa. We observe a shift of $T'_V$ to lower temperatures with increasing Ge-concentration, reaching $T'_V$ $\approx$ 90 K for x = 0.1, while still showing a non-magnetic ground state. Remarkably, on further increasing x to 0.2 we find a stable Eu$^{(2+δ)+}$ valence with long-range antiferromagnetic order below $T_N$ = (47.5 $\pm$ 0.1) K, reflecting a close competition between two energy scales in this system. In fact, by the application of hydrostatic pressure as small as 0.1 GPa, the ground state of this system can be changed from long-range antiferromagnetic order for $p$ $<$ 0.1 GPa to an intermediate-valence state for $p$ $\geq$ 0.1 GPa. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: 9 pages, 3 figures

SIR-HUXt -- a particle filter data assimilation scheme for assimilating CME time-elongation profiles

Authors: Luke Barnard, Mathew Owens, Chris Scott, Matthew Lang, Mike Lockwood

Abstract: We present the development of SIR-HUXt, the integration of a sequential importance resampling (SIR) data assimilation scheme with the HUXt solar wind model. SIR-HUXt is designed to assimilate the time-elongation profiles of CME fronts in the low heliosphere, such as those typically extracted from heliospheric imager data returned by the STEREO, Parker Solar Probe, and Solar Orbiter missions. We us… ▽ More We present the development of SIR-HUXt, the integration of a sequential importance resampling (SIR) data assimilation scheme with the HUXt solar wind model. SIR-HUXt is designed to assimilate the time-elongation profiles of CME fronts in the low heliosphere, such as those typically extracted from heliospheric imager data returned by the STEREO, Parker Solar Probe, and Solar Orbiter missions. We use Observing System Simulation Experiments to explore the performance of SIR-HUXt for a simple synthetic CME scenario of a fully Earth directed CME flowing through a uniform ambient solar wind, where the CME is initialised with the average observed CME speed and width. These experiments are performed for a range of observer locations, from 20 deg to 90 deg behind Earth, spanning the L5 point where ESA's future Vigil space weather monitor will return heliospheric imager data for operational space weather forecasting. We show that SIR-HUXt performs well at constraining the CME speed, and has some success at constraining the CME longitude. The CME width is largely unconstrained by the SIR-HUXt assimilations, and more experiments are required to determine if this is due to this specific CME scenario, or is a general feature of assimilating time-elongation profiles. Rank-histograms suggest that the SIR-HUXt ensembles are well calibrated, with no clear indications of bias or under/over dispersion. Improved constraints on the initial CME speed lead directly to improvements in the CME transit time to Earth and arrival speed. For an observer in the L5 region, SIR-HUXt returned a 69% reduction in the CME transit time uncertainty, and a 63% reduction in the arrival speed uncertainty. This suggests SIR-HUXt has potential to improve the real-world representivity of HUXt simulations, and therefore has potential to reduce the uncertainty of CME arrival time hindcasts and forecasts. △ Less

Submitted 5 October, 2022; originally announced October 2022.

Comments: 29 pages, 10 figures, 1 table

Optimal independent generating system for the congruence subgroups $Γ_0(p)$ and $Γ_0(p^2)$

Authors: Nhat Minh Doan, Sang-hyun Kim, Mong Lung Lang, Ser Peow Tan

Abstract: Let $n$ be a prime or its square. We prove that the congruence subgroup $Γ_0(n)$ admits a free product decomposition into cyclic factors in such a way that the $(2,1)$-component of each cyclic generator is either $n$ or $0$, answering a conjecture of Kulkarni. We can also require that the Frobenius norm of each generator is less than $2n-1$. A crucial observation is that if $P$ denotes the convex… ▽ More Let $n$ be a prime or its square. We prove that the congruence subgroup $Γ_0(n)$ admits a free product decomposition into cyclic factors in such a way that the $(2,1)$-component of each cyclic generator is either $n$ or $0$, answering a conjecture of Kulkarni. We can also require that the Frobenius norm of each generator is less than $2n-1$. A crucial observation is that if $P$ denotes the convex hull of the extended Farey sequence of order $\lfloor \sqrt{n} \rfloor$ in the hyperbolic plane $\mathbb{H}^2$, then the projection $π: \mathbb{H}^2\to \mathbb{H}^2/Γ_0(n)$ is injective on the interior of $P$ and each connected component of $π(\mathbb{H}^2)\setminusπ(P)$ is either an order-three cone of area $π/3$ or an ideal triangle. Denoting by $m(Γ_0(n))$ the minimum of the largest denominator in the cusp set of $Q$ where $Q$ ranges over all possible special (fundamental) polygons for $Γ_0(n)$, we establish the inequality $ \lfloor \sqrt{n} \rfloor \le m(Γ_0(n))\le \lfloor \sqrt{4n/3} \rfloor$, and completely characterize the cases in which the bounds are achieved. We also prove analogous results when $n$ is the multiplication of two sufficiently close odd primes. △ Less

Submitted 9 February, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: With some minor revisions, the multiplications of two nearby primes are also considered

MSC Class: 11F06 (Primary); 11B57; 30F35 (Secondary)

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

High-dimensional Imputation for the Social Sciences: a Comparison of State-of-the-art Methods

Authors: Edoardo Costantini, Kyle M. Lang, Tim Reeskens, Klaas Sijtsma

Abstract: Including a large number of predictors in the imputation model underlying a multiple imputation (MI) procedure is one of the most challenging tasks imputers face. A variety of high-dimensional MI techniques can help, but there has been limited research on their relative performance. In this study, we investigated a wide range of extant high-dimensional MI techniques that can handle a large number… ▽ More Including a large number of predictors in the imputation model underlying a multiple imputation (MI) procedure is one of the most challenging tasks imputers face. A variety of high-dimensional MI techniques can help, but there has been limited research on their relative performance. In this study, we investigated a wide range of extant high-dimensional MI techniques that can handle a large number of predictors in the imputation models and general missing data patterns. We assessed the relative performance of seven high-dimensional MI methods with a Monte Carlo simulation study and a resampling study based on real survey data. The performance of the methods was defined by the degree to which they facilitate unbiased and confidencevalid estimates of the parameters of complete data analysis models. We found that using lasso penalty or forward selection to select the predictors used in the MI model and using principal component analysis to reduce the dimensionality of auxiliary data produce the best results. △ Less

Submitted 14 August, 2023; v1 submitted 29 August, 2022; originally announced August 2022.

Awaiting for Godot: Stateless Model Checking that Avoids Executions where Nothing Happens

Authors: Bengt Jonsson, Magnus Lång, Konstantinos Sagonas

Abstract: Stateless Model Checking (SMC) is a verification technique for concurrent programs that checks for safety violations by exploring all possible thread schedulings. It is highly effective when coupled with Dynamic Partial Order Reduction (DPOR), which introduces an equivalence on schedulings and need explore only one in each equivalence class. Even with DPOR, SMC often spends unnecessary effort in e… ▽ More Stateless Model Checking (SMC) is a verification technique for concurrent programs that checks for safety violations by exploring all possible thread schedulings. It is highly effective when coupled with Dynamic Partial Order Reduction (DPOR), which introduces an equivalence on schedulings and need explore only one in each equivalence class. Even with DPOR, SMC often spends unnecessary effort in exploring loop iterations that are pure, i.e., have no effect on the program state. We present techniques for making SMC with DPOR more effective on programs with pure loop iterations. The first is a static program analysis to detect loop purity and an associated program transformation, called Partial Loop Purity Elimination, that inserts assume statements to block pure loop iterations. Subsequently, some of these assumes are turned into await statements that completely remove many assume-blocked executions. Finally, we present an extension of the standard DPOR equivalence, obtained by weakening the conflict relation between events. All these techniques are incorporated into a new DPOR algorithm, Optimal-DPOR-Await, which can handle both awaits and the weaker conflict relation, is optimal in the sense that it explores exactly one execution in each equivalence class, and can also diagnose livelocks. Our implementation in Nidhugg shows that these techniques can significantly speed up the analysis of concurrent programs that are currently challenging for SMC tools, both for exploring their complete set of interleavings, but even for detecting concurrency errors in them. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Characterization of electroless nickel-phosphorus plating for ultracold-neutron storage

Authors: H. Akatsuka, T. Andalib, B. Bell, J. Berean-Dutcher, N. Bernier, C. P. Bidinosti, C. Cude-Woods, S. A. Currie, C. A. Davis, B. Franke, R. Gaur, P. Giampa, S. Hansen-Romu, M. T. Hassan, K. Hatanaka, T. Higuchi, C. Gibson, G. Ichikawa, I. Ide, S. Imajo, T. M. Ito, B. Jamieson, S. Kawasaki, M. Kitaguchi, W. Klassen , et al. (29 additional authors not shown)

Abstract: Electroless nickel plating is an established industrial process that provides a robust and relatively low-cost coating suitable for transporting and storing ultracold neutrons (UCN). Using roughness measurements and UCN-storage experiments we characterized UCN guides made from polished aluminum or stainless-steel tubes plated by several vendors. All electroless nickel platings were similarly suite… ▽ More Electroless nickel plating is an established industrial process that provides a robust and relatively low-cost coating suitable for transporting and storing ultracold neutrons (UCN). Using roughness measurements and UCN-storage experiments we characterized UCN guides made from polished aluminum or stainless-steel tubes plated by several vendors. All electroless nickel platings were similarly suited for UCN storage with an average loss probability per wall bounce of $2.8\cdot10^{-4}$ to $4.1\cdot10^{-4}$ for energies between 90 neV and 190 neV, or a ratio of imaginary to real Fermi potential $η$ of $1.7\cdot10^{-4}$ to $3.3\cdot10^{-4}$. Measurements at different elevations indicate that the energy dependence of UCN losses is well described by the imaginary Fermi potential. Some special considerations are required to avoid an increase in surface roughness during the plating process and hence a reduction in UCN transmission. Increased roughness had only a minor impact on storage properties. Based on these findings we chose a vendor to plate the UCN-production vessel that will contain the superfluid-helium converter for the new TRIUMF UltraCold Advanced Neutron (TUCAN) source, achieving acceptable UCN-storage properties with ${η=3.5(5)\cdot10^{-4}}$. △ Less

Submitted 7 February, 2023; v1 submitted 10 August, 2022; originally announced August 2022.

Comments: 8 pages, 7 figures

Journal ref: Nucl. Instrum. Methods Phys. Res. A 1049 (2023) 168106