Showing 1–50 of 198 results for author: Russell, R

Quantum decoherence by magnetic fluctuations in a candidate axion insulator

Authors: Ruben Saatjian, Kohtaro Yamakawa, Ryan S. Russell, James G. Analytis, John W. Harter

Abstract: In magnetic topological insulators, spontaneous time-reversal symmetry breaking by intrinsic magnetic order can open an energy gap in the topological surface spectrum. In the resulting state, exotic properties like axion electrodynamics, the quantum anomalous Hall effect, and other topological magnetoelectric responses are expected to emerge. A detailed understanding of the magnetic order and its… ▽ More In magnetic topological insulators, spontaneous time-reversal symmetry breaking by intrinsic magnetic order can open an energy gap in the topological surface spectrum. In the resulting state, exotic properties like axion electrodynamics, the quantum anomalous Hall effect, and other topological magnetoelectric responses are expected to emerge. A detailed understanding of the magnetic order and its coupling to the topological surface states is essential to harness and tune these properties. Here, we leverage near-resonant electric quadrupole optical second harmonic generation to probe magnetic fluctuations in the candidate axion insulator EuSn$_2$(As,P)$_2$ across its antiferromagnetic phase boundary. We observe a pronounced dimensional crossover in the quantum decoherence induced by magnetic fluctuations, where two-dimensional in-plane ferromagnetic correlations at high temperatures give way to three-dimensional long-range order at the Néel temperature. We also observe the breaking of rotational symmetry within the long-range-ordered antiferromagnetic state and map out the resulting spatial domain structure. More generally, we demonstrate the unique capabilities of nonlinear optical spectroscopy to study quantum coherence and fluctuations in magnetic quantum materials. △ Less

Submitted 3 July, 2024; originally announced July 2024.

A framework for developing a knowledge management platform

Authors: Marie Lisandra Zepeda Mendoza, Sonali Agarwal, James A. Blackshaw, Vanesa Bol, Audrey Fazzi, Filippo Fiorini, Amy Louise Foreman, Nancy George, Brett R. Johnson, Brian Martin, Dave McComb, Euphemia Mutasa-Gottgens, Helen Parkinson, Martin Romacker, Rolf Russell, Valérien Ségard, Shawn Zheng Kai Tan, Wei Kheng Teh, F. P. Winstanley, Benedict Wong, Adrian M. Smith

Abstract: Knowledge management (KM) involves collecting, organizing, storing, and disseminating information to improve decision-making, innovation, and performance. Implementing KM at scale has become essential for organizations to effectively leverage vast accessible data. This paper is a compilation of concepts that emerged from KM workshops hosted by EMBL-EBI, attended by SMEs and industry. We provide gu… ▽ More Knowledge management (KM) involves collecting, organizing, storing, and disseminating information to improve decision-making, innovation, and performance. Implementing KM at scale has become essential for organizations to effectively leverage vast accessible data. This paper is a compilation of concepts that emerged from KM workshops hosted by EMBL-EBI, attended by SMEs and industry. We provide guidance on envisioning, executing, evaluating, and evolving knowledge management platforms. We emphasize essential considerations such as setting knowledge domain boundaries and measuring success, as well as the importance of making knowledge accessible for downstream applications and non-computational users and highlights necessary personal and organizational skills for success. We stress the importance of collaboration and the need for convergence on shared principles and commitment to provide or seek resources to advance KM. The community is invited to join the journey of KM and contribute to the advancement of the field by applying and improving on the guidelines described. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 18 pages, 1 figure

Deep Modeling of Non-Gaussian Aleatoric Uncertainty

Authors: Aastha Acharya, Caleb Lee, Marissa D'Alonzo, Jared Shamwell, Nisar R. Ahmed, Rebecca Russell

Abstract: Deep learning offers promising new ways to accurately model aleatoric uncertainty in robotic estimation systems, particularly when the uncertainty distributions do not conform to traditional assumptions of being fixed and Gaussian. In this study, we formulate and evaluate three fundamental deep learning approaches for conditional probability density modeling to quantify non-Gaussian aleatoric unce… ▽ More Deep learning offers promising new ways to accurately model aleatoric uncertainty in robotic estimation systems, particularly when the uncertainty distributions do not conform to traditional assumptions of being fixed and Gaussian. In this study, we formulate and evaluate three fundamental deep learning approaches for conditional probability density modeling to quantify non-Gaussian aleatoric uncertainty: parametric, discretized, and generative modeling. We systematically compare the respective strengths and weaknesses of these three methods on simulated non-Gaussian densities as well as on real-world terrain-relative navigation data. Our results show that these deep learning methods can accurately capture complex uncertainty patterns, highlighting their potential for improving the reliability and robustness of estimation systems. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 8 pages, 7 figures

Evaluating experiences in a digital nutrition education program for people with multiple sclerosis: a qualitative study

Authors: RD Russell, J He, LJ Black, A Begley

Abstract: Background Multiple sclerosis (MS) is a complex immune-mediated disease with no currently known cure. There is growing evidence to support the role of diet in reducing some of the symptoms and disease progression in MS, and we previously developed and tested the feasibility of a digital nutrition education program for people with MS. Objective The aim of this study was to explore factors that infl… ▽ More Background Multiple sclerosis (MS) is a complex immune-mediated disease with no currently known cure. There is growing evidence to support the role of diet in reducing some of the symptoms and disease progression in MS, and we previously developed and tested the feasibility of a digital nutrition education program for people with MS. Objective The aim of this study was to explore factors that influenced engagement in the digital nutrition education program, including features influencing capability, opportunity, and motivation to change their dietary behaviours. Methods Semi-structured interviews were conducted with people who MS who completed some or all of the program, until data saturation was reached. Interviews were analysed inductively using thematic analysis. Themes were deductively mapped against the COM-B behaviour change model. Results 16 interviews were conducted with participants who completed all (n=10) or some of the program (n=6). Four themes emerged: 1) Acquiring and validating nutrition knowledge; 2) Influence of time and social support; 3) Getting in early to improve health; and 4) Accounting for food literacy experiences. Discussion This is the first online nutrition program with suitable behavioural supports for people with MS. It highlights the importance of disease-specific and evidence-based nutrition education to support people with MS to make dietary changes. Acquiring nutrition knowledge, coupled with practical support mechanisms such as recipe booklets and goal-setting, emerged as crucial for facilitating engagement with the program. Conclusions When designing education programs for people with MS and other neurological conditions, healthcare professionals and program designers should consider flexible delivery and building peer support to address the needs and challenges faced by participants. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Feasibility of a co-designed online nutrition education program for people with multiple sclerosis

Authors: Rebecca D. Russell, Andrea Begley, Alison Daly, Eleanor Dunlop, Minh N. Pham, Lisa Grech, Lucinda J. Black

Abstract: Objective: Diet quality is important for people with multiple sclerosis (MS), but conflicting online information causes them confusion. People with MS want evidence-based MS-specific information to help them make healthy dietary changes, and we co-designed an asynchronous, online nutrition education program (Eating Well with MS) with the MS community. Our aim was to determine the feasibility of Ea… ▽ More Objective: Diet quality is important for people with multiple sclerosis (MS), but conflicting online information causes them confusion. People with MS want evidence-based MS-specific information to help them make healthy dietary changes, and we co-designed an asynchronous, online nutrition education program (Eating Well with MS) with the MS community. Our aim was to determine the feasibility of Eating Well with MS. Methods: We used a single-arm pre-post design. The feasibility trial was a nine-week intervention with adults with confirmed MS. Feasibility outcomes: 1) demand (recruitment); 2) practicality (completion); 3) acceptability (Intrinsic Motivation Inventory: interest/enjoyment and value/usefulness subscales); and 4) limited efficacy testing (Diet Habits Questionnaire (DHQ); Critical Nutrition Literacy Tool (CNLT); Food Literacy Behaviour Checklist (FLBC)). Results: The recruitment target (n=70) was reached. 87% completed at least one module and 57% completed the full program (five modules). The median interest/enjoyment rating was 5 out of 7 and median value/usefulness rating was 6 out of 7 (where 7 = very true). Compared to pre-program, participants who completed any of the program had statistically significantly improved DHQ, CNLT, and FLBC scores. Conclusion: Eating Well with MS was well received by the MS community and improved their dietary behaviours; demonstrating feasibility. Our findings support the use of co-design methods when developing resources to improve dietary behaviours. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Complex Velocity Structure of Nebular Gas in Active Galaxies Centred in Cooling X-ray Atmospheres

Authors: Marie-Joëlle Gingras, Alison L. Coil, B. R. McNamara, Serena Perrotta, Fabrizio Brighenti, H. R. Russell, S. Peng Oh

Abstract: [OII] emission maps obtained with the Keck Cosmic Web Imager (KCWI) are presented for four galaxies lying at the centers of cooling X-ray cluster atmospheres. Nebular emission reaching altitudes of tens of kpc is found in systems covering a broad range of atmospheric cooling rates, cluster masses, and dynamical states. The central galaxy in Abell 262 hosts high angular momentum gas in a kpc-scale… ▽ More [OII] emission maps obtained with the Keck Cosmic Web Imager (KCWI) are presented for four galaxies lying at the centers of cooling X-ray cluster atmospheres. Nebular emission reaching altitudes of tens of kpc is found in systems covering a broad range of atmospheric cooling rates, cluster masses, and dynamical states. The central galaxy in Abell 262 hosts high angular momentum gas in a kpc-scale disk. The nebular gas in RXJ0820.9+0752 is offset and redshifted with respect to the central galaxy by $10-20$ kpc and $150$ km s$^{-1}$, respectively. The nebular gas in PKS 0745-191 and Abell 1835, both experiencing strong radio-mechanical feedback, is being churned to higher velocity dispersion by the buoyantly rising bubbles and jets. Churning gas flows, likely outflows behind the rising radio bubbles, are likely driven by buoyancy and ram pressure due to the galaxies' motion with respect to the gas. The churned gas is enveloped by larger scale, lower velocity dispersion quiescent nebular emission. The mean radial speeds of the churned gas, quiescent gas, and the central galaxy each differ by up to $\sim 150$ km s$^{-1}$, although speeds upward of $800$ km s$^{-1}$ are found. Nebular gas is dynamically complex due to feedback, motion of the central galaxy, and perhaps relative motion of the hot atmosphere from which it presumably condensed. These motions will affect thermally unstable cooling models, the dispersal of jet energy, and the angular momentum of gas accreting onto the galaxy and its nuclear black hole. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 37 pages, 26 figures, submitted to ApJ

