Showing 1–46 of 46 results for author: Wheeler, A

The impact of incorrect dissociation energies on inferred photospheric abundances

Authors: Sarah E Aquilina, Andrew R Casey, Adam J Wheeler

Abstract: Spectral synthesis codes are essential for inferring stellar parameters and detailed chemical abundances. These codes require many physical inputs to predict an emergent spectrum. Developers adopt the best measurements of those inputs at the time they release their code, but those measurements usually improve over time faster than the software is updated. In general, the impact of using incorrect… ▽ More Spectral synthesis codes are essential for inferring stellar parameters and detailed chemical abundances. These codes require many physical inputs to predict an emergent spectrum. Developers adopt the best measurements of those inputs at the time they release their code, but those measurements usually improve over time faster than the software is updated. In general, the impact of using incorrect or uncertain dissociation energies are largely unknown. Here we evaluate how incorrect dissociation energies impact abundances measured from C2, CN, CH, TiO, and MgO features. For each molecule we synthesised optical spectra of FGKM-type main-sequence and giant stars using the literature dissociation energy, and an incorrect (perturbed) dissociation energy. We find that the uncertainties in the dissociation energies adopted by spectral synthesis codes for CN, CH, TiO, and MgO lead to negligible differences in flux or abundances. C2 is the only diatomic molecule where the uncertainty of the inputted dissociation energy translates to a significant difference in flux, and carbon abundance differences of up to 0.2 dex. For Solar-like stars, the impact on carbon abundance is up to 0.09 dex. These large abundance differences demonstrate the importance of updating the inputs adopted by spectral synthesis codes, as well as a consensus on appropriate values between different codes. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 7 pages, 4 figures

Journal ref: MNRAS 531(4), 2024, 4538-4544

Linear stability analysis for a system of singular amplitude equations arising in biomorphology

Authors: Aric Wheeler, Kevin Zumbrun

Abstract: We study linear stability of exponential periodic solutions of a system of singular amplitude equations associated with convective Turing bifurcation in the presence of conservation laws, as arises in modern biomorphology models, binary fluids, and elsewhere. Consisting of a complex Ginzburg-Landau equation coupled with a singular convection-diffusion equation in "mean modes" associated with conse… ▽ More We study linear stability of exponential periodic solutions of a system of singular amplitude equations associated with convective Turing bifurcation in the presence of conservation laws, as arises in modern biomorphology models, binary fluids, and elsewhere. Consisting of a complex Ginzburg-Landau equation coupled with a singular convection-diffusion equation in "mean modes" associated with conservation laws, these were shown previously by the authors to admit a constant-coefficient linearized stability analysis as in the classical Ginzburg-Landau case -- albeit now singular in wave amplitude epsilon -- yielding useful necessary conditions for stability, both of the exponential functions as solutions of the amplitude equations, and of the associated periodic pattern solving the underlying PDE. Here, we show by a delicate two-parameter matrix perturbation analysis that (strict) satisfaction of these necessary conditions is also sufficient for diffusive stability in the sense of Schneider, yielding a corresponding result, and nonlinear stability, for the underlying PDE. Moreover, we show that they may be interpreted as stability along a non-normally hyperbolic slow manifold approximated by Darcy-type reduction, together with attraction along transverse mean modes, connecting with finite-time approximation theorems of Hacker-Schneider-Zimmerman. △ Less

Submitted 7 June, 2024; originally announced June 2024.

MSC Class: 35

Reducing the Barriers to Entry for Foundation Model Training

Authors: Paolo Faraboschi, Ellis Giles, Justin Hotard, Konstanty Owczarek, Andrew Wheeler

Abstract: The world has recently witnessed an unprecedented acceleration in demands for Machine Learning and Artificial Intelligence applications. This spike in demand has imposed tremendous strain on the underlying technology stack in supply chain, GPU-accelerated hardware, software, datacenter power density, and energy consumption. If left on the current technological trajectory, future demands show insur… ▽ More The world has recently witnessed an unprecedented acceleration in demands for Machine Learning and Artificial Intelligence applications. This spike in demand has imposed tremendous strain on the underlying technology stack in supply chain, GPU-accelerated hardware, software, datacenter power density, and energy consumption. If left on the current technological trajectory, future demands show insurmountable spending trends, further limiting market players, stifling innovation, and widening the technology gap. To address these challenges, we propose a fundamental change in the AI training infrastructure throughout the technology ecosystem. The changes require advancements in supercomputing and novel AI training approaches, from high-end software to low-level hardware, microprocessor, and chip design, while advancing the energy efficiency required by a sustainable infrastructure. This paper presents the analytical framework that quantitatively highlights the challenges and points to the opportunities to reduce the barriers to entry for training large language models. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Chemical Cartography with APOGEE: Two-process Parameters and Residual Abundances for 288,789 Stars from Data Release 17

Authors: Tawny Sit, David H. Weinberg, Adam Wheeler, Christian R. Hayes, Sten Hasselquist, Thomas Masseron, Jennifer Sobeck

Abstract: Stellar abundance measurements are subject to systematic errors that induce extra scatter and artificial correlations in elemental abundance patterns. We derive empirical calibration offsets to remove systematic trends with surface gravity $\log(g)$ in 17 elemental abundances of 288,789 evolved stars from the SDSS APOGEE survey. We fit these corrected abundances as the sum of a prompt process trac… ▽ More Stellar abundance measurements are subject to systematic errors that induce extra scatter and artificial correlations in elemental abundance patterns. We derive empirical calibration offsets to remove systematic trends with surface gravity $\log(g)$ in 17 elemental abundances of 288,789 evolved stars from the SDSS APOGEE survey. We fit these corrected abundances as the sum of a prompt process tracing core-collapse supernovae and a delayed process tracing Type Ia supernovae, thus recasting each star's measurements into the amplitudes $A_{\text{cc}}$ and $A_{\text{Ia}}$ and the element-by-element residuals from this two-parameter fit. As a first application of this catalog, which is $8\times$ larger than that of previous analyses that used a restricted $\log(g)$ range, we examine the median residual abundances of 14 open clusters, nine globular clusters, and four dwarf satellite galaxies. Relative to field Milky Way disk stars, the open clusters younger than 2 Gyr show $\approx 0.1-0.2$ dex enhancements of the neutron-capture element Ce, and the two clusters younger than 0.5 Gyr also show elevated levels of C+N, Na, S, and Cu. Globular clusters show elevated median abundances of C+N, Na, Al, and Ce, and correlated abundance residuals that follow previously known trends. The four dwarf satellites show similar residual abundance patterns despite their different star formation histories, with $\approx 0.2-0.3$ dex depletions in C+N, Na, and Al and $\approx 0.1$ dex depletions in Ni, V, Mn, and Co. We provide our catalog of corrected APOGEE abundances, two-process amplitudes, and residual abundances, which will be valuable for future studies of abundance patterns in different stellar populations and of additional enrichment processes that affect galactic chemical evolution. △ Less

Submitted 24 May, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Updated to ApJ accepted version. 35 pages, 18 figures, 5 tables. Catalog available at https://doi.org/10.5281/zenodo.10659205

Scalable Agent-Based Modeling for Complex Financial Market Simulations

Authors: Aaron Wheeler, Jeffrey D. Varner

Abstract: In this study, we developed a computational framework for simulating large-scale agent-based financial markets. Our platform supports trading multiple simultaneous assets and leverages distributed computing to scale the number and complexity of simulated agents. Heterogeneous agents make decisions in parallel, and their orders are processed through a realistic, continuous double auction matching e… ▽ More In this study, we developed a computational framework for simulating large-scale agent-based financial markets. Our platform supports trading multiple simultaneous assets and leverages distributed computing to scale the number and complexity of simulated agents. Heterogeneous agents make decisions in parallel, and their orders are processed through a realistic, continuous double auction matching engine. We present a baseline model implementation and show that it captures several known statistical properties of real financial markets (i.e., stylized facts). Further, we demonstrate these results without fitting models to historical financial data. Thus, this framework could be used for direct applications such as human-in-the-loop machine learning or to explore theoretically exciting questions about market microstructure's role in forming the statistical regularities of real markets. To the best of our knowledge, this study is the first to implement multiple assets, parallel agent decision-making, a continuous double auction mechanism, and intelligent agent types in a scalable real-time environment. △ Less

Submitted 31 January, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Ensemble variational Monte Carlo for optimization of correlated excited state wave functions

Authors: William A. Wheeler, Kevin G. Kleiner, Lucas K. Wagner

