Showing 1–42 of 42 results for author: Bailey, L

Energy, strength, and alpha width measurements of $E_{\rm{c.m.}} = 1323$ and $1487$ keV resonances in $^{15}$N($α,γ$)$^{19}$F

Authors: R. Fang, J. Görres, R. J. deBoer, S. Moylan, A. Sanchez, T. L. Bailey, S. Carmichael, J. Koros, K. Lee, K. Manukyan, M. Matney, J. P. McDonaugh, D. Robertson, J. Rufino, E. Stech, M. Couder

Abstract: The $^{15}$N($α,γ$)$^{19}$F reaction produces $^{19}$F in asymptotic giant branch (AGB) stars, where the low energy tails of two resonances at $E_{\rm{c.m.}} = 1323 \pm 2$ and $1487 \pm 1.7$ keV are estimated to contribute about $30\%$ of the total reaction rate in these environments. However, recent measurements have shown discrepancies in the energies, the strengths, and the corresponding alpha… ▽ More The $^{15}$N($α,γ$)$^{19}$F reaction produces $^{19}$F in asymptotic giant branch (AGB) stars, where the low energy tails of two resonances at $E_{\rm{c.m.}} = 1323 \pm 2$ and $1487 \pm 1.7$ keV are estimated to contribute about $30\%$ of the total reaction rate in these environments. However, recent measurements have shown discrepancies in the energies, the strengths, and the corresponding alpha widths of these two resonances, resulting in an increase in the systematic uncertainty of the extrapolated cross section to helium burning energies. With this motivation, we have undertaken new measurements of the $^{15}$N$(α,γ)^{19}$F at the University of Notre Dame Nuclear Science Laboratory. The setup consisted of an alpha particle beam impinged on a solid Ti$^{15}$N target with gamma-ray spectroscopy accomplished using a high purity germanium detector. Using the Doppler corrected gamma-ray energies, we confirmed the lower resonance energy to be $1321.6 \pm 0.6$ keV and found a value for the higher one of $1479.4 \pm 0.6$ keV that is more consistent with those found from previous elastic scattering studies. We found that the resonance strengths for both were consistent with most values found in the literature, but a larger alpha width has been recommended for the $E_{\rm{c.m.}} = 1487$ keV resonance. The larger alpha width suggests a reaction rate increase of about $15\%$ at temperatures $T < 0.1$ GK relevant to low mass AGB stars. The impact of the increased reaction rate requires further investigations. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 12 pages, 7 figures, 3 tables. Manuscript submitted to Phys. Rev. C on March 21, 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.

Comments: 48 pages, 29 figures; 4 pages of Supplement Material; version accepted for publication in Physical Review D

Report number: FERMILAB-PUB-24-0084-AD-CSAID-PPD

Pairwise Rearrangement is Fixed-Parameter Tractable in the Single Cut-and-Join Model

Authors: Lora Bailey, Heather Smith Blake, Garner Cochran, Nathan Fox, Michael Levet, Reem Mahmoud, Inne Singgih, Grace Stadnyk, Alexander Wiedemann

Abstract: Genome rearrangement is a common model for molecular evolution. In this paper, we consider the Pairwise Rearrangement problem, which takes as input two genomes and asks for the number of minimum-length sequences of permissible operations transforming the first genome into the second. In the Single Cut-and-Join model (Bergeron, Medvedev, & Stoye, J. Comput. Biol. 2010), Pairwise Rearrangement is… ▽ More Genome rearrangement is a common model for molecular evolution. In this paper, we consider the Pairwise Rearrangement problem, which takes as input two genomes and asks for the number of minimum-length sequences of permissible operations transforming the first genome into the second. In the Single Cut-and-Join model (Bergeron, Medvedev, & Stoye, J. Comput. Biol. 2010), Pairwise Rearrangement is $\#\textsf{P}$-complete (Bailey, et. al., COCOON 2023), which implies that exact sampling is intractable. In order to cope with this intractability, we investigate the parameterized complexity of this problem. We exhibit a fixed-parameter tractable algorithm with respect to the number of components in the adjacency graph that are not cycles of length $2$ or paths of length $1$. As a consequence, we obtain that Pairwise Rearrangement in the Single Cut-and-Join model is fixed-parameter tractable by distance. Our results suggest that the number of nontrivial components in the adjacency graph serves as the key obstacle for efficient sampling. △ Less

Submitted 23 April, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: Full version of paper to appear in SWAT 2024; arXiv admin note: text overlap with arXiv:2305.01851

MSC Class: 92-08; 92D10; 92D20; 68Q17 ACM Class: F.2.2

Tensor Trust: Interpretable Prompt Injection Attacks from an Online Game

Authors: Sam Toyer, Olivia Watkins, Ethan Adrian Mendes, Justin Svegliato, Luke Bailey, Tiffany Wang, Isaac Ong, Karim Elmaaroufi, Pieter Abbeel, Trevor Darrell, Alan Ritter, Stuart Russell

Abstract: While Large Language Models (LLMs) are increasingly being used in real-world applications, they remain vulnerable to prompt injection attacks: malicious third party prompts that subvert the intent of the system designer. To help researchers study this problem, we present a dataset of over 126,000 prompt injection attacks and 46,000 prompt-based "defenses" against prompt injection, all created by p… ▽ More While Large Language Models (LLMs) are increasingly being used in real-world applications, they remain vulnerable to prompt injection attacks: malicious third party prompts that subvert the intent of the system designer. To help researchers study this problem, we present a dataset of over 126,000 prompt injection attacks and 46,000 prompt-based "defenses" against prompt injection, all created by players of an online game called Tensor Trust. To the best of our knowledge, this is currently the largest dataset of human-generated adversarial examples for instruction-following LLMs. The attacks in our dataset have a lot of easily interpretable stucture, and shed light on the weaknesses of LLMs. We also use the dataset to create a benchmark for resistance to two types of prompt injection, which we refer to as prompt extraction and prompt hijacking. Our benchmark results show that many models are vulnerable to the attack strategies in the Tensor Trust dataset. Furthermore, we show that some attack strategies from the dataset generalize to deployed LLM-based applications, even though they have a very different set of constraints to the game. We release all data and source code at https://tensortrust.ai/paper △ Less

Submitted 2 November, 2023; originally announced November 2023.

Checking and Automating Confidentiality Theory in Isabelle/UTP

Authors: Lex Bailey, Jim Woodcock, Simon Foster, Roberto Metere

Abstract: The severity of recent vulnerabilities discovered on modern CPUs, e.g., Spectre [1], highlights how information leakage can have devas-tating effects to the security of computer systems. At the same time, it suggests that confidentiality should be promoted as a normal part of program verification, to discover and mitigate such vulnerabili-ties early in development. The theory we propose is primari… ▽ More The severity of recent vulnerabilities discovered on modern CPUs, e.g., Spectre [1], highlights how information leakage can have devas-tating effects to the security of computer systems. At the same time, it suggests that confidentiality should be promoted as a normal part of program verification, to discover and mitigate such vulnerabili-ties early in development. The theory we propose is primarily based on Bank's theory [2], a framework for reasoning about confidentiali-ty properties formalised in the Unifying Theories of Programming (UTP) [3]. We mechanised our encoding in the current implementa-tion of UTP in the Isabelle theorem prover, Isabelle/UTP [4]. We have identified some theoretical issues in Bank's original framework. Finally, we demonstrate how our mechanisation can be used to for-mally verify of some of the examples from Bank's work. △ Less

Submitted 7 September, 2023; originally announced October 2023.

Image Hijacks: Adversarial Images can Control Generative Models at Runtime

Authors: Luke Bailey, Euan Ong, Stuart Russell, Scott Emmons

Abstract: Are foundation models secure against malicious actors? In this work, we focus on the image input to a vision-language model (VLM). We discover image hijacks, adversarial images that control the behaviour of VLMs at inference time, and introduce the general Behaviour Matching algorithm for training image hijacks. From this, we derive the Prompt Matching method, allowing us to train hijacks matching… ▽ More Are foundation models secure against malicious actors? In this work, we focus on the image input to a vision-language model (VLM). We discover image hijacks, adversarial images that control the behaviour of VLMs at inference time, and introduce the general Behaviour Matching algorithm for training image hijacks. From this, we derive the Prompt Matching method, allowing us to train hijacks matching the behaviour of an arbitrary user-defined text prompt (e.g. 'the Eiffel Tower is now located in Rome') using a generic, off-the-shelf dataset unrelated to our choice of prompt. We use Behaviour Matching to craft hijacks for four types of attack, forcing VLMs to generate outputs of the adversary's choice, leak information from their context window, override their safety training, and believe false statements. We study these attacks against LLaVA, a state-of-the-art VLM based on CLIP and LLaMA-2, and find that all attack types achieve a success rate of over 80%. Moreover, our attacks are automated and require only small image perturbations. △ Less

Submitted 22 April, 2024; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: Project page at https://image-hijacks.github.io

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.

Comments: 8 pages, 3 figures

Report number: FERMILAB-PUB-23-385-AD-CSAID-PPD

Journal ref: Phys. Rev. Lett. 131, 161802 (2023)

Complexity and Enumeration in Models of Genome Rearrangement

Authors: Lora Bailey, Heather Smith Blake, Garner Cochran, Nathan Fox, Michael Levet, Reem Mahmoud, Elizabeth Matson, Inne Singgih, Grace Stadnyk, Xinyi Wang, Alexander Wiedemann

Abstract: In this paper, we examine the computational complexity of enumeration in certain genome rearrangement models. We first show that the Pairwise Rearrangement problem in the Single Cut-and-Join model (Bergeron, Medvedev, & Stoye, J. Comput. Biol. 2010) is $\#\textsf{P}$-complete under polynomial-time Turing reductions. Next, we show that in the Single Cut or Join model (Feijao & Meidanis, IEEE ACM Tr… ▽ More In this paper, we examine the computational complexity of enumeration in certain genome rearrangement models. We first show that the Pairwise Rearrangement problem in the Single Cut-and-Join model (Bergeron, Medvedev, & Stoye, J. Comput. Biol. 2010) is $\#\textsf{P}$-complete under polynomial-time Turing reductions. Next, we show that in the Single Cut or Join model (Feijao & Meidanis, IEEE ACM Trans. Comp. Biol. Bioinf. 2011), the problem of enumerating all medians ($\#$Median) is logspace-computable ($\textsf{FL}$), improving upon the previous polynomial-time ($\textsf{FP}$) bound of Miklós & Smith (RECOMB 2015). △ Less

