Neoclassical tearing mode seeding by Alfvén eigenmode coupling in MAST-U

Nuclear Fusion 66:12 (2026)

Authors:

K Thomas, L Bardoczi, K Gibson, J Ruiz Ruiz, M Dreval, D Ryan, K McClements, C Michael

Abstract:

We present evidence of tearing mode seeding by nonlinear three-wave coupling of Alfvén eigenmodes in a tokamak plasma. Here, a m,n = 2,1 magnetic island appears to be seeded in a MAST-U discharge by coupling of (Formula presented) (Formula presented) and (Formula presented) (Formula presented) Alfvén eigenmodes, which is demonstrated by increased squared normalised bicoherence calculated from magnetic probes preceding tearing mode onset. This observation offers a new explanation for tearing modes that would be otherwise addressed as spontaneous, and indicates three-wave coupling can seed 2,1 magnetic islands in plasmas without pre-existing (Formula presented) (Formula presented) tearing modes. Further, this mechanism has the potential to be more prevalent in burning plasmas expected in future fusion devices.

Overview of the MAST Upgrade physics programme: testing novel concepts at low aspect ratio to inform future devices

Nuclear Fusion IOP Publishing 66:11 (2026)

Authors:

JR Harrison, A Aboutaleb, M Aljunid, SY Allan, R Allan, A Alli, H Anand, Y Andrew, LC Appel, A Ash, J Ashton, F Auriemma, O Bachmann, S Bakes, P Balazs, O Bardsley, M Barnes, B Barrett, D Baver, C Beckley, J Bennett, J Bentley, J Berkery, M Bernert, N Bertelli, P Bilkova, S Blackmore, A Bock, W Boeglin, P Bohm, J Booth, A Bosnjak, C Bowman, J Bradley, L Bramucci, D Brida, M Brix, PK Browning, D Brunetti, P Bryant, J Bryant, J Buchanan, M Bull, M Bulman, N Bulmer, D Burke, P Cahill, P Cano-Megías, M Carpita, A Carruthers, L Casali, A Cassidy, F Causa, M Cecconello, A Celora, M Cha, B Chamberlain, B Chapman, B Chapman-Oplopoiou, R Chazal, X Chen, J Clark, M Clark, T Clayton, K Collie, D Collishaw-Schepman, JW Connor, M Coy, S Cramp, N Crocker, D Cruse, G Cunningham, M Czarski, I Cziegler, Y Damizia, P Davies, K Davies, I Day, E Delabie, GL Derks, S Dixon, J Dobrashian, M Dreval, X Du, H Dudding, D Dunai, M Dunne, F Ebrahimi, J Edmond, J Edwards, H El-Haroun, S Elmore, Y Enters, M Faitsch, F Federici, N Fedorczak, F Felici, N Ferraro, AR Field, P Figueiredo, I FitzGerald, M Fitzgerald, R Fitzpatrick, C Fitzpatrick, S Frankel, L Frassinetti, D Frattolillo, S Freethy, W Fuller, R Futtersack, S Gabriellini, K Gage, J Galdon, J Galdon-Quiroga, M Gambrioli, C Garner, L Garzotti, TE Gheorghiu, KJ Gibson, C Gibson, E Giovannozzi, C Giroud, J Goatley, A Goodyear, M Gorelenkova, S Gosden, G Grassler, JP Graves, D Greenhouse, R Griffiths, VH Hall-Chen, CJ Ham, E Harrington, R Harrison, A Haupt, J Hawes, S Hegedus, SS Henderson, C Heo, C Hickling, M Hill, B Hnat, C Hogben, B Honey, L Howlett, Z Huang, A Hudoba, J Hughes, R Hussain, K Imada, P Ivanov, A Jackson, P Jacquet, F Jaulmes, P Jepson, T Jones, P Jones, M Juvonen, V Kachkanov, B Kandan, I Katramados, S Kaye, YO Kazakov, D Keeling, D Kennedy, A Kenny, H-T Kim, D King, R King, A King, V Kiptily, A Kirjasuo, A Kirk, A Kit, A Kleiner, M Knolker, S Kobussen, M Kochan, L Kogan, B Kool, D Kos, M Kotschenreuther, M Lampert, A Lawson, K Lawson, K-W Lee, G Lee, J Lee, M Lees, S Leigh, AW Leonard, G Liddiard, B Lipschultz, E Litherland-Smith, YQ Liu, BA Lomanowski, J Lombardo, N Lonigro, J Lore, J Lovell, R Lucock, T Luong, A Lvovskiy, J Macdonald, T Macwan, S Mahajan, F Maiden, R Maingi, C Man-Friel, F Mansfield, M Markl, S Marsden, R Martin, R Mathew, R Maurizio, U Mazzarese, S Mazzi, R McAdams, G McArdle, J McBride, K McClements, J McClenaghan, D McConville, K McKay, C McKnight, P McKnight, A McLean, BF McMillan, A McShee, J Measures, N Mehay, S Menmuir, HF Meyer, CA Michael, F Militello, IG Miron, R Mishra, J Mitchell, D Moiraf, P Monaghan, R Mooney, N Mooring, R Morales Gomes, D Morbey, S Mordijck, C Morgan, J Morris, D Moulton, S Munaretto, A Munasinghe, A Muraro, O Myatra, Y-S Na, TF Neiser, AO Nelson, SL Newton, M Nicassio, MG O’Mullane, C Olde, HJ Oliver, P Ollus, J Omotani, M Ono, FP Orsitto, R Osawa, N Osborne, T Osborne, R Otin, E Ozturk, F Palermo, A Pankin, I Paradela Pérez, J Parisi, E Parr, B Parry, BS Patel, E Pawelec, D Payne, C Paz-Soldan, A Phelps, L Piron, C Piron, G Pokol, R Preece, M Price, B Pritchard, R Proudfoot, G Pucella, T Pumfrett, D-Y Pyo, H Reimerdes, T Rhodes, E Ribeiro, D Rigamonti, J Riquezes, JF Rivero-Rodriguez, J Roberts, M Robson, K Ronald, E Rose, D Ryan, P Ryan, S Saarelma, S Sabbagh, A Salmi, R Sarwar, P Saunders, O Sauter, R Scannell, R Sealey, R Seath, S Sharapov, R Sharma, H Sheikh, S Shiraiwa, B Sieglin, SA Silburn, M Simmonds, J Simpson, A Sladkomedova, J Smith, P Smith, M Sos, VA Soukhanovskii, D Speirs, C Srinivasan, G Staebler, R Stephen, P Stevenson, J Stobbs, C Stroud, H Sun, H Sun, G Szepesi, DM Takács, T Tala, C Tame, C Theiler, B Thomas, S Thomas, S Thomas, N Thomas-Davies, AJ Thornton, A Tilley, I Tirkova, M Tobin, E Tomasina, A Tonel, P Tonner, A Tookey, G Tvalashvili, M Vallar, M Valovic, RGL Vann, L Velarde, L Velarde, K Verhaegh, E Viezzer, C Vincent, M Walsh, T Walsh, M Warr, S Wiesen, TA Wijkamp, D Wilkins, J Willis, T Wilson, HR Wilson, N Winston, G Withenshaw, H Wong, M Wood, R Worrall, Q Xia, G Xia, L Xiang, G Xiang, T Xu, JH Yu, V Zamkovska, M Zerbini, VK Zotta, M Zurita, LE di Grazia