Structure and Dynamics of Magneto-Inertial, Differentially Rotating Laboratory Plasmas

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, D. R. Russell, S. Merlini, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise-time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set… ▽ More We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise-time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum, and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at $\approx 100$ km/s. The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$. As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value $\approx 20$ km/s. We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements. △ Less

Submitted 29 March, 2024; originally announced March 2024.

Comments: Submitted to the Journal of Plasma Physics. 22 pages, 9 figures

Dracula's Chivito: discovery of a large edge-on protoplanetary disk with Pan-STARRS

Authors: Ciprian T. Berghea, Ammar Bayyari, Michael L. Sitko, Jeremy J. Drake, Ana Mosquera, Cecilia Garraffo, Thomas Petit, Ray W. Russell, Korash D. Assani

Abstract: We report the serendipitous discovery of a large edge-on protoplanetary disk associated with the infrared source IRAS 23077+6707. The disk's apparent size in the Pan-STARRS (PS1) images is ~11", making this one of the largest known disks on the sky. It is likely a young system, still surrounded by the envelope which is very faint but still visible in the PS1 images in the northern part (alternativ… ▽ More We report the serendipitous discovery of a large edge-on protoplanetary disk associated with the infrared source IRAS 23077+6707. The disk's apparent size in the Pan-STARRS (PS1) images is ~11", making this one of the largest known disks on the sky. It is likely a young system, still surrounded by the envelope which is very faint but still visible in the PS1 images in the northern part (alternatively this structure could be filaments from the disk itself). We use the PS1 magnitudes and other available photometric data to construct the spectral energy distribution (SED) of the disk. An optical spectrum indicates that the obscured star is hot, most likely late A. We adopt a distance of 300 pc for this object based on Gaia DR3 extinctions. We model the system using the HOCHUNK3D radiative transfer software and find that the system is consistent with a hot star of effective temperature 8000 K surrounded by a disk of size 1650 AU and mass 0.2 M_solar at inclination 82 degrees. △ Less

Submitted 23 April, 2024; v1 submitted 1 February, 2024; originally announced February 2024.

Comments: 11 pages, 4 figures, accepted in ApJL

A cooling flow around the low-redshift quasar H1821+643

Authors: H. R. Russell, P. E. J. Nulsen, A. C. Fabian, T. E. Braben, W. N. Brandt, L. Clews, M. McDonald, C. S. Reynolds, J. S. Sanders, S. Veilleux

Abstract: H1821+643 is the nearest quasar hosted by a galaxy cluster. The energy output by the quasar, in the form of intense radiation and radio jets, is captured by the surrounding hot atmosphere. Here we present a new deep Chandra observation of H1821+643 and extract the hot gas properties into the region where Compton cooling by the quasar radiation is expected to dominate. Using detailed simulations to… ▽ More H1821+643 is the nearest quasar hosted by a galaxy cluster. The energy output by the quasar, in the form of intense radiation and radio jets, is captured by the surrounding hot atmosphere. Here we present a new deep Chandra observation of H1821+643 and extract the hot gas properties into the region where Compton cooling by the quasar radiation is expected to dominate. Using detailed simulations to subtract the quasar light, we show that the soft-band surface brightness of the hot atmosphere increases rapidly by a factor of ~ 30 within the central ~ 10 kpc. The gas temperature drops precipitously to < 0.4 keV and the density increases by over an order of magnitude. The remarkably low metallicity here is likely due to photo-ionization by the quasar emission. The variations in temperature and density are consistent with hydrostatic compression of the hot atmosphere. The extended soft-band peak cannot be explained by an undersubtraction of the quasar or scattered quasar light and is instead due to thermal ISM. The radiative cooling time of the gas falls to only 12 +/- 1 Myr, below the free fall time, and we resolve the sonic radius. H1821+643 is therefore embedded in a cooling flow with a mass deposition rate of up to 3000 Msolar/yr. Multi-wavelength observations probing the star formation rate and cold gas mass are consistent with a cooling flow. We show that the cooling flow extends to much larger radii than can be explained by Compton cooling. Instead, the AGN appears to be underheating the core of this cluster. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 16 pages, 14 figures, accepted to MNRAS

On The Structure of Plasma Jets in the Rotating Plasma Experiment

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, D. R. Russell, S. Merlini, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: Recent pulsed-power experiments have demonstrated the formation of astrophysically-relevant, differentially rotating plasmas [1]. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly-collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication we disc… ▽ More Recent pulsed-power experiments have demonstrated the formation of astrophysically-relevant, differentially rotating plasmas [1]. Key features of the plasma flows are the discovery of a quasi-Keplerian rotation curve, the launching of highly-collimated angular-momentum-transporting axial jets, and a hollow density structure sustained by the centrifugal barrier effect. In this communication we discuss several features of the plasma structure in these experiments through order-of-magnitude models. First, we show that the observed rotation velocity would produce a centrifugal force strong enough to support the hollow density profile. Second, we show that the axial jet should diverge much faster than what was observed, were it not for a magnetized halo with 3T which surrounds the jet and exerts pressure on the interface.Finally, we discuss the temperature structure in the axial jet and plasma halo.We show that a 3T magnetic field would also suppress electron heat conduction,leading to the flat profile observed experimentally. We also find that the axial jet is efficiently radiatively cooled,whereas the halo is not, which would explain the thermal decoupling between the two regions. △ Less

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

Comments: Second version with clarifications on text and figures, and explicit calculations regarding the omega-effect. Manuscript re-submitted to IEEE Transactions On Plasma Science. 4 Figures

Identifying the spin trapped character of the $^{32}$Si isomeric state

Authors: J. Williams, G. Hackman, K. Starosta, R. S. Lubna, Priyanka Choudhary, P. C. Srivastava, C. Andreoiu, D. Annen, H. Asch, M. D. H. K. G. Badanage, G. C. Ball, M. Beuschlein, H. Bidaman, V. Bildstein, R. Coleman, A. B. Garnsworthy, B. Greaves, G. Leckenby, V. Karayonchev, M. S. Martin, C. Natzke, C. M. Petrache, A. Radich, E. Raleigh-Smith, D. Rhodes , et al. (8 additional authors not shown)

Abstract: The properties of a nanosecond isomer in $^{32}$Si, disputed in previous studies, depend on the evolution of proton and neutron shell gaps near the `island of inversion'. We have placed the isomer at 5505.2(2) keV with $J^π = 5^-$, decaying primarily via an $E3$ transition to the $2^+_1$ state. The $E3$ strength of 0.0841(10) W.u. is unusually small and suggests that this isomer is dominated by th… ▽ More The properties of a nanosecond isomer in $^{32}$Si, disputed in previous studies, depend on the evolution of proton and neutron shell gaps near the `island of inversion'. We have placed the isomer at 5505.2(2) keV with $J^π = 5^-$, decaying primarily via an $E3$ transition to the $2^+_1$ state. The $E3$ strength of 0.0841(10) W.u. is unusually small and suggests that this isomer is dominated by the $(νd_{3/2})^{-1} \otimes (νf_{7/2})^{1}$ configuration, which is sensitive to the $N=20$ shell gap. A newly observed $4^+_1$ state is placed at 5881.4(13) keV; its energy is enhanced by the $Z=14$ subshell closure. This indicates that the isomer is located in a `yrast trap', a feature rarely seen at low mass numbers. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: Accepted, Physical Review C

Journal ref: Phys. Rev. C 108, L051305 (2023)

Cooling and Instabilities in Colliding Radiative Flows with Toroidal Magnetic Fields

Authors: R. N. Markwick, A. Frank, E. G. Blackman, J. Carroll-Nellenback, S. V. Lebedev, D. R. Russell, J. W. D. Halliday, L. G. Suttle, P. M. Hartigan

Abstract: We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calcul… ▽ More We report on the results of a simulation based study of colliding magnetized plasma flows. Our set-up mimics pulsed power laboratory astrophysical experiments but, with an appropriate frame change, are relevant to astrophysical jets with internal velocity variations. We track the evolution of the interaction region where the two flows collide. Cooling via radiative loses are included in the calculation. We systematically vary plasma beta ($β_m$) in the flows, the strength of the cooling ($Λ_0$) and the exponent ($α$) of temperature-dependence of the cooling function. We find that for strong magnetic fields a counter-propagating jet called a "spine" is driven by pressure from shocked toroidal fields. The spines eventually become unstable and break apart. We demonstrate how formation and evolution of the spines depends on initial flow parameters and provide a simple analytic model that captures the basic features of the flow. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 14 pages, 16 figures. Submitted to MNRAS

The evolution of galaxies and clusters at high spatial resolution with AXIS

Authors: H. R. Russell, L. A. Lopez, S. W. Allen, G. Chartas, P. P. Choudhury, R. A. Dupke, A. C. Fabian, A. M. Flores, K. Garofali, E. Hodges-Kluck, M. J. Koss, L. Lanz, B. D. Lehmer, J. -T. Li, W. P. Maksym, A. B. Mantz, M. McDonald, E. D. Miller, R. F. Mushotzky, Y. Qiu, C. S. Reynolds, F. Tombesi, P. Tozzi, A. Trindade-Falcao, S. A. Walker , et al. (3 additional authors not shown)

Abstract: Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks and producing a dynamic interstellar medium. On large scales bordering the cosmic web, feedback drives enriched gas out of galaxies and groups, seeding the intergalactic medium with heavy elements. In this way, feedback shapes galaxy evolution by shutting down sta… ▽ More Stellar and black hole feedback heat and disperse surrounding cold gas clouds, launching gas flows off circumnuclear and galactic disks and producing a dynamic interstellar medium. On large scales bordering the cosmic web, feedback drives enriched gas out of galaxies and groups, seeding the intergalactic medium with heavy elements. In this way, feedback shapes galaxy evolution by shutting down star formation and ultimately curtailing the growth of structure after the peak at redshift 2-3. To understand the complex interplay between gravity and feedback, we must resolve both the key physics within galaxies and map the impact of these processes over large scales, out into the cosmic web. The Advanced X-ray Imaging Satellite (AXIS) is a proposed X-ray probe mission for the 2030s with arcsecond spatial resolution, large effective area, and low background. AXIS will untangle the interactions of winds, radiation, jets, and supernovae with the surrounding ISM across the wide range of mass scales and large volumes driving galaxy evolution and trace the establishment of feedback back to the main event at cosmic noon. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 29 pages, 18 figures; this white paper is part of a series commissioned for the AXIS Probe mission concept

