Showing 1–50 of 142 results for author: Lloyd, J

Application of first- and second-order adjoint methods to glacial isostatic adjustment incorporating rotational feedbacks

Authors: Ziheng Yu, David Al-Attar, Frank Syvret, Andrew J. Lloyd

Abstract: This paper revisits and extends the adjoint theory for glacial isostatic adjustment (GIA) of Crawford et al. (2018). Rotational feedbacks are now incorporated, and the application of the second-order adjoint method is described for the first time. The first-order adjoint method provides an efficient means for computing sensitivity kernels for a chosen objective functional, while the second-order a… ▽ More This paper revisits and extends the adjoint theory for glacial isostatic adjustment (GIA) of Crawford et al. (2018). Rotational feedbacks are now incorporated, and the application of the second-order adjoint method is described for the first time. The first-order adjoint method provides an efficient means for computing sensitivity kernels for a chosen objective functional, while the second-order adjoint method provides second-derivative information in the form of Hessian kernels. These latter kernels are required by efficient Newton-type optimisation schemes and within methods for quantifying uncertainty for non-linear inverse problems. Most importantly, the entire theory has been reformulated so as to simplify its implementation by others within the GIA community. In particular, the rate-formulation for the GIA forward problem introduced by Crawford et al. (2018) has been replaced with the conventional equations for modelling GIA in laterally heterogeneous earth models. The implementation of the first- and second-order adjoint problems should be relatively easy within both existing and new GIA codes, with only the inclusions of more general force terms being required. △ Less

Submitted 24 August, 2024; originally announced August 2024.

Stochastic Continuation of Trajectories in the Circular Restricted Three-Body Problem via Differential Algebra

Authors: Giacomo Acciarini, Nicola Baresi, David J. B. Lloyd, Dario Izzo

Abstract: Numerical continuation techniques are powerful tools that have been extensively used to identify particular solutions of nonlinear dynamical systems and enable trajectory design in chaotic astrodynamics problems such as the Circular Restricted Three-Body Problem. However, the applicability of equilibrium points and periodic orbits may be questionable in real-world applications where the uncertaint… ▽ More Numerical continuation techniques are powerful tools that have been extensively used to identify particular solutions of nonlinear dynamical systems and enable trajectory design in chaotic astrodynamics problems such as the Circular Restricted Three-Body Problem. However, the applicability of equilibrium points and periodic orbits may be questionable in real-world applications where the uncertainties of the initial conditions of the spacecraft and dynamical parameters of the problem (e.g., mass ratio parameter) are taken into consideration. Due to uncertain parameters and initial conditions, the spacecraft might not follow the reference periodic orbit owing to growing uncertainties that cause the satellite to deviate from its nominal path. Hence, it is crucial to keep track of the probability of finding the spacecraft in a given region. Building on previous work, we extend numerical continuation to moments of the distribution (i.e., stochastic continuation) by directly continuing moments of the probability density function of the spacecraft state. Only assuming normality of the initial conditions, and leveraging moment-generating functions, Isserlis' theorem, and the algebra of truncated polynomials, we propagate the distribution of the spacecraft state at consecutive surface of section crossings while retaining a symbolic map of the final moments of the distribution that depend on the initial mean and covariance matrix only. The goal of the work is to offer a differential algebra-based general framework to continue 3D periodic orbits in the presence of uncertain dynamical systems. The proposed approach is compared against traditional Monte Carlo simulations to validate the uncertainty propagation approach and demonstrate the advantages of the proposed in terms of uncertainty propagation computational burden and access to higher-dimensional problems. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: Paper presented and published in the proceedings of the 2024 ISSFD Conference. (Darmstadt, Germany)

Tracking and forecasting oscillatory data streams using Koopman autoencoders and Kalman filtering

Authors: Stephen A Falconer, David J. B. Lloyd, Naratip Santitissadeekorn

Abstract: Data-driven modelling techniques provide a method for deriving models of dynamical systems directly from complicated data streams. However, tracking and forecasting such data streams poses a significant challenge to most methods, as they assume the underlying process and model does not change over time. In this paper, we apply one such data-driven method, the Koopman autoencoder (KAE), to high-dim… ▽ More Data-driven modelling techniques provide a method for deriving models of dynamical systems directly from complicated data streams. However, tracking and forecasting such data streams poses a significant challenge to most methods, as they assume the underlying process and model does not change over time. In this paper, we apply one such data-driven method, the Koopman autoencoder (KAE), to high-dimensional oscillatory data to generate a low-dimensional latent space and model, where the system's dynamics appear linear. This allows one to accurately track and forecast systems where the underlying model may change over time. States and the model in the reduced order latent space can then be efficiently updated as new data becomes available, using data assimilation techniques such as the ensemble Kalman filter (EnKF), in a technique we call the KAE EnKF. We demonstrate that this approach is able to effectively track and forecast time-varying, nonlinear dynamical systems in synthetic examples. We then apply the KAE EnKF to a video of a physical pendulum, and achieve a significant improvement over current state-of-the-art methods. By generating effective latent space reconstructions, we find that we are able to construct accurate short-term forecasts and efficient adaptations to externally forced changes to the pendulum's frequency. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 44 pages, 23 figures, 1 table

Localized Multi-Dimensional Patterns

Authors: Jason J. Bramburger, Dan J. Hill, David J. B. Lloyd

Abstract: Localized patterns are coherent structures embedded in a quiescent state and occur in both discrete and continuous media across a wide range of applications. While it is well-understood how domain covering patterns (for example stripes and hexagons) emerge from a pattern-forming/Turing instability, analyzing the emergence of their localized counterparts remains a significant challenge. There has b… ▽ More Localized patterns are coherent structures embedded in a quiescent state and occur in both discrete and continuous media across a wide range of applications. While it is well-understood how domain covering patterns (for example stripes and hexagons) emerge from a pattern-forming/Turing instability, analyzing the emergence of their localized counterparts remains a significant challenge. There has been considerable progress in studying localized patterns over the past few decades, often by employing innovative mathematical tools and techniques. In particular, the study of localized pattern formation has benefited greatly from numerical techniques; the continuing advancement in computational power has helped to both identify new types patterns and further our understanding of their behavior. We review recent advances regarding the complex behavior of localized patterns and the mathematical tools that have been developed to understand them, covering various topics from spatial dynamics, exponential asymptotics, and numerical methods. We observe that the mathematical understanding of localized patterns decreases as the spatial dimension increases, thus providing significant open problems that will form the basis for future investigations. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: Review paper, 76 pages, 51 figures

The James Webb Interferometer: Space-based interferometric detections of PDS 70 b and c at 4.8 $μ$m

Authors: Dori Blakely, Doug Johnstone, Gabriele Cugno, Anand Sivaramakrishnan, Peter Tuthill, Ruobing Dong, Benjamin J. S. Pope, Loïc Albert, Max Charles, Rachel A. Cooper, Matthew De Furio, Louis Desdoigts, René Doyon, Logan Francis, Alexandra Z. Greenbaum, David Lafrenière, James P. Lloyd, Michael R. Meyer, Laurent Pueyo, Shrishmoy Ray, Joel Sánchez-Bermúdez, Anthony Soulain, Deepashri Thatte, Thomas Vandal

Abstract: We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of… ▽ More We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of both PDS 70 b and c provide evidence for circumplanetary disk emission through fitting SED models to these new measurements and those found in the literature. We also newly detect emission within the disk gap at an SNR of $\sim$4, at a position angle of $207^{+11}_{-10}$ degrees, and an unconstrained separation within $\sim$200 mas. Follow-up observations will be needed to determine the nature of this emission. We place a 5$σ$ upper limit of $Δ$mag = 7.56 on the contrast of the candidate PDS 70 d at 4.8 $μ$m, which indicates that if the previously observed emission at shorter wavelengths is due to a planet, this putative planet has a different atmospheric composition than PDS 70 b or c. Finally, we place upper limits on emission from any additional planets in the disk gap. We find an azimuthally averaged 5$σ$ upper limit of $Δ$mag $\approx$ 7.5 at separations greater than 125 mas. These are the deepest limits to date within $\sim$250 mas at 4.8 $μ$m and the first space-based interferometric observations of this system. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

Quasiparticle cooling algorithms for quantum many-body state preparation

Authors: Jerome Lloyd, Alexios Michailidis, Xiao Mi, Vadim Smelyanskiy, Dmitry A. Abanin

Abstract: Probing correlated states of many-body systems is one of the central tasks for quantum simulators and processors. A promising approach to state preparation is to realize desired correlated states as steady states of engineered dissipative evolution. A recent experiment with a Google superconducting quantum processor [X. Mi et al., Science 383, 1332 (2024)] demonstrated a cooling algorithm utilizin… ▽ More Probing correlated states of many-body systems is one of the central tasks for quantum simulators and processors. A promising approach to state preparation is to realize desired correlated states as steady states of engineered dissipative evolution. A recent experiment with a Google superconducting quantum processor [X. Mi et al., Science 383, 1332 (2024)] demonstrated a cooling algorithm utilizing auxiliary degrees of freedom that are periodically reset to remove quasiparticles from the system, thereby driving it towards the ground state. We develop a kinetic theory framework to describe quasiparticle cooling dynamics, and employ it to compare the efficiency of different cooling algorithms. In particular, we introduce a protocol where coupling to auxiliaries is modulated in time to minimize heating processes, and demonstrate that it allows a high-fidelity preparation of ground states in different quantum phases. We verify the validity of the kinetic theory description by an extensive comparison with numerical simulations of a 1d transverse-field Ising model using a solvable model and tensor-network techniques. Further, the effect of noise, which limits efficiency of variational quantum algorithms in near-term quantum processors, can be naturally described within the kinetic theory. We investigate the steady state quasiparticle population as a function of noise strength, and establish maximum noise values for achieving high-fidelity ground states. This work establishes quasiparticle cooling algorithms as a practical, robust method for many-body state preparation on near-term quantum processors. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 25 pages, 14 figures

The Effect of Pulsar Geometry on the Observed Gamma-ray Spectrum of Millisecond Pulsars

Authors: Sheridan J. Lloyd, Paula M. Chadwick, Anthony M. Brown

