Showing 1–50 of 122 results for author: Forbes, J

Nonlinear Interactions of Planetary-Scale Waves in Mesospheric Winds Observed at 52°N Latitude and Two Longitudes

Authors: Maosheng He, Jeffrey M. Forbes, Gunter Stober, Christoph Jacobi, Guozhu Li, Libo Liu, Jiyao Xu

Abstract: Nine years of mesospheric wind data from two meteor radars at 52°N latitude were analyzed to investigate planetary waves (PWs) and tides by estimating their zonal wavenumber through longitudinal phase differences. Our results reveal that PW normal modes (NMs) primarily drive multi-day oscillations, showing seasonal variability and statistical associations with Sudden Stratospheric Warming (SSW) ev… ▽ More Nine years of mesospheric wind data from two meteor radars at 52°N latitude were analyzed to investigate planetary waves (PWs) and tides by estimating their zonal wavenumber through longitudinal phase differences. Our results reveal that PW normal modes (NMs) primarily drive multi-day oscillations, showing seasonal variability and statistical associations with Sudden Stratospheric Warming (SSW) events. Specifically, a significant 6-day NM emerges in April, followed by predominant 4- and 2-day NMs until June, with peaks of 2-, 4-, and 6-day NMs spanning July to October. Furthermore, our study provides the first observational verification of frequency and zonal wavenumber of over ten secondary waves from nonlinear interactions among planetary-scale waves. One notable finding is the prevalence of non-migrating components in winter 24-hour and summer 8-hour tides, attributed to these nonlinear interactions. Our findings underscore the diverse nonlinear dynamics of planetary-scale waves, triggering a variety of periodic oscillations. △ Less

Submitted 11 June, 2024; v1 submitted 9 June, 2024; originally announced June 2024.

Validating Automated Resonance Evaluation with Synthetic Data

Authors: Oleksii Zivenko, Noah A. W. Walton, William Fritsch, Jacob Forbes, Amanda M. Lewis, Aaron Clark, Jesse M. Brown, Vladimir Sobes

Abstract: The integrity and precision of nuclear data are crucial for a broad spectrum of applications, from national security and nuclear reactor design to medical diagnostics, where the associated uncertainties can significantly impact outcomes. A substantial portion of uncertainty in nuclear data originates from the subjective biases in the evaluation process, a crucial phase in the nuclear data producti… ▽ More The integrity and precision of nuclear data are crucial for a broad spectrum of applications, from national security and nuclear reactor design to medical diagnostics, where the associated uncertainties can significantly impact outcomes. A substantial portion of uncertainty in nuclear data originates from the subjective biases in the evaluation process, a crucial phase in the nuclear data production pipeline. Recent advancements indicate that automation of certain routines can mitigate these biases, thereby standardizing the evaluation process, reducing uncertainty and enhancing reproducibility. This article contributes to developing a framework for automated evaluation techniques testing, emphasizing automated fitting methods that do not require the user to provide any prior information. This approach simplifies the process and reduces the manual effort needed in the initial evaluation stage. It highlights the capability of the framework to validate and optimize subroutines, targeting the performance analysis and optimization of the fitting procedure using high-fidelity synthetic data (labeled experimental data) and the concept of a fully controlled computational experiment. An error metric is introduced to provide a clear and intuitive measure of the fitting quality by quantifying the accuracy and performance across the specified energy. This metric sets a scale for comparison and optimization of routines or hyperparameter selection, improving the entire evaluation process methodology and increasing reproducibility and objectivity. △ Less

Submitted 18 June, 2024; v1 submitted 3 June, 2024; originally announced June 2024.

Comments: 39 pages, 12 figures; As a follow-up to arXiv:2303.09698 and arXiv:2402.14122

MSC Class: 65K10; 62H12 ACM Class: J.2; G.1.6

A Hessian for Gaussian Mixture Likelihoods in Nonlinear Least Squares

Authors: Vassili Korotkine, Mitchell Cohen, James Richard Forbes

Abstract: This paper proposes a novel Hessian approximation for Maximum a Posteriori estimation problems in robotics involving Gaussian mixture likelihoods. The proposed Hessian leads to better convergence properties. Previous approaches manipulate the Gaussian mixture likelihood into a form that allows the problem to be represented as a nonlinear least squares (NLS) problem. However, they result in an inac… ▽ More This paper proposes a novel Hessian approximation for Maximum a Posteriori estimation problems in robotics involving Gaussian mixture likelihoods. The proposed Hessian leads to better convergence properties. Previous approaches manipulate the Gaussian mixture likelihood into a form that allows the problem to be represented as a nonlinear least squares (NLS) problem. However, they result in an inaccurate Hessian approximation due to additional nonlinearities that are not accounted for in NLS solvers. The proposed Hessian approximation is derived by setting the Hessians of the Gaussian mixture component errors to zero, which is the same starting point as for the Gauss-Newton Hessian approximation for NLS, and using the chain rule to account for additional nonlinearities. The proposed Hessian approximation is more accurate, resulting in improved convergence properties that are demonstrated on simulated and real-world experiments. A method to maintain compatibility with existing solvers, such as ceres, is also presented. Accompanying software and supplementary material can be found at https://github.com/decargroup/hessian_sum_mixtures. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 8 pages, 2 figures. Submitted to IEEE Robotics and Automation Letters

Forward-Backward Extended DMD with an Asymptotic Stability Constraint

Authors: Louis Lortie, Steven Dahdah, James Richard Forbes

Abstract: This paper presents a data-driven method to identify an asymptotically stable Koopman system from noisy data. In particular, the proposed approach combines approximations of the system's forward- and backward-in-time dynamics to reduce bias caused by noisy data while enforcing asymptotic stability. A Koopman model of an inherently asymptotically stable system can be unstable due to noisy data and… ▽ More This paper presents a data-driven method to identify an asymptotically stable Koopman system from noisy data. In particular, the proposed approach combines approximations of the system's forward- and backward-in-time dynamics to reduce bias caused by noisy data while enforcing asymptotic stability. A Koopman model of an inherently asymptotically stable system can be unstable due to noisy data and a poor choice of lifting functions. To prevent identifying an unstable model, the proposed approach imposes an asymptotic stability constraint on the Koopman model. The proposed method is formulated as a semidefinite program and its performance is compared to state-of-the-art methods with a simulated Duffing oscillator dataset and experimental soft robot dataset. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 20 pages, 7 figures

Optimal Robot Formations: Balancing Range-Based Observability and User-Defined Configurations

Authors: Syed Shabbir Ahmed, Mohammed Ayman Shalaby, Jerome Le Ny, James Richard Forbes

Abstract: This paper introduces a set of customizable and novel cost functions that enable the user to easily specify desirable robot formations, such as a ``high-coverage'' infrastructure-inspection formation, while maintaining high relative pose estimation accuracy. The overall cost function balances the need for the robots to be close together for good ranging-based relative localization accuracy and the… ▽ More This paper introduces a set of customizable and novel cost functions that enable the user to easily specify desirable robot formations, such as a ``high-coverage'' infrastructure-inspection formation, while maintaining high relative pose estimation accuracy. The overall cost function balances the need for the robots to be close together for good ranging-based relative localization accuracy and the need for the robots to achieve specific tasks, such as minimizing the time taken to inspect a given area. The formations found by minimizing the aggregated cost function are evaluated in a coverage path planning task in simulation and experiment, where the robots localize themselves and unknown landmarks using a simultaneous localization and mapping algorithm based on the extended Kalman filter. Compared to an optimal formation that maximizes ranging-based relative localization accuracy, these formations significantly reduce the time to cover a given area with minimal impact on relative pose estimation accuracy. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 8 pages, 9 figures, submitted to IEEE International Conference on Intelligent Robots and Systems 2024

Navigation and Control of Unconventional VTOL UAVs in Forward-Flight with Explicit Wind Velocity Estimation

Authors: Mitchell Cohen, James Richard Forbes

Abstract: This paper presents a solution for the state estimation and control problems for a class of unconventional vertical takeoff and landing (VTOL) UAVs operating in forward-flight conditions. A tightly-coupled state estimation approach is used to estimate the aircraft navigation states, sensor biases, and the wind velocity. State estimation is done within a matrix Lie group framework using the Invaria… ▽ More This paper presents a solution for the state estimation and control problems for a class of unconventional vertical takeoff and landing (VTOL) UAVs operating in forward-flight conditions. A tightly-coupled state estimation approach is used to estimate the aircraft navigation states, sensor biases, and the wind velocity. State estimation is done within a matrix Lie group framework using the Invariant Extended Kalman Filter (IEKF), which offers several advantages compared to standard multiplicative EKFs traditionally used in aerospace and robotics problems. An SO(3)- based attitude controller is employed, leading to a single attitude control law without a separate sideslip control loop. A control allocator is used to determine how to use multiple, possibly redundant, actuators to produce the desired control moments. The wind velocity estimates are used in the attitude controller and the control allocator to improve performance. A numerical example is considered using a sample VTOL tailsitter-type UAV with four control surfaces. Monte-Carlo simulations demonstrate robustness of the proposed control and estimation scheme to various initial conditions, noise levels, and flight trajectories. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 8 pages, 7 figures, published in Robotics and Automation Letters

Journal ref: IEEE Robotics and Automation Letters, vol. 5, no. 2, pp. 1151-1158, June 2020

The Invariant Rauch-Tung-Striebel Smoother

Authors: Niels van der Laan, Mitchell Cohen, Jonathan Arsenault, James Richard Forbes

