Search

Traversing rough terrain requires dynamic bipeds to stabilize themselves through foot placement without stepping in unsafe areas. Planning these footsteps online is challenging given non-convexity of the safe terrain, and imperfect perception and state estimation. This paper addresses these challenges with a full-stack perception and control system for achieving underactuated walking on discontinuous terrain. First, we develop model-predictive footstep control (MPFC), a single mixed-integer quadratic program which assumes a convex polygon terrain decomposition to optimize over discrete foothold choice, footstep position, ankle torque, template dynamics, and footstep timing at over 100 Hz. We then propose a novel approach for generating convex polygon terrain decompositions online. Our perception stack decouples safe-terrain classification from fitting planar polygons, generating a temporally consistent terrain segmentation in real time using a single CPU thread. We demonstrate the performance of our perception and control stack through outdoor experiments with the underactuated biped Cassie, achieving state of the art perceptive bipedal walking on discontinuous terrain. Supplemental Video: https://youtu.be/eCOD1bMi638

Successful operation of large particle detectors like the Compact Muon Solenoid (CMS) at the CERN Large Hadron Collider requires rapid, in-depth assessment of data quality. We introduce the ``AutoDQM'' system for Automated Data Quality Monitoring using advanced statistical techniques and unsupervised machine learning. Anomaly detection algorithms based on the beta-binomial probability function, principal component analysis, and neural network autoencoder image evaluation are tested on the full set of proton-proton collision data collected by CMS in 2022. AutoDQM identifies anomalous ``bad'' data affected by significant detector malfunction at a rate 4 -- 6 times higher than ``good'' data, demonstrating its effectiveness as a general data quality monitoring tool., Comment: 16 pages, 14 figures

We study a chemotaxis-consumption mechanism, in which some chemical signal and cells density interact each other. In order to control the concentration of such a population, sources involving gradient nonlinearities, which introduce a dampening effect on the model, are considered. Moreover, the system is characterized by nonlinear diffusion and sensitivity terms. We derive conditions on some data of the problem so to ensure the boundedness of related solutions.

In this work we present a novel linear and positivity preserving upwind discontinuous Galerkin (DG) approximation of a class of chemotaxis models with damping gradient nonlinearities. In particular, both a local and a nonlocal model including nonlinear diffusion, chemoattraction, chemorepulsion and logistic growth are considered. Some numerical experiments in the context of chemotactic collapse are presented, whose results are in accordance with the previous analysis of the approximation and show how the blow-up can be prevented by means of the damping gradient term., Comment: 20 pages, 7 figures

At shipping ports, some repetitive maneuvering tasks such as entering/leaving port, transporting goods inside it or just making surveillance activities, can be efficiently and quickly carried out by a domestic pilot according to his experience. This know-how can be seized by Unmanned Surface Vehicles (USV) in order to autonomously replicate the same tasks. However, the inherent nonlinearity of ship trajectories and environmental perturbations as wind or marine currents make it difficult to learn a model and its respective control. We therefore present a data-driven learning and control methodology for USV, which is based on Gaussian Mixture Model, Gaussian Mixture Regression and the Sontag's universal formula. Our approach is capable to learn the nonlinear dynamics as well as guarantee the convergence toward the target with a robust controller. Real data have been collected through experiments with a vessel at the port of Ceuta. The complex trajectories followed by an expert have been learned including the robust controller. The effect of the controller over noise/perturbations are presented, a measure of error is used to compare estimates and real data trajectories, and finally, an analysis of computational complexity is performed., Comment: Preprint submitted to journal (under review). 25 pages, 13 figures, 3 tables

We present the discovery of 11 new transiting brown dwarfs and low-mass M-dwarfs from NASA's TESS mission: TOI-2844, TOI-3122, TOI-3577, TOI-3755, TOI-4462, TOI-4635, TOI-4737, TOI-4759, TOI-5240, TOI-5467, and TOI-5882. They consist of 5 brown dwarf companions and 6 very low mass stellar companions ranging in mass from $25 M_{\rm J}$ to $128 M_{\rm J}$. We used a combination of photometric time-series, spectroscopic, and high resolution imaging follow-up as a part of the TESS Follow-up Observing Program (TFOP) in order to characterize each system. With over 50 transiting brown dwarfs confirmed, we now have a large enough sample to directly test different formation and evolutionary scenarios. We provide a renewed perspective on the transiting brown dwarf desert and its role in differentiating between planetary and stellar formation mechanisms. Our analysis of the eccentricity distribution for the transiting brown dwarf sample does not support previous claims of a transition between planetary and stellar formation at $\sim42$ $M_{\rm J}$. We also contribute a first look into the metallicity distribution of transiting companions in the range $7 - 150$ $M_{\rm J}$, showing that this too does not support a $\sim42$ $M_{\rm J}$ transition. Finally, we also detect a significant lithium absorption feature in one of the brown dwarf hosts (TOI-5882) but determine that the host star is likely old based on rotation, kinematic, and photometric measurements. We therefore claim that TOI-5882 may be a candidate for planetary engulfment., Comment: Submitted, 32 pages, 16 figures

