Search

The Euclid mission is generating a vast amount of imaging data in four broadband filters at high angular resolution. This will allow the detailed study of mass, metallicity, and stellar populations across galaxies, which will constrain their formation and evolutionary pathways. Transforming the Euclid imaging for large samples of galaxies into maps of physical parameters in an efficient and reliable manner is an outstanding challenge. We investigate the power and reliability of machine learning techniques to extract the distribution of physical parameters within well-resolved galaxies. We focus on estimating stellar mass surface density, mass-averaged stellar metallicity and age. We generate noise-free, synthetic high-resolution imaging data in the Euclid photometric bands for a set of 1154 galaxies from the TNG50 cosmological simulation. The images are generated with the SKIRT radiative transfer code, taking into account the complex 3D distribution of stellar populations and interstellar dust attenuation. We use a machine learning framework to map the idealised mock observational data to the physical parameters on a pixel-by-pixel basis. We find that stellar mass surface density can be accurately recovered with a $\leq 0.130 {\rm \,dex}$ scatter. Conversely, stellar metallicity and age estimates are, as expected, less robust, but still contain significant information which originates from underlying correlations at a sub-kpc scale between stellar mass surface density and stellar population properties.

Data-driven machine learning has made significant strides in medical image analysis. However, most existing methods are tailored to specific modalities and assume a particular resolution (often isotropic). This limits their generalizability in clinical settings, where variations in scan appearance arise from differences in sequence parameters, resolution, and orientation. Furthermore, most general-purpose models are designed for healthy subjects and suffer from performance degradation when pathology is present. We introduce UNA (Unraveling Normal Anatomy), the first modality-agnostic learning approach for normal brain anatomy reconstruction that can handle both healthy scans and cases with pathology. We propose a fluid-driven anomaly randomization method that generates an unlimited number of realistic pathology profiles on-the-fly. UNA is trained on a combination of synthetic and real data, and can be applied directly to real images with potential pathology without the need for fine-tuning. We demonstrate UNA's effectiveness in reconstructing healthy brain anatomy and showcase its direct application to anomaly detection, using both simulated and real images from 3D healthy and stroke datasets, including CT and MRI scans. By bridging the gap between healthy and diseased images, UNA enables the use of general-purpose models on diseased images, opening up new opportunities for large-scale analysis of uncurated clinical images in the presence of pathology. Code is available at https://github.com/peirong26/UNA., Comment: 16 pages, 6 figures

The physical and orbital parameters of Trans-Neptunian Objects (TNOs) provide valuable information about the Solar System's formation and evolution. In particular, the characterization of binaries provides insights into the formation mechanisms that may be playing a role at such large distances from the Sun. Studies show two distinct populations, and (38628) Huya occupies an intermediate position between the unequal-size binaries and those with components of roughly equal sizes. In this work, we predicted and observed three stellar occultation events by Huya. Huya and its satellite - S/2012 (38628) 1 - were detected during occultations in March 2021 and again in June 2023. Additionally, an attempt to detect Huya in February 2023 resulted in an additional single-chord detection of the secondary. A spherical body with a minimum diameter of D = 165 km can explain the three single-chord observations and provide a lower limit for the satellite size. The astrometry of Huya's system, as derived from the occultations and supplemented by observations from the Hubble Space Telescope and Keck Observatory, provided constraints on the satellite orbit and the mass of the system. Therefore, assuming the secondary is in an equatorial orbit around the primary, the limb fitting was constrained by the satellite orbit position angle. The system density, calculated by summing the most precise measurement of Huya's volume to the spherical satellite average volume, is $\rho_{1}$ = 1073 $\pm$ 66 kg m$^{-3}$. The density that the object would have assuming a Maclaurin equilibrium shape with a rotational period of 6.725 $\pm$ 0.01 hours is $\rho_{2}$ = 768 $\pm$ 42 kg m$^{-3}$. This difference rules out the Maclaurin equilibrium assumption for the main body shape., Comment: Manuscript accepted by Planetary Science Journal (PSJ)