Submitted 23 April, 2024; v1 submitted 2 May, 2023; originally announced May 2023.

Comments: Full version of paper that appeared in COCOON 2023: https://doi.org/10.1007/978-3-031-49190-0_1

MSC Class: 92-08; 92D10; 92D20; 68Q17 ACM Class: F.2.2

The effect of loops on the mean square displacement of Rouse-model chromatin

Authors: Tianyu Yuan, Hao Yan, Mary Lou P. Bailey, Jessica F. Williams, Ivan Surovtsev, Megan C. King, Simon G. J. Mochrie

Abstract: Many researchers have been encouraged to describe the dynamics of chromosomal loci in chromatin using the classical Rouse model of polymer dynamics by the agreement between the measured mean square displacement (MSD) versus time of fluorescently-labelled loci and the Rouse-model predictions. However, the discovery of intermediate-scale chromatin organization, known as topologically associating dom… ▽ More Many researchers have been encouraged to describe the dynamics of chromosomal loci in chromatin using the classical Rouse model of polymer dynamics by the agreement between the measured mean square displacement (MSD) versus time of fluorescently-labelled loci and the Rouse-model predictions. However, the discovery of intermediate-scale chromatin organization, known as topologically associating domains (TADs), together with the proposed explanation of TADs in terms of chromatin loops and loop extrusion, is at odds with the classical Rouse model, which does not contain loops. Accordingly, we introduce an extended Rouse model that incorporates chromatin loop configurations from loop-extrusion-factor-model simulations. Specifically, we extend the classical Rouse model by modifying the polymer's dynamical matrix to incorporate extra springs that represent loop bases. We also theoretically generalize the friction coefficient matrix so that the Rouse beads with non-uniform friction coefficients are compatible with our Rouse model simulation method. This extended Rouse model allowes us to investigate the impact of loops and loop extrusion on the dynamics of chromatin. We show that loops significantly suppress the averaged MSD of a chromosomal locus, consistent with recent experiments that track fluorescently-labelled chromatin loci in fission yeast [M. L. P. Bailey, I. Surovtsev, J. F. Williams, H. Yan, T. Yuan, S. G. Mochrie, and M. C. King, Mol. Biol. Cell (in press)]. We also find that loops slightly reduce the MSD's stretching exponent from the classical Rouse-model value of 0.5 to a loop-density-dependent value in the 0.45-0.40 range. Remarkably, stretching exponent values in this range have also been reported in recent experiments [S. C. Weber, A. J. Spakowitz, and J. A. Theriot, Phys. Rev. Lett. 104, 238102 (2010) and Bailey et al., Mol. Biol. Cell (in press)]. △ Less

Submitted 21 April, 2023; originally announced April 2023.

SpeedLimit: Neural Architecture Search for Quantized Transformer Models

Authors: Yuji Chai, Luke Bailey, Yunho Jin, Matthew Karle, Glenn G. Ko, David Brooks, Gu-Yeon Wei, H. T. Kung

Abstract: While research in the field of transformer models has primarily focused on enhancing performance metrics such as accuracy and perplexity, practical applications in industry often necessitate a rigorous consideration of inference latency constraints. Addressing this challenge, we introduce SpeedLimit, a novel Neural Architecture Search (NAS) technique that optimizes accuracy whilst adhering to an u… ▽ More While research in the field of transformer models has primarily focused on enhancing performance metrics such as accuracy and perplexity, practical applications in industry often necessitate a rigorous consideration of inference latency constraints. Addressing this challenge, we introduce SpeedLimit, a novel Neural Architecture Search (NAS) technique that optimizes accuracy whilst adhering to an upper-bound latency constraint. Our method incorporates 8-bit integer quantization in the search process to outperform the current state-of-the-art technique. Our results underline the feasibility and efficacy of seeking an optimal balance between performance and latency, providing new avenues for deploying state-of-the-art transformer models in latency-sensitive environments. △ Less

Submitted 13 October, 2023; v1 submitted 24 September, 2022; originally announced September 2022.

Common Patterns in Block-Based Robot Programs

Authors: Florian Obermüller, Robert Pernerstorfer, Lisa Bailey, Ute Heuer, Gordon Fraser