In this work, we study the $\Xi_c$ baryon electromagnetic decay widths within the constituent quark model formalism through an analysis of the transitions from $P-$wave states to ground states. We use the non-relativistic limit of the Hamiltonian of the electromagnetic interaction on keeping all the terms up to the order $m_j^{-1}$. We calculate the electromagnetic decay widths analytically, for the first time, without any further approximation. Specifically, our theoretical results for the $\Xi_c(2790)^{+/0}$ and $\Xi_c(2815)^{+/0}$ radiative decay widths, without the introduction of any additional parameters, display a significant agreement with the recent experimental values obtained by the Belle experiment. The agreement is due to the fact that we have not introduced further approximation to simplify the calculation of the difficult convective term, unlike what was done in previous studies. Our predictions may be useful for future experiments at the Belle, BABAR, and LHC experiments., Comment: 5 pages

Context. The supernova remnant (SNR) W44 and its surroundings are a prime target for studying the acceleration of cosmic rays (CRs). Several previous studies established an extended gamma-ray emission that is set apart from the radio shell of W44. This emission is thought to originate from escaped high-energy CRs that interact with a surrounding dense molecular cloud complex. Aims. We present a detailed analysis of Fermi-LAT data with an emphasis on the spatial and spectral properties of W44 and its surroundings. We also report the results of the observations performed with the MAGIC telescopes of the northwestern region of W44. Finally, we present an interpretation model to explain the gamma-ray emission of the SNR and its surroundings. Methods. We first performed a detailed spatial analysis of 12 years of Fermi-LAT data at energies above 1 GeV, in order to exploit the better angular resolution, while we set a threshold of 100MeV for the spectral analysis. We performed a likelihood analysis of 174 hours of MAGIC data above 130 GeV using the spatial information obtained with Fermi-LAT. Results. The combined spectra of Fermi-LAT and MAGIC, extending from 100MeV to several TeV, were used to derive constraints on the escape of CRs. Using a time-dependent model to describe the particle acceleration and escape from the SNR, we show that the maximum energy of the accelerated particles has to be ' 40 GeV. However, our gamma-ray data suggest that a small number of lower-energy particles also needs to escape. We propose a novel model, the broken-shock scenario, to account for this effect and explain the gamma-ray emission.

This paper represents one contribution to a larger Roadmap article reviewing the current status of the FHI-aims code. In this contribution, the implementation of polarization, Born-effective charges and topological invariants using a Berry-phase approach in a all-electron, numeric atom-centered orbitals framework is summarized. Guidelines on usage and links to tutorials are provided.

We present results from a high cadence multi-wavelength observational campaign of the enigmatic changing look AGN 1ES 1927+654 from May 2022- April 2024, coincident with an unprecedented radio flare (an increase in flux by a factor of $\sim 60$ over a few months) and the emergence of a spatially resolved jet at $0.1-0.3$ pc scales (Meyer et al. 2024). Companion work has also detected a recurrent quasi-periodic oscillation (QPO) in the $2-10$ keV band with an increasing frequency ($1-2$ mHz) over the same period (Masterson et al., 2025). During this time, the soft X-rays ($0.3-2$ keV) monotonically increased by a factor of $\sim 8$, while the UV emission remained near-steady with $<30\%$ variation and the $2-10$ keV flux showed variation by a factor $\lesssim 2$. The weak variation of the $2-10$ keV X-ray emission and the stability of the UV emission suggest that the magnetic energy density and accretion rate are relatively unchanged, and that the jet could be launched due to a reconfiguration of the magnetic field (toroidal to poloidal) close to the black hole. Advecting poloidal flux onto the event horizon would trigger the Blandford-Znajek (BZ) mechanism, leading to the onset of the jet. The concurrent softening of the coronal slope (from $\Gamma= 2.70\pm 0.04$ to $\Gamma=3.27\pm 0.04$), the appearance of a QPO, and low coronal temperature ($kT_{e}=8_{-3}^{+8}$ keV) during the radio outburst suggest that the poloidal field reconfiguration can significantly impact coronal properties and thus influence jet dynamics. These extraordinary findings in real time are crucial for coronal and jet plasma studies, particularly as our results are independent of coronal geometry., Comment: Submitted to ApJ after minor referee comments