With users demanding seamless connectivity, handovers (HOs) have become a fundamental element of cellular networks. However, optimizing HOs is a challenging problem, further exacerbated by the growing complexity of mobile networks. This paper presents the first countrywide study of HO optimization, through the prism of Smoothed Online Learning (SOL). We first analyze an extensive dataset from a commercial mobile network operator (MNO) in Europe with more than 40M users, to understand and reveal important features and performance impacts on HOs. Our findings highlight a correlation between HO failures/delays, and the characteristics of radio cells and end-user devices, showcasing the impact of heterogeneity in mobile networks nowadays. We subsequently model UE-cell associations as dynamic decisions and propose a realistic system model for smooth and accurate HOs that extends existing approaches by (i) incorporating device and cell features on HO optimization, and (ii) eliminating (prior) strong assumptions about requiring future signal measurements and knowledge of end-user mobility. Our algorithm, aligned with the O-RAN paradigm, provides robust dynamic regret guarantees, even in challenging environments, and shows superior performance in multiple scenarios with real-world and synthetic data.

The formation time-scales of quiescent galaxies can be estimated in two different ways, by their star formation history and by their chemistry. Previously, both methods yielded conflicting results, especially when considering $\alpha$-enhanced objects. This is primarily due to the time resolution limitations of very old stellar populations that prevent us from accurately constraining their star formation histories. We analyse the JWST observations of the extremely massive galaxy ZF-UDS-7329 at z$\sim$3.2 and show that we can achieve the higher time resolution necessary to match the chemical formation time-scales using stellar population synthesis by studying galaxies at high redshift. We compare it to the well known relic galaxy NGC 1277, arguing that ZF-UDS-7329 is an early Universe example of the cores of present day massive elliptical galaxies or, if left untouched, a relic galaxy., Comment: Accepted for publication in A&A, 5 pages, 4 figures

Generative machine learning models are increasingly being used to design novel proteins for therapeutic and biotechnological applications. However, the current methods mostly focus on the design of proteins with a fixed backbone structure, which leads to their limited ability to account for protein flexibility, one of the crucial properties for protein function. Learning to engineer protein flexibility is problematic because the available data are scarce, heterogeneous, and costly to obtain using computational as well as experimental methods. Our contributions to address this problem are three-fold. First, we comprehensively compare methods for quantifying protein flexibility and identify data relevant to learning. Second, we design and train flexibility predictors utilizing sequential or both sequential and structural information on the input. We overcome the data scarcity issue by leveraging a pre-trained protein language model. Third, we introduce a method for fine-tuning a protein inverse folding model to steer it toward desired flexibility in specified regions. We demonstrate that our method Flexpert-Design enables guidance of inverse folding models toward increased flexibility. This opens up new possibilities for protein flexibility engineering and the development of proteins with enhanced biological activities.

We report on the crystallographic structure of penta-twinned gold nanoparticles. Although gold typically exhibits a face-centered cubic (fcc) lattice, other phases have been reported in some nanoscale systems. We show that the crystallographic system and the lattice parameters of the gold unit cell strongly depend on the nanoparticle geometry, for a wide size range. Specifically, we show that decahedra exhibit a body-centered tetragonal structure (I4/mmm), whereas rods and bipyramids exhibit a body-centered orthorhombic structure (Immm). These changes in the crystallographic structure are explained by the elastic lattice distortions required to close the mismatch gap in penta-twinned nanoparticles, with respect to fcc single-crystal gold nanoparticles. The effects of nanoparticle shape and size on the surface pressure and the subsequent distortions are additionally discussed., Comment: 13 pages, 5 figures

Stellar disk truncations, also referred to as galaxy edges, are key indicators of galactic size, determined by the radial location of the gas density threshold for star formation. Accurately measuring galaxy sizes for millions of galaxies is essential for understanding the physical processes driving galaxy evolution over cosmic time. In this study, we aim to explore the potential of the Segment Anything Model (SAM), a foundation model designed for image segmentation, to automatically identify disk truncations in galaxy images. With the Euclid Wide Survey poised to deliver vast datasets, our goal is to assess SAM's capability to measure galaxy sizes in a fully automated manner. SAM was applied to a labeled dataset of 1,047 disk-like galaxies with $M_* > 10^{10} M_{\odot}$ at redshifts up to $z \sim 1$, sourced from the HST CANDELS fields. We 'euclidized' the HST galaxy images by creating composite RGB images, using the F160W (H-band), F125W (J-band), and F814W + F606W (I-band + V-band) HST filters, respectively. Using these processed images as input for SAM, we retrieved various truncation masks for each galaxy image under different configurations of the input data. We find excellent agreement between the galaxy sizes identified by SAM and those measured manually (i.e., by using the radial positions of the stellar disk edges in galaxy light profiles), with an average deviation of approximately $3\%$. This error reduces to about $1\%$ when excluding problematic cases. Our results highlight the strong potential of SAM for detecting disk truncations and measuring galaxy sizes across large datasets in an automated way. SAM performs well without requiring extensive image preprocessing, labeled training datasets for truncations (used only for validation), fine-tuning, or additional domain-specific adaptations such as transfer learning., Comment: Accepted for publication in Astronomy & Astrophysics (11 pages, 10 figures)

