Showing 1–50 of 94 results for author: Johnston, H

Determination of the stably free cancellation property for orders

Authors: Werner Bley, Tommy Hofmann, Henri Johnston

Abstract: Let $K$ be a number field, let $A$ be a finite-dimensional semisimple $K$-algebra, and let $Λ$ be an $\mathcal{O}_{K}$-order in $A$. We give practical algorithms that determine whether or not $Λ$ has the stably free cancellation property (SFC). As an application, we determine all finite groups $G$ of order at most $383$ such that the integral group ring $\mathbb{Z}[G]$ has SFC. Let $K$ be a number field, let $A$ be a finite-dimensional semisimple $K$-algebra, and let $Λ$ be an $\mathcal{O}_{K}$-order in $A$. We give practical algorithms that determine whether or not $Λ$ has the stably free cancellation property (SFC). As an application, we determine all finite groups $G$ of order at most $383$ such that the integral group ring $\mathbb{Z}[G]$ has SFC. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 40 pages, comments welcome!

MSC Class: 16H10; 16H20; 16Z05; 20C05; 20C10; 11R33; 11Y40

Formation of giant planets around intermediate-mass stars

Authors: Heather Johnston, Olja Panic, Beibei Liu

Abstract: To understand giant planet formation, we need to focus on host stars close to $1.7\ \rm M_{\odot}$, where the occurrence rate of these planets is the highest. In this initial study, we carry out pebble-driven core accretion planet formation modelling to investigate the trends and optimal conditions for the formation of giant planets around host stars in the range of $1{-}2.4\ \rm M_{\odot}$. We fi… ▽ More To understand giant planet formation, we need to focus on host stars close to $1.7\ \rm M_{\odot}$, where the occurrence rate of these planets is the highest. In this initial study, we carry out pebble-driven core accretion planet formation modelling to investigate the trends and optimal conditions for the formation of giant planets around host stars in the range of $1{-}2.4\ \rm M_{\odot}$. We find that giant planets are more likely to form in systems with a larger initial disk radius; higher disk gas accretion rate; pebbles of $\sim$ millimeter in size; and birth location of the embryo at a moderate radial distance of $\sim 10$ AU. We also conduct a population synthesis study of our model and find that the frequency of giant planets and super-Earths decreases with increasing stellar mass. This contrasts the observational peak at $1.7\ \rm M_{\odot}$, stressing the need for strong assumptions on stellar mass dependencies in this range. Investigating the combined effect of stellar mass dependent disk masses, sizes, and lifetimes in the context of planet population synthesis studies is a promising avenue to alleviate this discrepancy. The hot-Jupiter occurrence rate in our models is $\sim 0.7{-}0.8\%$ around $1\ \rm M_{\odot}$ - similar to RV observations around Sun-like stars, but drastically decreases for higher mass stars. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 21 pages, 11 figures

On the existence of free sublattices of bounded index and arithmetic applications

Authors: Henri Johnston, Alex Torzewski

Abstract: Let $\mathcal{O}$ be a Dedekind domain whose field of fractions $K$ is a global field. Let $A$ be a finite-dimensional separable $K$-algebra and let $Λ$ be an $\mathcal{O}$-order in $A$. Let $n$ be a positive integer and suppose that $X$ is a $Λ$-lattice such that $K \otimes_{\mathcal{O}} X$ is free of rank $n$ over $A$. Then $X$ contains a (non-unique) free $Λ$-sublattice of rank $n$. The main re… ▽ More Let $\mathcal{O}$ be a Dedekind domain whose field of fractions $K$ is a global field. Let $A$ be a finite-dimensional separable $K$-algebra and let $Λ$ be an $\mathcal{O}$-order in $A$. Let $n$ be a positive integer and suppose that $X$ is a $Λ$-lattice such that $K \otimes_{\mathcal{O}} X$ is free of rank $n$ over $A$. Then $X$ contains a (non-unique) free $Λ$-sublattice of rank $n$. The main result of the present article is to show there exists such a sublattice $Y$ such that the generalised module index $[X : Y]_{\mathcal{O}}$ has explicit upper bounds with respect to division that are independent of $X$ and can be chosen to satisfy certain conditions. We give examples of applications to the approximation of normal integral bases and strong Minkowski units, and to the Galois module structure of rational points over abelian varieties. △ Less

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

Comments: 22 pages; changes and corrections following referee reports

MSC Class: 16H20; 11R33; 11R27; 11G05; 11G10

Journal ref: Journal of Algebra 657 (2024) 81-108

Design Considerations for 3RRR Parallel Robots with Lightweight, Approximate Static-Balancing

Authors: Giuseppe Del Giudice, Garrison L. H. Johnston, Nabil Simaan

Abstract: Balancing parallel robots throughout their workspace while avoiding the use of balancing masses and respecting design practicality constraints is difficult. Medical robots demand such compact and lightweight designs. This paper considers the difficult task of achieving optimal approximate balancing of a parallel robot throughout a desired task-based dexterous workspace using balancing springs only… ▽ More Balancing parallel robots throughout their workspace while avoiding the use of balancing masses and respecting design practicality constraints is difficult. Medical robots demand such compact and lightweight designs. This paper considers the difficult task of achieving optimal approximate balancing of a parallel robot throughout a desired task-based dexterous workspace using balancing springs only. While it is possible to achieve perfect balancing in a path, only approximate balancing may be achieved without the addition of balancing masses. Design considerations for optimal robot base placement and the effects of placement of torsional balancing springs are presented. Using a modal representation for the balancing torque requirements, we use recent results on the design of wire-wrapped cam mechanisms to achieve balancing throughout a task-based workspace. A simulation study shows that robot base placement can have a detrimental effect on the attainability of a practical design solution for static balancing. We also show that optimal balancing using torsional springs is best achieved when all springs are at the actuated joints and that the wire-wrapped cam design can significantly improve the performance of static balancing. The methodology presented in this paper provides practical design solutions that yield simple, lightweight and compact designs suitable for medical applications where such traits are paramount. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Journal ref: 2023 ASME International Design Engineering Technical Conferences & Computers and Information in Engineering Conference (IDTEC-CIE)

Task and Configuration Space Compliance of Continuum Robots via Lie Group and Modal Shape Formulations

Authors: Andrew L. Orekhov, Garrison L. H. Johnston, Nabil Simaan

Abstract: Continuum robots suffer large deflections due to internal and external forces. Accurate modeling of their passive compliance is necessary for accurate environmental interaction, especially in scenarios where direct force sensing is not practical. This paper focuses on deriving analytic formulations for the compliance of continuum robots that can be modeled as Kirchhoff rods. Compared to prior work… ▽ More Continuum robots suffer large deflections due to internal and external forces. Accurate modeling of their passive compliance is necessary for accurate environmental interaction, especially in scenarios where direct force sensing is not practical. This paper focuses on deriving analytic formulations for the compliance of continuum robots that can be modeled as Kirchhoff rods. Compared to prior works, the approach presented herein is not subject to the constant-curvature assumptions to derive the configuration space compliance, and we do not rely on computationally-expensive finite difference approximations to obtain the task space compliance. Using modal approximations over curvature space and Lie group integration, we obtain closed-form expressions for the task and configuration space compliance matrices of continuum robots, thereby bridging the gap between constant-curvature analytic formulations of configuration space compliance and variable curvature task space compliance. We first present an analytic expression for the compliance of a single Kirchhoff rod. We then extend this formulation for computing both the task space and configuration space compliance of a tendon-actuated continuum robot. We then use our formulation to study the tradeoffs between computation cost and modeling accuracy as well as the loss in accuracy from neglecting the Jacobian derivative term in the compliance model. Finally, we experimentally validate the model on a tendon-actuated continuum segment, demonstrating the model's ability to predict passive deflections with error below 11.5\% percent of total arc length. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Journal ref: 2023 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS)

Design Considerations and Robustness to Parameter Uncertainty in Wire-Wrapped Cam Mechanisms

Authors: Garrison L. H. Johnston, Andrew L. Orekhov, Nabil Simaan

Abstract: Collaborative robots must simultaneously be safe enough to operate in close proximity to human operators and powerful enough to assist users in industrial tasks such as lifting heavy equipment. The requirement for safety necessitates that collaborative robots are designed with low-powered actuators. However, some industrial tasks may require the robot to have high payload capacity and/or long reac… ▽ More Collaborative robots must simultaneously be safe enough to operate in close proximity to human operators and powerful enough to assist users in industrial tasks such as lifting heavy equipment. The requirement for safety necessitates that collaborative robots are designed with low-powered actuators. However, some industrial tasks may require the robot to have high payload capacity and/or long reach. For collaborative robot designs to be successful, they must find ways of addressing these conflicting design requirements. One promising strategy for navigating this tradeoff is through the use of static balancing mechanisms to offset the robot's self weight, thus enabling the selection of lower-powered actuators. In this paper, we introduce a novel, 2 degree of freedom static balancing mechanism based on spring-loaded, wire-wrapped cams. We also present an optimization-based cam design method that guarantees the cams stay convex, ensures the springs stay below their extensions limits, and minimizes sensitivity to unmodeled deviations from the nominal spring constant. Additionally, we present a model of the effect of friction between the wire and the cam. Lastly, we show experimentally that the torque generated by the cam mechanism matches the torque predicted in our modeling approach. Our results also suggest that the effects of wire-cam friction are significant for non-circular cams. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Journal ref: ASME Journal of Mechanisms and Robotics. February 2024. 16(2): 021001

The Intrinsic Alignment of Galaxy Clusters and Impact of Projection Effects

Authors: Jingjing Shi, Tomomi Sunayama, Toshiki Kurita, Masahiro Takada, Sunao Sugiyama, Rachel Mandelbaum, Hironao Miyatake, Surhud More, Takahiro Nishimichi, Harry Johnston

Abstract: Galaxy clusters, being the most massive objects in the Universe, exhibit the strongest alignment with the large-scale structure. However, mis-identification of members due to projection effects from the large scale structure can occur. We studied the impact of projection effects on the measurement of the intrinsic alignment of galaxy clusters, using galaxy cluster mock catalogs. Our findings showe… ▽ More Galaxy clusters, being the most massive objects in the Universe, exhibit the strongest alignment with the large-scale structure. However, mis-identification of members due to projection effects from the large scale structure can occur. We studied the impact of projection effects on the measurement of the intrinsic alignment of galaxy clusters, using galaxy cluster mock catalogs. Our findings showed that projection effects result in a decrease of the large scale intrinsic alignment signal of the cluster and produce a bump at $r_p\sim 1h^{-1}/Mpc$, most likely due to interlopers and missed member galaxies. This decrease in signal explains the observed similar alignment strength between bright central galaxies and clusters in the SDSS redMaPPer cluster catalog. The projection effect and cluster intrinsic alignment signal are coupled, with clusters having lower fractions of missing members or having higher fraction of interlopers exhibiting higher alignment signals in their projected shapes. We aim to use these findings to determine the impact of projection effects on galaxy cluster cosmology in future studies. △ Less

Submitted 10 January, 2024; v1 submitted 16 June, 2023; originally announced June 2023.

Comments: 15 pages, 14 figures, Accepted by MNRAS

Applications of representation theory and of explicit units to Leopoldt's conjecture