Despite the vast amount of information encoded in Knowledge Graphs (KGs), information about the class affiliation of entities remains often incomplete. Graph Convolutional Networks (GCNs) have been shown to be effective predictors of complete information about the class affiliation of entities in KGs. However, these models do not learn the class affiliation of entities in KGs incorporating the complexity of the task, which negatively affects the models prediction capabilities. To address this problem, we introduce a Markov process-based architecture into well-known GCN architectures. This end-to-end network learns the prediction of class affiliation of entities in KGs within a Markov process. The number of computational steps is learned during training using a geometric distribution. At the same time, the loss function combines insights from the field of evidential learning. The experiments show a performance improvement over existing models in several studied architectures and datasets. Based on the chosen hyperparameters for the geometric distribution, the expected number of computation steps can be adjusted to improve efficiency and accuracy during training.

The BL Lacertae object VER J0521+211 underwent a notable flaring episode in February 2020. A short-term monitoring campaign, led by the MAGIC (Major Atmospheric Gamma Imaging Cherenkov) collaboration, covering a wide energy range from radio to very-high-energy (VHE, 100 GeV < E < 100 TeV) gamma rays was organised to study its evolution. These observations resulted in a consistent detection of the source over six consecutive nights in the VHE gamma-ray domain. Combining these nightly observations with an extensive set of multiwavelength data made modelling of the blazar's spectral energy distribution (SED) possible during the flare. This modelling was performed with a focus on two plausible emission mechanisms: i) a leptonic two-zone synchrotron-self-Compton scenario, and ii) a lepto-hadronic one-zone scenario. Both models effectively replicated the observed SED from radio to the VHE gamma-ray band. Furthermore, by introducing a set of evolving parameters, both models were successful in reproducing the evolution of the fluxes measured in different bands throughout the observing campaign. Notably, the lepto-hadronic model predicts enhanced photon and neutrino fluxes at ultra-high energies (E > 100 TeV). While the photon component, generated via decay of neutral pions, is not directly observable as it is subject to intense pair production (and therefore extinction) through interactions with the cosmic microwave background photons, neutrino detectors (e.g. IceCube) can probe the predicted neutrino component. Finally, the analysis of the gamma-ray spectra, as observed by MAGIC and the Fermi-LAT telescopes, yielded a conservative 95\% confidence upper limit of z \leq 0.244 for the redshift of this blazar., Comment: Accepted for publication in A&A

In deep-inelastic positron-proton scattering, the lepton-jet azimuthal angular asymmetry is measured using data collected with the H1 detector at HERA. When the average transverse momentum of the lepton-jet system, $\lvert \vec{P}_\perp \rvert $, is much larger than the total transverse momentum of the system, $\lvert \vec{q}_\perp \rvert$, the asymmetry between parallel and antiparallel configurations, $\vec{P}_\perp$ and $\vec{q}_\perp$, is expected to be generated by initial and final state soft gluon radiation and can be predicted using perturbation theory. Quantifying the angular properties of the asymmetry therefore provides an additional test of the strong force. Studying the asymmetry is important for future measurements of intrinsic asymmetries generated by the proton's constituents through Transverse Momentum Dependent (TMD) Parton Distribution Functions (PDFs), where this asymmetry constitutes a dominant background. Moments of the azimuthal asymmetries are measured using a machine learning method for unfolding that does not require binning., Comment: 15 pages, 4 figures, 2 tables, updated affiliations and acknowledgements

We present a comprehensive evaluation of a2z-1, an artificial intelligence (AI) model designed to analyze abdomen-pelvis CT scans for 21 time-sensitive and actionable findings. Our study focuses on rigorous assessment of the model's performance and generalizability. Large-scale retrospective analysis demonstrates an average AUC of 0.931 across 21 conditions. External validation across two distinct health systems confirms consistent performance (AUC 0.923), establishing generalizability to different evaluation scenarios, with notable performance in critical findings such as small bowel obstruction (AUC 0.958) and acute pancreatitis (AUC 0.961). Subgroup analysis shows consistent accuracy across patient sex, age groups, and varied imaging protocols, including different slice thicknesses and contrast administration types. Comparison of high-confidence model outputs to radiologist reports reveals instances where a2z-1 identified overlooked findings, suggesting potential for quality assurance applications.