This letter presents the measurement of the energy-dependent neutrino-nucleon cross section in tungsten and the differential flux of muon neutrinos and anti-neutrinos. The analysis is performed using proton-proton collision data at a center-of-mass energy of $13.6 \, {\rm TeV}$ and corresponding to an integrated luminosity of $(65.6 \pm 1.4) \, \mathrm{fb^{-1}}$. Using the active electronic components of the FASER detector, $338.1 \pm 21.0$ charged current muon neutrino interaction events are identified, with backgrounds from other processes subtracted. We unfold the neutrino events into a fiducial volume corresponding to the sensitive regions of the FASER detector and interpret the results in two ways: We use the expected neutrino flux to measure the cross section, and we use the predicted cross section to measure the neutrino flux. Both results are presented in six bins of neutrino energy, achieving the first differential measurement in the TeV range. The observed distributions align with Standard Model predictions. Using this differential data, we extract the contributions of neutrinos from pion and kaon decays.

Cellular networks rely on handovers (HOs) as a fundamental element to enable seamless connectivity for mobile users. A comprehensive analysis of HOs can be achieved through data from Mobile Network Operators (MNOs); however, the vast majority of studies employ data from measurement campaigns within confined areas and with limited end-user devices, thereby providing only a partial view of HOs. This paper presents the first countrywide analysis of HO performance, from the perspective of a top-tier MNO in a European country. We collect traffic from approximately 40M users for 4 weeks and study the impact of the radio access technologies (RATs), device types, and manufacturers on HOs across the country. We characterize the geo-temporal dynamics of horizontal (intra-RAT) and vertical (inter-RATs) HOs, at the district level and at millisecond granularity, and leverage open datasets from the country's official census office to associate our findings with the population. We further delve into the frequency, duration, and causes of HO failures, and model them using statistical tools. Our study offers unique insights into mobility management, highlighting the heterogeneity of the network and devices, and their effect on HOs.

Designing resilient Internet of Things (IoT) systems requires i) identification of IoT Critical Objects (ICOs) such as services, devices, and resources, ii) threat analysis, and iii) mitigation strategy selection. However, the traditional process for designing resilient IoT systems is still manual, leading to inefficiencies and increased risks. In addition, while tools such as ChatGPT could support this manual and highly error-prone process, their use raises concerns over data privacy, inconsistent outputs, and internet dependence. Therefore, we propose RITA, an automated, open-source framework that uses a fine-tuned RoBERTa-based Named Entity Recognition (NER) model to identify ICOs from IoT requirement documents, correlate threats, and recommend countermeasures. RITA operates entirely offline and can be deployed on-site, safeguarding sensitive information and delivering consistent outputs that enhance standardization. In our empirical evaluation, RITA outperformed ChatGPT in four of seven ICO categories, particularly in actuator, sensor, network resource, and service identification, using both human-annotated and ChatGPT-generated test data. These findings indicate that RITA can improve resilient IoT design by effectively supporting key security operations, offering a practical solution for developing robust IoT architectures.

