Showing 1–50 of 228 results for author: Ball, J

A nonlinear elasticity model in computer vision

Authors: John M. Ball, Christopher L. Horner

Abstract: The purpose of this paper is to analyze a nonlinear elasticity model previously introduced by the authors for comparing two images, regarded as bounded open subsets of $\R^n$ together with associated vector-valued intensity maps. Optimal transformations between the images are sought as minimisers of an integral functional among orientation-preserving homeomorphisms. The existence of minimisers is… ▽ More The purpose of this paper is to analyze a nonlinear elasticity model previously introduced by the authors for comparing two images, regarded as bounded open subsets of $\R^n$ together with associated vector-valued intensity maps. Optimal transformations between the images are sought as minimisers of an integral functional among orientation-preserving homeomorphisms. The existence of minimisers is proved under natural coercivity and polyconvexity conditions, assuming only that the intensity functions are bounded measurable. Variants of the existence theorem are also proved, first under the constraint that finite sets of landmark points in the two images are mapped one to the other, and second when one image is to be compared to an unknown part of another. The question is studied as to whether for images related by a linear mapping the unique minimizer is given by that linear mapping. For a natural class of functional integrands an example is given guaranteeing that this property holds for pairs of images in which the second is a scaling of the first by a constant factor. However for the property to hold for arbitrary pairs of linearly related images it is shown that the integrand has to depend on the gradient of the transformation as a convex function of its determinant alone. This suggests a new model in which the integrand depends also on second derivatives of the transformation, and an example is given for which both existence of minimizers is assured and the above property holds for all pairs of linearly related images. △ Less

Submitted 30 August, 2024; originally announced August 2024.

MSC Class: 94A08; 74B20

Physics of the low momentum diffusivity regime in tokamaks and its experimental applicability

Authors: Haomin Sun, Justin Ball, Stephan Brunner, Anthony Field, Bhavin Patel, Alessandro Balestri, Daniel Kennedy, Colin Roach, Diego Jose Cruz-Zabala, Fernando Puentes Del Pozo, Eleonora Viezzer, Manuel Garcia Munoz

Abstract: Strong $E\times B$ plasma flow shear is beneficial for reducing turbulent transport. However, traditional methods of driving flow shear do not scale well to large devices such as future fusion power plants. In this paper, we use a large number of nonlinear gyrokinetic simulations to study a novel approach to increase flow shear: decreasing the momentum diffusivity to make the plasma ``easier to pu… ▽ More Strong $E\times B$ plasma flow shear is beneficial for reducing turbulent transport. However, traditional methods of driving flow shear do not scale well to large devices such as future fusion power plants. In this paper, we use a large number of nonlinear gyrokinetic simulations to study a novel approach to increase flow shear: decreasing the momentum diffusivity to make the plasma ``easier to push''. We first use an idealized circular geometry and find that one can obtain low momentum diffusivity at tight aspect ratio, low safety factor, high magnetic shear and low temperature gradient. This is the so-called Low Momentum Diffusivity (LMD) regime. To drive intrinsic momentum flux, we then tilt the flux surface, making it up-down asymmetric. In the LMD regime, this intrinsic momentum flux drives strong flow shear that can significantly reduce the heat flux and increase the critical temperature gradient. We also consider the actual experimental geometry of the MAST tokamak to illustrate that this strategy can be practical and create experimentally significant flow shear. Lastly, a preliminary prediction for the SMART tokamak is made. △ Less

Submitted 22 August, 2024; originally announced August 2024.

Comments: 43 pages, 23 figures

URCDM: Ultra-Resolution Image Synthesis in Histopathology

Authors: Sarah Cechnicka, James Ball, Matthew Baugh, Hadrien Reynaud, Naomi Simmonds, Andrew P. T. Smith, Catherine Horsfield, Candice Roufosse, Bernhard Kainz

Abstract: Diagnosing medical conditions from histopathology data requires a thorough analysis across the various resolutions of Whole Slide Images (WSI). However, existing generative methods fail to consistently represent the hierarchical structure of WSIs due to a focus on high-fidelity patches. To tackle this, we propose Ultra-Resolution Cascaded Diffusion Models (URCDMs) which are capable of synthesising… ▽ More Diagnosing medical conditions from histopathology data requires a thorough analysis across the various resolutions of Whole Slide Images (WSI). However, existing generative methods fail to consistently represent the hierarchical structure of WSIs due to a focus on high-fidelity patches. To tackle this, we propose Ultra-Resolution Cascaded Diffusion Models (URCDMs) which are capable of synthesising entire histopathology images at high resolutions whilst authentically capturing the details of both the underlying anatomy and pathology at all magnification levels. We evaluate our method on three separate datasets, consisting of brain, breast and kidney tissue, and surpass existing state-of-the-art multi-resolution models. Furthermore, an expert evaluation study was conducted, demonstrating that URCDMs consistently generate outputs across various resolutions that trained evaluators cannot distinguish from real images. All code and additional examples can be found on GitHub. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Comments: arXiv admin note: text overlap with arXiv:2312.01152

A zero dimensional study of Ohmically heated negative triangularity tokamaks

Authors: Alessandro Balestri, Justin Ball, Stefano Coda

Abstract: Because negative triangularity plasma scenarios remain in L-mode, they do not require external heating systems that exceed the H-mode power threshold. Operating with less heating has the potential to improve performance as heating generally degrades confinement in tokamaks. Using simple zero dimensional power balance and standard empirical scaling laws for confinement, we analyze the impact of ext… ▽ More Because negative triangularity plasma scenarios remain in L-mode, they do not require external heating systems that exceed the H-mode power threshold. Operating with less heating has the potential to improve performance as heating generally degrades confinement in tokamaks. Using simple zero dimensional power balance and standard empirical scaling laws for confinement, we analyze the impact of external heating on several different reactor-relevant devices (i.e. SPARC, MANTA, ITER and DEMO). We compare the nominal externally heated scenarios against equivalent negative tringularity cases without external heating. For most of these devices, the Ohmically heated negative triangularity versions achieve better performance, particularly for devices with high magnetic field and/or high fusion gain. △ Less

Submitted 8 July, 2024; originally announced July 2024.

A green solvent system for precursor phase-engineered sequential deposition of stable formamidinium lead triiodide for perovskite solar cells

Authors: Benjamin M. Gallant, Philippe Holzhey, Joel A. Smith, Saqlain Choudhary, Karim A. Elmestekawy, Pietro Caprioglio, Igal Levine, Alex Sheader, Fengning Yang, Daniel T. W. Toolan, Rachel C. Kilbride, Augustin K. A. Zaininger, James M. Ball, M. Greyson Christoforo, Nakita Noel, Laura M. Herz, Dominik J. Kubicki, Henry J. Snaith

Abstract: Perovskite solar cells (PSCs) offer an efficient, inexpensive alternative to current photovoltaic technologies, with the potential for manufacture via high-throughput coating methods. However, challenges for commercial-scale solution-processing of metal-halide perovskites include the use of harmful solvents, the expense of maintaining controlled atmospheric conditions, and the inherent instabiliti… ▽ More Perovskite solar cells (PSCs) offer an efficient, inexpensive alternative to current photovoltaic technologies, with the potential for manufacture via high-throughput coating methods. However, challenges for commercial-scale solution-processing of metal-halide perovskites include the use of harmful solvents, the expense of maintaining controlled atmospheric conditions, and the inherent instabilities of PSCs under operation. Here, we address these challenges by introducing a high volatility, low toxicity, biorenewable solvent system to fabricate a range of 2D perovskites, which highly effective precursor phases for subsequent transformation to alpha-formamidinium lead triiodide (FAPbI3), fully processed under ambient conditions. PSCs utilising our FAPbI3 reproducibly show remarkable stability under illumination and elevated temperature (ISOS-L-2) and "damp heat" (ISOS-D-3) stressing, surpassing other state-of-the-art perovskite compositions. We determine that this enhancement is a consequence of the 2D precursor phase crystallisation route, which simultaneously avoids retention of residual low-volatility solvents (such as DMF and DMSO) and reduces the rate of degradation of FA+ in the material. Our findings highlight both the critical role of the initial crystallisation process in determining the operational stability of perovskite materials, and that neat FA+-based perovskites can be competitively stable despite the inherent metastability of the alpha-phase. △ Less

Submitted 14 June, 2024; v1 submitted 12 June, 2024; originally announced June 2024.

MANTA: A Negative-Triangularity NASEM-Compliant Fusion Pilot Plant

Authors: MANTA Collaboration, G. Rutherford, H. S. Wilson, A. Saltzman, D. Arnold, J. L. Ball, S. Benjamin, R. Bielajew, N. de Boucaud, M. Calvo-Carrera, R. Chandra, H. Choudhury, C. Cummings, L. Corsaro, N. DaSilva, R. Diab, A. R. Devitre, S. Ferry, S. J. Frank, C. J. Hansen, J. Jerkins, J. D. Johnson, P. Lunia, J. van de Lindt, S. Mackie , et al. (16 additional authors not shown)

Abstract: The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicin… ▽ More The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicine report ``Bringing Fusion to the U.S. Grid". A self-consistent integrated modeling workflow predicts a fusion power of 450 MW and a plasma gain of 11.5 with only 23.5 MW of power to the scrape-off layer (SOL). This low $P_\text{SOL}$ together with impurity seeding and high density at the separatrix results in a peak heat flux of just 2.8 MW/m$^{2}$. MANTA's high aspect ratio provides space for a large central solenoid (CS), resulting in ${\sim}$15 minute inductive pulses. In spite of the high B fields on the CS and the other REBCO-based magnets, the electromagnetic stresses remain below structural and critical current density limits. Iterative optimization of neutron shielding and tritium breeding blanket yield tritium self-sufficiency with a breeding ratio of 1.15, a blanket power multiplication factor of 1.11, toroidal field coil lifetimes of $3100 \pm 400$ MW-yr, and poloidal field coil lifetimes of at least $890 \pm 40$ MW-yr. Following balance of plant modeling, MANTA is projected to generate 90 MW of net electricity at an electricity gain factor of ${\sim}2.4$. Systems-level economic analysis estimates an overnight cost of US\$3.4 billion, meeting the NASEM FPP requirement that this first-of-a-kind be less than US\$5 billion. The toroidal field coil cost and replacement time are the most critical upfront and lifetime cost drivers, respectively. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Self-interaction of turbulent eddies in tokamaks with low magnetic shear

Authors: Arnas Volčokas, Justin Ball, Stephan Brunner

