Search

We analyse continuum and molecular emission, observed with ALMA, from the dust-enshrouded intermediate-mass AGB star OH 30.1 -0.7. We find a secondary peak in the continuum maps, "feature B", separated by 4.6" from the AGB star, which corresponds to a projected separation of $1.8 \times 10^4$ au, placing a lower limit on the physical separation. This feature is most likely composed of cold dust and is likely to be ejecta associated with the AGB star, though we cannot rule out that it is a background object. The molecular emission we detect includes lines of CO, SiS, CS, SO$_2$, NS, NaCl, and KCl. We find that the NS emission is off centre and arranged along an axis perpendicular to the direction of feature B, indicative of a UV-emitting binary companion (e.g. a G-type main sequence star or hotter), perhaps on an eccentric orbit, contributing to its formation. However, the NaCl and KCl emission constrain the nature of that companion to not be hotter than a late B-type main sequence star. We find relatively warm emission arising from the inner wind and detect several vibrationally excited lines of SiS (3 = 1), NaCl (up to 3 = 4) and KCl (up to 3 = 2), and emission from low energy levels in the mid to outer envelope, as traced by SO$_2$. The CO emission is abruptly truncated around 3.5" or 14,000 au from the continuum peak, suggesting that mass loss at a high rate may have commenced as little as 2800 years ago., Comment: Accepted for publication in MNRAS

We study a flaring activity of the HSP Mrk421 that was characterized from radio to very-high-energy (VHE; E $>0.1$TeV) gamma rays with MAGIC, Fermi-LAT, Swift, XMM-Newton and several optical and radio telescopes. These observations included, for the first time for a gamma-ray flare of a blazar, simultaneous X-ray polarization measurements with IXPE. We find substantial variability in both X-rays and VHE gamma rays throughout the campaign, with the highest VHE flux above 0.2 TeV occurring during the IXPE observing window, and exceeding twice the flux of the Crab Nebula. However, the VHE and X-ray spectra are on average softer, and the correlation between these two bands weaker that those reported in previous flares of Mrk421. IXPE reveals an X-ray polarization degree significantly higher than that at radio and optical frequencies. The X-ray polarization angle varies by $\sim$100$^\circ$ on timescales of days, and the polarization degree changes by more than a factor 4. The highest X-ray polarization degree reaches 26%, around which a X-ray counter-clockwise hysteresis loop is measured with XMM-Newton. It suggests that the X-ray emission comes from particles close to the high-energy cutoff, hence possibly probing an extreme case of the Turbulent Extreme Multi-Zone model. We model the broadband emission with a simplified stratified jet model throughout the flare. The polarization measurements imply an electron distribution in the X-ray emitting region with a very high minimum Lorentz factor, which is expected in electron-ion plasma, as well as a variation of the emitting region size up to a factor of three during the flaring activity. We find no correlation between the fluxes and the evolution of the model parameters, which indicates a stochastic nature of the underlying physical mechanism. Such behaviour would be expected in a highly turbulent electron-ion plasma crossing a shock front., Comment: Submitted to Astronomy and Astrophysics. Corresponding authors: Axel Arbet-Engels, Lea Heckmann, David Paneque

OT 081 is a well-known, luminous blazar that is remarkably variable in many energy bands. We present the first broadband study of the source which includes very-high-energy (VHE, $E>$100\,GeV) $\gamma$-ray data taken by the MAGIC and H.E.S.S. imaging Cherenkov telescopes. The discovery of VHE $\gamma$-ray emission happened during a high state of $\gamma$-ray activity in July 2016, observed by many instruments from radio to VHE $\gamma$-rays. We identify four states of activity of the source, one of which includes VHE $\gamma$-ray emission. Variability in the VHE domain is found on daily timescales. The intrinsic VHE spectrum can be described by a power-law with index $3.27\pm0.44_{\rm stat}\pm0.15_{\rm sys}$ (MAGIC) and $3.39\pm0.58_{\rm stat}\pm0.64_{\rm sys}$ (H.E.S.S.) in the energy range of 55--300\,GeV and 120--500\,GeV, respectively. The broadband emission cannot be sucessfully reproduced by a simple one-zone synchrotron self-Compton model. Instead, an additional external Compton component is required. We test a lepto-hadronic model that reproduces the dataset well and a proton-synchrotron dominated model that requires an extreme proton luminosity. Emission models that are able to successfully represent the data place the emitting region well outside of the Broad Line Region (BLR) to a location at which the radiative environment is dominated by the infrared thermal radiation field of the dusty torus. In the scenario described by this flaring activity, the source appears to be an FSRQ, in contrast with past categorizations. This suggests that the source can be considered to be a transitional blazar, intermediate between BL~Lac and FSRQ objects., Comment: Accepted on MNRAS Corresponding authors: M. Manganaro, J. Becerra Gonz\'alez, M. Seglar-Arroyo, D. A. Sanchez