Radiologists face increasing workload pressures amid growing imaging volumes, creating risks of burnout and delayed reporting times. While artificial intelligence (AI) based automated radiology report generation shows promise for reporting workflow optimization, evidence of its real-world impact on clinical accuracy and efficiency remains limited. This study evaluated the effect of draft reports on radiology reporting workflows by conducting a three reader multi-case study comparing standard versus AI-assisted reporting workflows. In both workflows, radiologists reviewed the cases and modified either a standard template (standard workflow) or an AI-generated draft report (AI-assisted workflow) to create the final report. For controlled evaluation, we used GPT-4 to generate simulated AI drafts and deliberately introduced 1-3 errors in half the cases to mimic real AI system performance. The AI-assisted workflow significantly reduced average reporting time from 573 to 435 seconds (p=0.003), without a statistically significant difference in clinically significant errors between workflows. These findings suggest that AI-generated drafts can meaningfully accelerate radiology reporting while maintaining diagnostic accuracy, offering a practical solution to address mounting workload challenges in clinical practice.

Due to its peculiar and highly variable nature, the blazar 3C 454.3 has been extensively monitored by the WEBT team. Here, we present for the first time these long-term optical flux and color variability results using data acquired in B, V, R, and I bands over a time span of $\sim$ 2 decades. We include data from WEBT collaborators and public archives such as SMARTS, Steward Observatory, and ZTF. The data are binned and segmented to study the source over this long term when more regular sampling was available. During our study, the long-term spectral variability reveals a redder when brighter (RWB) trend, which, however, stabilizes at a particular brightness cutoff $\sim$ 14.5 mag in the I-band, after which it saturates and evolves into a complex state. This trend indicates increasing jet emission dominance over accretion disk emission until jet emission completely dominates. Plots of the spectral index variation (following $F_{\nu} \propto \nu^{-\alpha}$) reveal a bimodal distribution using a one-day binning. These correlate with two extreme phases of 3C 454.3, an outburst or high flux state and quiescent or low flux state, which are respectively jet and accretion disk dominated. We have also conducted intra-day variability studies of nine light curves and found that six of them are variable. Discrete Correlation Function (DCF) analysis between different optical waveband pairs peak at zero lags, indicating co-spatial emission in different optical bands., Comment: 18 pages, 6 figures, 5 tables

Atmospheric models demand a lot of computational power and solving the chemical processes is one of its most computationally intensive components. This work shows how to improve the computational performance of the Multiscale Online Nonhydrostatic AtmospheRe CHemistry model (MONARCH), a chemical weather prediction system developed by the Barcelona Supercomputing Center. The model implements the new flexible external package Chemistry Across Multiple Phases (CAMP) for the solving of gas- and aerosol-phase chemical processes, that allows multiple chemical processes to be solved simultaneously as a single system. We introduce a novel strategy to simultaneously solve multiple instances of a chemical mechanism, represented in the model as grid-cells, obtaining a speedup up to 9x using thousands of cells. In addition, we present a GPU strategy for the most time-consuming function of CAMP. The GPU version achieves up to 1.2x speedup compared to CPU. Also, we optimize the memory access in the GPU to increase its speedup up to 1.7x.

Two substantial technological advances have reshaped the public square in recent decades: first with the advent of the internet and second with the recent introduction of large language models (LLMs). LLMs offer opportunities for a paradigm shift towards more decentralized, participatory online spaces that can be used to facilitate deliberative dialogues at scale, but also create risks of exacerbating societal schisms. Here, we explore four applications of LLMs to improve digital public squares: collective dialogue systems, bridging systems, community moderation, and proof-of-humanity systems. Building on the input from over 70 civil society experts and technologists, we argue that LLMs both afford promising opportunities to shift the paradigm for conversations at scale and pose distinct risks for digital public squares. We lay out an agenda for future research and investments in AI that will strengthen digital public squares and safeguard against potential misuses of AI., Comment: 40 pages, 5 figures