Abstract: Using local nonlinear gyrokinetic simulations, we demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns in tokamaks with weak or zero magnetic shear $\hat{s}$. We observe that this parallel eddy length scales inversely with magnetic shear and at $\hat{s}=0$ is limited by the thermal speed of electrons $v_{th,e}$. We examine the consequences of these… ▽ More Using local nonlinear gyrokinetic simulations, we demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns in tokamaks with weak or zero magnetic shear $\hat{s}$. We observe that this parallel eddy length scales inversely with magnetic shear and at $\hat{s}=0$ is limited by the thermal speed of electrons $v_{th,e}$. We examine the consequences of these "ultra long" eddies on turbulent transport, in particular, how field line topology mediates strong parallel self-interaction. Our investigation reveals that, through this process, field line topology can strongly affect transport. It can cause transitions between different turbulent instabilities and in some cases triple the logarithmic gradient needed to drive a given amount of heat flux. We also identify a novel "eddy squeezing" effect, which reduces the perpendicular size of eddies and their ability to transport energy, thus representing a novel approach to improve confinement. Finally, we investigate the triggering mechanism of Internal Transport Barriers (ITBs) using low magnetic shear simulations, shedding light on why ITBs are often easier to trigger where the safety factor has a low-order rational value. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Assessing the Risk of Proliferation via Fissile Material Breeding in ARC-class Fusion Power Plants

Authors: J. L. Ball, E. E. Peterson, R. S. Kemp, S. E. Ferry

Abstract: Construction of a nuclear weapon requires access to kilogram-scale quantities of fissile material, which can be bred from fertile material like U-238 and Th-232 via neutron capture. Future fusion power plants, with total neutron source rates in excess of $10^{20}$ n/s, could breed weapons-relevant quantities of fissile material on short timescales, posing a breakout proliferation risk. The ARC-cla… ▽ More Construction of a nuclear weapon requires access to kilogram-scale quantities of fissile material, which can be bred from fertile material like U-238 and Th-232 via neutron capture. Future fusion power plants, with total neutron source rates in excess of $10^{20}$ n/s, could breed weapons-relevant quantities of fissile material on short timescales, posing a breakout proliferation risk. The ARC-class fusion reactor design is characterized by demountable high temperature superconducting magnets, a FLiBe liquid immersion blanket, and a relatively small size ($\sim$ 4 m major radius, $\sim$ 1 m minor radius). We use the open-source Monte Carlo neutronics code OpenMC to perform self-consistent time-dependent simulations of a representative ARC-class blanket to assess the feasibility of a fissile breeding breakout scenario. We find that a significant quantity of fissile material can be bred in less than six months of full power operation for initial fertile inventories ranging from 5 to 50 metric tons, representing a non-negligible proliferation risk. We further study the feasibility of this scenario by examining other consequences of fissile breeding such as reduced tritium breeding ratio, extra heat from fission and decay heat, isotopic purity of bred material, and self-protection time of irradiated blanket material. We also examine the impact of Li-6 enrichment on fissile breeding and find that it substantially reduces breeding rate, motivating its use as a proliferation resistance tool. △ Less

Submitted 5 June, 2024; v1 submitted 18 April, 2024; originally announced April 2024.

Comments: Submitted to Nuclear Fusion

Geometric Burn Control For Tokamaks

Authors: J. F. Parisi, J. W. Berkery, A. Sladkomedova, S. Guizzo, M. R. Hardman, J. R. Ball, A. O. Nelson, S. M. Kaye, M. Anastopoulos-Tzanis, S. A. M. McNamara, J. Dominski, S. Janhunen, M. Romanelli, D. Dickinson, A. Diallo, A. Dnestrovskii, W. Guttenfelder, C. Hansen, O. Myatra, H. R. Wilson

Abstract: A new burn control scheme for tokamaks is described where the total fusion power is controlled by adjusting the plasma volume fraction that is packed into power dense regions. In an example spherical tokamak burning plasma, by modifying the plasma edge squareness the total fusion power is doubled at almost constant total plasma volume and fusion power density. Therefore, increased plasma squarenes… ▽ More A new burn control scheme for tokamaks is described where the total fusion power is controlled by adjusting the plasma volume fraction that is packed into power dense regions. In an example spherical tokamak burning plasma, by modifying the plasma edge squareness the total fusion power is doubled at almost constant total plasma volume and fusion power density. Therefore, increased plasma squareness could be extremely beneficial to a fusion reactor and squareness control could be desirable for power load balancing. Experiments have observed the impact of increased edge squareness on modified core plasma volume, highlighting the practical relevance of this approach. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: 7 pages, 7 figures

Physical insights from the aspect ratio dependence of turbulence in negative triangularity plasmas

Authors: Alessandro Balestri, Justin Ball, Stefano Coda, Diego Jose Cruz-Zabala, Manuel Garcia-Munoz, Eleonora Viezzer

Abstract: In this work, we study the impact of aspect ratio A = R0 /r (the ratio of major radius R0 to minor radius r) on the confinement benefits of Negative Triangularity (NT) plasma shaping. We use high-fidelity flux tube gyrokinetic GENE simulations and consider several different scenarios: four of them inspired by TCV experimental data, a scenario inspired by DIII-D experimental data and a scenario exp… ▽ More In this work, we study the impact of aspect ratio A = R0 /r (the ratio of major radius R0 to minor radius r) on the confinement benefits of Negative Triangularity (NT) plasma shaping. We use high-fidelity flux tube gyrokinetic GENE simulations and consider several different scenarios: four of them inspired by TCV experimental data, a scenario inspired by DIII-D experimental data and a scenario expected in the new SMART spherical tokamak. The present study reveals a surprising and non-trivial dependence. NT improves confinement at any value of A for ITG turbulence, while for TEM turbulence confinement is improved only in the case of large and conventional aspect ratios. Additionally, through a detailed study of a large aspect ratio case with pure ITG drive, we develop an intuitive physical picture that explains the beneficial effect of NT at large and conventional aspect ratios. This picture does not hold in TEM-dominated regimes, where a complex synergistic effect of many factors is found. Finally, we performed the first linear gyrokinetic simulations of SMART, finding that both NT and PT scenarios are dominated by micro-tearing-mode (MTM) turbulence and that NT is more susceptible to MTMs at tight aspect ratio. However, we found that a regime where ITG dominates in SMART can be found, and in this regime NT is more linearly stable. △ Less

Submitted 22 May, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Risk assessment and observation of driver with pedestrian using instantaneous heart rate and HRV

Authors: Riku Kikuta, Daniel Carruth, John Ball, Reuben Burch, Ichiro Kageyama

Abstract: Currently, human drivers outperform self-driving vehicles in many conditions such as collision avoidance. Therefore, understanding human driver behaviour in these conditions will provide insight for future autonomous vehicles. For understanding driver behaviour, risk assessment is applied so far as one of the approaches by using both subjective and objective measurement. Subjective measurement met… ▽ More Currently, human drivers outperform self-driving vehicles in many conditions such as collision avoidance. Therefore, understanding human driver behaviour in these conditions will provide insight for future autonomous vehicles. For understanding driver behaviour, risk assessment is applied so far as one of the approaches by using both subjective and objective measurement. Subjective measurement methods such as questionnaires may provide insight into driver risk assessment but there is often significant variability between drivers.Physiological measurements such as heart rate (HR), electroencephalogram (EEG), and electromyogram (EMG) provide more objective measurements of driver risk assessment. HR is often used for measuring driver risk assessment based on observed correlations between HR and risk perception. Previous work has used HR to measure driver risk assessment in self-driving systems, but pedestrian dynamics is not considered for the research. In this study, we observed driver behaviour in certain scenarios which have pedestrian on driving simulator. The scenarios have safe/unsafe situations (i.e., pedestrian crosses road and vehicle may hit pedestrian in one scenario), HR analysis in time/frequency domain is processed for risk assessment. As a result, HR analysis in frequency domain shows certain reasonability for driver risk assessment when driver has pedestrian in its traffic. △ Less

Submitted 10 February, 2024; originally announced February 2024.

Journal ref: 2023 AHFE Open Access, vol 95

Voice Activity Detection (VAD) in Noisy Environments

Authors: Joshua Ball

Abstract: In the realm of digital audio processing, Voice Activity Detection (VAD) plays a pivotal role in distinguishing speech from non-speech elements, a task that becomes increasingly complex in noisy environments. This paper details the development and implementation of a VAD system, specifically engineered to maintain high accuracy in the presence of various ambient noises. We introduce a novel algori… ▽ More In the realm of digital audio processing, Voice Activity Detection (VAD) plays a pivotal role in distinguishing speech from non-speech elements, a task that becomes increasingly complex in noisy environments. This paper details the development and implementation of a VAD system, specifically engineered to maintain high accuracy in the presence of various ambient noises. We introduce a novel algorithm enhanced with a specially designed filtering technique, effectively isolating speech even amidst diverse background sounds. Our comprehensive testing and validation demonstrate the system's robustness, highlighting its capability to discern speech from noise with remarkable precision. The exploration delves into: (1) the core principles underpinning VAD and its crucial role in modern audio processing; (2) the methodologies we employed to filter ambient noise; and (3) a presentation of evidence affirming our system's superior performance in noisy conditions. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 7 pages

Ultra-Resolution Cascaded Diffusion Model for Gigapixel Image Synthesis in Histopathology

Authors: Sarah Cechnicka, Hadrien Reynaud, James Ball, Naomi Simmonds, Catherine Horsfield, Andrew Smith, Candice Roufosse, Bernhard Kainz

Abstract: Diagnoses from histopathology images rely on information from both high and low resolutions of Whole Slide Images. Ultra-Resolution Cascaded Diffusion Models (URCDMs) allow for the synthesis of high-resolution images that are realistic at all magnification levels, focusing not only on fidelity but also on long-distance spatial coherency. Our model beats existing methods, improving the pFID-50k [2]… ▽ More Diagnoses from histopathology images rely on information from both high and low resolutions of Whole Slide Images. Ultra-Resolution Cascaded Diffusion Models (URCDMs) allow for the synthesis of high-resolution images that are realistic at all magnification levels, focusing not only on fidelity but also on long-distance spatial coherency. Our model beats existing methods, improving the pFID-50k [2] score by 110.63 to 39.52 pFID-50k. Additionally, a human expert evaluation study was performed, reaching a weighted Mean Absolute Error (MAE) of 0.11 for the Lower Resolution Diffusion Models and a weighted MAE of 0.22 for the URCDM. △ Less

Submitted 2 December, 2023; originally announced December 2023.