Imaging atmospheric Cherenkov telescopes (IACTs) are used to observe very high-energy photons from the ground. Gamma rays are indirectly detected through the Cherenkov light emitted by the air showers they induce. The new generation of experiments, in particular the Cherenkov Telescope Array Observatory (CTAO), sets ambitious goals for discoveries of new gamma-ray sources and precise measurements of the already discovered ones. To achieve these goals, both hardware and data analysis must employ cutting-edge techniques. This also applies to the LST-1, the first IACT built for the CTAO, which is currently taking data on the Canary island of La Palma. This paper introduces a new event reconstruction technique for IACT data, aiming to improve the image reconstruction quality and the discrimination between the signal and the background from misidentified hadrons and electrons. The technique models the development of the extensive air shower signal, recorded as a waveform per pixel, seen by CTAO telescopes' cameras. Model parameters are subsequently passed to random forest regressors and classifiers to extract information on the primary particle. The new reconstruction was applied to simulated data and to data from observations of the Crab Nebula performed by the LST-1. The event reconstruction method presented here shows promising performance improvements. The angular and energy resolution, and the sensitivity, are improved by 10 to 20% over most of the energy range. At low energy, improvements reach up to 22%, 47%, and 50%, respectively. A future extension of the method to stereoscopic analysis for telescope arrays will be the next important step., Comment: Accepted in A&A

Sparse autoencoders (SAEs) are a promising technique for decomposing language model activations into interpretable linear features. However, current SAEs fall short of completely explaining model performance, resulting in "dark matter": unexplained variance in activations. This work investigates dark matter as an object of study in its own right. Surprisingly, we find that much of SAE dark matter--about half of the error vector itself and >90% of its norm--can be linearly predicted from the initial activation vector. Additionally, we find that the scaling behavior of SAE error norms at a per token level is remarkably predictable: larger SAEs mostly struggle to reconstruct the same contexts as smaller SAEs. We build on the linear representation hypothesis to propose models of activations that might lead to these observations, including postulating a new type of "introduced error"; these insights imply that the part of the SAE error vector that cannot be linearly predicted ("nonlinear" error) might be fundamentally different from the linearly predictable component. To validate this hypothesis, we empirically analyze nonlinear SAE error and show that 1) it contains fewer not yet learned features, 2) SAEs trained on it are quantitatively worse, 3) it helps predict SAE per-token scaling behavior, and 4) it is responsible for a proportional amount of the downstream increase in cross entropy loss when SAE activations are inserted into the model. Finally, we examine two methods to reduce nonlinear SAE error at a fixed sparsity: inference time gradient pursuit, which leads to a very slight decrease in nonlinear error, and linear transformations from earlier layer SAE outputs, which leads to a larger reduction., Comment: Code at https://github.com/JoshEngels/SAE-Dark-Matter

Sparse autoencoders have recently produced dictionaries of high-dimensional vectors corresponding to the universe of concepts represented by large language models. We find that this concept universe has interesting structure at three levels: 1) The "atomic" small-scale structure contains "crystals" whose faces are parallelograms or trapezoids, generalizing well-known examples such as (man-woman-king-queen). We find that the quality of such parallelograms and associated function vectors improves greatly when projecting out global distractor directions such as word length, which is efficiently done with linear discriminant analysis. 2) The "brain" intermediate-scale structure has significant spatial modularity; for example, math and code features form a "lobe" akin to functional lobes seen in neural fMRI images. We quantify the spatial locality of these lobes with multiple metrics and find that clusters of co-occurring features, at coarse enough scale, also cluster together spatially far more than one would expect if feature geometry were random. 3) The "galaxy" scale large-scale structure of the feature point cloud is not isotropic, but instead has a power law of eigenvalues with steepest slope in middle layers. We also quantify how the clustering entropy depends on the layer., Comment: 13 pages, 12 figures

Sparse autoencoders (SAEs) are a recent technique for decomposing neural network activations into human-interpretable features. However, in order for SAEs to identify all features represented in frontier models, it will be necessary to scale them up to very high width, posing a computational challenge. In this work, we introduce Switch Sparse Autoencoders, a novel SAE architecture aimed at reducing the compute cost of training SAEs. Inspired by sparse mixture of experts models, Switch SAEs route activation vectors between smaller "expert" SAEs, enabling SAEs to efficiently scale to many more features. We present experiments comparing Switch SAEs with other SAE architectures, and find that Switch SAEs deliver a substantial Pareto improvement in the reconstruction vs. sparsity frontier for a given fixed training compute budget. We also study the geometry of features across experts, analyze features duplicated across experts, and verify that Switch SAE features are as interpretable as features found by other SAE architectures., Comment: Code available at https://github.com/amudide/switch_sae

