Showing 1–50 of 83 results for author: Curtis, A

Towards Practical Finite Sample Bounds for Motion Planning in TAMP

Authors: Seiji Shaw, Aidan Curtis, Leslie Pack Kaelbling, Tomás Lozano-Pérez, Nicholas Roy

Abstract: When using sampling-based motion planners, such as PRMs, in configuration spaces, it is difficult to determine how many samples are required for the PRM to find a solution consistently. This is relevant in Task and Motion Planning (TAMP), where many motion planning problems must be solved in sequence. We attempt to solve this problem by proving an upper bound on the number of samples that are suff… ▽ More When using sampling-based motion planners, such as PRMs, in configuration spaces, it is difficult to determine how many samples are required for the PRM to find a solution consistently. This is relevant in Task and Motion Planning (TAMP), where many motion planning problems must be solved in sequence. We attempt to solve this problem by proving an upper bound on the number of samples that are sufficient, with high probability, to find a solution by drawing on prior work in deterministic sampling and sample complexity theory. We also introduce a numerical algorithm to compute a tighter number of samples based on the proof of the sample complexity theorem we apply to derive our bound. Our experiments show that our numerical bounding algorithm is tight within two orders of magnitude on planar planning problems and becomes looser as the problem's dimensionality increases. When deployed as a heuristic to schedule samples in a TAMP planner, we also observe planning time improvements in planar problems. While our experiments show much work remains to tighten our bounds, the ideas presented in this paper are a step towards a practical sample bound. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Survey and Analysis of IoT Operating Systems: A Comparative Study on the Effectiveness and Acquisition Time of Open Source Digital Forensics Tools

Authors: Jeffrey Fairbanks, Md Mashrur Arifin, Sadia Afreen, Alex Curtis

Abstract: The main goal of this research project is to evaluate the effectiveness and speed of open-source forensic tools for digital evidence collecting from various Internet-of-Things (IoT) devices. The project will create and configure many IoT environments, across popular IoT operating systems, and run common forensics tasks in order to accomplish this goal. To validate these forensic analysis operation… ▽ More The main goal of this research project is to evaluate the effectiveness and speed of open-source forensic tools for digital evidence collecting from various Internet-of-Things (IoT) devices. The project will create and configure many IoT environments, across popular IoT operating systems, and run common forensics tasks in order to accomplish this goal. To validate these forensic analysis operations, a variety of open-source forensic tools covering four standard digital forensics tasks. These tasks will be utilized across each sample IoT operating system and will have its time spent on record carefully tracked down and examined, allowing for a thorough evaluation of the effectiveness and speed for performing forensics on each type of IoT device. The research also aims to offer recommendations to IoT security experts and digital forensic practitioners about the most efficient open-source tools for forensic investigations with IoT devices while maintaining the integrity of gathered evidence and identifying challenges that exist with these new device types. The results will be shared widely and well-documented in order to provide significant contributions to the field of internet-of-things device makers and digital forensics. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Trust the PRoC3S: Solving Long-Horizon Robotics Problems with LLMs and Constraint Satisfaction

Authors: Aidan Curtis, Nishanth Kumar, Jing Cao, Tomás Lozano-Pérez, Leslie Pack Kaelbling

Abstract: Recent developments in pretrained large language models (LLMs) applied to robotics have demonstrated their capacity for sequencing a set of discrete skills to achieve open-ended goals in simple robotic tasks. In this paper, we examine the topic of LLM planning for a set of continuously parameterized skills whose execution must avoid violations of a set of kinematic, geometric, and physical constra… ▽ More Recent developments in pretrained large language models (LLMs) applied to robotics have demonstrated their capacity for sequencing a set of discrete skills to achieve open-ended goals in simple robotic tasks. In this paper, we examine the topic of LLM planning for a set of continuously parameterized skills whose execution must avoid violations of a set of kinematic, geometric, and physical constraints. We prompt the LLM to output code for a function with open parameters, which, together with environmental constraints, can be viewed as a Continuous Constraint Satisfaction Problem (CCSP). This CCSP can be solved through sampling or optimization to find a skill sequence and continuous parameter settings that achieve the goal while avoiding constraint violations. Additionally, we consider cases where the LLM proposes unsatisfiable CCSPs, such as those that are kinematically infeasible, dynamically unstable, or lead to collisions, and re-prompt the LLM to form a new CCSP accordingly. Experiments across three different simulated 3D domains demonstrate that our proposed strategy, PRoC3S, is capable of solving a wide range of complex manipulation tasks with realistic constraints on continuous parameters much more efficiently and effectively than existing baselines. △ Less

Submitted 8 June, 2024; originally announced June 2024.

TAR reconfiguration for vertex set parameters

Authors: Bryan A. Curtis, Mary Flagg, Leslie Hogben