Comments: MedNeurIPS 2023 poster

A new quasilinear model for turbulent momentum transport in tokamaks with flow shear and plasma shaping

Authors: Haomin Sun, Justin Ball, Stephan Brunner, Arnas Volčokas

Abstract: In tokamak experiments, sufficiently strong $E\times B$ flow shear reduces turbulent transport, thereby improving the prospects for fusion power plants. It is therefore of great importance to efficiently explore parameter space to find where strong plasma flow can be achieved. To this end, we propose a new, physically motivated quasi-linear model for estimating momentum transport from turbulence i… ▽ More In tokamak experiments, sufficiently strong $E\times B$ flow shear reduces turbulent transport, thereby improving the prospects for fusion power plants. It is therefore of great importance to efficiently explore parameter space to find where strong plasma flow can be achieved. To this end, we propose a new, physically motivated quasi-linear model for estimating momentum transport from turbulence in the presence of toroidal flow shear and plasma shaping. The method gives good estimates of momentum transport for up-down asymmetric geometries as well as low magnetic shear and tight aspect ratio. The results are benchmarked with high-fidelity nonlinear GENE simulations, demonstrating that it provides a fast and accurate estimate of momentum transport. △ Less

Submitted 29 January, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 34 pages, 17 figures

Grain-level effects on in-situ deformation-induced phase transformations in a complex-phase steel using 3DXRD and EBSD

Authors: James A. D. Ball, Claire Davis, Carl Slater, Himanshu Vashishtha, Mohammed Said, Louis Hébrard, Florian Steinhilber, Jonathan P. Wright, Thomas Connolley, Stefan Michalik, David M. Collins

Abstract: A novel complex-phase steel alloy is conceived with a deliberately unstable austenite, $γ$, phase that enables the deformation-induced martensitic transformations (DIMT) to be explored at low levels of plastic strain. The DIMT was thus explored, in-situ and non-destructively, using both far-field Three-Dimensional X-Ray Diffraction (3DXRD) and Electron Back-Scatter Diffraction (EBSD). Substantial… ▽ More A novel complex-phase steel alloy is conceived with a deliberately unstable austenite, $γ$, phase that enables the deformation-induced martensitic transformations (DIMT) to be explored at low levels of plastic strain. The DIMT was thus explored, in-situ and non-destructively, using both far-field Three-Dimensional X-Ray Diffraction (3DXRD) and Electron Back-Scatter Diffraction (EBSD). Substantial $α'$ martensite formation was observed under 10% applied strain with EBSD, and many $\varepsilon$ grain formation events were captured with 3DXRD, indicative of the indirect transformation of martensite via the reaction $γ\rightarrow \varepsilon \rightarrow α'$. Using $\varepsilon$ grain formation as a direct measurement of $γ$ grain stability, the influence of several microstructural properties, such as grain size, orientation and neighbourhood configuration, on $γ$ stability have been identified. Larger $γ$ grains were found to be less stable than smaller grains. Any $γ$ grains oriented with {100} parallel to the loading direction preferentially transformed with lower stresses. Parent $\varepsilon$-forming $γ$ grains possessed a neighbourhood with increased ferritic/martensitic volume fraction. This finding shows, unambiguously, that $α$/$α'$ promotes $\varepsilon$ formation in neighbouring grains. The minimum strain work criterion model for $\varepsilon$ variant prediction was also evaluated, which worked well for most grains. However, $\varepsilon$-forming grains with a lower stress were less well predicted by the model, indicating crystal-level behaviour must be considered for accurate $\varepsilon$ formation. The findings from this work are considered key for the future design of alloys where the deformation response can be controlled by tailoring microstructure and local or macroscopic crystal orientations. △ Less

Submitted 7 November, 2023; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: 20 pages, 3 supplementary pages, 14 figures, 3 supplementary figures. Preprint submitted to Acta Materialia. Updated to correct Figure 14

Journal ref: Acta Materialia 265 (2024) 119608

Extension of convex functions from a hyperplane to a half-space

Authors: John M. Ball, Christopher L. Horner

Abstract: It is shown that a possibly infinite-valued proper lower semicontinuous convex function on ${\mathbb R}^n$ has an extension to a convex function on the half-space ${\mathbb R}^n\times[0,\infty)$ which is finite and smooth on the open half-space ${\mathbb R}^n\times(0,\infty)$. The result is applied to nonlinear elasticity, where it clarifies how the condition of polyconvexity of the free-energy de… ▽ More It is shown that a possibly infinite-valued proper lower semicontinuous convex function on ${\mathbb R}^n$ has an extension to a convex function on the half-space ${\mathbb R}^n\times[0,\infty)$ which is finite and smooth on the open half-space ${\mathbb R}^n\times(0,\infty)$. The result is applied to nonlinear elasticity, where it clarifies how the condition of polyconvexity of the free-energy density $ψ(Dy)$ is best expressed when $ψ(A)\to\infty$ as $\det A\to 0+$. △ Less

Submitted 29 March, 2024; v1 submitted 5 November, 2023; originally announced November 2023.

Comments: To appear in Calc Var PDE

MSC Class: 26B25 49J45 74B20

Dilation and Model Theory for Pairs of Commuting Contractions

Authors: Joseph A. Ball, Haripada Sau

Abstract: This manuscript is an effort to extend the Sz.-Nagy--Foias dilation and model theory for a single contraction to the case of commuting pair of contractions. Fundamental to the Sz.-Nagy--Foias model theory is the functional model for the minimal isometric dilation. The first step in our approach for the pair case is to obtain further information, beyond that in the original paper of Ando, concernin… ▽ More This manuscript is an effort to extend the Sz.-Nagy--Foias dilation and model theory for a single contraction to the case of commuting pair of contractions. Fundamental to the Sz.-Nagy--Foias model theory is the functional model for the minimal isometric dilation. The first step in our approach for the pair case is to obtain further information, beyond that in the original paper of Ando, concerning the structure of the plethora of minimal commuting isometric lifts. We exhibit an explicit simple example of two minimal isometric lifts of a commuting contractive pair that are not unitarily equivalent -- see Chapter 5. We provide two constructive new proofs of Ando's Dilation Theorem, each of which leads to a new functional-model representation for such a lift -- see Theorem 4.3.8 and Remark 4.5.7. The construction leads to the identification of a set of additional free parameters which serves to classify the distinct unitary-equivalence classes of minimal Ando lifts. However this lack of uniqueness limits the utility of such minimal Ando lifts for the construction of a functional model for a commuting contractive pair. We identify an intermediate type of lift, called pseudo-commuting contractive lift, which paves the way for a functional model. In the model form, the Sz.-Nagy--Foias characteristic function is augmented by what is called the fundamental operator pair, together with a canonical pair of commuting unitary operators, so that the augmented collection, called the characteristic triple, is a complete unitary invariant for a commuting contractive pair. There is also a notion of admissible triple as the substitute for a purely contractive analytic function in the Sz.-Nagy--Foias theory, from which one can construct a functional model commuting contractive pair having its characteristic triple coinciding with the original admissible triple in an appropriate sense. △ Less

Submitted 15 August, 2023; originally announced August 2023.

Comments: 150 pages

Remarks on the linear wave equation

Authors: John M. Ball

Abstract: We make some remarks on the linear wave equation concerning the existence and uniqueness of weak solutions, satisfaction of the energy equation, growth properties of solutions, the passage from bounded to unbounded domains, and reconciliation of different representations of solutions. We make some remarks on the linear wave equation concerning the existence and uniqueness of weak solutions, satisfaction of the energy equation, growth properties of solutions, the passage from bounded to unbounded domains, and reconciliation of different representations of solutions. △ Less

Submitted 14 September, 2023; v1 submitted 31 July, 2023; originally announced July 2023.

Comments: Some changes to section on growth of L^2 norm, plus minor corrections. To appear in Quarterly of Applied Mathematics

MSC Class: 35L05 47D06 42B37

13 New Light Curves and Updated Mid-Transit Time and Period for Hot Jupiter WASP-104 b with EXOTIC

Authors: Heather B. Hewitt, Federico Noguer, Suber Corley, James Ball, Claudia Chastain, Richard Cochran-White, Kendall Collins, Kris Ganzel, Kimberly Merriam Gray, Mike Logan, Steve Marquez-Perez, Chyna Merchant, Matthew Pedone, Gina Plumey, Matthew Rice, Zachary Ruybal, Molly N. Simon, Isabela Huckabee, Robert T. Zellem, Kyle A. Pearson

Abstract: Using the EXOplanet Transit Interpretation Code (EXOTIC), we reduced 52 sets of images of WASP-104 b, a Hot Jupiter-class exoplanet orbiting WASP-104, in order to obtain an updated mid-transit time (ephemeris) and orbital period for the planet. We performed this reduction on images taken with a 6-inch telescope of the Center for Astrophysics | Harvard & Smithsonian MicroObservatory. Of the reduced… ▽ More Using the EXOplanet Transit Interpretation Code (EXOTIC), we reduced 52 sets of images of WASP-104 b, a Hot Jupiter-class exoplanet orbiting WASP-104, in order to obtain an updated mid-transit time (ephemeris) and orbital period for the planet. We performed this reduction on images taken with a 6-inch telescope of the Center for Astrophysics | Harvard & Smithsonian MicroObservatory. Of the reduced light curves, 13 were of sufficient accuracy to be used in updating the ephemerides for WASP-104 b, meeting or exceeding the three-sigma standard for determining a significant detection. Our final mid-transit value was 2457805.170208 +/- 0.000036 BJD_TBD and the final period value was 1.75540644 +/- 0.00000016 days. The true significance of our results is in their derivation from image sets gathered over time by a small, ground-based telescope as part of the Exoplanet Watch citizen science initiative, and their competitive results to an ephemeris generated from data gathered by the TESS telescope. We use these results to further show how such techniques can be employed by amateur astronomers and citizen scientists to maximize the efficacy of larger telescopes by reducing the use of expensive observation time. The work done in the paper was accomplished as part of the first fully online Course-Based Undergraduate Research Experience (CURE) for astronomy majors in the only online Bachelor of Science program in Astronomical and Planetary Sciences. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 6 pages, 5 figures, published in JAAVSO

Journal ref: Journal of the American Association of Variable Star Observers, 51(1) (2023)

Registration between DCT and EBSD datasets for multiphase microstructures

Authors: James A. D. Ball, Jette Oddershede, Claire Davis, Carl Slater, Mohammed Said, Himanshu Vashishtha, Stefan Michalik, David M. Collins