The Flat Spectrum Radio Quasar OP~313 entered an enhanced activity phase in November 2023 and has undergone multiple flares since then which have motivated the organisation of several large multi-wavelength campaigns, including two deep observations from the hard X-ray telescope NuSTAR. The broadband emission from OP~313 during these two observations is investigated under a new unified analysis framework, with data spanning from optical to gamma rays. Traditional methods for analyzing blazar emissions often rely on proprietary software tailored to specific instruments, making it challenging to integrate and interpret data from multi-wavelength campaigns comprehensively. This study demonstrates the feasibility of utilizing gammapy, an open-source Python package, and common data formats originally developed for gamma-ray instrumentation, to perform a consistent multi-instrument analysis. This enables a forward folding approach that fully incorporates source observations, detector responses, and various instrumental and astrophysical backgrounds. The methodology is applied as an example to recent data collected from the distant quasar OP~313. We present a comprehensive data reconstruction and analysis for instruments including the Liverpool Telescope's IO:O detector, Swift-UVOT, Swift-XRT, NuSTAR, and Fermi-LAT. The resulting spectral analysis is validated against the native tools for each instrument. Additionally, a multi-wavelength phenomenological model of the source emission, encompassing optical to gamma-ray bands, is developed, incorporating absorption components across different energy regimes. We introduce and validate a new unified framework for multi-wavelength forward folding data analysis based on gammapy and open data formats, demonstrating its application to spectral data from the quasar OP~313., Comment: 20 pages, 10 figures github: https://github.com/mireianievas/gammapy_mwl_workflow zenodo: https://zenodo.org/records/13837637

Instruments for gamma-ray astronomy at Very High Energies ($E>100\,{\rm GeV}$) have traditionally derived their scientific results through proprietary data and software. Data standardisation has become a prominent issue in this field both as a requirement for the dissemination of data from the next generation of gamma-ray observatories and as an effective solution to realise public data legacies of current-generation instruments. Specifications for a standardised gamma-ray data format have been proposed as a community effort and have already been successfully adopted by several instruments. We present the first production of standardised data from the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes. We converted $166\,{\rm h}$ of observations from different sources and validated their analysis with the open-source software Gammapy. We consider six data sets representing different scientific and technical analysis cases and compare the results obtained analysing the standardised data with open-source software against those produced with the MAGIC proprietary data and software. Aiming at a systematic production of MAGIC data in this standardised format, we also present the implementation of a database-driven pipeline automatically performing the MAGIC data reduction from the calibrated down to the standardised data level. In all the cases selected for the validation, we obtain results compatible with the MAGIC proprietary software, both for the manual and for the automatic data productions. Part of the validation data set is also made publicly available, thus representing the first large public release of MAGIC data. This effort and this first data release represent a technical milestone toward the realisation of a public MAGIC data legacy., Comment: Accepted for publication in the Journal of High Energy Astrophysics

