-
High contrast at short separation with VLTI/GRAVITY: Bringing Gaia companions to light
Authors:
N. Pourré,
T. O. Winterhalder,
J. -B. Le Bouquin,
S. Lacour,
A. Bidot,
M. Nowak,
A. -L. Maire,
D. Mouillet,
C. Babusiaux,
J. Woillez,
R. Abuter,
A. Amorim,
R. Asensio-Torres,
W. O. Balmer,
M. Benisty,
J. -P. Berger,
H. Beust,
S. Blunt,
A. Boccaletti,
M. Bonnefoy,
H. Bonnet,
M. S. Bordoni,
G. Bourdarot,
W. Brandner,
F. Cantalloube
, et al. (151 additional authors not shown)
Abstract:
Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working…
▽ More
Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($Δ\mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($Δ\mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
Zero-echo-time sequences in highly inhomogeneous fields
Authors:
Jose Borreguero,
Fernando Galve,
José Miguel Algarín,
Joseba Alonso
Abstract:
Zero echo time (ZTE) sequences have proven a powerful tool for Magnetic Resonance Imaging (MRI) of ultrashort T2 tissues, but they fail to produce useful images in the presence of strong field inhomogeneities. Here we present a method to correct artifacts induced by strong B0 inhomogeneities in non-Cartesian sequences, based on magnetic field maps obtained from the phase difference between two fas…
▽ More
Zero echo time (ZTE) sequences have proven a powerful tool for Magnetic Resonance Imaging (MRI) of ultrashort T2 tissues, but they fail to produce useful images in the presence of strong field inhomogeneities. Here we present a method to correct artifacts induced by strong B0 inhomogeneities in non-Cartesian sequences, based on magnetic field maps obtained from the phase difference between two fast pointwise acquisitions. These are free of geometric distortions and can be used for model based image reconstruction with iterative algebraic techniques. In this paper, we show that distortions and artifacts coming from inhomogeneites in B0 are largely reverted by our method, as opposed to standard Conjugate Phase reconstructions. We base this technique on a Single-Point Double-Shot (SPDS) approach, and we have shown it to work even for intra-voxel bandwidths (determined by B0 inhomogeneities) comparable to the encoding bandwidth (determined by the gradient fields). We benchmark the performance of SPDS for ZTE acquisitions, which can be exploited for e.g. dental imaging in affordable low-field MRI systems. Furthermore, SPDS can be used for arbitrary pulse sequences and it may prove useful for extreme magnet geometries, as in e.g. single-sided MRI.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
IsoDAR@Yemilab: Preliminary Design Report -- Volume I: Cyclotron Driver
Authors:
Daniel Winklehner,
Joshua Spitz,
Jose R. Alonso,
Janet M. Conrad,
Jarrett Moon,
Michel Abs,
Alexander Herrod,
Sébastien De Neuter,
Eric Forton,
Denis Joassin,
Erik Van der Kraaij,
Gil Wéry
Abstract:
This Preliminary Design Report (PDR) describes the IsoDAR electron-antineutrino source. Volumes I and II are site-independent and describe the cyclotron driver providing a 10~mA proton beam, and the medium energy beam transport line and target, respectively. Volume III describes the installation at the Yemilab underground laboratory in South Korea. The IsoDAR driver and target will produce a mole…
▽ More
This Preliminary Design Report (PDR) describes the IsoDAR electron-antineutrino source. Volumes I and II are site-independent and describe the cyclotron driver providing a 10~mA proton beam, and the medium energy beam transport line and target, respectively. Volume III describes the installation at the Yemilab underground laboratory in South Korea. The IsoDAR driver and target will produce a mole of electron-antineutrinos over the course of five years. Paired with a kton-scale liquid scintillator detector, it will enable an impressive particle physics program including searches for new symmetries, new interactions and new particles. Here in Volume I, we describe the driver, which includes the ion source, low energy beam transport, and cyclotron. The latter features radiofrequency quadrupole (RFQ) direct axial injection and represents the first accelerator purpose-built to make use of vortex motion.
△ Less
Submitted 11 April, 2024; v1 submitted 9 April, 2024;
originally announced April 2024.
-
Waveguide QED at the onset of spin-spin correlations
Authors:
Sebastián Roca-Jerat,
Marcos Rubín-Osanz,
Mark D. Jenkins,
Agustín Camón,
Pablo J. Alonso,
David Zueco,
Fernando Luis
Abstract:
We experimentally explore the competition between light-mediated and direct matter-matter interactions in waveguide quantum electrodynamics. For this, we couple a superconducting line to a model magnetic material, made of organic free radical DPPH molecules with a spin $S=1/2$ and a $g_{S}$ factor very close to that of a free electron. The microwave transmission has been measured in a wide range o…
▽ More
We experimentally explore the competition between light-mediated and direct matter-matter interactions in waveguide quantum electrodynamics. For this, we couple a superconducting line to a model magnetic material, made of organic free radical DPPH molecules with a spin $S=1/2$ and a $g_{S}$ factor very close to that of a free electron. The microwave transmission has been measured in a wide range of temperatures ($0.013$ K $\leq T \leq 2$ K), magnetic fields ($0\leq B \leq 0.5$ T) and frequencies ($0 \leq ω/2 π\leq 14$ GHz). We find that molecules belonging to the crystal sublattice B form one-dimensional spin chains. Temperature then controls intrinsic spin correlations along the chain in a continuous and monotonic way. In the paramagnetic region ($T > 0.7$ K), the microwave transmission shows evidences for the collective coupling of quasi-identical spins to the propagating photons, with coupling strengths that reach values close to the dissipation rates. As $T$ decreases, the growth of intrinsic spin correlations, combined with the anisotropy in the spin-spin exchange constants, break down the collective spin-photon coupling. In this regime, the temperature dependence of the spin resonance visibility reflects the change in the nature of the dominant spin excitations, from single spin flips to bosonic magnons, which is brought about by the magnetic correlation growth.
△ Less
Submitted 4 April, 2024;
originally announced April 2024.
-
Conditional integrability and stability for the homogeneous Boltzmann equation with very soft potentials
Authors:
Ricardo J. Alonso,
Pierre Gervais,
Bertrand Lods
Abstract:
We introduce a practical criterion that justifies the propagation and appearance of $L^{p}$-norms for the solutions to the spatially homogeneous Boltzmann equation with very soft potentials without cutoff. Such criterion also provides a new conditional stability result for classical solutions to the equation. All results are quantitative. Our approach is inspired by a recent analogous result for t…
▽ More
We introduce a practical criterion that justifies the propagation and appearance of $L^{p}$-norms for the solutions to the spatially homogeneous Boltzmann equation with very soft potentials without cutoff. Such criterion also provides a new conditional stability result for classical solutions to the equation. All results are quantitative. Our approach is inspired by a recent analogous result for the Landau equation derived in arXiv:2306.15729 and generalises existing conditional results related to higher integrability properties and stability of solutions to the Boltzmann equation with very soft potentials.
△ Less
Submitted 22 March, 2024;
originally announced March 2024.
-
Exotic magnetism in perovskite KOsO3
Authors:
Jie Chen,
Hongze Li,
Javier Gainza,
Angel Munoz,
Jose A. Alonso,
Jue Liu,
Yu-Sheng Chen,
Alexei A. Belik,
Kazunari Yamaura,
Jiaming He,
Xinyu Li,
John B. Goodenough,
J. -S. Zhou
Abstract:
A new perovskite KOsO3 has been stabilized under high-pressure and high temperature conditions. It is cubic at 500 K (Pm-3m) and undergoes subsequent phase transitions to tetragonal at 320 K (P4/mmm) and rhombohedral (R-3m) at 230 K as shown from refining synchrotron X-ray powder diffraction (SXRD) data. The larger orbital overlap integral and the extended wavefunction of 5d electrons in the perov…
▽ More
A new perovskite KOsO3 has been stabilized under high-pressure and high temperature conditions. It is cubic at 500 K (Pm-3m) and undergoes subsequent phase transitions to tetragonal at 320 K (P4/mmm) and rhombohedral (R-3m) at 230 K as shown from refining synchrotron X-ray powder diffraction (SXRD) data. The larger orbital overlap integral and the extended wavefunction of 5d electrons in the perovskite KOsO3 allow to explore physics from the regime where Mott and Hund's rule couplings dominate to the state where the multiple interactions are on equal footing. We demonstrate an exotic magnetic ordering phase found by neutron powder diffraction along with physical properties via a suite of measurements including magnetic and transport properties, differential scanning calorimetry, and specific heat, which provide comprehensive information for a system at the crossover from localized to itinerant electronic behavior.
△ Less
Submitted 19 March, 2024;
originally announced March 2024.
-
Probabilistic Analysis of Multiparameter Persistence Decompositions
Authors:
Ángel Javier Alonso,
Michael Kerber,
Primoz Skraba
Abstract:
Multiparameter persistence modules can be uniquely decomposed into indecomposable summands. Among these indecomposables, intervals stand out for their simplicity, making them preferable for their ease of interpretation in practical applications and their computational efficiency. Empirical observations indicate that modules that decompose into only intervals are rare. To support this observation,…
▽ More
Multiparameter persistence modules can be uniquely decomposed into indecomposable summands. Among these indecomposables, intervals stand out for their simplicity, making them preferable for their ease of interpretation in practical applications and their computational efficiency. Empirical observations indicate that modules that decompose into only intervals are rare. To support this observation, we show that for numerous common multiparameter constructions, such as density- or degree-Rips bifiltrations, and across a general category of point samples, the probability of the homology-induced persistence module decomposing into intervals goes to zero as the sample size goes to infinity.
△ Less
Submitted 18 March, 2024;
originally announced March 2024.
-
Discrete Painlevé equations and pencils of quadrics in $\mathbb P^3$
Authors:
Jaume Alonso,
Yuri B. Suris,
Kangning Wei
Abstract:
Discrete Painlevé equations constitute a famous class of integrable non-autonomous second order difference equations. A classification scheme proposed by Sakai interprets a discrete Painlevé equation as a birational map between generalized Halphen surfaces (surfaces obtained from $\mathbb P^1\times\mathbb P^1$ by blowing up at eight points). Sakai's classification is thus based on the classificati…
▽ More
Discrete Painlevé equations constitute a famous class of integrable non-autonomous second order difference equations. A classification scheme proposed by Sakai interprets a discrete Painlevé equation as a birational map between generalized Halphen surfaces (surfaces obtained from $\mathbb P^1\times\mathbb P^1$ by blowing up at eight points). Sakai's classification is thus based on the classification of generalized Halphen surfaces. We propose a novel geometric interpretation of discrete Painlevé equations, where the family of generalized Halphen surfaces is replaced by a pencil of quadrics in $\mathbb P^3$. A discrete Painlevé equation is viewed as a birational transformation of $\mathbb P^3$ that preserves the pencil and maps each quadric of the pencil to a different one, according to a Möbius transformation of the pencil parameter. Thus, our scheme is based on the classification of pencils of quadrics in $\mathbb P^3$.
While historically discrete Painlevé equations appeared as de-autonomizations of QRT maps, in our scheme they are viewed as deformations of 3D QRT maps introduced in our previous paper. A 3D QRT map is defined geometrically as a composition of involutions along generators of quadrics of a pencil, preserving the intersection curves with a second pencil of quadrics. The base set of the net of quadrics spanned by both pencils consists of eight points (which play the role of the eight blow-up points of generalized Halphen surfaces). A Painlevé deformation of a 3D QRT map is obtained by composing involutions along generators with a birational (often linear) transformation of $\mathbb P^3$ under which the pencil remains invariant, but the individual quadrics are mapped according to a Möbius transformation of the pencil parameter.
△ Less
Submitted 17 March, 2024;
originally announced March 2024.
-
Elliptical Halbach magnet and gradient modules for low-field portable MRI
Authors:
Fernando Galve,
Eduardo Pallás,
Teresa Guallart-Naval,
Pablo García-Cristóbal,
Pablo Martínez,
José M. Algarín,
José Borreguero,
Rubén Bosch,
Francisco Juan-Lloris,
José M. Benlloch,
Joseba Alonso
Abstract:
Objective. To develop methods to design the complete magnetic system for a truly portable MRI scanner for neurological and musculoskeletal (MSK) applications, optimized for field homogeneity, field of view (FoV) and gradient performance compared to existing low-weight configurations. Approach. We explore optimal elliptic-bore Halbach configurations based on discrete arrays of permanent magnets. In…
▽ More
Objective. To develop methods to design the complete magnetic system for a truly portable MRI scanner for neurological and musculoskeletal (MSK) applications, optimized for field homogeneity, field of view (FoV) and gradient performance compared to existing low-weight configurations. Approach. We explore optimal elliptic-bore Halbach configurations based on discrete arrays of permanent magnets. In this way, we seek to improve the field homogeneity and remove constraints to the extent of the gradient coils typical of Halbach magnets. Specifically, we have optimized a tightly-packed distribution of magnetic Nd$_2$Fe$_14$B cubes with differential evolution algorithms, and a second array of shimming magnets with interior point and differential evolution methods. We have also designed and constructed an elliptical set of gradient coils that extend over the whole magnet length, maximizing the distance between the lobe centers. These are optimized with a target field method minimizing a cost function that considers also heat dissipation. Main result. We have employed the new toolbox to build the main magnet and gradient modules for a portable MRI scanner designed for point-of-care and residential use. The elliptical Halbach bore has semi-axes of 10 & 14 cm and the magnet generates a field of 87 mT homogeneous down to 5,700 ppm (parts per million) in a 20 cm diameter FoV, it weighs 216 kg and has a width of 65 cm and a height of 72 cm. Gradient efficiencies go up to around 0.8 mT/m/A, for a maximum of 12 mT/m with in 0.5 ms with 15 A & 15 V amplifier. The distance between lobes is 28 cm, significantly increased with respect to other Halbach-based scanners. Heat dissipation is around 25 W at maximum power, and gradient deviations from linearity are below 20% in a 20 cm sphere.
△ Less
Submitted 12 March, 2024;
originally announced March 2024.
-
Performance of a modular ton-scale pixel-readout liquid argon time projection chamber
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
B. Aimard,
F. Akbar,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
A. Alton,
R. Alvarez,
T. Alves,
H. Amar,
P. Amedo,
J. Anderson,
D. A. Andrade
, et al. (1340 additional authors not shown)
Abstract:
The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmi…
▽ More
The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmic ray events collected in the spring of 2021. We use this sample to demonstrate the imaging performance of the charge and light readout systems as well as the signal correlations between the two. We also report argon purity and detector uniformity measurements, and provide comparisons to detector simulations.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Geometric flavours of Quantum Field theory on a Cauchy hypersurface. Part II: Methods of quantization and evolution
Authors:
José Luis Alonso,
Carlos Bouthelier-Madre,
Jesús Clemente-Gallardo,
David Martínez-Crespo
Abstract:
In this series of papers we aim to provide a mathematically comprehensive framework to the Hamiltonian pictures of quantum field theory in curved spacetimes. Our final goal is to study the kinematics and the dynamics of the theory from the point of differential geometry in infinite dimensions. In this second part we use the tools of Gaussian analysis in infinite dimensional spaces introduced in th…
▽ More
In this series of papers we aim to provide a mathematically comprehensive framework to the Hamiltonian pictures of quantum field theory in curved spacetimes. Our final goal is to study the kinematics and the dynamics of the theory from the point of differential geometry in infinite dimensions. In this second part we use the tools of Gaussian analysis in infinite dimensional spaces introduced in the first part to describe rigorously the procedures of geometric quantization in the space of Cauchy data of a scalar theory. This leads us to discuss and establish relations between different pictures of QFT. We also apply these tools to describe the geometrization of the space of pure states of quantum field theory as a Khäler manifold. We use this to derive an evolution equation that preserves the geometric structure and avoids norm losses in the evolution. This leads us to a modification of the Schrödinger equation via a quantum connection that we discuss and exemplify in a simple case.
△ Less
Submitted 10 June, 2024; v1 submitted 10 February, 2024;
originally announced February 2024.
-
Memristor variability and stochastic physical properties modeling from a multivariate time series approach
Authors:
Francisco J. Alonso,
David Maldonado,
Ana M. Aguilera,
Juan B. Roldán
Abstract:
A powerful time series analysis modeling technique is presented to describe cycle-to-cycle variability in memristors. These devices show variability linked to the inherent stochasticity of device operation and it needs to be accurately modeled to build compact models for circuit simulation and design purposes. A new multivariate approach is proposed for the reset and set voltages that accurately d…
▽ More
A powerful time series analysis modeling technique is presented to describe cycle-to-cycle variability in memristors. These devices show variability linked to the inherent stochasticity of device operation and it needs to be accurately modeled to build compact models for circuit simulation and design purposes. A new multivariate approach is proposed for the reset and set voltages that accurately describes the statistical data structure of a resistive switching series. Experimental data were measured from advanced hafnium oxide based devices. The models reproduce the experiments correctly and a comparison of the multivariate and univariate approaches is shown for comparison.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Doping Liquid Argon with Xenon in ProtoDUNE Single-Phase: Effects on Scintillation Light
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
B. Aimard,
F. Akbar,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
A. Alton,
R. Alvarez,
H. Amar Es-sghir,
P. Amedo,
J. Anderson,
D. A. Andrade,
C. Andreopoulos
, et al. (1300 additional authors not shown)
Abstract:
Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUN…
▽ More
Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUNE-SP) at CERN, featuring 770 t of total liquid argon mass with 410 t of fiducial mass. The goal of the run was to measure the light and charge response of the detector to the addition of xenon, up to a concentration of 18.8 ppm. The main purpose was to test the possibility for reduction of non-uniformities in light collection, caused by deployment of photon detectors only within the anode planes. Light collection was analysed as a function of the xenon concentration, by using the pre-existing photon detection system (PDS) of ProtoDUNE-SP and an additional smaller set-up installed specifically for this run. In this paper we first summarize our current understanding of the argon-xenon energy transfer process and the impact of the presence of nitrogen in argon with and without xenon dopant. We then describe the key elements of ProtoDUNE-SP and the injection method deployed. Two dedicated photon detectors were able to collect the light produced by xenon and the total light. The ratio of these components was measured to be about 0.65 as 18.8 ppm of xenon were injected. We performed studies of the collection efficiency as a function of the distance between tracks and light detectors, demonstrating enhanced uniformity of response for the anode-mounted PDS. We also show that xenon doping can substantially recover light losses due to contamination of the liquid argon by nitrogen.
△ Less
Submitted 9 February, 2024; v1 submitted 2 February, 2024;
originally announced February 2024.
-
Multiobjective Optimization Analysis for Finding Infrastructure-as-Code Deployment Configurations
Authors:
Eneko Osaba,
Josu Diaz-de-Arcaya,
Juncal Alonso,
Jesus L. Lobo,
Gorka Benguria,
Iñaki Etxaniz
Abstract:
Multiobjective optimization is a hot topic in the artificial intelligence and operations research communities. The design and development of multiobjective methods is a frequent task for researchers and practitioners. As a result of this vibrant activity, a myriad of techniques have been proposed in the literature to date, demonstrating a significant effectiveness for dealing with situations comin…
▽ More
Multiobjective optimization is a hot topic in the artificial intelligence and operations research communities. The design and development of multiobjective methods is a frequent task for researchers and practitioners. As a result of this vibrant activity, a myriad of techniques have been proposed in the literature to date, demonstrating a significant effectiveness for dealing with situations coming from a wide range of real-world areas. This paper is focused on a multiobjective problem related to optimizing Infrastructure-as-Code deployment configurations. The system implemented for solving this problem has been coined as IaC Optimizer Platform (IOP). Despite the fact that a prototypical version of the IOP has been introduced in the literature before, a deeper analysis focused on the resolution of the problem is needed, in order to determine which is the most appropriate multiobjective method for embedding in the IOP. The main motivation behind the analysis conducted in this work is to enhance the IOP performance as much as possible. This is a crucial aspect of this system, deeming that it will be deployed in a real environment, as it is being developed as part of a H2020 European project. Going deeper, we resort in this paper to nine different evolutionary computation-based multiobjective algorithms. For assessing the quality of the considered solvers, 12 different problem instances have been generated based on real-world settings. Results obtained by each method after 10 independent runs have been compared using Friedman's non-parametric tests. Findings reached from the tests carried out lad to the creation of a multi-algorithm system, capable of applying different techniques according to the user's needs.
△ Less
Submitted 18 January, 2024;
originally announced January 2024.
-
MaRGA: a Graphical and Application Interface for MaRCoS
Authors:
J. M. Algarín,
T. Guallart-Naval,
J. Borreguero,
F. Galve,
J. Alonso
Abstract:
The open-source console MaRCoS, which stands for "Magnetic Resonance Control System", combines hardware, firmware and software elements for integral control of Magnetic Resonance Imaging (MRI) scanners. Previous developments have focused on making the system robust and reliable, rather than on users, who have been somewhat overlooked. This work describes a Graphical User Interface (GUI) designed f…
▽ More
The open-source console MaRCoS, which stands for "Magnetic Resonance Control System", combines hardware, firmware and software elements for integral control of Magnetic Resonance Imaging (MRI) scanners. Previous developments have focused on making the system robust and reliable, rather than on users, who have been somewhat overlooked. This work describes a Graphical User Interface (GUI) designed for intuitive control of MaRCoS, as well as compatibility with clinical environments. The GUI is based on an uncomplicated set of panels and a renewed Application Program Interface (API). Compared to the previous versions, the MaRGA package ("MaRCoS Graphical Application") includes new functionalities such as the possibility to export images to standard DICOM formats, create and manage clinical protocols, or display and process image reconstructions, among other features conceived to simplify the operation of MRI scanners. All prototypes in our facilities are commanded by MaRCoS and operated with the new GUI. Here we report on its performance on an experimental 0.2 T scanner designed for hard-tissue imaging, as well as a 72~mT portable scanner presently installed in the radiology department of a large hospital. The possibility to customize, adapt and streamline processes has substantially improved our workflows and overall experience.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
Pedestrian and Passenger Interaction with Autonomous Vehicles: Field Study in a Crosswalk Scenario
Authors:
Rubén Izquierdo,
Javier Alonso,
Ola Benderius,
Miguel Ángel Sotelo,
David Fernández Llorca
Abstract:
This study presents the outcomes of empirical investigations pertaining to human-vehicle interactions involving an autonomous vehicle equipped with both internal and external Human Machine Interfaces (HMIs) within a crosswalk scenario. The internal and external HMIs were integrated with implicit communication techniques, incorporating a combination of gentle and aggressive braking maneuvers within…
▽ More
This study presents the outcomes of empirical investigations pertaining to human-vehicle interactions involving an autonomous vehicle equipped with both internal and external Human Machine Interfaces (HMIs) within a crosswalk scenario. The internal and external HMIs were integrated with implicit communication techniques, incorporating a combination of gentle and aggressive braking maneuvers within the crosswalk. Data were collected through a combination of questionnaires and quantifiable metrics, including pedestrian decision to cross related to the vehicle distance and speed. The questionnaire responses reveal that pedestrians experience enhanced safety perceptions when the external HMI and gentle braking maneuvers are used in tandem. In contrast, the measured variables demonstrate that the external HMI proves effective when complemented by the gentle braking maneuver. Furthermore, the questionnaire results highlight that the internal HMI enhances passenger confidence only when paired with the aggressive braking maneuver.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
The DUNE Far Detector Vertical Drift Technology, Technical Design Report
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
B. Aimard,
F. Akbar,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
A. Alton,
R. Alvarez,
H. Amar,
P. Amedo,
J. Anderson,
D. A. Andrade,
C. Andreopoulos
, et al. (1304 additional authors not shown)
Abstract:
DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precisi…
▽ More
DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precision measurements of the PMNS matrix parameters, including the CP-violating phase. It will also stand ready to observe supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model.
The DUNE far detector implements liquid argon time-projection chamber (LArTPC) technology, and combines the many tens-of-kiloton fiducial mass necessary for rare event searches with the sub-centimeter spatial resolution required to image those events with high precision. The addition of a photon detection system enhances physics capabilities for all DUNE physics drivers and opens prospects for further physics explorations. Given its size, the far detector will be implemented as a set of modules, with LArTPC designs that differ from one another as newer technologies arise.
In the vertical drift LArTPC design, a horizontal cathode bisects the detector, creating two stacked drift volumes in which ionization charges drift towards anodes at either the top or bottom. The anodes are composed of perforated PCB layers with conductive strips, enabling reconstruction in 3D. Light-trap-style photon detection modules are placed both on the cryostat's side walls and on the central cathode where they are optically powered.
This Technical Design Report describes in detail the technical implementations of each subsystem of this LArTPC that, together with the other far detector modules and the near detector, will enable DUNE to achieve its physics goals.
△ Less
Submitted 5 December, 2023;
originally announced December 2023.
-
Optimizing IaC Configurations: a Case Study Using Nature-inspired Computing
Authors:
Eneko Osaba,
Gorka Benguria,
Jesus L. Lobo,
Josu Diaz-de-Arcaya,
Juncal Alonso,
Iñaki Etxaniz
Abstract:
In the last years, one of the fields of artificial intelligence that has been investigated the most is nature-inspired computing. The research done on this specific topic showcases the interest that sparks in researchers and practitioners, who put their focus on this paradigm because of the adaptability and ability of nature-inspired algorithms to reach high-quality outcomes on a wide range of pro…
▽ More
In the last years, one of the fields of artificial intelligence that has been investigated the most is nature-inspired computing. The research done on this specific topic showcases the interest that sparks in researchers and practitioners, who put their focus on this paradigm because of the adaptability and ability of nature-inspired algorithms to reach high-quality outcomes on a wide range of problems. In fact, this kind of methods has been successfully applied to solve real-world problems in heterogeneous fields such as medicine, transportation, industry, or software engineering. Our main objective with this paper is to describe a tool based on nature-inspired computing for solving a specific software engineering problem. The problem faced consists of optimizing Infrastructure as Code deployment configurations. For this reason, the name of the system is IaC Optimizer Platform. A prototypical version of the IOP was described in previous works, in which the functionality of this platform was introduced. With this paper, we take a step forward by describing the final release of the IOP, highlighting its main contribution regarding the current state-of-the-art, and justifying the decisions made on its implementation. Also, we contextualize the IOP within the complete platform in which it is embedded, describing how a user can benefit from its use. To do that, we also present and solve a real-world use case.
△ Less
Submitted 15 November, 2023;
originally announced November 2023.
-
A New Family of High-Current Cyclotrons for Isotope Production
Authors:
D. Winklehner,
J. R. Alonso,
J. M. Conrad
Abstract:
We are developing a high-current cyclotron as a driver for the IsoDAR neutrino experiment. It accelerates 5 mA H2+ to 60 MeV/amu, after which the electron is removed to produce a 10 mA, 60 MeV proton beam. The enabling innovations that offset space-charge effects occur at injection and in the first few turns, allowing one to construct cyclotrons with energies ranging from below 5 MeV up to 60 MeV/…
▽ More
We are developing a high-current cyclotron as a driver for the IsoDAR neutrino experiment. It accelerates 5 mA H2+ to 60 MeV/amu, after which the electron is removed to produce a 10 mA, 60 MeV proton beam. The enabling innovations that offset space-charge effects occur at injection and in the first few turns, allowing one to construct cyclotrons with energies ranging from below 5 MeV up to 60 MeV/amu, or possibly higher, with the same performance for accelerated ions with Q/A = 0.5 (H2+, D+, He++, ...). In this paper, we discuss the possible uses of such cyclotrons for isotope production, including production of long-lived generator parents (68Ga, 44Ti, 82Sr,...), as well as intense fast neutron beams from deuteron breakup for (n,2n) production of isotopes like 225Ac.
△ Less
Submitted 26 April, 2024; v1 submitted 29 October, 2023;
originally announced October 2023.
-
Delaunay Bifiltrations of Functions on Point Clouds
Authors:
Ángel Javier Alonso,
Michael Kerber,
Tung Lam,
Michael Lesnick
Abstract:
The Delaunay filtration $\mathcal{D}_{\bullet}(X)$ of a point cloud $X\subset \mathbb{R}^d$ is a central tool of computational topology. Its use is justified by the topological equivalence of $\mathcal{D}_{\bullet}(X)$ and the offset (i.e., union-of-balls) filtration of $X$. Given a function $γ: X \to \mathbb{R}$, we introduce a Delaunay bifiltration $\mathcal{DC}_{\bullet}(γ)$ that satisfies an a…
▽ More
The Delaunay filtration $\mathcal{D}_{\bullet}(X)$ of a point cloud $X\subset \mathbb{R}^d$ is a central tool of computational topology. Its use is justified by the topological equivalence of $\mathcal{D}_{\bullet}(X)$ and the offset (i.e., union-of-balls) filtration of $X$. Given a function $γ: X \to \mathbb{R}$, we introduce a Delaunay bifiltration $\mathcal{DC}_{\bullet}(γ)$ that satisfies an analogous topological equivalence, ensuring that $\mathcal{DC}_{\bullet}(γ)$ topologically encodes the offset filtrations of all sublevel sets of $γ$, as well as the topological relations between them. $\mathcal{DC}_{\bullet}(γ)$ is of size $O(|X|^{\lceil\frac{d+1}{2}\rceil})$, which for $d$ odd matches the worst-case size of $\mathcal{D}_{\bullet}(X)$. Adapting the Bowyer-Watson algorithm for computing Delaunay triangulations, we give a simple, practical algorithm to compute $\mathcal{DC}_{\bullet}(γ)$ in time $O(|X|^{\lceil \frac{d}{2}\rceil +1})$. Our implementation, based on CGAL, computes $\mathcal{DC}_{\bullet}(γ)$ with modest overhead compared to computing $\mathcal{D}_{\bullet}(X)$, and handles tens of thousands of points in $\mathbb{R}^3$ within seconds.
△ Less
Submitted 24 October, 2023;
originally announced October 2023.
-
The twisting index in semitoric systems
Authors:
Jaume Alonso,
Sonja Hohloch,
Joseph Palmer
Abstract:
Semitoric integrable systems were symplectically classified by Pelayo and Vu Ngoc in 2009-2011 in terms of five invariants. Four of these invariants were already well-understood prior to the classification, but the fifth invariant, the so-called twisting index invariant, came as a surprise. Intuitively, the twisting index encodes how the structure in a neighborhood of a focus-focus fiber compares…
▽ More
Semitoric integrable systems were symplectically classified by Pelayo and Vu Ngoc in 2009-2011 in terms of five invariants. Four of these invariants were already well-understood prior to the classification, but the fifth invariant, the so-called twisting index invariant, came as a surprise. Intuitively, the twisting index encodes how the structure in a neighborhood of a focus-focus fiber compares to the large-scale structure of the semitoric system and it was originally defined by comparing certain momentum maps.
In the first half of the present paper, we produce several new formulations of the twisting index which give rise to dynamical, geometric, and topological interpretations. More specifically, we describe it in terms of differences of action variables, Taylor series, and homology cycles. In the second half of the paper, we compute the twisting index invariant of a specific family of systems with two focus-focus singular points (the so-called generalized coupled angular momenta), which is the first time that the twisting index has been computed for a system with more than one focus-focus point. Moreover, we also compute the terms of the Taylor series invariant up to second order. Since the other invariants of this family were already computed, this becomes the third family of semitoric systems for which all invariants are known, after the coupled spin oscillators and the coupled angular momenta.
△ Less
Submitted 8 January, 2024; v1 submitted 28 September, 2023;
originally announced September 2023.
-
Physics Potential of a Few Kiloton Scale Neutrino Detector at a Deep Underground Lab in Korea
Authors:
Seon-Hee Seo,
Jose Alonso,
Pouya Bakhti,
Janet Conrad,
Steve Dye,
Doojin Kim,
Jost Migenda,
Marco Pallavicini,
Jong-Chul Park,
Meshkat Rajaee,
Mike Shaevitz,
Seodong Shin,
Joshua Spitz,
Daniel Winklehner,
Slawomir Wronka,
Michael Wurm,
Minfang Yeh
Abstract:
The demand for underground labs for neutrino and rare event search experiments has been increasing over the last few decades. Yemilab, constructed in October 2022, is the first deep ($\sim$1~km) underground lab dedicated to science in Korea, where a large cylindrical cavern (D: 20~m, H: 20~m) was excavated in addition to the main caverns and halls. The large cavern could be utilized for a low back…
▽ More
The demand for underground labs for neutrino and rare event search experiments has been increasing over the last few decades. Yemilab, constructed in October 2022, is the first deep ($\sim$1~km) underground lab dedicated to science in Korea, where a large cylindrical cavern (D: 20~m, H: 20~m) was excavated in addition to the main caverns and halls. The large cavern could be utilized for a low background neutrino experiment by a liquid scintillator-based detector (LSC) where a 2.26 kiloton LS target would be filled. It's timely to have such a large but ultra-pure LS detector after the shutdown of the Borexino experiment so that solar neutrinos can be measured much more precisely. Interesting BSM physics searches can be also pursued with this detector when it's combined with an electron linac, a proton cyclotron (IsoDAR source), or a radioactive source. This article discusses the concept of a candidate detector and the physics potential of a large liquid scintillator detector.
△ Less
Submitted 23 September, 2023;
originally announced September 2023.
-
Effective nonlinear Ehrenfest hybrid quantum-classical dynamics
Authors:
J. L. Alonso,
C. Bouthelier-Madre,
J. Clemente-Gallardo,
D. Martínez-Crespo,
J. Pomar
Abstract:
The definition of a consistent evolution equation for statistical hybrid quantum-classical systems is still an open problem. In this paper we analyze the case of Ehrenfest dynamics on systems defined by a probability density and identify the relations of the non-linearity of the dynamics with the obstructions to define a consistent dynamics for the first quantum moment of the distribution. This fi…
▽ More
The definition of a consistent evolution equation for statistical hybrid quantum-classical systems is still an open problem. In this paper we analyze the case of Ehrenfest dynamics on systems defined by a probability density and identify the relations of the non-linearity of the dynamics with the obstructions to define a consistent dynamics for the first quantum moment of the distribution. This first quantum moment represents the physical states as a family of classically-parametrized density matrices $\hat ρ(ξ)$, for $ξ$ a classical point; and it is the most common representation of hybrid systems in the literature. Due to this obstruction, we consider higher order quantum moments, and argue that only a finite number of them are physically measurable. Because of this, we propose an effective solution for the hybrid dynamics problem based on approximating the distribution by those moments and representing the states by them.
△ Less
Submitted 29 August, 2023; v1 submitted 28 August, 2023;
originally announced August 2023.
-
The Landau equation with moderate soft potentials: An approach using $\varepsilon$-Poincaré inequality and Lorentz spaces
Authors:
R. J. Alonso,
V. Bagland,
B. Lods
Abstract:
This document presents an elementary approach using $\varepsilon$-Poincaré inequality to prove generation of $L^{p}$-bounds, $p\in(1,\infty)$, for the homogeneous Landau equation with moderate soft potentials $γ\in[-2,0)$. The critical case $γ=-2$ uses an interpolation approach in the realm of Lorentz spaces and entropy. Alternatively, a direct approach using the Hardy-Littlewood-Sobolev (HLS) ine…
▽ More
This document presents an elementary approach using $\varepsilon$-Poincaré inequality to prove generation of $L^{p}$-bounds, $p\in(1,\infty)$, for the homogeneous Landau equation with moderate soft potentials $γ\in[-2,0)$. The critical case $γ=-2$ uses an interpolation approach in the realm of Lorentz spaces and entropy. Alternatively, a direct approach using the Hardy-Littlewood-Sobolev (HLS) inequality and entropy is also presented. On this basis, the generation of pointwise bounds $p=\infty$ is deduced from a De Giorgi argument.
△ Less
Submitted 8 August, 2023;
originally announced August 2023.
-
Human-vehicle interaction for autonomous vehicles in crosswalk scenarios: Field experiments with pedestrians and passengers
Authors:
R. Izquierdo,
S. Martín,
J. Alonso,
I. Parra,
M. A.,
D. Fernández-Llorca
Abstract:
This paper presents the results of real-world testing of human-vehicle interactions with an autonomous vehicle equipped with internal and external Human Machine Interfaces (HMIs) in a crosswalk scenario. The internal and external HMIs were combined with implicit communication techniques using gentle and aggressive braking maneuvers in the crosswalk. Results have been collected in the form of quest…
▽ More
This paper presents the results of real-world testing of human-vehicle interactions with an autonomous vehicle equipped with internal and external Human Machine Interfaces (HMIs) in a crosswalk scenario. The internal and external HMIs were combined with implicit communication techniques using gentle and aggressive braking maneuvers in the crosswalk. Results have been collected in the form of questionnaires and measurable variables such as distance or speed when the pedestrian decides to cross. The questionnaires show that pedestrians feel safer when the external HMI or the gentle braking maneuver is used interchangeably, while the measured variables show that external HMI only helps in combination with the gentle braking maneuver. The questionnaires also show that internal HMI only improves passenger confidence in combination with the aggressive braking maneuver.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
First glycine isomer detected in the interstellar medium: glycolamide (NH$_2$C(O)CH$_2$OH)
Authors:
Víctor M. Rivilla,
Miguel Sanz-Novo,
Izaskun Jiménez-Serra,
Jesús Martín-Pintado,
Laura Colzi,
Shaoshan Zeng,
Andrés Megías,
Álvaro López-Gallifa,
Antonio Martínez-Henares,
Sarah Massalkhi,
Belén Tercero,
Pablo de Vicente,
Sergio Martín,
David San Andrés,
Miguel A. Requena-Torres,
José Luis Alonso
Abstract:
We report the first detection in the interstellar medium of a C$_2$H$_5$O$_2$N isomer: $syn$-glycolamide (NH$_2$C(O)CH$_2$OH). The exquisite sensitivity at sub-mK levels of an ultra-deep spectral survey carried out with the Yebes 40m and IRAM 30m telescopes towards the G+0.693-0.027 molecular cloud have allowed us to unambiguously identify multiple transitions of this species. We derived a column…
▽ More
We report the first detection in the interstellar medium of a C$_2$H$_5$O$_2$N isomer: $syn$-glycolamide (NH$_2$C(O)CH$_2$OH). The exquisite sensitivity at sub-mK levels of an ultra-deep spectral survey carried out with the Yebes 40m and IRAM 30m telescopes towards the G+0.693-0.027 molecular cloud have allowed us to unambiguously identify multiple transitions of this species. We derived a column density of (7.4 $\pm$ 0.7)$\times$10$^{12}$ cm$^{-2}$, which implies a molecular abundance with respect to H$_2$ of 5.5$\times$10$^{-11}$. The other C$_2$H$_5$O$_2$N isomers, including the higher-energy $anti$ conformer of glycolamide, and two conformers of glycine, were not detected. The upper limit derived for the abundance of glycine indicates that this amino acid is surely less abundant than its isomer glycolamide in the ISM. The abundances of the C$_2$H$_5$O$_2$N isomers cannot be explained in terms of thermodynamic equilibrium, and thus chemical kinetics need to be invoked. While the low abundance of glycine might not be surprising, based on the relative low abundances of acids in the ISM compared to other compounds (e.g. alcohols, aldehydes or amines), several chemical pathways can favour the formation of its isomer glycolamide. It can be formed through radical-radical reactions on the surface of dust grains. The abundances of these radicals can be significantly boosted in an environment affected by a strong ultraviolet field induced by cosmic rays, such as that expected in G+0.693-0.027. Therefore, as shown by several recent molecular detections towards this molecular cloud, it stands out as the best target to discover new species with carbon, oxygen and nitrogen with increasing chemical complexity.
△ Less
Submitted 11 August, 2023; v1 submitted 21 July, 2023;
originally announced July 2023.
-
Role of positional disorder in fully textured ensembles of Ising-like dipoles
Authors:
J. J. Alonso,
B. Alles,
J. G. Malherbe,
V. Russier
Abstract:
We study by numerical simulation the magnetic order in ensembles of randomly packed magnetic spherical particles which, induced by their uniaxial anisotropy in the strong coupling limit, behave as Ising dipoles. We explore the role of the frozen disorder in the positions of the particles assuming a common fixed direction for the easy axes of all spheres. We look at two types of spatially disordere…
▽ More
We study by numerical simulation the magnetic order in ensembles of randomly packed magnetic spherical particles which, induced by their uniaxial anisotropy in the strong coupling limit, behave as Ising dipoles. We explore the role of the frozen disorder in the positions of the particles assuming a common fixed direction for the easy axes of all spheres. We look at two types of spatially disordered configurations. In the first one we consider isotropic positional distributions which can be obtained from the liquid state of the hard sphere fluid. We derive the phase diagram in the $T$-$Φ$ plane where $T$ is the temperature and $Φ$ the volume fraction. This diagram exhibits long-range ferromagnetic order at low $T$ for volume fractions above the threshold $Φ_{c} = 0.157$ predicted by mean-field calculations. For $Φ\lesssim Φ_{c}$ a spin-glass phase forms with the same marginal behavior found for other strongly disordered dipolar systems. The second type of spatial configurations we study are anisotropic distributions that can be obtained by freezing a dipolar hard sphere liquid in its polarized state at low temperature. This structural anisotropy enhances the ferromagnetic order present in isotropic hard sphere configurations.
△ Less
Submitted 4 October, 2023; v1 submitted 19 July, 2023;
originally announced July 2023.
-
Dynamical degrees of birational maps from indices of polynomials with respect to blow-ups II. 3D examples
Authors:
Jaume Alonso,
Yuri B. Suris,
Kangning Wei
Abstract:
The goal of this paper is the exact computation of the degrees $\text{deg}(f^n)$ of the iterates of birational maps $f: \mathbb{P}^N \dashrightarrow \mathbb{P}^N$. In the preceding companion paper, a new method has been proposed based on the use of indices of polynomials associated to the local blow-ups used to resolve contractions of hypersurfaces by $f$, and on the control of the factorization o…
▽ More
The goal of this paper is the exact computation of the degrees $\text{deg}(f^n)$ of the iterates of birational maps $f: \mathbb{P}^N \dashrightarrow \mathbb{P}^N$. In the preceding companion paper, a new method has been proposed based on the use of indices of polynomials associated to the local blow-ups used to resolve contractions of hypersurfaces by $f$, and on the control of the factorization of pull-backs of polynomials. This leads to recurrence relations for the degrees and the indices. We apply this method to several illustrative examples in three dimensions. These examples demonstrate the flexibility of the method which, in particular, does not require the construction of an algebraically stable lift of $f$, unlike the previously known methods based on the Picard group.
△ Less
Submitted 21 July, 2023; v1 submitted 19 July, 2023;
originally announced July 2023.
-
Hybrid Geometrodynamics: A Hamiltonian description of classical gravity coupled to quantum matter
Authors:
J. L. Alonso,
C. Bouthelier-Madre,
J. Clemente-Gallardo,
D. Martínez-Crespo
Abstract:
We generalize the Hamiltonian picture of General Relativity coupled to classical matter, known as geometrodynamics, to the case where such matter is described by a Quantum Field Theory in Curved Spacetime, but gravity is still described by a classical metric tensor field over a spatial hypersurface and its associated momentum. Thus, in our approach there is no non-dynamic background structure, apa…
▽ More
We generalize the Hamiltonian picture of General Relativity coupled to classical matter, known as geometrodynamics, to the case where such matter is described by a Quantum Field Theory in Curved Spacetime, but gravity is still described by a classical metric tensor field over a spatial hypersurface and its associated momentum. Thus, in our approach there is no non-dynamic background structure, apart from the manifold of events, and the gravitational and quantum degrees of freedom have their dynamics inextricably coupled. Given the Hamiltonian natureof the framework, we work with the generators of hypersurface deformations over the manifold of quantum states. The construction relies heavily on the differential geometry of a fibration of the set of quantum states over the set of gravitational variables. An important feature of this work is the use of Gaussian measures over the space of matter fields and of Hida distributions to define a common superspace to all possible Hilbert spaces with different measures, to properly characterize the Schrodinger wave functional picture of QFT in curved spacetime. This allows us to relate states within different Hilbert spaces in the case of vacuum states or measures that depend on the gravitational degrees of freedom, as the ones associated to Ashtekar's complex structure. This is achieved through the inclusion of a quantum Hermitian connection for the fibration, which will have profound physical implications. The most remarkable physical features of the construction are norm conservation of the quantum state (even if the total dynamics are non-unitary), the clear identification of the hybrid conserved quantities and the description of a dynamical backreaction of quantum matter on geometry and vice versa, which shall modify the physical properties the gravitational field would have in the absence of backreaction.
△ Less
Submitted 14 February, 2024; v1 submitted 3 July, 2023;
originally announced July 2023.
-
Geometric flavours of Quantum Field theory on a Cauchy hypersurface. Part I: Gaussian analysis and other Mathematical aspects
Authors:
José Luis Alonso,
Carlos Bouthelier-Madre,
Jesús Clemente-Gallardo,
David Martínez-Crespo
Abstract:
In this series of papers we aim to provide a mathematically comprehensive framework to the Hamiltonian pictures of quantum field theory in curved spacetimes. Our final goal is to study the kinematics and the dynamics of the theory from the point of differential geometry in infinite dimensions.In this first part we introduce the tools of Gaussian analysis in infinite dimensional spaces of distribut…
▽ More
In this series of papers we aim to provide a mathematically comprehensive framework to the Hamiltonian pictures of quantum field theory in curved spacetimes. Our final goal is to study the kinematics and the dynamics of the theory from the point of differential geometry in infinite dimensions.In this first part we introduce the tools of Gaussian analysis in infinite dimensional spaces of distributions. These spaces will serve the basis to understand the Schrödinger and Holomorphic pictures, over arbitrary Cauchy hypersurfaces, using tools of Hida-Malliavin calculus. Here the Wiener-Ito decomposition theorem provides the QFT particle interpretation. Special emphasis is done in the applications to quantization of these tools in the second part of this paper. We devote a section to introduce Hida test functions as a notion of second quantized test functions. We also analyze of the ingredients of classical field theory modeled as distributions paving the way for quantization procedures that will be analyzed in the second part of this series.
△ Less
Submitted 10 June, 2024; v1 submitted 26 June, 2023;
originally announced June 2023.
-
The Cauchy problem for Boltzmann bi-linear systems: The mixing of monatomic and polyatomic gases
Authors:
Ricardo J. Alonso,
Irene M. Gamba,
Milana Pavic-Colic
Abstract:
From a unified vision of vector valued solutions in weighted Banach spaces, this manuscript establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such mod…
▽ More
From a unified vision of vector valued solutions in weighted Banach spaces, this manuscript establishes the existence and uniqueness for space homogeneous Boltzmann bi-linear systems with conservative collisional forms arising in complex gas dynamical structures. This broader vision is directly applied to dilute multi-component gas mixtures composed of both monatomic and polyatomic gases. Such models can be viewed as extensions of scalar Boltzmann binary elastic flows, as much as monatomic gas mixtures with disparate masses and single polyatomic gases, providing a unify approach for vector valued solutions in weighted vector Banach spaces.
Novel aspects of this work include developing the extension of a general ODE theory in vector valued weighted Banach spaces, precise lower bounds for the collision frequency in terms of the weighted Banach norm, energy identities and the consequently angular or compact manifold averaging lemmas which provide coerciveness resulting into global in time stability, a new combinatorics estimate for $p$-binomial forms producing sharper estimates for the $k$-moments of bi-linear collisional forms. These techniques enable the Cauchy problem improvement that resolves the model with initial data corresponding to strictly positive and bounded initial vector valued mass and total energy, in addition to only $2^+$ moment determined by the hard potential rates discrepancy, a result comparable in generality to the classical Cauchy theory of the scalar homogeneous Boltzmann equation.
△ Less
Submitted 12 April, 2023;
originally announced April 2023.
-
Impact of cross-section uncertainties on supernova neutrino spectral parameter fitting in the Deep Underground Neutrino Experiment
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
Z. Ahmad,
J. Ahmed,
B. Aimard,
F. Akbar,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
A. Alton,
R. Alvarez,
P. Amedo,
J. Anderson,
D. A. Andrade
, et al. (1294 additional authors not shown)
Abstract:
A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the $\mathcal{O}(10)$ MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the $ν_e$ component of the supernova flux, enabling a wide variety of physics…
▽ More
A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the $\mathcal{O}(10)$ MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the $ν_e$ component of the supernova flux, enabling a wide variety of physics and astrophysics measurements. A key requirement for a correct interpretation of these measurements is a good understanding of the energy-dependent total cross section $σ(E_ν)$ for charged-current $ν_e$ absorption on argon. In the context of a simulated extraction of supernova $ν_e$ spectral parameters from a toy analysis, we investigate the impact of $σ(E_ν)$ modeling uncertainties on DUNE's supernova neutrino physics sensitivity for the first time. We find that the currently large theoretical uncertainties on $σ(E_ν)$ must be substantially reduced before the $ν_e$ flux parameters can be extracted reliably: in the absence of external constraints, a measurement of the integrated neutrino luminosity with less than 10\% bias with DUNE requires $σ(E_ν)$ to be known to about 5%. The neutrino spectral shape parameters can be known to better than 10% for a 20% uncertainty on the cross-section scale, although they will be sensitive to uncertainties on the shape of $σ(E_ν)$. A direct measurement of low-energy $ν_e$-argon scattering would be invaluable for improving the theoretical precision to the needed level.
△ Less
Submitted 7 July, 2023; v1 submitted 29 March, 2023;
originally announced March 2023.
-
Dynamical degrees of birational maps from indices of polynomials with respect to blow-ups I. General theory and 2D examples
Authors:
Jaume Alonso,
Yuri B. Suris,
Kangning Wei
Abstract:
In this paper we address the problem of computing $\text{deg}(f^n)$, the degrees of iterates of a birational map $f:\mathbb{P}^N\rightarrow\mathbb{P}^N$. For this goal, we develop a method based on two main ingredients: the factorization of a polynomial under pull-back of $f$, based on local indices of a polynomial associated to blow-ups used to resolve the contraction of hypersurfaces by $f$, and…
▽ More
In this paper we address the problem of computing $\text{deg}(f^n)$, the degrees of iterates of a birational map $f:\mathbb{P}^N\rightarrow\mathbb{P}^N$. For this goal, we develop a method based on two main ingredients: the factorization of a polynomial under pull-back of $f$, based on local indices of a polynomial associated to blow-ups used to resolve the contraction of hypersurfaces by $f$, and the propagation of these indices along orbits of $f$.
For maps admitting algebraically stable modifications $f_X:X\rightarrow X$, where $X$ is a variety obtained from $\mathbb P^N$ by a finite number of blow-ups, this method leads to an algorithm producing a finite system of recurrent equations relating the degrees and indices of iterated pull-backs of linear polynomials. We illustrate the method by three representative two-dimensional examples. It is actually applicable in any dimension, and we will provide a number of three-dimensional examples as a separate companion paper.
△ Less
Submitted 30 March, 2023; v1 submitted 28 March, 2023;
originally announced March 2023.
-
Portable MRI for major sporting events -- a case study on the MotoGP World Championship
Authors:
J. M. Algarín,
T. Guallart-Naval,
E. Gastaldi-Orquín,
R. Bosch,
F. J. Lloris,
E. Pallás,
J. P. Rigla,
P. Martínez,
J. Borreguero,
R. Alamar,
L. Martí-Bonmatí,
J. M. Benlloch,
F. Galve,
J. Alonso
Abstract:
The goal of this work is to showcase the clinical value that portable MRI can provide in crowded events and major sports competitions. We temporarily installed a low-field and low-cost portable MRI system for extremity imaging in the medical facilities of the Ricardo Tormo Motor Racing Circuit during the four days of the Motorcycle Grand Prix held in Valencia (Spain), which closed the 2022 season…
▽ More
The goal of this work is to showcase the clinical value that portable MRI can provide in crowded events and major sports competitions. We temporarily installed a low-field and low-cost portable MRI system for extremity imaging in the medical facilities of the Ricardo Tormo Motor Racing Circuit during the four days of the Motorcycle Grand Prix held in Valencia (Spain), which closed the 2022 season of the MotoGP. During this time, we scanned 14 subjects, running a total of 21 protocols for wrist, knee and ankle imaging. Each protocol included a minimum of one T1-weighted 3D-RARE sequence for general anatomical information, and one 3D-STIR sequence to highlight fluid accumulation and inflammation. The circuit medical staff were able to visualize a number of lesions and conditions in the low-field reconstructions, including gonarthrosis, effusion, or Haglund's syndrome, as well as metallic implants and tissue changes due to surgical interventions. Out of eight low-field acquisitions on previously diagnosed lesions, only two (a meniscus tear and a Baker cyst) were not detected by the experts that evaluated our images. The main highlight was that a low-field MRI scan on a subject reporting pain in a wrist revealed a traumatic arthritis which an X-ray radiograph and visual inspection had missed. We have operated in a scenario where high-field MRI is unlikely to play a role but where a low-field system can lead to improved medical attention. In the case reported here, system transport, installation in the circuit facilities and calibration were all uncomplicated. The images presented to the medical staff were mostly unprocessed and there is thus room for improvement. In conclusion, this work supports the claim that low-field MRI can likely provide added value whenever concepts such as accessibility, portability and low-cost outweigh exquisite detail in images.
△ Less
Submitted 17 March, 2023; v1 submitted 16 March, 2023;
originally announced March 2023.
-
A numerically stable communication-avoiding s-step GMRES algorithm
Authors:
Zan Xu,
Juan J. Alonso,
Eric Darve
Abstract:
Krylov subspace methods are extensively used in scientific computing to solve large-scale linear systems. However, the performance of these iterative Krylov solvers on modern supercomputers is limited by expensive communication costs. The $s$-step strategy generates a series of $s$ Krylov vectors at a time to avoid communication. Asymptotically, the $s$-step approach can reduce communication laten…
▽ More
Krylov subspace methods are extensively used in scientific computing to solve large-scale linear systems. However, the performance of these iterative Krylov solvers on modern supercomputers is limited by expensive communication costs. The $s$-step strategy generates a series of $s$ Krylov vectors at a time to avoid communication. Asymptotically, the $s$-step approach can reduce communication latency by a factor of $s$. Unfortunately, due to finite-precision implementation, the step size has to be kept small for stability. In this work, we tackle the numerical instabilities encountered in the $s$-step GMRES algorithm. By choosing an appropriate polynomial basis and block orthogonalization schemes, we construct a communication avoiding $s$-step GMRES algorithm that automatically selects the optimal step size to ensure numerical stability. To further maximize communication savings, we introduce scaled Newton polynomials that can increase the step size $s$ to a few hundreds for many problems. An initial step size estimator is also developed to efficiently choose the optimal step size for stability. The guaranteed stability of the proposed algorithm is demonstrated using numerical experiments. In the process, we also evaluate how the choice of polynomial and preconditioning affects the stability limit of the algorithm. Finally, we show parallel scalability on more than 14,000 cores in a distributed-memory setting. Perfectly linear scaling has been observed in both strong and weak scaling studies with negligible communication costs.
△ Less
Submitted 17 March, 2023; v1 submitted 15 March, 2023;
originally announced March 2023.
-
Decomposition of zero-dimensional persistence modules via rooted subsets
Authors:
Ángel Javier Alonso,
Michael Kerber
Abstract:
We study the decomposition of zero-dimensional persistence modules, viewed as functors valued in the category of vector spaces factorizing through sets. Instead of working directly at the level of vector spaces, we take a step back and first study the decomposition problem at the level of sets.
This approach allows us to define the combinatorial notion of rooted subsets. In the case of a filtere…
▽ More
We study the decomposition of zero-dimensional persistence modules, viewed as functors valued in the category of vector spaces factorizing through sets. Instead of working directly at the level of vector spaces, we take a step back and first study the decomposition problem at the level of sets.
This approach allows us to define the combinatorial notion of rooted subsets. In the case of a filtered metric space $M$, rooted subsets relate the clustering behavior of the points of $M$ with the decomposition of the associated persistence module. In particular, we can identify intervals in such a decomposition quickly. In addition, rooted subsets can be understood as a generalization of the elder rule, and are also related to the notion of constant conqueror of Cai, Kim, Mémoli and Wang. As an application, we give a lower bound on the number of intervals that we can expect in the decomposition of zero-dimensional persistence modules of a density-Rips filtration in Euclidean space: in the limit, and under very general circumstances, we can expect that at least 25% of the indecomposable summands are interval modules.
△ Less
Submitted 10 March, 2023;
originally announced March 2023.
-
Low temperature Terahertz Spectroscopy of LaFeO$_3$, PrFeO$_3$, ErFeO$_3$, and LuFeO$_3$: Quasimagnon resonances and ground multiplet transitions
Authors:
Nestor E. Massa,
Leire del Campo,
Vinh Ta Puoc,
Paula Kaiser,
Jose Antonio Alonso
Abstract:
We report on zone center THz excitations of non-Jahn Teller LaFeO$_3$, PrFeO$_3$, ErFeO$_3$, and LuFeO$_3$ distorted perovskites under external magnetic fields up 7 T. Low temperature-low energy absorptions of LaFeO$_3$ show antiferromagnetic and ferromagnetic quasimagnons at $ω$AFM ~31.4 and $ω$FM ~26.7 cm$^{-1}$ in the $Γ$4 (Gx, Ay, Fz) representation. LuFeO$_3$ is characterized by zero field ma…
▽ More
We report on zone center THz excitations of non-Jahn Teller LaFeO$_3$, PrFeO$_3$, ErFeO$_3$, and LuFeO$_3$ distorted perovskites under external magnetic fields up 7 T. Low temperature-low energy absorptions of LaFeO$_3$ show antiferromagnetic and ferromagnetic quasimagnons at $ω$AFM ~31.4 and $ω$FM ~26.7 cm$^{-1}$ in the $Γ$4 (Gx, Ay, Fz) representation. LuFeO$_3$ is characterized by zero field magnetic resonances at $ω$AFM ~26.3 cm$^{-1}$ and $ω$FM ~22.4 cm$^{-1}$ in addition to Fe$^{3+}$ Zeeman-split crystal field (CF) 6A$_1$ ground transitions at ~10.4 cm$^{-1}$ triggered by structural deviations induced by smaller Lu 4f$^{14}$. This local non-centrosymmetric departure is also found in ErFeO$_3$ (Kramers 4f$^{11}$ Er$^{3+}$ (4I15/2); Γ2 (Fx, Cy, Gz) <TSR ~93 K), but with the ~4 cm$^{-1}$ Fe$^{3+}$ Zeeman branching strongly biased toward higher energies. Quasimagnons at $ω$AFM ~31.5 cm$^{-1}$ and $ω$FM ~21.5 cm$^{-1}$ in ErFeO$_3$ do not undergo field induced band splits but a 13-fold increase in the antiferro ($ω$AMF) /ferro($ω$AFM) intensity ratio. There is a remarkable field-dependent CF matching population balance between Fe$^{3+}$ higher and Er$^{3+}$ lower Zeeman branches. Antiferro- and ferro- resonances in PrFeO$_3$ turn much broader as non-Kramers Pr$^3$ introduces ligand changes at the A site leading into near degeneracy the antiferromagnetic mode and the lowest Pr$^{3+}$ CF transition. We conclude that low energy excitations in RFeO$_3$ (R=rare earth) strongly depend on the lanthanide ionic size. Minute lattice displacements also underlie considering non-centrosymmetric the most distorted RFeO$_3$ (R=rare earth). Changes triggered by the smaller rare earth and the nonlinear intrinsic oxygen ion polarizability provide grounds for interplay of ionic and electronic interactions yielding ferroelectric spontaneous polarization.
△ Less
Submitted 10 September, 2023; v1 submitted 5 March, 2023;
originally announced March 2023.
-
Data augmentation for battery materials using lattice scaling
Authors:
Eduardo Abenza,
César Alonso,
Isabel Sobrados,
José M. Amarilla,
Javier L. Rodríguez,
José A. Alonso,
Roberto G. E. Martín,
Maria C. Asensio
Abstract:
A significant step forward in Lithium-ion batteries (LIBs) developments can only be achieved by proposing mold-breaking research based on selecting the best materials for the cell components, optimizing cell manufacture, anticipating the degradation mechanisms of the LIBs, and consolidating the regeneration processes of damaged batteries. LIBs with longer recycling life, better safety, and the abi…
▽ More
A significant step forward in Lithium-ion batteries (LIBs) developments can only be achieved by proposing mold-breaking research based on selecting the best materials for the cell components, optimizing cell manufacture, anticipating the degradation mechanisms of the LIBs, and consolidating the regeneration processes of damaged batteries. LIBs with longer recycling life, better safety, and the ability to be reused will establish sustainable state-of-the-art batteries with maximum energy efficiency, low costs, and minimal CO2 emissions within a circular economy, promoting sustainability in areas as relevant as electromobility and portable electronics. Recently, there has been increasing interest in applying Artificial Intelligence (AI) techniques and their subclasses to better predict novel materials with designed properties. This collection of methods has already obtained considerable success, having been used to predict numerous physical properties of materials. However, compared to other fields, the materials data are typically much smaller and sometimes more diverse, which undoubtedly affects the construction and effectiveness of AI models. At present, several Data Augmentation (DA) methods have been proposed in materials science based on flipping, rotating, and distorting the unit cells of materials, which have been demonstrated to be very efficient in increasing the size and quality of data. Here we present an even more effective new method of Data Augmentation based on the lattice scaling of crystal structures. In the lattice scaling DA method, the unit cell is perturbed, undergoing an increase and decrease of the unit cell volume in an isotropic or anisotropic way. This transformation is particularly pertinent for battery components since volume changes of up to 5% have been reported for the insertion-based LIBs during cycling.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
Highly-parallelized simulation of a pixelated LArTPC on a GPU
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
Z. Ahmad,
J. Ahmed,
B. Aimard,
F. Akbar,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
C. Alt,
A. Alton,
R. Alvarez,
P. Amedo,
J. Anderson
, et al. (1282 additional authors not shown)
Abstract:
The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we pr…
▽ More
The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we present the first implementation of a full microphysical simulator of a liquid argon time projection chamber (LArTPC) equipped with light readout and pixelated charge readout, developed for the DUNE Near Detector. The software is implemented with an end-to-end set of GPU-optimized algorithms. The algorithms have been written in Python and translated into CUDA kernels using Numba, a just-in-time compiler for a subset of Python and NumPy instructions. The GPU implementation achieves a speed up of four orders of magnitude compared with the equivalent CPU version. The simulation of the current induced on $10^3$ pixels takes around 1 ms on the GPU, compared with approximately 10 s on the CPU. The results of the simulation are compared against data from a pixel-readout LArTPC prototype.
△ Less
Submitted 28 February, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Ferromagnetic frozen structures from the dipolar hard spheres fluid at moderate and small volume fractions
Authors:
Jean-Guillaume Malherbe,
Vincent Russier,
Juan-Jose Alonso
Abstract:
We study the magnetic phase diagram of an ensemble of dipolar hard spheres (DHS) with or without uniaxial anisotropy and frozen in position on a disordered structure by tempered Monte Carlo simulations. The crucial point is to consider an anisotropic structure, obtained from the liquid state of the dipolar hard spheres fluid, frozen in its polarized state at low temperature. The freezing inverse t…
▽ More
We study the magnetic phase diagram of an ensemble of dipolar hard spheres (DHS) with or without uniaxial anisotropy and frozen in position on a disordered structure by tempered Monte Carlo simulations. The crucial point is to consider an anisotropic structure, obtained from the liquid state of the dipolar hard spheres fluid, frozen in its polarized state at low temperature. The freezing inverse temperature $β_f$ determines the degree of anisotropy of the structure which is quantified through a structural nematic order parameter, $λ_s$. The case of the non zero uniaxial anisotropy is considered only in its infinitely strong strength limit where the system transforms in a dipolar Ising model (DIM). The important finding of this work is that both the DHS and the DIM with a frozen structure build in this way present a ferromagnetic phase at volume fractions below the threshold value where the corresponding isotropic DHS systems exhibit a spin glass phase at low temperature.
△ Less
Submitted 25 April, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
The Boltzmann equation for hard potentials with integrable angular transition: Coerciveness, exponential tails rates, and Lebesgue integrability
Authors:
Ricardo J. Alonso,
Irene M. Gamba
Abstract:
This manuscript focus on an extensive survey with new techniques on the problem of solving the Boltzmann flow by bringing a unified approach to the Cauchy problem to homogeneous kinetic equations with Boltzmann-like collision operators under integrability assumption of the scattering profile in the particle-particle interaction mechanism. The work focuses on the relevant hard potential case where…
▽ More
This manuscript focus on an extensive survey with new techniques on the problem of solving the Boltzmann flow by bringing a unified approach to the Cauchy problem to homogeneous kinetic equations with Boltzmann-like collision operators under integrability assumption of the scattering profile in the particle-particle interaction mechanism. The work focuses on the relevant hard potential case where the solution properties are studied with a modern take. While many of the discussed results can be found the literature spread over several papers along the years, we bring a complete program that includes a new approach to the existence and uniqueness theorem securing the well-posedness theory, to moments estimates, and integrability propagation for the homogeneous Boltzmann flow. In particular, a detailed calculation of classical polynomial moments upper bounds as function of the coerciveness is described, which characterized the rate of exponential moments obtained by summability of the polynomial ones. In addition, a proof of uniform propagation of $L^\infty$ regularity under general integrable scattering kernels is performed. Along the way, constants appearing in estimates are carefully calculated, improving most of previous exiting results in the literature. For the non expert reader we also include a general discussion of the basic elements of the Boltzmann model and important key results for the understanding of the mathematical discussion of the equation and include an extensive set of references that enrich and motivate further discussions on Boltzmann flows for broader gas modeling configuration such as gas mixtures systems, polyatomic gases, multilinear collisional forms such as, ternary or quartic, to those derived from symmetry braking quantum mean field theories or weak turbulence models from spectral energy waves in classical fluid.
△ Less
Submitted 10 January, 2023; v1 submitted 16 November, 2022;
originally announced November 2022.
-
Filtration-Domination in Bifiltered Graphs
Authors:
Ángel Javier Alonso,
Michael Kerber,
Siddharth Pritam
Abstract:
Bifiltered graphs are a versatile tool for modelling relations between data points across multiple grades of a two-dimensional scale. They are especially popular in topological data analysis, where the homological properties of the induced clique complexes are studied. To reduce the large size of these clique complexes, we identify filtration-dominated edges of the graph, whose removal preserves t…
▽ More
Bifiltered graphs are a versatile tool for modelling relations between data points across multiple grades of a two-dimensional scale. They are especially popular in topological data analysis, where the homological properties of the induced clique complexes are studied. To reduce the large size of these clique complexes, we identify filtration-dominated edges of the graph, whose removal preserves the relevant topological properties. We give two algorithms to detect filtration-dominated edges in a bifiltered graph and analyze their complexity. These two algorithms work directly on the bifiltered graph, without first extracting the clique complexes, which are generally much bigger. We present extensive experimental evaluation which shows that in most cases, more than 90% of the edges can be removed. In turn, we demonstrate that this often leads to a substantial speedup, and reduction in the memory usage, of the computational pipeline of multiparameter topological data analysis.
△ Less
Submitted 10 November, 2022;
originally announced November 2022.
-
One-dimensional inelastic Boltzmann equation: Regularity \& uniqueness of self-similar profiles for moderately hard potentials
Authors:
Ricardo J. Alonso,
Véronique Bagland,
José A. Cañizo,
Bertrand Lods,
Sebastian Throm
Abstract:
We prove uniqueness of self-similar profiles for the one-dimensional inelastic Boltzmann equation with moderately hard potentials, that is with collision kernel of the form | $\bullet$ | $γ$ for $γ$ > 0 small enough (explicitly quantified). Our result provides the first uniqueness statement for self-similar profiles of inelastic Boltzmann models allowing for strong inelasticity besides the explici…
▽ More
We prove uniqueness of self-similar profiles for the one-dimensional inelastic Boltzmann equation with moderately hard potentials, that is with collision kernel of the form | $\bullet$ | $γ$ for $γ$ > 0 small enough (explicitly quantified). Our result provides the first uniqueness statement for self-similar profiles of inelastic Boltzmann models allowing for strong inelasticity besides the explicitly solvable case of Maxwell interactions (corresponding to $γ$ = 0). Our approach relies on a perturbation argument from the corresponding Maxwell model through a careful study of the associated linearised operator. In particular, a part of the paper is devoted to the trend to equilibrium for the Maxwell model in suitable weighted Sobolev spaces, an extension of results which are known to hold in weaker topologies. Our results can be seen as a first step towards a full proof, in the one-dimensional setting, of a conjecture in Ernst \& Brito (2002) regarding the determination of the long-time behaviour of solutions to inelastic Boltzmann equation.
△ Less
Submitted 7 November, 2022;
originally announced November 2022.
-
Identification and reconstruction of low-energy electrons in the ProtoDUNE-SP detector
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
Z. Ahmad,
J. Ahmed,
B. Aimard,
F. Akbar,
K. Allison,
S. Alonso Monsalve,
M. Alrashed,
C. Alt,
A. Alton,
R. Alvarez,
P. Amedo,
J. Anderson
, et al. (1235 additional authors not shown)
Abstract:
Measurements of electrons from $ν_e$ interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is…
▽ More
Measurements of electrons from $ν_e$ interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is one of the prototypes for the DUNE far detector, built and operated at CERN as a charged particle test beam experiment. A sample of low-energy electrons produced by the decay of cosmic muons is selected with a purity of 95%. This sample is used to calibrate the low-energy electron energy scale with two techniques. An electron energy calibration based on a cosmic ray muon sample uses calibration constants derived from measured and simulated cosmic ray muon events. Another calibration technique makes use of the theoretically well-understood Michel electron energy spectrum to convert reconstructed charge to electron energy. In addition, the effects of detector response to low-energy electron energy scale and its resolution including readout electronics threshold effects are quantified. Finally, the relation between the theoretical and reconstructed low-energy electron energy spectrum is derived and the energy resolution is characterized. The low-energy electron selection presented here accounts for about 75% of the total electron deposited energy. After the addition of lost energy using a Monte Carlo simulation, the energy resolution improves from about 40% to 25% at 50~MeV. These results are used to validate the expected capabilities of the DUNE far detector to reconstruct low-energy electrons.
△ Less
Submitted 31 May, 2023; v1 submitted 2 November, 2022;
originally announced November 2022.
-
A sufficient condition for confinement in QCD
Authors:
J. L. Alonso,
C. Bouthelier-Madre,
J. Clemente-Gallardo,
D. Martínez-Crespo
Abstract:
This letter is about confinement in QCD. At the moment we have pictures of confinement to complete our understanding of the physics of strongly interacting particles, interaction which asks for confinement. As it is said in [1] : " In principle it should be possible to derive the confinement hypothesis from the QCD Lagrangian. At this time, no rigorous derivation exists, so it is not absolutely cl…
▽ More
This letter is about confinement in QCD. At the moment we have pictures of confinement to complete our understanding of the physics of strongly interacting particles, interaction which asks for confinement. As it is said in [1] : " In principle it should be possible to derive the confinement hypothesis from the QCD Lagrangian. At this time, no rigorous derivation exists, so it is not absolutely clear that the confinement hypothesis is a bone fide prediction of QCD" . In this letter we show that a sufficient (of course not necessary) condition for confinement is that topological structure of vacuum in Nature does not correspond to the $θ$-vacuum. Therefore, if different vacua with nontrivial winding number cannot be connected by tunneling, we obtain confinement as a consequence.
△ Less
Submitted 7 December, 2022; v1 submitted 2 November, 2022;
originally announced November 2022.
-
Prevention of core particle depletion in stellarators by turbulence
Authors:
H. Thienpondt,
J. M. García-Regaña,
I. Calvo,
J. A. Alonso,
J. L. Velasco,
A. González-Jerez,
M. Barnes,
K. Brunner,
O. Ford,
G. Fuchert,
J. Knauer,
E. Pasch,
L. Vanó,
the Wendelstein 7-X team
Abstract:
In reactor-relevant plasmas, neoclassical transport drives an outward particle flux in the core of large stellarators and predicts strongly hollow density profiles. However, this theoretical prediction is contradicted by experiments. In particular, in Wendelstein 7-X, the first large optimized stellarator, flat or weakly peaked density profiles are generally measured, indicating that neoclassical…
▽ More
In reactor-relevant plasmas, neoclassical transport drives an outward particle flux in the core of large stellarators and predicts strongly hollow density profiles. However, this theoretical prediction is contradicted by experiments. In particular, in Wendelstein 7-X, the first large optimized stellarator, flat or weakly peaked density profiles are generally measured, indicating that neoclassical theory is not sufficient and that an inward contribution to the particle flux is missing in the core. In this Research Letter, it is shown that the turbulent contribution to the particle flux can explain the difference between experimental measurements and neoclassical predictions. The results of this Research Letter also prove that theoretical and numerical tools are approaching the level of maturity needed for the prediction of equilibrium density profiles in stellarator plasmas, which is a fundamental requirement for the design of operation scenarios of present devices and future reactors.
△ Less
Submitted 8 November, 2023; v1 submitted 9 September, 2022;
originally announced September 2022.
-
MaRCoS, an open-source electronic control system for low-field MRI
Authors:
Vlad Negnevitsky,
Yolanda Vives-Gilabert,
José M. Algarín,
Lincoln Craven-Brightman,
Rubén Pellicer-Guridi,
Thomas O'Reilly,
Jason P. Stockmann,
Andrew Webb,
Joseba Alonso,
Benjamin Menküc
Abstract:
Every magnetic resonance imaging (MRI) device requires an electronic control system that handles pulse sequences and signal detection and processing. Here we provide details on the architecture and performance of MaRCoS, a MAgnetic Resonance COntrol System developed by an open international community of low-field MRI researchers. MaRCoS is inexpensive and can handle cycle-accurate sequences withou…
▽ More
Every magnetic resonance imaging (MRI) device requires an electronic control system that handles pulse sequences and signal detection and processing. Here we provide details on the architecture and performance of MaRCoS, a MAgnetic Resonance COntrol System developed by an open international community of low-field MRI researchers. MaRCoS is inexpensive and can handle cycle-accurate sequences without hard length limitations, rapid bursts of events, and arbitrary waveforms. It can also be easily adapted to meet further specifications required by the various academic and private institutions participating in its development. We describe the MaRCoS hardware, firmware and software that enable all of the above, including a Python-based graphical user interface for pulse sequence implementation, data processing and image reconstruction.
△ Less
Submitted 2 August, 2022;
originally announced August 2022.
-
A three-dimensional generalization of QRT maps
Authors:
Jaume Alonso,
Yuri B. Suris,
Kangning Wei
Abstract:
We propose a geometric construction of three-dimensional birational maps that preserve two pencils of quadrics. The maps act as compositions of involutions, which, in turn, act along the straight line generators of the quadrics of the first pencil and are defined by the intersections with quadrics of the second pencil. On each quadric of the first pencil, the maps act as two-dimensional QRT maps.…
▽ More
We propose a geometric construction of three-dimensional birational maps that preserve two pencils of quadrics. The maps act as compositions of involutions, which, in turn, act along the straight line generators of the quadrics of the first pencil and are defined by the intersections with quadrics of the second pencil. On each quadric of the first pencil, the maps act as two-dimensional QRT maps.
While these maps are of a pretty high degree in general, we find geometric conditions which guarantee that the degree is reduced to 3. The resulting degree 3 maps are illustrated by two known and two novel Kahan-type discretizations of three-dimensional Nambu systems, including the Euler top and the Zhukovski-Volterra gyrostat with two non-vanishing components of the gyrostatic momentum.
△ Less
Submitted 11 June, 2023; v1 submitted 13 July, 2022;
originally announced July 2022.
-
Reconstruction of interactions in the ProtoDUNE-SP detector with Pandora
Authors:
DUNE Collaboration,
A. Abed Abud,
B. Abi,
R. Acciarri,
M. A. Acero,
M. R. Adames,
G. Adamov,
M. Adamowski,
D. Adams,
M. Adinolfi,
C. Adriano,
A. Aduszkiewicz,
J. Aguilar,
Z. Ahmad,
J. Ahmed,
B. Aimard,
F. Akbar,
B. Ali-Mohammadzadeh,
K. Allison,
S. Alonso Monsalve,
M. AlRashed,
C. Alt,
A. Alton,
R. Alvarez,
P. Amedo
, et al. (1203 additional authors not shown)
Abstract:
The Pandora Software Development Kit and algorithm libraries provide pattern-recognition logic essential to the reconstruction of particle interactions in liquid argon time projection chamber detectors. Pandora is the primary event reconstruction software used at ProtoDUNE-SP, a prototype for the Deep Underground Neutrino Experiment far detector. ProtoDUNE-SP, located at CERN, is exposed to a char…
▽ More
The Pandora Software Development Kit and algorithm libraries provide pattern-recognition logic essential to the reconstruction of particle interactions in liquid argon time projection chamber detectors. Pandora is the primary event reconstruction software used at ProtoDUNE-SP, a prototype for the Deep Underground Neutrino Experiment far detector. ProtoDUNE-SP, located at CERN, is exposed to a charged-particle test beam. This paper gives an overview of the Pandora reconstruction algorithms and how they have been tailored for use at ProtoDUNE-SP. In complex events with numerous cosmic-ray and beam background particles, the simulated reconstruction and identification efficiency for triggered test-beam particles is above 80% for the majority of particle type and beam momentum combinations. Specifically, simulated 1 GeV/$c$ charged pions and protons are correctly reconstructed and identified with efficiencies of 86.1$\pm0.6$% and 84.1$\pm0.6$%, respectively. The efficiencies measured for test-beam data are shown to be within 5% of those predicted by the simulation.
△ Less
Submitted 17 July, 2023; v1 submitted 29 June, 2022;
originally announced June 2022.
-
Bayesian additive regression trees for probabilistic programming
Authors:
Miriana Quiroga,
Pablo G Garay,
Juan M. Alonso,
Juan Martin Loyola,
Osvaldo A Martin
Abstract:
Bayesian additive regression trees (BART) is a non-parametric method to approximate functions. It is a black-box method based on the sum of many trees where priors are used to regularize inference, mainly by restricting trees' learning capacity so that no individual tree is able to explain the data, but rather the sum of trees. We discuss BART in the context of probabilistic programming languages…
▽ More
Bayesian additive regression trees (BART) is a non-parametric method to approximate functions. It is a black-box method based on the sum of many trees where priors are used to regularize inference, mainly by restricting trees' learning capacity so that no individual tree is able to explain the data, but rather the sum of trees. We discuss BART in the context of probabilistic programming languages (PPL), i.e., we present BART as a primitive that can be used as a component of a probabilistic model rather than as a standalone model. Specifically, we introduce the Python library PyMC-BART, which works by extending PyMC, a library for probabilistic programming. We showcase a few examples of models that can be built using PyMC-BART, discuss recommendations for the selection of hyperparameters, and finally, we close with limitations of our implementation and future directions for improvement.
△ Less
Submitted 15 August, 2023; v1 submitted 7 June, 2022;
originally announced June 2022.