-
Experimental timing and control using microcontrollers
Authors:
Philip T Starkey,
Carter Turnbaugh,
Patrick Miller,
Kermit-James LeBlanc,
David H Meyer
Abstract:
Modern physics experiments rely on precise timing provided by programmable digital pulse generators. In many experimental control systems, this role is filled by custom devices built on Field Programmable Gate Arrays (FPGAs). While highly flexible and performant, these systems can be difficult to scale to very large systems due to cost and complexity. Recent advances in microcontroller systems all…
▽ More
Modern physics experiments rely on precise timing provided by programmable digital pulse generators. In many experimental control systems, this role is filled by custom devices built on Field Programmable Gate Arrays (FPGAs). While highly flexible and performant, these systems can be difficult to scale to very large systems due to cost and complexity. Recent advances in microcontroller systems allows these much simpler systems to fill the role of digital pulse generators. Here we demonstrate one such alternative based on the Raspberry Pi Pico microcontroller which allows for timing resolution down to 7.5 ns with a minimum pulse width of of 37.5 ns.
△ Less
Submitted 25 June, 2024;
originally announced June 2024.
-
Differential fuzz testing to detect tampering in sensor systems and its application to arms control authentication
Authors:
Jayson R Vavrek,
Luozhong Zhou,
Joshua Boverhof,
Elisa R Heymann,
Barton P Miller,
Sean Peisert
Abstract:
In future nuclear arms control treaties, it will be necessary to authenticate the hardware and software components of verification measurement systems, i.e., to ensure these systems are functioning as intended and have not been tampered with by malicious actors. While methods such as source code hashing and static analysis can help verify the integrity of software components, they may not be capab…
▽ More
In future nuclear arms control treaties, it will be necessary to authenticate the hardware and software components of verification measurement systems, i.e., to ensure these systems are functioning as intended and have not been tampered with by malicious actors. While methods such as source code hashing and static analysis can help verify the integrity of software components, they may not be capable of detecting tampering with environment variables, external libraries, or the firmware and hardware of radiation measurement systems. In this article, we introduce the concept of physical differential fuzz testing as a challenge-response-style tamper indicator that can holistically and simultaneously test all the above components in a cyber-physical system. In essence, we randomly sample (or "fuzz") the untampered system's parameter space, including both normal and off-normal parameter values, and consider the time series of outputs as the baseline signature of the system. Re-running the same input sequence on a untampered system will produce an output sequence consistent with this baseline, while running the same input sequence on a tampered system will produce a modified output sequence and raise an alarm. We then apply this concept to authenticating the radiation measurement equipment in nuclear weapon verification systems and conduct demonstration fuzz testing measurements with a sodium iodide (NaI) gamma ray spectrometer. Because there is Poisson noise in the measured output spectra, we also use a mechanism for comparing inherently noisy or stochastic fuzzing sequences. We show that physical differential fuzz testing can detect two types of tamper attempts, and conclude that it is a promising framework for authenticating future cyber-physical systems in nuclear arms control, safeguards, and beyond.
△ Less
Submitted 8 April, 2024;
originally announced April 2024.
-
DROID: A Large-Scale In-The-Wild Robot Manipulation Dataset
Authors:
Alexander Khazatsky,
Karl Pertsch,
Suraj Nair,
Ashwin Balakrishna,
Sudeep Dasari,
Siddharth Karamcheti,
Soroush Nasiriany,
Mohan Kumar Srirama,
Lawrence Yunliang Chen,
Kirsty Ellis,
Peter David Fagan,
Joey Hejna,
Masha Itkina,
Marion Lepert,
Yecheng Jason Ma,
Patrick Tree Miller,
Jimmy Wu,
Suneel Belkhale,
Shivin Dass,
Huy Ha,
Arhan Jain,
Abraham Lee,
Youngwoon Lee,
Marius Memmel,
Sungjae Park
, et al. (74 additional authors not shown)
Abstract:
The creation of large, diverse, high-quality robot manipulation datasets is an important stepping stone on the path toward more capable and robust robotic manipulation policies. However, creating such datasets is challenging: collecting robot manipulation data in diverse environments poses logistical and safety challenges and requires substantial investments in hardware and human labour. As a resu…
▽ More
The creation of large, diverse, high-quality robot manipulation datasets is an important stepping stone on the path toward more capable and robust robotic manipulation policies. However, creating such datasets is challenging: collecting robot manipulation data in diverse environments poses logistical and safety challenges and requires substantial investments in hardware and human labour. As a result, even the most general robot manipulation policies today are mostly trained on data collected in a small number of environments with limited scene and task diversity. In this work, we introduce DROID (Distributed Robot Interaction Dataset), a diverse robot manipulation dataset with 76k demonstration trajectories or 350 hours of interaction data, collected across 564 scenes and 84 tasks by 50 data collectors in North America, Asia, and Europe over the course of 12 months. We demonstrate that training with DROID leads to policies with higher performance and improved generalization ability. We open source the full dataset, policy learning code, and a detailed guide for reproducing our robot hardware setup.
△ Less
Submitted 19 March, 2024;
originally announced March 2024.
-
Hard-label based Small Query Black-box Adversarial Attack
Authors:
Jeonghwan Park,
Paul Miller,
Niall McLaughlin
Abstract:
We consider the hard label based black box adversarial attack setting which solely observes predicted classes from the target model. Most of the attack methods in this setting suffer from impractical number of queries required to achieve a successful attack. One approach to tackle this drawback is utilising the adversarial transferability between white box surrogate models and black box target mod…
▽ More
We consider the hard label based black box adversarial attack setting which solely observes predicted classes from the target model. Most of the attack methods in this setting suffer from impractical number of queries required to achieve a successful attack. One approach to tackle this drawback is utilising the adversarial transferability between white box surrogate models and black box target model. However, the majority of the methods adopting this approach are soft label based to take the full advantage of zeroth order optimisation. Unlike mainstream methods, we propose a new practical setting of hard label based attack with an optimisation process guided by a pretrained surrogate model. Experiments show the proposed method significantly improves the query efficiency of the hard label based black-box attack across various target model architectures. We find the proposed method achieves approximately 5 times higher attack success rate compared to the benchmarks, especially at the small query budgets as 100 and 250.
△ Less
Submitted 9 March, 2024;
originally announced March 2024.
-
Detailed Report on the Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm
Authors:
D. P. Aguillard,
T. Albahri,
D. Allspach,
A. Anisenkov,
K. Badgley,
S. Baeßler,
I. Bailey,
L. Bailey,
V. A. Baranov,
E. Barlas-Yucel,
T. Barrett,
E. Barzi,
F. Bedeschi,
M. Berz,
M. Bhattacharya,
H. P. Binney,
P. Bloom,
J. Bono,
E. Bottalico,
T. Bowcock,
S. Braun,
M. Bressler,
G. Cantatore,
R. M. Carey,
B. C. K. Casey
, et al. (168 additional authors not shown)
Abstract:
We present details on a new measurement of the muon magnetic anomaly, $a_μ= (g_μ-2)/2$. The result is based on positive muon data taken at Fermilab's Muon Campus during the 2019 and 2020 accelerator runs. The measurement uses $3.1$ GeV$/c$ polarized muons stored in a $7.1$-m-radius storage ring with a $1.45$ T uniform magnetic field. The value of $ a_μ$ is determined from the measured difference b…
▽ More
We present details on a new measurement of the muon magnetic anomaly, $a_μ= (g_μ-2)/2$. The result is based on positive muon data taken at Fermilab's Muon Campus during the 2019 and 2020 accelerator runs. The measurement uses $3.1$ GeV$/c$ polarized muons stored in a $7.1$-m-radius storage ring with a $1.45$ T uniform magnetic field. The value of $ a_μ$ is determined from the measured difference between the muon spin precession frequency and its cyclotron frequency. This difference is normalized to the strength of the magnetic field, measured using Nuclear Magnetic Resonance (NMR). The ratio is then corrected for small contributions from beam motion, beam dispersion, and transient magnetic fields. We measure $a_μ= 116 592 057 (25) \times 10^{-11}$ (0.21 ppm). This is the world's most precise measurement of this quantity and represents a factor of $2.2$ improvement over our previous result based on the 2018 dataset. In combination, the two datasets yield $a_μ(\text{FNAL}) = 116 592 055 (24) \times 10^{-11}$ (0.20 ppm). Combining this with the measurements from Brookhaven National Laboratory for both positive and negative muons, the new world average is $a_μ$(exp) $ = 116 592 059 (22) \times 10^{-11}$ (0.19 ppm).
△ Less
Submitted 22 May, 2024; v1 submitted 23 February, 2024;
originally announced February 2024.
-
On Strong Zero-Dispersion Asymptotics for Benjamin-Ono Soliton Ensembles
Authors:
Elliot Blackstone,
Louise Gassot,
Peter D. Miller
Abstract:
A soliton ensemble is a particular kind of approximation of the solution of an initial-value problem for an integrable equation by a reflectionless potential that is well adapted to singular asymptotics like the small-dispersion limit. We show how soliton ensembles for the Benjamin-Ono equation can be analyzed in this limit via the construction of local approximations that capture highly oscillato…
▽ More
A soliton ensemble is a particular kind of approximation of the solution of an initial-value problem for an integrable equation by a reflectionless potential that is well adapted to singular asymptotics like the small-dispersion limit. We show how soliton ensembles for the Benjamin-Ono equation can be analyzed in this limit via the construction of local approximations that capture highly oscillatory features of the solution and hence provide more information than weak convergence results that are easier to obtain. These local approximations are deduced from the distributions of eigenvalues of two related matrices, one Hermitian and another non-Hermitian. We perform careful numerical experiments to deduce the asymptotic behavior of the eigenvalues of these matrices in the small-dispersion limit, and formulate conjectures reflecting our observations. Then we apply the conjectures to construct the local approximations of slowly varying profiles and rapidly oscillating profiles as well. We show that the latter profiles are consistent with the predictions of Whitham modulation theory as originally developed for the Benjamin-Ono equation by Dobrokhotov and Krichever.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
Open X-Embodiment: Robotic Learning Datasets and RT-X Models
Authors:
Open X-Embodiment Collaboration,
Abby O'Neill,
Abdul Rehman,
Abhinav Gupta,
Abhiram Maddukuri,
Abhishek Gupta,
Abhishek Padalkar,
Abraham Lee,
Acorn Pooley,
Agrim Gupta,
Ajay Mandlekar,
Ajinkya Jain,
Albert Tung,
Alex Bewley,
Alex Herzog,
Alex Irpan,
Alexander Khazatsky,
Anant Rai,
Anchit Gupta,
Andrew Wang,
Andrey Kolobov,
Anikait Singh,
Animesh Garg,
Aniruddha Kembhavi,
Annie Xie
, et al. (267 additional authors not shown)
Abstract:
Large, high-capacity models trained on diverse datasets have shown remarkable successes on efficiently tackling downstream applications. In domains from NLP to Computer Vision, this has led to a consolidation of pretrained models, with general pretrained backbones serving as a starting point for many applications. Can such a consolidation happen in robotics? Conventionally, robotic learning method…
▽ More
Large, high-capacity models trained on diverse datasets have shown remarkable successes on efficiently tackling downstream applications. In domains from NLP to Computer Vision, this has led to a consolidation of pretrained models, with general pretrained backbones serving as a starting point for many applications. Can such a consolidation happen in robotics? Conventionally, robotic learning methods train a separate model for every application, every robot, and even every environment. Can we instead train generalist X-robot policy that can be adapted efficiently to new robots, tasks, and environments? In this paper, we provide datasets in standardized data formats and models to make it possible to explore this possibility in the context of robotic manipulation, alongside experimental results that provide an example of effective X-robot policies. We assemble a dataset from 22 different robots collected through a collaboration between 21 institutions, demonstrating 527 skills (160266 tasks). We show that a high-capacity model trained on this data, which we call RT-X, exhibits positive transfer and improves the capabilities of multiple robots by leveraging experience from other platforms. More details can be found on the project website https://robotics-transformer-x.github.io.
△ Less
Submitted 1 June, 2024; v1 submitted 13 October, 2023;
originally announced October 2023.
-
Large Eddy Simulation of Non-stationary Hurricane Boundary Layer Winds
Authors:
Tianqi Ma,
Chao Sun,
Paul Miller
Abstract:
Recent extreme tropical cyclones caused extensive damages to critical civil infrastructure globally. To better capture the unique hurricane wind characteristics, a large eddy simulation (LES) Hurricane Boundary Layer (HBL) model is developed by considering the variation of meso-scale kinematic and thermodynamic conditions. An asymmetric model is adopted to obtain the gradient wind velocity using t…
▽ More
Recent extreme tropical cyclones caused extensive damages to critical civil infrastructure globally. To better capture the unique hurricane wind characteristics, a large eddy simulation (LES) Hurricane Boundary Layer (HBL) model is developed by considering the variation of meso-scale kinematic and thermodynamic conditions. An asymmetric model is adopted to obtain the gradient wind velocity using the National Hurricane Center data. The meso-scale thermal conditions are obtained by extracting the hourly air temperature and relative humidity profiles from generated proxy soundings. Measurements recorded at the Aransas County airport during Hurricane Harvey and that at the City of Naples during Irma are used to validate the developed LES model. Research results show that the simulated 10-minute average wind speed and direction are consistent with the observations. The developed model can well predict the high wind turbulences, which are around 20% in Hurricane Harvey and 26% in Hurricane Irma. The 3-s gust wind speeds reach 62.4 m/s at 10-m elevation during Hurricane Harvey and 53.5 m/s at 15-m elevation during Hurricane Irma, close to the field observed data of 61.3 m/s and 54.2 m/s, respectively. The simulated 3-s gust factors are close to the observation except at some moments with significant variations because of the poorly understood physical phenomena. The simulated wind spectrum in longitudinal and lateral directions agrees well with the observed results. In summary, the developed LES-based HBL model can capture the main characteristics of hurricane structure and turbulence characteristics and is applicable for modeling civil infrastructure exposed to hurricanes.
△ Less
Submitted 5 October, 2023;
originally announced October 2023.
-
Differential Equations for Approximate Solutions of Painlevé Equations: Application to the Algebraic Solutions of the Painlevé-III $({\rm D}_7)$ Equation
Authors:
Robert J. Buckingham,
Peter D. Miller
Abstract:
It is well known that the Painlevé equations can formally degenerate to autonomous differential equations with elliptic function solutions in suitable scaling limits. A way to make this degeneration rigorous is to apply Deift-Zhou steepest-descent techniques to a Riemann-Hilbert representation of a family of solutions. This method leads to an explicit approximation formula in terms of theta functi…
▽ More
It is well known that the Painlevé equations can formally degenerate to autonomous differential equations with elliptic function solutions in suitable scaling limits. A way to make this degeneration rigorous is to apply Deift-Zhou steepest-descent techniques to a Riemann-Hilbert representation of a family of solutions. This method leads to an explicit approximation formula in terms of theta functions and related algebro-geometric ingredients that is difficult to directly link to the expected limiting differential equation. However, the approximation arises from an outer parametrix that satisfies relatively simple conditions. By applying a method that we learned from Alexander Its, it is possible to use these simple conditions to directly obtain the limiting differential equation, bypassing the details of the algebro-geometric solution of the outer parametrix problem. In this paper, we illustrate the use of this method to relate an approximation of the algebraic solutions of the Painlevé-III (D$_7$) equation valid in the part of the complex plane where the poles and zeros of the solutions asymptotically reside to a form of the Weierstrass equation.
△ Less
Submitted 20 January, 2024; v1 submitted 30 August, 2023;
originally announced August 2023.
-
Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm
Authors:
D. P. Aguillard,
T. Albahri,
D. Allspach,
A. Anisenkov,
K. Badgley,
S. Baeßler,
I. Bailey,
L. Bailey,
V. A. Baranov,
E. Barlas-Yucel,
T. Barrett,
E. Barzi,
F. Bedeschi,
M. Berz,
M. Bhattacharya,
H. P. Binney,
P. Bloom,
J. Bono,
E. Bottalico,
T. Bowcock,
S. Braun,
M. Bressler,
G. Cantatore,
R. M. Carey,
B. C. K. Casey
, et al. (166 additional authors not shown)
Abstract:
We present a new measurement of the positive muon magnetic anomaly, $a_μ\equiv (g_μ- 2)/2$, from the Fermilab Muon $g\!-\!2$ Experiment using data collected in 2019 and 2020. We have analyzed more than 4 times the number of positrons from muon decay than in our previous result from 2018 data. The systematic error is reduced by more than a factor of 2 due to better running conditions, a more stable…
▽ More
We present a new measurement of the positive muon magnetic anomaly, $a_μ\equiv (g_μ- 2)/2$, from the Fermilab Muon $g\!-\!2$ Experiment using data collected in 2019 and 2020. We have analyzed more than 4 times the number of positrons from muon decay than in our previous result from 2018 data. The systematic error is reduced by more than a factor of 2 due to better running conditions, a more stable beam, and improved knowledge of the magnetic field weighted by the muon distribution, $\tildeω'^{}_p$, and of the anomalous precession frequency corrected for beam dynamics effects, $ω_a$. From the ratio $ω_a / \tildeω'^{}_p$, together with precisely determined external parameters, we determine $a_μ= 116\,592\,057(25) \times 10^{-11}$ (0.21 ppm). Combining this result with our previous result from the 2018 data, we obtain $a_μ\text{(FNAL)} = 116\,592\,055(24) \times 10^{-11}$ (0.20 ppm). The new experimental world average is $a_μ(\text{Exp}) = 116\,592\,059(22)\times 10^{-11}$ (0.19 ppm), which represents a factor of 2 improvement in precision.
△ Less
Submitted 4 October, 2023; v1 submitted 11 August, 2023;
originally announced August 2023.
-
Efficient Model Adaptation for Continual Learning at the Edge
Authors:
Zachary A. Daniels,
Jun Hu,
Michael Lomnitz,
Phil Miller,
Aswin Raghavan,
Joe Zhang,
Michael Piacentino,
David Zhang
Abstract:
Most machine learning (ML) systems assume stationary and matching data distributions during training and deployment. This is often a false assumption. When ML models are deployed on real devices, data distributions often shift over time due to changes in environmental factors, sensor characteristics, and task-of-interest. While it is possible to have a human-in-the-loop to monitor for distribution…
▽ More
Most machine learning (ML) systems assume stationary and matching data distributions during training and deployment. This is often a false assumption. When ML models are deployed on real devices, data distributions often shift over time due to changes in environmental factors, sensor characteristics, and task-of-interest. While it is possible to have a human-in-the-loop to monitor for distribution shifts and engineer new architectures in response to these shifts, such a setup is not cost-effective. Instead, non-stationary automated ML (AutoML) models are needed. This paper presents the Encoder-Adaptor-Reconfigurator (EAR) framework for efficient continual learning under domain shifts. The EAR framework uses a fixed deep neural network (DNN) feature encoder and trains shallow networks on top of the encoder to handle novel data. The EAR framework is capable of 1) detecting when new data is out-of-distribution (OOD) by combining DNNs with hyperdimensional computing (HDC), 2) identifying low-parameter neural adaptors to adapt the model to the OOD data using zero-shot neural architecture search (ZS-NAS), and 3) minimizing catastrophic forgetting on previous tasks by progressively growing the neural architecture as needed and dynamically routing data through the appropriate adaptors and reconfigurators for handling domain-incremental and class-incremental continual learning. We systematically evaluate our approach on several benchmark datasets for domain adaptation and demonstrate strong performance compared to state-of-the-art algorithms for OOD detection and few-/zero-shot NAS.
△ Less
Submitted 13 October, 2023; v1 submitted 3 August, 2023;
originally announced August 2023.
-
Painlevé-III Monodromy Maps Under the $D_6\to D_8$ Confluence and Applications to the Large-Parameter Asymptotics of Rational Solutions
Authors:
Ahmad Barhoumi,
Oleg Lisovyy,
Peter D. Miller,
Andrei Prokhorov
Abstract:
The third Painlevé equation in its generic form, often referred to as Painlevé-III($D_6$), is given by $$ \frac{{\rm d}^2u}{{\rm d}x^2} =\frac{1}{u}\left(\frac{{\rm d}u}{{\rm d}x}\right)^2-\frac{1}{x}\frac{{\rm d}u}{{\rm d}x}+\frac{αu^2+β}{x}+4u^3-\frac{4}{u}, \qquad α,β\in \mathbb C.$$ Starting from a generic initial solution $u_0(x)$ corresponding to parameters $α$, $β$, denoted as the triple…
▽ More
The third Painlevé equation in its generic form, often referred to as Painlevé-III($D_6$), is given by $$ \frac{{\rm d}^2u}{{\rm d}x^2} =\frac{1}{u}\left(\frac{{\rm d}u}{{\rm d}x}\right)^2-\frac{1}{x}\frac{{\rm d}u}{{\rm d}x}+\frac{αu^2+β}{x}+4u^3-\frac{4}{u}, \qquad α,β\in \mathbb C.$$ Starting from a generic initial solution $u_0(x)$ corresponding to parameters $α$, $β$, denoted as the triple $(u_0(x),α,β)$, we apply an explicit Bäcklund transformation to generate a family of solutions $(u_n(x),α+4n,β+4n)$ indexed by $n \in \mathbb N$. We study the large $n$ behavior of the solutions $(u_n(x),α+4n,β+4n)$ under the scaling $x=z/n$ in two different ways: (a) analyzing the convergence properties of series solutions to the equation, and (b) using a Riemann-Hilbert representation of the solution $u_n(z/n)$. Our main result is a proof that the limit of solutions $u_n(z/n)$ exists and is given by a solution of the degenerate Painlevé-III equation, known as Painlevé-III($D_8$), $$ \frac{{\rm d}^2U}{{\rm d}z^2} =\frac{1}{U}\left(\frac{{\rm d}U}{{\rm d}z}\right)^2-\frac{1}{z}\frac{{\rm d}U}{{\rm d}z}+\frac{4U^2+4}{z}.$$ A notable application of our result is to rational solutions of Painlevé-III($D_6$), which are constructed using the seed solution $(1,4m,-4m)$ where $m \in \mathbb C \setminus \big(\mathbb Z +\frac{1}{2}\big)$ and can be written as a particular ratio of Umemura polynomials. We identify the limiting solution in terms of both its initial condition at $z=0$ when it is well defined, and by its monodromy data in the general case. Furthermore, as a consequence of our analysis, we deduce the asymptotic behavior of generic solutions of Painlevé-III, both $D_6$ and $D_8$ at $z=0$. We also deduce the large $n$ behavior of the Umemura polynomials in a neighborhood of $z=0$.
△ Less
Submitted 9 March, 2024; v1 submitted 20 July, 2023;
originally announced July 2023.
-
An Open-Source Multi-functional Testing Platform for Optical Phase Change Materials
Authors:
Cosmin-Constantin Popescu,
Khoi Phuong Dao,
Luigi Ranno,
Brian Mills,
Louis Martin,
Yifei Zhang,
David Bono. Brian Neltner,
Tian Gu,
Juejun Hu,
Kiumars Aryana,
William M. Humphreys,
Hyun Jung Kim,
Steven Vitale,
Paul Miller,
Christopher Roberts,
Sarah Geiger,
Dennis Callahan,
Michael Moebius,
Myungkoo Kang,
Kathleen Richardson,
Carlos A. Ríos Ocampo
Abstract:
Owing to their unique tunable optical properties, chalcogenide phase change materials are increasingly being investigated for optics and photonics applications. However, in situ characterization of their phase transition characteristics is a capability that remains inaccessible to many researchers. In this article, we introduce a multi-functional silicon microheater platform capable of in situ mea…
▽ More
Owing to their unique tunable optical properties, chalcogenide phase change materials are increasingly being investigated for optics and photonics applications. However, in situ characterization of their phase transition characteristics is a capability that remains inaccessible to many researchers. In this article, we introduce a multi-functional silicon microheater platform capable of in situ measurement of structural, kinetic, optical, and thermal properties of these materials. The platform can be fabricated leveraging industry-standard silicon foundry manufacturing processes. We fully open-sourced this platform, including complete hardware design and associated software codes.
△ Less
Submitted 12 July, 2023;
originally announced July 2023.
-
A Virtual Diffractometer For Creating Synthetic HEDM Images of Tessellated and Meshed Finite Element Polycrystals
Authors:
Paul R. Dawson,
Matthew P. Miller
Abstract:
To assist in the planning of in situ loading, HEDM experiments by generating synthetic diffraction images of virtual samples in loaded and unloaded states. The user designates a target grain in the virtual sample and specifies the set of reflections for which images are to be generated. The code generates several intermediate images, including: (1) the points of intersection between the diffracted…
▽ More
To assist in the planning of in situ loading, HEDM experiments by generating synthetic diffraction images of virtual samples in loaded and unloaded states. The user designates a target grain in the virtual sample and specifies the set of reflections for which images are to be generated. The code generates several intermediate images, including: (1) the points of intersection between the diffracted beam direction vector the detector plane; (2) reflection-dependent frequency distributions of the diffraction angle, $ω$; and (3) plots of the diffraction volume-weighted intensity distributions for the specified set of reflections. The final sets of plots are facsimiles of pixelated detector images which take into account characteristics of specific detectors, including its pixel size and point spread behavior.
△ Less
Submitted 30 March, 2023;
originally announced March 2023.
-
Prognostic Covariate Adjustment for Binary Outcomes Using Stratification
Authors:
Alyssa M. Vanderbeek,
Jessica L. Ross,
David P. Miller,
Alejandro Schuler
Abstract:
Covariate adjustment and methods of incorporating historical data in randomized clinical trials (RCTs) each provide opportunities to increase trial power. We unite these approaches for the analysis of RCTs with binary outcomes based on the Cochran-Mantel-Haenszel (CMH) test for marginal risk ratio (RR). In PROCOVA-CMH, subjects are stratified on a single prognostic covariate reflective of their pr…
▽ More
Covariate adjustment and methods of incorporating historical data in randomized clinical trials (RCTs) each provide opportunities to increase trial power. We unite these approaches for the analysis of RCTs with binary outcomes based on the Cochran-Mantel-Haenszel (CMH) test for marginal risk ratio (RR). In PROCOVA-CMH, subjects are stratified on a single prognostic covariate reflective of their predicted outcome on the control treatment (e.g. placebo). This prognostic score is generated based on baseline covariates through a model trained on historical data. We propose two closed-form prospective estimators for the asymptotic sampling variance of the log RR that rely only on values obtainable from observed historical outcomes and the prognostic model. Importantly, these estimators can be used to inform sample size during trial planning. PROCOVA-CMH demonstrates type I error control and appropriate asymptotic coverage for valid inference. Like other covariate adjustment methods, PROCOVA-CMH can reduce the variance of the treatment effect estimate when compared to an unadjusted (unstratified) CMH analysis. In addition to statistical methods, simulations and a case study in Alzheimer's Disease are given to demonstrate performance. Results show that PROCOVA-CMH can provide a gain in power, which can be used to conduct smaller trials.
△ Less
Submitted 19 December, 2022;
originally announced December 2022.
-
Letter of Intent: Muonium R&D/Physics Program at the MTA
Authors:
C. Gatto,
C. Izzo,
C. J. Johnstone,
D. M. Kaplan,
K. R. Lynch,
D. C. Mancini,
A. Mazzacane,
B. McMorran,
J. P. Miller,
J. D. Phillips,
T. J. Phillips,
R. D. Reasenberg,
T. J. Roberts,
J. Terry
Abstract:
With the planned turn-on of the PIP-II 800 MeV superconducting proton linac, Fermilab will potentially become the world's best laboratory at which to carry out fundamental muon measurements, sensitive searches for symmetry violation, and precision tests of theory. In preparation, we propose to develop the techniques that will be needed. An R&D and physics program is proposed at the Fermilab MeV Te…
▽ More
With the planned turn-on of the PIP-II 800 MeV superconducting proton linac, Fermilab will potentially become the world's best laboratory at which to carry out fundamental muon measurements, sensitive searches for symmetry violation, and precision tests of theory. In preparation, we propose to develop the techniques that will be needed. An R&D and physics program is proposed at the Fermilab MeV Test Area to use the existing 400 MeV Linac to demonstrate the efficient production of a slow muonium beam using $μ^+$ stopped in a ~100-$μ$m-thick layer of superfluid helium, and to use that beam to measure muonium gravity.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Energetic electron precipitation driven by electromagnetic ion cyclotron waves from ELFIN's low altitude perspective
Authors:
V. Angelopoulos,
X. -J. Zhang,
A. V. Artemyev,
D. Mourenas,
E. Tsai,
C. Wilkins,
A. Runov,
J. Liu,
D. L. Turner,
W. Li,
K. Khurana,
R. E. Wirz,
V. A. Sergeev,
X. Meng,
J. Wu,
M. D. Hartinger,
T. Raita,
Y. Shen,
X. An,
X. Shi,
M. F. Bashir,
X. Shen,
L. Gan,
M. Qin,
L. Capannolo
, et al. (61 additional authors not shown)
Abstract:
We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data from the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibi…
▽ More
We review comprehensive observations of electromagnetic ion cyclotron (EMIC) wave-driven energetic electron precipitation using data from the energetic electron detector on the Electron Losses and Fields InvestigatioN (ELFIN) mission, two polar-orbiting low-altitude spinning CubeSats, measuring 50-5000 keV electrons with good pitch-angle and energy resolution. EMIC wave-driven precipitation exhibits a distinct signature in energy-spectrograms of the precipitating-to-trapped flux ratio: peaks at 0.5 MeV which are abrupt (bursty) with significant substructure (occasionally down to sub-second timescale). Multiple ELFIN passes over the same MLT sector allow us to study the spatial and temporal evolution of the EMIC wave - electron interaction region. Using two years of ELFIN data, we assemble a statistical database of 50 events of strong EMIC wave-driven precipitation. Most reside at L=5-7 at dusk, while a smaller subset exists at L=8-12 at post-midnight. The energies of the peak-precipitation ratio and of the half-peak precipitation ratio (our proxy for the minimum resonance energy) exhibit an L-shell dependence in good agreement with theoretical estimates based on prior statistical observations of EMIC wave power spectra. The precipitation ratio's spectral shape for the most intense events has an exponential falloff away from the peak (i.e., on either side of 1.45 MeV). It too agrees well with quasi-linear diffusion theory based on prior statistics of wave spectra. Sub-MeV electron precipitation observed concurrently with strong EMIC wave-driven 1MeV precipitation has a spectral shape that is consistent with efficient pitch-angle scattering down to 200-300 keV by much less intense higher frequency EMIC waves. These results confirm the critical role of EMIC waves in driving relativistic electron losses. Nonlinear effects may abound and require further investigation.
△ Less
Submitted 28 November, 2022;
originally announced November 2022.
-
Generation and motion of interfaces in a mass-conserving reaction-diffusion system
Authors:
Pearson W. Miller,
Daniel Fortunato,
Matteo Novaga,
Stanislav Y. Shvartsman,
Cyrill B. Muratov
Abstract:
Reaction-diffusion models with nonlocal constraints naturally arise as limiting cases of coupled bulk-surface models of intracellular signalling. In this paper, a minimal, mass-conserving model of cell-polarization on a curved membrane is analyzed in the limit of slow surface diffusion. Using the tools of formal asymptotics and calculus of variations, we study the characteristic wave-pinning behav…
▽ More
Reaction-diffusion models with nonlocal constraints naturally arise as limiting cases of coupled bulk-surface models of intracellular signalling. In this paper, a minimal, mass-conserving model of cell-polarization on a curved membrane is analyzed in the limit of slow surface diffusion. Using the tools of formal asymptotics and calculus of variations, we study the characteristic wave-pinning behavior of this system on three dynamical timescales. On the short timescale, generation of an interface separating high- and low-concentration domains is established under suitable conditions. Intermediate timescale dynamics is shown to lead to a uniform growth or shrinking of these domains to sizes which are fixed by global parameters. Finally, the long time dynamics reduces to area-preserving geodesic curvature flow that may lead to multi-interface steady state solutions. These results provide a foundation for studying cell polarization and related phenomena in biologically relevant geometries.
△ Less
Submitted 2 October, 2022;
originally announced October 2022.
-
Restricted mean survival time estimate using covariate adjusted pseudovalue regression to improve precision
Authors:
Yunfan Li,
Jessica L. Ross,
Aaron M. Smith,
David P. Miller
Abstract:
Covariate adjustment is desired by both practitioners and regulators of randomized clinical trials because it improves precision for estimating treatment effects. However, covariate adjustment presents a particular challenge in time-to-event analysis. We propose to apply covariate adjusted pseudovalue regression to estimate between-treatment difference in restricted mean survival times (RMST). Our…
▽ More
Covariate adjustment is desired by both practitioners and regulators of randomized clinical trials because it improves precision for estimating treatment effects. However, covariate adjustment presents a particular challenge in time-to-event analysis. We propose to apply covariate adjusted pseudovalue regression to estimate between-treatment difference in restricted mean survival times (RMST). Our proposed method incorporates a prognostic covariate to increase precision of treatment effect estimate, maintaining strict type I error control without introducing bias. In addition, the amount of increase in precision can be quantified and taken into account in sample size calculation at the study design stage. Consequently, our proposed method provides the ability to design smaller randomized studies at no expense to statistical power.
△ Less
Submitted 18 July, 2023; v1 submitted 8 August, 2022;
originally announced August 2022.
-
The storage ring proton EDM experiment
Authors:
Jim Alexander,
Vassilis Anastassopoulos,
Rick Baartman,
Stefan Baeßler,
Franco Bedeschi,
Martin Berz,
Michael Blaskiewicz,
Themis Bowcock,
Kevin Brown,
Dmitry Budker,
Sergey Burdin,
Brendan C. Casey,
Gianluigi Casse,
Giovanni Cantatore,
Timothy Chupp,
Hooman Davoudiasl,
Dmitri Denisov,
Milind V. Diwan,
George Fanourakis,
Antonios Gardikiotis,
Claudio Gatti,
James Gooding,
Renee Fatemi,
Wolfram Fischer,
Peter Graham
, et al. (52 additional authors not shown)
Abstract:
We describe a proposal to search for an intrinsic electric dipole moment (EDM) of the proton with a sensitivity of \targetsens, based on the vertical rotation of the polarization of a stored proton beam. The New Physics reach is of order $10^~3$TeV mass scale. Observation of the proton EDM provides the best probe of CP-violation in the Higgs sector, at a level of sensitivity that may be inaccessib…
▽ More
We describe a proposal to search for an intrinsic electric dipole moment (EDM) of the proton with a sensitivity of \targetsens, based on the vertical rotation of the polarization of a stored proton beam. The New Physics reach is of order $10^~3$TeV mass scale. Observation of the proton EDM provides the best probe of CP-violation in the Higgs sector, at a level of sensitivity that may be inaccessible to electron-EDM experiments. The improvement in the sensitivity to $θ_{QCD}$, a parameter crucial in axion and axion dark matter physics, is about three orders of magnitude.
△ Less
Submitted 25 April, 2022;
originally announced May 2022.
-
Weak topologies for unbounded nets in CAT($0$) spaces
Authors:
Philip Miller,
Arian Berdellima,
Max Wardetzky
Abstract:
Weak topologies that yield weak convergence for bounded sequences and nets in CAT($0$) spaces have been studied in the past. We are here concerned with weak topologies that yield weak convergence of unbounded sequences and nets. We analyze two such topologies that generalize the weak topology on Hilbert spaces and that agree with the strong topology on a CAT($0$) space if and only if the space is…
▽ More
Weak topologies that yield weak convergence for bounded sequences and nets in CAT($0$) spaces have been studied in the past. We are here concerned with weak topologies that yield weak convergence of unbounded sequences and nets. We analyze two such topologies that generalize the weak topology on Hilbert spaces and that agree with the strong topology on a CAT($0$) space if and only if the space is locally compact.
△ Less
Submitted 4 April, 2022;
originally announced April 2022.
-
Electric dipole moments and the search for new physics
Authors:
Ricardo Alarcon,
Jim Alexander,
Vassilis Anastassopoulos,
Takatoshi Aoki,
Rick Baartman,
Stefan Baeßler,
Larry Bartoszek,
Douglas H. Beck,
Franco Bedeschi,
Robert Berger,
Martin Berz,
Hendrick L. Bethlem,
Tanmoy Bhattacharya,
Michael Blaskiewicz,
Thomas Blum,
Themis Bowcock,
Anastasia Borschevsky,
Kevin Brown,
Dmitry Budker,
Sergey Burdin,
Brendan C. Casey,
Gianluigi Casse,
Giovanni Cantatore,
Lan Cheng,
Timothy Chupp
, et al. (118 additional authors not shown)
Abstract:
Static electric dipole moments of nondegenerate systems probe mass scales for physics beyond the Standard Model well beyond those reached directly at high energy colliders. Discrimination between different physics models, however, requires complementary searches in atomic-molecular-and-optical, nuclear and particle physics. In this report, we discuss the current status and prospects in the near fu…
▽ More
Static electric dipole moments of nondegenerate systems probe mass scales for physics beyond the Standard Model well beyond those reached directly at high energy colliders. Discrimination between different physics models, however, requires complementary searches in atomic-molecular-and-optical, nuclear and particle physics. In this report, we discuss the current status and prospects in the near future for a compelling suite of such experiments, along with developments needed in the encompassing theoretical framework.
△ Less
Submitted 4 April, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
Mu2e-II: Muon to electron conversion with PIP-II
Authors:
K. Byrum,
S. Corrodi,
Y. Oksuzian,
P. Winter,
L. Xia,
A. W. J. Edmonds,
J. P. Miller,
J. Mott,
W. J. Marciano,
R. Szafron,
R. Bonventre,
D. N. Brown,
Yu. G. Kolomensky,
O. Ning,
V. Singh,
E. Prebys,
L. Borrel,
B. Echenard,
D. G. Hitlin,
C. Hu,
D. X. Lin,
S. Middleton,
F. C. Porter,
L. Zhang,
R. -Y. Zhu
, et al. (83 additional authors not shown)
Abstract:
An observation of Charged Lepton Flavor Violation (CLFV) would be unambiguous evidence for physics beyond the Standard Model. The Mu2e and COMET experiments, under construction, are designed to push the sensitivity to CLFV in the mu to e conversion process to unprecedented levels. Whether conversion is observed or not, there is a strong case to be made for further improving sensitivity, or for exa…
▽ More
An observation of Charged Lepton Flavor Violation (CLFV) would be unambiguous evidence for physics beyond the Standard Model. The Mu2e and COMET experiments, under construction, are designed to push the sensitivity to CLFV in the mu to e conversion process to unprecedented levels. Whether conversion is observed or not, there is a strong case to be made for further improving sensitivity, or for examining the process on additional target materials. Mu2e-II is a proposed upgrade to Mu2e, with at least an additional order of magnitude in sensitivity to the conversion rate over Mu2e. The approach and challenges for this proposal are summarized. Mu2e-II may be regarded as the next logical step in a continued high-intensity muon program at FNAL.
△ Less
Submitted 16 March, 2022; v1 submitted 14 March, 2022;
originally announced March 2022.
-
Convergence Rates for Oversmoothing Banach Space Regularization
Authors:
Philip Miller,
Thorsten Hohage
Abstract:
This paper studies Tikhonov regularization for finitely smoothing operators in Banach spaces when the penalization enforces too much smoothness in the sense that the penalty term is not finite at the true solution. In a Hilbert space setting, Natterer (1984) showed with the help of spectral theory that optimal rates can be achieved in this situation. ('Oversmoothing does not harm.') For oversmooth…
▽ More
This paper studies Tikhonov regularization for finitely smoothing operators in Banach spaces when the penalization enforces too much smoothness in the sense that the penalty term is not finite at the true solution. In a Hilbert space setting, Natterer (1984) showed with the help of spectral theory that optimal rates can be achieved in this situation. ('Oversmoothing does not harm.') For oversmoothing variational regularization in Banach spaces only very recently progress has been achieved in several papers on different settings, all of which construct families of smooth approximations to the true solution. In this paper we propose to construct such a family of smooth approximations based on $K$-interpolation theory. We demonstrate that this leads to simple, self-contained proofs and to rather general results. In particular, we obtain optimal convergence rates for bounded variation regularization, general Besov penalty terms and $\ell^p$ wavelet penalization with $p<1$ which cannot be treated by previous approaches. We also derive minimax optimal rates for white noise models. Our theoretical results are confirmed in numerical experiments.
△ Less
Submitted 4 March, 2022;
originally announced March 2022.
-
On the algebraic solutions of the Painleve-III (D7) equation
Authors:
Robert J. Buckingham,
Peter D. Miller
Abstract:
The D7 degeneration of the Painleve-III equation has solutions that are rational functions of $x^{1/3}$ for certain parameter values. We apply the isomonodromy method to obtain a Riemann-Hilbert representation of these solutions. We demonstrate the utility of this representation by analyzing rigorously the behavior of the solutions in the large parameter limit.
The D7 degeneration of the Painleve-III equation has solutions that are rational functions of $x^{1/3}$ for certain parameter values. We apply the isomonodromy method to obtain a Riemann-Hilbert representation of these solutions. We demonstrate the utility of this representation by analyzing rigorously the behavior of the solutions in the large parameter limit.
△ Less
Submitted 8 February, 2022;
originally announced February 2022.
-
Broader Universality of Rogue Waves of Infinite Order
Authors:
Deniz Bilman,
Peter D. Miller
Abstract:
We show that the same special solution of the focusing nonlinear Schrödinger equation that has been shown to arise in a certain near-field/large-order limit from soliton and Peregrine-like rogue wave solutions actually arises universally from an arbitrary background solution when subjected to a sequence of iterated Bäcklund transformations.
We show that the same special solution of the focusing nonlinear Schrödinger equation that has been shown to arise in a certain near-field/large-order limit from soliton and Peregrine-like rogue wave solutions actually arises universally from an arbitrary background solution when subjected to a sequence of iterated Bäcklund transformations.
△ Less
Submitted 30 December, 2021;
originally announced December 2021.
-
Heroin addiction hijacks the Nucleus Accumbens: craving and reactivity to naturalistic stimuli
Authors:
Greg Kronberg,
Ahmet O. Ceceli,
Yuefeng Huang,
Pierre-Olivier Gaudreault,
Sarah King,
Natalie McClain,
Pazia Miller,
Lily Gabay,
Devarshi Vasa,
Pias Malaker,
Defne Ekin,
Nelly Alia-Klein,
Rita Z. Goldstein
Abstract:
Drug-related cues hijack attention away from alternative reinforcers in drug addiction, inducing craving and motivating drug-seeking. However, the neural correlates underlying this biased processing, its expression in the real-world, and its relationship to cue-induced craving are not fully established, especially in opioid addiction. Here we tracked inter-brain synchronization in the Nucleus Accu…
▽ More
Drug-related cues hijack attention away from alternative reinforcers in drug addiction, inducing craving and motivating drug-seeking. However, the neural correlates underlying this biased processing, its expression in the real-world, and its relationship to cue-induced craving are not fully established, especially in opioid addiction. Here we tracked inter-brain synchronization in the Nucleus Accumbens (NAc), a hub of motivational salience, while heroin-addicted individuals and healthy control subjects watched the same engaging heroin-related movie. Strikingly, the left NAc was synchronized during drug scenes in the addicted individuals and non-drug scenes in controls, predicting scene- and movie-induced heroin craving in the former. Our results open a window into the neurobiology underlying shared drug-biased processing of naturalistic stimuli and cue-induced craving in opiate addiction as they unfold in the real world.
△ Less
Submitted 15 December, 2021;
originally announced December 2021.
-
Evaluation of Static Vulnerability Detection Tools with Java Cryptographic API Benchmarks
Authors:
Sharmin Afrose,
Ya Xiao,
Sazzadur Rahaman,
Barton P. Miller,
Danfeng,
Yao
Abstract:
Several studies showed that misuses of cryptographic APIs are common in real-world code (e.g., Apache projects and Android apps). There exist several open-sourced and commercial security tools that automatically screen Java programs to detect misuses. To compare their accuracy and security guarantees, we develop two comprehensive benchmarks named CryptoAPI-Bench and ApacheCryptoAPI-Bench. CryptoAP…
▽ More
Several studies showed that misuses of cryptographic APIs are common in real-world code (e.g., Apache projects and Android apps). There exist several open-sourced and commercial security tools that automatically screen Java programs to detect misuses. To compare their accuracy and security guarantees, we develop two comprehensive benchmarks named CryptoAPI-Bench and ApacheCryptoAPI-Bench. CryptoAPI-Bench consists of 181 unit test cases that cover basic cases, as well as complex cases, including interprocedural, field sensitive, multiple class test cases, and path sensitive data flow of misuse cases. The benchmark also includes correct cases for testing false-positive rates. The ApacheCryptoAPI-Bench consists of 121 cryptographic cases from 10 Apache projects. We evaluate four tools, namely, SpotBugs, CryptoGuard, CrySL, and Coverity using both benchmarks. We present their performance and comparative analysis. The ApacheCryptoAPI-Bench also examines the scalability of the tools. Our benchmarks are useful for advancing state-of-the-art solutions in the space of misuse detection.
△ Less
Submitted 7 December, 2021;
originally announced December 2021.
-
A Measurement of Proton, Deuteron, Triton and Alpha Particle Emission after Nuclear Muon Capture on Al, Si and Ti with the AlCap Experiment
Authors:
AlCap Collaboration,
Andrew Edmonds,
John Quirk,
Ming-Liang Wong,
Damien Alexander,
Robert H. Bernstein,
Aji Daniel,
Eleonora Diociaiuti,
Raffaella Donghia,
Ewen L. Gillies,
Ed V. Hungerford,
Peter Kammel,
Benjamin E. Krikler,
Yoshitaka Kuno,
Mark Lancaster,
R. Phillip Litchfield,
James P. Miller,
Anthony Palladino,
Jose Repond,
Akira Sato,
Ivano Sarra,
Stefano Roberto Soleti,
Vladimir Tishchenko,
Nam H. Tran,
Yoshi Uchida
, et al. (2 additional authors not shown)
Abstract:
Heavy charged particles after nuclear muon capture are an important nuclear physics background to the muon-to-electron conversion experiments Mu2e and COMET, which will search for charged lepton flavor violation at an unprecedented level of sensitivity. The AlCap experiment measured the yield and energy spectra of protons, deuterons, tritons, and alpha particles emitted after the nuclear capture o…
▽ More
Heavy charged particles after nuclear muon capture are an important nuclear physics background to the muon-to-electron conversion experiments Mu2e and COMET, which will search for charged lepton flavor violation at an unprecedented level of sensitivity. The AlCap experiment measured the yield and energy spectra of protons, deuterons, tritons, and alpha particles emitted after the nuclear capture of muons stopped in Al, Si, and Ti in the low energy range relevant for the muon-to-electron conversion experiments. Individual charged particle types were identified in layered silicon detector packages and their initial energy distributions were unfolded from the observed energy spectra. Detailed information on yields and energy spectra for all observed nuclei are presented in the paper.
△ Less
Submitted 1 April, 2022; v1 submitted 19 October, 2021;
originally announced October 2021.
-
Using Discrete Harmonic Expansions and Equilibrium Conditions to Estimate Intragrain Stress Distributions in Polycrystals from Grain-Averaged Data
Authors:
Paul R. Dawson,
Matthew. P. Miller
Abstract:
The application of harmonic expansions to estimate intra-grain stress distributions from grain-averaged stress data is presented that extends the capabilities of the open source code, MechMonics. The method is based on using an optimization algorithm to determine the harmonic expansion weights that reduce the violation of equilibrium while maintaining prescribed grain-averages. The method is demon…
▽ More
The application of harmonic expansions to estimate intra-grain stress distributions from grain-averaged stress data is presented that extends the capabilities of the open source code, MechMonics. The method is based on using an optimization algorithm to determine the harmonic expansion weights that reduce the violation of equilibrium while maintaining prescribed grain-averages. The method is demonstrated using synthetic data generated for uniaxial extension of a virtual polycrystal with the mechMet code.
△ Less
Submitted 29 June, 2021;
originally announced June 2021.
-
On the Maxwell-Bloch System in the Sharp-Line Limit Without Solitons
Authors:
Sitai Li,
Peter D. Miller
Abstract:
We study the (characteristic) Cauchy problem for the Maxwell-Bloch equations of light-matter interaction via asymptotics, under assumptions that prevent the generation of solitons. Our analysis clarifies some features of the sense in which physically-motivated initial/boundary conditions are satisfied. In particular, we present a proper Riemann-Hilbert problem that generates the unique causal solu…
▽ More
We study the (characteristic) Cauchy problem for the Maxwell-Bloch equations of light-matter interaction via asymptotics, under assumptions that prevent the generation of solitons. Our analysis clarifies some features of the sense in which physically-motivated initial/boundary conditions are satisfied. In particular, we present a proper Riemann-Hilbert problem that generates the unique causal solution to the Cauchy problem, that is, the solution vanishes outside of the light cone. Inside the light cone, we relate the leading-order asymptotics to self-similar solutions that satisfy a system of ordinary differential equations related to the Painlevé-III (PIII) equation. We identify these solutions and show that they are related to a family of PIII solutions recently discovered in connection with several limiting processes involving the focusing nonlinear Schrödinger equation. We fully explain a resulting boundary layer phenomenon in which, even for smooth initial data (an incident pulse), the solution makes a sudden transition over an infinitesimally small propagation distance. At a formal level, this phenomenon has been described by other authors in terms of the PIII self-similar solutions. We make this observation precise and for the first time we relate the PIII self-similar solutions to the Cauchy problem. Our analysis of the asymptotic behavior satisfied by the optical field and medium density matrix reveals slow decay of the optical field in one direction that is actually inconsistent with the simplest version of scattering theory. Our results identify a precise generic condition on an optical pulse incident on an initially-unstable medium sufficient for the pulse to stimulate the decay of the medium to its stable state.
△ Less
Submitted 22 May, 2021;
originally announced May 2021.
-
Ultra-compact nonvolatile phase shifter based on electrically reprogrammable transparent phase change materials
Authors:
Carlos Ríos,
Qingyang Du,
Yifei Zhang,
Cosmin-Constantin Popescu,
Mikhail Y. Shalaginov,
Paul Miller,
Christopher Roberts,
Myungkoo Kang,
Kathleen A. Richardson,
Tian Gu,
Steven A. Vitale,
Juejun Hu
Abstract:
Energy-efficient programmable photonic integrated circuits (PICs) are the cornerstone of on-chip classical and quantum optical technologies. Optical phase shifters constitute the fundamental building blocks which enable these programmable PICs. Thus far, carrier modulation and thermo-optical effect are the chosen phenomena for ultrafast and low-loss phase shifters, respectively; however, the state…
▽ More
Energy-efficient programmable photonic integrated circuits (PICs) are the cornerstone of on-chip classical and quantum optical technologies. Optical phase shifters constitute the fundamental building blocks which enable these programmable PICs. Thus far, carrier modulation and thermo-optical effect are the chosen phenomena for ultrafast and low-loss phase shifters, respectively; however, the state and information they carry are lost once the power is turned off-they are volatile. The volatility not only compromises energy efficiency due to their demand for constant power supply, but also precludes them from emerging applications such as in-memory computing. To circumvent this limitation, we introduce a novel phase shifting mechanism that exploits the nonvolatile refractive index modulation upon structural phase transition of Sb$_{2}$Se$_{3}$, a bi-stable transparent phase change material. A zero-static power and electrically-driven phase shifter was realized on a foundry-processed silicon-on-insulator platform, featuring record phase modulation up to 0.09 $π$/$μ$m and a low insertion loss of 0.3 dB/$π$, which can be further improved upon streamlined design. We also pioneered a one-step partial amorphization scheme to enhance the speed and energy efficiency of PCM devices. A diverse cohort of programmable photonic devices were demonstrated based on the ultra-compact PCM phase shifter.
△ Less
Submitted 21 March, 2022; v1 submitted 12 May, 2021;
originally announced May 2021.
-
Analysis of a Three-Dimensional Slip Field in a Hexagonal Ti Alloy from in-situ High-Energy X-ray Diffraction Microscopy Data
Authors:
Darren C. Pagan,
Kelly E. Nygren,
Matthew P. Miller
Abstract:
Here we analyze a three-dimensional distribution of crystallographic slip measured in-situ during the uniaxial deformation of hexagonal Ti-7Al. The slip field is reconstructed using a novel methodology that combines spatially resolved lattice orientation fields and grain-averaged stresses measured using high-energy X-ray diffraction microscopy (HEDM) with crystal plasticity. Analysis is performed…
▽ More
Here we analyze a three-dimensional distribution of crystallographic slip measured in-situ during the uniaxial deformation of hexagonal Ti-7Al. The slip field is reconstructed using a novel methodology that combines spatially resolved lattice orientation fields and grain-averaged stresses measured using high-energy X-ray diffraction microscopy (HEDM) with crystal plasticity. Analysis is performed to explore lattice orientation dependence, stress dependence, and connectivity (network relationships) of grains experiencing elevated amounts of slip. Elevated slip is found to be primarily associated with a single large network of connected grains, and within this network, a clustered group of grains oriented favorably for slip are found to have outsized structural importance. The effect of different rate sensitivities of families of slip systems on reconstructed slip activity is also discussed.
△ Less
Submitted 25 April, 2021;
originally announced April 2021.
-
Measurement of the Positive Muon Anomalous Magnetic Moment to 0.46 ppm
Authors:
B. Abi,
T. Albahri,
S. Al-Kilani,
D. Allspach,
L. P. Alonzi,
A. Anastasi,
A. Anisenkov,
F. Azfar,
K. Badgley,
S. Baeßler,
I. Bailey,
V. A. Baranov,
E. Barlas-Yucel,
T. Barrett,
E. Barzi,
A. Basti,
F. Bedeschi,
A. Behnke,
M. Berz,
M. Bhattacharya,
H. P. Binney,
R. Bjorkquist,
P. Bloom,
J. Bono,
E. Bottalico
, et al. (212 additional authors not shown)
Abstract:
We present the first results of the Fermilab Muon g-2 Experiment for the positive muon magnetic anomaly $a_μ\equiv (g_μ-2)/2$. The anomaly is determined from the precision measurements of two angular frequencies. Intensity variation of high-energy positrons from muon decays directly encodes the difference frequency $ω_a$ between the spin-precession and cyclotron frequencies for polarized muons in…
▽ More
We present the first results of the Fermilab Muon g-2 Experiment for the positive muon magnetic anomaly $a_μ\equiv (g_μ-2)/2$. The anomaly is determined from the precision measurements of two angular frequencies. Intensity variation of high-energy positrons from muon decays directly encodes the difference frequency $ω_a$ between the spin-precession and cyclotron frequencies for polarized muons in a magnetic storage ring. The storage ring magnetic field is measured using nuclear magnetic resonance probes calibrated in terms of the equivalent proton spin precession frequency ${\tildeω'^{}_p}$ in a spherical water sample at 34.7$^{\circ}$C. The ratio $ω_a / {\tildeω'^{}_p}$, together with known fundamental constants, determines $a_μ({\rm FNAL}) = 116\,592\,040(54)\times 10^{-11}$ (0.46\,ppm). The result is 3.3 standard deviations greater than the standard model prediction and is in excellent agreement with the previous Brookhaven National Laboratory (BNL) E821 measurement. After combination with previous measurements of both $μ^+$ and $μ^-$, the new experimental average of $a_μ({\rm Exp}) = 116\,592\,061(41)\times 10^{-11}$ (0.35\,ppm) increases the tension between experiment and theory to 4.2 standard deviations
△ Less
Submitted 7 April, 2021;
originally announced April 2021.
-
Measurement of the anomalous precession frequency of the muon in the Fermilab Muon g-2 experiment
Authors:
T. Albahri,
A. Anastasi,
A. Anisenkov,
K. Badgley,
S. Baeßler,
I. Bailey,
V. A. Baranov,
E. Barlas-Yucel,
T. Barrett,
A. Basti,
F. Bedeschi,
M. Berz,
M. Bhattacharya,
H. P. Binney,
P. Bloom,
J. Bono,
E. Bottalico,
T. Bowcock,
G. Cantatore,
R. M. Carey,
B. C. K. Casey,
D. Cauz,
R. Chakraborty,
S. P. Chang,
A. Chapelain
, et al. (153 additional authors not shown)
Abstract:
The Muon g-2 Experiment at Fermi National Accelerator Laboratory (FNAL) has measured the muon anomalous precession frequency $ω_a$ to an uncertainty of 434 parts per billion (ppb), statistical, and 56 ppb, systematic, with data collected in four storage ring configurations during its first physics run in 2018. When combined with a precision measurement of the magnetic field of the experiment's muo…
▽ More
The Muon g-2 Experiment at Fermi National Accelerator Laboratory (FNAL) has measured the muon anomalous precession frequency $ω_a$ to an uncertainty of 434 parts per billion (ppb), statistical, and 56 ppb, systematic, with data collected in four storage ring configurations during its first physics run in 2018. When combined with a precision measurement of the magnetic field of the experiment's muon storage ring, the precession frequency measurement determines a muon magnetic anomaly of $a_μ({\rm FNAL}) = 116\,592\,040(54) \times 10^{-11}$ (0.46 ppm). This article describes the multiple techniques employed in the reconstruction, analysis and fitting of the data to measure the precession frequency. It also presents the averaging of the results from the eleven separate determinations of ω_a, and the systematic uncertainties on the result.
△ Less
Submitted 7 April, 2021;
originally announced April 2021.
-
Beam dynamics corrections to the Run-1 measurement of the muon anomalous magnetic moment at Fermilab
Authors:
T. Albahri,
A. Anastasi,
K. Badgley,
S. Baeßler,
I. Bailey,
V. A. Baranov,
E. Barlas-Yucel,
T. Barrett,
F. Bedeschi,
M. Berz,
M. Bhattacharya,
H. P. Binney,
P. Bloom,
J. Bono,
E. Bottalico,
T. Bowcock,
G. Cantatore,
R. M. Carey,
B. C. K. Casey,
D. Cauz,
R. Chakraborty,
S. P. Chang,
A. Chapelain,
S. Charity,
R. Chislett
, et al. (152 additional authors not shown)
Abstract:
This paper presents the beam dynamics systematic corrections and their uncertainties for the Run-1 data set of the Fermilab Muon g-2 Experiment. Two corrections to the measured muon precession frequency $ω_a^m$ are associated with well-known effects owing to the use of electrostatic quadrupole (ESQ) vertical focusing in the storage ring. An average vertically oriented motional magnetic field is fe…
▽ More
This paper presents the beam dynamics systematic corrections and their uncertainties for the Run-1 data set of the Fermilab Muon g-2 Experiment. Two corrections to the measured muon precession frequency $ω_a^m$ are associated with well-known effects owing to the use of electrostatic quadrupole (ESQ) vertical focusing in the storage ring. An average vertically oriented motional magnetic field is felt by relativistic muons passing transversely through the radial electric field components created by the ESQ system. The correction depends on the stored momentum distribution and the tunes of the ring, which has relatively weak vertical focusing. Vertical betatron motions imply that the muons do not orbit the ring in a plane exactly orthogonal to the vertical magnetic field direction. A correction is necessary to account for an average pitch angle associated with their trajectories. A third small correction is necessary because muons that escape the ring during the storage time are slightly biased in initial spin phase compared to the parent distribution. Finally, because two high-voltage resistors in the ESQ network had longer than designed RC time constants, the vertical and horizontal centroids and envelopes of the stored muon beam drifted slightly, but coherently, during each storage ring fill. This led to the discovery of an important phase-acceptance relationship that requires a correction. The sum of the corrections to $ω_a^m$ is 0.50 $\pm$ 0.09 ppm; the uncertainty is small compared to the 0.43 ppm statistical precision of $ω_a^m$.
△ Less
Submitted 23 April, 2021; v1 submitted 7 April, 2021;
originally announced April 2021.
-
Magnetic Field Measurement and Analysis for the Muon g-2 Experiment at Fermilab
Authors:
T. Albahri,
A. Anastasi,
K. Badgley,
S. Baeßler,
I. Bailey,
V. A. Baranov,
E. Barlas-Yucel,
T. Barrett,
F. Bedeschi,
M. Berz,
M. Bhattacharya,
H. P. Binney,
P. Bloom,
J. Bono,
E. Bottalico,
T. Bowcock,
G. Cantatore,
R. M. Carey,
B. C. K. Casey,
D. Cauz,
R. Chakraborty,
S. P. Chang,
A. Chapelain,
S. Charity,
R. Chislett
, et al. (148 additional authors not shown)
Abstract:
The Fermi National Accelerator Laboratory has measured the anomalous precession frequency $a^{}_μ= (g^{}_μ-2)/2$ of the muon to a combined precision of 0.46 parts per million with data collected during its first physics run in 2018. This paper documents the measurement of the magnetic field in the muon storage ring. The magnetic field is monitored by nuclear magnetic resonance systems and calibrat…
▽ More
The Fermi National Accelerator Laboratory has measured the anomalous precession frequency $a^{}_μ= (g^{}_μ-2)/2$ of the muon to a combined precision of 0.46 parts per million with data collected during its first physics run in 2018. This paper documents the measurement of the magnetic field in the muon storage ring. The magnetic field is monitored by nuclear magnetic resonance systems and calibrated in terms of the equivalent proton spin precession frequency in a spherical water sample at 34.7$^\circ$C. The measured field is weighted by the muon distribution resulting in $\tildeω'^{}_p$, the denominator in the ratio $ω^{}_a$/$\tildeω'^{}_p$ that together with known fundamental constants yields $a^{}_μ$. The reported uncertainty on $\tildeω'^{}_p$ for the Run-1 data set is 114 ppb consisting of uncertainty contributions from frequency extraction, calibration, mapping, tracking, and averaging of 56 ppb, and contributions from fast transient fields of 99 ppb.
△ Less
Submitted 17 June, 2022; v1 submitted 7 April, 2021;
originally announced April 2021.
-
Mechanical Analyses using Discrete Harmonic Expansions of Grain-Based Field Data
Authors:
Paul R. Dawson,
Matthew P. Miller
Abstract:
A methodology for computing expansion basis functions using discrete harmonic modes is presented. The discrete harmonic modes are determined grain-by-grain for virtual polycrystals for which finite element meshes are available. The expansion weights associated with representing field variables over grain domains are determined by exploiting the orthogonality of the harmonic modes. The methodology…
▽ More
A methodology for computing expansion basis functions using discrete harmonic modes is presented. The discrete harmonic modes are determined grain-by-grain for virtual polycrystals for which finite element meshes are available. The expansion weights associated with representing field variables over grain domains are determined by exploiting the orthogonality of the harmonic modes. The methodology is demonstrated with the representation of the axial stress distributions during tensile loading of a polycrystalline sample. An open source code, MechMonics, is available to researchers wishing to use the methodology to analyze data.
△ Less
Submitted 9 March, 2021;
originally announced March 2021.
-
Extreme Superposition: High-Order Fundamental Rogue Waves in the Far-Field Regime
Authors:
Deniz Bilman,
Peter D. Miller
Abstract:
We study fundamental rogue-wave solutions of the focusing nonlinear Schrödinger equation in the limit that the order of the rogue wave is large and the independent variables $(x,t)$ are proportional to the order (the far-field limit). We first formulate a Riemann-Hilbert representation of these solutions that allows the order to vary continuously rather than by integer increments. The intermediate…
▽ More
We study fundamental rogue-wave solutions of the focusing nonlinear Schrödinger equation in the limit that the order of the rogue wave is large and the independent variables $(x,t)$ are proportional to the order (the far-field limit). We first formulate a Riemann-Hilbert representation of these solutions that allows the order to vary continuously rather than by integer increments. The intermediate solutions in this continuous family include also soliton solutions for zero boundary conditions spectrally encoded by a single complex-conjugate pair of poles of arbitrary order, as well as other solutions having nonzero boundary conditions matching those of the rogue waves albeit with far slower decay as $x\to\pm\infty$. The large-order far-field asymptotic behavior of the solution depends on which of three disjoint regions $\mathcal{C}$, $\mathcal{S}$, and $\mathcal{E}$ contains the rescaled variables. On the regions $\mathcal{C}$ and $\mathcal{S}$ we show that the asymptotic behavior is the same for all continuous orders, while in the region $\mathcal{E}$ the discrete sequence of rogue-wave orders produces distinctive asymptotic behavior that is different from other cases.
△ Less
Submitted 27 February, 2021;
originally announced March 2021.
-
Variational Regularization Theory Based on Image Space Approximation Rates
Authors:
Philip Miller
Abstract:
We present a new approach to convergence rate results for variational regularization. Avoiding Bregman distances and using image space approximation rates as source conditions we prove a nearly minimax theorem showing that the modulus of continuity is an upper bound on the reconstruction error up to a constant. Applied to Besov space regularization we obtain convergence rate results for $0,2,q$- a…
▽ More
We present a new approach to convergence rate results for variational regularization. Avoiding Bregman distances and using image space approximation rates as source conditions we prove a nearly minimax theorem showing that the modulus of continuity is an upper bound on the reconstruction error up to a constant. Applied to Besov space regularization we obtain convergence rate results for $0,2,q$- and $0,p,p$-penalties without restrictions on $p,q\in (1,\infty).$ Finally we prove equivalence of Hölder-type variational source conditions, bounds on the defect of the Tikhonov functional, and image space approximation rates.
△ Less
Submitted 29 March, 2021; v1 submitted 1 September, 2020;
originally announced September 2020.
-
The Relevance of Classic Fuzz Testing: Have We Solved This One?
Authors:
Barton P. Miller,
Mengxiao Zhang,
Elisa R. Heymann
Abstract:
As fuzz testing has passed its 30th anniversary, and in the face of the incredible progress in fuzz testing techniques and tools, the question arises if the classic, basic fuzz technique is still useful and applicable? In that tradition, we have updated the basic fuzz tools and testing scripts and applied them to a large collection of Unix utilities on Linux, FreeBSD, and MacOS. As before, our fai…
▽ More
As fuzz testing has passed its 30th anniversary, and in the face of the incredible progress in fuzz testing techniques and tools, the question arises if the classic, basic fuzz technique is still useful and applicable? In that tradition, we have updated the basic fuzz tools and testing scripts and applied them to a large collection of Unix utilities on Linux, FreeBSD, and MacOS. As before, our failure criteria was whether the program crashed or hung. We found that 9 crash or hang out of 74 utilities on Linux, 15 out of 78 utilities on FreeBSD, and 12 out of 76 utilities on MacOS. A total of 24 different utilities failed across the three platforms. We note that these failure rates are somewhat higher than our in previous 1995, 2000, and 2006 studies of the reliability of command line utilities. In the basic fuzz tradition, we debugged each failed utility and categorized the causes the failures. Classic categories of failures, such as pointer and array errors and not checking return codes, were still broadly present in the current results. In addition, we found a couple of new categories of failures appearing. We present examples of these failures to illustrate the programming practices that allowed them to happen. As a side note, we tested the limited number of utilities available in a modern programming language (Rust) and found them to be of no better reliability than the standard ones.
△ Less
Submitted 4 January, 2021; v1 submitted 14 August, 2020;
originally announced August 2020.
-
Large-degree asymptotics of rational Painleve-IV solutions by the isomonodromy method
Authors:
Robert J. Buckingham,
Peter D. Miller
Abstract:
The Painleve-IV equation has two families of rational solutions generated respectively by the generalized Hermite polynomials and the generalized Okamoto polynomials. We apply the isomonodromy method to represent all of these rational solutions by means of two related Riemann-Hilbert problems, each of which involves two integer-valued parameters related to the two parameters in the Painleve-IV equ…
▽ More
The Painleve-IV equation has two families of rational solutions generated respectively by the generalized Hermite polynomials and the generalized Okamoto polynomials. We apply the isomonodromy method to represent all of these rational solutions by means of two related Riemann-Hilbert problems, each of which involves two integer-valued parameters related to the two parameters in the Painleve-IV equation. We then use the steepest-descent method to analyze the rational solutions in the limit that at least one of the parameters is large. Our analysis provides rigorous justification for formal asymptotic arguments that suggest that in general solutions of Painleve-IV with large parameters behave either as an algebraic function or an elliptic function. Moreover, the results show that the elliptic approximation holds on the union of a curvilinear rectangle and, in the case of the generalized Okamoto rational solutions, four curvilinear triangles each of which shares an edge with the rectangle; the algebraic approximation is valid in the complementary unbounded domain. We compare the theoretical predictions for the locations of the poles and zeros with numerical plots of the actual poles and zeros obtained from the generating polynomials, and find excellent agreement.
△ Less
Submitted 2 August, 2020;
originally announced August 2020.
-
Using Machine Learning to Identify the Most At-Risk Students in Physics Classes
Authors:
Jie Yang,
Seth DeVore,
Dona Hewagallage,
Paul Miller,
Qing X. Ryan,
John Stewart
Abstract:
Machine learning algorithms have recently been used to predict students' performance in an introductory physics class. The prediction model classified students as those likely to receive an A or B or students likely to receive a grade of C, D, F or withdraw from the class. Early prediction could better allow the direction of educational interventions and the allocation of educational resources. Ho…
▽ More
Machine learning algorithms have recently been used to predict students' performance in an introductory physics class. The prediction model classified students as those likely to receive an A or B or students likely to receive a grade of C, D, F or withdraw from the class. Early prediction could better allow the direction of educational interventions and the allocation of educational resources. However, the performance metrics used in that study become unreliable when used to classify whether a student would receive an A, B or C (the ABC outcome) or if they would receive a D, F or withdraw (W) from the class (the DFW outcome) because the outcome is substantially unbalanced with between 10\% to 20\% of the students receiving a D, F, or W. This work presents techniques to adjust the prediction models and alternate model performance metrics more appropriate for unbalanced outcome variables. These techniques were applied to three samples drawn from introductory mechanics classes at two institutions ($N=7184$, $1683$, and $926$). Applying the same methods as the earlier study produced a classifier that was very inaccurate, classifying only 16\% of the DFW cases correctly; tuning the model increased the DFW classification accuracy to 43\%. Using a combination of institutional and in-class data improved DFW accuracy to 53\% by the second week of class. As in the prior study, demographic variables such as gender, underrepresented minority status, first-generation college student status, and low socioeconomic status were not important variables in the final prediction models.
△ Less
Submitted 23 July, 2020;
originally announced July 2020.
-
An In-Depth Security Assessment of Maritime Container Terminal Software Systems
Authors:
Joseph O. Eichenhofer,
Elisa Heymann,
Barton P. Miller,
Arnold Kang
Abstract:
Attacks on software systems occur world-wide on a daily basis targeting individuals, corporations, and governments alike. The systems that facilitate maritime shipping are at risk of serious disruptions, and these disruptions can stem from vulnerabilities in the software and processes used in these systems. These vulnerabilities leave such systems open to cyber-attack. Assessments of the security…
▽ More
Attacks on software systems occur world-wide on a daily basis targeting individuals, corporations, and governments alike. The systems that facilitate maritime shipping are at risk of serious disruptions, and these disruptions can stem from vulnerabilities in the software and processes used in these systems. These vulnerabilities leave such systems open to cyber-attack. Assessments of the security of maritime shipping systems have focused on identifying risks but have not taken the critical (and expensive) next step of actually identifying vulnerabilities present in these systems. While such risk assessments are important, they have not provided the detailed identification of security issues in the systems that control these ports and their terminals. In response, we formed a key collaboration between an experienced academic cybersecurity team and a well-known commercial software provider that manages maritime shipping. We performed an analysis of the information flow involved in the maritime shipping process, and then executed an in-depth vulnerability assessment of the software that manages freight systems. In this paper, we show the flow of information involved in the freight shipping process and explain how we performed the in-depth assessment, summarizing our findings. Like every large software system, maritime shipping systems have vulnerabilities.
△ Less
Submitted 22 June, 2020;
originally announced June 2020.
-
The ELFIN Mission
Authors:
V. Angelopoulos,
E. Tsai,
L. Bingley,
C. Shaffer,
D. L. Turner,
A. Runov,
W. Li,
J. Liu,
A. V. Artemyev,
X. -J. Zhang,
R. J. Strangeway,
R. E. Wirz,
Y. Y. Shprits,
V. A. Sergeev,
R. P. Caron,
M. Chung,
P. Cruce,
W. Greer,
E. Grimes,
K. Hector,
M. J. Lawson,
D. Leneman,
E. V. Masongsong,
C. L. Russell,
C. Wilkins
, et al. (57 additional authors not shown)
Abstract:
The Electron Loss and Fields Investigation with a Spatio-Temporal Ambiguity-Resolving option (ELFIN-STAR, or simply: ELFIN) mission comprises two identical 3-Unit (3U) CubeSats on a polar (~93deg inclination), nearly circular, low-Earth (~450 km altitude) orbit. Launched on September 15, 2018, ELFIN is expected to have a >2.5 year lifetime. Its primary science objective is to resolve the mechanism…
▽ More
The Electron Loss and Fields Investigation with a Spatio-Temporal Ambiguity-Resolving option (ELFIN-STAR, or simply: ELFIN) mission comprises two identical 3-Unit (3U) CubeSats on a polar (~93deg inclination), nearly circular, low-Earth (~450 km altitude) orbit. Launched on September 15, 2018, ELFIN is expected to have a >2.5 year lifetime. Its primary science objective is to resolve the mechanism of storm-time relativistic electron precipitation, for which electromagnetic ion cyclotron (EMIC) waves are a prime candidate. From its ionospheric vantage point, ELFIN uses its unique pitch-angle-resolving capability to determine whether measured relativistic electron pitch-angle and energy spectra within the loss cone bear the characteristic signatures of scattering by EMIC waves or whether such scattering may be due to other processes. Pairing identical ELFIN satellites with slowly-variable along-track separation allows disambiguation of spatial and temporal evolution of the precipitation over minutes-to-tens-of-minutes timescales, faster than the orbit period of a single low-altitude satellite (~90min). Each satellite carries an energetic particle detector for electrons (EPDE) that measures 50keV to 5MeV electrons with deltaE/E<40% and a fluxgate magnetometer (FGM) on a ~72cm boom that measures magnetic field waves (e.g., EMIC waves) in the range from DC to 5Hz Nyquist (nominally) with <0.3nT/sqrt(Hz) noise at 1Hz. The spinning satellites (T_spin~3s) are equipped with magnetorquers that permit spin-up/down and reorientation maneuvers. The spin axis is placed normal to the orbit plane, allowing full pitch-angle resolution twice per spin. An energetic particle detector for ions (EPDI) measures 250keV-5MeV ions, addressing secondary science. Funded initially by CalSpace and the University Nanosat Program, ELFIN was selected for flight with joint support from NSF and NASA between 2014 and 2018.
△ Less
Submitted 16 June, 2020; v1 submitted 13 June, 2020;
originally announced June 2020.
-
Maximal Spaces for Approximation Rates in $\ell^1$-regularization
Authors:
Philip Miller,
Thorsten Hohage
Abstract:
We study Tikhonov regularization for possibly nonlinear inverse problems with weighted $\ell^1$-penalization. The forward operator, mapping from a sequence space to an arbitrary Banach space, typically an $L^2$-space, is assumed to satisfy a two-sided Lipschitz condition with respect to a weighted $\ell^2$-norm and the norm of the image space. We show that in this setting approximation rates of ar…
▽ More
We study Tikhonov regularization for possibly nonlinear inverse problems with weighted $\ell^1$-penalization. The forward operator, mapping from a sequence space to an arbitrary Banach space, typically an $L^2$-space, is assumed to satisfy a two-sided Lipschitz condition with respect to a weighted $\ell^2$-norm and the norm of the image space. We show that in this setting approximation rates of arbitrarily high Hölder-type order in the regularization parameter can be achieved, and we characterize maximal subspaces of sequences on which these rates are attained. On these subspaces the method also converges with optimal rates in terms of the noise level with the discrepancy principle as parameter choice rule. Our analysis includes the case that the penalty term is not finite at the exact solution ('oversmoothing'). As a standard example we discuss wavelet regularization in Besov spaces $B^r_{1,1}$. In this setting we demonstrate in numerical simulations for a parameter identification problem in a differential equation that our theoretical results correctly predict improved rates of convergence for piecewise smooth unknown coefficients.
△ Less
Submitted 18 October, 2021; v1 submitted 29 May, 2020;
originally announced May 2020.
-
Parallel Binary Code Analysis
Authors:
Xiaozhu Meng,
Jonathon M. Anderson,
John Mellor-Crummey,
Mark W. Krentel,
Barton P. Miller,
Srđan Milaković
Abstract:
Binary code analysis is widely used to assess a program's correctness, performance, and provenance. Binary analysis applications often construct control flow graphs, analyze data flow, and use debugging information to understand how machine code relates to source lines, inlined functions, and data types. To date, binary analysis has been single-threaded, which is too slow for applications such as…
▽ More
Binary code analysis is widely used to assess a program's correctness, performance, and provenance. Binary analysis applications often construct control flow graphs, analyze data flow, and use debugging information to understand how machine code relates to source lines, inlined functions, and data types. To date, binary analysis has been single-threaded, which is too slow for applications such as performance analysis and software forensics, where it is becoming common to analyze binaries that are gigabytes in size and in large batches that contain thousands of binaries.
This paper describes our design and implementation for accelerating the task of constructing control flow graphs (CFGs) from binaries with multithreading. Existing research focuses on addressing challenging code constructs encountered during constructing CFGs, including functions sharing code, jump table analysis, non-returning functions, and tail calls. However, existing analyses do not consider the complex interactions between concurrent analysis of shared code, making it difficult to extend existing serial algorithms to be parallel. A systematic methodology to guide the design of parallel algorithms is essential. We abstract the task of constructing CFGs as repeated applications of several core CFG operations regarding to creating functions, basic blocks, and edges. We then derive properties among CFG operations, including operation dependency, commutativity, monotonicity. These operation properties guide our design of a new parallel analysis for constructing CFGs. We achieved as much as 25$\times$ speedup for constructing CFGs on 64 hardware threads. Binary analysis applications are significantly accelerated with the new parallel analysis: we achieve 8$\times$ for a performance analysis tool and 7$\times$ for a software forensic tool with 16 hardware threads.
△ Less
Submitted 16 May, 2020; v1 submitted 28 January, 2020;
originally announced January 2020.
-
Exploring the Structure of Misconceptions in the Force and Motion Conceptual Evaluation with Modified Module Analysis
Authors:
James Wells,
Rachel Henderson,
Adrienne Traxler,
Paul Miller,
John Stewart
Abstract:
Investigating student learning and understanding of conceptual physics is a primary research area within Physics Education Research (PER). Multiple quantitative methods have been employed to analyze commonly used mechanics conceptual inventories: the Force Concept Inventory (FCI) and the Force and Motion Conceptual Evaluation (FMCE). Recently, researchers have applied network analytic techniques t…
▽ More
Investigating student learning and understanding of conceptual physics is a primary research area within Physics Education Research (PER). Multiple quantitative methods have been employed to analyze commonly used mechanics conceptual inventories: the Force Concept Inventory (FCI) and the Force and Motion Conceptual Evaluation (FMCE). Recently, researchers have applied network analytic techniques to explore the structure of the incorrect responses to the FCI identifying communities of incorrect responses which could be mapped on to common misconceptions. In this study, the method used to analyze the FCI, Modified Module Analysis (MMA), was applied to a large sample of FMCE pretest and post-test responses ($N_{pre}=3956$, $N_{post}=3719$). The communities of incorrect responses identified were consistent with the item groups described in previous works. As in the work with the FCI, the network was simplified by only retaining nodes selected by a substantial number of students. Retaining nodes selected by 20\% of the students produced communities associated with only four misconceptions. The incorrect response communities identified for men and women were substantially different, as was the change in these communities from pretest to post-test. The 20% threshold was far more restrictive than the 4% threshold applied to the FCI in the prior work which generated similar structures. Retaining nodes selected by 5% or 10% of students generated a large number of complex communities. The communities identified at the 10\% threshold were generally associated with common misconceptions producing a far richer set of incorrect communities than the FCI; this may indicate that the FMCE is a superior instrument for characterizing the breadth of student misconceptions about Newtonian mechanics.
△ Less
Submitted 3 January, 2020;
originally announced January 2020.
-
Universality near the gradient catastrophe point in the semiclassical sine-Gordon equation
Authors:
Bing-Ying Lu,
Peter D. Miller
Abstract:
We study the semiclassical limit of the sine-Gordon (sG) equation with below threshold pure impulse initial data of Klaus-Shaw type. The Whitham averaged approximation of this system exhibits a gradient catastrophe in finite time. In accordance with a conjecture of Dubrovin, Grava and Klein, we found that in a $\mathcal{O}(ε^{4/5})$ neighborhood near the gradient catastrophe point, the asymptotics…
▽ More
We study the semiclassical limit of the sine-Gordon (sG) equation with below threshold pure impulse initial data of Klaus-Shaw type. The Whitham averaged approximation of this system exhibits a gradient catastrophe in finite time. In accordance with a conjecture of Dubrovin, Grava and Klein, we found that in a $\mathcal{O}(ε^{4/5})$ neighborhood near the gradient catastrophe point, the asymptotics of the sG solution are universally described by the Painlevé I tritronquée solution. A linear map can be explicitly made from the tritronquée solution to this neighborhood. Under this map: away from the tritronquée poles, the first correction of sG is universally given by the real part of the Hamiltonian of the tritronquée solution; localized defects appear at locations mapped from the poles of tritronquée solution; the defects are proved universally to be a two parameter family of special localized solutions on a periodic background for the sG equation. We are able to characterize the solution in detail. Our approach is the rigorous steepest descent method for matrix Riemann--Hilbert problems, substantially generalizing Bertola and Tovbis's results on the nonlinear Schrödinger equation to establish universality beyond the context of solutions of a single equation.
△ Less
Submitted 19 December, 2019;
originally announced December 2019.
-
Feature2Vec: Distributional semantic modelling of human property knowledge
Authors:
Steven Derby,
Paul Miller,
Barry Devereux
Abstract:
Feature norm datasets of human conceptual knowledge, collected in surveys of human volunteers, yield highly interpretable models of word meaning and play an important role in neurolinguistic research on semantic cognition. However, these datasets are limited in size due to practical obstacles associated with exhaustively listing properties for a large number of words. In contrast, the development…
▽ More
Feature norm datasets of human conceptual knowledge, collected in surveys of human volunteers, yield highly interpretable models of word meaning and play an important role in neurolinguistic research on semantic cognition. However, these datasets are limited in size due to practical obstacles associated with exhaustively listing properties for a large number of words. In contrast, the development of distributional modelling techniques and the availability of vast text corpora have allowed researchers to construct effective vector space models of word meaning over large lexicons. However, this comes at the cost of interpretable, human-like information about word meaning. We propose a method for mapping human property knowledge onto a distributional semantic space, which adapts the word2vec architecture to the task of modelling concept features. Our approach gives a measure of concept and feature affinity in a single semantic space, which makes for easy and efficient ranking of candidate human-derived semantic properties for arbitrary words. We compare our model with a previous approach, and show that it performs better on several evaluation tasks. Finally, we discuss how our method could be used to develop efficient sampling techniques to extend existing feature norm datasets in a reliable way.
△ Less
Submitted 29 August, 2019;
originally announced August 2019.