Since the early days of quantum mechanics hydrogen, as the simplest of all atoms, has been studied or used to investigate new physics. In parallel, this knowledge leads to different applications, e.g. a spin filter to separate metastable hydrogen atoms in single hyperfine substates with electron spin m_s=1/2. Subsequently, this work provides the necessary theory as well as experimental conditions to build a new generation of spin filter which permits the separation of all four individual metastable hydrogen hyperfine states as well as for its isotopes in a corresponding beam.

Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$\sigma$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$\sigma$), as well as the so-called bridge emission (5.7$\sigma$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources., Comment: Accepted by A&A

The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This energy dependence, although suppressed by the massive QG energy scale $E_\mathrm{QG}$, expected to be on the level of the Planck energy $1.22 \times 10^{19}$ GeV, is potentially detectable in astrophysical observations. In this scenario, the cosmological distances traversed by photons act as an amplifier for this effect. By leveraging the observation of a remarkable flare from the blazar Mrk\,421, recorded at energies above 100 GeV by the MAGIC telescopes on the night of April 25 to 26, 2014, we look for time delays scaling linearly and quadratically with the photon energies. Using for the first time in LIV studies a binned-likelihood approach we set constraints on the QG energy scale. For the linear scenario, we set $95\%$ lower limits $E_\mathrm{QG}>2.7\times10^{17}$ GeV for the subluminal case and $E_\mathrm{QG}> 3.6 \times10^{17}$ GeV for the superluminal case. For the quadratic scenario, the $95\%$ lower limits for the subluminal and superluminal cases are $E_\mathrm{QG}>2.6 \times10^{10}$ GeV and $E_\mathrm{QG}>2.5\times10^{10}$ GeV, respectively.

Railroad transportation plays a vital role in the future of sustainable mobility. Besides building new infrastructure, capacity can be improved by modern train control systems, e.g., based on moving blocks. At the same time, there is only limited work on how to optimally route trains using the potential gained by these systems. Recently, an initial approach for train routing with moving block control has been proposed to address this demand. However, detailed evaluations on so-called lazy constraints are missing, and no publicly available implementation exists. In this work, we close this gap by providing an extended approach as well as a flexible open-source implementation that can use different solving strategies. Using that, we experimentally evaluate what choices should be made when implementing a lazy constraint approach. The corresponding implementation and benchmarks are publicly available as part of the Munich Train Control Toolkit (MTCT) at https://github.com/cda-tum/mtct., Comment: Accepted Version: 19th Conference on Computer Science and Intelligence Systems (FedCSIS 2024)

The main science aim of the BlackGEM array is to detect optical counterparts to gravitational wave mergers. Additionally, the array will perform a set of synoptic surveys to detect Local Universe transients and short time-scale variability in stars and binaries, as well as a six-filter all-sky survey down to ~22nd mag. The BlackGEM Phase-I array consists of three optical wide-field unit telescopes. Each unit uses an f/5.5 modified Dall-Kirkham (Harmer-Wynne) design with a triplet corrector lens, and a 65cm primary mirror, coupled with a 110Mpix CCD detector, that provides an instantaneous field-of-view of 2.7~square degrees, sampled at 0.564\arcsec/pixel. The total field-of-view for the array is 8.2 square degrees. Each telescope is equipped with a six-slot filter wheel containing an optimised Sloan set (BG-u, BG-g, BG-r, BG-i, BG-z) and a wider-band 440-720 nm (BG-q) filter. Each unit telescope is independent from the others. Cloud-based data processing is done in real time, and includes a transient-detection routine as well as a full-source optimal-photometry module. BlackGEM has been installed at the ESO La Silla observatory as of October 2019. After a prolonged COVID-19 hiatus, science operations started on April 1, 2023 and will run for five years. Aside from its core scientific program, BlackGEM will give rise to a multitude of additional science cases in multi-colour time-domain astronomy, to the benefit of a variety of topics in astrophysics, such as infant supernovae, luminous red novae, asteroseismology of post-main-sequence objects, (ultracompact) binary stars, and the relation between gravitational wave counterparts and other classes of transients, Comment: Published in PASP

Recent work has proposed that language models perform computation by manipulating one-dimensional representations of concepts ("features") in activation space. In contrast, we explore whether some language model representations may be inherently multi-dimensional. We begin by developing a rigorous definition of irreducible multi-dimensional features based on whether they can be decomposed into either independent or non-co-occurring lower-dimensional features. Motivated by these definitions, we design a scalable method that uses sparse autoencoders to automatically find multi-dimensional features in GPT-2 and Mistral 7B. These auto-discovered features include strikingly interpretable examples, e.g. circular features representing days of the week and months of the year. We identify tasks where these exact circles are used to solve computational problems involving modular arithmetic in days of the week and months of the year. Next, we provide evidence that these circular features are indeed the fundamental unit of computation in these tasks with intervention experiments on Mistral 7B and Llama 3 8B. Finally, we find further circular representations by breaking down the hidden states for these tasks into interpretable components, and we examine the continuity of the days of the week feature in Mistral 7B., Comment: Code and data at https://github.com/JoshEngels/MultiDimensionalFeatures

Automated volumetric segmentation of the pancreas on cross-sectional imaging is needed for diagnosis and follow-up of pancreatic diseases. While CT-based pancreatic segmentation is more established, MRI-based segmentation methods are understudied, largely due to a lack of publicly available datasets, benchmarking research efforts, and domain-specific deep learning methods. In this retrospective study, we collected a large dataset (767 scans from 499 participants) of T1-weighted (T1W) and T2-weighted (T2W) abdominal MRI series from five centers between March 2004 and November 2022. We also collected CT scans of 1,350 patients from publicly available sources for benchmarking purposes. We developed a new pancreas segmentation method, called PanSegNet, combining the strengths of nnUNet and a Transformer network with a new linear attention module enabling volumetric computation. We tested PanSegNet's accuracy in cross-modality (a total of 2,117 scans) and cross-center settings with Dice and Hausdorff distance (HD95) evaluation metrics. We used Cohen's kappa statistics for intra and inter-rater agreement evaluation and paired t-tests for volume and Dice comparisons, respectively. For segmentation accuracy, we achieved Dice coefficients of 88.3% (std: 7.2%, at case level) with CT, 85.0% (std: 7.9%) with T1W MRI, and 86.3% (std: 6.4%) with T2W MRI. There was a high correlation for pancreas volume prediction with R^2 of 0.91, 0.84, and 0.85 for CT, T1W, and T2W, respectively. We found moderate inter-observer (0.624 and 0.638 for T1W and T2W MRI, respectively) and high intra-observer agreement scores. All MRI data is made available at https://osf.io/kysnj/. Our source code is available at https://github.com/NUBagciLab/PaNSegNet., Comment: Peer-reviewer version

The momentum space Josephson effect describes the supercurrent flow between weakly coupled Bose-Einstein condensates (BECs) at two discrete momentum states. Here, we experimentally observe this exotic phenomenon using a BEC with Raman-induced spin-orbit coupling, where the tunneling between two local band minima is implemented by the momentum kick of an additional optical lattice. A sudden quench of the Raman detuning induces coherent spin-momentum oscillations of the BEC, which is analogous to the a.c. Josephson effect. We observe both plasma and regular Josephson oscillations in different parameter regimes. The experimental results agree well with the theoretical model and numerical simulation, and showcase the important role of nonlinear interactions. We also show that the measurement of the Josephson plasma frequency gives the Bogoliubov zero quasimomentum gap, which determines the mass of the corresponding pseudo-Goldstone mode, a long-sought phenomenon in particle physics. The observation of momentum space Josephson physics offers an exciting platform for quantum simulation and sensing utilizing momentum states as a synthetic degree., Comment: 12 pages, 11 figures

We developed an open-source scalar wave transport model to estimate the generalized scattering matrix (S matrix) of a disordered medium in the diffusion regime. Here, the term generalization refers to the incorporation of evanescent wave field modes in addition to propagating modes while estimating the S matrix. For that we used the scalar Kirchhoff-Helmholtz boundary integral formulation together with the Green's function perturbation method to generalize the conventional Fisher-Lee relations to include evanescent modes as well. The estimated S matrix, which satisfies generalized unitarity and reciprocity conditions, is modeled for a 2D disordered waveguide. The generalized transmission matrix contained in the S matrix is used to estimate the optimal phase-conjugate wavefront for focusing onto an evanescent mode. The phenomena of universal transmission value of 2/3 for such an optimal phase conjugate wavefront is also shown in the context of evanescent wave mode focusing through a diffusive disorder. The presented code framework may be of interest to wavefront shaping researchers for visualizing and estimating wave transport properties in general.

We developed a 2D Finite-Difference Time-Domain (FDTD) method for modeling a space-time modulated guidestar targeting wavefront shaping applications in disordered media. Space-time modulation in general (a particular example being the acousto-optic effect) is used here as a guidestar for the transverse confinement of light around the tagged region surrounded by disorder. Together with the guidestar, the iterative optical phase conjugation (IOPC) method is used to overcome the diffusion of light due to multiple scattering. A phase sensitive lock-in detection technique is utilized to estimate the steady-state amplitude and phase of the modulated wavefronts emerging from the guidestar region continuously operating in the Raman-Nath regime. As the IOPC scheme naturally converges to the maximally transmitting eigenchannel profile, one could use the position of the guidestar within the disorder to channelize the maximal transmission through the tagged region. The associated code developed in MATLAB is provided as an open source (MIT license) package. The code package is referred by the acronym STAR-FDTD where STAR stands for Space-Time modulated Acousto-optic guidestaR., Comment: Accepted in Journal of Physics : Photonics, Focus Issue on Foundational Skills and Tools for Building Wavefront Shaping Systems, Original submission on Nov 30, 2023 and revised submission on Jul 31, 2024

PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

Within the 'Medicina/Effelsberg H2O maser monitoring program' we observed U Her and RR Aql at 22-GHz for about two decades between 1990 and 2011, with a gap between 1997 and 2000 in the case of RR Aql. In addition, maps were obtained in the period 1990-1992 of U Her with the Very Large Array. We find that the strongest emission in U Her is located in a shell with boundaries 11-25 AU. The gas crossing time is 8.5 years. We derive lifetimes for individual maser clouds of less than 4 years, based on the absence of detectable line-of-sight velocity drifts of the maser emission. The shell is not evenly filled, and its structure is maintained on timescales much longer than those of individual maser clouds. Both stars show brightness variability on several timescales. The prevalent variation is periodic, following the optical variability of the stars with a lag of 2-3 months. Superposed are irregular fluctuations, of a few months' duration, of increased or decreased excitation at particular locations, and long-term systematic variations on timescales of a decade or more. The properties of the maser emission are governed by those of the stellar wind while traversing the water maser shell. Inhomogeneities in the wind affecting the excitation conditions and prevalent beaming directions likely cause the variations seen on timescales longer than the stellar pulsation period. We propose the existence of long-living regions in the shells, which maintain favourable excitation conditions on timescales of the wind crossing times through the shells or orbital periods of (sub-)stellar companions., Comment: Main paper: 26 pages, 5 tables, 22 figures. Appendices: 27 pages (with figures of all spectra, and tables). Accepted for publication in Astronomy & Astrophysics

Coloring a graph is a well known problem and used in many different contexts. Here we want to assign $k \geq 1$ colors to each vertex of a graph $G$ such that each edge has two different colors at each endpoint. Such a vertex-coloring, if exists, is called a feasible coloring of $G$. \textsc{Distance Coloring} is an extension to the standard \textsc{Coloring} problem. Here we want to enforce that every pair of distinct vertices of distance less than or equal to $d$ have different colors, for integers $d \geq 1$ and $k \geq d+1$. Reconfiguration problems ask if two given configurations can be transformed into each other with certain rules. For example, the well-known \textsc{Coloring Reconfiguration} asks if there is a way to change one vertex's color at a time, starting from a feasible given coloring $\alpha$ of a graph $G$ to reach another feasible given coloring $\beta$ of $G$, such that all intermediate colorings are also feasible. In this paper, we study the reconfiguration of distance colorings on certain graph classes. We show that even for planar, bipartite, and $2$-degenerate graphs, reconfiguring distance colorings is $\mathsf{PSPACE}$-complete for $d \geq 2$ and $k = \Omega(d^2)$ via a reduction from the well-known \textsc{Sliding Tokens} problem. Additionally, we show that the problem on split graphs remains $\mathsf{PSPACE}$-complete when $d = 2$ and large $k$ but can be solved in polynomial time when $d \geq 3$ and $k \geq d+1$, and design a quadratic-time algorithm to solve the problem on paths for any $d \geq 2$ and $k \geq d+1$., Comment: 28 pages, 8 figures, v2: minor revision of v1

In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC's standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory's northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths., Comment: 18 pages, 13 figures, submitted to MNRAS

We present our past and current long-term monitoring program of water masers in the circumstellar envelopes of evolved stars, augmented by occasional interferometric observations. Using as example the Mira-variable U Her, we identify three types of variability: periodic (following the optical variation), long-term (years-decades) and short-term irregular (weeks-months). We show there are regions in the maser shell where excitation conditions are favourable, which remain stable for many years. Lifetimes of maser clouds in the wind-acceleration zone are of the order of up to a few years. Much longer lifetimes are found for the peculiar case of a maser cloud outside that zone (as in RT Vir), or in some cases where the motion of spectral features can be followed for the entire 2 decade monitoring period (as in red supergiant VX Sgr)., Comment: 4 pages, 8 figures. To appear in the proceedings of IAU Symposium 380 (on "Cosmic Masers")

We define and investigate the problem of $\textit{c-approximate window search}$: approximate nearest neighbor search where each point in the dataset has a numeric label, and the goal is to find nearest neighbors to queries within arbitrary label ranges. Many semantic search problems, such as image and document search with timestamp filters, or product search with cost filters, are natural examples of this problem. We propose and theoretically analyze a modular tree-based framework for transforming an index that solves the traditional c-approximate nearest neighbor problem into a data structure that solves window search. On standard nearest neighbor benchmark datasets equipped with random label values, adversarially constructed embeddings, and image search embeddings with real timestamps, we obtain up to a $75\times$ speedup over existing solutions at the same level of recall., Comment: Code available: https://github.com/JoshEngels/RangeFilteredANN

Detecting cosmic gamma rays at high rates is the key to time-resolve the acceleration of particles within some of the most powerful events in the universe. Time-resolving the emission of gamma rays from merging celestial bodies, apparently random bursts of gamma rays, recurring novas in binary systems, flaring jets from active galactic nuclei, clocking pulsars, and many more became a critical contribution to astronomy. For good timing on account of high rates, we would ideally collect the naturally more abundant, low energetic gamma rays in the domain of one giga electronvolt in large areas. Satellites detect low energetic gamma rays but only in small collecting areas. Cherenkov telescopes have large collecting areas but can only detect the rare, high energetic gamma rays. To detect gamma rays with lower energies, Cherenkov-telescopes need to increase in precision and size. But when we push the concept of the -- far/tele -- seeing Cherenkov telescope accordingly, the telescope's physical limits show more clearly. The narrower depth-of-field of larger mirrors, the aberrations of mirrors, and the deformations of mirrors and mechanics all blur the telescope's image. To overcome these limits, we propose to record the -- full/plenum -- Cherenkov-light field of an atmospheric shower, i.e. recording the directions and impacts of each individual Cherenkov photon simultaneously, with a novel class of instrument. This novel Cherenkov plenoscope can turn a narrow depth-of-field into the perception of depth, can compensate aberrations, and can tolerate deformations. We design a Cherenkov plenoscope to explore timing by detecting low energetic gamma rays in large areas., Comment: An interdisciplinary effort investigated by physicists and civil engineers. For simulations, see https://github.com/cherenkov-plenoscope. 28 pages, 25 figures, 50 references, 100% motivation. Log:(2024 February 09) Add journal reference, update email and affiliations

Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the radio galaxy in the center of the cluster: NGC 1275. By modeling the magnetic field surrounding this target, we searched for spectral indications of ALP presence. Despite finding no statistical evidence of ALP signatures, we were able to exclude ALP models in the sub-micro electronvolt range. Our analysis improved upon previous work by calculating the full likelihood and statistical coverage for all considered models across the parameter space. Consequently, we achieved the most stringent limits to date for ALP masses around 50 neV, with cross sections down to $g_{a\gamma} = 3 \times 10^{-12}$ GeV$^{-1}$., Comment: 25 pages, 10 figures, accepted for publication in Physics of the Dark Universe

We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured 0.2-1 TeV flux is close to the quiescent state and ranges from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the VHE and X-ray emission are positively correlated at a $4\sigma$ significance level. The IXPE measurements unveil a X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring unveils an X-ray flux increase with a clear spectral hardening. It suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counter-clockwise), implying important changes in the particle acceleration efficiency on $\sim$hour timescales., Comment: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