Domain randomization through synthesis is a powerful strategy to train networks that are unbiased with respect to the domain of the input images. Randomization allows networks to see a virtually infinite range of intensities and artifacts during training, thereby minimizing overfitting to appearance and maximizing generalization to unseen data. While powerful, this approach relies on the accurate tuning of a large set of hyper-parameters governing the probabilistic distribution of the synthesized images. Instead of manually tuning these parameters, we introduce Learn2Synth, a novel procedure in which synthesis parameters are learned using a small set of real labeled data. Unlike methods that impose constraints to align synthetic data with real data (e.g., contrastive or adversarial techniques), which risk misaligning the image and its label map, we tune an augmentation engine such that a segmentation network trained on synthetic data has optimal accuracy when applied to real data. This approach allows the training procedure to benefit from real labeled examples, without ever using these real examples to train the segmentation network, which avoids biasing the network towards the properties of the training set. Specifically, we develop both parametric and nonparametric strategies to augment the synthetic images, enhancing the segmentation network's performance. Experimental results on both synthetic and real-world datasets demonstrate the effectiveness of this learning strategy. Code is available at: https://github.com/HuXiaoling/Learn2Synth., Comment: 14 pages, 5 figures

The claustrum is a band-like gray matter structure located between putamen and insula whose exact functions are still actively researched. Its sheet-like structure makes it barely visible in in vivo Magnetic Resonance Imaging (MRI) scans at typical resolutions and neuroimaging tools for its study, including methods for automatic segmentation, are currently very limited. In this paper, we propose a contrast- and resolution-agnostic method for claustrum segmentation at ultra-high resolution (0.35 mm isotropic); the method is based on the SynthSeg segmentation framework (Billot et al., 2023), which leverages the use of synthetic training intensity images to achieve excellent generalization. In particular, SynthSeg requires only label maps to be trained, since corresponding intensity images are synthesized on the fly with random contrast and resolution. We trained a deep learning network for automatic claustrum segmentation, using claustrum manual labels obtained from 18 ultra-high resolution MRI scans (mostly ex vivo). We demonstrated the method to work on these 18 high resolution cases (Dice score = 0.632, mean surface distance = 0.458 mm, and volumetric similarity = 0.867 using 6-fold Cross Validation (CV)), and also on in vivo T1-weighted MRI scans at typical resolutions (~1 mm isotropic). We also demonstrated that the method is robust in a test-retest setting and when applied to multimodal imaging (T2-weighted, Proton Density and quantitative T1 scans). To the best of our knowledge this is the first accurate method for automatic ultra-high resolution claustrum segmentation, which is robust against changes in contrast and resolution. The method is released at https://github.com/chiara-mauri/claustrum_segmentation and as part of the neuroimaging package Freesurfer (Fischl, 2012)., Comment: 14 pages, 10 figures, 3 tables

BabyIAXO is the intermediate stage of the International Axion Observatory (IAXO) to be hosted at DESY. Its primary goal is the detection of solar axions following the axion helioscope technique. Axions are converted into photons in a large magnet that is pointing to the sun. The resulting X-rays are focused by appropriate X-ray optics and detected by sensitive low-background detectors placed at the focal spot. The aim of this article is to provide an accurate quantitative description of the different components (such as the magnet, optics, and X-ray detectors) involved in the detection of axions. Our efforts have focused on developing robust and integrated software tools to model these helioscope components, enabling future assessments of modifications or upgrades to any part of the IAXO axion helioscope and evaluating the potential impact on the experiment's sensitivity. In this manuscript, we demonstrate the application of these tools by presenting a precise signal calculation and response analysis of BabyIAXO's sensitivity to the axion-photon coupling. Though focusing on the Primakoff solar flux component, our virtual helioscope model can be used to test different production mechanisms, allowing for direct comparisons within a unified framework., Comment: 36 pages, 18 figures, 4 tables, Submitted to JHEP

This paper introduces and examines numerical approximation schemes for computing risk budgeting portfolios associated to positive homogeneous and sub-additive risk measures. We employ Mirror Descent algorithms to determine the optimal risk budgeting weights in both deterministic and stochastic settings, establishing convergence along with an explicit non-asymptotic quantitative rate for the averaged algorithm. A comprehensive numerical analysis follows, illustrating our theoretical findings across various risk measures -- including standard deviation, Expected Shortfall, deviation measures, and Variantiles -- and comparing the performance with that of the standard stochastic gradient descent method recently proposed in the literature.

In 1986, Albert proposed a Marsden-Weinstein reduction process for cosymplectic structures. In this paper, we present the limitations of this theory in the application of the reduction of symmetric time-dependent Hamiltonian systems. As a consequence, we conclude that cosymplectic geometry is not appropriate for this reduction. Motived for this fact, we replace cosymplectic structures by more general structures: mechanical presymplectic structures. Then, we develop Marsden-Weinstein reduction for this kind of structures and we apply this theory to interesting examples of time-dependent Hamiltonian systems for which Albert's reduction method doesn't work., Comment: 23 pages, 2 figures