M dwarfs have become increasingly important in the detection of exoplanets and the study of Earth-sized planets and their habitability. However, 20-30% of M dwarfs have companions that can impact the formation and evolution of planetary systems. We use high-resolution imaging and Gaia astrometry to detect stellar companions around M dwarf exoplanet hosts discovered by TESS and determine the projected separation and estimated stellar masses for each system. We find 47 companions around 216 M dwarfs and a multiplicity rate of $19.4\pm2.7$% that is consistent with field M dwarfs. The binary projected separation distribution is shifted to larger separations, confirming the lack of close binaries hosting transiting exoplanets seen in previous studies. We correct the radii of planets with nearby companions and examine the properties of planets in M dwarf multi-star systems. We also note three multi-planet systems that occur in close binaries ($\lesssim 50$ au) where planet formation is expected to be suppressed., Comment: 28 pages, 13 figures, 5 tables, accepted for publication in The Astronomical Journal

The previously-unknown nucleus 20Al has been observed for the first time by detecting its in-flight decays. Tracking trajectories of all decay products with silicon micro-strip detectors allowed for a conclusion that 20Al is unbound with respect to three-proton (3p) emission. The 3p-decay energy of 20Al ground state has been determined to be 1.93(+0.11,-0.09) MeV through a detailed study of angular correlations of its decay products, 17Ne+p+p+p. This value is much smaller in comparison with the predictions inferred from the isospin symmetry by using the known energy of its mirror nucleus 20N, which indicates a possible mirror symmetry violation in the structure of 3p emitters. Such an isospin symmetry breaking is supported by the calculations of the continuum embedded theoretical frameworks, describing the observed 20Al ground state as an 1p s-wave state with a spin-parity of 1-, which contradicts to the spin-parity (2-) of the 20N ground state. The 20Al ground state decays by sequential 1p-2p emission via intermediate ground state of 19Mg, which is the first observed case of daughter two-proton radioactivity following 1p decay of the parent state., Comment: 9 pages, 6 figures, to be submitted to APS

The BL Lac object 3C 371 is one of the targets that are regularly monitored by the Whole Earth Blazar Telescope (WEBT) Collaboration to study blazar variability on both short and long timescales. We aim to evaluate the long-term multiwavelength (MWL) behaviour of 3C 371, comparing it with the results derived for its optical emission in our previous study. For this, we make use of the multi-band campaigns organized by the WEBT Collaboration in optical and radio between January 2018 and December 2020, and of public data from Swift and Fermi satellites and the MOJAVE Very Large Interferometry programme. We evaluate the variability shown by the source in each band with the amplitude variability quantification, as well as possible interband correlation using the z-Discrete Correlation Function. We also present a deep analysis of the optical-UV, X-ray and $\gamma$-ray spectral variability. With the MOJAVE data we perform a kinematics analysis, looking for components propagating along the jet, calculating its kinematics parameters. This set of parameters is later used for the interpretation of the source MWL behaviour, modelling the broadband spectral energy distribution (SED) of the source with theoretical blazar emission scenarios., Comment: Accepted in A&A, 23 pages, 15 figures

This paper introduces a novel methodology for the mathematical modelling of first and second order phase transformations. It will be shown that this methodology can be related to certain limiting cases of the Cahn-Hilliard equation, specifically the cases of having (i) a convex molar free energy function and (ii) a convex molar free energy function with no regularization. The latter case is commonly regarded as unstable; however, by modifying the variational approach and solving for rate-dependent variables, we obtain a stabilized method capable of handling the missing regularization. While the specific numerical method used to solve the equations (a mixed finite element approach) has been previously employed in related contexts (e.g., to stabilize solutions of the Laplace equation), its application to diffusion and diffusional phase transformations is novel. We prove the thermodynamic consistency of the derived method and discuss several use cases. Our work contributes to the development of new mathematical tools for modeling complex phase transformations in materials science., Comment: 21 pages, 12 figures, preprint

AI-driven models have demonstrated significant potential in automating radiology report generation for chest X-rays. However, there is no standardized benchmark for objectively evaluating their performance. To address this, we present ReXrank, https://rexrank.ai, a public leaderboard and challenge for assessing AI-powered radiology report generation. Our framework incorporates ReXGradient, the largest test dataset consisting of 10,000 studies, and three public datasets (MIMIC-CXR, IU-Xray, CheXpert Plus) for report generation assessment. ReXrank employs 8 evaluation metrics and separately assesses models capable of generating only findings sections and those providing both findings and impressions sections. By providing this standardized evaluation framework, ReXrank enables meaningful comparisons of model performance and offers crucial insights into their robustness across diverse clinical settings. Beyond its current focus on chest X-rays, ReXrank's framework sets the stage for comprehensive evaluation of automated reporting across the full spectrum of medical imaging.