Abstract: The ability to characterise the three-dimensional microstructure of multiphase materials is essential for understanding the interaction between phases and associated materials properties. Here, laboratory-based diffraction-contrast tomography (DCT), a recently-established materials characterization technique that can determine grain phases, morphologies, positions and orientations in a voxel-based… ▽ More The ability to characterise the three-dimensional microstructure of multiphase materials is essential for understanding the interaction between phases and associated materials properties. Here, laboratory-based diffraction-contrast tomography (DCT), a recently-established materials characterization technique that can determine grain phases, morphologies, positions and orientations in a voxel-based reconstruction method, was used to map part of a dual-phase steel alloy sample. To assess the resulting microstructures that were produced by the DCT technique, an EBSD map was collected within the same sample volume. To identify the 2D slice of the 3D DCT reconstruction that best corresponded to the EBSD map, a novel registration technique based solely on grain-averaged orientations was developed -- this registration technique requires very little a priori knowledge of dataset alignment and can be extended to other techniques that only recover grain-averaged orientation data such as far-field 3D X-ray diffraction microscopy. Once the corresponding 2D slice was identified in the DCT dataset, comparisons of phase balance, grain size, shape and texture were performed between DCT and EBSD techniques. More complicated aspects of the microstructural morphology such as grain boundary shape and grains less than a critical size were poorly reproduced by the DCT reconstruction, primarily due to the difference in resolutions of the technique compared with EBSD. However, lab-based DCT is shown to accurately determine the centre-of-mass position, orientation, and size of the large grains for each phase present, austenite and martensitic ferrite. The results reveals a complex ferrite grain network of similar crystal orientations that are absent from the EBSD dataset. Such detail demonstrates that lab-based DCT, as a technique, shows great promise in the field of multi-phase material characterization. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: 15 pages, 11 figures. Preprint submitted to Materials Characterization

Journal ref: Materials Characterization, October 2023, Volume 204, Page 113228

Realistic Data Enrichment for Robust Image Segmentation in Histopathology

Authors: Sarah Cechnicka, James Ball, Hadrien Reynaud, Callum Arthurs, Candice Roufosse, Bernhard Kainz

Abstract: Poor performance of quantitative analysis in histopathological Whole Slide Images (WSI) has been a significant obstacle in clinical practice. Annotating large-scale WSIs manually is a demanding and time-consuming task, unlikely to yield the expected results when used for fully supervised learning systems. Rarely observed disease patterns and large differences in object scales are difficult to mode… ▽ More Poor performance of quantitative analysis in histopathological Whole Slide Images (WSI) has been a significant obstacle in clinical practice. Annotating large-scale WSIs manually is a demanding and time-consuming task, unlikely to yield the expected results when used for fully supervised learning systems. Rarely observed disease patterns and large differences in object scales are difficult to model through conventional patient intake. Prior methods either fall back to direct disease classification, which only requires learning a few factors per image, or report on average image segmentation performance, which is highly biased towards majority observations. Geometric image augmentation is commonly used to improve robustness for average case predictions and to enrich limited datasets. So far no method provided sampling of a realistic posterior distribution to improve stability, e.g. for the segmentation of imbalanced objects within images. Therefore, we propose a new approach, based on diffusion models, which can enrich an imbalanced dataset with plausible examples from underrepresented groups by conditioning on segmentation maps. Our method can simply expand limited clinical datasets making them suitable to train machine learning pipelines, and provides an interpretable and human-controllable way of generating histopathology images that are indistinguishable from real ones to human experts. We validate our findings on two datasets, one from the public domain and one from a Kidney Transplant study. △ Less

Submitted 7 August, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: 11 pages, 2 figures, 1 table

The two rings of (50000) Quaoar

Authors: C. L. Pereira, B. Sicardy, B. E. Morgado, F. Braga-Ribas, E. Fernández-Valenzuela, D. Souami, B. J. Holler, R. C. Boufleur, G. Margoti, M. Assafin, J. L. Ortiz, P. Santos-Sanz, B. Epinat, P. Kervella, J. Desmars, R. Vieira-Martins, Y. Kilic, A. R. Gomes-Júnior, J. I. B. Camargo, M. Emilio, M. Vara-Lubiano, M. Kretlow, L. Albert, C. Alcock, J. G. Ball , et al. (44 additional authors not shown)