We study, for plane complex branches of genus one, the topological type of its generic polar curve, as a function of the semigroup of values and the Zariski invariant of the branch. We improve some results given by Casas-Alvero in 2023, since we filter the topological type fixed for the branch by the possible values of Zariski invariants.

The recent discovery of superconductivity in infinite-layer (IL, ABO$_2$) nickelates has opened a new avenue to deepen the understanding of high-temperature superconductivity. However, progress in this field is slowed by significant challenges in material synthesis and the scarcity of research groups capable of producing high quality superconducting samples. IL nickelates are obtained from a reduction of the perovskite ABO$_3$ phase, typically achieved by annealing using CaH$_2$ as a reducing agent. Here, we present a new method to synthesize superconducting infinite-layer nickelate Pr$_{0.8}$Sr$_{0.2}$NiO$_2$ thin films using an aluminum overlayer deposited by sputtering as a reducing agent. We systematically optimized the aluminum deposition parameters and obtained superconducting samples reduced either in situ or ex situ (after air exposure of the precursor ABO$_3$ films). A comparison of their crystalline quality and transport properties shows that in situ Al reduction enhances the quality of the superconducting Pr$_{0.8}$Sr$_{0.2}$NiO$_2$ thin films, achieving a maximum superconducting transition temperature $T_{c}^{onset}$ of 17 K, in agreement with the optimum value reported for this compound. This simple synthesis route, much more accessible than existing methods, offers better control and reproducibility over the topotactic transformation, opening new opportunities to gain insights into the physics of superconductivity in nickelates.

Layered perovskite-based compounds offer a range of unconventional properties enabled by their naturally anisotropic structure. While most renowned for the superconductivity observed in the Ruddlesden-Popper phases, many of these layered compounds are also ferroelectric and exhibit a sizeable in-plane polarization. Among these, the Carpy-Galy phases (A${_n}$B${_n}$O$_{3n+2}$), characterized by 110-oriented perovskite planes interleaved with additional oxygen layers, have been debated as platforms for hosting not only a robust polarization but also multiferroicity and polar metallicity. However, the challenges associated with the synthesis of ultrathin Carpy-Galy films and understanding the impact of strain on their properties limit their integration into devices. Addressing this issue, our study focuses on La$_2$Ti$_2$O$_7$, an $n$=4 (A$_2$B$_2$O$_7$) representative of the Carpy-Galy family, exploring its growth and concurrent phase stability on various substrates under different strain conditions. Remarkably, we demonstrate that a 3% tensile strain from DyScO$_3$ (100) substrates promotes a controlled layer-by-layer growth mode, while SrTiO$_3$ (110) and LaAlO$_3$-Sr$_2$TaAlO$_6$ (110), that exert negligible and compressive strains respectively, require post-deposition annealing to achieve similar results. Using scanning probe microscopy, X-ray diffraction, scanning transmission electron microscopy, and polarization switching experiments, we confirm that these films possess exceptional ferroelectric properties, including a polarization of 18 $\mu$C/cm$^2$ - more than three times higher than previously reported - as well as persistence of ferroelectricity down to a single-unit-cell thickness. This study not only advances our understanding of Carpy-Galy phases in thin films but also lays a foundation for their application in advanced ferroelectric device architectures., Comment: 19 pages, 5 main-text figures, 6 figures in Supporting Information

Viruses like influenza have long coevolved with host immune systems, gradually shaping the evolutionary trajectory of these pathogens. Host immune systems develop immunity against circulating strains, which in turn avoid extinction by exploiting antigenic escape mutations that render new strains immune from existing antibodies in the host population. Infected hosts are also mobile, which can spread the virus to regions without developed host immunity, offering additional reservoirs for viral growth. While the effects of migration on long term stability have been investigated, we know little about how antigenic escape coupled with migration changes the survival and spread of emerging viruses. By considering the two processes on equal footing, we show that on short timescales an intermediate host mobility rate increases the survival probability of the virus through antigenic escape. We show that more strongly connected migratory networks decrease the survival probability of the virus. Using data from high traffic airports we argue that current human migration rates are beneficial for viral survival.