The disk of the eruptive protostar V900 Mon; a MATISSE/VLTI and MUSE/VLT perspective

Authors: F. Lykou, P. Ábrahám, F. Cruz-Sáenz de Miera, J. Varga, Á. Kóspál, J. Bouwman, L. Chen, S. Kraus, M. L. Sitko, R. W. Russell, M. Pikhartova

Abstract: In this work, we study the silicate dust content in the disk of one of the youngest eruptive stars, V900 Mon, at the highest angular resolution probing down to the inner 10 au of said disk, and study the historical evolution of the system traced in part by a newly discovered emission clump. We performed high-angular resolution mid-infrared interferometric observations of V900 Mon with MATISSE/VLTI… ▽ More In this work, we study the silicate dust content in the disk of one of the youngest eruptive stars, V900 Mon, at the highest angular resolution probing down to the inner 10 au of said disk, and study the historical evolution of the system traced in part by a newly discovered emission clump. We performed high-angular resolution mid-infrared interferometric observations of V900 Mon with MATISSE/VLTI with a spatial coverage ranging from 38-m to 130-m baselines, and compared them to archival MIDI/VLTI data. We also mined and re-analyzed archival optical and infrared photometry of the star to study its long-term evolution since its eruption in the 1990s. We complemented our findings with integral field spectroscopy data from MUSE/VLT. The MATISSE/VLTI data suggest a radial variation of the silicate feature in the dusty disk, whereby at large spatial scales ($\geq10$ au) the protostellar disk's emission is dominated by large-sized ($\geq1\,μm$) silicate grains, while at smaller spatial scales and closer to the star ($\leq5$ au), silicate emission is absent suggesting self-shielding. We propose that the self-shielding may be the result of small dust grains at the base of the collimated CO outflow previously detected by ALMA. A newly discovered knot in the MUSE/VLT data, located at a projected distance approximately 27,000 au from the star, is co-aligned with the molecular gas outflow at a P.A. of $250^o$ ($\pm5^o$) consistent with the position angle and inclination of the disk. The knot is seen in emission in H$α$, [N II], and the [S II] doublet and its kinematic age is about 5150 years. This ejected material could originate from a previous eruption. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 18 pages, 16 figures, accepted for publication by Astronomy & Astrophysics

Overview of the Advanced X-ray Imaging Satellite (AXIS)

Authors: Christopher S. Reynolds, Erin A. Kara, Richard F. Mushotzky, Andrew Ptak, Michael J. Koss, Brian J. Williams, Steven W. Allen, Franz E. Bauer, Marshall Bautz, Arash Bodaghee, Kevin B. Burdge, Nico Cappelluti, Brad Cenko, George Chartas, Kai-Wing Chan, Lía Corrales, Tansu Daylan, Abraham D. Falcone, Adi Foord, Catherine E. Grant, Mélanie Habouzit, Daryl Haggard, Sven Herrmann, Edmund Hodges-Kluck, Oleg Kargaltsev , et al. (18 additional authors not shown)

Abstract: The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics usin… ▽ More The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics using single-crystal silicon, and developments in the fabrication of large-format, small-pixel, high readout rate CCD detectors with good spectral resolution, allowing a robust and cost-effective design. Further, AXIS will be responsive to target-of-opportunity alerts and, with onboard transient detection, will be a powerful facility for studying the time-varying X-ray universe, following on from the legacy of the Neil Gehrels (Swift) X-ray observatory that revolutionized studies of the transient X-ray Universe. In this paper, we present an overview of AXIS, highlighting the prime science objectives driving the AXIS concept and how the observatory design will achieve these objectives. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: Published in Proceedings of SPIE Optics & Photonics 2023, San Diego

Constraints on thermal conductivity in the merging cluster Abell 2146

Authors: A. Richard-Laferrière, H. R. Russell, A. C. Fabian, U. Chadayammuri, C. S. Reynolds, R. E. A. Canning, A. C. Edge, J. Hlavacek-Larrondo, L. J. King, B. R. McNamara, P. E. J. Nulsen, J. S. Sanders

Abstract: The cluster of galaxies Abell 2146 is undergoing a major merger and is an ideal cluster to study ICM physics, as it has a simple geometry with the merger axis in the plane of the sky, its distance allows us to resolve features across the relevant scales and its temperature lies within Chandra's sensitivity. Gas from the cool core of the subcluster has been partially stripped into a tail of gas, wh… ▽ More The cluster of galaxies Abell 2146 is undergoing a major merger and is an ideal cluster to study ICM physics, as it has a simple geometry with the merger axis in the plane of the sky, its distance allows us to resolve features across the relevant scales and its temperature lies within Chandra's sensitivity. Gas from the cool core of the subcluster has been partially stripped into a tail of gas, which gives a unique opportunity to look at the survival of such gas and determine the rate of conduction in the ICM. We use deep 2.4 Ms Chandra observations of Abell 2146 to produce a high spatial resolution map of the temperature structure along a plume in the ram-pressure stripped tail, described by a partial cone, which is distinguishable from the hot ambient gas. Previous studies of conduction in the ICM typically rely on estimates of the survival time for key structures, such as cold fronts. Here we use detailed hydrodynamical simulations of Abell 2146 to determine the flow velocities along the stripped plume and measure the timescale of the temperature increase along its length. We find that conduction must be highly suppressed by multiple orders of magnitude compared to the Spitzer rate, as the energy used is about 1% of the energy available. We discuss magnetic draping around the core as a possible mechanism for suppressing conduction. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 10 pages, 3 figures, 3 tables, accepted for publication in MNRAS

Wide-Area Geolocalization with a Limited Field of View Camera in Challenging Urban Environments

Authors: Lena M. Downes, Ted J. Steiner, Rebecca L. Russell, Jonathan P. How

Abstract: Cross-view geolocalization, a supplement or replacement for GPS, localizes an agent within a search area by matching ground-view images to overhead images. Significant progress has been made assuming a panoramic ground camera. Panoramic cameras' high complexity and cost make non-panoramic cameras more widely applicable, but also more challenging since they yield less scene overlap between ground a… ▽ More Cross-view geolocalization, a supplement or replacement for GPS, localizes an agent within a search area by matching ground-view images to overhead images. Significant progress has been made assuming a panoramic ground camera. Panoramic cameras' high complexity and cost make non-panoramic cameras more widely applicable, but also more challenging since they yield less scene overlap between ground and overhead images. This paper presents Restricted FOV Wide-Area Geolocalization (ReWAG), a cross-view geolocalization approach that combines a neural network and particle filter to globally localize a mobile agent with only odometry and a non-panoramic camera. ReWAG creates pose-aware embeddings and provides a strategy to incorporate particle pose into the Siamese network, improving localization accuracy by a factor of 100 compared to a vision transformer baseline. This extended work also presents ReWAG*, which improves upon ReWAG's generalization ability in previously unseen environments. ReWAG* repeatedly converges accurately on a dataset of images we have collected in Boston with a 72 degree field of view (FOV) camera, a location and FOV that ReWAG* was not trained on. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Comments: 10 pages, 16 figures. Extension of ICRA 2023 paper arXiv:2209.11854

Multi-Shooting Differential Dynamic Programming for Hybrid Systems using Analytical Derivatives

Authors: Shubham Singh, Ryan P. Russell, Patrick M. Wensing

Abstract: Differential Dynamic Programming (DDP) is a popular technique used to generate motion for dynamic-legged robots in the recent past. However, in most cases, only the first-order partial derivatives of the underlying dynamics are used, resulting in the iLQR approach. Neglecting the second-order terms often slows down the convergence rate compared to full DDP. Multi-Shooting is another popular techni… ▽ More Differential Dynamic Programming (DDP) is a popular technique used to generate motion for dynamic-legged robots in the recent past. However, in most cases, only the first-order partial derivatives of the underlying dynamics are used, resulting in the iLQR approach. Neglecting the second-order terms often slows down the convergence rate compared to full DDP. Multi-Shooting is another popular technique to improve robustness, especially if the dynamics are highly non-linear. In this work, we consider Multi-Shooting DDP for trajectory optimization of a bounding gait for a simplified quadruped model. As the main contribution, we develop Second-Order analytical partial derivatives of the rigid-body contact dynamics, extending our previous results for fixed/floating base models with multi-DoF joints. Finally, we show the benefits of a novel Quasi-Newton method for approximating second-order derivatives of the dynamics, leading to order-of-magnitude speedups in the convergence compared to the full DDP method. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: https://www.youtube.com/watch?v=C0h6mEpcnAE

Direct images and spectroscopy of a giant protoplanet driving spiral arms in MWC 758

Authors: Kevin Wagner, Jordan Stone, Andrew Skemer, Steve Ertel, Ruobing Dong, Dániel Apai, Eckhart Spalding, Jarron Leisenring, Michael Sitko, Kaitlin Kratter, Travis Barman, Mark Marley, Brittany Miles, Anthony Boccaletti, Korash Assani, Ammar Bayyari, Taichi Uyama, Charles E. Woodward, Phil Hinz, Zackery Briesemeister, Kellen Lawson, François Ménard, Eric Pantin, Ray W. Russell, Michael Skrutskie , et al. (1 additional authors not shown)

Abstract: Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets t… ▽ More Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets than past observations at shorter wavelengths. We detect a giant protoplanet, MWC 758c, at a projected separation of ~100 au from the star. The spectrum of MWC 758c is distinct from the rest of the disk and consistent with emission from a planetary atmosphere with Teff = 500 +/- 100 K for a low level of extinction (AV<30), or a hotter object with a higher level of extinction. Both scenarios are commensurate with the predicted properties of the companion responsible for driving the spiral arms. MWC 758c provides evidence that spiral arms in protoplanetary disks can be caused by cold giant planets or by those whose optical emission is highly attenuated. MWC 758c stands out both as one of the youngest giant planets known, and also as one of the coldest and/or most attenuated. Furthermore, MWC 758c is among the first planets to be observed within a system hosting a protoplanetary disk. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Comments: Published in Nature Astronomy

Modeling the Multiwavelength Evolution of the V960 Mon System

Authors: Adolfo S. Carvalho, Lynne A. Hillenbrand, Franz-Josef Hambsch, Shawn Dvorak, Michael Sitko, Ray W. Russell, Victoria Hammond, Michael Connelley, Michael C. B. Ashley, Matthew J. Hankins