Abstract: This paper presents an invariant Rauch-Tung- Striebel (IRTS) smoother applicable to systems with states that are an element of a matrix Lie group. In particular, the extended Rauch-Tung-Striebel (RTS) smoother is adapted to work within a matrix Lie group framework. The main advantage of the invariant RTS (IRTS) smoother is that the linearization of the process and measurement models is independent… ▽ More This paper presents an invariant Rauch-Tung- Striebel (IRTS) smoother applicable to systems with states that are an element of a matrix Lie group. In particular, the extended Rauch-Tung-Striebel (RTS) smoother is adapted to work within a matrix Lie group framework. The main advantage of the invariant RTS (IRTS) smoother is that the linearization of the process and measurement models is independent of the state estimate resulting in state-estimate-independent Jacobians when certain technical requirements are met. A sample problem is considered that involves estimation of the three dimensional pose of a rigid body on SE(3), along with sensor biases. The multiplicative RTS (MRTS) smoother is also reviewed and is used as a direct comparison to the proposed IRTS smoother using experimental data. Both smoothing methods are also compared to invariant and multiplicative versions of the Gauss-Newton approach to solving the batch state estimation problem. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures, published in Robotics and Automation Letters

Journal ref: IEEE Robotics and Automation Letters, vol. 5, no. 4, pp 5067-5074, June 2020

Laser-to-Vehicle Extrinsic Calibration in Low-Observability Scenarios for Subsea Mapping

Authors: Thomas Hitchcox, James Richard Forbes

Abstract: Laser line scanners are increasingly being used in the subsea industry for high-resolution mapping and infrastructure inspection. However, calibrating the 3D pose of the scanner relative to the vehicle is a perennial source of confusion and frustration for industrial surveyors. This work describes three novel algorithms for laser-to-vehicle extrinsic calibration using naturally occurring features.… ▽ More Laser line scanners are increasingly being used in the subsea industry for high-resolution mapping and infrastructure inspection. However, calibrating the 3D pose of the scanner relative to the vehicle is a perennial source of confusion and frustration for industrial surveyors. This work describes three novel algorithms for laser-to-vehicle extrinsic calibration using naturally occurring features. Each algorithm makes a different assumption on the quality of the vehicle trajectory estimate, enabling good calibration results in a wide range of situations. A regularization technique is used to address low-observability scenarios frequently encountered in practice with large, rotationally stable subsea vehicles. Experimental results are provided for two field datasets, including the recently discovered wreck of the Endurance. △ Less

Submitted 8 March, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 8 pages, 8 figures. Accepted for publication in IEEE Robotics and Automation Letters. Updated to include vol., no., and page nos

Automated Resonance Identification in Nuclear Data Evaluation

Authors: Noah A. W. Walton, Oleksii Zivenko, William Fritsch, Jacob Forbes, Amanda Lewis, Jesse Brown, Vlad Sobes

Abstract: Global and national efforts to deliver high-quality nuclear data to users have a broad impact across applications such as national security, reactor operation, basic science, medical fields, and more. Cross section evaluation is a large part this effort as it combines theory and experiment to produce suggested values and uncertainty for reaction probabilities. In most isotopes, the cross section e… ▽ More Global and national efforts to deliver high-quality nuclear data to users have a broad impact across applications such as national security, reactor operation, basic science, medical fields, and more. Cross section evaluation is a large part this effort as it combines theory and experiment to produce suggested values and uncertainty for reaction probabilities. In most isotopes, the cross section exhibits resonant behavior in what is called the resonance region of incident neutron energy. Resonance region evaluation is a specialized type of nuclear data evaluation that can require significant, manual effort and months of time from expert scientists. In this article, non-convex, non-linear optimization methods are combined with concepts of inferential statistics to infer a set of optimized resonance models from experimental data in an automated manner that is not dependent on prior evaluation(s). This methodology aims to enhance the workflow of a resonance evaluator by minimizing time, effort, and prior biases while improving reproducibility and document-ability, addressing widely recognized challenges in the field. △ Less

Submitted 21 February, 2024; originally announced February 2024.

Gaussian-Sum Filter for Range-based 3D Relative Pose Estimation in the Presence of Ambiguities

Authors: Syed S. Ahmed, Mohammed A. Shalaby, Charles C. Cossette, Jerome Le Ny, James R. Forbes

Abstract: Multi-robot systems must have the ability to accurately estimate relative states between robots in order to perform collaborative tasks, possibly with no external aiding. Three-dimensional relative pose estimation using range measurements oftentimes suffers from a finite number of non-unique solutions, or ambiguities. This paper: 1) identifies and accurately estimates all possible ambiguities in 2… ▽ More Multi-robot systems must have the ability to accurately estimate relative states between robots in order to perform collaborative tasks, possibly with no external aiding. Three-dimensional relative pose estimation using range measurements oftentimes suffers from a finite number of non-unique solutions, or ambiguities. This paper: 1) identifies and accurately estimates all possible ambiguities in 2D; 2) treats them as components of a Gaussian mixture model; and 3) presents a computationally-efficient estimator, in the form of a Gaussian-sum filter (GSF), to realize range-based relative pose estimation in an infrastructure-free, 3D, setup. This estimator is evaluated in simulation and experiment and is shown to avoid divergence to local minima induced by the ambiguous poses. Furthermore, the proposed GSF outperforms an extended Kalman filter, demonstrates similar performance to the computationally-demanding particle filter, and is shown to be consistent. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 7 pages, 9 figures, submitted to IEEE Conference on Control Technology and Applications 2024

Monitoring the evolution of relative product populations at early times during a photochemical reaction

Authors: Joao Pedro Figueira Nunes, Lea Maria Ibele, Shashank Pathak, Andrew R. Attar, Surjendu Bhattacharyya, Rebecca Boll, Kurtis Borne, Martin Centurion, Benjamin Erk, Ming-Fu Lin, Ruaridh J. G. Forbes, Nate Goff, Christopher S. Hansen, Matthias Hoffmann, David M. P. Holland, Rebecca A. Ingle, Duan Luo, Sri Bhavya Muvva, Alex Reid, Arnaud Rouzée, Artem Rudenko, Sajib Kumar Saha, Xiaozhe Shen, Anbu Selvam Venkatachalam, Xijie Wang , et al. (9 additional authors not shown)

Abstract: Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps towards understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species amongst the emerging re… ▽ More Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps towards understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species amongst the emerging reaction products. Here, we show that mega-electronvolt ultrafast electron diffraction in combination with ab initio molecular dynamics calculations offer a powerful route to determining time-resolved populations of the various isomeric products formed after UV (266 nm) excitation of the five-membered heterocyclic molecule 2(5H)-thiophenone. This strategy provides experimental validation of the predicted high (~50%) yield of an episulfide isomer containing a strained 3-membered ring within ~1 ps of photoexcitation and highlights the rapidity of interconversion between the rival highly vibrationally excited photoproducts in their ground electronic state. △ Less

Submitted 21 November, 2023; originally announced November 2023.

navlie: A Python Package for State Estimation on Lie Groups

Authors: Charles Champagne Cossette, Mitchell Cohen, Vassili Korotkine, Arturo del Castillo Bernal, Mohammed Ayman Shalaby, James Richard Forbes

Abstract: The ability to rapidly test a variety of algorithms for an arbitrary state estimation task is valuable in the prototyping phase of navigation systems. Lie group theory is now mainstream in the robotics community, and hence estimation prototyping tools should allow state definitions that belong to manifolds. A new package, called navlie, provides a framework that allows a user to model a large clas… ▽ More The ability to rapidly test a variety of algorithms for an arbitrary state estimation task is valuable in the prototyping phase of navigation systems. Lie group theory is now mainstream in the robotics community, and hence estimation prototyping tools should allow state definitions that belong to manifolds. A new package, called navlie, provides a framework that allows a user to model a large class of problems by implementing a set of classes complying with a generic interface. Once accomplished, navlie provides a variety of on-manifold estimation algorithms that can run directly on these classes. The package also provides a built-in library of common models, as well as many useful utilities. The open-source project can be found at https://github.com/decargroup/navlie. △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: 6 pages, 8 figures, presented at the 2023 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS)

Why do semi-analytic models predict higher scatter in the stellar mass-halo mass relation than cosmological hydrodynamic simulations?

Authors: Antonio J. Porras-Valverde, John C. Forbes, Rachel S. Somerville, Adam R. H. Stevens, Kelly Holley-Bockelmann, Andreas A. Berlind, Shy Genel

Abstract: Semi-analytic models (SAMs) systematically predict higher stellar-mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAM Dark Sage. We design two black hole formation models that approximate results from the IllustrisTNG 300-1 hydrodynamical simu… ▽ More Semi-analytic models (SAMs) systematically predict higher stellar-mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAM Dark Sage. We design two black hole formation models that approximate results from the IllustrisTNG 300-1 hydrodynamical simulation. In the first model, we assign a fixed black hole mass of $10^{6}\, \mathrm{M}_{\odot}$ to every halo that reaches $10^{10.5}\, \mathrm{M}_{\odot}$. In the second model, we disregard any black hole growth as implemented in the standard Dark Sage model. Instead, we force all black hole masses to follow the median black hole mass-halo mass relation in IllustrisTNG 300-1 with a fixed scatter. We find that each model on its own does not significantly reduce the scatter in stellar mass. To do this, we replace the native Dark Sage AGN feedback model with a simple model where we turn off cooling for galaxies with black hole masses above $10^{8}\, \mathrm{M}_{\odot}$. With this additional modification, the SMBH seeding and fixed conditional distribution models find a significant reduction in the scatter in stellar mass at halo masses between $10^{11-14}\, \mathrm{M}_{\odot}$. These results suggest that AGN feedback in SAMs acts in a qualitatively different way than feedback implemented in cosmological simulations. Either or both may require substantial modification to match the empirically inferred scatter in the Stellar Mass Halo Mass Relation (SMHMR). △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 21 pages, 16 figures

