-
A Comprehensive Guide to Combining R and Python code for Data Science, Machine Learning and Reinforcement Learning
Authors:
Alejandro L. García Navarro,
Nataliia Koneva,
Alfonso Sánchez-Macián,
José Alberto Hernández
Abstract:
Python has gained widespread popularity in the fields of machine learning, artificial intelligence, and data engineering due to its effectiveness and extensive libraries. R, on its side, remains a dominant language for statistical analysis and visualization. However, certain libraries have become outdated, limiting their functionality and performance. Users can use Python's advanced machine learni…
▽ More
Python has gained widespread popularity in the fields of machine learning, artificial intelligence, and data engineering due to its effectiveness and extensive libraries. R, on its side, remains a dominant language for statistical analysis and visualization. However, certain libraries have become outdated, limiting their functionality and performance. Users can use Python's advanced machine learning and AI capabilities alongside R's robust statistical packages by combining these two programming languages. This paper explores using R's reticulate package to call Python from R, providing practical examples and highlighting scenarios where this integration enhances productivity and analytical capabilities. With a few hello-world code snippets, we demonstrate how to run Python's scikit-learn, pytorch and OpenAI gym libraries for building Machine Learning, Deep Learning, and Reinforcement Learning projects easily.
△ Less
Submitted 19 July, 2024;
originally announced July 2024.
-
Non-linear gauge-Higgs CP violation
Authors:
Akanksha Bhardwaj,
Christoph Englert,
Dorival Gonçalves,
Alberto Navarro
Abstract:
A critical element of the LHC physics program is the search for an additional source of CP violation. This is largely unexplored in the context of non-linear Higgs physics, which is naturally described in Higgs Effective Field Theory (HEFT). Relevant new higher-dimensional operators modify the production rate and branching ratios of the Higgs boson, de-correlating different Higgs multiplicities. I…
▽ More
A critical element of the LHC physics program is the search for an additional source of CP violation. This is largely unexplored in the context of non-linear Higgs physics, which is naturally described in Higgs Effective Field Theory (HEFT). Relevant new higher-dimensional operators modify the production rate and branching ratios of the Higgs boson, de-correlating different Higgs multiplicities. In this work, we consider single Higgs and Higgs pair production via weak boson fusion from the perspective of gauge-Higgs CP violation through the lens of Higgs non-linearity. This generalizes existing rate-based searches and analyses by the ATLAS and CMS experiments. Particular focus is given to the phenomenological differences in the expected BSM sensitivity pattern when comparing HEFT constraints with Standard Model Effective Field Theory (SMEFT) limits.
△ Less
Submitted 19 July, 2024;
originally announced July 2024.
-
Hardware Neural Control of CartPole and F1TENTH Race Car
Authors:
Marcin Paluch,
Florian Bolli,
Xiang Deng,
Antonio Rios Navarro,
Chang Gao,
Tobi Delbruck
Abstract:
Nonlinear model predictive control (NMPC) has proven to be an effective control method, but it is expensive to compute. This work demonstrates the use of hardware FPGA neural network controllers trained to imitate NMPC with supervised learning. We use these Neural Controllers (NCs) implemented on inexpensive embedded FPGA hardware for high frequency control on physical cartpole and F1TENTH race ca…
▽ More
Nonlinear model predictive control (NMPC) has proven to be an effective control method, but it is expensive to compute. This work demonstrates the use of hardware FPGA neural network controllers trained to imitate NMPC with supervised learning. We use these Neural Controllers (NCs) implemented on inexpensive embedded FPGA hardware for high frequency control on physical cartpole and F1TENTH race car. Our results show that the NCs match the control performance of the NMPCs in simulation and outperform it in reality, due to the faster control rate that is afforded by the quick FPGA NC inference. We demonstrate kHz control rates for a physical cartpole and offloading control to the FPGA hardware on the F1TENTH car. Code and hardware implementation for this paper are available at https:// github.com/SensorsINI/Neural-Control-Tools.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
On the impact of VR/AR applications on optical transport networks: First experiments with Meta Quest 3 gaming and conferencing application
Authors:
C. de Quinto,
A. Navarro,
G. Otero,
N. Koneva,
J. A. Hernández,
M. Quagliotti,
A. Sánchez-Macian,
F. Arpanaei,
P. Reviriego,
Ó. González de Dios,
J. M. Rivas-Moscoso,
E. Riccardi,
D. Larrabeiti
Abstract:
With the advent of next-generation AR/VR headsets, many of them with affordable prices, telecom operators have forecasted an explosive growth of traffic in their networks. Penetration of AR/VR services and applications is estimated to grow exponentially in the next few years. This work attempts to shed light on the bandwidth capacity requirements and latency of popular AR/VR applications with four…
▽ More
With the advent of next-generation AR/VR headsets, many of them with affordable prices, telecom operators have forecasted an explosive growth of traffic in their networks. Penetration of AR/VR services and applications is estimated to grow exponentially in the next few years. This work attempts to shed light on the bandwidth capacity requirements and latency of popular AR/VR applications with four different real experimental settings on the Meta Quest 3 headsets, and their potential impact on the network.
△ Less
Submitted 29 July, 2024; v1 submitted 10 July, 2024;
originally announced July 2024.
-
Analytical Insights on Hadronic Top Quark Polarimetry
Authors:
Zhongtian Dong,
Dorival Gonçalves,
Kyoungchul Kong,
Andrew J. Larkoski,
Alberto Navarro
Abstract:
Top quark polarization provides an important tool for studying its production mechanisms, spin correlations, top quark properties, and new physics searches. Unlike lighter quarks, the top quark's polarization remains intact until its decay, enabling precise spin measurements. While the down-type fermions from $W$ boson decay are known to be effective spin analyzers, charged leptons have typically…
▽ More
Top quark polarization provides an important tool for studying its production mechanisms, spin correlations, top quark properties, and new physics searches. Unlike lighter quarks, the top quark's polarization remains intact until its decay, enabling precise spin measurements. While the down-type fermions from $W$ boson decay are known to be effective spin analyzers, charged leptons have typically been the main target for most analyses. In this paper, we investigate the relevance of global jet dynamics -- considering kinematics, jet charges, and particle multiplicity -- for hadronic top quark polarimetry. The formalism used allows for analytical derivations obtained throughout the manuscript, offering deeper insights into the corresponding phenomenology.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Segment-Based Interactive Machine Translation for Pre-trained Models
Authors:
Angel Navarro,
Francisco Casacuberta
Abstract:
Pre-trained large language models (LLM) are starting to be widely used in many applications. In this work, we explore the use of these models in interactive machine translation (IMT) environments. In particular, we have chosen mBART (multilingual Bidirectional and Auto-Regressive Transformer) and mT5 (multilingual Text-to-Text Transfer Transformer) as the LLMs to perform our experiments. The syste…
▽ More
Pre-trained large language models (LLM) are starting to be widely used in many applications. In this work, we explore the use of these models in interactive machine translation (IMT) environments. In particular, we have chosen mBART (multilingual Bidirectional and Auto-Regressive Transformer) and mT5 (multilingual Text-to-Text Transfer Transformer) as the LLMs to perform our experiments. The system generates perfect translations interactively using the feedback provided by the user at each iteration. The Neural Machine Translation (NMT) model generates a preliminary hypothesis with the feedback, and the user validates new correct segments and performs a word correction--repeating the process until the sentence is correctly translated. We compared the performance of mBART, mT5, and a state-of-the-art (SoTA) machine translation model on a benchmark dataset regarding user effort, Word Stroke Ratio (WSR), Key Stroke Ratio (KSR), and Mouse Action Ratio (MAR). The experimental results indicate that mBART performed comparably with SoTA models, suggesting that it is a viable option for this field of IMT. The implications of this finding extend to the development of new machine translation models for interactive environments, as it indicates that some novel pre-trained models exhibit SoTA performance in this domain, highlighting the potential benefits of adapting these models to specific needs.
△ Less
Submitted 9 July, 2024;
originally announced July 2024.
-
Overeducation under different macroeconomic conditions: The case of Spanish university graduates
Authors:
Maite Blázquez Cuesta,
Marco A. Pérez Navarro,
Rocío Sánchez-Mangas
Abstract:
This paper examines the incidence and persistence of overeducation in the early careers of Spanish university graduates. We investigate the role played by the business cycle and field of study and their interaction in shaping both phenomena. We also analyse the relevance of specific types of knowledge and skills as driving factors in reducing overeducation risk. We use data from the Survey on the…
▽ More
This paper examines the incidence and persistence of overeducation in the early careers of Spanish university graduates. We investigate the role played by the business cycle and field of study and their interaction in shaping both phenomena. We also analyse the relevance of specific types of knowledge and skills as driving factors in reducing overeducation risk. We use data from the Survey on the Labour Insertion of University Graduates (EILU) conducted by the Spanish National Statistics Institute in 2014 and 2019. The survey collects rich information on cohorts that graduated in the 2009/2010 and 2014/2015 academic years during the Great Recession and the subsequent economic recovery, respectively. Our results show, first, the relevance of the economic scenario when graduates enter the labour market. Graduation during a recession increased overeducation risk and persistence. Second, a clear heterogeneous pattern occurs across fields of study, with health sciences graduates displaying better performance in terms of both overeducation incidence and persistence and less impact of the business cycle. Third, we find evidence that some transversal skills (language, IT, management) can help to reduce overeducation risk in the absence of specific knowledge required for the job, thus indicating some kind of compensatory role. Finally, our findings have important policy implications. Overeducation, and more importantly overeducation persistence, imply a non-neglectable misallocation of resources. Therefore, policymakers need to address this issue in the design of education and labour market policies.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Hadronic Top Quark Polarimetry with ParticleNet
Authors:
Zhongtian Dong,
Dorival Gonçalves,
Kyoungchul Kong,
Andrew J. Larkoski,
Alberto Navarro
Abstract:
Precision studies for top quark physics are a cornerstone of the Large Hadron Collider program. Polarization, probed through decay kinematics, provides a unique tool to scrutinize the top quark across its various production modes and to explore potential new physics effects. However, the top quark most often decays hadronically, for which unambiguous identification of its decay products sensitive…
▽ More
Precision studies for top quark physics are a cornerstone of the Large Hadron Collider program. Polarization, probed through decay kinematics, provides a unique tool to scrutinize the top quark across its various production modes and to explore potential new physics effects. However, the top quark most often decays hadronically, for which unambiguous identification of its decay products sensitive to top quark polarization is not possible. In this Letter, we introduce a jet flavor tagging method to significantly improve spin analyzing power in hadronic decays, going beyond exclusive kinematic information employed in previous studies. We provide parametric estimates of the improvement from flavor tagging with any set of measured observables and demonstrate this in practice on simulated data using a Graph Neural Network (GNN). We find that the spin analyzing power in hadronic decays can improve by approximately 20% (40%) compared to the kinematic approach, assuming an efficiency of 0.5 (0.2) for the network.
△ Less
Submitted 1 July, 2024;
originally announced July 2024.
-
CAT: Cellular Automata on Tensor cores
Authors:
Cristóbal A. Navarro,
Felipe A. Quezada,
Enzo Meneses,
Héctor Ferrada,
Nancy Hitschfeld
Abstract:
Cellular automata (CA) are simulation models that can produce complex emergent behaviors from simple local rules. Although state-of-the-art GPU solutions are already fast due to their data-parallel nature, their performance can rapidly degrade in CA with a large neighborhood radius. With the inclusion of tensor cores across the entire GPU ecosystem, interest has grown in finding ways to leverage t…
▽ More
Cellular automata (CA) are simulation models that can produce complex emergent behaviors from simple local rules. Although state-of-the-art GPU solutions are already fast due to their data-parallel nature, their performance can rapidly degrade in CA with a large neighborhood radius. With the inclusion of tensor cores across the entire GPU ecosystem, interest has grown in finding ways to leverage these fast units outside the field of artificial intelligence, which was their original purpose. In this work, we present CAT, a GPU tensor core approach that can accelerate CA in which the cell transition function acts on a weighted summation of its neighborhood. CAT is evaluated theoretically, using an extended PRAM cost model, as well as empirically using the Larger Than Life (LTL) family of CA as case studies. The results confirm that the cost model is accurate, showing that CAT exhibits constant time throughout the entire radius range $1 \le r \le 16$, and its theoretical speedups agree with the empirical results. At low radius $r=1,2$, CAT is competitive and is only surpassed by the fastest state-of-the-art GPU solution. Starting from $r=3$, CAT progressively outperforms all other approaches, reaching speedups of up to $101\times$ over a GPU baseline and up to $\sim 14\times$ over the fastest state-of-the-art GPU approach. In terms of energy efficiency, CAT is competitive in the range $1 \le r \le 4$ and from $r \ge 5$ it is the most energy efficient approach. As for performance scaling across GPU architectures, CAT shows a promising trend that if continues for future generations, it would increase its performance at a higher rate than classical GPU solutions. The results obtained in this work put CAT as an attractive GPU approach for scientists that need to study emerging phenomena on CA with large neighborhood radius.
△ Less
Submitted 25 June, 2024;
originally announced June 2024.
-
von Mises Quasi-Processes for Bayesian Circular Regression
Authors:
Yarden Cohen,
Alexandre Khae Wu Navarro,
Jes Frellsen,
Richard E. Turner,
Raziel Riemer,
Ari Pakman
Abstract:
The need for regression models to predict circular values arises in many scientific fields. In this work we explore a family of expressive and interpretable distributions over circle-valued random functions related to Gaussian processes targeting two Euclidean dimensions conditioned on the unit circle. The resulting probability model has connections with continuous spin models in statistical physi…
▽ More
The need for regression models to predict circular values arises in many scientific fields. In this work we explore a family of expressive and interpretable distributions over circle-valued random functions related to Gaussian processes targeting two Euclidean dimensions conditioned on the unit circle. The resulting probability model has connections with continuous spin models in statistical physics. Moreover, its density is very simple and has maximum-entropy, unlike previous Gaussian process-based approaches, which use wrapping or radial marginalization. For posterior inference, we introduce a new Stratonovich-like augmentation that lends itself to fast Markov Chain Monte Carlo sampling. We argue that transductive learning in these models favors a Bayesian approach to the parameters. We present experiments applying this model to the prediction of (i) wind directions and (ii) the percentage of the running gait cycle as a function of joint angles.
△ Less
Submitted 18 June, 2024;
originally announced June 2024.
-
Reinforcement-Learning based routing for packet-optical networks with hybrid telemetry
Authors:
A. L. García Navarro,
Nataliia Koneva,
Alfonso Sánchez-Macián,
José Alberto Hernández,
Óscar González de Dios,
J. M. Rivas-Moscoso
Abstract:
This article provides a methodology and open-source implementation of Reinforcement Learning algorithms for finding optimal routes in a packet-optical network scenario. The algorithm uses measurements provided by the physical layer (pre-FEC bit error rate and propagation delay) and the link layer (link load) to configure a set of latency-based rewards and penalties based on such measurements. Then…
▽ More
This article provides a methodology and open-source implementation of Reinforcement Learning algorithms for finding optimal routes in a packet-optical network scenario. The algorithm uses measurements provided by the physical layer (pre-FEC bit error rate and propagation delay) and the link layer (link load) to configure a set of latency-based rewards and penalties based on such measurements. Then, the algorithm executes Q-learning based on this set of rewards for finding the optimal routing strategies. It is further shown that the algorithm dynamically adapts to changing network conditions by re-calculating optimal policies upon either link load changes or link degradation as measured by pre-FEC BER.
△ Less
Submitted 21 June, 2024; v1 submitted 18 June, 2024;
originally announced June 2024.
-
Assessing the Impact of Alpha Particles on Thermal Confinement in JET D-T Plasmas through Global GENE-Tango Simulations
Authors:
A. Di Siena,
J. Garcia,
R. Bilato,
K. Kirov,
J. Varela A. Banon Navarro,
Hyun-Tae Kim,
C. Challis,
J. Hobirk,
A. Kappatou,
E. Lerche,
D. Spong,
C. Angioni,
T. Gorler,
E. Poli,
M. Bergmann,
F. Jenko,
JET contributors
Abstract:
The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially…
▽ More
The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge - despite this being the shot with the highest fusion output - is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
Comparing bulge RR Lyrae stars with bulge giants -- Insight from 3D kinematics
Authors:
J. Olivares Carvajal,
M. Zoccali,
M. De Leo,
R. Contreras Ramos,
C. Quezada,
Á. Rojas-Arriagada,
E. Valenti,
R. Albarracín,
Á. Valenzuela Navarro
Abstract:
The structure and kinematics of the old component of the Galactic bulge are still a matter of debate. The bulk of the bulge as traced by red clump stars includes two main components, which are usually identified as the metal-rich and metal-poor components. They have different shapes, kinematics, mean metallicities, and alpha-element abundances. It is our current understanding that they are associa…
▽ More
The structure and kinematics of the old component of the Galactic bulge are still a matter of debate. The bulk of the bulge as traced by red clump stars includes two main components, which are usually identified as the metal-rich and metal-poor components. They have different shapes, kinematics, mean metallicities, and alpha-element abundances. It is our current understanding that they are associated with a bar and a spheroid, respectively. On the other hand, RR Lyrae variables trace the oldest population of the bulge. While it would be natural to think that they follow the structure and kinematics of the metal-poor component, the data analysed in the literature show conflicting results. We aim to derive a rotation curve for bulge RR Lyrae stars in order to determine that the old component traced by these stars is distinct from the two main components observed in the Galactic bulge. This paper combines APOGEE-2S spectra with OGLE-IV light curves, near-IR photometry, and proper motions from the VISTA Variables in the Vía Láctea survey for 4197 RR Lyrae stars. Six-dimensional phase-space coordinates were used to calculate orbits within an updated Galactic potential and to isolate the stars. The stars that stay confined within the bulge represent 57% of our sample. Our results show that bulge RR Lyrae variables rotate more slowly than metal-rich red clump stars and have a lower velocity dispersion. Their kinematics is compatible with them being the low-metallicity tail of the metal-poor component. We confirm that a rather large fraction of halo RR Lyrae stars pass by the bulge within their orbits, increasing the velocity dispersion. A proper orbital analysis is therefore critical to isolate bona fide bulge variables. Finally, bulge RR Lyrae seem to trace a spheroidal component, although the current data do now allow us to reach a firm conclusion about the spatial distribution.
△ Less
Submitted 16 May, 2024; v1 submitted 14 May, 2024;
originally announced May 2024.
-
Radiation Surveys in Active Nuclear Facilities with Heterogeneous Collaborative Mobile Robots
Authors:
Mitchell Pryor,
Alex Navarro,
Janak Panthi,
Kevin Torres,
Mary Tebben,
Daniel Meza,
Caleb Horan,
Alex Macris
Abstract:
Nuclear facilities must routinely survey their infrastructure for radiation contamination. Generally, this is done by trained professionals, wearing personal protective equipment (PPE) that swipe potentially contaminated surfaces and test the wipes under detectors. This approach leaves personnel vulnerable to radiation exposure and is not comprehensive. Robots address these inadequacies, offering…
▽ More
Nuclear facilities must routinely survey their infrastructure for radiation contamination. Generally, this is done by trained professionals, wearing personal protective equipment (PPE) that swipe potentially contaminated surfaces and test the wipes under detectors. This approach leaves personnel vulnerable to radiation exposure and is not comprehensive. Robots address these inadequacies, offering a cost-effective solution with negligible downtime. We present a Robot Radiation Survey System (RRSS): a heterogeneous robot team to perform comprehensive alpha/beta/gamma radiation surveys. The RRSS system members, core capabilities, and comprehensive survey plan are addresses in this paper.
△ Less
Submitted 22 February, 2024;
originally announced February 2024.
-
Global gyrokinetic analysis of Wendelstein 7-X discharge: unveiling the importance of trapped-electron-mode and electron-temperature-gradient turbulence
Authors:
Felix Wilms,
Alejandro Bañón Navarro,
Thomas Windisch,
Sergey Bozhenkov,
Felix Warmer,
Golo Fuchert,
Oliver Ford,
Daihong Zhang,
Torsten Stange,
Frank Jenko,
the W7-X Team
Abstract:
We present the first nonlinear, gyrokinetic, radially global simulation of a discharge of the Wendelstein 7-X-like stellarator (W7-X), including kinetic electrons, an equilibrium radial electric field, as well as electromagnetic and collisional effects. By comparison against flux-tube and full-flux-surface simulations, we assess the impact of the equilibrium ExB-flow and flow shear on the stabilis…
▽ More
We present the first nonlinear, gyrokinetic, radially global simulation of a discharge of the Wendelstein 7-X-like stellarator (W7-X), including kinetic electrons, an equilibrium radial electric field, as well as electromagnetic and collisional effects. By comparison against flux-tube and full-flux-surface simulations, we assess the impact of the equilibrium ExB-flow and flow shear on the stabilisation of turbulence. In contrast to the existing literature, we further provide substantial evidence for the turbulent electron heat flux being driven by trapped-electron-mode (TEM) and electron-temperature-gradient (ETG) turbulence in the core of the plasma. The former manifests as a hybrid together with ion-temperature-gradient (ITG) turbulence and is primarily driven by the finite electron temperature gradient, which has largely been neglected in nonlinear stellarator simulations presented in the existing literature.
△ Less
Submitted 22 February, 2024;
originally announced February 2024.
-
New candidate hypervelocity red clump stars in the inner Galactic bulge
Authors:
A. Luna,
T. Marchetti,
M. Rejkuba,
N. W. C. Leigh,
J. Alonso-García,
A. Valenzuela Navarro,
D. Minniti,
L. C. Smith
Abstract:
We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the de…
▽ More
We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the detection of HVSs, ejection mechanism models can be constrained by exploring the stellar dynamics in the Galactic centre through a recent stellar interaction with the SMBH. Based on a previously developed methodology by our group, we searched with a sample of preliminary data from version 2 of the Vista Variables in the Via Lactea (VVV) Infrared Astrometric Catalogue (VIRAC2) and Gaia DR3 data, including accurate optical and NIR proper motions. This search resulted in a sample of 46 stars with transverse velocities larger than the local escape velocity within the Galactic bulge, of which 4 are prime candidate HVSs with high-proper motions consistent with being ejections from the Galactic centre. Adding to that, we studied a sample of reddened stars without a Gaia DR3 counterpart and found 481 stars with transverse velocities larger than the local escape velocity, from which 65 stars have proper motions pointing out of the Galactic centre and are candidate HVSs. In total, we found 69 candidate HVSs pointing away from the Galactic centre with transverse velocities larger than the local escape velocity.
△ Less
Submitted 10 January, 2024;
originally announced January 2024.
-
The OBDT-theta board: time digitization for the theta view of Drift Tubes chambers
Authors:
J. Sastre,
C. F. Bedoya,
S. Cuadrado,
J. Cuchillo,
D. Francia,
C. de Lara,
A. Navarro,
R. Paz,
I. Redondo,
D. Redondo
Abstract:
We present the design and performance of the On-Board electronics for the Drift Tubes (OBDT) for the superlayer theta along the direction parallel to the beam-line, the new board built to substitute part of the CMS DT Muon on-detector electronics. The OBDT-theta is responsible for the time digitization of the DT chamber signals for the theta view, allowing further tracking and triggering of the ba…
▽ More
We present the design and performance of the On-Board electronics for the Drift Tubes (OBDT) for the superlayer theta along the direction parallel to the beam-line, the new board built to substitute part of the CMS DT Muon on-detector electronics. The OBDT-theta is responsible for the time digitization of the DT chamber signals for the theta view, allowing further tracking and triggering of the barrel muons. It is also in charge of part of the slow-control of the DT chamber inner electronics in the theta view. Prototypes of the OBDT-theta board are under validation in different laboratories in CERN, as well as in demonstrator chambers installed in the CMS experiment. This allows evaluation of the full functionality of the boards in real conditions, showing very satisfactory results.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Mancha3D code: Multi-purpose Advanced Non-ideal MHD Code for High resolution simulations in Astrophysics
Authors:
M. Modestov,
E. Khomenko,
N. Vitas,
A. de Vicente,
A. Navarro,
P. A. Gonzalez-Morales,
M. Collados,
T. Felipe,
D. Martinez-Gomez,
P. Hunana,
M. Luna,
M. Koll Pistarini,
B. Popescu Braileanu,
A. Perdomo Garcia,
V. Liakh,
I. Santamaria,
M. M. Gomez Miguez
Abstract:
The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, suc…
▽ More
The Mancha3D code is a versatile tool for numerical simulations of magnetohydrodynamic processes in solar/stellar atmospheres. The code includes non-ideal physics derived from plasma partial ionization, a realistic equation of state and radiative transfer, which allows performing high quality realistic simulations of magneto-convection, as well as idealized simulations of particular processes, such as wave propagation, instabilities or energetic events. The paper summarizes the equations and methods used in the Mancha3D code. It also describes its numerical stability and parallel performance and efficiency. The code is based on a finite difference discretization and memory-saving Runge-Kutta (RK) scheme. It handles non-ideal effects through super-time stepping and Hall diffusion schemes, and takes into account thermal conduction by solving an additional hyperbolic equation for the heat flux. The code is easily configurable to perform different kinds of simulations. Several examples of the code usage are given. It is demonstrated that splitting variables into equilibrium and perturbation parts is essential for simulations of wave propagation in a static background. A perfectly matched layer (PML) boundary condition built into the code greatly facilitates a non-reflective open boundary implementation. Spatial filtering is an important numerical remedy to eliminate grid-size perturbations enhancing the code stability. Parallel performance analysis reveals that the code is strongly memory bound, which is a natural consequence of the numerical techniques used, such as split variables and PML boundary conditions. Both strong and weak scalings show adequate performance up till several thousands of CPUs.
△ Less
Submitted 7 December, 2023;
originally announced December 2023.
-
Assessing global ion thermal confinement in critical-gradient-optimized stellarators
Authors:
A. Bañón Navarro,
G. T. Roberg-Clark,
G. G. Plunk,
D. Fernando,
A. Di Siena,
F. Wilms,
F. Jenko
Abstract:
We investigate the confinement properties of two recently devised quasi-helically symmetric stellarator configurations, HSK and QSTK. Both have been optimized for large critical gradients of the ion temperature gradient mode, which is an important driver of turbulent transport in magnetic confinement fusion devices. To predict the resulting core plasma profiles, we utilize an advanced theoretical…
▽ More
We investigate the confinement properties of two recently devised quasi-helically symmetric stellarator configurations, HSK and QSTK. Both have been optimized for large critical gradients of the ion temperature gradient mode, which is an important driver of turbulent transport in magnetic confinement fusion devices. To predict the resulting core plasma profiles, we utilize an advanced theoretical framework based on the gyrokinetic codes GENE and GENE-3D, coupled to the transport code TANGO. Compared to the HSX stellarator, both HSK and QSTK achieve significantly higher core-to-edge temperature ratios, partly thanks to their smaller aspect ratios, with the other part due to more detailed shaping of the magnetic geometry achieved during optimization. The computed confinement time, however, is less sensitive to core temperature than edge temperature, simply due to the disproportionate influence the edge has on stored plasma energy. We therefore emphasize the possible benefits of further optimizing turbulence in the outer core region, and the need to include accurate modelling of confinement in the edge region in order to assess overall plasma performance of turbulence optimized stellarators.
△ Less
Submitted 28 October, 2023;
originally announced October 2023.
-
FUTURE-AI: International consensus guideline for trustworthy and deployable artificial intelligence in healthcare
Authors:
Karim Lekadir,
Aasa Feragen,
Abdul Joseph Fofanah,
Alejandro F Frangi,
Alena Buyx,
Anais Emelie,
Andrea Lara,
Antonio R Porras,
An-Wen Chan,
Arcadi Navarro,
Ben Glocker,
Benard O Botwe,
Bishesh Khanal,
Brigit Beger,
Carol C Wu,
Celia Cintas,
Curtis P Langlotz,
Daniel Rueckert,
Deogratias Mzurikwao,
Dimitrios I Fotiadis,
Doszhan Zhussupov,
Enzo Ferrante,
Erik Meijering,
Eva Weicken,
Fabio A González
, et al. (95 additional authors not shown)
Abstract:
Despite major advances in artificial intelligence (AI) for medicine and healthcare, the deployment and adoption of AI technologies remain limited in real-world clinical practice. In recent years, concerns have been raised about the technical, clinical, ethical and legal risks associated with medical AI. To increase real world adoption, it is essential that medical AI tools are trusted and accepted…
▽ More
Despite major advances in artificial intelligence (AI) for medicine and healthcare, the deployment and adoption of AI technologies remain limited in real-world clinical practice. In recent years, concerns have been raised about the technical, clinical, ethical and legal risks associated with medical AI. To increase real world adoption, it is essential that medical AI tools are trusted and accepted by patients, clinicians, health organisations and authorities. This work describes the FUTURE-AI guideline as the first international consensus framework for guiding the development and deployment of trustworthy AI tools in healthcare. The FUTURE-AI consortium was founded in 2021 and currently comprises 118 inter-disciplinary experts from 51 countries representing all continents, including AI scientists, clinicians, ethicists, and social scientists. Over a two-year period, the consortium defined guiding principles and best practices for trustworthy AI through an iterative process comprising an in-depth literature review, a modified Delphi survey, and online consensus meetings. The FUTURE-AI framework was established based on 6 guiding principles for trustworthy AI in healthcare, i.e. Fairness, Universality, Traceability, Usability, Robustness and Explainability. Through consensus, a set of 28 best practices were defined, addressing technical, clinical, legal and socio-ethical dimensions. The recommendations cover the entire lifecycle of medical AI, from design, development and validation to regulation, deployment, and monitoring. FUTURE-AI is a risk-informed, assumption-free guideline which provides a structured approach for constructing medical AI tools that will be trusted, deployed and adopted in real-world practice. Researchers are encouraged to take the recommendations into account in proof-of-concept stages to facilitate future translation towards clinical practice of medical AI.
△ Less
Submitted 8 July, 2024; v1 submitted 11 August, 2023;
originally announced September 2023.
-
Global Optimization-Based Calibration Algorithm for a 2D Distributed Hydrologic-Hydrodynamic and Water Quality Model
Authors:
Marcus N. Gomes Jr.,
Marcio H. Giacomoni,
Fabricio A. R. Navarro,
Eduardo M. Mendiondo
Abstract:
Hydrodynamic models with rain-on-the-grid capabilities are usually computationally expensive. This makes the use of automatic calibration algorithms hard to apply due to the large number of model runs. However, with the recent advances in parallel processing, computational resources, and increasing high-resolution climatologic and GIS data, high-resolution hydrodynamic models can be used for optim…
▽ More
Hydrodynamic models with rain-on-the-grid capabilities are usually computationally expensive. This makes the use of automatic calibration algorithms hard to apply due to the large number of model runs. However, with the recent advances in parallel processing, computational resources, and increasing high-resolution climatologic and GIS data, high-resolution hydrodynamic models can be used for optimization-based calibration. This paper presents a global optimization-based algorithm to calibrate a fully distributed hydrologic-hydrodynamic and water quality model (HydroPol2D) using observed data (i.e., discharge, or pollutant concentration) as input. The algorithm can find a near-optimal set of parameters to explain observed gauged data. The modeling framework presented here, although applied in a poorly-gauged catchment, can be adapted for catchments with more detailed observations. We applied the algorithm in different cases of the V-Tilted Catchment, the Wooden-Board catchment, and in an existing urban catchment with heterogeneous data. The results of automatic calibration indicate $\mathrm{NSE} = 0.99$ for the V-Tilted catchment, $\mathrm{RMSE} = 830~\mathrm{mgL^{-1}}$ for salt concentration pollutographs (i.e., 8.3% of the event mean concentration), and $\mathrm{NSE} = 0.89$ for the urban catchment case study. This paper also explores the issue of equifinality in modeling calibration (EqMC). Equifinality is defined as the set of different parameter combinations that can provide equally good or accepted results, within the physical parameter ranges. EqMC decreases with the number of events and increases with the choice of partially or nonproducing runoff ones. Furthermore, results indicate that providing more accurate parameter ranges based on a priori knowledge of the catchment is fundamental to reduce the chances of finding a set of parameters with equifinality.
△ Less
Submitted 31 August, 2023;
originally announced August 2023.
-
Non-linear top-Higgs CP violation
Authors:
Akanksha Bhardwaj,
Christoph Englert,
Dorival Gonçalves,
Alberto Navarro
Abstract:
Searches for additional sources of CP violation at the Large Hadron Collider are a central part of the Higgs physics programme beyond the Standard Model. Studies employing so-called signed observables that track CP violation through purpose-built asymmetries bolster efforts based on Higgs boson rate analyses under clear assumptions. A possibility, which is so far unexplored at the LHC, is a signif…
▽ More
Searches for additional sources of CP violation at the Large Hadron Collider are a central part of the Higgs physics programme beyond the Standard Model. Studies employing so-called signed observables that track CP violation through purpose-built asymmetries bolster efforts based on Higgs boson rate analyses under clear assumptions. A possibility, which is so far unexplored at the LHC, is a significant non-linear realisation of CP-violation, which is naturally described in non-linear Higgs Effective Field Theory (HEFT). We perform an analysis of the HL-LHC potential to constrain such interactions considering a large range of single and double Higgs production processes, including differential information where this is statistically and theoretically possible. A particular emphasis of our work is distinguishing expected correlations in the Standard Model Effective Field Theory from those attainable in HEFT.
△ Less
Submitted 22 August, 2023;
originally announced August 2023.
-
Universal stability of coherently diffusive 1D systems with respect to decoherence
Authors:
F. S. Lozano-Negro,
E. Alvarez Navarro,
N. C. Chávez,
F. Mattiotti,
F. Borgonovi,
H. M. Pastawski,
G. L. Celardo
Abstract:
Static disorder in a 3D crystal degrades the ideal ballistic dynamics until it produces a localized regime. This Metal-Insulator Transition is often preceded by coherent diffusion. By studying three paradigmatic 1D models, namely the Harper-Hofstadter-Aubry-André and Fibonacci tight-binding chains, along with the power-banded random matrix model, we show that whenever coherent diffusion is present…
▽ More
Static disorder in a 3D crystal degrades the ideal ballistic dynamics until it produces a localized regime. This Metal-Insulator Transition is often preceded by coherent diffusion. By studying three paradigmatic 1D models, namely the Harper-Hofstadter-Aubry-André and Fibonacci tight-binding chains, along with the power-banded random matrix model, we show that whenever coherent diffusion is present, transport is exceptionally stable against decoherent noise. This is completely at odds with what happens for coherently ballistic and localized dynamics, where the diffusion coefficient strongly depends on the environmental decoherence. A universal dependence of the diffusion coefficient on the decoherence strength is analytically derived: the diffusion coefficient remains almost decoherence-independent until the coherence time becomes comparable with the mean elastic scattering time.
Thus, systems with a quantum diffusive regime could be used to design robust quantum wires. Moreover our results might shed new light on the functionality of many biological systems, which often operate at the border between the ballistic and localized regimes.
△ Less
Submitted 27 March, 2024; v1 submitted 11 July, 2023;
originally announced July 2023.
-
RaViTT: Random Vision Transformer Tokens
Authors:
Felipe A. Quezada,
Carlos F. Navarro,
Cristian Muñoz,
Manuel Zamorano,
Jorge Jara-Wilde,
Violeta Chang,
Cristóbal A. Navarro,
Mauricio Cerda
Abstract:
Vision Transformers (ViTs) have successfully been applied to image classification problems where large annotated datasets are available. On the other hand, when fewer annotations are available, such as in biomedical applications, image augmentation techniques like introducing image variations or combinations have been proposed. However, regarding ViT patch sampling, less has been explored outside…
▽ More
Vision Transformers (ViTs) have successfully been applied to image classification problems where large annotated datasets are available. On the other hand, when fewer annotations are available, such as in biomedical applications, image augmentation techniques like introducing image variations or combinations have been proposed. However, regarding ViT patch sampling, less has been explored outside grid-based strategies. In this work, we propose Random Vision Transformer Tokens (RaViTT), a random patch sampling strategy that can be incorporated into existing ViTs. We experimentally evaluated RaViTT for image classification, comparing it with a baseline ViT and state-of-the-art (SOTA) augmentation techniques in 4 datasets, including ImageNet-1k and CIFAR-100. Results show that RaViTT increases the accuracy of the baseline in all datasets and outperforms the SOTA augmentation techniques in 3 out of 4 datasets by a significant margin +1.23% to +4.32%. Interestingly, RaViTT accuracy improvements can be achieved even with fewer tokens, thus reducing the computational load of any ViT model for a given accuracy value.
△ Less
Submitted 19 June, 2023;
originally announced June 2023.
-
Reassessing the impact of megaelectronvolt ions in fusion plasmas via gyrokinetic simulations
Authors:
Alessandro Di Siena,
Gabriele Merlo,
Alejandro Banon Navarro,
Tobias Goerler,
Frank Jenko
Abstract:
Gyrokinetic simulations conducted by Mazzi et al. reveal the suppression of turbulence in fusion plasmas through the destabilization of Toroidal Alfvén Eigenmodes (TAEs) by megaelectronvolt ions. Our analysis demonstrates that the authors' numerical findings are strongly influenced by the selected simulation settings, calling into question their claim that the resulting heat conductivity aligns wi…
▽ More
Gyrokinetic simulations conducted by Mazzi et al. reveal the suppression of turbulence in fusion plasmas through the destabilization of Toroidal Alfvén Eigenmodes (TAEs) by megaelectronvolt ions. Our analysis demonstrates that the authors' numerical findings are strongly influenced by the selected simulation settings, calling into question their claim that the resulting heat conductivity aligns with the TRANSP power balance. Specifically, we assert that the numerical results presented by Mazzi et al. are a direct consequence of the inadequacy of resolution employed in their numerical simulations. Notably, there are three primary factors contributing to this issue: (i) the employed radial box size is insufficient, leading to an undesirable impact of the boundary conditions on the simulations; (ii) the adoption of a higher minimum toroidal mode number fails to accurately resolve the entire range of TAEs, resulting in an underestimation of TAE drive and introducing spurious effects on wave-particle resonances; (iii) an insufficient resolution in the magnetic moment direction exacerbates these challenges. Upon addressing the aforementioned numerical issues, a significant increase in heat conductivity for each plasma species was observed by more than tenfold, diverging from the expected values derived from the TRANSP power balance calculations. Consequently, our results raise serious doubts on the authors' assertion of enhanced performance in the presence of strongly unstable TAEs, emphasizing the need for a reevaluation of their claims.
△ Less
Submitted 19 June, 2023;
originally announced June 2023.
-
Accelerating Range Minimum Queries with Ray Tracing Cores
Authors:
Enzo Meneses,
Cristóbal A. Navarro,
Héctor Ferrada,
Felipe A. Quezada
Abstract:
During the last decade GPU technology has shifted from pure general purpose computation to the inclusion of application specific integrated circuits (ASICs), such as Tensor Cores and Ray Tracing (RT) cores. Although these special purpose GPU cores were designed to further accelerate specific fields such as AI and real-time rendering, recent research has managed to exploit them to further accelerat…
▽ More
During the last decade GPU technology has shifted from pure general purpose computation to the inclusion of application specific integrated circuits (ASICs), such as Tensor Cores and Ray Tracing (RT) cores. Although these special purpose GPU cores were designed to further accelerate specific fields such as AI and real-time rendering, recent research has managed to exploit them to further accelerate other tasks that typically used regular GPU computing. In this work we present RTXRMQ, a new approach that can compute range minimum queries (RMQs) with RT cores. The main contribution is the proposal of a geometric solution for RMQ, where elements become triangles that are placed and shaped according to the element's value and position in the array, respectively, such that the closest hit of a ray launched from a point given by the query parameters corresponds to the result of that query. Experimental results show that RTXRMQ is currently best suited for small query ranges relative to the problem size, achieving up to $5\times$ and $2.3\times$ of speedup over state of the art CPU (HRMQ) and GPU (LCA) approaches, respectively. Although for medium and large query ranges RTXRMQ is currently surpassed by LCA, it is still competitive by being $2.5\times$ and $4\times$ faster than HRMQ which is a highly parallel CPU approach. Furthermore, performance scaling experiments across the latest RTX GPU architectures show that if the current RT scaling trend continues, then RTXRMQ's performance would scale at a higher rate than HRMQ and LCA, making the approach even more relevant for future high performance applications that employ batches of RMQs.
△ Less
Submitted 5 June, 2023;
originally announced June 2023.
-
When the Machine Chimes the Bell: Entanglement and Bell Inequalities with Boosted $t\bar{t}$
Authors:
Zhongtian Dong,
Dorival Gonçalves,
Kyoungchul Kong,
Alberto Navarro
Abstract:
The Large Hadron Collider provides a unique opportunity to study quantum entanglement and violation of Bell inequalities at the highest energy available today. In this paper, we will investigate these quantum correlations with top quark pair production, which represents a system of two-qubits. The spacelike separation requirement for the two causally disconnected top quarks requires they fly relat…
▽ More
The Large Hadron Collider provides a unique opportunity to study quantum entanglement and violation of Bell inequalities at the highest energy available today. In this paper, we will investigate these quantum correlations with top quark pair production, which represents a system of two-qubits. The spacelike separation requirement for the two causally disconnected top quarks requires they fly relativistically away from each other, which motivates the use of the boosted top-tagging with the semi-leptonic top pair channel. Although measuring the spin polarization of the hadronic top quark is known to be challenging, our study indicates that it is feasible to reconstruct the spin density matrix of the two-qubit system using an optimal hadronic polarimeter. This is achieved with the aid of jet substructure techniques and NN-inspired reconstruction methods, which improve the mapping between subjets and quarks. We find that entanglement can already be observed at more than $5σ$ level with existing data, and violation of Bell inequalities may be probed above 4$σ$ level at the HL-LHC with 3 ab$^{-1}$ of data.
△ Less
Submitted 5 July, 2024; v1 submitted 11 May, 2023;
originally announced May 2023.
-
An Evaluation of GPU Filters for Accelerating the 2D Convex Hull
Authors:
Roberto Carrasco,
Héctor Ferrada,
Cristóbal A. Navarro,
Nancy Hitschfeld
Abstract:
The Convex Hull algorithm is one of the most important algorithms in computational geometry, with many applications such as in computer graphics, robotics, and data mining. Despite the advances in the new algorithms in this area, it is often needed to improve the performance to solve more significant problems quickly or in real-time processing. This work presents an experimental evaluation of GPU…
▽ More
The Convex Hull algorithm is one of the most important algorithms in computational geometry, with many applications such as in computer graphics, robotics, and data mining. Despite the advances in the new algorithms in this area, it is often needed to improve the performance to solve more significant problems quickly or in real-time processing. This work presents an experimental evaluation of GPU filters to reduce the cost of computing the 2D convex hull. The technique first performs a preprocessing of the input set, filtering all points within an eight-vertex polygon in logarithmic time, to obtain a reduced set of candidate points. We use parallel computation and the use of the Manhattan distance as a metric to find the vertices of the polygon and perform the point filtering. For the filtering stage we study different approaches; from custom CUDA kernels to libraries such as Thrust and CUB. Three types of point distributions are tested: a normal distribution (favorable case), circumference (the worst case), and a case where points are shifted randomly from the circumference (intermediate case). Experimental evaluation shows that the GPU filtering algorithm can be up to 23x faster than a sequential CPU implementation, and the whole convex hull computation can be up to 30x faster than the fastest implementation provided by the CGAL library.
△ Less
Submitted 19 March, 2023;
originally announced March 2023.
-
HERMES: qualification of High pErformance pRogrammable Microprocessor and dEvelopment of Software ecosystem
Authors:
Nadia Ibellaatti,
Edouard Lepape,
Alp Kilic,
Kaya Akyel,
Kassem Chouayakh,
Fabrizio Ferrandi,
Claudio Barone,
Serena Curzel,
Michele Fiorito,
Giovanni Gozzi,
Miguel Masmano,
Ana Risquez Navarro,
Manuel Muñoz,
Vicente Nicolau Gallego,
Patricia Lopez Cueva,
Jean-noel Letrillard,
Franck Wartel
Abstract:
European efforts to boost competitiveness in the sector of space services promote the research and development of advanced software and hardware solutions. The EU-funded HERMES project contributes to the effort by qualifying radiation-hardened, high-performance programmable microprocessors, and by developing a software ecosystem that facilitates the deployment of complex applications on such platf…
▽ More
European efforts to boost competitiveness in the sector of space services promote the research and development of advanced software and hardware solutions. The EU-funded HERMES project contributes to the effort by qualifying radiation-hardened, high-performance programmable microprocessors, and by developing a software ecosystem that facilitates the deployment of complex applications on such platforms. The main objectives of the project include reaching a technology readiness level of 6 (i.e., validated and demonstrated in relevant environment) for the rad-hard NG-ULTRA FPGA with its ceramic hermetic package CGA 1752, developed within projects of the European Space Agency, French National Centre for Space Studies and the European Union. An equally important share of the project is dedicated to the development and validation of tools that support multicore software programming and FPGA acceleration, including Bambu for High-Level Synthesis and the XtratuM hypervisor with a level one boot loader for virtualization.
△ Less
Submitted 9 February, 2023;
originally announced February 2023.
-
A Learning Model Applied to the Calculation of the Velocities of 20 Stars Relative to the Sun
Authors:
Rafael Edgardo Carlos_Reyes,
Atilio Buendia Giribaldi,
Felipe Americo Reyes Navarro
Abstract:
We aim to explain the paradigm of a learning model, as well as to validate it in an applied case of an astronomy problem where the data used are declination, parallax, radial velocity of a star, as well as its annual variation in right ascension and declination. This study is based on a socio critical and positivist paradigm in the context of basic and applied science; algorithms and astronomical…
▽ More
We aim to explain the paradigm of a learning model, as well as to validate it in an applied case of an astronomy problem where the data used are declination, parallax, radial velocity of a star, as well as its annual variation in right ascension and declination. This study is based on a socio critical and positivist paradigm in the context of basic and applied science; algorithms and astronomical models were used as an instrument, which allowed us to address a specific case such as the calculation of the velocity of a star relative to the Sun.
△ Less
Submitted 29 December, 2022;
originally announced December 2022.
-
The Riemann-Roch theorem for the Adams operations
Authors:
A. Navarro,
J. Navarro
Abstract:
We prove the classical Riemann-Roch theorems for the Adams operations $\,ψ^j\,$ on $K$-theory: a statement with coefficients on $\mathbb{Z}[j^{-1}]$, that holds for arbitrary projective morphisms, as well as another one with integral coefficients, that is valid for closed immersions. In presence of rational coefficients, we also analyze the relation between the corresponding Riemann-Roch formula f…
▽ More
We prove the classical Riemann-Roch theorems for the Adams operations $\,ψ^j\,$ on $K$-theory: a statement with coefficients on $\mathbb{Z}[j^{-1}]$, that holds for arbitrary projective morphisms, as well as another one with integral coefficients, that is valid for closed immersions. In presence of rational coefficients, we also analyze the relation between the corresponding Riemann-Roch formula for one Adams operation and the analogous formula for the Chern character.
To do so, we complete the elementary exposition of the work of Panin-Smirnov that was initiated by the first author in a previous work. Their notion of oriented cohomology theory of algebraic varieties allows to use classical arguments to prove general and neat statements, which imply all the aforementioned results as particular cases.
△ Less
Submitted 18 November, 2022;
originally announced November 2022.
-
Predictions of improved confinement in SPARC via energetic particle turbulence stabilization
Authors:
A. Di Siena,
P. Rodriguez-Fernandez,
N. T. Howard,
A. Banon Navarro,
R. Bilato,
T. Goerler,
1 E. Poli,
G. Merlo,
J. Wrigh,
M. Greenwald,
F. Jenko
Abstract:
The recent progress in high-temperature superconductor technologies has led to the design and construction of SPARC, a compact tokamak device expected to reach plasma breakeven with up to $25$MW of external ion cyclotron resonant heating (ICRH) power. This manuscript presents local (flux-tube) and radially global gyrokinetic GENE (Jenko et al 2000 Phys. Plasmas {\bf 7} 1904) simulations for a redu…
▽ More
The recent progress in high-temperature superconductor technologies has led to the design and construction of SPARC, a compact tokamak device expected to reach plasma breakeven with up to $25$MW of external ion cyclotron resonant heating (ICRH) power. This manuscript presents local (flux-tube) and radially global gyrokinetic GENE (Jenko et al 2000 Phys. Plasmas {\bf 7} 1904) simulations for a reduced-field and current H-mode SPARC scenario showing that supra-thermal particles - generated via ICRH - strongly suppress ion-scale turbulent transport by triggering a fast ion-induced anomalous transport barrier (F-ATB). The trigger mechanism is identified as a wave-particle resonant interaction between the fast particle population and plasma micro-instabilities (Di Siena et al 2021 Phys. Rev. Lett. {\bf 125} 025002). By performing a series of global simulations employing different profiles for the thermal ions, we show that the fusion gain of this SPARC scenario could be substantially enhanced up to $\sim 80\%$ by exploiting this fast ion stabilizing mechanism. A study is also presented to further optimize the energetic particle profiles, thus possibly leading experimentally to an even more significant fusion gain.
△ Less
Submitted 21 October, 2022;
originally announced October 2022.
-
Instabilities and turbulence in stellarators from the perspective of global codes
Authors:
E. Sánchez,
A. Bañón Navarro,
F. Wilms,
M. Borchardt,
R. Kleiber,
F. Jenko
Abstract:
In this work, a comparison of the global gyrokinetic codes EUTERPE and GENE-3D in stellarator configurations of LHD and W7-X is carried out. In linear simulations with adiabatic electrons, excellent agreement is found in the mode numbers, growth rate and frequency, mode structure, and spatial localization of the most unstable mode in LHD. In W7-X, the dependence of the growth rate and frequency wi…
▽ More
In this work, a comparison of the global gyrokinetic codes EUTERPE and GENE-3D in stellarator configurations of LHD and W7-X is carried out. In linear simulations with adiabatic electrons, excellent agreement is found in the mode numbers, growth rate and frequency, mode structure, and spatial localization of the most unstable mode in LHD. In W7-X, the dependence of the growth rate and frequency with the mode number is well reproduced by both codes. The codes are also compared in linear simulations with kinetic ions and electrons in W7-X using model profiles, and reasonable agreement is found in the wavenumber of the most unstable modes. A stabilization of small-scale modes in kinetic-electron simulations with respect to the adiabatic-electron case is consistently found in both codes. Nonlinear simulations using adiabatic electrons and model profiles are also studied and the heat fluxes are compared. Very good agreement is found in the turbulent ion heat fluxes in both LHD and W7-X. Two problems that cannot be properly accounted for in local flux tube codes are studied: the localization of instabilities and turbulence over the flux surface and the influence of a background long-wavelength electric field. Good agreement between codes is found with respect to the spatial localization of instabilities and turbulence over the flux surface. The localization of saturated turbulence is found in both codes to be much smaller than that of the linear instabilities and smaller than previously reported in full-surface radially-local simulations. The influence of the electric field on the localization is also found to be smaller in the developed turbulent state than in the linear phase, and smaller than in previous works. A stabilizing effect of a constant electric field on the linearly unstable modes is found in both codes. A moderate reduction of turbulent transport by the radial electric field...
△ Less
Submitted 11 October, 2022;
originally announced October 2022.
-
First-principles based plasma profile predictions for optimized stellarators
Authors:
A. Bañón Navarro,
A. Di Siena,
J. L. Velasco,
F. Wilms,
G. Merlo,
T. Windisch,
L. L. LoDestro,
J. B. Parker,
F. Jenko
Abstract:
In the present Letter, first-of-its-kind computer simulations predicting plasma profiles for modern optimized stellarators -- while self-consistently retaining neoclassical transport, turbulent transport with 3D effects, and external physical sources -- are presented. These simulations exploit a newly developed coupling framework involving the global gyrokinetic turbulence code GENE-3D, the neocla…
▽ More
In the present Letter, first-of-its-kind computer simulations predicting plasma profiles for modern optimized stellarators -- while self-consistently retaining neoclassical transport, turbulent transport with 3D effects, and external physical sources -- are presented. These simulations exploit a newly developed coupling framework involving the global gyrokinetic turbulence code GENE-3D, the neoclassical transport code KNOSOS, and the 1D transport solver TANGO. This framework is used to analyze the recently observed degradation of energy confinement in electron-heated plasmas in the Wendelstein 7-X stellarator, where the central ion temperature was "clamped" to $T_i \approx 1.5$ keV regardless of the external heating power. By performing first-principles based simulations, we provide key evidence to understand this effect, namely the inefficient thermal coupling between electrons and ions in a turbulence-dominated regime, which is exacerbated by the large $T_e/T_i$ ratios, and show that a more efficient ion heat source, such as direct ion heating, will increase the on-axis ion temperature. This work paves the way towards the use of high-fidelity models for the development of the next generation of stellarators, in which neoclassical and turbulent transport are optimized simultaneously.
△ Less
Submitted 24 March, 2023; v1 submitted 4 October, 2022;
originally announced October 2022.
-
Accelerating the Convex Hull Computation with a Parallel GPU Algorithm
Authors:
Alan Keith,
Héctor Ferrada,
Cristóbal A. Navarro
Abstract:
The convex hull is a fundamental geometrical structure for many applications where groups of points must be enclosed or represented by a convex polygon. Although efficient sequential convex hull algorithms exist, and are constantly being used in applications, their computation time is often considered an issue for time-sensitive tasks such as real-time collision detection, clustering or image proc…
▽ More
The convex hull is a fundamental geometrical structure for many applications where groups of points must be enclosed or represented by a convex polygon. Although efficient sequential convex hull algorithms exist, and are constantly being used in applications, their computation time is often considered an issue for time-sensitive tasks such as real-time collision detection, clustering or image processing for virtual reality, among others, where fast response times are required. In this work we propose a parallel GPU-based adaptation of heaphull, which is a state of the art CPU algorithm that computes the convex hull by first doing a efficient filtering stage followed by the actual convex hull computation. More specifically, this work parallelizes the filtering stage, adapting it to the GPU programming model as a series of parallel reductions. Experimental evaluation shows that the proposed implementation significantly improves the performance of the convex hull computation, reaching up to $4\times$ of speedup over the sequential CPU-based heaphull and between $3\times \sim 4\times$ over existing GPU based approaches.
△ Less
Submitted 25 September, 2022;
originally announced September 2022.
-
GPU Voronoi Diagrams for Random Moving Seeds
Authors:
Rodrigo Stevenson,
Cristóbal A. Navarro
Abstract:
The Voronoi Diagram is a geometrical structure that is widely used in scientific or technological applications where proximity is a relevant aspect to consider, and it also resembles natural phenomena such as cellular banks, rock formations or bee hives, among others. Typically, computing the Voronoi Diagram is done in a static context, that is, the location of the input seeds is defined once and…
▽ More
The Voronoi Diagram is a geometrical structure that is widely used in scientific or technological applications where proximity is a relevant aspect to consider, and it also resembles natural phenomena such as cellular banks, rock formations or bee hives, among others. Typically, computing the Voronoi Diagram is done in a static context, that is, the location of the input seeds is defined once and does not change. In this work we study the dynamic case where seeds move, which leads to a dynamic Voronoi Diagram that changes over time. In particular, we consider uniform random moving seeds, for which we propose the \textit{dynamic Jump Flooding Algorithm} (dJFA), a variant of JFA that uses less iterations than the standard JFA. An experimental evaluation shows that dJFA achieves a speedup of up to $\sim 5.3 \times$ over JFA, while maintaining a similarity of at least $88\%$ and close to $100\%$ in many cases. These results contribute with a step towards the achievement of real-time GPU-based computation of dynamic Voronoi diagrams for any particle simulation.
△ Less
Submitted 31 August, 2022;
originally announced September 2022.
-
GGArray: A Dynamically Growable GPU Array
Authors:
Enzo Meneses,
Cristóbal A. Navarro,
Héctor Ferrada
Abstract:
We present a dynamically Growable GPU array (GGArray) fully implemented in GPU that does not require synchronization with the host. The idea is to improve the programming of GPU applications that require dynamic memory, by offering a structure that does not require pre-allocating GPU VRAM for the worst case scenario. The GGArray is based on the LFVector, by utilizing an array of them in order to t…
▽ More
We present a dynamically Growable GPU array (GGArray) fully implemented in GPU that does not require synchronization with the host. The idea is to improve the programming of GPU applications that require dynamic memory, by offering a structure that does not require pre-allocating GPU VRAM for the worst case scenario. The GGArray is based on the LFVector, by utilizing an array of them in order to take advantage of the GPU architecture and the synchronization offered by thread blocks. This structure is compared to other state of the art ones such as a pre-allocated static array and a semi-static array that needs to be resized through communication with the host. Experimental evaluation shows that the GGArray has a competitive insertion and resize performance, but it is slower for regular parallel memory accesses. Given the results, the GGArray is a potentially useful structure for applications with high uncertainty on the memory usage as well as applications that have phases, such as an insertion phase followed by a regular GPU phase. In such cases, the GGArray can be used for the first phase and then data can be flattened for the second phase in order to allow the classical GPU memory accesses which are faster. These results constitute a step towards achieving a parallel efficient C++ like vector for modern GPU architectures.
△ Less
Submitted 7 September, 2022; v1 submitted 31 August, 2022;
originally announced September 2022.
-
A Scalable and Energy Efficient GPU Thread Map for m-Simplex Domains
Authors:
Cristóbal A. Navarro,
Felipe A. Quezada,
Benjamin Bustos,
Nancy Hitschfeld,
Rolando Kindelan
Abstract:
This work proposes a new GPU thread map for $m$-simplex domains, that scales its speedup with dimension and is energy efficient compared to other state of the art approaches. The main contributions of this work are i) the formulation of the new block-space map $\mathcal{H}: \mathbb{Z}^m \mapsto \mathbb{Z}^m$ for regular orthogonal simplex domains, which is analyzed in terms of resource usage, and…
▽ More
This work proposes a new GPU thread map for $m$-simplex domains, that scales its speedup with dimension and is energy efficient compared to other state of the art approaches. The main contributions of this work are i) the formulation of the new block-space map $\mathcal{H}: \mathbb{Z}^m \mapsto \mathbb{Z}^m$ for regular orthogonal simplex domains, which is analyzed in terms of resource usage, and ii) the experimental evaluation in terms of speedup over a bounding box approach and energy efficiency as elements per second per Watt. Results from the analysis show that $\mathcal{H}$ has a potential speedup of up to $2\times$ and $6\times$ for $2$ and $3$-simplices, respectively. Experimental evaluation shows that $\mathcal{H}$ is competitive for $2$-simplices, reaching $1.2\times \sim 2.0\times$ of speedup for different tests, which is on par with the fastest state of the art approaches. For $3$-simplices $\mathcal{H}$ reaches up to $1.3\times \sim 6.0\times$ of speedup making it the fastest of all. The extension of $\mathcal{H}$ to higher dimensional $m$-simplices is feasible and has a potential speedup that scales as $m!$ given a proper selection of parameters $r, β$ which are the scaling and replication factors, respectively. In terms of energy consumption, although $\mathcal{H}$ is among the highest in power consumption, it compensates by its short duration, making it one of the most energy efficient approaches. Lastly, further improvements with Tensor and Ray Tracing Cores are analyzed, giving insights to leverage each one of them. The results obtained in this work show that $\mathcal{H}$ is a scalable and energy efficient map that can contribute to the efficiency of GPU applications when they need to process $m$-simplex domains, such as Cellular Automata or PDE simulations.
△ Less
Submitted 12 September, 2022; v1 submitted 24 August, 2022;
originally announced August 2022.
-
Modeling GPU Dynamic Parallelism for Self Similar Density Workloads
Authors:
Felipe A. Quezada,
Cristóbal A. Navarro,
Miguel Romero,
Cristhian Aguilera
Abstract:
Dynamic Parallelism (DP) is a runtime feature of the GPU programming model that allows GPU threads to execute additional GPU kernels, recursively. Apart from making the programming of parallel hierarchical patterns easier, DP can also speedup problems that exhibit a heterogeneous data layout by focusing, through a subdivision process, the finite GPU resources on the sub-regions that exhibit more p…
▽ More
Dynamic Parallelism (DP) is a runtime feature of the GPU programming model that allows GPU threads to execute additional GPU kernels, recursively. Apart from making the programming of parallel hierarchical patterns easier, DP can also speedup problems that exhibit a heterogeneous data layout by focusing, through a subdivision process, the finite GPU resources on the sub-regions that exhibit more parallelism. However, doing an optimal subdivision process is not trivial, as there are different parameters that play an important role in the final performance of DP. Moreover, the current programming abstraction for DP also introduces an overhead that can penalize the final performance. In this work we present a subdivision cost model for problems that exhibit self similar density (SSD) workloads (such as fractals), in order understand what parameters provide the fastest subdivision approach. Also, we introduce a new subdivision implementation, named \textit{Adaptive Serial Kernels} (ASK), as a smaller overhead alternative to CUDA's Dynamic Parallelism. Using the cost model on the Mandelbrot Set as a case study shows that the optimal scheme is to start with an initial subdivision between $g=[2,16]$, then keep subdividing in regions of $r=2,4$, and stop when regions reach a size of $B \sim 32$. The experimental results agree with the theoretical parameters, confirming the usability of the cost model. In terms of performance, the proposed ASK approach runs up to $\sim 60\%$ faster than Dynamic Parallelism in the Mandelbrot set, and up to $12\times$ faster than a basic exhaustive implementation, whereas DP is up to $7.5\times$.
△ Less
Submitted 5 June, 2022;
originally announced June 2022.
-
Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D
Authors:
Anamaría Navarro,
E. Khomenko,
M. Modestov,
N. Vitas
Abstract:
Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper con…
▽ More
Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper convection zone into the solar corona. Two different schemes to model heat conduction are implemented: (1) a standard scheme where a parabolic term is added to the energy equation, and (2) a scheme where the hyperbolic heat flux equation is solved. The first scheme limits the time step due to the explicit integration of a parabolic term, which makes the simulations computationally expensive. The second scheme solves the limitations on the time step by artificially limiting the heat conduction speed to computationally manageable values. The validation of both schemes is carried out with standard tests in one, two, and three spatial dimensions. Furthermore, we implement the model for heat flux derived by Braginskii (1965) in its most general form, when the expression for the heat flux depends on the ratio of the collisional to cyclotron frequencies of the plasma, and, therefore on the magnetic field strength. Additionally, our implementation takes into account the heat conduction in parallel, perpendicular, and transverse directions, and provides the contributions from ions and electrons separately. The model also transitions smoothly between field-aligned conductivity and isotropic conductivity for regions with a low or null magnetic field. Finally, we present a two-dimensional test for heat conduction using realistic values of the solar atmosphere where we prove the robustness of the two schemes implemented.
△ Less
Submitted 20 May, 2022; v1 submitted 18 May, 2022;
originally announced May 2022.
-
Global gyrokinetic simulations of ASDEX Upgrade up to the transport time-scale with GENE-Tango
Authors:
A. Di Siena,
A. Banon Navarro,
T. Luda,
G. Merlo,
M. Bergmann,
L. Leppin,
T. Goerler,
J. B. Parker,
L. LoDestro,
J. Hittinger,
B. Dorland,
G. Hammett,
F. Jenko,
the ASDEX Upgrade Team,
the EUROfusion MST1 Team
Abstract:
An accurate description of turbulence up to the transport time scale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport time scales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collision…
▽ More
An accurate description of turbulence up to the transport time scale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport time scales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collisions, realistic geometries, toroidal rotation and electromagnetic effects for the first time. This approach overcomes gyrokinetic codes' limitations and enables high-fidelity profile calculations in experimentally relevant plasma conditions, significantly reducing the computational cost.
We present numerical results of GENE-Tango for two ASDEX Upgrade discharges, one of which exhibits a pronounced peaking of the ion temperature profile not reproduced by TGLF-ASTRA. We show that GENE-Tango can correctly capture the ion temperature peaking observed in the experiment. By retaining different physical effects in the GENE simulations, e.g., collisions, toroidal rotation and electromagnetic effects, we demonstrate that the ion temperature profile's peaking is due to electromagnetic effects of submarginal MHD instability. Based on these results, the expected GENE-Tango speedup for the ITER standard scenario is larger than two orders of magnitude compared to a single gyrokinetic simulation up to the transport time scale, possibly making first-principles ITER simulations feasible on current computing resources.
△ Less
Submitted 12 April, 2022;
originally announced April 2022.
-
Some results for semi-stable radial solutions of $k$-Hessian equations
Authors:
Miguel Angel Navarro,
Justino Sanchez
Abstract:
We devote this paper to study semi-stable nonconstant radial solutions of $S_k(D^2u) = w(|x|)g(u)$ on the Euclidean space $R^n$. We establish pointwise estimates and necessary conditions for the existence of such solutions (not necessarily bounded) for this equation. For bounded solutions we estimate their asymptotic behavior at infinity. All the estimates are given in terms of the spatial dimensi…
▽ More
We devote this paper to study semi-stable nonconstant radial solutions of $S_k(D^2u) = w(|x|)g(u)$ on the Euclidean space $R^n$. We establish pointwise estimates and necessary conditions for the existence of such solutions (not necessarily bounded) for this equation. For bounded solutions we estimate their asymptotic behavior at infinity. All the estimates are given in terms of the spatial dimension $n$, the values of $k$ and the behavior at infinity of the growth rate function of $w$.
△ Less
Submitted 18 February, 2022;
originally announced February 2022.
-
Squeeze: Efficient Compact Fractals for Tensor Core GPUs
Authors:
Felipe A. Quezada,
Cristóbal A. Navarro,
Nancy Hitschfeld,
Benjamin Bustos
Abstract:
This work presents Squeeze, an efficient compact fractal processing scheme for tensor core GPUs. By combining discrete-space transformations between compact and expanded forms, one can do data-parallel computation on a fractal with neighborhood access without needing to expand the fractal in memory. The space transformations are formulated as two GPU tensor-core accelerated thread maps, $λ(ω)$ and…
▽ More
This work presents Squeeze, an efficient compact fractal processing scheme for tensor core GPUs. By combining discrete-space transformations between compact and expanded forms, one can do data-parallel computation on a fractal with neighborhood access without needing to expand the fractal in memory. The space transformations are formulated as two GPU tensor-core accelerated thread maps, $λ(ω)$ and $ν(ω)$, which act as compact-to-expanded and expanded-to-compact space functions, respectively. The cost of the maps is $\mathcal{O}(\log_2 \log_s(n))$ time, with $n$ being the side of a $n \times n$ embedding for the fractal in its expanded form, and $s$ the linear scaling factor. The proposed approach works for any fractal that belongs to the Non-overlapping-Bounding-Boxes (NBB) class of discrete fractals, and can be extended to three dimensions as well. Experimental results using a discrete Sierpinski Triangle as a case study shows up to $\sim12\times$ of speedup and a memory reduction factor of up to $\sim 315\times$ with respect to a GPU-based expanded-space bounding box approach. These results show that the proposed compact approach will allow the scientific community to efficiently tackle problems that up to now could not fit into GPU memory.
△ Less
Submitted 3 January, 2022;
originally announced January 2022.
-
Accelerating Compact Fractals with Tensor Core GPUs
Authors:
Felipe A. Quezada,
Cristóbal A. Navarro
Abstract:
This work presents a GPU thread mapping approach that allows doing fast parallel stencil-like computations on discrete fractals using their compact representation. The intuition behind is to employ two GPU tensor-core accelerated thread maps, $λ(ω)$ and $ν(ω)$, which act as threadspace-to-dataspace and dataspace-to-threadspace functions, respectively. By combining these maps, threads can access co…
▽ More
This work presents a GPU thread mapping approach that allows doing fast parallel stencil-like computations on discrete fractals using their compact representation. The intuition behind is to employ two GPU tensor-core accelerated thread maps, $λ(ω)$ and $ν(ω)$, which act as threadspace-to-dataspace and dataspace-to-threadspace functions, respectively. By combining these maps, threads can access compact space and interact with their neighbors. The cost of the maps is $\mathcal{O}(\log \log(n))$ time, with $n$ being the side of a $n \times n$ embedding for the fractal in its expanded form. The technique works on any fractal that belongs to the Non-overlapping-Bounding-Boxes (NBB) class of discrete fractals, and can be extended to three dimensions as well. Results using an A100 GPU on the Sierpinski Triangle as a case study show up to $\sim11\times$ of speedup and a memory usage reduction of $234\times$ with respect to a Bounding Box approach. These results show that the proposed compact approach can allow the scientific community to tackle larger problems that did not fit in GPU memory before, and run even faster than a bounding box approach.
△ Less
Submitted 25 October, 2021;
originally announced October 2021.
-
Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps
Authors:
M. Gatti,
B. Jain,
C. Chang,
M. Raveri,
D. Zürcher,
L. Secco,
L. Whiteway,
N. Jeffrey,
C. Doux,
T. Kacprzak,
D. Bacon,
P. Fosalba,
A. Alarcon,
A. Amon,
K. Bechtol,
M. Becker,
G. Bernstein,
J. Blazek,
A. Campos,
A. Choi,
C. Davis,
J. Derose,
S. Dodelson,
F. Elsner,
J. Elvin-Poole
, et al. (85 additional authors not shown)
Abstract:
We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing…
▽ More
We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analysed in the context of the $Λ$CDM model, varying 5 cosmological parameters and 19 nuisance parameters modelling astrophysical and measurement systematics. Our modelling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7\% measurement of the amplitude of fluctuations parameter $S_8\equiv σ_8 (Ω_m/0.3)^{0.5} = 0.784\pm 0.013$. The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the $Λ$CDM model. The addition of the third moment improves the constraints on $S_8$ and $Ω_{\rm m}$ by $\sim$15\% and $\sim$25\% compared to an analysis that only uses second moments. We compare our results with {\it Planck} constraints from the Cosmic Microwave Background (CMB), finding a $2.2$ \textendash $2.8σ$ tension in the full parameter space, depending on the combination of moments considered. The third moment independently is in $2.8σ$ tension with {\it Planck}, and thus provides a cross-check on analyses of 2-point correlations.
△ Less
Submitted 9 September, 2022; v1 submitted 19 October, 2021;
originally announced October 2021.
-
Core transport barriers induced by fast ions in global gyrokinetic GENE simulations
Authors:
A. Di Siena,
R. Bilato,
T. Goerler,
E. Poli,
A. Banon Navarro,
D. Jarema,
F. Jenko
Abstract:
A novel type of internal transport barrier (ITB) called F-ATB (fast ion-induced anomalous transport barrier) has been recently observed in state-of-the-art global gyrokinetic simulations on a properly optimized ASDEX Upgrade experiment [A. Di Siena et al. Phys. Rev. Lett. {\bf 127} 025002 (2021)]. Unlike the transport barriers previously reported in literature, the trigger mechanism for the F-ATB…
▽ More
A novel type of internal transport barrier (ITB) called F-ATB (fast ion-induced anomalous transport barrier) has been recently observed in state-of-the-art global gyrokinetic simulations on a properly optimized ASDEX Upgrade experiment [A. Di Siena et al. Phys. Rev. Lett. {\bf 127} 025002 (2021)]. Unlike the transport barriers previously reported in literature, the trigger mechanism for the F-ATB is a basically electrostatic wave-particle resonant interaction between supra-thermal particles - generated via ion cyclotron resonance heating (ICRH) - and ion scale plasma turbulence. This resonant effect strongly depends on the particular shape of the fast ion temperature and density profiles. Therefore, to further improve our theoretical understanding of this transport barrier, we present results exploring the parameter space and physical conditions for the F-ATB generation by performing a systematic study with global GENE simulations. Particular emphasis is given to the transport barrier width and its localization by scanning over different energetic particle temperature profiles. The latter are varied in amplitude, half-width, and radial localization of an ad-hoc Gaussian-like energetic particle logarithmic temperature gradient profile. For the reference parameters at hand, a threshold in the amplitude of the fast ion logarithmic temperature gradient is identified to trigger the transport barrier effectively.
△ Less
Submitted 18 October, 2021;
originally announced October 2021.
-
Upgrade of the CMS Resistive Plate Chambers for the High Luminosity LHC
Authors:
A. Samalan,
M. Tytgat,
G. A. Alves,
F. Marujo,
F. Torres Da Silva De Araujo,
E. M. DaCosta,
D. De Jesus Damiao,
H. Nogima,
A. Santoro,
S. Fonseca De Souza,
A. Aleksandrov,
R. Hadjiiska,
P. Iaydjiev,
M. Rodozov,
M. Shopova,
G. Soultanov,
M. Bonchev,
A. Dimitrov,
L. Litov,
B. Pavlov,
P. Petkov,
A. Petrov,
S. J. Qian,
C. Bernal,
A. Cabrera
, et al. (86 additional authors not shown)
Abstract:
During the upcoming High Luminosity phase of the Large Hadron Collider (HL-LHC), the integrated luminosity of the accelerator will increase to 3000 fb$^{-1}$. The expected experimental conditions in that period in terms of background rates, event pileup, and the probable aging of the current detectors present a challenge for all the existing experiments at the LHC, including the Compact Muon Solen…
▽ More
During the upcoming High Luminosity phase of the Large Hadron Collider (HL-LHC), the integrated luminosity of the accelerator will increase to 3000 fb$^{-1}$. The expected experimental conditions in that period in terms of background rates, event pileup, and the probable aging of the current detectors present a challenge for all the existing experiments at the LHC, including the Compact Muon Solenoid (CMS) experiment. To ensure a highly performing muon system for this period, several upgrades of the Resistive Plate Chamber (RPC) system of the CMS are currently being implemented. These include the replacement of the readout system for the present system, and the installation of two new RPC stations with improved chamber and front-end electronics designs. The current overall status of this CMS RPC upgrade project is presented.
△ Less
Submitted 2 November, 2021; v1 submitted 29 September, 2021;
originally announced September 2021.
-
Constant-roll inflation in the generalized SU(2) Proca theory
Authors:
Juan C. Garnica,
L. Gabriel Gomez,
Andres A. Navarro,
Yeinzon Rodriguez
Abstract:
The generalized SU(2) Proca theory (GSU2P) is a variant of the well known generalized Proca theory where the vector field belongs to the Lie algebra of the SU(2) group of global transformations under which the action is made invariant. New interesting possibilities arise in this framework because of the existence of new interactions of purely non-Abelian character and new configurations of the vec…
▽ More
The generalized SU(2) Proca theory (GSU2P) is a variant of the well known generalized Proca theory where the vector field belongs to the Lie algebra of the SU(2) group of global transformations under which the action is made invariant. New interesting possibilities arise in this framework because of the existence of new interactions of purely non-Abelian character and new configurations of the vector field resulting in spatial spherical symmetry and the cosmological dynamics being driven by the propagating degrees of freedom. We study the two-dimensional phase space of the system that results when the cosmic triad configuration is employed in the Friedmann-Lemaitre-Robertson-Walker background and find an attractor curve whose attraction basin both covers almost all the allowed region and does not include a Big-Bang singularity. Such an attractor curve corresponds to a primordial inflationary solution that has the following characteristic properties: 1.) it is a de Sitter solution whose Hubble parameter is regulated by a generalized version of the SU(2) group coupling constant, 2.) it is constant-roll including, as a limiting case, the slow-roll variety, 3.) a number of e-folds $N > 60$ is easily reached, 4.) it has a graceful exit into a radiation dominated period powered by the canonical kinetic term of the vector field and the Einstein-Hilbert term. The free parameters of the action are chosen such that the tensor sector of the theory is the same as that of general relativity at least up to second-order perturbations, thereby avoiding the presence of ghost and Laplacian instabilities in the tensor sector as well as making the gravity waves propagate at light speed. This is a proof of concept of the interesting properties we could find in this scenario when the coupling constants be replaced by general coupling functions and more terms be discovered in the GSU2P.
△ Less
Submitted 6 January, 2022; v1 submitted 21 September, 2021;
originally announced September 2021.
-
Gyrokinetic simulations in stellarators using different computational domains
Authors:
E. Sánchez,
J. M. García-Regaña,
A. Bañón Navarro,
J. H. E. Proll,
C. Mora Moreno,
A. González-Jerez,
I. Calvo,
R. Kleiber,
J. Riemann,
J. Smoniewski,
M. Barnes,
F. I. Parra
Abstract:
In this work, we compare gyrokinetic simulations in stellarators using different computational domains, namely, flux tube, full-flux-surface, and radially global domains. Two problems are studied: the linear relaxation of zonal flows and the linear stability of ion temperature gradient (ITG) modes. Simulations are carried out with the codes EUTERPE, GENE, GENE-3D, and stella in magnetic configurat…
▽ More
In this work, we compare gyrokinetic simulations in stellarators using different computational domains, namely, flux tube, full-flux-surface, and radially global domains. Two problems are studied: the linear relaxation of zonal flows and the linear stability of ion temperature gradient (ITG) modes. Simulations are carried out with the codes EUTERPE, GENE, GENE-3D, and stella in magnetic configurations of LHD and W7-X using adiabatic electrons. The zonal flow relaxation properties obtained in different flux tubes are found to differ with each other and with the radially global result, except for sufficiently long flux tubes, in general. The flux tube length required for convergence is configuration-dependent. Similarly, for ITG instabilities, different flux tubes provide different results, but the discrepancy between them diminishes with increasing flux tube length. Full-flux-surface and flux tube simulations show good agreement in the calculation of the growth rate and frequency of the most unstable modes in LHD, while for W7-X differences in the growth rates are found between the flux tube and the full-flux-surface domains. Radially global simulations provide results close to the full-flux-surface ones. The radial scale of unstable ITG modes is studied in global and flux tube simulations finding that in W7-X, the radial scale of the most unstable modes depends on the binormal wavenumber, while in LHD no clear dependency is found.
△ Less
Submitted 5 June, 2021;
originally announced June 2021.
-
Analysis of zero inflated dichotomous variables from a Bayesian perspective: Application to occupational health
Authors:
David Moriña,
Pedro Puig,
Albert Navarro
Abstract:
This work proposes a new methodology to fit zero inflated Bernoulli data from a Bayesian approach, able to distinguish between two potential sources of zeros (structurals and non-structurals). Its usage is illustrated by means of a real example from the field of occupational health as the phenomenon of sickness presenteeism, in which it is reasonable to think that some individuals will never be at…
▽ More
This work proposes a new methodology to fit zero inflated Bernoulli data from a Bayesian approach, able to distinguish between two potential sources of zeros (structurals and non-structurals). Its usage is illustrated by means of a real example from the field of occupational health as the phenomenon of sickness presenteeism, in which it is reasonable to think that some individuals will never be at risk of suffering it because they have not been sick in the period of study (structural zeros). Without separating structural and non-structural zeros one would one would be studying jointly the general health status and the presenteeism itself, and therefore obtaining potentially biased estimates as the phenomenon is being implicitly underestimated by diluting it into the general health status. The proposed methodology performance has been evaluated through a comprehensive simulation study, and it has been compiled as an R package freely available to the community.
△ Less
Submitted 3 May, 2021;
originally announced May 2021.