Authors: Fabio Ferri, Henri Johnston

Abstract: Let $L/K$ be a Galois extension of number fields and let $G=\mathrm{Gal}(L/K)$. We show that under certain hypotheses on $G$, for a fixed prime number $p$, Leopoldt's conjecture at $p$ for certain proper intermediate fields of $L/K$ implies Leopoldt's conjecture at $p$ for $L$. We also obtain relations between the Leopoldt defects of intermediate extensions of $L/K$. By applying a result of Buchma… ▽ More Let $L/K$ be a Galois extension of number fields and let $G=\mathrm{Gal}(L/K)$. We show that under certain hypotheses on $G$, for a fixed prime number $p$, Leopoldt's conjecture at $p$ for certain proper intermediate fields of $L/K$ implies Leopoldt's conjecture at $p$ for $L$. We also obtain relations between the Leopoldt defects of intermediate extensions of $L/K$. By applying a result of Buchmann and Sands together with an explicit description of units and a special case of the above results, we show that given any finite set of prime numbers $\mathcal{P}$, there exists an infinite family $\mathcal{F}$ of totally real $S_{3}$-extensions of $\mathbb{Q}$ such that Leopoldt's conjecture for $F$ at $p$ holds for every $F \in \mathcal{F}$ and $p \in \mathcal{P}$. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 23 pages; preliminary version; comments welcome

MSC Class: 11R23 (Primary) 11R27; 19A22; 20C15 (Secondary)

Long Horizon Planning through Contact using Discrete Search and Continuous Optimization

Authors: Ramkumar Natarajan, Garrison L. H. Johnston, Nabil Simaan, Maxim Likhachev, Howie Choset

Abstract: Robots often have to perform manipulation tasks in close proximity to people. As such, it is desirable to use a robot arm that has limited joint torques to not injure the nearby person and interacts with the environment to explore new possibilities for completing a task. By bracing against the environment, robots can expand their reachable workspace, which would otherwise be inaccessible due to ex… ▽ More Robots often have to perform manipulation tasks in close proximity to people. As such, it is desirable to use a robot arm that has limited joint torques to not injure the nearby person and interacts with the environment to explore new possibilities for completing a task. By bracing against the environment, robots can expand their reachable workspace, which would otherwise be inaccessible due to exceeding actuator torque limits, and accomplish tasks beyond their design specifications. However, motion planning for complex contact-rich tasks requires reasoning through the permutations of different possible contact modes and bracing locations, which grow exponentially with the number of contact points and links in the robot. To address this combinatorial problem, we developed INSAT, which interleaves graph search to explore the manipulator joint configuration and the contact mode space with incremental trajectory optimizations seeded by neighborhood solutions to find a dynamically feasible trajectory through contact. In this paper, we present recent additions to the INSAT algorithm that improve its runtime performance. In particular, we propose Lazy INSAT with reduced optimization rejection that systematically procrastinates its calls to trajectory optimization while reusing feasible solutions that violate boundary constraints. The algorithm is evaluated on a heavy payload transportation task in simulation and on physical hardware. In simulation, we show that Lazy INSAT can discover solutions for tasks that cannot be accomplished within its design limits and without interacting with the environment. In comparison to executing the same trajectory without environment support, we demonstrate that the utilization of bracing contacts reduces the overall torque required to execute the trajectory. △ Less

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

Comments: Updated journal version under review

KiDS-1000: Combined halo-model cosmology constraints from galaxy abundance, galaxy clustering and galaxy-galaxy lensing

Authors: Andrej Dvornik, Catherine Heymans, Marika Asgari, Constance Mahony, Benjamin Joachimi, Maciej Bilicki, Elisa Chisari, Hendrik Hildebrandt, Henk Hoekstra, Harry Johnston, Konrad Kuijken, Alexander Mead, Hironao Miyatake, Takahiro Nishimichi, Robert Reischke, Sandra Unruh, Angus H. Wright

Abstract: We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matt… ▽ More We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matter field. Our magnitude-limited bright galaxy sample combines 9-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter $S_8=σ_8\sqrt{Ω_{\mathrm{m}}/0.3}=0.773^{+0.028}_{-0.030}$, and the matter density parameter $Ω_{\mathrm{m}}=0.290^{+0.021}_{-0.017}$. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint $3\times2{\mathrm{pt}}$ clustering-lensing analyses that additionally include a cosmic shear observable. This analysis therefore brings attention to the significant constraining power in the often-excluded non-linear scales for galaxy clustering and galaxy-galaxy lensing observables. By adopting a theoretical model that accounts for non-linear halo bias, halo exclusion, scale-dependent galaxy bias and the impact of baryon feedback, this work demonstrates the potential and a way forward to include non-linear scales in cosmological analyses. Varying the width of the satellite galaxy distribution with an additional parameter yields a strong preference for sub-Poissonian variance, improving the goodness of fit by 0.18 in reduced $χ^{2}$ value compared to a fixed Poisson distribution. △ Less

Submitted 13 June, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 25 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 675, A189 (2023)

Intrinsic correlations of galaxy sizes in a hydrodynamical cosmological simulation

Authors: Harry Johnston, Dana Sophia Westbeek, Sjoerd Weide, Nora Elisa Chisari, Yohan Dubois, Julien Devriendt, Christophe Pichon

Abstract: Residuals between measured galactic radii and those predicted by the Fundamental Plane (FP) are possible tracers of weak lensing magnification. However, observations have shown these to be systematically correlated with the large-scale structure. We use the Horizon-AGN hydrodynamical cosmological simulation to analyse these intrinsic size correlations (ISCs) for both elliptical (early-type) and sp… ▽ More Residuals between measured galactic radii and those predicted by the Fundamental Plane (FP) are possible tracers of weak lensing magnification. However, observations have shown these to be systematically correlated with the large-scale structure. We use the Horizon-AGN hydrodynamical cosmological simulation to analyse these intrinsic size correlations (ISCs) for both elliptical (early-type) and spiral (late-type) galaxies at $z=0.06$. We fit separate FPs to each sample, finding similarly distributed radius residuals, $λ$, in each case. We find persistent $λλ$ correlations over three-dimensional separations $0.5-17\,h^{-1}{\rm{Mpc}}$ in the case of spiral galaxies, at $>3σ$ significance. When relaxing a mass-selection, applied for better agreement with galaxy clustering constraints, the spiral $λλ$ detection strengthens to $9σ$; we detect a $5σ$ density-$λ$ correlation; and we observe intrinsically-large spirals to cluster more strongly than small spirals over scales $\lesssim10\,h^{-1}{\rm{Mpc}}$, at $>5σ$ significance. Conversely, and in agreement with the literature, we observe lower-mass, intrinsically-small ellipticals to cluster more strongly than their large counterparts over scales $0.5-17\,h^{-1}{\rm{Mpc}}$, at $>5σ$ significance. We model $λλ$ correlations using a phenomenological non-linear size model, and predict the level of contamination for cosmic convergence analyses. We find the systematic contribution to be of similar order to, or dominant over the cosmological signal. We make a mock measurement of an intrinsic, systematic contribution to the projected surface mass density $Σ(r)$ and find statistically significant, low-amplitude, positive (negative) contributions from lower-mass spirals (ellipticals), which may be of concern for large-scale ($\gtrsim\,7\,h^{-1}$ Mpc) measurements. △ Less

Submitted 17 January, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 26 pages, 11 figures, matches MNRAS accepted version

Euclid: Calibrating photometric redshifts with spectroscopic cross-correlations

Authors: K. Naidoo, H. Johnston, B. Joachimi, J. L. van den Busch, H. Hildebrandt, O. Ilbert, O. Lahav, N. Aghanim, B. Altieri, A. Amara, M. Baldi, R. Bender, C. Bodendorf, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou , et al. (104 additional authors not shown)

Abstract: Cosmological constraints from key probes of the Euclid imaging survey rely critically on the accurate determination of the true redshift distributions, $n(z)$, of tomographic redshift bins. We determine whether the mean redshift, $<z>$, of ten Euclid tomographic redshift bins can be calibrated to the Euclid target uncertainties of $σ(<z>)<0.002\,(1+z)$ via cross-correlation, with spectroscopic sam… ▽ More Cosmological constraints from key probes of the Euclid imaging survey rely critically on the accurate determination of the true redshift distributions, $n(z)$, of tomographic redshift bins. We determine whether the mean redshift, $<z>$, of ten Euclid tomographic redshift bins can be calibrated to the Euclid target uncertainties of $σ(<z>)<0.002\,(1+z)$ via cross-correlation, with spectroscopic samples akin to those from the Baryon Oscillation Spectroscopic Survey (BOSS), Dark Energy Spectroscopic Instrument (DESI), and Euclid's NISP spectroscopic survey. We construct mock Euclid and spectroscopic galaxy samples from the Flagship simulation and measure small-scale clustering redshifts up to redshift $z<1.8$ with an algorithm that performs well on current galaxy survey data. The clustering measurements are then fitted to two $n(z)$ models: one is the true $n(z)$ with a free mean; the other a Gaussian Process modified to be restricted to non-negative values. We show that $<z>$ is measured in each tomographic redshift bin to an accuracy of order 0.01 or better. By measuring the clustering redshifts on subsets of the full Flagship area, we construct scaling relations that allow us to extrapolate the method performance to larger sky areas than are currently available in the mock. For the full expected Euclid, BOSS, and DESI overlap region of approximately 6000 deg$^{2}$, the uncertainties attainable by clustering redshifts exceeds the Euclid requirement by at least a factor of three for both $n(z)$ models considered, although systematic biases limit the accuracy. Clustering redshifts are an extremely effective method for redshift calibration for Euclid if the sources of systematic biases can be determined and removed, or calibrated-out with sufficiently realistic simulations. We outline possible future work, in particular an extension to higher redshifts with quasar reference samples. △ Less

Submitted 24 February, 2023; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: 14 pages, 8 figures, accepted for publication in A&A

Journal ref: A&A 670, A149 (2023)

Computation of lattice isomorphisms and the integral matrix similarity problem

Authors: Werner Bley, Tommy Hofmann, Henri Johnston

Abstract: Let $K$ be a number field, let $A$ be a finite-dimensional $K$-algebra, let $\mathrm{J}(A)$ denote the Jacobson radical of $A$, and let $Λ$ be an $\mathcal{O}_{K}$-order in $A$. Suppose that each simple component of the semisimple $K$-algebra $A/{\mathrm{J}(A)}$ is isomorphic to a matrix ring over a field. Under this hypothesis on $A$, we give an algorithm that given two $Λ$-lattices $X$ and $Y$,… ▽ More Let $K$ be a number field, let $A$ be a finite-dimensional $K$-algebra, let $\mathrm{J}(A)$ denote the Jacobson radical of $A$, and let $Λ$ be an $\mathcal{O}_{K}$-order in $A$. Suppose that each simple component of the semisimple $K$-algebra $A/{\mathrm{J}(A)}$ is isomorphic to a matrix ring over a field. Under this hypothesis on $A$, we give an algorithm that given two $Λ$-lattices $X$ and $Y$, determines whether $X$ and $Y$ are isomorphic, and if so, computes an explicit isomorphism $X \rightarrow Y$. This algorithm reduces the problem to standard problems in computational algebra and algorithmic algebraic number theory in polynomial time. As an application, we give an algorithm for the following long-standing problem: given a number field $K$, a positive integer $n$ and two matrices $A,B \in \mathrm{Mat}_{n}(\mathcal{O}_{K})$, determine whether $A$ and $B$ are similar over $\mathcal{O}_{K}$, and if so, return a matrix $C \in \mathrm{GL}_{n}(\mathcal{O}_{K})$ such that $B= CAC^{-1}$. We give explicit examples that show that the implementation of the latter algorithm for $\mathcal{O}_{K}=\mathbb{Z}$ vastly outperforms implementations of all previous algorithms, as predicted by our complexity analysis. △ Less