Abstract: Quaoar is a classical Trans-Neptunian Object (TNO) with an area equivalent diameter of 1,100 km and an orbital semi-major axis of 43.3 astronomical units. Based on stellar occultations observed between 2018 and 2021, an inhomogeneous ring (Q1R, Quaoar's first ring) was detected around this body. Aims. A new stellar occultation by Quaoar was observed on August 9th, 2022 aiming to improve Quaoar's s… ▽ More Quaoar is a classical Trans-Neptunian Object (TNO) with an area equivalent diameter of 1,100 km and an orbital semi-major axis of 43.3 astronomical units. Based on stellar occultations observed between 2018 and 2021, an inhomogeneous ring (Q1R, Quaoar's first ring) was detected around this body. Aims. A new stellar occultation by Quaoar was observed on August 9th, 2022 aiming to improve Quaoar's shape models and the physical parameters of Q1R while searching for additional material around the body. Methods. The occultation provided nine effective chords across Quaoar, pinning down its size, shape, and astrometric position. Large facilities, such as Gemini North and the Canada-France-Hawaii Telescope (CFHT), were used to obtain high acquisition rates and signal-to-noise ratios. The light curves were also used to characterize the Q1R ring (radial profiles and orbital elements). Results. Quaoar's elliptical fit to the occultation chords yields the limb with an apparent semi-major axis of $579.5\pm4.0$ km, apparent oblateness of $0.12\pm0.01$, and area-equivalent radius of $543\pm2$ km. Quaoar's limb orientation is consistent with Q1R and Weywot orbiting in Quaoar's equatorial plane. The orbital radius of Q1R is refined to a value of $4,057\pm6$ km. The radial opacity profile of the more opaque ring profile follows a Lorentzian shape that extends over 60 km, with a full width at half maximum (FWHM) of $\sim5$ km and a peak normal optical depth of 0.4. Besides the secondary events related to the already reported rings, new secondary events detected during the August 2022 occultation in three different data sets are consistent with another ring around Quaoar with a radius of $2,520\pm20$ km, assuming the ring is circular and co-planar with Q1R. This new ring has a typical width of 10 km and a normal optical depth of $\sim$0.004. Like Q1R, it also lies outside Quaoar's classical Roche limit. △ Less

Submitted 20 April, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: Accepted for publication in Astronomy & Astrophysics (17-April-2023). 18 pages, 12 figures

Journal ref: A&A 673, L4 (2023)

Image comparison and scaling via nonlinear elasticity

Authors: John M. Ball, Christopher L. Horner

Abstract: A nonlinear elasticity model for comparing images is formulated and analyzed, in which optimal transformations between images are sought as minimizers of an integral functional. The existence of minimizers in a suitable class of homeomorphisms between image domains is established under natural hypotheses. We investigate whether for linearly related images the minimization algorithm delivers the li… ▽ More A nonlinear elasticity model for comparing images is formulated and analyzed, in which optimal transformations between images are sought as minimizers of an integral functional. The existence of minimizers in a suitable class of homeomorphisms between image domains is established under natural hypotheses. We investigate whether for linearly related images the minimization algorithm delivers the linear transformation as the unique minimizer. △ Less

Submitted 20 March, 2023; v1 submitted 17 March, 2023; originally announced March 2023.

Comments: SSVM2023 Proceedings to appear. New references added plus related minor changes

MSC Class: 65D19 74B20 49J10

Synthetic Experience Replay

Authors: Cong Lu, Philip J. Ball, Yee Whye Teh, Jack Parker-Holder

Abstract: A key theme in the past decade has been that when large neural networks and large datasets combine they can produce remarkable results. In deep reinforcement learning (RL), this paradigm is commonly made possible through experience replay, whereby a dataset of past experiences is used to train a policy or value function. However, unlike in supervised or self-supervised learning, an RL agent has to… ▽ More A key theme in the past decade has been that when large neural networks and large datasets combine they can produce remarkable results. In deep reinforcement learning (RL), this paradigm is commonly made possible through experience replay, whereby a dataset of past experiences is used to train a policy or value function. However, unlike in supervised or self-supervised learning, an RL agent has to collect its own data, which is often limited. Thus, it is challenging to reap the benefits of deep learning, and even small neural networks can overfit at the start of training. In this work, we leverage the tremendous recent progress in generative modeling and propose Synthetic Experience Replay (SynthER), a diffusion-based approach to flexibly upsample an agent's collected experience. We show that SynthER is an effective method for training RL agents across offline and online settings, in both proprioceptive and pixel-based environments. In offline settings, we observe drastic improvements when upsampling small offline datasets and see that additional synthetic data also allows us to effectively train larger networks. Furthermore, SynthER enables online agents to train with a much higher update-to-data ratio than before, leading to a significant increase in sample efficiency, without any algorithmic changes. We believe that synthetic training data could open the door to realizing the full potential of deep learning for replay-based RL algorithms from limited data. Finally, we open-source our code at https://github.com/conglu1997/SynthER. △ Less

Submitted 26 October, 2023; v1 submitted 12 March, 2023; originally announced March 2023.

Comments: Published at NeurIPS, 2023

Slip and twinning in Bravais lattices

Authors: John M. Ball

Abstract: A unified treatment of slip and twinning in Bravais lattices is given, focussing on the case of cubic symmetry, and using the Ericksen energy well formulation, so that interfaces correspond to rank-one connections between the infinitely many crystallographically equivalent energy wells. Twins are defined to be such rank-one connections involving a nontrivial reflection of the lattice across some p… ▽ More A unified treatment of slip and twinning in Bravais lattices is given, focussing on the case of cubic symmetry, and using the Ericksen energy well formulation, so that interfaces correspond to rank-one connections between the infinitely many crystallographically equivalent energy wells. Twins are defined to be such rank-one connections involving a nontrivial reflection of the lattice across some plane. The slips and twins minimizing shear magnitude for cubic lattices are rigorously calculated, and the conjugates of these and other slips analyzed. It is observed that all rank-one connections between the energy wells for the dual of a Bravais lattice can be obtained explicitly from those for the original lattice, so that in particular the rank-one connections for fcc can be obtained explicitly from those for bcc. △ Less

Submitted 17 August, 2023; v1 submitted 8 March, 2023; originally announced March 2023.

Comments: Various corrections and clarifications made. The result that the rank-one connections for fcc can be obtained from those for bcc (and vice versa) now generalized to apply to any pair of dual lattices. To appear in Journal of Elasticity

MSC Class: 74N15; 74G65

Per-grain and neighbourhood stress interactions during deformation of a ferritic steel obtained using three-dimensional X-ray diffraction

Authors: James A. D. Ball, Anna Kareer, Oxana V. Magdysyuk, Stefan Michalik, Thomas Connolley, David M. Collins

Abstract: Three-dimensional X-ray diffraction (3DXRD) has been used to measure, in-situ, the evolution of $\sim 1800$ grains in a single phase low carbon ferritic steel sample during uniaxial deformation. The distribution of initial residual grain stresses in the material was observed to prevail as plasticity builds, though became less pronounced, and therefore less influential as strain increased. The init… ▽ More Three-dimensional X-ray diffraction (3DXRD) has been used to measure, in-situ, the evolution of $\sim 1800$ grains in a single phase low carbon ferritic steel sample during uniaxial deformation. The distribution of initial residual grain stresses in the material was observed to prevail as plasticity builds, though became less pronounced, and therefore less influential as strain increased. The initial Schmid factor of a grain was found to be strongly correlated to the intergranular stress change and the range of stresses that are permissible; a grain well aligned for easy slip is more likely to exhibit a range of stresses than those orientated poorly for dislocation motion. The orientation path of a grain, however, is not only dependent on its initial orientation, but hypothesised to be influenced by its stress state and the stress state of its grain environment. A grain neighbourhood effect is observed: the Schmid factor of serial adjoining grains influences the stress state of a grain of interest, whereas parallel neighbours are much less influential. This phenomenon is strongest at low plastic strains only, with the effect diminishing as plasticity builds. The influence of initial residual stresses becomes less evident, and grains rotate to eliminate any orientation dependent load shedding. The ability of the BCC ferrite to exhaust such neighbourhood interactions, which would otherwise be detrimental in crystal structures with lower symmetric and fewer slip systems, is considered key to the high ductility possessed by these materials. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 17 pages (+2 supplementary pages), 13 figures (+3 supplementary figures). Preprint submitted to Acta Materialia

Journal ref: Commun Mater 5, 27 (2024)

Efficient Online Reinforcement Learning with Offline Data

Authors: Philip J. Ball, Laura Smith, Ilya Kostrikov, Sergey Levine

Abstract: Sample efficiency and exploration remain major challenges in online reinforcement learning (RL). A powerful approach that can be applied to address these issues is the inclusion of offline data, such as prior trajectories from a human expert or a sub-optimal exploration policy. Previous methods have relied on extensive modifications and additional complexity to ensure the effective use of this dat… ▽ More Sample efficiency and exploration remain major challenges in online reinforcement learning (RL). A powerful approach that can be applied to address these issues is the inclusion of offline data, such as prior trajectories from a human expert or a sub-optimal exploration policy. Previous methods have relied on extensive modifications and additional complexity to ensure the effective use of this data. Instead, we ask: can we simply apply existing off-policy methods to leverage offline data when learning online? In this work, we demonstrate that the answer is yes; however, a set of minimal but important changes to existing off-policy RL algorithms are required to achieve reliable performance. We extensively ablate these design choices, demonstrating the key factors that most affect performance, and arrive at a set of recommendations that practitioners can readily apply, whether their data comprise a small number of expert demonstrations or large volumes of sub-optimal trajectories. We see that correct application of these simple recommendations can provide a $\mathbf{2.5\times}$ improvement over existing approaches across a diverse set of competitive benchmarks, with no additional computational overhead. We have released our code at https://github.com/ikostrikov/rlpd. △ Less

Submitted 31 May, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: Short Presentation at ICML 2023; to reproduce our results and use our codebase, see https://github.com/ikostrikov/rlpd

Simulations of Precessing Jets and the Formation of X-shaped Radio Galaxies

Authors: Chris Nolting, Jay Ball, Tri M. Nguyen

Abstract: Jet precession is sometimes invoked to explain asymmetries in radio galaxy (RG) jets and "X/S/Z-shape" radio galaxies, caused by the presence of a binary black hole companion to the source active galactic nucleus (AGN) or by accretion instabilities. We present a series of simulations of radio galaxy jet precession to examine how these sources would evolve over time, including a passive distributio… ▽ More Jet precession is sometimes invoked to explain asymmetries in radio galaxy (RG) jets and "X/S/Z-shape" radio galaxies, caused by the presence of a binary black hole companion to the source active galactic nucleus (AGN) or by accretion instabilities. We present a series of simulations of radio galaxy jet precession to examine how these sources would evolve over time, including a passive distribution of cosmic ray electrons (CRe) so we can model radio synchrotron emissions and create synthetic radio maps of the sources. We find that a single source viewed from different angles can result in differing RG morphological classifications, confusing physical implications of these classifications. Additionally, the jet trajectories can become unstable due to their own self-interactions and lead to "reorientation events" that may look like the effects of external dynamics such as shocks, winds, or cold fronts in the medium. Finally, something akin to an "Odd Radio Circle" may be observed in the case of viewing the radio remnant of such a precessing source from a line of sight near the precession axis. △ Less

Submitted 5 May, 2023; v1 submitted 11 January, 2023; originally announced January 2023.

Comments: 12 pages, 5 figures, published in ApJ. Animated versions of the figures can be found at: https://www.youtube.com/playlist?list=PLNl1xrpDmrjb0pcTUwSrTcgmpZrSYx81I

Journal ref: ApJ 948 25 (2023)

A Generalist, Automated ALFALFA Baryonic Tully-Fisher Relation

Authors: Catie J. Ball, Martha P. Haynes, Michael G. Jones, Bo Peng, Adriana Durbala, Rebecca A. Koopmann, Joseph Ribaudo, Aileen O'Donoghue

Abstract: The Baryonic Tully-Fisher Relation (BTFR) has applications in galaxy evolution as a testbed for the galaxy-halo connection and in observational cosmology as a redshift-independent secondary distance indicator. We use the 31,000+ galaxy ALFALFA sample -- which provides redshifts, velocity widths, and HI content for a large number of gas-bearing galaxies in the local universe -- to fit and test an e… ▽ More The Baryonic Tully-Fisher Relation (BTFR) has applications in galaxy evolution as a testbed for the galaxy-halo connection and in observational cosmology as a redshift-independent secondary distance indicator. We use the 31,000+ galaxy ALFALFA sample -- which provides redshifts, velocity widths, and HI content for a large number of gas-bearing galaxies in the local universe -- to fit and test an extensive local universe BTFR. This BTFR is designed to be as inclusive of ALFALFA and comparable samples as possible. Velocity widths measured via an automated method and $M_{b}$ proxies extracted from survey data can be uniformly and efficiently measured for other samples, giving this analysis broad applicability. We also investigate the role of sample demographics in determining the best-fit relation. We find that the best-fit relations are changed significantly by changes to the sample mass range and to second order, mass sampling, gas fraction, different stellar mass and velocity width measurements. We use a subset of ALFALFA with demographics that reflect the full sample to measure a robust BTFR slope of $3.30\pm0.06$. We apply this relation and estimate source distances, finding general agreement with flow-model distances as well as average distance uncertainties of $\sim0.17$ dex for the full ALFALFA sample. We demonstrate the utility of these distance estimates by applying them to a sample of sources in the Virgo vicinity, recovering signatures of infall consistent with previous work. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 23 pages, 10 figures, submitted to ApJ

Parametrized Asymmetric Neutral hydrogen Disk Integrated Spectrum Characterization (PANDISC) I: Introduction to A Physically Motivated H I Model

Authors: Bo Peng, Martha P. Haynes, Catie J. Ball, Michael G. Jones

Abstract: Modelling the integrated H I spectra of galaxies has been a difficult task due to their diverse shapes, but more dynamical information is waiting to be explored in Hi line profiles. Based on simple assumptions, we construct a physically motivated model for the integrated Hi spectra: Parametrized Asymmetric Neutral hydrogen Disk Integrated Spectrum Characterization (PANDISC). The model shows great… ▽ More Modelling the integrated H I spectra of galaxies has been a difficult task due to their diverse shapes, but more dynamical information is waiting to be explored in Hi line profiles. Based on simple assumptions, we construct a physically motivated model for the integrated Hi spectra: Parametrized Asymmetric Neutral hydrogen Disk Integrated Spectrum Characterization (PANDISC). The model shows great flexibility in reproducing the diverse Hi profiles. We use Monte-Carlo Markov Chain (MCMC) for fitting the model to global H I profiles, producing statistically robust quantitative results. Comparing with several samples of H I data available in the literature , we find the model-fitted results agree with catalogued velocity widths (e.g., W50) down to the lowest S/N. The model is also shown to be useful for applications like the baryonic Tully-Fisher relation (BTFR) and profile-based sample control. By comparing v_r to v_flat , we uncover how the H I width is affected by the structure of the rotation curve, following a trend consistent with the difference in the BTFR slope. We also select a sample of spectra with broad wing-like features suggestive of a population of galaxies with unusual gas dynamics. The PANDISC model bears both promise and limitations for potential use beyond H I lines. Further application on the whole ALFALFA sample will enable us to perform large scale ensemble studies of the H I properties and dynamics in nearby galaxies. △ Less

Submitted 15 May, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 28 pages, 18 figures, accepted by ApJ

A free discontinuity model for smectic thin films

Authors: John M. Ball, Giacomo Canevari, Bianca Stroffolini

Abstract: We attempt to describe surface defects in smectic A thin films by formulating a free discontinuity problem - that is, a variational problem in which the order parameter is allowed to have jump discontinuities on some (unknown) set. The free energy functional contains an interfacial energy which penalizes dislocations of the smectic layers at the jump. We discuss mathematical issues related to the… ▽ More We attempt to describe surface defects in smectic A thin films by formulating a free discontinuity problem - that is, a variational problem in which the order parameter is allowed to have jump discontinuities on some (unknown) set. The free energy functional contains an interfacial energy which penalizes dislocations of the smectic layers at the jump. We discuss mathematical issues related to the existence of minimizers and provide examples of minimizers in some simplified settings. △ Less

Submitted 9 March, 2023; v1 submitted 20 November, 2022; originally announced November 2022.

Learning General World Models in a Handful of Reward-Free Deployments

Authors: Yingchen Xu, Jack Parker-Holder, Aldo Pacchiano, Philip J. Ball, Oleh Rybkin, Stephen J. Roberts, Tim Rocktäschel, Edward Grefenstette

Abstract: Building generally capable agents is a grand challenge for deep reinforcement learning (RL). To approach this challenge practically, we outline two key desiderata: 1) to facilitate generalization, exploration should be task agnostic; 2) to facilitate scalability, exploration policies should collect large quantities of data without costly centralized retraining. Combining these two properties, we i… ▽ More Building generally capable agents is a grand challenge for deep reinforcement learning (RL). To approach this challenge practically, we outline two key desiderata: 1) to facilitate generalization, exploration should be task agnostic; 2) to facilitate scalability, exploration policies should collect large quantities of data without costly centralized retraining. Combining these two properties, we introduce the reward-free deployment efficiency setting, a new paradigm for RL research. We then present CASCADE, a novel approach for self-supervised exploration in this new setting. CASCADE seeks to learn a world model by collecting data with a population of agents, using an information theoretic objective inspired by Bayesian Active Learning. CASCADE achieves this by specifically maximizing the diversity of trajectories sampled by the population through a novel cascading objective. We provide theoretical intuition for CASCADE which we show in a tabular setting improves upon naïve approaches that do not account for population diversity. We then demonstrate that CASCADE collects diverse task-agnostic datasets and learns agents that generalize zero-shot to novel, unseen downstream tasks on Atari, MiniGrid, Crafter and the DM Control Suite. Code and videos are available at https://ycxuyingchen.github.io/cascade/ △ Less