Bayesian Filtering for Homography Estimation

Authors: Arturo Del Castillo Bernal, Philippe Decoste, James Richard Forbes

Abstract: This paper considers homography estimation in a Bayesian filtering framework using rate gyro and camera measurements. The use of rate gyro measurements facilitates a more reliable estimate of homography in the presence of occlusions, while a Bayesian filtering approach generates both a homography estimate along with an uncertainty. Uncertainty information opens the door to adaptive filtering appro… ▽ More This paper considers homography estimation in a Bayesian filtering framework using rate gyro and camera measurements. The use of rate gyro measurements facilitates a more reliable estimate of homography in the presence of occlusions, while a Bayesian filtering approach generates both a homography estimate along with an uncertainty. Uncertainty information opens the door to adaptive filtering approaches, post-processing procedures, and safety protocols. In particular, herein an iterative extended Kalman filter and an interacting multiple model (IMM) filter are tested using both simulated and experimental datasets. The IMM is shown to have good consistency properties and better overall performance when compared to the state-of-the-art homography nonlinear deterministic observer in both simulations and experiments. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 8 pages and 6 figures, to be published in IEEE Robotics and Automation Letters

IMU Preintegration for Multi-Robot Systems in the Presence of Bias and Communication Constraints

Authors: Mohammed Ayman Shalaby, Charles Champagne Cossette, Jerome Le Ny, James Richard Forbes

Abstract: This document is in supplement to the paper titled "Multi-Robot Relative Pose Estimation and IMU Preintegration Using Passive UWB Transceivers", available at [1]. The purpose of this document is to show how IMU biases can be incorporated into the framework presented in [1], while maintaining the differential Sylvester equation form of the process model. This document is in supplement to the paper titled "Multi-Robot Relative Pose Estimation and IMU Preintegration Using Passive UWB Transceivers", available at [1]. The purpose of this document is to show how IMU biases can be incorporated into the framework presented in [1], while maintaining the differential Sylvester equation form of the process model. △ Less

Submitted 16 October, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

STAR-loc: Dataset for STereo And Range-based localization

Authors: Frederike Dümbgen, Mohammed A. Shalaby, Connor Holmes, Charles C. Cossette, James R. Forbes, Jerome Le Ny, Timothy D. Barfoot