Abstract: We analyse 13 yrs of $\textit{Fermi}$-LAT PASS 8 events from 127 gamma-ray emitting millisecond pulsars (MSPs) in the energy range 0.1$-$100 GeV and significantly detect 118 MSPs. We fit the stacked emission with a log parabola (LP) spectral model which we show is preferred to two previously published models. We consider the influence of pulsar properties and observer geometric effects on spectral… ▽ More We analyse 13 yrs of $\textit{Fermi}$-LAT PASS 8 events from 127 gamma-ray emitting millisecond pulsars (MSPs) in the energy range 0.1$-$100 GeV and significantly detect 118 MSPs. We fit the stacked emission with a log parabola (LP) spectral model which we show is preferred to two previously published models. We consider the influence of pulsar properties and observer geometric effects on spectral features by defining energy flux colours for both the individual MSPs, and our stacked model as a baseline. There is no correlation of colours with pulsar luminosity, $\dot{E}$, surface magnetic field or magnetic impact angle. We also find that pulsar geometry has little effect on the observed gamma-ray spectrum which is in tension with previous modelling of gamma-ray emission with respect to pulsar geometry. Our LP MSP model is applicable to problems where an ensemble of gamma-ray MSPs is considered, such as that of the Galactic centre excess or in the case of emission from globular clusters. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. Published by Oxford University Press on behalf of the Royal Astronomical Society

Radial amplitude equations for fully localised planar patterns

Authors: Dan J. Hill, David J. B. Lloyd

Abstract: Isolated patches of spatially oscillating pattern have been found to emerge near a pattern-forming instability in a wide variety of experiments and mathematical models. However, there is currently no mathematical theory to explain this emergence or characterise the structure of these patches. We provide a method for formally deriving radial amplitude equations to planar patterns via non-autonomous… ▽ More Isolated patches of spatially oscillating pattern have been found to emerge near a pattern-forming instability in a wide variety of experiments and mathematical models. However, there is currently no mathematical theory to explain this emergence or characterise the structure of these patches. We provide a method for formally deriving radial amplitude equations to planar patterns via non-autonomous multiple-scale analysis and convolutional sums of products of Bessel functions. Our novel approach introduces nonautonomous differential operators, which allow for the systematic manipulation of Bessel functions, as well as previously unseen identities involving infinite sums of Bessel functions. Solutions of the amplitude equations describe fully localised patterns with non-trivial angular dependence, where localisation occurs in a purely radial direction. Amplitude equations are derived for multiple examples of patterns with dihedral symmetry, including fully localised hexagons and quasipatterns with twelve-fold rotational symmetry. In particular, we show how to apply the asymptotic method to the Swift--Hohenberg equation and general reaction-diffusion systems. △ Less

Submitted 16 August, 2024; v1 submitted 5 March, 2024; originally announced March 2024.

Comments: 19 pages, 4 figures

Pose and shear-based tactile servoing

Authors: John Lloyd, Nathan F. Lepora

Abstract: Tactile servoing is an important technique because it enables robots to manipulate objects with precision and accuracy while adapting to changes in their environments in real-time. One approach for tactile servo control with high-resolution soft tactile sensors is to estimate the contact pose relative to an object surface using a convolutional neural network (CNN) for use as a feedback signal. In… ▽ More Tactile servoing is an important technique because it enables robots to manipulate objects with precision and accuracy while adapting to changes in their environments in real-time. One approach for tactile servo control with high-resolution soft tactile sensors is to estimate the contact pose relative to an object surface using a convolutional neural network (CNN) for use as a feedback signal. In this paper, we investigate how the surface pose estimation model can be extended to include shear, and utilize these combined pose-and-shear models to develop a tactile robotic system that can be programmed for diverse non-prehensile manipulation tasks, such as object tracking, surface following, single-arm object pushing and dual-arm object pushing. In doing this, two technical challenges had to be overcome. Firstly, the use of tactile data that includes shear-induced slippage can lead to error-prone estimates unsuitable for accurate control, and so we modified the CNN into a Gaussian-density neural network and used a discriminative Bayesian filter to improve the predictions with a state dynamics model that utilizes the robot kinematics. Secondly, to achieve smooth robot motion in 3D space while interacting with objects, we used SE(3) velocity-based servo control, which required re-deriving the Bayesian filter update equations using Lie group theory, as many standard assumptions do not hold for state variables defined on non-Euclidean manifolds. In future, we believe that pose and shear-based tactile servoing will enable many object manipulation tasks and the fully-dexterous utilization of multi-fingered tactile robot hands. Video: https://www.youtube.com/watch?v=xVs4hd34ek0 △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted in International Journal of Robotics Research (IJRR). 29 pages, 20 figures. Related technical report: arXiv:2306.08560. Video: https://www.youtube.com/watch?v=xVs4hd34ek0

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems V: Do Self-Consistent Atmospheric Models Represent JWST Spectra? A Showcase With VHS 1256 b

Authors: Simon Petrus, Niall Whiteford, Polychronis Patapis, Beth A. Biller, Andrew Skemer, Sasha Hinkley, Genaro Suárez, Anna Lueber, Paulina Palma-Bifani, Jordan M. Stone, Johanna M. Vos, Caroline V. Morley, Pascal Tremblin, Benjamin Charnay, Christiane Helling, Brittany E. Miles, Aarynn L. Carter, Jason J. Wang, Markus Janson, Eileen C. Gonzales, Ben Sutlieff, Kielan K. W. Hoch, Mickaël Bonnefoy, Gaël Chauvin, Olivier Absil , et al. (97 additional authors not shown)

Abstract: The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. W… ▽ More The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models. △ Less

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

Comments: 32 pages, 16 figures, 6 tables, 2 appendices

The ballistic to diffusive crossover in a weakly-interacting Fermi gas

Authors: Jerome Lloyd, Tibor Rakovszky, Frank Pollmann, Curt von Keyserlingk

Abstract: Charge and energy are expected to diffuse in interacting systems of fermions at finite temperatures, even in the absence of disorder, with the interactions inducing a crossover from the coherent and ballistic streaming of quasi-particles at early times, to incoherent diffusive behavior at late times. The relevant crossover timescales and the transport coefficients are both controlled by the streng… ▽ More Charge and energy are expected to diffuse in interacting systems of fermions at finite temperatures, even in the absence of disorder, with the interactions inducing a crossover from the coherent and ballistic streaming of quasi-particles at early times, to incoherent diffusive behavior at late times. The relevant crossover timescales and the transport coefficients are both controlled by the strength of interactions. In this work we develop a numerical method to simulate such systems at high temperatures, applicable in a wide range of interaction strengths, by adapting Dissipation-assisted Operator Evolution (DAOE) to fermions. Our fermion DAOE, which approximates the exact dynamics by systematically discarding information from high $n$-point functions, is tailored to capture non-interacting dynamics exactly, thus providing a good starting point for the weakly interacting problem. Applying our method to a microscopic model of weakly interacting fermions, we numerically demonstrate that the crossover from ballistic to diffusive transport happens at a time $t_D\sim1/Δ^{2}$ and that the diffusion constant similarly scales as $D \sim 1/Δ^2$, where $Δ$ is the interaction strength. We substantiate this scaling with a Fermi's golden rule calculation in the operator spreading picture, interpreting $t_D$ as the fermion-fermion scattering time and lifetime of the single-particle Green's function. △ Less

Submitted 24 October, 2023; originally announced October 2023.

The \textit{JWST} Early Release Science Program for Direct Observations of Exoplanetary Systems III: Aperture Masking Interferometric Observations of the star HIP\,65426 at $\boldsymbol{3.8\,\rm{μm}}$

Authors: Shrishmoy Ray, Steph Sallum, Sasha Hinkley, Anand Sivamarakrishnan, Rachel Cooper, Jens Kammerer, Alexandra Z. Greebaum, Deepashri Thatte, Cecilia Lazzoni, Andrei Tokovinin, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan , et al. (98 additional authors not shown)

Abstract: We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{JWST} Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of ${}0.5λ/D$ f… ▽ More We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{JWST} Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of ${}0.5λ/D$ for an interferometer), which are inaccessible with the classical inner working angles of the \textit{JWST} coronagraphs. When combined with \textit{JWST}'s unprecedented infrared sensitivity, this mode has the potential to probe a new portion of parameter space across a wide array of astronomical observations. Using this mode, we are able to achieve a contrast of $Δm_{F380M}{\sim }7.8$\,mag relative to the host star at a separation of ${\sim}0.07\arcsec$ but detect no additional companions interior to the known companion HIP\,65426\,b. Our observations thus rule out companions more massive than $10{-}12\,\rm{M\textsubscript{Jup}}$ at separations ${\sim}10{-}20\,\rm{au}$ from HIP\,65426, a region out of reach of ground or space-based coronagraphic imaging. These observations confirm that the AMI mode on \textit{JWST} is sensitive to planetary mass companions orbiting at the water frost line, even for more distant stars at $\sim$100\,pc. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening essentially unexplored parameter space. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 15 pages, 9 figures, submitted to ApJ Letters

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems IV: NIRISS Aperture Masking Interferometry Performance and Lessons Learned

Authors: Steph Sallum, Shrishmoy Ray, Jens Kammerer, Anand Sivaramakrishnan, Rachel Cooper, Alexandra Z. Greebaum, Deepashri Thatte, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Sasha Hinkley, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan, Sarah K. Betti, Anthony Boccaletti , et al. (98 additional authors not shown)

Abstract: We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early… ▽ More We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early Release Science (ERS) 1386 program with a deep search for close-in companions in the HIP 65426 exoplanetary system. As part of ERS 1386, we use the same data set to explore the random, static, and calibration errors of NIRISS AMI observables. We compare the observed noise properties and achievable contrast to theoretical predictions. We explore possible sources of calibration errors and show that differences in charge migration between the observations of HIP 65426 and point-spread function calibration stars can account for the achieved contrast curves. Lastly, we use self-calibration tests to demonstrate that with adequate calibration NIRISS F380M AMI can reach contrast levels of $\sim9-10$ mag at $\gtrsim λ/D$. These tests lead us to observation planning recommendations and strongly motivate future studies aimed at producing sophisticated calibration strategies taking these systematic effects into account. This will unlock the unprecedented capabilities of JWST/NIRISS AMI, with sensitivity to significantly colder, lower-mass exoplanets than lower-contrast ground-based AMI setups, at orbital separations inaccessible to JWST coronagraphy. △ Less

Submitted 11 March, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

The Simons Observatory: A fully remote controlled calibration system with a sparse wire grid for cosmic microwave background telescopes

Authors: Masaaki Murata, Hironobu Nakata, Kengo Iijima, Shunsuke Adachi, Yudai Seino, Kenji Kiuchi, Frederick Matsuda, Michael J. Randall, Kam Arnold, Nicholas Galitzki, Bradley R. Johnson, Brian Keating, Akito Kusaka, John B. Lloyd, Joseph Seibert, Maximiliano Silva-Feaver, Osamu Tajima, Tomoki Terasaki, Kyohei Yamada