Abstract: Programmable robots are engaging and fun to play with, interact with the real world, and are therefore well suited to introduce young learners to programming. Introductory robot programming languages often extend existing block-based languages such as Scratch. While teaching programming with such languages is well established, the interaction with the real world in robot programs leads to specific… ▽ More Programmable robots are engaging and fun to play with, interact with the real world, and are therefore well suited to introduce young learners to programming. Introductory robot programming languages often extend existing block-based languages such as Scratch. While teaching programming with such languages is well established, the interaction with the real world in robot programs leads to specific challenges, for which learners and educators may require assistance and feedback. A practical approach to provide this feedback is by identifying and pointing out patterns in the code that are indicative of good or bad solutions. While such patterns have been defined for regular block-based programs, robot-specific programming aspects have not been considered so far. The aim of this paper is therefore to identify patterns specific to robot programming for the Scratch-based mBlock programming language, which is used for the popular mBot and Codey Rocky robots. We identify: (1) 26 bug patterns, which indicate erroneous code; (2) three code smells, which indicate code that may work but is written in a confusing or difficult to understand way; and (3) 18 code perfumes, which indicate aspects of code that are likely good. We extend the LitterBox analysis framework to automatically identify these patterns in mBlock programs. Evaluated on a dataset of 3,540 mBlock programs, we find a total of 6,129 instances of bug patterns, 592 code smells and 14,495 code perfumes. This demonstrates the potential of our approach to provide feedback and assistance to learners and educators alike for their mBlock robot programs. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: To be published in the Proceedings of the 17th Workshop in Primary and Secondary Computing Education (WIPSCE'22)

MSC Class: 97P50 ACM Class: D.2.5; K.3.2

Automatic Detection of Interplanetary Coronal Mass Ejections in Solar Wind In Situ Data

Authors: Hannah T. Rüdisser, Andreas Windisch, Ute V. Amerstorfer, Christian Möstl, Tanja Amerstorfer, Rachel L. Bailey, Martin A. Reiss

Abstract: Interplanetary coronal mass ejections (ICMEs) are one of the main drivers for space weather disturbances. In the past, different approaches have been used to automatically detect events in existing time series resulting from solar wind in situ observations. However, accurate and fast detection still remains a challenge when facing the large amount of data from different instruments. For the automa… ▽ More Interplanetary coronal mass ejections (ICMEs) are one of the main drivers for space weather disturbances. In the past, different approaches have been used to automatically detect events in existing time series resulting from solar wind in situ observations. However, accurate and fast detection still remains a challenge when facing the large amount of data from different instruments. For the automatic detection of ICMEs we propose a pipeline using a method that has recently proven successful in medical image segmentation. Comparing it to an existing method, we find that while achieving similar results, our model outperforms the baseline regarding training time by a factor of approximately 20, thus making it more applicable for other datasets. The method has been tested on in situ data from the Wind spacecraft between 1997 and 2015 with a True Skill Statistic (TSS) of 0.64. Out of the 640 ICMEs, 466 were detected correctly by our algorithm, producing a total of 254 False Positives. Additionally, it produced reasonable results on datasets with fewer features and smaller training sets from Wind, STEREO-A and STEREO-B with True Skill Statistics of 0.56, 0.57 and 0.53, respectively. Our pipeline manages to find the start of an ICME with a mean absolute error (MAE) of around 2 hours and 56 minutes, and the end time with a MAE of 3 hours and 20 minutes. The relatively fast training allows straightforward tuning of hyperparameters and could therefore easily be used to detect other structures and phenomena in solar wind data, such as corotating interaction regions. △ Less

Submitted 29 October, 2022; v1 submitted 7 May, 2022; originally announced May 2022.

Journal ref: Space Weather, 20, e2022SW003149

Writhed Analytical Magnetic Flux Rope Model

Authors: Andreas J. Weiss, Teresa Nieves-Chinchilla, Christian Möstl, Martin A. Reiss, Tanja Amerstorfer, Rachel L. Bailey

Abstract: Observations of magnetic clouds, within interplanetary coronal mass ejections (ICMEs), are often well described by flux rope models. Most of these assume either a cylindrical or toroidal geometry. In some cases, these models are also capable of accounting for non-axisymmetric crosssections but they generally all assume axial invariance. It can be expected that any ICME, and its flux rope, will be… ▽ More Observations of magnetic clouds, within interplanetary coronal mass ejections (ICMEs), are often well described by flux rope models. Most of these assume either a cylindrical or toroidal geometry. In some cases, these models are also capable of accounting for non-axisymmetric crosssections but they generally all assume axial invariance. It can be expected that any ICME, and its flux rope, will be deformed along its axis due to influences such as the solar wind. In this work, we aim to develop a writhed analytical magnetic flux rope model which would allow us to analytically describe a flux rope structure with varying curvature and torsion so that we are no longer constrained to a cylindrical or toroidal geometry. In this first iteration of our model we will solely focus on a circular cross-section of constant size. We describe our flux rope geometry in terms of a parametrized flux rope axis and a parallel transport frame. We derive expressions for the axial and poloidal magnetic field components under the assumption that the total axial magnetic flux is conserved. We find an entire class of possible solutions, which differ by the choice of integration constants, and present the results for a specific example. In general, we find that the twist of the magnetic field locally changes when the geometry deviates from a cylinder or torus. This new approach also allows us to generate completely new types of in situ magnetic field profiles which strongly deviate from those generated by cylindrical or toroidal models. △ Less

Submitted 2 December, 2022; v1 submitted 21 February, 2022; originally announced February 2022.

Comments: 19 pages, 4 figures

Uncovering diffusive states of the yeast proton pump, Pma1, and how labeling method can change diffusive behavior

Authors: Mary Lou P. Bailey, Susan E. Pratt, Yongdeng Zhang, Michael Hinrichsen, Joerg Bewersdorf, Lynne J. Regan, Simon G. J. Mochrie

Abstract: We present and analyze video-microscopy-based single-particle-tracking measurements of the budding yeast (Saccharomyces cerevisiae) membrane protein, Pma1, fluorescently-labeled either by direct fusion to the switchable fluorescent protein, mEos3.2, or by a novel, light-touch, labeling scheme, in which a 5 amino acid tag is directly fused to the C-terminus of Pma1, which then binds mEos3.2. The di… ▽ More We present and analyze video-microscopy-based single-particle-tracking measurements of the budding yeast (Saccharomyces cerevisiae) membrane protein, Pma1, fluorescently-labeled either by direct fusion to the switchable fluorescent protein, mEos3.2, or by a novel, light-touch, labeling scheme, in which a 5 amino acid tag is directly fused to the C-terminus of Pma1, which then binds mEos3.2. The diffusivity distributions of these two populations of single particle tracks differ significantly, demonstrating that labeling method can be an important determinant of diffusive behavior. We also applied perturbation expectation maximization (pEMv2) [Physical Review E 94, 052412 (2016)], which sorts trajectories into the statistically-optimum number of diffusive states. For both TRAP-labeled Pma1 and Pma1-mEos3.2, pEMv2 sorts the tracks into two diffusive states: an essentially immobile state and a more mobile state. However, the mobile fraction of Pma1-mEos3.2 tracks is much smaller (0.1) than the mobile fraction of TRAP-labeled Pma1 tracks (0.5). In addition, the diffusivity of Pma1-mEos3.2's mobile state is several times smaller than the diffusivity of TRAP-labeled Pma1's mobile state. To critically assess pEMv2's performance, we compare the diffusivity and covariance distributions of the experimental pEMv2-sorted populations to corresponding theoretical distributions, assuming that Pma1 displacements realize a Gaussian random process. The experiment-theory comparisons for both the TRAP-labeled Pma1 and Pma1-mEos3.2 reveal good agreement, bolstering the pEMv2 approach. △ Less

Submitted 30 June, 2022; v1 submitted 23 January, 2022; originally announced January 2022.

Comments: 19 pages, 14 figures

Forecasting GICs and geoelectric fields from solar wind data using LSTMs: application in Austria

Authors: R. L. Bailey, R. Leonhardt, C. Möstl, C. Beggan, M. A. Reiss, A. Bhaskar, A. J. Weiss

Abstract: The forecasting of local GIC effects has largely relied on the forecasting of dB/dt as a proxy and, to date, little attention has been paid to directly forecasting the geoelectric field or GICs themselves. We approach this problem with machine learning tools, specifically recurrent neural networks or LSTMs by taking solar wind observations as input and training the models to predict two different… ▽ More The forecasting of local GIC effects has largely relied on the forecasting of dB/dt as a proxy and, to date, little attention has been paid to directly forecasting the geoelectric field or GICs themselves. We approach this problem with machine learning tools, specifically recurrent neural networks or LSTMs by taking solar wind observations as input and training the models to predict two different kinds of output: first, the geoelectric field components Ex and Ey; and second, the GICs in specific substations in Austria. The training is carried out on the geoelectric field and GICs modelled from 26 years of one-minute geomagnetic field measurements, and results are compared to GIC measurements from recent years. The GICs are generally predicted better by an LSTM trained on values from a specific substation, but only a fraction of the largest GICs are correctly predicted. This model had a correlation with measurements of around 0.6, and a root-mean-square error of 0.7 A. The probability of detecting mild activity in GICs is around 50%, and 15% for larger GICs. △ Less

Submitted 16 February, 2022; v1 submitted 14 September, 2021; originally announced September 2021.

Multipoint interplanetary coronal mass ejections observed with Solar Orbiter, BepiColombo, Parker Solar Probe, Wind and STEREO-A

Authors: C. Möstl, A. J. Weiss, M. A. Reiss, T. Amerstorfer, R. L. Bailey, J. Hinterreiter, M. Bauer, D. Barnes, J. A. Davies, R. A. Harrison, J. L. Freiherr von Forstner, E. E. Davies, D. Heyner, T. Horbury, S. D. Bale

Abstract: We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP),… ▽ More We report the result of the first search for multipoint in situ and imaging observations of interplanetary coronal mass ejections (ICMEs) starting with the first Solar Orbiter (SolO) data in 2020 April - 2021 April. A data exploration analysis is performed including visualizations of the magnetic field and plasma observations made by the five spacecraft SolO, BepiColombo, Parker Solar Probe (PSP), Wind and STEREO-A, in connection with coronagraph and heliospheric imaging observations from STEREO-A/SECCHI and SOHO/LASCO. We identify ICME events that could be unambiguously followed with the STEREO-A heliospheric imagers during their interplanetary propagation to their impact at the aforementioned spacecraft, and look for events where the same ICME is seen in situ by widely separated spacecraft. We highlight two events: (1) a small streamer blowout CME on 2020 June 23 observed with a triple lineup by PSP, BepiColombo and Wind, guided by imaging with STEREO-A, and (2) the first fast CME of solar cycle 25 ($ \approx 1600$ km s$^{-1}$) on 2020 November 29 observed in situ by PSP and STEREO-A. These results are useful for modeling the magnetic structure of ICMEs and the interplanetary evolution and global shape of their flux ropes and shocks, and for studying the propagation of solar energetic particles. The combined data from these missions are already turning out to be a treasure trove for space weather research and are expected to become even more valuable with an increasing number of ICME events expected during the rise and maximum of solar cycle 25. △ Less

Submitted 4 January, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: in press at ApJ Letters (submitted 2021 September 15, revised 2021 November 24, accepted 2021 December 9). 11 pages, 3 figures, 1 table

Machine learning for predicting the Bz magnetic field component from upstream in situ observations of solar coronal mass ejections

Authors: Martin A. Reiss, Christian Möstl, Rachel L. Bailey, Hannah T. Rüdisser, Ute V. Amerstorfer, Tanja Amerstorfer, Andreas J. Weiss, Jürgen Hinterreiter, Andreas Windisch

Abstract: Predicting the Bz magnetic field embedded within ICMEs, also known as the Bz problem, is a key challenge in space weather forecasting. We study the hypothesis that upstream in situ measurements of the sheath region and the first few hours of the magnetic obstacle provide sufficient information for predicting the downstream Bz component. To do so, we develop a predictive tool based on machine learn… ▽ More Predicting the Bz magnetic field embedded within ICMEs, also known as the Bz problem, is a key challenge in space weather forecasting. We study the hypothesis that upstream in situ measurements of the sheath region and the first few hours of the magnetic obstacle provide sufficient information for predicting the downstream Bz component. To do so, we develop a predictive tool based on machine learning that is trained and tested on 348 ICMEs from Wind, STEREO-A, and STEREO-B measurements. We train the machine learning models to predict the minimum value of the Bz component and the maximum value of the total magnetic field Bz in the magnetic obstacle. To validate the tool, we let the ICMEs sweep over the spacecraft and assess how continually feeding in situ measurements into the tool improves the Bz prediction. Because the application of the tool in operations needs an automated detection of ICMEs, we implement an existing automated ICME detection algorithm and test its robustness for the time intervals under scrutiny. We find that the predictive tool can predict the minimum value of the Bz component in the magnetic obstacle with a mean absolute error of 3.12 nT and a Pearson correlation coefficient of 0.71 when the sheath region and the first 4 hours of the magnetic obstacle are observed. While the underlying hypothesis is unlikely to solve the Bz problem, the tool shows promise for ICMEs that have a recognizable magnetic flux rope signature. Transitioning the tool to operations could lead to improved space weather forecasting. △ Less

Submitted 5 November, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Data Enforced: An Exploratory Impact Analysis of Automated Speed Enforcement in the District of Columbia

Authors: Awad Abdelhalim, Linda Bailey, Emily Dalphy, Kelli Raboy

Abstract: In 2015, the District of Columbia framed a Vision Zero mission and action plan, with a target of achieving zero traffic fatalities by 2024. This study examines the impacts of Automated Speed Enforcement (ASE) and its role in achieving the goals of Vision Zero. Independent datasets containing detailed information about traffic crashes, ASE camera locations, and citation records, and driving speeds… ▽ More In 2015, the District of Columbia framed a Vision Zero mission and action plan, with a target of achieving zero traffic fatalities by 2024. This study examines the impacts of Automated Speed Enforcement (ASE) and its role in achieving the goals of Vision Zero. Independent datasets containing detailed information about traffic crashes, ASE camera locations, and citation records, and driving speeds across the District's streets were collected, combined, and analyzed to identify patterns and trends in crashes, speed limit violations, and speeding behavior before and after the ASE camera installation. The results of this exploratory analysis confirm the safety benefits of ASE systems in Washington, D.C. The study also provides a blueprint for the different means of evaluating the short-term impact of ASE systems using different data sources which can aid practitioners in better evaluating existing systems and support the decision-making process regarding future installations. △ Less

Submitted 18 June, 2021; originally announced June 2021.

Multi point analysis of coronal mass ejection flux ropes using combined data from Solar Orbiter, BepiColombo and Wind

Authors: A. J. Weiss, C. Moestl, E. E. Davies, T. Amerstorfer, M. Bauer, J. Hinterreiter, M. Reiss, R. L. Bailey, T. S. Horbury, H. O'Brien, V. Evans, V. Angelini, D. Heiner, I. Richter, H-U. Auster, W. Magnes, D. Fischer, W. Baumjohann