Abstract:

The research programme performed on the Mega Amp Spherical Tokamak (MAST) Upgrade experiment has made significant advances in developing the physics understanding of low aspect ratio tokamaks in support of the operation of ITER and design of fusion powerplants. High performance plasma scenarios have been developed to facilitate a broad programme of experiments, in which confinement is constrained by the presence of m/n = 2/1 modes that cause substantial losses of fast ions. The onset of these modes coincides with the q = 2 surface residing in a local minimum in the toroidal current density profile. The maximum electron temperature at the pedestal top, Te,ped is limited with gas fuelling to ∼350 eV to maintain regular ELMs; higher Te,ped results in a transition to a non-stationary ELM-free regime. The operational space of spherical tokamaks has been expanded into small and ELM-free regimes. Strong shaping of the last closed flux surface can induce a transition from large to small ELMs, and ELM suppression with resonant magnetic perturbations has been observed for the first time in a low aspect ratio tokamak. Negative triangularity shaping has induced a transition from ELMy H-mode to a high-performance L-mode regime for the first time in a low aspect ratio tokamak. In studies of fast ion confinement, losses of fast particles due to Global Alfvén Eigenmodes have been identified. Interactions between fast ions generated by off-axis neutral beam injection and thermal neutrals can result in significant losses of fast ions. Experiments with on- and off-axis neutral beam injection exhibit a flux pumping mechanism, where the central safety factor is held to ∼1 in the absence of sawteeth. In studies of pedestal physics, it has been found that elevated main chamber neutral pressures result in an increase in the electron density and reduction in the temperature at the pedestal top. Advances in understanding plasma exhaust include the integration of a high-performance plasma core with detached outer divertors in the X-point target configuration. A newly commissioned lower divertor cryopump reduces the lower divertor neutral pressure by up to 50%, with minimal effect on the main chamber or upper divertor. New measurements and SOLPS-ITER simulations emphasise the importance of plasma–neutral interactions on divertor detachment in the conditions accessible in experiments. Real-time control of the ionisation front location in both divertor chambers independently has been demonstrated in double null experiments, enabled by the tightly baffled divertor chambers.

