-
FEBio FINESSE: An open-source finite element simulation approach to estimate in vivo heart valve strains using shape enforcement
Authors:
Devin W. Laurence,
Patricia M. Sabin,
Analise M. Sulentic,
Matthew Daemer,
Steve A. Maas,
Jeffrey A. Weiss,
Matthew A. Jolley
Abstract:
Finite element simulations are an enticing tool to evaluate heart valve function in healthy and diseased patients; however, patient-specific simulations derived from 3D echocardiography are hampered by several technical challenges. In this work, we present an open-source method to enforce matching between finite element simulations and in vivo image-derived heart valve geometry in the absence of p…
▽ More
Finite element simulations are an enticing tool to evaluate heart valve function in healthy and diseased patients; however, patient-specific simulations derived from 3D echocardiography are hampered by several technical challenges. In this work, we present an open-source method to enforce matching between finite element simulations and in vivo image-derived heart valve geometry in the absence of patient-specific material properties, leaflet thickness, and chordae tendineae structures. We evaluate FEBio Finite Element Simulations with Shape Enforcement (FINESSE) using three synthetic test cases covering a wide range of model complexity. Our results suggest that FINESSE can be used to not only enforce finite element simulations to match an image-derived surface, but to also estimate the first principal leaflet strains within +/- 0.03 strain. Key FINESSE considerations include: (i) appropriately defining the user-defined penalty, (ii) omitting the leaflet commissures to improve simulation convergence, and (iii) emulating the chordae tendineae behavior via prescribed leaflet free edge motion or a chordae emulating force. We then use FINESSE to estimate the in vivo valve behavior and leaflet strains for three pediatric patients. In all three cases, FINESSE successfully matched the target surface with median errors similar to or less than the smallest voxel dimension. Further analysis revealed valve-specific findings, such as the tricuspid valve leaflet strains of a 2-day old patient with HLHS being larger than those of two 13-year old patients. The development of this open source pipeline will enable future studies to begin linking in vivo leaflet mechanics with patient outcomes
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
A Perspective on Foundation Models for the Electric Power Grid
Authors:
Hendrik F. Hamann,
Thomas Brunschwiler,
Blazhe Gjorgiev,
Leonardo S. A. Martins,
Alban Puech,
Anna Varbella,
Jonas Weiss,
Juan Bernabe-Moreno,
Alexandre Blondin Massé,
Seong Choi,
Ian Foster,
Bri-Mathias Hodge,
Rishabh Jain,
Kibaek Kim,
Vincent Mai,
François Mirallès,
Martin De Montigny,
Octavio Ramos-Leaños,
Hussein Suprême,
Le Xie,
El-Nasser S. Youssef,
Arnaud Zinflou,
Alexander J. Belvi,
Ricardo J. Bessa,
Bishnu Prasad Bhattari
, et al. (2 additional authors not shown)
Abstract:
Foundation models (FMs) currently dominate news headlines. They employ advanced deep learning architectures to extract structural information autonomously from vast datasets through self-supervision. The resulting rich representations of complex systems and dynamics can be applied to many downstream applications. Therefore, FMs can find uses in electric power grids, challenged by the energy transi…
▽ More
Foundation models (FMs) currently dominate news headlines. They employ advanced deep learning architectures to extract structural information autonomously from vast datasets through self-supervision. The resulting rich representations of complex systems and dynamics can be applied to many downstream applications. Therefore, FMs can find uses in electric power grids, challenged by the energy transition and climate change. In this paper, we call for the development of, and state why we believe in, the potential of FMs for electric grids. We highlight their strengths and weaknesses amidst the challenges of a changing grid. We argue that an FM learning from diverse grid data and topologies could unlock transformative capabilities, pioneering a new approach in leveraging AI to redefine how we manage complexity and uncertainty in the electric grid. Finally, we discuss a power grid FM concept, namely GridFM, based on graph neural networks and show how different downstream tasks benefit.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Distorted Magnetic Flux Ropes within Interplanetary Coronal Mass Ejections
Authors:
Andreas J. Weiss,
Teresa Nieves-Chinchilla,
Christian Möstl
Abstract:
Magnetic flux ropes within interplanetary coronal mass ejections are often characterized as simplistic cylindrical or toroidal tubes with field lines that twist around the cylinder or torus axis. Recent multi-point observations suggest that the overall geometry of these large-scale structures may be significantly more complex, so that the contemporary modeling approaches would be, in some cases, i…
▽ More
Magnetic flux ropes within interplanetary coronal mass ejections are often characterized as simplistic cylindrical or toroidal tubes with field lines that twist around the cylinder or torus axis. Recent multi-point observations suggest that the overall geometry of these large-scale structures may be significantly more complex, so that the contemporary modeling approaches would be, in some cases, insufficient to properly understand the global structure of any interplanetary coronal mass ejection. In an attempt to partially rectify this issue, we have developed a novel magnetic flux rope model that allows for the description of arbitrary distortions of the cross-section or deformation of the magnetic axis. The distorted magnetic flux rope model is a fully analytic flux rope model, that can be used to describe significantly more complex geometries and is numerically efficient enough to be used for large ensemble simulations. To demonstrate the usefulness of our model, we focus on a specific implementation of our model and apply it to an ICME event that was observed \textit{in situ} on 2023 April 23 at the L1 point by the Wind spacecraft and also by the STEREO-A spacecraft that was $10.2^\circ$ further east and $0.9^\circ$ south in heliographic coordinates. We demonstrate that our model can accurately reconstruct each observation individual and also gives a fair reconstruction of both events simultaneously using a multi-point reconstruction approach, which results in a geometry that is not fully constistent with a cylindrical or toroidal approximation.
△ Less
Submitted 18 June, 2024;
originally announced June 2024.
-
Automatic detection of large-scale flux ropes and their geoeffectiveness with a machine learning approach
Authors:
Sanchita Pal,
Luiz F. G. dos Santos,
Andreas J. Weiss,
Thomas Narock,
Ayris Narock,
Teresa Nieves-Chinchilla,
Lan K. Jian,
Simon W. Good
Abstract:
Detecting large-scale flux ropes (FRs) embedded in interplanetary coronal mass ejections (ICMEs) and assessing their geoeffectiveness are essential since they can drive severe space weather. At 1 au, these FRs have an average duration of 1 day. Their most common magnetic features are large, smoothly rotating magnetic fields. Their manual detection has become a relatively common practice over decad…
▽ More
Detecting large-scale flux ropes (FRs) embedded in interplanetary coronal mass ejections (ICMEs) and assessing their geoeffectiveness are essential since they can drive severe space weather. At 1 au, these FRs have an average duration of 1 day. Their most common magnetic features are large, smoothly rotating magnetic fields. Their manual detection has become a relatively common practice over decades, although visual detection can be time-consuming and subject to observer bias. Our study proposes a pipeline that utilizes two supervised binary-classification machine learning (ML) models trained with solar wind magnetic properties to automatically detect large-scale FRs and additionally determine their geoeffectiveness. The first model is used to generate a list of auto-detected FRs. Using the properties of southward magnetic field the second model determines the geoeffectiveness of FRs. Our method identifies 88.6\% and 80\% large-scale ICMEs (duration $\ge 1$ day) observed at 1 au by Wind and Sun Earth Connection Coronal and Heliospheric Investigation (STEREO) mission, respectively. While testing with a continuous solar wind data obtained from Wind, our pipeline detected 56 of the 64 large-scale ICMEs during 2008 - 2014 period (recall= 0.875) but many false positives (precision= 0.56) as we do not take into account any additional solar wind properties than the magnetic properties. We found an accuracy of 0.88 when estimating the geoeffectiveness of the auto-detected FRs using our method. Thus, in space weather now-casting and forecasting at L1 or any planetary missions, our pipeline can be utilized to offer a first-order detection of large-scale FRs and geoeffectiveness.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
TotalVibeSegmentator: Full Torso Segmentation for the NAKO and UK Biobank in Volumetric Interpolated Breath-hold Examination Body Images
Authors:
Robert Graf,
Paul-Sören Platzek,
Evamaria Olga Riedel,
Constanze Ramschütz,
Sophie Starck,
Hendrik Kristian Möller,
Matan Atad,
Henry Völzke,
Robin Bülow,
Carsten Oliver Schmidt,
Julia Rüdebusch,
Matthias Jung,
Marco Reisert,
Jakob Weiss,
Maximilian Löffler,
Fabian Bamberg,
Bene Wiestler,
Johannes C. Paetzold,
Daniel Rueckert,
Jan Stefan Kirschke
Abstract:
Objectives: To present a publicly available torso segmentation network for large epidemiology datasets on volumetric interpolated breath-hold examination (VIBE) images. Materials & Methods: We extracted preliminary segmentations from TotalSegmentator, spine, and body composition networks for VIBE images, then improved them iteratively and retrained a nnUNet network. Using subsets of NAKO (85 subje…
▽ More
Objectives: To present a publicly available torso segmentation network for large epidemiology datasets on volumetric interpolated breath-hold examination (VIBE) images. Materials & Methods: We extracted preliminary segmentations from TotalSegmentator, spine, and body composition networks for VIBE images, then improved them iteratively and retrained a nnUNet network. Using subsets of NAKO (85 subjects) and UK Biobank (16 subjects), we evaluated with Dice-score on a holdout set (12 subjects) and existing organ segmentation approach (1000 subjects), generating 71 semantic segmentation types for VIBE images. We provide an additional network for the vertebra segments 22 individual vertebra types. Results: We achieved an average Dice score of 0.89 +- 0.07 overall 71 segmentation labels. We scored > 0.90 Dice-score on the abdominal organs except for the pancreas with a Dice of 0.70. Conclusion: Our work offers a detailed and refined publicly available full torso segmentation on VIBE images.
△ Less
Submitted 31 May, 2024;
originally announced June 2024.
-
TotalSegmentator MRI: Sequence-Independent Segmentation of 59 Anatomical Structures in MR images
Authors:
Tugba Akinci D'Antonoli,
Lucas K. Berger,
Ashraya K. Indrakanti,
Nathan Vishwanathan,
Jakob Weiß,
Matthias Jung,
Zeynep Berkarda,
Alexander Rau,
Marco Reisert,
Thomas Küstner,
Alexandra Walter,
Elmar M. Merkle,
Martin Segeroth,
Joshy Cyriac,
Shan Yang,
Jakob Wasserthal
Abstract:
Purpose: To develop an open-source and easy-to-use segmentation model that can automatically and robustly segment most major anatomical structures in MR images independently of the MR sequence.
Materials and Methods: In this study we extended the capabilities of TotalSegmentator to MR images. 298 MR scans and 227 CT scans were used to segment 59 anatomical structures (20 organs, 18 bones, 11 mus…
▽ More
Purpose: To develop an open-source and easy-to-use segmentation model that can automatically and robustly segment most major anatomical structures in MR images independently of the MR sequence.
Materials and Methods: In this study we extended the capabilities of TotalSegmentator to MR images. 298 MR scans and 227 CT scans were used to segment 59 anatomical structures (20 organs, 18 bones, 11 muscles, 7 vessels, 3 tissue types) relevant for use cases such as organ volumetry, disease characterization, and surgical planning. The MR and CT images were randomly sampled from routine clinical studies and thus represent a real-world dataset (different ages, pathologies, scanners, body parts, sequences, contrasts, echo times, repetition times, field strengths, slice thicknesses and sites). We trained an nnU-Net segmentation algorithm on this dataset and calculated Dice similarity coefficients (Dice) to evaluate the model's performance.
Results: The model showed a Dice score of 0.824 (CI: 0.801, 0.842) on the test set, which included a wide range of clinical data with major pathologies. The model significantly outperformed two other publicly available segmentation models (Dice score, 0.824 versus 0.762; p<0.001 and 0.762 versus 0.542; p<0.001). On the CT image test set of the original TotalSegmentator paper it almost matches the performance of the original TotalSegmentator (Dice score, 0.960 versus 0.970; p<0.001).
Conclusion: Our proposed model extends the capabilities of TotalSegmentator to MR images. The annotated dataset (https://zenodo.org/doi/10.5281/zenodo.11367004) and open-source toolkit (https://www.github.com/wasserth/TotalSegmentator) are publicly available.
△ Less
Submitted 29 May, 2024;
originally announced May 2024.
-
Flux rope modeling of the 2022 Sep 5 CME observed by Parker Solar Probe and Solar Orbiter from 0.07 to 0.69 au
Authors:
Emma E. Davies,
Hannah T. Rüdisser,
Ute V. Amerstorfer,
Christian Möstl,
Maike Bauer,
Eva Weiler,
Tanja Amerstorfer,
Satabdwa Majumdar,
Phillip Hess,
Andreas J. Weiss,
Martin A. Reiss,
Lucie M. Green,
David M. Long,
Teresa Nieves-Chinchilla,
Domenico Trotta,
Timothy S. Horbury,
Helen O'Brien,
Edward Fauchon-Jones,
Jean Morris,
Christopher J. Owen,
Stuart D. Bale,
Justin C. Kasper
Abstract:
As both Parker Solar Probe (PSP) and Solar Orbiter (SolO) reach heliocentric distances closer to the Sun, they present an exciting opportunity to study the structure of CMEs in the inner heliosphere. We present an analysis of the global flux rope structure of the 2022 September 5 CME event that impacted PSP at a heliocentric distance of only 0.07 au and SolO at 0.69 au. We compare in situ measurem…
▽ More
As both Parker Solar Probe (PSP) and Solar Orbiter (SolO) reach heliocentric distances closer to the Sun, they present an exciting opportunity to study the structure of CMEs in the inner heliosphere. We present an analysis of the global flux rope structure of the 2022 September 5 CME event that impacted PSP at a heliocentric distance of only 0.07 au and SolO at 0.69 au. We compare in situ measurements at PSP and SolO to determine global and local expansion measures, finding a good agreement between magnetic field relationships with heliocentric distance, but significant differences with respect to flux rope size. We use PSP/WISPR images as input to the ELEvoHI model, providing a direct link between remote and in situ observations; we find a large discrepancy between the resulting modeled arrival times, suggesting that the underlying model assumptions may not be suitable when using data obtained close to the Sun, where the drag regime is markedly different in comparison to larger heliocentric distances. Finally, we fit the SolO/MAG and PSP/FIELDS data independently with the 3DCORE model and find that many parameters are consistent between spacecraft, however, challenges are apparent when reconstructing a global 3D structure that aligns with arrival times at PSP and Solar Orbiter, likely due to the large radial and longitudinal separations between spacecraft. From our model results, it is clear the solar wind background speed and drag regime strongly affect the modeled expansion and propagation of CMEs and need to be taken into consideration.
△ Less
Submitted 19 July, 2024; v1 submitted 17 May, 2024;
originally announced May 2024.
-
Adolescent sports participation and health in early adulthood: An observational study
Authors:
Ajinkya H. Kokandakar,
Yuzhou Lin,
Steven Jin,
Jordan Weiss,
Amanda R. Rabinowitz,
Reuben A. Buford May,
Dylan Small,
Sameer K. Deshpande
Abstract:
We study the impact of teenage sports participation on early-adulthood health using longitudinal data from the National Study of Youth and Religion. We focus on two primary outcomes measured at ages 23--28 -- self-rated health and total score on the PHQ9 Patient Depression Questionnaire -- and control for several potential confounders related to demographics and family socioeconomic status. To pro…
▽ More
We study the impact of teenage sports participation on early-adulthood health using longitudinal data from the National Study of Youth and Religion. We focus on two primary outcomes measured at ages 23--28 -- self-rated health and total score on the PHQ9 Patient Depression Questionnaire -- and control for several potential confounders related to demographics and family socioeconomic status. To probe the possibility that certain types of sports participation may have larger effects on health than others, we conduct a matched observational study at each level within a hierarchy of exposures. Our hierarchy ranges from broadly defined exposures (e.g., participation in any organized after-school activity) to narrow (e.g., participation in collision sports). We deployed an ordered testing approach that exploits the hierarchical relationships between our exposure definitions to perform our analyses while maintaining a fixed family-wise error rate. Compared to teenagers who did not participate in any after-school activities, those who participated in sports had statistically significantly better self-rated and mental health outcomes in early adulthood.
△ Less
Submitted 6 May, 2024;
originally announced May 2024.
-
Understanding the effects of spacecraft trajectories through solar coronal mass ejection flux ropes using 3DCOREweb
Authors:
Hannah Theresa Rüdisser,
Andreas Jeffrey Weiss,
Justin Le Louëdec,
Ute V. Amerstorfer,
Christian Möstl,
Emma E. Davies,
Helmut Lammer
Abstract:
This study investigates the impact of spacecraft positioning and trajectory on in situ signatures of coronal mass ejections (CMEs). Employing the 3DCORE model, a 3D flux rope model that can generate in situ profiles for any given point in space and time, we conduct forward modeling to analyze such signatures for various latitudinal and longitudinal positions, with respect to the flux rope apex, at…
▽ More
This study investigates the impact of spacecraft positioning and trajectory on in situ signatures of coronal mass ejections (CMEs). Employing the 3DCORE model, a 3D flux rope model that can generate in situ profiles for any given point in space and time, we conduct forward modeling to analyze such signatures for various latitudinal and longitudinal positions, with respect to the flux rope apex, at 0.8~au. Using this approach, we explore the appearance of the resulting in situ profiles for different flux rope types, with different handedness and inclination angles, for both high and low twist CMEs. Our findings reveal that CMEs exhibit distinct differences in signatures between apex hits and flank encounters, with the latter displaying elongated profiles with reduced rotation. However, constant, non-rotating in situ signatures are only observed for flank encounters of low twist CMEs, suggesting the existence of untwisted magnetic field lines within CME legs. Additionally, our study confirms the unambiguous appearance of different flux rope types in in situ signatures, contributing to the broader understanding and interpretation of observational data. Given the model assumptions, this may refute trajectory effects to be the cause for mismatching flux rope types as identified in solar signatures. While acknowledging limitations inherent in our model, such as the assumption of constant twist and non-deformable torus-like shape, we still draw relevant conclusions within the context of global magnetic field structures of CMEs and the potential for distinguishing flux rope types based on in situ observations.
△ Less
Submitted 17 July, 2024; v1 submitted 6 May, 2024;
originally announced May 2024.
-
Aftershocks as a time independant phenomenon
Authors:
A. Mathey,
J. Crassous,
D. Marsan,
J. Weiss,
A. Amon
Abstract:
Sequences of aftershocks following Omori's empirical law are observed after most major earthquakes, as well as in laboratory-scale fault-mimicking experiments. Nevertheless, the origin of this memory effect is still unclear. In this letter, we present an analytical framework for treating labquake and earthquake catalogs on an equal footing. Using this analysis method, we show that when memory is c…
▽ More
Sequences of aftershocks following Omori's empirical law are observed after most major earthquakes, as well as in laboratory-scale fault-mimicking experiments. Nevertheless, the origin of this memory effect is still unclear. In this letter, we present an analytical framework for treating labquake and earthquake catalogs on an equal footing. Using this analysis method, we show that when memory is considered to be in deformation and not in time, all data collapse onto a single master curve, showing that the timescale is entirely fixed by the inverse of the strain rate.
△ Less
Submitted 27 March, 2024;
originally announced March 2024.
-
On the Mesoscale Structure of CMEs at Mercury's Orbit: BepiColombo and Parker Solar Probe Observations
Authors:
Erika Palmerio,
Fernando Carcaboso,
Leng Ying Khoo,
Tarik M. Salman,
Beatriz Sánchez-Cano,
Benjamin J. Lynch,
Yeimy J. Rivera,
Sanchita Pal,
Teresa Nieves-Chinchilla,
Andreas J. Weiss,
David Lario,
Johannes Z. D. Mieth,
Daniel Heyner,
Michael L. Stevens,
Orlando M. Romeo,
Andrei N. Zhukov,
Luciano Rodriguez,
Christina O. Lee,
Christina M. S. Cohen,
Laura Rodríguez-García,
Phyllis L. Whittlesey,
Nina Dresing,
Philipp Oleynik,
Immanuel C. Jebaraj,
David Fischer
, et al. (5 additional authors not shown)
Abstract:
On 2022 February 15, an impressive filament eruption was observed off the solar eastern limb from three remote-sensing viewpoints, namely Earth, STEREO-A, and Solar Orbiter. In addition to representing the most-distant observed filament at extreme ultraviolet wavelengths -- captured by Solar Orbiter's field of view extending to above 6 $R_{\odot}$ -- this event was also associated with the release…
▽ More
On 2022 February 15, an impressive filament eruption was observed off the solar eastern limb from three remote-sensing viewpoints, namely Earth, STEREO-A, and Solar Orbiter. In addition to representing the most-distant observed filament at extreme ultraviolet wavelengths -- captured by Solar Orbiter's field of view extending to above 6 $R_{\odot}$ -- this event was also associated with the release of a fast ($\sim$2200 km$\cdot$s$^{-1}$) coronal mass ejection (CME) that was directed towards BepiColombo and Parker Solar Probe. These two probes were separated by 2$^{\circ}$ in latitude, 4$^{\circ}$ in longitude, and 0.03 au in radial distance around the time of the CME-driven shock arrival in situ. The relative proximity of the two probes to each other and to the Sun ($\sim$0.35 au) allows us to study the mesoscale structure of CMEs at Mercury's orbit for the first time. We analyse similarities and differences in the main CME-related structures measured at the two locations, namely the interplanetary shock, the sheath region, and the magnetic ejecta. We find that, despite the separation between the two spacecraft being well within the typical uncertainties associated with determination of CME geometric parameters from remote-sensing observations, the two sets of in-situ measurements display some profound differences that make understanding of the overall 3D CME structure particularly challenging. Finally, we discuss our findings within the context of space weather at Mercury's distances and in terms of the need to investigate solar transients via spacecraft constellations with small separations, which has been gaining significant attention during recent years.
△ Less
Submitted 3 January, 2024;
originally announced January 2024.
-
Temporal Supervised Contrastive Learning for Modeling Patient Risk Progression
Authors:
Shahriar Noroozizadeh,
Jeremy C. Weiss,
George H. Chen
Abstract:
We consider the problem of predicting how the likelihood of an outcome of interest for a patient changes over time as we observe more of the patient data. To solve this problem, we propose a supervised contrastive learning framework that learns an embedding representation for each time step of a patient time series. Our framework learns the embedding space to have the following properties: (1) nea…
▽ More
We consider the problem of predicting how the likelihood of an outcome of interest for a patient changes over time as we observe more of the patient data. To solve this problem, we propose a supervised contrastive learning framework that learns an embedding representation for each time step of a patient time series. Our framework learns the embedding space to have the following properties: (1) nearby points in the embedding space have similar predicted class probabilities, (2) adjacent time steps of the same time series map to nearby points in the embedding space, and (3) time steps with very different raw feature vectors map to far apart regions of the embedding space. To achieve property (3), we employ a nearest neighbor pairing mechanism in the raw feature space. This mechanism also serves as an alternative to data augmentation, a key ingredient of contrastive learning, which lacks a standard procedure that is adequately realistic for clinical tabular data, to our knowledge. We demonstrate that our approach outperforms state-of-the-art baselines in predicting mortality of septic patients (MIMIC-III dataset) and tracking progression of cognitive impairment (ADNI dataset). Our method also consistently recovers the correct synthetic dataset embedding structure across experiments, a feat not achieved by baselines. Our ablation experiments show the pivotal role of our nearest neighbor pairing.
△ Less
Submitted 10 December, 2023;
originally announced December 2023.
-
On the Low-Frequency Dynamics of Turbulent Separation Bubbles
Authors:
C. Cura,
A. Hanifi,
A. V. G. Cavalieri,
J. Weiss
Abstract:
The low-frequency modal and non-modal stability characteristics of an incompressible, pressure-gradient-induced turbulent separation bubble (TSB) are investigated with the objective of studying the mechanism responsible for the low-frequency contraction and expansion (breathing) commonly observed in experimental studies. The configuration of interest is a TSB generated on a flat test surface by a…
▽ More
The low-frequency modal and non-modal stability characteristics of an incompressible, pressure-gradient-induced turbulent separation bubble (TSB) are investigated with the objective of studying the mechanism responsible for the low-frequency contraction and expansion (breathing) commonly observed in experimental studies. The configuration of interest is a TSB generated on a flat test surface by a succession of adverse and favourable pressure gradients. The base flow selected for the analysis is the average TSB from the direct numerical simulation of Coleman et al. (J. Fluid Mech., vol. 847, 2018). Global linear stability analysis reveals that the flow is globally stable for wavenumbers. The mode closest to the stability threshold appears to occur at zero frequency and low, non-zero spanwise wavenumber. Resolvent analysis is then employed to examine the forced dynamics of the flow. At low frequency, a region of low, non-zero spanwise wavenumber is also discernible, where the receptivity appears to be driven by the identified weakly damped global mode. The results from resolvent analysis are compared to the unsteady experimental database of Le Floc'h et al. (J. Fluid Mech., vol. 902, 2020) in a similar TSB flow. The alignment between the optimal response and the first spectral proper orthogonal decomposition mode computed from the experiments is shown to exceed 95 %, while the spanwise wavenumber of the optimal response is consistent with that of the low-frequency breathing motion captured experimentally. This indicates that the fluctuations observed experimentally at low frequency closely match the response computed from resolvent analysis. Based on these results, we propose that the forced dynamics of the flow, driven by the weakly damped global mode, serve as a plausible mechanism for the origin of the low-frequency breathing motion commonly observed in experimental studies of TSBs.
△ Less
Submitted 22 November, 2023;
originally announced November 2023.
-
The Future of Astronomical Data Infrastructure: Meeting Report
Authors:
Michael R. Blanton,
Janet D. Evans,
Dara Norman,
William O'Mullane,
Adrian Price-Whelan,
Luca Rizzi,
Alberto Accomazzi,
Megan Ansdell,
Stephen Bailey,
Paul Barrett,
Steven Berukoff,
Adam Bolton,
Julian Borrill,
Kelle Cruz,
Julianne Dalcanton,
Vandana Desai,
Gregory P. Dubois-Felsmann,
Frossie Economou,
Henry Ferguson,
Bryan Field,
Dan Foreman-Mackey,
Jaime Forero-Romero,
Niall Gaffney,
Kim Gillies,
Matthew J. Graham
, et al. (47 additional authors not shown)
Abstract:
The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and productio…
▽ More
The astronomical community is grappling with the increasing volume and complexity of data produced by modern telescopes, due to difficulties in reducing, accessing, analyzing, and combining archives of data. To address this challenge, we propose the establishment of a coordinating body, an "entity," with the specific mission of enhancing the interoperability, archiving, distribution, and production of both astronomical data and software. This report is the culmination of a workshop held in February 2023 on the Future of Astronomical Data Infrastructure. Attended by 70 scientists and software professionals from ground-based and space-based missions and archives spanning the entire spectrum of astronomical research, the group deliberated on the prevailing state of software and data infrastructure in astronomy, identified pressing issues, and explored potential solutions. In this report, we describe the ecosystem of astronomical data, its existing flaws, and the many gaps, duplication, inconsistencies, barriers to access, drags on productivity, missed opportunities, and risks to the long-term integrity of essential data sets. We also highlight the successes and failures in a set of deep dives into several different illustrative components of the ecosystem, included as an appendix.
△ Less
Submitted 7 November, 2023;
originally announced November 2023.
-
A Leptonic ALP Portal to the Dark Sector
Authors:
Giovanni Armando,
Paolo Panci,
Joachim Weiss,
Robert Ziegler
Abstract:
We discuss the leptonic ALP portal as a simple scenario that connects observed discrepancies in anomalous magnetic moments to the Dark Matter relic abundance. In this framework an axion-like particle in the multi-MeV range couples to SM leptons and a DM fermion, with mass above the ALP mass but below a GeV. The ALP contributes to $(g-2)_μ$ and $(g-2)_e$ dominantly through 2-loop Barr-Zee diagrams,…
▽ More
We discuss the leptonic ALP portal as a simple scenario that connects observed discrepancies in anomalous magnetic moments to the Dark Matter relic abundance. In this framework an axion-like particle in the multi-MeV range couples to SM leptons and a DM fermion, with mass above the ALP mass but below a GeV. The ALP contributes to $(g-2)_μ$ and $(g-2)_e$ dominantly through 2-loop Barr-Zee diagrams, while the DM abundance is generated by $p$-wave annihilation to ALP pairs. Constraints from beam-dump experiments, colliders and CMB probes are very stringent, and restrict the viable parameter space to a rather narrow region that will be tested in the near future.
△ Less
Submitted 27 March, 2024; v1 submitted 9 October, 2023;
originally announced October 2023.
-
LHC tau-pair production constraints on $a_τ$ and $d_τ$
Authors:
Ulrich Haisch,
Luc Schnell,
Joachim Weiss
Abstract:
We point out that relevant constraints on the anomalous magnetic ($a_τ$) and electric ($d_τ$) moment of the tau lepton can be derived from tau-pair production measurements performed at the LHC. Our conclusion is based on the observation that the leading relative deviations from the Standard Model prediction for $pp \to τ^+ τ^-$ due to $a_τ$ and $d_τ$ are enhanced at high energies. Less precise mea…
▽ More
We point out that relevant constraints on the anomalous magnetic ($a_τ$) and electric ($d_τ$) moment of the tau lepton can be derived from tau-pair production measurements performed at the LHC. Our conclusion is based on the observation that the leading relative deviations from the Standard Model prediction for $pp \to τ^+ τ^-$ due to $a_τ$ and $d_τ$ are enhanced at high energies. Less precise measurements at hadron colliders can therefore offer the same or better sensitivity to new physics with respect to high-precision low-energy measurements performed at lepton machines. We derive bounds on $a_τ$ and $d_τ$ using the full LHC Run II data set on tau-pair production and compare our findings with the current best limits on the tau anomalous moments.
△ Less
Submitted 15 February, 2024; v1 submitted 26 July, 2023;
originally announced July 2023.
-
Observation of high-energy neutrinos from the Galactic plane
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Athanasiadou,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus
, et al. (364 additional authors not shown)
Abstract:
The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth's atmosphere, has been a mystery for over a century. Due to deflection in interstellar magnetic fields, cosmic rays from the Milky Way arrive at Earth from random directions. However, near their sources and during propagation, cosmic rays interact with matter and produce high-energy neutrinos. We search for neutrin…
▽ More
The origin of high-energy cosmic rays, atomic nuclei that continuously impact Earth's atmosphere, has been a mystery for over a century. Due to deflection in interstellar magnetic fields, cosmic rays from the Milky Way arrive at Earth from random directions. However, near their sources and during propagation, cosmic rays interact with matter and produce high-energy neutrinos. We search for neutrino emission using machine learning techniques applied to ten years of data from the IceCube Neutrino Observatory. We identify neutrino emission from the Galactic plane at the 4.5$σ$ level of significance, by comparing diffuse emission models to a background-only hypothesis. The signal is consistent with modeled diffuse emission from the Galactic plane, but could also arise from a population of unresolved point sources.
△ Less
Submitted 10 July, 2023;
originally announced July 2023.
-
EZClone: Improving DNN Model Extraction Attack via Shape Distillation from GPU Execution Profiles
Authors:
Jonah O'Brien Weiss,
Tiago Alves,
Sandip Kundu
Abstract:
Deep Neural Networks (DNNs) have become ubiquitous due to their performance on prediction and classification problems. However, they face a variety of threats as their usage spreads. Model extraction attacks, which steal DNNs, endanger intellectual property, data privacy, and security. Previous research has shown that system-level side-channels can be used to leak the architecture of a victim DNN,…
▽ More
Deep Neural Networks (DNNs) have become ubiquitous due to their performance on prediction and classification problems. However, they face a variety of threats as their usage spreads. Model extraction attacks, which steal DNNs, endanger intellectual property, data privacy, and security. Previous research has shown that system-level side-channels can be used to leak the architecture of a victim DNN, exacerbating these risks. We propose two DNN architecture extraction techniques catering to various threat models. The first technique uses a malicious, dynamically linked version of PyTorch to expose a victim DNN architecture through the PyTorch profiler. The second, called EZClone, exploits aggregate (rather than time-series) GPU profiles as a side-channel to predict DNN architecture, employing a simple approach and assuming little adversary capability as compared to previous work. We investigate the effectiveness of EZClone when minimizing the complexity of the attack, when applied to pruned models, and when applied across GPUs. We find that EZClone correctly predicts DNN architectures for the entire set of PyTorch vision architectures with 100% accuracy. No other work has shown this degree of architecture prediction accuracy with the same adversarial constraints or using aggregate side-channel information. Prior work has shown that, once a DNN has been successfully cloned, further attacks such as model evasion or model inversion can be accelerated significantly.
△ Less
Submitted 6 April, 2023;
originally announced April 2023.
-
Multi-label classification of open-ended questions with BERT
Authors:
Matthias Schonlau,
Julia Weiß,
Jan Marquardt
Abstract:
Open-ended questions in surveys are valuable because they do not constrain the respondent's answer, thereby avoiding biases. However, answers to open-ended questions are text data which are harder to analyze. Traditionally, answers were manually classified as specified in the coding manual. Most of the effort to automate coding has gone into the easier problem of single label prediction, where ans…
▽ More
Open-ended questions in surveys are valuable because they do not constrain the respondent's answer, thereby avoiding biases. However, answers to open-ended questions are text data which are harder to analyze. Traditionally, answers were manually classified as specified in the coding manual. Most of the effort to automate coding has gone into the easier problem of single label prediction, where answers are classified into a single code. However, open-ends that require multi-label classification, i.e., that are assigned multiple codes, occur frequently. This paper focuses on multi-label classification of text answers to open-ended survey questions in social science surveys. We evaluate the performance of the transformer-based architecture BERT for the German language in comparison to traditional multi-label algorithms (Binary Relevance, Label Powerset, ECC) in a German social science survey, the GLES Panel (N=17,584, 55 labels). We find that classification with BERT (forcing at least one label) has the smallest 0/1 loss (13.1%) among methods considered (18.9%-21.6%). As expected, it is much easier to correctly predict answer texts that correspond to a single label (7.1% loss) than those that correspond to multiple labels ($\sim$50% loss). Because BERT predicts zero labels for only 1.5% of the answers, forcing at least one label, while recommended, ultimately does not lower the 0/1 loss by much. Our work has important implications for social scientists: 1) We have shown multi-label classification with BERT works in the German language for open-ends. 2) For mildly multi-label classification tasks, the loss now appears small enough to allow for fully automatic classification (as compared to semi-automatic approaches). 3) Multi-label classification with BERT requires only a single model. The leading competitor, ECC, iterates through individual single label predictions.
△ Less
Submitted 6 April, 2023;
originally announced April 2023.
-
Critical Current Measurement of REBCO Cables by Using a Superconducting Transformer
Authors:
H. Yu,
J. Lu,
J. D. Weiss,
D. C. van der Laan
Abstract:
Development of REBCO cables that carry high electrical current in high magnetic field is crucial for future large-scale magnet applications. This experimental work presents the critical current measurements of two different REBCO cables by a test facility at the National High Magnetic Field Laboratory (NHMFL). The simple-stacked cable is made by the NHMFL by stacking 21 REBCO tapes without solderi…
▽ More
Development of REBCO cables that carry high electrical current in high magnetic field is crucial for future large-scale magnet applications. This experimental work presents the critical current measurements of two different REBCO cables by a test facility at the National High Magnetic Field Laboratory (NHMFL). The simple-stacked cable is made by the NHMFL by stacking 21 REBCO tapes without soldering. The Conductor-on-Round-Core (CORC) cable provided by Advanced Conductor Technologies has 21 layers of REBCO tapes with 2 tapes/layer. The test facility consists of a 12 T split solenoid magnet with 15 cm bore providing transverse field to the samples, a superconducting transformer (SCT) as a current source providing up to 45 kA current. Special attentions were paid to fabrication of solder joints between REBCO cables and the SCT output. The voltage-current traces were measured as a function of magnetic field at 4.2 K, from which the critical currents are determined. The details of this measurement are discussed.
△ Less
Submitted 13 February, 2023;
originally announced February 2023.
-
Individual Fairness under Uncertainty
Authors:
Wenbin Zhang,
Zichong Wang,
Juyong Kim,
Cheng Cheng,
Thomas Oommen,
Pradeep Ravikumar,
Jeremy Weiss
Abstract:
Algorithmic fairness, the research field of making machine learning (ML) algorithms fair, is an established area in ML. As ML technologies expand their application domains, including ones with high societal impact, it becomes essential to take fairness into consideration during the building of ML systems. Yet, despite its wide range of socially sensitive applications, most work treats the issue of…
▽ More
Algorithmic fairness, the research field of making machine learning (ML) algorithms fair, is an established area in ML. As ML technologies expand their application domains, including ones with high societal impact, it becomes essential to take fairness into consideration during the building of ML systems. Yet, despite its wide range of socially sensitive applications, most work treats the issue of algorithmic bias as an intrinsic property of supervised learning, i.e., the class label is given as a precondition. Unlike prior studies in fairness, we propose an individual fairness measure and a corresponding algorithm that deal with the challenges of uncertainty arising from censorship in class labels, while enforcing similar individuals to be treated similarly from a ranking perspective, free of the Lipschitz condition in the conventional individual fairness definition. We argue that this perspective represents a more realistic model of fairness research for real-world application deployment and show how learning with such a relaxed precondition draws new insights that better explains algorithmic fairness. We conducted experiments on four real-world datasets to evaluate our proposed method compared to other fairness models, demonstrating its superiority in minimizing discrimination while maintaining predictive performance with uncertainty present.
△ Less
Submitted 11 December, 2023; v1 submitted 15 February, 2023;
originally announced February 2023.
-
The effects of leaflet material properties on the simulated function of regurgitant mitral valves
Authors:
Wensi Wu,
Stephen Ching,
Patricia Sabin,
Devin W. Laurence,
Steve A. Maas,
Andras Lasso,
Jeffrey A. Weiss,
Matthew A. Jolley
Abstract:
Advances in three-dimensional imaging provide the ability to construct and analyze finite element (FE) models to evaluate the biomechanical behavior and function of atrioventricular valves. However, while obtaining patient-specific valve geometry is now possible, non-invasive measurement of patient-specific leaflet material properties remains nearly impossible. Both valve geometry and tissue prope…
▽ More
Advances in three-dimensional imaging provide the ability to construct and analyze finite element (FE) models to evaluate the biomechanical behavior and function of atrioventricular valves. However, while obtaining patient-specific valve geometry is now possible, non-invasive measurement of patient-specific leaflet material properties remains nearly impossible. Both valve geometry and tissue properties play a significant role in governing valve dynamics, leading to the central question of whether clinically relevant insights can be attained from FE analysis of atrioventricular valves without precise knowledge of tissue properties. As such we investigated 1) the influence of tissue extensibility and 2) the effects of constitutive model parameters and leaflet thickness on simulated valve function and mechanics. We compared metrics of valve function (e.g., leaflet coaptation and regurgitant orifice area) and mechanics (e.g., stress and strain) across one normal and three regurgitant mitral valve (MV) models with common mechanisms of regurgitation (annular dilation, leaflet prolapse, leaflet tethering) of both moderate and severe degree. We developed a novel fully-automated approach to accurately quantify regurgitant orifice areas of complex valve geometries. We found that the relative ordering of the mechanical and functional metrics was maintained across a group of valves using material properties up to 15% softer than the representative adult mitral constitutive model. Our findings suggest that FE simulations can be used to qualitatively compare how differences and alterations in valve structure affect relative atrioventricular valve function even in populations where material properties are not precisely known.
△ Less
Submitted 25 April, 2023; v1 submitted 9 February, 2023;
originally announced February 2023.
-
Fair Decision-making Under Uncertainty
Authors:
Wenbin Zhang,
Jeremy C. Weiss
Abstract:
There has been concern within the artificial intelligence (AI) community and the broader society regarding the potential lack of fairness of AI-based decision-making systems. Surprisingly, there is little work quantifying and guaranteeing fairness in the presence of uncertainty which is prevalent in many socially sensitive applications, ranging from marketing analytics to actuarial analysis and re…
▽ More
There has been concern within the artificial intelligence (AI) community and the broader society regarding the potential lack of fairness of AI-based decision-making systems. Surprisingly, there is little work quantifying and guaranteeing fairness in the presence of uncertainty which is prevalent in many socially sensitive applications, ranging from marketing analytics to actuarial analysis and recidivism prediction instruments. To this end, we study a longitudinal censored learning problem subject to fairness constraints, where we require that algorithmic decisions made do not affect certain individuals or social groups negatively in the presence of uncertainty on class label due to censorship. We argue that this formulation has a broader applicability to practical scenarios concerning fairness. We show how the newly devised fairness notions involving censored information and the general framework for fair predictions in the presence of censorship allow us to measure and mitigate discrimination under uncertainty that bridges the gap with real-world applications. Empirical evaluations on real-world discriminated datasets with censorship demonstrate the practicality of our approach.
△ Less
Submitted 29 January, 2023;
originally announced January 2023.
-
Evidence for neutrino emission from the nearby active galaxy NGC 1068
Authors:
IceCube Collaboration,
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
C. Alispach,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
A. Barbano,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay
, et al. (361 additional authors not shown)
Abstract:
We report three searches for high energy neutrino emission from astrophysical objects using data recorded with IceCube between 2011 and 2020. Improvements over previous work include new neutrino reconstruction and data calibration methods. In one search, the positions of 110 a priori selected gamma-ray sources were analyzed individually for a possible surplus of neutrinos over atmospheric and cosm…
▽ More
We report three searches for high energy neutrino emission from astrophysical objects using data recorded with IceCube between 2011 and 2020. Improvements over previous work include new neutrino reconstruction and data calibration methods. In one search, the positions of 110 a priori selected gamma-ray sources were analyzed individually for a possible surplus of neutrinos over atmospheric and cosmic background expectations. We found an excess of $79_{-20}^{+22}$ neutrinos associated with the nearby active galaxy NGC 1068 at a significance of 4.2$\,σ$. The excess, which is spatially consistent with the direction of the strongest clustering of neutrinos in the Northern Sky, is interpreted as direct evidence of TeV neutrino emission from a nearby active galaxy. The inferred flux exceeds the potential TeV gamma-ray flux by at least one order of magnitude.
△ Less
Submitted 8 February, 2024; v1 submitted 17 November, 2022;
originally announced November 2022.
-
Pre-analysis protocol for an observational study on the effects of adolescent sports participation on health in early adulthood
Authors:
Ajinkya H Kokandakar,
Yuzhou Lin,
Steven Jin,
Jordan Weiss,
Amanda R Rabinowitz,
Reuben A Buford May,
Dylan Small,
Sameer K Deshpande
Abstract:
We will study the impact of adolescent sports participation on early-adulthood health using longitudinal data from the National Study of Youth and Religion. We focus on two primary outcomes measured at ages 23--28 -- self-rated health and total score on the PHQ9 Patient Depression Questionnaire -- and control for several potential confounders related to demographics and family socioeconomic status…
▽ More
We will study the impact of adolescent sports participation on early-adulthood health using longitudinal data from the National Study of Youth and Religion. We focus on two primary outcomes measured at ages 23--28 -- self-rated health and total score on the PHQ9 Patient Depression Questionnaire -- and control for several potential confounders related to demographics and family socioeconomic status. Comparing outcomes between sports participants and matched non-sports participants with similar confounders is straightforward. Unfortunately, an analysis based on such a broad exposure cannot probe the possibility that participation in certain types of sports (e.g., collision sports like football or soccer) may have larger effects on health than others.
In this study, we introduce a hierarchy of exposure definitions, ranging from broad (participation in any after-school organized activity) to narrow (e.g., participation in limited-contact sports). We will perform separate matched observational studies, one for each definition, to estimate the health effects of several levels of sports participation. In order to conduct these studies while maintaining a fixed family-wise error rate, we deployed an ordered testing approach that exploits the logical relationships between exposure definitions. Our study will also consider several secondary outcomes including body mass index, life satisfaction, and problematic drinking behavior.
△ Less
Submitted 30 November, 2023; v1 submitted 3 November, 2022;
originally announced November 2022.
-
G-Augment: Searching for the Meta-Structure of Data Augmentation Policies for ASR
Authors:
Gary Wang,
Ekin D. Cubuk,
Andrew Rosenberg,
Shuyang Cheng,
Ron J. Weiss,
Bhuvana Ramabhadran,
Pedro J. Moreno,
Quoc V. Le,
Daniel S. Park
Abstract:
Data augmentation is a ubiquitous technique used to provide robustness to automatic speech recognition (ASR) training. However, even as so much of the ASR training process has become automated and more "end-to-end", the data augmentation policy (what augmentation functions to use, and how to apply them) remains hand-crafted. We present Graph-Augment, a technique to define the augmentation space as…
▽ More
Data augmentation is a ubiquitous technique used to provide robustness to automatic speech recognition (ASR) training. However, even as so much of the ASR training process has become automated and more "end-to-end", the data augmentation policy (what augmentation functions to use, and how to apply them) remains hand-crafted. We present Graph-Augment, a technique to define the augmentation space as directed acyclic graphs (DAGs) and search over this space to optimize the augmentation policy itself. We show that given the same computational budget, policies produced by G-Augment are able to perform better than SpecAugment policies obtained by random search on fine-tuning tasks on CHiME-6 and AMI. G-Augment is also able to establish a new state-of-the-art ASR performance on the CHiME-6 evaluation set (30.7% WER). We further demonstrate that G-Augment policies show better transfer properties across warm-start to cold-start training and model size compared to random-searched SpecAugment policies.
△ Less
Submitted 24 October, 2022; v1 submitted 19 October, 2022;
originally announced October 2022.
-
Learning Clinical Concepts for Predicting Risk of Progression to Severe COVID-19
Authors:
Helen Zhou,
Cheng Cheng,
Kelly J. Shields,
Gursimran Kochhar,
Tariq Cheema,
Zachary C. Lipton,
Jeremy C. Weiss
Abstract:
With COVID-19 now pervasive, identification of high-risk individuals is crucial. Using data from a major healthcare provider in Southwestern Pennsylvania, we develop survival models predicting severe COVID-19 progression. In this endeavor, we face a tradeoff between more accurate models relying on many features and less accurate models relying on a few features aligned with clinician intuition. Co…
▽ More
With COVID-19 now pervasive, identification of high-risk individuals is crucial. Using data from a major healthcare provider in Southwestern Pennsylvania, we develop survival models predicting severe COVID-19 progression. In this endeavor, we face a tradeoff between more accurate models relying on many features and less accurate models relying on a few features aligned with clinician intuition. Complicating matters, many EHR features tend to be under-coded, degrading the accuracy of smaller models. In this study, we develop two sets of high-performance risk scores: (i) an unconstrained model built from all available features; and (ii) a pipeline that learns a small set of clinical concepts before training a risk predictor. Learned concepts boost performance over the corresponding features (C-index 0.858 vs. 0.844) and demonstrate improvements over (i) when evaluated out-of-sample (subsequent time periods). Our models outperform previous works (C-index 0.844-0.872 vs. 0.598-0.810).
△ Less
Submitted 27 August, 2022;
originally announced August 2022.
-
Spin-wave theory in a randomly disordered lattice: A Heisenberg ferromagnet
Authors:
Z. J. Weiss,
A. R. Massih
Abstract:
Starting from the hamiltonian for the Heisenberg ferromagnet which comprise randomly distributed nonmagnetic ions as impurities in a Bravais lattice, we express the spin operators by means of the Dyson-Maleev transformation in terms of the Bose operators of the second quantization. Then by using methods of quantum statistical field theory, we derive the partition function and the free energy for t…
▽ More
Starting from the hamiltonian for the Heisenberg ferromagnet which comprise randomly distributed nonmagnetic ions as impurities in a Bravais lattice, we express the spin operators by means of the Dyson-Maleev transformation in terms of the Bose operators of the second quantization. Then by using methods of quantum statistical field theory, we derive the partition function and the free energy for the system. We adopt the Matsubara thermal perturbation method to a portion of the hamiltonian which describes the interaction between magnons and the stationary field of nonmagnetic ions. Upon averaging over all possible distributions of impurities, we express the free energy of the system as a function of the mean impurity concentration. Subsequently, we set up the double-time single particle Green function at temperature T in the momentum space in terms of magnon operators and derive the equation of motion for the Green function through the Heisenberg equation of motion and then solve the resulting equation. From this, we calculate the self-energy and then the spectral density for the system. We apply the formalism to the case of the simple cubic lattice and compute the density of states, the spectral density function and the lifetime of the magnons as a function of energy for several values of the mean concentration of nonmagnetic ions in the lattice. We calculate magnon energy spectrum as a function of average impurity concentration fraction c, which shows that for low lying states, the excitation energy increases continuously with c in the studied range 0.1 < c < 0.7. We use the spectral density function to compute some thermal quantities. We have obtained closed form expressions for the configurationally averaged physical quantities of interest in a unified fashion as functions of c to any order of c applicable below a critical percolation concentration.
△ Less
Submitted 14 August, 2022;
originally announced August 2022.
-
Rapid Production of Accurate Embedded-Atom Method Potentials for Metal Alloys
Authors:
Elan J. Weiss,
Logan Ward,
Christian Oberdorfer,
Travis Withrow,
David C. Riegner,
Anupriya Agrawal,
Wolfgang Windl
Abstract:
A critical limitation to the wide-scale use of classical molecular dynamics for alloy design is the limited availability of suitable interatomic potentials. Here, we introduce the Rapid Alloy Method for Producing Accurate General Empirical Potentials or RAMPAGE, a computationally economical procedure to generate binary embedded-atom model potentials from already-existing single-element potentials…
▽ More
A critical limitation to the wide-scale use of classical molecular dynamics for alloy design is the limited availability of suitable interatomic potentials. Here, we introduce the Rapid Alloy Method for Producing Accurate General Empirical Potentials or RAMPAGE, a computationally economical procedure to generate binary embedded-atom model potentials from already-existing single-element potentials that can be further combined into multi-component alloy potentials. We present the quality of RAMPAGE calibrated Finnis-Sinclair type EAM potentials using binary Ag-Al and ternary Ag-Au-Cu as case studies. We demonstrate that RAMPAGE potentials can reproduce bulk properties and forces with greater accuracy than that of other alloy potentials. In some simulations, it is observed the quality of the optimized cross interactions can exceed that of the original off-the-shelf elemental potential inputs.
△ Less
Submitted 3 August, 2022;
originally announced August 2022.
-
History effects in the creep of a disordered brittle material
Authors:
Tero Mäkinen,
Jérôme Weiss,
David Amitrano,
Philippe Roux
Abstract:
We study the creep behavior of a disordered brittle material (concrete) under successive loading steps, using acoustic emission and ultrasonic sensing to track internal damage. The primary creep rate is observed to follow a (Omori-type) power-law decay in the strain rate, the number of acoustic emission events as well as the amplitudes of the ultrasonic beams, supporting a brittle-creep mechanism.…
▽ More
We study the creep behavior of a disordered brittle material (concrete) under successive loading steps, using acoustic emission and ultrasonic sensing to track internal damage. The primary creep rate is observed to follow a (Omori-type) power-law decay in the strain rate, the number of acoustic emission events as well as the amplitudes of the ultrasonic beams, supporting a brittle-creep mechanism. The main outcome is however the discovery of unexpected history effects that make the material less prone to creep when it has been previously deformed and damaged under primary creep at a lower applied stress. With the help of a progressive damage model implementing thermal activation, we interpret this as an aging-under-stress phenomenon: during an initial creep step at relatively low applied stress, the easy-to-damage sites are exhausted first, depleting the excitation spectrum at low stress gap values. Consequently, upon reloading under a larger applied stress, although previously damaged, the material creeps (and damages) less than it would under the same stress but without pre-creeping. Besides shedding a new light on the fundamental physics of creep of disordered brittle materials, this has important practical consequences in the interpretation of some experimental procedures, such as stress-stepping experiments.
△ Less
Submitted 22 June, 2022;
originally announced June 2022.
-
Hardening DNNs against Transfer Attacks during Network Compression using Greedy Adversarial Pruning
Authors:
Jonah O'Brien Weiss,
Tiago Alves,
Sandip Kundu
Abstract:
The prevalence and success of Deep Neural Network (DNN) applications in recent years have motivated research on DNN compression, such as pruning and quantization. These techniques accelerate model inference, reduce power consumption, and reduce the size and complexity of the hardware necessary to run DNNs, all with little to no loss in accuracy. However, since DNNs are vulnerable to adversarial in…
▽ More
The prevalence and success of Deep Neural Network (DNN) applications in recent years have motivated research on DNN compression, such as pruning and quantization. These techniques accelerate model inference, reduce power consumption, and reduce the size and complexity of the hardware necessary to run DNNs, all with little to no loss in accuracy. However, since DNNs are vulnerable to adversarial inputs, it is important to consider the relationship between compression and adversarial robustness. In this work, we investigate the adversarial robustness of models produced by several irregular pruning schemes and by 8-bit quantization. Additionally, while conventional pruning removes the least important parameters in a DNN, we investigate the effect of an unconventional pruning method: removing the most important model parameters based on the gradient on adversarial inputs. We call this method Greedy Adversarial Pruning (GAP) and we find that this pruning method results in models that are resistant to transfer attacks from their uncompressed counterparts.
△ Less
Submitted 15 June, 2022;
originally announced June 2022.
-
Logarithmic vs Andrade's transient creep: the role of elastic stress redistribution
Authors:
Jerome Weiss,
David Amitrano
Abstract:
Creep is defined as time-dependent deformation and rupture processes taking place within a material subjected to a constant applied stress smaller than its athermal, time-independent strength. This time-dependence is classically attributed to thermal activation of local deformation events. The phenomenology of creep is characterized by several ubiquitous but empirical rheological and scaling laws.…
▽ More
Creep is defined as time-dependent deformation and rupture processes taking place within a material subjected to a constant applied stress smaller than its athermal, time-independent strength. This time-dependence is classically attributed to thermal activation of local deformation events. The phenomenology of creep is characterized by several ubiquitous but empirical rheological and scaling laws. We focus here on primary creep following the onset of loading, for which a power law decay of the strain-rate is observed, with the exponent p varying between '0.4 and 1, this upper bound defining the so-called logarithmic creep. Although this phenomenology is known for more than a century, the physical origin of Andrade-like (p <1) creep remains unclear and debated. Here we show that p <1 values arise from the interplay between thermal activation and elastic stress redistribution. The latter stimulates creep dynamics from a shortening of waiting times between successive events, is associated to material damage and possibly, at high temperature and/or stresses, gives rise to avalanches of deformation events.
△ Less
Submitted 23 February, 2023; v1 submitted 14 June, 2022;
originally announced June 2022.
-
FedNorm: Modality-Based Normalization in Federated Learning for Multi-Modal Liver Segmentation
Authors:
Tobias Bernecker,
Annette Peters,
Christopher L. Schlett,
Fabian Bamberg,
Fabian Theis,
Daniel Rueckert,
Jakob Weiß,
Shadi Albarqouni
Abstract:
Given the high incidence and effective treatment options for liver diseases, they are of great socioeconomic importance. One of the most common methods for analyzing CT and MRI images for diagnosis and follow-up treatment is liver segmentation. Recent advances in deep learning have demonstrated encouraging results for automatic liver segmentation. Despite this, their success depends primarily on t…
▽ More
Given the high incidence and effective treatment options for liver diseases, they are of great socioeconomic importance. One of the most common methods for analyzing CT and MRI images for diagnosis and follow-up treatment is liver segmentation. Recent advances in deep learning have demonstrated encouraging results for automatic liver segmentation. Despite this, their success depends primarily on the availability of an annotated database, which is often not available because of privacy concerns. Federated Learning has been recently proposed as a solution to alleviate these challenges by training a shared global model on distributed clients without access to their local databases. Nevertheless, Federated Learning does not perform well when it is trained on a high degree of heterogeneity of image data due to multi-modal imaging, such as CT and MRI, and multiple scanner types. To this end, we propose Fednorm and its extension \fednormp, two Federated Learning algorithms that use a modality-based normalization technique. Specifically, Fednorm normalizes the features on a client-level, while Fednorm+ employs the modality information of single slices in the feature normalization. Our methods were validated using 428 patients from six publicly available databases and compared to state-of-the-art Federated Learning algorithms and baseline models in heterogeneous settings (multi-institutional, multi-modal data). The experimental results demonstrate that our methods show an overall acceptable performance, achieve Dice per patient scores up to 0.961, consistently outperform locally trained models, and are on par or slightly better than centralized models.
△ Less
Submitted 23 May, 2022;
originally announced May 2022.
-
Framework and Tools for the Simulation and Analysis of the Radio Emission from Air Showers at IceCube
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus
, et al. (361 additional authors not shown)
Abstract:
The Surface Enhancement of the IceTop air-shower array will include the addition of radio antennas and scintillator panels, co-located with the existing ice-Cherenkov tanks and covering an area of about 1 km$^2$. Together, these will increase the sensitivity of the IceCube Neutrino Observatory to the electromagnetic and muonic components of cosmic-ray-induced air showers at the South Pole. The inc…
▽ More
The Surface Enhancement of the IceTop air-shower array will include the addition of radio antennas and scintillator panels, co-located with the existing ice-Cherenkov tanks and covering an area of about 1 km$^2$. Together, these will increase the sensitivity of the IceCube Neutrino Observatory to the electromagnetic and muonic components of cosmic-ray-induced air showers at the South Pole. The inclusion of the radio technique necessitates an expanded set of simulation and analysis tools to explore the radio-frequency emission from air showers in the 70 MHz to 350 MHz band. In this paper we describe the software modules that have been developed to work with time- and frequency-domain information within IceCube's existing software framework, IceTray, which is used by the entire IceCube collaboration. The software includes a method by which air-shower simulation, generated using CoREAS, can be reused via waveform interpolation, thus overcoming a significant computational hurdle in the field.
△ Less
Submitted 4 May, 2022;
originally announced May 2022.
-
First Search for Unstable Sterile Neutrinos with the IceCube Neutrino Observatory
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus
, et al. (359 additional authors not shown)
Abstract:
We present a search for an unstable sterile neutrino by looking for a matter-induced signal in eight years of atmospheric $ν_μ$ data collected from 2011 to 2019 at the IceCube Neutrino Observatory. Both the (stable) three-neutrino and the 3+1 sterile neutrino models are disfavored relative to the unstable sterile neutrino model, though with $p$-values of 2.5\% and 0.81\%, respectively, we do not o…
▽ More
We present a search for an unstable sterile neutrino by looking for a matter-induced signal in eight years of atmospheric $ν_μ$ data collected from 2011 to 2019 at the IceCube Neutrino Observatory. Both the (stable) three-neutrino and the 3+1 sterile neutrino models are disfavored relative to the unstable sterile neutrino model, though with $p$-values of 2.5\% and 0.81\%, respectively, we do not observe evidence for 3+1 neutrinos with neutrino decay. The best-fit parameters for the sterile neutrino with decay model from this study are $Δm_{41}^2=6.7^{+3.9}_{-2.5}\,\rm{eV}^2$, $\sin^2 2θ_{24}=0.33^{+0.20}_{-0.17}$, and $g^2=2.5π\pm1.5π$, where $g$ is the decay-mediating coupling. The preferred regions from short-baseline oscillation searches are excluded at 90\% C.L.
△ Less
Submitted 1 April, 2022;
originally announced April 2022.
-
CMEs and SEPs During November-December 2020: A Challenge for Real-Time Space Weather Forecasting
Authors:
Erika Palmerio,
Christina O. Lee,
M. Leila Mays,
Janet G. Luhmann,
David Lario,
Beatriz Sánchez-Cano,
Ian G. Richardson,
Rami Vainio,
Michael L. Stevens,
Christina M. S. Cohen,
Konrad Steinvall,
Christian Möstl,
Andreas J. Weiss,
Teresa Nieves-Chinchilla,
Yan Li,
Davin E. Larson,
Daniel Heyner,
Stuart D. Bale,
Antoinette B. Galvin,
Mats Holmström,
Yuri V. Khotyaintsev,
Milan Maksimovic,
Igor G. Mitrofanov
Abstract:
Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric sy…
▽ More
Predictions of coronal mass ejections (CMEs) and solar energetic particles (SEPs) are a central issue in space weather forecasting. In recent years, interest in space weather predictions has expanded to include impacts at other planets beyond Earth as well as spacecraft scattered throughout the heliosphere. In this sense, the scope of space weather science now encompasses the whole heliospheric system, and multi-point measurements of solar transients can provide useful insights and validations for prediction models. In this work, we aim to analyse the whole inner heliospheric context between two eruptive flares that took place in late 2020, i.e. the M4.4 flare of November 29 and the C7.4 flare of December 7. This period is especially interesting because the STEREO-A spacecraft was located ~60° east of the Sun-Earth line, giving us the opportunity to test the capabilities of "predictions at 360°" using remote-sensing observations from the Lagrange L1 and L5 points as input. We simulate the CMEs that were ejected during our period of interest and the SEPs accelerated by their shocks using the WSA-Enlil-SEPMOD modelling chain and four sets of input parameters, forming a "mini-ensemble". We validate our results using in-situ observations at six locations, including Earth and Mars. We find that, despite some limitations arising from the models' architecture and assumptions, CMEs and shock-accelerated SEPs can be reasonably studied and forecast in real time at least out to several tens of degrees away from the eruption site using the prediction tools employed here.
△ Less
Submitted 5 May, 2022; v1 submitted 30 March, 2022;
originally announced March 2022.
-
Longitudinal Fairness with Censorship
Authors:
Wenbin Zhang,
Jeremy C. Weiss
Abstract:
Recent works in artificial intelligence fairness attempt to mitigate discrimination by proposing constrained optimization programs that achieve parity for some fairness statistic. Most assume availability of the class label, which is impractical in many real-world applications such as precision medicine, actuarial analysis and recidivism prediction. Here we consider fairness in longitudinal right-…
▽ More
Recent works in artificial intelligence fairness attempt to mitigate discrimination by proposing constrained optimization programs that achieve parity for some fairness statistic. Most assume availability of the class label, which is impractical in many real-world applications such as precision medicine, actuarial analysis and recidivism prediction. Here we consider fairness in longitudinal right-censored environments, where the time to event might be unknown, resulting in censorship of the class label and inapplicability of existing fairness studies. We devise applicable fairness measures, propose a debiasing algorithm, and provide necessary theoretical constructs to bridge fairness with and without censorship for these important and socially-sensitive tasks. Our experiments on four censored datasets confirm the utility of our approach.
△ Less
Submitted 30 March, 2022; v1 submitted 29 March, 2022;
originally announced March 2022.
-
Common Coil Dipole for High Field Magnet Design and R&D
Authors:
Ramesh Gupta,
Kathleen Amm,
Julien Avronsart,
Michael Anerella,
Anis Ben Yahia,
John Cozzolino,
Piyush Joshi,
Mithlesh Kumar,
Febin Kurian,
Chris Runyan,
William Sampson,
Jesse Schmalzle,
Stephan Kahn,
Ronald Scanlan,
Robert Weggel,
Erich Willen,
Qingjin Xu,
Javier Munilla,
Fernando Toral,
Paolo Ferracin,
Steve Gourlay,
GianLuca Sabbi,
Xiaorong Wang,
Danko van der Laan,
Jeremy Weiss
Abstract:
The common coil geometry provides an alternate design to the conventional cosine theta dipoles. It allows a wider range of conductor and magnet technologies. It also facilitates a low-cost, rapid-turn-around design and R&D program. Recent studies carried out as a part of the US Magnet Development Program revealed that at high fields (20 T with 15% operating margin or more), the common coil design…
▽ More
The common coil geometry provides an alternate design to the conventional cosine theta dipoles. It allows a wider range of conductor and magnet technologies. It also facilitates a low-cost, rapid-turn-around design and R&D program. Recent studies carried out as a part of the US Magnet Development Program revealed that at high fields (20 T with 15% operating margin or more), the common coil design also uses significantly less conductor (particularly much less HTS), as compared to that in the other designs.
△ Less
Submitted 16 March, 2022;
originally announced March 2022.
-
Low Energy Event Reconstruction in IceCube DeepCore
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus
, et al. (360 additional authors not shown)
Abstract:
The reconstruction of event-level information, such as the direction or energy of a neutrino interacting in IceCube DeepCore, is a crucial ingredient to many physics analyses. Algorithms to extract this high level information from the detector's raw data have been successfully developed and used for high energy events. In this work, we address unique challenges associated with the reconstruction o…
▽ More
The reconstruction of event-level information, such as the direction or energy of a neutrino interacting in IceCube DeepCore, is a crucial ingredient to many physics analyses. Algorithms to extract this high level information from the detector's raw data have been successfully developed and used for high energy events. In this work, we address unique challenges associated with the reconstruction of lower energy events in the range of a few to hundreds of GeV and present two separate, state-of-the-art algorithms. One algorithm focuses on the fast directional reconstruction of events based on unscattered light. The second algorithm is a likelihood-based multipurpose reconstruction offering superior resolutions, at the expense of larger computational cost.
△ Less
Submitted 4 March, 2022;
originally announced March 2022.
-
Search for High-Energy Neutrino Emission from Galactic X-ray Binaries with IceCube
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus
, et al. (358 additional authors not shown)
Abstract:
We present the first comprehensive search for high-energy neutrino emission from high- and low-mass X-ray binaries conducted by IceCube. Galactic X-ray binaries are long-standing candidates for the source of Galactic hadronic cosmic rays and neutrinos. The compact object in these systems can be the site of cosmic-ray acceleration, and neutrinos can be produced by interactions of cosmic rays with r…
▽ More
We present the first comprehensive search for high-energy neutrino emission from high- and low-mass X-ray binaries conducted by IceCube. Galactic X-ray binaries are long-standing candidates for the source of Galactic hadronic cosmic rays and neutrinos. The compact object in these systems can be the site of cosmic-ray acceleration, and neutrinos can be produced by interactions of cosmic rays with radiation or gas, in the jet of a microquasar, in the stellar wind, or in the atmosphere of the companion star. We study X-ray binaries using 7.5 years of IceCube data with three separate analyses. In the first, we search for periodic neutrino emission from 55 binaries in the Northern Sky with known orbital periods. In the second, the X-ray light curves of 102 binaries across the entire sky are used as templates to search for time-dependent neutrino emission. Finally, we search for time-integrated emission of neutrinos for a list of 4 notable binaries identified as microquasars. In the absence of a significant excess, we place upper limits on the neutrino flux for each hypothesis and compare our results with theoretical predictions for several binaries. In addition, we evaluate the sensitivity of the next generation neutrino telescope at the South Pole, IceCube-Gen2, and demonstrate its power to identify potential neutrino emission from these binary sources in the Galaxy.
△ Less
Submitted 23 February, 2022;
originally announced February 2022.
-
Writhed Analytical Magnetic Flux Rope Model
Authors:
Andreas J. Weiss,
Teresa Nieves-Chinchilla,
Christian Möstl,
Martin A. Reiss,
Tanja Amerstorfer,
Rachel L. Bailey
Abstract:
Observations of magnetic clouds, within interplanetary coronal mass ejections (ICMEs), are often well described by flux rope models. Most of these assume either a cylindrical or toroidal geometry. In some cases, these models are also capable of accounting for non-axisymmetric crosssections but they generally all assume axial invariance. It can be expected that any ICME, and its flux rope, will be…
▽ More
Observations of magnetic clouds, within interplanetary coronal mass ejections (ICMEs), are often well described by flux rope models. Most of these assume either a cylindrical or toroidal geometry. In some cases, these models are also capable of accounting for non-axisymmetric crosssections but they generally all assume axial invariance. It can be expected that any ICME, and its flux rope, will be deformed along its axis due to influences such as the solar wind. In this work, we aim to develop a writhed analytical magnetic flux rope model which would allow us to analytically describe a flux rope structure with varying curvature and torsion so that we are no longer constrained to a cylindrical or toroidal geometry. In this first iteration of our model we will solely focus on a circular cross-section of constant size. We describe our flux rope geometry in terms of a parametrized flux rope axis and a parallel transport frame. We derive expressions for the axial and poloidal magnetic field components under the assumption that the total axial magnetic flux is conserved. We find an entire class of possible solutions, which differ by the choice of integration constants, and present the results for a specific example. In general, we find that the twist of the magnetic field locally changes when the geometry deviates from a cylinder or torus. This new approach also allows us to generate completely new types of in situ magnetic field profiles which strongly deviate from those generated by cylindrical or toroidal models.
△ Less
Submitted 2 December, 2022; v1 submitted 21 February, 2022;
originally announced February 2022.
-
Fairness Amidst Non-IID Graph Data: Current Achievements and Future Directions
Authors:
Wenbin Zhang,
Shimei Pan,
Shuigeng Zhou,
Toby Walsh,
Jeremy C. Weiss
Abstract:
The importance of understanding and correcting algorithmic bias in machine learning (ML) has led to an increase in research on fairness in ML, which typically assumes that the underlying data is independent and identically distributed (IID). However, in reality, data is often represented using non-IID graph structures that capture connections among individual units. To address bias in ML systems,…
▽ More
The importance of understanding and correcting algorithmic bias in machine learning (ML) has led to an increase in research on fairness in ML, which typically assumes that the underlying data is independent and identically distributed (IID). However, in reality, data is often represented using non-IID graph structures that capture connections among individual units. To address bias in ML systems, it is crucial to bridge the gap between the traditional fairness literature designed for IID data and the ubiquity of non-IID graph data. In this survey, we review such recent advance in fairness amidst non-IID graph data and identify datasets and evaluation metrics available for future research. We also point out the limitations of existing work as well as promising future directions.
△ Less
Submitted 9 February, 2023; v1 submitted 14 February, 2022;
originally announced February 2022.
-
Deep Learning-based Assessment of Hepatic Steatosis on chest CT
Authors:
Zhongyi Zhang,
Jakob Weiss,
Jana Taron,
Roman Zeleznik,
Michael T. Lu,
Hugo J. W. L. Aerts
Abstract:
Purpose: Automatic methods are required for the early detection of hepatic steatosis to avoid progression to cirrhosis and cancer. Here, we developed a fully automated deep learning pipeline to quantify hepatic steatosis on non-contrast enhanced chest computed tomography (CT) scans. Materials and Methods: We developed and evaluated our pipeline on chest CT images of 1,431 randomly selected Nationa…
▽ More
Purpose: Automatic methods are required for the early detection of hepatic steatosis to avoid progression to cirrhosis and cancer. Here, we developed a fully automated deep learning pipeline to quantify hepatic steatosis on non-contrast enhanced chest computed tomography (CT) scans. Materials and Methods: We developed and evaluated our pipeline on chest CT images of 1,431 randomly selected National Lung Screening Trial (NLST) participants. A dataset of 451 CT scans with volumetric liver segmentations of expert readers was used for training a deep learning model. For testing, in an independent dataset of 980 CT scans hepatic attenuation was manually measured by an expert reader on three cross-sectional images at different hepatic levels by selecting three circular regions of interest. Additionally, 100 randomly selected cases of the test set were volumetrically segmented by expert readers. Hepatic steatosis on the test set was defined as mean hepatic attenuation of < 40 Hounsfield unit. Spearman correlation was conducted to analyze liver fat quantification accuracy and the Cohen's Kappa coefficient was calculated for hepatic steatosis prediction reliability. Results: Our pipeline demonstrated strong performance and achieved a mean dice score of 0.970 for the volumetric liver segmentation. The spearman correlation of the liver fat quantification was 0.954 (P <0.0001) between the automated and expert reader measurements. The cohen's kappa coefficient was 0.875 for automatic assessment of hepatic steatosis. Conclusion: We developed a fully automatic deep learning-based pipeline for the assessment of hepatic steatosis in chest CT images. With the fast and cheap screening of hepatic steatosis, our pipeline has the potential to help initiate preventive measures to avoid progression to cirrhosis and cancer.
△ Less
Submitted 4 February, 2022;
originally announced February 2022.
-
A Computational Framework for Atrioventricular Valve Modeling using Open-Source Software
Authors:
Wensi Wu,
Stephen Ching,
Steve A. Maas,
Andras Lasso,
Patricia Sabin,
Jeffrey A. Weiss,
Matthew A. Jolley
Abstract:
Atrioventricular valve regurgitation is a significant cause of morbidity and mortality in patients with acquired and congenital cardiac valve disease. Image-derived computational modeling of atrioventricular valves has advanced substantially over the last decade and holds particular promise to inform valve repair in small and heterogeneous populations which are less likely to be optimized through…
▽ More
Atrioventricular valve regurgitation is a significant cause of morbidity and mortality in patients with acquired and congenital cardiac valve disease. Image-derived computational modeling of atrioventricular valves has advanced substantially over the last decade and holds particular promise to inform valve repair in small and heterogeneous populations which are less likely to be optimized through empiric clinical application. While an abundance of computational biomechanics studies have investigated mitral and tricuspid valve disease in adults, few studies have investigated application to vulnerable pediatric and congenital heart populations. Further, to date, investigators have primarily relied upon a series of commercial applications that are neither designed for image-derived modeling of cardiac valves, nor freely available to facilitate transparent and reproducible valve science. To address this deficiency, we aimed to build an open-source computational framework for the image-derived biomechanical analysis of atrioventricular valves. In the present work, we integrated an open-source valve modeling platform, SlicerHeart, and an open-source biomechanics finite element modeling software, FEBio, to facilitate image-derived atrioventricular valve model creation and finite element analysis. We present a detailed verification and sensitivity analysis to demonstrate the fidelity of this modeling in application to 3D echocardiography-derived pediatric mitral and tricuspid valve models. Our analyses achieved excellent agreement with those reported in the literature. As such, this evolving computational framework offers a promising initial foundation for future development and investigation of valve mechanics, in particular collaborative efforts targeting the development of improved repairs for children with congenital heart disease.
△ Less
Submitted 31 January, 2022;
originally announced January 2022.
-
Density of GeV muons in air showers measured with IceTop
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker,
J. Becker Tjus
, et al. (355 additional authors not shown)
Abstract:
We present a measurement of the density of GeV muons in near-vertical air showers using three years of data recorded by the IceTop array at the South Pole. Depending on the shower size, the muon densities have been measured at lateral distances between 200 m and 1000 m. From these lateral distributions, we derive the muon densities as functions of energy at reference distances of 600 m and 800 m f…
▽ More
We present a measurement of the density of GeV muons in near-vertical air showers using three years of data recorded by the IceTop array at the South Pole. Depending on the shower size, the muon densities have been measured at lateral distances between 200 m and 1000 m. From these lateral distributions, we derive the muon densities as functions of energy at reference distances of 600 m and 800 m for primary energies between 2.5 PeV and 40 PeV and between 9 PeV and 120 PeV, respectively. The muon densities are determined using, as a baseline, the hadronic interaction model Sibyll 2.1 together with various composition models. The measurements are consistent with the predicted muon densities within these baseline interaction and composition models. The measured muon densities have also been compared to simulations using the post-LHC models EPOS-LHC and QGSJet-II.04. The result of this comparison is that the post-LHC models together with any given composition model yield higher muon densities than observed. This is in contrast to the observations above 1 EeV where all model simulations yield for any mass composition lower muon densities than the measured ones. The post-LHC models in general feature higher muon densities so that the agreement with experimental data at the highest energies is improved but the muon densities are not correct in the energy range between 2.5 PeV and about 100 PeV.
△ Less
Submitted 18 May, 2022; v1 submitted 29 January, 2022;
originally announced January 2022.
-
Volume Transport by a 3D Quasigeostrophic Heton
Authors:
Adhithiya Sivakumar,
Jeffrey B. Weiss
Abstract:
Oceanic flows self-organize into coherent vortices which strongly influence their transport and mixing properties. Counter-rotating vortex pairs can travel long distances and carry trapped fluid as they move. These structures are often modeled as hetons, viz. counter-rotating quasigeostrophic point vortex pairs with equal circulations. Here, we investigate the structure of the transport induced by…
▽ More
Oceanic flows self-organize into coherent vortices which strongly influence their transport and mixing properties. Counter-rotating vortex pairs can travel long distances and carry trapped fluid as they move. These structures are often modeled as hetons, viz. counter-rotating quasigeostrophic point vortex pairs with equal circulations. Here, we investigate the structure of the transport induced by a single three-dimensional heton. The transport is determined by the Hamiltonian structure of the velocity field induced by the heton's component vortices. The dynamics displays a sequence of bifurcations as one moves through the heton-induced velocity field in height. These bifurcations create and destroy unstable fixed points whose associated invariant manifolds bound the trapped volume. Heton configurations fall into three categories. Vertically aligned hetons do not move and do not transport fluid. Horizontally aligned hetons have a single parameter, the horizontal vortex half-separation $Y$, and simple scaling shows the dimensional trapped volume scales as $Y^3$. Tilted hetons are described by two parameters, $Y$ and the vertical vortex half-separation $Z$, rendering the scaling analysis more complex. A scaling theory is developed for the trapped volume of tilted hetons showing that it scales as $Z^4/Y$ for large $Z$. Numerical calculations illustrate the structure of the trapped volume and verify the scaling theory.
△ Less
Submitted 20 January, 2022;
originally announced January 2022.
-
Search for Spatial Correlations of Neutrinos with Ultra-High-Energy Cosmic Rays
Authors:
The ANTARES collaboration,
A. Albert,
S. Alves,
M. André,
M. Anghinolfi,
M. Ardid,
S. Ardid,
J. -J. Aubert,
J. Aublin,
B. Baret,
S. Basa,
B. Belhorma,
M. Bendahman,
V. Bertin,
S. Biagi,
M. Bissinger,
J. Boumaaza,
M. Bouta,
M. C. Bouwhuis,
H. Brânzaş,
R. Bruijn,
J. Brunner,
J. Busto,
B. Caiffi,
D. Calvo
, et al. (1025 additional authors not shown)
Abstract:
For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for corre…
▽ More
For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrinos clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.
△ Less
Submitted 23 August, 2022; v1 submitted 18 January, 2022;
originally announced January 2022.
-
Strong constraints on neutrino nonstandard interactions from TeV-scale $ν_μ$ disappearance at IceCube
Authors:
IceCube Collaboration,
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty,
K. -H. Becker
, et al. (359 additional authors not shown)
Abstract:
We report a search for nonstandard neutrino interactions (NSI) using eight years of TeV-scale atmospheric muon neutrino data from the IceCube Neutrino Observatory. By reconstructing incident energies and zenith angles for atmospheric neutrino events, this analysis presents unified confidence intervals for the NSI parameter $ε_{μτ}$. The best-fit value is consistent with no NSI at a p-value of 25.2…
▽ More
We report a search for nonstandard neutrino interactions (NSI) using eight years of TeV-scale atmospheric muon neutrino data from the IceCube Neutrino Observatory. By reconstructing incident energies and zenith angles for atmospheric neutrino events, this analysis presents unified confidence intervals for the NSI parameter $ε_{μτ}$. The best-fit value is consistent with no NSI at a p-value of 25.2%. With a 90% confidence interval of $-0.0041 \leq ε_{μτ} \leq 0.0031$ along the real axis and similar strength in the complex plane, this result is the strongest constraint on any NSI parameter from any oscillation channel to date.
△ Less
Submitted 5 June, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
Improved Characterization of the Astrophysical Muon-Neutrino Flux with 9.5 Years of IceCube Data
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
C. Alispach,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
A. Barbano,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty
, et al. (355 additional authors not shown)
Abstract:
We present a measurement of the high-energy astrophysical muon-neutrino flux with the IceCube Neutrino Observatory. The measurement uses a high-purity selection of ~650k neutrino-induced muon tracks from the Northern celestial hemisphere, corresponding to 9.5 years of experimental data. With respect to previous publications, the measurement is improved by the increased size of the event sample and…
▽ More
We present a measurement of the high-energy astrophysical muon-neutrino flux with the IceCube Neutrino Observatory. The measurement uses a high-purity selection of ~650k neutrino-induced muon tracks from the Northern celestial hemisphere, corresponding to 9.5 years of experimental data. With respect to previous publications, the measurement is improved by the increased size of the event sample and the extended model testing beyond simple power-law hypotheses. An updated treatment of systematic uncertainties and atmospheric background fluxes has been implemented based on recent models. The best-fit single power-law parameterization for the astrophysical energy spectrum results in a normalization of $φ_{\mathrm{@100TeV}}^{ν_μ+\barν_μ} = 1.44_{-0.26}^{+0.25} \times 10^{-18}\,\mathrm{GeV}^{-1}\mathrm{cm}^{-2}\mathrm{s}^{-1}\mathrm{sr}^{-1}$ and a spectral index $γ_{\mathrm{SPL}} = 2.37_{-0.09}^{+0.09}$, constrained in the energy range from $15\,\mathrm{TeV}$ to $5\,\mathrm{PeV}$. The model tests include a single power law with a spectral cutoff at high energies, a log-parabola model, several source-class specific flux predictions from the literature and a model-independent spectral unfolding. The data is well consistent with a single power law hypothesis, however, spectra with softening above one PeV are statistically more favorable at a two sigma level.
△ Less
Submitted 19 November, 2021;
originally announced November 2021.
-
A search for neutrino emission from cores of Active Galactic Nuclei
Authors:
R. Abbasi,
M. Ackermann,
J. Adams,
J. A. Aguilar,
M. Ahlers,
M. Ahrens,
J. M. Alameddine,
C. Alispach,
A. A. Alves Jr.,
N. M. Amin,
K. Andeen,
T. Anderson,
G. Anton,
C. Argüelles,
Y. Ashida,
S. Axani,
X. Bai,
A. Balagopal V.,
A. Barbano,
S. W. Barwick,
B. Bastian,
V. Basu,
S. Baur,
R. Bay,
J. J. Beatty
, et al. (355 additional authors not shown)
Abstract:
The sources of the majority of the high-energy astrophysical neutrinos observed with the IceCube neutrino telescope at the South Pole are unknown. So far, only a gamma-ray blazar was compellingly associated with the emission of high-energy neutrinos. In addition, several studies suggest that the neutrino emission from the gamma-ray blazar population only accounts for a small fraction of the total…
▽ More
The sources of the majority of the high-energy astrophysical neutrinos observed with the IceCube neutrino telescope at the South Pole are unknown. So far, only a gamma-ray blazar was compellingly associated with the emission of high-energy neutrinos. In addition, several studies suggest that the neutrino emission from the gamma-ray blazar population only accounts for a small fraction of the total astrophysical neutrino flux. In this work we probe the production of high-energy neutrinos in the cores of Active Galactic Nuclei (AGN), induced by accelerated cosmic rays in the accretion disk region. We present a likelihood analysis based on eight years of IceCube data, searching for a cumulative neutrino signal from three AGN samples created for this work. The neutrino emission is assumed to be proportional to the accretion disk luminosity estimated from the soft X-ray flux. Next to the observed soft X-ray flux, the objects for the three samples have been selected based on their radio emission and infrared color properties. For the largest sample in this search, an excess of high-energy neutrino events with respect to an isotropic background of atmospheric and astrophysical neutrinos is found, corresponding to a post-trial significance of 2.60$σ$. Assuming a power-law spectrum, the best-fit spectral index is $2.03^{+0.14}_{-0.11}$, consistent with expectations from particle acceleration in astrophysical sources. If interpreted as a genuine signal with the assumptions of a proportionality of X-ray and neutrino fluxes and a model for the sub-threshold flux distribution, this observation implies that at 100 TeV, 27$\%$ - 100$\%$ of the observed neutrinos arise from particle acceleration in the core of AGN.
△ Less
Submitted 19 November, 2021;
originally announced November 2021.