Abstract: The recent launch of Solar Orbiter and BepiColombo opened a brief window in which these two spacecraft were positioned in a constellation that allows for the detailed sampling of any Earth-directed CMEs. Fortunately, two such events occurred with in situ detections of an ICME by Solar Orbiter on the 19th of April and the 28th of May 2020. These two events were subsequently also observed in situ by… ▽ More The recent launch of Solar Orbiter and BepiColombo opened a brief window in which these two spacecraft were positioned in a constellation that allows for the detailed sampling of any Earth-directed CMEs. Fortunately, two such events occurred with in situ detections of an ICME by Solar Orbiter on the 19th of April and the 28th of May 2020. These two events were subsequently also observed in situ by BepiColombo and Wind around a day later. We attempt to reconstruct the observed in situ magnetic field measurements for all three spacecraft simultaneously using an empirical magnetic flux rope model. This allows us to test the validity of our flux rope model on a larger and more global scale and allows for cross-validation of the analysis with different spacecraft combinations. Finally, we can also compare the results from the in situ modeling to remote observations obtained from the STEREO-A heliospheric imagers. We make use of the 3D coronal rope ejection model in order to simulate the ICME evolution. We adapt a previously developed ABC-SMC fitting algorithm for the application to multi point scenarios. We show that we are able to generally reconstruct the flux ropes signatures at three different spacecraft positions simultaneously using our model in combination with the flux rope fitting algorithm. For the well-behaved 19th of April ICME our approach works very well. The 28th of May ICME, on the other hand, shows the limitations of our approach. Unfortunately, the usage of multi-point observations for these events does not appear to solve inherent issues, such as the estimation of the magnetic field twist or flux rope aspect-ratios due to the specific constellation of the spacecraft positions. As our general approach can be used for any fast forward simulation-based model we give a blueprint for future studies using more advanced ICME models. △ Less

Submitted 17 May, 2021; v1 submitted 30 March, 2021; originally announced March 2021.

Comments: 11 pages, 7 figures

Journal ref: A&A 656, A13 (2021)

Measurement of the generalized spin polarizabilities of the neutron in the low $Q^2$ region

Authors: V. Sulkosky, C. Peng, J. -P. Chen, A. Deur, S. Abrahamyan, K. A. Aniol, D. S. Armstrong, T. Averett, S. L. Bailey, A. Beck, P. Bertin, F. Butaru, W. Boeglin, A. Camsonne, G. D. Cates, C. C. Chang, Seonho Choi, E. Chudakov, L. Coman, J. C Cornejo, B. Craver, F. Cusanno, R. De Leo, C. W. de Jager, J. D. Denton , et al. (84 additional authors not shown)

Abstract: Understanding the nucleon spin structure in the regime where the strong interaction becomes truly strong poses a challenge to both experiment and theory. At energy scales below the nucleon mass of about 1 GeV, the intense interaction among the quarks and gluons inside the nucleon makes them highly correlated. Their coherent behaviour causes the emergence of effective degrees of freedom, requiring… ▽ More Understanding the nucleon spin structure in the regime where the strong interaction becomes truly strong poses a challenge to both experiment and theory. At energy scales below the nucleon mass of about 1 GeV, the intense interaction among the quarks and gluons inside the nucleon makes them highly correlated. Their coherent behaviour causes the emergence of effective degrees of freedom, requiring the application of non-perturbative techniques, such as chiral effective field theory. Here, we present measurements of the neutron's generalized spin-polarizabilities that quantify the neutron's spin precession under electromagnetic fields at very low energy-momentum transfer squared down to 0.035 GeV$^2$. In this regime, chiral effective field theory calculations are expected to be applicable. Our data, however, show a strong discrepancy with these predictions, presenting a challenge to the current description of the neutron's spin properties. △ Less

Submitted 23 February, 2022; v1 submitted 4 March, 2021; originally announced March 2021.

Comments: V1: initial version submitted to Nature Physics. V2: Published version. 16 pages, 7 figures. Additional material: 4 data tables (18 pages) V3: Typo corrected in author list. Paper content unchanged

Report number: JLAB-PHY-21-3312, DOE/OR/23177-5119

Journal ref: Nature Physics, Vol. 17 687-692 (2021)

Extrusion of chromatin loops by a composite loop extrusion factor

Authors: Hao Yan, Ivan Surovtsev, Jessica F Williams, Mary Lou P Bailey, Megan C King, Simon G J Mochrie

Abstract: Chromatin loop extrusion by Structural Maintenance of Chromosome (SMC) complexes is thought to underlie intermediate-scale chromatin organization inside cells. Motivated by a number of experiments suggesting that nucleosomes may block loop extrusion by SMCs, such as cohesin and condensin complexes, we introduce and characterize theoretically a composite loop extrusion factor (composite LEF) model.… ▽ More Chromatin loop extrusion by Structural Maintenance of Chromosome (SMC) complexes is thought to underlie intermediate-scale chromatin organization inside cells. Motivated by a number of experiments suggesting that nucleosomes may block loop extrusion by SMCs, such as cohesin and condensin complexes, we introduce and characterize theoretically a composite loop extrusion factor (composite LEF) model. In addition to an SMC complex that creates a chromatin loop by encircling two threads of DNA, this model includes a remodeling complex that relocates or removes nucleosomes as it progresses along the chromatin, and nucleosomes that block SMC translocation along the DNA. Loop extrusion is enabled by SMC motion along nucleosome-free DNA, created in the wake of the remodeling complex, while nucleosome re-binding behind the SMC acts as a ratchet, holding the SMC close to the remodeling complex. We show that, for a wide range of parameter values, this collection of factors constitutes a composite LEF that extrudes loops with a velocity, comparable to the velocity of remodeling complex translocation on chromatin in the absence of SMC, and much faster than loop extrusion by an isolated SMC that is blocked by nucleosomes. △ Less

Submitted 1 March, 2021; originally announced March 2021.

Comments: 12 pages, 7 figures

Journal ref: Phys. Rev. E 104, 024414 (2021)

Why are ELEvoHI CME arrival predictions different if based on STEREO-A or STEREO-B heliospheric imager observations?

Authors: Jürgen Hinterreiter, Tanja Amerstorfer, Martin A. Reiss, Christian Möstl, Manuela Temmer, Maike Bauer, Ute V. Amerstorfer, Rachel L. Bailey, Andreas J. Weiss, Jackie A. Davies, Luke A. Barnard, Mathew J. Owens

Abstract: Accurate forecasting of the arrival time and arrival speed of coronal mass ejections (CMEs) is a unsolved problem in space weather research. In this study, a comparison of the predicted arrival times and speeds for each CME based, independently, on the inputs from the two STEREO vantage points is carried out. We perform hindcasts using ELlipse Evolution model based on Heliospheric Imager observati… ▽ More Accurate forecasting of the arrival time and arrival speed of coronal mass ejections (CMEs) is a unsolved problem in space weather research. In this study, a comparison of the predicted arrival times and speeds for each CME based, independently, on the inputs from the two STEREO vantage points is carried out. We perform hindcasts using ELlipse Evolution model based on Heliospheric Imager observations (ELEvoHI) ensemble modelling. An estimate of the ambient solar wind conditions is obtained by the Wang-Sheeley-Arge/Heliospheric Upwind eXtrapolation (WSA/HUX) model combination that serves as input to ELEvoHI. We carefully select 12 CMEs between February 2010 and July 2012 that show clear signatures in both STEREO-A and STEREO-B HI time-elongation maps, that propagate close to the ecliptic plane, and that have corresponding in situ signatures at Earth. We find a mean arrival time difference of 6.5 hrs between predictions from the two different viewpoints, which can reach up to 9.5 hrs for individual CMEs, while the mean arrival speed difference is 63 km s$^{-1}$. An ambient solar wind with a large speed variance leads to larger differences in the STEREO-A and STEREO-B CME arrival time predictions ($cc~=~0.92$). Additionally, we compare the predicted arrivals, from both spacecraft, to the actual in situ arrivals at Earth and find a mean absolute error of 7.5 $\pm$ 9.5 hrs for the arrival time and 87 $\pm$ 111 km s$^{-1}$ for the arrival speed. There is no tendency for one spacecraft to provide more accurate arrival predictions than the other. △ Less

Submitted 15 February, 2021; originally announced February 2021.

In situ multi-spacecraft and remote imaging observations of the first CME detected by Solar Orbiter and BepiColombo

Authors: E. E. Davies, C. Möstl, M. J. Owens, A. J. Weiss, T. Amerstorfer, J. Hinterreiter, M. Bauer, R. L. Bailey, M. A. Reiss, R. J. Forsyth, T. S. Horbury, H. O'Brien, V. Evans, V. Angelini, D. Heyner, I. Richter, H-U. Auster, W. Magnes, W. Baumjohann, D. Fischer, D. Barnes, J. A. Davies, R. A. Harrison

Abstract: On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ i… ▽ More On 2020 April 19 a coronal mass ejection (CME) was detected in situ by Solar Orbiter at a heliocentric distance of about 0.8 AU. The CME was later observed in situ on April 20th by the Wind and BepiColombo spacecraft whilst BepiColombo was located very close to Earth. This CME presents a good opportunity for a triple radial alignment study, as the spacecraft were separated by less than 5$^\circ$ in longitude. The source of the CME, which was launched on April 15th, was an almost entirely isolated streamer blowout. STEREO-A observed the event remotely from -75.1$^\circ$ longitude, which is an exceptionally well suited viewpoint for heliospheric imaging of an Earth directed CME. The configuration of the four spacecraft has provided an exceptionally clean link between remote imaging and in situ observations of the CME. We have used the in situ observations of the CME at Solar Orbiter, Wind, and BepiColombo, and the remote observations of the CME at STEREO-A in combination with flux rope models to determine the global shape of the CME and its evolution as it propagated through the inner heliosphere. A clear flattening of the CME cross-section has been observed by STEREO-A, and further confirmed by comparing profiles of the flux rope models to the in situ data, where the distorted flux rope cross-section qualitatively agrees most with in situ observations of the magnetic field at Solar Orbiter. Comparing in situ observations of the magnetic field between spacecraft, we find that the dependence of the maximum (mean) magnetic field strength decreases with heliocentric distance as $r^{-1.24 \pm 0.50}$ ($r^{-1.12 \pm 0.14}$), in disagreement with previous studies. Further assessment of the axial and poloidal magnetic field strength dependencies suggests that the expansion of the CME is likely neither self-similar nor cylindrically symmetric. △ Less