Submitted 31 August, 2022; v1 submitted 7 February, 2022; originally announced February 2022.

Comments: 37 pages; v3 revised and accepted version to appear in Forum of Mathematics, Sigma; v2 has many minor corrections and changes

MSC Class: 11R33; 11Y40; 16Z05; 20G30

KiDS-1000: Constraints on the intrinsic alignment of luminous red galaxies

Authors: Maria Cristina Fortuna, Henk Hoekstra, Harry Johnston, Mohammadjavad Vakili, Arun Kannawadi, Christos Georgiou, Benjamin Joachimi, Angus H. Wright, Marika Asgari, Maciej Bilicki, Catherine Heymans, Hendrik Hildebrandt, Konrad Kuijken, Maximilian Von Wietersheim-Kramsta

Abstract: We constrain the luminosity and redshift dependence of the intrinsic alignment (IA) of a nearly volume-limited sample of luminous red galaxies selected from the fourth public data release of the Kilo-Degree Survey (KiDS-1000). To measure the shapes of the galaxies, we used two complementary algorithms, finding consistent IA measurements for the overlapping galaxy sample. The global significance of… ▽ More We constrain the luminosity and redshift dependence of the intrinsic alignment (IA) of a nearly volume-limited sample of luminous red galaxies selected from the fourth public data release of the Kilo-Degree Survey (KiDS-1000). To measure the shapes of the galaxies, we used two complementary algorithms, finding consistent IA measurements for the overlapping galaxy sample. The global significance of IA detection across our two independent luminous red galaxy samples, with our favoured method of shape estimation, is $\sim10.7σ$. We find no significant dependence with redshift of the IA signal in the range $0.2<z<0.8$, nor a dependence with luminosity below $L_r\lesssim 2.9 \times 10^{10} h^{-2} L_{r,\odot}$. Above this luminosity, however, we find that the IA signal increases as a power law, although our results are also compatible with linear growth within the current uncertainties. This behaviour motivates the use of a broken power law model when accounting for the luminosity dependence of IA contamination in cosmic shear studies. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: 21 pages, 14 figures. Accepted for publication in A&A

The Theory of Functional Connections: A journey from theory to application

Authors: Hunter Johnston

Abstract: The Theory of Functional Connections (TFC) is a general methodology for functional interpolation that can embed a set of user-specified linear constraints. The functionals derived from this method, called \emph{constrained expressions}, analytically satisfy the imposed constraints and can be leveraged to transform constrained optimization problems to unconstrained ones. By simplifying the optimiza… ▽ More The Theory of Functional Connections (TFC) is a general methodology for functional interpolation that can embed a set of user-specified linear constraints. The functionals derived from this method, called \emph{constrained expressions}, analytically satisfy the imposed constraints and can be leveraged to transform constrained optimization problems to unconstrained ones. By simplifying the optimization problem, this technique has been shown to produce a numerical scheme that is faster, more accurate, and robust to poor initialization. The content of this dissertation details the complete development of the Theory of Functional Connections. First, the seminal paper on the Theory of Functional Connections is discussed and motivates the discovery of a more general formulation of the constrained expressions. Leveraging this formulation, a rigorous structure of the constrained expression is produced with associated mathematical definitions, claims, and proofs. Furthermore, the second part of this dissertation explains how this technique can be used to solve ordinary differential equations providing a wide variety of examples compared to the state-of-the-art. The final part of this work focuses on unitizing the techniques and algorithms produced in the prior sections to explore the feasibility of using the Theory of Functional Connections to solve real-time optimal control problems, namely optimal landing problems. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 283 pages

Fast 2-impulse non-Keplerian orbit-transfer using the Theory of Functional Connections

Authors: Allan K. de Almeida Junior, Hunter Johnston, Carl Leake, Daniele Mortari

Abstract: This study applies a new approach, the Theory of Functional Connections (TFC), to solve the two-point boundary-value problem (TPBVP) in non-Keplerian orbit transfer. The perturbations considered are drag, solar radiation pressure, higher-order gravitational potential harmonic terms, and multiple bodies. The proposed approach is applied to Earth-to-Moon transfers, and obtains exact boundary conditi… ▽ More This study applies a new approach, the Theory of Functional Connections (TFC), to solve the two-point boundary-value problem (TPBVP) in non-Keplerian orbit transfer. The perturbations considered are drag, solar radiation pressure, higher-order gravitational potential harmonic terms, and multiple bodies. The proposed approach is applied to Earth-to-Moon transfers, and obtains exact boundary condition satisfaction and with very fast convergence. Thanks to this highly efficient approach, perturbed pork-chop plots of Earth-to-Moon transfers are generated, and individual analyses on the transfers' parameters are easily done at low computational costs. The minimum fuel analysis is provided in terms of the time of flight, thrust application points, and relative geometry of the Moon and Sun. The transfer costs obtained are in agreement with the literature's best solutions, and in some cases are even slightly better. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Journal ref: Eur. Phys. J. Plus 136, 223 (2021)

Bright galaxy sample in the Kilo-Degree Survey Data Release 4: selection, photometric redshifts, and physical properties

Authors: M. Bilicki, A. Dvornik, H. Hoekstra, A. H. Wright, N. E. Chisari, M. Vakili, M. Asgari, B. Giblin, C. Heymans, H. Hildebrandt, B. W. Holwerda, A. Hopkins, H. Johnston, A. Kannawadi, K. Kuijken, S. J. Nakoneczny, H. Y. Shan, A. Sonnenfeld, E. Valentijn

Abstract: We present a bright galaxy sample with accurate and precise photometric redshifts (photo-zs), selected using $ugriZYJHK_\mathrm{s}$ photometry from the Kilo-Degree Survey (KiDS) Data Release 4 (DR4). The highly pure and complete dataset is flux-limited at $r<20$ mag, covers $\sim1000$ deg$^2$, and contains about 1 million galaxies after artifact masking. We exploit the overlap with Galaxy And Mass… ▽ More We present a bright galaxy sample with accurate and precise photometric redshifts (photo-zs), selected using $ugriZYJHK_\mathrm{s}$ photometry from the Kilo-Degree Survey (KiDS) Data Release 4 (DR4). The highly pure and complete dataset is flux-limited at $r<20$ mag, covers $\sim1000$ deg$^2$, and contains about 1 million galaxies after artifact masking. We exploit the overlap with Galaxy And Mass Assembly (GAMA) spectroscopy as calibration to determine photo-zs with the supervised machine learning neural network algorithm implemented in the ANNz2 software. The photo-zs have mean error of $|\langle δz \rangle| \sim 5 \times 10^{-4}$ and low scatter (scaled mean absolute deviation of $\sim 0.018(1+z)$), both practically independent of the $r$-band magnitude and photo-z at $0.05 < z_\mathrm{phot} < 0.5$. Combined with the 9-band photometry, these allow us to estimate robust absolute magnitudes and stellar masses for the full sample. As a demonstration of the usefulness of these data we split the dataset into red and blue galaxies, use them as lenses and measure the weak gravitational lensing signal around them for five stellar mass bins. We fit a halo model to these high-precision measurements to constrain the stellar-mass--halo-mass relations for blue and red galaxies. We find that for high stellar mass ($M_\star>5\times 10^{11} M_\odot$), the red galaxies occupy dark matter halos that are much more massive than those occupied by blue galaxies with the same stellar mass. The data presented here are publicly released via the KiDS webpage at http://kids.strw.leidenuniv.nl/DR4/brightsample.php. △ Less

Submitted 21 July, 2021; v1 submitted 15 January, 2021; originally announced January 2021.

Comments: Matches the published version. Data available at http://kids.strw.leidenuniv.nl/DR4/brightsample.php

Journal ref: A&A 653, A82 (2021)

Organised Randoms: Learning and correcting for systematic galaxy clustering patterns in KiDS using self-organising maps

Authors: Harry Johnston, Angus H. Wright, Benjamin Joachimi, Maciej Bilicki, Nora Elisa Chisari, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Shahab Joudaki, Mohammadjavad Vakili

Abstract: We present a new method for the mitigation of observational systematic effects in angular galaxy clustering via corrective random galaxy catalogues. Real and synthetic galaxy data, from the Kilo Degree Survey's (KiDS) 4$^{\rm{th}}$ Data Release (KiDS-$1000$) and the Full-sky Lognormal Astro-fields Simulation Kit (FLASK) package respectively, are used to train self-organising maps (SOMs) to learn t… ▽ More We present a new method for the mitigation of observational systematic effects in angular galaxy clustering via corrective random galaxy catalogues. Real and synthetic galaxy data, from the Kilo Degree Survey's (KiDS) 4$^{\rm{th}}$ Data Release (KiDS-$1000$) and the Full-sky Lognormal Astro-fields Simulation Kit (FLASK) package respectively, are used to train self-organising maps (SOMs) to learn the multivariate relationships between observed galaxy number density and up to six systematic-tracer variables, including seeing, Galactic dust extinction, and Galactic stellar density. We then create `organised' randoms, i.e. random galaxy catalogues with spatially variable number densities, mimicking the learnt systematic density modes in the data. Using realistically biased mock data, we show that these organised randoms consistently subtract spurious density modes from the two-point angular correlation function $w(\vartheta)$, correcting biases of up to $12σ$ in the mean clustering amplitude to as low as $0.1σ$, over a high signal-to-noise angular range of 7-100 arcmin. Their performance is also validated for angular clustering cross-correlations in a bright, flux-limited subset of KiDS-$1000$, comparing against an analogous sample constructed from highly-complete spectroscopic redshift data. Each organised random catalogue object is a `clone' carrying the properties of a real galaxy, and is distributed throughout the survey footprint according to the parent galaxy's position in systematics-space. Thus, sub-sample randoms are readily derived from a single master random catalogue via the same selection as applied to the real galaxies. Our method is expected to improve in performance with increased survey area, galaxy number density, and systematic contamination, making organised randoms extremely promising for current and future clustering analyses of faint samples. △ Less

Submitted 21 April, 2021; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: 19 pages (7 appendix pages), 12 figures (8 appendix figures), accepted to A&A

Journal ref: A&A 648, A98 (2021)

A Flare-Type IV Burst Event from Proxima Centauri and Implications for Space Weather

Authors: Andrew Zic, Tara Murphy, Christene Lynch, George Heald, Emil Lenc, David L. Kaplan, Iver H. Cairns, David Coward, Bruce Gendre, Helen Johnston, Meredith MacGregor, Danny C. Price, Michael S. Wheatland

