-
How to Build the Virtual Cell with Artificial Intelligence: Priorities and Opportunities
Authors:
Charlotte Bunne,
Yusuf Roohani,
Yanay Rosen,
Ankit Gupta,
Xikun Zhang,
Marcel Roed,
Theo Alexandrov,
Mohammed AlQuraishi,
Patricia Brennan,
Daniel B. Burkhardt,
Andrea Califano,
Jonah Cool,
Abby F. Dernburg,
Kirsty Ewing,
Emily B. Fox,
Matthias Haury,
Amy E. Herr,
Eric Horvitz,
Patrick D. Hsu,
Viren Jain,
Gregory R. Johnson,
Thomas Kalil,
David R. Kelley,
Shana O. Kelley,
Anna Kreshuk
, et al. (17 additional authors not shown)
Abstract:
The cell is arguably the most fundamental unit of life and is central to understanding biology. Accurate modeling of cells is important for this understanding as well as for determining the root causes of disease. Recent advances in artificial intelligence (AI), combined with the ability to generate large-scale experimental data, present novel opportunities to model cells. Here we propose a vision…
▽ More
The cell is arguably the most fundamental unit of life and is central to understanding biology. Accurate modeling of cells is important for this understanding as well as for determining the root causes of disease. Recent advances in artificial intelligence (AI), combined with the ability to generate large-scale experimental data, present novel opportunities to model cells. Here we propose a vision of leveraging advances in AI to construct virtual cells, high-fidelity simulations of cells and cellular systems under different conditions that are directly learned from biological data across measurements and scales. We discuss desired capabilities of such AI Virtual Cells, including generating universal representations of biological entities across scales, and facilitating interpretable in silico experiments to predict and understand their behavior using virtual instruments. We further address the challenges, opportunities and requirements to realize this vision including data needs, evaluation strategies, and community standards and engagement to ensure biological accuracy and broad utility. We envision a future where AI Virtual Cells help identify new drug targets, predict cellular responses to perturbations, as well as scale hypothesis exploration. With open science collaborations across the biomedical ecosystem that includes academia, philanthropy, and the biopharma and AI industries, a comprehensive predictive understanding of cell mechanisms and interactions has come into reach.
△ Less
Submitted 14 October, 2024; v1 submitted 17 September, 2024;
originally announced September 2024.
-
Thermal ion kinetic effects and Landau damping in fishbone modes
Authors:
Chang Liu,
Stephen C. Jardin,
Jian Bao,
Nikolai Gorelenkov,
Dylan P. Brennan,
James Yang,
Mario Podesta
Abstract:
The kinetic-MHD hybrid simulation approach for macroscopic instabilities in plasmas can be extended to include the kinetic effects of both thermal ions and energetic ions. The new coupling scheme includes synchronization of density and parallel velocity between thermal ions and MHD, in addition to pressure coupling, to ensure the quasineutrality condition and avoid numerical errors. The new approa…
▽ More
The kinetic-MHD hybrid simulation approach for macroscopic instabilities in plasmas can be extended to include the kinetic effects of both thermal ions and energetic ions. The new coupling scheme includes synchronization of density and parallel velocity between thermal ions and MHD, in addition to pressure coupling, to ensure the quasineutrality condition and avoid numerical errors. The new approach has been implemented in the kinetic-MHD code M3D-C1-K, and was used to study the thermal ion kinetic effects and Landau damping in fishbone modes in both DIII-D and NSTX. It is found that the thermal ion kinetic effects can cause an increase of the frequencies of the non-resonant $n=1$ fishbone modes driven by energetic particles for $q_\mathrm{min}>1$, and Landau damping can provide additional stabilization effects. A nonlinear simulation for $n=1$ fishbone mode in NSTX is also performed, and the perturbation on magnetic flux surfaces and the transport of energetic particles are calculated.
△ Less
Submitted 7 June, 2022;
originally announced June 2022.
-
Interpretable Machine Learning Classifiers for Brain Tumour Survival Prediction
Authors:
Colleen E. Charlton,
Michael Tin Chung Poon,
Paul M. Brennan,
Jacques D. Fleuriot
Abstract:
Prediction of survival in patients diagnosed with a brain tumour is challenging because of heterogeneous tumour behaviours and responses to treatment. Better estimations of prognosis would support treatment planning and patient support. Advances in machine learning have informed development of clinical predictive models, but their integration into clinical practice is almost non-existent. One reas…
▽ More
Prediction of survival in patients diagnosed with a brain tumour is challenging because of heterogeneous tumour behaviours and responses to treatment. Better estimations of prognosis would support treatment planning and patient support. Advances in machine learning have informed development of clinical predictive models, but their integration into clinical practice is almost non-existent. One reasons for this is the lack of interpretability of models. In this paper, we use a novel brain tumour dataset to compare two interpretable rule list models against popular machine learning approaches for brain tumour survival prediction. All models are quantitatively evaluated using standard performance metrics. The rule lists are also qualitatively assessed for their interpretability and clinical utility. The interpretability of the black box machine learning models is evaluated using two post-hoc explanation techniques, LIME and SHAP. Our results show that the rule lists were only slightly outperformed by the black box models. We demonstrate that rule list algorithms produced simple decision lists that align with clinical expertise. By comparison, post-hoc interpretability methods applied to black box models may produce unreliable explanations of local model predictions. Model interpretability is essential for understanding differences in predictive performance and for integration into clinical practice.
△ Less
Submitted 17 June, 2021;
originally announced June 2021.
-
Exascale Landau collision operator in the Cuda programming model applied to thermal quench plasmas
Authors:
M. F. Adams,
D. P. Brennan,
M. G. Knepley,
P. Wang
Abstract:
Collisional processes are critical in the understanding of non-Maxwellian plasmas. The Landau form of the Fokker-Planck equation is the gold standard for modeling collisions in most plasmas, however O(N^2) work complexity inhibits its widespread use. We show that with advanced numerical methods and GPU hardware this cost can be effectively mitigated. This paper extends previous work on a conservat…
▽ More
Collisional processes are critical in the understanding of non-Maxwellian plasmas. The Landau form of the Fokker-Planck equation is the gold standard for modeling collisions in most plasmas, however O(N^2) work complexity inhibits its widespread use. We show that with advanced numerical methods and GPU hardware this cost can be effectively mitigated. This paper extends previous work on a conservative, high order accurate, finite element discretization with adaptive mesh refinement of the Landau operator, with extensions to GPU hardware and implementations in both the CUDA and Kokkos programming languages. This work focuses on the Landau kernels and on NVIDIA hardware, however preliminary results on AMD and Fujitsu/ARM hardware, as well as end-to-end performance of a velocity space model of a plasma thermal quench, are also presented. Both the fully implicit Landau time integrator and the plasma thermal quench model are publicly available in PETSc (Portable, Extensible, Toolkit for Scientific computing).
△ Less
Submitted 18 May, 2022; v1 submitted 7 April, 2021;
originally announced April 2021.
-
Compressional Alfvén eigenmodes excited by runaway electrons
Authors:
Chang Liu,
Dylan P. Brennan,
Andrey Lvovskiy,
Carlos Paz-Soldan,
Eric D. Fredrickson,
Amitava Bhattacharjee
Abstract:
Compressional Alfvén eigenmodes (CAE) driven by energetic ions have been observed in magnetic fusion experiments. In this paper, we show that the modes can also be driven by runaway electrons formed in post-disruption plasma, which may explain kinetic instabilities observed in DIII-D disruption experiments with massive gas injection. The mode-structure is calculated, as are the frequencies which a…
▽ More
Compressional Alfvén eigenmodes (CAE) driven by energetic ions have been observed in magnetic fusion experiments. In this paper, we show that the modes can also be driven by runaway electrons formed in post-disruption plasma, which may explain kinetic instabilities observed in DIII-D disruption experiments with massive gas injection. The mode-structure is calculated, as are the frequencies which are in agreement with experimental observations. Using a runaway electron distribution function obtained from a kinetic simulation, the mode growth rates are calculated and found to exceed the collisional damping rate when the runaway electron density exceeds a threshold value. The excitation of CAE poses a new possible approach to mitigate seed runaway electrons during the current quench and surpassing the avalanche.
△ Less
Submitted 1 October, 2020; v1 submitted 24 September, 2020;
originally announced September 2020.
-
Canonical Degrees of Cohen-Macaulay Rings and Modules: a Survey
Authors:
J. P. Brennan,
L. Ghezzi,
J. Hong,
L. Hutson,
W. V. Vasconcelos
Abstract:
The aim of this survey is to discuss invariants of Cohen-Macaulay local rings that admit a canonical module. Attached to each such ring R with a canonical ideal C, there are integers--the type of R, the reduction number of C--that provide valuable metrics to express the deviation of R from being a Gorenstein ring. We enlarge this list with other integers--the roots of R and several canonical degre…
▽ More
The aim of this survey is to discuss invariants of Cohen-Macaulay local rings that admit a canonical module. Attached to each such ring R with a canonical ideal C, there are integers--the type of R, the reduction number of C--that provide valuable metrics to express the deviation of R from being a Gorenstein ring. We enlarge this list with other integers--the roots of R and several canonical degrees. The latter are multiplicity based functions of the Rees algebra of C. We give a uniform presentation of three degrees arising from common roots. Finally we experiment with ways to extend one of these degrees to rings where C is not necessarily an ideal.
△ Less
Submitted 24 June, 2020;
originally announced June 2020.
-
Structure and overstability of resistive modes with runaway electrons
Authors:
Chang Liu,
Chen Zhao,
Stephen C. Jardin,
Amitava Bhattacharjee,
Dylan P. Brennan,
Nathanial M. Ferraro
Abstract:
We investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge wh…
▽ More
We investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge whose properties depend on the ratio of the Alfvén velocity to the runaway electron convection speed. Due to the contribution from this layer, both the tearing mode and kink mode will have a real frequency in addition to a growth rate. The derived dispersion relation has been compared with numerical results using both a simplified eigenvalue calculation and a M3D-C1 linear simulation, and good agreement is found in both cases.
△ Less
Submitted 9 September, 2020; v1 submitted 15 June, 2020;
originally announced June 2020.
-
Dark-field signal extraction in propagation-based phase-contrast imaging
Authors:
T. E. Gureyev,
D. M. Paganin,
B. Arhatari,
S. T. Taba,
S. Lewis,
P. C. Brennan,
H. M. Quiney
Abstract:
A method for extracting the dark-field signal in propagation-based phase-contrast imaging is proposed. In the case of objects consisting predominantly of a single material, or several different materials with similar ratios of the real decrement to the imaginary part of the complex refractive index, the proposed method requires a single image for extraction of the dark-field signal in two-dimensio…
▽ More
A method for extracting the dark-field signal in propagation-based phase-contrast imaging is proposed. In the case of objects consisting predominantly of a single material, or several different materials with similar ratios of the real decrement to the imaginary part of the complex refractive index, the proposed method requires a single image for extraction of the dark-field signal in two-dimensional projection imaging. In the case of three-dimensional tomographic imaging, the method needs only one image to be collected at each projection angle. A preliminary example demonstrates that this method can improve the visualization of microcalcifications in propagation-based X-ray breast cancer imaging. It is suggested that the proposed approach may be useful in other forms of biomedical imaging, where it can help one to obtain additional small-angle scattering information without increasing the radiation dose to the sample.
△ Less
Submitted 2 April, 2020; v1 submitted 27 March, 2020;
originally announced March 2020.
-
Nonlinear error field response in the presence of plasma rotation and real frequencies due to favorable curvature
Authors:
Cihan Akcay,
John M. Finn,
Andrew C. Cole,
Dylan P. Brennan
Abstract:
We present nonlinear resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TM's), stabilized by pressure gradient and favorable curvature. Favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index $Δ'$. A cylinder with hollow pressure is used…
▽ More
We present nonlinear resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TM's), stabilized by pressure gradient and favorable curvature. Favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index $Δ'$. A cylinder with hollow pressure is used to model the toroidal favorable curvature. Linear simulations with rotation and an error field $\tildeψ_w$ show, in agreement with analytic results, that the peak reconnected flux occurs for a rotation rate near the TM phase velocity. Nonlinear simulations with small $\tildeψ_w$ show that the real frequency and stabilization by favorable average curvature are masked by a nonlinear effect that occurs for very thin islands: flattening of the pressure across the island, mainly due to sound wave propagation. This flattening causes the disappearance of real frequencies destabilization of the mode, allowing it to grow to large amplitude similar to a $β=0$ unstable TM. The flattening of the current for larger islands saturates the mode nonlinearly. In the post-saturation phase, the interaction of the error field with the destabilized spontaneous tearing mode, which rotates with the plasma, leads to oscillations in the Maxwell torque and therefore modulations in the plasma rotation. The islands also rotate with modulated phase velocity, undergoing small-amplitude oscillations due to these modulations. We also present a quasilinear model with an unstable spontaneous TM and error fields, showing that the superposition of these fields results in similar oscillations.
△ Less
Submitted 6 January, 2020;
originally announced January 2020.
-
Binospec: A Wide-field Imaging Spectrograph for the MMT
Authors:
Daniel Fabricant,
Robert Fata,
Harland Epps,
Thomas Gauron,
Mark Mueller,
Joseph Zajac,
Stephen Amato,
Jack Barberis,
Henry Bergner,
Patricia Brennan,
Warren Brown,
Igor Chilingarian,
John Geary,
Vladimir Kradinov,
Brian McLeod,
Matthew Smith,
Deborah Woods
Abstract:
Binospec is a high throughput, 370 to 1000 nm, imaging spectrograph that addresses two adjacent 8' by 15' fields of view. Binospec was commissioned in late 2017 at the f/5 focus of the 6.5m MMT and is now available to all MMT observers. Aperture masks cut from stainless steel with a laser cutter are used to define the entrance apertures that range from 15' long slits to hundreds of 2" slitlets. Sy…
▽ More
Binospec is a high throughput, 370 to 1000 nm, imaging spectrograph that addresses two adjacent 8' by 15' fields of view. Binospec was commissioned in late 2017 at the f/5 focus of the 6.5m MMT and is now available to all MMT observers. Aperture masks cut from stainless steel with a laser cutter are used to define the entrance apertures that range from 15' long slits to hundreds of 2" slitlets. System throughputs, including the MMT's mirrors and the f/5 wide-field corrector peak at ~30%. Three reflection gratings, duplicated for the two beams, provide resolutions ($λ$/$Δλ$) between 1300 and $>$5000 with a 1" wide slit. Two through-the-mask guiders are used for target acquisition, mask alignment, guiding, and precision offsets. A full-time Shack-Hartmann wave front sensor allows continuous adjustment of primary mirror support forces, telescope collimation and focus. Active flexure control maintains spectrograph alignment and focus under varying gravity and thermal conditions.
△ Less
Submitted 8 May, 2019;
originally announced May 2019.
-
Linkage and Intermediate C-Gorenstein Dimensions
Authors:
Joseph P. Brennan,
Alexander York
Abstract:
This paper brings together two theories in algebra that have had been extensively developed in recent years. First is the study of various homological dimensions and what information such invariants can give about a ring and its modules. A collection of intermediate C-Gorenstein dimensions is defined and this allows generalizations of results concerning C-Gorenstein dimension and certain Serre-lik…
▽ More
This paper brings together two theories in algebra that have had been extensively developed in recent years. First is the study of various homological dimensions and what information such invariants can give about a ring and its modules. A collection of intermediate C-Gorenstein dimensions is defined and this allows generalizations of results concerning C-Gorenstein dimension and certain Serre-like conditions. Second is the theory of linkage first introduced by Peskine and Szpiro and generalized to modules by Martinskovsky and Strooker. Using the further generalization of module linkage of Nagel, results are proven connecting linkage with these homological dimensions and Serre-like conditions.
△ Less
Submitted 8 October, 2018; v1 submitted 2 October, 2018;
originally announced October 2018.
-
An Extension of a Theorem of Frobenius and Stickelberger to Modules of Projective Dimension One over a Factorial Domain
Authors:
Joseph P. Brennan,
Alexander York
Abstract:
Let $R$ be a commutative ring. A quasi-Gorenstein $R$-module is an $R$-module such that the grade of the module and the projective dimension of the module are equal and the canonical module of the module is isomorphic to the module itself. After discussing properties of finitely generated quasi-Gorenstein modules, it is shown that this definition allows for a characterization of diagonal matrices…
▽ More
Let $R$ be a commutative ring. A quasi-Gorenstein $R$-module is an $R$-module such that the grade of the module and the projective dimension of the module are equal and the canonical module of the module is isomorphic to the module itself. After discussing properties of finitely generated quasi-Gorenstein modules, it is shown that this definition allows for a characterization of diagonal matrices of maximal rank over a Noetherian factorial domain $R$ extending a theorem of Frobenius and Stickelberger to modules of projective dimension 1 over a commutative Noetherian factorial domain.
△ Less
Submitted 26 June, 2018;
originally announced June 2018.
-
The effects of kinetic instabilities on the electron cyclotron emission from runaway electrons
Authors:
Chang Liu,
Lei Shi,
Eero Hirvijoki,
Dylan P. Brennan,
Amitava Bhattacharjee,
Carlos Paz-Soldan,
Max E. Austin
Abstract:
In this paper we show that the kinetic instabilities associated with runaway electron beams play an essential role for the production of high-level non-thermal electron-cyclotron-emission (ECE) radiation. Most of the non-thermal ECE comes from runaway electrons in the low-energy regime with large pitch angle, which are strongly scattered by the excited whistler waves. The power of ECE from runaway…
▽ More
In this paper we show that the kinetic instabilities associated with runaway electron beams play an essential role for the production of high-level non-thermal electron-cyclotron-emission (ECE) radiation. Most of the non-thermal ECE comes from runaway electrons in the low-energy regime with large pitch angle, which are strongly scattered by the excited whistler waves. The power of ECE from runaway electrons is obtained using a synthetic diagnostic model based on the reciprocity method. The electron distribution function is calculated using a kinetic simulation model including the whistler wave instabilities and the quasilinear diffusion effects. Simulations based on DIII-D low-density discharge reproduces the rapid growth of the ECE signals observed in DIII-D experiments. Unlike the thermal ECE where radiation for a certain frequency is strongly localized inside the resonance region, the non-thermal ECE radiation from runaway electrons is nonlocal, and the emission-absorption ratio is higher than that of thermal electrons. The runaway electron tail is more significant for ECE with higher frequencies, and the ECE spectrum becomes flatter as RE population grows. The nonlinear behavior of the kinetic instabilities is illustrated in the osculations of the ECE waves. The good agreement with the DIII-D experimental observations after including the kinetic instabilities clearly illustrate the significance of the scattering effects from wave-particle interactions, which can also be important for runaway electrons produced in disruptions.
△ Less
Submitted 27 March, 2018;
originally announced March 2018.
-
Role of kinetic instability in runaway electron avalanche and elevated critical electric fields
Authors:
Chang Liu,
Eero Hirvijoki,
Guo-yong Fu,
Dylan P. Brennan,
Amitava Bhattacharjee,
Carlos Paz-Soldan
Abstract:
The effects of kinetic whistler-wave instabilities on the runaway-electron (RE) avalanche is investigated. With parameters from DIII-D experiments, we show that RE scattering from excited whistler waves can explain several poorly understood experimental results seen in a variety of tokamaks. We find an increase of the avalanche growth rate and threshold electric field, bringing the present model m…
▽ More
The effects of kinetic whistler-wave instabilities on the runaway-electron (RE) avalanche is investigated. With parameters from DIII-D experiments, we show that RE scattering from excited whistler waves can explain several poorly understood experimental results seen in a variety of tokamaks. We find an increase of the avalanche growth rate and threshold electric field, bringing the present model much closer to observations than previous results. The excitation of kinetic instabilities and the scattering of resonant electrons are calculated self-consistently using a quasilinear model. We also explain the observed fast growth of electron cyclotron emission (ECE) signals and excitation of very low-frequency whistler modes observed in the quiescent RE experiments at DIII-D [D. A. Spong et al., submitted to Phys. Rev. Lett.]. These results indicate that by controlling the background thermal plasma temperature, the plasma wave can be excited spontaneously in tokamak disruptions and the avalanche generation of runaway electrons may be suppressed.
△ Less
Submitted 5 January, 2018;
originally announced January 2018.
-
Real frequency tearing layers with parallel dynamics and the effect on locking and resistive wall modes
Authors:
John M. Finn,
Andrew J. Cole,
Dylan P. Brennan
Abstract:
Tearing modes with real frequencies in the plasma frame (i.e. in addition to the Doppler shift due to $E\times B$ rotation) are of potential importance because of their effect on the locking process. In particular, it has recently been shown [J. M. Finn, A. J. Cole, and D. P. Brennan, Phys. Plasmas (Letters), 22:120701, 2015] that the Maxwell torque on the plasma in the presence of an applied erro…
▽ More
Tearing modes with real frequencies in the plasma frame (i.e. in addition to the Doppler shift due to $E\times B$ rotation) are of potential importance because of their effect on the locking process. In particular, it has recently been shown [J. M. Finn, A. J. Cole, and D. P. Brennan, Phys. Plasmas (Letters), 22:120701, 2015] that the Maxwell torque on the plasma in the presence of an applied error field is modified significantly for tearing modes having real frequencies near marginal stability. In addition, it is known [J. M. Finn and R. A. Gerwin, Phys. Plasmas, 3:2344, 1996] that resistive wall tearing modes can be destabilized below their no-wall limits by rotation, if the tearing modes have real frequencies near marginal stability. In this paper we first derive the tearing mode dispersion relation with pressure gradient, field line curvature and parallel dynamics in the resistive-inertial (RI) regime, neglecting the divergence of the $E\times B$ drift and perpendicular resistivity. The results show that the usual Glasser effect, a toroidal effect which involves real frequencies, occurs in this simplified model, which ignores perpendicular resistivity and the divergence of the $E\times B$ drift. We also find, using a similar simple model, the surprising result that in the viscoresistive regime with pressure gradient, favorable curvature due to toroidal effects, and parallel dynamics, a similar Glasser-like effect is found. We show that in both regimes the existence of tearing modes with complex frequencies is related to nearby electrostatic resistive interchange modes with complex frequencies. We discuss the effect on locking to an error field and the significant lowering of the threshold for destabilization of resistive wall tearing modes, which can be much more pronounced than the weak effect observed for RI tearing modes without pressure-curvature drive in Finn and Gerwin, 1996.
△ Less
Submitted 15 August, 2017;
originally announced August 2017.
-
A Model of Energetic Ion Effects on Pressure Driven Tearing Modes in Tokamaks
Authors:
Michael R. Halfmoon,
Dylan P. Brennan
Abstract:
The effects that energetic trapped ions have on linear resistive magnetohydrodynamic (MHD) instabilities are studied in a reduced model that captures the essential physics driving or damping the modes through variations in the magnetic shear. The drift-kinetic orbital interaction of a slowing down distribution of trapped energetic ions with a resistive MHD instability is integrated to a scalar con…
▽ More
The effects that energetic trapped ions have on linear resistive magnetohydrodynamic (MHD) instabilities are studied in a reduced model that captures the essential physics driving or damping the modes through variations in the magnetic shear. The drift-kinetic orbital interaction of a slowing down distribution of trapped energetic ions with a resistive MHD instability is integrated to a scalar contribution to the perturbed pressure, and entered into an asymptotic matching formalism for the resistive MHD dispersion relation. Toroidal magnetic field line curvature is included to model trapping in the particle distribution, in an otherwise cylindrical model. The focus is on a configuration that is driven unstable to the m/n = 2/1 mode by increasing pressure, where m is the poloidal mode number and n the toroidal. The particles and pressure can affect the mode both in the core region where there can be low and reversed shear and outside the resonant surface in significant positive shear. The results show that the energetic ions damp and stabilize the mode when orbiting in significant positive shear, increasing the marginal stability boundary. However, the inner core region contribution with low and reversed shear can drive the mode unstable. This effect of shear on the energetic ion pressure contribution is found to be consistent with the literature. These results explain the observation that the 2/1 mode was found to be damped and stabilized by energetic ions in δf - MHD simulations of tokamak experiments with positive shear throughout, while the 2/1 mode was found to be driven unstable in simulations of experiments with weakly reversed shear in the core. This is also found to be consistent with related experimental observations of the stability of the 2/1 mode changing significantly with core shear.
△ Less
Submitted 2 March, 2017; v1 submitted 22 February, 2017;
originally announced February 2017.
-
Adjoint Fokker-Planck equation and runaway electron dynamics
Authors:
Chang Liu,
Dylan P. Brennan,
Amitava Bhattacharjee,
Allen H. Boozer
Abstract:
The adjoint Fokker-Planck equation method is applied to study the runaway probability function and the expected slowing-down time for highly relativistic runaway electrons, including the loss of energy due to synchrotron radiation. In direct correspondence to Monte Carlo simulation methods, the runaway probability function has a smooth transition across the runaway separatrix, which can be attribu…
▽ More
The adjoint Fokker-Planck equation method is applied to study the runaway probability function and the expected slowing-down time for highly relativistic runaway electrons, including the loss of energy due to synchrotron radiation. In direct correspondence to Monte Carlo simulation methods, the runaway probability function has a smooth transition across the runaway separatrix, which can be attributed to effect of the pitch angle scattering term in the kinetic equation. However, for the same numerical accuracy, the adjoint method is more efficient than the Monte Carlo method. The expected slowing-down time gives a novel method to estimate the runaway current decay time in experiments. A new result from this work is that the decay rate of high energyelectrons is very slow when E is close to the critical electric field. This effect contributes further to a hysteresis previously found in the runaway electron population.
△ Less
Submitted 21 January, 2016; v1 submitted 15 September, 2015;
originally announced September 2015.
-
Error field penetration and locking to the backward propagating wave
Authors:
John M. Finn,
Andrew J. Cole,
Dylan P. Brennan
Abstract:
Resonant field amplification or error field penetration involves driving a weakly stable tearing perturbation in a rotating toroidal plasma. In this paper it is shown that the locking characteristics for modes with finite real frequencies $ω_{r}$ are quite different from the conventional results. A calculation of the tearing mode amplitude assuming modes with frequencies $\pmω_{r}$ in the plasma f…
▽ More
Resonant field amplification or error field penetration involves driving a weakly stable tearing perturbation in a rotating toroidal plasma. In this paper it is shown that the locking characteristics for modes with finite real frequencies $ω_{r}$ are quite different from the conventional results. A calculation of the tearing mode amplitude assuming modes with frequencies $\pmω_{r}$ in the plasma frame shows that it is maximized when the frequency of the stable backward propagating mode ($-ω_{r}$) in the lab frame is zero, i.e. when $v=ω_{r}/k$. Even more importantly, the locking torque is exactly zero at the mode phase velocity, with a pronounced peak at just higher rotation, leading to a locked state with plasma velocity $v$ just above the mode phase velocity in the lab frame. Real frequencies $\pmω_{r}$, leading to a $v\rightarrow-v$ symmetry, are known to occur due to the Glasser effect [A.H. Glasser, J.M. Greene, and J.M. Johnson, Phys. Fluids {\bf 19}, 567 (1976).] for modes in the resistive-inertial (RI) regime. This therefore leads to locking of the plasma velocity to just above the phase velocity. It is also shown that similar real frequencies occur over a range of parameters in the visco-resistive (VR) regime with pressure, and the locking torque is similar to that in the RI regime. Real frequencies occur due to diamagnetic effects in other tearing mode regimes and also show this effect, but without the $v\rightarrow-v$ symmetry. Nonlinear effects on the mode amplitude and torque for weakly stable modes or large error fields are discussed. Also, the possibility of applying external fields of different helicities to drive sheared flows in toroidal plasmas is discussed.
△ Less
Submitted 10 August, 2015; v1 submitted 14 July, 2015;
originally announced July 2015.
-
Dynamics of Ion Temperature Gradient Turbulence and Transport with a Static Magnetic Island
Authors:
Olivier Izacard,
Christopher Holland,
Spencer D. James,
Dylan P. Brennan
Abstract:
Understanding the interaction mechanisms between large-scale magnetohydrodynamic instabilities and small-scale drift-wave microturbulence is essential for predicting and optimizing the performance of magnetic confinement based fusion energy experiments. We report progress on understanding these interactions using both analytic theory and numerical simulations performed with the BOUT++ [B. Dudson e…
▽ More
Understanding the interaction mechanisms between large-scale magnetohydrodynamic instabilities and small-scale drift-wave microturbulence is essential for predicting and optimizing the performance of magnetic confinement based fusion energy experiments. We report progress on understanding these interactions using both analytic theory and numerical simulations performed with the BOUT++ [B. Dudson et al., Comput. Phys. Comm. 180, 1467 (2009)] framework. This work focuses upon the dynamics of the ion temperature gradient instability in the presence of a background static magnetic island, using a weakly electromagnetic two-dimensional five-field fluid model. It is found that the island width must exceed a threshold size (comparable to the turbulent correlation length in the no-island limit) to significantly impact the turbulence dynamics, with the primary impact being an increase in turbulent fluctuation and heat flux amplitudes. The turbulent radial ion energy flux is shown to localize near the X-point, but does so asymmetrically in the poloidal dimension. An effective turbulent resistivity which acts upon the island outer layer is also calculated, and shown to always be significantly (10x - 100x) greater than the collisional resistivity used in the simulations.
△ Less
Submitted 12 September, 2015; v1 submitted 24 June, 2015;
originally announced June 2015.
-
Control of resistive wall modes in a cylindrical tokamak with plasma rotation and complex gain
Authors:
D. P. Brennan,
J. M. Finn
Abstract:
Feedback stabilization of magnetohydrodynamic (MHD) modes is studied in a cylindrical model for a tokamak with resistivity, viscosity and toroidal rotation. The control is based on a linear combination of the normal and tangential components of the magnetic field just inside the resistive wall. The feedback includes complex gain, for both the normal and for the tangential components, and the imagi…
▽ More
Feedback stabilization of magnetohydrodynamic (MHD) modes is studied in a cylindrical model for a tokamak with resistivity, viscosity and toroidal rotation. The control is based on a linear combination of the normal and tangential components of the magnetic field just inside the resistive wall. The feedback includes complex gain, for both the normal and for the tangential components, and the imaginary part of the feedback for the former is equivalent to plasma rotation. The work includes (1) analysis with a reduced resistive MHD model for a tokamak with finite βand with stepfunction current density and pressure profiles, and (2) computations with full compressible visco-resistive MHD and smooth decreasing profiles of current density and pressure. The equilibria are stable for β=0 and the marginal stability values $β_{rp,rw}<β_{rp,iw}<β_{ip,rw}<β_{ip,iw}$ (resistive plasma, resistive wall; resistive plasma, ideal wall; ideal plasma, resistive wall; ideal plasma, ideal wall) are computed for both cases. The main results are: (a) imaginary gain with normal sensors or plasma rotation stabilizes below $β_{rp,iw}$ because rotation supresses the diffusion of flux from the plasma out through the wall and, more surprisingly, (b) rotation or imaginary gain with normal sensors destabilizes above $β_{rp,iw}$ because it prevents the feedback flux from entering the plasma through the resistive wall to form a virtual wall. The effect of imaginary gain with tangential sensors is more complicated but essentially destabilizes above and below $β_{rp,iw}$. A method of using complex gain to optimize in the presence of rotation in the $β>β_{rp,iw}$ regime is presented.
△ Less
Submitted 19 June, 2014;
originally announced June 2014.
-
Noether normalizations, reductions of ideals, and matroids
Authors:
Joseph P. Brennan,
Neil Epstein
Abstract:
We show that given a finitely generated standard-graded algebra of dimension $d$ over an infinite field, its graded Noether normalizations obey a certain kind of `generic exchange', allowing one to pass between any two of them in at most $d$ steps. We prove analogous generic exchange theorems for minimal reductions of an ideal, minimal complete reductions of a set of ideals, and minimal complete r…
▽ More
We show that given a finitely generated standard-graded algebra of dimension $d$ over an infinite field, its graded Noether normalizations obey a certain kind of `generic exchange', allowing one to pass between any two of them in at most $d$ steps. We prove analogous generic exchange theorems for minimal reductions of an ideal, minimal complete reductions of a set of ideals, and minimal complete reductions of multigraded $k$-algebras. Finally, we unify all these results into a common axiomatic framework by introducing a new topological-combinatorial structure we call a generic matroid, which is a common generalization of a topological space and a matroid.
△ Less
Submitted 1 August, 2010;
originally announced August 2010.
-
Efficient quantum processing of ideals in finite rings
Authors:
Pawel M. Wocjan,
Stephen P. Jordan,
Hamed Ahmadi,
Joseph P. Brennan
Abstract:
Suppose we are given black-box access to a finite ring R, and a list of generators for an ideal I in R. We show how to find an additive basis representation for I in poly(log |R|) time. This generalizes a quantum algorithm of Arvind et al. which finds a basis representation for R itself. We then show that our algorithm is a useful primitive allowing quantum computers to rapidly solve a wide variet…
▽ More
Suppose we are given black-box access to a finite ring R, and a list of generators for an ideal I in R. We show how to find an additive basis representation for I in poly(log |R|) time. This generalizes a quantum algorithm of Arvind et al. which finds a basis representation for R itself. We then show that our algorithm is a useful primitive allowing quantum computers to rapidly solve a wide variety of problems regarding finite rings. In particular we show how to test whether two ideals are identical, find their intersection, find their quotient, prove whether a given ring element belongs to a given ideal, prove whether a given element is a unit, and if so find its inverse, find the additive and multiplicative identities, compute the order of an ideal, solve linear equations over rings, decide whether an ideal is maximal, find annihilators, and test the injectivity and surjectivity of ring homomorphisms. These problems appear to be hard classically.
△ Less
Submitted 5 July, 2023; v1 submitted 31 July, 2009;
originally announced August 2009.
-
Duality, a-invariants and canonical modules of rings arising from linear optimization problems
Authors:
Joseph P. Brennan,
Luis A. Dupont,
Rafael H. Villarreal
Abstract:
The aim of this paper is to study integer rounding properties of various systems of linear inequalities to gain insight about the algebraic properties of Rees algebras of monomial ideals and monomial subrings. We study the normality and Gorenstein property--as well as the canonical module and the a-invariant--of Rees algebras and subrings arising from systems with the integer rounding property.…
▽ More
The aim of this paper is to study integer rounding properties of various systems of linear inequalities to gain insight about the algebraic properties of Rees algebras of monomial ideals and monomial subrings. We study the normality and Gorenstein property--as well as the canonical module and the a-invariant--of Rees algebras and subrings arising from systems with the integer rounding property. We relate the algebraic properties of Rees algebras and monomial subrings with integer rounding properties and present a duality theorem.
△ Less
Submitted 3 December, 2008;
originally announced December 2008.
-
ELM triggering conditions for the integrated modeling of H-mode plasmas
Authors:
A. Y. Pankin,
G. Bateman,
D. P. Brennan,
D. D. Schnack,
P. B. Snyder,
I. Voitsekhovitch,
A. H. Kritz,
G. Janeschitz,
S. Kruger,
T. Onjun,
G. W. Pacher,
H. D. Pacher
Abstract:
Recent advances in the integrated modeling of ELMy H-mode plasmas are presented. A model for the H-mode pedestal and for the triggering of ELMs predicts the height, width, and shape of the H-mode pedestal and the frequency and width of ELMs. Formation of the pedestal and the L-H transition is the direct result of ExB flow shear suppression of anomalous transport. The periodic ELM crashes are tri…
▽ More
Recent advances in the integrated modeling of ELMy H-mode plasmas are presented. A model for the H-mode pedestal and for the triggering of ELMs predicts the height, width, and shape of the H-mode pedestal and the frequency and width of ELMs. Formation of the pedestal and the L-H transition is the direct result of ExB flow shear suppression of anomalous transport. The periodic ELM crashes are triggered by either the ballooning or peeling MHD instabilities. The BALOO, DCON, and ELITE ideal MHD stability codes are used to derive a new parametric expression for the peeling-ballooning threshold. The new dependence for the peeling-ballooning threshold is implemented in the ASTRA transport code. Results of integrated modeling of DIII-D like discharges are presented and compared with experimental observations. The results from the ideal MHD stability codes are compared with results from the resistive MHD stability code NIMROD.
△ Less
Submitted 8 November, 2004;
originally announced November 2004.