Submitted 23 October, 2022; originally announced October 2022.

Comments: To be published at NeurIPS 2022. Code and videos available at https://ycxuyingchen.github.io/cascade/

Gemini North Adaptive Optics (GNAO) facility overview and status updates

Authors: Gaetano Sivo, Julia Scharwächter, Manuel Lazo, Célia Blain, Stephen Goodsell, Marcos van Dam, Martin Tschimmel, Henry Roe, Jennifer Lotz, Kim Tomassino-Reed, William Rambold, Courtney Raich, Ricardo Cardenes, Angelic Ebbers, Tim Gaggstatter, Pedro Gigoux, Thomas Schneider, Charles Cavedoni, Stacy Kang, Stanislas Karewicz, Heather Carr, Jesse Ball, Paul Hirst, Emmanuel Chirre, John White , et al. (32 additional authors not shown)

Abstract: The Gemini North Adaptive Optics (GNAO) facility is the upcoming AO facility for Gemini North providing a state-of-the-art AO system for surveys and time domain science in the era of JWST and Rubin operations. GNAO will be optimized to feed the Gemini infrared Multi Object Spectrograph (GIRMOS). While GIRMOS is the primary science driver for defining the capabilities of GNAO, any instrument oper… ▽ More The Gemini North Adaptive Optics (GNAO) facility is the upcoming AO facility for Gemini North providing a state-of-the-art AO system for surveys and time domain science in the era of JWST and Rubin operations. GNAO will be optimized to feed the Gemini infrared Multi Object Spectrograph (GIRMOS). While GIRMOS is the primary science driver for defining the capabilities of GNAO, any instrument operating with an f/32 beam can be deployed using GNAO. The GNAO project includes the development of a new laser guide star facility which will consist of four side-launched laser beams supporting the two primary AO modes of GNAO: a wide-field mode providing an improved image quality over natural seeing for a 2-arcminute circular field-of-view and a narrow-field mode providing near diffraction-limited performance over a 20x20 arcsecond square field-of-view. The GNAO wide field mode will enable GIRMOS's multi-IFU configuration in which the science beam to each individual IFU will be additionally corrected using multi-object AO within GIRMOS. The GNAO narrow field mode will feed the GIRMOS tiled IFU configuration in which all IFUs are combined into a "super"-IFU in the center of the field. GNAO also includes the development of a new Real Time Controller, a new GNAO Facility System Controller and finally the development of a new AO Bench. We present in this paper an overview of the GNAO facility and provide a status update of each product. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: SPIE conference 2022 Montreal

Local gyrokinetic simulations of tokamaks with non-uniform magnetic shear

Authors: Justin Ball, Stephan Brunner

Abstract: In this work, we modify the standard flux tube simulation domain to include arbitrary ion gyroradius-scale variation in the radial profile of the safety factor. To determine how to appropriately include such a modification, we add a strong ion gyroradius-scale source (inspired by electron cyclotron current drive) to the Fokker-Planck equation, then perform a multi-scale analysis that distinguishes… ▽ More In this work, we modify the standard flux tube simulation domain to include arbitrary ion gyroradius-scale variation in the radial profile of the safety factor. To determine how to appropriately include such a modification, we add a strong ion gyroradius-scale source (inspired by electron cyclotron current drive) to the Fokker-Planck equation, then perform a multi-scale analysis that distinguishes the fast electrons driven by the source from the slow bulk thermal electrons. This allows us to systematically derive the needed changes to the gyrokinetic model. We find new terms that adjust the ion and electron parallel streaming to be along the modified field lines. These terms have been successfully implemented in a gyrokinetic code (while retaining the typical Fourier representation), which enables flux tube studies of non-monotonic safety factor profiles and the associated profile shearing. As an illustrative example, we investigate tokamaks with positive versus negative triangularity plasma shaping and find that the importance of profile shearing is not significantly affected by the change in shape. △ Less

Submitted 13 December, 2022; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 33 pages, 4 figures

Journal ref: Plasma Phys. Control. Fusion 65 (2023) 014004

Tree species classification from hyperspectral data using graph-regularized neural networks

Authors: Debmita Bandyopadhyay, Subhadip Mukherjee, James Ball, Grégoire Vincent, David A. Coomes, Carola-Bibiane Schönlieb

Abstract: We propose a novel graph-regularized neural network (GRNN) algorithm for tree species classification. The proposed algorithm encompasses superpixel-based segmentation for graph construction, a pixel-wise neural network classifier, and the label propagation technique to generate an accurate and realistic (emulating tree crowns) classification map on a sparsely annotated data set. GRNN outperforms s… ▽ More We propose a novel graph-regularized neural network (GRNN) algorithm for tree species classification. The proposed algorithm encompasses superpixel-based segmentation for graph construction, a pixel-wise neural network classifier, and the label propagation technique to generate an accurate and realistic (emulating tree crowns) classification map on a sparsely annotated data set. GRNN outperforms several state-of-the-art techniques not only for the standard Indian Pines HSI but also achieves a high classification accuracy (approx. 92%) on a new HSI data set collected over the heterogeneous forests of French Guiana (FG) when less than 1% of the pixels are labeled. We further show that GRNN is competitive with the state-of-the-art semi-supervised methods and exhibits a small deviation in accuracy for different numbers of training samples and over repeated trials with randomly sampled labeled pixels for training. △ Less

Submitted 5 May, 2023; v1 submitted 18 August, 2022; originally announced August 2022.

Ultra Long Turbulent Eddies, Magnetic Topology, and the Triggering of Internal Transport Barriers in Tokamaks

Authors: A. Volčokas, J. Ball, S. Brunner

Abstract: Local nonlinear gyrokinetic simulations of tokamak plasmas demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns when the magnetic shear is very small. By accurately modeling different field line topologies (e.g. low-order, almost rational, or irrational safety factor), we show that the parallel self-interaction of such "ultra long" eddies can drama… ▽ More Local nonlinear gyrokinetic simulations of tokamak plasmas demonstrate that turbulent eddies can extend along magnetic field lines for hundreds of poloidal turns when the magnetic shear is very small. By accurately modeling different field line topologies (e.g. low-order, almost rational, or irrational safety factor), we show that the parallel self-interaction of such "ultra long" eddies can dramatically reduce heat transport. This reveals novel strategies to improve confinement, constitutes experimentally testable predictions, and illuminates past observations of internal transport barriers. △ Less

Submitted 12 August, 2022; originally announced August 2022.

Bayesian Generational Population-Based Training

Authors: Xingchen Wan, Cong Lu, Jack Parker-Holder, Philip J. Ball, Vu Nguyen, Binxin Ru, Michael A. Osborne

Abstract: Reinforcement learning (RL) offers the potential for training generally capable agents that can interact autonomously in the real world. However, one key limitation is the brittleness of RL algorithms to core hyperparameters and network architecture choice. Furthermore, non-stationarities such as evolving training data and increased agent complexity mean that different hyperparameters and architec… ▽ More Reinforcement learning (RL) offers the potential for training generally capable agents that can interact autonomously in the real world. However, one key limitation is the brittleness of RL algorithms to core hyperparameters and network architecture choice. Furthermore, non-stationarities such as evolving training data and increased agent complexity mean that different hyperparameters and architectures may be optimal at different points of training. This motivates AutoRL, a class of methods seeking to automate these design choices. One prominent class of AutoRL methods is Population-Based Training (PBT), which have led to impressive performance in several large scale settings. In this paper, we introduce two new innovations in PBT-style methods. First, we employ trust-region based Bayesian Optimization, enabling full coverage of the high-dimensional mixed hyperparameter search space. Second, we show that using a generational approach, we can also learn both architectures and hyperparameters jointly on-the-fly in a single training run. Leveraging the new highly parallelizable Brax physics engine, we show that these innovations lead to large performance gains, significantly outperforming the tuned baseline while learning entire configurations on the fly. Code is available at https://github.com/xingchenwan/bgpbt. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: AutoML Conference 2022. 10 pages, 4 figure, 3 tables (28 pages, 10 figures, 7 tables including references and appendices)

Functional Models for Commuting Hilbert-space Contractions

Authors: Joseph A. Ball, Haripada Sau

Abstract: We develop a Sz.-Nagy--Foias-type functional model for a commutative contractive operator tuple $\underline{T} = (T_1, \dots, T_d)$ having $T = T_1 \cdots T_d$ equal to a completely nonunitary contraction. We identify additional invariants ${\mathbb G}_\sharp, {\mathbb W}_\sharp$ in addition to the Sz.-Nagy--Foias characteristic function $Θ_T$ for the product operator $T$ so that the combined trip… ▽ More We develop a Sz.-Nagy--Foias-type functional model for a commutative contractive operator tuple $\underline{T} = (T_1, \dots, T_d)$ having $T = T_1 \cdots T_d$ equal to a completely nonunitary contraction. We identify additional invariants ${\mathbb G}_\sharp, {\mathbb W}_\sharp$ in addition to the Sz.-Nagy--Foias characteristic function $Θ_T$ for the product operator $T$ so that the combined triple $({\mathbb G}_\sharp, {\mathbb W}_\sharp, Θ_T)$ becomes a complete unitary invariant for the original operator tuple $\underline{T}$. For the case $d \ge 3$ in general there is no commutative isometric lift of $\underline{T}$; however there is a (not necessarily commutative) isometric lift having some additional structure so that, when compressed to the minimal isometric-lift space for the product operator $T$, generates a special kind of lift of $\underline{T}$, herein called a {\em pseudo-commutative contractive lift} of $\underline{T}$, which in turn leads to the functional model for $\underline{T}$. This work has many parallels with recently developed model theories for symmetrized-bidisk contractions (commutative operator pairs $(S,P)$ having the symmetrized bidisk $Γ$ as a spectral set) and for tetrablock contractions (commutative operator triples $(A, B, P)$ having the tetrablock domain ${\mathbb E}$ as a spectral set). △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: Dedicated to the memory of Ron Douglas. It has appeared in Operator Theory: Advances and Applications (Ronald G. Douglas Memorial Volume)