Evaluation of beam properties using a checkerboard masking configuration at the NNBI test stand BATMAN Upgrade

Fusion Engineering and Design Elsevier 232 (2026) 116017

Authors:

Araceli Navarro, Eleni Nanou, Michael Barnes, Christian Wimmer, Dirk Wünderlich, Ursel Fantz

Abstract:

Negative ion based neutral beam injection (NBI) systems will be used for heating and current drive at the ITER tokamak and require stringent beam quality. Including a beamlet core divergence of less than 7 mrad and high uniformity across individual beamlets (better than 90%), at extracted current densities of 329 A/m 2 for hydrogen and 286 A/m 2 for deuterium. The BATMAN Upgrade (BUG) test facility, a 1/8 scale ITER source, is used to investigate and optimize these beam properties. Beamlet divergence and uniformity at BUG test facility are measured using a retractable Carbon Fiber Composite (CFC) calorimeter. However, under standard operation (with all 70 apertures open), strong beamlet overlap prevents the extraction of single-beamlet divergences from the CFC data. To overcome this limitation, a checkerboard masking configuration was implemented. This masking, combined with a multi-Gaussian model allows to obtain beamlet properties, including size, shape, and intensity along the whole extraction area. The experimental methodology and data analysis techniques are presented in detail. Their capabilities are showcased by presenting measurements performed at BUG under different configurations of the magnetic filter field (drift-up and drift-down conditions) and for both hydrogen and deuterium operation.

Double-adiabatic equations of state for relativistic plasmas

Journal of Plasma Physics Cambridge University Press (CUP) 92:5 (2026) e132

Authors:

A Wierzchucka, PJ Bilbao, AGR Thomas, DA Uzdensky, AA Schekochihin

Abstract:

The adiabatic equation of state upper P proportional to n Superscript upper Gamma P ∝ n Γ $P \propto n^{\varGamma }$ describes the pressure evolution of highly collisional, isotropic plasmas in terms of their density, providing a possible closure of the fluid moment hierarchy in the absence of heat fluxes and dissipation. An analogous closure exists for collisionless, magnetised plasmas, whose pressure tensor is anisotropic with respect to the magnetic field, and the closure is therefore double adiabatic, prescribing the evolution of the parallel and perpendicular pressures in terms of the magnetic-field strength and density. Here, we present a general first-principles formalism to derive adiabatic laws using symmetries and phase-space conservation properties of the distribution function. With this theory we recover the adiabatic equation of state upper P proportional to n Superscript upper Gamma P ∝ n Γ $P \propto n^{\varGamma }$ for isotropic plasmas and the double-adiabatic equations of state for collisionless, magnetised plasmas. We extend the latter to the relativistic regime, finding that their exact functional form depends on the pressure anisotropy and is not a simple power law. Our double-adiabatic equations of state describe simple geometries, like magnetic mirrors or compressed homogeneous plasmas, as well as complex high-energy astrophysical processes, such as the evolution of plasmoid structures formed during magnetic reconnection.

High-beta runaway transitions in a fluid model of electromagnetic ion-temperature-gradient turbulence

Plasma Physics and Controlled Fusion IOP Publishing (2026)

Authors:

Yujia Zhang, Michael Barnes, Alex Schekochihin, Plamen G Ivanov, Toby Adkins

Abstract:

Abstract Gyrokinetic simulations of tokamak turbulence indicate that fluctuation levels increase abruptly and dramatically when the plasma beta exceeds a certain critical value. This increase in fluctuation levels coincides with a transition from a state dominated by zonal flow to one in which turbulent eddies form radially-elongated 'streamers'. Here we derive from gyrokinetics a minimal fluid model for electromagnetic ion-temperature-gradient (ITG) turbulence that captures the key features of this transition. Due to the relative simplicity of the model, we are able to conduct a detailed numerical study of the interplay between the turbulence and the zonal flow across a broad range of values of the plasma beta and the ITG. We find that the transition occurs when the Reynolds stress, which tends to strengthen zonal flows, is overwhelmed by the Maxwell and diamagnetic stresses, which tend to weaken them. Power-law scalings of the stress ratios with plasma beta and ITG are obtained, indicating a possible means by which the location of the transition could be predicted with minimal computational cost.