Abstract: For cosmic microwave background (CMB) polarization observations, calibration of detector polarization angles is essential. We have developed a fully remote controlled calibration system with a sparse wire grid that reflects linearly polarized light along the wire direction. The new feature is a remote-controlled system for regular calibration, which has not been possible in sparse wire grid calibr… ▽ More For cosmic microwave background (CMB) polarization observations, calibration of detector polarization angles is essential. We have developed a fully remote controlled calibration system with a sparse wire grid that reflects linearly polarized light along the wire direction. The new feature is a remote-controlled system for regular calibration, which has not been possible in sparse wire grid calibrators in past experiments. The remote control can be achieved by two electric linear actuators that load or unload the sparse wire grid into a position centered on the optical axis of a telescope between the calibration time and CMB observation. Furthermore, the sparse wire grid can be rotated by a motor. A rotary encoder and a gravity sensor are installed on the sparse wire grid to monitor the wire direction. They allow us to achieve detector angle calibration with expected systematic error of $0.08^{\circ}$. The calibration system will be installed in small-aperture telescopes at Simons Observatory. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Journal ref: Rev. Sci. Instrum. 94, 124502 (2023)

Sim-to-Real Model-Based and Model-Free Deep Reinforcement Learning for Tactile Pushing

Authors: Max Yang, Yijiong Lin, Alex Church, John Lloyd, Dandan Zhang, David A. W. Barton, Nathan F. Lepora

Abstract: Object pushing presents a key non-prehensile manipulation problem that is illustrative of more complex robotic manipulation tasks. While deep reinforcement learning (RL) methods have demonstrated impressive learning capabilities using visual input, a lack of tactile sensing limits their capability for fine and reliable control during manipulation. Here we propose a deep RL approach to object pushi… ▽ More Object pushing presents a key non-prehensile manipulation problem that is illustrative of more complex robotic manipulation tasks. While deep reinforcement learning (RL) methods have demonstrated impressive learning capabilities using visual input, a lack of tactile sensing limits their capability for fine and reliable control during manipulation. Here we propose a deep RL approach to object pushing using tactile sensing without visual input, namely tactile pushing. We present a goal-conditioned formulation that allows both model-free and model-based RL to obtain accurate policies for pushing an object to a goal. To achieve real-world performance, we adopt a sim-to-real approach. Our results demonstrate that it is possible to train on a single object and a limited sample of goals to produce precise and reliable policies that can generalize to a variety of unseen objects and pushing scenarios without domain randomization. We experiment with the trained agents in harsh pushing conditions, and show that with significantly more training samples, a model-free policy can outperform a model-based planner, generating shorter and more reliable pushing trajectories despite large disturbances. The simplicity of our training environment and effective real-world performance highlights the value of rich tactile information for fine manipulation. Code and videos are available at https://sites.google.com/view/tactile-rl-pushing/. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: Accepted by IEEE Robotics and Automation Letters (RA-L)

Reciprocity and sensitivity kernels for sea level fingerprints

Authors: David Al-Attar, Frank Syvret, Ophelia Crawford, Jerry X. Mitrovica, Andrew J. Lloyd

Abstract: Reciprocity theorems are established for the elastic sea level fingerprint problem including rotational feedbacks. In their simplest form, these results show that the sea level change at a location x due to melting a unit point mass of ice at x' is equal to the sea level change at x' due to melting a unit point mass of ice at x. This identity holds irrespective of the shoreline geometry or of late… ▽ More Reciprocity theorems are established for the elastic sea level fingerprint problem including rotational feedbacks. In their simplest form, these results show that the sea level change at a location x due to melting a unit point mass of ice at x' is equal to the sea level change at x' due to melting a unit point mass of ice at x. This identity holds irrespective of the shoreline geometry or of lateral variations in elastic Earth structure. Using the reciprocity theorems, sensitivity kernels for sea level and related observables with respect to the ice load can be readily derived. It is notable that calculation of the sensitivity kernels is possible using standard fingerprint codes, though for some types of observable a slight generalisation to the fingerprint problem must be considered. These results are of use within coastal hazard assessment and have a range of potential applications within studies of modern-day sea level change. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: Paper submitted to Geophysical Journal International

Bi-Touch: Bimanual Tactile Manipulation with Sim-to-Real Deep Reinforcement Learning

Authors: Yijiong Lin, Alex Church, Max Yang, Haoran Li, John Lloyd, Dandan Zhang, Nathan F. Lepora

Abstract: Bimanual manipulation with tactile feedback will be key to human-level robot dexterity. However, this topic is less explored than single-arm settings, partly due to the availability of suitable hardware along with the complexity of designing effective controllers for tasks with relatively large state-action spaces. Here we introduce a dual-arm tactile robotic system (Bi-Touch) based on the Tactile… ▽ More Bimanual manipulation with tactile feedback will be key to human-level robot dexterity. However, this topic is less explored than single-arm settings, partly due to the availability of suitable hardware along with the complexity of designing effective controllers for tasks with relatively large state-action spaces. Here we introduce a dual-arm tactile robotic system (Bi-Touch) based on the Tactile Gym 2.0 setup that integrates two affordable industrial-level robot arms with low-cost high-resolution tactile sensors (TacTips). We present a suite of bimanual manipulation tasks tailored towards tactile feedback: bi-pushing, bi-reorienting and bi-gathering. To learn effective policies, we introduce appropriate reward functions for these tasks and propose a novel goal-update mechanism with deep reinforcement learning. We also apply these policies to real-world settings with a tactile sim-to-real approach. Our analysis highlights and addresses some challenges met during the sim-to-real application, e.g. the learned policy tended to squeeze an object in the bi-reorienting task due to the sim-to-real gap. Finally, we demonstrate the generalizability and robustness of this system by experimenting with different unseen objects with applied perturbations in the real world. Code and videos are available at https://sites.google.com/view/bi-touch/. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: Accepted by IEEE Robotics and Automation Letters (RA-L)

A pose and shear-based tactile robotic system for object tracking, surface following and object pushing

Authors: John Lloyd, Nathan Lepora

Abstract: Tactile perception is a crucial sensing modality in robotics, particularly in scenarios that require precise manipulation and safe interaction with other objects. Previous research in this area has focused extensively on tactile perception of contact poses as this is an important capability needed for tasks such as traversing an object's surface or edge, manipulating an object, or pushing an objec… ▽ More Tactile perception is a crucial sensing modality in robotics, particularly in scenarios that require precise manipulation and safe interaction with other objects. Previous research in this area has focused extensively on tactile perception of contact poses as this is an important capability needed for tasks such as traversing an object's surface or edge, manipulating an object, or pushing an object along a predetermined path. Another important capability needed for tasks such as object tracking and manipulation is estimation of post-contact shear but this has received much less attention. Indeed, post-contact shear has often been considered a "nuisance variable" and is removed if possible because it can have an adverse effect on other types of tactile perception such as contact pose estimation. This paper proposes a tactile robotic system that can simultaneously estimate both the contact pose and post-contact shear, and use this information to control its interaction with other objects. Moreover, our new system is capable of interacting with other objects in a smooth and continuous manner, unlike the stepwise, position-controlled systems we have used in the past. We demonstrate the capabilities of our new system using several different controller configurations, on tasks including object tracking, surface following, single-arm object pushing, and dual-arm object pushing. △ Less

Submitted 26 June, 2023; v1 submitted 14 June, 2023; originally announced June 2023.

Comments: A video demonstrating the methods described in this paper is available at https://www.youtube.com/watch?v=xVs4hd34ek0

Stable Quantum-Correlated Many Body States through Engineered Dissipation

Authors: X. Mi, A. A. Michailidis, S. Shabani, K. C. Miao, P. V. Klimov, J. Lloyd, E. Rosenberg, R. Acharya, I. Aleiner, T. I. Andersen, M. Ansmann, F. Arute, K. Arya, A. Asfaw, J. Atalaya, J. C. Bardin, A. Bengtsson, G. Bortoli, A. Bourassa, J. Bovaird, L. Brill, M. Broughton, B. B. Buckley, D. A. Buell, T. Burger , et al. (142 additional authors not shown)

Abstract: Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-… ▽ More Engineered dissipative reservoirs have the potential to steer many-body quantum systems toward correlated steady states useful for quantum simulation of high-temperature superconductivity or quantum magnetism. Using up to 49 superconducting qubits, we prepared low-energy states of the transverse-field Ising model through coupling to dissipative auxiliary qubits. In one dimension, we observed long-range quantum correlations and a ground-state fidelity of 0.86 for 18 qubits at the critical point. In two dimensions, we found mutual information that extends beyond nearest neighbors. Lastly, by coupling the system to auxiliaries emulating reservoirs with different chemical potentials, we explored transport in the quantum Heisenberg model. Our results establish engineered dissipation as a scalable alternative to unitary evolution for preparing entangled many-body states on noisy quantum processors. △ Less

Submitted 5 April, 2024; v1 submitted 26 April, 2023; originally announced April 2023.

Journal ref: Science 383, 1332-1337 (2024)

Tactile-Driven Gentle Grasping for Human-Robot Collaborative Tasks

Authors: Christopher J. Ford, Haoran Li, John Lloyd, Manuel G. Catalano, Matteo Bianchi, Efi Psomopoulou, Nathan F. Lepora

Abstract: This paper presents a control scheme for force sensitive, gentle grasping with a Pisa/IIT anthropomorphic SoftHand equipped with a miniaturised version of the TacTip optical tactile sensor on all five fingertips. The tactile sensors provide high-resolution information about a grasp and how the fingers interact with held objects. We first describe a series of hardware developments for performing as… ▽ More This paper presents a control scheme for force sensitive, gentle grasping with a Pisa/IIT anthropomorphic SoftHand equipped with a miniaturised version of the TacTip optical tactile sensor on all five fingertips. The tactile sensors provide high-resolution information about a grasp and how the fingers interact with held objects. We first describe a series of hardware developments for performing asynchronous sensor data acquisition and processing, resulting in a fast control loop sufficient for real-time grasp control. We then develop a novel grasp controller that uses tactile feedback from all five fingertip sensors simultaneously to gently and stably grasp 43 objects of varying geometry and stiffness, which is then applied to a human-to-robot handover task. These developments open the door to more advanced manipulation with underactuated hands via fast reflexive control using high-resolution tactile sensing. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: Manuscript accepted to ICRA 2023. 6+n pages, 7 figures

Hamiltonian Cycles on Ammann-Beenker Tilings

Authors: Shobhna Singh, Jerome Lloyd, Felix Flicker