Abstract: Variational Monte Carlo methods have recently been applied to the calculation of excited states; however, it is still an open question what objective function is most effective. A promising approach is to optimize excited states using a penalty to minimize overlap with lower eigenstates, which has the drawback that states must be computed one at a time. We derive a general framework for constructi… ▽ More Variational Monte Carlo methods have recently been applied to the calculation of excited states; however, it is still an open question what objective function is most effective. A promising approach is to optimize excited states using a penalty to minimize overlap with lower eigenstates, which has the drawback that states must be computed one at a time. We derive a general framework for constructing objective functions with minima at the the lowest $N$ eigenstates of a many-body Hamiltonian. The objective function uses a weighted average of the energies and an overlap penalty, which must satisfy several conditions. We show this objective function has a minimum at the exact eigenstates for a finite penalty, and provide a few strategies to minimize the objective function. The method is demonstrated using ab initio variational Monte Carlo to calculate the degenerate first excited state of a CO molecule. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 7 pages, 4 figures

Korg: fitting, model atmosphere interpolation, and Brackett lines

Authors: Adam J Wheeler, Andrew R Casey, Matthew W Abruzzo

Abstract: We describe several updates to Korg, a package for 1D LTE spectral synthesis of FGKM stars. Built-in functions to fit observed spectra via synthesis or equivalent widths make it easy to take advantage of Korg's automatic differentiation. Comparison to a past analysis of 18 Sco shows that we obtain significantly reduced line-to-line abundance scatter with Korg. Fitting and synthesis are facilitated… ▽ More We describe several updates to Korg, a package for 1D LTE spectral synthesis of FGKM stars. Built-in functions to fit observed spectra via synthesis or equivalent widths make it easy to take advantage of Korg's automatic differentiation. Comparison to a past analysis of 18 Sco shows that we obtain significantly reduced line-to-line abundance scatter with Korg. Fitting and synthesis are facilitated by a rigorously-tested model atmosphere interpolation method, which introduces negligible error to synthesized spectra for stars with $T_\mathrm{eff} \gtrsim 4000\,\mathrm{K}$. For cooler stars, atmosphere interpolation is complicated by the presence of molecules, though we demonstrate an adequate method for cool dwarfs. The chemical equilibrium solver has been extended to include polyatomic and charged molecules, extending Korg's regime of applicability to M stars. We also discuss a common oversight regarding the synthesis of hydrogen lines in the infrared, and show that Korg's Brackett line profiles are a much closer match to observations than others available. Documentation, installation instructions, and tutorials are available at https://github.com/ajwheeler/Korg.jl. △ Less

Submitted 9 April, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: updated to reflect published version

Convective Turing bifurcation with conservation laws

Authors: Aric Wheeler, Kevin Zumbrun

Abstract: Generalizing results of \cite{MC,S} and \cite{HSZ} for certain model reaction-diffusion and reaction-convection-diffusion equations, we derive and rigorously justify weakly nonlinear amplitude equations governing general Turing bifurcation in the presence of conservation laws. In the nonconvective, reaction-diffusion case, this is seen similarly as in \cite{MC,S} to be a real Ginsburg-Landau equat… ▽ More Generalizing results of \cite{MC,S} and \cite{HSZ} for certain model reaction-diffusion and reaction-convection-diffusion equations, we derive and rigorously justify weakly nonlinear amplitude equations governing general Turing bifurcation in the presence of conservation laws. In the nonconvective, reaction-diffusion case, this is seen similarly as in \cite{MC,S} to be a real Ginsburg-Landau equation coupled with a diffusion equation in a large-scale mean-mode vector comprising variables associated with conservation laws. In the general, convective case, by contrast, the amplitude equations as noted in \cite{HSZ} consist of a complex Ginsburg-Landau equation coupled with a singular convection-diffusion equation featuring rapidly-propagating modes with speed $\sim 1/\eps$ where $\eps$ measures amplitude of the wave as a disturbance from a background steady state. Different from the partially coupled case considered in \cite{HSZ} in the context of Bénard-Marangoni convection/inclined flow, the Ginzburg Landau and mean-mode equations are here fully coupled, leading to substantial new difficulties in the analysis. Applications are to biological morphogenesis, in particular vasculogenesis, as described by the Murray-Oster and other mechanochemical/hydrodynamical models △ Less

Submitted 25 May, 2023; originally announced May 2023.

Surfing, sweeping, and assembly of particles by a moving liquid crystal phase boundary

Authors: Tom Shneer, Jocelyn Ochoa, Alauna C. Wheeler, Isabella C. Reyes, Chaitanya Joshi, Benjamin J. Stokes, Linda S. Hirst, Timothy J. Atherton

Abstract: Non-equilibrium transport of particles embedded in a liquid crystal host can, by cooling through a phase transition, be exploited to create a remarkable variety of structures including shells, foams, and gels. Due to the complexity of the multicomponent system and protocol-dependent experimental results, the physical mechanisms behind structure selection remain only partially understood. Here we f… ▽ More Non-equilibrium transport of particles embedded in a liquid crystal host can, by cooling through a phase transition, be exploited to create a remarkable variety of structures including shells, foams, and gels. Due to the complexity of the multicomponent system and protocol-dependent experimental results, the physical mechanisms behind structure selection remain only partially understood. Here we formulate a new model coupling LC physics to a Fokker-Planck equation as is commonly used in studies of transport. The resulting model allows us to draw an analogy between the LC-nanocomposite system and chemotaxis, enriching the space of possible target structures that could be produced. We study the model in one dimension both analytically and numerically to identify different parameter regimes where soliton-like pulses of particles ``surf'' the phase boundary or where the interface ``sweeps'' particles from one domain to another. We also consider an extended model that includes agglomeration of the particles and observe formation of periodic structures as a prototypical example of hierarchical self assembly. Results are compared with experimental observations of transport by isolated phase boundaries. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 11 pages, 6 figures

The Eighteenth Data Release of the Sloan Digital Sky Surveys: Targeting and First Spectra from SDSS-V

Authors: Andrés Almeida, Scott F. Anderson, Maria Argudo-Fernández, Carles Badenes, Kat Barger, Jorge K. Barrera-Ballesteros, Chad F. Bender, Erika Benitez, Felipe Besser, Dmitry Bizyaev, Michael R. Blanton, John Bochanski, Jo Bovy, William Nielsen Brandt, Joel R. Brownstein, Johannes Buchner, Esra Bulbul, Joseph N. Burchett, Mariana Cano Díaz, Joleen K. Carlberg, Andrew R. Casey, Vedant Chandra, Brian Cherinka, Cristina Chiappini, Abigail A. Coker , et al. (129 additional authors not shown)

Abstract: The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM),… ▽ More The eighteenth data release of the Sloan Digital Sky Surveys (SDSS) is the first one for SDSS-V, the fifth generation of the survey. SDSS-V comprises three primary scientific programs, or "Mappers": Milky Way Mapper (MWM), Black Hole Mapper (BHM), and Local Volume Mapper (LVM). This data release contains extensive targeting information for the two multi-object spectroscopy programs (MWM and BHM), including input catalogs and selection functions for their numerous scientific objectives. We describe the production of the targeting databases and their calibration- and scientifically-focused components. DR18 also includes ~25,000 new SDSS spectra and supplemental information for X-ray sources identified by eROSITA in its eFEDS field. We present updates to some of the SDSS software pipelines and preview changes anticipated for DR19. We also describe three value-added catalogs (VACs) based on SDSS-IV data that have been published since DR17, and one VAC based on the SDSS-V data in the eFEDS field. △ Less

Submitted 6 July, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: Accepted to ApJS

PyQMC: an all-Python real-space quantum Monte Carlo module in PySCF

Authors: William A. Wheeler, Shivesh Pathak, Kevin Kleiner, Shunyue Yuan, João N. B. Rodrigues, Cooper Lorsung, Kittithat Krongchon, Yueqing Chang, Yiqing Zhou, Brian Busemeyer, Kiel T. Williams, Alexander Muñoz, Chun Yu Chow, Lucas K. Wagner

Abstract: We describe a new open-source Python-based package for high accuracy correlated electron calculations using quantum Monte Carlo (QMC) in real space: PyQMC. PyQMC implements modern versions of QMC algorithms in an accessible format, enabling algorithmic development and easy implementation of complex workflows. Tight integration with the PySCF environment allows for simple comparison between QMC cal… ▽ More We describe a new open-source Python-based package for high accuracy correlated electron calculations using quantum Monte Carlo (QMC) in real space: PyQMC. PyQMC implements modern versions of QMC algorithms in an accessible format, enabling algorithmic development and easy implementation of complex workflows. Tight integration with the PySCF environment allows for simple comparison between QMC calculations and other many-body wave function techniques, as well as access to high accuracy trial wave functions. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Korg: a modern 1D LTE spectral synthesis package