Abstract: Studies of solar radio bursts play an important role in understanding the dynamics and acceleration processes behind solar space weather events, and the influence of solar magnetic activity on solar system planets. Similar low-frequency bursts detected from active M-dwarfs are expected to probe their space weather environments and therefore the habitability of their planetary companions. Active M-… ▽ More Studies of solar radio bursts play an important role in understanding the dynamics and acceleration processes behind solar space weather events, and the influence of solar magnetic activity on solar system planets. Similar low-frequency bursts detected from active M-dwarfs are expected to probe their space weather environments and therefore the habitability of their planetary companions. Active M-dwarfs produce frequent, powerful flares which, along with radio emission, reveal conditions within their atmospheres. However, to date, only one candidate solar-like coherent radio burst has been identified from these stars, preventing robust observational constraints on their space weather environment. During simultaneous optical and radio monitoring of the nearby dM5.5e star Proxima Centauri, we detected a bright, long-duration optical flare, accompanied by a series of intense, coherent radio bursts. These detections include the first example of an interferometrically detected coherent stellar radio burst temporally coincident with a flare, strongly indicating a causal relationship between these transient events. The polarization and temporal structure of the trailing long-duration burst enable us to identify it as a type IV burst. This represents the most compelling detection of a solar-like radio burst from another star to date. Solar type IV bursts are strongly associated with space weather events such as coronal mass ejections and solar energetic particle events, suggesting that stellar type IV bursts may be used as a tracer of stellar coronal mass ejections. We discuss the implications of this event for the occurrence of coronal mass ejections from Proxima Cen and other active M-dwarfs. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 18 pages, 7 figures, 1 table. Published in the Astrophysical Journal

The PAU Survey: Intrinsic alignments and clustering of narrow-band photometric galaxies

Authors: Harry Johnston, Benjamin Joachimi, Peder Norberg, Henk Hoekstra, Martin Eriksen, Maria Cristina Fortuna, Giorgio Manzoni, Santiago Serrano, Malgorzata Siudek, Luca Tortorelli, Laura Cabayol, Jorge Carretero, Ricard Casas, Francisco Castander, Enrique Fernandez, Juan García-Bellido, Enrique Gaztanaga, Hendrik Hildebrandt, Ramon Miquel, Cristobal Padilla, Eusebio Sanchez, Ignacio Sevilla-Noarbe, Pau Tallada-Crespí