Submitted 23 February, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Journal ref: A&A 656, A2 (2021)

Lessons Learned from the 1st ARIEL Machine Learning Challenge: Correcting Transiting Exoplanet Light Curves for Stellar Spots

Authors: Nikolaos Nikolaou, Ingo P. Waldmann, Angelos Tsiaras, Mario Morvan, Billy Edwards, Kai Hou Yip, Giovanna Tinetti, Subhajit Sarkar, James M. Dawson, Vadim Borisov, Gjergji Kasneci, Matej Petkovic, Tomaz Stepisnik, Tarek Al-Ubaidi, Rachel Louise Bailey, Michael Granitzer, Sahib Julka, Roman Kern, Patrick Ofner, Stefan Wagner, Lukas Heppe, Mirko Bunse, Katharina Morik

Abstract: The last decade has witnessed a rapid growth of the field of exoplanet discovery and characterisation. However, several big challenges remain, many of which could be addressed using machine learning methodology. For instance, the most prolific method for detecting exoplanets and inferring several of their characteristics, transit photometry, is very sensitive to the presence of stellar spots. The… ▽ More The last decade has witnessed a rapid growth of the field of exoplanet discovery and characterisation. However, several big challenges remain, many of which could be addressed using machine learning methodology. For instance, the most prolific method for detecting exoplanets and inferring several of their characteristics, transit photometry, is very sensitive to the presence of stellar spots. The current practice in the literature is to identify the effects of spots visually and correct for them manually or discard the affected data. This paper explores a first step towards fully automating the efficient and precise derivation of transit depths from transit light curves in the presence of stellar spots. The methods and results we present were obtained in the context of the 1st Machine Learning Challenge organized for the European Space Agency's upcoming Ariel mission. We first present the problem, the simulated Ariel-like data and outline the Challenge while identifying best practices for organizing similar challenges in the future. Finally, we present the solutions obtained by the top-5 winning teams, provide their code and discuss their implications. Successful solutions either construct highly non-linear (w.r.t. the raw data) models with minimal preprocessing -deep neural networks and ensemble methods- or amount to obtaining meaningful statistics from the light curves, constructing linear models on which yields comparably good predictive performance. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: 20 pages, 7 figures, 2 tables, Submitted to The Astrophysics Journal (ApJ)

Covariance Distributions in Single Particle Tracking

Authors: Mary Lou P Bailey, Hao Yan, Ivan Surovtsev, Jessica F Williams, Megan C King, Simon G J Mochrie

Abstract: Several recent experiments, including our own in the fission yeast, S. pombe, have characterized the motions of gene loci within living nuclei by measuring the locus position over time, then proceeding to obtain the statistical properties of this motion. To address the question of whether a population of single particle tracks, obtained from many different cells, corresponds to a single mode of di… ▽ More Several recent experiments, including our own in the fission yeast, S. pombe, have characterized the motions of gene loci within living nuclei by measuring the locus position over time, then proceeding to obtain the statistical properties of this motion. To address the question of whether a population of single particle tracks, obtained from many different cells, corresponds to a single mode of diffusion, we derive theoretical equations describing the probability distribution of the displacement covariance, assuming the displacement is a zero-mean multivariate Gaussian random variable. We also determine the corresponding theoretical means, variances, and third central moments. Bolstering the theory is good agreement between its predictions and the results obtained for various simulated and measured data sets, including simulated particle trajectories of simple and anomalous diffusion, and the measured trajectories of an optically-trapped bead in water, and in a viscoelastic solution. We also show that, for sufficiently long tracks, each covariance distribution in these examples is well-described by a skew-normal distribution with mean, variance, and skewness given by theory. For experimental S. pombe gene locus data, however, we find that the first two covariance distributions are wider than predicted, although the third and subsequent covariances are well-described by theory. This suggests that the origin of the theory-experiment discrepancy is associated with localization noise, which influences only the first two covariances. Thus, we hypothesize that the discrepancy is caused by locus-to-locus heterogeneity in the localization noise. Further simulations reveal excess covariance widths can be largely recreated on the basis of heterogeneous noise. We conclude that the motion of gene loci in fission yeast is consistent with a single mode of diffusion. △ Less

Submitted 11 January, 2021; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: Main text 19 pages, 22 figures (28 pages total)

Journal ref: Phys. Rev. E 103, 032405 (2021)

Analysis of coronal mass ejection flux rope signatures using 3DCORE and approximate Bayesian Computation

Authors: Andreas J. Weiss, Christian Möstl, Tanja Amerstorfer, Rachel L. Bailey, Martin A. Reiss, Jürgen Hinterreiter, Ute A. Amerstorfer, Maike Bauer

Abstract: We present a major update to the 3D coronal rope ejection (3DCORE) technique for modeling coronal mass ejection flux ropes in conjunction with an Approximate Bayesian Computation (ABC) algorithm that is used for fitting the model to in situ magnetic field measurements. The model assumes an empirically motivated torus-like flux rope structure that expands self-similarly within the heliosphere, is i… ▽ More We present a major update to the 3D coronal rope ejection (3DCORE) technique for modeling coronal mass ejection flux ropes in conjunction with an Approximate Bayesian Computation (ABC) algorithm that is used for fitting the model to in situ magnetic field measurements. The model assumes an empirically motivated torus-like flux rope structure that expands self-similarly within the heliosphere, is influenced by a simplified interaction with the solar wind environment, and carries along an embedded analytical magnetic field. The improved 3DCORE implementation allows us to generate extremely large ensemble simulations which we then use to find global best-fit model parameters using an ABC sequential Monte Carlo (SMC) algorithm. The usage of this algorithm, under some basic assumptions on the uncertainty of the magnetic field measurements, allows us to furthermore generate estimates on the uncertainty of model parameters using only a single in situ observation. We apply our model to synthetically generated measurements to prove the validity of our implementation for the fitting procedure. We also present a brief analysis, within the scope of our model, of an event captured by Parker Solar Probe (PSP) shortly after its first fly-by of the Sun on 2018 November 12 at 0.25 AU. The presented toolset is also easily extendable to the analysis of events captured by multiple spacecraft and will therefore facilitate future multi-point studies. △ Less

Submitted 17 November, 2020; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: 18 pages, 11 figures, in press at ApJS, accepted 10 November 2020

An Experience of Introducing Primary School Children to Programming using Ozobots (Practical Report)

Authors: Nina Körber, Lisa Bailey, Luisa Greifenstein, Gordon Fraser, Barbara Sabitzer, Marina Rottenhofer

Abstract: Algorithmic thinking is a central concept in the context of computational thinking, and it is commonly taught by computer programming. A recent trend is to introduce basic programming concepts already very early on at primary school level. There are, however, several challenges in teaching programming at this level: Schools and teachers are often neither equipped nor trained appropriately, and the… ▽ More Algorithmic thinking is a central concept in the context of computational thinking, and it is commonly taught by computer programming. A recent trend is to introduce basic programming concepts already very early on at primary school level. There are, however, several challenges in teaching programming at this level: Schools and teachers are often neither equipped nor trained appropriately, and the best way to move from initial "unplugged" activities to creating programs on a computer are still a matter of open debate. In this paper, we describe our experience of a small INTERREG-project aiming at supporting local primary schools in introducing children to programming concepts using Ozobot robots. These robots have two distinct advantages: First, they can be programmed with and without computers, thus helping the transition from unplugged programming to programming with a computer. Second, they are small and easy to transport, even when used together with tablet computers. Although we learned in our outreach events that the use of Ozobots is not without challenges, our overall experience is positive and can hopefully support others in setting up first encounters with programming at primary schools. △ Less

Submitted 16 August, 2021; v1 submitted 28 August, 2020; originally announced August 2020.

Comments: To be published in the proceedings of the 16th Workshop in Primary and Secondary Computing Education (WIPSCE)

MSC Class: 97P50 ACM Class: K.3.2; D.2.6

Evaluation of CME Arrival Prediction Using Ensemble Modeling Based on Heliospheric Imaging Observations

Authors: Tanja Amerstorfer, Jürgen Hinterreiter, Martin A. Reiss, Christian Möstl, Jackie A. Davies, Rachel L. Bailey, Andreas J. Weiss, Mateja Dumbović, Maike Bauer, Ute V. Amerstorfer, Richard A. Harrison

Abstract: In this study, we evaluate a coronal mass ejection (CME) arrival prediction tool that utilizes the wide-angle observations made by STEREO's heliospheric imagers (HI). The unsurpassable advantage of these imagers is the possibility to observe the evolution and propagation of a CME from close to the Sun out to 1 AU and beyond. We believe that by exploiting this capability, instead of relying on coro… ▽ More In this study, we evaluate a coronal mass ejection (CME) arrival prediction tool that utilizes the wide-angle observations made by STEREO's heliospheric imagers (HI). The unsurpassable advantage of these imagers is the possibility to observe the evolution and propagation of a CME from close to the Sun out to 1 AU and beyond. We believe that by exploiting this capability, instead of relying on coronagraph observations only, it is possible to improve today's CME arrival time predictions. The ELlipse Evolution model based on HI observations (ELEvoHI) assumes that the CME frontal shape within the ecliptic plane is an ellipse, and allows the CME to adjust to the ambient solar wind speed, i.e. it is drag-based. ELEvoHI is used to perform ensemble simulations by varying the CME frontal shape within given boundary conditions that are consistent with the observations made by HI. In this work, we evaluate different set-ups of the model by performing hindcasts for 15 well-defined isolated CMEs that occurred when STEREO was near L4/5, between the end of 2008 and the beginning of 2011. In this way, we find a mean absolute error of between $6.2\pm7.9$ h and $9.9\pm13$ h depending on the model set-up used. ELEvoHI is specified for using data from future space weather missions carrying HIs located at L5 or L1. It can also be used with near real-time STEREO-A HI beacon data to provide CME arrival predictions during the next $\sim7$ years when STEREO-A is observing the Sun-Earth space. △ Less