Abstract: We study the evolution of the FU Ori object V960 Mon since its outburst, using available multi-wavelength photometric time series over 8 years, complemented by several epochs of moderate-dispersion spectrophotometry. We find that the source fading can be well-described by a decrease in the temperature of the inner disk, which results from a combination of decreasing accretion rate and increasing i… ▽ More We study the evolution of the FU Ori object V960 Mon since its outburst, using available multi-wavelength photometric time series over 8 years, complemented by several epochs of moderate-dispersion spectrophotometry. We find that the source fading can be well-described by a decrease in the temperature of the inner disk, which results from a combination of decreasing accretion rate and increasing inner disk radius. We model the system with a disk atmosphere model that produces the observed variations in multi-band photometry (this paper) and high resolution spectral lines (a companion paper). △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 15 pages, 13 figures, 2 tables, Accepted to ApJ

Radiative cooling effects on reverse shocks formed by magnetised supersonic plasma flows

Authors: S. Merlini, J. D. Hare, G. C. Burdiak, J. W. D. Halliday, A. Ciardi, J. P. Chittenden, T. Clayson, A. J. Crilly, S. J. Eardley, K. E. Marrow, D. R. Russell, R. A. Smith, N. Stuart, L. G. Suttle, E. R. Tubman, V. Valenzuela-Villaseca, T. W. O. Varnish, S. V. Lebedev

Abstract: We study the structure of reverse shocks formed by the collision of supersonic, magnetised plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminium wire arrays, the shock is sharp and well… ▽ More We study the structure of reverse shocks formed by the collision of supersonic, magnetised plasma flows driven by an inverse (or exploding) wire array with a planar conducting obstacle. We observe that the structure of these reverse shocks varies dramatically with wire material, despite the similar upstream flow velocities and mass densities. For aluminium wire arrays, the shock is sharp and well defined, consistent with magneto-hydrodynamic theory. In contrast, we do not observe a well-defined shock using tungsten wires, instead, we see a broad region dominated by density fluctuations on a wide range of spatial scales. We diagnose these two very different interactions using interferometry, Thomson scattering, shadowgraphy, and a newly developed imaging refractometer which is sensitive to small deflections of the probing laser corresponding to small-scale density perturbations. We conclude that the differences in shock structure are most likely due to radiative cooling instabilities which create small-scale density perturbations elongated along magnetic field lines in the tungsten plasma. These instabilities grow more slowly and are smoothed by thermal conduction in the aluminium plasma. △ Less

Submitted 7 August, 2023; v1 submitted 2 June, 2023; originally announced June 2023.

Journal ref: Phys. Plasmas 30, 092102 (2023)

Wavelength-Dependent Extinction and Grain Sizes in Dippers

Authors: Michael L. Sitko, Ray W. Russell, Zachary C. Long, Korash Assani, Monika Pikhartova, Ammar Bayyari, Carol A. Grady, Carey M. Lisse, Massimo Marengo, John P. Wisniewski, William Danchi

Abstract: We have examined inter-night variability of K2-discovered Dippers that are not close to being viewed edge-on, as determined from previously-reported ALMA images, using the SpeX spectrograph and the NASA Infrared Telescope facility (IRTF). The three objects observed were EPIC 203850058, EPIC 205151387, and EPIC 204638512 (2MASS J16042165-2130284). Using the ratio of the fluxes between two successiv… ▽ More We have examined inter-night variability of K2-discovered Dippers that are not close to being viewed edge-on, as determined from previously-reported ALMA images, using the SpeX spectrograph and the NASA Infrared Telescope facility (IRTF). The three objects observed were EPIC 203850058, EPIC 205151387, and EPIC 204638512 (2MASS J16042165-2130284). Using the ratio of the fluxes between two successive nights, we find that for EPIC 204638512 and EPIC 205151387, we find that the properties of the dust differ from that seen in the diffuse interstellar medium and denser molecular clouds. However, the grain properties needed to explain the extinction does resemble those used to model the disks of many young stellar objects. The wavelength-dependent extinction models of both EPIC 204638512 and EPIC 205151387 includes grains at least 500 microns in size, but lacks grains smaller than 0.25 microns. The change in extinction during the dips, and the timescale for these variations to occur, imply obscuration by the surface layers of the inner disks. The recent discovery of a highly mis-inclined inner disk in EPIC 204638512 is suggests that the variations in this disk system may point to due to rapid changes in obscuration by the surface layers of its inner disk, and that other face-on Dippers might have similar geometries. The He I line at 1.083 microns in EPIC 205151387 and EPIC 20463851 were seen to change from night to night, suggesting that we are seeing He I gas mixed in with the surface dust. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: 13 pages, 6 figures, 2 tables

Electronic nematic order in the normal state of strontium ruthenate

Authors: Ryan Russell, Hari P. Nair, Kyle M. Shen, Darrell G. Schlom, John W. Harter

Abstract: Despite significant achievements in characterizing the properties of Sr$_2$RuO$_4$ over the last three decades, the precise nature of its electronic ground state is still unresolved. In this work, we provide a missing piece of the puzzle by uncovering evidence of electronic nematic order in the normal state of Sr$_2$RuO$_4$, revealed by ultrafast time-resolved optical dichroism measurements of uni… ▽ More Despite significant achievements in characterizing the properties of Sr$_2$RuO$_4$ over the last three decades, the precise nature of its electronic ground state is still unresolved. In this work, we provide a missing piece of the puzzle by uncovering evidence of electronic nematic order in the normal state of Sr$_2$RuO$_4$, revealed by ultrafast time-resolved optical dichroism measurements of uniaxially strained thin films. This nematic order, whose domains are aligned by the strain, spontaneously breaks the four-fold rotational symmetry of the crystal. The temperature dependence of the dichroism resembles an Ising-like order parameter, and optical pumping induces a coherent oscillation of its amplitude mode. A microscopic model of intra-unit-cell nematic order is presented, highlighting the importance of Coulomb repulsion between neighboring oxygen $p$-orbitals. The existence of electronic nematic order in the normal state of Sr$_2$RuO$_4$ may have consequences for the form and mechanism of superconductivity in this material. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Journal ref: Physical Review B 108, L081105 (2023)

Picosecond volume expansion drives a later-time insulator-metal transition in a nano-textured Mott Insulator

Authors: Anita Verma, Denis Golež, Oleg Yu. Gorobtsov, Kelson Kaj, Ryan Russell, Jeffrey Z. Kaaret, Erik Lamb, Guru Khalsa, Hari P Nair, Yifei Sun, Ryan Bouck, Nathaniel Schreiber, Jacob P. Ruf, Varun Ramaprasad, Yuya Kubota, Tadashi Togashi, Vladimir A. Stoica, Hari Padmanabhan, John W. Freeland, Nicole A. Benedek, Oleg Shpyrko, John W. Harter, Richard D. Averitt, Darrell G. Schlom, Kyle M. Shen , et al. (2 additional authors not shown)

Abstract: Technology moves towards ever faster switching between different electronic and magnetic states of matter. Manipulating properties at terahertz rates requires accessing the intrinsic timescales of electrons (femtoseconds) and associated phonons (10s of femtoseconds to few picoseconds), which is possible with short-pulse photoexcitation. Yet, in many Mott insulators, the electronic transition is ac… ▽ More Technology moves towards ever faster switching between different electronic and magnetic states of matter. Manipulating properties at terahertz rates requires accessing the intrinsic timescales of electrons (femtoseconds) and associated phonons (10s of femtoseconds to few picoseconds), which is possible with short-pulse photoexcitation. Yet, in many Mott insulators, the electronic transition is accompanied by the nucleation and growth of percolating domains of the changed lattice structure, leading to empirical time scales dominated by slow coarsening dynamics. Here, we use time-resolved X-ray diffraction and reflectivity measurements to investigate the photoinduced insulator-to-metal transition in an epitaxially strained thin film Mott insulator Ca2RuO4. The dynamical transition occurs without observable domain formation and coarsening effects, allowing the study of the intrinsic electronic and lattice dynamics. Above a fluence threshold, the initial electronic excitation drives a fast lattice rearrangement, followed by a slower electronic evolution into a metastable non-equilibrium state. Microscopic calculations based on time-dependent dynamical mean-field theory and semiclassical lattice dynamics within a recently published equilibrium energy landscape picture explain the threshold-behavior and elucidate the delayed onset of the electronic phase transition in terms of kinematic constraints on recombination. Analysis of satellite scattering peaks indicates the persistence of a strain-induced nano-texture in the photoexcited film. This work highlights the importance of combined electronic and structural studies to unravel the physics of dynamic transitions and elucidates the role of strain in tuning the timescales of photoinduced processes. △ Less

Submitted 6 April, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Surrogate Neural Networks for Efficient Simulation-based Trajectory Planning Optimization

Authors: Evelyn Ruff, Rebecca Russell, Matthew Stoeckle, Piero Miotto, Jonathan P. How

Abstract: This paper presents a novel methodology that uses surrogate models in the form of neural networks to reduce the computation time of simulation-based optimization of a reference trajectory. Simulation-based optimization is necessary when there is no analytical form of the system accessible, only input-output data that can be used to create a surrogate model of the simulation. Like many high-fidelit… ▽ More This paper presents a novel methodology that uses surrogate models in the form of neural networks to reduce the computation time of simulation-based optimization of a reference trajectory. Simulation-based optimization is necessary when there is no analytical form of the system accessible, only input-output data that can be used to create a surrogate model of the simulation. Like many high-fidelity simulations, this trajectory planning simulation is very nonlinear and computationally expensive, making it challenging to optimize iteratively. Through gradient descent optimization, our approach finds the optimal reference trajectory for landing a hypersonic vehicle. In contrast to the large datasets used to create the surrogate models in prior literature, our methodology is specifically designed to minimize the number of simulation executions required by the gradient descent optimizer. We demonstrated this methodology to be more efficient than the standard practice of hand-tuning the inputs through trial-and-error or randomly sampling the input parameter space. Due to the intelligently selected input values to the simulation, our approach yields better simulation outcomes that are achieved more rapidly and to a higher degree of accuracy. Optimizing the hypersonic vehicle's reference trajectory is very challenging due to the simulation's extreme nonlinearity, but even so, this novel approach found a 74% better-performing reference trajectory compared to nominal, and the numerical results clearly show a substantial reduction in computation time for designing future trajectories. △ Less

Submitted 30 March, 2023; originally announced March 2023.