This chapter delves into the emerging field of neuro-symbolic query optimization for knowledge graphs (KGs), presenting a comprehensive exploration of how neural and symbolic techniques can be integrated to enhance query processing. Traditional query optimizers in knowledge graphs rely heavily on symbolic methods, utilizing dataset summaries, statistics, and cost models to select efficient execution plans. However, these approaches often suffer from misestimations and inaccuracies, particularly when dealing with complex queries or large-scale datasets. Recent advancements have introduced neural models, which capture non-linear aspects of query optimization, offering promising alternatives to purely symbolic methods. In this chapter, we introduce neuro-symbolic query optimizers, a novel approach that combines the strengths of symbolic reasoning with the adaptability of neural computation. We discuss the architecture of these hybrid systems, highlighting the interplay between neural and symbolic components to improve the optimizer's ability to navigate the search space and produce efficient execution plans. Additionally, the chapter reviews existing neural components tailored for optimizing queries over knowledge graphs and examines the limitations and challenges in deploying neuro-symbolic query optimizers in real-world environments.

Superparamagnetic iron-oxide nanoparticles (SPIONs) are promising probes for biomedical imaging, but the heterogeneity of their magnetic properties is difficult to characterize with existing methods. Here, we perform widefield imaging of the stray magnetic fields produced by hundreds of isolated ~30-nm SPIONs using a magnetic microscope based on nitrogen-vacancy centers in diamond. By analyzing the SPION magnetic field patterns as a function of applied magnetic field, we observe substantial field-dependent transverse magnetization components that are typically obscured with ensemble characterization methods. We find negligible hysteresis in each of the three magnetization components for nearly all SPIONs in our sample. Most SPIONs exhibit a sharp Langevin saturation curve, enumerated by a characteristic polarizing applied field, B_c. The B_c distribution is highly asymmetric, with a standard deviation (1.4 mT) that is larger than the median (0.6 mT). Using time-resolved magnetic microscopy, we directly record SPION N\'eel relaxation, after switching off a 31 mT applied field, with a temporal resolution of ~60 ms that is limited by the ring-down time of the electromagnet coils. For small bias fields B_{hold}=1.5-3.5 mT, we observe a broad range of SPION N\'eel relaxation times--from milliseconds to seconds--that are consistent with an exponential dependence on B_{hold}. Our time-resolved diamond magnetic microscopy study reveals rich SPION sample heterogeneity and may be extended to other fundamental studies of nanomagnetism., Comment: Main text: 8 pages, 5 figures. Entire manuscript including Appendices: 28 pages, 22 figures

Future quantum networks will be composed of both terrestrial links for metropolitan and continent-scale connections and space-based links for global coverage and infrastructure resilience. However, the propagation of quantum signals through the atmosphere is severely impacted by the effects of turbulence. This is even more the case for entanglement-based quantum communication protocols requiring two free-space channels to be considered simultaneously. In this work, we assess the advantage of turbulence mitigation by adaptive optics, in particular during daytime link operation, so as to increase the coupling of the received signal into an optical fiber. We show in particular that this improves the performance of entanglement-based quantum key distribution by up to a few hundred bits per second when compared with the uncorrected scenario

We analyze the dynamics of a gravitational bounded cloud of generic bosons or a boson star in the non-relativistic regime, assuming from the very beginning that the boson star is in a Bose-Einstein condensed state. Also, we have taken into account multi-particle interactions within the boson star, which are described through the so-called logarithmic potential. Assuming energy conservation conditions at any time, we are able to describe the boson star as a some kind of gravitational quantum liquid, or more specifically, as a macroscopic gravitational quantum droplet. The assumption of logarithmic interactions, lead us to establish several scenarios in which the system is able to oscillates around the equilibrium radius. The approach analyzed in the present report shows that the interpretation of a boson star as a macroscopic gravitational quantum droplet can be applied to the dynamics analysis of a large number of systems on different size scales., Comment: 15 pages, 10 figures

Generating accurate radiology reports from medical images is a clinically important but challenging task. While current Vision Language Models (VLMs) show promise, they are prone to generating hallucinations, potentially compromising patient care. We introduce RadFlag, a black-box method to enhance the accuracy of radiology report generation. Our method uses a sampling-based flagging technique to find hallucinatory generations that should be removed. We first sample multiple reports at varying temperatures and then use a Large Language Model (LLM) to identify claims that are not consistently supported across samples, indicating that the model has low confidence in those claims. Using a calibrated threshold, we flag a fraction of these claims as likely hallucinations, which should undergo extra review or be automatically rejected. Our method achieves high precision when identifying both individual hallucinatory sentences and reports that contain hallucinations. As an easy-to-use, black-box system that only requires access to a model's temperature parameter, RadFlag is compatible with a wide range of radiology report generation models and has the potential to broadly improve the quality of automated radiology reporting., Comment: 17 pages, 6 figures