Dilation theory and functional models for tetrablock contractions

Authors: Joseph A. Ball, Haripada Sau

Abstract: A classical result of Sz.-Nagy asserts that a Hilbert space contraction operator $T$ can be dilated to a unitary $\cU$. A more general multivariable setting for these ideas is the setup where (i) the unit disk is replaced by a domain $Ω$ contained in ${\mathbb C}^d$, (ii) the contraction operator $T$ is replaced by a commuting tuple $\bfT = (T_1, \dots, T_d)$ such that… ▽ More A classical result of Sz.-Nagy asserts that a Hilbert space contraction operator $T$ can be dilated to a unitary $\cU$. A more general multivariable setting for these ideas is the setup where (i) the unit disk is replaced by a domain $Ω$ contained in ${\mathbb C}^d$, (ii) the contraction operator $T$ is replaced by a commuting tuple $\bfT = (T_1, \dots, T_d)$ such that $\| r(T_1, \dots, T_d) \|_{\cL(\cH)} \le \sup_{\lam \in Ω} | r(\lam) |$ for all rational functions with no singularities in $\overlineΩ$ and the unitary operator $\cU$ is replaced by an $Ω$-unitary operator tuple, i.e., a commutative operator $d$-tuple $\bfU = (U_1, \dots, U_d)$ of commuting normal operators with joint spectrum contained in the distinguished boundary $bΩ$ of $Ω$. For a given domain $Ω\subset {\mathbb C}^d$, the {\em rational dilation question} asks: given an $Ω$-contraction $\bfT$ on $\cH$, is it always possible to find an $Ω$-unitary $\bfU$ on a larger Hilbert space $\cK \supset \cH$ so that, for any $d$-variable rational function without singularities in $\overlineΩ$, one can recover $r(T)$ as $r(T) = P_\cH r(\bfU)|_\cH$. We focus here on the case where $Ω$ is the {\em tetrablock}. (i) We identify a complete set of unitary invariants for a ${\mathbb E}$-contraction $(A,B,T)$ which can then be used to write down a functional model for $(A,B,T)$, thereby extending earlier results only done for a special case, (ii) we identify the class of {\em pseudo-commutative ${\mathbb E}$-isometries} (a priori slightly larger than the class of ${\mathbb E}$-isometries) to which any ${\mathbb E}$-contraction can be lifted, and (iii) we use our functional model to recover an earlier result on the existence and uniqueness of a ${\mathbb E}$-isometric lift $(V_1, V_2, V_3)$ of a special type for a ${\mathbb E}$-contraction $(A,B,T)$. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: Dedicated to the memory of Jörg Eschmeier

Models for q-commutative tuples of isometries

Authors: Joseph A. Ball, Haripada Sau

Abstract: A pair of Hilbert space linear operators $(V_1,V_2)$ is said to be $q$-commutative, for a unimodular complex number $q$, if $V_1V_2=qV_2V_1$. A concrete functional model for $q$-commutative pairs of isometries is obtained. The functional model is parametrized by a collection of Hilbert spaces and operators acting on them. As a consequence, the collection serves as a complete unitary invariance for… ▽ More A pair of Hilbert space linear operators $(V_1,V_2)$ is said to be $q$-commutative, for a unimodular complex number $q$, if $V_1V_2=qV_2V_1$. A concrete functional model for $q$-commutative pairs of isometries is obtained. The functional model is parametrized by a collection of Hilbert spaces and operators acting on them. As a consequence, the collection serves as a complete unitary invariance for $q$-commutative pairs of isometries. A $q$-commutative operator pair $(V_1,V_2)$ is said to be doubly $q$-commutative, if in addition, it satisfies $V_2V_1^*=qV_1^*V_2$. Doubly $q$-commutative pairs of isometries are also characterized. Special attention is given to doubly $q$-commutative pairs of shift operators. The notion of $q$-commutativity is then naturally extended to the case of general tuples of operators to obtain a similar model for tuples of $q$-commutative isometries. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Stabilizing Off-Policy Deep Reinforcement Learning from Pixels

Authors: Edoardo Cetin, Philip J. Ball, Steve Roberts, Oya Celiktutan

Abstract: Off-policy reinforcement learning (RL) from pixel observations is notoriously unstable. As a result, many successful algorithms must combine different domain-specific practices and auxiliary losses to learn meaningful behaviors in complex environments. In this work, we provide novel analysis demonstrating that these instabilities arise from performing temporal-difference learning with a convolutio… ▽ More Off-policy reinforcement learning (RL) from pixel observations is notoriously unstable. As a result, many successful algorithms must combine different domain-specific practices and auxiliary losses to learn meaningful behaviors in complex environments. In this work, we provide novel analysis demonstrating that these instabilities arise from performing temporal-difference learning with a convolutional encoder and low-magnitude rewards. We show that this new visual deadly triad causes unstable training and premature convergence to degenerate solutions, a phenomenon we name catastrophic self-overfitting. Based on our analysis, we propose A-LIX, a method providing adaptive regularization to the encoder's gradients that explicitly prevents the occurrence of catastrophic self-overfitting using a dual objective. By applying A-LIX, we significantly outperform the prior state-of-the-art on the DeepMind Control and Atari 100k benchmarks without any data augmentation or auxiliary losses. △ Less

Submitted 3 July, 2022; originally announced July 2022.

Comments: Short presentation at ICML 2022

Challenges and Opportunities in Offline Reinforcement Learning from Visual Observations

Authors: Cong Lu, Philip J. Ball, Tim G. J. Rudner, Jack Parker-Holder, Michael A. Osborne, Yee Whye Teh

Abstract: Offline reinforcement learning has shown great promise in leveraging large pre-collected datasets for policy learning, allowing agents to forgo often-expensive online data collection. However, offline reinforcement learning from visual observations with continuous action spaces remains under-explored, with a limited understanding of the key challenges in this complex domain. In this paper, we esta… ▽ More Offline reinforcement learning has shown great promise in leveraging large pre-collected datasets for policy learning, allowing agents to forgo often-expensive online data collection. However, offline reinforcement learning from visual observations with continuous action spaces remains under-explored, with a limited understanding of the key challenges in this complex domain. In this paper, we establish simple baselines for continuous control in the visual domain and introduce a suite of benchmarking tasks for offline reinforcement learning from visual observations designed to better represent the data distributions present in real-world offline RL problems and guided by a set of desiderata for offline RL from visual observations, including robustness to visual distractions and visually identifiable changes in dynamics. Using this suite of benchmarking tasks, we show that simple modifications to two popular vision-based online reinforcement learning algorithms, DreamerV2 and DrQ-v2, suffice to outperform existing offline RL methods and establish competitive baselines for continuous control in the visual domain. We rigorously evaluate these algorithms and perform an empirical evaluation of the differences between state-of-the-art model-based and model-free offline RL methods for continuous control from visual observations. All code and data used in this evaluation are open-sourced to facilitate progress in this domain. △ Less

Submitted 6 July, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: Published at TMLR, 2023

Interpolation in multivariable de Branges-Rovnyak spaces

Authors: Joseph A. Ball, Vladimir Bolotnikov, Sanne ter Horst

Abstract: We study a general metric constrained interpolation problem in a de Branges-Rovnyak space $\mathcal{H}(K_S)$ associated with a contractive multiplier $S$ between two Fock spaces along with its commutative counterpart, a de Branges-Rovnyak space associated with a Schur multiplier on the Drury-Arveson space of the unit ball of $\mathbb{C}^n$. We study a general metric constrained interpolation problem in a de Branges-Rovnyak space $\mathcal{H}(K_S)$ associated with a contractive multiplier $S$ between two Fock spaces along with its commutative counterpart, a de Branges-Rovnyak space associated with a Schur multiplier on the Drury-Arveson space of the unit ball of $\mathbb{C}^n$. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 47 pages

MSC Class: 46E22; 47A57; 30E05

Three-Dimensional Inhomogeneity of Electron-Temperature-Gradient Turbulence in the Edge of Tokamak Plasmas

Authors: J. F. Parisi, F. I. Parra, C. M. Roach, M. R. Hardman, A. A. Schekochihin, I. G. Abel, N. Aiba, J. Ball, M. Barnes, B. Chapman-Oplopoiou, D. Dickinson, W. Dorland, C. Giroud, D. R. Hatch, J. C. Hillesheim, J. Ruiz Ruiz, S. Saarelma, D. St-Onge

Abstract: Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a… ▽ More Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport -- in our case by $\sim$40 \% -- through multiscale interactions. △ Less

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

Comments: 17 pages, 14 figures

Free noncommutative hereditary kernels: Jordan decomposition, Arveson extension, kernel domination

Authors: Joseph A. Ball, Gregory Marx, Victor Vinnikov

