Showing 1–50 of 114 results for author: Anderson, K

Online regularization of Poincaré map of storage rings with Shannon entropy

Authors: Yongjun Li, Kelly Anderson, Derong Xu, Yue Hao, Kiman Ha, Yoshiteru Hidaka, Minghao Song, Robert Rainer, Victor Smaluk, Timur Shaftan

Abstract: Shannon Entropy is adopted to quantify the chaos of measured Poincaré maps in the National Synchrotron Light Source-II (NSLS-II) storage ring. The recurrent Poincaré maps, constructed from beam position monitor's turn-by-turn readings, are commonly used to observe the nonlinearity in ring-based accelerators. However, these observations typically only provide a qualitative observation. With some ca… ▽ More Shannon Entropy is adopted to quantify the chaos of measured Poincaré maps in the National Synchrotron Light Source-II (NSLS-II) storage ring. The recurrent Poincaré maps, constructed from beam position monitor's turn-by-turn readings, are commonly used to observe the nonlinearity in ring-based accelerators. However, these observations typically only provide a qualitative observation. With some canonical transformations on Poincaré maps, not only can the commonly used nonlinear characterizations be extracted, but more importantly, the chaos can be quantitatively measured with entropy. Entropy, therefore as a chaos indicator, is used for online Poincaré map regularization and dynamic aperture optimization in the NSLS-II ring. △ Less

Submitted 26 August, 2024; originally announced August 2024.

Comments: 5 pages, 6 figures

Cooling of gold cluster anions, Au$_N^-$, $N=2-13,15$, in a cryogenic ion-beam storage ring

Authors: Klavs Hansen, Tian Weihao, Emma K. Anderson, Mikael Björkhage, Henrik Cederquist, Ji MingChao, Stefan Rosén, Alice Schmidt-May, Mark H. Stockett, Henning Zettergren, Vitali Zhaunerchyk, Henning T. Schmidt

Abstract: We have measured the spontaneous and photo-induced decays of anionic gold clusters, Au$_N^-$, with sizes ranging from $N = 2$ to 13, and 15. After production in a sputter ion source, the size-selected clusters were stored in the cryogenic electrostatic ion-beam storage ring DESIREE and their neutralization decays were measured for storage times between 0.1 and 100 s. The dimer was observed to deca… ▽ More We have measured the spontaneous and photo-induced decays of anionic gold clusters, Au$_N^-$, with sizes ranging from $N = 2$ to 13, and 15. After production in a sputter ion source, the size-selected clusters were stored in the cryogenic electrostatic ion-beam storage ring DESIREE and their neutralization decays were measured for storage times between 0.1 and 100 s. The dimer was observed to decay by electron emission in parallel to neutral atom emission at long times, analogously to the behavior of copper and silver dimers, implying a breakdown of the Born-Oppenheimer approximation. Radiative cooling is observed for all cluster sizes except for the dimer. The decay rates of clusters $N=3,6,8-13,15$ show only a single radiative cooling time. For $N=6-13$ the cooling times have a strong odd-even oscillation with an amplitude that decrease with cluster size, and with the even $N$ having the fast cooling. We compare our results with previous measurements of radiative cooling rates of the corresponding cationic gold clusters, Au$_N^+$, which also show an odd-even effect with a similar oscillation amplitude but at orders of magnitude shorter time scales, and out of phase with the anions. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 12 pages, 12 figures

Optimization of a Welding Procedure for Making Critical Aluminum Welds on the LBNF Absorber Core Block

Authors: K. E. Anderson, A. Deshpande, V. I. Sidorov, J. Zahurones

Abstract: The LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. A weld quality optimization was performed to produce National Aeronautical Standard 1514 Class I quality welds on the aluminum core blocks. This was… ▽ More The LBNF Absorber consists of thirteen 6061-T6 aluminum core blocks. The core blocks are water cooled with de-ionized (DI) water which becomes radioactive during beam operations. The cooling water flows through gun-drilled channels in the core blocks. A weld quality optimization was performed to produce National Aeronautical Standard 1514 Class I quality welds on the aluminum core blocks. This was not successful in all cases. An existing Gas Tungsten Arc Welding Procedure Specification was fine tuned to minimize, in most cases, and eliminate detect-able tungsten inclusions in the welds. All the weld coupons, however passed welding inspection as per the piping code: ASME B31.3 Normal Fluid Service. Tungsten electrode diameter, type, and manufacturer were varied. Some of the samples were pre-heated and others were not. It was observed that larger diameter electrodes, 5/32 in., with pre-heated joints resulted in welds with the least number of tungsten inclusions. It is hypothesized that thinner electrodes breakdown easily and get lodged into the weld pool during the welding process. This breakdown is further enhanced by the large temperature differential between the un-preheated sample and the hot electrode. △ Less

Submitted 29 May, 2024; originally announced June 2024.

Comments: IPAC'24 - 15th International Particle Accelerator Conference

Report number: FERMILAB-CONF-23-798-AD-LBNF

Journal ref: JACoW IPAC2024 (2024) TUPS43

Thermalization and Criticality on an Analog-Digital Quantum Simulator

Authors: Trond I. Andersen, Nikita Astrakhantsev, Amir H. Karamlou, Julia Berndtsson, Johannes Motruk, Aaron Szasz, Jonathan A. Gross, Alexander Schuckert, Tom Westerhout, Yaxing Zhang, Ebrahim Forati, Dario Rossi, Bryce Kobrin, Agustin Di Paolo, Andrey R. Klots, Ilya Drozdov, Vladislav D. Kurilovich, Andre Petukhov, Lev B. Ioffe, Andreas Elben, Aniket Rath, Vittorio Vitale, Benoit Vermersch, Rajeev Acharya, Laleh Aghababaie Beni , et al. (202 additional authors not shown)

Abstract: Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators. Unlocking the full potential of such systems toward this goal requires flexible initial state preparation, precise time evolution, and extensive probes for final state characterization. We present a quantum simulator comprising 69 superconducting qubits which supports both universal qua… ▽ More Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators. Unlocking the full potential of such systems toward this goal requires flexible initial state preparation, precise time evolution, and extensive probes for final state characterization. We present a quantum simulator comprising 69 superconducting qubits which supports both universal quantum gates and high-fidelity analog evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. Emulating a two-dimensional (2D) XY quantum magnet, we leverage a wide range of measurement techniques to study quantum states after ramps from an antiferromagnetic initial state. We observe signatures of the classical Kosterlitz-Thouless phase transition, as well as strong deviations from Kibble-Zurek scaling predictions attributed to the interplay between quantum and classical coarsening of the correlated domains. This interpretation is corroborated by injecting variable energy density into the initial state, which enables studying the effects of the eigenstate thermalization hypothesis (ETH) in targeted parts of the eigenspectrum. Finally, we digitally prepare the system in pairwise-entangled dimer states and image the transport of energy and vorticity during thermalization. These results establish the efficacy of superconducting analog-digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics. △ Less

Submitted 8 July, 2024; v1 submitted 27 May, 2024; originally announced May 2024.

Growth of GeO2 on R-plane and C-plane Sapphires by MOCVD

Authors: Imteaz Rahaman, Hunter D. Ellis, Kathy Anderson, Michael A. Scarpulla, Kai Fu

Abstract: Rutile Germanium Dioxide (GeO2) has been recently theoretically identified as an ultrawide bandgap (UWBG) semiconductor with bandgap 4.68 eV similar to Ga2O3 but having bipolar dopability and ~2x higher electron mobility, Baliga figure of merit (BFOM) and thermal conductivity than Ga2O3. Bulk crystal growth is rapidly moving towards making large sized native substrates available. These outstanding… ▽ More Rutile Germanium Dioxide (GeO2) has been recently theoretically identified as an ultrawide bandgap (UWBG) semiconductor with bandgap 4.68 eV similar to Ga2O3 but having bipolar dopability and ~2x higher electron mobility, Baliga figure of merit (BFOM) and thermal conductivity than Ga2O3. Bulk crystal growth is rapidly moving towards making large sized native substrates available. These outstanding material properties position GeO2 as a highly attractive UWBG semiconductor for various applications. However, the epitaxial growth in the most advantageous polymorph (rutile), ensuring controlled phase, pristine surface/interface quality, precise microstructure, and optimal functional properties, is still in its infancy. In this work, we explored growth of GeO2 by metal-organic chemical vapor deposition (MOCVD) on both C- and R-plane sapphire. Utilizing tetramethylgermane (TMGe) as a precursor, we have investigated the influences of different parameters on the film properties, including growth temperature, chamber pressure, TMGe flow rate, oxygen flow rate, shroud gas flow rate, and rotation speed. The total pressure emerged as a crucial parameter while growth attempts at low total pressure resulted in no films for a wide range of temperatures, precursor flow rate, argon flow rates, and susceptor rotation rate. A phase diagram, derived from our experimental findings, delineates the growth windows for GeO2 films on sapphire substrates. This study serves as a pioneering guide for the MOCVD growth of GeO2 films. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 22 pages, 14 Figures

Prior Effective Sample Size When Borrowing on the Treatment Effect Scale

Authors: Hongtao Zhang, Keaven M Anderson, Zachary Zimmer, Gregory Golm, Aditi Sapre, Joseph G Ibrahim

Abstract: With the robust uptick in the applications of Bayesian external data borrowing, eliciting a prior distribution with the proper amount of information becomes increasingly critical. The prior effective sample size (ESS) is an intuitive and efficient measure for this purpose. The majority of ESS definitions have been proposed in the context of borrowing control information. While many Bayesian models… ▽ More With the robust uptick in the applications of Bayesian external data borrowing, eliciting a prior distribution with the proper amount of information becomes increasingly critical. The prior effective sample size (ESS) is an intuitive and efficient measure for this purpose. The majority of ESS definitions have been proposed in the context of borrowing control information. While many Bayesian models can be naturally extended to leveraging external information on the treatment effect scale, very little attention has been directed to computing the prior ESS in this setting. In this research, we bridge this methodological gap by extending the popular ELIR ESS definition. We lay out the general framework, and derive the prior ESS for various types of endpoints and treatment effect measures. The posterior distribution and the predictive consistency property of ESS are also examined. The methods are implemented in R programs available on GitHub: https://github.com/squallteo/TrtEffESS. △ Less

Submitted 20 April, 2024; originally announced April 2024.

Probing Electrical Properties of A Silicon Nanocrystal Thin Film Using X-ray Photoelectron Spectroscopy

Authors: Amrit Laudari, Sameera Pathiranage, Salim A. Thomas, Reed J. Petersen, Kenneth J. Anderson, Todd A. Pringle, Erik K. Hobbie, Nuri Oncel

Abstract: We performed X-ray photoelectron spectroscopy (XPS) measurements on a thin film of Si nanocrystals (SiNCs) while applying DC or AC external biases to extract the resistance and the capacitance of the thin film. The measurement consists of the application of 10 V DC or square wave pulses of 10 V amplitude to the sample at various frequencies ranging from 0.01 Hz to 1 MHz while recording X-ray photo… ▽ More We performed X-ray photoelectron spectroscopy (XPS) measurements on a thin film of Si nanocrystals (SiNCs) while applying DC or AC external biases to extract the resistance and the capacitance of the thin film. The measurement consists of the application of 10 V DC or square wave pulses of 10 V amplitude to the sample at various frequencies ranging from 0.01 Hz to 1 MHz while recording X-ray photoemission data. To analyze the data, we propose three different models with varying degrees of accuracy. The calculated capacitance of SiNCs agrees with the experimental value in the literature. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Journal ref: Rev. Sci. Instrum. 93, 083906 (2022)

Interspecific dispersal constraints suppress pattern formation in metacommunities

Authors: Patrick Lawton, Ashkaan K. Fahimipour, Kurt E. Anderson

Abstract: Decisions to disperse from a habitat stand out among organismal behaviors as pivotal drivers of ecosystem dynamics across scales. Encounters with other species are an important component of adaptive decision-making in dispersal, resulting in widespread behaviors like tracking resources or avoiding consumers in space. Despite this, metacommunity models often treat dispersal as a function of intrasp… ▽ More Decisions to disperse from a habitat stand out among organismal behaviors as pivotal drivers of ecosystem dynamics across scales. Encounters with other species are an important component of adaptive decision-making in dispersal, resulting in widespread behaviors like tracking resources or avoiding consumers in space. Despite this, metacommunity models often treat dispersal as a function of intraspecific density alone. We show, focusing initially on three-species network motifs, that interspecific dispersal rules generally drive a transition in metacommunities from homogeneous steady states to self-organized heterogeneous spatial patterns. However, when ecologically realistic constraints reflecting adaptive behaviors are imposed -- prey tracking and predator avoidance -- a pronounced homogenizing effect emerges where spatial pattern formation is suppressed. We demonstrate this effect for each motif by computing master stability functions that separate the contributions of local and spatial interactions to pattern formation. We extend this result to species rich food webs using a random matrix approach, where we find that eventually webs become large enough to override the homogenizing effect of adaptive dispersal behaviors, leading once again to predominately pattern forming dynamics. Our results emphasize the critical role of interspecific dispersal rules in shaping spatial patterns across landscapes, highlighting the need to incorporate adaptive behavioral constraints in efforts to link local species interactions and metacommunity structure. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Gemini 1.5: Unlocking multimodal understanding across millions of tokens of context

Authors: Gemini Team, Petko Georgiev, Ving Ian Lei, Ryan Burnell, Libin Bai, Anmol Gulati, Garrett Tanzer, Damien Vincent, Zhufeng Pan, Shibo Wang, Soroosh Mariooryad, Yifan Ding, Xinyang Geng, Fred Alcober, Roy Frostig, Mark Omernick, Lexi Walker, Cosmin Paduraru, Christina Sorokin, Andrea Tacchetti, Colin Gaffney, Samira Daruki, Olcan Sercinoglu, Zach Gleicher, Juliette Love , et al. (1110 additional authors not shown)

Abstract: In this report, we introduce the Gemini 1.5 family of models, representing the next generation of highly compute-efficient multimodal models capable of recalling and reasoning over fine-grained information from millions of tokens of context, including multiple long documents and hours of video and audio. The family includes two new models: (1) an updated Gemini 1.5 Pro, which exceeds the February… ▽ More In this report, we introduce the Gemini 1.5 family of models, representing the next generation of highly compute-efficient multimodal models capable of recalling and reasoning over fine-grained information from millions of tokens of context, including multiple long documents and hours of video and audio. The family includes two new models: (1) an updated Gemini 1.5 Pro, which exceeds the February version on the great majority of capabilities and benchmarks; (2) Gemini 1.5 Flash, a more lightweight variant designed for efficiency with minimal regression in quality. Gemini 1.5 models achieve near-perfect recall on long-context retrieval tasks across modalities, improve the state-of-the-art in long-document QA, long-video QA and long-context ASR, and match or surpass Gemini 1.0 Ultra's state-of-the-art performance across a broad set of benchmarks. Studying the limits of Gemini 1.5's long-context ability, we find continued improvement in next-token prediction and near-perfect retrieval (>99%) up to at least 10M tokens, a generational leap over existing models such as Claude 3.0 (200k) and GPT-4 Turbo (128k). Finally, we highlight real-world use cases, such as Gemini 1.5 collaborating with professionals on completing their tasks achieving 26 to 75% time savings across 10 different job categories, as well as surprising new capabilities of large language models at the frontier; when given a grammar manual for Kalamang, a language with fewer than 200 speakers worldwide, the model learns to translate English to Kalamang at a similar level to a person who learned from the same content. △ Less

Submitted 8 August, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Group Sequential Design Under Non-proportional Hazards

Authors: Yujie Zhao, Yilong Zhang, Larry Leon, Keaven M. Anderson

Abstract: Non-proportional hazards (NPH) are often observed in clinical trials with time-to-event endpoints. A common example is a long-term clinical trial with a delayed treatment effect in immunotherapy for cancer. When designing clinical trials with time-to-event endpoints, it is crucial to consider NPH scenarios to gain a complete understanding of design operating characteristics. In this paper, we focu… ▽ More Non-proportional hazards (NPH) are often observed in clinical trials with time-to-event endpoints. A common example is a long-term clinical trial with a delayed treatment effect in immunotherapy for cancer. When designing clinical trials with time-to-event endpoints, it is crucial to consider NPH scenarios to gain a complete understanding of design operating characteristics. In this paper, we focus on group sequential design for three NPH methods: the average hazard ratio, the weighted logrank test, and the MaxCombo combination test. For each of these approaches, we provide analytic forms of design characteristics that facilitate sample size calculation and bound derivation for group sequential designs. Examples are provided to illustrate the proposed methods. To facilitate statisticians in designing and comparing group sequential designs under NPH, we have implemented the group sequential design methodology in the gsDesign2 R package at https://cran.r-project.org/web/packages/gsDesign2/. △ Less

Submitted 4 December, 2023; originally announced December 2023.

AI in Pharma for Personalized Sequential Decision-Making: Methods, Applications and Opportunities

Authors: Yuhan Li, Hongtao Zhang, Keaven Anderson, Songzi Li, Ruoqing Zhu

Abstract: In the pharmaceutical industry, the use of artificial intelligence (AI) has seen consistent growth over the past decade. This rise is attributed to major advancements in statistical machine learning methodologies, computational capabilities and the increased availability of large datasets. AI techniques are applied throughout different stages of drug development, ranging from drug discovery to pos… ▽ More In the pharmaceutical industry, the use of artificial intelligence (AI) has seen consistent growth over the past decade. This rise is attributed to major advancements in statistical machine learning methodologies, computational capabilities and the increased availability of large datasets. AI techniques are applied throughout different stages of drug development, ranging from drug discovery to post-marketing benefit-risk assessment. Kolluri et al. provided a review of several case studies that span these stages, featuring key applications such as protein structure prediction, success probability estimation, subgroup identification, and AI-assisted clinical trial monitoring. From a regulatory standpoint, there was a notable uptick in submissions incorporating AI components in 2021. The most prevalent therapeutic areas leveraging AI were oncology (27%), psychiatry (15%), gastroenterology (12%), and neurology (11%). The paradigm of personalized or precision medicine has gained significant traction in recent research, partly due to advancements in AI techniques \cite{hamburg2010path}. This shift has had a transformative impact on the pharmaceutical industry. Departing from the traditional "one-size-fits-all" model, personalized medicine incorporates various individual factors, such as environmental conditions, lifestyle choices, and health histories, to formulate customized treatment plans. By utilizing sophisticated machine learning algorithms, clinicians and researchers are better equipped to make informed decisions in areas such as disease prevention, diagnosis, and treatment selection, thereby optimizing health outcomes for each individual. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Adjusted inference for multiple testing procedure in group sequential designs

Authors: Yujie Zhao, Qi Liu, Linda Z. Sun, Keaven M. Anderson

Abstract: Adjustment of statistical significance levels for repeated analysis in group sequential trials has been understood for some time. Similarly, methods for adjustment accounting for testing multiple hypotheses are common. There is limited research on simultaneously adjusting for both multiple hypothesis testing and multiple analyses of one or more hypotheses. We address this gap by proposing adjusted… ▽ More Adjustment of statistical significance levels for repeated analysis in group sequential trials has been understood for some time. Similarly, methods for adjustment accounting for testing multiple hypotheses are common. There is limited research on simultaneously adjusting for both multiple hypothesis testing and multiple analyses of one or more hypotheses. We address this gap by proposing adjusted-sequential p-values that reject an elementary hypothesis when its adjusted-sequential p-values are less than or equal to the family-wise Type I error rate (FWER) in a group sequential design. We also propose sequential p-values for intersection hypotheses as a tool to compute adjusted sequential p-values for elementary hypotheses. We demonstrate the application using weighted Bonferroni tests and weighted parametric tests, comparing adjusted sequential p-values to a desired FWER for inference on each elementary hypothesis tested. △ Less

Submitted 26 November, 2023; originally announced November 2023.

Dynamics of magnetization at infinite temperature in a Heisenberg spin chain

Authors: Eliott Rosenberg, Trond Andersen, Rhine Samajdar, Andre Petukhov, Jesse Hoke, Dmitry Abanin, Andreas Bengtsson, Ilya Drozdov, Catherine Erickson, Paul Klimov, Xiao Mi, Alexis Morvan, Matthew Neeley, Charles Neill, Rajeev Acharya, Richard Allen, Kyle Anderson, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Juan Atalaya, Joseph Bardin, A. Bilmes, Gina Bortoli , et al. (156 additional authors not shown)

Abstract: Understanding universal aspects of quantum dynamics is an unresolved problem in statistical mechanics. In particular, the spin dynamics of the 1D Heisenberg model were conjectured to belong to the Kardar-Parisi-Zhang (KPZ) universality class based on the scaling of the infinite-temperature spin-spin correlation function. In a chain of 46 superconducting qubits, we study the probability distributio… ▽ More Understanding universal aspects of quantum dynamics is an unresolved problem in statistical mechanics. In particular, the spin dynamics of the 1D Heisenberg model were conjectured to belong to the Kardar-Parisi-Zhang (KPZ) universality class based on the scaling of the infinite-temperature spin-spin correlation function. In a chain of 46 superconducting qubits, we study the probability distribution, $P(\mathcal{M})$, of the magnetization transferred across the chain's center. The first two moments of $P(\mathcal{M})$ show superdiffusive behavior, a hallmark of KPZ universality. However, the third and fourth moments rule out the KPZ conjecture and allow for evaluating other theories. Our results highlight the importance of studying higher moments in determining dynamic universality classes and provide key insights into universal behavior in quantum systems. △ Less

Submitted 4 April, 2024; v1 submitted 15 June, 2023; originally announced June 2023.

Journal ref: Science 384, 48-53 (2024)

The Spatial Correlation of High Mass X-ray Binaries and Young Star Clusters in Nearby Star-Forming Galaxies

Authors: Breanna A. Binder, Ashley K. Anderson, Kristen Garofali, Margaret Lazzarini, Benjamin F. Williams

Abstract: We present an analysis of the two-point spatial correlation functions of high-mass X-ray binary (HMXB) and young star cluster (YSC) populations in M31 and M33. We find evidence that HMXBs are spatially correlated with YSCs to a higher degree than would be expected from random chance in both galaxies. When supplemented with similar studies in the Milky Way, Small Magellanic Cloud, and NGC 4449, we… ▽ More We present an analysis of the two-point spatial correlation functions of high-mass X-ray binary (HMXB) and young star cluster (YSC) populations in M31 and M33. We find evidence that HMXBs are spatially correlated with YSCs to a higher degree than would be expected from random chance in both galaxies. When supplemented with similar studies in the Milky Way, Small Magellanic Cloud, and NGC 4449, we find that the peak value of the spatial correlation function correlates strongly with the specific star formation rate of the host galaxy. We additionally perform an X-ray stacking analysis of 211 non-X-ray detected YSCs in M31 and 463 YSCs in M33. We do not detect excess X-ray emission at the stacked cluster locations down to 3$σ$ upper limits of $\sim10^{33}$ erg s$^{-1}$ (0.35-8 keV) in both galaxies, which strongly suggests that dynamical formation within YSCs is not a major HMXB formation channel. We interpret our results in the context of (1) the recent star formation histories of the galaxies, which may produce differences in the demographics of compact objects powering the HMXBs, and (2) the differences in natal kicks experienced by compact objects during formation, which can eject newly-formed HMXB from their birth clusters. △ Less

Submitted 2 May, 2023; originally announced May 2023.

Comments: Accepted to MNRAS; 11 pages, 7 figures

Phase transition in Random Circuit Sampling

Authors: A. Morvan, B. Villalonga, X. Mi, S. Mandrà, A. Bengtsson, P. V. Klimov, Z. Chen, S. Hong, C. Erickson, I. K. Drozdov, J. Chau, G. Laun, R. Movassagh, A. Asfaw, L. T. A. N. Brandão, R. Peralta, D. Abanin, R. Acharya, R. Allen, T. I. Andersen, K. Anderson, M. Ansmann, F. Arute, K. Arya, J. Atalaya , et al. (160 additional authors not shown)

Abstract: Undesired coupling to the surrounding environment destroys long-range correlations on quantum processors and hinders the coherent evolution in the nominally available computational space. This incoherent noise is an outstanding challenge to fully leverage the computation power of near-term quantum processors. It has been shown that benchmarking Random Circuit Sampling (RCS) with Cross-Entropy Benc… ▽ More Undesired coupling to the surrounding environment destroys long-range correlations on quantum processors and hinders the coherent evolution in the nominally available computational space. This incoherent noise is an outstanding challenge to fully leverage the computation power of near-term quantum processors. It has been shown that benchmarking Random Circuit Sampling (RCS) with Cross-Entropy Benchmarking (XEB) can provide a reliable estimate of the effective size of the Hilbert space coherently available. The extent to which the presence of noise can trivialize the outputs of a given quantum algorithm, i.e. making it spoofable by a classical computation, is an unanswered question. Here, by implementing an RCS algorithm we demonstrate experimentally that there are two phase transitions observable with XEB, which we explain theoretically with a statistical model. The first is a dynamical transition as a function of the number of cycles and is the continuation of the anti-concentration point in the noiseless case. The second is a quantum phase transition controlled by the error per cycle; to identify it analytically and experimentally, we create a weak link model which allows varying the strength of noise versus coherent evolution. Furthermore, by presenting an RCS experiment with 67 qubits at 32 cycles, we demonstrate that the computational cost of our experiment is beyond the capabilities of existing classical supercomputers, even when accounting for the inevitable presence of noise. Our experimental and theoretical work establishes the existence of transitions to a stable computationally complex phase that is reachable with current quantum processors. △ Less

Submitted 21 December, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

Overcoming leakage in scalable quantum error correction

Authors: Kevin C. Miao, Matt McEwen, Juan Atalaya, Dvir Kafri, Leonid P. Pryadko, Andreas Bengtsson, Alex Opremcak, Kevin J. Satzinger, Zijun Chen, Paul V. Klimov, Chris Quintana, Rajeev Acharya, Kyle Anderson, Markus Ansmann, Frank Arute, Kunal Arya, Abraham Asfaw, Joseph C. Bardin, Alexandre Bourassa, Jenna Bovaird, Leon Brill, Bob B. Buckley, David A. Buell, Tim Burger, Brian Burkett , et al. (92 additional authors not shown)

Abstract: Leakage of quantum information out of computational states into higher energy states represents a major challenge in the pursuit of quantum error correction (QEC). In a QEC circuit, leakage builds over time and spreads through multi-qubit interactions. This leads to correlated errors that degrade the exponential suppression of logical error with scale, challenging the feasibility of QEC as a path… ▽ More Leakage of quantum information out of computational states into higher energy states represents a major challenge in the pursuit of quantum error correction (QEC). In a QEC circuit, leakage builds over time and spreads through multi-qubit interactions. This leads to correlated errors that degrade the exponential suppression of logical error with scale, challenging the feasibility of QEC as a path towards fault-tolerant quantum computation. Here, we demonstrate the execution of a distance-3 surface code and distance-21 bit-flip code on a Sycamore quantum processor where leakage is removed from all qubits in each cycle. This shortens the lifetime of leakage and curtails its ability to spread and induce correlated errors. We report a ten-fold reduction in steady-state leakage population on the data qubits encoding the logical state and an average leakage population of less than $1 \times 10^{-3}$ throughout the entire device. The leakage removal process itself efficiently returns leakage population back to the computational basis, and adding it to a code circuit prevents leakage from inducing correlated error across cycles, restoring a fundamental assumption of QEC. With this demonstration that leakage can be contained, we resolve a key challenge for practical QEC at scale. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: Main text: 7 pages, 5 figures

Deep Learning Object Detection Approaches to Signal Identification

Authors: Luke Wood, Kevin Anderson, Peter Gerstoft, Richard Bell, Raghab Subbaraman, Dinesh Bharadia

Abstract: Traditionally source identification is solved using threshold based energy detection algorithms. These algorithms frequently sum up the activity in regions, and consider regions above a specific activity threshold to be sources. While these algorithms work for the majority of cases, they often fail to detect signals that occupy small frequency bands, fail to distinguish sources with overlapping fr… ▽ More Traditionally source identification is solved using threshold based energy detection algorithms. These algorithms frequently sum up the activity in regions, and consider regions above a specific activity threshold to be sources. While these algorithms work for the majority of cases, they often fail to detect signals that occupy small frequency bands, fail to distinguish sources with overlapping frequency bands, and cannot detect any signals under a specified signal to noise ratio. Through the conversion of raw signal data to spectrogram, source identification can be framed as an object detection problem. By leveraging modern advancements in deep learning based object detection, we propose a system that manages to alleviate the failure cases encountered when using traditional source identification algorithms. Our contributions include framing source identification as an object detection problem, the publication of a spectrogram object detection dataset, and evaluation of the RetinaNet and YOLOv5 object detection models trained on the dataset. Our final models achieve Mean Average Precisions of up to 0.906. With such a high Mean Average Precision, these models are sufficiently robust for use in real world applications. △ Less

Submitted 1 November, 2022; v1 submitted 26 October, 2022; originally announced October 2022.

Baby Physical Safety Monitoring in Smart Home Using Action Recognition System

Authors: Victor Adewopo, Nelly Elsayed, Kelly Anderson

Abstract: Humans are able to intuitively deduce actions that took place between two states in observations via deductive reasoning. This is because the brain operates on a bidirectional communication model, which has radically improved the accuracy of recognition and prediction based on features connected to previous experiences. During the past decade, deep learning models for action recognition have signi… ▽ More Humans are able to intuitively deduce actions that took place between two states in observations via deductive reasoning. This is because the brain operates on a bidirectional communication model, which has radically improved the accuracy of recognition and prediction based on features connected to previous experiences. During the past decade, deep learning models for action recognition have significantly improved. However, deep neural networks struggle with these tasks on a smaller dataset for specific Action Recognition (AR) tasks. As with most action recognition tasks, the ambiguity of accurately describing activities in spatial-temporal data is a drawback that can be overcome by curating suitable datasets, including careful annotations and preprocessing of video data for analyzing various recognition tasks. In this study, we present a novel lightweight framework combining transfer learning techniques with a Conv2D LSTM layer to extract features from the pre-trained I3D model on the Kinetics dataset for a new AR task (Smart Baby Care) that requires a smaller dataset and less computational resources. Furthermore, we developed a benchmark dataset and an automated model that uses LSTM convolution with I3D (ConvLSTM-I3D) for recognizing and predicting baby activities in a smart baby room. Finally, we implemented video augmentation to improve model performance on the smart baby care task. Compared to other benchmark models, our experimental framework achieved better performance with less computational resources. △ Less

Submitted 29 April, 2023; v1 submitted 22 October, 2022; originally announced October 2022.

Purification-based quantum error mitigation of pair-correlated electron simulations

Authors: T. E. O'Brien, G. Anselmetti, F. Gkritsis, V. E. Elfving, S. Polla, W. J. Huggins, O. Oumarou, K. Kechedzhi, D. Abanin, R. Acharya, I. Aleiner, R. Allen, T. I. Andersen, K. Anderson, M. Ansmann, F. Arute, K. Arya, A. Asfaw, J. Atalaya, D. Bacon, J. C. Bardin, A. Bengtsson, S. Boixo, G. Bortoli, A. Bourassa , et al. (151 additional authors not shown)

Abstract: An important measure of the development of quantum computing platforms has been the simulation of increasingly complex physical systems. Prior to fault-tolerant quantum computing, robust error mitigation strategies are necessary to continue this growth. Here, we study physical simulation within the seniority-zero electron pairing subspace, which affords both a computational stepping stone to a ful… ▽ More An important measure of the development of quantum computing platforms has been the simulation of increasingly complex physical systems. Prior to fault-tolerant quantum computing, robust error mitigation strategies are necessary to continue this growth. Here, we study physical simulation within the seniority-zero electron pairing subspace, which affords both a computational stepping stone to a fully correlated model, and an opportunity to validate recently introduced ``purification-based'' error-mitigation strategies. We compare the performance of error mitigation based on doubling quantum resources in time (echo verification) or in space (virtual distillation), on up to $20$ qubits of a superconducting qubit quantum processor. We observe a reduction of error by one to two orders of magnitude below less sophisticated techniques (e.g. post-selection); the gain from error mitigation is seen to increase with the system size. Employing these error mitigation strategies enables the implementation of the largest variational algorithm for a correlated chemistry system to-date. Extrapolating performance from these results allows us to estimate minimum requirements for a beyond-classical simulation of electronic structure. We find that, despite the impressive gains from purification-based error mitigation, significant hardware improvements will be required for classically intractable variational chemistry simulations. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 10 pages, 13 page supplementary material, 12 figures. Experimental data available at https://doi.org/10.5281/zenodo.7225821

Journal ref: Nat. Phys. (2023)

Clustering Egocentric Images in Passive Dietary Monitoring with Self-Supervised Learning

Authors: Jiachuan Peng, Peilun Shi, Jianing Qiu, Xinwei Ju, Frank P. -W. Lo, Xiao Gu, Wenyan Jia, Tom Baranowski, Matilda Steiner-Asiedu, Alex K. Anderson, Megan A McCrory, Edward Sazonov, Mingui Sun, Gary Frost, Benny Lo

Abstract: In our recent dietary assessment field studies on passive dietary monitoring in Ghana, we have collected over 250k in-the-wild images. The dataset is an ongoing effort to facilitate accurate measurement of individual food and nutrient intake in low and middle income countries with passive monitoring camera technologies. The current dataset involves 20 households (74 subjects) from both the rural a… ▽ More In our recent dietary assessment field studies on passive dietary monitoring in Ghana, we have collected over 250k in-the-wild images. The dataset is an ongoing effort to facilitate accurate measurement of individual food and nutrient intake in low and middle income countries with passive monitoring camera technologies. The current dataset involves 20 households (74 subjects) from both the rural and urban regions of Ghana, and two different types of wearable cameras were used in the studies. Once initiated, wearable cameras continuously capture subjects' activities, which yield massive amounts of data to be cleaned and annotated before analysis is conducted. To ease the data post-processing and annotation tasks, we propose a novel self-supervised learning framework to cluster the large volume of egocentric images into separate events. Each event consists of a sequence of temporally continuous and contextually similar images. By clustering images into separate events, annotators and dietitians can examine and analyze the data more efficiently and facilitate the subsequent dietary assessment processes. Validated on a held-out test set with ground truth labels, the proposed framework outperforms baselines in terms of clustering quality and classification accuracy. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: accepted to BHI 2022

A gated group sequential design for seamless Phase II/III trial with subpopulation selection

Authors: Guanhong Miao, Jason J. Z. Liao, Jing Yang, Keaven Anderson

Abstract: Due to the high cost and high failure rate of Phase III trials, seamless Phase II/III designs are more and more popular to trial efficiency. A potential attraction of Phase II/III design is to allow a randomized proof-of-concept stage prior to committing to the full cost of the Phase III trial. Population selection during the trial allows a trial to adapt and focus investment where it is most like… ▽ More Due to the high cost and high failure rate of Phase III trials, seamless Phase II/III designs are more and more popular to trial efficiency. A potential attraction of Phase II/III design is to allow a randomized proof-of-concept stage prior to committing to the full cost of the Phase III trial. Population selection during the trial allows a trial to adapt and focus investment where it is most likely to provide patient benefit. Motivated by a clinical trial to find the population that potential benefits with dual-primary endpoints progression free survival (PFS) and overall survival (OS), we propose a gated group sequential design for a seamless Phase II/III trial design with population selection. The investigated design controls the familywise error rate and allows multiple interim analyses to enable early stopping for efficacy or futility. Simulations and an illustrative example suggest that the proposed gated group sequential design can have more power than the commonly used classical group sequential design, and reduces the patient's exposure to less effective treatment if the complementary sub-group has less significant treatment effect. The proposed design has the potential to save drug development cost and more quickly fulfill unmet medical needs. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 3 figures, 5 tables

MSC Class: 62P10

Dynamics of a Rotated Orthogonal Gravitational Wedge Billiard

Authors: K. D. Anderson

Abstract: We investigate a rotated, orthogonal gravitational wedge billiard - a special case of the asymmetric wedge billiard - in which the dynamics are integrable. We derive equations and conditions under which periodic orbits may be constructed for this model, and show that any other trajectory will be dense in the configuration space. We investigate a rotated, orthogonal gravitational wedge billiard - a special case of the asymmetric wedge billiard - in which the dynamics are integrable. We derive equations and conditions under which periodic orbits may be constructed for this model, and show that any other trajectory will be dense in the configuration space. △ Less

Submitted 10 June, 2022; originally announced June 2022.

Comments: 13 pages, 7 figures

MSC Class: 70F99

Journal ref: International Journal of Bifurcation and Chaos, vol 33, no 12 (2023)

Were the Obliquities in DI Herculis Excited by an Unseen Tertiary Companion?

Authors: Kassandra R. Anderson, Joshua N. Winn

Abstract: The eclipsing binary DI Herculis garnered interest for several decades because of an apparent disagreement between the observed and calculated values of the apsidal precession rate. The problem was resolved when both stars were found to have high obliquities, but the reason for the high obliquities is unknown. Here we investigate the possibility that the obliquities are (or were) excited by an uns… ▽ More The eclipsing binary DI Herculis garnered interest for several decades because of an apparent disagreement between the observed and calculated values of the apsidal precession rate. The problem was resolved when both stars were found to have high obliquities, but the reason for the high obliquities is unknown. Here we investigate the possibility that the obliquities are (or were) excited by an unseen tertiary star. Obliquity excitation in the current orbital configuration can be ruled out with existing data; any tertiary star that is sufficiently close or massive to overcome the strong spin-orbit coupling of the binary would have been detected through various dynamical effects. It remains possible that the orbit of DI Herculis was initially wider and the obliquity was excited during high-eccentricity tidal migration driven by a tertiary companion, but in this scenario it would be difficult to explain why the observed spin rates are much faster than the pseudo-synchronous rate. In addition, inward migration is most likely to arise when the mass of the perturbing star is comparable to the binary mass, and such a bright tertiary would have been detected in imaging or spectroscopic data. Alternative explanations that do not invoke a tertiary star should be sought for the large obliquities in DI Herculis. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: 13 pages, 6 figures. Accepted for publication in ApJ

Sloppy model analysis identifies bifurcation parameters without Normal Form analysis

Authors: Christian N. K. Anderson, Mark K. Transtrum

Abstract: Bifurcation phenomena are common in multi-dimensional multi-parameter dynamical systems. Normal form theory suggests that the bifurcations themselves are driven by relatively few parameters; however, these are often nonlinear combinations of the bare parameters in which the equations are expressed. Discovering reparameterizations to transform such complex original equations into normal-form is oft… ▽ More Bifurcation phenomena are common in multi-dimensional multi-parameter dynamical systems. Normal form theory suggests that the bifurcations themselves are driven by relatively few parameters; however, these are often nonlinear combinations of the bare parameters in which the equations are expressed. Discovering reparameterizations to transform such complex original equations into normal-form is often very difficult, and the reparameterization may not even exist in a closed-form. Recent advancements have tied both information geometry and bifurcations to the Renormalization Group. Here, we show that sloppy model analysis (a method of information geometry) can be used directly on bifurcations of increasing time scales to rapidly characterize the system's topological inhomogeneities, whether the system is in normal form or not. We anticipate that this novel analytical method, which we call time-widening information geometry (TWIG), will be useful in applied network analysis. △ Less

Submitted 27 November, 2023; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: 18 pages, 11 postscript figures

MSC Class: 34C23; 37G05; 53B12 ACM Class: J.2.8

Signatures of a Distant Planet on the Inclination Distribution of the Detached Kuiper Belt

Authors: Kalee E. Anderson, Nathan A. Kaib

Abstract: A distant, massive planet in the outer solar system has recently been proposed to explain some observed features of extreme trans-Neptunian objects (TNOs). Here we use N-body simulations of the formation of the Kuiper belt and Oort cloud as well as a survey simulator to compare models of the solar system with and without a 9th planet to one another as well as to observations. The main mechanism fo… ▽ More A distant, massive planet in the outer solar system has recently been proposed to explain some observed features of extreme trans-Neptunian objects (TNOs). Here we use N-body simulations of the formation of the Kuiper belt and Oort cloud as well as a survey simulator to compare models of the solar system with and without a 9th planet to one another as well as to observations. The main mechanism for TNOs to be deposited into the distant ($a$ > au), detached ($q$ > au) region of the Kuiper Belt in the 8-planet model is Kozai-Lidov oscillation of objects in mean motion resonances (MMR) with Neptune. This effect does not deposit low-inclination ($i \lesssim$ 20°) objects into this region. However, we find that the 9th planet generates a group of distant, detached TNOs at low inclinations that are not present in the 8-planet model. This disparity between the 8-planet and the 9-planet models could provide a strong constraint on a possible planet 9 with further detections of TNOs in the distant, detached region of the Kuiper Belt. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 8 pages, 3 figures, accepted to ApJL

Challenges in Detoxifying Language Models

Authors: Johannes Welbl, Amelia Glaese, Jonathan Uesato, Sumanth Dathathri, John Mellor, Lisa Anne Hendricks, Kirsty Anderson, Pushmeet Kohli, Ben Coppin, Po-Sen Huang

Abstract: Large language models (LM) generate remarkably fluent text and can be efficiently adapted across NLP tasks. Measuring and guaranteeing the quality of generated text in terms of safety is imperative for deploying LMs in the real world; to this end, prior work often relies on automatic evaluation of LM toxicity. We critically discuss this approach, evaluate several toxicity mitigation strategies wit… ▽ More Large language models (LM) generate remarkably fluent text and can be efficiently adapted across NLP tasks. Measuring and guaranteeing the quality of generated text in terms of safety is imperative for deploying LMs in the real world; to this end, prior work often relies on automatic evaluation of LM toxicity. We critically discuss this approach, evaluate several toxicity mitigation strategies with respect to both automatic and human evaluation, and analyze consequences of toxicity mitigation in terms of model bias and LM quality. We demonstrate that while basic intervention strategies can effectively optimize previously established automatic metrics on the RealToxicityPrompts dataset, this comes at the cost of reduced LM coverage for both texts about, and dialects of, marginalized groups. Additionally, we find that human raters often disagree with high automatic toxicity scores after strong toxicity reduction interventions -- highlighting further the nuances involved in careful evaluation of LM toxicity. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 23 pages, 6 figures, published in Findings of EMNLP 2021

ACM Class: I.2.6; I.2.7

Further Evidence for Tidal Spin-Up of Hot Jupiter Host Stars

Authors: Roberto A. Tejada Arevalo, Joshua N. Winn, Kassandra R. Anderson

Abstract: For most hot Jupiters around main-sequence Sun-like stars, tidal torques are expected to transfer angular momentum from the planet's orbit to the star's rotation. The timescale for this process is difficult to calculate, leading to uncertainties in the history of orbital evolution of hot Jupiters. We present evidence for tidal spin-up by taking advantage of recent advances in planet detection and… ▽ More For most hot Jupiters around main-sequence Sun-like stars, tidal torques are expected to transfer angular momentum from the planet's orbit to the star's rotation. The timescale for this process is difficult to calculate, leading to uncertainties in the history of orbital evolution of hot Jupiters. We present evidence for tidal spin-up by taking advantage of recent advances in planet detection and host-star characterization. We compared the projected rotation velocities and rotation periods of Sun-like stars with hot Jupiters and spectroscopically similar stars with (i) wider-orbiting giant planets, and (ii) less massive planets. The hot Jupiter hosts tend to spin faster than the stars in either of the control samples. Reinforcing earlier studies, the results imply that hot Jupiters alter the spins of their host stars while they are on the main sequence, and that the ages of hot-Jupiter hosts cannot be reliably determined using gyrochronology. △ Less

Submitted 12 July, 2021; originally announced July 2021.

Egocentric Image Captioning for Privacy-Preserved Passive Dietary Intake Monitoring

Authors: Jianing Qiu, Frank P. -W. Lo, Xiao Gu, Modou L. Jobarteh, Wenyan Jia, Tom Baranowski, Matilda Steiner-Asiedu, Alex K. Anderson, Megan A McCrory, Edward Sazonov, Mingui Sun, Gary Frost, Benny Lo

Abstract: Camera-based passive dietary intake monitoring is able to continuously capture the eating episodes of a subject, recording rich visual information, such as the type and volume of food being consumed, as well as the eating behaviours of the subject. However, there currently is no method that is able to incorporate these visual clues and provide a comprehensive context of dietary intake from passive… ▽ More Camera-based passive dietary intake monitoring is able to continuously capture the eating episodes of a subject, recording rich visual information, such as the type and volume of food being consumed, as well as the eating behaviours of the subject. However, there currently is no method that is able to incorporate these visual clues and provide a comprehensive context of dietary intake from passive recording (e.g., is the subject sharing food with others, what food the subject is eating, and how much food is left in the bowl). On the other hand, privacy is a major concern while egocentric wearable cameras are used for capturing. In this paper, we propose a privacy-preserved secure solution (i.e., egocentric image captioning) for dietary assessment with passive monitoring, which unifies food recognition, volume estimation, and scene understanding. By converting images into rich text descriptions, nutritionists can assess individual dietary intake based on the captions instead of the original images, reducing the risk of privacy leakage from images. To this end, an egocentric dietary image captioning dataset has been built, which consists of in-the-wild images captured by head-worn and chest-worn cameras in field studies in Ghana. A novel transformer-based architecture is designed to caption egocentric dietary images. Comprehensive experiments have been conducted to evaluate the effectiveness and to justify the design of the proposed architecture for egocentric dietary image captioning. To the best of our knowledge, this is the first work that applies image captioning for dietary intake assessment in real life settings. △ Less

Submitted 1 March, 2023; v1 submitted 1 July, 2021; originally announced July 2021.

Journal ref: IEEE Transactions on Cybernetics, 2023

Microgrid Resilience: A Holistic and Context-Aware Resilience Metric

Authors: Sakshi Mishra, Ted Kwasnik, Kate Anderson

Abstract: Microgrids present an effective solution for the coordinated deployment of various distributed energy resources and furthermore provide myriad additional benefits such as resilience, decreased carbon footprint, and reliability to energy consumers and the energy system as a whole. Boosting the resilience of distribution systems is another major benefit of microgrids. This is because they can also s… ▽ More Microgrids present an effective solution for the coordinated deployment of various distributed energy resources and furthermore provide myriad additional benefits such as resilience, decreased carbon footprint, and reliability to energy consumers and the energy system as a whole. Boosting the resilience of distribution systems is another major benefit of microgrids. This is because they can also serve as a backup power source when the utility grid operations are interrupted due to either high-probability low-impact events like a component failure or low-probability high-impact events - be it a natural disaster or a planned cyberattack. However, the degree to which any particular system can defend, adapt, and restore normal operation depends on various factors including the type and severity of events to which a microgrid is subjected. These factors, in turn, are dependent on the geographical location of the deployed microgrid as well as the cyber risk profile of the site where the microgrid is operating. Therefore, in this work, we attempt to capture this multi-dimensional interplay of various factors in quantifying the ability of the microgrid to be resilient in these varying aspects. This paper, thus, proposes a customized site-specific quantification of the resilience strength for the individual microgrid capability to absorb, restore, and adapt to the changing circumstances for sustaining the critical load when a low-probability high-impact event occurs - termed as - context-aware resilience metric. We also present a case study to illustrate the key elements of our integrated analytical approach. △ Less

Submitted 31 July, 2021; v1 submitted 17 June, 2021; originally announced June 2021.

Comments: 14 pages, 5 figures, in journal review

A unified framework for weighted parametric group sequential design (WPGSD)

Authors: Keaven M. Anderson, Zifang Guo, Jing Zhao, Linda Z. Sun

Abstract: Group sequential design (GSD) is widely used in clinical trials in which correlated tests of multiple hypotheses are used. Multiple primary objectives resulting in tests with known correlations include evaluating 1) multiple experimental treatment arms, 2) multiple populations, 3) the combination of multiple arms and multiple populations, or 4) any asymptotically multivariate normal tests. In this… ▽ More Group sequential design (GSD) is widely used in clinical trials in which correlated tests of multiple hypotheses are used. Multiple primary objectives resulting in tests with known correlations include evaluating 1) multiple experimental treatment arms, 2) multiple populations, 3) the combination of multiple arms and multiple populations, or 4) any asymptotically multivariate normal tests. In this paper, we focus on the first 3 of these and extend the framework of the weighted parametric multiple test procedure from fixed designs with a single analysis per objective to a GSD setting where different objectives may be assessed at the same or different times, each in a group sequential fashion. Pragmatic methods for design and analysis of weighted parametric group sequential design(WPGSD) under closed testing procedures are proposed to maintain the strong control of familywise Type I error rate (FWER) when correlations between tests are incorporated. This results in the ability to relax testing bounds compared to designs not fully adjusting for known correlations, increasing power or allowing decreased sample size. We illustrate the proposed methods using clinical trial examples and conduct a simulation study to evaluate the operating characteristics. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Comments: 24 pages (37 with appendices); 2 figures (+ 2 in appendices); submitted for publication

MSC Class: 62N03

A Far Ultraviolet Spectroscopic Explorer Survey of Interstellar Molecular Hydrogen in the Galactic Disk

Authors: J. Michael Shull, Charles W. Danforth, Katherine L. Anderson

Abstract: We report results from a FUSE survey of interstellar molecular hydrogen (H2) in the Galactic disk toward 139 O-type and early B-type stars at Galactic latitudes $|b| < 10^{\circ}$, with updated photometric and parallax distances. The H2 absorption is measured using the far-ultraviolet Lyman and Werner bands, including strong R(0), R(1), and P(1) lines from rotational levels $J = 0$ and $J = 1$ and… ▽ More We report results from a FUSE survey of interstellar molecular hydrogen (H2) in the Galactic disk toward 139 O-type and early B-type stars at Galactic latitudes $|b| < 10^{\circ}$, with updated photometric and parallax distances. The H2 absorption is measured using the far-ultraviolet Lyman and Werner bands, including strong R(0), R(1), and P(1) lines from rotational levels $J = 0$ and $J = 1$ and excited states up to $J = 5$ (sometimes $J = 6$ and 7). For each sight line, we report column densities $N_{H2}$, $N_{HI}$, $N(J)$, $N_H = N_{HI} + 2N_{H2}$, and molecular fraction, $f_{H2} = 2N_{H2}/N_H$. Our survey extends the 1977 Copernicus H2 survey up to $N_H \sim 5\times10^{21}$ cm$^{-2}$. The lowest rotational states have mean excitation temperatures and rms dispersions, $T_{01} = 88\pm 20$ K and $T_{02} = 77\pm18$ K, suggesting that J = 0,1,2 are coupled to the gas kinetic temperature. Populations of higher-J states exhibit mean excitation temperatures, $T_{24} = 237\pm91$ K and $T_{35} = 304\pm108$ K, produced primarily by UV radiative pumping. Correlations of $f_{H2}$ with E(B-V) and N_H show a transition to $f_{H2} \geq 0.1$ at $N_ H \geq 10^{21}$ cm$^{-2}$ and $E(B-V) > 0.2$, interpreted with an analytic model of H2 formation-dissociation equilibrium and attenuation of the far-UV radiation field by self-shielding and dust opacity. Results of this disk survey are compared to previous FUSE studies of H2 in translucent clouds, at high Galactic latitudes, and in the Magellanic Clouds. Using updated distances to the target stars, we find average sight-line values $\langle f_{H2} \rangle \geq 0.20$ and $\langle N_H/E(B-V) \rangle = (6.07\pm1.01)\times10^{21}$ cm$^{-2}$ mag$^{-1}$. △ Less

Submitted 5 March, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: Accepted to ApJ. 40 pages, 11 figures, 5 tables (four will be machine readable in journal). Revised text in Sections 3.4 and 3.5

On a Possible Solution to the Tidal Realignment Problem for Hot Jupiters

Authors: Kassandra R. Anderson, Joshua N. Winn, Kaloyan Penev

Abstract: Hot stars with hot Jupiters have a wide range of obliquities, while cool stars with hot Jupiters tend to have low obliquities. An enticing explanation for this pattern is tidal realignment of the cool host stars, although this explanation assumes that obliquity damping occurs faster than orbital decay, an assumption that needs further exploration. Here we revisit this tidal realignment problem, bu… ▽ More Hot stars with hot Jupiters have a wide range of obliquities, while cool stars with hot Jupiters tend to have low obliquities. An enticing explanation for this pattern is tidal realignment of the cool host stars, although this explanation assumes that obliquity damping occurs faster than orbital decay, an assumption that needs further exploration. Here we revisit this tidal realignment problem, building on previous work identifying a low-frequency component of the time-variable tidal potential that affects the obliquity but not the orbital separation. We adopt a recent empirically-based model for the stellar tidal quality factor and its sharp increase with forcing frequency. This leads to enhanced dissipation at low frequencies, and efficient obliquity damping. We model the tidal evolution of 46 observed hot Jupiters orbiting cool stars. A key parameter is the stellar age, which we determine in a homogeneous manner for the sample, taking advantage of Gaia DR2 data. We explore a variety of tidal histories and futures for each system, finding in most cases that the stellar obliquity is successfully damped before the planet is destroyed. A testable prediction of our model is that hot-Jupiter hosts with orbital periods shorter than 2--3 days should have obliquities much smaller than $1^\circ$. With the possible exception of WASP-19b, the predicted future lifetimes of the planets range from $10^8$\,yr to more than $10^{10}$\,yr. Thus, our model implies that these hot Jupiters are probably not in immediate danger of being devoured by their host stars while they are on the main sequence. △ Less

Submitted 15 April, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 24 pages, 14 figures, 1 table, accepted for publication in ApJ. Revised version following the referee report

Using Unity to Help Solve Intelligence

Authors: Tom Ward, Andrew Bolt, Nik Hemmings, Simon Carter, Manuel Sanchez, Ricardo Barreira, Seb Noury, Keith Anderson, Jay Lemmon, Jonathan Coe, Piotr Trochim, Tom Handley, Adrian Bolton

Abstract: In the pursuit of artificial general intelligence, our most significant measurement of progress is an agent's ability to achieve goals in a wide range of environments. Existing platforms for constructing such environments are typically constrained by the technologies they are founded on, and are therefore only able to provide a subset of scenarios necessary to evaluate progress. To overcome these… ▽ More In the pursuit of artificial general intelligence, our most significant measurement of progress is an agent's ability to achieve goals in a wide range of environments. Existing platforms for constructing such environments are typically constrained by the technologies they are founded on, and are therefore only able to provide a subset of scenarios necessary to evaluate progress. To overcome these shortcomings, we present our use of Unity, a widely recognized and comprehensive game engine, to create more diverse, complex, virtual simulations. We describe the concepts and components developed to simplify the authoring of these environments, intended for use predominantly in the field of reinforcement learning. We also introduce a practical approach to packaging and re-distributing environments in a way that attempts to improve the robustness and reproducibility of experiment results. To illustrate the versatility of our use of Unity compared to other solutions, we highlight environments already created using our approach from published papers. We hope that others can draw inspiration from how we adapted Unity to our needs, and anticipate increasingly varied and complex environments to emerge from our approach as familiarity grows. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Excitation of Spin-Orbit Misalignments in Stellar Binaries with Circumbinary Disks: Application to DI Herculis

Authors: Kassandra R. Anderson, Dong Lai

Abstract: The large spin-orbit misalignments in the DI Herculis stellar binary system have resolved the decades-long puzzle of the anomalously slow apsidal precession rate, but raise new questions regarding the origin of the obliquities. This paper investigates obliquity evolution in stellar binaries hosting modestly-inclined circumbinary disks. As the disk and binary axes undergo mutual precession, each ob… ▽ More The large spin-orbit misalignments in the DI Herculis stellar binary system have resolved the decades-long puzzle of the anomalously slow apsidal precession rate, but raise new questions regarding the origin of the obliquities. This paper investigates obliquity evolution in stellar binaries hosting modestly-inclined circumbinary disks. As the disk and binary axes undergo mutual precession, each oblate star experiences a torque from its companion star, so that the spin and orbital axes undergo mutual precession. As the disk loses mass through a combination of winds and accretion, the system may be captured into a high-obliquity Cassini state (a spin-orbit resonance). The final obliquity depends on the details of the disk dispersal. We construct a simple disk model to emulate disk dispersal due to viscous accretion and photoevaporation, and identify the necessary disk properties for producing the observed obliquities in DI Herculis. The disk must be massive (at least $10 \%$ of the binary mass). If accretion onto the binary is suppressed, the observed high stellar obliquities are reproduced with a binary-disk inclination of $\sim 5^\circ - 10^\circ$, but if substantial accretion occurs, the inclination must be larger, $\sim 20^\circ - 30^\circ$. If moderate accretion occurs, initially the disk must lose its mass slowly, but eventually lose its remaining mass abruptly, analogous to the observed two-timescale behavior for disks around T-Tauri stars. The spin feedback on the binary orbit causes the binary-disk inclination to decay as the obliquity evolves, a feature that is absent from the standard Cassini state treatment. △ Less

Submitted 24 November, 2020; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 13 pages, 6 figures, accepted for publication in ApJ

Computational Framework for Behind-The-Meter DER Techno-Economic Modeling and Optimization -- REopt Lite

Authors: Sakshi Mishra, Josiah Pohl, Nick Laws, Dylan Cutler, Ted Kwasnik, William Becker, Alex Zolan, Kate Anderson, Dan Olis, Emma Elgqvist

Abstract: The global energy system is undergoing a major transformation. Renewable energy generation is growing and is projected to accelerate further with the global emphasis on decarbonization. Furthermore, distributed generation is projected to play a significant role in the new energy system, and energy models are playing a key role in understanding how distributed generation can be integrated reliably… ▽ More The global energy system is undergoing a major transformation. Renewable energy generation is growing and is projected to accelerate further with the global emphasis on decarbonization. Furthermore, distributed generation is projected to play a significant role in the new energy system, and energy models are playing a key role in understanding how distributed generation can be integrated reliably and economically. The deployment of massive amounts of distributed generation requires understanding the interface of technology, economics, and policy in the energy modeling process. In this work, we present an end-to-end computational framework for distributed energy resource (DER) modeling, REopt Lite which addresses this need effectively. We describe the problem space, the building blocks of the model, the scaling capabilities of the design, the optimization formulation, and the accessibility of the model. We present a framework for accelerating the techno-economic analysis of behind-the-meter distributed energy resources to enable rapid planning and decision-making, thereby significantly boosting the rate the renewable energy deployment. Lastly, but equally importantly, this computation framework is open-sourced to facilitate transparency, flexibility, and wider collaboration opportunities within the worldwide energy modeling community. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 18 pages, 6 figures, under journal review

Probabilistic Resilience of DER systems -- A Simulation Assisted Optimization Approach

Authors: Sakshi Mishra, Kate Anderson

Abstract: Energy systems resilience is becoming increasingly important as the frequency of major grid outages increases. In this work, we present a methodology to optimize a behind-the-meter distributed energy resource system to sustain a site's critical loads during a pre-defined utility outage period. With the fixed system design, we then propose an outage simulation approach to estimate the resilience po… ▽ More Energy systems resilience is becoming increasingly important as the frequency of major grid outages increases. In this work, we present a methodology to optimize a behind-the-meter distributed energy resource system to sustain a site's critical loads during a pre-defined utility outage period. With the fixed system design, we then propose an outage simulation approach to estimate the resilience potential of the DER system to sustain loads beyond the fixed outage period - a yearlong resilience performance analysis. We apply statistical analysis to assess the system's resilience performance over a broader parametric problem space on an hourly, monthly, and yearly basis. We present a case study to demonstrate the impact of the pre-defined outage period on the resilience performance of the system. The case study shows that the probability of surviving a random outage of a given duration changes from 20 percent to 95 percent when the outage is modeled for a weekday instead of a weekend for the case study building. △ Less

Submitted 1 December, 2020; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: 5 pages, 6 figures, conference paper under review

Giant Planet Scatterings and Collisions: Hydrodynamics, Merger-Ejection Branching Ratio, and Properties of the Remnants

Authors: Jiaru Li, Dong Lai, Kassandra R. Anderson, Bonan Pu