Comments: 8 pages, 11 figures, submitted to the IEEE Conference of Decision and Control 2023

SN2017egm: A Helium-rich Superluminous Supernova with Multiple Bumps in the Light Curves

Authors: Jiazheng Zhu, Ning Jiang, Subo Dong, Alexei V. Filippenko, Richard J. Rudy, A. Pastorello, Christopher Ashall, Subhash Bose, R. S. Post, D. Bersier, Stefano Benetti, Thomas G. Brink, Ping Chen, Liming Dou, N. Elias-Rosa, Peter Lundqvist, Seppo Mattila, Ray W. Russell, Michael L. Sitko, Auni Somero, M. D. Stritzinger, Tinggui Wang, Peter J. Brown, E. Cappellaro, Morgan Fraser , et al. (6 additional authors not shown)

Abstract: When discovered, SN~2017egm was the closest (redshift $z=0.03$) hydrogen-poor superluminous supernova (SLSN-I) and a rare case that exploded in a massive and metal-rich galaxy. Thus, it has since been extensively observed and studied. We report spectroscopic data showing strong emission at around He~I $λ$10,830 and four He~I absorption lines in the optical. Consequently, we classify SN~2017egm as… ▽ More When discovered, SN~2017egm was the closest (redshift $z=0.03$) hydrogen-poor superluminous supernova (SLSN-I) and a rare case that exploded in a massive and metal-rich galaxy. Thus, it has since been extensively observed and studied. We report spectroscopic data showing strong emission at around He~I $λ$10,830 and four He~I absorption lines in the optical. Consequently, we classify SN~2017egm as a member of an emerging population of helium-rich SLSNe-I (i.e., SLSNe-Ib). We also present our late-time photometric observations. By combining them with archival data, we analyze high-cadence ultra-violet, optical, and near-infrared light curves spanning from early pre-peak ($\sim -20\,d$) to late phases ($\sim +300\,d$). We obtain its most complete bolometric light curve, in which multiple bumps are identified. None of the previously proposed models can satisfactorily explain all main light-curve features, while multiple interactions between the ejecta and circumstellar material (CSM) may explain the undulating features. The prominent infrared excess with a blackbody luminosity of $10^7$--$10^8\,L_{sun}$ detected in SN~2017egm could originate from the emission of either an echo of a pre-existing dust shell, or newly-formed dust, offering an additional piece of evidence supporting the ejecta-CSM interaction model. Moreover, our analysis of deep $Chandra$ observations yields the tightest-ever constraint on the X-ray emission of an SLSN-I, amounting to an X-ray-to-optical luminosity ratio $\lesssim 10^{-3}$ at late phases ($\sim100-200\,d$), which could help explore its close environment and central engine. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 25 pages, 14 Figures, 4 Tables; accepted for publication in ApJ (Mar. 2023)

Learning to Forecast Aleatoric and Epistemic Uncertainties over Long Horizon Trajectories

Authors: Aastha Acharya, Rebecca Russell, Nisar R. Ahmed

Abstract: Giving autonomous agents the ability to forecast their own outcomes and uncertainty will allow them to communicate their competencies and be used more safely. We accomplish this by using a learned world model of the agent system to forecast full agent trajectories over long time horizons. Real world systems involve significant sources of both aleatoric and epistemic uncertainty that compound and i… ▽ More Giving autonomous agents the ability to forecast their own outcomes and uncertainty will allow them to communicate their competencies and be used more safely. We accomplish this by using a learned world model of the agent system to forecast full agent trajectories over long time horizons. Real world systems involve significant sources of both aleatoric and epistemic uncertainty that compound and interact over time in the trajectory forecasts. We develop a deep generative world model that quantifies aleatoric uncertainty while incorporating the effects of epistemic uncertainty during the learning process. We show on two reinforcement learning problems that our uncertainty model produces calibrated outcome uncertainty estimates over the full trajectory horizon. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: Accepted to ICRA 2023

On Second-Order Derivatives of Rigid-Body Dynamics: Theory & Implementation

Authors: Shubham Singh, Ryan P. Russell, Patrick M. Wensing

Abstract: Model-based control for robots has increasingly been dependent on optimization-based methods like Differential Dynamic Programming and iterative LQR (iLQR). These methods can form the basis of Model-Predictive Control (MPC), which is commonly used for controlling legged robots. Computing the partial derivatives of the dynamics is often the most expensive part of these algorithms, regardless of whe… ▽ More Model-based control for robots has increasingly been dependent on optimization-based methods like Differential Dynamic Programming and iterative LQR (iLQR). These methods can form the basis of Model-Predictive Control (MPC), which is commonly used for controlling legged robots. Computing the partial derivatives of the dynamics is often the most expensive part of these algorithms, regardless of whether analytical methods, Finite Difference, Automatic Differentiation (AD), or Chain-Rule accumulation is used. Since the second-order derivatives of dynamics result in tensor computations, they are often ignored, leading to the use of iLQR, instead of the full second-order DDP method. In this paper, we present analytical methods to compute the second-order derivatives of inverse and forward dynamics for open-chain rigid-body systems with multi-DoF joints and fixed/floating bases. An extensive comparison of accuracy and run-time performance with AD and other methods is provided, including the consideration of code-generation techniques in C/C++ to speed up the computations. For the 36 DoF ATLAS humanoid, the second-order Inverse, and the Forward dynamics derivatives take approx 200 mu s, and approx 2.1 ms respectively, resulting in a 3x speedup over the AD approach. △ Less

Submitted 12 February, 2023; originally announced February 2023.

Comments: Evolved work from https://ieeexplore.ieee.org/document/9981356. arXiv admin note: text overlap with arXiv:2203.01497

AGN Feeding and Feedback in M84: From Kiloparsec Scales to the Bondi Radius

Authors: C. J. Bambic, H. R. Russell, C. S. Reynolds, A. C. Fabian, B. R. McNamara, P. E. J. Nulsen

Abstract: We present the deepest Chandra observation to date of the galaxy M84 in the Virgo Cluster, with over 840 kiloseconds of data provided by legacy observations and a recent 730 kilosecond campaign. The increased signal-to-noise allows us to study the origins of the accretion flow feeding the supermassive black hole in the center of M84 from the kiloparsec scales of the X-ray halo to the Bondi radius,… ▽ More We present the deepest Chandra observation to date of the galaxy M84 in the Virgo Cluster, with over 840 kiloseconds of data provided by legacy observations and a recent 730 kilosecond campaign. The increased signal-to-noise allows us to study the origins of the accretion flow feeding the supermassive black hole in the center of M84 from the kiloparsec scales of the X-ray halo to the Bondi radius, $R_{\rm B}$. Temperature, metallicity, and deprojected density profiles are obtained in four sectors about M84's AGN, extending into the Bondi radius. Rather than being dictated by the potential of the black hole, the accretion flow is strongly influenced by the AGN's bipolar radio jets. Along the jet axis, the density profile is consistent with $n_e \propto r^{-1}$; however, the profiles flatten perpendicular to the jet. Radio jets produce a significant asymmetry in the flow, violating a key assumption of Bondi accretion. Temperature in the inner kiloparsec is approximately constant, with only a slight increase from 0.6 to 0.7 keV approaching $R_{\rm B}$, and there is no evidence for a temperature rise imposed by the black hole. The Bondi accretion rate $\dot{M}_{\rm B}$ exceeds the rate inferred from AGN luminosity and jet power by over four orders of magnitude. In sectors perpendicular to the jet, $\dot{M}_{\rm B}$ measurements agree; however, the accretion rate is $> 4 σ$ lower in the North sector along the jet, likely due to cavities in the X-ray gas. Our measurements provide unique insight into the fueling of AGN responsible for radio mode feedback in galaxy clusters. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 17 pages, 9 figures, 2 tables, submitted to MNRAS

Radio jet-ISM interaction and positive radio-mechanical feedback in Abell 1795

Authors: Prathamesh D. Tamhane, Brian R. McNamara, Helen R. Russell, Francoise Combes, Yu Qiu, Alastair C. Edge, Roberto Maiolino, Andrew C. Fabian, Paul E. J. Nulsen, R. Johnstone, Stefano Carniani

Abstract: We present XSHOOTER observations with previous ALMA, MUSE and $HST$ observations to study the nature of radio-jet triggered star formation and the interaction of radio jets with the interstellar medium in the brightest cluster galaxy (BCG) in the Abell 1795 cluster. Using $HST$ UV data we determined an ongoing star formation rate of 9.3 M$_\odot$ yr$^{-1}$. The star formation follows the global Ke… ▽ More We present XSHOOTER observations with previous ALMA, MUSE and $HST$ observations to study the nature of radio-jet triggered star formation and the interaction of radio jets with the interstellar medium in the brightest cluster galaxy (BCG) in the Abell 1795 cluster. Using $HST$ UV data we determined an ongoing star formation rate of 9.3 M$_\odot$ yr$^{-1}$. The star formation follows the global Kennicutt-Schmidt law, however, it has a low efficiency compared to circumnuclear starbursts in nearby galaxies with an average depletion time of $\sim$1 Gyr. The star formation and molecular gas are offset by $\sim1$ kpc indicating that stars have decoupled from the gas. We detected an arc of high linewidth in ionized gas where electron densities are elevated by a factor of $\sim$4 suggesting a shock front driven by radio jets or peculiar motion of the BCG. An analysis of nebular emission line flux ratios suggests that the gas is predominantly ionized by star formation with a small contribution from shocks. We also calculated the velocity structure function (VSF) of the ionized and molecular gases using velocity maps to characterize turbulent motion in the gas. The ionized gas VSF suggests that the radio jets are driving supersonic turbulence in the gas. Thus radio jets can not only heat the atmosphere on large scales and may quench star formation on longer timescales while triggering star formation in positive feedback on short timescales of a few million years. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Accepted for publication in MNRAS. 22 pages, 13 figures 1 table

Symmetry Detection in Trajectory Data for More Meaningful Reinforcement Learning Representations

Authors: Marissa D'Alonzo, Rebecca Russell