This paper studies density-based clustering of point sets. These methods use dense regions of points to detect clusters of arbitrary shapes. In particular, we study variants of density peaks clustering, a popular type of algorithm that has been shown to work well in practice. Our goal is to cluster large high-dimensional datasets, which are prevalent in practice. Prior solutions are either sequential, and cannot scale to large data, or are specialized for low-dimensional data. This paper unifies the different variants of density peaks clustering into a single framework, PECANN, by abstracting out several key steps common to this class of algorithms. One such key step is to find nearest neighbors that satisfy a predicate function, and one of the main contributions of this paper is an efficient way to do this predicate search using graph-based approximate nearest neighbor search (ANNS). To provide ample parallelism, we propose a doubling search technique that enables points to find an approximate nearest neighbor satisfying the predicate in a small number of rounds. Our technique can be applied to many existing graph-based ANNS algorithms, which can all be plugged into PECANN. We implement five clustering algorithms with PECANN and evaluate them on synthetic and real-world datasets with up to 1.28 million points and up to 1024 dimensions on a 30-core machine with two-way hyper-threading. Compared to the state-of-the-art FASTDP algorithm for high-dimensional density peaks clustering, which is sequential, our best algorithm is 45x-734x faster while achieving competitive ARI scores. Compared to the state-of-the-art parallel DPC-based algorithm, which is optimized for low dimensions, we show that PECANN is two orders of magnitude faster. As far as we know, our work is the first to evaluate DPC variants on large high-dimensional real-world image and text embedding datasets.