Abstract: We present the first measurements of the projected clustering and intrinsic alignments (IA) of galaxies observed by the Physics of the Accelerating Universe Survey (PAUS). With photometry in 40 narrow optical passbands ($450\rm{nm}-850\rm{nm}$), the quality of photometric redshift estimation is $σ_{z} \sim 0.01(1 + z)$ for galaxies in the $19\,\rm{deg}^{2}$ Canada-France-Hawaii Telescope Legacy Su… ▽ More We present the first measurements of the projected clustering and intrinsic alignments (IA) of galaxies observed by the Physics of the Accelerating Universe Survey (PAUS). With photometry in 40 narrow optical passbands ($450\rm{nm}-850\rm{nm}$), the quality of photometric redshift estimation is $σ_{z} \sim 0.01(1 + z)$ for galaxies in the $19\,\rm{deg}^{2}$ Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) W3 field, allowing us to measure the projected 3D clustering and IA for flux-limited, faint galaxies ($i < 22.5$) out to $z\sim0.8$. To measure two-point statistics, we developed, and tested with mock photometric redshift samples, `cloned' random galaxy catalogues which can reproduce data selection functions in 3D and account for photometric redshift errors. In our fiducial colour-split analysis, we made robust null detections of IA for blue galaxies and tentative detections of radial alignments for red galaxies ($\sim1-3σ$), over scales of $0.1-18\,h^{-1}\rm{Mpc}$. The galaxy clustering correlation functions in the PAUS samples are comparable to their counterparts in a spectroscopic population from the Galaxy and Mass Assembly survey, modulo the impact of photometric redshift uncertainty which tends to flatten the blue galaxy correlation function, whilst steepening that of red galaxies. We investigate the sensitivity of our correlation function measurements to choices in the random catalogue creation and the galaxy pair-binning along the line of sight, in preparation for an optimised analysis over the full PAUS area. △ Less

Submitted 1 March, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: 18 pages (7 appendix pages), 14 figures (6 appendix figures), accepted in A&A

Journal ref: A&A 646, A147 (2021)

An unconditional proof of the abelian equivariant Iwasawa main conjecture and applications

Authors: Henri Johnston, Andreas Nickel

Abstract: Let $p$ be an odd prime. We give an unconditional proof of the equivariant Iwasawa main conjecture for totally real fields for every admissible one-dimensional $p$-adic Lie extension whose Galois group has an abelian Sylow $p$-subgroup. Crucially, this result does not depend on the vanishing of any $μ$-invariant. As applications, we deduce the Coates-Sinnott conjecture away from its $2$-primary pa… ▽ More Let $p$ be an odd prime. We give an unconditional proof of the equivariant Iwasawa main conjecture for totally real fields for every admissible one-dimensional $p$-adic Lie extension whose Galois group has an abelian Sylow $p$-subgroup. Crucially, this result does not depend on the vanishing of any $μ$-invariant. As applications, we deduce the Coates-Sinnott conjecture away from its $2$-primary part and new cases of the equivariant Tamagawa number conjecture for Tate motives. △ Less

Submitted 30 November, 2021; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: 43 pages; large parts of article rewritten; proof of main result simplified and new appendix included

MSC Class: 11R23; 11R42

Searching for quicksand ideals in partially ordered sets

Authors: Alexas Iams, Hannah Johnston, Robert Muth

Abstract: We consider a combinatorial question about searching for an unknown ideal $μ$ within a known poset $λ$. Elements of $λ$ may be queried for membership in $μ$, but at most $k$ positive query results are permitted. The goal is to find a search strategy which guarantees a solution in a minimal total number $m_k(λ)$ of queries. We provide tight bounds for $m_k(λ)$, and construct optimal search strategi… ▽ More We consider a combinatorial question about searching for an unknown ideal $μ$ within a known poset $λ$. Elements of $λ$ may be queried for membership in $μ$, but at most $k$ positive query results are permitted. The goal is to find a search strategy which guarantees a solution in a minimal total number $m_k(λ)$ of queries. We provide tight bounds for $m_k(λ)$, and construct optimal search strategies for the case where $k=2$ and $λ$ is the product poset of totally ordered finite sets, one of which has cardinality not more than six. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Clustering of red-sequence galaxies in the fourth data release ofthe Kilo-Degree Survey

Authors: Mohammadjavad Vakili, Henk Hoekstra, Maciej Bilicki, Maria-Cristina Fortuna, Konrad Kuijken, Angus H. Wright, Marika Asgari, Michael Brown, Elisabeth Dombrovskij, Thomas Erben, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Harry Johnston, Shahab Joudaki, Arun Kannawadi

Abstract: We present a sample of luminous red-sequence galaxies to study the large-scale structure in the fourth data release of the Kilo-Degree Survey. The selected galaxies are defined by a red-sequence template, in the form of a data-driven model of the colour-magnitude relation conditioned on redshift. In this work, the red-sequence template is built using the broad-band optical+near infrared photometry… ▽ More We present a sample of luminous red-sequence galaxies to study the large-scale structure in the fourth data release of the Kilo-Degree Survey. The selected galaxies are defined by a red-sequence template, in the form of a data-driven model of the colour-magnitude relation conditioned on redshift. In this work, the red-sequence template is built using the broad-band optical+near infrared photometry of KiDS-VIKING and the overlapping spectroscopic data sets. The selection process involves estimating the red-sequence redshifts, assessing the purity of the sample, and estimating the underlying redshift distributions of redshift bins. After performing the selection, we mitigate the impact of survey properties on the observed number density of galaxies by assigning photometric weights to the galaxies. We measure the angular two-point correlation function of the red galaxies in four redshift bins, and constrain the large scale bias of our red-sequence sample assuming a fixed $Λ$CDM cosmology. We find consistent linear biases for two luminosity-threshold samples (dense and luminous). We find that our constraints are well characterized by the passive evolution model. △ Less

Submitted 30 August, 2020; originally announced August 2020.

Comments: submitted to A&A

Journal ref: A&A 675, A202 (2023)

Kinematic Modeling and Compliance Modulation of Redundant Manipulators Under Bracing Constraints

Authors: Garrison L. H. Johnston, Andrew L. Orekhov, Nabil Simaan

Abstract: Collaborative robots should ideally use low torque actuators for passive safety reasons. However, some applications require these collaborative robots to reach deep into confined spaces while assisting a human operator in physically demanding tasks. In this paper, we consider the use of in-situ collaborative robots (ISCRs) that balance the conflicting demands of passive safety dictating low torque… ▽ More Collaborative robots should ideally use low torque actuators for passive safety reasons. However, some applications require these collaborative robots to reach deep into confined spaces while assisting a human operator in physically demanding tasks. In this paper, we consider the use of in-situ collaborative robots (ISCRs) that balance the conflicting demands of passive safety dictating low torque actuation and the need to reach into deep confined spaces. We consider the judicious use of bracing as a possible solution to these conflicting demands and present a modeling framework that takes into account the constrained kinematics and the effect of bracing on the end-effector compliance. We then define a redundancy resolution framework that minimizes the directional compliance of the end-effector while maximizing end-effector dexterity. Kinematic simulation results show that the redundancy resolution strategy successfully decreases compliance and improves kinematic conditioning while satisfying the constraints imposed by the bracing task. Applications of this modeling framework can support future research on the choice of bracing locations and support the formation of an admittance control framework for collaborative control of ISCRs under bracing constraints. Such robots can benefit workers in the future by reducing the physiological burdens that contribute to musculoskeletal injury. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: Accepted for publication at the 2020 IEEE International Conference on Robotics and Automation (ICRA)

A successful search for intervening 21 cm HI absorption in galaxies at 0.4 < z <1.0 with the Australian Square Kilometre Array Pathfinder (ASKAP)

Authors: Elaine M. Sadler, Vanessa A. Moss, James R. Allison, Elizabeth K. Mahony, Matthew T. Whiting, Helen M. Johnston, Sara L. Ellison, Claudia del P. Lagos, Bärbel S. Koribalski

Abstract: We have used the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope to search for intervening 21 cm neutral hydrogen (HI) absorption along the line of sight to 53 bright radio continuum sources. Our observations are sensitive to HI column densities typical of Damped Lyman Alpha absorbers (DLAs) in cool gas with an HI spin temperature below about 300-500 K. The six-dish Boolardy E… ▽ More We have used the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope to search for intervening 21 cm neutral hydrogen (HI) absorption along the line of sight to 53 bright radio continuum sources. Our observations are sensitive to HI column densities typical of Damped Lyman Alpha absorbers (DLAs) in cool gas with an HI spin temperature below about 300-500 K. The six-dish Boolardy Engineering Test Array (BETA) and twelve-antenna Early Science array (ASKAP-12) covered a frequency range corresponding to redshift $0.4<z<1.0$ and $0.37<z<0.77$ respectively for the HI line. Fifty of the 53 radio sources observed have reliable optical redshifts, giving a total redshift path $Δz$ = 21.37. This was a spectroscopically-untargeted survey, with no prior assumptions about the location of the lines in redshift space. Four intervening HI lines were detected, two of them new. In each case, the estimated HI column density lies above the DLA limit for HI spin temperatures above 50-80 K, and we estimate a DLA number density at redshift $z\sim0.6$ of $n(z)=0.19\substack{+0.15 \\ -0.09}$. This value lies somewhat above the general trend of $n(z)$ with redshift seen in optical DLA studies. Although the current sample is small, it represents an important proof of concept for the much larger 21cm First Large Absorption Survey in HI (FLASH) project to be carried out with the full 36-antenna ASKAP telescope, probing a total redshift path $Δz\sim\,50,000$. △ Less

Submitted 7 August, 2020; v1 submitted 10 July, 2020; originally announced July 2020.

Comments: 21 pages, 11 figures. Accepted for publication in MNRAS

The Multivariate Theory of Functional Connections: Theory, Proofs, and Application in Partial Differential Equations

Authors: Carl Leake, Hunter Johnston, Daniele Mortari

Abstract: This article presents a reformulation of the Theory of Functional Connections: a general methodology for functional interpolation that can embed a set of user-specified linear constraints. The reformulation presented in this paper exploits the underlying functional structure presented in the seminal paper on the Theory of Functional Connections to ease the derivation of these interpolating functio… ▽ More This article presents a reformulation of the Theory of Functional Connections: a general methodology for functional interpolation that can embed a set of user-specified linear constraints. The reformulation presented in this paper exploits the underlying functional structure presented in the seminal paper on the Theory of Functional Connections to ease the derivation of these interpolating functionals--called constrained expressions--and provides rigorous terminology that lends itself to straightforward derivations of mathematical proofs regarding the properties of these constrained expressions. Furthermore, the extension of the technique to and proofs in n-dimensions is immediate through a recursive application of the univariate formulation. In all, the results of this reformulation are compared to prior work to highlight the novelty and mathematical convenience of using this approach. Finally, the methodology presented in this paper is applied to two partial differential equations with different boundary conditions, and, when data is available, the results are compared to state-of-the-art methods. △ Less

Submitted 6 August, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: 33 pages, 6 figures, 14 tables

MSC Class: 35D99; 65N35; 65M70

Journal ref: Mathematics 2020, 8, 1303

Arakelov class groups of random number fields

Authors: Alex Bartel, Henri Johnston, Hendrik W. Lenstra Jr

Abstract: The main purpose of the paper is to formulate a probabilistic model for Arakelov class groups in families of number fields, offering a correction to the Cohen--Lenstra--Martinet heuristic on ideal class groups. To that end, we show that Chinburg's Omega(3) conjecture implies tight restrictions on the Galois module structure of oriented Arakelov class groups. As a consequence, we construct a new in… ▽ More The main purpose of the paper is to formulate a probabilistic model for Arakelov class groups in families of number fields, offering a correction to the Cohen--Lenstra--Martinet heuristic on ideal class groups. To that end, we show that Chinburg's Omega(3) conjecture implies tight restrictions on the Galois module structure of oriented Arakelov class groups. As a consequence, we construct a new infinite series of counterexamples to the Cohen--Lenstra--Martinet heuristic, which have the novel feature that their Galois groups are non-abelian. △ Less

Submitted 27 March, 2024; v1 submitted 23 May, 2020; originally announced May 2020.

Comments: 22 pages; minor expository changes; to appear in Math. Ann

MSC Class: 11R29; 11R33; 11R65

Extreme Theory of Functional Connections: A Physics-Informed Neural Network Method for Solving Parametric Differential Equations

Authors: Enrico Schiassi, Carl Leake, Mario De Florio, Hunter Johnston, Roberto Furfaro, Daniele Mortari

Abstract: In this work we present a novel, accurate, and robust physics-informed method for solving problems involving parametric differential equations (DEs) called the Extreme Theory of Functional Connections, or X-TFC. The proposed method is a synergy of two recently developed frameworks for solving problems involving parametric DEs, 1) the Theory of Functional Connections, TFC, and the Physics-Informed… ▽ More In this work we present a novel, accurate, and robust physics-informed method for solving problems involving parametric differential equations (DEs) called the Extreme Theory of Functional Connections, or X-TFC. The proposed method is a synergy of two recently developed frameworks for solving problems involving parametric DEs, 1) the Theory of Functional Connections, TFC, and the Physics-Informed Neural Networks, PINN. Although this paper focuses on the solution of exact problems involving parametric DEs (i.e. problems where the modeling error is negligible) with known parameters, X-TFC can also be used for data-driven solutions and data-driven discovery of parametric DEs. In the proposed method, the latent solution of the parametric DEs is approximated by a TFC constrained expression that uses a Neural Network (NN) as the free-function. This approximate solution form always analytically satisfies the constraints of the DE, while maintaining a NN with unconstrained parameters, like the Deep-TFC method. X-TFC differs from PINN and Deep-TFC; whereas PINN and Deep-TFC use a deep-NN, X-TFC uses a single-layer NN, or more precisely, an Extreme Learning Machine, ELM. This choice is based on the properties of the ELM algorithm. In order to numerically validate the method, it was tested over a range of problems including the approximation of solutions to linear and non-linear ordinary DEs (ODEs), systems of ODEs (SODEs), and partial DEs (PDEs). Furthermore, a few of these problems are of interest in physics and engineering such as the Classic Emden-Fowler equation, the Radiative Transfer (RT) equation, and the Heat-Transfer (HT) equation. The results show that X-TFC achieves high accuracy with low computational time and thus it is comparable with the other state-of-the-art methods. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Comments: 28 pages, 12 figures, 6 tables

Recursive Star-Identification Algorithm using an Adaptive SVD-based Angular Velocity Estimator

Authors: Hunter Johnston, Carl Leake, Marcelino M. de Almeida, Daniele Mortari

Abstract: This paper describes an algorithm obtained by merging a recursive star identification algorithm with a recently developed adaptive SVD-based estimator of the angular velocity vector (QuateRA). In a recursive algorithm, the more accurate the angular velocity estimate, the quicker and more robust to noise the resultant recursive algorithm is. Hence, combining these two techniques produces an algorit… ▽ More This paper describes an algorithm obtained by merging a recursive star identification algorithm with a recently developed adaptive SVD-based estimator of the angular velocity vector (QuateRA). In a recursive algorithm, the more accurate the angular velocity estimate, the quicker and more robust to noise the resultant recursive algorithm is. Hence, combining these two techniques produces an algorithm capable of handling a variety of dynamics scenarios. The speed and robustness of the algorithm are highlighted in a selection of simulated scenarios. First, a speed comparison is made with the state-of-the-art lost-in-space star identification algorithm, Pyramid. This test shows that in the best case the algorithm is on average an order of magnitude faster than Pyramid. Next, the recursive algorithm is validated for a variety of dynamic cases including a ground-based "Stellar Compass" scenario, a satellite in geosynchronous orbit, a satellite during a re-orientation maneuver, and a satellite undergoing non-pure-spin dynamics. △ Less

Submitted 29 May, 2020; v1 submitted 25 March, 2020; originally announced March 2020.

Comments: 15 pages, 11 figures, 6 tables

Journal ref: Journal of Spacecraft and Rockets, 2021

The halo model as a versatile tool to predict intrinsic alignments

Authors: Maria Cristina Fortuna, Henk Hoekstra, Benjamin Joachimi, Harry Johnston, Nora Elisa Chisari, Christos Georgiou, Constance Mahony

Abstract: Intrinsic alignments (IAs) of galaxies are an important contaminant for cosmic shear studies, but the modelling is complicated by the dependence of the signal on the source galaxy sample. In this paper, we use the halo model formalism to capture this diversity and examine its implications for Stage-III and Stage-IV cosmic shear surveys. We account for the different IA signatures at large and small… ▽ More Intrinsic alignments (IAs) of galaxies are an important contaminant for cosmic shear studies, but the modelling is complicated by the dependence of the signal on the source galaxy sample. In this paper, we use the halo model formalism to capture this diversity and examine its implications for Stage-III and Stage-IV cosmic shear surveys. We account for the different IA signatures at large and small scales, as well for the different contributions from central/satellite and red/blue galaxies, and we use realistic mocks to account for the characteristics of the galaxy populations as a function of redshift. We inform our model using the most recent observational findings: we include a luminosity dependence at both large and small scales and a radial dependence of the signal within the halo. We predict the impact of the total IA signal on the lensing angular power spectra, including the current uncertainties from the IA best-fits to illustrate the range of possible impact on the lensing signal: the lack of constraints for fainter galaxies is the main source of uncertainty for our predictions of the IA signal. We investigate how well effective models with limited degrees of freedom can account for the complexity of the IA signal. Although these lead to negligible biases for Stage-III surveys, we find that, for Stage-IV surveys, it is essential to at least include an additional parameter to capture the redshift dependence. △ Less

Submitted 14 December, 2020; v1 submitted 5 March, 2020; originally announced March 2020.

Comments: 21 pages, 15 figure, accepted for publication in MNRAS

Least-squares Solutions of Eighth-order Boundary Value Problems using the Theory of Functional Connections

Authors: Hunter Johnston, Carl Leake, Daniele Mortari

Abstract: This paper shows how to obtain highly accurate solutions of eighth-order boundary-value problems of linear and nonlinear ordinary differential equations. The presented method is based on the Theory of Functional Connections, and is solved in two steps. First, the Theory of Functional Connections analytically embeds the differential equation constraints into a candidate function (called a… ▽ More This paper shows how to obtain highly accurate solutions of eighth-order boundary-value problems of linear and nonlinear ordinary differential equations. The presented method is based on the Theory of Functional Connections, and is solved in two steps. First, the Theory of Functional Connections analytically embeds the differential equation constraints into a candidate function (called a $constrained \, expression$) that contains a function that the user is free to choose. This expression always satisfies the constraints, no matter what the free function is. Second, the free-function is expanded as a linear combination of orthogonal basis functions with unknown coefficients. The constrained expression (and its derivatives) are then substituted into the eighth-order differential equation, transforming the problem into an unconstrained optimization problem where the coefficients in the linear combination of orthogonal basis functions are the optimization parameters. These parameters are then found by linear/nonlinear least-squares. The solution obtained from this method is a highly accurate analytical approximation of the true solution. Comparisons with alternative methods appearing in literature validate the proposed approach. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: 14 pages, 1 figure, 8 tables

MSC Class: 34K10; 34K28; 65D05; 65L10; 65L60

Journal ref: Mathematics 2020, 8(3), 397

Fuel-Efficient Powered Descent Guidance on Large Planetary Bodies via Theory of Functional Connections

Authors: Hunter Johnston, Enrico Schiassi, Roberto Furfaro, Daniele Mortari

Abstract: In this paper we present a new approach to solve the fuel-efficient powered descent guidance problem on large planetary bodies with no atmosphere (e.g. the Moon or Mars) using the recently developed Theory of Functional Connections. The problem is formulated using the indirect method which casts the optimal guidance problem as a system of nonlinear two-point boundary value problems. Using the Theo… ▽ More In this paper we present a new approach to solve the fuel-efficient powered descent guidance problem on large planetary bodies with no atmosphere (e.g. the Moon or Mars) using the recently developed Theory of Functional Connections. The problem is formulated using the indirect method which casts the optimal guidance problem as a system of nonlinear two-point boundary value problems. Using the Theory of Functional Connections, the problem constraints are analytically embedded into a "constrained expression," which maintains a free-function that is expanded using orthogonal polynomials with unknown coefficients. The constraints are satisfied regardless of the values of the unknown coefficients which convert the two-point boundary value problem into an unconstrained optimization problem. This process casts the solution into the admissible subspace of the problem and therefore simple numerical techniques can be used (i.e. in this paper a nonlinear least-squares method is used). In addition to the derivation of this technique, the method is validated in two scenarios and the results are compared to those obtained by the general purpose optimal control software, GPOPS-II. In general, the proposed technique produces solutions of $\mathcal{O}(10^{-10})$. Additionally, for the proposed test cases, it is reported that each individual TFC-based inner-loop iteration converges within 6 iterations, each iteration exhibiting a computational time between 72 and 81 milliseconds within the MATLAB legacy implementation. Consequently, the proposed methodology is potentially suitable for on-board generation of optimal trajectories in real-time. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: 17 pages, 10 figures, 6 tables

MSC Class: 49M99

Journal ref: The Journal of the Astronautical Sciences (2020)

Molybdenum dichalcogenide cathodes for aluminium-ion batteries

Authors: Shalini Divya, James H. Johnston, Thomas Nann

Abstract: Many successful battery electrodes are based on 2D-layered materials. We have studied aluminium-ion batteries using molybdenum dichalcogenides: \ce{MoS2}, \ce{MoSe2} and MoSSe as active cathode materials. The batteries showed clear discharge voltage plateaus in the ranges 1.6 - 1.4 V for \ce{MoS2} and \ce{MoSe2}, and 0.6 - 0.5 V for MoSSe. \ce{MoS2} and \ce{MoSe2} have similar crystal structures,… ▽ More Many successful battery electrodes are based on 2D-layered materials. We have studied aluminium-ion batteries using molybdenum dichalcogenides: \ce{MoS2}, \ce{MoSe2} and MoSSe as active cathode materials. The batteries showed clear discharge voltage plateaus in the ranges 1.6 - 1.4 V for \ce{MoS2} and \ce{MoSe2}, and 0.6 - 0.5 V for MoSSe. \ce{MoS2} and \ce{MoSe2} have similar crystal structures, interestingly we found that \ce{MoSe2} performed better than \ce{MoS2}. MoSSe exhibited a higher specific capacity over \ce{MoS2} and \ce{MoSe2}, but the energy density was lower than \ce{MoSe2} at a current rate of 40 mA g$^{-1}$. \ce{MoSe2} cells recorded a discharge capacity of $\sim$ 110 mAh g$^{-1}$ with an average potential in the range of 2.0 - 1.8 V and 1.5 - 0.8 V during discharge. The cells were stable at 100 mA g$^{-1}$ for over 200 cycles with 90\% coulombic efficiency. △ Less

Submitted 30 March, 2020; v1 submitted 22 December, 2019; originally announced December 2019.

Comments: Energy Technology, 2020

Least-squares solutions of boundary-value problems in hybrid systems

Authors: Hunter Johnston, Daniele Mortari

Abstract: This paper looks to apply the mathematical framework of the Theory of Functional Connections to the solution of boundary-value problems arising from hybrid systems. The Theory of Functional Connections is a technique to derive constrained expressions which are analytical expressions with embedded constraints. These expressions are particularly suitable to transform a large class of constrained opt… ▽ More This paper looks to apply the mathematical framework of the Theory of Functional Connections to the solution of boundary-value problems arising from hybrid systems. The Theory of Functional Connections is a technique to derive constrained expressions which are analytical expressions with embedded constraints. These expressions are particularly suitable to transform a large class of constrained optimization problems into unconstrained problems. The initial and most useful application of this technique is in the solution of differential equations where the problem can be posed as an unconstrained optimization problem and solved with simple numerical techniques (i.e. least-squares). A hybrid system is simply a sequence of different differential equations. The approach developed in this work derives an analytical constrained expression for the entire range of a hybrid system, enforcing both the boundary conditions as well as the continuity conditions across the sequence of differential equations. This reduces the searched solution space of the hybrid system to only admissible solutions. The transformation allows for a least-squares solution of the sequence for linear differential equations and a iterative least-squares solution for nonlinear differential equations. Lastly, the general formulation for "n" segments is developed and validation is provided through numerical tests for three differential equation sequences: linear/linear, linear/nonlinear, and nonlinear/nonlinear. The accuracy level obtained are all at machine-error, which is consistent with the accuracy experienced in past studies on the application of the Theory of Functional Connections to solve single ordinary differential equations. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Comments: 15 pages, 7 figures, 5 tables

MSC Class: 65L10

Journal ref: Journal of Computational and Applied Mathematics (2021)

Ionisation of the atomic gas in redshifted radio sources

Authors: S. J. Curran, R. W. Hunstead, H. M. Johnston, M. T. Whiting, E. M. Sadler, J. R. Allison, R. Athreya

Abstract: We report the results of a survey for HI 21-cm absorption at z < 0.4 in a new sample of radio sources with the Giant Metrewave Radio Telescope. Of the 11 sources for which there are good data, we obtain zero detections, where four are expected upon accounting for the ionising photon rates and sensitivity. Adding these to the previously published values, we confirm that the non-detection of 21-cm a… ▽ More We report the results of a survey for HI 21-cm absorption at z < 0.4 in a new sample of radio sources with the Giant Metrewave Radio Telescope. Of the 11 sources for which there are good data, we obtain zero detections, where four are expected upon accounting for the ionising photon rates and sensitivity. Adding these to the previously published values, we confirm that the non-detection of 21-cm absorption in active sources at high redshift is due to photo-ionisation of the gas rather than excitation by 21-cm photons (significant at 6.09 sigma and 2.90 sigma, respectively). We also confirm a strong correlation between the absorption strength and the reddening of the source, suggesting that dust plays a significant role in shielding the gas from the ambient ultra-violet field. An anti-correlation between the 21-cm detection rate and the radio turnover frequency is also found, which runs contrary to what is expected on the basis that the higher the turnover frequency, the more compact the source. It is, however, consistent with the hypothesis that the turnover frequency is related to the electron density, supported by a correlation between the turnover frequency and ionising photon rate. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: Accepted by MNRAS

Analytically Embedding Differential Equation Constraints into Least Squares Support Vector Machines using the Theory of Functional Connections

Authors: Carl Leake, Hunter Johnston, Lidia Smith, Daniele Mortari

Abstract: Differential equations (DEs) are used as numerical models to describe physical phenomena throughout the field of engineering and science, including heat and fluid flow, structural bending, and systems dynamics. While there are many other techniques for finding approximate solutions to these equations, this paper looks to compare the application of the Theory of Functional Connections (TFC) with on… ▽ More Differential equations (DEs) are used as numerical models to describe physical phenomena throughout the field of engineering and science, including heat and fluid flow, structural bending, and systems dynamics. While there are many other techniques for finding approximate solutions to these equations, this paper looks to compare the application of the Theory of Functional Connections (TFC) with one based on least-squares support vector machines (LS-SVM). The TFC method uses a constrained expression, an expression that always satisfies the DE constraints, which transforms the process of solving a DE into solving an unconstrained optimization problem that is ultimately solved via least-squares (LS). In addition to individual analysis, the two methods are merged into a new methodology, called constrained SVMs (CSVM), by incorporating the LS-SVM method into the TFC framework to solve unconstrained problems. Numerical tests are conducted on four sample problems: One first order linear ordinary differential equation (ODE), one first order nonlinear ODE, one second order linear ODE, and one two-dimensional linear partial differential equation (PDE). Using the LS-SVM method as a benchmark, a speed comparison is made for all the problems by timing the training period, and an accuracy comparison is made using the maximum error and mean squared error on the training and test sets. In general, TFC is shown to be slightly faster (by an order of magnitude or less) and more accurate (by multiple orders of magnitude) than the LS-SVM and CSVM approaches. △ Less

Submitted 9 October, 2019; v1 submitted 13 December, 2018; originally announced December 2018.

Comments: 22 pages, 8 Figures, 12 tables

Journal ref: Mach. Learn. Knowl. Extr. 2019, 1(4), 1058-1083

KiDS+GAMA: Intrinsic alignment model constraints for current and future weak lensing cosmology

Authors: Harry Johnston, Christos Georgiou, Benjamin Joachimi, Henk Hoekstra, Nora Elisa Chisari, Daniel Farrow, Maria Cristina Fortuna, Catherine Heymans, Shahab Joudaki, Konrad Kuijken, Angus Wright

Abstract: We directly constrain the non-linear alignment (NLA) model of intrinsic galaxy alignments, analysing the most representative and complete flux-limited sample of spectroscopic galaxies available for cosmic shear surveys. We measure the projected galaxy position-intrinsic shear correlations and the projected galaxy clustering signal using high-resolution imaging from the Kilo Degree Survey (KiDS) ov… ▽ More We directly constrain the non-linear alignment (NLA) model of intrinsic galaxy alignments, analysing the most representative and complete flux-limited sample of spectroscopic galaxies available for cosmic shear surveys. We measure the projected galaxy position-intrinsic shear correlations and the projected galaxy clustering signal using high-resolution imaging from the Kilo Degree Survey (KiDS) overlapping with the GAMA spectroscopic survey, and data from the Sloan Digital Sky Survey. Separating samples by colour, we make no significant detection of blue galaxy alignments, constraining the blue galaxy NLA amplitude $A_{\textrm{IA}}^{\textrm{B}}=0.21^{+0.37}_{-0.36}$ to be consistent with zero. We make robust detections ($\sim9σ$) for red galaxies, with $A_{\textrm{IA}}^{\textrm{R}}=3.18^{+0.47}_{-0.46}$, corresponding to a net radial alignment with the galaxy density field, and we find no evidence for any scaling of alignments with galaxy luminosity. We provide informative priors for current and future weak lensing surveys, an improvement over de facto wide priors that allow for unrealistic levels of intrinsic alignment contamination. For a colour-split cosmic shear analysis of the final KiDS survey area, we forecast that our priors will improve the constraining power on $S_{8}$ and the dark energy equation of state $w_{0}$, by up to $62\%$ and $51\%$, respectively. Our results indicate, however, that the modelling of red/blue-split galaxy alignments may be insufficient to describe samples with variable central/satellite galaxy fractions. △ Less

Submitted 20 February, 2019; v1 submitted 23 November, 2018; originally announced November 2018.

Comments: 27 pages (incl. 7 appendix pages), 10 figures, accepted by A&A

Journal ref: A&A 624, A30 (2019)

The dependence of intrinsic alignment of galaxies on wavelength using KiDS and GAMA

Authors: Christos Georgiou, Harry Johnston, Henk Hoekstra, Massimo Viola, Konrad Kuijken, Benjamin Joachimi, Nora Elisa Chisari, Hendrik Hildebrandt, Arun Kannawadi

Abstract: The outer regions of galaxies are more susceptible to the tidal interactions that lead to intrinsic alignments of galaxies. The resulting alignment signal may therefore depend on the passband if the colours of galaxies vary spatially. To quantify this, we measured the shapes of galaxies with spectroscopic redshifts from the GAMA survey using deep gri imaging data from the KiloDegree Survey. The pe… ▽ More The outer regions of galaxies are more susceptible to the tidal interactions that lead to intrinsic alignments of galaxies. The resulting alignment signal may therefore depend on the passband if the colours of galaxies vary spatially. To quantify this, we measured the shapes of galaxies with spectroscopic redshifts from the GAMA survey using deep gri imaging data from the KiloDegree Survey. The performance of the moment-based shape measurement algorithm DEIMOS was assessed using dedicated image simulations, which showed that the ellipticities could be determined with an accuracy better than 1% in all bands. Additional tests for potential systematic errors did not reveal any issues. We measure a significant difference of the alignment signal between the g,r and i-band observations. This difference exceeds the amplitude of the linear alignment model on scales below 2 Mpc/h. Separating the sample into central/satellite and red/blue galaxies, we find that that the difference is dominated by red satellite galaxies. △ Less

Submitted 21 January, 2019; v1 submitted 10 September, 2018; originally announced September 2018.

Comments: 16 pages, 13 figures, accepted, to appear in A&A

Journal ref: A&A 622, A90 (2019)

Computing isomorphisms between lattices

Authors: Tommy Hofmann, Henri Johnston

Abstract: Let K be a number field, let A be a finite dimensional semisimple K-algebra and let Lambda be an O_K-order in A. It was shown in previous work that, under certain hypotheses on A, there exists an algorithm that for a given (left) Lambda-lattice X either computes a free basis of X over Lambda or shows that X is not free over Lambda. In the present article, we generalise this by showing that, under… ▽ More Let K be a number field, let A be a finite dimensional semisimple K-algebra and let Lambda be an O_K-order in A. It was shown in previous work that, under certain hypotheses on A, there exists an algorithm that for a given (left) Lambda-lattice X either computes a free basis of X over Lambda or shows that X is not free over Lambda. In the present article, we generalise this by showing that, under weaker hypotheses on A, there exists an algorithm that for two given Lambda-lattices X and Y either computes an isomorphism X -> Y or determines that X and Y are not isomorphic. The algorithm is implemented in Magma for A=Q[G], Lambda=Z[G] and Lambda-lattices X and Y contained in Q[G], where G is a finite group satisfying certain hypotheses. This is used to investigate the Galois module structure of rings of integers and ambiguous ideals of tamely ramified Galois extensions of Q with Galois group isomorphic to Q_8 x C_2, the direct product of the quaternion group of order 8 and the cyclic group of order 2. △ Less

Submitted 2 March, 2020; v1 submitted 22 June, 2018; originally announced June 2018.

Comments: 30 pages; v3 revised and accepted version to appear in Mathematics of Computation; v2 has many minor corrections with additional explanation in section 10

MSC Class: 11R33; 11Y40; 16Z05

Multiband counterparts of two eclipsing ultraluminous X-ray sources in M 51

Authors: R. Urquhart, R. Soria, H. M. Johnston, M. W. Pakull, C. Motch, A. Schwope, J. C. A. Miller-Jones, G. E. Anderson

Abstract: We present the discovery and interpretation of ionized nebulae around two ultraluminous X-ray sources in M 51; both sources share the rare property of showing X-ray eclipses by their companion stars, and are therefore prime targets for follow-up studies. Using archival Hubble Space Telescope images, we found an elongated, 100-pc-long emission-line structure associated with one X-ray source (CXOM51… ▽ More We present the discovery and interpretation of ionized nebulae around two ultraluminous X-ray sources in M 51; both sources share the rare property of showing X-ray eclipses by their companion stars, and are therefore prime targets for follow-up studies. Using archival Hubble Space Telescope images, we found an elongated, 100-pc-long emission-line structure associated with one X-ray source (CXOM51 J132940.0+471237; ULX-1 for simplicity), and a more circular, ionized nebula at the location of the second source (CXOM51 J132939.5+471244; ULX-2 for simplicity). We observed both nebulae with the Large Binocular Telescope's Multi-Object Double Spectrograph. From our analysis of the optical spectra, we argue that the gas in the ULX-1 bubble is shock-ionized, consistent with the effect of a jet with a kinetic power of $\approx$2 $\times 10^{39}$ erg s$^{-1}$. Additional X-ray photo-ionization may also be present, to explain the strength of high-ionization lines such as He II $λ4686$ and [Ne V] $λ3426$. On the other hand, the emission lines from the ULX-2 bubble are typical for photoionization by normal O stars suggesting that the nebula is actually an H II region not physically related to the ULX but is simply a chance alignment. From archival Very Large Array data, we also detect spatially extended, steep-spectrum radio emission at the location of the ULX-1 bubble (consistent with its jet origin), but no radio counterpart for ULX-2 (consistent with the lack of shock-ionized gas around that source). △ Less

Submitted 5 January, 2018; originally announced January 2018.

Comments: 18 pages, 10 figures, 2 tables. Accepted for publication in MNRAS

Assimilation of nearly turbulent Rayleigh-Bénard flow through vorticity or local circulation measurements: a computational study

Authors: Aseel Farhat, Hans Johnston, Michael S. Jolly, Edriss S. Titi

Abstract: We introduce a continuous (downscaling) data assimilation algorithm for the 2D Bénard convection problem using vorticity or local circulation measurements only. In this algorithm, a nudging term is added to the vorticity equation to constrain the model. Our numerical results indicate that the approximate solution of the algorithm is converging to the unknown reference solution (vorticity and tempe… ▽ More We introduce a continuous (downscaling) data assimilation algorithm for the 2D Bénard convection problem using vorticity or local circulation measurements only. In this algorithm, a nudging term is added to the vorticity equation to constrain the model. Our numerical results indicate that the approximate solution of the algorithm is converging to the unknown reference solution (vorticity and temperature) corresponding to the measurements of the 2D Bénard convection problem when only spatial coarse-grain measurements of vorticity are assimilated. Moreover, this convergence is realized using data which is much more coarse than the resolution needed to satisfy rigorous analytical estimates. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Further Observational Evidence for a Critical Ionising Luminosity in Active Galaxies

Authors: S. J. Curran, R. W. Hunstead, H. M. Johnston, M. T. Whiting, E. M. Sadler, J. R. Allison, C. Bignell

Abstract: We report the results of a survey for HI 21-cm absorption at redshifts of z > 2.6 in a new sample of radio sources with the Green Bank and Giant Metrewave Radio Telescopes. From a total of 25 targets, we report zero detections in the 16 for which optical depth limits could be obtained. Based upon the detection rate for z > 0.1 associated absorption, we would expect approximately four detections. O… ▽ More We report the results of a survey for HI 21-cm absorption at redshifts of z > 2.6 in a new sample of radio sources with the Green Bank and Giant Metrewave Radio Telescopes. From a total of 25 targets, we report zero detections in the 16 for which optical depth limits could be obtained. Based upon the detection rate for z > 0.1 associated absorption, we would expect approximately four detections. Of the 11 which have previously not been searched, there is sufficient source-frame optical/ultra-violet photometry to determine the ionising photon rate for four. Adding these to the literature, the hypothesis that there is a critical rate of logQ = 56 ionising photons per second is now significant at ~7 sigma. This reaffirms our assertion that searching z > 3 active galaxies for which optical redshifts are available selects sources in which the ultra-violet luminosity is sufficient to ionise all of the neutral gas in the host galaxy. △ Less

Submitted 7 July, 2017; v1 submitted 20 June, 2017; originally announced June 2017.

Comments: Accepted by MNRAS

Galaxy And Mass Assembly (GAMA): The environments of high- and low- excitation radio galaxies

Authors: J. H. Y. Ching, S. M. Croom, E. M. Sadler, A. S. G. Robotham, S. Brough, I. K. Baldry, J. Bland-Hawthorn, M. Colless, S. P. Driver, B. W. Holwerda, A. M. Hopkins, M. J. Jarvis, H. M. Johnston, L. S. Kelvin, J. Liske, J. Loveday, P. Norberg, M. B. Pracy, O. Steele, D. Thomas, L. Wang

Abstract: We study the environments of low- and high- excitation radio galaxies (LERGs and HERGs respectively) in the redshift range $0.01 < z < 0.4$, using a sample of 399 radio galaxies and environmental measurements from the Galaxy And Mass Assembly (GAMA) survey. In our analysis we use the fifth nearest neighbour density ($Σ_{5}$) and the GAMA galaxy groups catalogue (G3Cv6) and construct control sample… ▽ More We study the environments of low- and high- excitation radio galaxies (LERGs and HERGs respectively) in the redshift range $0.01 < z < 0.4$, using a sample of 399 radio galaxies and environmental measurements from the Galaxy And Mass Assembly (GAMA) survey. In our analysis we use the fifth nearest neighbour density ($Σ_{5}$) and the GAMA galaxy groups catalogue (G3Cv6) and construct control samples of galaxies matched in {\update stellar mass and colour} to the radio-detected sample. We find that LERGs and HERGs exist in different environments and that this difference is dependent on radio luminosity. High-luminosity LERGs ($L_{\rm NVSS} \gtrsim 10^{24}$ W Hz$^{-1}$) lie in much denser environments than a matched radio-quiet control sample (about three times as dense, as measured by $Σ_{5}$), and are more likely to be members of galaxy groups ($82^{+5}_{-7}$ percent of LERGs are in GAMA groups, compared to $58^{+3}_{-3}$ percent of the control sample). In contrast, the environments of the HERGs and lower luminosity LERGs are indistinguishable from that of a matched control sample. Our results imply that high-luminosity LERGs lie in more massive haloes than non-radio galaxies of similar stellar mass and colour, in agreement with earlier studies (Wake et al. 2008; Donoso et al. 2010). When we control for the preference of LERGs to be found in groups, both high- and low- luminosity LERGs are found in higher-mass haloes ($\sim 0.2$ dex; at least 97 percent significant) than the non-radio control sample. △ Less

Submitted 12 May, 2017; originally announced May 2017.

Comments: 18 pages, 11 figures. Accepted for publication by MNRAS

On the p-adic Stark conjecture at s=1 and applications

Authors: Henri Johnston, Andreas Nickel

Abstract: Let E/F be a finite Galois extension of totally real number fields and let p be a prime. The `p-adic Stark conjecture at s=1' relates the leading terms at s=1 of p-adic Artin L-functions to those of the complex Artin L-functions attached to E/F. We prove this conjecture unconditionally when E/Q is abelian. We also show that for certain non-abelian extensions E/F the p-adic Stark conjecture at s=1… ▽ More Let E/F be a finite Galois extension of totally real number fields and let p be a prime. The `p-adic Stark conjecture at s=1' relates the leading terms at s=1 of p-adic Artin L-functions to those of the complex Artin L-functions attached to E/F. We prove this conjecture unconditionally when E/Q is abelian. We also show that for certain non-abelian extensions E/F the p-adic Stark conjecture at s=1 is implied by Leopoldt's conjecture for E at p. Moreover, we prove that for a fixed prime p, the p-adic Stark conjecture at s=1 for E/F implies Stark's conjecture at s=1 for E/F. This leads to a `prime-by-prime' descent theorem for the `equivariant Tamagawa number conjecture' (ETNC) for Tate motives at s=1. As an application of these results, we provide strong new evidence for special cases of the ETNC for Tate motives and the closely related `leading term conjectures' at s=0 and s=1. △ Less

