-
Volume enclosed by a flux surface
Authors:
Robert S. MacKay
Abstract:
The paper describes ways that the computation of the volume enclosed by an invariant torus (flux surface) for a magnetic field can be reduced from a 3D integral to a 2D integral.
The paper describes ways that the computation of the volume enclosed by an invariant torus (flux surface) for a magnetic field can be reduced from a 3D integral to a 2D integral.
△ Less
Submitted 18 June, 2024;
originally announced June 2024.
-
Characterization of admissible quasisymmetries
Authors:
J. W. Burby,
N. Kallinikos,
R. S. MacKay,
D. Perrella,
D. Pfefferlé
Abstract:
We solve "half" the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined n…
▽ More
We solve "half" the problem of finding three-dimensional quasisymmetric magnetic fields that do not necessarily satisfy force balance. This involves determining which hidden symmetries are admissible as quasisymmetries, and then showing explicitly how to construct quasisymmetric magnetic fields given an admissible symmetry. The admissibility conditions take the form of a system of overdetermined nonlinear partial differential equations involving second derivatives of the symmetry's infinitesimal generator.
△ Less
Submitted 3 July, 2024; v1 submitted 5 March, 2024;
originally announced March 2024.
-
Regions without flux surfaces of given class for magnetic fields in toroidal geometry
Authors:
Nikos Kallinikos,
Robert S. MacKay,
David Martinez-del-Rio
Abstract:
A Converse KAM method for 3D vector fields, establishing regions through which pass no invariant 2-tori transverse to a given direction field, is tested on some helical perturbations of an axisymmetric magnetic field in toroidal geometry. It finds regions corresponding to magnetic islands and chaos for the fieldline flow. Minimization of these regions is proposed as a tool to help in the design of…
▽ More
A Converse KAM method for 3D vector fields, establishing regions through which pass no invariant 2-tori transverse to a given direction field, is tested on some helical perturbations of an axisymmetric magnetic field in toroidal geometry. It finds regions corresponding to magnetic islands and chaos for the fieldline flow. Minimization of these regions is proposed as a tool to help in the design of plasma confinement devices of tokamak and stellarator type.
△ Less
Submitted 6 July, 2023; v1 submitted 19 April, 2023;
originally announced April 2023.
-
Isodrastic Magnetic fields for suppressing transitions in guiding-centre motion
Authors:
J. W. Burby,
R. S. MacKay,
S. Naik
Abstract:
In a magnetic field, transitions between classes of guiding-centre motion can lead to cross-field diffusion and escape. We say a magnetic field is isodrastic if guiding centres make no transitions between classes of motion. Therefore, this is an important ideal for enhancing confinement. First, we present a weak formulation, based on the longitudinal adiabatic invariant, generalising omnigenity. T…
▽ More
In a magnetic field, transitions between classes of guiding-centre motion can lead to cross-field diffusion and escape. We say a magnetic field is isodrastic if guiding centres make no transitions between classes of motion. Therefore, this is an important ideal for enhancing confinement. First, we present a weak formulation, based on the longitudinal adiabatic invariant, generalising omnigenity. To demonstrate that isodrasticity is strictly more general than omnigenity, we construct weakly isodrastic mirror fields that are not omnigenous. Then we present a strong formulation that is exact for guiding-centre motion. We develop a first-order treatment of the strong version via a Melnikov function and show that it recovers the weak version. The theory provides quantification of deviations from isodrasticity that can be used as objective functions in optimal design. The theory is illustrated on some simple examples.
△ Less
Submitted 8 March, 2023; v1 submitted 23 November, 2022;
originally announced November 2022.
-
A new stochastic framework for ship capsizing
Authors:
Manuela L. Bujorianu,
Robert S. MacKay,
Tobias Grafke,
Shibabrat Naik,
Evangelos Boulougouris
Abstract:
We present a new stochastic framework for studying ship capsize. It is a synthesis of two strands of transition state theory. The first is an extension of deterministic transition state theory to dissipative non-autonomous systems, together with a probability distribution over the forcing functions. The second is stochastic reachability and large deviation theory for transition paths in Markovian…
▽ More
We present a new stochastic framework for studying ship capsize. It is a synthesis of two strands of transition state theory. The first is an extension of deterministic transition state theory to dissipative non-autonomous systems, together with a probability distribution over the forcing functions. The second is stochastic reachability and large deviation theory for transition paths in Markovian systems. In future work we aim to bring these together to make a tool for predicting capsize rate in different stochastic sea states, suggesting control strategies and improving designs.
△ Less
Submitted 12 May, 2021;
originally announced May 2021.
-
Approximate symmetries of guiding-centre motion
Authors:
Joshua W. Burby,
Nikos Kallinikos,
Robert S. MacKay
Abstract:
Quasisymmetry builds a third invariant for charged-particle motion besides energy and magnetic moment. We address quasisymmetry at the level of approximate symmetries of first-order guiding-centre motion. We find that the conditions to leading order are the same as for exact quasisymmetry if one insists that the symmetry is purely spatial. We also generalise to allow for approximate phase-space sy…
▽ More
Quasisymmetry builds a third invariant for charged-particle motion besides energy and magnetic moment. We address quasisymmetry at the level of approximate symmetries of first-order guiding-centre motion. We find that the conditions to leading order are the same as for exact quasisymmetry if one insists that the symmetry is purely spatial. We also generalise to allow for approximate phase-space symmetries, and derive weaker conditions. The latter recover "weak quasisymmetry" as a subcase, thus we prove it is spatial only to leading order, but also that it implies the existence of a wider class of independent approximate conserved quantities. Finally, we demonstrate that magnetohydrostatics imposes quasisymmetry to leading order.
△ Less
Submitted 18 October, 2020;
originally announced October 2020.
-
Generalized Grad-Shafranov equation for non-axisymmetric MHD equilibria
Authors:
J. W. Burby,
N. Kallinikos,
R. S. MacKay
Abstract:
The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad-Shafranov equation, which is a single elliptic PDE in two space dimensions. By revealing a hidden symmetry, we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away fro…
▽ More
The structure of static MHD equilibria that admit continuous families of Euclidean symmetries is well understood. Such field configurations are governed by the classical Grad-Shafranov equation, which is a single elliptic PDE in two space dimensions. By revealing a hidden symmetry, we show that in fact all smooth solutions of the equilibrium equations with non-vanishing pressure gradients away from the magnetic axis satisfy a generalization of the Grad-Shafranov equation. In contrast to solutions of the classical Grad-Shafranov equation, solutions of he generalized equation are not automatically equilibria, but instead only satisfy force balance averaged over the one-parameter hidden symmetry. We then explain how the generalized Grad-Shafranov equation can be used to reformulate the problem of finding exact three-dimensional smooth solutions of the equilibrium equations as finding an optimal volume-preserving symmetry.
△ Less
Submitted 30 September, 2020; v1 submitted 27 May, 2020;
originally announced May 2020.
-
How directed is a directed network?
Authors:
R. S. MacKay,
S. Johnson,
B. Sansom
Abstract:
The trophic levels of nodes in directed networks can reveal their functional properties. Moreover, the trophic coherence of a network, defined in terms of trophic levels, is related to properties such as cycle structure, stability and percolation. The standard definition of trophic levels, however, borrowed from ecology, suffers from drawbacks such as requiring source nodes, which limit its applic…
▽ More
The trophic levels of nodes in directed networks can reveal their functional properties. Moreover, the trophic coherence of a network, defined in terms of trophic levels, is related to properties such as cycle structure, stability and percolation. The standard definition of trophic levels, however, borrowed from ecology, suffers from drawbacks such as requiring source nodes, which limit its applicability. Here we propose a simple new definition of trophic level that can be computed on any directed network. We demonstrate how the method can identify node function in examples including ecosystems, supply chain networks, gene expression, and global language networks. We also explore how trophic levels and coherence relate to other topological properties, such as non-normality and cycle structure, and show that our method reveals the extent to which the edges in a directed network are aligned in a global direction.
△ Less
Submitted 15 January, 2020;
originally announced January 2020.
-
Differential forms for plasma physics
Authors:
Robert S. MacKay
Abstract:
Differential forms provide a coordinate-free way to express many quantities and relations in mathematical physics. In particular, they are useful in plasma physics. This tutorial gives a guide so that you can read the plasma physics literature that uses them and apply them yourself.
Differential forms provide a coordinate-free way to express many quantities and relations in mathematical physics. In particular, they are useful in plasma physics. This tutorial gives a guide so that you can read the plasma physics literature that uses them and apply them yourself.
△ Less
Submitted 17 December, 2019;
originally announced December 2019.
-
Calibration with confidence: A principled method for panel assessment
Authors:
Robert S. MacKay,
Ralph Kenna,
Robert J. Low,
Sarah Parker
Abstract:
Frequently, a set of objects has to be evaluated by a panel of assessors, but not every object is assessed by every assessor. A problem facing such panels is how to take into account different standards amongst panel members and varying levels of confidence in their scores. Here, a mathematically-based algorithm is developed to calibrate the scores of such assessors, addressing both of these issue…
▽ More
Frequently, a set of objects has to be evaluated by a panel of assessors, but not every object is assessed by every assessor. A problem facing such panels is how to take into account different standards amongst panel members and varying levels of confidence in their scores. Here, a mathematically-based algorithm is developed to calibrate the scores of such assessors, addressing both of these issues. The algorithm is based on the connectivity of the graph of assessors and objects evaluated, incorporating declared confidences as weights on its edges. If the graph is sufficiently well connected, relative standards can be inferred by comparing how assessors rate objects they assess in common, weighted by the levels of confidence of each assessment. By removing these biases, "true" values are inferred for all the objects. Reliability estimates for the resulting values are obtained. The algorithm is tested in two case studies, one by computer simulation and another based on realistic evaluation data. The process is compared to the simple averaging procedure in widespread use, and to Fisher's additive incomplete block analysis. It is anticipated that the algorithm will prove useful in a wide variety of situations such as evaluation of the quality of research submitted to national assessment exercises; appraisal of grant proposals submitted to funding panels; ranking of job applicants; and judgement of performances on degree courses wherein candidates can choose from lists of options.
△ Less
Submitted 13 February, 2017; v1 submitted 28 March, 2015;
originally announced April 2015.
-
Morse bifurcations of transition states in bimolecular reactions
Authors:
Robert S. MacKay,
Dayal C. Strub
Abstract:
The transition states and dividing surfaces used to find rate constants for bimolecular reactions are shown to undergo qualitative changes, known as Morse bifurcations, and to exist for a large range of energies, not just immediately above the critical energy for first connection between reactants and products. Specifically, we consider capture between two molecules and the associated transition s…
▽ More
The transition states and dividing surfaces used to find rate constants for bimolecular reactions are shown to undergo qualitative changes, known as Morse bifurcations, and to exist for a large range of energies, not just immediately above the critical energy for first connection between reactants and products. Specifically, we consider capture between two molecules and the associated transition states for the case of non-zero angular momentum and general attitudes. The capture between an atom and a diatom, and then a general molecule are presented, providing concrete examples of Morse bifurcations of transition states and dividing surfaces. The reduction of the $n$-body systems representing the reactions is discussed and reviewed with comments on the difficulties associated with choosing appropriate charts and the global geometry of the reduced spaces.
△ Less
Submitted 11 November, 2015; v1 submitted 1 January, 2015;
originally announced January 2015.
-
Opportunities at the Mathematics/Future Cities Interface
Authors:
Peter Grindrod,
Desmond J. Higham,
Robert S. MacKay
Abstract:
We make the case for mathematicians and statisticians to stake their claim in the fast-moving and high-impact research field that is becoming known as Future Cities. After assessing the Future Cities arena, we provide some illustrative challenges where mathematical scientists can make an impact.
We make the case for mathematicians and statisticians to stake their claim in the fast-moving and high-impact research field that is becoming known as Future Cities. After assessing the Future Cities arena, we provide some illustrative challenges where mathematical scientists can make an impact.
△ Less
Submitted 5 September, 2014;
originally announced September 2014.
-
Frontiers of chaotic advection
Authors:
Hassan Aref,
John R. Blake,
Marko Budišić,
Silvana S. S. Cardoso,
Julyan H. E. Cartwright,
Herman J. H. Clercx,
Kamal El Omari,
Ulrike Feudel,
Ramin Golestanian,
Emmanuelle Gouillart,
GertJan F. van Heijst,
Tatyana S. Krasnopolskaya,
Yves Le Guer,
Robert S. MacKay,
Vyacheslav V. Meleshko,
Guy Metcalfe,
Igor Mezić,
Alessandro P. S. de Moura,
Oreste Piro,
Michel F. M. Speetjens,
Rob Sturman,
Jean-Luc Thiffeault,
Idan Tuval
Abstract:
This work reviews the present position of and surveys future perspectives in the physics of chaotic advection: the field that emerged three decades ago at the intersection of fluid mechanics and nonlinear dynamics, which encompasses a range of applications with length scales ranging from micrometers to hundreds of kilometers, including systems as diverse as mixing and thermal processing of viscous…
▽ More
This work reviews the present position of and surveys future perspectives in the physics of chaotic advection: the field that emerged three decades ago at the intersection of fluid mechanics and nonlinear dynamics, which encompasses a range of applications with length scales ranging from micrometers to hundreds of kilometers, including systems as diverse as mixing and thermal processing of viscous fluids, microfluidics, biological flows, and oceanographic and atmospheric flows.
△ Less
Submitted 14 June, 2017; v1 submitted 12 March, 2014;
originally announced March 2014.
-
Bifurcations of transition states: Morse bifurcations
Authors:
Robert S MacKay,
Dayal C Strub
Abstract:
A transition state for a Hamiltonian system is a closed, invariant, oriented, codimension-2 submanifold of an energy-level that can be spanned by two compact codimension-1 surfaces of unidirectional flux whose union, called a dividing surface, locally separates the energy-level into two components and has no local recrossings. For this to happen robustly to all smooth perturbations, the transition…
▽ More
A transition state for a Hamiltonian system is a closed, invariant, oriented, codimension-2 submanifold of an energy-level that can be spanned by two compact codimension-1 surfaces of unidirectional flux whose union, called a dividing surface, locally separates the energy-level into two components and has no local recrossings. For this to happen robustly to all smooth perturbations, the transition state must be normally hyperbolic. The dividing surface then has locally minimal geometric flux through it, giving an upper bound on the rate of transport in either direction. Transition states diffeomorphic to $\mathbb S^{2m-3}$ are known to exist for energies just above any index-1 critical point of a Hamiltonian of $m$ degrees of freedom, with dividing surfaces $\mathbb S^{2m-2}$. The question addressed here is what qualitative changes in the transition state, and consequently the dividing surface, may occur as the energy or other parameters are varied? We find that there is a class of systems for which the transition state becomes singular and then regains normal hyperbolicity with a change in diffeomorphism class. These are Morse bifurcations. Various examples are considered. Firstly, some simple examples in which transition states connect or disconnect, and the dividing surface may become a torus or other. Then, we show how sequences of Morse bifurcations producing various interesting forms of transition state and dividing surface are present in reacting systems, by considering a hypothetical class of bimolecular reactions in gas phase.
△ Less
Submitted 17 April, 2014; v1 submitted 8 May, 2013;
originally announced May 2013.
-
Challenges in Complex Systems Science
Authors:
Maxi San Miguel,
Jeffrey H. Johnson,
Janos Kertesz,
Kimmo Kaski,
Albert Díaz-Guilera,
Robert S. MacKay,
Vittorio Loreto,
Peter Erdi,
Dirk Helbing
Abstract:
FuturICT foundations are social science, complex systems science, and ICT. The main concerns and challenges in the science of complex systems in the context of FuturICT are laid out in this paper with special emphasis on the Complex Systems route to Social Sciences. This include complex systems having: many heterogeneous interacting parts; multiple scales; complicated transition laws; unexpected o…
▽ More
FuturICT foundations are social science, complex systems science, and ICT. The main concerns and challenges in the science of complex systems in the context of FuturICT are laid out in this paper with special emphasis on the Complex Systems route to Social Sciences. This include complex systems having: many heterogeneous interacting parts; multiple scales; complicated transition laws; unexpected or unpredicted emergence; sensitive dependence on initial conditions; path-dependent dynamics; networked hierarchical connectivities; interaction of autonomous agents; self-organisation; non-equilibrium dynamics; combinatorial explosion; adaptivity to changing environments; co-evolving subsystems; ill-defined boundaries; and multilevel dynamics. In this context, science is seen as the process of abstracting the dynamics of systems from data. This presents many challenges including: data gathering by large-scale experiment, participatory sensing and social computation, managing huge distributed dynamic and heterogeneous databases; moving from data to dynamical models, going beyond correlations to cause-effect relationships, understanding the relationship between simple and comprehensive models with appropriate choices of variables, ensemble modeling and data assimilation, modeling systems of systems of systems with many levels between micro and macro; and formulating new approaches to prediction, forecasting, and risk, especially in systems that can reflect on and change their behaviour in response to predictions, and systems whose apparently predictable behaviour is disrupted by apparently unpredictable rare or extreme events. These challenges are part of the FuturICT agenda.
△ Less
Submitted 22 April, 2012;
originally announced April 2012.
-
Nonlinear structures and thermodynamic instabilities in a one-dimensional lattice system
Authors:
Nikos Theodorakopoulos,
Michel Peyrard,
Robert S. MacKay
Abstract:
The equilibrium states of the discrete Peyrard-Bishop Hamiltonian with one end fixed are computed exactly from the two-dimensional nonlinear Morse map. These exact nonlinear structures are interpreted as domain walls (DW), interpolating between bound and unbound segments of the chain. The free energy of the DWs is calculated to leading order beyond the Gaussian approximation. Thermodynamic insta…
▽ More
The equilibrium states of the discrete Peyrard-Bishop Hamiltonian with one end fixed are computed exactly from the two-dimensional nonlinear Morse map. These exact nonlinear structures are interpreted as domain walls (DW), interpolating between bound and unbound segments of the chain. The free energy of the DWs is calculated to leading order beyond the Gaussian approximation. Thermodynamic instabilities (e.g. DNA unzipping and/or thermal denaturation) can be understood in terms of DW formation.
△ Less
Submitted 8 November, 2004;
originally announced November 2004.