Abstract: This document contains a detailed description of the STAR-loc dataset. For a quick starting guide please refer to the associated Github repository (https://github.com/utiasASRL/starloc). The dataset consists of stereo camera data (rectified/raw images and inertial measurement unit measurements) and ultra-wideband (UWB) data (range measurements) collected on a sensor rig in a Vicon motion capture a… ▽ More This document contains a detailed description of the STAR-loc dataset. For a quick starting guide please refer to the associated Github repository (https://github.com/utiasASRL/starloc). The dataset consists of stereo camera data (rectified/raw images and inertial measurement unit measurements) and ultra-wideband (UWB) data (range measurements) collected on a sensor rig in a Vicon motion capture arena. The UWB anchors and visual landmarks (Apriltags) are of known position, so the dataset can be used for both localization and Simultaneous Localization and Mapping (SLAM). △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 15 pages, 15 figures

The Global Structure of Molecular Clouds: I. Trends with Mass and Star Formation Rate

Authors: Nia Imara, John C. Forbes

Abstract: We introduce a model for the large-scale, global 3D structure of molecular clouds. Motivated by the morphological appearance of clouds in surface density maps, we model clouds as cylinders, with the aim of backing out information about the volume density distribution of gas and its relationship to star formation. We test our model by applying it to surface density maps for a sample of nearby cloud… ▽ More We introduce a model for the large-scale, global 3D structure of molecular clouds. Motivated by the morphological appearance of clouds in surface density maps, we model clouds as cylinders, with the aim of backing out information about the volume density distribution of gas and its relationship to star formation. We test our model by applying it to surface density maps for a sample of nearby clouds and find solutions that fit each of the observed radial surface density profiles remarkably well. Our most salient findings are that clouds with higher central volume densities are more compact and also have lower total mass. These same lower-mass clouds tend to have shorter gas depletion times, regardless of whether we consider their total mass or dense mass. Our analyses lead us to conclude that cylindrical clouds can be characterized by a universal structure that sets the timescale on which they form stars. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: Accepted to ApJ. 12 pages, 9 figures

Data-Driven Batch Localization and SLAM Using Koopman Linearization

Authors: Zi Cong Guo, Frederike Dümbgen, James R. Forbes, Timothy D. Barfoot

Abstract: We present a framework for model-free batch localization and SLAM. We use lifting functions to map a control-affine system into a high-dimensional space, where both the process model and the measurement model are rendered bilinear. During training, we solve a least-squares problem using groundtruth data to compute the high-dimensional model matrices associated with the lifted system purely from da… ▽ More We present a framework for model-free batch localization and SLAM. We use lifting functions to map a control-affine system into a high-dimensional space, where both the process model and the measurement model are rendered bilinear. During training, we solve a least-squares problem using groundtruth data to compute the high-dimensional model matrices associated with the lifted system purely from data. At inference time, we solve for the unknown robot trajectory and landmarks through an optimization problem, where constraints are introduced to keep the solution on the manifold of the lifting functions. The problem is efficiently solved using a sequential quadratic program (SQP), where the complexity of an SQP iteration scales linearly with the number of timesteps. Our algorithms, called Reduced Constrained Koopman Linearization Localization (RCKL-Loc) and Reduced Constrained Koopman Linearization SLAM (RCKL-SLAM), are validated experimentally in simulation and on two datasets: one with an indoor mobile robot equipped with a laser rangefinder that measures range to cylindrical landmarks, and one on a golf cart equipped with RFID range sensors. We compare RCKL-Loc and RCKL-SLAM with classic model-based nonlinear batch estimation. While RCKL-Loc and RCKL-SLAM have similar performance compared to their model-based counterparts, they outperform the model-based approaches when the prior model is imperfect, showing the potential benefit of the proposed data-driven technique. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: Submitted to IEEE T-RO. 18 pages, 9 figures, 1 table

Magnetic Navigation using Attitude-Invariant Magnetic Field Information for Loop Closure Detection

Authors: Natalia Pavlasek, Charles Champagne Cossette, David Roy-Guay, James Richard Forbes

Abstract: Indoor magnetic fields are a combination of Earth's magnetic field and disruptions induced by ferromagnetic objects, such as steel structural components in buildings. As a result of these disruptions, pervasive in indoor spaces, magnetic field data is often omitted from navigation algorithms in indoor environments. This paper leverages the spatially-varying disruptions to Earth's magnetic field to… ▽ More Indoor magnetic fields are a combination of Earth's magnetic field and disruptions induced by ferromagnetic objects, such as steel structural components in buildings. As a result of these disruptions, pervasive in indoor spaces, magnetic field data is often omitted from navigation algorithms in indoor environments. This paper leverages the spatially-varying disruptions to Earth's magnetic field to extract positional information for use in indoor navigation algorithms. The algorithm uses a rate gyro and an array of four magnetometers to estimate the robot's pose. Additionally, the magnetometer array is used to compute attitude-invariant measurements associated with the magnetic field and its gradient. These measurements are used to detect loop closure points. Experimental results indicate that the proposed approach can estimate the pose of a ground robot in an indoor environment within meter accuracy. △ Less

Submitted 5 September, 2023; originally announced September 2023.

The Galactic Interstellar Object Population: A Framework for Prediction and Inference

Authors: Matthew J. Hopkins, Chris Lintott, Michele T. Bannister, J. Ted Mackereth, John C. Forbes

Abstract: The Milky Way is thought to host a huge population of interstellar objects (ISOs), numbering approximately $10^{15}\mathrm{pc}^{-3}$ around the Sun, which are formed and shaped by a diverse set of processes ranging from planet formation to galactic dynamics. We define a novel framework: firstly to predict the properties of this Galactic ISO population by combining models of processes across planet… ▽ More The Milky Way is thought to host a huge population of interstellar objects (ISOs), numbering approximately $10^{15}\mathrm{pc}^{-3}$ around the Sun, which are formed and shaped by a diverse set of processes ranging from planet formation to galactic dynamics. We define a novel framework: firstly to predict the properties of this Galactic ISO population by combining models of processes across planetary and galactic scales, and secondly to make inferences about the processes modelled, by comparing the predicted population to what is observed. We predict the spatial and compositional distribution of the Galaxy's population of ISOs by modelling the Galactic stellar population with data from the APOGEE survey and combining this with a protoplanetary disk chemistry model. Selecting ISO water mass fraction as an example observable quantity, we evaluate its distribution both at the position of the Sun and averaged over the Galactic disk; our prediction for the Solar neighbourhood is compatible with the inferred water mass fraction of 2I/Borisov. We show that the well-studied Galactic stellar metallicity gradient has a corresponding ISO compositional gradient. We also demonstrate the inference part of the framework by using the current observed ISO composition distribution to constrain the parent star metallicity dependence of the ISO production rate. This constraint, and other inferences made with this framework, will improve dramatically as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) progresses and more ISOs are observed. Finally, we explore generalisations of this framework to other Galactic populations, such as that of exoplanets. △ Less

Submitted 15 November, 2023; v1 submitted 10 August, 2023; originally announced August 2023.

Comments: Accepted to AJ

Journal ref: AJ 166 241 (2023)

Combining DVL-INS and Laser-Based Loop Closures in a Batch Estimation Framework for Underwater Positioning

Authors: Amro Al-Baali, Thomas Hitchcox, James Richard Forbes

Abstract: Correcting gradual position drift is a challenge in long-term subsea navigation. Though highly accurate, modern inertial navigation system (INS) estimates will drift over time due to the accumulated effects of sensor noise and biases, even with acoustic aiding from a Doppler velocity log (DVL). The raw sensor measurements and estimation algorithms used by the DVL-aided INS are often proprietary, w… ▽ More Correcting gradual position drift is a challenge in long-term subsea navigation. Though highly accurate, modern inertial navigation system (INS) estimates will drift over time due to the accumulated effects of sensor noise and biases, even with acoustic aiding from a Doppler velocity log (DVL). The raw sensor measurements and estimation algorithms used by the DVL-aided INS are often proprietary, which restricts the fusion of additional sensors that could bound navigation drift over time. In this letter, the raw sensor measurements and their respective covariances are estimated from the DVL-aided INS output using semidefinite programming tools. The estimated measurements are then augmented with laser-based loop-closure measurements in a batch state estimation framework to correct planar position errors. The heading uncertainty from the DVL-aided INS is also considered in the estimation of the updated positions. The pipeline is tested in simulation and on experimental field data. The proposed methodology reduces the long-term navigation drift by more than 30 times compared to the DVL-aided INS estimate. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 17 pages, 13 figures. Accepted for publication at the IEEE Journal of Oceanic Engineering

Synthesizing Control Laws from Data using Sum-of-Squares Optimization

Authors: Jason J. Bramburger, Steven Dahdah, James Richard Forbes

Abstract: The control Lyapunov function (CLF) approach to nonlinear control design is well established. Moreover, when the plant is control affine and polynomial, sum-of-squares (SOS) optimization can be used to find a polynomial controller as a solution to a semidefinite program. This letter considers the use of data-driven methods to design a polynomial controller by leveraging Koopman operator theory, CL… ▽ More The control Lyapunov function (CLF) approach to nonlinear control design is well established. Moreover, when the plant is control affine and polynomial, sum-of-squares (SOS) optimization can be used to find a polynomial controller as a solution to a semidefinite program. This letter considers the use of data-driven methods to design a polynomial controller by leveraging Koopman operator theory, CLFs, and SOS optimization. First, Extended Dynamic Mode Decomposition (EDMD) is used to approximate the Lie derivative of a given CLF candidate with polynomial lifting functions. Then, the polynomial Koopman model of the Lie derivative is used to synthesize a polynomial controller via SOS optimization. The result is a flexible data-driven method that skips the intermediary process of system identification and can be applied widely to control problems. The proposed approach is used to successfully synthesize a controller to stabilize an inverted pendulum on a cart. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Cascaded Model Predictive Control of a Tandem-Rotor Helicopter

Authors: Faraaz Ahmed, Ludwik Sobiesiak, James Richard Forbes

Abstract: This letter considers cascaded model predictive control (MPC) as a computationally lightweight method for controlling a tandem-rotor helicopter. A traditional single MPC structure is split into separate outer and inner-loops. The outer-loop MPC uses an $SE_2(3)$ error to linearize the translational dynamics about a reference trajectory. The inner-loop MPC uses the optimal angular velocity sequence… ▽ More This letter considers cascaded model predictive control (MPC) as a computationally lightweight method for controlling a tandem-rotor helicopter. A traditional single MPC structure is split into separate outer and inner-loops. The outer-loop MPC uses an $SE_2(3)$ error to linearize the translational dynamics about a reference trajectory. The inner-loop MPC uses the optimal angular velocity sequence of the outer-loop MPC to linearize the rotational dynamics. The outer-loop MPC is run at a slower rate than the inner-loop allowing for longer prediction time and improved performance. Monte-Carlo simulations demonstrate robustness to model uncertainty and environmental disturbances. The proposed control structure is benchmarked against a single MPC algorithm where it shows significant improvements in position and velocity tracking while using significantly less computational resources. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 6 pages, 3 figures

MSC Class: 9306 (Primary); 4906 (Secondary)

Journal ref: in IEEE Control Systems Letters, vol. 7, pp. 1345-1350, 2023

An Adaptive Graduated Nonconvexity Loss Function for Robust Nonlinear Least Squares Solutions

Authors: Kyungmin Jung, Thomas Hitchcox, James Richard Forbes

Abstract: Many problems in robotics, such as estimating the state from noisy sensor data or aligning two point clouds, can be posed and solved as least-squares problems. Unfortunately, vanilla nonminimal solvers for least-squares problems are notoriously sensitive to outliers. As such, various robust loss functions have been proposed to reduce the sensitivity to outliers. Examples of loss functions include… ▽ More Many problems in robotics, such as estimating the state from noisy sensor data or aligning two point clouds, can be posed and solved as least-squares problems. Unfortunately, vanilla nonminimal solvers for least-squares problems are notoriously sensitive to outliers. As such, various robust loss functions have been proposed to reduce the sensitivity to outliers. Examples of loss functions include pseudo-Huber, Cauchy, and Geman-McClure. Recently, these loss functions have been generalized into a single loss function that enables the best loss function to be found adaptively based on the distribution of the residuals. However, even with the generalized robust loss function, most nonminimal solvers can only be solved locally given a prior state estimate due to the nonconvexity of the problem. The first contribution of this paper is to combine graduated nonconvexity (GNC) with the generalized robust loss function to solve least-squares problems without a prior state estimate and without the need to specify a loss function. Moreover, existing loss functions, including the generalized loss function, are based on Gaussian-like distribution. However, residuals are often defined as the squared norm of a multivariate error and distributed in a Chi-like fashion. The second contribution of this paper is to apply a norm-aware adaptive robust loss function within a GNC framework. The proposed approach enables a GNC formulation of a generalized loss function such that GNC can be readily applied to a wider family of loss functions. Furthermore, simulations and experiments demonstrate that the proposed method is more robust compared to non-GNC counterparts, and yields faster convergence times compared to other GNC formulations. △ Less

Submitted 10 May, 2024; v1 submitted 11 May, 2023; originally announced May 2023.

Decentralized State Estimation: An Approach using Pseudomeasurements and Preintegration

Authors: Charles Champagne Cossette, Mohammed Ayman Shalaby, David Saussié, James Richard Forbes

Abstract: This paper addresses the problem of decentralized, collaborative state estimation in robotic teams. In particular, this paper considers problems where individual robots estimate similar physical quantities, such as each other's position relative to themselves. The use of pseudomeasurements is introduced as a means of modelling such relationships between robots' state estimates, and is shown to be… ▽ More This paper addresses the problem of decentralized, collaborative state estimation in robotic teams. In particular, this paper considers problems where individual robots estimate similar physical quantities, such as each other's position relative to themselves. The use of pseudomeasurements is introduced as a means of modelling such relationships between robots' state estimates, and is shown to be a tractable way to approach the decentralized state estimation problem. Moreover, this formulation easily leads to a general-purpose observability test that simultaneously accounts for measurements that robots collect from their own sensors, as well as the communication structure within the team. Finally, input preintegration is proposed as a communication-efficient way of sharing odometry information between robots, and the entire theory is appropriate for both vector-space and Lie-group state definitions. To overcome the need for communicating preintegrated-covariance information, a deep autoencoder is proposed that reconstructs the covariance information from the inputs, hence further reducing the communication requirements. The proposed framework is evaluated on three different simulated problems, and one experiment involving three quadcopters. △ Less

Submitted 22 January, 2024; v1 submitted 8 April, 2023; originally announced April 2023.

Comments: 15 pages, 13 figures, submitted to IJRR

Multi-Robot Relative Pose Estimation and IMU Preintegration Using Passive UWB Transceivers

Authors: Mohammed Ayman Shalaby, Charles Champagne Cossette, Jerome Le Ny, James Richard Forbes

Abstract: Ultra-wideband (UWB) systems are becoming increasingly popular as a means of inter-robot ranging and communication. A major constraint associated with UWB is that only one pair of UWB transceivers can range at a time to avoid interference, hence hindering the scalability of UWB-based localization. In this paper, a ranging protocol is proposed that allows all robots to passively listen on neighbour… ▽ More Ultra-wideband (UWB) systems are becoming increasingly popular as a means of inter-robot ranging and communication. A major constraint associated with UWB is that only one pair of UWB transceivers can range at a time to avoid interference, hence hindering the scalability of UWB-based localization. In this paper, a ranging protocol is proposed that allows all robots to passively listen on neighbouring communicating robots without any hierarchical restrictions on the role of the robots. This is utilized to allow each robot to obtain more range measurements and to broadcast preintegrated inertial measurement unit (IMU) measurements for relative extended pose state estimation directly on SE_2(3). Consequently, a simultaneous clock-synchronization and relative-pose estimator (CSRPE) is formulated using an on-manifold extended Kalman filter (EKF) and is evaluated in simulation using Monte-Carlo runs for up to 7 robots. The ranging protocol is implemented in C on custom-made UWB boards fitted to 3 quadcopters, and the proposed filter is evaluated over multiple experimental trials, yielding up to 56% improvement in localization accuracy. △ Less

Submitted 11 April, 2023; v1 submitted 7 April, 2023; originally announced April 2023.

The interplay between feedback, accretion, transport and winds in setting gas-phase metal distribution in galaxies

Authors: Piyush Sharda, Omri Ginzburg, Mark R. Krumholz, John C. Forbes, Emily Wisnioski, Matilde Mingozzi, Henry R. M. Zovaro, Avishai Dekel

Abstract: The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic win… ▽ More The recent decade has seen an exponential growth in spatially-resolved metallicity measurements in the interstellar medium (ISM) of galaxies. To first order, these measurements are characterised by the slope of the radial metallicity profile, known as the metallicity gradient. In this work, we model the relative role of star formation feedback, gas transport, cosmic gas accretion, and galactic winds in driving radial metallicity profiles and setting the mass-metallicity gradient relation (MZGR). We include a comprehensive treatment of these processes by including them as sources that supply mass, metals, and energy to marginally unstable galactic discs in pressure and energy balance. We show that both feedback and accretion that can drive turbulence and enhance metal-mixing via diffusion are crucial to reproduce the observed MZGR in local galaxies. Metal transport also contributes to setting metallicity profiles, but it is sensitive to the strength of radial gas flows in galaxies. While the mass loading of galactic winds is important to reproduce the mass metallicity relation (MZR), we find that metal mass loading is more important to reproducing the MZGR. Specifically, our model predicts preferential metal enrichment of galactic winds in low-mass galaxies. This conclusion is robust against our adopted scaling of the wind mass-loading factor, uncertainties in measured wind metallicities, and systematics due to metallicity calibrations. Overall, we find that at $z \sim 0$, galactic winds and metal transport are more important in setting metallicity gradients in low-mass galaxies whereas star formation feedback and gas accretion dominate setting metallicity gradients in massive galaxies. △ Less

Submitted 11 January, 2024; v1 submitted 28 March, 2023; originally announced March 2023.

Comments: 21 pages (+ appendix), 11 figures, 4 tables. Accepted by MNRAS

Closed-Loop Koopman Operator Approximation

Authors: Steven Dahdah, James Richard Forbes

Abstract: This paper proposes a method to identify a Koopman model of a feedback-controlled system given a known controller. The Koopman operator allows a nonlinear system to be rewritten as an infinite-dimensional linear system by viewing it in terms of an infinite set of lifting functions. A finite-dimensional approximation of the Koopman operator can be identified from data by choosing a finite subset of… ▽ More This paper proposes a method to identify a Koopman model of a feedback-controlled system given a known controller. The Koopman operator allows a nonlinear system to be rewritten as an infinite-dimensional linear system by viewing it in terms of an infinite set of lifting functions. A finite-dimensional approximation of the Koopman operator can be identified from data by choosing a finite subset of lifting functions and solving a regression problem in the lifted space. Existing methods are designed to identify open-loop systems. However, it is impractical or impossible to run experiments on some systems, such as unstable systems, in an open-loop fashion. The proposed method leverages the linearity of the Koopman operator, along with knowledge of the controller and the structure of the closed-loop system, to simultaneously identify the closed-loop and plant systems. The advantages of the proposed closed-loop Koopman operator approximation method are demonstrated in simulation using a Duffing oscillator and experimentally using a rotary inverted pendulum system. An open-source software implementation of the proposed method is publicly available, along with the experimental dataset generated for this paper. △ Less

Submitted 1 May, 2024; v1 submitted 27 March, 2023; originally announced March 2023.

Comments: 13 pages, 11 figures, 3 tables, accepted for accepted for publication in Machine Learning: Science and Technology

Journal ref: Mach. Learn.: Sci. Technol. 5 025038 (2024)

Characterizing the Conditional Galaxy Property Distribution using Gaussian Mixture Models

Authors: Yucheng Zhang, Anthony R. Pullen, Rachel S. Somerville, Patrick C. Breysse, John C. Forbes, Shengqi Yang, Yin Li, Abhishek S. Maniyar

Abstract: Line-intensity mapping (LIM) is a promising technique to constrain the global distribution of galaxy properties. To combine LIM experiments probing different tracers with traditional galaxy surveys and fully exploit the scientific potential of these observations, it is necessary to have a physically motivated modeling framework. As part of developing such a framework, in this work we introduce and… ▽ More Line-intensity mapping (LIM) is a promising technique to constrain the global distribution of galaxy properties. To combine LIM experiments probing different tracers with traditional galaxy surveys and fully exploit the scientific potential of these observations, it is necessary to have a physically motivated modeling framework. As part of developing such a framework, in this work we introduce and model the conditional galaxy property distribution (CGPD), i.e. the distribution of galaxy properties conditioned on the host halo mass and redshift. We consider five galaxy properties, including the galaxy stellar mass, molecular gas mass, galaxy radius, gas phase metallicity and star formation rate (SFR), which are important for predicting the emission lines of interest. The CGPD represents the full distribution of galaxies in the five dimensional property space; many important galaxy distribution functions and scaling relations, such as the stellar mass function and SFR main sequence, can be derived from integrating and projecting it. We utilize two different kinds of cosmological galaxy simulations, a semi-analytic model and the IllustrisTNG hydrodynamic simulation, to characterize the CGPD and explore how well it can be represented using a Gaussian mixture model (GMM). We find that with just a few ($\sim 3$) Gaussian components, a GMM can describe the CGPD of the simulated galaxies to high accuracy for both simulations. The CGPD can be mapped to LIM or other observables by constructing the appropriate relationship between galaxy properties and the relevant observable tracers. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 17 pages, 10 figures

The star formation history of the Sco-Cen association. Coherent star formation patterns in space and time

Authors: Sebastian Ratzenböck, Josefa E. Großschedl, João Alves, Núria Miret-Roig, Immanuel Bomze, John Forbes, Alyssa Goodman, Alvaro Hacar, Doug Lin, Stefan Meingast, Torsten Möller, Martin Piecka, Laura Posch, Alena Rottensteiner, Cameren Swiggum, Catherine Zucker

Abstract: We reconstruct the star formation history of the Sco-Cen OB association using a novel high-resolution age map of the region. We develop an approach to produce robust ages for Sco-Cen's recently identified 37 stellar clusters using the SigMA algorithm. The Sco-Cen star formation timeline reveals four periods of enhanced star formation activity, or bursts, remarkably separated by about 5 Myr. Of the… ▽ More We reconstruct the star formation history of the Sco-Cen OB association using a novel high-resolution age map of the region. We develop an approach to produce robust ages for Sco-Cen's recently identified 37 stellar clusters using the SigMA algorithm. The Sco-Cen star formation timeline reveals four periods of enhanced star formation activity, or bursts, remarkably separated by about 5 Myr. Of these, the second burst, which occurred about 15 million years ago, is by far the dominant, and most of Sco-Cen's stars and clusters were in place by the end of this burst. The formation of stars and clusters in Sco-Cen is correlated but not linearly, implying that more stars were formed per cluster during the peak of the star formation rate. Most of the clusters that are large enough to have supernova precursors were formed during the 15 Myr period. Star and cluster formation activity has been continuously declining since then. We have clear evidence that Sco-Cen formed from the inside out and contains 100-pc long chains of contiguous clusters exhibiting well-defined age gradients, from massive older clusters to smaller young clusters. These observables suggest an important role for feedback in forming about half of Sco-Cen stars, although follow-up work is needed to quantify this statement. Finally, we confirm that the Upper-Sco age controversy discussed in the literature during the last decades is solved: the nine clusters previously lumped together as Upper-Sco, a benchmark region for planet formation studies, exhibit a wide range of ages from 3 to 19 Myr. △ Less

Submitted 14 May, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: 20 pages, 15 figures, 3 tables, submitted to A&A

Journal ref: A&A 678, A71 (2023)

How Low Can Q Go?

Authors: John C. Forbes

Abstract: Gravitational instability plays a substantial role in the evolution of galaxies. Various schemes to include it in galaxy evolution models exist, generally assuming that the Toomre $Q$ parameter is self-regulated to $Q_\mathrm{crit}$, the critical $Q$ dividing stable from unstable conditions in a linear stability analysis. This assumption is in tension with observational estimates of $Q$ that find… ▽ More Gravitational instability plays a substantial role in the evolution of galaxies. Various schemes to include it in galaxy evolution models exist, generally assuming that the Toomre $Q$ parameter is self-regulated to $Q_\mathrm{crit}$, the critical $Q$ dividing stable from unstable conditions in a linear stability analysis. This assumption is in tension with observational estimates of $Q$ that find values far below any plausible value of $Q_\mathrm{crit}$. While the observations are subject to some uncertainty, this tension can more easily be relieved on the theoretical side by relaxing the common assumption that $Q\ge Q_\mathrm{crit}$. Based on observations of both $z\sim 2$ disks and local face-on galaxies, we estimate the effect of gravitational instability necessary to balance out every other physical process that affects $Q$. In particular we find that the disk's response to low $Q$ values can be described by simple functions that depend only on $Q$. These response functions allow galaxies to maintain $Q$ values below $Q_\mathrm{crit}$ in equilibrium over a wide range of parameters. Extremely low values of $Q$ are predicted when the gas surface density is greater than $\sim 10^3$ M$_\odot$ pc$^{-2}$, the rotation curve provides minimal shear, the orbital time becomes long, and/or when the gas is much more unstable than the stellar component. We suggest that these response functions should be used in place of the $Q\ge Q_\mathrm{crit}$ ansatz. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: Submitted to AAS Journals

The Neon Gap: Probing Ionization with Dwarf Galaxies at z~1

Authors: John Pharo, Yicheng Guo, David C. Koo, John C. Forbes, Puragra Guhathakurta

Abstract: We present measurements of [NeIII]λ3869 emission in z~1 low-mass galaxies taken from the Keck/DEIMOS spectroscopic surveys HALO7D and DEEPWinds. We identify 167 individual galaxies with significant [NeIII] emission lines, including 112 "dwarf" galaxies with log(M_{\star}/M_{\odot}) < 9.5, with 0.3 < z < 1.4. We also measure [NeIII] emission from composite spectra derived from all [OII]λλ3727,3729… ▽ More We present measurements of [NeIII]λ3869 emission in z~1 low-mass galaxies taken from the Keck/DEIMOS spectroscopic surveys HALO7D and DEEPWinds. We identify 167 individual galaxies with significant [NeIII] emission lines, including 112 "dwarf" galaxies with log(M_{\star}/M_{\odot}) < 9.5, with 0.3 < z < 1.4. We also measure [NeIII] emission from composite spectra derived from all [OII]λλ3727,3729 line emitters in this range. This provides a unique sample of [NeIII]-emitters in the gap between well-studied emitters at z = 0 and 2 < z < 3. To study evolution in ionization conditions in the ISM over this time, we analyze the log([NeIII]λ3869/[OII]λλ3727,3729) ratio (Ne3O2) as a function of the stellar mass and of the log([OIII]λλ4959,5007/[OII]λλ3727,3729) ratio (O32). We find that the typical star-forming dwarf galaxy at this redshift, as measured from the composite spectra, shares the Ne3O2-M_{\star} relation with local galaxies, but have higher O32 at given Ne3O2. This finding implies that the ionization and metallicity characteristics of the z~1 dwarf population do not evolve substantially from z~1 to z=0, suggesting that the known evolution in those parameter from z~2 has largely taken place by z~1. Individual [NeIII]-detected galaxies have emission characteristics situated between local and z~2 galaxies, with elevated Ne3O2 and O32 emission potentially explained by variations in stellar and nebular metallicity. We also compare our dwarf sample to similarly low-mass z > 7 galaxies identified in JWST Early Release Observations, finding four HALO7D dwarfs with similar size, metallicity, and star formation properties. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: Accepted to ApJL, 16.01.2023

Arkenstone I: A Novel Method for Robustly Capturing High Specific Energy Outflows In Cosmological Simulations

Authors: Matthew C. Smith, Drummond B. Fielding, Greg L. Bryan, Chang-Goo Kim, Eve C. Ostriker, Rachel S. Somerville, Jonathan Stern, Kung-Yi Su, Rainer Weinberger, Chia-Yu Hu, John C. Forbes, Lars Hernquist, Blakesley Burkhart, Yuan Li

Abstract: Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arep… ▽ More Arkenstone is a new model for multiphase, stellar feedback driven galactic winds designed for inclusion in coarse resolution cosmological simulations. In this first paper of a series, we describe the features that allow Arkenstone to properly treat high specific energy wind components and demonstrate them using idealised non-cosmological simulations of a galaxy with a realistic CGM, using the Arepo code. Hot, fast gas phases with low mass loadings are predicted to dominate the energy content of multiphase outflows. In order to treat the huge dynamic range of spatial scales involved in cosmological galaxy formation at feasible computational expense, cosmological volume simulations typically employ a Lagrangian code or else use adaptive mesh refinement with a quasi-Lagrangian refinement strategy. However, it is difficult to inject a high specific energy wind in a Lagrangian scheme without incurring artificial burstiness. Additionally, the low densities inherent to this type of flow result in poor spatial resolution. Arkenstone addresses these issues with a novel scheme for coupling energy into the ISM/CGM transition region which also provides the necessary level of refinement at the base of the wind. In the absence of our improvements, we show that poor spatial resolution near the sonic point of a hot, fast outflow leads to an underestimation of gas acceleration as the wind propagates. We explore the different mechanisms by which low and high specific energy winds can regulate the SFR of galaxies. In future work, we will demonstrate other aspects of the Arkenstone model. △ Less

Submitted 7 March, 2024; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: Published MNRAS, 27 pages, 17 figures; updated for consistency with journal version

Journal ref: MNRAS, 527, 1216 (2024)

Improving Self-Consistency in Underwater Mapping Through Laser-Based Loop Closure (Extended)

Authors: Thomas Hitchcox, James Richard Forbes

Abstract: Accurate, self-consistent bathymetric maps are needed to monitor changes in subsea environments and infrastructure. These maps are increasingly collected by underwater vehicles, and mapping requires an accurate vehicle navigation solution. Commercial off-the-shelf (COTS) navigation solutions for underwater vehicles often rely on external acoustic sensors for localization, however survey-grade acou… ▽ More Accurate, self-consistent bathymetric maps are needed to monitor changes in subsea environments and infrastructure. These maps are increasingly collected by underwater vehicles, and mapping requires an accurate vehicle navigation solution. Commercial off-the-shelf (COTS) navigation solutions for underwater vehicles often rely on external acoustic sensors for localization, however survey-grade acoustic sensors are expensive to deploy and limit the range of the vehicle. Techniques from the field of simultaneous localization and mapping, particularly loop closures, can improve the quality of the navigation solution over dead-reckoning, but are difficult to integrate into COTS navigation systems. This work presents a method to improve the self-consistency of bathymetric maps by smoothly integrating loop-closure measurements into the state estimate produced by a commercial subsea navigation system. Integration is done using a white-noise-on-acceleration motion prior, without access to raw sensor measurements or proprietary models. Improvements in map self-consistency are shown for both simulated and experimental datasets, including a 3D scan of an underwater shipwreck in Wiarton, Ontario, Canada. △ Less

Submitted 9 June, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 26 pages, 18 figures. V2 correct Table III x2 parameter values, Table VIII 'INS' values, and equation A.28

Journal ref: IEEE Transactions on Robotics, vol. 39, no. 3, pp. 1873-1892, 2023

The Earliest Stage of Galactic Star Formation

Authors: Charles L. Steinhardt, Vadim Rusakov, Thomas H. Clark, Andrei Diaconu, Conor McPartland, John Forbes, Albert Sneppen, John Weaver

Abstract: Using a recently-developed technique to estimate gas temperatures ($T_\textrm{SF}$) in star-forming regions from large photometric surveys, we propose a diagram, analogous to the Hertzsprung-Russell diagram for individual stars, to probe the evolution of individual galaxies. On this $T_\textrm{SF}$-sSFR (specific star formation rate) diagram, a small fraction of star-forming galaxies appear to be… ▽ More Using a recently-developed technique to estimate gas temperatures ($T_\textrm{SF}$) in star-forming regions from large photometric surveys, we propose a diagram, analogous to the Hertzsprung-Russell diagram for individual stars, to probe the evolution of individual galaxies. On this $T_\textrm{SF}$-sSFR (specific star formation rate) diagram, a small fraction of star-forming galaxies appear to be dominated by different feedback mechanisms than typical star-forming galaxies. These galaxies generically have younger stellar populations, lower stellar masses and increase in relative abundance towards higher redshifts, so we argue that these objects are in an earlier stage of galactic star formation. Further, Hubble observations find that these "core-forming" galaxies also exhibit distinct morphology, and that tracks on the $T_\textrm{SF}$-sSFR diagram are also a morphological sequence. Thus, unlike starburst phases which can be triggered environmentally, these earliest, core-forming galaxies, appear to be a stage that typical galaxies go through early in their star formation history. We therefore argue that most galaxies first go through a core formation stage, then subsequently disk formation, and finally become quiescent. △ Less

Submitted 22 June, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: ApJL 949, L38

Know What You Don't Know: Consistency in Sliding Window Filtering with Unobservable States Applied to Visual-Inertial SLAM (Extended Version)

Authors: Daniil Lisus, Mitchell Cohen, James Richard Forbes

Abstract: Estimation algorithms, such as the sliding window filter, produce an estimate and uncertainty of desired states. This task becomes challenging when the problem involves unobservable states. In these situations, it is critical for the algorithm to ``know what it doesn't know'', meaning that it must maintain the unobservable states as unobservable during algorithm deployment. This letter presents ge… ▽ More Estimation algorithms, such as the sliding window filter, produce an estimate and uncertainty of desired states. This task becomes challenging when the problem involves unobservable states. In these situations, it is critical for the algorithm to ``know what it doesn't know'', meaning that it must maintain the unobservable states as unobservable during algorithm deployment. This letter presents general requirements for maintaining consistency in sliding window filters involving unobservable states. The value of these requirements for designing navigation solutions is experimentally shown within the context of visual-inertial SLAM making use of IMU preintegration. △ Less

Submitted 8 May, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: Main paper accepted to Robotics and Automation Letters. Main paper has 8 pages, 3 figures. Supplemental materials are 6 pages, 0 figures after the main paper

Journal ref: IEEE Robotics and Automation Letters, vol. 8, no. 6, pp. 3382-3389, June 2023

A 3D View of Orion: I. Barnard's Loop

Authors: Michael M. Foley, Alyssa Goodman, Catherine Zucker, John C. Forbes, Ralf Konietzka, Cameren Swiggum, João Alves, John Bally, Juan D. Soler, Josefa E. Großschedl, Shmuel Bialy, Michael Y. Grudić, Reimar Leike, Torsten Ensslin

Abstract: Barnard's Loop is a famous arc of H$α$ emission located in the Orion star-forming region. Here, we provide evidence of a possible formation mechanism for Barnard's Loop and compare our results with recent work suggesting a major feedback event occurred in the region around 6 Myr ago. We present a 3D model of the large-scale Orion region, indicating coherent, radial, 3D expansion of the OBP-Near/Br… ▽ More Barnard's Loop is a famous arc of H$α$ emission located in the Orion star-forming region. Here, we provide evidence of a possible formation mechanism for Barnard's Loop and compare our results with recent work suggesting a major feedback event occurred in the region around 6 Myr ago. We present a 3D model of the large-scale Orion region, indicating coherent, radial, 3D expansion of the OBP-Near/Briceño-1 (OBP-B1) cluster in the middle of a large dust cavity. The large-scale gas in the region also appears to be expanding from a central point, originally proposed to be Orion X. OBP-B1 appears to serve as another possible center, and we evaluate whether Orion X or OBP-B1 is more likely to be the cause of the expansion. We find that neither cluster served as the single expansion center, but rather a combination of feedback from both likely propelled the expansion. Recent 3D dust maps are used to characterize the 3D topology of the entire region, which shows Barnard's Loop's correspondence with a large dust cavity around the OPB-B1 cluster. The molecular clouds Orion A, Orion B, and Orion $λ$ reside on the shell of this cavity. Simple estimates of gravitational effects from both stars and gas indicate that the expansion of this asymmetric cavity likely induced anisotropy in the kinematics of OBP-B1. We conclude that feedback from OBP-B1 has affected the structure of the Orion A, Orion B, and Orion $λ$ molecular clouds and may have played a major role in the formation of Barnard's Loop. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: 25 pages, 10 figures. Submitted to ApJ

A unified model for the co-evolution of galaxies and their circumgalactic medium: the relative roles of turbulence and atomic cooling physics

Authors: Viraj Pandya, Drummond B. Fielding, Greg L. Bryan, Christopher Carr, Rachel S. Somerville, Jonathan Stern, Claude-Andre Faucher-Giguere, Zachary Hafen, Daniel Angles-Alcazar, John C. Forbes

Abstract: The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and… ▽ More The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (``thermalization'') from a cool, turbulence-supported phase to a virial-temperature, thermally-supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high-redshift in ultrafaint progenitors and at low redshift in classical $M_{\rm vir}\sim10^{11}M_{\odot}$ dwarfs, while Milky Way-mass halos undergo the transition at $z\approx0.5$. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multi-phase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions. △ Less

Submitted 27 August, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: Updated to reflect version accepted by ApJ; minor text edits to clarify limitations and caveats of model

Reducing Two-Way Ranging Variance by Signal-Timing Optimization

Authors: Mohammed Ayman Shalaby, Charles Champagne Cossette, James Richard Forbes, Jerome Le Ny

Abstract: Time-of-flight-based ranging among transceivers with different clocks requires protocols that accommodate varying rates of the clocks. Double-sided two-way ranging (DS-TWR) is widely adopted as a standard protocol due to its accuracy; however, the precision of DS-TWR has not been clearly addressed. In this paper, an analytical model of the variance of DS-TWR is derived as a function of the user-pr… ▽ More Time-of-flight-based ranging among transceivers with different clocks requires protocols that accommodate varying rates of the clocks. Double-sided two-way ranging (DS-TWR) is widely adopted as a standard protocol due to its accuracy; however, the precision of DS-TWR has not been clearly addressed. In this paper, an analytical model of the variance of DS-TWR is derived as a function of the user-programmed response delays, which is then compared to the Cramer-Rao Lower Bound (CRLB). This is then used to formulate an optimization problem over the response delays in order to maximize the information gained from range measurements. The derived analytical variance model and optimized protocol are validated experimentally with 2 ranging UWB transceivers, where 29 million range measurements are collected. △ Less

Submitted 7 February, 2024; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: 7 pages, 5 figures, accepted for publication at IEEE Transactions on Aerospace and Electronic Systems

Calibration and Uncertainty Characterization for Ultra-Wideband Two-Way-Ranging Measurements

Authors: Mohammed Ayman Shalaby, Charles Champagne Cossette, James Richard Forbes, Jerome Le Ny

Abstract: Ultra-Wideband (UWB) systems are becoming increasingly popular for indoor localization, where range measurements are obtained by measuring the time-of-flight of radio signals. However, the range measurements typically suffer from a systematic error or bias that must be corrected for high-accuracy localization. In this paper, a ranging protocol is proposed alongside a robust and scalable antenna-de… ▽ More Ultra-Wideband (UWB) systems are becoming increasingly popular for indoor localization, where range measurements are obtained by measuring the time-of-flight of radio signals. However, the range measurements typically suffer from a systematic error or bias that must be corrected for high-accuracy localization. In this paper, a ranging protocol is proposed alongside a robust and scalable antenna-delay calibration procedure to accurately and efficiently calibrate antenna delays for many UWB tags. Additionally, the bias and uncertainty of the measurements are modelled as a function of the received-signal power. The full calibration procedure is presented using experimental training data of 3 aerial robots fitted with 2 UWB tags each, and then evaluated on 2 test experiments. A localization problem is then formulated on the experimental test data, and the calibrated measurements and their modelled uncertainty are fed into an extended Kalman filter (EKF). The proposed calibration is shown to yield an average of 46% improvement in localization accuracy. Lastly, the paper is accompanied by an open-source UWB-calibration Python library, which can be found at https://github.com/decargroup/uwb_calibration. △ Less

Submitted 16 February, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 7 pages, 7 figures, to appear at the 2023 International Conference on Robotics and Automation (ICRA)

Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching

Authors: Thomas K. Waters, Colton Peterson, Razieh Emami, Xuejian Shen, Lars Hernquist, Randall Smith, Mark Vogelsberger, Charles Alcock, Grant Tremblay, Matthew Liska, John C. Forbes, Jorge Moreno

Abstract: We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner… ▽ More We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback. △ Less

Submitted 5 January, 2024; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 12 figures

Performance Evaluation of 3D Keypoint Detectors and Descriptors on Coloured Point Clouds in Subsea Environments

Authors: Kyungmin Jung, Thomas Hitchcox, James Richard Forbes

Abstract: The recent development of high-precision subsea optical scanners allows for 3D keypoint detectors and feature descriptors to be leveraged on point cloud scans from subsea environments. However, the literature lacks a comprehensive survey to identify the best combination of detectors and descriptors to be used in these challenging and novel environments. This paper aims to identify the best detecto… ▽ More The recent development of high-precision subsea optical scanners allows for 3D keypoint detectors and feature descriptors to be leveraged on point cloud scans from subsea environments. However, the literature lacks a comprehensive survey to identify the best combination of detectors and descriptors to be used in these challenging and novel environments. This paper aims to identify the best detector/descriptor pair using a challenging field dataset collected using a commercial underwater laser scanner. Furthermore, studies have shown that incorporating texture information to extend geometric features adds robustness to feature matching on synthetic datasets. This paper also proposes a novel method of fusing images with underwater laser scans to produce coloured point clouds, which are used to study the effectiveness of 6D point cloud descriptors. △ Less

Submitted 26 February, 2024; v1 submitted 26 September, 2022; originally announced September 2022.

Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods

Authors: Chia-Yu Hu, Matthew C. Smith, Romain Teyssier, Greg L. Bryan, Robbert Verbeke, Andrew Emerick, Rachel S. Somerville, Blakesley Burkhart, Yuan Li, John C. Forbes, Tjitske Starkenburg

Abstract: We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid mo… ▽ More We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density-temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger supernova-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: gas in Lagrangian codes collapses to much higher densities than in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed. △ Less

Submitted 3 May, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: accepted version by ApJ (including a new simulation in Appendix B suggested by the referee)

Stochastic Modelling of Star Formation Histories III. Constraints from Physically-Motivated Gaussian Processes

Authors: Kartheik G. Iyer, Joshua S. Speagle, Neven Caplar, John C. Forbes, Eric Gawiser, Joel Leja, Sandro Tacchella

Abstract: Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow,… ▽ More Galaxy formation and evolution involves a variety of effectively stochastic processes that operate over different timescales. The Extended Regulator model provides an analytic framework for the resulting variability (or `burstiness') in galaxy-wide star formation due to these processes. It does this by relating the variability in Fourier space to the effective timescales of stochastic gas inflow, equilibrium, and dynamical processes influencing GMC creation and destruction using the power spectral density (PSD) formalism. We use the connection between the PSD and auto-covariance function (ACF) for general stochastic processes to reformulate this model as an auto-covariance function, which we use to model variability in galaxy star formation histories (SFHs) using physically-motivated Gaussian Processes in log SFR space. Using stellar population synthesis models, we then explore how changes in model stochasticity can affect spectral signatures across galaxy populations with properties similar to the Milky Way and present-day dwarfs as well as at higher redshifts. We find that, even at fixed scatter, perturbations to the stochasticity model (changing timescales vs overall variability) leave unique spectral signatures across both idealized and more realistic galaxy populations. Distributions of spectral features including H$α$ and UV-based SFR indicators, H$δ$ and Ca-H,K absorption line strengths, D$_n$(4000) and broadband colors provide testable predictions for galaxy populations from present and upcoming surveys with Hubble, Webb \& Roman. The Gaussian process SFH framework provides a fast, flexible implementation of physical covariance models for the next generation of SED modeling tools. Code to reproduce our results can be found at https://github.com/kartheikiyer/GP-SFH △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 33 pages, 18 figures. Submitted to ApJ, comments are welcome! Code available at https://github.com/kartheikiyer/GP-SFH

System Identification and Two-Degree-of-Freedom Control of Nonlinear, Viscoelastic Tissues

Authors: Amanda Bianco, Raphael Zonis, Anne-Marie Lauzon, James Richard Forbes, Gijs Ijpma

Abstract: Objective: This paper presents a force control scheme for brief isotonic holds in an isometrically contracted muscle tissue, with minimal overshoot and settling time to measure its shortening velocity, a key parameter of muscle function. Methods: A two-degree-of-freedom control configuration, formed by a feedback controller and a feedforward controller, is explored. The feedback controller is a pr… ▽ More Objective: This paper presents a force control scheme for brief isotonic holds in an isometrically contracted muscle tissue, with minimal overshoot and settling time to measure its shortening velocity, a key parameter of muscle function. Methods: A two-degree-of-freedom control configuration, formed by a feedback controller and a feedforward controller, is explored. The feedback controller is a proportional-integral controller and the feedforward controller is designed using the inverse of a control-oriented model of muscle tissue. A generalized linear model and a nonlinear model of muscle tissue are explored using input-output data and system identification techniques. The force control scheme is tested on equine airway smooth muscle and its robustness confirmed with murine flexor digitorum brevis muscle. Results: Performance and repeatability of the force control scheme as well as the number of inputs and level of supervision required from the user were assessed with a series of experiments. The force control scheme was able to fulfill the stated control objectives in most cases, including the requirements for settling time and overshoot. Conclusion: The proposed control scheme is shown to enable automation of force control for characterizing muscle mechanics with minimal user input required. Significance: This paper leverages an inversion-based feedforward controller based on a nonlinear physiological model in a system identification context that is superior to classic linear system identification. The control scheme can be used as a steppingstone for generalized control of nonlinear, viscoelastic materials. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 9 pages, 11 figures, accepted for publication in IEEE Transactions on Biomedical Engineering

Relative Position Estimation in Multi-Agent Systems Using Attitude-Coupled Range Measurements

Authors: Mohammed Shalaby, Charles Champagne Cossette, James Richard Forbes, Jerome Le Ny

Abstract: The ability to accurately estimate the position of robotic agents relative to one another, in possibly GPS-denied environments, is crucial to execute collaborative tasks. Inter-agent range measurements are available at a low cost, due to technologies such as ultra-wideband radio. However, the task of three-dimensional relative position estimation using range measurements in multi-agent systems suf… ▽ More The ability to accurately estimate the position of robotic agents relative to one another, in possibly GPS-denied environments, is crucial to execute collaborative tasks. Inter-agent range measurements are available at a low cost, due to technologies such as ultra-wideband radio. However, the task of three-dimensional relative position estimation using range measurements in multi-agent systems suffers from unobservabilities. This letter presents a sufficient condition for the observability of the relative positions, and satisfies the condition using a simple framework with only range measurements, an accelerometer, a rate gyro, and a magnetometer. The framework has been tested in simulation and in experiments, where 40-50 cm positioning accuracy is achieved using inexpensive off-the-shelf hardware. △ Less

Submitted 16 July, 2022; originally announced July 2022.

Comments: 7 pages, 10 figures, published in IEEE Robotics and Automation Letters

Journal ref: in IEEE Robotics and Automation Letters, vol. 6, no. 3, pp. 4955 - 4961, July 2021

Mind the Gap: Norm-Aware Adaptive Robust Loss for Multivariate Least-Squares Problems

Authors: Thomas Hitchcox, James Richard Forbes

Abstract: Measurement outliers are unavoidable when solving real-world robot state estimation problems. A large family of robust loss functions (RLFs) exists to mitigate the effects of outliers, including newly developed adaptive methods that do not require parameter tuning. All of these methods assume that residuals follow a zero-mean Gaussian-like distribution. However, in multivariate problems the residu… ▽ More Measurement outliers are unavoidable when solving real-world robot state estimation problems. A large family of robust loss functions (RLFs) exists to mitigate the effects of outliers, including newly developed adaptive methods that do not require parameter tuning. All of these methods assume that residuals follow a zero-mean Gaussian-like distribution. However, in multivariate problems the residual is often defined as a norm, and norms follow a Chi-like distribution with a non-zero mode value. This produces a "mode gap" that impacts the convergence rate and accuracy of existing RLFs. The proposed approach, "Adaptive MB," accounts for this gap by first estimating the mode of the residuals using an adaptive Chi-like distribution. Applying an existing adaptive weighting scheme only to residuals greater than the mode leads to more robust performance and faster convergence times in two fundamental state estimation problems, point cloud alignment and pose averaging. △ Less

Submitted 10 June, 2023; v1 submitted 18 June, 2022; originally announced June 2022.

Comments: 8 pages, 4 figures. This paper has been accepted for publication in IEEE Robotics and Automation Letters. V2: Update weighting in (13), (28) and re-run results. Hypothesis, methodology, and general findings remain unchanged. Update Sec. II-A to reference IRLS, and update citation [11] accordingly. Include acknowledgement to Mitchell Cohen. V3: Update Section II-A to least-squares

Journal ref: IEEE Robotics and Automation Letters, vol. 7, no. 3, pp. 7116-7123, 2022

Model Predictive Control of a Tandem-Rotor Helicopter With a Nonuniformly Spaced Prediction Horizon

Authors: Faraaz Ahmed, Ludwik Sobiesiak, James Richard Forbes

Abstract: This letter considers model predictive control of a tandem-rotor helicopter. The error is formulated using the matrix Lie group $SE_2(3)$. A reference trajectory to a target is calculated using a quartic guidance law, leveraging the differentially flat properties of the system, and refined using a finite-horizon linear quadratic regulator. The nonlinear system is linearized about the reference tra… ▽ More This letter considers model predictive control of a tandem-rotor helicopter. The error is formulated using the matrix Lie group $SE_2(3)$. A reference trajectory to a target is calculated using a quartic guidance law, leveraging the differentially flat properties of the system, and refined using a finite-horizon linear quadratic regulator. The nonlinear system is linearized about the reference trajectory enabling the formulation of a quadratic program with control input, attitude keep-in zone, and attitude error constraints. A non-uniformly spaced prediction horizon is leveraged to capture the multi-timescale dynamics while keeping the problem size tractable. Monte-Carlo simulations demonstrate robustness of the proposed control structure to initial conditions, model uncertainty, and environmental disturbances. △ Less

Submitted 14 June, 2022; originally announced June 2022.

Comments: 6 pages, 5 figures

MSC Class: 9306 (Primary); 4906 (Secondary)

Journal ref: in IEEE Control Systems Letters, vol. 6, pp. 2828-2833, 2022

Optimal Multi-robot Formations for Relative Pose Estimation Using Range Measurements

Authors: Charles Champagne Cossette, Mohammed Ayman Shalaby, David Saussie, Jerome Le Ny, James Richard Forbes

Abstract: In multi-robot missions, relative position and attitude information between agents is valuable for a variety of tasks such as mapping, planning, and formation control. In this paper, the problem of estimating relative poses from a set of inter-agent range measurements is investigated. Specifically, it is shown that the estimation accuracy is highly dependent on the true relative poses themselves,… ▽ More In multi-robot missions, relative position and attitude information between agents is valuable for a variety of tasks such as mapping, planning, and formation control. In this paper, the problem of estimating relative poses from a set of inter-agent range measurements is investigated. Specifically, it is shown that the estimation accuracy is highly dependent on the true relative poses themselves, which prompts the desire to find multi-agent formations that provide the best estimation performance. By direct maximization of Fischer information, it is shown in simulation and experiment that large improvements in estimation accuracy can be obtained by optimizing the formation geometry of a team of robots. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: 7 pages, 8 figures, submitted to International Conference on Intelligent Robots and Systems

Empirical Dust Attenuation Model Leads to More Realistic UVJ Diagram for TNG100 Galaxies

Authors: Gautam Nagaraj, John C. Forbes, Joel Leja, Dan Foreman-Mackey, Christopher C. Hayward

Abstract: Dust attenuation varies substantially from galaxy to galaxy and as of yet cannot be reproduced from first principles in theoretical models. In Nagaraj et al. (2022), we developed the first Bayesian population model of dust attenuation as a function of stellar population properties and projected galaxy shape, built on spectral energy distribution (SED) fits of nearly 30,000 galaxies in the 3D-HST g… ▽ More Dust attenuation varies substantially from galaxy to galaxy and as of yet cannot be reproduced from first principles in theoretical models. In Nagaraj et al. (2022), we developed the first Bayesian population model of dust attenuation as a function of stellar population properties and projected galaxy shape, built on spectral energy distribution (SED) fits of nearly 30,000 galaxies in the 3D-HST grism survey with broadband photometric coverage from the rest-frame UV to IR. In this paper, we apply the model to galaxies from the large-volume cosmological simulation TNG100. We produce a UVJ diagram and compare it with one obtained in previous work by applying approximate radiative transfer to the simulated galaxies. We find that the UVJ diagram based on our empirical model is in better agreement with observations than the previous effort, especially in the number density of dusty star forming galaxies. We also construct the intrinsic dust-free UVJ diagram for TNG and 3D-HST galaxies at z ~ 1, finding qualitative agreement but residual differences at the 10-20% level. These differences can be largely attributed to the finding that TNG galaxies have, on average, 29% younger stellar populations and 0.28 dex higher metallicities than observed galaxies. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 12 pages, 5 figures, submitted to ApJ