Submitted 14 October, 2019; v1 submitted 20 March, 2017; originally announced March 2017.

Comments: 35 pages; v5 accepted version to appear in Journal of the London Mathematical Society; v4 minor revisions and appendix by Tommy Hofmann, Henri Johnston and Andreas Nickel added; v3 minor revisions; v2 minor revisions; comments welcome

MSC Class: 11R23; 11R42

Optical Spectra of Candidate International Celestial Reference Frame (ICRF) Flat-Spectrum Radio Sources. III

Authors: O. Titov, T. Pursimo, H. M. Johnston, L. M. Stanford, R. W. Hunstead, D. L. Jauncey, K. A. Zenere

Abstract: In extending our spectroscopic program, which targets sources drawn from the International Celestial Reference Frame (ICRF) Catalog, we have obtained spectra for about 160 compact flat-spectrum radio sources and determined redshifts for 112 quasars and radio galaxies. A further 14 sources with featureless spectra have been classified as BL Lac objects. Spectra were obtained at three telescopes: th… ▽ More In extending our spectroscopic program, which targets sources drawn from the International Celestial Reference Frame (ICRF) Catalog, we have obtained spectra for about 160 compact flat-spectrum radio sources and determined redshifts for 112 quasars and radio galaxies. A further 14 sources with featureless spectra have been classified as BL Lac objects. Spectra were obtained at three telescopes: the 3.58-m ESO New Technology Telescope (NTT), and the two 8.2-m Gemini telescopes in Hawaii and Chile. While most of the sources are powerful quasars, a significant fraction of radio galaxies is also included from the list of non-defining ICRF radio sources. △ Less