Submitted 16 February, 2021; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 25 pages, 9 figures, 2 tables; This article also appears in: Heliophysics and Space Weather Studies from the Sun-Earth Lagrange Points

Journal ref: AGU Space Weather, 2021, Volume19, Issue 1

Prediction of the in situ coronal mass ejection rate for solar cycle 25: Implications for Parker Solar Probe in situ observations

Authors: Christian Möstl, Andreas J. Weiss, Rachel L. Bailey, Martin A. Reiss, Tanja Amerstorfer, Jürgen Hinterreiter, Maike Bauer, Scott W. McIntosh, Noé Lugaz, David Stansby

Abstract: The Parker Solar Probe (PSP) and Solar Orbiter missions are designed to make groundbreaking observations of the Sun and interplanetary space within this decade. We show that a particularly interesting in situ observation of an interplanetary coronal mass ejection (ICME) by PSP may arise during close solar flybys ($< 0.1$~AU). During these times, the same magnetic flux rope inside an ICME could be… ▽ More The Parker Solar Probe (PSP) and Solar Orbiter missions are designed to make groundbreaking observations of the Sun and interplanetary space within this decade. We show that a particularly interesting in situ observation of an interplanetary coronal mass ejection (ICME) by PSP may arise during close solar flybys ($< 0.1$~AU). During these times, the same magnetic flux rope inside an ICME could be observed in situ by PSP twice, by impacting its frontal part as well as its leg. Investigating the odds of this situation, we forecast the ICME rate in solar cycle 25 based on 2 models for the sunspot number (SSN): (1) the forecast of an expert panel in 2019 (maximum SSN = 115), and (2) a prediction by McIntosh et al. (2020, maximum SSN = 232). We link the SSN to the observed ICME rates in solar cycles 23 and 24 with the Richardson and Cane list and our own ICME catalog, and calculate that between 1 and 7 ICMEs will be observed by PSP at heliocentric distances $< 0.1$ AU until 2025, including 1$σ$ uncertainties. We then model the potential flux rope signatures of such a double-crossing event with the semi-empirical 3DCORE flux rope model, showing a telltale elevation of the radial magnetic field component $B_R$, and a sign reversal in the component $B_N$ normal to the solar equator compared to field rotation in the first encounter. This holds considerable promise to determine the structure of CMEs close to their origin in the solar corona. △ Less

Submitted 17 September, 2020; v1 submitted 29 July, 2020; originally announced July 2020.

Comments: 11 pages, 6 figures, accepted for publication in the Astrophysical Journal on 2020 September 17

Journal ref: ApJ, 903, 92 (2020)

Using gradient boosting regression to improve ambient solar wind model predictions

Authors: R. L. Bailey, M. A. Reiss, C. N. Arge, C. Möstl, M. J. Owens, U. V. Amerstorfer, C. J. Henney, T. Amerstorfer, A. J. Weiss, J. Hinterreiter

Abstract: Studying the ambient solar wind, a continuous pressure-driven plasma flow emanating from our Sun, is an important component of space weather research. The ambient solar wind flows in interplanetary space determine how solar storms evolve through the heliosphere before reaching Earth, and especially during solar minimum are themselves a driver of activity in the Earth's magnetic field. Accurately f… ▽ More Studying the ambient solar wind, a continuous pressure-driven plasma flow emanating from our Sun, is an important component of space weather research. The ambient solar wind flows in interplanetary space determine how solar storms evolve through the heliosphere before reaching Earth, and especially during solar minimum are themselves a driver of activity in the Earth's magnetic field. Accurately forecasting the ambient solar wind flow is therefore imperative to space weather awareness. Here we present a machine learning approach in which solutions from magnetic models of the solar corona are used to output the solar wind conditions near the Earth. The results are compared to observations and existing models in a comprehensive validation analysis, and the new model outperforms existing models in almost all measures. In addition, this approach offers a new perspective to discuss the role of different input data to ambient solar wind modeling, and what this tells us about the underlying physical processes. The final model discussed here represents an extremely fast, well-validated and open-source approach to the forecasting of ambient solar wind at Earth. △ Less

Submitted 23 March, 2021; v1 submitted 23 June, 2020; originally announced June 2020.

Comments: 15 pages, 6 figures, 3 tables

Prediction of Dst during solar minimum using in situ measurements at L5

Authors: R. L. Bailey, C. Möstl, M. A. Reiss, A. J. Weiss, U. V. Amerstorfer, T. Amerstorfer, J. Hinterreiter, W. Magnes, R. Leonhardt

Abstract: Geomagnetic storms resulting from high-speed streams can have significant negative impacts on modern infrastructure due to complex interactions between the solar wind and geomagnetic field. One measure of the extent of this effect is the Kyoto $Dst$ index. We present a method to predict $Dst$ from data measured at the Lagrange 5 (L5) point, which allows for forecasts of solar wind development 4.5… ▽ More Geomagnetic storms resulting from high-speed streams can have significant negative impacts on modern infrastructure due to complex interactions between the solar wind and geomagnetic field. One measure of the extent of this effect is the Kyoto $Dst$ index. We present a method to predict $Dst$ from data measured at the Lagrange 5 (L5) point, which allows for forecasts of solar wind development 4.5 days in advance of the stream reaching the Earth. Using the STEREO-B satellite as a proxy, we map data measured near L5 to the near-Earth environment and make a prediction of the $Dst$ from this point using the Temerin-Li $Dst$ model enhanced from the original using a machine learning approach. We evaluate the method accuracy with both traditional point-to-point error measures and an event-based validation approach. The results show that predictions using L5 data outperform a 27-day solar wind persistence model in all validation measures but do not achieve a level similar to an L1 monitor. Offsets in timing and the rapidly-changing development of $B_z$ in comparison to $B_x$ and $B_y$ reduce the accuracy. Predictions of $Dst$ from L5 have an RMSE of $9$ nT, which is double the error of $4$ nT using measurements conducted near the Earth. The most useful application of L5 measurements is shown to be in predicting the minimum $Dst$ for the next four days. This method is being implemented in a real-time forecast setting using STEREO-A as an L5 proxy, and has implications for the usefulness of future L5 missions. △ Less

Submitted 1 May, 2020; originally announced May 2020.

Comments: 20 pages, 6 figures, in press at AGU Space Weather

Measurement of the 3He Spin-Structure Functions and of Neutron (3He) Spin-Dependent Sum Rules at 0.035<Q^2<0.24 GeV^2

Authors: V. Sulkosky, J. T. Singh, C. Peng, J. -P. Chen, A. Deur, S. Abrahamyan, K. A. Aniol, D. S. Armstrong, T. Averett, S. L. Bailey, A. Beck, P. Bertin, F. Butaru, W. Boeglin, A. Camsonne, G. D. Cates, C. C. Chang, Seonho Choi, E. Chudakov, L. Coman, J. C Cornejo, B. Craver, F. Cusanno, R. De Leo, C. W. de Jager , et al. (84 additional authors not shown)

Abstract: The spin-structure functions $g_1$ and $g_2$, and the spin-dependent partial cross-section $σ_\mathrm{TT}$ have been extracted from the polarized cross-sections differences, $Δσ_{\parallel}\hspace{-0.06cm}\left(ν,Q^{2}\right)$ and $Δσ_{\perp}\hspace{-0.06cm}\left(ν,Q^{2}\right)$ measured for the $\vec{^\textrm{3}\textrm{He}}(\vec{\textrm{e}},\textrm{e}')\textrm{X}$ reaction, in the E97-110 experim… ▽ More The spin-structure functions $g_1$ and $g_2$, and the spin-dependent partial cross-section $σ_\mathrm{TT}$ have been extracted from the polarized cross-sections differences, $Δσ_{\parallel}\hspace{-0.06cm}\left(ν,Q^{2}\right)$ and $Δσ_{\perp}\hspace{-0.06cm}\left(ν,Q^{2}\right)$ measured for the $\vec{^\textrm{3}\textrm{He}}(\vec{\textrm{e}},\textrm{e}')\textrm{X}$ reaction, in the E97-110 experiment at Jefferson Lab. Polarized electrons with energies from 1.147 to 4.404 GeV were scattered at angles of 6$^{\circ}$ and 9$^{\circ}$ from a longitudinally or transversely polarized $^{3}$He target. The data cover the kinematic regions of the quasi-elastic, resonance production and beyond. From the extracted spin-structure functions, the first moments $\overline{Γ_1}\hspace{-0.06cm}\left(Q^{2}\right)$, $Γ_2\hspace{-0.06cm}\left(Q^{2}\right)$ and $I_{\mathrm{TT}}\hspace{-0.06cm}\left(Q^{2}\right)$ are evaluated with high precision for the neutron in the $Q^2$ range from 0.035 to 0.24~GeV$^{2}$. The comparison of the data and the chiral effective field theory predictions reveals the importance of proper treatment of the $Δ$ degree of freedom for spin observables. △ Less

Submitted 23 April, 2020; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 7 pages, 5 figures, 3 tables. Version published in Phys. Lett. B

Report number: JLAB-PHY-19-3015

Journal ref: Phys.Lett.B 805 135428 (2020)

Using Nab to determine correlations in unpolarized neutron decay

Authors: L. J. Broussard, S. Baeßler, T. L. Bailey, N. Birge, J. D. Bowman, C. B. Crawford, C. Cude-Woods, D. E. Fellers, N. Fomin, E. Frlež, M. T. W. Gericke, L. Hayen, A. P. Jezghani, H. Li, N. Macsai, M. F. Makela, R. R. Mammei, D. Mathews, P. L. McGaughey, P. E. Mueller, D. Počanić, C. A. Royse, A. Salas-Bacci, S. K. L. Sjue, J. C. Ramsey , et al. (6 additional authors not shown)

Abstract: The Nab experiment will measure the ratio of the weak axial-vector and vector coupling constants $λ=g_A/g_V$ with precision $δλ/λ\sim3\times10^{-4}$ and search for a Fierz term $b_F$ at a level $Δb_F<10^{-3}$. The Nab detection system uses thick, large area, segmented silicon detectors to very precisely determine the decay proton's time of flight and the decay electron's energy in coincidence and… ▽ More The Nab experiment will measure the ratio of the weak axial-vector and vector coupling constants $λ=g_A/g_V$ with precision $δλ/λ\sim3\times10^{-4}$ and search for a Fierz term $b_F$ at a level $Δb_F<10^{-3}$. The Nab detection system uses thick, large area, segmented silicon detectors to very precisely determine the decay proton's time of flight and the decay electron's energy in coincidence and reconstruct the correlation between the antineutrino and electron momenta. Excellent understanding of systematic effects affecting timing and energy reconstruction using this detection system are required. To explore these effects, a series of ex situ studies have been undertaken, including a search for a Fierz term at a less sensitive level of $Δb_F<10^{-2}$ in the beta decay of $^{45}$Ca using the UCNA spectrometer. △ Less

Submitted 19 December, 2018; originally announced December 2018.

Comments: Proceedings of the 7th International Syposium on Symmetries in Subatomic Physics SSP2018, Aachen (Germany), 10 - 15 Jun 2018. This is a pre-print of an article published in Hyperfine Interactions. The final authenticated version is available online at: https://doi.org/10.1007/s10751-018-1538-7

Journal ref: Hyperfine Interact. 240 (2019) no.1, 1

Solid deuterium surface degradation at ultracold neutron sources

Authors: A. Anghel, T. L. Bailey, G. Bison, B. Blau, L. J. Broussard, S. M. Clayton, C. Cude-Woods, M. Daum, A. Hawari, N. Hild, P. Huffman, T. M. Ito, K. Kirch, E. Korobkina, B. Lauss, K. Leung, E. M. Lutz, M. Makela, G. Medlin, C. L. Morris, R. W. Pattie, D. Ries, A. Saunders, P. Schmidt-Wellenburg, V. Talanov , et al. (5 additional authors not shown)

