Source author record

Eve Armstrong

Eve Armstrong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Inference of bipolar neutrino flavor oscillations near a core-collapse supernova, based on multiple measurements at Earth

Neutrinos in compact-object environments, such as core-collapse supernovae, can experience various kinds of collective effects in flavor space, engendered by neutrino-neutrino interactions. These include "bipolar" collective oscillations, which are exhibited by neutrino ensembles where different flavors dominate at different energies. Considering the importance of neutrinos in the dynamics and nucleosynthesis in these environments, it is desirable to ascertain whether an Earth-based detection could contain signatures of bipolar oscillations that occurred within a supernova envelope. To that end, we continue examining a cost-function formulation of statistical data assimilation (SDA) to infer solutions to a small-scale model of neutrino flavor transformation. SDA is an inference paradigm designed to optimize a model with sparse data. Our model consists of two mono-energetic neutrino beams with different energies emanating from a source and coherently interacting with each other and with a matter background, with time-varying interaction strengths. We attempt to infer flavor transformation histories of these beams using simulated measurements of the flavor content at locations in vacuum (that is, far from the source), which could in principle correspond to earth-based detectors. Within the scope of this small-scale model, we found that: (i) based on such measurements, the SDA procedure is able to infer \textit{whether} bipolar oscillations had occurred within the protoneutron star envelope, and (ii) if the measurements are able to sample the full amplitude of the neutrino oscillations in vacuum, then the amplitude of the prior bipolar oscillations is also well predicted. This result intimates that the inference paradigm can well complement numerical integration codes, via its ability to infer flavor evolution at physically inaccessible locations.

preprint2022arXiv

The Three Little Pigs and the Big Bad Wolf: Case Studies of Peer Review

I present for your appraisal three independent cases of the manuscript referee process conducted by a venerable peer-reviewed scientific journal. Each case involves a little pig, who submitted for consideration a theoretical plan for a house to be constructed presently, in a faraway land. An anonymous big bad wolf was assigned by the journal to assess the merit of these manuscripts. The pigs proposed three distinct construction frameworks, which varied in physical and mathematical sophistication. The first little pig submitted a model of straw, based on the numerical method of toe-counting. His design included odd features, such as spilled millet and cloven-hoofprints on the window sill -- possibly a ploy to distract the wolf from the manuscript's facile mathematical foundation. The second little pig used a more advanced approach, employing Newton's classical laws of motion, to propose a house of sticks. This pig included in her manuscript copious citations by a specific wolf, possibly aiming to ensure acceptance by flattering the wolf whom she anticipated would be the referee. The third little pig described an ostentatious house of bricks based on an elaborate dynamical systems and stability analysis, possibly scheming to dazzle and impress. The big bad wolf did not appear moved by any of the pigs' tactics. His recommendations were, for straw: the minor revision of water-proofing; for sticks: the major revision of fire-proofing, given concerns surrounding climate change; for bricks: unequivocal rejection, accompanied by multiple derogatory comments regarding "high-and-mighty theorists." I describe each case in detail, and suggest that the wolf's reports might be driven as much by self interest as the manuscripts themselves -- namely, that at the time the wolf wrote the reviews, he was rather hungry. Finally, I examine morals learned, if any.

preprint2021arXiv

Inference solves a boundary-value collision problem, with relevance to neutrino flavor transformation

Understanding neutrino flavor transformation in dense environments such as core-collapse supernovae (CCSN) is critical for inferring nucleosynthesis and interpreting a detected neutrino signal. The role of direction-changing collisions in shaping the neutrino flavor field in these environments is important and poorly understood; it has not been treated self-consistently. There has been progress, via numerical integration, to include the effects of collisions in the dynamics of the neutrino flavor field. While this has led to important insights, integration is limited by its requirement that full initial conditions must be assumed known. On the contrary, feedback from collisions to the neutrino field is a boundary value problem. Numerical integration techniques are poorly equipped to handle that formulation. This paper demonstrates that an inference formulation of the problem can solve a simple collision-only model representing a CCSN core -- without full knowledge of initial conditions. Rather, the procedure solves a two-point boundary value problem with partial information at the bounds. The model is sufficiently simple that physical reasoning may be used as a confidence check on the inference-based solution, and the procedure recovers the expected model dynamics. This result demonstrates that inference can solve a problem that is artificially hidden from integration techniques -- a problem that is an important feature of flavor evolution in dense environments. Thus, it is worthwhile to explore means of augmenting the existing powerful integration tools with inference-based approaches.

preprint2020arXiv

An Artificially-intelligent Means to Escape Discreetly from the Departmental Holiday Party; guide for the socially awkward