Submitted 21 February, 2017; originally announced February 2017.

Comments: Accepted for publication in AJ

The Large Area Radio Galaxy Evolution Spectroscopic Survey (LARGESS): Survey design, data catalogue and GAMA/WiggleZ spectroscopy

Authors: John H. Y. Ching, Elaine M. Sadler, Scott M. Croom, Helen M. Johnston, Michael B. Pracy, Warrick J. Couch, A. M. Hopkins, Russell J. Jurek, K. A. Pimbblet

Abstract: We present the Large Area Radio Galaxy Evolution Spectroscopic Survey (LARGESS), a spectroscopic catalogue of radio sources designed to include the full range of radio AGN populations out to redshift z = 0.8. The catalogue covers roughly 800 square degrees of sky, and provides optical identifications for 19,179 radio sources from the 1.4 GHz Faint Images of the Radio Sky at Twenty-cm (FIRST) surve… ▽ More We present the Large Area Radio Galaxy Evolution Spectroscopic Survey (LARGESS), a spectroscopic catalogue of radio sources designed to include the full range of radio AGN populations out to redshift z = 0.8. The catalogue covers roughly 800 square degrees of sky, and provides optical identifications for 19,179 radio sources from the 1.4 GHz Faint Images of the Radio Sky at Twenty-cm (FIRST) survey down to an optical magnitude limit of i_mod < 20.5 in Sloan Digital Sky Survey (SDSS) images. Both galaxies and point-like objects are included, and no colour cuts are applied. In collaboration with the WiggleZ and Galaxy And Mass Assembly (GAMA) spectroscopic survey teams, we have obtained new spectra for over 5,000 objects in the LARGESS sample. Combining these new spectra with data from earlier surveys provides spectroscopic data for 12,329 radio sources in the survey area, of which 10,856 have reliable redshifts. 85% of the LARGESS spectroscopic sample are radio AGN (median redshift z = 0.44), and 15% are nearby star-forming galaxies (median z = 0.08). Low-excitation radio galaxies (LERGs) comprise the majority (83%) of LARGESS radio AGN at z < 0.8, with 12% being high-excitation radio galaxies (HERGs) and 5% radio-loud QSOs. Unlike the more homogeneous LERG and QSO sub-populations, HERGs are a heterogeneous class of objects with relatively blue optical colours and a wide dispersion in mid-infrared colours. This is consistent with a picture in which most HERGs are hosted by galaxies with recent or ongoing star formation as well as a classical accretion disk. △ Less