Abstract: Solid deuterium (sD_2) is used as an efficient converter to produce ultracold neutrons (UCN). It is known that the sD_2 must be sufficiently cold, of high purity and mostly in its ortho-state in order to guarantee long lifetimes of UCN in the solid from which they are extracted into vacuum. Also the UCN transparency of the bulk sD_2 material must be high because crystal inhomogeneities limit the m… ▽ More Solid deuterium (sD_2) is used as an efficient converter to produce ultracold neutrons (UCN). It is known that the sD_2 must be sufficiently cold, of high purity and mostly in its ortho-state in order to guarantee long lifetimes of UCN in the solid from which they are extracted into vacuum. Also the UCN transparency of the bulk sD_2 material must be high because crystal inhomogeneities limit the mean free path for elastic scattering and reduce the extraction efficiency. Observations at the UCN sources at Paul Scherrer Institute and at Los Alamos National Laboratory consistently show a decrease of the UCN yield with time of operation after initial preparation or later treatment (`conditioning') of the sD_2. We show that, in addition to the quality of the bulk sD_2, the quality of its surface is essential. Our observations and simulations support the view that the surface is deteriorating due to a build-up of D_2 frost-layers under pulsed operation which leads to strong albedo reflections of UCN and subsequent loss. We report results of UCN yield measurements, temperature and pressure behavior of deuterium during source operation and conditioning, and UCN transport simulations. This, together with optical observations of sD_2 frost formation on initially transparent sD_2 in offline studies with pulsed heat input at the North Carolina State University UCN source results in a consistent description of the UCN yield decrease. △ Less

Submitted 28 August, 2018; v1 submitted 23 April, 2018; originally announced April 2018.

Comments: 15 pages, 22 figures, accepted by EPJ-A

Probit models for capture-recapture data subject to imperfect detection, individual heterogeneity and misidentification

Authors: Brett T. McClintock, Larissa L. Bailey, Brian P. Dreher, William A. Link

Abstract: As noninvasive sampling techniques for animal populations have become more popular, there has been increasing interest in the development of capture-recapture models that can accommodate both imperfect detection and misidentification of individuals (e.g., due to genotyping error). However, current methods do not allow for individual variation in parameters, such as detection or survival probabilit… ▽ More As noninvasive sampling techniques for animal populations have become more popular, there has been increasing interest in the development of capture-recapture models that can accommodate both imperfect detection and misidentification of individuals (e.g., due to genotyping error). However, current methods do not allow for individual variation in parameters, such as detection or survival probability. Here we develop misidentification models for capture-recapture data that can simultaneously account for temporal variation, behavioral effects and individual heterogeneity in parameters. To facilitate Bayesian inference using our approach, we extend standard probit regression techniques to latent multinomial models where the dimension and zeros of the response cannot be observed. We also present a novel Metropolis-Hastings within Gibbs algorithm for fitting these models using Markov chain Monte Carlo. Using closed population abundance models for illustration, we re-visit a DNA capture-recapture population study of black bears in Michigan, USA and find evidence of misidentification due to genotyping error, as well as temporal, behavioral and individual variation in detection probability. We also estimate a salamander population of known size from laboratory experiments evaluating the effectiveness of a marking technique commonly used for amphibians and fish. Our model was able to reliably estimate the size of this population and provided evidence of individual heterogeneity in misidentification probability that is attributable to variable mark quality. Our approach is more computationally demanding than previously proposed methods, but it provides the flexibility necessary for a much broader suite of models to be explored while properly accounting for uncertainty introduced by misidentification and imperfect detection. In the absence of misidentification, our probit formulation also provides a convenient and efficient Gibbs sampler for Bayesian analysis of traditional closed population capture-recapture data. △ Less

Submitted 3 February, 2015; v1 submitted 14 January, 2014; originally announced January 2014.

Comments: Published in at http://dx.doi.org/10.1214/14-AOAS783 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org)

Report number: IMS-AOAS-AOAS783

Journal ref: Annals of Applied Statistics 2014, Vol. 8, No. 4, 2461-2484

Ionisation in atmospheres of brown dwarfs and extrasolar planets VI: Properties of large-scale discharge events

Authors: R. L. Bailey, Ch. Helling, G. Hodosán, C. Bilger, C. R. Stark

Abstract: Mineral clouds in substellar atmospheres play a special role as a catalyst for a variety of charge processes. If clouds are charged, the surrounding environment becomes electrically activated, and ensembles of charged grains are electrically discharging (e.g. by lightning), which significantly infuences the local chemistry creating conditions similar to those thought responsible for life in early… ▽ More Mineral clouds in substellar atmospheres play a special role as a catalyst for a variety of charge processes. If clouds are charged, the surrounding environment becomes electrically activated, and ensembles of charged grains are electrically discharging (e.g. by lightning), which significantly infuences the local chemistry creating conditions similar to those thought responsible for life in early planetary atmospheres. We note that such lightning discharges contribute also to the ionisation state of the atmosphere. We apply scaling laws for electrical discharge processes from laboratory measurements and numerical experiments to Drift-Phoenix model atmosphere results to model the discharge's propagation downwards (as lightning) and upwards (as sprites) through the atmospheric clouds. We evaluate the spatial extent and energetics of lightning discharges. The atmospheric volume affected (e.g. by increase of temperature or electron number) is larger in a brown dwarf atmosphere ($10^8 -~10^{10}$m$^3$) than in a giant gas planet's ($10^4 -~10^{6}$m$^3$). Our results suggest that the total dissipated energy in one event is $<10^{12}$ J for all models of initial solar metallicity. First attempts to show the infuence of lightning on the local gas phase indicate an increase of small carbohydrate molecules like CH and CH$_2$ at the expense of CO and CH$_4$. Dust forming molecules are destroyed and the cloud particle properties are frozen-in unless enough time is available for complete evaporation. We summarise instruments potentially suitable to observe lightning on extrasolar objects. △ Less

Submitted 7 January, 2014; v1 submitted 24 December, 2013; originally announced December 2013.

First Measurement of the Neutral Current Excitation of the Delta Resonance on a Proton Target

Authors: G0 Collaboration, D. Androic, D. S. Armstrong, J. Arvieux, S. L. Bailey, D. H. Beck, E. J. Beise, J. Benesch, F. Benmokhtar, L. Bimbot, J. Birchall, P. Bosted, H. Breuer, C. L. Capuano, Y. -C. Chao, A. Coppens, C. A. Davis, C. Ellis, G. Flores, G. Franklin, C. Furget, D. Gaskell, J. Grames, M. T. W. Gericke, G. Guillard , et al. (43 additional authors not shown)

Abstract: The parity-violating asymmetry arising from inelastic electron-nucleon scattering at backward angle (~95 degrees) near the Delta(1232) resonance has been measured using a hydrogen target. From this asymmetry, we extracted the axial transition form factor G^A_{NΔ}, a function of the axial Adler form factors C^A_i. Though G^A_{NΔ} has been previously studied using charged current reactions, this is… ▽ More The parity-violating asymmetry arising from inelastic electron-nucleon scattering at backward angle (~95 degrees) near the Delta(1232) resonance has been measured using a hydrogen target. From this asymmetry, we extracted the axial transition form factor G^A_{NΔ}, a function of the axial Adler form factors C^A_i. Though G^A_{NΔ} has been previously studied using charged current reactions, this is the first measurement of the weak neutral current excitation of the Delta using a proton target. For Q^2 = 0.34 (GeV/c)^2 and W = 1.18 GeV, the asymmetry was measured to be -33.4 \pm (5.3)_{stat} \pm (5.1)_{sys} ppm. The value of G^A_{NΔ} determined from the hydrogen asymmetry was -0.05 \pm (0.35)_{stat} \pm (0.34)_{sys} \pm (0.06)_{theory}. These findings agree within errors with theoretical predictions for both the total asymmetry and the form factor. In addition to the hydrogen measurement, the asymmetry was measured at the same kinematics using a deuterium target. The asymmetry for deuterium was determined to be -43.6 \pm (14.6)_{stat} \pm (6.2)_{sys} ppm. △ Less

Submitted 7 December, 2012; originally announced December 2012.

Comments: 8 pages, 3 figures, submitted to Physics Letters B

New Measurements of the Transverse Beam Asymmetry for Elastic Electron Scattering from Selected Nuclei

Authors: The HAPPEX, PREX Collaborations, :, S. Abrahamyan, A. Acha, A. Afanasev, Z. Ahmed, H. Albataineh, K. Aniol, D. S. Armstrong, W. Armstrong, J. Arrington, T. Averett, B. Babineau, S. L. Bailey, J. Barber, A. Barbieri, A. Beck, V. Bellini, R. Beminiwattha, H. Benaoum, J. Benesch, F. Benmokhtar, P. Bertin, T. Bielarski , et al. (173 additional authors not shown)

Abstract: We have measured the beam-normal single-spin asymmetry $A_n$ in the elastic scattering of 1-3 GeV transversely polarized electrons from $^1$H and for the first time from $^4$He, $^{12}$C, and $^{208}$Pb. For $^1$H, $^4$He and $^{12}$C, the measurements are in agreement with calculations that relate $A_n$ to the imaginary part of the two-photon exchange amplitude including inelastic intermediate st… ▽ More We have measured the beam-normal single-spin asymmetry $A_n$ in the elastic scattering of 1-3 GeV transversely polarized electrons from $^1$H and for the first time from $^4$He, $^{12}$C, and $^{208}$Pb. For $^1$H, $^4$He and $^{12}$C, the measurements are in agreement with calculations that relate $A_n$ to the imaginary part of the two-photon exchange amplitude including inelastic intermediate states. Surprisingly, the $^{208}$Pb result is significantly smaller than the corresponding prediction using the same formalism. These results suggest that a systematic set of new $A_n$ measurements might emerge as a new and sensitive probe of the structure of heavy nuclei. △ Less

Submitted 12 October, 2012; v1 submitted 30 August, 2012; originally announced August 2012.

Comments: 5 pages, 3 figures, accepted by PRL. v3: fixed one author name and affiliation, otherwise no change

Journal ref: Phys. Rev. Lett. 109, (2012) 192501

Measurement of the parity-violating asymmetry in inclusive electroproduction of $π^-$ near the $Δ^0$ resonance

Authors: G0 Collaboration, D. Androić, D. S. Armstrong, J. Arvieux, S. L. Bailey, D. H. Beck, E. J. Beise, J. Benesch, F. Benmokhtar, L. Bimbot, J. Birchall, P. Bosted, H. Breuer, C. L. Capuano, Y. -C. Chao, A. Coppens, C. A. Davis, C. Ellis, G. Flores, G. Franklin, C. Furget, D. Gaskell, M. T. W. Gericke, J. Grames, G. Guillard , et al. (43 additional authors not shown)

Abstract: The parity-violating (PV) asymmetry of inclusive $π^-$ production in electron scattering from a liquid deuterium target was measured at backward angles. The measurement was conducted as a part of the G0 experiment, at a beam energy of 360 MeV. The physics process dominating pion production for these kinematics is quasi-free photoproduction off the neutron via the $Δ^0$ resonance. In the context of… ▽ More The parity-violating (PV) asymmetry of inclusive $π^-$ production in electron scattering from a liquid deuterium target was measured at backward angles. The measurement was conducted as a part of the G0 experiment, at a beam energy of 360 MeV. The physics process dominating pion production for these kinematics is quasi-free photoproduction off the neutron via the $Δ^0$ resonance. In the context of heavy-baryon chiral perturbation theory (HB$χ$PT), this asymmetry is related to a low energy constant $d_Δ^-$ that characterizes the parity-violating $γ$N$Δ$ coupling. Zhu et al. calculated $d_Δ^-$ in a model benchmarked by the large asymmetries seen in hyperon weak radiative decays, and predicted potentially large asymmetries for this process, ranging from $A_γ^-=-5.2$ to $+5.2$ ppm. The measurement performed in this work leads to $A_γ^-=-0.36\pm 1.06\pm 0.37\pm 0.03$ ppm (where sources of statistical, systematic and theoretical uncertainties are included), which would disfavor enchancements considered by Zhu et al. proportional to $V_{ud}/V_{us}$. The measurement is part of a program of inelastic scattering measurements that were conducted by the G0 experiment, seeking to determine the $N-Δ$ axial transition form-factors using PV electron scattering. △ Less

Submitted 2 February, 2012; v1 submitted 7 December, 2011; originally announced December 2011.

Comments: 5 pages, 3 figures

Transverse Beam Spin Asymmetries at Backward Angles in Elastic Electron-Proton and Quasi-elastic Electron-Deuteron Scattering

Authors: The G0 Collaboration, D. Androić, D. S. Armstrong, J. Arvieux, S. L. Bailey, D. H. Beck, E. J. Beise, J. Benesch, F. Benmokhtar, L. Bimbot, J. Birchall, P. Bosted, H. Breuer, C. L. Capuano, Y. -C. Chao, A. Coppens, C. A. Davis, C. Ellis, G. Flores, G. Franklin, C. Furget, D. Gaskell, M. T. W. Gericke, J. Grames, G. Guillard , et al. (45 additional authors not shown)

Abstract: We have measured the beam-normal single-spin asymmetries in elastic scattering of transversely polarized electrons from the proton, and performed the first measurement in quasi-elastic scattering on the deuteron, at backward angles (lab scattering angle of 108 degrees) for Q2 = 0.22 GeV^2/c^2 and 0.63 GeV^2/c^2 at beam energies of 362 MeV and 687 MeV, respectively. The asymmetry arises due to the… ▽ More We have measured the beam-normal single-spin asymmetries in elastic scattering of transversely polarized electrons from the proton, and performed the first measurement in quasi-elastic scattering on the deuteron, at backward angles (lab scattering angle of 108 degrees) for Q2 = 0.22 GeV^2/c^2 and 0.63 GeV^2/c^2 at beam energies of 362 MeV and 687 MeV, respectively. The asymmetry arises due to the imaginary part of the interference of the two-photon exchange amplitude with that of single photon exchange. Results for the proton are consistent with a model calculation which includes inelastic intermediate hadronic (piN) states. An estimate of the beam-normal single-spin asymmetry for the scattering from the neutron is made using a quasi-static deuterium approximation, and is also in agreement with theory. △ Less

Submitted 16 June, 2011; v1 submitted 18 March, 2011; originally announced March 2011.

Journal ref: Phys.Rev.Lett.107:022501,2011

The G0 Experiment: Apparatus for Parity-Violating Electron Scattering Measurements at Forward and Backward Angles

Authors: G0 Collaboration, D. Androic, D. S. Armstrong, J. Arvieux, R. Asaturyan, T. D. Averett, S. L. Bailey, G. Batigne, D. H. Beck, E. J. Beise, J. Benesch, F. Benmokhtar, L. Bimbot, J. Birchall, A. Biselli, P. Bosted, H. Breuer, P. Brindza, C. L. Capuano, R. D. Carlini, R. Carr, N. Chant, Y. -C. Chao, R. Clark, A. Coppens , et al. (105 additional authors not shown)

Abstract: In the G0 experiment, performed at Jefferson Lab, the parity-violating elastic scattering of electrons from protons and quasi-elastic scattering from deuterons is measured in order to determine the neutral weak currents of the nucleon. Asymmetries as small as 1 part per million in the scattering of a polarized electron beam are determined using a dedicated apparatus. It consists of specialized bea… ▽ More In the G0 experiment, performed at Jefferson Lab, the parity-violating elastic scattering of electrons from protons and quasi-elastic scattering from deuterons is measured in order to determine the neutral weak currents of the nucleon. Asymmetries as small as 1 part per million in the scattering of a polarized electron beam are determined using a dedicated apparatus. It consists of specialized beam-monitoring and control systems, a cryogenic hydrogen (or deuterium) target, and a superconducting, toroidal magnetic spectrometer equipped with plastic scintillation and aerogel Cerenkov detectors, as well as fast readout electronics for the measurement of individual events. The overall design and performance of this experimental system is discussed. △ Less

Submitted 3 March, 2011; originally announced March 2011.

Comments: Submitted to Nuclear Instruments and Methods

Journal ref: Nucl. Instrum. Methods A 646 (2011)59

Two Dynamical Classes of Centaurs

Authors: B. L. Bailey, Renu Malhotra

Abstract: The Centaurs are a transient population of small bodies in the outer solar system whose orbits are strongly chaotic. These objects typically suffer significant changes of orbital parameters on timescales of a few thousand years, and their orbital evolution exhibits two types of behaviors described qualitatively as random-walk and resonance-sticking. We have analyzed the chaotic behavior of the k… ▽ More The Centaurs are a transient population of small bodies in the outer solar system whose orbits are strongly chaotic. These objects typically suffer significant changes of orbital parameters on timescales of a few thousand years, and their orbital evolution exhibits two types of behaviors described qualitatively as random-walk and resonance-sticking. We have analyzed the chaotic behavior of the known Centaurs. Our analysis has revealed that the two types of chaotic evolution are quantitatively distinguishable: (1) the random walk-type behavior is well described by so-called generalized diffusion in which the rms deviation of the semimajor axis grows with time t as ~t^H, with Hurst exponent H in the range 0.22--0.95, however (2) orbital evolution dominated by intermittent resonance sticking, with sudden jumps from one mean motion resonance to another, has poorly defined H. We further find that these two types of behavior are correlated with Centaur dynamical lifetime: most Centaurs whose dynamical lifetime is less than ~22 Myr exhibit generalized diffusion, whereas most Centaurs of longer dynamical lifetimes exhibit intermittent resonance sticking. We also find that Centaurs in the diffusing class are likely to evolve into Jupiter-family comets during their dynamical lifetimes, while those in the resonance-hopping class do not. △ Less

Submitted 25 June, 2009; originally announced June 2009.

Comments: 31 pages, including 12 figures and 2 tables. Accepted for publication in Icarus