Authors: Adam J. Wheeler, Matthew W. Abruzzo, Andrew R. Casey, Melissa K. Ness

Abstract: We present Korg, a new package for 1D LTE (local thermodynamic equilibrium) spectral synthesis of FGK stars, which computes theoretical spectra from the near-ultraviolet to the near-infrared, and implements both plane-parallel and spherical radiative transfer. We outline the inputs and internals of Korg, and compare synthetic spectra from Korg, MOOG, Turbospectrum, and SME. The disagreements betwe… ▽ More We present Korg, a new package for 1D LTE (local thermodynamic equilibrium) spectral synthesis of FGK stars, which computes theoretical spectra from the near-ultraviolet to the near-infrared, and implements both plane-parallel and spherical radiative transfer. We outline the inputs and internals of Korg, and compare synthetic spectra from Korg, MOOG, Turbospectrum, and SME. The disagreements between Korg and the other codes are no larger than those between the other codes, although disagreement between codes is substantial. We examine the case of a C$_2$ band in detail, finding that uncertainties on physical inputs to spectral synthesis account for a significant fraction of the disagreement. Korg is 1-100 times faster than other codes in typical use, compatible with automatic differentiation libraries, and easily extensible, making it ideal for statistical inference and parameter estimation applied to large data sets. Documentation and installation instructions are available at https://ajwheeler.github.io/Korg.jl/stable/. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: to appear in ApJ; comments welcome

Majda and ZND models for detonation: nonlinear stability vs. formation of singularities

Authors: Paul Blochas, Aric Wheeler

Abstract: For the ZND model, we show also singularity formation on the downstream side for arbitrary exponentially-growing weighted norms. For the Majda model on the other hand, we establish for appropriate such weighted norms a set of energy estimates implying not only non-formation of singularities near waves of arbitrary amplitude but also full asymptotic orbital stability for small-amplitude ones. For the ZND model, we show also singularity formation on the downstream side for arbitrary exponentially-growing weighted norms. For the Majda model on the other hand, we establish for appropriate such weighted norms a set of energy estimates implying not only non-formation of singularities near waves of arbitrary amplitude but also full asymptotic orbital stability for small-amplitude ones. △ Less

Submitted 28 October, 2022; originally announced October 2022.

Comments: 52 pages

MSC Class: 35B35 35B40 35B44 ACM Class: G.1.8

Continuous guts poker and numerical optimization of generalized recursive games

Authors: Kevin Buck, Jae Hwan Lee, Jacob Platnick, Aric Wheeler, Kevin Zumbrun

Abstract: We study a type of generalized recursive game introduced by Castronova, Chen, and Zumbrun featuring increasing stakes, with an emphasis on continuous guts poker and $1$ v. $n$ coalitions. Our main results are to develop practical numerical algorithms with rigorous underlying theory for the approximation of optimal mutiplayer strategies, and to use these to obtain a number of interesting observatio… ▽ More We study a type of generalized recursive game introduced by Castronova, Chen, and Zumbrun featuring increasing stakes, with an emphasis on continuous guts poker and $1$ v. $n$ coalitions. Our main results are to develop practical numerical algorithms with rigorous underlying theory for the approximation of optimal mutiplayer strategies, and to use these to obtain a number of interesting observations about guts. Outcomes are a striking 2-strategy optimum for $n$-player coalitions, with asymptotic advantage approximately $16\%$; convergence of Fictitious Play to symmetric Nash equilibrium; and a malevolent interactive $n$-player "bot" for demonstration. △ Less

Submitted 11 December, 2022; v1 submitted 4 August, 2022; originally announced August 2022.

Toric and tropical Bertini theorems in positive characteristic

Authors: Francesca Gandini, Milena Hering, Diane Maclagan, Fatemeh Mohammadi, Jenna Rajchgot, Ashley K. Wheeler, Josephine Yu

Abstract: We generalize the toric Bertini theorem of Fuchs, Mantova, and Zannier to positive characteristic. A key part of the proof is a new algebraically closed field containing the field \kk(t_1,\dots,t_d) of rational functions over an algebraically closed field \kk of prime characteristic. As a corollary, we extend the tropical Bertini theorem of Maclagan and Yu to arbitrary characteristic, which remove… ▽ More We generalize the toric Bertini theorem of Fuchs, Mantova, and Zannier to positive characteristic. A key part of the proof is a new algebraically closed field containing the field \kk(t_1,\dots,t_d) of rational functions over an algebraically closed field \kk of prime characteristic. As a corollary, we extend the tropical Bertini theorem of Maclagan and Yu to arbitrary characteristic, which removes the characteristic dependence from the d-connectivity result for tropical varieties from that paper. △ Less

Submitted 25 November, 2021; originally announced November 2021.

MSC Class: 14M25 (Primary); 12J25; 13A35; 14T15 (Secondary)

Structure-Property-Performance Relationships of Cuprous Oxide Nanostructures for Dielectric Mie Resonance-Enhanced Photocatalysis

Authors: Ravi Teja A. Tirumala, Sunil Gyawali, Aaron Wheeler, Sundaram Bhardwaj Ramakrishnan, Rishmali Sooriyagoda, Farshid Mohammadparast, Susheng Tan, A. Kaan Kalkan, Alan D. Bristow, Marimuthu Andiappan

Abstract: Nanostructured metal oxides, such as Cu2O, CeO2, α-Fe2O3, and TiO2 can efficiently mediate photocatalysis for solar-to-chemical energy conversion and pollution remediation. In this contribution, we report a novel approach, dielectric Mie resonance-enhanced photocatalysis, to enhance the catalytic activity of metal oxide photocatalysts. Specifically, we demonstrate that Cu2O nanostructures exhibiti… ▽ More Nanostructured metal oxides, such as Cu2O, CeO2, α-Fe2O3, and TiO2 can efficiently mediate photocatalysis for solar-to-chemical energy conversion and pollution remediation. In this contribution, we report a novel approach, dielectric Mie resonance-enhanced photocatalysis, to enhance the catalytic activity of metal oxide photocatalysts. Specifically, we demonstrate that Cu2O nanostructures exhibiting dielectric Mie resonances can exhibit up to an order of magnitude higher photocatalytic rate as compared to Cu2O nanostructures not exhibiting dielectric Mie resonances. Our finite-difference time-domain (FDTD) simulation and experimental results predict a volcano-type relationship between the photocatalytic rate and the size of Cu2O nanospheres and nanocubes. Using transient absorption measurements, we reveal that a coherent electronic process associated with dielectric Mie resonance-mediated charge carrier generation is dominant in Cu2O nanostructures that exhibit higher photocatalytic rates. Although we experimentally demonstrate dielectric Mie resonance-enhanced photocatalysis using Cu2O particles here, based on our FDTD simulations, we anticipate the same can be achieved with other metal oxide photocatalysts, including CeO2, α-Fe2O3, and TiO2. △ Less

Submitted 17 January, 2022; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: Combined draft. Both Main draft & SI. Total of 62 pages. arXiv admin note: text overlap with arXiv:2109.13285

The homogeneity of the star forming environment of the Milky Way disk over time

Authors: Melissa K. Ness, Adam J. Wheeler, Kevin McKinnon, Danny Horta, Andrew R. Casey, Emily C. Cunningham, Adrian M. Price-Whelan

Abstract: Stellar abundances and ages afford the means to link chemical enrichment to galactic formation. In the Milky Way, individual element abundances show tight correlations with age, which vary in slope across ([Fe/H]-[$α$/Fe]). Here, we step from characterising abundances as measures of age, to understanding how abundances trace properties of stellar birth-environment in the disk over time. Using meas… ▽ More Stellar abundances and ages afford the means to link chemical enrichment to galactic formation. In the Milky Way, individual element abundances show tight correlations with age, which vary in slope across ([Fe/H]-[$α$/Fe]). Here, we step from characterising abundances as measures of age, to understanding how abundances trace properties of stellar birth-environment in the disk over time. Using measurements from $\sim$27,000 APOGEE stars (R=22,500, SNR$>$200), we build simple local linear models to predict a sample of elements (X = Si, O, Ca, Ti, Ni, Al, Mn, Cr) using (Fe, Mg) abundances alone, as fiducial tracers of supernovae production channels. Given [Fe/H] and [Mg/H], we predict these elements, [X/H], to about double the uncertainty of their measurements. The intrinsic dispersion, after subtracting measurement errors in quadrature is $\approx 0.015-0.04$~dex. The residuals of the prediction (measurement $-$ model) for each element demonstrate that each element has an individual link to birth properties at fixed (Fe, Mg). Residuals from primarily massive-star supernovae (i.e. Si, O, Al) partially correlate with guiding radius. Residuals from primarily supernovae Ia (i.e. Mn, Ni) partially correlate with age. A fraction of the intrinsic scatter that persists at fixed (Fe, Mg), however, after accounting for correlations, does not appear to further discriminate between birth properties that can be traced with present-day measurements. Presumably, this is because the residuals are also, in part, a measure of the typical (in)-homogeneity of the disk's stellar birth environments, previously inferred only using open-cluster systems. Our study implies at fixed birth radius and time, there is a median scatter of $\approx 0.01-0.015$ dex in elements generated in supernovae sources. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: submitted to ApJ, 16 pages, 12 figures