Abstract: Knowledge of the symmetries of reinforcement learning (RL) systems can be used to create compressed and semantically meaningful representations of a low-level state space. We present a method of automatically detecting RL symmetries directly from raw trajectory data without requiring active control of the system. Our method generates candidate symmetries and trains a recurrent neural network (RNN)… ▽ More Knowledge of the symmetries of reinforcement learning (RL) systems can be used to create compressed and semantically meaningful representations of a low-level state space. We present a method of automatically detecting RL symmetries directly from raw trajectory data without requiring active control of the system. Our method generates candidate symmetries and trains a recurrent neural network (RNN) to discriminate between the original trajectories and the transformed trajectories for each candidate symmetry. The RNN discriminator's accuracy for each candidate reveals how symmetric the system is under that transformation. This information can be used to create high-level representations that are invariant to all symmetries on a dataset level and to communicate properties of the RL behavior to users. We show in experiments on two simulated RL use cases (a pusher robot and a UAV flying in wind) that our method can determine the symmetries underlying both the environment physics and the trained RL policy. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Appears in Proceedings of AAAI FSS-22 Symposium "Lessons Learned for Autonomous Assessment of Machine Abilities (LLAAMA)"

Learning and Understanding a Disentangled Feature Representation for Hidden Parameters in Reinforcement Learning

Authors: Christopher Reale, Rebecca Russell

Abstract: Hidden parameters are latent variables in reinforcement learning (RL) environments that are constant over the course of a trajectory. Understanding what, if any, hidden parameters affect a particular environment can aid both the development and appropriate usage of RL systems. We present an unsupervised method to map RL trajectories into a feature space where distance represents the relative diffe… ▽ More Hidden parameters are latent variables in reinforcement learning (RL) environments that are constant over the course of a trajectory. Understanding what, if any, hidden parameters affect a particular environment can aid both the development and appropriate usage of RL systems. We present an unsupervised method to map RL trajectories into a feature space where distance represents the relative difference in system behavior due to hidden parameters. Our approach disentangles the effects of hidden parameters by leveraging a recurrent neural network (RNN) world model as used in model-based RL. First, we alter the standard world model training algorithm to isolate the hidden parameter information in the world model memory. Then, we use a metric learning approach to map the RNN memory into a space with a distance metric approximating a bisimulation metric with respect to the hidden parameters. The resulting disentangled feature space can be used to meaningfully relate trajectories to each other and analyze the hidden parameter. We demonstrate our approach on four hidden parameters across three RL environments. Finally we present two methods to help identify and understand the effects of hidden parameters on systems. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Appears in Proceedings of AAAI FSS-22 Symposium "Lessons Learned for Autonomous Assessment of Machine Abilities (LLAAMA)"

Wide-Area Geolocalization with a Limited Field of View Camera

Authors: Lena M. Downes, Ted J. Steiner, Rebecca L. Russell, Jonathan P. How

Abstract: Cross-view geolocalization, a supplement or replacement for GPS, localizes an agent within a search area by matching images taken from a ground-view camera to overhead images taken from satellites or aircraft. Although the viewpoint disparity between ground and overhead images makes cross-view geolocalization challenging, significant progress has been made assuming that the ground agent has access… ▽ More Cross-view geolocalization, a supplement or replacement for GPS, localizes an agent within a search area by matching images taken from a ground-view camera to overhead images taken from satellites or aircraft. Although the viewpoint disparity between ground and overhead images makes cross-view geolocalization challenging, significant progress has been made assuming that the ground agent has access to a panoramic camera. For example, our prior work (WAG) introduced changes in search area discretization, training loss, and particle filter weighting that enabled city-scale panoramic cross-view geolocalization. However, panoramic cameras are not widely used in existing robotic platforms due to their complexity and cost. Non-panoramic cross-view geolocalization is more applicable for robotics, but is also more challenging. This paper presents Restricted FOV Wide-Area Geolocalization (ReWAG), a cross-view geolocalization approach that generalizes WAG for use with standard, non-panoramic ground cameras by creating pose-aware embeddings and providing a strategy to incorporate particle pose into the Siamese network. ReWAG is a neural network and particle filter system that is able to globally localize a mobile agent in a GPS-denied environment with only odometry and a 90 degree FOV camera, achieving similar localization accuracy as what WAG achieved with a panoramic camera and improving localization accuracy by a factor of 100 compared to a baseline vision transformer (ViT) approach. A video highlight that demonstrates ReWAG's convergence on a test path of several dozen kilometers is available at https://youtu.be/U_OBQrt8qCE. △ Less

Submitted 18 May, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 7 pages, 10 figures. Accepted to ICRA 2023

The structure of 3D collisional magnetized bow shocks in pulsed-power-driven plasma flow

Authors: Rishabh Datta, Danny R. Russell, Iek Tang, Thomas Clayson, Lee G. Suttle, Jeremy P. Chittenden, Sergey V. Lebedev, Jack D. Hare

Abstract: We investigate 3D bow shocks in a highly collisional magnetized aluminum plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ($M_S \sim 7$), super-Alfvénic ($M_A > 1$) magnetized flows with frozen-in magnetic flux ($R_M \gg 1$). These flows collide with an ind… ▽ More We investigate 3D bow shocks in a highly collisional magnetized aluminum plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ($M_S \sim 7$), super-Alfvénic ($M_A > 1$) magnetized flows with frozen-in magnetic flux ($R_M \gg 1$). These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. 3D simulations performed using the resistive magnetohydrodynamic (MHD) code GORGON confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Molecular Flows in Contemporary Active Galaxies and the Efficacy of Radio-Mechanical Feedback

Authors: Prathamesh D. Tamhane, Brian R. McNamara, Helen R. Russell, Alastair C. Edge, Andrew C. Fabian, Paul E. J. Nulsen, Iurii V. Babyk

Abstract: Molecular gas flows are analyzed in 14 cluster galaxies (BCGs) centered in cooling hot atmospheres. The BCGs contain $10^{9}-10^{11}~\rm M_\odot$ of molecular gas, much of which is being moved by radio jets and lobes. The molecular flows and radio jet powers are compared to molecular outflows in 45 active galaxies within $z<0.2$. We seek to understand the relative efficacy of radio, quasar, and st… ▽ More Molecular gas flows are analyzed in 14 cluster galaxies (BCGs) centered in cooling hot atmospheres. The BCGs contain $10^{9}-10^{11}~\rm M_\odot$ of molecular gas, much of which is being moved by radio jets and lobes. The molecular flows and radio jet powers are compared to molecular outflows in 45 active galaxies within $z<0.2$. We seek to understand the relative efficacy of radio, quasar, and starburst feedback over a range of active galaxy types. Molecular flows powered by radio feedback in BCGs are $\sim$10--1000 times larger in extent compared to contemporary galaxies hosting quasar nuclei and starbursts. Radio feedback yields lower flow velocities but higher momenta compared to quasar nuclei, as the molecular gas flows in BCGs are usually $\sim$10--100 times more massive. The product of the molecular gas mass and lifting altitude divided by the AGN or starburst power -- a parameter referred to as the lifting factor -- exceeds starbursts and quasar nuclei by two to three orders of magnitude, respectively. When active, radio feedback is generally more effective at lifting gas in galaxies compared to quasars and starburst winds. The kinetic energy flux of molecular clouds generally lies below and often substantially below a few percent of the driving power. We find tentatively that star formation is suppressed in BCGs relative to other active galaxies, perhaps because these systems rarely form molecular disks that are more impervious to feedback and are better able to promote star formation. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: Accepted for Publication in MNRAS

Journal ref: MNRAS 516 October 2022 Pages 861 882

Testing the Limits of AGN Feedback and the Onset of Thermal Instability in the Most Rapidly Star Forming Brightest Cluster Galaxies

Authors: Michael S. Calzadilla, Michael McDonald, Megan Donahue, Brian R. McNamara, Kevin Fogarty, Massimo Gaspari, Myriam Gitti, Helen R. Russell, Grant R. Tremblay, G. Mark Voit, Francesco Ubertosi

Abstract: We present new, deep, narrow- and broad-band Hubble Space Telescope observations of seven of the most star-forming brightest cluster galaxies (BCGs). Continuum-subtracted [O II] maps reveal the detailed, complex structure of warm ($T \sim 10^4$ K) ionized gas filaments in these BCGs, allowing us to measure spatially-resolved star formation rates (SFRs) of ~60-600 Msun/yr. We compare the SFRs in th… ▽ More We present new, deep, narrow- and broad-band Hubble Space Telescope observations of seven of the most star-forming brightest cluster galaxies (BCGs). Continuum-subtracted [O II] maps reveal the detailed, complex structure of warm ($T \sim 10^4$ K) ionized gas filaments in these BCGs, allowing us to measure spatially-resolved star formation rates (SFRs) of ~60-600 Msun/yr. We compare the SFRs in these systems and others from the literature to their intracluster medium (ICM) cooling rates (dM/dt), measured from archival Chandra X-ray data, finding a best-fit relation of log(SFR) = (1.67+/-0.17) log(dM/dt) + (-3.25+/-0.38) with an intrinsic scatter of 0.39+/-0.09 dex. This steeper-than-unity slope implies an increasingly efficient conversion of hot ($T \sim 10^7$ K) gas into young stars with increasing dM/dt, or conversely a gradual decrease in the effectiveness of AGN feedback in the strongest cool cores. We also seek to understand the physical extent of these multiphase filaments that we observe in cluster cores. We show, for the first time, that the average extent of the multiphase gas is always smaller than the radii at which the cooling time reaches 1 Gyr, the tcool/tff profile flattens, and that X-ray cavities are observed. This implies a close connection between the multiphase filaments, the thermodynamics of the cooling core, and the dynamics of X-ray bubbles. Interestingly, we find a one-to-one correlation between the average extent of cool multiphase filaments and the radius at which the cooling time reaches 0.5 Gyr, which may be indicative of a universal condensation timescale in cluster cores. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 26 pages, 5 figures, 5 tables. Submitted to ApJ

Uncertainty Quantification for Competency Assessment of Autonomous Agents

Authors: Aastha Acharya, Rebecca Russell, Nisar R. Ahmed

Abstract: For safe and reliable deployment in the real world, autonomous agents must elicit appropriate levels of trust from human users. One method to build trust is to have agents assess and communicate their own competencies for performing given tasks. Competency depends on the uncertainties affecting the agent, making accurate uncertainty quantification vital for competency assessment. In this work, we… ▽ More For safe and reliable deployment in the real world, autonomous agents must elicit appropriate levels of trust from human users. One method to build trust is to have agents assess and communicate their own competencies for performing given tasks. Competency depends on the uncertainties affecting the agent, making accurate uncertainty quantification vital for competency assessment. In this work, we show how ensembles of deep generative models can be used to quantify the agent's aleatoric and epistemic uncertainties when forecasting task outcomes as part of competency assessment. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: Accepted at the Workshop on Safe and Reliable Robot Autonomy under Uncertainty at ICRA 2022, Philadelphia, USA