We revisit the mean field parametrization of shallow neural networks, using signed measures on unbounded parameter spaces and duality pairings that take into account the regularity and growth of activation functions. This setting directly leads to the use of unbalanced Kantorovich-Rubinstein norms defined by duality with Lipschitz functions, and of spaces of measures dual to those of continuous functions with controlled growth. These allow to make transparent the need for total variation and moment bounds or penalization to obtain existence of minimizers of variational formulations, under which we prove a compactness result in strong Kantorovich-Rubinstein norm, and in the absence of which we show several examples demonstrating undesirable behavior. Further, the Kantorovich-Rubinstein setting enables us to combine the advantages of a completely linear parametrization and ensuing reproducing kernel Banach space framework with optimal transport insights. We showcase this synergy with representer theorems and uniform large data limits for empirical risk minimization, and in proposed formulations for distillation and fusion applications., Comment: 39 pages, 1 table

Characterizing the dimensionality of the superconducting state in the infinite-layer (IL) nickelates is crucial to understand its nature. Most studies have addressed the problem by studying the anisotropy of the upper critical fields. Yet, the dominance of Pauli-paramagnetism effects over orbital ones makes it challenging to interpret the experiments in terms of dimensionality. Here we address the question from a different perspective, by investigating the vortex phase diagram in the mixed-state. We demonstrate that superconducting Pr0.8Sr0.2NiO2 thin films present a vortex liquid-to-glass transition of a purely two-dimensional nature. The obtained results suggest that bidimensionality is an intrinsic property, and that superconductivity resides in fully-decoupled NiO2 planes. In this scenario, the coherence length along the c-axis must be shorter than the distance between those planes, while Josephson and magnetostatic coupling between them must be negligible. We believe that these conclusions are relevant for theories on the origin of superconductivity in the IL-nickelates.

The first FASER search for a light, long-lived particle decaying into a pair of photons is reported. The search uses LHC proton-proton collision data at $\sqrt{s}=13.6~\text{TeV}$ collected in 2022 and 2023, corresponding to an integrated luminosity of $57.7\text{fb}^{-1}$. A model with axion-like particles (ALPs) dominantly coupled to weak gauge bosons is the primary target. Signal events are characterised by high-energy deposits in the electromagnetic calorimeter and no signal in the veto scintillators. One event is observed, compared to a background expectation of $0.44 \pm 0.39$ events, which is entirely dominated by neutrino interactions. World-leading constraints on ALPs are obtained for masses up to $300~\text{MeV}$ and couplings to the Standard Model W gauge boson, $g_{aWW}$, around $10^{-4}$ GeV$^{-1}$, testing a previously unexplored region of parameter space. Other new particle models that lead to the same experimental signature, including ALPs coupled to gluons or photons, U(1)$_B$ gauge bosons, up-philic scalars, and a Type-I two-Higgs doublet model, are also considered for interpretation, and new constraints on previously viable parameter space are presented in this paper., Comment: 37 pages, 22 figures

In this work, we propose to use the Reduced-Basis Method (RBM) as a model order reduction approach to solve Maxwell's equations in electromagnetic (EM) scatterers based on plasma to build a metasurface, taking into account a parameter, namely, the plasma frequency. We build up the reduced-order model in an adaptive fashion following a greedy algorithm. This method enables a fast sweep over a wide range of plasma frequencies, thus providing an efficient way to characterize electromagnetic structures based on Drude-like plasma scatterers. We validate and test the proposed technique on several plasma metasurfaces and compare it with the finite element method (FEM) approach.

The notion of defects in crystalline phases of matter has been extremely powerful for understanding crystal growth, deformation and melting. These discontinuities in the periodic order of crystals are mathematically described by the Burgers vector, derived from the particle displacements, which encapsulates the direction and magnitude of slip relative to the undeformed state. Since the reference structure of the crystal is known a priori, the Burgers vector can be determined experimentally using both imaging and diffraction methods to measure the lattice distortion, and thus infer the particle displacements. Glasses have structures that lack the periodicity of crystals, and thus a well-defined reference state. Yet, measurable structural parameters can still be obtained from diffraction from a glass. Here we examine the usefulness of these parameters to probe deformation in glasses. We find that coordinated transformations in the centrosymmetry of local particle arrangements are a strong marker of plastic events. Moreover, we investigate two geometric indicators that can be derived from distortions in local diffraction patterns, namely the continuous Burgers vector and the quadrupolar strain. We find that the Burgers vector again emerges as a robust and sensitive metric for understanding local structural transformations due to mechanical deformation, even in structural glasses., Comment: 13 pages 5 figures

Over recent years, deep learning based image registration has achieved impressive accuracy in many domains, including medical imaging and, specifically, human neuroimaging with magnetic resonance imaging (MRI). However, the uncertainty estimation associated with these methods has been largely limited to the application of generic techniques (e.g., Monte Carlo dropout) that do not exploit the peculiarities of the problem domain, particularly spatial modeling. Here, we propose a principled way to propagate uncertainties (epistemic or aleatoric) estimated at the level of spatial location by these methods, to the level of global transformation models, and further to downstream tasks. Specifically, we justify the choice of a Gaussian distribution for the local uncertainty modeling, and then propose a framework where uncertainties spread across hierarchical levels, depending on the choice of transformation model. Experiments on publicly available data sets show that Monte Carlo dropout correlates very poorly with the reference registration error, whereas our uncertainty estimates correlate much better. % with the reference registration error. Crucially, the results also show that uncertainty-aware fitting of transformations improves the registration accuracy of brain MRI scans. Finally, we illustrate how sampling from the posterior distribution of the transformations can be used to propagate uncertainties to downstream neuroimaging tasks. Code is available at: https://github.com/HuXiaoling/Regre4Regis., Comment: 15 pages, 6 figures

While functional Magnetic Resonance Imaging (fMRI) offers valuable insights into cognitive processes, its inherent spatial limitations pose challenges for detailed analysis of the fine-grained functional architecture of the brain. More specifically, MRI scanner and sequence specifications impose a trade-off between temporal resolution, spatial resolution, signal-to-noise ratio, and scan time. Deep Learning (DL) Super-Resolution (SR) methods have emerged as a promising solution to enhance fMRI resolution, generating high-resolution (HR) images from low-resolution (LR) images typically acquired with lower scanning times. However, most existing SR approaches depend on supervised DL techniques, which require training ground truth (GT) HR data, which is often difficult to acquire and simultaneously sets a bound for how far SR can go. In this paper, we introduce a novel self-supervised DL SR model that combines a DL network with an analytical approach and Total Variation (TV) regularization. Our method eliminates the need for external GT images, achieving competitive performance compared to supervised DL techniques and preserving the functional maps., Comment: Preprint Submitted to ISBI 2025

While the physics of disordered packing in non-growing systems is well understood, unexplored phenomena can emerge when packing takes place in growing domains. We study the arrangements of pigment cells (chromatophores) on squid skin as a biological example of a packed system on an expanding surface. We find that relative density fluctuations in cell numbers grow with spatial scale. We term this behavior ''hyperdisordered'', in contrast with hyperuniform behavior in which relative fluctuations tend to zero at large scale. We find that hyperdisordered scaling, akin to that of a critical system, is quantitatively reproduced by a model in which hard disks are randomly inserted in a homogeneously growing surface. In addition, we find that chromatophores increase in size during animal development, but maintain a stationary size distribution. The physical mechanisms described in our work may apply to a broad class of growing dense systems.

The Optional Public Goods Game is a three-strategy game in which an individual can play as a cooperator or defector or decide not to participate. Despite its simplicity, this model can effectively represent many human social dilemmas, such as those found in the use of public services, environmental concerns, or other activities related to society. In this contribution, we present a comprehensive analysis of the conditions under which spontaneous, sustained cooperation emerges and the characteristics of these cooperative states. Through simulations, we demonstrate the conditions leading to the coexistence of the three strategies in a steady equilibrium or the alternate dominance of each strategy in a rock-paper-scissors fashion. The results identify each of the possible scenarios in terms of two key parameters: the multiplication rate of the public good game (reward) and the size of the group of potential players. We also discuss other details of the game that may influence the appearance of cycles, along with relevant characteristics of these cycles, such as the prevalence of cooperation., Comment: 8 pages, 6 figures