Chemodynamical signatures of bar resonances in the Galactic disk: current data and future prospects

Authors: Adam Wheeler, Irene Abril-Cabezas, Wilma H. Trick, Francesca Fragkoudi, Melissa Ness

Abstract: The Galactic disk exhibits complex chemical and dynamical substructure thought to be induced by the bar, spiral arms, and satellites. Here, we explore the chemical signatures of bar resonances in action and velocity space and characterize the differences between the signatures of corotation and higher-order resonances using test particle simulations. Thanks to recent surveys, we now have large dat… ▽ More The Galactic disk exhibits complex chemical and dynamical substructure thought to be induced by the bar, spiral arms, and satellites. Here, we explore the chemical signatures of bar resonances in action and velocity space and characterize the differences between the signatures of corotation and higher-order resonances using test particle simulations. Thanks to recent surveys, we now have large datasets containing metallicities and kinematics of stars outside the solar neighborhood. We compare the simulations to the observational data from Gaia EDR3 and LAMOST DR5 and find weak evidence for a slow bar with the "hat" moving group ($250~\text{km/s} \lesssim v_φ\lesssim 270~\text{km/s}$) associated with its outer Lindblad resonance and "Hercules" ($170~\textrm{km/s} \lesssim v_φ\lesssim 195~\text{km/s}$) with corotation. While constraints from current data are limited by their spatial footprint, stars closer in azimuth than the Sun to the bar's minor axis show much stronger signatures of the bar's outer Lindblad and corotation resonances in test particle simulations. Future datasets with greater azimuthal coverage, including the final Gaia data release, will allow reliable chemodynamical identification of bar resonances. △ Less

Submitted 10 August, 2022; v1 submitted 11 May, 2021; originally announced May 2021.

Comments: Published in ApJ

Disentangled Representation Learning for Astronomical Chemical Tagging

Authors: Damien de Mijolla, Melissa Ness, Serena Viti, Adam Wheeler

Abstract: Modern astronomical surveys are observing spectral data for millions of stars. These spectra contain chemical information that can be used to trace the Galaxy's formation and chemical enrichment history. However, extracting the information from spectra, and making precise and accurate chemical abundance measurements are challenging. Here, we present a data-driven method for isolating the chemical… ▽ More Modern astronomical surveys are observing spectral data for millions of stars. These spectra contain chemical information that can be used to trace the Galaxy's formation and chemical enrichment history. However, extracting the information from spectra, and making precise and accurate chemical abundance measurements are challenging. Here, we present a data-driven method for isolating the chemical factors of variation in stellar spectra from those of other parameters (i.e. \teff, \logg, \feh). This enables us to build a spectral projection for each star with these parameters removed. We do this with no ab initio knowledge of elemental abundances themselves, and hence bypass the uncertainties and systematics associated with modeling that rely on synthetic stellar spectra. To remove known non-chemical factors of variation, we develop and implement a neural network architecture that learns a disentangled spectral representation. We simulate our recovery of chemically identical stars using the disentangled spectra in a synthetic APOGEE-like dataset. We show that this recovery declines as a function of the signal to noise ratio, but that our neural network architecture outperforms simpler modeling choices. Our work demonstrates the feasibility of data-driven abundance-free chemical tagging. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: Accepted by ApJ. code available at github.com/drd13/tagging-package

Autonomous object harvesting using synchronized optoelectronic microrobots

Authors: Christopher Bendkowski, Laurent Mennillo, Tao Xu, Mohamed Elsayed, Filip Stojic, Harrison Edwards, Shuailong Zhang, Cindi Morshead, Vijay Pawar, Aaron R. Wheeler, Danail Stoyanov, Michael Shaw

Abstract: Optoelectronic tweezer-driven microrobots (OETdMs) are a versatile micromanipulation technology based on the use of light induced dielectrophoresis to move small dielectric structures (microrobots) across a photoconductive substrate. The microrobots in turn can be used to exert forces on secondary objects and carry out a wide range of micromanipulation operations, including collecting, transportin… ▽ More Optoelectronic tweezer-driven microrobots (OETdMs) are a versatile micromanipulation technology based on the use of light induced dielectrophoresis to move small dielectric structures (microrobots) across a photoconductive substrate. The microrobots in turn can be used to exert forces on secondary objects and carry out a wide range of micromanipulation operations, including collecting, transporting and depositing microscopic cargos. In contrast to alternative (direct) micromanipulation techniques, OETdMs are relatively gentle, making them particularly well suited to interacting with sensitive objects such as biological cells. However, at present such systems are used exclusively under manual control by a human operator. This limits the capacity for simultaneous control of multiple microrobots, reducing both experimental throughput and the possibility of cooperative multi-robot operations. In this article, we describe an approach to automated targeting and path planning to enable open-loop control of multiple microrobots. We demonstrate the performance of the method in practice, using microrobots to simultaneously collect, transport and deposit silica microspheres. Using computational simulations based on real microscopic image data, we investigate the capacity of microrobots to collect target cells from within a dissociated tissue culture. Our results indicate the feasibility of using OETdMs to autonomously carry out micromanipulation tasks within complex, unstructured environments. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 7 pages, 5 figures

Diffusive stability of convective Turing patterns

Authors: Aric Wheeler, Kevin Zumbrun

Abstract: Following the approach of [E1, M1, M2, S1, S2, SZJV] for reaction diffusion systems, we justify rigorously the Eckhaus stability criterion for stability of convective Turing patterns, as derived formally by complex Ginzburg-Landau approximation [SS, NW, WZ]. Notably, our analysis includes also higher-order, nonlocal, and even certain semilinear hyperbolic systems. Following the approach of [E1, M1, M2, S1, S2, SZJV] for reaction diffusion systems, we justify rigorously the Eckhaus stability criterion for stability of convective Turing patterns, as derived formally by complex Ginzburg-Landau approximation [SS, NW, WZ]. Notably, our analysis includes also higher-order, nonlocal, and even certain semilinear hyperbolic systems. △ Less

Submitted 17 June, 2021; v1 submitted 20 January, 2021; originally announced January 2021.

Convective Turing Bifurcation

Authors: Aric Wheeler, Kevin Zumbrun

Abstract: Following the approach pioneered by Eckhaus, Mielke, Schneider, and others for reaction diffusion systems [E, M1, M2, S1, S2, SZJV], we systematically derive formally by multiscale expansion and justify rigorously by Lyapunov-Schmidt reduction amplitude equations describing Turing-type bifurcations of general reaction diffusion convection systems. Notably, our analysis includes also higher-order,… ▽ More Following the approach pioneered by Eckhaus, Mielke, Schneider, and others for reaction diffusion systems [E, M1, M2, S1, S2, SZJV], we systematically derive formally by multiscale expansion and justify rigorously by Lyapunov-Schmidt reduction amplitude equations describing Turing-type bifurcations of general reaction diffusion convection systems. Notably, our analysis includes also higher-order, nonlocal, and even certain semilinear hyperbolic systems. △ Less

Submitted 18 January, 2021; originally announced January 2021.

The Swan: Data-Driven Inference of Stellar Surface Gravities for Cool Stars from Photometric Light Curves

Authors: Maryum Sayeed, Daniel Huber, Adam Wheeler, Melissa Ness