Postdoctoral researchers are important contributors to European higher education, playing significant roles in research, but also in teaching and societal engagement. Despite their valuable contribution, the situation of postdoctoral researchers remains challenging in Europe as issues such as temporary or short-term contracts, high workloads and financial insecurity impact their work and mental health. In this context, many European universities are stepping up the support they provide at institutional level to shape highly skilled and well-rounded postdoctoral researchers, ready to explore a broad variety of job opportunities both inside and beyond academia. This report outlines the conclusions of the EUA-CDE Thematic Peer Group on "Profiles and institutional support for postdoctoral researchers", established in 2023 to provide an overview of the current services and resources available at the postdoctoral level. With findings gathered from 28 universities in 16 European countries, this publication reflects the group's discussions concerning the current state of play and the most pressing issues faced by postdoctoral researchers. This report addresses existing institutional support in place for the postdoctoral community, discusses new strategies, provides examples of good practices and suggests possible ways forward for European universities on how to improve the postdoctoral experience. It also serves as a source of reflection for institutions considering creating a more sustainable environment for this population and fostering their postdoctoral services and support.

Sizable hyperpolarisation, i.e. an imbalance of the occupation numbers of nuclear spins in a sample deviating from thermal equilibrium, is needed in various fields of science. For example, hyperpolarised tracers are utilised in magnetic resonance imaging in medicine (MRI) and polarised beams and targets are employed in nuclear physics to study the spin dependence of nuclear forces. Here we show that the quantum interference of transitions induced by radio-wave pumping with longitudinal and radial pulses are able to produce large polarisations at small magnetic fields. This method is easier than established methods, theoretically understood and experimentally proven for beams of metastable hydrogen atoms in the keV energy range. It should also work for a variety of samples at rest. Thus, this technique opens the door for a new generation of polarised tracers, possibly low-field MRI with better spatial resolution or the production of polarised fuel to increase the efficiency of fusion reactors by manipulating the involved cross sections., Comment: 22 pages, 9 figures, splitted into two parts: main paper and methods