Abstract: We provide a simple algorithm for constructing Hamiltonian graph cycles (visiting every vertex exactly once) on a set of arbitrarily large finite subgraphs of aperiodic two-dimensional Ammann-Beenker (AB) tilings. Using this result, and the discrete scale symmetry of AB tilings, we find exact solutions to a range of other problems which lie in the complexity class NP-complete for general graphs. T… ▽ More We provide a simple algorithm for constructing Hamiltonian graph cycles (visiting every vertex exactly once) on a set of arbitrarily large finite subgraphs of aperiodic two-dimensional Ammann-Beenker (AB) tilings. Using this result, and the discrete scale symmetry of AB tilings, we find exact solutions to a range of other problems which lie in the complexity class NP-complete for general graphs. These include the equal-weight traveling salesperson problem, providing, for example, the most efficient route a scanning tunneling microscope tip could take to image the atoms of physical quasicrystals with AB symmetries; the longest path problem, whose solution demonstrates that collections of flexible molecules of any length can adsorb onto AB quasicrystal surfaces at density one, with possible applications to catalysis; and the three-coloring problem, giving ground states for the $q$-state Potts model ($q \ge 3$) of magnetic interactions defined on the planar dual to AB, which may provide useful models for protein folding. △ Less

Submitted 30 June, 2024; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 25 pages, 16 figures

Analysing transitions from a Turing instability to large periodic patterns in a reaction-diffusion system

Authors: Christopher Brown, Gianne Derks, Peter van Heijster, David J. B. Lloyd

Abstract: Analytically tracking patterns emerging from a small amplitude Turing instability to large amplitude remains a challenge as no general theory exists. In this paper, we consider a three component reaction-diffusion system with one of its components singularly perturbed, this component is known as the fast variable. We develop an analytical theory describing the periodic patterns emerging from a Tur… ▽ More Analytically tracking patterns emerging from a small amplitude Turing instability to large amplitude remains a challenge as no general theory exists. In this paper, we consider a three component reaction-diffusion system with one of its components singularly perturbed, this component is known as the fast variable. We develop an analytical theory describing the periodic patterns emerging from a Turing instability using geometric singular perturbation theory. We show analytically that after the initial Turing instability, spatially periodic patterns evolve into a small amplitude spike in the fast variable whose amplitude grows as one moves away from onset. This is followed by a secondary transition where the spike in the fast variable widens, its periodic pattern develops two sharp transitions between two flat states and the amplitudes of the other variables grow. The final type of transition we uncover analytically is where the flat states of the fast variable develop structure in the periodic pattern. The analysis is illustrated and motivated by a numerical investigation. We conclude with a preliminary numerical investigation where we uncover more complicated periodic patterns and snaking-like behaviour that are driven by the three transitions analysed in this paper. This paper provides a crucial step towards understanding how periodic patterns transition from a Turing instability to large amplitude. △ Less

Submitted 3 November, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: 38 pages, 13 figures

MSC Class: 34D15; 34E15; 35K40; 35K57; 37J46

Combining Dynamic Mode Decomposition with Ensemble Kalman Filtering for Tracking and Forecasting

Authors: Stephen A Falconer, David J. B. Lloyd, Naratip Santitissadeekorn

Abstract: Data assimilation techniques, such as ensemble Kalman filtering, have been shown to be a highly effective and efficient way to combine noisy data with a mathematical model to track and forecast dynamical systems. However, when dealing with high-dimensional data, in many situations one does not have a model, so data assimilation techniques cannot be applied. In this paper, we use dynamic mode decom… ▽ More Data assimilation techniques, such as ensemble Kalman filtering, have been shown to be a highly effective and efficient way to combine noisy data with a mathematical model to track and forecast dynamical systems. However, when dealing with high-dimensional data, in many situations one does not have a model, so data assimilation techniques cannot be applied. In this paper, we use dynamic mode decomposition to generate a low-dimensional, linear model of a dynamical system directly from high-dimensional data, which is defined by temporal and spatial modes, that we can then use with data assimilation techniques such as the ensemble Kalman filter. We show how the dynamic mode decomposition can be combined with the ensemble Kalman filter (which we call the DMDEnKF) to iteratively update the current state and temporal modes as new data becomes available. We demonstrate that this approach is able to track time varying dynamical systems in synthetic examples, and experiment with the use of time-delay embeddings. We then apply the DMDEnKF to real world seasonal influenza-like illness data from the USA Centers for Disease Control and Prevention, and find that for short term forecasting, the DMDEnKF is comparable to the best mechanistic models in the ILINet competition. △ Less

Submitted 13 January, 2023; originally announced January 2023.

The Near Infrared Imager and Slitless Spectrograph for JWST -- V. Kernel Phase Imaging and Data Analysis

Authors: Jens Kammerer, Rachel A. Cooper, Thomas Vandal, Deepashri Thatte, Frantz Martinache, Anand Sivaramakrishnan, Alexander Chaushev, Tomas Stolker, James P. Lloyd, Loïc Albert, René Doyon, Steph Sallum, Marshall D. Perrin, Laurent Pueyo, Antoine Mérand, Alexandre Gallenne, Alexandra Greenbaum, Joel Sanchez-Bermudez, Dori Blakely, Doug Johnstone, Kevin Volk, Andre Martel, Paul Goudfrooij, Michael R. Meyer, Chris J. Willott , et al. (4 additional authors not shown)

Abstract: Kernel phase imaging (KPI) enables the direct detection of substellar companions and circumstellar dust close to and below the classical (Rayleigh) diffraction limit. We present a kernel phase analysis of JWST NIRISS full pupil images taken during the instrument commissioning and compare the performance to closely related NIRISS aperture masking interferometry (AMI) observations. For this purpose,… ▽ More Kernel phase imaging (KPI) enables the direct detection of substellar companions and circumstellar dust close to and below the classical (Rayleigh) diffraction limit. We present a kernel phase analysis of JWST NIRISS full pupil images taken during the instrument commissioning and compare the performance to closely related NIRISS aperture masking interferometry (AMI) observations. For this purpose, we develop and make publicly available the custom "Kpi3Pipeline" enabling the extraction of kernel phase observables from JWST images. The extracted observables are saved into a new and versatile kernel phase FITS file (KPFITS) data exchange format. Furthermore, we present our new and publicly available "fouriever" toolkit which can be used to search for companions and derive detection limits from KPI, AMI, and long-baseline interferometry observations while accounting for correlated uncertainties in the model fitting process. Among the four KPI targets that were observed during NIRISS instrument commissioning, we discover a low-contrast (~1:5) close-in (~1 $λ/D$) companion candidate around CPD-66~562 and a new high-contrast (~1:170) detection separated by ~1.5 $λ/D$ from 2MASS~J062802.01-663738.0. The 5-$σ$ companion detection limits around the other two targets reach ~6.5 mag at ~200 mas and ~7 mag at ~400 mas. Comparing these limits to those obtained from the NIRISS AMI commissioning observations, we find that KPI and AMI perform similar in the same amount of observing time. Due to its 5.6 times higher throughput if compared to AMI, KPI is beneficial for observing faint targets and superior to AMI at separations >325 mas. At very small separations (<100 mas) and between ~250-325 mas, AMI slightly outperforms KPI which suffers from increased photon noise from the core and the first Airy ring of the point-spread function. △ Less

Submitted 3 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 34 pages, 17 figures, accepted for publication in PASP

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- IV. Aperture Masking Interferometry

Authors: Anand Sivaramakrishnan, Peter Tuthill, James P. Lloyd, Alexandra Z. Greenbaum, Deepashri Thatte, Rachel A. Cooper, Thomas Vandal, Jens Kammerer, Joel Sanchez-Bermudez, Benjamin J. S. Pope, Dori Blakely, Loïc Albert, Neil J. Cook, Doug Johnstone, André R. Martel, Kevin Volk, Anthony Soulain, Étienne Artigau, David Lafrenière, Chris J. Willott, Sébastien Parmentier, K. E. Saavik Ford, Barry McKernan, M. Begoña Vila, Neil Rowlands , et al. (14 additional authors not shown)

Abstract: The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, pres… ▽ More The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, present some sample science cases, explain its operational observing strategies, indicate how AMI proposals can be developed with data simulations, and how AMI data can be analyzed. We also present key results from commissioning AMI. Since the allied Kernel Phase Imaging (KPI) technique benefits from AMI operational strategies, we also cover NIRISS KPI methods and analysis techniques, including a new user-friendly KPI pipeline. The NIRISS KPI bright limit is $\simeq 8$ W2 magnitudes. AMI (and KPI) achieve an inner working angle of $\sim 70$ mas that is well inside the $\sim 400$ mas NIRCam inner working angle for its circular occulter coronagraphs at comparable wavelengths. △ Less

Submitted 7 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 30 pages, 10 figures

Dihedral rings of patterns emerging from a Turing bifurcation

Authors: Dan J. Hill, Jason J. Bramburger, David J. B. Lloyd

Abstract: Collective organisation of patterns into ring-like configurations has been well-studied when patterns are subject to either weak or semi-strong interactions. However, little is known numerically or analytically about their formation when the patterns are strongly interacting. We prove that approximate strongly interacting patterns can emerge in various ring-like dihedral configurations, bifurcatin… ▽ More Collective organisation of patterns into ring-like configurations has been well-studied when patterns are subject to either weak or semi-strong interactions. However, little is known numerically or analytically about their formation when the patterns are strongly interacting. We prove that approximate strongly interacting patterns can emerge in various ring-like dihedral configurations, bifurcating from quiescence near a Turing instability in generic two-component reaction-diffusion systems. The methods used are constructive and provide accurate initial conditions for numerical continuation methods to path-follow these ring-like patterns in parameter space. Our analysis is complemented by numerical investigations that illustrate our findings. △ Less

Submitted 7 February, 2024; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 35 pages, 11 figures

Journal ref: Nonlinearity 37, 035015 (2024)

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems II: A 1 to 20 Micron Spectrum of the Planetary-Mass Companion VHS 1256-1257 b

Authors: Brittany E. Miles, Beth A. Biller, Polychronis Patapis, Kadin Worthen, Emily Rickman, Kielan K. W. Hoch, Andrew Skemer, Marshall D. Perrin, Niall Whiteford, Christine H. Chen, B. Sargent, Sagnick Mukherjee, Caroline V. Morley, Sarah E. Moran, Mickael Bonnefoy, Simon Petrus, Aarynn L. Carter, Elodie Choquet, Sasha Hinkley, Kimberly Ward-Duong, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Laurent Pueyo, Shrishmoy Ray, Karl R. Stapelfeldt , et al. (79 additional authors not shown)