Abstract: Stellar light curves are well known to encode physical stellar properties. Precise, automated and computationally inexpensive methods to derive physical parameters from light curves are needed to cope with the large influx of these data from space-based missions such as Kepler and TESS. Here we present a new methodology which we call The Swan, a fast, generalizable and effective approach for deriv… ▽ More Stellar light curves are well known to encode physical stellar properties. Precise, automated and computationally inexpensive methods to derive physical parameters from light curves are needed to cope with the large influx of these data from space-based missions such as Kepler and TESS. Here we present a new methodology which we call The Swan, a fast, generalizable and effective approach for deriving stellar surface gravity ($\log g$) for main sequence, subgiant and red giant stars from Kepler light curves using local linear regression on the full frequency content of Kepler long cadence power spectra. With this inexpensive data-driven approach, we recover $\log g$ to a precision of $\sim$0.02 dex for 13,822 stars with seismic $\log g$ values between 0.2-4.4 dex, and $\sim$0.11 dex for 4,646 stars with Gaia derived $\log g$ values between 2.3-4.6 dex. We further develop a signal-to-noise metric and find that granulation is difficult to detect in many cool main sequence stars ($T_{\text{eff}}$ $\lesssim$ 5500 K), in particular K dwarfs. By combining our $\log g$ measurements with Gaia radii, we derive empirical masses for 4,646 subgiant and main sequence stars with a median precision of $\sim$7%. Finally, we demonstrate that our method can be used to recover $\log g$ to a similar mean absolute deviation precision for TESS-baseline of 27 days. Our methodology can be readily applied to photometric time-series observations to infer stellar surface gravities to high precision across evolutionary states. △ Less

Submitted 15 February, 2021; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: 21 pages, 14 figures, 4 tables. Accepted for publishing in the Astronomical Journal. The electronic versions of Tables 1, 3 and 4 are available as an ancillary file (sidebar on the right)

Hundreds of new periodic signals detected in the first year of $\it{TESS}$ with the ${\tt weirddetector}$

Authors: Joheen Chakraborty, Adam Wheeler, David Kipping

Abstract: We apply the ${\tt weirddetector}$, a nonparametric signal detection algorithm based on phase dispersion minimization, in a search for low duty-cycle periodic signals in the Transiting Exoplanet Survey Satellite (${\it TESS}$) photometry. Our approach, in contrast to commonly used model-based approaches specifically for flagging transits, eclipsing binaries, or other similarly periodic events, mak… ▽ More We apply the ${\tt weirddetector}$, a nonparametric signal detection algorithm based on phase dispersion minimization, in a search for low duty-cycle periodic signals in the Transiting Exoplanet Survey Satellite (${\it TESS}$) photometry. Our approach, in contrast to commonly used model-based approaches specifically for flagging transits, eclipsing binaries, or other similarly periodic events, makes minimal assumptions about the shape of a periodic signal, with the goal of finding "weird" signals of unexpected or arbitrary shape. In total, 248,301 ${\it TESS}$ sources from the first-year Southern sky survey are run through the ${\tt weirddetector}$, of which we manually inspect the top 21,500 for periodicity. To minimize false-positives, we here only report on the upper decile in terms of signal score, a sample for which we obtain 97% recall of ${\it TESS}$ eclipsing binaries and 62% of the TOIs. In our sample, we find 377 previously unreported periodic signals, for which we make a first-pass assignment that 26 are ultra-short periods ($<0.3$ d), 313 are likely eclipsing binaries, 28 appear planet-like, and 10 are miscellaneous signals. △ Less

Submitted 24 September, 2020; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: Accepted to MNRAS, ref. MN-20-3245-MJ.R1. Data available at https://github.com/joheenc/woi_data . 09/23/20: edited Table 1 and Section 2.1; eliminated several previously known signals from ASAS-SN and VSX catalogs

An unsupervised method for identifying $X$-enriched stars directly from spectra: Li in LAMOST

Authors: Adam Wheeler, Melissa Ness, David W. Hogg

Abstract: Stars with peculiar element abundances are important markers of chemical enrichment mechanisms. We present a simple method, tangent space projection (TSP), for the detection of $X$-enriched stars, for arbitrary elements $X$, even from blended lines. Our method does not require stellar labels, but instead directly estimates the counterfactual unrenriched spectrum from other unlabelled spectra. As a… ▽ More Stars with peculiar element abundances are important markers of chemical enrichment mechanisms. We present a simple method, tangent space projection (TSP), for the detection of $X$-enriched stars, for arbitrary elements $X$, even from blended lines. Our method does not require stellar labels, but instead directly estimates the counterfactual unrenriched spectrum from other unlabelled spectra. As a case study, we apply this method to the $6708~$Å Li doublet in LAMOST DR5, identifying 8,428 Li-enriched stars seamlessly across evolutionary state. We comment on the explanation for Li-enrichement for different subpopulations, including planet accretion, nonstandard mixing, and youth. △ Less

Submitted 8 March, 2021; v1 submitted 8 September, 2020; originally announced September 2020.

Comments: Published in ApJ; Data available at https://doi.org/10.7916/d8-3ap9-qe35

Journal ref: ApJ 908 247 (2021)

Abundances in the Milky Way across five nucleosynthetic channels from 4 million LAMOST stars

Authors: Adam Wheeler, Melissa Ness, Sven Buder, Joss Bland-Hawthorn, Gayandhi De Silva, Michael Hayden, Janez Kos, Geraint F. Lewis, Sarah Martell, Sanjib Sharma, Jeffrey D. Simpson, D. B. Zucker, Thomaž Zwitter

Abstract: Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. W… ▽ More Large stellar surveys are revealing the chemodynamical structure of the Galaxy across a vast spatial extent. However, the many millions of low-resolution spectra observed to date are yet to be fully exploited. We employ The Cannon, a data-driven approach to estimating abundances, to obtain detailed abundances from low-resolution (R = 1800) LAMOST spectra, using the GALAH survey as our reference. We deliver five (for dwarfs) or six (for giants) estimated abundances representing five different nucleosynthetic channels, for 3.9 million stars, to a precision of 0.05 - 0.23 dex. Using wide binary pairs, we demonstrate that our abundance estimates provide chemical discriminating power beyond metallicity alone. We show the coverage of our catalogue with radial, azimuthal and dynamical abundance maps, and examine the neutron capture abundances across the disk and halo, which indicate different origins for the in-situ and accreted halo populations. LAMOST has near-complete Gaia coverage and provides an unprecedented perspective on chemistry across the Milky Way. △ Less

Submitted 11 June, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: 17 pages + refs and appendices, to appear in ApJ. Data: https://doi.org/10.7910/DVN/5VWKMC

Judge, Jury & Encryptioner: Exceptional Device Access with a Social Cost

Authors: Sacha Servan-Schreiber, Archer Wheeler

Abstract: We present Judge, Jury and Encryptioner (JJE) an exceptional access scheme for unlocking devices that does not give unilateral power to any single authority. JJE achieves this by placing final approval to unlock a device in the hands of peer devices. JJE distributes maintenance of the protocol across a network of "custodians" such as courts, government agencies, civil rights watchdogs, and academi… ▽ More We present Judge, Jury and Encryptioner (JJE) an exceptional access scheme for unlocking devices that does not give unilateral power to any single authority. JJE achieves this by placing final approval to unlock a device in the hands of peer devices. JJE distributes maintenance of the protocol across a network of "custodians" such as courts, government agencies, civil rights watchdogs, and academic institutions. Unlock requests, however, can only be approved by a randomly selected set of recently active peer devices that must be physically located by law enforcement in order to gain access to the locked device. This requires that law enforcement expend both human and monetary resources and pay a "social cost" in order to find and request the participation of random device owners in the unlock process. Compared to other proposed exceptional access schemes, we believe that JJE mitigates the risk of mass surveillance, law enforcement abuse, and vulnerability to unlawful attackers. While we propose a concrete construction, our primary goal with JJE is to spur discussion on ethical exceptional access schemes that balance privacy of individuals and the desires for law enforcement. JJE transparently reveals the use of exceptional access to the public and enforces a fixed social cost that, we believe, can be an effective deterrent to mass surveillance and abuse. △ Less

Submitted 6 March, 2020; v1 submitted 11 December, 2019; originally announced December 2019.

Geometric Equations for Matroid Varieties

Authors: Jessica Sidman, Will Traves, Ashley Wheeler

Abstract: Each point $x$ in Gr$(r,n)$ corresponds to an $r \times n$ matrix $A_x$ which gives rise to a matroid $M_x$ on its columns. Gel'fand, Goresky, MacPherson, and Serganova showed that the sets $\{y \in \mathrm{Gr}(r,n) | M_y = M_x\}$ form a stratification of Gr$(r,n)$ with many beautiful properties. However, results of Mnëv and Sturmfels show that these strata can be quite complicated, and in particu… ▽ More Each point $x$ in Gr$(r,n)$ corresponds to an $r \times n$ matrix $A_x$ which gives rise to a matroid $M_x$ on its columns. Gel'fand, Goresky, MacPherson, and Serganova showed that the sets $\{y \in \mathrm{Gr}(r,n) | M_y = M_x\}$ form a stratification of Gr$(r,n)$ with many beautiful properties. However, results of Mnëv and Sturmfels show that these strata can be quite complicated, and in particular may have arbitrary singularities. We study the ideals $I_x$ of matroid varieties, the Zariski closures of these strata. We construct several classes of examples based on theorems from projective geometry and describe how the Grassmann-Cayley algebra may be used to derive non-trivial elements of $I_x$ geometrically when the combinatorics of the matroid is sufficiently rich. △ Less