We explore the use of glide symmetry (GS) and electromagnetic bandgap (EBG) technology in a glide-symmetric double corrugated gap waveguide (GSDC-GW) slow wave structure (SWS) for traveling wave tube (TWT) applications. Notably, this GS structure provides the advantage of wide-band operation and the EBG eliminates the need for a conductive connection between the top and bottom waveguide plates. The TWT performance is evaluated via particle-in-cell simulations that reveal a 3-dB bandwidth of approximately 12 GHz spanning from 54.5 GHz to 66.3 GHz, accompanied by a maximum gain of 23 dB. Because of GS, the backward wave in the first spatial harmonic is not longitudinally polarized, leading to a low risk of backward wave oscillations in the TWT. This work places the GSDC-EBG structure within the arena of potential SWS topologies for TWTs operating under similar conditions., Comment: 7 pages, 12 figures

Surface-based analysis of the cerebral cortex is ubiquitous in human neuroimaging with MRI. It is crucial for cortical registration, parcellation, and thickness estimation. Traditionally, these analyses require high-resolution, isotropic scans with good gray-white matter contrast, typically a 1mm T1-weighted scan. This excludes most clinical MRI scans, which are often anisotropic and lack the necessary T1 contrast. To enable large-scale neuroimaging studies using vast clinical data, we introduce recon-all-clinical, a novel method for cortical reconstruction, registration, parcellation, and thickness estimation in brain MRI scans of any resolution and contrast. Our approach employs a hybrid analysis method that combines a convolutional neural network (CNN) trained with domain randomization to predict signed distance functions (SDFs) and classical geometry processing for accurate surface placement while maintaining topological and geometric constraints. The method does not require retraining for different acquisitions, thus simplifying the analysis of heterogeneous clinical datasets. We tested recon-all-clinical on multiple datasets, including over 19,000 clinical scans. The method consistently produced precise cortical reconstructions and high parcellation accuracy across varied MRI contrasts and resolutions. Cortical thickness estimates are precise enough to capture aging effects independently of MRI contrast, although accuracy varies with slice thickness. Our method is publicly available at https://surfer.nmr.mgh.harvard.edu/fswiki/recon-all-clinical, enabling researchers to perform detailed cortical analysis on the huge amounts of already existing clinical MRI scans. This advancement may be particularly valuable for studying rare diseases and underrepresented populations where research-grade MRI data is scarce., Comment: 16 pages in the manuscript with 11 page supplementary material

In an attempt to further characterise the near-Earth asteroid (NEA) population we present 38 new light-curves acquired between September 2020 and November 2023 for NEAs (7335) 1989 JA, (7822) 1991 CS, (154244) 2002 KL6 and (159402) 1999 AP10, obtained from observations taken at the Teide Observatory (Tenerife, Spain). With these new observations along with archival data, we computed their first shape models and spin solutions by applying the light curve inversion method. The obtained rotation periods are in good agreement with those reported in previous works, with improved uncertainties. Additionally, besides the constant period models for (7335) 1989 JA, (7822) 1991 CS and (159402) 1999 AP10, our results for (154244) 2002 KL6 suggest that it could be affected by a Yarkovsky-O'Keefe-Radzievskii-Paddack acceleration with a value of $\upsilon \simeq -7 \times 10^{-9}$ rad d$^{-2}$. This would be one of the first detections of this effect slowing down an asteroid.

Streaming data analysis is increasingly required in applications, e.g., IoT, cybersecurity, robotics, mechatronics or cyber-physical systems. Despite its relevance, it is still an emerging field with open challenges. SDO is a recent anomaly detection method designed to meet requirements of speed, interpretability and intuitive parameterization. In this work, we present SDOoop, which extends the capabilities of SDO's streaming version to retain temporal information of data structures. SDOoop spots contextual anomalies undetectable by traditional algorithms, while enabling the inspection of data geometries, clusters and temporal patterns. We used SDOoop to model real network communications in critical infrastructures and extract patterns that disclose their dynamics. Moreover, we evaluated SDOoop with data from intrusion detection and natural science domains and obtained performances equivalent or superior to state-of-the-art approaches. Our results show the high potential of new model-based methods to analyze and explain streaming data. Since SDOoop operates with constant per-sample space and time complexity, it is ideal for big data, being able to instantly process large volumes of information. SDOoop conforms to next-generation machine learning, which, in addition to accuracy and speed, is expected to provide highly interpretable and informative models.