We spot a hole in the area of succinct data structures for graph classes from a universe of size at most $n^n$. Very often, the input graph is labeled by the user in an arbitrary and easy-to-use way, and the data structure for the graph relabels the input graph in some way. For any access, the user needs to store these labels or compute the new labels in an online manner. This might require more bits than the information-theoretic minimum of the original graph class, hence, defeating the purpose of succinctness. Given this, the data structure designer must allow the user to access the data structure with the original labels, i.e., relabeling is not allowed. We call such a graph data structure ``accessible''. In this paper, we study the complexity of such accessible data structures for interval graphs, a graph class with information-theoretic minimum less than $n\log n$ bits. - We formalize the concept of "accessibility" (which was implicitly assumed), and propose the "universal interval representation", for interval graphs. - Any data structure for interval graphs in universal interval representation, which supports both adjacency and degree query simultaneously with time cost $t_1$ and $t_2$ respectively, must consume at least $\log_2(n!)+n/(\log n)^{O(t_1+t_2)}$ bits of space. This is also the first lower bound for graph classes with information-theoretic minimum less than $n\log_2n$ bits. - We provide efficient succinct data structures for interval graphs in universal interval representation supporting adjacency query and degree query individually in constant time and space costs. Therefore, two upper bounds together with the lower bound show that the two elementary queries for interval graphs are incompatible with each other in the context of succinct data structure. To the best of our knowledge, this is the first proof of such incompatibility phenomenon.

Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C is modelled together with multi-wavelength data within a synchrotron and synchrotron-self Compton (SSC) scenario. We find that SSC can explain the broadband data well from the optical to the very-high-energy band. For the late-time radio data, a different component is needed to account for the observed emission. Differently from previous GRBs detected in the very-high-energy range, the model for GRB~201216C strongly favors a wind-like medium. The model parameters have values similar to those found in past studies of the afterglows of GRBs detected up to GeV energies., Comment: 13 pages, 6 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Water fountains (WFs) are thought to be objects in the morphological evolution of the circumstellar envelopes of low- and intermediate-mass evolved stars, transitioning from spherically symmetric to asymmetric shapes. We used databases of circumstellar 1612 MHz OH and 22.235 GHz H$_2$O maser sources to search for new WF candidates using the criterion of a larger velocity range of the H$_2$O maser emission compared to that of the OH maser emission. Thus, it is in principle possible to identify WFs with an H$_2$O velocity ranges smaller than those for the previously known WFs. For the OH maser line, we analyzed database entries of 8,474 observations from 2,195 sources, and 6,085 observations from 3,642 sources for H$_2$O maser line. After a close examination of the velocity ranges and line profiles, we identified 11 sources that meet the criterion mentioned above. We examined the IRAS colors of the selected sources and found that two of them (IRAS 19069+0916 and IRAS 19319+2214) are in the color region for post-AGB stars. We find that the maser velocity criterion can discover other astrophysically interesting objects than just WFs. Such objects may include peculiar planetary nebulae with maser emissions and stellar merger remnants., Comment: 34 pages, 10 figures, Accepted for publication in ApJS

The preservation of nuclear spin alignment in plasmas is a prerequisite for important applications, such as energy production through polarized fusion or the acceleration of polarized particle beams. Although this conservation property has been the basis of numerous theoretical papers, it has never been experimentally confirmed. Here, we report on first experimental data from a polarized $^3\mkern-2mu$He target heated by a PW laser pulse, showing evidence for persistence of the nuclear polarization after acceleration to MeV energies. The finding also validates the concept of using pre-polarized targets for experiments at high-power laser facilities., Comment: 5 pages, 2 figures