At high-energy collider experiments, generative models can be used for a wide range of tasks, including fast detector simulations, unfolding, searches of physics beyond the Standard Model, and inference tasks. In particular, it has been demonstrated that score-based diffusion models can generate high-fidelity and accurate samples of jets or collider events. This work expands on previous generative models in three distinct ways. First, our model is trained to generate entire collider events, including all particle species with complete kinematic information. We quantify how well the model learns event-wide constraints such as the conservation of momentum and discrete quantum numbers. We focus on the events at the future Electron-Ion Collider, but we expect that our results can be extended to proton-proton and heavy-ion collisions. Second, previous generative models often relied on image-based techniques. The sparsity of the data can negatively affect the fidelity and sampling time of the model. We address these issues using point clouds and a novel architecture combining edge creation with transformer modules called Point Edge Transformers. Third, we adapt the foundation model OmniLearn, to generate full collider events. This approach may indicate a transition toward adapting and fine-tuning foundation models for downstream tasks instead of training new models from scratch.

Blazars are beamed active galactic nuclei known for their strong multi-wavelength variability on timescales from years down to minutes. We aim to investigate the suitability of the twisting jet model presented in previous works to explain the multi-wavelength behaviour of BL Lacertae, the prototype of one of the blazar classes. According to this model, the jet is inhomogeneous, curved, and twisting, and the long-term variability is due to changes in the Doppler factor due to variations in the orientation of the jet-emitting regions. We analysed optical data of the source obtained during monitoring campaigns organised by the Whole Earth Blazar Telescope (WEBT) in 2019-2022, together with radio data from the WEBT and other teams, and gamma-ray data from the Fermi satellite. In this period, BL Lacertae underwent an extraordinary activity phase, reaching its historical optical and gamma-ray brightness maxima. The application of the twisting jet model to the source light curves allows us to infer the wiggling motion of the optical, radio, and gamma-ray jet-emitting regions. The optical-radio correlation shows that the changes in the radio viewing angle follow those in the optical viewing angle by about 120 days, and it suggests that the jet is composed of plasma filaments, which is in agreement with some radio high-resolution observations of other sources. The gamma-ray emitting region is found to be co-spatial with the optical one, and the analysis of the gamma-optical correlation is consistent with both the geometric interpretation and a synchrotron self-Compton (SSC) origin of the high-energy photons. We propose a geometric scenario where the jet is made up of a pair of emitting plasma filaments in a sort of double-helix curved rotating structure, whose wiggling motion produces changes in the Doppler beaming and can thus explain the observed multi-wavelength long-term variability., Comment: In press for A&A

This paper presents a new hybrid Quantum Machine Learning (QML) model composed of three elements: a classical computer in charge of the data preparation and interpretation; a Gate-based Quantum Computer running the Variational Quantum Algorithm (VQA) representing the Quantum Neural Network (QNN); and an adiabatic Quantum Computer where the optimization function is executed to find the best parameters for the VQA. As of the moment of this writing, the majority of QNNs are being trained using gradient-based classical optimizers having to deal with the barren-plateau effect. Some gradient-free classical approaches such as Evolutionary Algorithms have also been proposed to overcome this effect. To the knowledge of the authors, adiabatic quantum models have not been used to train VQAs. The paper compares the results of gradient-based classical algorithms against adiabatic optimizers showing the feasibility of integration for gate-based and adiabatic quantum computing models, opening the door to modern hybrid QML approaches for High Performance Computing., Comment: 12 pages, 6 figures, Euro PAR 2024 EuroQHPC Workshop

It is known that the Bessel--Fourier coefficients $f_m$ of a function $f$ such that $\sqrt{x}f(x)$ is integrable over $[0,1]$ satisfy $f_m/\sqrt{m}\to 0$. We show a partial converse, namely that for $0\leq \alpha<1/2$ and any non-negative $a_m\to 0$, there is a function $f$ such that $x^{\alpha+1}f(x)$ is integrable and its Bessel--Fourier coefficients $f_m$ satisfy $m^{-\alpha}f_m\geq a_m$ and $m^{-\alpha}f_m\to 0$. We conjecture that the same should be true when $\alpha=\frac{1}{2}$, and discuss some consequences of this conjecture., Comment: 21 pages, 3 figures

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