Abstract: We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude… ▽ More We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude diagram where substellar atmospheres transition from cloudy to clear. We observed VHS 1256~b with \textit{JWST}'s NIRSpec IFU and MIRI MRS modes for coverage from 1 $μ$m to 20 $μ$m at resolutions of $\sim$1,000 - 3,700. Water, methane, carbon monoxide, carbon dioxide, sodium, and potassium are observed in several portions of the \textit{JWST} spectrum based on comparisons from template brown dwarf spectra, molecular opacities, and atmospheric models. The spectral shape of VHS 1256 b is influenced by disequilibrium chemistry and clouds. We directly detect silicate clouds, the first such detection reported for a planetary-mass companion. △ Less

Submitted 4 July, 2024; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: Accepted ApJL. Iterations of spectra reduced by the ERS team are hosted at this link: https://github.com/bemiles/JWST_VHS1256b_Reduction/tree/main/reduced_spectra

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems I: High Contrast Imaging of the Exoplanet HIP 65426 b from 2-16 $μ$m

Authors: Aarynn L. Carter, Sasha Hinkley, Jens Kammerer, Andrew Skemer, Beth A. Biller, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Simon Petrus, Jordan M. Stone, Kimberly Ward-Duong, Jason J. Wang, Julien H. Girard, Dean C. Hines, Marshall D. Perrin, Laurent Pueyo, William O. Balmer, Mariangela Bonavita, Mickael Bonnefoy, Gael Chauvin, Elodie Choquet, Valentin Christiaens, Camilla Danielski, Grant M. Kennedy, Elisabeth C. Matthews, Brittany E. Miles , et al. (86 additional authors not shown)

Abstract: We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exo… ▽ More We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exoplanet to be obtained by JWST, and the first ever direct detection of an exoplanet beyond 5 $μ$m. These observations demonstrate that JWST is exceeding its nominal predicted performance by up to a factor of 10, depending on separation and subtraction method, with measured 5$σ$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $μ$m and MIRI at 11.3 $μ$m, respectively. These contrast limits provide sensitivity to sub-Jupiter companions with masses as low as 0.3$M_\mathrm{Jup}$ beyond separations of $\sim$100 au. Together with existing ground-based near-infrared data, the JWST photometry are well fit by a BT-SETTL atmospheric model from 1-16 $μ$m, and span $\sim$97% of HIP 65426 b's luminous range. Independent of the choice of model atmosphere we measure an empirical bolometric luminosity that is tightly constrained between $\mathrm{log}\!\left(L_\mathrm{bol}/L_{\odot}\right)$=-4.31 to $-$4.14, which in turn provides a robust mass constraint of 7.1$\pm$1.2 $M_\mathrm{Jup}$. In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterise the population of exoplanets amenable to high-contrast imaging in greater detail. △ Less

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

Comments: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

The dynamics of interacting multi-pulses in the one-dimensional quintic complex Ginzburg-Landau equation

Authors: T. Rossides, D. J. B. Lloyd, S. Zelik, M. R. Turner

Abstract: We formulate an effective numerical scheme that can readily, and accurately, calculate the dynamics of weakly interacting multi-pulse solutions of the quintic complex Ginzburg-Landau equation (QCGLE) in one space dimension. The scheme is based on a global centre-manifold reduction where one considers the solution of the QCGLE as the composition of individual pulses plus a remainder function, which… ▽ More We formulate an effective numerical scheme that can readily, and accurately, calculate the dynamics of weakly interacting multi-pulse solutions of the quintic complex Ginzburg-Landau equation (QCGLE) in one space dimension. The scheme is based on a global centre-manifold reduction where one considers the solution of the QCGLE as the composition of individual pulses plus a remainder function, which is orthogonal to the adjoint eigenfunctions of the linearised operator about a single pulse. This centre-manifold projection overcomes the difficulties of other, more orthodox, numerical schemes, by yielding a fast-slow system describing 'slow' ordinary differential equations for the locations and phases of the individual pulses, and a 'fast' partial differential equation for the remainder function. With small parameter $ε=e^{-λ_r d}$ where $λ_r$ is a constant and $d>0$ is the pulse separation distance, we write the fast-slow system in terms of first-order and second-order correction terms only, a formulation which is solved more efficiently than the full system. This fast-slow system is integrated numerically using adaptive time-stepping. Results are presented here for two- and three-pulse interactions. For the two-pulse problem, cells of periodic behaviour, separated by an infinite set of heteroclinic orbits, are shown to 'split' under perturbation creating complex spiral behaviour. For the case of three pulse interaction a range of dynamics, including chaotic pulse interaction, are found. While results are presented for pulse interaction in the QCGLE, the numerical scheme can also be applied to a wider class of parabolic PDEs. △ Less

Submitted 21 April, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 33 pages

Tactile Gym 2.0: Sim-to-real Deep Reinforcement Learning for Comparing Low-cost High-Resolution Robot Touch

Authors: Yijiong Lin, John Lloyd, Alex Church, Nathan F. Lepora

Abstract: High-resolution optical tactile sensors are increasingly used in robotic learning environments due to their ability to capture large amounts of data directly relating to agent-environment interaction. However, there is a high barrier of entry to research in this area due to the high cost of tactile robot platforms, specialised simulation software, and sim-to-real methods that lack generality acros… ▽ More High-resolution optical tactile sensors are increasingly used in robotic learning environments due to their ability to capture large amounts of data directly relating to agent-environment interaction. However, there is a high barrier of entry to research in this area due to the high cost of tactile robot platforms, specialised simulation software, and sim-to-real methods that lack generality across different sensors. In this letter we extend the Tactile Gym simulator to include three new optical tactile sensors (TacTip, DIGIT and DigiTac) of the two most popular types, Gelsight-style (image-shading based) and TacTip-style (marker based). We demonstrate that a single sim-to-real approach can be used with these three different sensors to achieve strong real-world performance despite the significant differences between real tactile images. Additionally, we lower the barrier of entry to the proposed tasks by adapting them to an inexpensive 4-DoF robot arm, further enabling the dissemination of this benchmark. We validate the extended environment on three physically-interactive tasks requiring a sense of touch: object pushing, edge following and surface following. The results of our experimental validation highlight some differences between these sensors, which may help future researchers select and customize the physical characteristics of tactile sensors for different manipulations scenarios. △ Less

Submitted 27 July, 2022; v1 submitted 21 July, 2022; originally announced July 2022.

DigiTac: A DIGIT-TacTip Hybrid Tactile Sensor for Comparing Low-Cost High-Resolution Robot Touch

Authors: Nathan F. Lepora, Yijiong Lin, Ben Money-Coomes, John Lloyd

Abstract: Deep learning combined with high-resolution tactile sensing could lead to highly capable dexterous robots. However, progress is slow because of the specialist equipment and expertise. The DIGIT tactile sensor offers low-cost entry to high-resolution touch using GelSight-type sensors. Here we customize the DIGIT to have a 3D-printed sensing surface based on the TacTip family of soft biomimetic opti… ▽ More Deep learning combined with high-resolution tactile sensing could lead to highly capable dexterous robots. However, progress is slow because of the specialist equipment and expertise. The DIGIT tactile sensor offers low-cost entry to high-resolution touch using GelSight-type sensors. Here we customize the DIGIT to have a 3D-printed sensing surface based on the TacTip family of soft biomimetic optical tactile sensors. The DIGIT-TacTip (DigiTac) enables direct comparison between these distinct tactile sensor types. For this comparison, we introduce a tactile robot system comprising a desktop arm, mounts and 3D-printed test objects. We use tactile servo control with a PoseNet deep learning model to compare the DIGIT, DigiTac and TacTip for edge- and surface-following over 3D-shapes. All three sensors performed similarly at pose prediction, but their constructions led to differing performances at servo control, offering guidance for researchers selecting or innovating tactile sensors. All hardware and software for reproducing this study will be openly released. Project website: www.lepora.com/digitac. Project repository: www.github.com/nlepora/digitac-design. △ Less

Submitted 18 July, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: 7 pages. Published in RA-L and accepted in IROS 2022

Journal ref: IEEE Robotics and Automation Letters, 2022

The Beyond the Fence Musical and Computer Says Show Documentary

Authors: Simon Colton, Maria Teresa Llano, Rose Hepworth, John Charnley, Catherine V. Gale, Archie Baron, Francois Pachet, Pierre Roy, Pablo Gervas, Nick Collins, Bob Sturm, Tillman Weyde, Daniel Wolff, James Robert Lloyd

Abstract: During 2015 and early 2016, the cultural application of Computational Creativity research and practice took a big leap forward, with a project where multiple computational systems were used to provide advice and material for a new musical theatre production. Billed as the world's first 'computer musical... conceived by computer and substantially crafted by computer', Beyond The Fence was staged in… ▽ More During 2015 and early 2016, the cultural application of Computational Creativity research and practice took a big leap forward, with a project where multiple computational systems were used to provide advice and material for a new musical theatre production. Billed as the world's first 'computer musical... conceived by computer and substantially crafted by computer', Beyond The Fence was staged in the Arts Theatre in London's West End during February and March of 2016. Various computational approaches to analytical and generative sub-projects were used to bring about the musical, and these efforts were recorded in two 1-hour documentary films made by Wingspan Productions, which were aired on SkyArts under the title Computer Says Show. We provide details here of the project conception and execution, including details of the systems which took on some of the creative responsibility in writing the musical, and the contributions they made. We also provide details of the impact of the project, including a perspective from the two (human) writers with overall control of the creative aspects the musical. △ Less

Submitted 11 May, 2022; originally announced June 2022.

Journal ref: The Seventh International Conference on Computational Creativity, {ICCC} 2016

Factored Conditional Filtering: Tracking States and Estimating Parameters in High-Dimensional Spaces

Authors: Dawei Chen, Samuel Yang-Zhao, John Lloyd, Kee Siong Ng

Abstract: This paper introduces factored conditional filters, new filtering algorithms for simultaneously tracking states and estimating parameters in high-dimensional state spaces. The conditional nature of the algorithms is used to estimate parameters and the factored nature is used to decompose the state space into low-dimensional subspaces in such a way that filtering on these subspaces gives distributi… ▽ More This paper introduces factored conditional filters, new filtering algorithms for simultaneously tracking states and estimating parameters in high-dimensional state spaces. The conditional nature of the algorithms is used to estimate parameters and the factored nature is used to decompose the state space into low-dimensional subspaces in such a way that filtering on these subspaces gives distributions whose product is a good approximation to the distribution on the entire state space. The conditions for successful application of the algorithms are that observations be available at the subspace level and that the transition model can be factored into local transition models that are approximately confined to the subspaces; these conditions are widely satisfied in computer science, engineering, and geophysical filtering applications. We give experimental results on tracking epidemics and estimating parameters in large contact networks that show the effectiveness of our approach. △ Less