Abstract: This paper surveys results about token addition and removal (TAR) reconfiguration for several well-known vertex set parameters including domination, power domination, standard zero forcing, and PSD zero forcing. We also expand the range of parameters to which universal $X$-set TAR graph results apply, for $X$-sets and their TAR graphs introduced in [B. Bjorkman, C. Bozeman, D. Ferrero, M. Flagg, C… ▽ More This paper surveys results about token addition and removal (TAR) reconfiguration for several well-known vertex set parameters including domination, power domination, standard zero forcing, and PSD zero forcing. We also expand the range of parameters to which universal $X$-set TAR graph results apply, for $X$-sets and their TAR graphs introduced in [B. Bjorkman, C. Bozeman, D. Ferrero, M. Flagg, C. Grood, L. Hogben, B. Jacob, C. Reinhart, Power domination reconfiguration, arXiv:2201.01798] and [N.H. Bong, J. Carlson, B. Curtis, R. Haas, L. Hogben, Isomorphisms and properties of TAR reconfiguration graphs for zero forcing and other $X$-set parameters, \emph{Graphs Combin.} {39} (2023), Paper No. 86]. Here we examine which of the $X$-set axioms are needed for which results. The main results apply to skew zero forcing and vertex covering, and results about TAR reconfiguration graphs of these parameters are presented. While $X$-sets are defined for parameters that take the minimum cardinality over the $X$-sets of a graph, and $X$-set results are restricted to such minimizing parameters, our expansion of the universal perspective allows these results to be applied to parameters that take the maximum value among relevant sets, called $Y$-sets. Maximizing parameters to which the main results apply include independence number, (upper) irredundance number, and (upper) zero forcing irredundance number, and failed zero forcing number; TAR reconfiguration results are presented for these parameters. We also show that the equivalence of connectedness in certain token jumping reconfiguration graphs and certain TAR reconfiguration graphs for independent sets established in [M. Kamiński, P. Medvedev, M. Milani\v c. Complexity of independent set reconfiguration problems. {\em J. Theoretical Computer Science} 439 (2012), 9--15.] extends to $X$-set and $Y$-set parameters. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Variational Prior Replacement in Bayesian Inference and Inversion

Authors: Xuebin Zhao, Andrew Curtis

Abstract: Many scientific investigations require that the values of a set of model parameters are estimated using recorded data. In Bayesian inference, information from both observed data and prior knowledge is combined to update model parameters probabilistically. Prior information represents our belief about the range of values that the variables can take, and their relative probabilities when considered… ▽ More Many scientific investigations require that the values of a set of model parameters are estimated using recorded data. In Bayesian inference, information from both observed data and prior knowledge is combined to update model parameters probabilistically. Prior information represents our belief about the range of values that the variables can take, and their relative probabilities when considered independently of recorded data. Situations arise in which we wish to change prior information: (i) the subjective nature of prior information, (ii) cases in which we wish to test different states of prior information as hypothesis tests, and (iii) information from new studies may emerge so prior information may evolve over time. Estimating the solution to any single inference problem is usually computationally costly, as it typically requires thousands of model samples and their forward simulations. Therefore, recalculating the Bayesian solution every time prior information changes can be extremely expensive. We develop a mathematical formulation that allows prior information to be changed in a solution using variational methods, without performing Bayesian inference on each occasion. In this method, existing prior information is removed from a previously obtained posterior distribution and is replaced by new prior information. We therefore call the methodology variational prior replacement (VPR). We demonstrate VPR using a 2D seismic full waveform inversion example, where VPR provides almost identical posterior solutions compared to those obtained by solving independent inference problems using different priors. The former can be completed within minutes even on a laptop whereas the latter requires days of computations using high-performance computing resources. We demonstrate the value of the method by comparing the posterior solutions obtained using three different types of prior information. △ Less

Submitted 6 June, 2024; originally announced June 2024.

A Single Motor Nano Aerial Vehicle with Novel Peer-to-Peer Communication and Sensing Mechanism

Authors: Jingxian Wang, Andrew G. Curtis, Mark Yim, Michael Rubenstein

Abstract: Communication and position sensing are among the most important capabilities for swarm robots to interact with their peers and perform tasks collaboratively. However, the hardware required to facilitate communication and position sensing is often too complicated, expensive, and bulky to be carried on swarm robots. Here we present Maneuverable Piccolissimo 3 (MP3), a minimalist, single motor drone… ▽ More Communication and position sensing are among the most important capabilities for swarm robots to interact with their peers and perform tasks collaboratively. However, the hardware required to facilitate communication and position sensing is often too complicated, expensive, and bulky to be carried on swarm robots. Here we present Maneuverable Piccolissimo 3 (MP3), a minimalist, single motor drone capable of executing inter-robot communication via infrared light and triangulation-based sensing of relative bearing, distance, and elevation using message arrival time. Thanks to its novel design, MP3 can communicate with peers and localize itself using simple components, keeping its size and mass small and making it inherently safe for human interaction. We present the hardware and software design of MP3 and demonstrate its capability to localize itself, fly stably, and maneuver in the environment using peer-to-peer communication and sensing. △ Less

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

Continuous Sculpting: Persistent Swarm Shape Formation Adaptable to Local Environmental Changes

Authors: Andrew G. Curtis, Mark Yim, Michael Rubenstein

Abstract: Despite their growing popularity, swarms of robots remain limited by the operating time of each individual. We present algorithms which allow a human to sculpt a swarm of robots into a shape that persists in space perpetually, independent of onboard energy constraints such as batteries. Robots generate a path through a shape such that robots cycle in and out of the shape. Robots inside the shape r… ▽ More Despite their growing popularity, swarms of robots remain limited by the operating time of each individual. We present algorithms which allow a human to sculpt a swarm of robots into a shape that persists in space perpetually, independent of onboard energy constraints such as batteries. Robots generate a path through a shape such that robots cycle in and out of the shape. Robots inside the shape react to human initiated changes and adapt the path through the shape accordingly. Robots outside the shape recharge and return to the shape so that the shape can persist indefinitely. The presented algorithms communicate shape changes throughout the swarm using message passing and robot motion. These algorithms enable the swarm to persist through any arbitrary changes to the shape. We describe these algorithms in detail and present their performance in simulation and on a swarm of mobile robots. The result is a swarm behavior more suitable for extended duration, dynamic shape-based tasks in applications such as agriculture and emergency response. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 20 pages, 17 figures

Wehnelt Photoemission in an Ultrafast Electron Microscope: Stability and Usability

Authors: Simon A. Willis, Wyatt A. Curtis, David J. Flannigan

Abstract: We tested and compared the stability and usability of three different cathode materials and configurations in a thermionic-based ultrafast electron microscope: (1) on-axis thermionic and photoemission from a 0.1-mm diameter LaB6 source with graphite guard ring, (2) off-axis photoemission from the Ni aperture surface of the Wehnelt electrode, and (3) on-axis thermionic and photoemission from a 0.2-… ▽ More We tested and compared the stability and usability of three different cathode materials and configurations in a thermionic-based ultrafast electron microscope: (1) on-axis thermionic and photoemission from a 0.1-mm diameter LaB6 source with graphite guard ring, (2) off-axis photoemission from the Ni aperture surface of the Wehnelt electrode, and (3) on-axis thermionic and photoemission from a 0.2-mm diameter polycrystalline Ta source. For each cathode type and configuration, we illustrate how the photoelectron beam-current stability is deleteriously impacted by simultaneous cooling of the source following thermionic heating. Further, we demonstrate usability via collection of parallel- and convergent-beam electron diffraction patterns and by formation of optimum probe size. We find that usability of the off-axis Ni Wehnelt-aperture photoemission is at least comparable to on-axis LaB6 thermionic emission, as well as to on-axis photoemission. However, the stability and achievable beam currents for off-axis photoemission from the Wehnelt aperture were superior to that of the other cathode types and configurations, regardless of the electron-emission mechanism. Beam-current stability for this configuration was found to be within 1% of the mean for 70 minutes (longest duration tested), and steady-state beam current was reached within the sampling-time resolution used here (~1 s) for 15 pA beam currents (i.e., 460 electrons per packet for a 200 kHz repetition rate). Repeatability and robustness of the steady-state condition was also found to be within 1% of the mean. We discuss the implications of these findings for UEM imaging and diffraction experiments, for pulsed-beam damage measurements, and for practical switching between optimum conventional TEM and UEM operation within the same instrument. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 25 pages, 6 figures

$^{171}$Yb$^+$ optical clock with $2.2\times 10^{-18}$ systematic uncertainty and absolute frequency measurements

Authors: Alexandra Tofful, Charles F. A. Baynham, E. Anne Curtis, Adam O. Parsons, Billy I. Robertson, Marco Schioppo, Jacob Tunesi, Helen S. Margolis, Richard J. Hendricks, Josh Whale, Richard C. Thompson, Rachel M. Godun

Abstract: A full evaluation of the uncertainty budget for the ytterbium ion optical clock at the National Physical Laboratory (NPL) was performed on the electric octupole (E3) $^2\mathrm{S}_{1/2}\,\rightarrow\, ^2\mathrm{F}_{7/2}$ transition. The total systematic frequency shift was measured with a fractional standard systematic uncertainty of $2.2\times 10^{-18}$. Furthermore, the absolute frequency of the… ▽ More A full evaluation of the uncertainty budget for the ytterbium ion optical clock at the National Physical Laboratory (NPL) was performed on the electric octupole (E3) $^2\mathrm{S}_{1/2}\,\rightarrow\, ^2\mathrm{F}_{7/2}$ transition. The total systematic frequency shift was measured with a fractional standard systematic uncertainty of $2.2\times 10^{-18}$. Furthermore, the absolute frequency of the E3 transition of the $^{171}$Yb$^+$ ion was measured between 2019 and 2023 via a link to International Atomic Time (TAI) and against the local caesium fountain NPL-CsF2. The absolute frequencies were measured with fractional standard uncertainties between $3.7 \times 10^{-16}$ and $1.1 \times 10^{-15}$, and all were in agreement with the 2021 BIPM recommended frequency. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 16 pages, 8 figures, submitted to Metrologia

Partially Observable Task and Motion Planning with Uncertainty and Risk Awareness

Authors: Aidan Curtis, George Matheos, Nishad Gothoskar, Vikash Mansinghka, Joshua Tenenbaum, Tomás Lozano-Pérez, Leslie Pack Kaelbling

Abstract: Integrated task and motion planning (TAMP) has proven to be a valuable approach to generalizable long-horizon robotic manipulation and navigation problems. However, the typical TAMP problem formulation assumes full observability and deterministic action effects. These assumptions limit the ability of the planner to gather information and make decisions that are risk-aware. We propose a strategy fo… ▽ More Integrated task and motion planning (TAMP) has proven to be a valuable approach to generalizable long-horizon robotic manipulation and navigation problems. However, the typical TAMP problem formulation assumes full observability and deterministic action effects. These assumptions limit the ability of the planner to gather information and make decisions that are risk-aware. We propose a strategy for TAMP with Uncertainty and Risk Awareness (TAMPURA) that is capable of efficiently solving long-horizon planning problems with initial-state and action outcome uncertainty, including problems that require information gathering and avoiding undesirable and irreversible outcomes. Our planner reasons under uncertainty at both the abstract task level and continuous controller level. Given a set of closed-loop goal-conditioned controllers operating in the primitive action space and a description of their preconditions and potential capabilities, we learn a high-level abstraction that can be solved efficiently and then refined to continuous actions for execution. We demonstrate our approach on several robotics problems where uncertainty is a crucial factor and show that reasoning under uncertainty in these problems outperforms previously proposed determinized planning, direct search, and reinforcement learning strategies. Lastly, we demonstrate our planner on two real-world robotics problems using recent advancements in probabilistic perception. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Zero forcing irredundant sets

Authors: Bryan A. Curtis, Leslie Hogben, Adriana Roux

Abstract: Irredundance has been studied in the context of dominating sets, via the concept of private neighbor. Here irredundance of zero forcing sets is introduced via the concept of a private fort and the upper and lower zero forcing irrdedundance numbers $\mbox{ZIR}(G)$ and $\mbox{zir}(G)$ are defined. Bounds on $\mbox{ZIR}(G)$ and $\mbox{zir}(G)$ are established and graphs having extreme values of… ▽ More Irredundance has been studied in the context of dominating sets, via the concept of private neighbor. Here irredundance of zero forcing sets is introduced via the concept of a private fort and the upper and lower zero forcing irrdedundance numbers $\mbox{ZIR}(G)$ and $\mbox{zir}(G)$ are defined. Bounds on $\mbox{ZIR}(G)$ and $\mbox{zir}(G)$ are established and graphs having extreme values of $\mbox{ZIR}(G)$ and $\mbox{zir}(G)$ are characterized. The effect of the join and corona operations is studied. As the concept of a zero forcing irrdedundant set is new, there are many questions for future research. △ Less

Submitted 17 April, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Bayesian Inversion, Uncertainty Analysis and Interrogation using Boosting Variational Inference

Authors: Xuebin Zhao, Andrew Curtis

Abstract: Geoscientists use observed data to estimate properties of the Earth's interior. This often requires non-linear inverse problems to be solved and uncertainties to be estimated. Bayesian inference solves inverse problems under a probabilistic framework, in which uncertainty is represented by a so-called posterior probability distribution. Recently, variational inference has emerged as an efficient m… ▽ More Geoscientists use observed data to estimate properties of the Earth's interior. This often requires non-linear inverse problems to be solved and uncertainties to be estimated. Bayesian inference solves inverse problems under a probabilistic framework, in which uncertainty is represented by a so-called posterior probability distribution. Recently, variational inference has emerged as an efficient method to estimate Bayesian solutions. By seeking the closest approximation to the posterior distribution within any chosen family of distributions, variational inference yields a fully probabilistic solution. It is important to define expressive variational families so that the posterior distribution can be represented accurately. We introduce boosting variational inference (BVI) as a computationally efficient means to construct a flexible approximating family comprising all possible finite mixtures of simpler component distributions. We use Gaussian mixture components due to their fully parametric nature and the ease with which they can be optimised. We apply BVI to seismic travel time tomography and full waveform inversion, comparing its performance with other methods of solution. The results demonstrate that BVI achieves reasonable efficiency and accuracy while enabling the construction of a fully analytic expression for the posterior distribution. Samples that represent major components of uncertainty in the solution can be obtained analytically from each mixture component. We demonstrate that these samples can be used to solve an interrogation problem: to assess the size of a subsurface target structure. To the best of our knowledge, this is the first method in geophysics that provides both analytic and reasonably accurate probabilistic solutions to fully non-linear, high-dimensional Bayesian full waveform inversion problems. △ Less

Submitted 29 December, 2023; originally announced December 2023.

Bayes3D: fast learning and inference in structured generative models of 3D objects and scenes

Authors: Nishad Gothoskar, Matin Ghavami, Eric Li, Aidan Curtis, Michael Noseworthy, Karen Chung, Brian Patton, William T. Freeman, Joshua B. Tenenbaum, Mirko Klukas, Vikash K. Mansinghka

Abstract: Robots cannot yet match humans' ability to rapidly learn the shapes of novel 3D objects and recognize them robustly despite clutter and occlusion. We present Bayes3D, an uncertainty-aware perception system for structured 3D scenes, that reports accurate posterior uncertainty over 3D object shape, pose, and scene composition in the presence of clutter and occlusion. Bayes3D delivers these capabilit… ▽ More Robots cannot yet match humans' ability to rapidly learn the shapes of novel 3D objects and recognize them robustly despite clutter and occlusion. We present Bayes3D, an uncertainty-aware perception system for structured 3D scenes, that reports accurate posterior uncertainty over 3D object shape, pose, and scene composition in the presence of clutter and occlusion. Bayes3D delivers these capabilities via a novel hierarchical Bayesian model for 3D scenes and a GPU-accelerated coarse-to-fine sequential Monte Carlo algorithm. Quantitative experiments show that Bayes3D can learn 3D models of novel objects from just a handful of views, recognizing them more robustly and with orders of magnitude less training data than neural baselines, and tracking 3D objects faster than real time on a single GPU. We also demonstrate that Bayes3D learns complex 3D object models and accurately infers 3D scene composition when used on a Panda robot in a tabletop scenario. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Ookami: An A64FX Computing Resource

Authors: A. C. Calder, E. Siegmann, C. Feldman, S. Chheda, D. C. Smolarski, F. D. Swesty, A. Curtis, J. Dey, D. Carlson, B. Michalowicz, R. J. Harrison

Abstract: We present a look at Ookami, a project providing community access to a testbed supercomputer with the ARM-based A64FX processors developed by a collaboration between RIKEN and Fujitsu and deployed in the Japanese supercomputer Fugaku. We describe the project, provide details about the user base and education/training program, and present highlights from performance studies of two astrophysical sim… ▽ More We present a look at Ookami, a project providing community access to a testbed supercomputer with the ARM-based A64FX processors developed by a collaboration between RIKEN and Fujitsu and deployed in the Japanese supercomputer Fugaku. We describe the project, provide details about the user base and education/training program, and present highlights from performance studies of two astrophysical simulation codes. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: 9 pages, 3 figures, submitted to the Proceedings of 15th International Conference on Numerical Modeling of Space Plasma Flows

The inverse nullity pair problem and the strong nullity interlacing property

Authors: Aida Abiad, Bryan A. Curtis, Mary Flagg, H. Tracy Hall, Jephian C. -H. Lin, Bryan Shader

Abstract: The inverse eigenvalue problem studies the possible spectra among matrices whose off-diagonal entries have their zero-nonzero patterns described by the adjacency of a graph $G$. In this paper, we refer to the $i$-nullity pair of a matrix $A$ as $(\operatorname{null}(A), \operatorname{null}(A(i))$, where $A(i)$ is the matrix obtained from $A$ by removing the $i$-th row and column. The inverse $i$-n… ▽ More The inverse eigenvalue problem studies the possible spectra among matrices whose off-diagonal entries have their zero-nonzero patterns described by the adjacency of a graph $G$. In this paper, we refer to the $i$-nullity pair of a matrix $A$ as $(\operatorname{null}(A), \operatorname{null}(A(i))$, where $A(i)$ is the matrix obtained from $A$ by removing the $i$-th row and column. The inverse $i$-nullity pair problem is considered for complete graphs, cycles, and trees. The strong nullity interlacing property is introduced, and the corresponding supergraph lemma and decontraction lemma are developed as new tools for constructing matrices with a given nullity pair. △ Less

Submitted 20 October, 2023; originally announced October 2023.

MSC Class: 05C50; 05C83; 15A03; 15B57; 65F18

VIP -- Variational Inversion Package with example implementations of Bayesian tomographic imaging

Authors: Xin Zhang, Andrew Curtis

Abstract: Bayesian inference has become an important tool to solve inverse problems and to quantify uncertainties in their solutions. Variational inference is a method that provides probabilistic, Bayesian solutions efficiently by using optimization. In this study we present a Python Variational Inversion Package (VIP), to solve inverse problems using variational inference methods. The package includes auto… ▽ More Bayesian inference has become an important tool to solve inverse problems and to quantify uncertainties in their solutions. Variational inference is a method that provides probabilistic, Bayesian solutions efficiently by using optimization. In this study we present a Python Variational Inversion Package (VIP), to solve inverse problems using variational inference methods. The package includes automatic differential variational inference (ADVI), Stein variational gradient descent (SVGD) and stochastic SVGD (sSVGD), and provides implementations of 2D travel time tomography and 2D full waveform inversion including test examples and solutions. Users can solve their own problems by supplying an appropriate forward function and a gradient calculation code. In addition, the package provides a scalable implementation which can be deployed easily on a desktop machine or using modern high performance computational facilities. The examples demonstrate that VIP is an efficient, scalable, extensible and user-friendly package, and can be used to solve a wide range of low or high dimensional inverse problems in practice. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Apportionable matrices and gracefully labelled graphs

Authors: Antwan Clark, Bryan A. Curtis, Edinah K. Gnang, Leslie Hogben

Abstract: To apportion a complex matrix means to apply a similarity so that all entries of the resulting matrix have the same magnitude. We initiate the study of apportionment, both by unitary matrix similarity and general matrix similarity. There are connections between apportionment and classical graph decomposition problems, including graceful labelings of graphs, Hadamard matrices, and equiangluar lines… ▽ More To apportion a complex matrix means to apply a similarity so that all entries of the resulting matrix have the same magnitude. We initiate the study of apportionment, both by unitary matrix similarity and general matrix similarity. There are connections between apportionment and classical graph decomposition problems, including graceful labelings of graphs, Hadamard matrices, and equiangluar lines, and potential applications to instantaneous uniform mixing in quantum walks. The connection between apportionment and graceful labelings allows the construction of apportionable matrices from trees. A generalization of the well-known Eigenvalue Interlacing Inequalities using graceful labelings is also presented. It is shown that every rank one matrix can be apportioned by a unitary similarity, but there are $2\x 2$ matrices that cannot be apportioned. A necessary condition for a matrix to be apportioned by unitary matrix is established. This condition is used to construct a set of matrices with nonzero Lebesgue measure that are not apportionable by a unitary matrix. △ Less

Submitted 1 June, 2024; v1 submitted 22 August, 2023; originally announced August 2023.

MSC Class: 05C50; 05C78; 15A04; 15A18

Bayesian Variational Time-lapse Full-waveform Inversion

Authors: Xin Zhang, Andrew Curtis

Abstract: Time-lapse seismic full-waveform inversion (FWI) provides estimates of dynamic changes in the subsurface by performing multiple seismic surveys at different times. Since FWI problems are highly non-linear and non-unique, it is important to quantify uncertainties in such estimates to allow robust decision making. Markov chain Monte Carlo (McMC) methods have been used for this purpose, but due to th… ▽ More Time-lapse seismic full-waveform inversion (FWI) provides estimates of dynamic changes in the subsurface by performing multiple seismic surveys at different times. Since FWI problems are highly non-linear and non-unique, it is important to quantify uncertainties in such estimates to allow robust decision making. Markov chain Monte Carlo (McMC) methods have been used for this purpose, but due to their high computational cost, those studies often require an accurate baseline model and estimates of the locations of potential velocity changes, and neglect uncertainty in the baseline velocity model. Such detailed and accurate prior information is not always available in practice. In this study we use an efficient optimization method called stochastic Stein variational gradient descent (sSVGD) to solve time-lapse FWI problems without assuming such prior knowledge, and to estimate uncertainty both in the baseline velocity model and the velocity change. We test two Bayesian strategies: separate Bayesian inversions for each seismic survey, and a single join inversion for baseline and repeat surveys, and compare the methods with the standard linearised double difference inversion. The results demonstrate that all three methods can produce accurate velocity change estimates in the case of having fixed (exactly repeatable) acquisition geometries, but that the two Bayesian methods generate more accurate results when the acquisition geometry changes between surveys. Furthermore the joint inversion provides the most accurate velocity change and uncertainty estimates in all cases. We therefore conclude that Bayesian time-lapse inversion, especially adopting a joint inversion strategy, may be useful to image and monitor the subsurface changes, in particular where uncertainty in the results might lead to significantly different decisions. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Variational Bayesian experimental design for geophysical applications: seismic source location, amplitude versus offset inversion, and estimating CO2 saturations in a subsurface reservoir

Authors: Dominik Strutz, Andrew Curtis

Abstract: This paper introduces variational design methods that are novel to Geophysics, and discusses their benefits and limitations in the context of geophysical applications and more established design methods. Variational methods rely on functional approximations to probability distributions and model-data relationships. They can be used to design experiments that best resolve either all model parameter… ▽ More This paper introduces variational design methods that are novel to Geophysics, and discusses their benefits and limitations in the context of geophysical applications and more established design methods. Variational methods rely on functional approximations to probability distributions and model-data relationships. They can be used to design experiments that best resolve either all model parameters, or the answer to specific questions about the system to be interrogated. The methods are tested in three schematic geophysical applications: (i) estimating a source location given arrival times at sensor locations, and (ii) estimating the contrast in seismic wavefield velocity across a stratal interface given measurements of the amplitudes of seismic wavefield reflections from that interface, and (iii) designing a survey to best constrain CO2 saturation in a subsurface storage scenario. Variational methods allow the value of an experiment to be calculated and optimised simultaneously, which results in substantial savings in computational cost. In the context of designing a survey to best constrain CO2 in a subsurface storage scenario, we show that optimal designs may change substantially depending on the questions of interest. Overall, this work indicates that optimal design methods should be used more widely in Geophysics, as they are in other scientifically advanced fields. △ Less

Submitted 23 January, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: Submitted to Geophysical Journal International in June 2023

Travel times and ray paths for acoustic and elastic waves in generally anisotropic media

Authors: James Ludlam, Katherine Tant, Victorita Dolean, Andrew Curtis

Abstract: Wavefield travel time tomography is used for a variety of purposes in acoustics, geophysics and non-destructive testing. Since the problem is non-linear, assessing uncertainty in the results requires many forward evaluations. It is therefore important that the forward evaluation of travel times and ray paths is efficient, which is challenging in anisotropic media. Given a computed travel time fiel… ▽ More Wavefield travel time tomography is used for a variety of purposes in acoustics, geophysics and non-destructive testing. Since the problem is non-linear, assessing uncertainty in the results requires many forward evaluations. It is therefore important that the forward evaluation of travel times and ray paths is efficient, which is challenging in anisotropic media. Given a computed travel time field, ray tracing can be performed to obtain the fastest ray path from any point in the medium to the source of the travel time field. These rays can then be used to speed up gradient based inversion methods. We present a forward modeller for calculating travel time fields by localised estimation of wavefronts, and a novel approach to ray tracing through travel time fields. These methods have been tested in a complex anisotropic weld and give travel times comparable to those obtained using finite element modelling while being computationally cheaper. △ Less

Submitted 18 August, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

MSC Class: 65N06 (Primary) 86-10; 00A69; 35L05 (Secondary)

Analysis of atomic-clock data to constrain variations of fundamental constants

Authors: Nathaniel Sherrill, Adam O. Parsons, Charles F. A. Baynham, William Bowden, E. Anne Curtis, Richard Hendricks, Ian R. Hill, Richard Hobson, Helen S. Margolis, Billy I. Robertson, Marco Schioppo, Krzysztof Szymaniec, Alexandra Tofful, Jacob Tunesi, Rachel M. Godun, Xavier Calmet

Abstract: We present a new framework to study the time variation of fundamental constants in a model-independent way. Model independence implies more free parameters than assumed in previous studies. Using data from atomic clocks based on $^{87}$Sr, $^{171}$Yb$^+$ and $^{133}$Cs, we set bounds on parameters controlling the variation of the fine-structure constant, $α$, and the electron-to-proton mass ratio,… ▽ More We present a new framework to study the time variation of fundamental constants in a model-independent way. Model independence implies more free parameters than assumed in previous studies. Using data from atomic clocks based on $^{87}$Sr, $^{171}$Yb$^+$ and $^{133}$Cs, we set bounds on parameters controlling the variation of the fine-structure constant, $α$, and the electron-to-proton mass ratio, $μ$. We consider variations on timescales ranging from a minute to almost a day. In addition, we use our results to derive some of the tightest limits to date on the parameter space of models of ultralight dark matter and axion-like particles. △ Less

Submitted 15 December, 2023; v1 submitted 9 February, 2023; originally announced February 2023.

Journal ref: New J.Phys. 25 (2023) 9, 093012

Orthogonal realizations of random sign patterns and other applications of the SIPP

Authors: Zachary Brennan, Christopher Cox, Bryan A. Curtis, Enrique Gomez-Leos, Kimberly P. Hadaway, Leslie Hogben, Conor Thompson

Abstract: A sign pattern is an array with entries in $\{+,-,0\}$. A matrix $Q$ is row orthogonal if $QQ^T = I$. The Strong Inner Product Property (SIPP), introduced in [B.A.~Curtis and B.L.~Shader, Sign patterns of orthogonal matrices and the strong inner product property, Linear Algebra Appl. 592: 228--259, 2020], is an important tool when determining whether a sign pattern allows row orthogonality because… ▽ More A sign pattern is an array with entries in $\{+,-,0\}$. A matrix $Q$ is row orthogonal if $QQ^T = I$. The Strong Inner Product Property (SIPP), introduced in [B.A.~Curtis and B.L.~Shader, Sign patterns of orthogonal matrices and the strong inner product property, Linear Algebra Appl. 592: 228--259, 2020], is an important tool when determining whether a sign pattern allows row orthogonality because it guarantees there is a nearby matrix with the same property, allowing zero entries to be perturbed to nonzero entries, while preserving the sign of every nonzero entry. This paper uses the SIPP to initiate the study of conditions under which random sign patterns allow row orthogonality with high probability. Building on prior work, $5\times n$ nowhere zero sign patterns that minimally allow orthogonality are determined. Conditions on zero entries in a sign pattern are established that guarantee any row orthogonal matrix with such a sign pattern has the SIPP. △ Less

Submitted 10 December, 2022; originally announced December 2022.

MSC Class: 15B10; 15B35; 15B52; 60B20

Task-Directed Exploration in Continuous POMDPs for Robotic Manipulation of Articulated Objects

Authors: Aidan Curtis, Leslie Kaelbling, Siddarth Jain

Abstract: Representing and reasoning about uncertainty is crucial for autonomous agents acting in partially observable environments with noisy sensors. Partially observable Markov decision processes (POMDPs) serve as a general framework for representing problems in which uncertainty is an important factor. Online sample-based POMDP methods have emerged as efficient approaches to solving large POMDPs and hav… ▽ More Representing and reasoning about uncertainty is crucial for autonomous agents acting in partially observable environments with noisy sensors. Partially observable Markov decision processes (POMDPs) serve as a general framework for representing problems in which uncertainty is an important factor. Online sample-based POMDP methods have emerged as efficient approaches to solving large POMDPs and have been shown to extend to continuous domains. However, these solutions struggle to find long-horizon plans in problems with significant uncertainty. Exploration heuristics can help guide planning, but many real-world settings contain significant task-irrelevant uncertainty that might distract from the task objective. In this paper, we propose STRUG, an online POMDP solver capable of handling domains that require long-horizon planning with significant task-relevant and task-irrelevant uncertainty. We demonstrate our solution on several temporally extended versions of toy POMDP problems as well as robotic manipulation of articulated objects using a neural perception frontend to construct a distribution of possible models. Our results show that STRUG outperforms the current sample-based online POMDP solvers on several tasks. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Visibility-Aware Navigation Among Movable Obstacles

Authors: Jose Muguira-Iturralde, Aidan Curtis, Yilun Du, Leslie Pack Kaelbling, Tomás Lozano-Pérez

Abstract: In this paper, we examine the problem of visibility-aware robot navigation among movable obstacles (VANAMO). A variant of the well-known NAMO robotic planning problem, VANAMO puts additional visibility constraints on robot motion and object movability. This new problem formulation lifts the restrictive assumption that the map is fully visible and the object positions are fully known. We provide a… ▽ More In this paper, we examine the problem of visibility-aware robot navigation among movable obstacles (VANAMO). A variant of the well-known NAMO robotic planning problem, VANAMO puts additional visibility constraints on robot motion and object movability. This new problem formulation lifts the restrictive assumption that the map is fully visible and the object positions are fully known. We provide a formal definition of the VANAMO problem and propose the Look and Manipulate Backchaining (LaMB) algorithm for solving such problems. LaMB has a simple vision-based API that makes it more easily transferable to real-world robot applications and scales to the large 3D environments. To evaluate LaMB, we construct a set of tasks that illustrate the complex interplay between visibility and object movability that can arise in mobile base manipulation problems in unknown environments. We show that LaMB outperforms NAMO and visibility-aware motion planning approaches as well as simple combinations of them on complex manipulation problems with partial observability. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Introducing Conceptual Geological Information into Bayesian Tomographic Imaging

Authors: Hugo Bloem, Andrew Curtis, Daniel Tetzlaff

Abstract: Geological process models simulate a range of dynamic processes to evolve a base topography into a final 2D cross-section or 3D geological scenario. In principle, process parameters may be updated to better align with observed geophysical or geological data; however, many realisations of process models that embody different conceptual models may provide similar consistency with observed data, and… ▽ More Geological process models simulate a range of dynamic processes to evolve a base topography into a final 2D cross-section or 3D geological scenario. In principle, process parameters may be updated to better align with observed geophysical or geological data; however, many realisations of process models that embody different conceptual models may provide similar consistency with observed data, and finding all such realisations may be infeasible due to the computational demands of the task. Alternatively, geophysical probabilistic tomographic methods may be used to estimate the family of models of a target subsurface structure that are consistent both with information obtained from previous experiments and with new data (the Bayesian posterior distribution). However, this family seldom embodies geologically reasonable images. We show that the posterior distribution of tomographic images can be enhanced by injecting geological prior information into Bayesian inference, and we can do this near-instantaneously by trained MDNs. We invoke two geological concepts as geological prior information: a braided river system, and a set of marine parasequences, each parameterised by a GAN. Data from a target structure can then be used to infer the image parameter values using either geological concept using MDNs. Our MDN solutions closely resemble those found using expensive McMC methods, and while the use of incorrect geological conceptual models provides less accurate results the mean structures still approximate the target. We then show that a classification network can infer the correct geological conceptual model. Thus, imposing even incorrect geological prior information may improve geophysical tomographic images compared to those obtained without geological prior information, and in principle geological conceptual models can be inferred directly from geophysical travel time data. △ Less

Submitted 14 October, 2022; originally announced October 2022.

Journal ref: BRE-092-2023

3D Bayesian Variational Full Waveform Inversion

Authors: Xin Zhang, Angus Lomas, Muhong Zhou, York Zheng, Andrew Curtis

Abstract: Seismic full-waveform inversion (FWI) provides high resolution images of the subsurface by exploiting information in the recorded seismic waveforms. This is achieved by solving a highly nonnlinear and nonunique inverse problem. Bayesian inference is therefore used to quantify uncertainties in the solution. Variational inference is a method that provides probabilistic, Bayesian solutions efficientl… ▽ More Seismic full-waveform inversion (FWI) provides high resolution images of the subsurface by exploiting information in the recorded seismic waveforms. This is achieved by solving a highly nonnlinear and nonunique inverse problem. Bayesian inference is therefore used to quantify uncertainties in the solution. Variational inference is a method that provides probabilistic, Bayesian solutions efficiently using optimization. The method has been applied to 2D FWI problems to produce full Bayesian posterior distributions. However, due to higher dimensionality and more expensive computational cost, the performance of the method in 3D FWI problems remains unknown. We apply three variational inference methods to 3D FWI and analyse their performance. Specifically we apply automatic differential variational inference (ADVI), Stein variational gradient descent (SVGD) and stochastic SVGD (sSVGD), to a 3D FWI problem, and compare their results and computational cost. The results show that ADVI is the most computationally efficient method but systematically underestimates the uncertainty. The method can therefore be used to provide relatively rapid but approximate insights into the subsurface together with a lower bound estimate of the uncertainty. SVGD demands the highest computational cost, and still produces biased results. In contrast, by including a randomized term in the SVGD dynamics, sSVGD becomes a Markov chain Monte Carlo method and provides the most accurate results at intermediate computational cost. We thus conclude that 3D variational full-waveform inversion is practically applicable, at least in small problems, and can be used to image the Earth's interior and to provide reasonable uncertainty estimates on those images. △ Less

Submitted 10 October, 2022; v1 submitted 7 October, 2022; originally announced October 2022.

Modular interface for managing cognitive bias in experts

Authors: Melody G Whitehead, Andrew Curtis

Abstract: Expert knowledge is required to interpret data across a range of fields. Experts bridge gaps that often exists in our knowledge about relationships between data and the parameters of interest. This is especially true in geoscientific applications, where knowledge of the Earth is derived from interpretations of observable features and relies on predominantly unproven but widely accepted theories. T… ▽ More Expert knowledge is required to interpret data across a range of fields. Experts bridge gaps that often exists in our knowledge about relationships between data and the parameters of interest. This is especially true in geoscientific applications, where knowledge of the Earth is derived from interpretations of observable features and relies on predominantly unproven but widely accepted theories. Thus, experts facilitate solutions to otherwise unsolvable problems. However, experts are inherently subjective, and susceptible to cognitive biases and adverse external effects. This work examines this problem within geoscience. Three compelling examples are provided of the prevalence of cognitive biases from previous work. The problem is then formally defined, and a set of design principles which ensure that any solution is sufficiently flexible to be readily applied to the range of geoscientific problems. No solutions exist that reliably capture and reduce cognitive bias in experts. However, formal expert elicitation methods can be used to assess expert variation, and a variety of approaches exist that may help to illuminate uncertainties, avoid misunderstandings, and reduce herding behaviours or single-expert over-dominance. This work combines existing and future approaches to reduce expert suboptimality through a flexible modular design where each module provides a specific function. The design centres around action modules that force a stop-and-perform step into interpretation tasks. A starter-pack of modules is provided as an example of the conceptual design. This simple bias-reduction system may readily be applied in organisations and during everyday interpretations through to tasks for major commercial ventures. △ Less

Submitted 8 August, 2022; originally announced August 2022.

On Using Linux Kernel Huge Pages with FLASH, an Astrophysical Simulation Code

Authors: Alan C. Calder, Catherine Feldman, Eva Siegmann, John Dey, Anthony Curtis, Smeet Chheda, Robert J. Harrison

Abstract: We present efforts at improving the performance of FLASH, a multi-scale, multi-physics simulation code principally for astrophysical applications, by using huge pages on Ookami, an HPE Apollo 80 A64FX platform. FLASH is written principally in modern Fortran and makes use of the PARAMESH library to manage a block-structured adaptive mesh. We explored options for enabling the use of huge pages with… ▽ More We present efforts at improving the performance of FLASH, a multi-scale, multi-physics simulation code principally for astrophysical applications, by using huge pages on Ookami, an HPE Apollo 80 A64FX platform. FLASH is written principally in modern Fortran and makes use of the PARAMESH library to manage a block-structured adaptive mesh. We explored options for enabling the use of huge pages with several compilers, but we were only able to successfully use huge pages when compiling with the Fujitsu compiler. The use of huge pages substantially reduced the number of translation lookaside buffer misses, but overall performance gains were marginal. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 6 pages, 1 figure, accepted to Embracing Arm for HPC, An IEEE Cluster 2022 Workshop

An introduction to variational inference in Geophysical inverse problems

Authors: Xin Zhang, Muhammad Atif Nawaz, Xuebin Zhao, Andrew Curtis

Abstract: In a variety of scientific applications we wish to characterize a physical system using measurements or observations. This often requires us to solve an inverse problem, which usually has non-unique solutions so uncertainty must be quantified in order to define the family of all possible solutions. Bayesian inference provides a powerful theoretical framework which defines the set of solutions to i… ▽ More In a variety of scientific applications we wish to characterize a physical system using measurements or observations. This often requires us to solve an inverse problem, which usually has non-unique solutions so uncertainty must be quantified in order to define the family of all possible solutions. Bayesian inference provides a powerful theoretical framework which defines the set of solutions to inverse problems, and variational inference is a method to solve Bayesian inference problems using optimization while still producing fully probabilistic solutions. This chapter provides an introduction to variational inference, and reviews its applications to a range of geophysical problems, including petrophysical inversion, travel time tomography and full-waveform inversion. We demonstrate that variational inference is an efficient and scalable method which can be deployed in many practical scenarios. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Journal ref: Advances in Geophysics, Vol. 62, pp. 73-140, 2021. Elsevier

PG3: Policy-Guided Planning for Generalized Policy Generation

Authors: Ryan Yang, Tom Silver, Aidan Curtis, Tomas Lozano-Perez, Leslie Pack Kaelbling

Abstract: A longstanding objective in classical planning is to synthesize policies that generalize across multiple problems from the same domain. In this work, we study generalized policy search-based methods with a focus on the score function used to guide the search over policies. We demonstrate limitations of two score functions and propose a new approach that overcomes these limitations. The main idea b… ▽ More A longstanding objective in classical planning is to synthesize policies that generalize across multiple problems from the same domain. In this work, we study generalized policy search-based methods with a focus on the score function used to guide the search over policies. We demonstrate limitations of two score functions and propose a new approach that overcomes these limitations. The main idea behind our approach, Policy-Guided Planning for Generalized Policy Generation (PG3), is that a candidate policy should be used to guide planning on training problems as a mechanism for evaluating that candidate. Theoretical results in a simplified setting give conditions under which PG3 is optimal or admissible. We then study a specific instantiation of policy search where planning problems are PDDL-based and policies are lifted decision lists. Empirical results in six domains confirm that PG3 learns generalized policies more efficiently and effectively than several baselines. Code: https://github.com/ryangpeixu/pg3 △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: IJCAI 2022

Surface wave dispersion inversion using an energy likelihood function

Authors: Xin Zhang, York Zheng, Andrew Curtis

Abstract: Seismic surface wave dispersion inversion is used widely to study the subsurface structure of the Earth. The dispersion property is usually measured by using frequency-phase velocity (f-c) analysis and by picking phase velocities from the obtained f-c spectrum. However, because of potential contamination the f-c spectrum often has multimodalities at each frequency for each mode. These introduce un… ▽ More Seismic surface wave dispersion inversion is used widely to study the subsurface structure of the Earth. The dispersion property is usually measured by using frequency-phase velocity (f-c) analysis and by picking phase velocities from the obtained f-c spectrum. However, because of potential contamination the f-c spectrum often has multimodalities at each frequency for each mode. These introduce uncertainty and errors in the picked phase velocities, and consequently the obtained shear velocity structure is biased. To overcome this issue, in this study we introduce a new method which directly uses the spectrum as data. We achieve this by solving the inverse problem in a Bayesian framework and define a new likelihood function, the energy likelihood function, which uses the spectrum energy to define data fit. We apply the new method to a land dataset recorded by a dense receiver array, and compare the results to those obtained using the traditional method. The results show that the new method produces more accurate results since they better match independent data from refraction tomography. This real-data application also shows that it can be applied efficiently since it removes the need to pick phase velocities, and with relatively little adjustment to current practice since it uses standard f-c panels to define the likelihood. We therefore recommend using the energy likelihood function rather than explicitly picking phase velocities in surface wave dispersion inversion. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Let's Handle It: Generalizable Manipulation of Articulated Objects

Authors: Zhutian Yang, Aidan Curtis

Abstract: In this project we present a framework for building generalizable manipulation controller policies that map from raw input point clouds and segmentation masks to joint velocities. We took a traditional robotics approach, using point cloud processing, end-effector trajectory calculation, inverse kinematics, closed-loop position controllers, and behavior trees. We demonstrate our framework on four m… ▽ More In this project we present a framework for building generalizable manipulation controller policies that map from raw input point clouds and segmentation masks to joint velocities. We took a traditional robotics approach, using point cloud processing, end-effector trajectory calculation, inverse kinematics, closed-loop position controllers, and behavior trees. We demonstrate our framework on four manipulation skills on common household objects that comprise the SAPIEN ManiSkill Manipulation challenge. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Measuring the stability of fundamental constants with a network of clocks

Authors: G. Barontini, L. Blackburn, V. Boyer, F. Butuc-Mayer, X. Calmet, J. R. Crespo Lopez-Urrutia, E. A. Curtis, B. Darquie, J. Dunningham, N. J. Fitch, E. M. Forgan, K. Georgiou, P. Gill, R. M. Godun, J. Goldwin, V. Guarrera, A. C. Harwood, I. R. Hill, R. J. Hendricks, M. Jeong, M. Y. H. Johnson, M. Keller, L. P. Kozhiparambil Sajith, F. Kuipers, H. S. Margolis , et al. (19 additional authors not shown)

Abstract: The detection of variations of fundamental constants of the Standard Model would provide us with compelling evidence of new physics, and could lift the veil on the nature of dark matter and dark energy. In this work, we discuss how a network of atomic and molecular clocks can be used to look for such variations with unprecedented sensitivity over a wide range of time scales. This is precisely the… ▽ More The detection of variations of fundamental constants of the Standard Model would provide us with compelling evidence of new physics, and could lift the veil on the nature of dark matter and dark energy. In this work, we discuss how a network of atomic and molecular clocks can be used to look for such variations with unprecedented sensitivity over a wide range of time scales. This is precisely the goal of the recently launched QSNET project: A network of clocks for measuring the stability of fundamental constants. QSNET will include state-of-the-art atomic clocks, but will also develop next-generation molecular and highly charged ion clocks with enhanced sensitivity to variations of fundamental constants. We describe the technological and scientific aims of QSNET and evaluate its expected performance. We show that in the range of parameters probed by QSNET, either we will discover new physics, or we will impose new constraints on violations of fundamental symmetries and a range of theories beyond the Standard Model, including dark matter and dark energy models. △ Less

Submitted 11 May, 2022; v1 submitted 20 December, 2021; originally announced December 2021.

Journal ref: EPJ Quantum Technology volume 9, Article number: 12 (2022)

Map Induction: Compositional spatial submap learning for efficient exploration in novel environments

Authors: Sugandha Sharma, Aidan Curtis, Marta Kryven, Josh Tenenbaum, Ila Fiete

Abstract: Humans are expert explorers. Understanding the computational cognitive mechanisms that support this efficiency can advance the study of the human mind and enable more efficient exploration algorithms. We hypothesize that humans explore new environments efficiently by inferring the structure of unobserved spaces using spatial information collected from previously explored spaces. This cognitive pro… ▽ More Humans are expert explorers. Understanding the computational cognitive mechanisms that support this efficiency can advance the study of the human mind and enable more efficient exploration algorithms. We hypothesize that humans explore new environments efficiently by inferring the structure of unobserved spaces using spatial information collected from previously explored spaces. This cognitive process can be modeled computationally using program induction in a Hierarchical Bayesian framework that explicitly reasons about uncertainty with strong spatial priors. Using a new behavioral Map Induction Task, we demonstrate that this computational framework explains human exploration behavior better than non-inductive models and outperforms state-of-the-art planning algorithms when applied to a realistic spatial navigation domain. △ Less

Submitted 17 March, 2022; v1 submitted 23 October, 2021; originally announced October 2021.

Discovering State and Action Abstractions for Generalized Task and Motion Planning

Authors: Aidan Curtis, Tom Silver, Joshua B. Tenenbaum, Tomas Lozano-Perez, Leslie Pack Kaelbling

Abstract: Generalized planning accelerates classical planning by finding an algorithm-like policy that solves multiple instances of a task. A generalized plan can be learned from a few training examples and applied to an entire domain of problems. Generalized planning approaches perform well in discrete AI planning problems that involve large numbers of objects and extended action sequences to achieve the g… ▽ More Generalized planning accelerates classical planning by finding an algorithm-like policy that solves multiple instances of a task. A generalized plan can be learned from a few training examples and applied to an entire domain of problems. Generalized planning approaches perform well in discrete AI planning problems that involve large numbers of objects and extended action sequences to achieve the goal. In this paper, we propose an algorithm for learning features, abstractions, and generalized plans for continuous robotic task and motion planning (TAMP) and examine the unique difficulties that arise when forced to consider geometric and physical constraints as a part of the generalized plan. Additionally, we show that these simple generalized plans learned from only a handful of examples can be used to improve the search efficiency of TAMP solvers. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Long-Horizon Manipulation of Unknown Objects via Task and Motion Planning with Estimated Affordances

Authors: Aidan Curtis, Xiaolin Fang, Leslie Pack Kaelbling, Tomás Lozano-Pérez, Caelan Reed Garrett

Abstract: We present a strategy for designing and building very general robot manipulation systems involving the integration of a general-purpose task-and-motion planner with engineered and learned perception modules that estimate properties and affordances of unknown objects. Such systems are closed-loop policies that map from RGB images, depth images, and robot joint encoder measurements to robot joint po… ▽ More We present a strategy for designing and building very general robot manipulation systems involving the integration of a general-purpose task-and-motion planner with engineered and learned perception modules that estimate properties and affordances of unknown objects. Such systems are closed-loop policies that map from RGB images, depth images, and robot joint encoder measurements to robot joint position commands. We show that following this strategy a task-and-motion planner can be used to plan intelligent behaviors even in the absence of a priori knowledge regarding the set of manipulable objects, their geometries, and their affordances. We explore several different ways of implementing such perceptual modules for segmentation, property detection, shape estimation, and grasp generation. We show how these modules are integrated within the PDDLStream task and motion planning framework. Finally, we demonstrate that this strategy can enable a single system to perform a wide variety of real-world multi-step manipulation tasks, generalizing over a broad class of objects, object arrangements, and goals, without any prior knowledge of the environment and without re-training. △ Less

Submitted 10 August, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: The first two authors contributed equally and are listed in alphabetical order

Principles for Evaluation of AI/ML Model Performance and Robustness

Authors: Olivia Brown, Andrew Curtis, Justin Goodwin

Abstract: The Department of Defense (DoD) has significantly increased its investment in the design, evaluation, and deployment of Artificial Intelligence and Machine Learning (AI/ML) capabilities to address national security needs. While there are numerous AI/ML successes in the academic and commercial sectors, many of these systems have also been shown to be brittle and nonrobust. In a complex and ever-cha… ▽ More The Department of Defense (DoD) has significantly increased its investment in the design, evaluation, and deployment of Artificial Intelligence and Machine Learning (AI/ML) capabilities to address national security needs. While there are numerous AI/ML successes in the academic and commercial sectors, many of these systems have also been shown to be brittle and nonrobust. In a complex and ever-changing national security environment, it is vital that the DoD establish a sound and methodical process to evaluate the performance and robustness of AI/ML models before these new capabilities are deployed to the field. This paper reviews the AI/ML development process, highlights common best practices for AI/ML model evaluation, and makes recommendations to DoD evaluators to ensure the deployment of robust AI/ML capabilities for national security needs. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Interrogating probabilistic inversion results for subsurface structural information

Authors: Xin Zhang, Andrew Curtis

Abstract: The goal of a scientific investigation is to find answers to specific questions. In geosciences this is typically achieved by solving an inference or inverse problem and interpreting the solution. However, the answer obtained is often biased because the solution to an inverse problem is nonunique and human interpretation is a biased process. Interrogation theory provides a systematic way to find o… ▽ More The goal of a scientific investigation is to find answers to specific questions. In geosciences this is typically achieved by solving an inference or inverse problem and interpreting the solution. However, the answer obtained is often biased because the solution to an inverse problem is nonunique and human interpretation is a biased process. Interrogation theory provides a systematic way to find optimal answers by considering their full uncertainty estimates, and by designing an objective function that defines desirable qualities in the answer. In this study we demonstrate interrogation theory by quantifying the size of a particular subsurface structure. The results show that interrogation theory provides an accurate estimate of the true answer, which cannot be obtained by direct, subjective interpretation of the solution mean and standard deviation. This demonstrates the value of interrogation theory. It also shows that fully nonlinear uncertainty assessments may be critical in order to address real-world scientific problems, which goes some way towards justifying their computational expense. △ Less

Submitted 6 July, 2021; originally announced July 2021.

Real-time super-resolution mapping of locally anisotropic grain orientations for ultrasonic non-destructive evaluation of crystalline material

Authors: Jonathan Singh, Katherine M. M. Tant, Andrew Curtis, Anthony J. Mulholland

Abstract: Estimating the spatially varying microstructures of heterogeneous and locally anisotropic media non-destructively is necessary for the accurate detection of flaws and reliable monitoring of manufacturing processes. Conventional algorithms used for solving this inverse problem come with significant computational cost, particularly in the case of high dimensional non-linear tomographic problems. In… ▽ More Estimating the spatially varying microstructures of heterogeneous and locally anisotropic media non-destructively is necessary for the accurate detection of flaws and reliable monitoring of manufacturing processes. Conventional algorithms used for solving this inverse problem come with significant computational cost, particularly in the case of high dimensional non-linear tomographic problems. In this paper, we propose a framework which uses deep neural networks (DNNs) with full aperture, pitch-catch and pulse-echo transducer configurations to reconstruct material maps of crystallographic orientation. We also present the first ever application of generative adversarial networks (GANs) to achieve super resolution of ultrasonic tomographic images, providing a factor-four increase in image resolution and up to a 50% increase in structural similarity. The importance of including appropriate prior knowledge in the GAN training dataset to increase inversion accuracy is highlighted; known information about the material's structure should be present in the training data. We show that after a computationally expensive training process, the DNNs and GANs can be used in less that one second (0.9 seconds on a standard desktop computer) to provide a high resolution map of the material's grain orientations. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Broadband acoustic invisibility and illusions

Authors: Theodor S. Becker, Dirk-Jan van Manen, Thomas Haag, Christoph Bärlocher, Xun Li, Nele Börsing, Andrew Curtis, Marc Serra-Garcia, Johan O. A. Robertsson

Abstract: Rendering objects invisible to impinging acoustic waves (cloaking) and creating acoustic illusions (holography) has been attempted using active and passive approaches. While passive methods are applicable only to narrow frequency bands, active approaches attempt to respond dynamically, interfering with broadband incident or scattered wavefields by emitting secondary waves. Without prior knowledge… ▽ More Rendering objects invisible to impinging acoustic waves (cloaking) and creating acoustic illusions (holography) has been attempted using active and passive approaches. While passive methods are applicable only to narrow frequency bands, active approaches attempt to respond dynamically, interfering with broadband incident or scattered wavefields by emitting secondary waves. Without prior knowledge of the primary wavefield, the signals for the secondary sources need to be estimated and adapted in real-time. This has thus far impeded active cloaking and holography for broadband wavefields. We present experimental results of active acoustic cloaking and holography without prior knowledge of the wavefield so that objects remain invisible and illusions intact even for broadband moving sources. This opens novel research directions and facilitates practical applications including architectural acoustics, and stealth. △ Less

Submitted 21 May, 2021; v1 submitted 17 May, 2021; originally announced May 2021.

Seismic interferometry from correlated noise sources

Authors: Daniella Ayala-Garcia, Andrew Curtis, Michal Branicki

Abstract: It is a well-established principle that cross-correlating seismic observations at different receiver locations can yield estimates of band-limited inter-receiver Green's functions. This principle, known as seismic interferometry, is a powerful technique that can transform noise into signals which allow us to remotely image and interrogate subsurface Earth structures. In practice it is often necess… ▽ More It is a well-established principle that cross-correlating seismic observations at different receiver locations can yield estimates of band-limited inter-receiver Green's functions. This principle, known as seismic interferometry, is a powerful technique that can transform noise into signals which allow us to remotely image and interrogate subsurface Earth structures. In practice it is often necessary and even desirable to rely on noise already present in the environment. Theory that underpins many applications of ambient noise interferometry makes an assumption that the noise sources are uncorrelated in space and time. However, many real-world noise sources such as trains, highway traffic and ocean waves are inherently correlated both in space and time, in direct contradiction to the current theoretical foundations. Applying standard interferometric techniques to recordings from correlated energy sources makes the Green's function liable to estimation errors that so far have not been fully accounted for theoretically nor in practice. We show that these errors are significant for common noise sources, always perturbing and sometimes obscuring the phase one wishes to retrieve. Our analysis explains why stacking may reduce the phase errors, but also shows that in commonly-encountered circumstances stacking will not remediate the problem. This analytical insight allowed us to develop a novel workflow that significantly mitigates effects arising from the use of correlated noise sources. Our methodology can be used in conjunction with already existing approaches, and improves results from both correlated and uncorrelated ambient noise. Hence, we expect it to be widely applicable in real life ambient noise studies. △ Less

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

Bayesian Full-waveform Inversion with Realistic Priors

Authors: Xin Zhang, Andrew Curtis

Abstract: Seismic full-waveform inversion (FWI) uses full seismic records to estimate subsurface velocity structure. This requires a highly nonlinear and nonunique inverse problem to be solved, and Bayesian methods have been used to quantify uncertainties in the solution. Variational Bayesian inference uses optimization to provide solutions efficiently. The method has been applied to solve a transmission FW… ▽ More Seismic full-waveform inversion (FWI) uses full seismic records to estimate subsurface velocity structure. This requires a highly nonlinear and nonunique inverse problem to be solved, and Bayesian methods have been used to quantify uncertainties in the solution. Variational Bayesian inference uses optimization to provide solutions efficiently. The method has been applied to solve a transmission FWI problem using data generated by known earthquake-like sources, with strong prior information imposed on the velocity. Unfortunately such prior information about velocity structure and earthquake sources is never available in practice. We present the first application of the method in a seismic reflection setting, and with realistically weak prior information. We thus demonstrate that the method can produce high resolution images and reliable uncertainties given practically reasonable prior information. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: 7 pages, 3 figures

Planning with Learned Object Importance in Large Problem Instances using Graph Neural Networks

Authors: Tom Silver, Rohan Chitnis, Aidan Curtis, Joshua Tenenbaum, Tomas Lozano-Perez, Leslie Pack Kaelbling

Abstract: Real-world planning problems often involve hundreds or even thousands of objects, straining the limits of modern planners. In this work, we address this challenge by learning to predict a small set of objects that, taken together, would be sufficient for finding a plan. We propose a graph neural network architecture for predicting object importance in a single inference pass, thus incurring little… ▽ More Real-world planning problems often involve hundreds or even thousands of objects, straining the limits of modern planners. In this work, we address this challenge by learning to predict a small set of objects that, taken together, would be sufficient for finding a plan. We propose a graph neural network architecture for predicting object importance in a single inference pass, thus incurring little overhead while greatly reducing the number of objects that must be considered by the planner. Our approach treats the planner and transition model as black boxes, and can be used with any off-the-shelf planner. Empirically, across classical planning, probabilistic planning, and robotic task and motion planning, we find that our method results in planning that is significantly faster than several baselines, including other partial grounding strategies and lifted planners. We conclude that learning to predict a sufficient set of objects for a planning problem is a simple, powerful, and general mechanism for planning in large instances. Video: https://youtu.be/FWsVJc2fvCE Code: https://git.io/JIsqX △ Less

Submitted 8 December, 2020; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: AAAI 2021

ThreeDWorld: A Platform for Interactive Multi-Modal Physical Simulation

Authors: Chuang Gan, Jeremy Schwartz, Seth Alter, Damian Mrowca, Martin Schrimpf, James Traer, Julian De Freitas, Jonas Kubilius, Abhishek Bhandwaldar, Nick Haber, Megumi Sano, Kuno Kim, Elias Wang, Michael Lingelbach, Aidan Curtis, Kevin Feigelis, Daniel M. Bear, Dan Gutfreund, David Cox, Antonio Torralba, James J. DiCarlo, Joshua B. Tenenbaum, Josh H. McDermott, Daniel L. K. Yamins

Abstract: We introduce ThreeDWorld (TDW), a platform for interactive multi-modal physical simulation. TDW enables simulation of high-fidelity sensory data and physical interactions between mobile agents and objects in rich 3D environments. Unique properties include: real-time near-photo-realistic image rendering; a library of objects and environments, and routines for their customization; generative procedu… ▽ More We introduce ThreeDWorld (TDW), a platform for interactive multi-modal physical simulation. TDW enables simulation of high-fidelity sensory data and physical interactions between mobile agents and objects in rich 3D environments. Unique properties include: real-time near-photo-realistic image rendering; a library of objects and environments, and routines for their customization; generative procedures for efficiently building classes of new environments; high-fidelity audio rendering; realistic physical interactions for a variety of material types, including cloths, liquid, and deformable objects; customizable agents that embody AI agents; and support for human interactions with VR devices. TDW's API enables multiple agents to interact within a simulation and returns a range of sensor and physics data representing the state of the world. We present initial experiments enabled by TDW in emerging research directions in computer vision, machine learning, and cognitive science, including multi-modal physical scene understanding, physical dynamics predictions, multi-agent interactions, models that learn like a child, and attention studies in humans and neural networks. △ Less

Submitted 28 December, 2021; v1 submitted 9 July, 2020; originally announced July 2020.

Comments: Oral Presentation at NeurIPS 21 Datasets and Benchmarks Track. Project page: http://www.threedworld.org

Flexible and Efficient Long-Range Planning Through Curious Exploration

Authors: Aidan Curtis, Minjian Xin, Dilip Arumugam, Kevin Feigelis, Daniel Yamins

Abstract: Identifying algorithms that flexibly and efficiently discover temporally-extended multi-phase plans is an essential step for the advancement of robotics and model-based reinforcement learning. The core problem of long-range planning is finding an efficient way to search through the tree of possible action sequences. Existing non-learned planning solutions from the Task and Motion Planning (TAMP) l… ▽ More Identifying algorithms that flexibly and efficiently discover temporally-extended multi-phase plans is an essential step for the advancement of robotics and model-based reinforcement learning. The core problem of long-range planning is finding an efficient way to search through the tree of possible action sequences. Existing non-learned planning solutions from the Task and Motion Planning (TAMP) literature rely on the existence of logical descriptions for the effects and preconditions for actions. This constraint allows TAMP methods to efficiently reduce the tree search problem but limits their ability to generalize to unseen and complex physical environments. In contrast, deep reinforcement learning (DRL) methods use flexible neural-network-based function approximators to discover policies that generalize naturally to unseen circumstances. However, DRL methods struggle to handle the very sparse reward landscapes inherent to long-range multi-step planning situations. Here, we propose the Curious Sample Planner (CSP), which fuses elements of TAMP and DRL by combining a curiosity-guided sampling strategy with imitation learning to accelerate planning. We show that CSP can efficiently discover interesting and complex temporally-extended plans for solving a wide range of physically realistic 3D tasks. In contrast, standard planning and learning methods often fail to solve these tasks at all or do so only with a huge and highly variable number of training samples. We explore the use of a variety of curiosity metrics with CSP and analyze the types of solutions that CSP discovers. Finally, we show that CSP supports task transfer so that the exploration policies learned during experience with one task can help improve efficiency on related tasks. △ Less

Submitted 8 July, 2020; v1 submitted 22 April, 2020; originally announced April 2020.

KELT-25b and KELT-26b: A Hot Jupiter and a Substellar Companion Transiting Young A-stars Observed by TESS

Authors: Romy Rodríguez Martínez, B. Scott Gaudi, Joseph E. Rodriguez, George Zhou, Jonathan Labadie-Bartz, Samuel N. Quinn, Kaloyan Minev Penev, Thiam-Guan Tan, David W. Latham, Leonardo A. Paredes, John Kielkopf, Brett C. Addison, Duncan J. Wright, Johanna K. Teske, Steve B. Howell, David R. Ciardi, Carl Ziegler, Keivan G. Stassun, Marshall C. Johnson, Jason D. Eastman, Robert J. Siverd, Thomas G. Beatty, Luke G. Bouma, Joshua Pepper, Michael B. Lund , et al. (67 additional authors not shown)

Abstract: We present the discoveries of KELT-25b (TIC 65412605, TOI-626.01) and KELT-26b (TIC 160708862, TOI-1337.01), two transiting companions orbiting relatively bright, early A-stars. The transit signals were initially detected by the KELT survey, and subsequently confirmed by \textit{TESS} photometry. KELT-25b is on a 4.40-day orbit around the V = 9.66 star CD-24 5016 (… ▽ More We present the discoveries of KELT-25b (TIC 65412605, TOI-626.01) and KELT-26b (TIC 160708862, TOI-1337.01), two transiting companions orbiting relatively bright, early A-stars. The transit signals were initially detected by the KELT survey, and subsequently confirmed by \textit{TESS} photometry. KELT-25b is on a 4.40-day orbit around the V = 9.66 star CD-24 5016 ($T_{\rm eff} = 8280^{+440}_{-180}$ K, $M_{\star}$ = $2.18^{+0.12}_{-0.11}$ $M_{\odot}$), while KELT-26b is on a 3.34-day orbit around the V = 9.95 star HD 134004 ($T_{\rm eff}$ =$8640^{+500}_{-240}$ K, $M_{\star}$ = $1.93^{+0.14}_{-0.16}$ $M_{\odot}$), which is likely an Am star. We have confirmed the sub-stellar nature of both companions through detailed characterization of each system using ground-based and \textit{TESS} photometry, radial velocity measurements, Doppler Tomography, and high-resolution imaging. For KELT-25, we determine a companion radius of $R_{\rm P}$ = $1.64^{+0.039}_{-0.043}$ $R_{\rm J}$, and a 3-sigma upper limit on the companion's mass of $\sim64~M_{\rm J}$. For KELT-26b, we infer a planetary mass and radius of $M_{\rm P}$ = $1.41^{+0.43}_{-0.51}$ $M_{\rm J}$ and $R_{\rm P}$ = $1.940^{+0.060}_{-0.058}$ $R_{\rm J}$. From Doppler Tomographic observations, we find KELT-26b to reside in a highly misaligned orbit. This conclusion is weakly corroborated by a subtle asymmetry in the transit light curve from the \textit{TESS} data. KELT-25b appears to be in a well-aligned, prograde orbit, and the system is likely a member of a cluster or moving group. △ Less

Submitted 2 December, 2019; originally announced December 2019.

Comments: 24 pages, 18 figures, 8 tables

Probabilistic Neural Network Tomography across Grane field (North Sea) from Surface Wave Dispersion Data

Authors: Stephanie Earp, Andrew Curtis, Xin Zhang, Fredrik Hansteen

Abstract: Surface wave tomography uses measured dispersion properties of surface waves to infer the spatial distribution of subsurface properties such as shear-wave velocities. These properties can be estimated vertically below any geographical location at which surface wave dispersion data are available. As the inversion is significantly non-linear, Monte Carlo methods are often used to invert dispersion c… ▽ More Surface wave tomography uses measured dispersion properties of surface waves to infer the spatial distribution of subsurface properties such as shear-wave velocities. These properties can be estimated vertically below any geographical location at which surface wave dispersion data are available. As the inversion is significantly non-linear, Monte Carlo methods are often used to invert dispersion curves for shear-wave velocity profiles with depth to give a probabilistic solution. Such methods provide uncertainty information but are computationally expensive. Neural network based inversion provides a more efficient way to obtain probabilistic solutions when those solutions are required beneath many geographical locations. Unlike Monte Carlo methods, once a network has been trained it can be applied rapidly to perform any number of inversions. We train a class of neural networks called mixture density networks, to invert dispersion curves for shear-wave velocity models and their non-linearised uncertainty. Mixture density networks are able to produce fully probabilistic solutions in the form of weighted sums of multivariate analytic kernels such as Gaussians, and we show that including data uncertainties in the mixture density network gives more reliable mean velocity estimates when data contains significant noise. The networks were applied to data from the Grane field in the Norwegian North sea to produce shear-wave velocity maps at several depth levels. Post-training we obtained probabilistic velocity profiles with depth beneath 26,772 locations to produce a 3D velocity model in 21 seconds on a standard desktop computer. This method is therefore ideally suited for rapid, repeated 3D subsurface imaging and monitoring. △ Less

Submitted 26 August, 2019; originally announced August 2019.

Comments: Submitted to Geophysical Journal International

Seismic tomography using variational inference methods

Authors: Xin Zhang, Andrew Curtis

Abstract: Seismic tomography is a methodology to image the interior of solid or fluid media, and is often used to map properties in the subsurface of the Earth. In order to better interpret the resulting images it is important to assess imaging uncertainties. Since tomography is significantly nonlinear, Monte Carlo sampling methods are often used for this purpose, but they are generally computationally intr… ▽ More Seismic tomography is a methodology to image the interior of solid or fluid media, and is often used to map properties in the subsurface of the Earth. In order to better interpret the resulting images it is important to assess imaging uncertainties. Since tomography is significantly nonlinear, Monte Carlo sampling methods are often used for this purpose, but they are generally computationally intractable for large datasets and high-dimensional parameter spaces. To extend uncertainty analysis to larger systems we use variational inference methods to conduct seismic tomography. In contrast to Monte Carlo sampling, variational methods solve the Bayesian inference problem as an optimization problem, yet still provide probabilistic results. In this study, we applied two variational methods, automatic differential variational inference (ADVI) and Stein variational gradient descent (SVGD), to 2D seismic tomography problems using both synthetic and real data and we compare the results to those from two different Monte Carlo sampling methods. The results show that variational inference methods can produce accurate approximations to the results of Monte Carlo sampling methods at significantly lower computational cost, provided that gradients of parameters with respect to data can be calculated efficiently. We expect that the methods can be applied fruitfully to many other types of geophysical inverse problems. △ Less

Submitted 27 August, 2019; v1 submitted 22 August, 2019; originally announced August 2019.

Comments: 26 pages, 14 figures

Sign Patterns of Orthogonal Matrices and the Strong Inner Product Property

Authors: Bryan A. Curtis, Bryan L. Shader

Abstract: A new condition, the strong inner product property, is introduced and used to construct sign patterns of row orthogonal matrices. Using this property, infinite families of sign patterns allowing row orthogonality are found. These provide insight into the underlying combinatorial structure of row orthogonal matrices. Algorithmic techniques for verifying that a matrix has the strong inner product pr… ▽ More A new condition, the strong inner product property, is introduced and used to construct sign patterns of row orthogonal matrices. Using this property, infinite families of sign patterns allowing row orthogonality are found. These provide insight into the underlying combinatorial structure of row orthogonal matrices. Algorithmic techniques for verifying that a matrix has the strong inner product property are also presented. These techniques lead to a generalization of the strong inner product property and can be easily implemented using various software. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Search for transient variations of the fine structure constant and dark matter using fiber-linked optical atomic clocks

Authors: B. M. Roberts, P. Delva, A. Al-Masoudi, A. Amy-Klein, C. Bærentsen, C. F. A. Baynham, E. Benkler, S. Bilicki, S. Bize, W. Bowden, J. Calvert, V. Cambier, E. Cantin, E. A. Curtis, S. Dörscher, M. Favier, F. Frank, P. Gill, R. M. Godun, G. Grosche, C. Guo, A. Hees, I. R. Hill, R. Hobson, N. Huntemann , et al. (29 additional authors not shown)

Abstract: We search for transient variations of the fine structure constant using data from a European network of fiber-linked optical atomic clocks. By searching for coherent variations in the recorded clock frequency comparisons across the network, we significantly improve the constraints on transient variations of the fine structure constant. For example, we constrain the variation in alpha to <5*10^-17… ▽ More We search for transient variations of the fine structure constant using data from a European network of fiber-linked optical atomic clocks. By searching for coherent variations in the recorded clock frequency comparisons across the network, we significantly improve the constraints on transient variations of the fine structure constant. For example, we constrain the variation in alpha to <5*10^-17 for transients of duration 10^3 s. This analysis also presents a possibility to search for dark matter, the mysterious substance hypothesised to explain galaxy dynamics and other astrophysical phenomena that is thought to dominate the matter density of the universe. At the current sensitivity level, we find no evidence for dark matter in the form of topological defects (or, more generally, any macroscopic objects), and we thus place constraints on certain potential couplings between the dark matter and standard model particles, substantially improving upon the existing constraints, particularly for large (>~10^4 km) objects. △ Less

Submitted 8 July, 2019; v1 submitted 4 July, 2019; originally announced July 2019.

Journal ref: New J. Phys. 22, 093010 (2020)