This paper explores a unique cave art found in southern Puerto Rico that depicts a comet over a tomb. Through interdisciplinary methods, including art interpretation, historical documentation, and demographic analysis, this study uncovered the artist's identity, the societal context of the period, and the potential motivations behind the creation of this art. The investigation revealed a connection to the passage of Halley's Comet in 1910 and the widespread panic it induced., Comment: 7 pages, 6 figures. Submitted to the Journal of Astronomical History and Heritage

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: 17 pages, 10 figures github: https://github.com/mireianievas/gammapy_mwl_workflow zenodo: https://zenodo.org/records/13837637

The Aitken method is a useful approach for decomposing mid-IR polarimetry of silicates into emissive and absorptive components along the line-of-sight in astronomical sources. Here we extend this method to include the effects of polarization caused by scattering from graphite or similar particles located along the same sightlines. To illustrate the application of the extended method, we apply it in the reduction and analysis of CanariCam multi-wavelength imaging polarimetry observations of the Egg Nebula, W3 IRS5, and W51 IRS2, and also spectropolarimetry of W3 IRS5. We compare these results to those obtained with the original Aitken method and show that the data are generally fit much better when this third component, i.e., scattering, is incorporated into the analysis. Including scattering in the decomposition of polarimetry produces better results for the emissive and absorptive components, and in particular for the position angle (PA) of those polarization components. The distribution of the difference between absorptive PA and emissive PA is then found to be more peaked at a fixed angle, nearly perpendicular. This supports the theory that mid-IR polarization arises from elongated dust grains aligned along magnetic field lines, since then the PA of emissive and absorptive polarization would be approximately perpendicular. When significant scattering is not present the extended Aitken method produces the same results as the original method., Comment: 27 pages, 24 figures, 1 table

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

This paper introduces a new pseudodifferential preconditioner for the Helmholtz equation in variable media with absorption. The pseudodifferential operator is associated with the multiplicative inverse to the symbol of the Helmholtz operator. This approach is well-suited for the intermediate and high-frequency regimes. The main novel idea for the fast evaluation of the preconditioner is to interpolate its symbol, not as a function of the (high-dimensional) phase-space variables, but as a function of the wave speed itself. Since the wave speed is a real-valued function, this approach allows us to interpolate in a univariate setting even when the original problem is posed in a multidimensional physical space. As a result, the needed number of interpolation points is small, and the interpolation coefficients can be computed using the fast Fourier transform. The overall computational complexity is log-linear with respect to the degrees of freedom as inherited from the fast Fourier transform. We present some numerical experiments to illustrate the effectiveness of the preconditioner to solve the discrete Helmholtz equation using the GMRES iterative method. The implementation of an absorbing layer for scattering problems using a complex-valued wave speed is also developed. Limitations and possible extensions are also discussed.

We present Head CT Ontology Normalized Evaluation (HeadCT-ONE), a metric for evaluating head CT report generation through ontology-normalized entity and relation extraction. HeadCT-ONE enhances current information extraction derived metrics (such as RadGraph F1) by implementing entity normalization through domain-specific ontologies, addressing radiological language variability. HeadCT-ONE compares normalized entities and relations, allowing for controllable weighting of different entity types or specific entities. Through experiments on head CT reports from three health systems, we show that HeadCT-ONE's normalization and weighting approach improves the capture of semantically equivalent reports, better distinguishes between normal and abnormal reports, and aligns with radiologists' assessment of clinically significant errors, while offering flexibility to prioritize specific aspects of report content. Our results demonstrate how HeadCT-ONE enables more flexible, controllable, and granular automated evaluation of head CT reports.

Accurately interpreting medical images and writing radiology reports is a critical but challenging task in healthcare. Both human-written and AI-generated reports can contain errors, ranging from clinical inaccuracies to linguistic mistakes. To address this, we introduce ReXErr, a methodology that leverages Large Language Models to generate representative errors within chest X-ray reports. Working with board-certified radiologists, we developed error categories that capture common mistakes in both human and AI-generated reports. Our approach uses a novel sampling scheme to inject diverse errors while maintaining clinical plausibility. ReXErr demonstrates consistency across error categories and produces errors that closely mimic those found in real-world scenarios. This method has the potential to aid in the development and evaluation of report correction algorithms, potentially enhancing the quality and reliability of radiology reporting.