We employ simulated annealing to identify the global solution of a dynamical model, to make a favorable impression upon colleagues at the departmental holiday party and then exit undetected as soon as possible. The procedure, Gradual Freeze-out of an Optimal Estimation via Optimization of Parameter Quantification - GFOOEOPQ, is designed for the socially awkward. The socially awkward among us possess little instinct for pulling off such a maneuver, and may benefit from a machine to do it for us. The method rests upon Bayes' Theorem, where the probability of a future model state depends on current knowledge of the model. Here, model state vectors are party attendees, and the future event of interest is their disposition toward us at times following the party. We want these dispositions to be favorable. To this end, we first interact so as to make favorable impressions, or at least ensure that these people remember having seen us there. Then we identify the exit that minimizes the chance that anyone notes how early we high-tailed it. Now, poorly-resolved estimates will correspond to degenerate solutions. As noted, we possess no instinct to identify a global optimum by ourselves. This can have disastrous consequences. For this reason, GFOOEOPQ employs annealing to iteratively home in on this optimum. The method is illustrated via a simulated event hosted by someone in the physics department (I am not sure who), in a two-bedroom apartment on the fifth floor of an elevator building in Manhattan, with viable Exit parameters: front door, side door to a stairwell, fire escape, and a bathroom window that opens onto the fire escape. Preliminary tests are reported at two real social celebrations. The procedure is generalizable to corporate events and family gatherings. Readers are encouraged to report novel applications of GFOOEOPQ, to expand the algorithm.

preprint2020arXiv

Identifying the measurements required to estimate rates of COVID-19 transmission, infection, and detection, using variational data assimilation

We demonstrate the ability of statistical data assimilation to identify the measurements required for accurate state and parameter estimation in an epidemiological model for the novel coronavirus disease COVID-19. Our context is an effort to inform policy regarding social behavior, to mitigate strain on hospital capacity. The model unknowns are taken to be: the time-varying transmission rate, the fraction of exposed cases that require hospitalization, and the time-varying detection probabilities of new asymptomatic and symptomatic cases. In simulations, we obtain accurate estimates of undetected (that is, unmeasured) infectious populations, by measuring the detected cases together with the recovered and dead - and without assumed knowledge of the detection rates. Given a noiseless measurement of the recovered population, excellent estimates of all quantities are obtained using a temporal baseline of 101 days, with the exception of the time-varying transmission rate at times prior to the implementation of social distancing. With low noise added to the recovered population, accurate state estimates require a lengthening of the temporal baseline of measurements. Estimates of all parameters are sensitive to the contamination, highlighting the need for accurate and uniform methods of reporting. The aim of this paper is to exemplify the power of SDA to determine what properties of measurements will yield estimates of unknown parameters to a desired precision, in a model with the complexity required to capture important features of the COVID-19 pandemic.

preprint2020arXiv

Inference offers a metric to constrain dynamical models of neutrino flavor transformation

The multi-messenger astrophysics of compact objects presents a vast range of environments where neutrino flavor transformation may occur and may be important for nucleosynthesis, dynamics, and a detected neutrino signal. Development of efficient techniques for surveying flavor evolution solution spaces in these diverse environments, which augment and complement existing sophisticated computational tools, could leverage progress in this field. To this end we continue our exploration of statistical data assimilation (SDA) to identify solutions to a small-scale model of neutrino flavor transformation. SDA is a machine learning (ML) formula wherein a dynamical model is assumed to generate any measured quantities. Specifically, we use an optimization formulation of SDA wherein a cost function is extremized via the variational method. Regions of state space in which the extremization identifies the global minimum of the cost function will correspond to parameter regimes in which a model solution can exist. Our example study seeks to infer the flavor transformation histories of two mono-energetic neutrino beams coherently interacting with each other and with a matter background. We require that the solution be consistent with measured neutrino flavor fluxes at the point of detection, and with constraints placed upon the flavor content at various locations along their trajectories, such as the point of emission, and the locations of the Mikheyev-Smirnov-Wolfenstein (MSW) resonances. We show how the procedure efficiently identifies solution regimes and rules out regimes where solutions are infeasible. Overall, results intimate the promise of this "variational annealing" methodology to efficiently probe an array of fundamental questions that traditional numerical simulation codes render difficult to access.

preprint2019arXiv

Statistical data assimilation for estimating electrophysiology simultaneously with connectivity within a biological neuronal network

A method of data assimilation (DA) is employed to estimate electrophysiological parameters of neurons simultaneously with their synaptic connectivity in a small model biological network. The DA procedure is cast as an optimization, with a cost function consisting of both a measurement error and a model error term. An iterative reweighting of these terms permits a systematic method to identify the lowest minimum, within a local region of state space, on the surface of a non-convex cost function. In the model, two sets of parameter values are associated with two particular functional modes of network activity: simultaneous firing of all neurons, and a pattern-generating mode wherein the neurons burst in sequence. The DA procedure is able to recover these modes if: i) the stimulating electrical currents have chaotic waveforms, and ii) the measurements consist of the membrane voltages of all neurons in the circuit. Further, this method is able to prune a model of unnecessarily high dimensionality to a representation that contains the maximum dimensionality required to reproduce the provided measurements. This paper offers a proof-of-concept that DA has the potential to inform laboratory designs for estimating properties in small and isolatable functional circuits.