Submitted 1 July, 2020; v1 submitted 3 August, 2019; originally announced August 2019.

Comments: Updated Proposition 2.1.3. Added Theorem 2.1.7 and Remark 3.0.3

MSC Class: 14M15; 05B35

The Weird Detector: Flagging periodic, coherent signals of arbitrary shape in time series photometry

Authors: Adam Wheeler, David Kipping

Abstract: By design, model-based approaches for flagging transiting exoplanets in light curves, such as boxed least squares, excel at detecting planets with low S/N at the expense of finding signals that are not well described by the assumed model, such as self-lensing binaries, disintegrating or evaporating planets, or planets with large rings. So far, such signals have typically been found through visual… ▽ More By design, model-based approaches for flagging transiting exoplanets in light curves, such as boxed least squares, excel at detecting planets with low S/N at the expense of finding signals that are not well described by the assumed model, such as self-lensing binaries, disintegrating or evaporating planets, or planets with large rings. So far, such signals have typically been found through visual searches by professional or citizen scientists, or by inspection of the photometric power-spectra. We present a nonparametric detection algorithm, for short duty-cycle periodic signals in photometric time series based on phase dispersion minimization. We apply our code to 161,786 Kepler sources and detect 18 new periodic signals consistent with heartbeat binaries/planets, 4 new singly-transiting systems, and 2 new doubly-transiting systems. We show that our code is able to recover the majority of known Kepler objects of interest (KOIs) to high confidence, as well as more unusual events such as Boyajian's star and a comet passing through the Kepler field. Nonparametric signal-flagging techniques, such as the one presented here, will become increasingly valuable with the coming data from TESS and future transit surveys as the volume of data available to us exceeds that which can be feasibly examined manually. △ Less

Submitted 26 March, 2019; v1 submitted 21 March, 2019; originally announced March 2019.

Comments: 13 pages, 11 figures, accepted to MNRAS

Many-Body Electric Multipole Operators in Extended Systems

Authors: William A. Wheeler, Lucas K. Wagner, Taylor L. Hughes

Abstract: The quantum mechanical position operators, and their products, are not well-defined in systems obeying periodic boundary conditions. Here we extend the work of Resta who developed a formalism to calculate the electronic polarization as an expectation value of a many-body operator, to include higher multipole moments, e.g., quadrupole and octupole. We define $n$-th order multipole operators whose e… ▽ More The quantum mechanical position operators, and their products, are not well-defined in systems obeying periodic boundary conditions. Here we extend the work of Resta who developed a formalism to calculate the electronic polarization as an expectation value of a many-body operator, to include higher multipole moments, e.g., quadrupole and octupole. We define $n$-th order multipole operators whose expectation values can be used to calculate the $n$-th multipole moment when all of the lower moments are vanishing (modulo a quantum). We show that changes in our operators are tied to flows of $n-1$-st multipole currents, and encode the adiabatic evolution of the system in the presence of an $n-1$-st gradient of the electric field. Finally, we test our operators on a set of tightbinding models to show that they correctly determine the phase diagrams of topological quadrupole and octupole models, capture an adiabatic quadrupole pump, and distinguish a bulk quadrupole moment from other mechanisms that generate corner charges. △ Less

Submitted 30 September, 2019; v1 submitted 17 December, 2018; originally announced December 2018.

Comments: 12 pages (6 figures), v2 contains updates for clarity and a new discussion linking this work to other recent related work

Journal ref: Phys. Rev. B 100, 245135 (2019)

Geometric optics for Rayleigh wavetrains in d-dimensional nonlinear elasticity

Authors: Aric Wheeler, Mark Williams

Abstract: A Rayleigh wave is a type of surface wave that propagates in the boundary of an elastic solid with traction (or Neumann) boundary conditions. Since the 1980s much work has been done on the problem of constructing a leading term in an \emph{approximate} solution to the rather complicated second-order quasilinear hyperbolic boundary value problem with fully nonlinear Neumann boundary conditions that… ▽ More A Rayleigh wave is a type of surface wave that propagates in the boundary of an elastic solid with traction (or Neumann) boundary conditions. Since the 1980s much work has been done on the problem of constructing a leading term in an \emph{approximate} solution to the rather complicated second-order quasilinear hyperbolic boundary value problem with fully nonlinear Neumann boundary conditions that governs the propagation of Rayleigh waves. The question has remained open whether or not this leading term approximate solution is really close in a precise sense to the exact solution of the governing equations. We prove a positive answer to this question for the case of Rayleigh wavetrains in any space dimension $d\geq 2$. The case of Rayleigh pulses in dimension $d=2$ has already been treated by Coulombel and Williams. For highly oscillatory Rayleigh wavetrains we are able to construct high-order approximate solutions consisting of the leading term plus an arbitrary number of correctors. Using those high-order solutions we then perform an error analysis which shows (among other things) that for small wavelengths the leading term is close to the exact solution in $L^\infty$ on a fixed time interval independent of the wavelength. The error analysis is carried out in a somewhat general setting and is applicable to other types of waves for which high order approximate solutions can be constructed. △ Less

Submitted 31 December, 2017; originally announced January 2018.

The Sandpile Group of a Thick Cycle Graph

Authors: Diane Christine Alar, Jonathan Celaya, Luis David García Puente, Micah Henson, Ashley K. Wheeler

Abstract: The majority of graphs whose sandpile groups are known are either regular or simple. We give an explicit formula for a family of non-regular multi-graphs called thick cycles. A thick cycle graph is a cycle where multi-edges are permitted. Its sandpile group is the direct sum of cyclic groups of orders given by quotients of greatest common divisors of minors of its Laplacian matrix. We show these g… ▽ More The majority of graphs whose sandpile groups are known are either regular or simple. We give an explicit formula for a family of non-regular multi-graphs called thick cycles. A thick cycle graph is a cycle where multi-edges are permitted. Its sandpile group is the direct sum of cyclic groups of orders given by quotients of greatest common divisors of minors of its Laplacian matrix. We show these greatest common divisors can be expressed in terms of monomials in the graph's edge multiplicities. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 12 pages

MSC Class: 05E15 ACM Class: G.2.1; G.2.2

Closed-Loop Policies for Operational Tests of Safety-Critical Systems

Authors: Jeremy Morton, Tim A. Wheeler, Mykel J. Kochenderfer

Abstract: Manufacturers of safety-critical systems must make the case that their product is sufficiently safe for public deployment. Much of this case often relies upon critical event outcomes from real-world testing, requiring manufacturers to be strategic about how they allocate testing resources in order to maximize their chances of demonstrating system safety. This work frames the partially observable a… ▽ More Manufacturers of safety-critical systems must make the case that their product is sufficiently safe for public deployment. Much of this case often relies upon critical event outcomes from real-world testing, requiring manufacturers to be strategic about how they allocate testing resources in order to maximize their chances of demonstrating system safety. This work frames the partially observable and belief-dependent problem of test scheduling as a Markov decision process, which can be solved efficiently to yield closed-loop manufacturer testing policies. By solving for policies over a wide range of problem formulations, we are able to provide high-level guidance for manufacturers and regulators on issues relating to the testing of safety-critical systems. This guidance spans an array of topics, including circumstances under which manufacturers should continue testing despite observed incidents, when manufacturers should test aggressively, and when regulators should increase or reduce the real-world testing requirements for an autonomous vehicle. △ Less

Submitted 19 May, 2018; v1 submitted 25 July, 2017; originally announced July 2017.

Comments: 12 pages, 5 figures, 5 tables

Deep Stochastic Radar Models

Authors: Tim Allan Wheeler, Martin Holder, Hermann Winner, Mykel Kochenderfer

Abstract: Accurate simulation and validation of advanced driver assistance systems requires accurate sensor models. Modeling automotive radar is complicated by effects such as multipath reflections, interference, reflective surfaces, discrete cells, and attenuation. Detailed radar simulations based on physical principles exist but are computationally intractable for realistic automotive scenes. This paper d… ▽ More Accurate simulation and validation of advanced driver assistance systems requires accurate sensor models. Modeling automotive radar is complicated by effects such as multipath reflections, interference, reflective surfaces, discrete cells, and attenuation. Detailed radar simulations based on physical principles exist but are computationally intractable for realistic automotive scenes. This paper describes a methodology for the construction of stochastic automotive radar models based on deep learning with adversarial loss connected to real-world data. The resulting model exhibits fundamental radar effects while remaining real-time capable. △ Less