Time-resolved velocity and ion sound speed measurements from simultaneous bow shock imaging and inductive probe measurements

Authors: R. Datta, D. R. Russell, T. Clayson, J. P. Chittenden, S. V. Lebedev, J. D. Hare

Abstract: We present a technique to measure the time-resolved velocity and ion sound speed in magnetized, supersonic high-energy-density plasmas. We place an inductive (`b-dot') probe in a supersonic pulsed-power-driven plasma flow and measure the magnetic field advected by the plasma. As the magnetic Reynolds number is large ($R_M > 10$), the plasma flow advects a magnetic field proportional to the current… ▽ More We present a technique to measure the time-resolved velocity and ion sound speed in magnetized, supersonic high-energy-density plasmas. We place an inductive (`b-dot') probe in a supersonic pulsed-power-driven plasma flow and measure the magnetic field advected by the plasma. As the magnetic Reynolds number is large ($R_M > 10$), the plasma flow advects a magnetic field proportional to the current at the load. This enables us to estimate the plasma flow velocity as a function of time from the delay between the current at the load and the signal at the probe. The supersonic flow also generates a detached hydrodynamic bow shock around the probe, the structure of which depends on the upstream sonic Mach number. By imaging the shock around the probe with a Mach-Zehnder interferometer, we determine the upstream Mach number from the shock Mach angle, which we then use to determine the ion sound speed from the known upstream velocity. We use the measured sound speed to infer the value of $\bar{Z}T_e$, where $\bar{Z}$ is the average ionization, and $T_e$ is the electron temperature. We use this diagnostic to measure the time-resolved velocity and sound speed of a supersonic $(M_S \sim 8)$, super-Alfvénic $(M_A \sim 2)$ aluminum plasma generated during the ablation stage of an exploding wire array on the MAGPIE generator (1.4 MA, 250 ns). Velocity and $\bar{Z}T_e$ measured using this technique agree well with optical Thompson scattering measurements reported in literature, and with 3D resistive MHD simulations in GORGON. △ Less

Submitted 12 May, 2022; originally announced May 2022.

The structure of cluster merger shocks: turbulent width and the electron heating timescale

Authors: H. R. Russell, P. E. J. Nulsen, D. Caprioli, U. Chadayammuri, A. C. Fabian, M. W. Kunz, B. R. McNamara, J. S. Sanders, A. Richard-Laferrière, M. Beleznay, R. E. A. Canning, J. Hlavacek-Larrondo, L. J. King

Abstract: We present a new 2 Ms Chandra observation of the cluster merger Abell 2146, which hosts two huge M~2 shock fronts each ~500 kpc across. For the first time, we resolve and measure the width of cluster merger shocks. The best-fit width for the bow shock is 17+/-1 kpc and for the upstream shock is 10.7+/-0.3 kpc. A narrow collisionless shock will appear broader in projection if its smooth shape is wa… ▽ More We present a new 2 Ms Chandra observation of the cluster merger Abell 2146, which hosts two huge M~2 shock fronts each ~500 kpc across. For the first time, we resolve and measure the width of cluster merger shocks. The best-fit width for the bow shock is 17+/-1 kpc and for the upstream shock is 10.7+/-0.3 kpc. A narrow collisionless shock will appear broader in projection if its smooth shape is warped by local gas motions. We show that both shock widths are consistent with collisionless shocks blurred by local gas motions of 290+/-30 km/s. The upstream shock forms later on in the merger than the bow shock and is therefore expected to be significantly narrower. From the electron temperature profile behind the bow shock, we measure the timescale for the electrons and ions to come back into thermal equilibrium. We rule out rapid thermal equilibration of the electrons with the shock-heated ions at the 6 sigma level. The observed temperature profile instead favours collisional equilibration. For these cluster merger shocks, which have low sonic Mach numbers and propagate through a high $β$ plasma, we find no evidence for electron heating over that produced by adiabatic compression. Our findings are expected to be valid for collisionless shocks with similar parameters in other environments and support the existing picture from the solar wind and supernova remnants. The upstream shock is consistent with this result but has a more complex structure, including a ~2 keV increase in temperature ~50 kpc ahead of the shock. △ Less

Submitted 9 May, 2022; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: 18 pages, 17 figures, 1 table, accepted to MNRAS

Competency Assessment for Autonomous Agents using Deep Generative Models

Authors: Aastha Acharya, Rebecca Russell, Nisar R. Ahmed

Abstract: For autonomous agents to act as trustworthy partners to human users, they must be able to reliably communicate their competency for the tasks they are asked to perform. Towards this objective, we develop probabilistic world models based on deep generative modelling that allow for the simulation of agent trajectories and accurate calculation of tasking outcome probabilities. By combining the streng… ▽ More For autonomous agents to act as trustworthy partners to human users, they must be able to reliably communicate their competency for the tasks they are asked to perform. Towards this objective, we develop probabilistic world models based on deep generative modelling that allow for the simulation of agent trajectories and accurate calculation of tasking outcome probabilities. By combining the strengths of conditional variational autoencoders with recurrent neural networks, the deep generative world model can probabilistically forecast trajectories over long horizons to task completion. We show how these forecasted trajectories can be used to calculate outcome probability distributions, which enable the precise assessment of agent competency for specific tasks and initial settings. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Investigating radiatively driven, magnetised plasmas with a university scale pulsed-power generator

Authors: Jack W. D. Halliday, Aidan Crilly, Jeremy Chittenden, Roberto C. Mancini, Stefano Merlini, Steven Rose, Danny R. Russell, Lee G. Suttle, Vicente Valenzuela-Villaseca, Simon N. Bland, Sergey V. Lebedev

Abstract: We present first results from a novel experimental platform which is able to access physics relevant to topics including indirect-drive magnetised ICF; laser energy deposition; various topics in atomic physics; and laboratory astrophysics (for example the penetration of B-fields into HED plasmas). This platform uses the X-Rays from a wire array Z-Pinch to irradiate a silicon target, producing an o… ▽ More We present first results from a novel experimental platform which is able to access physics relevant to topics including indirect-drive magnetised ICF; laser energy deposition; various topics in atomic physics; and laboratory astrophysics (for example the penetration of B-fields into HED plasmas). This platform uses the X-Rays from a wire array Z-Pinch to irradiate a silicon target, producing an outflow of ablated plasma. The ablated plasma expands into ambient, dynamically significant B-fields (~5 T) which are supported by the current flowing through the Z-Pinch. The outflows have a well-defined (quasi-1D) morphology, enabling the study of fundamental processes typically only available in more complex, integrated schemes. Experiments were fielded on the MAGPIE pulsed-power generator (1.4 MA, 240 ns rise time). On this machine a wire array Z-Pinch produces an X-Ray pulse carrying a total energy of ~15 kJ over ~30 ns. This equates to an average brightness temperature of around 10 eV on-target. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Analytical Second-Order Partial Derivatives of Rigid-Body Inverse Dynamics

Authors: Shubham Singh, Ryan P. Russell, Patrick M. Wensing

Abstract: Optimization-based robot control strategies often rely on first-order dynamics approximation methods, as in iLQR. Using second-order approximations of the dynamics is expensive due to the costly second-order partial derivatives of the dynamics with respect to the state and control. Current approaches for calculating these derivatives typically use automatic differentiation (AD) and chain-rule accu… ▽ More Optimization-based robot control strategies often rely on first-order dynamics approximation methods, as in iLQR. Using second-order approximations of the dynamics is expensive due to the costly second-order partial derivatives of the dynamics with respect to the state and control. Current approaches for calculating these derivatives typically use automatic differentiation (AD) and chain-rule accumulation or finite-difference. In this paper, for the first time, we present analytical expressions for the second-order partial derivatives of inverse dynamics for open-chain rigid-body systems with floating base and multi-DoF joints. A new extension of spatial vector algebra is proposed that enables the analysis. A recursive algorithm with complexity of $\mathcal{O}(Nd^2)$ is also provided where $N$ is the number of bodies and $d$ is the depth of the kinematic tree. A comparison with AD in CasADi shows speedups of 1.5-3$\times$ for serial kinematic trees with $N> 5$, and a C++ implementation shows runtimes of $\approx$51$μs$ for a quadruped. △ Less

Submitted 14 August, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: Accepted for IROS 2022 (Oct 23-27, 2022 Kyoto)

Details of Second-Order Partial Derivatives of Rigid-Body Inverse Dynamics

Authors: Shubham Singh, Ryan P. Russell, Patrick M. Wensing

Abstract: The details of second-order partial derivatives of rigid-body Inverse/Forward dynamics are provided. Several properties and identities using Spatial Vector Algebra are listed, along with their detailed derivations. The expressions build upon previous work by the author on first-order partial derivatives of inverse dynamics. The first/second-order derivatives are also extended for systems with exte… ▽ More The details of second-order partial derivatives of rigid-body Inverse/Forward dynamics are provided. Several properties and identities using Spatial Vector Algebra are listed, along with their detailed derivations. The expressions build upon previous work by the author on first-order partial derivatives of inverse dynamics. The first/second-order derivatives are also extended for systems with external forces. Finally, the KKT Forward dynamics and Impact dynamics derivatives are derived. △ Less

Submitted 29 July, 2023; v1 submitted 1 March, 2022; originally announced March 2022.

Characterization of Quasi-Keplerian, Differentially Rotating, Free-Boundary Laboratory Plasmas

Authors: V. Valenzuela-Villaseca, L. G. Suttle, F. Suzuki-Vidal, J. W. D. Halliday, S. Merlini, D. R. Russell, E. R. Tubman, J. D. Hare, J. P. Chittenden, M. E. Koepke, E. G. Blackman, S. V. Lebedev

Abstract: We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gra… ▽ More We present results from pulsed-power driven differentially rotating plasma experiments designed to simulate physics relevant to astrophysical disks and jets. In these experiments, angular momentum is injected by the ram pressure of the ablation flows from a wire array Z pinch. In contrast to previous liquid metal and plasma experiments, rotation is not driven by boundary forces. Axial pressure gradients launch a rotating plasma jet upwards, which is confined by a combination of ram, thermal, and magnetic pressure of a surrounding plasma halo. The jet has subsonic rotation, with a maximum rotation velocity $23 \pm 3$ km/s. The rotational velocity profile is quasi-Keplerian with a positive Rayleigh discriminant $κ^2 \propto r^{-2.8\pm0.8}$ rad$^2$/s$^2$. The plasma completes $0.5 - 2$ full rotations in the experimental time frame ($\sim 150$ ns). △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 4 figures

Perpendicular subcritical shock structure in a collisional plasma experiment

Authors: D. R. Russell, G. C. Burdiak, J. J. Carroll-Nellenback, J. W. D. Halliday, J. D. Hare, S. Merlini, L. G. Suttle, V. Valenzuela-Villaseca, S. J. Eardley, J. A. Fullalove, G. C. Rowland, R. A. Smith, A. Frank, P. Hartigan, A. L. Velikovich, S. V. Lebedev

Abstract: We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the super-magnetosonic outflow from an inverse wire array z-pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absenc… ▽ More We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the super-magnetosonic outflow from an inverse wire array z-pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absence of a hydrodynamic jump. We calculate the classical (Spitzer) resistive diffusion length and show that it is approximately equal to the shock width. We measure little heating across the shock (< 10 % of the ion kinetic energy) which is consistent with an absence of viscous dissipation. △ Less

Submitted 22 January, 2022; originally announced January 2022.

Comments: 6 pages, 5 figures

A Time-Resolved Imaging System for the Diagnosis of X-ray Self-Emission in High Energy Density Physics Experiments

Authors: Jack W. D. Halliday, Simon N. Bland, Jack D. Hare, Susan Parker, Lee G. Suttle, Danny R. Russell, Sergey V. Lebedev

Abstract: A diagnostic capable of recording spatially and temporally resolved X-ray self emission data was developed to characterise experiments on the MAGPIE pulsed-power generator. The diagnostic used two separate imaging systems: A pinhole imaging system with two dimensional spatial resolution and a slit imaging system with one dimensional spatial resolution. The two dimensional imaging system imaged lig… ▽ More A diagnostic capable of recording spatially and temporally resolved X-ray self emission data was developed to characterise experiments on the MAGPIE pulsed-power generator. The diagnostic used two separate imaging systems: A pinhole imaging system with two dimensional spatial resolution and a slit imaging system with one dimensional spatial resolution. The two dimensional imaging system imaged light onto image plate. The one dimensional imaging system imaged light onto the same piece of image plate and a linear array of silicon photodiodes. This design allowed the cross-comparison of different images, allowing a picture of the spatial and temporal distribution of X-ray self emission to be established. The design was tested in a series of pulsed-power driven magnetic-reconnection experiments. △ Less

Submitted 22 November, 2021; originally announced November 2021.

New constraints on light Axion-Like Particles using Chandra Transmission Grating Spectroscopy of the powerful cluster-hosted quasar H1821+643

Authors: Júlia Sisk Reynés, James H. Matthews, Christopher S. Reynolds, Helen R. Russell, Robyn N. Smith, M. C. David Marsh

Abstract: Axion-Like Particles (ALPs) are predicted by several Beyond the Standard Model theories, in particular, string theory. In the presence of an external magnetic field perpendicular to the direction of propagation, ALPs can couple to photons. Therefore, if an X-ray source is viewed through a magnetised plasma, such as a luminous quasar in a galaxy cluster, we may expect spectral distortions that are… ▽ More Axion-Like Particles (ALPs) are predicted by several Beyond the Standard Model theories, in particular, string theory. In the presence of an external magnetic field perpendicular to the direction of propagation, ALPs can couple to photons. Therefore, if an X-ray source is viewed through a magnetised plasma, such as a luminous quasar in a galaxy cluster, we may expect spectral distortions that are well described by photon-ALP oscillations. We present a $571 \ \mathrm{ks}$ combined High and Low Energy Transmission Grating (HETG/LETG) Chandra observation of the powerful radio-quiet quasar H1821+643, hosted by a cool-core cluster at redshift $0.3$. The spectrum is well described by a double power-law continuum and broad$+$narrow iron line emission typical of type-1 Active Galactic Nuclei (AGN), with remaining spectral features $< 2.5\%$. Using a cell-based approach to describe the turbulent cluster magnetic field, we compare our spectrum with photon-ALP mixing curves for 500 field realisations assuming that the thermal-to-magnetic pressure ratio remains constant up to the virial radius. At $99.7\%$ credibility and taking $β= 100$, we exclude all couplings $g_\mathrm{aγ} > 6.3 \times 10^{-13} \ {\mathrm{GeV}}^{-1}$ for most ALP masses $< 10^{-12} \ \mathrm{eV}$. Our results are moderately more sensitive to constraining ALPs than the best previous result from Chandra observations of the Perseus cluster, albeit with a less constrained field model. We reflect on the promising future of ALP studies with bright AGN embedded in rich clusters, especially with the upcoming Athena mission. △ Less

Submitted 25 November, 2021; v1 submitted 7 September, 2021; originally announced September 2021.

Comments: 14 pages, 7 figures. Accepted for publication in MNRAS on 24th November 2021. Comments welcome

Manipulating thermal fields with inhomogeneous heat spreaders

Authors: Eleanor R. Russell, Raphaël C. Assier, William J. Parnell

Abstract: We design a class of spatially inhomogeneous heat spreaders in the context of steady-state thermal conduction leading to spatially uniform thermal fields across a large convective surface. Each spreader has a funnel-shaped design, either in the form of a trapezoidal prism or truncated cone, and is forced by a thermal source at its base. We employ transformation-based techniques, commonly used to s… ▽ More We design a class of spatially inhomogeneous heat spreaders in the context of steady-state thermal conduction leading to spatially uniform thermal fields across a large convective surface. Each spreader has a funnel-shaped design, either in the form of a trapezoidal prism or truncated cone, and is forced by a thermal source at its base. We employ transformation-based techniques, commonly used to study metamaterials, to determine the require thermal conductivity for the spreaders. The obtained materials, although strongly anisotropic and inhomogeneous, can be accurately approximated by assembling isotropic, homogeneous layers, rendering them realisable. An alternative approach is then considered for the conical and trapezoidal spreaders by dividing them into two or three isotropic, homogeneous components respectively. We refer to these simple configurations as neutral layers. All designs are validated numerically both with and without the effects of thermal contact resistance between interfaces. Such novel designs pave the way for future materials that can manipulate and control the flow of heat, helping to solve traditional heat transfer problems such as controlling the temperature of an object and energy harvesting. △ Less

Submitted 1 February, 2022; v1 submitted 6 August, 2021; originally announced August 2021.

Efficient Analytical Derivatives of Rigid-Body Dynamics using Spatial Vector Algebra

Authors: Shubham Singh, Ryan P. Russell, Patrick M. Wensing

Abstract: An essential need for many model-based robot control algorithms is the ability to quickly and accurately compute partial derivatives of the equations of motion. State of the art approaches to this problem often use analytical methods based on the chain rule applied to existing dynamics algorithms. Although these methods are an improvement over finite differences in terms of accuracy, they are not… ▽ More An essential need for many model-based robot control algorithms is the ability to quickly and accurately compute partial derivatives of the equations of motion. State of the art approaches to this problem often use analytical methods based on the chain rule applied to existing dynamics algorithms. Although these methods are an improvement over finite differences in terms of accuracy, they are not always the most efficient. In this paper, we contribute new closed-form expressions for the first-order partial derivatives of inverse dynamics, leading to a recursive algorithm. The algorithm is benchmarked against chain-rule approaches in Fortran and against an existing algorithm from the Pinocchio library in C++. Tests consider computing the partial derivatives of inverse and forward dynamics for robots ranging from kinematic chains to humanoids and quadrupeds. Compared to the previous open-source Pinocchio implementation, our new analytical results uncover a key computational restructuring that enables efficiency gains. Speedups of up to 1.4x are reported for calculating the partial derivatives of inverse dynamics for the 50-dof Talos humanoid. △ Less

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

Journal ref: IEEE Robotics and Automation Letters

Suppressed cooling and turbulent heating in the core of X-ray luminous clusters RXCJ1504.1-0248 and Abell 1664

Authors: Haonan Liu, Andrew C. Fabian, Ciro Pinto, Helen R. Russell, Jeremy S. Sanders, Brian R. McNamara

Abstract: We present the analysis of XMM-Newton observations of two X-ray luminous cool core clusters, RXCJ1504.1-0248 and Abell 1664. The Reflection Grating Spectrometer reveals a radiative cooling rate of $180\pm 40\, \rm M_{\odot}\rm\,yr^{-1}$ and $34\pm 6\, \rm M_{\odot}\rm\,yr^{-1}$ in RXCJ1504.1-0248 and Abell 1664 for gas above 0.7 keV, respectively. These cooling rates are higher than the star forma… ▽ More We present the analysis of XMM-Newton observations of two X-ray luminous cool core clusters, RXCJ1504.1-0248 and Abell 1664. The Reflection Grating Spectrometer reveals a radiative cooling rate of $180\pm 40\, \rm M_{\odot}\rm\,yr^{-1}$ and $34\pm 6\, \rm M_{\odot}\rm\,yr^{-1}$ in RXCJ1504.1-0248 and Abell 1664 for gas above 0.7 keV, respectively. These cooling rates are higher than the star formation rates observed in the clusters, and support simultaneous star formation and molecular gas mass growth on a timescale of 3$\times 10^8$ yr or longer. At these rates, the energy of the X-ray cooling gas is inadequate to power the observed UV/optical line-emitting nebulae, which suggests additional strong heating. No significant residual cooling is detected below 0.7 keV in RXCJ1504.1-0248. By simultaneously fitting the first and second order spectra, we place an upper limit on turbulent velocity of 300 km$\rm s^{-1}$ at 90 per cent confidence level for the soft X-ray emitting gas in both clusters. The turbulent energy density is considered to be less than 8.9 and 27 per cent of the thermal energy density in RXCJ1504.1-0248 and Abell 1664, respectively. This means it is insufficient for AGN heating to fully propagate throughout the cool core via turbulence. We find the cool X-ray component of Abell 1664 ($\sim$0.8 keV) is blueshifted from the systemic velocity by 750$^{+800}_{-280}$ km$\rm s^{-1}$. This is consistent with one component of the molecular gas in the core and suggests a similar dynamical structure for the two phases. We find that an intrinsic absorption model allows the cooling rate to increase to $520\pm 30\, \rm M_{\odot}\rm\,yr^{-1}$ in RXCJ1504.1-0248. △ Less

Submitted 9 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS