-
Active Healing of Microtubule-Motor Networks
Authors:
Fan Yang,
Shichen Liu,
Heun Jin Lee,
Rob Phillips,
Matt Thomson
Abstract:
Cytoskeletal networks have a self-healing property where networks can repair defects to maintain structural integrity. However, both the mechanisms and dynamics of healing remain largely unknown. Here we report an unexplored healing mechanism in microtubule-motor networks by active crosslinking. We directly generate network cracks using a light-controlled microtubule-motor system, and observe that…
▽ More
Cytoskeletal networks have a self-healing property where networks can repair defects to maintain structural integrity. However, both the mechanisms and dynamics of healing remain largely unknown. Here we report an unexplored healing mechanism in microtubule-motor networks by active crosslinking. We directly generate network cracks using a light-controlled microtubule-motor system, and observe that the cracks can self-heal. Combining theory and experiment, we find that the networks must overcome internal elastic resistance in order to heal cracks, giving rise to a bifurcation of dynamics dependent on the initial opening angle of the crack: the crack heals below a critical angle and opens up at larger angles. Simulation of a continuum model reproduces the bifurcation dynamics, revealing the importance of a boundary layer where free motors and microtubules can actively crosslink and thereby heal the crack. We also formulate a simple elastic-rod model that can qualitatively predict the critical angle, which is found to be tunable by two dimensionless geometric parameters, the ratio of the boundary layer and network width, and the aspect ratio of the network. Our results provide a new framework for understanding healing in cytoskeletal networks and designing self-healable biomaterials.
△ Less
Submitted 30 June, 2024;
originally announced July 2024.
-
A minimal scenario for the origin of non-equilibrium order
Authors:
Riccardo Ravasio,
Kabir Husain,
Constantine G. Evans,
Rob Phillips,
Marco Ribezzi,
Jack W. Szostak,
Arvind Murugan
Abstract:
Extant life contains numerous non-equilibrium mechanisms to create order not achievable at equilibrium; it is generally assumed that these mechanisms evolved because the resulting order was sufficiently beneficial to overcome associated costs of time and energy. Here, we identify a broad range of conditions under which non-equilibrium order-creating mechanisms will evolve as an inevitable conseque…
▽ More
Extant life contains numerous non-equilibrium mechanisms to create order not achievable at equilibrium; it is generally assumed that these mechanisms evolved because the resulting order was sufficiently beneficial to overcome associated costs of time and energy. Here, we identify a broad range of conditions under which non-equilibrium order-creating mechanisms will evolve as an inevitable consequence of self-replication, even if the order is not directly functional. We show that models of polymerases, when expanded to include known stalling effects, can evolve kinetic proofreading through selection for fast replication alone, consistent with data from recent mutational screens. Similarly, replication contingent on fast self-assembly can select for non-equilibrium instabilities and result in more ordered structures without any direct selection for order. We abstract these results into a framework that predicts that self-replication intrinsically amplifies dissipative order-enhancing mechanisms if the distribution of replication times is wide enough. Our work suggests the intriguing possibility that non-equilibrium order can arise more easily than assumed, even before that order is directly functional, with consequences impacting mutation rate evolution and kinetic traps in self-assembly to the origin of life.
△ Less
Submitted 17 May, 2024;
originally announced May 2024.
-
The Shape of Money Laundering: Subgraph Representation Learning on the Blockchain with the Elliptic2 Dataset
Authors:
Claudio Bellei,
Muhua Xu,
Ross Phillips,
Tom Robinson,
Mark Weber,
Tim Kaler,
Charles E. Leiserson,
Arvind,
Jie Chen
Abstract:
Subgraph representation learning is a technique for analyzing local structures (or shapes) within complex networks. Enabled by recent developments in scalable Graph Neural Networks (GNNs), this approach encodes relational information at a subgroup level (multiple connected nodes) rather than at a node level of abstraction. We posit that certain domain applications, such as anti-money laundering (A…
▽ More
Subgraph representation learning is a technique for analyzing local structures (or shapes) within complex networks. Enabled by recent developments in scalable Graph Neural Networks (GNNs), this approach encodes relational information at a subgroup level (multiple connected nodes) rather than at a node level of abstraction. We posit that certain domain applications, such as anti-money laundering (AML), are inherently subgraph problems and mainstream graph techniques have been operating at a suboptimal level of abstraction. This is due in part to the scarcity of annotated datasets of real-world size and complexity, as well as the lack of software tools for managing subgraph GNN workflows at scale. To enable work in fundamental algorithms as well as domain applications in AML and beyond, we introduce Elliptic2, a large graph dataset containing 122K labeled subgraphs of Bitcoin clusters within a background graph consisting of 49M node clusters and 196M edge transactions. The dataset provides subgraphs known to be linked to illicit activity for learning the set of "shapes" that money laundering exhibits in cryptocurrency and accurately classifying new criminal activity. Along with the dataset we share our graph techniques, software tooling, promising early experimental results, and new domain insights already gleaned from this approach. Taken together, we find immediate practical value in this approach and the potential for a new standard in anti-money laundering and forensic analytics in cryptocurrencies and other financial networks.
△ Less
Submitted 1 May, 2024; v1 submitted 29 April, 2024;
originally announced April 2024.
-
From chiral EFT to perturbative QCD: a Bayesian model mixing approach to symmetric nuclear matter
Authors:
A. C. Semposki,
C. Drischler,
R. J. Furnstahl,
J. A. Melendez,
D. R. Phillips
Abstract:
Constraining the equation of state (EOS) of strongly interacting, dense matter is the focus of intense experimental, observational, and theoretical effort. Chiral effective field theory ($χ$EFT) can describe the EOS between the typical densities of nuclei and those in the outer cores of neutron stars while perturbative QCD (pQCD) can be applied to properties of deconfined quark matter, both with q…
▽ More
Constraining the equation of state (EOS) of strongly interacting, dense matter is the focus of intense experimental, observational, and theoretical effort. Chiral effective field theory ($χ$EFT) can describe the EOS between the typical densities of nuclei and those in the outer cores of neutron stars while perturbative QCD (pQCD) can be applied to properties of deconfined quark matter, both with quantified theoretical uncertainties. However, describing the complete range of densities between nuclear saturation and an almost-free quark gas with a single EOS that has well-quantified uncertainties is a challenging problem. In this work, we argue that Bayesian multi-model inference from $χ$EFT and pQCD can help bridge the gap between the two theories: we combine the Gaussian random variables that constitute the theories' predictions for the pressure as a function of the density in symmetric nuclear matter. We do this using two Bayesian model mixing procedures: a pointwise approach, and a correlated approach implemented via a Gaussian process (GP), and present results for the pressure and speed of sound in each. The second method produces a smooth $χ$EFT-to-pQCD EOS. Without input data in the intermediate region, the choice of prior on the EOS, encoded through the GP kernel, as the prior on the EOS function space significantly affects the result in that region. We also discuss future extensions and applications to neutron star matter guided by recent EOS constraints from nuclear theory, nuclear experiment, and multi-messenger astronomy.
△ Less
Submitted 9 April, 2024;
originally announced April 2024.
-
Assessing Correlated Truncation Errors in Modern Nucleon-Nucleon Potentials
Authors:
P. J. Millican,
R. J. Furnstahl,
J. A. Melendez,
D. R. Phillips,
M. T. Pratola
Abstract:
We test the BUQEYE model of correlated effective field theory (EFT) truncation errors on Reinert, Krebs, and Epelbaum's semi-local momentum-space implementation of the chiral EFT ($χ$EFT) expansion of the nucleon-nucleon (NN) potential. This Bayesian model hypothesizes that dimensionless coefficient functions extracted from the order-by-order corrections to NN observables can be treated as draws f…
▽ More
We test the BUQEYE model of correlated effective field theory (EFT) truncation errors on Reinert, Krebs, and Epelbaum's semi-local momentum-space implementation of the chiral EFT ($χ$EFT) expansion of the nucleon-nucleon (NN) potential. This Bayesian model hypothesizes that dimensionless coefficient functions extracted from the order-by-order corrections to NN observables can be treated as draws from a Gaussian process (GP). We combine a variety of graphical and statistical diagnostics to assess when predicted observables have a $χ$EFT convergence pattern consistent with the hypothesized GP statistical model. Our conclusions are: First, the BUQEYE model is generally applicable to the potential investigated here, which enables statistically principled estimates of the impact of higher EFT orders on observables. Second, parameters defining the extracted coefficients such as the expansion parameter $Q$ must be well chosen for the coefficients to exhibit a regular convergence pattern -- a property we exploit to obtain posterior distributions for such quantities. Third, the assumption of GP stationarity across lab energy and scattering angle is not generally met; this necessitates adjustments in future work. We provide a workflow and interpretive guide for our analysis framework, and show what can be inferred about probability distributions for $Q$, the EFT breakdown scale $Λ_b$, the scale associated with soft physics in the $χ$EFT potential $m_{\rm eff}$, and the GP hyperparameters. All our results can be reproduced using a publicly available Jupyter notebook, which can be straightforwardly modified to analyze other $χ$EFT NN potentials.
△ Less
Submitted 20 February, 2024;
originally announced February 2024.
-
Dueling Over Dessert, Mastering the Art of Repeated Cake Cutting
Authors:
Simina Brânzei,
MohammadTaghi Hajiaghayi,
Reed Phillips,
Suho Shin,
Kun Wang
Abstract:
We consider the setting of repeated fair division between two players, denoted Alice and Bob, with private valuations over a cake. In each round, a new cake arrives, which is identical to the ones in previous rounds. Alice cuts the cake at a point of her choice, while Bob chooses the left piece or the right piece, leaving the remainder for Alice. We consider two versions: sequential, where Bob obs…
▽ More
We consider the setting of repeated fair division between two players, denoted Alice and Bob, with private valuations over a cake. In each round, a new cake arrives, which is identical to the ones in previous rounds. Alice cuts the cake at a point of her choice, while Bob chooses the left piece or the right piece, leaving the remainder for Alice. We consider two versions: sequential, where Bob observes Alice's cut point before choosing left/right, and simultaneous, where he only observes her cut point after making his choice. The simultaneous version was first considered by Aumann and Maschler (1995).
We observe that if Bob is almost myopic and chooses his favorite piece too often, then he can be systematically exploited by Alice through a strategy akin to a binary search. This strategy allows Alice to approximate Bob's preferences with increasing precision, thereby securing a disproportionate share of the resource over time.
We analyze the limits of how much a player can exploit the other one and show that fair utility profiles are in fact achievable. Specifically, the players can enforce the equitable utility profile of $(1/2, 1/2)$ in the limit on every trajectory of play, by keeping the other player's utility to approximately $1/2$ on average while guaranteeing they themselves get at least approximately $1/2$ on average. We show this theorem using a connection with Blackwell approachability.
Finally, we analyze a natural dynamic known as fictitious play, where players best respond to the empirical distribution of the other player. We show that fictitious play converges to the equitable utility profile of $(1/2, 1/2)$ at a rate of $O(1/\sqrt{T})$.
△ Less
Submitted 18 February, 2024; v1 submitted 13 February, 2024;
originally announced February 2024.
-
Compton Scattering on 4He with Nuclear One- and Two-Body Densities
Authors:
Harald W. Griesshammer,
Junjie Liao,
Judith A. McGovern,
Andreas Nogga,
Daniel R. Phillips
Abstract:
We present the first \emph{ab initio} calculation of elastic Compton scattering from 4He. It is carried out to $\mathcal{O}(e^2 δ^3)$ [N3LO] in the $δ$ expansion of $χ$EFT. At this order and for this target, the only free parameters are the scalar-isoscalar electric and magnetic dipole polarisabilities of the nucleon. Adopting current values for these yields a parameter-free prediction. This compa…
▽ More
We present the first \emph{ab initio} calculation of elastic Compton scattering from 4He. It is carried out to $\mathcal{O}(e^2 δ^3)$ [N3LO] in the $δ$ expansion of $χ$EFT. At this order and for this target, the only free parameters are the scalar-isoscalar electric and magnetic dipole polarisabilities of the nucleon. Adopting current values for these yields a parameter-free prediction. This compares favourably with the world data from HI$γ$S, Illinois and Lund for photon energies $50\;\mathrm{MeV}\lesssimω\lesssim120\;\mathrm{MeV}$ within our theoretical uncertainties of $\pm10\%$. We predict a cross section up to 7 times that for deuterium. As in 3He, this emphasises and tests the key role of meson-exchange currents between np pairs in Compton scattering on light nuclei. We assess the sensitivity of the cross section and beam asymmetry to the nucleon polarisabilities, providing clear guidance to future experiments seeking to further constrain them. The calculation becomes tractable by use of the Transition Density Method. The one- and two-body densities generated from 5 chiral potentials and the AV18$+$UIX potential are available using the python package provided at \url{https://pypi.org/project/nucdens/}.
△ Less
Submitted 20 June, 2024; v1 submitted 30 January, 2024;
originally announced January 2024.
-
Deciphering regulatory architectures from synthetic single-cell expression patterns
Authors:
Rosalind Wenshan Pan,
Tom Roeschinger,
Kian Faizi,
Hernan Garcia,
Rob Phillips
Abstract:
For the vast majority of genes in sequenced genomes, there is limited understanding of how they are regulated. Without such knowledge, it is not possible to perform a quantitative theory-experiment dialogue on how such genes give rise to physiological and evolutionary adaptation. One category of high-throughput experiments used to understand the sequence-phenotype relationship of the transcriptome…
▽ More
For the vast majority of genes in sequenced genomes, there is limited understanding of how they are regulated. Without such knowledge, it is not possible to perform a quantitative theory-experiment dialogue on how such genes give rise to physiological and evolutionary adaptation. One category of high-throughput experiments used to understand the sequence-phenotype relationship of the transcriptome is massively parallel reporter assays (MPRAs). However, to improve the versatility and scalability of MPRA pipelines, we need a "theory of the experiment" to help us better understand the impact of various biological and experimental parameters on the interpretation of experimental data. To that end, in this paper we create tens of thousands of synthetic single-cell gene expression outputs using both equilibrium and out-of-equilibrium models. These models make it possible to imitate the summary statistics (information footprints and expression shift matrices) used to characterize the output of MPRAs and from this summary statistic to infer the underlying regulatory architecture. Specifically, we use a more refined implementation of the so-called thermodynamic models in which the binding energies of each sequence variant are derived from energy matrices. Our simulations reveal important effects of the parameters on MPRA data and we demonstrate our ability to optimize MPRA experimental designs with the goal of generating thermodynamic models of the transcriptome with base-pair specificity. Further, this approach makes it possible to carefully examine the mapping between mutations in binding sites and their corresponding expression profiles, a tool useful not only for better designing MPRAs, but also for exploring regulatory evolution.
△ Less
Submitted 5 June, 2024; v1 submitted 28 January, 2024;
originally announced January 2024.
-
Using the Transition-Density Formalism in the First Computation of 4He Compton Scattering
Authors:
Harald W. Griesshammer,
Junjie Liao,
Judith A. McGovern,
Andreas Nogga,
Daniel R. Phillips
Abstract:
The method and results of the first theory description of 4He Compton scattering at nuclear energies is presented, with a focus on figures. An upcoming publication [1] contains details and a comprehensive list of references.
The method and results of the first theory description of 4He Compton scattering at nuclear energies is presented, with a focus on figures. An upcoming publication [1] contains details and a comprehensive list of references.
△ Less
Submitted 28 January, 2024;
originally announced January 2024.
-
Force Propagation in Active Cytoskeletal Networks
Authors:
Shichen Liu,
Rosalind Wenshan Pan,
Heun Jin Lee,
Shahriar Shadkhoo,
Fan Yang,
Chunhe Li,
Zijie Qu,
Rob Phillips,
Matt Thomson
Abstract:
In biological systems, molecular-scale forces and motions are pivotal for enabling processes like motility, shape change, and replication. These forces and motions are organized, amplified, and transmitted across macroscopic scales by active materials such as the cytoskeleton, which drives micron-scale cellular movement and re-organization. Despite the integral role of active materials, understand…
▽ More
In biological systems, molecular-scale forces and motions are pivotal for enabling processes like motility, shape change, and replication. These forces and motions are organized, amplified, and transmitted across macroscopic scales by active materials such as the cytoskeleton, which drives micron-scale cellular movement and re-organization. Despite the integral role of active materials, understanding how molecular-scale interactions alter macroscopic structure and force propagation remains elusive. This knowledge gap presents challenges to the harnessing and regulation of such dynamics across diverse length scales. Here, we demonstrate how mediating the bundling of microtubules can shift active matter between a global force-transmitting phase and a local force-dissipating phase. A fivefold increase in microtubule effective length results in the transition from local to global phase with a hundredfold increase in velocity autocorrelation. Through theory and simulation, we identify signatures of a percolation-driven transition between the two phases. This provides evidence for how force propagation can be generated when local molecular interactions reach a sufficient length scale. We show that force propagation in the active matter system enables material transport. Consequently, we demonstrate that the global phase is capable of facilitating millimeter-scale human cell transport and manipulation, as well as powering the movement of aqueous droplets. These findings underscore the potential for designing active materials capable of force organization and transmission. Our results lay the foundation for further exploration into the organization and propagation of forces/stresses in biological systems, thereby paving the way for the engineering of active materials in synthetic biology and soft robotics.
△ Less
Submitted 12 April, 2024; v1 submitted 8 January, 2024;
originally announced January 2024.
-
Universality of $nn$ distributions of $s$-wave $2n$ halos and the unitary limit
Authors:
Matthias Göbel,
Hans-Werner Hammer,
Daniel R. Phillips
Abstract:
We calculate neutron-neutron relative-energy distributions of $s$-wave two-neutron ($2n$) halo nuclei using Halo Effective Field Theory (Halo EFT) at leading order. At this order these systems are described by the $2n$ separation energy, the neutron-core ($nc$) virtual-state energy and the neutron-neutron ($nn$) scattering length. We focus on knockout reactions where the removal of the core is sud…
▽ More
We calculate neutron-neutron relative-energy distributions of $s$-wave two-neutron ($2n$) halo nuclei using Halo Effective Field Theory (Halo EFT) at leading order. At this order these systems are described by the $2n$ separation energy, the neutron-core ($nc$) virtual-state energy and the neutron-neutron ($nn$) scattering length. We focus on knockout reactions where the removal of the core is sudden, such that the final-state interactions are dominated by the $nn$ interaction. We consider the neutron relative-energy distribution for the nuclei $^{11}$Li, $^{14}$Be, $^{17}$B, $^{19}$B, and $^{22}$C. We show that the ground-state neutron momentum distributions of all these nuclei stem from a single curve, which can be obtained by taking both the neutron-core and neutron-neutron interaction to the unitary limit. This universal description can be extended to the final distribution measured in experiment by including $nn$ final-state interactions via the approximate technique of enhancement factors. For all the nuclei considered we find good agreement between the full leading-order Halo EFT calculation and the universal prediction obtained in this way. The universality of the ground-state momentum distribution in two-neutron Borromean halos can thus be tested by dividing the experimental results from sudden core knockout by the enhancement factor and comparing to the unitary-limit prediction.
△ Less
Submitted 7 September, 2023;
originally announced September 2023.
-
Effective Field Theory for the Bound States and Scattering of a Heavy Charged Particle and a Neutral Atom
Authors:
Daniel Odell,
Daniel R. Phillips,
Ubirajara van Kolck
Abstract:
We show the system of a heavy charged particle and a neutral atom can be described by a low-energy effective field theory where the attractive $1/r^4$ induced dipole potential determines the long-distance/low-energy wave functions. The $1/r^4$ interaction is renormalized by a contact interaction at leading order. Derivative corrections to that contact interaction give rise to higher-order terms. W…
▽ More
We show the system of a heavy charged particle and a neutral atom can be described by a low-energy effective field theory where the attractive $1/r^4$ induced dipole potential determines the long-distance/low-energy wave functions. The $1/r^4$ interaction is renormalized by a contact interaction at leading order. Derivative corrections to that contact interaction give rise to higher-order terms. We show that this ``Induced-dipole EFT'' (ID-EFT) reproduces the $π^+$-hydrogen phase shifts of a more microscopic potential, the Temkin-Lamkin potential, over a wide range of energies. Already at leading order it also describes the highest-lying excited bound states of the pionic-hydrogen ion. Lower-lying bound states receive substantial corrections at next-to-leading order, with the size of the correction proportional to their distance from the scattering threshold. Our next-to-leading order calculation shows that the three highest-lying bound states of the Temkin-Lamkin potential are well-described in ID-EFT.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
High-sensitivity dual-comb and cross-comb spectroscopy across the infrared using a widely-tunable and free-running optical parametric oscillator
Authors:
Carolin P. Bauer,
Zofia A. Bejm,
Michelle K. Bollier,
Justinas Pupeikis,
Benjamin Willenberg,
Ursula Keller,
Christopher R. Phillips
Abstract:
Coherent dual-comb spectroscopy (DCS) enables high-resolution measurements at high speeds without the trade-off between resolution and update rate inherent to mechanical delay scanning approaches. However, high system complexity and limited measurement sensitivity remain major challenges for DCS. Here, we address these challenges via a wavelength-tunable dual-comb optical parametric oscillator (OP…
▽ More
Coherent dual-comb spectroscopy (DCS) enables high-resolution measurements at high speeds without the trade-off between resolution and update rate inherent to mechanical delay scanning approaches. However, high system complexity and limited measurement sensitivity remain major challenges for DCS. Here, we address these challenges via a wavelength-tunable dual-comb optical parametric oscillator (OPO) combined with an up-conversion detection method. The OPO is tunable in the short-wave infrared (1300-1670 nm range) and mid-infrared (2700- 5000 nm range) where many molecules have strong absorption bands. Both OPO pump beams are generated in a single spatially-multiplexed laser cavity, while both signal and idler beams are generated in a single spatially-multiplexed OPO cavity. The near-common path of the combs in this new configuration enables comb-line-resolved and aliasing-free measurements in free-running operation. By limiting the instantaneous idler bandwidth to below 1 THz, we obtain a high power per comb line in the mid-infrared of up to 160 $μ$W. With a novel intra-cavity nonlinear up-conversion scheme based on cross-comb spectroscopy, we leverage these power levels while overcoming the sensitivity limitations of direct mid-infrared detection, leading to a high signal-to-noise ratio (50.2 dB Hz$^{1/2}$) and record-level dual-comb figure of merit (3.5\times 10^8 Hz$^{1/2}$). As a proof of concept, we demonstrate the detection of methane with 2-ppm concentration over 3-m path length. Our results demonstrate a new paradigm for DCS compatible with high-sensitivity and high-resolution measurements over a wide spectral range.
△ Less
Submitted 13 March, 2024; v1 submitted 4 May, 2023;
originally announced May 2023.
-
Application of targeted maximum likelihood estimation in public health and epidemiological studies: a systematic review
Authors:
Matthew J. Smith,
Rachael V. Phillips,
Miguel Angel Luque-Fernandez,
Camille Maringe
Abstract:
The Targeted Maximum Likelihood Estimation (TMLE) statistical data analysis framework integrates machine learning, statistical theory, and statistical inference to provide a least biased, efficient and robust strategy for estimation and inference of a variety of statistical and causal parameters. We describe and evaluate the epidemiological applications that have benefited from recent methodologic…
▽ More
The Targeted Maximum Likelihood Estimation (TMLE) statistical data analysis framework integrates machine learning, statistical theory, and statistical inference to provide a least biased, efficient and robust strategy for estimation and inference of a variety of statistical and causal parameters. We describe and evaluate the epidemiological applications that have benefited from recent methodological developments. We conducted a systematic literature review in PubMed for articles that applied any form of TMLE in observational studies. We summarised the epidemiological discipline, geographical location, expertise of the authors, and TMLE methods over time. We used the Roadmap of Targeted Learning and Causal Inference to extract key methodological aspects of the publications. We showcase the contributions to the literature of these TMLE results. Of the 81 publications included, 25% originated from the University of California at Berkeley, where the framework was first developed by Professor Mark van der Laan. By the first half of 2022, 70% of the publications originated from outside the United States and explored up to 7 different epidemiological disciplines in 2021-22. Double-robustness, bias reduction and model misspecification were the main motivations that drew researchers towards the TMLE framework. Through time, a wide variety of methodological, tutorial and software-specific articles were cited, owing to the constant growth of methodological developments around TMLE. There is a clear dissemination trend of the TMLE framework to various epidemiological disciplines and to increasing numbers of geographical areas. The availability of R packages, publication of tutorial papers, and involvement of methodological experts in applied publications have contributed to an exponential increase in the number of studies that understood the benefits, and adoption, of TMLE.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
Gigahertz Single-cavity Dual-comb Laser for Rapid Time-domain Spectroscopy: from Few Terahertz to Optical Frequencies
Authors:
Benjamin Willenberg,
Christopher R. Phillips,
Justinas Pupeikis,
Sandro L. Camenzind,
Lars Liebermeister,
Robert B. Kohlhass,
Björn Globisch,
Ursula Keller
Abstract:
Dual-comb generation from a single laser cavity provides a simple and high-performance solution to time sampling applications. We demonstrate a compact single-cavity dual-comb laser operating at gigahertz repetition rates and high repetition rate differences up to more than 100 kHz with sub-100 fs pulse duration. The single cavity approach leads to passive common noise suppression resulting in ult…
▽ More
Dual-comb generation from a single laser cavity provides a simple and high-performance solution to time sampling applications. We demonstrate a compact single-cavity dual-comb laser operating at gigahertz repetition rates and high repetition rate differences up to more than 100 kHz with sub-100 fs pulse duration. The single cavity approach leads to passive common noise suppression resulting in ultra-low relative timing jitter and fully resolvable comb lines in free-running operation. We showcase the laser performance with two application demonstrations: (a) time-domain spectroscopy of acetylene in the near-infrared via computational comb line tracking and (b) free-space THz time domain spectroscopy and thickness-measurements via adaptive sampling. For (b) we use efficient state-of the art iron-doped InGaAs photoconductive antennas to generate and detect the THz light. Here we operate these devices with an efficient Yb-based gigahertz repetition rate laser for the first time. One optical comb generates the THz light, while the other probes it via equivalent time sampling. We obtain signal strengths comparable to reference measurements with MHz repetition rate Er-based laser systems while achieving close to 1 GHz spectral resolution (defined by the comb line spacing) and generating THz frequencies up to 3 THz. By carrying out a careful investigation of the noise properties of the laser we confirm that the free-running gigahertz dual-comb oscillator provides a rapid yet highly precise optical delay sweep from a simple setup. Therefore, our approach will be beneficial for high-update rate time sampling and time-domain spectroscopy applications.
△ Less
Submitted 21 February, 2023;
originally announced February 2023.
-
Multi-task Highly Adaptive Lasso
Authors:
Ivana Malenica,
Rachael V. Phillips,
Daniel Lazzareschi,
Jeremy R. Coyle,
Romain Pirracchio,
Mark J. van der Laan
Abstract:
We propose a novel, fully nonparametric approach for the multi-task learning, the Multi-task Highly Adaptive Lasso (MT-HAL). MT-HAL simultaneously learns features, samples and task associations important for the common model, while imposing a shared sparse structure among similar tasks. Given multiple tasks, our approach automatically finds a sparse sharing structure. The proposed MTL algorithm at…
▽ More
We propose a novel, fully nonparametric approach for the multi-task learning, the Multi-task Highly Adaptive Lasso (MT-HAL). MT-HAL simultaneously learns features, samples and task associations important for the common model, while imposing a shared sparse structure among similar tasks. Given multiple tasks, our approach automatically finds a sparse sharing structure. The proposed MTL algorithm attains a powerful dimension-free convergence rate of $o_p(n^{-1/4})$ or better. We show that MT-HAL outperforms sparsity-based MTL competitors across a wide range of simulation studies, including settings with nonlinear and linear relationships, varying levels of sparsity and task correlations, and different numbers of covariates and sample size.
△ Less
Submitted 27 January, 2023;
originally announced January 2023.
-
Solving the Discretised Boltzmann Transport Equations using Neural Networks: Applications in Neutron Transport
Authors:
T. R. F. Phillips,
C. E. Heaney,
C. Boyang,
A. G. Buchan,
C. C. Pain
Abstract:
In this paper we solve the Boltzmann transport equation using AI libraries. The reason why this is attractive is because it enables one to use the highly optimised software within AI libraries, enabling one to run on different computer architectures and enables one to tap into the vast quantity of community based software that has been developed for AI and ML applications e.g. mixed arithmetic pre…
▽ More
In this paper we solve the Boltzmann transport equation using AI libraries. The reason why this is attractive is because it enables one to use the highly optimised software within AI libraries, enabling one to run on different computer architectures and enables one to tap into the vast quantity of community based software that has been developed for AI and ML applications e.g. mixed arithmetic precision or model parallelism. Here we take the first steps towards developing this approach for the Boltzmann transport equation and develop the necessary methods in order to do that effectively. This includes: 1) A space-angle multigrid solution method that can extract the level of parallelism necessary to run efficiently on GPUs or new AI computers. 2) A new Convolutional Finite Element Method (ConvFEM) that greatly simplifies the implementation of high order finite elements (quadratic to quintic, say). 3) A new non-linear Petrov-Galerkin method that introduces dissipation anisotropically.
△ Less
Submitted 25 January, 2023; v1 submitted 24 January, 2023;
originally announced January 2023.
-
Solving the Discretised Neutron Diffusion Equations using Neural Networks
Authors:
T. R. F. Phillips,
C. E. Heaney,
C. Boyang,
A. G. Buchan,
C. C. Pain
Abstract:
This paper presents a new approach which uses the tools within Artificial Intelligence (AI) software libraries as an alternative way of solving partial differential equations (PDEs) that have been discretised using standard numerical methods. In particular, we describe how to represent numerical discretisations arising from the finite volume and finite element methods by pre-determining the weight…
▽ More
This paper presents a new approach which uses the tools within Artificial Intelligence (AI) software libraries as an alternative way of solving partial differential equations (PDEs) that have been discretised using standard numerical methods. In particular, we describe how to represent numerical discretisations arising from the finite volume and finite element methods by pre-determining the weights of convolutional layers within a neural network. As the weights are defined by the discretisation scheme, no training of the network is required and the solutions obtained are identical (accounting for solver tolerances) to those obtained with standard codes often written in Fortran or C++. We also explain how to implement the Jacobi method and a multigrid solver using the functions available in AI libraries. For the latter, we use a U-Net architecture which is able to represent a sawtooth multigrid method. A benefit of using AI libraries in this way is that one can exploit their power and their built-in technologies. For example, their executions are already optimised for different computer architectures, whether it be CPUs, GPUs or new-generation AI processors. In this article, we apply the proposed approach to eigenvalue problems in reactor physics where neutron transport is described by diffusion theory. For a fuel assembly benchmark, we demonstrate that the solution obtained from our new approach is the same (accounting for solver tolerances) as that obtained from the same discretisation coded in a standard way using Fortran. We then proceed to solve a reactor core benchmark using the new approach.
△ Less
Submitted 24 January, 2023;
originally announced January 2023.
-
Effective field theory analysis of the Coulomb breakup of the one-neutron halo nucleus 19C
Authors:
Pierre Capel,
Daniel R. Phillips,
Andrew Andis,
Mirko Bagnarol,
Behnaz Behzadmoghaddam,
Francesca Bonaiti,
Rishabh Bubna,
Ylenia Capitani,
Pierre-Yves Duerinck,
Victoria Durant,
Niklas Döpper,
Aya El Boustani,
Roland Farrell,
Maurus Geiger,
Michael Gennari,
Nitzan Goldberg,
Jakub Herko,
Tanja Kirchner,
Live-Palm Kubushishi,
Zhen Li,
Simone S. Li Muli,
Alexander Long,
Brady Martin,
Kamyar Mohseni,
Imane Moumene
, et al. (7 additional authors not shown)
Abstract:
We analyse the Coulomb breakup of 19C measured at 67A MeV at RIKEN. We use the Coulomb-Corrected Eikonal (CCE) approximation to model the reaction and describe the one-neutron halo nucleus 19C within Halo Effective Field Theory (EFT). At leading order we obtain a fair reproduction of the measured cross section as a function of energy and angle. The description is insensitive to the choice of optic…
▽ More
We analyse the Coulomb breakup of 19C measured at 67A MeV at RIKEN. We use the Coulomb-Corrected Eikonal (CCE) approximation to model the reaction and describe the one-neutron halo nucleus 19C within Halo Effective Field Theory (EFT). At leading order we obtain a fair reproduction of the measured cross section as a function of energy and angle. The description is insensitive to the choice of optical potential, as long as it accurately represents the size of 18C. It is also insensitive to the interior of the 19C wave function. Comparison between theory and experiment thus enables us to infer asymptotic properties of the ground state of 19C: these data put constraints on the one-neutron separation energy of this nucleus and, for a given binding energy, can be used to extract an asymptotic normalisation coefficient (ANC). These results are confirmed by CCE calculations employing next-to-leading order Halo EFT descriptions of 19C: at this order the results for the Coulomb breakup cross section are completely insensitive to the choice of the regulator. Accordingly, this reaction can be used to constrain the one-neutron separation energy and ANC of 19C.
△ Less
Submitted 28 October, 2023; v1 submitted 16 January, 2023;
originally announced January 2023.
-
Forward Orthogonal Deviations GMM and the Absence of Large Sample Bias
Authors:
Robert F. Phillips
Abstract:
It is well-known that generalized method of moments (GMM) estimators of dynamic panel data models can have asymptotic bias if the number of time periods (T) and the number of cross-sectional units (n) are both large (Alvarez and Arellano, 2003). This conclusion, however, follows when all available instruments are used. This paper provides results supporting a more optimistic conclusion when less t…
▽ More
It is well-known that generalized method of moments (GMM) estimators of dynamic panel data models can have asymptotic bias if the number of time periods (T) and the number of cross-sectional units (n) are both large (Alvarez and Arellano, 2003). This conclusion, however, follows when all available instruments are used. This paper provides results supporting a more optimistic conclusion when less than all available instrumental variables are used. If the number of instruments used per period increases with T sufficiently slowly, the bias of GMM estimators based on the forward orthogonal deviations transformation (FOD-GMM) disappears as T and n increase, regardless of the relative rate of increase in T and n. Monte Carlo evidence is provided that corroborates this claim. Moreover, a large-n, large-T distribution result is provided for FOD-GMM.
△ Less
Submitted 28 December, 2022;
originally announced December 2022.
-
Elastic scattering of 3He+4He with SONIK
Authors:
S. N. Paneru,
C. R. Brune,
D. Connolly,
D. Odell,
M. Poudel,
D. R. Phillips,
J. Karpesky,
B. Davids,
C. Ruiz,
A. Lennarz,
U. Greife,
M. Alcorta,
R. Giri,
M. Lovely,
M. Bowry,
M. Delgado,
N. E. Esker,
A. Garnsworthy,
C. Seeman,
P. Machule,
J. Fallis,
A. A. Chen,
F. Laddaran,
A. Firmino,
C. Weinerman
Abstract:
Measurements of the elastic scattering cross section of 3He and 4He are important in order to improve constraints on theoretical models of 4He(3He,g)7Be, a key reaction in Big Bang nucleosynthesis and solar neutrino production. The astrophysical S-factor for this reaction is a significant source of uncertainty in the standard solar-model prediction of the 7Be and 8B solar neutrino fluxes. The elas…
▽ More
Measurements of the elastic scattering cross section of 3He and 4He are important in order to improve constraints on theoretical models of 4He(3He,g)7Be, a key reaction in Big Bang nucleosynthesis and solar neutrino production. The astrophysical S-factor for this reaction is a significant source of uncertainty in the standard solar-model prediction of the 7Be and 8B solar neutrino fluxes. The elastic scattering measurements reported in the literature do not extend to low energies and lack proper uncertainty quantification. A new measurement of the 4He(3He,3He)4He reaction has been made at center-of-mass energies Ec.m. = 0.38-3.13 MeV using the Scattering of Nuclei in Inverse Kinematics (SONIK) scattering chamber: a windowless, extended gas target surrounded by an array of 30 collimated silicon charged particle detectors situated at TRIUMF. This is the first elastic scattering measurement of 3He+4He made below 500 keV and it has greater angular range and better precision than previous measurements. The elastic scattering data were analyzed using both R-matrix and Halo Effective Field Theory (Halo EFT) frameworks, and values of the s-wave scattering length and effective range were extracted. The resulting improvement in knowledge of the s-wave effective-range function at low energies will reduce the overall uncertainty in S34 at solar energies.
△ Less
Submitted 26 November, 2022;
originally announced November 2022.
-
Long Range Plan: Dense matter theory for heavy-ion collisions and neutron stars
Authors:
Alessandro Lovato,
Travis Dore,
Robert D. Pisarski,
Bjoern Schenke,
Katerina Chatziioannou,
Jocelyn S. Read,
Philippe Landry,
Pawel Danielewicz,
Dean Lee,
Scott Pratt,
Fabian Rennecke,
Hannah Elfner,
Veronica Dexheimer,
Rajesh Kumar,
Michael Strickland,
Johannes Jahan,
Claudia Ratti,
Volodymyr Vovchenko,
Mikhail Stephanov,
Dekrayat Almaalol,
Gordon Baym,
Mauricio Hippert,
Jacquelyn Noronha-Hostler,
Jorge Noronha,
Enrico Speranza
, et al. (39 additional authors not shown)
Abstract:
Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theo…
▽ More
Since the release of the 2015 Long Range Plan in Nuclear Physics, major events have occurred that reshaped our understanding of quantum chromodynamics (QCD) and nuclear matter at large densities, in and out of equilibrium. The US nuclear community has an opportunity to capitalize on advances in astrophysical observations and nuclear experiments and engage in an interdisciplinary effort in the theory of dense baryonic matter that connects low- and high-energy nuclear physics, astrophysics, gravitational waves physics, and data science
△ Less
Submitted 7 November, 2022; v1 submitted 3 November, 2022;
originally announced November 2022.
-
Coherently averaged dual-comb spectroscopy with a low-noise and high-power free-running gigahertz dual-comb laser
Authors:
C. R. Phillips,
B. Willenberg,
A. Nussbaum-Lapping,
F. Callegari,
S. L. Camenzind,
J. Pupeikis,
U. Keller
Abstract:
We present a new type of dual optical frequency comb source capable of scaling applications to high measurement speeds while combining high average power, ultra-low noise operation, and a compact setup. Our approach is based on a diode-pumped solid-state laser cavity which includes an intracavity biprism operated at Brewster angle to generate two spatially-separated modes with highly correlated pr…
▽ More
We present a new type of dual optical frequency comb source capable of scaling applications to high measurement speeds while combining high average power, ultra-low noise operation, and a compact setup. Our approach is based on a diode-pumped solid-state laser cavity which includes an intracavity biprism operated at Brewster angle to generate two spatially-separated modes with highly correlated properties. The 15-cm-long cavity uses an Yb:CALGO crystal and a SESAM as an end mirror to generate more than 3 W average power per comb, below 80 fs pulse duration, a repetition rate of 1.03 GHz, and a continuously tunable repetition rate difference up to 27 kHz. We carefully investigate the coherence properties of the dual-comb by a series of heterodyne measurements, revealing several important features: (1) ultra-low jitter on the uncorrelated part of the timing noise; (2) the radio frequency comb lines of the interferograms are fully resolved in free-running operation; (3) we validate that through a simple measurement of the interferograms we can determine the fluctuations of the phase of all the radio frequency comb lines; (4) this phase information is used in a post-processing routine to perform coherently averaged dual-comb spectroscopy of acetylene (C2H2) over long timescales. Our results represent a powerful and general approach to dual-comb applications by combining low noise and high power operation directly from a highly compact laser oscillator.
△ Less
Submitted 17 April, 2023; v1 submitted 27 October, 2022;
originally announced November 2022.
-
Dynamic Flow Control Through Active Matter Programming Language
Authors:
Fan Yang,
Shichen Liu,
Heun Jin Lee,
Rob Phillips,
Matt Thomson
Abstract:
Cells control fluid flows with a spatial and temporal precision that far exceeds the capabilities of current microfluidic technologies. Cells achieve this superior spatio-temporal control by harnessing dynamic networks of cytoskeleton and motor proteins. Thus, engineering systems to mimic cytoskeletal protein networks could lead to the development of a new, active-matter-powered microfluidic devic…
▽ More
Cells control fluid flows with a spatial and temporal precision that far exceeds the capabilities of current microfluidic technologies. Cells achieve this superior spatio-temporal control by harnessing dynamic networks of cytoskeleton and motor proteins. Thus, engineering systems to mimic cytoskeletal protein networks could lead to the development of a new, active-matter-powered microfluidic device with improved performance over the existing technologies. However, reconstituted motor-microtubule systems conventionally generate chaotic flows and cannot perform useful tasks. Here, we develop an all-optical platform for programming flow fields for transport, separation and mixing of cells and particles using networks of microtubules and motor proteins reconstituted in vitro. We employ mathematical modeling for design optimization, which enables the construction of flow fields that achieves micron-scale transport. We use the platform to demonstrate that active-matter-generated flow fields can probe the extensional rheology of polymers, such as DNA, achieve transport and mixing of beads and human cells, and isolation of human cell clusters. Our findings provide a bio-inspired pathway for programmatically engineering dynamic micron-scale flows and demonstrate the vast potential of active matter systems as an engineering technology.
△ Less
Submitted 5 October, 2023; v1 submitted 26 August, 2022;
originally announced August 2022.
-
Evaluating and improving real-world evidence with Targeted Learning
Authors:
Susan Gruber,
Rachael V. Phillips,
Hana Lee,
John Concato,
Mark van der Laan
Abstract:
Purpose: The Targeted Learning roadmap provides a systematic guide for generating and evaluating real-world evidence (RWE). From a regulatory perspective, RWE arises from diverse sources such as randomized controlled trials that make use of real-world data, observational studies, and other study designs. This paper illustrates a principled approach to assessing the validity and interpretability of…
▽ More
Purpose: The Targeted Learning roadmap provides a systematic guide for generating and evaluating real-world evidence (RWE). From a regulatory perspective, RWE arises from diverse sources such as randomized controlled trials that make use of real-world data, observational studies, and other study designs. This paper illustrates a principled approach to assessing the validity and interpretability of RWE.
Methods: We applied the roadmap to a published observational study of the dose-response association between ritodrine hydrochloride and pulmonary edema among women pregnant with twins in Japan. The goal was to identify barriers to causal effect estimation beyond unmeasured confounding reported by the study's authors, and to explore potential options for overcoming the barriers that robustify results.
Results: Following the roadmap raised issues that led us to formulate alternative causal questions that produced more reliable, interpretable RWE. The process revealed a lack of information in the available data to identify a causal dose-response curve. However, under explicit assumptions the effect of treatment with any amount of ritodrine versus none, albeit a less ambitious parameter, can be estimated from data.
Conclusion: Before RWE can be used in support of clinical and regulatory decision-making, its quality and reliability must be systematically evaluated. The TL roadmap prescribes how to carry out a thorough, transparent, and realistic assessment of RWE. We recommend this approach be a routine part of any decision-making process.
△ Less
Submitted 15 August, 2022;
originally announced August 2022.
-
Final-state interactions and spin structure in $E1$ breakup of $^11$Li in Halo EFT
Authors:
Matthias Göbel,
Bijaya Acharya,
Hans-Werner Hammer,
Daniel R. Phillips
Abstract:
We calculate the $E1$ breakup of the $2n$ halo nucleus $^{11}$Li in Halo Effective Field Theory (Halo EFT) at leading order. In Halo EFT, $^{11}$Li is treated as a three-body system of a $^{9}$Li core and two neutrons. We present a detailed investigation of final-state interactions (FSI) in the neutron-neutron $(nn)$ and neutron-core $(nc)$ channels. We employ Moller operators to formulate an expa…
▽ More
We calculate the $E1$ breakup of the $2n$ halo nucleus $^{11}$Li in Halo Effective Field Theory (Halo EFT) at leading order. In Halo EFT, $^{11}$Li is treated as a three-body system of a $^{9}$Li core and two neutrons. We present a detailed investigation of final-state interactions (FSI) in the neutron-neutron $(nn)$ and neutron-core $(nc)$ channels. We employ Moller operators to formulate an expansion scheme that satisfies the non-energy-weighted cluster sum rule and successively includes higher-order terms in the multiple-scattering series for the FSI. Computing the $E1$ strength up to third order in this scheme, we observe apparent convergence and good agreement with experiment. The neutron-neutron FSI is by far the most important contribution and largely determines the maximum value of the $E1$ distribution. However, inclusion of $nc$ FSI does shift the peak position to slightly lower energies. Moreover, we investigate the sensitivity of the $E1$ response to the spin structure of the neutron-${}^9$Li interaction. We contrast results for an interaction that is the same in the spin-1 and spin-2 channels with one that is only operative in the spin-2 channel, and find that good agreement with experimental data is only obtained if the interaction is present in both spin channels. The latter case is shown to be equivalent to a calculation in which the spin of $^9$Li is neglected.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Rubisco function, evolution, and engineering
Authors:
Noam Prywes,
Naiya R Phillips,
Owen T Tuck,
Luis E Valentin-Alvarado,
David F Savage
Abstract:
Carbon fixation is the process by which CO2 is converted from a gas into biomass. The Calvin Benson Bassham (CBB) cycle is the dominant carbon fixation pathway on earth, driving >99.5% of the ~120 billion tons of carbon that are "fixed" as sugar, by plants, algae and cyanobacteria. The carboxylase enzyme in the CBB, ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco), fixes one CO2 molecule…
▽ More
Carbon fixation is the process by which CO2 is converted from a gas into biomass. The Calvin Benson Bassham (CBB) cycle is the dominant carbon fixation pathway on earth, driving >99.5% of the ~120 billion tons of carbon that are "fixed" as sugar, by plants, algae and cyanobacteria. The carboxylase enzyme in the CBB, ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisco), fixes one CO2 molecule per turn of the cycle. Despite being critical to the assimilation of carbon, rubisco's kinetic rate is not very fast and it is a bottleneck in flux through the pathway. This presents a paradox - why hasn't rubisco evolved to be a better catalyst? Many hypothesize that the catalytic mechanism of rubisco is subject to one or more trade-offs, and that rubisco variants have been optimized for their native physiological environment. Here we review the evolution and biochemistry of rubisco through the lens of structure and mechanism in order to understand what trade-offs limit its improvement. We also review the many attempts to improve rubisco itself and, thereby, promote plant growth.
△ Less
Submitted 21 July, 2022;
originally announced July 2022.
-
Towards Precise and Accurate Calculations of Neutrinoless Double-Beta Decay: Project Scoping Workshop Report
Authors:
V. Cirigliano,
Z. Davoudi,
J. Engel,
R. J. Furnstahl,
G. Hagen,
U. Heinz,
H. Hergert,
M. Horoi,
C. W. Johnson,
A. Lovato,
E. Mereghetti,
W. Nazarewicz,
A. Nicholson,
T. Papenbrock,
S. Pastore,
M. Plumlee,
D. R. Phillips,
P. E. Shanahan,
S. R. Stroberg,
F. Viens,
A. Walker-Loud,
K. A. Wendt,
S. M. Wild
Abstract:
We present the results of a National Science Foundation (NSF) Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, a…
▽ More
We present the results of a National Science Foundation (NSF) Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a road map for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The road map contains specific projects in theoretical and computational physics as well as an uncertainty-quantification plan that employs Bayesian Model Mixing and an analysis of correlations between double-beta-decay rates and other observables. The goal of this program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play.
△ Less
Submitted 3 July, 2022;
originally announced July 2022.
-
Interpolating between small- and large-$g$ expansions using Bayesian Model Mixing
Authors:
A. C. Semposki,
R. J. Furnstahl,
D. R. Phillips
Abstract:
Bayesian Model Mixing (BMM) is a statistical technique that can be used to combine models that are predictive in different input domains into a composite distribution that has improved predictive power over the entire input space. We explore the application of BMM to the mixing of two expansions of a function of a coupling constant $g$ that are valid at small and large values of $g$ respectively.…
▽ More
Bayesian Model Mixing (BMM) is a statistical technique that can be used to combine models that are predictive in different input domains into a composite distribution that has improved predictive power over the entire input space. We explore the application of BMM to the mixing of two expansions of a function of a coupling constant $g$ that are valid at small and large values of $g$ respectively. This type of problem is quite common in nuclear physics, where physical properties are straightforwardly calculable in strong and weak interaction limits or at low and high densities or momentum transfers, but difficult to calculate in between. Interpolation between these limits is often accomplished by a suitable interpolating function, e.g., Padé approximants, but it is then unclear how to quantify the uncertainty of the interpolant. We address this problem in the simple context of the partition function of zero-dimensional $φ^4$ theory, for which the (asymptotic) expansion at small $g$ and the (convergent) expansion at large $g$ are both known. We consider three mixing methods: linear mixture BMM, localized bivariate BMM, and localized multivariate BMM with Gaussian processes. We find that employing a Gaussian process in the intermediate region between the two predictive models leads to the best results of the three methods. The methods and validation strategies we present here should be generalizable to other nuclear physics settings.
△ Less
Submitted 31 October, 2023; v1 submitted 8 June, 2022;
originally announced June 2022.
-
Wildebeest Herds on Rolling Hills: Flocking on Arbitrary Curved Surfaces
Authors:
Christina L. Hueschen,
Alexander R. Dunn,
Rob Phillips
Abstract:
The collective behavior of active agents, whether herds of wildebeest or microscopic actin filaments propelled by molecular motors, is an exciting frontier in biological and soft matter physics. Almost three decades ago, Toner and Tu developed a continuum theory of the collective action of flocks, or herds, that helped launch the modern field of active matter. One challenge faced when applying con…
▽ More
The collective behavior of active agents, whether herds of wildebeest or microscopic actin filaments propelled by molecular motors, is an exciting frontier in biological and soft matter physics. Almost three decades ago, Toner and Tu developed a continuum theory of the collective action of flocks, or herds, that helped launch the modern field of active matter. One challenge faced when applying continuum active matter theories to living phenomena is the complex geometric structure of biological environments. Both macroscopic and microscopic herds move on asymmetric curved surfaces, like undulating grass plains or the surface layers of cells or embryos, which can render problems analytically intractable. In this work, we present a formulation of the Toner-Tu flocking theory that uses the finite element method to solve the governing equations on arbitrary curved surfaces. First, we test the developed formalism and its numerical implementation in channel flow with scattering obstacles and on cylindrical and spherical surfaces, comparing our results to analytical solutions. We then progress to surfaces with arbitrary curvature, moving beyond previously accessible problems to explore herding behavior on a variety of landscapes. Our approach allows the investigation of transients and dynamic solutions not revealed by analytic methods. It also enables versatile incorporation of new geometries and boundary conditions and efficient sweeps of parameter space. Looking forward, the work presented here lays the groundwork for a dialogue between Toner-Tu theory and data on collective motion in biologically-relevant geometries, from drone footage of migrating animal herds to movies of microscopic cytoskeletal flows within cells.
△ Less
Submitted 24 August, 2023; v1 submitted 2 June, 2022;
originally announced June 2022.
-
Targeted learning: Towards a future informed by real-world evidence
Authors:
Susan Gruber,
Rachael V. Phillips,
Hana Lee,
Martin Ho,
John Concato,
Mark J. van der Laan
Abstract:
The 21st Century Cures Act of 2016 includes a provision for the U.S. Food and Drug Administration (FDA) to evaluate the potential use of real-world evidence (RWE) to support new indications for use for previously approved drugs, and to satisfy post-approval study requirements. Extracting reliable evidence from real-world data (RWD) is often complicated by a lack of treatment randomization, potenti…
▽ More
The 21st Century Cures Act of 2016 includes a provision for the U.S. Food and Drug Administration (FDA) to evaluate the potential use of real-world evidence (RWE) to support new indications for use for previously approved drugs, and to satisfy post-approval study requirements. Extracting reliable evidence from real-world data (RWD) is often complicated by a lack of treatment randomization, potential intercurrent events, and informative loss to follow up. Targeted Learning (TL) is a sub-field of statistics that provides a rigorous framework to help address these challenges. The TL Roadmap offers a step-by-step guide to generating valid evidence and assessing its reliability. Following these steps produces an extensive amount of information for assessing whether the study provides reliable scientific evidence in support regulatory decision making. This paper presents two case studies that illustrate the utility of following the roadmap. We use targeted minimum loss-based estimation combined with super learning to estimate causal effects. We also compared these findings with those obtained from an unadjusted analysis, propensity score matching, and inverse probability weighting. Non-parametric sensitivity analyses illuminate how departures from (untestable) causal assumptions would affect point estimates and confidence interval bounds that would impact the substantive conclusion drawn from the study. TL's thorough approach to learning from data provides transparency, allowing trust in RWE to be earned whenever it is warranted.
△ Less
Submitted 13 June, 2022; v1 submitted 17 May, 2022;
originally announced May 2022.
-
Practical considerations for specifying a super learner
Authors:
Rachael V. Phillips,
Mark J. van der Laan,
Hana Lee,
Susan Gruber
Abstract:
Common tasks encountered in epidemiology, including disease incidence estimation and causal inference, rely on predictive modeling. Constructing a predictive model can be thought of as learning a prediction function, i.e., a function that takes as input covariate data and outputs a predicted value. Many strategies for learning these functions from data are available, from parametric regressions to…
▽ More
Common tasks encountered in epidemiology, including disease incidence estimation and causal inference, rely on predictive modeling. Constructing a predictive model can be thought of as learning a prediction function, i.e., a function that takes as input covariate data and outputs a predicted value. Many strategies for learning these functions from data are available, from parametric regressions to machine learning algorithms. It can be challenging to choose an approach, as it is impossible to know in advance which one is the most suitable for a particular dataset and prediction task at hand. The super learner (SL) is an algorithm that alleviates concerns over selecting the one "right" strategy while providing the freedom to consider many of them, such as those recommended by collaborators, used in related research, or specified by subject-matter experts. It is an entirely pre-specified and data-adaptive strategy for predictive modeling. To ensure the SL is well-specified for learning the prediction function, the analyst does need to make a few important choices. In this Education Corner article, we provide step-by-step guidelines for making these choices, walking the reader through each of them and providing intuition along the way. In doing so, we aim to empower the analyst to tailor the SL specification to their prediction task, thereby ensuring their SL performs as well as possible. A flowchart provides a concise, easy-to-follow summary of key suggestions and heuristics, based on our accumulated experience, and guided by theory.
△ Less
Submitted 14 March, 2023; v1 submitted 12 April, 2022;
originally announced April 2022.
-
Modeling and mechanical perturbations reveal how spatially regulated anchorage gives rise to spatially distinct mechanics across the mammalian spindle
Authors:
Pooja Suresh,
Vahe Galstyan,
Rob Phillips,
Sophie Dumont
Abstract:
During cell division, the spindle generates force to move chromosomes. In mammals, microtubule bundles called kinetochore-fibers (k-fibers) attach to and segregate chromosomes. To do so, k-fibers must be robustly anchored to the dynamic spindle. We previously developed microneedle manipulation to mechanically challenge k-fiber anchorage, and observed spatially distinct response features revealing…
▽ More
During cell division, the spindle generates force to move chromosomes. In mammals, microtubule bundles called kinetochore-fibers (k-fibers) attach to and segregate chromosomes. To do so, k-fibers must be robustly anchored to the dynamic spindle. We previously developed microneedle manipulation to mechanically challenge k-fiber anchorage, and observed spatially distinct response features revealing the presence of heterogeneous anchorage (Suresh et al. 2020). How anchorage is precisely spatially regulated, and what forces are necessary and sufficient to recapitulate the k-fiber's response to force remain unclear. Here, we develop a coarse-grained k-fiber model and combine with manipulation experiments to infer underlying anchorage using shape analysis. By systematically testing different anchorage schemes, we find that forces solely at k-fiber ends are sufficient to recapitulate unmanipulated k-fiber shapes, but not manipulated ones for which lateral anchorage over a 3 $μ$m length scale near chromosomes is also essential. Such anchorage robustly preserves k-fiber orientation near chromosomes while allowing pivoting around poles. Anchorage over a shorter length scale cannot robustly restrict pivoting near chromosomes, while anchorage throughout the spindle obstructs pivoting at poles. Together, this work reveals how spatially regulated anchorage gives rise to spatially distinct mechanics in the mammalian spindle, which we propose are key for function.
△ Less
Submitted 9 April, 2022;
originally announced April 2022.
-
Spatially multiplexed single-cavity dual-comb laser
Authors:
J. Pupeikis,
B. Willenberg,
S. L. Camenzind,
A. Benayad,
P. Camy,
C. R. Phillips,
U. Keller
Abstract:
Single-cavity dual-comb lasers are a new class of ultrafast lasers which have a wide possible application space including pump-probe sampling, optical ranging, and gas absorption spectroscopy. However, to this date laser cavity multiplexing usually came to the trade-off in laser performance or relative timing noise suppression. We present a new method for multiplexing a single laser cavity to supp…
▽ More
Single-cavity dual-comb lasers are a new class of ultrafast lasers which have a wide possible application space including pump-probe sampling, optical ranging, and gas absorption spectroscopy. However, to this date laser cavity multiplexing usually came to the trade-off in laser performance or relative timing noise suppression. We present a new method for multiplexing a single laser cavity to support a pair of noise-correlated modes. These modes share all intracavity components and take a near-common path, but do not overlap on any active elements. We implement the method with an 80-MHz laser delivering more than 2.4 Watts of average power per comb with sub-140 fs pulses. We reach sub-cycle relative timing jitter of 2.2 fs [20 Hz, 100 kHz]. With this new multiplexing technique, we could implement slow feedback on the repetition rate difference Δfrep, enabling this quantity to be drift-free, low-jitter, and adjustable - a key combination for practical applications that was lacking in prior single-cavity dual-comb systems.
△ Less
Submitted 27 June, 2022; v1 submitted 16 March, 2022;
originally announced March 2022.
-
Analyzing rotational bands in odd-mass nuclei using Effective Field Theory and Bayesian methods
Authors:
I. K. Alnamlah,
E. A. Coello Pérez,
D. R. Phillips
Abstract:
We recently developed an Effective Field Theory (EFT) for rotational bands in odd-mass nuclei. Here we use EFT expressions to perform a Bayesian analysis of data on the rotational energy levels of $^{99}$Tc, ${}^{155,157}$Gd, ${}^{159}$Dy, ${}^{167, 169}$Er, ${}^{167, 169}$Tm, ${}^{183}$W, ${}^{235}$U and ${}^{239}$Pu. The error model in our Bayesian analysis includes both experimental and EFT tru…
▽ More
We recently developed an Effective Field Theory (EFT) for rotational bands in odd-mass nuclei. Here we use EFT expressions to perform a Bayesian analysis of data on the rotational energy levels of $^{99}$Tc, ${}^{155,157}$Gd, ${}^{159}$Dy, ${}^{167, 169}$Er, ${}^{167, 169}$Tm, ${}^{183}$W, ${}^{235}$U and ${}^{239}$Pu. The error model in our Bayesian analysis includes both experimental and EFT truncation uncertainties. It also accounts for the fact that low-energy constants (LECs) at even and odd orders are expected to have different sizes. We use Markov Chain Monte Carlo (MCMC) sampling to explore the joint posterior of the EFT and error-model parameters and show both the LECs and the expansion parameter, $Q$, can be reliably determined. We extract the LECs up to fourth order in the EFT and find that, provided we correctly account for EFT truncation errors in our likelihood, results for lower-order LECs are stable as we go to higher orders. LEC results are also stable with respect to the addition of higher-energy data. We extract the expansion parameter for all the nuclei listed above and find a clear correlation between the extracted and the expected $Q$ based on the single-particle and vibrational energy scales. However, the $Q$ that actually determines the convergence of the EFT expansion is markedly smaller than would be naively expected based on those scales.
△ Less
Submitted 3 March, 2022;
originally announced March 2022.
-
Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives
Authors:
Ingo Tews,
Zohreh Davoudi,
Andreas Ekström,
Jason D. Holt,
Kevin Becker,
Raúl Briceño,
David J. Dean,
William Detmold,
Christian Drischler,
Thomas Duguet,
Evgeny Epelbaum,
Ashot Gasparyan,
Jambul Gegelia,
Jeremy R. Green,
Harald W. Grießhammer,
Andrew D. Hanlon,
Matthias Heinz,
Heiko Hergert,
Martin Hoferichter,
Marc Illa,
David Kekejian,
Alejandro Kievsky,
Sebastian König,
Hermann Krebs,
Kristina D. Launey
, et al. (20 additional authors not shown)
Abstract:
This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges…
▽ More
This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program.
△ Less
Submitted 2 February, 2022;
originally announced February 2022.
-
Performing Bayesian analyses with AZURE2 using BRICK: an application to the ${}^7$Be system
Authors:
Daniel Odell,
Carl R. Brune,
Daniel R. Phillips,
Richard James deBoer,
Som Nath Paneru
Abstract:
Phenomenological $R$-matrix has been a standard framework for the evaluation of resolved resonance cross section data in nuclear physics for many years. It is a powerful method for comparing different types of experimental nuclear data and combining the results of many different experimental measurements in order to gain a better estimation of the true underlying cross sections. Yet a practical ch…
▽ More
Phenomenological $R$-matrix has been a standard framework for the evaluation of resolved resonance cross section data in nuclear physics for many years. It is a powerful method for comparing different types of experimental nuclear data and combining the results of many different experimental measurements in order to gain a better estimation of the true underlying cross sections. Yet a practical challenge has always been the estimation of the uncertainty on both the cross sections at the energies of interest and the fit parameters, which can take the form of standard level parameters. Frequentist ($χ^2$-based) estimation has been the norm. In this work, a Markov Chain Monte Carlo sampler, \texttt{emcee}, has been implemented for the $R$-matrix code \texttt{AZURE2}, creating the Bayesian $R$-matrix Inference Code Kit (\texttt{BRICK}). Bayesian uncertainty estimation has then been carried out for a simultaneous $R$-matrix fit of the $^3$He$(α,γ)^7$Be and $^3$He$(α,α)^3$He reactions in order to gain further insight into the fitting of capture and scattering data. Both data sets constrain the values of the bound state $α$-particle asymptotic normalization coefficients in $^7$Be. The analysis highlights the need for low-energy scattering data with well-documented uncertainty information and shows how misleading results can be obtained in its absence.
△ Less
Submitted 6 October, 2023; v1 submitted 23 December, 2021;
originally announced December 2021.
-
Timing jitter characterization of free-running dual-comb laser with sub-attosecond resolution using optical heterodyne detection
Authors:
Sandro L. Camenzind,
Daniel Koenen,
Benjamin Willenberg,
Justinas Pupeikis,
Christopher R. Phillips,
Ursula Keller
Abstract:
Pulse trains emitted from dual-comb systems are designed to have low relative timing jitter, making them useful for many optical measurement techniques such as optical ranging and spectroscopy. However, the characterization of low-jitter dual-comb systems is challenging because it requires measurement techniques with high sensitivity. Motivated by this challenge, we developed a technique based on…
▽ More
Pulse trains emitted from dual-comb systems are designed to have low relative timing jitter, making them useful for many optical measurement techniques such as optical ranging and spectroscopy. However, the characterization of low-jitter dual-comb systems is challenging because it requires measurement techniques with high sensitivity. Motivated by this challenge, we developed a technique based on an optical heterodyne detection approach for measuring the relative timing jitter of two pulse trains. The method is suitable for dual-comb systems with essentially any repetition rate difference. Furthermore, the proposed approach allows for continuous and precise tracking of the sampling rate. To demonstrate the technique, we perform a detailed characterization of a single-mode-diode pumped $\mathrm{Yb:CaF_2}$ dual-comb laser from a free-running polarization-multiplexed cavity. This new laser produces 115 fs pulses at 160 MHz repetition rate, with 130 mW of average power in each comb. The detection noise floor for the relative timing jitter between the two pulse trains reaches $8.0 \times 10^{-7} \, \mathrm{fs}^2/\mathrm{Hz} \; ( \sim 896 \: \mathrm{zs}/\sqrt{\mathrm{Hz}} )$, and the relative root mean square (rms) timing jitter is 13 fs when integrating from 100 Hz to 1 MHz. This performance indicates that the demonstrated laser is highly compatible with practical dual-comb spectroscopy, ranging, and sampling applications. Furthermore, our results show that the relative timing noise measurement technique can characterize dual-comb systems operating in free-running mode or with finite repetition rate differences while providing a sub-attosecond resolution, which was not feasible with any other approach before.
△ Less
Submitted 21 January, 2022; v1 submitted 11 November, 2021;
originally announced November 2021.
-
Effective Field Theory analysis of ${}^3$He-$α$ scattering data
Authors:
Maheshwor Poudel,
Daniel R. Phillips
Abstract:
We treat low-energy $^3$He-$α$ elastic scattering in an Effective Field Theory (EFT) that exploits the separation of scales in this reaction. We compute the amplitude up to Next-to-Next-to-Leading Order (NNLO), developing a hierarchy of the effective-range parameters that contribute at various orders. We use the resulting formalism to analyze data for recent measurements at center-of-mass energies…
▽ More
We treat low-energy $^3$He-$α$ elastic scattering in an Effective Field Theory (EFT) that exploits the separation of scales in this reaction. We compute the amplitude up to Next-to-Next-to-Leading Order (NNLO), developing a hierarchy of the effective-range parameters that contribute at various orders. We use the resulting formalism to analyze data for recent measurements at center-of-mass energies of 0.38-3.12 MeV using the SONIK gas target at TRIUMF as well as older data in this energy regime. We employ a likelihood function that incorporates the theoretical uncertainty due to truncation of the EFT and use Markov Chain Monte Carlo sampling to obtain the resulting posterior probability distribution. We find that the inclusion of a small amount of data on the analysing power $A_y$ is crucial to determine the sign of the p-wave splitting in such an analysis. The combination of $A_y$ and SONIK data constrains all effective-range parameters up to $O(p^4)$ in both s- and p-waves quite well. The ANCs and s-wave scattering length are consistent with a recent EFT analysis of the capture reaction ${}^3$He($α$,$γ$)${}^7$Be.
△ Less
Submitted 4 July, 2022; v1 submitted 4 October, 2021;
originally announced October 2021.
-
Evaluating the Robustness of Targeted Maximum Likelihood Estimators via Realistic Simulations in Nutrition Intervention Trials
Authors:
Haodong Li,
Sonali Rosete,
Jeremy Coyle,
Rachael V. Phillips,
Nima S. Hejazi,
Ivana Malenica,
Benjamin F. Arnold,
Jade Benjamin-Chung,
Andrew Mertens,
John M. Colford Jr,
Mark J. van der Laan,
Alan E. Hubbard
Abstract:
Several recently developed methods have the potential to harness machine learning in the pursuit of target quantities inspired by causal inference, including inverse weighting, doubly robust estimating equations and substitution estimators like targeted maximum likelihood estimation. There are even more recent augmentations of these procedures that can increase robustness, by adding a layer of cro…
▽ More
Several recently developed methods have the potential to harness machine learning in the pursuit of target quantities inspired by causal inference, including inverse weighting, doubly robust estimating equations and substitution estimators like targeted maximum likelihood estimation. There are even more recent augmentations of these procedures that can increase robustness, by adding a layer of cross-validation (cross-validated targeted maximum likelihood estimation and double machine learning, as applied to substitution and estimating equation approaches, respectively). While these methods have been evaluated individually on simulated and experimental data sets, a comprehensive analysis of their performance across ``real-world'' simulations have yet to be conducted.
In this work, we benchmark multiple widely used methods for estimation of the average treatment effect using ten different nutrition intervention studies data. A realistic set of simulations, based on a novel method, highly adaptive lasso, for estimating the data-generating distribution that guarantees a certain level of complexity (undersmoothing) is used to better mimic the complexity of the true data-generating distribution. We have applied this novel method for estimating the data-generating distribution by individual study and to subsequently use these fits to simulate data and estimate treatment effects parameters as well as their standard errors and resulting confidence intervals. Based on the analytic results, a general recommendation is put forth for use of the cross-validated variants of both substitution and estimating equation estimators. We conclude that the additional layer of cross-validation helps in avoiding unintentional over-fitting of nuisance parameter functionals and leads to more robust inferences.
△ Less
Submitted 28 September, 2021;
originally announced September 2021.
-
Personalized Online Machine Learning
Authors:
Ivana Malenica,
Rachael V. Phillips,
Romain Pirracchio,
Antoine Chambaz,
Alan Hubbard,
Mark J. van der Laan
Abstract:
In this work, we introduce the Personalized Online Super Learner (POSL) -- an online ensembling algorithm for streaming data whose optimization procedure accommodates varying degrees of personalization. Namely, POSL optimizes predictions with respect to baseline covariates, so personalization can vary from completely individualized (i.e., optimization with respect to baseline covariate subject ID)…
▽ More
In this work, we introduce the Personalized Online Super Learner (POSL) -- an online ensembling algorithm for streaming data whose optimization procedure accommodates varying degrees of personalization. Namely, POSL optimizes predictions with respect to baseline covariates, so personalization can vary from completely individualized (i.e., optimization with respect to baseline covariate subject ID) to many individuals (i.e., optimization with respect to common baseline covariates). As an online algorithm, POSL learns in real-time. POSL can leverage a diversity of candidate algorithms, including online algorithms with different training and update times, fixed algorithms that are never updated during the procedure, pooled algorithms that learn from many individuals' time-series, and individualized algorithms that learn from within a single time-series. POSL's ensembling of this hybrid of base learning strategies depends on the amount of data collected, the stationarity of the time-series, and the mutual characteristics of a group of time-series. In essence, POSL decides whether to learn across samples, through time, or both, based on the underlying (unknown) structure in the data. For a wide range of simulations that reflect realistic forecasting scenarios, and in a medical data application, we examine the performance of POSL relative to other current ensembling and online learning methods. We show that POSL is able to provide reliable predictions for time-series data and adjust to changing data-generating environments. We further cultivate POSL's practicality by extending it to settings where time-series enter/exit dynamically over chronological time.
△ Less
Submitted 21 September, 2021;
originally announced September 2021.
-
Picosecond ultrasonics with a free-running dual-comb laser
Authors:
Justinas Pupeikis,
Benjamin Willenberg,
Francois Bruno,
Mike Hettich,
Alexander Nussbaum-Lapping,
Matthias Golling,
Carolin P. Bauer,
Sandro L. Camenzind,
Abdelmjid Benayad,
Patrice Camy,
Bertrand Audoin,
Christopher R. Phillips,
Ursula Keller
Abstract:
We present a free-running 80-MHz dual-comb polarization-multiplexed solid-state laser which delivers 1.8 W of average power with 110-fs pulse duration per comb. With a high-sensitivity pump-probe setup, we apply this free-running dual-comb laser to picosecond ultrasonic measurements. The ultrasonic signatures in a semiconductor multi-quantum-well structure originating from the quantum wells and su…
▽ More
We present a free-running 80-MHz dual-comb polarization-multiplexed solid-state laser which delivers 1.8 W of average power with 110-fs pulse duration per comb. With a high-sensitivity pump-probe setup, we apply this free-running dual-comb laser to picosecond ultrasonic measurements. The ultrasonic signatures in a semiconductor multi-quantum-well structure originating from the quantum wells and superlattice regions are revealed and discussed. We further demonstrate ultrasonic measurements on a thin-film metalized sample and compare these measurements to ones obtained with a pair of locked femtosecond lasers. Our data show that a free-running dual-comb laser is well-suited for picosecond ultrasonic measurements and thus it offers a significant reduction in complexity and cost for this widely adopted non-destructive testing technique.
△ Less
Submitted 1 September, 2021;
originally announced September 2021.
-
The Number of Optimal Strategies in the Penney-Ante Game
Authors:
Reed Phillips,
A. J. Hildebrand
Abstract:
In the Penney-Ante game, Player I chooses a head/tail string of a predetermined length $n\ge3$. Player II, upon seeing Player I's choice, chooses another head/tail string of the same length. A coin is then tossed repeatedly and the player whose string appears first in the resulting head/tail sequence wins the game. The Penney-Ante game has gained notoriety as a source of counterintuitive probabili…
▽ More
In the Penney-Ante game, Player I chooses a head/tail string of a predetermined length $n\ge3$. Player II, upon seeing Player I's choice, chooses another head/tail string of the same length. A coin is then tossed repeatedly and the player whose string appears first in the resulting head/tail sequence wins the game. The Penney-Ante game has gained notoriety as a source of counterintuitive probabilities and nontransitivity phenomena. For example, Player II can always choose a string that beats the choice of Player I in the sense of being more likely to appear first in a random head/tail sequence.
It is known that Player II has a unique optimal strategy that maximizes her winning chances in this game. On the other hand, for Player I there exist multiple equivalent optimal strategies. In this paper we investigate the number, $c_n$, of optimal strategies for Player I, i.e., the number of head/tail strings of length $n$ that maximize the winning probability for Player I assuming optimal play by Player II. We derive a recurrence relation for $c_n$ and use this to obtain a sharp asymptotic estimate for $c_n$. In particular, we show that, as $n\to\infty$, a fixed proportion $α\approx 0.04062\dots$ of the $2^n$ head/tail strings of length $n$ are optimal from Player I's perspective.
△ Less
Submitted 14 July, 2021;
originally announced July 2021.
-
What hath Weinberg wrought? Reflections on what Weinberg's papers on 'Nuclear Forces from Chiral Lagrangians' did and did not accomplish
Authors:
Daniel R. Phillips
Abstract:
I discuss selected legacies of Weinberg's application of chiral Lagrangians to nuclear physics: (1) the use of the chiral expansion to organize the interaction of pions and photons with a nucleus; (2) the much-debated question of why and how the potential derived from a chiral Lagrangian should be inserted in the Schrodinger equation; (3) the emergence of "pionless EFT" as a tool for diagnosing un…
▽ More
I discuss selected legacies of Weinberg's application of chiral Lagrangians to nuclear physics: (1) the use of the chiral expansion to organize the interaction of pions and photons with a nucleus; (2) the much-debated question of why and how the potential derived from a chiral Lagrangian should be inserted in the Schrodinger equation; (3) the emergence of "pionless EFT" as a tool for diagnosing universal correlations that are present in quantum few-body systems of very different sizes, and, perhaps most important of all, (4) an epistemological shift in what is expected of a nuclear-physics calculation.
△ Less
Submitted 11 October, 2021; v1 submitted 7 July, 2021;
originally announced July 2021.
-
Does Bayesian Model Averaging improve polynomial extrapolations? Two toy problems as tests
Authors:
M. A. Connell,
I. Billig,
D. R. Phillips
Abstract:
We assess the accuracy of Bayesian polynomial extrapolations from small parameter values, x, to large values of x. We consider a set of polynomials of fixed order, intended as a proxy for a fixed-order effective field theory (EFT) description of data. We employ Bayesian Model Averaging (BMA) to combine results from different order polynomials (EFT orders). Our study considers two "toy problems" wh…
▽ More
We assess the accuracy of Bayesian polynomial extrapolations from small parameter values, x, to large values of x. We consider a set of polynomials of fixed order, intended as a proxy for a fixed-order effective field theory (EFT) description of data. We employ Bayesian Model Averaging (BMA) to combine results from different order polynomials (EFT orders). Our study considers two "toy problems" where the underlying function used to generate data sets is known. We use Bayesian parameter estimation to extract the polynomial coefficients that describe these data at low x. A "naturalness" prior is imposed on the coefficients, so that they are O(1). We Bayesian-Model-Average different polynomial degrees by weighting each according to its Bayesian evidence and compare the predictive performance of this Bayesian Model Average with that of the individual polynomials. The credibility intervals on the BMA forecast have the stated coverage properties more consistently than does the highest evidence polynomial, though BMA does not necessarily outperform every polynomial.
△ Less
Submitted 15 June, 2022; v1 submitted 10 June, 2021;
originally announced June 2021.
-
Dynamics of ramping bursts in a respiratory neuron model
Authors:
Muhammad U. Abdulla,
Ryan S. Phillips,
Jonathan E. Rubin
Abstract:
Intensive computational and theoretical work has led to the development of mutliple mathematical models for bursting in respiratory neurons in the pre-Bötzinger Complex (pre-BötC) of the mammalian brainstem. Nonetheless, these previous models have not captured the pre-inspiratory ramping aspects of these neurons' activity patterns, in which relatively slow tonic spiking gradually progresses to fas…
▽ More
Intensive computational and theoretical work has led to the development of mutliple mathematical models for bursting in respiratory neurons in the pre-Bötzinger Complex (pre-BötC) of the mammalian brainstem. Nonetheless, these previous models have not captured the pre-inspiratory ramping aspects of these neurons' activity patterns, in which relatively slow tonic spiking gradually progresses to faster spiking and a full blown burst, with a corresponding gradual development of an underlying plateau potential. In this work, we show that the incorporation of the dynamics of the extracellular potassium ion concentration into an existing model for pre-BötC neuron bursting, along with some parameter updates, suffices to induce this ramping behavior. Using fast-slow decomposition, we show that this activity can be considered as a form of parabolic bursting, but with burst termination at a homoclinic bifurcation rather than as a SNIC bifurcation. We also investigate the parameter-dependence of these solutions and show that the proposed model yields a greater dynamic range of burst frequencies, durations, and duty cycles than those produced by other models in the literature.
△ Less
Submitted 7 May, 2021;
originally announced May 2021.
-
Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables
Authors:
S. Wesolowski,
I. Svensson,
A. Ekström,
C. Forssén,
R. J. Furnstahl,
J. A. Melendez,
D. R. Phillips
Abstract:
We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory ($χ$EFT) that are provided by bound-state observables in the $A=3$ and $A=4$ sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the $χ$EFT expansion at next-to-next-to-leading order. A consistent solution…
▽ More
We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory ($χ$EFT) that are provided by bound-state observables in the $A=3$ and $A=4$ sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the $χ$EFT expansion at next-to-next-to-leading order. A consistent solution for the ${}^3$H binding energy, the ${}^4$He binding energy and radius, and the ${}^3$H $β$-decay rate can only be obtained if $χ$EFT truncation errors are included in the analysis. All of these except the $β$-decay rate give essentially degenerate constraints on the 3NF low-energy constants, so it is crucial for estimating these parameters. We use eigenvector continuation for fast and accurate emulation of No-Core Shell Model calculations of the considered few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the hyperparameters that quantify the truncation error. We find a $χ$EFT expansion parameter of $Q=0.33 \pm 0.06$ for these observables.
△ Less
Submitted 11 January, 2022; v1 submitted 9 April, 2021;
originally announced April 2021.
-
One for One, or All for All: Equilibria and Optimality of Collaboration in Federated Learning
Authors:
Avrim Blum,
Nika Haghtalab,
Richard Lanas Phillips,
Han Shao
Abstract:
In recent years, federated learning has been embraced as an approach for bringing about collaboration across large populations of learning agents. However, little is known about how collaboration protocols should take agents' incentives into account when allocating individual resources for communal learning in order to maintain such collaborations. Inspired by game theoretic notions, this paper in…
▽ More
In recent years, federated learning has been embraced as an approach for bringing about collaboration across large populations of learning agents. However, little is known about how collaboration protocols should take agents' incentives into account when allocating individual resources for communal learning in order to maintain such collaborations. Inspired by game theoretic notions, this paper introduces a framework for incentive-aware learning and data sharing in federated learning. Our stable and envy-free equilibria capture notions of collaboration in the presence of agents interested in meeting their learning objectives while keeping their own sample collection burden low. For example, in an envy-free equilibrium, no agent would wish to swap their sampling burden with any other agent and in a stable equilibrium, no agent would wish to unilaterally reduce their sampling burden.
In addition to formalizing this framework, our contributions include characterizing the structural properties of such equilibria, proving when they exist, and showing how they can be computed. Furthermore, we compare the sample complexity of incentive-aware collaboration with that of optimal collaboration when one ignores agents' incentives.
△ Less
Submitted 4 March, 2021;
originally announced March 2021.
-
Neutron-neutron scattering length from the $^6$He$(p,pα)nn$ reaction
Authors:
Matthias Göbel,
Thomas Aumann,
Carlos A. Bertulani,
Tobias Frederico,
Hans-Werner Hammer,
Daniel R. Phillips
Abstract:
We propose a novel method to measure the neutron-neutron scattering length using the $^{6}$He$(p,pα)nn$ reaction in inverse kinematics at high energies. The method is based on the final state interaction (FSI) between the neutrons after the sudden knockout of the $α$ particle. We show that the details of the neutron-neutron relative energy distribution allow for a precise extraction of the $s$-wav…
▽ More
We propose a novel method to measure the neutron-neutron scattering length using the $^{6}$He$(p,pα)nn$ reaction in inverse kinematics at high energies. The method is based on the final state interaction (FSI) between the neutrons after the sudden knockout of the $α$ particle. We show that the details of the neutron-neutron relative energy distribution allow for a precise extraction of the $s$-wave scattering length. We present the state-of-the-art in regard to the theory of this distribution. The distribution is calculated in two steps. First, we calculate the ground-state wave function of $^6$He as a $αn n$ three-body system. For this purpose we use Halo effective field theory (Halo EFT), which also provides uncertainty estimates for the results. We compare our results at this stage to model calculations done with the computer code FaCE. In a second step we determine the effects of the $nn$ FSI using the $nn$ t-matrix. We compare these FSI results to approximate FSI approaches based on standard FSI enhancement factors. While the final distribution is sensitive to the $nn$ scattering length, it depends only weakly on the effective range. Throughout we emphasize the impact of theoretical uncertainties on the neutron-neutron relative energy distribution, and discuss the extent to which those uncertainties limit the extraction of the neutron-neutron scattering length from the reaction $^{6}$He$(p,pα)nn$.
△ Less
Submitted 12 July, 2021; v1 submitted 4 March, 2021;
originally announced March 2021.
-
Generation of Radio Frequency Radiation by Femtosecond Filaments
Authors:
Travis Garrett,
Jennifer Elle,
Michael White,
Remington Reid,
Alexander Englesbe,
Ryan Phillips,
Peter Mardahl,
Erin Thornton,
James Wymer,
Anna Janicek,
Oliver Sale,
Andreas Schmitt-Sody
Abstract:
Recent experiments have shown that femtosecond filamentation plasmas generate ultra-broadband radio frequency radiation (RF). We show that a combination of plasma dynamics is responsible for the RF: a plasma wake field develops behind the laser pulse, and this wake excites (and copropagates with) a surface wave on the plasma column. The surface wave proceeds to detach from the end of the plasma an…
▽ More
Recent experiments have shown that femtosecond filamentation plasmas generate ultra-broadband radio frequency radiation (RF). We show that a combination of plasma dynamics is responsible for the RF: a plasma wake field develops behind the laser pulse, and this wake excites (and copropagates with) a surface wave on the plasma column. The surface wave proceeds to detach from the end of the plasma and propagates forward as the RF pulse. We have developed a four stage model of these plasma wake surface waves and find that it accurately predicts the RF from a wide range of experiments, including both 800 nm and 3.9 $μ$m laser systems.
△ Less
Submitted 17 December, 2021; v1 submitted 24 February, 2021;
originally announced February 2021.