Submitted 16 June, 2017; v1 submitted 31 January, 2017; originally announced January 2017.

Comments: IEEE Intelligent Vehicles Symposium 2017. Accepted for plenary presentation

Ultra-High Gradient Channeling Acceleration in Nanostructures: Design/Progress of Proof-of-Concept (POC) Experiments

Authors: Y. M. Shin, A. Green, A. H. Lumpkin, R. M. Thurman-Keup, V. Shiltsev, X. Zhang, D. M. -A. Farinella, P. Taborek, T. Tajima, J. A. Wheeler, G. Mourou

Abstract: This paper describes simulation analyses on beam and laser (X-ray)-driven accelerations in effective nanotube models obtained from Vsim and EPOCH codes. Experimental setups to detect wakefields are also outlined with accelerator facilities at Fermilab and NIU. In the FAST facility, the electron beamline was successfully commissioned at 50 MeV and it is being upgraded toward higher energies for ele… ▽ More This paper describes simulation analyses on beam and laser (X-ray)-driven accelerations in effective nanotube models obtained from Vsim and EPOCH codes. Experimental setups to detect wakefields are also outlined with accelerator facilities at Fermilab and NIU. In the FAST facility, the electron beamline was successfully commissioned at 50 MeV and it is being upgraded toward higher energies for electron accelerator R&D. The 50 MeV injector beamline of the facility is used for X-ray crystal-channeling radiation with a diamond target. It has been proposed to utilize the same diamond crystal for a channeling acceleration POC test. Another POC experiment is also designed for the NIU accelerator lab with time-resolved electron diffraction. Recently, a stable generation of single-cycle laser pulses with tens of Petawatt power based on thin film compression (TFC) technique has been investigated for target normal sheath acceleration (TNSA) and radiation pressure acceleration (RPA). The experimental plan with a nanometer foil is discussed with an available test facility such as Extreme Light Infrastructure - Nuclear Physics (ELI-NP). △ Less

Submitted 27 December, 2016; originally announced December 2016.

Comments: 4 pp

Report number: Fermilab-Conf-16-381-TD

Finiteness of Associated Primes of Local Cohomology Modules Over Stanley-Reisner Rings

Authors: Roberto Barrera, Jeffrey Madsen, Ashley K. Wheeler

Abstract: Local cohomology modules, even over a Noetherian ring $R$, are typically unwieldly. As such, it is of interest whether or not they have finitely many associated primes. We prove the affirmative in the case where $R$ is a Stanley-Reisner ring over a field, whose associated simplicial complex is a $T$-space. Local cohomology modules, even over a Noetherian ring $R$, are typically unwieldly. As such, it is of interest whether or not they have finitely many associated primes. We prove the affirmative in the case where $R$ is a Stanley-Reisner ring over a field, whose associated simplicial complex is a $T$-space. △ Less

Submitted 6 September, 2017; v1 submitted 17 September, 2016; originally announced September 2016.

Learning Discrete Bayesian Networks from Continuous Data

Authors: Yi-Chun Chen, Tim Allan Wheeler, Mykel John Kochenderfer

Abstract: Learning Bayesian networks from raw data can help provide insights into the relationships between variables. While real data often contains a mixture of discrete and continuous-valued variables, many Bayesian network structure learning algorithms assume all random variables are discrete. Thus, continuous variables are often discretized when learning a Bayesian network. However, the choice of discr… ▽ More Learning Bayesian networks from raw data can help provide insights into the relationships between variables. While real data often contains a mixture of discrete and continuous-valued variables, many Bayesian network structure learning algorithms assume all random variables are discrete. Thus, continuous variables are often discretized when learning a Bayesian network. However, the choice of discretization policy has significant impact on the accuracy, speed, and interpretability of the resulting models. This paper introduces a principled Bayesian discretization method for continuous variables in Bayesian networks with quadratic complexity instead of the cubic complexity of other standard techniques. Empirical demonstrations show that the proposed method is superior to the established minimum description length algorithm. In addition, this paper shows how to incorporate existing methods into the structure learning process to discretize all continuous variables and simultaneously learn Bayesian network structures. △ Less

Submitted 17 September, 2018; v1 submitted 8 December, 2015; originally announced December 2015.

Journal ref: Journal of Artificial Intelligence Research 59 (2017) 103-132

Spectroscopic super-resolution fluorescence cell imaging using ultra-small Ge quantum dots

Authors: Mingying Song, Ali Karatutlu, Osman Ersoy, Yun Zhou, Yongxin Yang, Yuanpeng Zhang, William R. Little, Ann P. Wheeler, Andrei V. Sapelkin

Abstract: In single molecule localisation super-resolution microscopy the need for repeated image capture limits the imaging speed, while the size of fluorescence probes limits the possible theoretical localisation resolution. Here, we demonstrated a spectral imaging based super-resolution approach by separating the overlapping diffraction spots into several detectors during a single scanning period and tak… ▽ More In single molecule localisation super-resolution microscopy the need for repeated image capture limits the imaging speed, while the size of fluorescence probes limits the possible theoretical localisation resolution. Here, we demonstrated a spectral imaging based super-resolution approach by separating the overlapping diffraction spots into several detectors during a single scanning period and taking advantage of the size-dependent emission wavelength in nanoparticles. This approach has been tested using off-the-shelf quantum dots (Qdot) and in-house novel ultra-small (~3 nm) Ge QDs. Furthermore, we developed a method-specific Gaussian fitting and maximum likelihood estimation based on a Matlab algorithm for fast QDs localisation. We demonstrate that this methodology results in ~ 40 nm localisation resolution using commercial QDs and ~12 nm localisation resolution using Ge QDs. Using a standard scanning confocal microscope we achieved data acquisition rate of 1.6 seconds/frame. However, we show that this approach has a potential to deliver data acquisition rates on ms scale thus providing super-resolution in live systems. △ Less

Submitted 9 April, 2015; v1 submitted 31 March, 2015; originally announced March 2015.

Comments: This paper has been withdrawn by the authors due to an author dispute. This paper will not be replaced or revised

Principal Minor Ideals and Rank Restrictions on their Vanishing Sets

Authors: Ashley K. Wheeler

Abstract: All matrices we consider have entries in a fixed algebraically closed field $K$. A minor of a square matrix is principal means it is defined by the same row and column indices. We study the ideal generated by size $t$ principal minors of a generic matrix, and restrict our attention to locally closed subsets of its vanishing set, given by matrices of a fixed rank. The main result is a computation o… ▽ More All matrices we consider have entries in a fixed algebraically closed field $K$. A minor of a square matrix is principal means it is defined by the same row and column indices. We study the ideal generated by size $t$ principal minors of a generic matrix, and restrict our attention to locally closed subsets of its vanishing set, given by matrices of a fixed rank. The main result is a computation of the dimension of the locally closed set of $n\times n$ rank $n-2$ matrices whose size $n-2$ principal minors vanish; this set has dimension $n^2-n-4$. △ Less

Submitted 3 August, 2015; v1 submitted 19 March, 2015; originally announced March 2015.

Ideals Generated by Principal Minors

Authors: Ashley K. Wheeler

Abstract: A minor is principal means it is defined by the same row and column indices. Let $X$ be a square generic matrix, $K[X]$ the polynomial ring in entries of $X$, over an algebraically closed field, $K$. For fixed $t\leq n$, let $\mathfrak P_t$ denote the ideal generated by the size $t$ principal minors of $X$. When $t=2$ the resulting quotient ring $K[X]/\mathfrak P_2$ is a normal complete intersecti… ▽ More A minor is principal means it is defined by the same row and column indices. Let $X$ be a square generic matrix, $K[X]$ the polynomial ring in entries of $X$, over an algebraically closed field, $K$. For fixed $t\leq n$, let $\mathfrak P_t$ denote the ideal generated by the size $t$ principal minors of $X$. When $t=2$ the resulting quotient ring $K[X]/\mathfrak P_2$ is a normal complete intersection domain. When $t>2$ we break the problem into cases depending on a fixed rank, $r$, of $X$. We show when $r=n$ for any $t$, the respective images of $\mathfrak P_t$ and $\mathfrak P_{n-t}$ in the localized polynomial ring, where we invert $\det X$, are isomorphic. From that we show the algebraic set given by $\mathfrak P_{n-1}$ has a codimension $n$ component, plus a codimension 4 component defined by the determinantal ideal (which is given by all the submaximal minors of $X$). When $n=4$ the two components are linked, and we prove some consequences. △ Less

Submitted 3 August, 2015; v1 submitted 7 October, 2014; originally announced October 2014.

Countering Trusting Trust through Diverse Double-Compiling

Authors: David A. Wheeler