Submitted 18 September, 2016; originally announced September 2016.

Comments: 35 pages, 14 tables, 21 figures, accepted for publication in MNRAS. The full data table (Table 10) will appear online. This table can be downloaded in advance of publication from http://tinyurl.com/jhyc-tab10

The nature of the companion star in Circinus X-1

Authors: Helen M. Johnston, Roberto Soria, Joel Gibson

Abstract: We present optical spectra and images of the X-ray binary Circinus X-1. The optical light curve of Cir X-1 is strongly variable, changing in brightness by 1.2 magnitudes in the space of four days. The shape of the light curve is consistent with that seen in the 1980s, when the X-ray and radio counterparts of the source were at least ten times as bright as they are currently. We detect strong, vari… ▽ More We present optical spectra and images of the X-ray binary Circinus X-1. The optical light curve of Cir X-1 is strongly variable, changing in brightness by 1.2 magnitudes in the space of four days. The shape of the light curve is consistent with that seen in the 1980s, when the X-ray and radio counterparts of the source were at least ten times as bright as they are currently. We detect strong, variable H-alpha emission lines, consisting of multiple components which vary with orbital phase. We estimate the extinction to the source from the strength of the diffuse interstellar bands and the Balmer decrement; the two methods give A_V = 7.6 +/- 0.6 mag and A_V > 9.1 mag respectively. The optical light curve can be modelled as arising from irradiation of the companion star by the central X-ray source, where a low-temperature star fills its Roche lobe in an orbit of moderate eccentricity (e ~ 0.4). We suggest that the companion star is over-luminous and under-dense, due to the impact of the supernova which occurred less than 5000 yr ago. △ Less

Submitted 11 November, 2015; originally announced November 2015.

Comments: MNRAS, accepted

Kinematic signatures of AGN feedback in moderately powerful radio galaxies at z~2 observed with SINFONI

Authors: C. Collet, N. P. H. Nesvadba, C. De Breuck, M. D. Lehnert, P. Best, J. J. Bryant, R. Hunstead, D. Dicken, H. Johnston

Abstract: Most successful galaxy formation scenarios now postulate that the intense star formation in massive, high-redshift galaxies during their major growth period was truncated when powerful AGNs launched galaxy-wide outflows of gas that removed large parts of the interstellar medium. The most powerful radio galaxies at z~2 show clear signatures of such winds, but are too rare to be good representatives… ▽ More Most successful galaxy formation scenarios now postulate that the intense star formation in massive, high-redshift galaxies during their major growth period was truncated when powerful AGNs launched galaxy-wide outflows of gas that removed large parts of the interstellar medium. The most powerful radio galaxies at z~2 show clear signatures of such winds, but are too rare to be good representatives of a generic phase in the evolution of all massive galaxies at high redshift. Here we present SINFONI imaging spectroscopy of 12 radio galaxies at z~2 that are intermediate between the most powerful radio and vigorous starburst galaxies in radio power, and common enough to represent a generic phase in the early evolution of massive galaxies. The kinematic properties are diverse, with regular velocity gradients with amplitudes of Delta v=200-400 km s^-1 as in rotating disks as well as irregular kinematics with multiple velocity jumps of a few 100 km s^-1. Line widths are generally high, typically around FWHM=800 km s^-1, consistent with wind velocities in hydrodynamic models. A broad H-alpha line in one target implies a black hole mass of a few 10^9 M$_sun. The ratio of line widths, sigma, to bulk velocity, v, is so large that even the gas in galaxies with regular velocity fields is unlikely to be gravitationally bound. It is unclear, however, whether the large line widths are due to turbulence or unresolved, local outflows as are sometimes observed at low redshifts. Comparison of the kinetic energy with the energy supply from the AGN through jet and radiation pressure suggests that the radio source still plays a dominant role for feedback, consistent with low-redshift radio-loud quasars. △ Less

Submitted 22 October, 2015; originally announced October 2015.

Comments: A&A in press

On the non-abelian Brumer-Stark conjecture and the equivariant Iwasawa main conjecture

Authors: Henri Johnston, Andreas Nickel

Abstract: We show that for an odd prime p, the p-primary parts of refinements of the (imprimitive) non-abelian Brumer and Brumer-Stark conjectures are implied by the equivariant Iwasawa main conjecture (EIMC) for totally real fields. Crucially, this result does not depend on the vanishing of the relevant Iwasawa mu-invariant. In combination with the authors' previous work on the EIMC, this leads to uncondit… ▽ More We show that for an odd prime p, the p-primary parts of refinements of the (imprimitive) non-abelian Brumer and Brumer-Stark conjectures are implied by the equivariant Iwasawa main conjecture (EIMC) for totally real fields. Crucially, this result does not depend on the vanishing of the relevant Iwasawa mu-invariant. In combination with the authors' previous work on the EIMC, this leads to unconditional proofs of the non-abelian Brumer and Brumer-Stark conjectures in many new cases. △ Less

Submitted 3 August, 2018; v1 submitted 1 September, 2015; originally announced September 2015.

Comments: 33 pages; to appear in Mathematische Zeitschrift; v3 many minor updates including new title; v2 some cohomological arguments simplified; v1 is a revised version of the second half of arXiv:1408.4934v2

MSC Class: 11R23; 11R42

Defying jet-gas alignment in two radio galaxies at z~2 with extended light profiles: Similarities to brightest cluster galaxies

Authors: C. Collet, N. P. H. Nesvadba, C. De Breuck, M. D. Lehnert, P. Best, J. J. Bryant, D. Dicken, H. Johnston, R. Hunstead, D. Wylezalek

Abstract: We report the detection of extended warm ionized gas in two powerful high-redshift radio galaxies, NVSS J210626-314003 at z=2.10 and TXS 2353-003 at z=1.49, that does not appear to be associated with the radio jets. This is contrary to what would be expected from the alignment effect, a characteristic feature of distant, powerful radio galaxies at z> 0.6. The gas also has smaller velocity gradient… ▽ More We report the detection of extended warm ionized gas in two powerful high-redshift radio galaxies, NVSS J210626-314003 at z=2.10 and TXS 2353-003 at z=1.49, that does not appear to be associated with the radio jets. This is contrary to what would be expected from the alignment effect, a characteristic feature of distant, powerful radio galaxies at z> 0.6. The gas also has smaller velocity gradients and line widths than most other high-z radio galaxies with similar data. Both galaxies are part of a systematic study of 50 high-redshift radio galaxies with SINFONI, and are the only two that are characterized by the presence of high surface-brightness gas not associated with the jet axis and by the absence of such gas aligned with the jet. Both galaxies are spatially resolved with ISAAC broadband imaging covering the rest-frame R band, and have extended wings that cannot be attributed to line contamination. We argue that the gas and stellar properties of these galaxies are more akin to gas-rich brightest cluster galaxies in cool-core clusters than the general population of high-redshift radio galaxies at z>2. In support of this interpretation, one of our sources, TXS 2353-003, for which we have Hα narrowband imaging, is associated with an overdensity of candidate Hα emitters by a factor of 8 relative to the field at z=1.5. We discuss possible scenarios of the evolutionary state of these galaxies and the nature of their emission line gas within the context of cyclical AGN feedback. △ Less

Submitted 11 June, 2015; originally announced June 2015.

Comments: A&A in press