Abstract: We discuss a (i) quantized version of the Jordan decomposition theorem for a complex Borel measure on a compact Hausdorff space, namely, the more general problem of decomposing a general noncommutative kernel (a quantization of the standard notion of kernel function) as a linear combination of completely positive noncommutative kernels (a quantization of the standard notion of positive definite ke… ▽ More We discuss a (i) quantized version of the Jordan decomposition theorem for a complex Borel measure on a compact Hausdorff space, namely, the more general problem of decomposing a general noncommutative kernel (a quantization of the standard notion of kernel function) as a linear combination of completely positive noncommutative kernels (a quantization of the standard notion of positive definite kernel). Other special cases of (i) include: the problem of decomposing a general operator-valued kernel function as a linear combination of positive kernels (not always possible), of decomposing a general bounded linear Hilbert-space operator as a linear combination of positive linear operators (always possible), of decomposing a completely bounded linear map from a $C^*$-algebra ${\mathcal A}$ to an injective $C^*$-algebra ${\mathcal L}({\mathcal Y})$ as a linear combination of completely positive maps from ${\mathcal A}$ to ${\mathcal L}({\mathcal Y})$ (always possible). We also discuss (ii) a noncommutative kernel generalization of the Arveson extension theorem (any completely positive map $φ$ from a operator system ${\mathbb S}$ to an injective $C^*$-algebra ${\mathcal L}({\mathcal Y})$ can be extended to a completely positive map $φ_e$ from a $C^*$-algebra containing ${\mathbb S}$ to ${\mathcal L}({\mathcal Y})$), and (iii) a noncommutative kernel version of a Positivstellensatz (i.e., finding a certificate to explain why one kernel is positive at points where another given kernel is positive). △ Less

Submitted 2 February, 2022; originally announced February 2022.

MSC Class: 47B32; 47A60

The mass distribution in the Galactic Centre from interferometric astrometry of multiple stellar orbits

Authors: GRAVITY Collaboration, R. Abuter, N. Aimar, A. Amorim, J. Ball, M. Bauböck, J. P. Berger, H. Bonnet, G. Bourdarot, W. Brandner, V. Cardoso, Y. Clénet, Y. Dallilar, R. Davies, P. T. de Zeeuw, J. Dexter, A. Drescher, F. Eisenhauer, N. M. Förster Schreiber, A. Foschi, P. Garcia, F. Gao, E. Gendron, R. Genzel, S. Gillessen , et al. (40 additional authors not shown)

Abstract: The stars orbiting the compact radio source Sgr A* in the Galactic Centre are precision probes of the gravitational field around the closest massive black hole. In addition to adaptive optics assisted astrometry (with NACO / VLT) and spectroscopy (with SINFONI / VLT, NIRC2 / Keck and GNIRS / Gemini) over three decades, since 2016/2017 we have obtained 30-100 mu-as astrometry with the four-telescop… ▽ More The stars orbiting the compact radio source Sgr A* in the Galactic Centre are precision probes of the gravitational field around the closest massive black hole. In addition to adaptive optics assisted astrometry (with NACO / VLT) and spectroscopy (with SINFONI / VLT, NIRC2 / Keck and GNIRS / Gemini) over three decades, since 2016/2017 we have obtained 30-100 mu-as astrometry with the four-telescope interferometric beam combiner GRAVITY / VLTI reaching a sensitivity of mK = 20 when combining data from one night. We present the simultaneous detection of several stars within the diffraction limit of a single telescope, illustrating the power of interferometry. The new data for the stars S2, S29, S38 and S55 yield significant accelerations between March and July 2021, as these stars pass the pericenters of their orbits between 2018 and 2023. This allows for a high-precision determination of the gravitational potential around Sgr A*. Our data are in excellent agreement with general relativity orbits around a single central point mass, M = 4.30 x 10^6 M_sun with a precision of about +-0.25%. We improve the significance of our detection of the Schwarzschild precession in the S2 orbit to 7 sigma. Assuming plausible density profiles, an extended mass component inside S2's apocentre (= 0.23" or 2.4 x 10^4 R_S) must be 3000 M_sun (1 sigma), or 0.1% of M. Adding the enclosed mass determinations from 13 stars orbiting Sgr A* at larger radii, the innermost radius at which the excess mass beyond Sgr A* tentatively is seen is r = 2.5" >= 10x the apocentre of S2. This is in full harmony with the stellar mass distribution (including stellar-mass black holes) obtained from the spatially resolved luminosity function. △ Less

Submitted 14 December, 2021; originally announced December 2021.

Comments: published in A&A

Journal ref: A&A, 2021,

Implementing and evaluating far-field 3D X-ray diffraction at the I12 JEEP beamline, Diamond Light Source

Authors: James A. D. Ball, Anna Kareer, Oxana V. Magdysyuk, Stefan Michalik, Anastasia Vrettou, Neal Parkes, Thomas Connolley, David M. Collins

Abstract: Three-dimensional X-ray diffraction (3DXRD) is shown to be feasible at the I12 Joint Engineering, Environmental and Processing (JEEP) beamline of Diamond Light Source. As a demonstration, a microstructually simple low-carbon ferritic steel was studied in a highly textured and annealed state. A processing pipeline suited to this beamline was created, using software already established in the 3DXRD… ▽ More Three-dimensional X-ray diffraction (3DXRD) is shown to be feasible at the I12 Joint Engineering, Environmental and Processing (JEEP) beamline of Diamond Light Source. As a demonstration, a microstructually simple low-carbon ferritic steel was studied in a highly textured and annealed state. A processing pipeline suited to this beamline was created, using software already established in the 3DXRD user community, enabling grain centre-of-mass positions, orientations and strain tensor elements to be determined. Orientations, with texture measurements independently validated from electron backscatter diffraction (EBSD) data, possessed a $\sim 0.1^\circ$ uncertainty, comparable with other 3DXRD instruments. The spatial resolution was limited by the far-field detector pixel size; the average of the grain centre of mass position errors was determined as $\pm \sim 80 μ$m. An average per-grain error of $\sim 1 \times 10^{-3}$ for the elastic strains was also measured; this could be reduced in future experiments by improving sample preparation, data collection and analysis techniques. Application of 3DXRD onto I12 shows great potential, where its implementation is highly desirable due the flexible, open architecture of the beamline. User-owned or designed sample environments can be used, thus 3DXRD could be applied to previously unexplored scientific areas. △ Less

Submitted 20 February, 2023; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 13 pages, 10 figures. Published in Journal of Synchrotron Radiation

Journal ref: Journal of Synchrotron Radiation, July 2022, Volume 29, Pages 1043-1053

Revisiting Design Choices in Offline Model-Based Reinforcement Learning

Authors: Cong Lu, Philip J. Ball, Jack Parker-Holder, Michael A. Osborne, Stephen J. Roberts

Abstract: Offline reinforcement learning enables agents to leverage large pre-collected datasets of environment transitions to learn control policies, circumventing the need for potentially expensive or unsafe online data collection. Significant progress has been made recently in offline model-based reinforcement learning, approaches which leverage a learned dynamics model. This typically involves construct… ▽ More Offline reinforcement learning enables agents to leverage large pre-collected datasets of environment transitions to learn control policies, circumventing the need for potentially expensive or unsafe online data collection. Significant progress has been made recently in offline model-based reinforcement learning, approaches which leverage a learned dynamics model. This typically involves constructing a probabilistic model, and using the model uncertainty to penalize rewards where there is insufficient data, solving for a pessimistic MDP that lower bounds the true MDP. Existing methods, however, exhibit a breakdown between theory and practice, whereby pessimistic return ought to be bounded by the total variation distance of the model from the true dynamics, but is instead implemented through a penalty based on estimated model uncertainty. This has spawned a variety of uncertainty heuristics, with little to no comparison between differing approaches. In this paper, we compare these heuristics, and design novel protocols to investigate their interaction with other hyperparameters, such as the number of models, or imaginary rollout horizon. Using these insights, we show that selecting these key hyperparameters using Bayesian Optimization produces superior configurations that are vastly different to those currently used in existing hand-tuned state-of-the-art methods, and result in drastically stronger performance. △ Less

Submitted 16 March, 2022; v1 submitted 8 October, 2021; originally announced October 2021.

Comments: Spotlight @ ICLR 2022; Spotlight @ RL4RealLife Workshop ICML2021

The Infinite-Dimensional Standard and Strict Bounded Real Lemmas in Continuous Time: The storage function approach

Authors: Joseph A. Ball, Sanne ter Horst, Mikael Kurula

Abstract: The bounded real lemma (BRL) is a classical result in systems theory, which provides a linear matrix inequality criterium for dissipativity, via the Kalman-Yakubovich-Popov (KYP) inequality. The BRL has many applications, among others in H-infinity control. Extensions to infinite dimensional systems, although already present in the work of Yakubovich, have only been studied systematically in the l… ▽ More The bounded real lemma (BRL) is a classical result in systems theory, which provides a linear matrix inequality criterium for dissipativity, via the Kalman-Yakubovich-Popov (KYP) inequality. The BRL has many applications, among others in H-infinity control. Extensions to infinite dimensional systems, although already present in the work of Yakubovich, have only been studied systematically in the last few decades. In this context various notions of stability, observability and controllability exist, and depending on the hypothesis one may have to allow the KYP-inequality to have unbounded solutions which forces one to consider the KYP-inequality in a spatial form. In the present paper we consider the BRL for continuous time, infinite dimensional, linear well-posed systems. Via an adaptation of Willems' storage function approach we present a unified way to address both the standard and strict forms of the BRL. We avoid making use of the Cayley transform and work only in continuous time. While for the standard bounded real lemma, we obtain analogous results as there exist for the discrete time case, when treating the strict case additional conditions are required, at least at this stage. This might be caused by the fact that the Cayley transform does not preserve exponential stability, an important property in the strict case, when transferring a continuous-time system to a discrete-time system. △ Less

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

Comments: 65 pages

Two-Orders-of-Magnitude Improvement in the Total Spin Angular Momentum of 131Xe Nuclei Using Spin Exchange Optical Pumping

Authors: Michael J. Molway, Liana Bales-Shaffer, Kaili Ranta, Dustin Basler, Megan Murphy, Bryce E. Kidd, Abdulbasit Tobi Gafar, Justin Porter, Kierstyn Albin, Boyd M. Goodson, Eduard Y. Chekmenev, Matthew S. Rosen, W. Michael Snow, James Ball, Eleanor Sparling, Mia Prince, Daniel Cocking, Michael J. Barlow

Abstract: We report on hyperpolarization of quadrupolar (I=3/2) 131Xe via spin-exchange optical pumping. Observations of the 131Xe polarization dynamics show that the effective alkali-metal/131Xe spin-exchange cross-sections are large enough to compete with 131Xe spin relaxation. 131Xe polarization up to 7.6 p/m 1.5 percent was achieved in ca. 8.5EE20 spins--a ca. 100-fold improvement in the total spin angu… ▽ More We report on hyperpolarization of quadrupolar (I=3/2) 131Xe via spin-exchange optical pumping. Observations of the 131Xe polarization dynamics show that the effective alkali-metal/131Xe spin-exchange cross-sections are large enough to compete with 131Xe spin relaxation. 131Xe polarization up to 7.6 p/m 1.5 percent was achieved in ca. 8.5EE20 spins--a ca. 100-fold improvement in the total spin angular momentum--enabling applications including measurement of spin-dependent neutron-131Xe s-wave scattering and sensitive searches for time-reversal violation in neutron-131Xe interactions beyond the Standard Model. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 6 page main document plus 23 page Supplemental Information (SI) document