Abstract: Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple "sticky-sphere" prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and $N$-body integrations to investi… ▽ More Planetary systems with sufficiently small orbital spacings can experience planetary mergers and ejections. The branching ratio of mergers vs ejections depends sensitively on the treatment of planetary close encounters. Previous works have adopted a simple "sticky-sphere" prescription, whose validity is questionable. We apply both smoothed particle hydrodynamics and $N$-body integrations to investigate the fluid effects in close encounters between gas giants and the long-term evolution of closely-packed planetary systems. Focusing on parabolic encounters between Jupiter-like planets with $M_J$ and $2M_J$, we find that quick mergers occur when the impact parameter $r_p$ (the pericenter separation between the planets) is less than $2R_J$, and the merger conserved 97% of the initial mass. Strong tidal effects can affect the "binary-planet" orbit when $r_p$ is between $2R_J$ and $4R_J$. We quantify these effects using a set of fitting formulae that can be implemented in $N$-body codes. We run a suite of $N$-body simulations with and without the formulae for systems of two giant planets initially in unstable, nearly circular coplanar orbits. The fluid (tidal) effects significantly increase the branching ratio of planetary mergers relative to ejections by doubling the effective collision radius. While the fluid effects do not change the distributions of semi-major axis and eccentricity of each type of remnant planets (mergers vs surviving planets in ejections), the overall orbital properties of planet scattering remnants are strongly affected due to the change in the branching ratio. We also find that the merger products have broad distributions of spin magnitudes and obliquities. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Comments: 13 pages, 14 Figures, 4 Tables, Submitted to MNRAS

New high-sensitivity searches for neutrons converting into antineutrons and/or sterile neutrons at the European Spallation Source

Authors: A. Addazi, K. Anderson, S. Ansell, K. Babu, J. Barrow, D. V. Baxter, P. M. Bentley, Z. Berezhiani, R. Bevilacqua, C. Bohm, G. Brooijmans, J. Broussard, R. Biondi, B. Dev, C. Crawford, A. Dolgov, K. Dunne, P. Fierlinger, M. R. Fitzsimmons, A. Fomin, M. Frost, S. Gardner, A. Galindo-Uribarri, E. Golubeva, S. Girmohanta , et al. (70 additional authors not shown)

Abstract: The violation of Baryon Number, $\mathcal{B}$, is an essential ingredient for the preferential creation of matter over antimatter needed to account for the observed baryon asymmetry in the universe. However, such a process has yet to be experimentally observed. The HIBEAM/NNBAR %experiment program is a proposed two-stage experiment at the European Spallation Source (ESS) to search for baryon numbe… ▽ More The violation of Baryon Number, $\mathcal{B}$, is an essential ingredient for the preferential creation of matter over antimatter needed to account for the observed baryon asymmetry in the universe. However, such a process has yet to be experimentally observed. The HIBEAM/NNBAR %experiment program is a proposed two-stage experiment at the European Spallation Source (ESS) to search for baryon number violation. The program will include high-sensitivity searches for processes that violate baryon number by one or two units: free neutron-antineutron oscillation ($n\rightarrow \bar{n}$) via mixing, neutron-antineutron oscillation via regeneration from a sterile neutron state ($n\rightarrow [n',\bar{n}'] \rightarrow \bar{n}$), and neutron disappearance ($n\rightarrow n'$); the effective $Δ\mathcal{B}=0$ process of neutron regeneration ($n\rightarrow [n',\bar{n}'] \rightarrow n$) is also possible. The program can be used to discover and characterise mixing in the neutron, antineutron, and sterile neutron sectors. The experiment addresses topical open questions such as the origins of baryogenesis, the nature of dark matter, and is sensitive to scales of new physics substantially in excess of those available at colliders. A goal of the program is to open a discovery window to neutron conversion probabilities (sensitivities) by up to three orders of magnitude compared with previous searches. The opportunity to make such a leap in sensitivity tests should not be squandered. The experiment pulls together a diverse international team of physicists from the particle (collider and low energy) and nuclear physics communities, while also including specialists in neutronics and magnetics. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Multiplicity for a Group Sequential Trial with Biomarker Subpopulations

Authors: Ting-Yu Chen, Jing Zhao, Linda Sun, Keaven Anderson

Abstract: Biomarker subpopulations have become increasingly important for drug development in targeted therapies. The use of biomarkers has the potential to facilitate more effective outcomes by guiding patient selection appropriately, thus enhancing the benefit-risk profile and improving trial power. Studying a broad population simultaneously with a more targeted one allows the trial to determine the popul… ▽ More Biomarker subpopulations have become increasingly important for drug development in targeted therapies. The use of biomarkers has the potential to facilitate more effective outcomes by guiding patient selection appropriately, thus enhancing the benefit-risk profile and improving trial power. Studying a broad population simultaneously with a more targeted one allows the trial to determine the population for which a treatment is effective and allows a goal of making approved regulatory labeling as inclusive as is appropriate. We examine new methods accounting for the complete correlation structure in group sequential designs with hypotheses in nested subgroups. The designs provide full control of family-wise Type I error rate. This extension of previous methods accounting for either group sequential design or correlation between subgroups improves efficiency (power or sample size) over a typical Bonferroni approach for testing nested populations. △ Less

Submitted 6 August, 2020; v1 submitted 17 April, 2020; originally announced April 2020.

Comments: 24 pages, 4 figures

Radial Star Formation Histories in 32 Nearby Galaxies

Authors: Daniel A. Dale, Kristin R. Anderson, Louis M. Bran, Isaiah S. Cox, Carolyn L. Drake, Nathan J. Lee, Jacob D. Pilawa, F. Alexander Slane, Susana Soto, Emily I. Jensen, Jessica S. Sutter, Jordan A. Turner, Henry A. Kobulnicky

Abstract: The spatially resolved star formation histories are studied for 32 normal star-forming galaxies drawn from the the Spitzer Extended Disk Galaxy Exploration Science survey. At surface brightness sensitivities fainter than 28 mag arcsec$^{-2}$, the new optical photometry is deep enough to complement archival ultraviolet and infrared imaging and to explore the properties of the emission well beyond t… ▽ More The spatially resolved star formation histories are studied for 32 normal star-forming galaxies drawn from the the Spitzer Extended Disk Galaxy Exploration Science survey. At surface brightness sensitivities fainter than 28 mag arcsec$^{-2}$, the new optical photometry is deep enough to complement archival ultraviolet and infrared imaging and to explore the properties of the emission well beyond the traditional optical extents of these nearby galaxies. Fits to the spectral energy distributions using a delayed star formation history model indicate a subtle but interesting average radial trend for the spiral galaxies: the inner stellar systems decrease in age with increasing radius, consistent with inside-out disk formation, but the trend reverses in the outermost regions with the stellar age nearly as old as the innermost stars. These results suggest an old stellar outer disk population formed through radial migration and/or the cumulative history of minor mergers and accretions of satellite dwarf galaxies. The subset of S0 galaxies studied here show the opposite trend compared to what is inferred for spirals: characteristic stellar ages that are increasingly older with radius for the inner portions of the galaxies, and increasingly younger stellar ages for the outer portions. This result suggests that either S0 galaxies are not well modeled by a delayed-$τ$ model, and/or that S0 galaxies have a more complicated formation history than spiral galaxies. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: Accepted for publication in the Astronomical Journal. arXiv admin note: text overlap with arXiv:1511.03285

Computational study of the dynamics of an asymmetric wedge billiard

Authors: K. D. Anderson, C. M. Villet

Abstract: We introduce the asymmetric wedge billiard as a generalization of the wedge billiard first introduced and studied by Lehtihet and Miller in 1986. We introduce the asymmetric wedge billiard as a generalization of the wedge billiard first introduced and studied by Lehtihet and Miller in 1986. △ Less

Submitted 23 February, 2020; originally announced February 2020.

Comments: 22 pages, 11 figures. International Journal of Bifurcation and Chaos (July 2020)

MSC Class: 70F99

Journal ref: International Journal of Bifurcation and Chaos vol 31, no 2 (2021)

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

Authors: Z. S. Qu, R. L. Dewar, F. Ebrahimi, J. K. Anderson, S. R. Hudson, M. J. Hole

Abstract: The Multi-region Relaxed MHD (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The Stepped Pressure Equilibrium Cod… ▽ More The Multi-region Relaxed MHD (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The Stepped Pressure Equilibrium Code (SPEC) [S.R. Hudson et al., Phys. Plasmas 19, 112502 (2012)] was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad-Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Finally, we provide an example equilibrium with a 3D helical field structure as the favoured lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated. △ Less

Submitted 20 January, 2020; originally announced January 2020.

Mining User Behaviour from Smartphone data: a literature review

Authors: Valentino Servizi, Francisco C. Pereira, Marie K. Anderson, Otto A. Nielsen

Abstract: To study users' travel behaviour and travel time between origin and destination, researchers employ travel surveys. Although there is consensus in the field about the potential, after over ten years of research and field experimentation, Smartphone-based travel surveys still did not take off to a large scale. Here, computer intelligence algorithms take the role that operators have in Traditional T… ▽ More To study users' travel behaviour and travel time between origin and destination, researchers employ travel surveys. Although there is consensus in the field about the potential, after over ten years of research and field experimentation, Smartphone-based travel surveys still did not take off to a large scale. Here, computer intelligence algorithms take the role that operators have in Traditional Travel Surveys; since we train each algorithm on data, performances rest on the data quality, thus on the ground truth. Inaccurate validations affect negatively: labels, algorithms' training, travel diaries precision, and therefore data validation, within a very critical loop. Interestingly, boundaries are proven burdensome to push even for Machine Learning methods. To support optimal investment decisions for practitioners, we expose the drivers they should consider when assessing what they need against what they get. This paper highlights and examines the critical aspects of the underlying research and provides some recommendations: (i) from the device perspective, on the main physical limitations; (ii) from the application perspective, the methodological framework deployed for the automatic generation of travel diaries; (iii)from the ground truth perspective, the relationship between user interaction, methods, and data. △ Less

Submitted 3 February, 2020; v1 submitted 24 December, 2019; originally announced December 2019.

Spontaneous electron emission from hot silver dimer anions: Breakdown of the Born-Oppenheimer approximation

Authors: E. K. Anderson, A. F. Schmidt-May, P. K. Najeeb, G. Eklund, K. C. Chartkunchand, S. Rosén, Å. Larson, K. Hansen, H. Cederquist, H. Zettergren, H. T. Schmidt

Abstract: We report the first experimental evidence of spontaneous electron emission from a homonuclear dimer anion through direct measurements of $\rm{Ag}_2^- \rightarrow \rm{Ag}_2 + \rm{e}^-$ decays on milliseconds and seconds time scales. This observation is very surprising as there is no avoided crossing between adiabatic energy curves to mediate such a process. The process is weak but yet dominates the… ▽ More We report the first experimental evidence of spontaneous electron emission from a homonuclear dimer anion through direct measurements of $\rm{Ag}_2^- \rightarrow \rm{Ag}_2 + \rm{e}^-$ decays on milliseconds and seconds time scales. This observation is very surprising as there is no avoided crossing between adiabatic energy curves to mediate such a process. The process is weak but yet dominates the decay signal after 100 ms when ensembles of internally hot Ag$_2^-$ ions are stored in the cryogenic ion-beam storage ring, DESIREE, for 10 seconds. The electron emission process is associated with an instantaneous, very large, reduction of the vibrational energy of the dimer system. This represents a dramatic deviation from a Born-Oppenheimer description of dimer dynamics. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Journal ref: Phys. Rev. Lett. 124, 173001 (2020)

PharML.Bind: Pharmacologic Machine Learning for Protein-Ligand Interactions

Authors: Aaron D. Vose, Jacob Balma, Damon Farnsworth, Kaylie Anderson, Yuri K. Peterson

Abstract: Is it feasible to create an analysis paradigm that can analyze and then accurately and quickly predict known drugs from experimental data? PharML.Bind is a machine learning toolkit which is able to accomplish this feat. Utilizing deep neural networks and big data, PharML.Bind correlates experimentally-derived drug affinities and protein-ligand X-ray structures to create novel predictions. The util… ▽ More Is it feasible to create an analysis paradigm that can analyze and then accurately and quickly predict known drugs from experimental data? PharML.Bind is a machine learning toolkit which is able to accomplish this feat. Utilizing deep neural networks and big data, PharML.Bind correlates experimentally-derived drug affinities and protein-ligand X-ray structures to create novel predictions. The utility of PharML.Bind is in its application as a rapid, accurate, and robust prediction platform for discovery and personalized medicine. This paper demonstrates that graph neural networks (GNNs) can be trained to screen hundreds of thousands of compounds against thousands of targets in minutes, a vastly shorter time than previous approaches. This manuscript presents results from training and testing using the entirety of BindingDB after cleaning; this includes a test set with 19,708 X-ray structures and 247,633 drugs, leading to 2,708,151 unique protein-ligand pairings. PharML.Bind achieves a prodigious 98.3% accuracy on this test set in under 25 minutes. PharML.Bind is premised on the following key principles: 1) speed and a high enrichment factor per unit compute time, provided by high-quality training data combined with a novel GNN architecture and use of high-performance computing resources, 2) the ability to generalize to proteins and drugs outside of the training set, including those with unknown active sites, through the use of an active-site-agnostic GNN mapping, and 3) the ability to be easily integrated as a component of increasingly-complex prediction and analysis pipelines. PharML.Bind represents a timely and practical approach to leverage the power of machine learning to efficiently analyze and predict drug action on any practical scale and will provide utility in a variety of discovery and medical applications. △ Less

