Showing 1–50 of 165 results for author: Robertson, D

Tilt-to-length coupling in LISA Pathfinder: long-term stability

Authors: M Armano, H Audley, J Baird, P Binetruy, M Born, D Bortoluzzi, E Castelli, A Cavalleri, A Cesarini, A M Cruise, K Danzmann, M de Deus Silva, I Diepholz, G Dixon, R Dolesi, L Ferraioli, V Ferroni, E D Fitzsimons, M Freschi, L Gesa, D Giardini, F Gibert, R Giusteri, C Grimani, J Grzymisch , et al. (53 additional authors not shown)

Abstract: The tilt-to-length coupling during the LISA Pathfinder mission has been numerically and analytically modeled for particular timespans. In this work, we investigate the long-term stability of the coupling coefficients of this noise. We show that they drifted slowly (by 1\,$μみゅー$m/rad and 6$\times10^{-6}$ in 100 days) and strongly correlated to temperature changes within the satellite (8\,$μみゅー$m/rad/K an… ▽ More The tilt-to-length coupling during the LISA Pathfinder mission has been numerically and analytically modeled for particular timespans. In this work, we investigate the long-term stability of the coupling coefficients of this noise. We show that they drifted slowly (by 1\,$μみゅー$m/rad and 6$\times10^{-6}$ in 100 days) and strongly correlated to temperature changes within the satellite (8\,$μみゅー$m/rad/K and 30$\times10^{-6}$/K). Based on analytical TTL coupling models, we attribute the temperature-driven coupling changes to rotations of the test masses and small distortions in the optical setup. Particularly, we show that LISA Pathfinder's optical baseplate was bent during the cooldown experiment, which started in late 2016 and lasted several months. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Precision measurements of the magnetic parameters of LISA Pathfinder test masses

Authors: M Armano, H Audley, J Baird, P Binetruy, M Born, D Bortoluzzi, E Castelli, A Cavalleri, A Cesarini, A M Cruise, K Danzmann, M De Deus Silva, I Diepholz, G Dixon, R Dolesi, L Ferraioli, V Ferroni, E D Fitzsimons, M Freschi, L Gesa, D Giardini, F Gibert, R Giusteri, C Grimani, J Grzymisch , et al. (53 additional authors not shown)

Abstract: A precise characterization of the magnetic properties of LISA Pathfinder free falling test-masses is of special interest for future gravitational wave observatory in space. Magnetic forces have an important impact on the instrument sensitivity in the low frequency regime below the millihertz. In this paper we report on the magnetic injection experiments performed throughout LISA Pathfinder operati… ▽ More A precise characterization of the magnetic properties of LISA Pathfinder free falling test-masses is of special interest for future gravitational wave observatory in space. Magnetic forces have an important impact on the instrument sensitivity in the low frequency regime below the millihertz. In this paper we report on the magnetic injection experiments performed throughout LISA Pathfinder operations. We show how these experiments allowed a high precision estimate of the instrument magnetic parameters. The remanent magnetic moment was found to have a modulus of $(0.245\pm0.081)\,\rm{nAm}^2$, the x-component of the background magnetic field within the test masses position was measured to be $(414 \pm 74)$ nT and its gradient had a value of $(-7.4\pm 2.1)\,μみゅー$T/m. Finally, we also measured the test mass magnetic susceptibility to be $(-3.35\pm0.15)\times$10$^{-5}$ in the low frequency regime. All results are in agreement with on-ground estimates. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Magnetic-induced force noise in LISA Pathfinder free-falling test masses

Authors: M Armano, H Audley, J Baird, P Binetruy, M Born, D Bortoluzzi, E Castelli, A Cavalleri, A Cesarini, A M Cruise, K Danzmann, M De Deus Silva, I Diepholz, G Dixon, R Dolesi, L Ferraioli, V Ferroni, E D Fitzsimons, M Freschi, L Gesa, D Giardini, F Gibert, R Giusteri, C Grimani, J Grzymisch , et al. (53 additional authors not shown)

Abstract: LISA Pathfinder was a mission designed to test key technologies required for gravitational wave detection in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime, which corresponds to the measurement band of interest for future space-borne gravitational wave observatories. Magnetic-induced forces couple to the test mass motion, introducing a c… ▽ More LISA Pathfinder was a mission designed to test key technologies required for gravitational wave detection in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime, which corresponds to the measurement band of interest for future space-borne gravitational wave observatories. Magnetic-induced forces couple to the test mass motion, introducing a contribution to the relative acceleration noise between the free falling test masses. In this Letter we present the first complete estimate of this term of the instrument performance model. Our results set the magnetic-induced acceleration noise during the February 2017 noise run of $\rm 0.25_{-0.08}^{+0.15}\,fm\,s^{-2}/\sqrt{Hz}$ at 1 mHz and $\rm 1.01_{-0.24}^{+0.73}\, fm\,s^{-2}/\sqrt{Hz}$ at 0.1 mHz. We also discuss how the non-stationarities of the interplanetary magnetic field can affect these values during extreme space weather conditions. △ Less

Submitted 5 July, 2024; originally announced July 2024.

The Blue Multi Unit Spectroscopic Explorer (BlueMUSE) on the VLT: science drivers and overview of instrument design

Authors: Johan Richard, Rémi Giroud, Florence Laurent, Davor Krajnović, Alexandre Jeanneau, Roland Bacon, Manuel Abreu, Angela Adamo, Ricardo Araujo, Nicolas Bouché, Jarle Brinchmann, Zhemin Cai, Norberto Castro, Ariadna Calcines, Diane Chapuis, Adélaïde Claeyssens, Luca Cortese, Emanuele Daddi, Christopher Davison, Michael Goodwin, Robert Harris, Matthew Hayes, Mathilde Jauzac, Andreas Kelz, Jean-Paul Kneib , et al. (25 additional authors not shown)

Abstract: BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its… ▽ More BlueMUSE is a blue-optimised, medium spectral resolution, panoramic integral field spectrograph under development for the Very Large Telescope (VLT). With an optimised transmission down to 350 nm, spectral resolution of R$\sim$3500 on average across the wavelength range, and a large FoV (1 arcmin$^2$), BlueMUSE will open up a new range of galactic and extragalactic science cases facilitated by its specific capabilities. The BlueMUSE consortium includes 9 institutes located in 7 countries and is led by the Centre de Recherche Astrophysique de Lyon (CRAL). The BlueMUSE project development is currently in Phase A, with an expected first light at the VLT in 2031. We introduce here the Top Level Requirements (TLRs) derived from the main science cases, and then present an overview of the BlueMUSE system and its subsystems fulfilling these TLRs. We specifically emphasize the tradeoffs that are made and the key distinctions compared to the MUSE instrument, upon which the system architecture is built. △ Less

Submitted 28 August, 2024; v1 submitted 19 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, proceedings of the SPIE astronomical telescopes and instrumentation conference, Yokohama, 16-21 June

In-depth analysis of LISA Pathfinder performance results: Time evolution, noise projection, physical models, and implications for LISA

Authors: M. Armano, H. Audley, J. Baird, P. Binetruy, M. Born, D. Bortoluzzi, E. Castelli, A. Cavalleri, A. Cesarini, V. Chiavegato, A. M. Cruise, D. Dal Bosco, K. Danzmann, M. De Deus Silva, I. Diepholz, G. Dixon, R. Dolesi, L. Ferraioli, V. Ferroni, E. D. Fitzsimons, M. Freschi, L. Gesa, D. Giardini, F. Gibert, R. Giusteri , et al. (55 additional authors not shown)

Abstract: We present an in-depth analysis of the LISA Pathfinder differential acceleration performance over the entire course of its science operations, spanning approximately 500 days. We find that: 1) the evolution of the Brownian noise that dominates the acceleration amplitude spectral density (ASD), for frequencies $f\gtrsim 1\,\text{mHz}$, is consistent with the decaying pressure due to the outgassing… ▽ More We present an in-depth analysis of the LISA Pathfinder differential acceleration performance over the entire course of its science operations, spanning approximately 500 days. We find that: 1) the evolution of the Brownian noise that dominates the acceleration amplitude spectral density (ASD), for frequencies $f\gtrsim 1\,\text{mHz}$, is consistent with the decaying pressure due to the outgassing of a single gaseous species. 2) between $f=36\,μみゅー\text{Hz}$ and $1\,\text{mHz}$, the acceleration ASD shows a $1/f$ tail in excess of the Brownian noise of almost constant amplitude, with $\simeq 20\%$ fluctuations over a period of a few days, with no particular time pattern over the course of the mission; 3) at the lowest considered frequency of $f=18\,μみゅー\text{Hz}$, the ASD significantly deviates from the $1/f$ behavior, because of temperature fluctuations that appear to modulate a quasi-static pressure gradient, sustained by the asymmetries of the outgassing pattern. We also present the results of a projection of the observed acceleration noise on the potential sources for which we had either a direct correlation measurement, or a quantitative estimate from dedicated experiments. These sources account for approximately $40\%$ of the noise power in the $1/f$ tail. Finally, we analyze the possible sources of the remaining unexplained fraction, and identify the possible measures that may be taken to keep those under control in LISA. △ Less

Submitted 4 September, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Journal ref: Phys. Rev. D 110, 042004 (2024)

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

Cyborgs for strategic communication on social media

Authors: Lynnette Hui Xian Ng, Dawn C. Robertson, Kathleen M. Carley

Abstract: Social media platforms are a key ground of information consumption and dissemination. Key figures like politicians, celebrities and activists have leveraged on its wide user base for strategic communication. Strategic communications, or StratCom, is the deliberate act of information creation and distribution. Its techniques are used by these key figures for establishing their brand and amplifying… ▽ More Social media platforms are a key ground of information consumption and dissemination. Key figures like politicians, celebrities and activists have leveraged on its wide user base for strategic communication. Strategic communications, or StratCom, is the deliberate act of information creation and distribution. Its techniques are used by these key figures for establishing their brand and amplifying their messages. Automated scripts are used on top of personal touches to quickly and effectively perform these tasks. The combination of automation and manual online posting creates a Cyborg social media profile, which is a hybrid between bot and human. In this study, we establish a quantitative definition for a Cyborg account, which is an account that are detected as bots in one time window, and identified as humans in another. This definition makes use of frequent changes of bot classification labels and large differences in bot likelihood scores to identify Cyborgs. We perform a large-scale analysis across over 3.1 million users from Twitter collected from two key events, the 2020 Coronavirus pandemic and 2020 US Elections. We extract Cyborgs from two datasets and employ tools from network science, natural language processing and manual annotation to characterize Cyborg accounts. Our analyses identify Cyborg accounts are mostly constructed for strategic communication uses, have a strong duality in their bot/human classification and are tactically positioned in the social media network, aiding these accounts to promote their desired content. Cyborgs are also discovered to have long online lives, indicating their ability to evade bot detectors, or the graciousness of platforms to allow their operations. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: To appear in Big Data and Society

NanoNewton electrostatic force actuators for femtoNewton-sensitive measurements: system performance test in the LISA Pathfinder mission

Authors: M Armano, H Audley, J Baird, M Bassan, P Binetruy, M Born, D Bortoluzzi, E Castelli, A Cavalleri, A Cesarini, V Chiavegato, A M Cruise, D Dal Bosco, K Danzmann, M De Deus Silva, R De Rosa, L Di Fiore, I Diepholz, G Dixon, R Dolesi, L Ferraioli V Ferroni, E D Fitzsimons, M Freschi, L Gesa, D Giardini , et al. (65 additional authors not shown)

Abstract: Electrostatic force actuation is a key component of the system of geodesic reference test masses (TM) for the LISA orbiting gravitational wave observatory and in particular for performance at low frequencies, below 1 mHz, where the observatory sensitivity is limited by stray force noise. The system needs to apply forces of order 10$^{-9}$ N while limiting fluctuations in the measurement band to le… ▽ More Electrostatic force actuation is a key component of the system of geodesic reference test masses (TM) for the LISA orbiting gravitational wave observatory and in particular for performance at low frequencies, below 1 mHz, where the observatory sensitivity is limited by stray force noise. The system needs to apply forces of order 10$^{-9}$ N while limiting fluctuations in the measurement band to levels approaching 10$^{-15}$ N/Hz$^{1/2}$. We present here the LISA actuation system design, based on audio-frequency voltage carrier signals, and results of its in-flight performance test with the LISA Pathfinder test mission. In LISA, TM force actuation is used to align the otherwise free-falling TM to the spacecraft-mounted optical metrology system, without any forcing along the critical gravitational wave-sensitive interferometry axes. In LISA Pathfinder, on the other hand, the actuation was used also to stabilize the TM along the critical $x$ axis joining the two TM, with the commanded actuation force entering directly into the mission's main differential acceleration science observable. The mission allowed demonstration of the full compatibility of the electrostatic actuation system with the LISA observatory requirements, including dedicated measurement campaigns to amplify, isolate, and quantify the two main force noise contributions from the actuation system, from actuator gain noise and from low frequency ``in band'' voltage fluctuations. These campaigns have shown actuation force noise to be a relevant, but not dominant, noise source in LISA Pathfinder and have allowed performance projections for the conditions expected in the LISA mission. △ Less

Submitted 30 December, 2023; originally announced January 2024.

Problems on infinite sumset configurations in the integers and beyond

Authors: Bryna Kra, Joel Moreira, Florian K. Richter, Donald Robertson

Abstract: In contrast to finite arithmetic configurations, relatively little is known about which infinite patterns can be found in every set of natural numbers with positive density. Building on recent advances showing infinite sumsets can be found, we explore numerous open problems and obstructions to finding other infinite configurations in every set of natural numbers with positive density. In contrast to finite arithmetic configurations, relatively little is known about which infinite patterns can be found in every set of natural numbers with positive density. Building on recent advances showing infinite sumsets can be found, we explore numerous open problems and obstructions to finding other infinite configurations in every set of natural numbers with positive density. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: 37 pages

MSC Class: 05D10; 11B13; 11B30; 05-02; 11-02; 00A27; 37A05

The use of restricted mean survival time to estimate treatment effect under model misspecification, a simulation study

Authors: Emily Alger, David S. Robertson, Abigail J. Burdon

Abstract: The use of the non-parametric Restricted Mean Survival Time endpoint (RMST) has grown in popularity as trialists look to analyse time-to-event outcomes without the restrictions of the proportional hazards assumption. In this paper, we evaluate the power and type I error rate of the parametric and non-parametric RMST estimators when treatment effect is explained by multiple covariates, including an… ▽ More The use of the non-parametric Restricted Mean Survival Time endpoint (RMST) has grown in popularity as trialists look to analyse time-to-event outcomes without the restrictions of the proportional hazards assumption. In this paper, we evaluate the power and type I error rate of the parametric and non-parametric RMST estimators when treatment effect is explained by multiple covariates, including an interaction term. Utilising the RMST estimator in this way allows the combined treatment effect to be summarised as a one-dimensional estimator, which is evaluated using a one-sided hypothesis Z-test. The estimators are either fully specified or misspecified, both in terms of unaccounted covariates or misspecified knot points (where trials exhibit crossing survival curves). A placebo-controlled trial of Gamma interferon is used as a motivating example to simulate associated survival times. When correctly specified, the parametric RMST estimator has the greatest power, regardless of the time of analysis. The misspecified RMST estimator generally performs similarly when covariates mirror those of the fitted case study dataset. However, as the magnitude of the unaccounted covariate increases, the associated power of the estimator decreases. In all cases, the non-parametric RMST estimator has the lowest power, and power remains very reliant on the time of analysis (with a later analysis time correlated with greater power). △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 26 pages, 15 figures, 9 tables

Giant anomalous Hall effect in epitaxial Mn$_{3.2}$Ge films with a cubic kagome structure

Authors: J. S. R. McCoombs, B. D. MacNeil, V. Askarpour, J. Myra, H. Herdin, M. Pula, M. D. Robertson, G. M. Luke, K. L. Kavanagh, J. Maassen, T. L. Monchesky

Abstract: We report on the first example of epitaxial Mn$_{3 + δでるた}$Ge thin films with a cubic $L1_2$ structure. The films are found to exhibit frustrated ferromagnetism with an average magnetization corresponding to 0.98$~\pm~$0.06$~μみゅー_B$/Mn, far larger than the parasitic ferromagnetism in hexagonal Mn$_3$Ge and the partially compensated ferrimagnetism in tetragonal Mn$_3$Ge. The Hall conductivity is the larg… ▽ More We report on the first example of epitaxial Mn$_{3 + δでるた}$Ge thin films with a cubic $L1_2$ structure. The films are found to exhibit frustrated ferromagnetism with an average magnetization corresponding to 0.98$~\pm~$0.06$~μみゅー_B$/Mn, far larger than the parasitic ferromagnetism in hexagonal Mn$_3$Ge and the partially compensated ferrimagnetism in tetragonal Mn$_3$Ge. The Hall conductivity is the largest reported for the kagome magnets with a low temperature value of $σしぐま_{xy} = 1587~$S/cm. Density functional calculations predict that a chiral antiferromagnetic structure is lower in energy than a ferromagnetic configuration in an ordered stoichiometric crystal. However, chemical disorder driven by the excess Mn in our films explains why a frustrated 120$^\circ$ spin structure is not observed. Comparisons between the magnetization and the Hall resistivity indicate that a non-coplanar spin structure contributes the Hall signal. Anisotropic magnetoresistance and planar Hall effect with hysteresis up to 14 T provides further insights into this material. △ Less

Submitted 17 July, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Asteroid 2023 NT1: A Cautionary Tale

Authors: Brin K. Bailey, Alexander N. Cohen, Dharv Patel, Philip Lubin, Mark Boslough, Darrel Robertson, Sasha Egan, Jeeya Khetia, Teagan Costa, Elizabeth Silber, Irina Sagert, Oleg Korobkin, Glenn Sjoden

Abstract: We investigate a variety of short warning time, terminal mitigation scenarios via fragmentation for a hypothetical impact of asteroid 2023 NT1, a Near-Earth Object (NEO) that was discovered on July 15, 2023, two days after its closest approach to Earth on July 13. The asteroid passed by Earth within ~0.25 lunar distances with a closest approach distance of ~10$^{5}$ km and speed of 11.27 km/s. Its… ▽ More We investigate a variety of short warning time, terminal mitigation scenarios via fragmentation for a hypothetical impact of asteroid 2023 NT1, a Near-Earth Object (NEO) that was discovered on July 15, 2023, two days after its closest approach to Earth on July 13. The asteroid passed by Earth within ~0.25 lunar distances with a closest approach distance of ~10$^{5}$ km and speed of 11.27 km/s. Its size remains largely uncertain, with an estimated diameter range of 26 - 58 m and probable diameter estimate (weighted by the NEO size frequency distribution) of 34 m (JPL Sentry, September 12, 2023). The asteroid approached Earth from the direction of the Sun, as did both the Chelyabinsk asteroid in 2013 and comet NEOWISE in 2021. As a result, 2023 NT1 remained undetected until after its closest approach. If it had been on a collision course, it would have had an impact energy of ~1.5 Mt (assuming a spherical asteroid with the probable diameter estimate of 34 m, 2.6 g/cm$^{3}$ uniform density, and impact speed of 15.59 km/s). 2023 NT1 represents a threat that could have caused significant local damage (~3x Chelyabinsk airburst energy). We utilize the PI ("Pulverize It") method for planetary defense to model potential mitigation scenarios of an object like 2023 NT1 through simulations of hypervelocity asteroid disruption and atmospheric ground effects for the case of a terminal defense mode. Simulations suggest that PI is an effective multimodal approach for planetary defense that can operate in extremely short interdiction modes (with intercepts as short as hours prior to impact), in addition to long interdiction time scales with months to years of warning. Our simulations support the proposition that threats like 2023 NT1 can be effectively mitigated with intercepts of one day (or less) prior to impact, yielding minimal to no ground damage, using modest resources and existing technologies. △ Less

Submitted 7 November, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

3D terrain mapping and filtering from coarse resolution data cubes extracted from real-aperture 94 GHz radar

Authors: William D. Harcourt, David G. Macfarlane, Duncan A. Robertson

Abstract: Accurate, high-resolution 3D mapping of environmental terrain is critical in a range of disciplines. In this study, we develop a new technique, called the PCFilt-94 algorithm, to extract 3D point clouds from coarse resolution millimetre-wave radar data cubes and quantify their associated uncertainties. A technique to non-coherently average neighbouring waveforms surrounding each AVTIS2 range profi… ▽ More Accurate, high-resolution 3D mapping of environmental terrain is critical in a range of disciplines. In this study, we develop a new technique, called the PCFilt-94 algorithm, to extract 3D point clouds from coarse resolution millimetre-wave radar data cubes and quantify their associated uncertainties. A technique to non-coherently average neighbouring waveforms surrounding each AVTIS2 range profile was developed in order to reduce speckle and was found to reduce point cloud uncertainty by 13% at long range and 20% at short range. Further, a Voronoi-based point cloud outlier removal algorithm was implemented which iteratively removes outliers in a point cloud until the process converges to the removal of 0 points. Taken together, the new processing methodology produces a stable point cloud, which means that: 1) it is repeatable even when using different point cloud extraction and filtering parameter values during pre-processing, and 2) is less sensitive to over-filtering through the point cloud processing workflow. Using an optimal number of Ground Control Points (GCPs) for georeferencing, which was determined to be 3 at close range (<1.5 km) and 5 at long range (>3 km), point cloud uncertainty was estimated to be approximately 1.5 m at 1.5 km to 3 m at 3 km and followed a Lorentzian distribution. These uncertainties are smaller than those reported for other close-range radar systems used for terrain mapping. The results of this study should be used as a benchmark for future application of millimetre-wave radar systems for 3D terrain mapping. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 19 pages, 22 figures, submitted to IEEE Transactions on Geoscience and Remote Sensing

Tilt-to-length coupling in LISA Pathfinder: a data analysis

Authors: M Armano, H Audley, J Baird, P Binetruy, M Born, D Bortoluzzi, E Castelli, A Cavalleri, A Cesarini, A M Cruise, K Danzmann, M de Deus Silva, I Diepholz, G Dixon, R Dolesi, L Ferraioli, V Ferroni, E D Fitzsimons, M Freschi, L Gesa, D Giardini, F Gibert, R Giusteri, C Grimani, J Grzymisch , et al. (54 additional authors not shown)

Abstract: We present a study of the tilt-to-length coupling noise during the LISA Pathfinder mission and how it depended on the system's alignment. Tilt-to-length coupling noise is the unwanted coupling of angular and lateral spacecraft or test mass motion into the primary interferometric displacement readout. It was one of the major noise sources in the LISA Pathfinder mission and is likewise expected to b… ▽ More We present a study of the tilt-to-length coupling noise during the LISA Pathfinder mission and how it depended on the system's alignment. Tilt-to-length coupling noise is the unwanted coupling of angular and lateral spacecraft or test mass motion into the primary interferometric displacement readout. It was one of the major noise sources in the LISA Pathfinder mission and is likewise expected to be a primary noise source in LISA. We demonstrate here that a recently derived and published analytical model describes the dependency of the LISA Pathfinder tilt-to-length coupling noise on the alignment of the two freely falling test masses. This was verified with the data taken before and after the realignments performed in March (engineering days) and June 2016, and during a two-day experiment in February 2017 (long cross-talk experiment). The latter was performed with the explicit goal of testing the tilt-to-length coupling noise dependency on the test mass alignment. Using the analytical model, we show that all realignments performed during the mission were only partially successful and explain the reasons why. In addition to the analytical model, we computed another physical tilt-to-length coupling model via a minimising routine making use of the long cross-talk experiment data. A similar approach could prove useful for the LISA mission. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Self-similar quantum groups

Authors: Nathan Brownlowe, David Robertson

Abstract: We introduce the notion of self-similarity for compact quantum groups. For a finite set $X$, we introduce a $C^*$-algebra $\mathbb{A}_X$, which is the quantum automorphism group of the infinite homogeneous rooted tree $X^*$. Self-similar quantum groups are then certain quantum subgroups of $\mathbb{A}_X$. Our main class of examples are called finitely-constrained self-similar quantum groups, and w… ▽ More We introduce the notion of self-similarity for compact quantum groups. For a finite set $X$, we introduce a $C^*$-algebra $\mathbb{A}_X$, which is the quantum automorphism group of the infinite homogeneous rooted tree $X^*$. Self-similar quantum groups are then certain quantum subgroups of $\mathbb{A}_X$. Our main class of examples are called finitely-constrained self-similar quantum groups, and we find a class of these examples that can be described as quantum wreath products by subgroups of the quantum permutation group. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Characterization of Fast Ion Transport via Position-Dependent Optical Deshelving

Authors: Craig R. Clark, Creston D. Herold, J. True Merrill, Holly N. Tinkey, Wade Rellergert, Robert Clark, Roger Brown, Wesley D. Robertson, Curtis Volin, Kara Maller, Chris Shappert, Brian J. McMahon, Brian C. Sawyer, Kenton R. Brown

Abstract: Ion transport is an essential operation in some models of quantum information processing, where fast ion shuttling with minimal motional excitation is necessary for efficient, high-fidelity quantum logic. While fast and cold ion shuttling has been demonstrated, the dynamics and specific trajectory of an ion during diabatic transport have not been studied in detail. Here we describe a position-depe… ▽ More Ion transport is an essential operation in some models of quantum information processing, where fast ion shuttling with minimal motional excitation is necessary for efficient, high-fidelity quantum logic. While fast and cold ion shuttling has been demonstrated, the dynamics and specific trajectory of an ion during diabatic transport have not been studied in detail. Here we describe a position-dependent optical deshelving technique useful for sampling an ion's position throughout its trajectory, and we demonstrate the technique on fast linear transport of a $^{40}\text{Ca}^+$ ion in a surface-electrode ion trap. At high speed, the trap's electrode filters strongly distort the transport potential waveform. With this technique, we observe deviations from the intended constant-velocity (100 m/s) transport: we measure an average speed of 83(2) m/s and a peak speed of 251(6) m/s over a distance of 120 $μみゅー$m △ Less

Submitted 28 April, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Needs, trends, and advances in scintillators for radiographic imaging and tomography

Authors: Zhehui Wang, Christophe Dujardin, Matthew S. Freeman, Amanda E. Gehring, James F. Hunter, Paul Lecoq, Wei Liu, Charles L. Melcher, C. L. Morris, Martin Nikl, Ghanshyam Pilania, Reeju Pokharel, Daniel G. Robertson, Daniel J. Rutstrom, Sky K. Sjue, Anton S. Tremsin, S. A. Watson, Brenden W. Wiggins, Nicola M. Winch, Mariya Zhuravleva

Abstract: Scintillators are important materials for radiographic imaging and tomography (RadIT), when ionizing radiations are used to reveal internal structures of materials. Since its invention by Röntgen, RadIT now come in many modalities such as absorption-based X-ray radiography, phase contrast X-ray imaging, coherent X-ray diffractive imaging, high-energy X- and $γがんま-$ray radiography at above 1 MeV, X-ra… ▽ More Scintillators are important materials for radiographic imaging and tomography (RadIT), when ionizing radiations are used to reveal internal structures of materials. Since its invention by Röntgen, RadIT now come in many modalities such as absorption-based X-ray radiography, phase contrast X-ray imaging, coherent X-ray diffractive imaging, high-energy X- and $γがんま-$ray radiography at above 1 MeV, X-ray computed tomography (CT), proton imaging and tomography (IT), neutron IT, positron emission tomography (PET), high-energy electron radiography, muon tomography, etc. Spatial, temporal resolution, sensitivity, and radiation hardness, among others, are common metrics for RadIT performance, which are enabled by, in addition to scintillators, advances in high-luminosity accelerators and high-power lasers, photodetectors especially CMOS pixelated sensor arrays, and lately data science. Medical imaging, nondestructive testing, nuclear safety and safeguards are traditional RadIT applications. Examples of growing or emerging applications include space, additive manufacturing, machine vision, and virtual reality or `metaverse'. Scintillator metrics such as light yield and decay time are correlated to RadIT metrics. More than 160 kinds of scintillators and applications are presented during the SCINT22 conference. New trends include inorganic and organic scintillator heterostructures, liquid phase synthesis of perovskites and $μみゅー$m-thick films, use of multiphysics models and data science to guide scintillator development, structural innovations such as photonic crystals, nanoscintillators enhanced by the Purcell effect, novel scintillator fibers, and multilayer configurations. Opportunities exist through optimization of RadIT with reduced radiation dose, data-driven measurements, photon/particle counting and tracking methods supplementing time-integrated measurements, and multimodal RadIT. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: 45 pages, 43 Figures, SCINT22 conference overview

Report number: Los Alamos report number LA-UR-22-32994

Journal ref: IEEE Transactions on Nuclear Science ( Volume: 70, Issue: 7, July 2023), pp. 1244 - 1280

Point estimation for adaptive trial designs II: practical considerations and guidance

Authors: David S. Robertson, Babak Choodari-Oskooei, Munya Dimairo, Laura Flight, Philip Pallmann, Thomas Jaki

Abstract: In adaptive clinical trials, the conventional end-of-trial point estimate of a treatment effect is prone to bias, that is, a systematic tendency to deviate from its true value. As stated in recent FDA guidance on adaptive designs, it is desirable to report estimates of treatment effects that reduce or remove this bias. However, it may be unclear which of the available estimators are preferable, an… ▽ More In adaptive clinical trials, the conventional end-of-trial point estimate of a treatment effect is prone to bias, that is, a systematic tendency to deviate from its true value. As stated in recent FDA guidance on adaptive designs, it is desirable to report estimates of treatment effects that reduce or remove this bias. However, it may be unclear which of the available estimators are preferable, and their use remains rare in practice. This paper is the second in a two-part series that studies the issue of bias in point estimation for adaptive trials. Part I provided a methodological review of approaches to remove or reduce the potential bias in point estimation for adaptive designs. In part II, we discuss how bias can affect standard estimators and assess the negative impact this can have. We review current practice for reporting point estimates and illustrate the computation of different estimators using a real adaptive trial example (including code), which we use as a basis for a simulation study. We show that while on average the values of these estimators can be similar, for a particular trial realisation they can give noticeably different values for the estimated treatment effect. Finally, we propose guidelines for researchers around the choice of estimators and the reporting of estimates following an adaptive design. The issue of bias should be considered throughout the whole lifecycle of an adaptive design, with the estimation strategy pre-specified in the statistical analysis plan. When available, unbiased or bias-reduced estimates are to be preferred. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: Part II of a two-part series, part I can be found at arXiv:2105.08836

MSC Class: 62F10

Charging of free-falling test masses in orbit due to cosmic rays: results from LISA Pathfinder

Authors: LISA Pathfinder Collaboration, M. Armano, H. Audley, J. Baird, P. Binetruy, M. Born, D. Bortoluzzi, E. Castelli, A. Cavalleri A. Cesarini, A. M Cruise, K. Danzmann, M. de Deus Silva, I. Diepholz, G. Dixon, R. Dolesi, L. Ferraioli, V. Ferroni, E. D. Fitzsimons, M. Freschi, L. Gesa, D. Giardini, F. Gibert, R. Giusteri, C. Grimani, J. Grzymisch , et al. (50 additional authors not shown)

Abstract: A comprehensive summary of the measurements made to characterize test mass charging due to the space environment during the LISA Pathfinder mission is presented. Measurements of the residual charge of the test mass after release by the grabbing and positioning mechanism, show that the initial charge of the test masses was negative after all releases, leaving the test mass with a potential in the r… ▽ More A comprehensive summary of the measurements made to characterize test mass charging due to the space environment during the LISA Pathfinder mission is presented. Measurements of the residual charge of the test mass after release by the grabbing and positioning mechanism, show that the initial charge of the test masses was negative after all releases, leaving the test mass with a potential in the range $-12$ mV to $-512$ mV. Variations in the neutral test mass charging rate between $21.7$ e s$^{-1}$ and $30.7$ e s$^{-1}$ were observed over the course of the 17-month science operations produced by cosmic ray flux changes including a Forbush decrease associated with a small solar energetic particle event. A dependence of the cosmic ray charging rate on the test mass potential between $-30.2$ e s$^{-1}$ V$^{-1}$ and $-40.3$ e s$^{-1}$ V$^{-1}$ was observed and this is attributed to a contribution to charging from low-energy electrons emitted from the gold surfaces of the gravitational reference sensor. Data from the on-board particle detector show a reliable correlation with the charging rate and with other environmental monitors of the cosmic ray flux. This correlation is exploited to extrapolate test mass charging rates to a 20-year period giving useful insight into the expected range of charging rate that may be observed in the LISA mission. △ Less

Submitted 23 March, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: 17 pages, 10 figures

Journal ref: Phys. Rev. D 107, 062007 (2023)

Galaxy And Mass Assembly: Galaxy Morphology in the Green Valley, Prominent rings and looser Spiral Arms

Authors: Dominic Smith, Lutz Haberzettl, L. E. Porter, Ren Porter-Temple, Christopher P. A. Henry, Benne Holwerda, A. R. Lopez-Sanchez, Steven Phillipps, Alister W. Graham, Sarah Brough, Kevin A. Pimbblet, Jochen Liske, Lee S. Kelvin, Clayton D. Robertson, Wade Roemer, Michael Walmsley, David O'Ryan, Tobias Geron

Abstract: Galaxies broadly fall into two categories: star-forming (blue) galaxies and quiescent (red) galaxies. In between, one finds the less populated ``green valley". Some of these galaxies are suspected to be in the process of ceasing their star-formation through a gradual exhaustion of gas supply or already dead and are experiencing a rejuvenation of star-formation through fuel injection. We use the Ga… ▽ More Galaxies broadly fall into two categories: star-forming (blue) galaxies and quiescent (red) galaxies. In between, one finds the less populated ``green valley". Some of these galaxies are suspected to be in the process of ceasing their star-formation through a gradual exhaustion of gas supply or already dead and are experiencing a rejuvenation of star-formation through fuel injection. We use the Galaxy And Mass Assembly database and the Galaxy Zoo citizen science morphological estimates to compare the morphology of galaxies in the green valley against those in the red sequence and blue cloud. Our goal is to examine the structural differences within galaxies that fall in the green valley, and what brings them there. Previous results found disc features such as rings and lenses are more prominently represented in the green valley population. We revisit this with a similar sized data set of galaxies with morphology labels provided by the Galaxy Zoo for the GAMA fields based on new KiDS images. Our aim is to compare qualitatively the results from expert classification to that of citizen science. We observe that ring structures are indeed found more commonly in green valley galaxies compared to their red and blue counterparts. We suggest that ring structures are a consequence of disc galaxies in the green valley actively exhibiting characteristics of fading discs and evolving disc morphology of galaxies. We note that the progression from blue to red correlates with loosening spiral arm structure. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 11 pages, 21 figures, accepted to MNRAS

Snowmass 2021 Topical Report on Synergies in Research at Underground Facilities

Authors: Catalina Curceanu, Derek Elsworth, Joe Formaggio, Jan Harms, Daniel Robertson, William Roggenthen, Herb Wang, Sijbrand de Jong, John L. Orrell

Abstract: This is a Snowmass 2021 Topical Report for the Underground Facilities and Infrastructure Frontier on Synergies in Research at Underground Facilities: A broad range of scientific and engineering research is possible in underground laboratories, beyond the physics-focused activities described in the other Underground Facilities and Infrastructure Topical Reports. These areas of research include nucl… ▽ More This is a Snowmass 2021 Topical Report for the Underground Facilities and Infrastructure Frontier on Synergies in Research at Underground Facilities: A broad range of scientific and engineering research is possible in underground laboratories, beyond the physics-focused activities described in the other Underground Facilities and Infrastructure Topical Reports. These areas of research include nuclear astrophysics, geology, geoengineering, gravitational wave detection, biology, and perhaps soon quantum information science. This UF Topical Report will survey those other scientific and engineering research activities that share interest in research-orientated Underground Facilities and Infrastructure. In most cases the breadth and depth of research aims is too large to cover in completeness and references to surveys or key documents for those fields are provided after introductory summaries. Additional attention is then given to shared, similar, and unique needs of each research area with respect to the broader underground research community's Underground Facilities and Infrastructure needs. Where potential conflicts of usage type, site, or duration might arise, these are identified. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: Snowmass 2021 Topical Report (UF5)

JWST's PEARLS: dust attenuation and gravitational lensing in the backlit-galaxy system VV 191

Authors: William C. Keel, Rogier A. Windhorst, Rolf A. Jansen, Seth H. Cohen, Jake Summers, Benne Holwerda, Sarah T. Bradford, Clayton D. Robertson, Giovanni Ferrami, Stuart Wyithe, Haojing Yan, Christopher J. Conselice, Simon P. Driver, Aaron Robotham, Norman A. Grogin, Christopher N. A. Willmer, Anton M. Koekemoer, Brenda L. Frye, Nimish P. Hathi, Russell E. Ryan, Jr., Nor Pirzkal, Madeline A. Marshall, Dan Coe, Jose M. Diego, Thomas J. Broadhurst , et al. (6 additional authors not shown)

Abstract: We derive the spatial and wavelength behavior of dust attenuation in the multiple-armed spiral galaxy VV191b using backlighting by the superimposed elliptical system VV191a in a pair with an exceptionally favorable geometry for this measurement. Imaging using JWST and HST spans the wavelength range 0.3-4.5 microns with high angular resolution, tracing the dust in detail from 0.6 to 1.5 microns. Di… ▽ More We derive the spatial and wavelength behavior of dust attenuation in the multiple-armed spiral galaxy VV191b using backlighting by the superimposed elliptical system VV191a in a pair with an exceptionally favorable geometry for this measurement. Imaging using JWST and HST spans the wavelength range 0.3-4.5 microns with high angular resolution, tracing the dust in detail from 0.6 to 1.5 microns. Distinct dust lanes continue well beyond the bright spiral arms, and trace a complex web, with a very sharp radial cutoff near 1.7 Petrosian radii. We present attenuation profiles and coverage statistics in each band at radii 14-21 kpc. We derive the attenuation law with wavelength; the data both within and between the dust lanes clearly favor a stronger reddening behavior (R ~ 2.0 between 0.6 and 0.9 microns, approaching unity by 1.5 microns) than found for starbursts and star-forming regions of galaxies. Power-law extinction behavior lambda^(-beta) gives beta=2.1 from 0.6-0.9 microns. R decreases at increasing wavelengths (R~1.1 between 0.9 and 1.5 microns), while beta steepens to 2.5. Mixing regions of different column density flattens the wavelength behavior, so these results suggest a different grain population than in our vicinity. The NIRCam images reveal a lens arc and counterimage from a background galaxy at z~1, spanning 90 degrees azimuthally at 2.8" from the foreground elliptical galaxy nucleus, and an additional weakly-lensed galaxy. The lens model and imaging data give a mass/light ratio 7.6 in solar units within the Einstein radius 2.0 kpc. △ Less

Submitted 21 February, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: Accepted by Astron. J. Analysis redone since submission, using updated JWST calibrations. Dust reddening behavior is steeper with wavelength and lensed galaxy redshift lower than we first derived

Online multiple hypothesis testing

Authors: David S. Robertson, James M. S. Wason, Aaditya Ramdas

Abstract: Modern data analysis frequently involves large-scale hypothesis testing, which naturally gives rise to the problem of maintaining control of a suitable type I error rate, such as the false discovery rate (FDR). In many biomedical and technological applications, an additional complexity is that hypotheses are tested in an online manner, one-by-one over time. However, traditional procedures that con… ▽ More Modern data analysis frequently involves large-scale hypothesis testing, which naturally gives rise to the problem of maintaining control of a suitable type I error rate, such as the false discovery rate (FDR). In many biomedical and technological applications, an additional complexity is that hypotheses are tested in an online manner, one-by-one over time. However, traditional procedures that control the FDR, such as the Benjamini-Hochberg procedure, assume that all p-values are available to be tested at a single time point. To address these challenges, a new field of methodology has developed over the past 15 years showing how to control error rates for online multiple hypothesis testing. In this framework, hypotheses arrive in a stream, and at each time point the analyst decides whether to reject the current hypothesis based both on the evidence against it, and on the previous rejection decisions. In this paper, we present a comprehensive exposition of the literature on online error rate control, with a review of key theory as well as a focus on applied examples. We also provide simulation results comparing different online testing algorithms and an up-to-date overview of the many methodological extensions that have been proposed. △ Less

Submitted 24 July, 2023; v1 submitted 24 August, 2022; originally announced August 2022.

Comments: Updated in response to reviewer comments

MSC Class: 62-02

Direct measurement of the low energy resonances in $^{22}\rm{Ne}(αあるふぁ,γがんま)^{26}\rm{Mg}$ reaction

Authors: S. Shahina, J. Gorres, D. Robertson, M. Couder, O. Gomez, A. Gula, M. Hanhardt, T. Kadlecek, R. Kelmar, P. Scholz, A. Simon, E. Stech, F. Strieder, M. Wiescher

Abstract: The $^{22}\rm{Ne}(αあるふぁ,γがんま)^{26}\rm{Mg}$ is an important reaction in stellar helium burning environments as it competes directly with one of the main neutron sources for the s-process, the $^{22}\rm{Ne}(αあるふぁ,n)^{25}\rm{Mg}$ reaction. The reaction rate of the $^{22}\rm{Ne}(αあるふぁ,γがんま)^{26}\rm{Mg}$ is dominated by the low energy resonances at $E_αあるふぁ^{lab}$ = 650 and 830 keV respectively. The $E_αあるふぁ^{lab}$ = 830 keV re… ▽ More The $^{22}\rm{Ne}(αあるふぁ,γがんま)^{26}\rm{Mg}$ is an important reaction in stellar helium burning environments as it competes directly with one of the main neutron sources for the s-process, the $^{22}\rm{Ne}(αあるふぁ,n)^{25}\rm{Mg}$ reaction. The reaction rate of the $^{22}\rm{Ne}(αあるふぁ,γがんま)^{26}\rm{Mg}$ is dominated by the low energy resonances at $E_αあるふぁ^{lab}$ = 650 and 830 keV respectively. The $E_αあるふぁ^{lab}$ = 830 keV resonance has been measured previously, but there are some uncertainties in the previous measurements. We confirmed the measurement of the $E_αあるふぁ^{lab}$ = 830 keV resonance using implanted $^{22}$Ne targets. We obtained a resonance strength of $ωおめがγがんま$ = 35 $\pm$ 4 $μみゅーeV$, and provide a weighted average of the present and previous measurements of $ωおめがγがんま$ = 35 $\pm$ 2 $μみゅーeV$ with reduced uncertainties compared to previous studies. We also attempted to measure the strength of the predicted resonance at $E_αあるふぁ^{lab}$ = 650 keV directly for the first time and found an upper limit of $ωおめがγがんま$ $\mathrm{<0.15}$ $μみゅーeV$ for the strength of this resonance. In addition, we also studied the $E_{P}^{lab}$= 851 keV resonance in $^{22}\rm{Ne}(p,γがんま)^{23}\rm{Na}$, and obtained a resonance strength of $ωおめがγがんま$ = 9.2 $\pm$ 0.7 eV with significantly lower uncertainties compared to previous measurements. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 10 pages, 6 figures

A proof of Erdős's $B+B+t$ conjecture

Authors: Bryna Kra, Joel Moreira, Florian K. Richter, Donald Robertson

Abstract: We show that every set $A$ of natural numbers with positive upper density can be shifted to contain the restricted sumset $\{b_1 + b_2 : b_1, b_2\in B \text{ and } b_1 \neq b_2 \}$ for some infinite set $B \subset A$. We show that every set $A$ of natural numbers with positive upper density can be shifted to contain the restricted sumset $\{b_1 + b_2 : b_1, b_2\in B \text{ and } b_1 \neq b_2 \}$ for some infinite set $B \subset A$. △ Less

Submitted 6 November, 2023; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 14 pages. Changes after referee comments, improved exposition

MSC Class: 05D10; 11B13; 37A05; 11B30

Online FDR Control for RNAseq Data

Authors: Lathan Liou, Milena Hornburg, David S. Robertson

Abstract: Motivation: While the analysis of a single RNA sequencing (RNAseq) dataset has been well described in the literature, modern research workflows often have additional complexity in that related RNAseq experiments are performed sequentially over time. The simplest and most widely used analysis strategy ignores the temporal aspects and analyses each dataset separately. However, this can lead to subst… ▽ More Motivation: While the analysis of a single RNA sequencing (RNAseq) dataset has been well described in the literature, modern research workflows often have additional complexity in that related RNAseq experiments are performed sequentially over time. The simplest and most widely used analysis strategy ignores the temporal aspects and analyses each dataset separately. However, this can lead to substantial inflation of the overall false discovery rate (FDR). We propose applying recently developed methodology for online hypothesis testing to analyse sequential RNAseq experiments in a principled way, guaranteeing FDR control at all times while never changing past decisions. Results: We show that standard offline approaches have variable control of FDR of related RNAseq experiments over time and a naively composed approach may improperly change historic decisions. We demonstrate that the online FDR algorithms are a principled way to guarantee control of FDR. Furthermore, in certain simulation scenarios, we observe empirically that online approaches have comparable power to offline approaches. Availability and Implementation: The onlineFDR package is freely available at http: //www.bioconductor.org/packages/onlineFDR. Additional code used for the simulation studies can be found at https://github.com/latlio/onlinefdr_rnaseq_simulation △ Less

Submitted 5 June, 2022; originally announced June 2022.

Infinite Sumsets in Sets with Positive Density

Authors: Bryna Kra, Joel Moreira, Florian K. Richter, Donald Robertson

Abstract: Motivated by questions asked by Erdos, we prove that any set $A\subset{\mathbb N}$ with positive upper density contains, for any $k\in{\mathbb N}$, a sumset $B_1+\cdots+B_k$, where $B_1,\dots,B_k\subset{\mathbb N}$ are infinite. Our proof uses ergodic theory and relies on structural results for measure preserving systems. Our techniques are new, even for the previously known case of $k=2$. Motivated by questions asked by Erdos, we prove that any set $A\subset{\mathbb N}$ with positive upper density contains, for any $k\in{\mathbb N}$, a sumset $B_1+\cdots+B_k$, where $B_1,\dots,B_k\subset{\mathbb N}$ are infinite. Our proof uses ergodic theory and relies on structural results for measure preserving systems. Our techniques are new, even for the previously known case of $k=2$. △ Less

Submitted 22 February, 2024; v1 submitted 3 June, 2022; originally announced June 2022.

Comments: 49 Pages

MSC Class: 05D10 11B13 37A05 11B30

Some performance considerations when using multi-armed bandit algorithms in the presence of missing data

Authors: Xijin Chen, Kim May Lee, Sofia S. Villar, David S. Robertson

Abstract: When comparing the performance of multi-armed bandit algorithms, the potential impact of missing data is often overlooked. In practice, it also affects their implementation where the simplest approach to overcome this is to continue to sample according to the original bandit algorithm, ignoring missing outcomes. We investigate the impact on performance of this approach to deal with missing data fo… ▽ More When comparing the performance of multi-armed bandit algorithms, the potential impact of missing data is often overlooked. In practice, it also affects their implementation where the simplest approach to overcome this is to continue to sample according to the original bandit algorithm, ignoring missing outcomes. We investigate the impact on performance of this approach to deal with missing data for several bandit algorithms through an extensive simulation study assuming the rewards are missing at random. We focus on two-armed bandit algorithms with binary outcomes in the context of patient allocation for clinical trials with relatively small sample sizes. However, our results apply to other applications of bandit algorithms where missing data is expected to occur. We assess the resulting operating characteristics, including the expected reward. Different probabilities of missingness in both arms are considered. The key finding of our work is that when using the simplest strategy of ignoring missing data, the impact on the expected performance of multi-armed bandit strategies varies according to the way these strategies balance the exploration-exploitation trade-off. Algorithms that are geared towards exploration continue to assign samples to the arm with more missing responses (which being perceived as the arm with less observed information is deemed more appealing by the algorithm than it would otherwise be). In contrast, algorithms that are geared towards exploitation would rapidly assign a high value to samples from the arms with a current high mean irrespective of the level observations per arm. Furthermore, for algorithms focusing more on exploration, we illustrate that the problem of missing responses can be alleviated using a simple mean imputation approach. △ Less

Submitted 7 July, 2022; v1 submitted 8 May, 2022; originally announced May 2022.

Comments: 30 pages, 6 figures

Familywise error rate control for block response-adaptive randomization

Authors: Ekkehard Glimm, David Robertson

Abstract: Response-adaptive randomization allows the probabilities of allocating patients to treatments in a clinical trial to change based on the previously observed response data, in order to achieve different experimental goals. One concern over the use of such designs in practice, particularly from a regulatory viewpoint, is controlling the type I error rate. To address this, Robertson and Wason (Biomet… ▽ More Response-adaptive randomization allows the probabilities of allocating patients to treatments in a clinical trial to change based on the previously observed response data, in order to achieve different experimental goals. One concern over the use of such designs in practice, particularly from a regulatory viewpoint, is controlling the type I error rate. To address this, Robertson and Wason (Biometrics, 2019) proposed methodology that guarantees familywise error rate control for a large class of response-adaptive designs. In this paper, we propose an improvement of their proposal that is conceptually simpler, in the specific context of block-randomised trials with a fixed allocation to the control arm. We show the modified method guarantees that there will never be negative weights for blocks of data, and can also provide a substantial power advantage in practice. △ Less

Submitted 12 April, 2022; originally announced April 2022.

MSC Class: 62L10

CUBES Phase A design overview -- The Cassegrain U-Band Efficient Spectrograph for the Very Large Telescope

Authors: Alessio Zanutta, Stefano Cristiani, David Atkinson, Veronica Baldini, Andrea Balestra, Beatriz Barbuy, Vanessa Bawden P. Macanhan, Ariadna Calcines, Giorgio Calderone, Scott Case, Bruno V. Castilho, Gabriele Cescutti, Roberto Cirami, Igor Coretti, Stefano Covino, Guido Cupani, Vincenzo De Caprio, Hans Dekker, Paolo Di Marcantonio, Valentina D'Odorico, Heitor Ernandes, Chris Evans, Tobias Feger, Carmen Feiz, Mariagrazia Franchini , et al. (29 additional authors not shown)

Abstract: We present the baseline conceptual design of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope. CUBES will provide unprecedented sensitivity for spectroscopy on a 8 - 10 m class telescope in the ground ultraviolet (UV), spanning a bandwidth of > 100 nm that starts at 300 nm, the shortest wavelength accessible from the ground. The design has been optimized for end-to… ▽ More We present the baseline conceptual design of the Cassegrain U-Band Efficient Spectrograph (CUBES) for the Very Large Telescope. CUBES will provide unprecedented sensitivity for spectroscopy on a 8 - 10 m class telescope in the ground ultraviolet (UV), spanning a bandwidth of > 100 nm that starts at 300 nm, the shortest wavelength accessible from the ground. The design has been optimized for end-to-end efficiency and provides a spectral resolving power of R > 20000, that will unlock a broad range of new topics across solar system, Galactic and extraglactic astronomy. The design also features a second, lower-resolution (R \sim 7000) mode and has the option of a fiberlink to the UVES instrument for simultaneous observations at longer wavelengths. Here we present the optical, mechanical and software design of the various subsystems of the instrument after the Phase A study of the project. We discuss the expected performances for the layout choices and highlight some of the performance trade-offs considered to best meet the instrument top-level requirements. We also introduce the model-based system engineering approach used to organize and manage the project activities and interfaces, in the context that it is increasingly necessary to integrate such tools in the development of complex astronomical projects. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: Accepted for publication in Experimental Astronomy

Online error control for platform trials

Authors: David S. Robertson, James M. S. Wason, Franz König, Martin Posch, Thomas Jaki

Abstract: Platform trials evaluate multiple experimental treatments under a single master protocol, where new treatment arms are added to the trial over time. Given the multiple treatment comparisons, there is the potential for inflation of the overall type I error rate, which is complicated by the fact that the hypotheses are tested at different times and are not all necessarily pre-specified. Online error… ▽ More Platform trials evaluate multiple experimental treatments under a single master protocol, where new treatment arms are added to the trial over time. Given the multiple treatment comparisons, there is the potential for inflation of the overall type I error rate, which is complicated by the fact that the hypotheses are tested at different times and are not all necessarily pre-specified. Online error control methodology provides a possible solution to the problem of multiplicity for platform trials where a relatively large number of hypotheses are expected to be tested over time. In the online testing framework, hypotheses are tested in a sequential manner, where at each time-step an analyst decides whether to reject the current null hypothesis without knowledge of future tests but based solely on past decisions. Methodology has recently been developed for online control of the false discovery rate as well as the familywise error rate (FWER). In this paper, we describe how to apply online error control to the platform trial setting, present extensive simulation results, and give some recommendations for the use of this new methodology in practice. We show that the algorithms for online error rate control can have a substantially lower FWER than uncorrected testing, while still achieving noticeable gains in power when compared with the use of a Bonferroni procedure. We also illustrate how online error control would have impacted a currently ongoing platform trial. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 26 pages, 13 figures

MSC Class: 62L10

The Status and Future of Direct Nuclear Reaction Measurements for Stellar Burning

Authors: M. Aliotta, R. Buompane, M. Couder, A. Couture, R. J. deBoer, A. Formicola, L. Gialanella, J. Glorius, G. Imbriani, M. Junker, C. Langer, A. Lennarz, Yu. A. Litvinov, W. -P. Liu, M. Lugaro, C. Matei, Z. Meisel, L. Piersanti, R. Reifarth, D. Robertson, A. Simon, O. Straniero, A. Tumino, M. Wiescher, Y. Xu

Abstract: The study of stellar burning began just over 100 years ago. Nonetheless, we do not yet have a detailed picture of the nucleosynthesis within stars and how nucleosynthesis impacts stellar structure and the remnants of stellar evolution. Achieving this understanding will require precise direct measurements of the nuclear reactions involved. This report summarizes the status of direct measurements fo… ▽ More The study of stellar burning began just over 100 years ago. Nonetheless, we do not yet have a detailed picture of the nucleosynthesis within stars and how nucleosynthesis impacts stellar structure and the remnants of stellar evolution. Achieving this understanding will require precise direct measurements of the nuclear reactions involved. This report summarizes the status of direct measurements for stellar burning, focusing on developments of the last couple of decades, and offering a prospectus of near-future developments. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: Accepted to Journal of Physics G as a Major Report. Corresponding author: Zach Meisel (meisel@ohio.edu)

Point estimation for adaptive trial designs I: a methodological review

Authors: David S. Robertson, Babak Choodari-Oskooei, Munya Dimairo, Laura Flight, Philip Pallmann, Thomas Jaki

Abstract: Recent FDA guidance on adaptive clinical trial designs defines bias as "a systematic tendency for the estimate of treatment effect to deviate from its true value", and states that it is desirable to obtain and report estimates of treatment effects that reduce or remove this bias. The conventional end-of-trial point estimates of the treatment effects are prone to bias in many adaptive designs, beca… ▽ More Recent FDA guidance on adaptive clinical trial designs defines bias as "a systematic tendency for the estimate of treatment effect to deviate from its true value", and states that it is desirable to obtain and report estimates of treatment effects that reduce or remove this bias. The conventional end-of-trial point estimates of the treatment effects are prone to bias in many adaptive designs, because they do not take into account the potential and realised trial adaptations. While much of the methodological developments on adaptive designs have tended to focus on control of type I error rates and power considerations, in contrast the question of biased estimation has received relatively less attention. This paper is the first in a two-part series that studies the issue of potential bias in point estimation for adaptive trials. Part I provides a comprehensive review of the methods to remove or reduce the potential bias in point estimation of treatment effects for adaptive designs, while part II illustrates how to implement these in practice and proposes a set of guidelines for trial statisticians. The methods reviewed in this paper can be broadly classified into unbiased and bias-reduced estimation, and we also provide a classification of estimators by the type of adaptive design. We compare the proposed methods, highlight available software and code, and discuss potential methodological gaps in the literature. △ Less

Submitted 28 November, 2022; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: Part I of a two-part series; accepted in Statistics in Medicine

MSC Class: 62F10

The impact and recovery of asteroid 2018 LA

Authors: Peter Jenniskens, Mohutsiwa Gabadirwe, Qing-Zhu Yin, Alexander Proyer, Oliver Moses, Tomas Kohout, Fulvio Franchi, Roger L. Gibson, Richard Kowalski, Eric J. Christensen, Alex R. Gibbs, Aren Heinze, Larry Denneau, Davide Farnocchia, Paul W. Chodas, William Gray, Marco Micheli, Nick Moskovitz, Christopher A. Onken, Christian Wolf, Hadrien A. R. Devillepoix, Quanzhi Ye, Darrel K. Robertson, Peter Brown, Esko Lyytinen , et al. (41 additional authors not shown)

Abstract: The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to… ▽ More The June 2, 2018, impact of asteroid 2018 LA over Botswana is only the second asteroid detected in space prior to impacting over land. Here, we report on the successful recovery of meteorites. Additional astrometric data refine the approach orbit and define the spin period and shape of the asteroid. Video observations of the fireball constrain the asteroid's position in its orbit and were used to triangulate the location of the fireball's main flare over the Central Kalahari Game Reserve. 23 meteorites were recovered. A consortium study of eight of these classifies Motopi Pan as a HED polymict breccia derived from howardite, cumulate and basaltic eucrite, and diogenite lithologies. Before impact, 2018 LA was a solid rock of about 156 cm diameter with high bulk density about 2.85 g/cm3, a relatively low albedo pV about 0.25, no significant opposition effect on the asteroid brightness, and an impact kinetic energy of about 0.2 kt. The orbit of 2018 LA is consistent with an origin at Vesta (or its Vestoids) and delivery into an Earth-impacting orbit via the nu_6 resonance. The impact that ejected 2018 LA in an orbit towards Earth occurred 22.8 +/- 3.8 Ma ago. Zircons record a concordant U-Pb age of 4563 +/- 11 Ma and a consistent 207Pb/206Pb age of 4563 +/- 6 Ma. A much younger Pb-Pb phosphate resetting age of 4234 +/- 41 Ma was found. From this impact chronology, we discuss what is the possible source crater of Motopi Pan and the age of Vesta's Veneneia impact basin. △ Less

Submitted 12 May, 2021; originally announced May 2021.

Comments: Meteoritics & Planetary Science (2021)

Determinantal formulas with major indices

Authors: Thomas McConville, Donald Robertson, Clifford Smyth

Abstract: We give a simple proof of a major index determinant formula in the symmetric group discovered by Krattenthaler and first proved by Thibon using noncommutative symmetric functions. We do so by proving a factorization of an element in the group ring of the symmetric group. By applying similar methods to the groups of signed permutations and colored permutations, we prove determinant formulas in thes… ▽ More We give a simple proof of a major index determinant formula in the symmetric group discovered by Krattenthaler and first proved by Thibon using noncommutative symmetric functions. We do so by proving a factorization of an element in the group ring of the symmetric group. By applying similar methods to the groups of signed permutations and colored permutations, we prove determinant formulas in these groups as conjectured by Krattenthaler. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 14 pages

MSC Class: 05A05

The Sariçiçek howardite fall in Turkey: Source crater of HED meteorites on Vesta and impact risk of Vestoids

Authors: Ozan Unsalan, Peter Jenniskens, Qing-Zhu Yin, Ersin Kaygisiz, Jim Albers, David L. Clark, Mikael Granvik, Iskender Demirkol, Ibrahim Y. Erdogan, Aydin S. Bengu, Mehmet E. Özel, Zahide Terzioglu, Nayeob GI, Peter Brown, Esref Yalcinkaya, Tuğba Temel, Dinesh K. Prabhu, Darrel K. Robertson, Mark Boslough, Daniel R. Ostrowski, Jamie Kimberley, Selman ER, Douglas J. Rowland, Kathryn L. Bryson, Cisem Altunayar-Unsalan , et al. (54 additional authors not shown)

Abstract: The Sariçiçek howardite meteorite shower consisting of 343 documented stones occurred on 2 September 2015 in Turkey and is the first documented howardite fall. Cosmogenic isotopes show that Sariçiçek experienced a complex cosmic ray exposure history, exposed during ~12-14 Ma in a regolith near the surface of a parent asteroid, and that an ca.1 m sized meteoroid was launched by an impact 22 +/- 2 M… ▽ More The Sariçiçek howardite meteorite shower consisting of 343 documented stones occurred on 2 September 2015 in Turkey and is the first documented howardite fall. Cosmogenic isotopes show that Sariçiçek experienced a complex cosmic ray exposure history, exposed during ~12-14 Ma in a regolith near the surface of a parent asteroid, and that an ca.1 m sized meteoroid was launched by an impact 22 +/- 2 Ma ago to Earth (as did one third of all HED meteorites). SIMS dating of zircon and baddeleyite yielded 4550.4 +/- 2.5 Ma and 4553 +/- 8.8 Ma crystallization ages for the basaltic magma clasts. The apatite U-Pb age of 4525 +/- 17 Ma, K-Ar age of ~3.9 Ga, and the U,Th-He ages of 1.8 +/- 0.7 and 2.6 +/- 0.3 Ga are interpreted to represent thermal metamorphic and impact-related resetting ages, respectively. Petrographic, geochemical and O-, Cr- and Ti- isotopic studies confirm that Sariçiçek belongs to the normal clan of HED meteorites. Petrographic observations and analysis of organic material indicate a small portion of carbonaceous chondrite material in the Sariçiçek regolith and organic contamination of the meteorite after a few days on soil. Video observations of the fall show an atmospheric entry at 17.3 +/- 0.8 kms-1 from NW, fragmentations at 37, 33, 31 and 27 km altitude, and provide a pre-atmospheric orbit that is the first dynamical link between the normal HED meteorite clan and the inner Main Belt. Spectral data indicate the similarity of Sariçiçek with the Vesta asteroid family spectra, a group of asteroids stretching to delivery resonances, which includes (4) Vesta. Dynamical modeling of meteoroid delivery to Earth shows that the disruption of a ca.1 km sized Vesta family asteroid or a ~10 km sized impact crater on Vesta is required to provide sufficient meteoroids <4 m in size to account for the influx of meteorites from this HED clan. △ Less

Submitted 7 February, 2021; originally announced February 2021.

Journal ref: Unsalan, Ozan, Peter Jenniskens,et al. "The Sariçiçek howardite fall in Turkey: Source crater of HED meteorites on Vesta and impact risk of Vestoids." Meteoritics & Planetary Science 54, no. 5 (2019): 953-1008

Conditional Power and Friends: The Why and How of (Un)planned, Unblinded Sample Size Recalculations in Confirmatory Trials

Authors: Kevin Kunzmann, Michael J. Grayling, Kim M. Lee, David S. Robertson, Kaspar Rufibach, James M. S. Wason

Abstract: Adapting the final sample size of a trial to the evidence accruing during the trial is a natural way to address planning uncertainty. Designs with adaptive sample size need to account for their optional stopping to guarantee strict type-I error-rate control. A variety of different methods to maintain type-I error-rate control after unplanned changes of the initial sample size have been proposed in… ▽ More Adapting the final sample size of a trial to the evidence accruing during the trial is a natural way to address planning uncertainty. Designs with adaptive sample size need to account for their optional stopping to guarantee strict type-I error-rate control. A variety of different methods to maintain type-I error-rate control after unplanned changes of the initial sample size have been proposed in the literature. This makes interim analyses for the purpose of sample size recalculation feasible in a regulatory context. Since the sample size is usually determined via an argument based on the power of the trial, an interim analysis raises the question of how the final sample size should be determined conditional on the accrued information. Conditional power is a concept often put forward in this context. Since it depends on the unknown effect size, we take a strict estimation perspective and compare assumed conditional power, observed conditional power, and predictive power with respect to their properties as estimators of the unknown conditional power. We then demonstrate that pre-planning an interim analysis using methodology for unplanned interim analyses is ineffective and naturally leads to the concept of optimal two-stage designs. We conclude that unplanned design adaptations should only be conducted as reaction to trial-external new evidence, operational needs to violate the originally chosen design, or post hoc changes in the objective criterion. Finally, we show that commonly discussed sample size recalculation rules can lead to paradoxical outcomes and propose two alternative ways of reacting to newly emerging trial-external evidence. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Optical Hemodynamic Imaging of Jugular Venous Dynamics During Altered Central Venous Pressure

Authors: Robert Amelard, Andrew D Robertson, Courtney A Patterson, Hannah Heigold, Essi Saarikoski, Richard L Hughson

Abstract: An optical imaging system is proposed for quantitatively assessing jugular venous response to altered central venous pressure. The proposed system assesses sub-surface optical absorption changes from jugular venous waveforms with a spatial calibration procedure to normalize incident tissue illumination. Widefield frames of the right lateral neck were captured and calibrated using a novel flexible… ▽ More An optical imaging system is proposed for quantitatively assessing jugular venous response to altered central venous pressure. The proposed system assesses sub-surface optical absorption changes from jugular venous waveforms with a spatial calibration procedure to normalize incident tissue illumination. Widefield frames of the right lateral neck were captured and calibrated using a novel flexible surface calibration method. A hemodynamic optical model was derived to quantify jugular venous optical attenuation (JVA) signals, and generate a spatial jugular venous pulsatility map. JVA was assessed in three cardiovascular protocols that altered central venous pressure: acute central hypovolemia (lower body negative pressure), venous congestion (head-down tilt), and impaired cardiac filling (Valsalva maneuver). JVA waveforms exhibited biphasic wave properties consistent with jugular venous pulse dynamics when time-aligned with an electrocardiogram. JVA correlated strongly (median, interquartile range) with invasive central venous pressure during graded central hypovolemia (r=0.85, [0.72, 0.95]), graded venous congestion (r=0.94, [0.84, 0.99]), and impaired cardiac filling (r=0.94, [0.85, 0.99]). Reduced JVA during graded acute hypovolemia was strongly correlated with reductions in stroke volume (SV) (r=0.85, [0.76, 0.92]) from baseline (SV: 79$\pm$15 mL, JVA: 0.56$\pm$0.10 a.u.) to -40 mmHg suction (SV: 59$\pm$18 mL, JVA: 0.47$\pm$0.05 a.u.; p$<$0.01). The proposed non-contact optical imaging system demonstrated jugular venous dynamics consistent with invasive central venous monitoring during three protocols that altered central venous pressure. This system provides non-invasive monitoring of pressure-induced jugular venous dynamics in clinically relevant conditions where catheterization is traditionally required, enabling monitoring in non-surgical environments. △ Less

Submitted 24 March, 2021; v1 submitted 22 July, 2020; originally announced July 2020.

A review of Bayesian perspectives on sample size derivation for confirmatory trials

Authors: Kevin Kunzmann, Michael J. Grayling, Kim May Lee, David S. Robertson, Kaspar Rufibach, James M. S. Wason

Abstract: Sample size derivation is a crucial element of the planning phase of any confirmatory trial. A sample size is typically derived based on constraints on the maximal acceptable type I error rate and a minimal desired power. Here, power depends on the unknown true effect size. In practice, power is typically calculated either for the smallest relevant effect size or a likely point alternative. The fo… ▽ More Sample size derivation is a crucial element of the planning phase of any confirmatory trial. A sample size is typically derived based on constraints on the maximal acceptable type I error rate and a minimal desired power. Here, power depends on the unknown true effect size. In practice, power is typically calculated either for the smallest relevant effect size or a likely point alternative. The former might be problematic if the minimal relevant effect is close to the null, thus requiring an excessively large sample size. The latter is dubious since it does not account for the a priori uncertainty about the likely alternative effect size. A Bayesian perspective on the sample size derivation for a frequentist trial naturally emerges as a way of reconciling arguments about the relative a priori plausibility of alternative effect sizes with ideas based on the relevance of effect sizes. Many suggestions as to how such `hybrid' approaches could be implemented in practice have been put forward in the literature. However, key quantities such as assurance, probability of success, or expected power are often defined in subtly different ways in the literature. Starting from the traditional and entirely frequentist approach to sample size derivation, we derive consistent definitions for the most commonly used `hybrid' quantities and highlight connections, before discussing and demonstrating their use in the context of sample size derivation for clinical trials. △ Less

Submitted 28 June, 2020; originally announced June 2020.

Journal ref: Am. Stat., 2021, 75(4), 424--432

Spacecraft and interplanetary contributions to the magnetic environment on-board LISA Pathfinder

Authors: M. Armano, H. Audley, J. Baird, P. Binetruy, M. Born, D. Bortoluzzi, E. Castelli, A. Cavalleri, A. Cesarini, A. M. Cruise, K. Danzmann, M. de Deus Silva, I. Diepholz, G. Dixon, R. Dolesi, L. Ferraioli, V. Ferroni, E. D. Fitzsimons, M. Freschi, L. Gesa, F. Gibert, D. Giardini, R. Giusteri, C. Grimani, J. Grzymisch , et al. (57 additional authors not shown)

Abstract: LISA Pathfinder (LPF) has been a space-based mission designed to test new technologies that will be required for a gravitational wave observatory in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime (mHz and below), the measurement band of interest for a space-based observatory. The magnetic field can couple to the magnetic susceptibility a… ▽ More LISA Pathfinder (LPF) has been a space-based mission designed to test new technologies that will be required for a gravitational wave observatory in space. Magnetically driven forces play a key role in the instrument sensitivity in the low-frequency regime (mHz and below), the measurement band of interest for a space-based observatory. The magnetic field can couple to the magnetic susceptibility and remanent magnetic moment from the test masses and disturb them from their geodesic movement. LISA Pathfinder carried on-board a dedicated magnetic measurement subsystem with noise levels of 10 $ \rm nT \ Hz^{-1/2}$ from 1 Hzへるつ down to 1 mHz. In this paper we report on the magnetic measurements throughout LISA Pathfinder operations. We characterise the magnetic environment within the spacecraft, study the time evolution of the magnetic field and its stability down to 20 $μみゅー$Hz, where we measure values around 200 $ \rm nT \ Hz^{-1/2}$ and identify two different frequency regimes, one related to the interplanetary magnetic field and the other to the magnetic field originating inside the spacecraft. Finally, we characterise the non-stationary component of the fluctuations of the magnetic field below the mHz and relate them to the dynamics of the solar wind. △ Less

Submitted 7 May, 2020; originally announced May 2020.

Comments: 16 pages, 17 figures. MNRAS LaTeX style file version 3.0

Journal ref: Monthly Notices of the Royal Astronomical Society, 2020, 494.2: 3014-3027

Response-adaptive randomization in clinical trials: from myths to practical considerations

Authors: David S. Robertson, Kim May Lee, Boryana C. Lopez-Kolkovska, Sofia S. Villar

Abstract: Response-Adaptive Randomization (RAR) is part of a wider class of data-dependent sampling algorithms, for which clinical trials are typically used as a motivating application. In that context, patient allocation to treatments is determined by randomization probabilities that change based on the accrued response data in order to achieve experimental goals. RAR has received abundant theoretical atte… ▽ More Response-Adaptive Randomization (RAR) is part of a wider class of data-dependent sampling algorithms, for which clinical trials are typically used as a motivating application. In that context, patient allocation to treatments is determined by randomization probabilities that change based on the accrued response data in order to achieve experimental goals. RAR has received abundant theoretical attention from the biostatistical literature since the 1930's and has been the subject of numerous debates. In the last decade, it has received renewed consideration from the applied and methodological communities, driven by well-known practical examples and its widespread use in machine learning. Papers on the subject present different views on its usefulness, and these are not easy to reconcile. This work aims to address this gap by providing a unified, broad and fresh review of methodological and practical issues to consider when debating the use of RAR in clinical trials. △ Less

Submitted 7 June, 2022; v1 submitted 1 May, 2020; originally announced May 2020.

Comments: Update in response to editor comments

MSC Class: 62-02

Graphical approaches for the control of generalised error rates

Authors: David S. Robertson, James M. S. Wason, Frank Bretz

Abstract: When simultaneously testing multiple hypotheses, the usual approach in the context of confirmatory clinical trials is to control the familywise error rate (FWER), which bounds the probability of making at least one false rejection. In many trial settings, these hypotheses will additionally have a hierarchical structure that reflects the relative importance and links between different clinical obje… ▽ More When simultaneously testing multiple hypotheses, the usual approach in the context of confirmatory clinical trials is to control the familywise error rate (FWER), which bounds the probability of making at least one false rejection. In many trial settings, these hypotheses will additionally have a hierarchical structure that reflects the relative importance and links between different clinical objectives. The graphical approach of Bretz et al. (2009) is a flexible and easily communicable way of controlling the FWER while respecting complex trial objectives and multiple structured hypotheses. However, the FWER can be a very stringent criterion that leads to procedures with low power, and may not be appropriate in exploratory trial settings. This motivates controlling generalised error rates, particularly when the number of hypotheses tested is no longer small. We consider the generalised familywise error rate (k-FWER), which is the probability of making k or more false rejections, as well as the tail probability of the false discovery proportion (FDP), which is the probability that the proportion of false rejections is greater than some threshold. We also consider asymptotic control of the false discovery rate (FDR), which is the expectation of the FDP. In this paper, we show how to control these generalised error rates when using the graphical approach and its extensions. We demonstrate the utility of the resulting graphical procedures on three clinical trial case studies. △ Less

Submitted 2 June, 2020; v1 submitted 3 April, 2020; originally announced April 2020.

Comments: Accepted for publication in Statistics in Medicine

MSC Class: 62-7

Optical suppression of tilt-to-length coupling in the LISA long-arm interferometer

Authors: M Chwalla, K Danzmann, M Dovale Álvarez, J J Esteban Delgado, G Fernández Barranco, E Fitzsimons, O Gerberding, G Heinzel, C J Killow, M Lieser, M Perreur-Lloyd, D I Robertson, J M Rohr, S Schuster, T S Schwarze, M Tröbs, G Wanner, H Ward

Abstract: The arm length and the isolation in space enable LISA to probe for signals unattainable on ground, opening a window to the sub-Hz gravitational-wave universe. The coupling of unavoidable angular spacecraft jitter into the longitudinal displacement measurement, an effect known as tilt-to-length (TTL) coupling, is critical for realizing the required sensitivity of picometer$/\sqrt{\rm{Hz}}$. An ultr… ▽ More The arm length and the isolation in space enable LISA to probe for signals unattainable on ground, opening a window to the sub-Hz gravitational-wave universe. The coupling of unavoidable angular spacecraft jitter into the longitudinal displacement measurement, an effect known as tilt-to-length (TTL) coupling, is critical for realizing the required sensitivity of picometer$/\sqrt{\rm{Hz}}$. An ultra-stable interferometer testbed has been developed in order to investigate this issue and validate mitigation strategies in a setup representative of LISA, and in this paper it is operated in the long-arm interferometer configuration. The testbed is fitted with a flat-top beam generator to simulate the beam received by a LISA spacecraft. We demonstrate a reduction of TTL coupling between this flat-top beam and a Gaussian reference beam via introducing two- and four-lens imaging systems. TTL coupling factors below $\pm 25\,μみゅー$m/rad for beam tilts within $\pm 300\,μみゅー$rad are obtained by careful optimization of the system. Moreover we show that the additional TTL coupling due to lateral alignment errors of elements of the imaging system can be compensated by introducing lateral shifts of the detector, and vice versa. These findings help validate the suitability of this noise-reduction technique for the LISA long-arm interferometer. △ Less

Submitted 5 June, 2020; v1 submitted 13 February, 2020; originally announced February 2020.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. Applied 14, 014030 (2020)

Field-based mechanical phenotyping of cereal crops to assess lodging resistance

Authors: Lindsay Erndwein, Douglas D. Cook, Daniel J. Robertson, Erin E. Sparks

Abstract: Plant mechanical failure, also known as lodging, is the cause of significant and unpredictable yield losses in cereal crops. Lodging occurs in two distinct failure modes - stalk lodging and root lodging. Despite the prevalence and detrimental impact of lodging on crop yields, there is little consensus on how to phenotype plants in the field for lodging resistance and thus breed for mechanically re… ▽ More Plant mechanical failure, also known as lodging, is the cause of significant and unpredictable yield losses in cereal crops. Lodging occurs in two distinct failure modes - stalk lodging and root lodging. Despite the prevalence and detrimental impact of lodging on crop yields, there is little consensus on how to phenotype plants in the field for lodging resistance and thus breed for mechanically resilient plants. This review provides an overview of field-based mechanical testing approaches to assess stalk and root lodging resistance. These approaches are placed in the context of future perspectives. Best practices and recommendations for acquiring field-based mechanical phenotypes of plants are also presented. △ Less

Submitted 20 February, 2020; v1 submitted 18 September, 2019; originally announced September 2019.

LISA Pathfinder Performance Confirmed in an Open-Loop Configuration: Results from the Free-Fall Actuation Mode

Authors: M. Armano, H. Audley, J. Baird, P. Binetruy, M. Born, D. Bortoluzzi, E. Castelli, A. Cavalleri, A. Cesarini, A. M. Cruise, K. Danzmann, M. de Deus Silva, I. Diepholz, G. Dixon, R. Dolesi, L. Ferraioli, V. Ferroni, E. D. Fitzsimons, M. Freschi, L. Gesa, F. Gibert, D. Giardini, R. Giusteri, C. Grimani, J. Grzymisch , et al. (53 additional authors not shown)

Abstract: We report on the results of the LISA Pathfinder (LPF) free-fall mode experiment, in which the control force needed to compensate the quasistatic differential force acting on two test masses is applied intermittently as a series of "impulse" forces lasting a few seconds and separated by roughly 350 s periods of true free fall. This represents an alternative to the normal LPF mode of operation in wh… ▽ More We report on the results of the LISA Pathfinder (LPF) free-fall mode experiment, in which the control force needed to compensate the quasistatic differential force acting on two test masses is applied intermittently as a series of "impulse" forces lasting a few seconds and separated by roughly 350 s periods of true free fall. This represents an alternative to the normal LPF mode of operation in which this balancing force is applied continuously, with the advantage that the acceleration noise during free fall is measured in the absence of the actuation force, thus eliminating associated noise and force calibration errors. The differential acceleration noise measurement presented here with the free-fall mode agrees with noise measured with the continuous actuation scheme, representing an important and independent confirmation of the LPF result. An additional measurement with larger actuation forces also shows that the technique can be used to eliminate actuation noise when this is a dominant factor. △ Less

Submitted 30 August, 2019; originally announced August 2019.

Standard Model parameters in the tadpole-free pure $\overline{\rm{MS}}$ scheme

Authors: Stephen P. Martin, David G. Robertson

Abstract: We present an implementation and numerical study of the Standard Model couplings, masses, and vacuum expectation value (VEV), using the pure $\overline{\rm{MS}}$ renormalization scheme based on dimensional regularization. Here, the $\overline{\rm{MS}}$ Lagrangian parameters are treated as the fundamental inputs, and the VEV is defined as the minimum of the Landau gauge effective potential, so that… ▽ More We present an implementation and numerical study of the Standard Model couplings, masses, and vacuum expectation value (VEV), using the pure $\overline{\rm{MS}}$ renormalization scheme based on dimensional regularization. Here, the $\overline{\rm{MS}}$ Lagrangian parameters are treated as the fundamental inputs, and the VEV is defined as the minimum of the Landau gauge effective potential, so that tadpole diagrams vanish, resulting in improved convergence of perturbation theory. State-of-the-art calculations relating the $\overline{\rm{MS}}$ inputs to on-shell observables are implemented in a consistent way within a public computer code library, SMDR (Standard Model in Dimensional Regularization), which can be run interactively or called by other programs. Included here for the first time are the full 2-loop contributions to the Fermi constant within this scheme and studies of the minimization condition for the VEV at 3-loop order with 4-loop QCD effects. We also implement, and study the scale dependence of, all known multi-loop contributions to the physical masses of the Higgs boson, the W and Z bosons, and the top quark, the fine structure constant and weak mixing angle, and the renormalization group equations and threshold matching relations for the gauge couplings, fermion masses, and Yukawa couplings. △ Less

Submitted 26 May, 2022; v1 submitted 4 July, 2019; originally announced July 2019.

Comments: 26 pages. The SMDR code can be obtained from https://davidgrobertson.github.io/SMDR/ v2: references added. v3: code updated with 3-loop QCD corrections to Higgs, W, Z masses. v4: The parameterization of the complex pole masses of W,Z has been changed in the code

Journal ref: Phys. Rev. D 100, 073004 (2019)

Temperature stability in the sub-milliHertz band with LISA Pathfinder

Authors: M. Armano, H. Audley, J. Baird, P. Binetruy, M. Born, D. Bortoluzzi, E. Castelli, A. Cavalleri, A. Cesarini, A. M. Cruise, K. Danzmann, M. de Deus Silva, I. Diepholz, G. Dixon, R. Dolesi, L. Ferraioli, V. Ferroni, E. D. Fitzsimons, M. Freschi, L. Gesa, F. Gibert, D. Giardini, R. Giusteri, C. Grimani, J. Grzymisch , et al. (57 additional authors not shown)

Abstract: LISA Pathfinder (LPF) was a technology pioneering mission designed to test key technologies required for gravitational wave detection in space. In the low frequency regime (milli-Hertz and below), where space-based gravitational wave observatories will operate, temperature fluctuations play a crucial role since they can couple into the interferometric measurement and the test masses' free-fall acc… ▽ More LISA Pathfinder (LPF) was a technology pioneering mission designed to test key technologies required for gravitational wave detection in space. In the low frequency regime (milli-Hertz and below), where space-based gravitational wave observatories will operate, temperature fluctuations play a crucial role since they can couple into the interferometric measurement and the test masses' free-fall accuracy in many ways. A dedicated temperature measurement subsystem, with noise levels in 10$\,μみゅー$K$\,$Hz$^{-1/2}$ down to $1\,$mHz was part of the diagnostics unit on board LPF. In this paper we report on the temperature measurements throughout mission operations, characterize the thermal environment, estimate transfer functions between different locations and report temperature stability (and its time evolution) at frequencies as low as 10$\,μみゅー$Hz, where typically values around $1\,$K$\,$Hz$^{-1/2}$ were measured. △ Less

Submitted 22 May, 2019; originally announced May 2019.

Comments: 13 pages, 16 figures. MNRAS LaTeX style file version 3.0

Journal ref: Mon.Not.Roy.Astron.Soc 486 (2019) no.3, 3368-3379

Micrometeoroid Events in LISA Pathfinder

Authors: James Ira Thorpe, Jacob Slutsky, John Baker, Tyson Littenberg, Sophie Hourihane, Nicole Pagane, Petr Pokorny, Diego Janches, Michele Armano, Heather Audley, G. Auger, Jonathan Baird, Massimo Bassan, Pierre Binetruy, Michael Born, D. Bortoluzzi, N. Brandt, M. Caleno, A Cavalleri, A Cesarini, A. M. Cruise, K. Danzmann, M. de Deus Silva, R. De Rosa, L. Di Fiore , et al. (82 additional authors not shown)

Abstract: The zodiacal dust complex, a population of dust and small particles that pervades the Solar System, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. Here we present a new set of data obtained using a novel technique: direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible b… ▽ More The zodiacal dust complex, a population of dust and small particles that pervades the Solar System, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. Here we present a new set of data obtained using a novel technique: direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible by the extreme precision of the instruments flown on the LISA Pathfinder spacecraft, a technology demonstrator for a future space-based gravitational wave observatory that operated near the first Sun-Earth Lagrange point from early 2016 through Summer of 2017. Using a simple model of the impacts and knowledge of the control system, we show that it is possible to detect impacts and measure properties such as the transferred momentum (related to the particle's mass and velocity), direction of travel, and location of impact on the spacecraft. In this paper, we present the results of a systematic search for impacts during 4348 hours of Pathfinder data. We report a total of 54 candidates with momenta ranging from 0.2$\,μみゅー\textrm{Ns}$ to 230$\,μみゅー\textrm{Ns}$. We furthermore make a comparison of these candidates with models of micrometeoroid populations in the inner solar system including those resulting from Jupiter-family comets, Oort-cloud comets, Hailey-type comets, and Asteroids. We find that our measured population is consistent with a population dominated by Jupiter-family comets with some evidence for a smaller contribution from Hailey-type comets. This is in agreement with consensus models of the zodiacal dust complex in the momentum range sampled by LISA Pathfinder. △ Less

Submitted 2 October, 2019; v1 submitted 7 May, 2019; originally announced May 2019.

Comments: 22 pages, 14 figures, accepted in ApJ

Journal ref: JApJ 883 53 (2019)

An Introduction to Topological Data Analysis for Physicists: From LGM to FRBs

Authors: Jeff Murugan, Duncan Robertson

Abstract: Topological Data Analysis (TDA) is a novel, and relatively new approach to analysing high-dimensional data sets. It does this by focussing on global properties like the shape and connectivity of the data giving it a significant advantage over more conventional tools based on cluster analysis, a localised property of the data. However, some of its mathematical foundations, like algebraic topology a… ▽ More Topological Data Analysis (TDA) is a novel, and relatively new approach to analysing high-dimensional data sets. It does this by focussing on global properties like the shape and connectivity of the data giving it a significant advantage over more conventional tools based on cluster analysis, a localised property of the data. However, some of its mathematical foundations, like algebraic topology and discrete Morse theory, are perceived as an intimidatingly steep upramp into the subject. Consequently, it has enjoyed much less popularity as a data-analysis tool than less abstract methods. This article aims to change this. By focusing on a small set of simple examples, chosen primarily for their pedagogical value, we introduce and explain TDA's two principle branches; persistent homology and the Mapper algorithm. We then illustrate the universality of the method by discussing its application to the intriguing data set of fast radio burst (FRB) observations. We close the article with a discussion of the resilience of topological data analysis to noise and some statistical and computational challenges faced by the method. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: 47 pages, nearly the same number of figures

Forbush decreases and $<$ 2-day GCR flux non-recurrent variations studied with LISA Pathfinder

Authors: C. Grimani, M. Armano, H. Audley, J. Baird, S. Benella, P. Binetruy, M. Born, D. Bortoluzzi, E. Castelli, A. Cavalleri, A. Cesarini, A. M. Cruise, K. Danzmann, M. de Deus Silva, I. Diepholz, G. Dixon, R. Dolesi, M. Fabi, L. Ferraioli, V. Ferroni, N. Finetti, E. D. Fitzsimons, M. Freschi, L. Gesa, F. Gibert , et al. (60 additional authors not shown)

Abstract: Non-recurrent short term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n$^{-1}$ were observed between 2016 February 18 and 2017 July 3 aboard the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5$\times$10$^6$ km from Earth. The energy dependence of three Forbush decreases (FDs) is studied and reported here. A comparison of these obser… ▽ More Non-recurrent short term variations of the galactic cosmic-ray (GCR) flux above 70 MeV n$^{-1}$ were observed between 2016 February 18 and 2017 July 3 aboard the European Space Agency LISA Pathfinder (LPF) mission orbiting around the Lagrange point L1 at 1.5$\times$10$^6$ km from Earth. The energy dependence of three Forbush decreases (FDs) is studied and reported here. A comparison of these observations with others carried out in space down to the energy of a few tens of MeV n$^{-1}$ shows that the same GCR flux parameterization applies to events of different intensity during the main phase. FD observations in L1 with LPF and geomagnetic storm occurrence is also presented. Finally, the characteristics of GCR flux non-recurrent variations (peaks and depressions) of duration $<$ 2 days and their association with interplanetary structures are investigated. It is found that, most likely, plasma compression regions between subsequent corotating high-speed streams cause peaks, while heliospheric current sheet crossing cause the majority of the depressions. △ Less

Submitted 9 April, 2019; originally announced April 2019.

Journal ref: M. Armano et al 2019 ApJ 874 167