Submitted 9 July, 2024; v1 submitted 5 June, 2022; originally announced June 2022.

Comments: 66 pages

MSC Class: 68T37 ACM Class: I.2.6

The JWST Early Release Science Program for the Direct Imaging & Spectroscopy of Exoplanetary Systems

Authors: Sasha Hinkley, Aarynn L. Carter, Shrishmoy Ray, Andrew Skemer, Beth Biller, Elodie Choquet, Maxwell A. Millar-Blanchaer, Stephanie Sallum, Brittany Miles, Niall Whiteford, Polychronis Patapis, Marshall D. Perrin, Laurent Pueyo, Glenn Schneider, Karl Stapelfeldt, Jason Wang, Kimberly Ward-Duong, Brendan P. Bowler, Anthony Boccaletti, Julien H. Girard, Dean Hines, Paul Kalas, Jens Kammerer, Pierre Kervella, Jarron Leisenring , et al. (61 additional authors not shown)

Abstract: The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe e… ▽ More The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe extrasolar planets at mid-infrared wavelengths beyond 5$μ$m, deliver detailed spectroscopy revealing much more precise chemical abundances and atmospheric conditions, and provide sensitivity to analogs of our solar system ice-giant planets at wide orbital separations, an entirely new class of exoplanet. However, in order to maximise the scientific output over the lifetime of the mission, an exquisite understanding of the instrumental performance of JWST is needed as early in the mission as possible. In this paper, we describe our 55-hour Early Release Science Program that will utilize all four JWST instruments to extend the characterisation of planetary mass companions to $\sim$15$μ$m as well as image a circumstellar disk in the mid-infrared with unprecedented sensitivity. Our program will also assess the performance of the observatory in the key modes expected to be commonly used for exoplanet direct imaging and spectroscopy, optimize data calibration and processing, and generate representative datasets that will enable a broad user base to effectively plan for general observing programs in future cycles. △ Less

Submitted 12 September, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 21 pages, 7 figures. Accepted for Publication in PASP

One-Shot Domain-Adaptive Imitation Learning via Progressive Learning

Authors: Dandan Zhang, Wen Fan, John Lloyd, Chenguang Yang, Nathan Lepora

Abstract: Traditional deep learning-based visual imitation learning techniques require a large amount of demonstration data for model training, and the pre-trained models are difficult to adapt to new scenarios. To address these limitations, we propose a unified framework using a novel progressive learning approach comprised of three phases: i) a coarse learning phase for concept representation, ii) a fine… ▽ More Traditional deep learning-based visual imitation learning techniques require a large amount of demonstration data for model training, and the pre-trained models are difficult to adapt to new scenarios. To address these limitations, we propose a unified framework using a novel progressive learning approach comprised of three phases: i) a coarse learning phase for concept representation, ii) a fine learning phase for action generation, and iii) an imaginary learning phase for domain adaptation. Overall, this approach leads to a one-shot domain-adaptive imitation learning framework. We use robotic pouring task as an example to evaluate its effectiveness. Our results show that the method has several advantages over contemporary end-to-end imitation learning approaches, including an improved success rate for task execution and more efficient training for deep imitation learning. In addition, the generalizability to new domains is improved, as demonstrated here with novel background, target container and granule combinations. We believe that the proposed method can be broadly applicable to different industrial or domestic applications that involve deep imitation learning for robotic manipulation, where the target scenarios have high diversity while the human demonstration data is limited. △ Less

Submitted 24 April, 2022; originally announced April 2022.

Comments: under review

Mobility Analysis Workflow (MAW): An accessible, interoperable, and reproducible container system for processing raw mobile data

Authors: Xiangyang Guan, Cynthia Chen, Ian Ren, Ka Yee Yeung, Ling-Hong Hung, Wes J. Lloyd

Abstract: Mobility analysis, or understanding and modeling of people's mobility patterns in terms of when, where, and how people move from one place to another, is fundamentally important as such information is the basis for large-scale investment decisions on the nation's multi-modal transportation infrastructure. Recent rise of using passively generated mobile data from mobile devices have raised question… ▽ More Mobility analysis, or understanding and modeling of people's mobility patterns in terms of when, where, and how people move from one place to another, is fundamentally important as such information is the basis for large-scale investment decisions on the nation's multi-modal transportation infrastructure. Recent rise of using passively generated mobile data from mobile devices have raised questions on using such data for capturing the mobility patterns of a population because: 1) there is a great variety of different kinds of mobile data and their respective properties are unknown; and 2) data pre-processing and analysis methods are often not explicitly reported. The high stakes involved with mobility analysis and issues associated with the passively generated mobile data call for mobility analysis (including data, methods and results) to be accessible to all, interoperable across different computing systems, reproducible and reusable by others. In this study, a container system named Mobility Analysis Workflow (MAW) that integrates data, methods and results, is developed. Built upon the containerization technology, MAW allows its users to easily create, configure, modify, execute and share their methods and results in the form of Docker containers. Tools for operationalizing MAW are also developed and made publicly available on GitHub. One use case of MAW is the comparative analysis for the impacts of different pre-processing and mobility analysis methods on inferred mobility patterns. This study finds that different pre-processing and analysis methods do have impacts on the resulting mobility patterns. The creation of MAW and a better understanding of the relationship between data, methods and resulting mobility patterns as facilitated by MAW represent an important first step toward promoting reproducibility and reusability in mobility analysis with passively-generated data. △ Less

Submitted 19 April, 2022; originally announced April 2022.

MSC Class: 91C20 ACM Class: J.6

Approximate localised dihedral patterns near a Turing instability

Authors: Dan J. Hill, Jason J. Bramburger, David J. B. Lloyd

Abstract: Fully localised patterns involving cellular hexagons or squares have been found experimentally and numerically in various continuum models. However, there is currently no mathematical theory for the emergence of these localised cellular patterns from a quiescent state. A key issue is that standard techniques for one-dimensional patterns have proven insufficient for understanding localisation in hi… ▽ More Fully localised patterns involving cellular hexagons or squares have been found experimentally and numerically in various continuum models. However, there is currently no mathematical theory for the emergence of these localised cellular patterns from a quiescent state. A key issue is that standard techniques for one-dimensional patterns have proven insufficient for understanding localisation in higher dimensions. In this work, we present a comprehensive approach to this problem by using techniques developed in the study of axisymmetric patterns. Our analysis covers localised patterns equipped with a wide range of dihedral symmetries, avoiding a restriction to solutions on a predetermined lattice. The context in this paper is a theory for the emergence of such patterns near a Turing instability for a general class of planar reaction-diffusion equations. Posing the reaction-diffusion system in polar coordinates, we carry out a finite-mode Fourier decomposition in the angular variable to yield a large system of coupled radial ordinary differential equations. We then utilise various radial spatial dynamics methods, such as invariant manifolds, rescaling charts, and normal form analysis, leading to an algebraic matching condition for localised patterns to exist in the finite-mode reduction. This algebraic matching condition is nontrivial, which we solve via a combination of by-hand calculations and Gröbner bases from polynomial algebra to reveal the existence of a plethora of localised dihedral patterns. These results capture the essence of the emergent localised hexagonal patterns witnessed in experiments. Moreover, we combine computer-assisted analysis and a Newton-Kantorovich procedure to prove the existence of localised patches with 6m-fold symmetry for arbitrarily large Fourier decompositions. This includes the localised hexagon patches that have been elusive to analytical treatment. △ Less

Submitted 13 January, 2023; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 60 pages, 15 figures

Journal ref: Nonlinearity 36, 2567 (2023)

In-beam $γ$-ray spectroscopy of $^{32}$Mg via direct reactions

Authors: N. Kitamura, K. Wimmer, T. Miyagi, A. Poves, N. Shimizu, J. A. Tostevin, V. M. Bader, C. Bancroft, D. Barofsky, T. Baugher, D. Bazin, J. S. Berryman, V. Bildstein, A. Gade, N. Imai, T. Kröll, C. Langer, J. Lloyd, E. Lunderberg, F. Nowacki, G. Perdikakis, F. Recchia, T. Redpath, S. Saenz, D. Smalley , et al. (4 additional authors not shown)

Abstract: Background: The nucleus $^{32}$Mg ($N=20$ and $Z=12$) plays a central role in the so-called "island of inversion" where in the ground states $sd$-shell neutrons are promoted to the $fp$-shell orbitals across the shell gap, resulting in the disappearance of the canonical neutron magic number $N=20$. Purpose: The primary goals of this work are to extend the level scheme of $^{32}$Mg, provide spin-pa… ▽ More Background: The nucleus $^{32}$Mg ($N=20$ and $Z=12$) plays a central role in the so-called "island of inversion" where in the ground states $sd$-shell neutrons are promoted to the $fp$-shell orbitals across the shell gap, resulting in the disappearance of the canonical neutron magic number $N=20$. Purpose: The primary goals of this work are to extend the level scheme of $^{32}$Mg, provide spin-parity assignments to excited states, and discuss the microscopic structure of each state through comparisons with theoretical calculations. Method: In-beam $γ$-ray spectroscopy of $^{32}$Mg was performed using two direct-reaction probes, one-neutron (two-proton) knockout reactions on $^{33}$Mg ($^{34}$Si). Final-state exclusive cross sections and parallel momentum distributions were extracted from the experimental data and compared with eikonal-based reaction model calculations combined with shell-model overlap functions. Results: Owing to the remarkable selectivity of the one-neutron and two-proton knockout reactions, a significantly updated level scheme for $^{32}$Mg, which exhibits negative-parity intruder and positive-parity normal states, was constructed. The experimental results were confronted with four different nuclear structure models. Conclusions: In some of these models, different aspects of $^{32}$Mg and the transition into the island of inversion are well described. However, unexplained discrepancies remain, and even with the help of these state-of-the-art theoretical approaches, the structure of this key nucleus is not yet fully captured. △ Less

Submitted 25 February, 2022; originally announced February 2022.

Constraints on the structure and seasonal variations of Triton's atmosphere from the 5 October 2017 stellar occultation and previous observations

Authors: J. Marques Oliveira, B. Sicardy, A. R. Gomes-Júnior, J. L. Ortiz, D. F. Strobel, T. Bertrand, F. Forget, E. Lellouch, J. Desmars, D. Bérard, A. Doressoundiram, J. Lecacheux, R. Leiva, E. Meza, F. Roques, D. Souami, T. Widemann, P. Santos-Sanz, N. Morales, R. Duffard, E. Fernández-Valenzuela, A. J. Castro-Tirado, F. Braga-Ribas, B. E. Morgado, M. Assafin , et al. (212 additional authors not shown)

Abstract: A stellar occultation by Neptune's main satellite, Triton, was observed on 5 October 2017 from Europe, North Africa, and the USA. We derived 90 light curves from this event, 42 of which yielded a central flash detection. We aimed at constraining Triton's atmospheric structure and the seasonal variations of its atmospheric pressure since the Voyager 2 epoch (1989). We also derived the shape of th… ▽ More A stellar occultation by Neptune's main satellite, Triton, was observed on 5 October 2017 from Europe, North Africa, and the USA. We derived 90 light curves from this event, 42 of which yielded a central flash detection. We aimed at constraining Triton's atmospheric structure and the seasonal variations of its atmospheric pressure since the Voyager 2 epoch (1989). We also derived the shape of the lower atmosphere from central flash analysis. We used Abel inversions and direct ray-tracing code to provide the density, pressure, and temperature profiles in the altitude range $\sim$8 km to $\sim$190 km, corresponding to pressure levels from 9 μbar down to a few nanobars. Results. (i) A pressure of 1.18$\pm$0.03 μbar is found at a reference radius of 1400 km (47 km altitude). (ii) A new analysis of the Voyager 2 radio science occultation shows that this is consistent with an extrapolation of pressure down to the surface pressure obtained in 1989. (iii) A survey of occultations obtained between 1989 and 2017 suggests that an enhancement in surface pressure as reported during the 1990s might be real, but debatable, due to very few high S/N light curves and data accessible for reanalysis. The volatile transport model analysed supports a moderate increase in surface pressure, with a maximum value around 2005-2015 no higher than 23 μbar. The pressures observed in 1995-1997 and 2017 appear mutually inconsistent with the volatile transport model presented here. (iv) The central flash structure does not show evidence of an atmospheric distortion. We find an upper limit of 0.0011 for the apparent oblateness of the atmosphere near the 8 km altitude. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 52 pages, 26 figures in the main paper, 2 figures in appendix B, 9 figures in appendix C, 1 long table over 5 pages

Journal ref: A&A 659, A136 (2022)

Computing the Morava K-theory of real Grassmanians using chromatic fixed point theory

Authors: Nicholas J. Kuhn, Christopher J. R. Lloyd

Abstract: We study K(n)*(Gr(d,m)) for all n - the 2-local Morava K-theories of the real Grassmanian Gr(d,m) of d-planes in R^m, about which very little has been previously computed. We conjecture that the Atiyah-Hirzebruch Spectral Sequences computing these all collapse after the first possible non-zero differential, and give much evidence that this is the case. Computational patterns for all n seem similar… ▽ More We study K(n)*(Gr(d,m)) for all n - the 2-local Morava K-theories of the real Grassmanian Gr(d,m) of d-planes in R^m, about which very little has been previously computed. We conjecture that the Atiyah-Hirzebruch Spectral Sequences computing these all collapse after the first possible non-zero differential, and give much evidence that this is the case. Computational patterns for all n seem similar to the known calculation of H*(Gr(d,m);Q), the n=0 case. We use a novel method to show that higher differentials can't occur: we get a lower bound on the size of K(n)*(Gr(d,m)) by constructing an action of C = the cyclic group of order 4, on our Grassmanians, and then applying the chromatic fixed point theory of the authors' previous paper. In essence, we bound the size of K(n)*(Gr(d,m)) from below by computing K(n-1)*(Gr(d,m)^C). Meanwhile, the AHSS after the first differential is determined by Q_n-homology, where Q_n is Milnor's nth primitive operation in mod 2 cohomology. Whenever we are able to calculate this, we have found that it agrees with our lower bound for the size of K(n)*(Gr(d,m)). We have two general families where we prove this: m at most 2^{n+1} and all d, and d=2 and all m and n. Computer calculations have allowed us to check many other examples with larger values of d. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 31 pages

MSC Class: 55M35

Journal ref: Algebr. Geom. Topol. 24 (2024) 919-950

Coexisting normal and intruder configurations in $^{32}$Mg

Authors: N. Kitamura, K. Wimmer, A. Poves, N. Shimizu, J. A. Tostevin, V. M. Bader, C. Bancroft, D. Barofsky, T. Baugher, D. Bazin, J. S. Berryman, V. Bildstein, A. Gade, N. Imai, T. Kröll, C. Langer, J. Lloyd, E. Lunderberg, F. Nowacki, G. Perdikakis, F. Recchia, T. Redpath, S. Saenz, D. Smalley, S. R. Stroberg , et al. (3 additional authors not shown)

Abstract: Situated in the so-called "island of inversion," the nucleus $^{32}$Mg is considered as an archetypal example of the disappearance of magicity at $N=20$. We report on high statistics in-beam spectroscopy of $^{32}$Mg with a unique approach, in that two direct reaction probes with different sensitivities to the underlying nuclear structure are employed at the same time. More specifically, states in… ▽ More Situated in the so-called "island of inversion," the nucleus $^{32}$Mg is considered as an archetypal example of the disappearance of magicity at $N=20$. We report on high statistics in-beam spectroscopy of $^{32}$Mg with a unique approach, in that two direct reaction probes with different sensitivities to the underlying nuclear structure are employed at the same time. More specifically, states in $^{32}$Mg were populated by knockout reactions starting from $^{33}$Mg and $^{34}$Si, lying inside and outside the island of inversion, respectively. The momentum distributions of the reaction residues and the cross sections leading to the individual final states were confronted with eikonal-based reaction calculations, yielding a significantly updated level scheme for $^{32}$Mg and spin-parity assignments. By fully exploiting observables obtained in this measurement, a variety of structures coexisting in 32Mg was unraveled. Comparisons with theoretical predictions based on shell-model overlaps allowed for clear discrimination between different structural models, revealing that the complete theoretical description of this key nucleus is yet to be achieved. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Moving window techniques to model shock wave propagation using the concurrent atomistic-continuum method

Authors: Alexander S. Davis, Jeffrey T. Lloyd, Vinamra Agrawal

Abstract: Atomistic methods have successfully modeled different aspects of shock wave propagation in mate-rials over the past several decades, but they suffer from limitations which restrict the total runtime and system size. Multiscale methods have been able to increase the length and time scales that can be modeled but employing such schemes to simulate wave propagation and evolution through engineering-s… ▽ More Atomistic methods have successfully modeled different aspects of shock wave propagation in mate-rials over the past several decades, but they suffer from limitations which restrict the total runtime and system size. Multiscale methods have been able to increase the length and time scales that can be modeled but employing such schemes to simulate wave propagation and evolution through engineering-scale domains is an active area of research. In this work, we develop two distinct moving window approaches within a Concurrent Atomistic-Continuum (CAC) framework to model shock wave propagation through a one-dimensional monatomic chain. In the first method, the entire CAC system travels with the shock in a conveyor fashion and maintains the shock front in the middle of the overall domain. In the second method, the atomistic region follows the shock by the simultaneous coarsening and refinement of the continuum regions. The CAC and moving window frameworks are verified through dispersion relation studies and phonon wave packet tests. We achieve good agreement between the simulated shock velocities and the values obtained from theory with the conveyor technique, while the coarsen-refine technique allows us to follow the propagating wave front through a large-scale domain. This work showcases the ability of the CAC method to accurately simulate propagating shocks and also demonstrates how a moving window technique can be used in a multiscale framework to study highly nonlinear, transient phenomena. △ Less

Submitted 15 November, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Tactile Sim-to-Real Policy Transfer via Real-to-Sim Image Translation

Authors: Alex Church, John Lloyd, Raia Hadsell, Nathan F. Lepora

Abstract: Simulation has recently become key for deep reinforcement learning to safely and efficiently acquire general and complex control policies from visual and proprioceptive inputs. Tactile information is not usually considered despite its direct relation to environment interaction. In this work, we present a suite of simulated environments tailored towards tactile robotics and reinforcement learning.… ▽ More Simulation has recently become key for deep reinforcement learning to safely and efficiently acquire general and complex control policies from visual and proprioceptive inputs. Tactile information is not usually considered despite its direct relation to environment interaction. In this work, we present a suite of simulated environments tailored towards tactile robotics and reinforcement learning. A simple and fast method of simulating optical tactile sensors is provided, where high-resolution contact geometry is represented as depth images. Proximal Policy Optimisation (PPO) is used to learn successful policies across all considered tasks. A data-driven approach enables translation of the current state of a real tactile sensor to corresponding simulated depth images. This policy is implemented within a real-time control loop on a physical robot to demonstrate zero-shot sim-to-real policy transfer on several physically-interactive tasks requiring a sense of touch. △ Less

Submitted 31 October, 2021; v1 submitted 16 June, 2021; originally announced June 2021.

Probabilistic Discriminative Models Address the Tactile Perceptual Aliasing Problem

Authors: John Lloyd, Yijiong Lin, Nathan F. Lepora

Abstract: In this paper, our aim is to highlight Tactile Perceptual Aliasing as a problem when using deep neural networks and other discriminative models. Perceptual aliasing will arise wherever a physical variable extracted from tactile data is subject to ambiguity between stimuli that are physically distinct. Here we address this problem using a probabilistic discriminative model implemented as a 5-compon… ▽ More In this paper, our aim is to highlight Tactile Perceptual Aliasing as a problem when using deep neural networks and other discriminative models. Perceptual aliasing will arise wherever a physical variable extracted from tactile data is subject to ambiguity between stimuli that are physically distinct. Here we address this problem using a probabilistic discriminative model implemented as a 5-component mixture density network comprised of a deep neural network that predicts the parameters of a Gaussian mixture model. We show that discriminative regression models such as deep neural networks and Gaussian process regression perform poorly on aliased data, only making accurate predictions when the sources of aliasing are removed. In contrast, the mixture density network identifies aliased data with improved prediction accuracy. The uncertain predictions of the model form patterns that are consistent with the various sources of perceptual ambiguity. In our view, perceptual aliasing will become an unavoidable issue for robot touch as the field progresses to training robots that act in uncertain and unstructured environments, such as with deep reinforcement learning. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: 11 pages. Accepted in Robotics: Science & Systems (RSS 2021)

A Robust Controller for Stable 3D Pinching using Tactile Sensing

Authors: Efi Psomopoulou, Nicholas Pestell, Fotios Papadopoulos, John Lloyd, Zoe Doulgeri, Nathan F. Lepora

Abstract: This paper proposes a controller for stable grasping of unknown-shaped objects by two robotic fingers with tactile fingertips. The grasp is stabilised by rolling the fingertips on the contact surface and applying a desired grasping force to reach an equilibrium state. The validation is both in simulation and on a fully-actuated robot hand (the Shadow Modular Grasper) fitted with custom-built optic… ▽ More This paper proposes a controller for stable grasping of unknown-shaped objects by two robotic fingers with tactile fingertips. The grasp is stabilised by rolling the fingertips on the contact surface and applying a desired grasping force to reach an equilibrium state. The validation is both in simulation and on a fully-actuated robot hand (the Shadow Modular Grasper) fitted with custom-built optical tactile sensors (based on the BRL TacTip). The controller requires the orientations of the contact surfaces, which are estimated by regressing a deep convolutional neural network over the tactile images. Overall, the grasp system is demonstrated to achieve stable equilibrium poses on various objects ranging in shape and softness, with the system being robust to perturbations and measurement errors. This approach also has promise to extend beyond grasping to stable in-hand object manipulation with multiple fingers. △ Less

Submitted 1 July, 2021; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 8 pages, 10 figures, 1 appendix. Accepted for publication in IEEE Robotics and Automation Letters and in IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2021). Supplemental video: https://youtu.be/rfQesw3FDA4

Neutron Star Quantum Death by Small Black Holes

Authors: Pierce Giffin, John Lloyd, Samuel D. McDermott, Stefano Profumo

Abstract: Neutron stars can be destroyed by black holes at their center accreting material and eventually swallowing the entire star. Here we note that the accretion model adopted in the literature, based on Bondi accretion or variations thereof, is inadequate for small black holes -- black holes whose Schwarzschild radius is comparable to, or smaller than, the neutron's de Broglie wavelength. In this case,… ▽ More Neutron stars can be destroyed by black holes at their center accreting material and eventually swallowing the entire star. Here we note that the accretion model adopted in the literature, based on Bondi accretion or variations thereof, is inadequate for small black holes -- black holes whose Schwarzschild radius is comparable to, or smaller than, the neutron's de Broglie wavelength. In this case, quantum mechanical aspects of the accretion process cannot be neglected, and give rise to a completely different accretion rate. We show that for the case of black holes seeded by the collapse of bosonic dark matter, this is the case for electroweak-scale dark matter particles. In the case of fermionic dark matter, typically the black holes that would form at the center of a neutron star are more massive, unless the dark matter particle mass is very large, larger than about 10$^{10}$ GeV. We calculate the lifetime of neutron stars harboring a ``small'' black hole, and find that black holes lighter than $\sim 10^{11}$ kg quickly evaporate, leaving no trace. More massive black holes destroy neutron stars via quantum accretion on time-scales much shorter than the age of observed neutron stars. We find that the range where seed black holes inside neutron stars are massive enough that they do not quickly evaporate away, but not so massive that a fluid accretion picture is warranted is limited to between $\sim10^{11}$ and $10^{12}$ kg, but our results are key to accurately determine the actual critical black hole mass corresponding to the onset of neutron star destruction △ Less

Submitted 15 June, 2022; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 7 pages, 2 figures; v2 matches version accepted for publication in PRD

Journal ref: Phys. Rev. D 105, 123030 (2022)

Statistical mechanics of dimers on quasiperiodic Ammann-Beenker tilings

Authors: Jerome Lloyd, Sounak Biswas, Steven H. Simon, S. A. Parameswaran, Felix Flicker

Abstract: We study classical dimers on two-dimensional quasiperiodic Ammann-Beenker (AB) tilings. Despite the lack of periodicity we prove that each infinite tiling admits 'perfect matchings' in which every vertex is touched by one dimer. We introduce an auxiliary 'AB$^*$' tiling obtained from the AB tiling by deleting all 8-fold coordinated vertices. The AB$^*$ tiling is again two-dimensional, infinite, an… ▽ More We study classical dimers on two-dimensional quasiperiodic Ammann-Beenker (AB) tilings. Despite the lack of periodicity we prove that each infinite tiling admits 'perfect matchings' in which every vertex is touched by one dimer. We introduce an auxiliary 'AB$^*$' tiling obtained from the AB tiling by deleting all 8-fold coordinated vertices. The AB$^*$ tiling is again two-dimensional, infinite, and quasiperiodic. The AB$^*$ tiling has a single connected component, which admits perfect matchings. We find that in all perfect matchings, dimers on the AB$^*$ tiling lie along disjoint one-dimensional loops and ladders, separated by 'membranes', sets of edges where dimers are absent. As a result, the dimer partition function of the AB$^*$ tiling factorizes into the product of dimer partition functions along these structures. We compute the partition function and free energy per edge on the AB$^*$ tiling using an analytic transfer matrix approach. Returning to the AB tiling, we find that membranes in the AB$^*$ tiling become 'pseudomembranes', sets of edges which collectively host at most one dimer. This leads to a remarkable discrete scale-invariance in the matching problem. The structure suggests that the AB tiling should exhibit highly inhomogenous and slowly decaying connected dimer correlations. Using Monte Carlo simulations, we find evidence supporting this supposition in the form of connected dimer correlations consistent with power law behaviour. Within the set of perfect matchings we find quasiperiodic analogues to the staggered and columnar phases observed in periodic systems. △ Less

Submitted 13 September, 2022; v1 submitted 1 March, 2021; originally announced March 2021.

Comments: 33 pages, 26 figures

Journal ref: Phys. Rev. B 106, 094202 (2022)

Towards integrated tactile sensorimotor control in anthropomorphic soft robotic hands

Authors: Nathan F. Lepora, Andrew Stinchcombe, Chris Ford, Alfred Brown, John Lloyd, Manuel G. Catalano, Matteo Bianchi, Benjamin Ward-Cherrier

Abstract: In this work, we report on the integrated sensorimotor control of the Pisa/IIT SoftHand, an anthropomorphic soft robot hand designed around the principle of adaptive synergies, with the BRL tactile fingertip (TacTip), a soft biomimetic optical tactile sensor based on the human sense of touch. Our focus is how a sense of touch can be used to control an anthropomorphic hand with one degree of actuat… ▽ More In this work, we report on the integrated sensorimotor control of the Pisa/IIT SoftHand, an anthropomorphic soft robot hand designed around the principle of adaptive synergies, with the BRL tactile fingertip (TacTip), a soft biomimetic optical tactile sensor based on the human sense of touch. Our focus is how a sense of touch can be used to control an anthropomorphic hand with one degree of actuation, based on an integration that respects the hand's mechanical functionality. We consider: (i) closed-loop tactile control to establish a light contact on an unknown held object, based on the structural similarity with an undeformed tactile image; and (ii) controlling the estimated pose of an edge feature of a held object, using a convolutional neural network approach developed for controlling other sensors in the TacTip family. Overall, this gives a foundation to endow soft robotic hands with human-like touch, with implications for autonomous grasping, manipulation, human-robot interaction and prosthetics. Supplemental video: https://youtu.be/ndsxj659bkQ △ Less

Submitted 5 February, 2021; originally announced February 2021.

Comments: 7 pages, 10 figures. Supplemental video: https://youtu.be/ndsxj659bkQ

First principles model for voids decorated by transmutation solutes: Short-range order effects and application to neutron irradiated tungsten

Authors: Duc Nguyen-Manh, Jan S. Wrobel, Michael Klimenkov, Matthew J. Lloyd, Luca Messina, Sergei L. Dudarev

Abstract: Understanding how properties of materials change due to nuclear transmutations is a major challenge for the design of structural components for a fusion power plant. In this study, by combining a first-principles matrix Hamiltonian approach with thermodynamic integration we investigate quasi-steady state configurations of multi-component alloys, containing defects, over a broad range of temperatur… ▽ More Understanding how properties of materials change due to nuclear transmutations is a major challenge for the design of structural components for a fusion power plant. In this study, by combining a first-principles matrix Hamiltonian approach with thermodynamic integration we investigate quasi-steady state configurations of multi-component alloys, containing defects, over a broad range of temperature and composition. The model enables simulating transmutation-induced segregation effects in materials, including tungsten where the phenomenon is strongly pronounced. Finite-temperature analysis shows that voids are decorated by Re and Os, but there is no decoration by tantalum (Ta). The difference between the elements is correlated with the sign of the short range order (SRO) parameter between impurity and vacancy species, in agreement with Atom Probe Tomography (APT) observations of irradiated W-Re, W-Os, W-Ta alloys in the solid solution limit. Statistical analyses of Re and Os impurities in vacancy-rich tungsten show that the SRO effects involving the two solutes are highly sensitive to the background concentration the species. In quaternary W-Re-Os-Vac alloys containing 1.5% Re and 0.1% Os, the SRO Re-Os parameter is negative at 1200K, driving the formation of concentrated Re and Os precipitates. Comparison with experimental Transmission Electron Microscopy (TEM) and APT data on W samples irradiated at the High Flux Reactor (HFR) shows that the model explains the origin of anomalous segregation of transmutation products (Re,Os) to vacancy clusters and voids in the high temperature limit pertinent to the operating conditions of a fusion power plant. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: 29 pages, 13 figures, 3 tables

Journal ref: Phys. Rev. Materials 5, 065401 (2021)

Pose-Based Tactile Servoing: Controlled Soft Touch using Deep Learning

Authors: Nathan F. Lepora, John Lloyd

Abstract: This article describes a new way of controlling robots using soft tactile sensors: pose-based tactile servo (PBTS) control. The basic idea is to embed a tactile perception model for estimating the sensor pose within a servo control loop that is applied to local object features such as edges and surfaces. PBTS control is implemented with a soft curved optical tactile sensor (the BRL TacTip) using a… ▽ More This article describes a new way of controlling robots using soft tactile sensors: pose-based tactile servo (PBTS) control. The basic idea is to embed a tactile perception model for estimating the sensor pose within a servo control loop that is applied to local object features such as edges and surfaces. PBTS control is implemented with a soft curved optical tactile sensor (the BRL TacTip) using a convolutional neural network trained to be insensitive to shear. In consequence, robust and accurate controlled motion over various complex 3D objects is attained. First, we review tactile servoing and its relation to visual servoing, before formalising PBTS control. Then, we assess tactile servoing over a range of regular and irregular objects. Finally, we reflect on the relation to visual servo control and discuss how controlled soft touch gives a route towards human-like dexterity in robots. △ Less

Submitted 2 July, 2021; v1 submitted 4 December, 2020; originally announced December 2020.

Comments: A summary video is available here https://youtu.be/12-DJeRcfn0 *NL and JL contributed equally to this work