Submitted 23 October, 2019; originally announced November 2019.

Microgrid Resilience: A holistic approach for assessing threats, identifying vulnerabilities, and designing corresponding mitigation strategies

Authors: Sakshi Mishra, Kate Anderson, Brian Miller, Kyle Boyer, Adam Warren

Abstract: Microgrids are being increasing deployed to improve the operational flexibility, resilience, coordinated-energy management capabilities, self-adequacy, and increased reliability of power systems. This strong market growth is also driven by advances in power electronics, improved control systems, and the rapidly falling price and increased adoption of distributed energy generation technologies, lik… ▽ More Microgrids are being increasing deployed to improve the operational flexibility, resilience, coordinated-energy management capabilities, self-adequacy, and increased reliability of power systems. This strong market growth is also driven by advances in power electronics, improved control systems, and the rapidly falling price and increased adoption of distributed energy generation technologies, like solar photovoltaics and storage. In the event of grid outages, microgrids can provide a backup source of power; providing resilience to the critical loads; however, this requires that the microgrid itself is resilient to physical and cyber threats. Building highly resilient microgrids requires a methodological assessment of potential threats, identification of vulnerabilities, and design of mitigation strategies. This paper provides a comprehensive review of threats, vulnerabilities, and mitigation strategies and develops a definition for microgrid resilience. The paper also develops a methodology for designing resilient microgrids by considering how microgrid designers and site owners evaluate threats, vulnerabilities, and consequences and choose the microgrid features required to address these threats under different situations. △ Less

Submitted 1 January, 2020; v1 submitted 2 October, 2019; originally announced October 2019.

Comments: 24 pages, 9 figures

Journal ref: https://doi.org/10.1016/j.apenergy.2020.114726

Alternative Analysis Methods for Time to Event Endpoints under Non-proportional Hazards: A Comparative Analysis

Authors: Ray S. Lin, Ji Lin, Satrajit Roychoudhury, Keaven M. Anderson, Tianle Hu, Bo Huang, Larry F Leon, Jason JZ Liao, Rong Liu, Xiaodong Luo, Pralay Mukhopadhyay, Rui Qin, Kay Tatsuoka, Xuejing Wang, Yang Wang, Jian Zhu, Tai-Tsang Chen, Renee Iacona, Cross-Pharma Non-proportional Hazards Working Group

Abstract: The log-rank test is most powerful under proportional hazards (PH). In practice, non-PH patterns are often observed in clinical trials, such as in immuno-oncology; therefore, alternative methods are needed to restore the efficiency of statistical testing. Three categories of testing methods were evaluated, including weighted log-rank tests, Kaplan-Meier curve-based tests (including weighted Kaplan… ▽ More The log-rank test is most powerful under proportional hazards (PH). In practice, non-PH patterns are often observed in clinical trials, such as in immuno-oncology; therefore, alternative methods are needed to restore the efficiency of statistical testing. Three categories of testing methods were evaluated, including weighted log-rank tests, Kaplan-Meier curve-based tests (including weighted Kaplan-Meier and Restricted Mean Survival Time, RMST), and combination tests (including Breslow test, Lee's combo test, and MaxCombo test). Nine scenarios representing the PH and various non-PH patterns were simulated. The power, type I error, and effect estimates of each method were compared. In general, all tests control type I error well. There is not a single most powerful test across all scenarios. In the absence of prior knowledge regarding the PH or non-PH patterns, the MaxCombo test is relatively robust across patterns. Since the treatment effect changes overtime under non-PH, the overall profile of the treatment effect may not be represented comprehensively based on a single measure. Thus, multiple measures of the treatment effect should be pre-specified as sensitivity analyses to evaluate the totality of the data. △ Less

Submitted 20 September, 2019; originally announced September 2019.

Report number: NPH12

Journal ref: Statistics in Biopharmaceutical Statistics 2020

Robust Design and Analysis of Clinical Trials With Non-proportional Hazards: A Straw Man Guidance from a Cross-pharma Working Group

Authors: Satrajit Roychoudhury, Keaven M Anderson, Jiabu Ye, Pralay Mukhopadhyay

Abstract: Loss of power and clear description of treatment differences are key issues in designing and analyzing a clinical trial where non-proportional hazard is a possibility. A log-rank test may be very inefficient and interpretation of the hazard ratio estimated using Cox regression is potentially problematic. In this case, the current ICH E9 (R1) addendum would suggest designing a trial with a clinical… ▽ More Loss of power and clear description of treatment differences are key issues in designing and analyzing a clinical trial where non-proportional hazard is a possibility. A log-rank test may be very inefficient and interpretation of the hazard ratio estimated using Cox regression is potentially problematic. In this case, the current ICH E9 (R1) addendum would suggest designing a trial with a clinically relevant estimand, e.g., expected life gain. This approach considers appropriate analysis methods for supporting the chosen estimand. However, such an approach is case specific and may suffer lack of power for important choices of the underlying alternate hypothesis distribution. On the other hand, there may be a desire to have robust power under different deviations from proportional hazards. Also, we would contend that no single number adequately describes treatment effect under non-proportional hazards scenarios. The cross-pharma working group has proposed a combination test to provide robust power under a variety of alternative hypotheses. These can be specified for primary analysis at the design stage and methods appropriately accounting for combination test correlations are efficient for a variety of scenarios. We have provided design and analysis considerations based on a combination test under different non-proportional hazard types and present a straw man proposal for practitioners. The proposals are illustrated with real life example and simulation. △ Less

Submitted 12 January, 2021; v1 submitted 19 August, 2019; originally announced August 2019.

Journal ref: Statistics in Biopharmaceutical Research 2021

In-Situ Scattering of Warm Jupiters and Implications for Dynamical Histories

Authors: Kassandra R. Anderson, Dong Lai, Bonan Pu

Abstract: Many warm Jupiters (WJs) have substantial eccentricities, which are linked to their formation and migration histories. This paper explores eccentricity excitation of WJs due to planet-planet scattering, beginning with 3-4 planets in unstable orbits, with the innermost planet placed in the range (0.1 - 1)AU. Such a setup is consistent with either in-situ formation or arrival at sub-AU orbits due to… ▽ More Many warm Jupiters (WJs) have substantial eccentricities, which are linked to their formation and migration histories. This paper explores eccentricity excitation of WJs due to planet-planet scattering, beginning with 3-4 planets in unstable orbits, with the innermost planet placed in the range (0.1 - 1)AU. Such a setup is consistent with either in-situ formation or arrival at sub-AU orbits due to disk migration. Most previous N-body experiments have focused on "cold" Jupiters at several AU, where scattering results in planet ejections, efficiently exciting the eccentricities of surviving planets. In contrast, scattering at sub-AU distances results in a mixture of collisions and ejections, and the final eccentricities of surviving planets are unclear. We conduct scattering experiments for a range of planet masses and initial spacings, including the effect of general relativistic apsidal precession, and systematically catalogue the scattering outcomes and properties of surviving planets. A comparable number of one-planet and two-planet systems are produced. Two-planet systems arise exclusively through planet-planet collisions, and tend to have low eccentricities/mutual inclinations and compact configurations. One-planet systems arise through a combination of ejections and collisions, resulting in higher eccentricities. The observed eccentricity distribution of solitary WJs (lacking detection of a giant planet companion) is consistent with roughly 60% of the systems having undergone in-situ scattering, and the remaining experiencing a quiescent history. △ Less

Submitted 5 November, 2019; v1 submitted 12 August, 2019; originally announced August 2019.

Comments: 16 pages, 12 figures, accepted for publication in MNRAS