Abstract: An Air Force evaluation of Multics, and Ken Thompson's famous Turing award lecture "Reflections on Trusting Trust," showed that compilers can be subverted to insert malicious Trojan horses into critical software, including themselves. If this attack goes undetected, even complete analysis of a system's source code will not find the malicious code that is running, and methods for detecting this pa… ▽ More An Air Force evaluation of Multics, and Ken Thompson's famous Turing award lecture "Reflections on Trusting Trust," showed that compilers can be subverted to insert malicious Trojan horses into critical software, including themselves. If this attack goes undetected, even complete analysis of a system's source code will not find the malicious code that is running, and methods for detecting this particular attack are not widely known. This paper describes a practical technique, termed diverse double-compiling (DDC), that detects this attack and some compiler defects as well. Simply recompile the source code twice: once with a second (trusted) compiler, and again using the result of the first compilation. If the result is bit-for-bit identical with the untrusted binary, then the source code accurately represents the binary. This technique has been mentioned informally, but its issues and ramifications have not been identified or discussed in a peer-reviewed work, nor has a public demonstration been made. This paper describes the technique, justifies it, describes how to overcome practical challenges, and demonstrates it. △ Less

Submitted 30 April, 2010; originally announced April 2010.

Comments: 13 pages.

Journal ref: Proc. of the 21st Annual Computer Security Applications Conference (ACSAC), December 5-9, 2005, Tucson, Arizona, pp. 28-40, Los Alamitos: IEEE Computer Society, ISBN 0-7695-2461-3, ISSN 1063-9527, IEEE Computer Society Order Number P2461.

Fully Countering Trusting Trust through Diverse Double-Compiling

Authors: David A. Wheeler

Abstract: An Air Force evaluation of Multics, and Ken Thompson's Turing award lecture ("Reflections on Trusting Trust"), showed that compilers can be subverted to insert malicious Trojan horses into critical software, including themselves. If this "trusting trust" attack goes undetected, even complete analysis of a system's source code will not find the malicious code that is running. Previously-known cou… ▽ More An Air Force evaluation of Multics, and Ken Thompson's Turing award lecture ("Reflections on Trusting Trust"), showed that compilers can be subverted to insert malicious Trojan horses into critical software, including themselves. If this "trusting trust" attack goes undetected, even complete analysis of a system's source code will not find the malicious code that is running. Previously-known countermeasures have been grossly inadequate. If this attack cannot be countered, attackers can quietly subvert entire classes of computer systems, gaining complete control over financial, infrastructure, military, and/or business systems worldwide. This dissertation's thesis is that the trusting trust attack can be detected and effectively countered using the "Diverse Double-Compiling" (DDC) technique, as demonstrated by (1) a formal proof that DDC can determine if source code and generated executable code correspond, (2) a demonstration of DDC with four compilers (a small C compiler, a small Lisp compiler, a small maliciously corrupted Lisp compiler, and a large industrial-strength C compiler, GCC), and (3) a description of approaches for applying DDC in various real-world scenarios. In the DDC technique, source code is compiled twice: the source code of the compiler's parent is compiled using a trusted compiler, and then the putative compiler source code is compiled using the result of the first compilation. If the DDC result is bit-for-bit identical with the original compiler-under-test's executable, and certain other assumptions hold, then the compiler-under-test's executable corresponds with its putative source code. △ Less

Submitted 30 April, 2010; originally announced April 2010.

Comments: PhD dissertation. Accepted by George Mason University, Fairfax, Virginia, USA's Volgenau School of Information Technology and Engineering in 2009. 199 single-side printed pages.

Retro-MACHOs: Pi in the sky?

Authors: Daniel E. Holz, John A. Wheeler

Abstract: Shine a flashlight on a black hole, and one is greeted with the return of a series of concentric rings of light. For a point source of light, and for perfect alignment of the lens, source, and observer, the rings are of infinite brightness (in the limit of geometric optics). In this manner, distant black holes can be revealed through their reflection of light from the Sun. Such retro-MACHO event… ▽ More Shine a flashlight on a black hole, and one is greeted with the return of a series of concentric rings of light. For a point source of light, and for perfect alignment of the lens, source, and observer, the rings are of infinite brightness (in the limit of geometric optics). In this manner, distant black holes can be revealed through their reflection of light from the Sun. Such retro-MACHO events involve photons leaving the Sun, making a Pi rotation about the black hole, and then returning to be detected at the Earth. Our calculations show that, although the light return is quite small, it may nonetheless be detectable for stellar-mass black holes at the edge of our solar system. For example, all (unobscured) black holes of mass M or greater will be observable to a limiting magnitude m, at a distance given by: 0.02 pc x \sqrt[3]{10^{(m-30)/2.5} (M/10 M_sun)^2}. Discovery of a Retro-MACHO offers a way to directly image the presence of a black hole, and would be a stunning confirmation of strong-field general relativity. △ Less

Submitted 2 September, 2002; originally announced September 2002.

Comments: 5 pages, 3 figures. ApJ in press

Journal ref: Astrophys.J. 578 (2002) 330-334

Quantum anti-centrifugal force

Authors: M. A. Cirone, K. Rzazewski, W. P. Schleich, F. Straub, J. A. Wheeler

Abstract: In a two-dimensional world a free quantum particle of vanishing angular momentum experiences an attractive force. This force originates from a modification of the classical centrifugal force due to the wave nature of the particle. For positive energies the quantum anti-centrifugal force manifests itself in a bunching of the nodes of the energy wave functions towards the origin. For negative ener… ▽ More In a two-dimensional world a free quantum particle of vanishing angular momentum experiences an attractive force. This force originates from a modification of the classical centrifugal force due to the wave nature of the particle. For positive energies the quantum anti-centrifugal force manifests itself in a bunching of the nodes of the energy wave functions towards the origin. For negative energies this force is sufficient to create a bound state in a two-dimensional delta function potential. In a counter-intuitive way the attractive force pushes the particle away from the location of the delta function potential. As a consequence, the particle is localized in a band-shaped domain around the origin △ Less

Submitted 16 August, 2001; v1 submitted 14 August, 2001; originally announced August 2001.

Comments: 8 pages, including three eps figures, submitted to Phys. Rev. A. Figures substituted

Journal ref: Phys.Rev.A65:022101,2002

100 Years of the Quantum

Authors: Max Tegmark, John Archibald Wheeler

Abstract: As quantum theory celebrates its 100th birthday, spectacular successes are mixed with outstanding puzzles and promises of new technologies. This article reviews both the successes of quantum theory and the ongoing debate about its consequences for issues ranging from quantum computation to consciousness, parallel universes and the nature of physical reality. We argue that modern experiments and… ▽ More As quantum theory celebrates its 100th birthday, spectacular successes are mixed with outstanding puzzles and promises of new technologies. This article reviews both the successes of quantum theory and the ongoing debate about its consequences for issues ranging from quantum computation to consciousness, parallel universes and the nature of physical reality. We argue that modern experiments and the discovery of decoherence have have shifted prevailing quantum interpretations away from wave function collapse towards unitary physics, and discuss quantum processes in the framework of a tripartite subject-object-environment decomposition. We conclude with some speculations on the bigger picture and the search for a unified theory of quantum gravity. △ Less

Submitted 17 January, 2001; originally announced January 2001.

Comments: This is the original (``director's cut'') version of the Scientific American article, with more text and inferior graphics. Sci. Am. links etc at http://www.hep.upenn.edu/~max/quantum.html

Journal ref: Sci.Am.284:68-75,2001; Spektrum Wiss.Dossier 2003N1:6-14,2003

Apparent horizons in simplicial Brill wave initial data

Authors: Adrian P. Gentle, Daniel E. Holz, Warner A. Miller, John A. Wheeler

Abstract: We construct initial data for a particular class of Brill wave metrics using Regge calculus, and compare the results to a corresponding continuum solution, finding excellent agreement. We then search for trapped surfaces in both sets of initial data, and provide an independent verification of the existence of an apparent horizon once a critical gravitational wave amplitude is passed. Our estimat… ▽ More We construct initial data for a particular class of Brill wave metrics using Regge calculus, and compare the results to a corresponding continuum solution, finding excellent agreement. We then search for trapped surfaces in both sets of initial data, and provide an independent verification of the existence of an apparent horizon once a critical gravitational wave amplitude is passed. Our estimate of this critical value, using both the Regge and continuum solutions, supports other recent findings. △ Less

Submitted 16 December, 1998; originally announced December 1998.

Comments: 7 pages, 6 EPS figures, LaTeX 2e. Submitted to Class. Quant. Grav

Report number: LAUR-98-5276

Journal ref: Class.Quant.Grav. 16 (1999) 1979-1985