preprint2016arXiv

HVC Interneuron Properties from Statistical Data Assimilation

Data assimilation (DA) solves the inverse problem of inferring initial conditions given data and a model. Here we use biophysically motivated Hodgkin-Huxley (HH) models of avian HVCI neurons, experimentally obtained recordings of these neurons, and our data assimilation algorithm to infer the full set of parameters and a minimal set of ionic currents precisely reproducing the observed waveform information. We find many distinct validated sets of parameters selected by our DA method and choice of model. We conclude exploring variations on the inverse problem applied to neurons producing accurate or inaccurate results; by manipulating data presented to the algorithm, varying sample rate and waveform; and by manipulating the model by adding and subtracting ionic currents.

preprint2016arXiv

Model of the Songbird Nucleus HVC as a Network of Central Pattern Generators

We propose a functional architecture of the adult songbird nucleus HVC in which the core element is a "functional syllable unit" (FSU). In this model, HVC is organized into FSUs, each of which provides the basis for the production of one syllable in vocalization. Within each FSU, the inhibitory neuron population takes one of two operational states: (A) simultaneous firing wherein all inhibitory neurons fire simultaneously, and (B) competitive firing of the inhibitory neurons. Switching between these basic modes of activity is accomplished via changes in the synaptic strengths among the inhibitory neurons. The inhibitory neurons connect to excitatory projection neurons such that during state (A) the activity of projection neurons is suppressed, while during state (B) patterns of sequential firing of projection neurons can occur. The latter state is stabilized by feedback from the projection to the inhibitory neurons. Song composition for specific species is distinguished by the manner in which different FSUs are functionally connected to each other. Ours is a computational model built with biophysically based neurons. We illustrate that many observations of HVC activity are explained by the dynamics of the proposed population of FSUs, and we identify aspects of the model that are currently testable experimentally. In addition, and standing apart from the core features of an FSU, we propose that the transition between modes may be governed by the biophysical mechanism of neuromodulation.

preprint2016arXiv

Nonlinear Statistical Data Assimilation for HVC$_{\text{RA}}$ Neurons in the Avian Song System

With the goal of building a model of the HVC nucleus in the avian song system, we discuss in detail a model of HVC$_{\text{RA}}$ projection neurons comprised of a somatic compartment with fast Na$^+$ and K$^+$ currents and a dendritic compartment with slower Ca$^{2+}$ dynamics. We show this model qualitatively exhibits many observed electrophysiological behaviors. We then show in numerical procedures how one can design and analyze feasible laboratory experiments that allow the estimation of all of the many parameters and unmeasured dynamical variables, given observations of the somatic voltage $V_s(t)$ alone. A key to this procedure is to initially estimate the slow dynamics associated with Ca, blocking the fast Na and K variations, and then with the Ca parameters fixed, estimate the fast Na and K dynamics. This separation of time scales provides a numerically robust method for completing the full neuron model, and the efficacy of the method is tested by prediction when observations are complete. The simulation provides a framework for the slice preparation experiments and illustrates the use of data assimilation methods for the design of those experiments.

preprint2016arXiv

Pipe-cleaner Model of Neuronal Network Dynamics

We present a functional model of neuronal network connectivity in which the single architectural element is the object commonly known in handicraft circles as a pipe cleaner. We argue that the dual nature of a neuronal circuit - that it be at times highly robust to external manipulation and yet sufficiently flexible to allow for learning and adaptation - is embodied in the pipe cleaner, and thus that a pipe cleaner framework serves as an instructive scaffold in which to examine network dynamics. Regarding the dynamics themselves: as pipe cleaners possess no intrinsic dynamics, in our model we attribute the emergent circuit dynamics to magic. Magic is a strategy that has been largely neglected in the neuroscience community, and may serve as an illuminating comparison to the common physics-based approaches. This model makes predictions that it would be really awesome to test experimentally. Moreover, the relative simplicity of the pipe cleaner - setting aside the fact that it comes in an overwhelming variety of colors - renders it an excellent theoretical building block with which to create simple network models. Also, they are incredibly cheap when bought wholesale on Amazon.

preprint2011arXiv

The New Eclipsing Cataclysmic Variable SDSS 154453+2553

The cataclysmic variable SDSS154453+2553 was recently identified in the Sloan Digital Sky Survey. We obtained spectra and photometry at the MDM Observatory, which revealed an eclipse with a 6.03 hour period. The Hα emission line exhibits a strong rotational disturbance during eclipse, indicating that it arises in an accretion disk. A contribution from an M-type companion is also observed. Time-series photometry during eclipse gives an ephemeris of 2454878.0062(15) + 0.251282(2)E. We present spectroscopy through the orbit and eclipse photometry. Our analysis of the secondary star indicates a distance of 800 {\pm} 180 pc.