## First Access to ELM-free Negative Triangularity at Low Aspect Ratio

A.O. Nelson, 1, ∗ C. Vincent, 2 H. Anand, 3 J. Lovell, 4 J.F. Parisi, 5 H. S. Wilson, 1 K. Imada, 6 W.P. Wehner, 3 M. Kochan, 2 S. Blackmore, 2 G. McArdle, 2 S. Guizzo, 1 L. Rondini, 1 S. Freiberger, 1 C. Paz-Soldan, 1 and the MAST-U Team 1 Columbia University, New York City, New York, USA 2 UKAEA-CCFE, Culham Science Centre, Abingdon, Oxon, UK 3 General Atomics, San Diego, California, USA 4 Oak Ridge National Laboratory, Oak Ridge, Tennessee, USA 5 Princeton Plasma Physics Laboratory, Princeton, New Jersey, USA

6 York Plasma Institute, Heslington, York, UK

(Dated: August 2, 2024)

A plasma scenario with negative triangularity (NT) shaping is achieved on MAST-U for the first time. While edge localized modes (ELMs) are eventually suppressed as the triangularity is decreased below δ ≲ -0 06, an extended period of H-mode operation with Type-III ELMs is sustained at less . negative δ even through access to the second stability region for ideal ballooning modes is closed. This documents a qualitative difference from the ELM-free access conditions documented in NT scenarios on conventional aspect ratio machines. The electron temperature at the pedestal top drops across the transition to ELM-free operation, but a steady rise in core temperature as δ is decreased allows for similar normalized β in the ELM-free NT and H-mode positive triangularity shapes.

## INTRODUCTION

Tokamak scenarios exploring negative triangularly (NT) shaping, where the plasma x-points are situated at a larger radius than the magnetic axis, have recently gathered significant interest from the international community due to their favorable confinement and powerhandling properties [1]. In machines with conventional aspect ratios A = R /a 0 minor ∼ 2 5 . -4, where R 0 and a minor are respectively the major and minor radii, NT experiments are often characterized by reduced electron heat flux throughout their volumes compared to similar discharges with positive triangularity (PT), which comes from a reduction in core turbulent transport [2, 3]. This leads to an enhancement of the core pressure gradient in NT plasmas that improves confinement over matched PT discharges [4, 5]. Notably, these same plasmas do not induce edge localized modes (ELMs), even at high heating powers [6]. This ELM-free behavior allows NT fusion power plant designs to avoid issues associated with plasma-wall interactions during ELM events [7-9]. The larger major radius of the x-point in NT scenarios also provides additional room for divertor heat flux spreading, highlighting the potential power-handling benefits of NT plasmas [10, 11].

While NT scenarios show tangible promise as an power plant candidate in conventional aspect ratio designs, they have not yet been fully explored in spherical tokamaks (STs) with lower aspect ratio ( A ≲ 2). Initial gyrokinetic simulations of NT scenarios at low A , both with gs2 [12] and GENE [13], predict that the confinement enhancement over PT discharges observed at conventional aspect ratio might be lost in NT ST scenarios, particularly in regimes dominated by trapped electron modes. In addition, NT

plasmas are predicted to be strongly unstable to longwavelength kinetic ballooning modes (KBMs) across the plasma in highβ , reactor-relevant STs [12]. This results from a closure of access to the so called 'second stability' region for ballooning modes, potentially limiting the maximum achievable β for NT scenarios in STs [7].

Notably, the loss of second stability access has also been strongly correlated with the prevention of ELMs in conventional aspect ratio NT discharges [6, 7, 14]. Access to second stability has long been associated with large H-mode pedestals [8, 15, 16], so it is initially unsurprising that closing the stability window for ideal ballooning modes would prevent H-mode access. On DIII-D, second stability access can be closed by directly scanning the triangularity δ below some critical value ( δ crit , DIII -D ∼ -0 15) . while holding other plasma and machine parameters (such as the injected power, plasma current, toroidal magnetic field and density) constant [6]. This shape change modifies the local shear near the plasma edge, closing access to the second stability regime and preventing the formation of an ELMy H-mode regime by severely limiting the allowable edge pressure gradient [14]. Negative triangularity plasmas on TCV follow a similar pattern, where closure of the second stability region is encountered in ELM-free NT plasmas [5, 7, 17].

As noted in the literature, however, closure of the second stability region is not itself always sufficient for Hmode prevention in PT plasmas, particularly at lower aspect ratios [8, 18, 19]. Further, experimentally measured gradients from DIII-D ELM-free plasmas with δ ∼ -0 5 . are not always indicative of being directly limited by the ideal ballooning mode [6]. These two facts suggest (1) that closure of the second ballooning stability region alone may not be sufficient to prevent ELMs, especially

FIG. 1: MAST-U discharge #49463 at 800 ms is diverted, nearly up-down symmetric, and has δ &lt; 0.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00180v1-figure-1d630100-c641-4566-99a6-9c43422e631c.jpg)

in regimes STs where first-stable H-modes are more common, and (2) that additional transport mechanisms that further limit the pressure gradient below the ballooning limit may be present even in conventional aspect ratio machines, especially when δ ≪ 0.

In this Letter, we present initial observations and modeling from the first ELM-free NT discharge achieved on MAST-U. Importantly, this work demonstrates a distinct decorrelation between the closure of the second stability region and the prevention of ELMy H-mode operation that has not yet been observed in conventional aspect ratio machines. Section details the access conditions for NT operation on MAST-U, leading to documentation of the stability modeling in section . In section , the potential consequences of this result are briefly discussed, particularly with respect to the construction of a predictive model for the NT pedestal height, which may ultimately set the performance of NT reactor scenarios.

## NT OPERATION ON THE MAST-U SPHERICAL TOKAMAK

The magnetic equilibrium for MAST-U discharge #49463, which features diverted negative triangularity

FIG. 2: (a) The plasma current I p , (b) injected beam power P inj (black) and β N (red), (c) triangularity δ , (d) line-integrated average density n e (black) and carbon intensity I C (red) from SPRED, (e) the plasma elongation κ (black) and normalized confinement H 98 y2 , (red), and (f) D α emission lines for MAST-U discharge #49463. The ELM-free period is marked in blue.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00180v1-figure-20234b31-e76c-4b64-8d04-abf5e668d157.jpg)

operation, is reconstructed by the TokaMaker code [20] in figure 1. This discharge was achieved during a snowflake divertor development experiment in which the GA-TokSys code [21] was employed to reduce the plasma squareness at high elongation [22, 23]. During this particular discharge, the lower x-point (subject to vertical position control) left the x-point control region as the plasma elongation was reduced, at which point the plasma control system continued to push the plasma in feedback. Luckily, this resulted in a smooth triangularity scan at constant power and near constant density, providing an excellent first NT discharge for the initial physics investigations presented in this work. Further details on this control scheme, including techniques that can be employed to create stable NT discharges on MAST-U, will be the subject of a subsequent publication.

Figure 2 shows major time traces for this discharge, highlighting the stable nature of the plasma as the shape is changed. The discharge has constant plasma current ( I p = 600kA) and beam heating power ( P inj = 3MW).

A shape scan starts around t = 200ms, as can be seen in the δ scan presented in figure 2c. We note that the strong dip in δ near t = 420ms is a result of the inner strikeline moving on to the outboard portion of the vacuum vessel. The line-integrated average density (figure 2d) grows until t = 400ms, at which point an ELMy H-mode appears. The density then remains roughly constant for the remainder of the discharge, although some slight rise in the density can be observed after the transition to the ELM-free state, presumably linked to reduced pumping in the absence of ELMs (as also occurs before H-mode is achieved.)

Evidence of the ELMs can be seen in the D α traces presented in figure 2f. Notably, the plasma remains in H-mode with rapid Type-III ELMs as the triangularity drops below δ = 0 until a transition to an ELM-free state is made at around t crit ∼ 710 ms when δ crit ≲ -0 06. . The observation of ELM-free behavior below a critical triangularity δ crit is qualitatively consistent with expectations concerning the behavior of NT scenarios from the DIIID and TCV tokamaks [5, 6]. However, it remains to be conclusively demonstrated that this phenomena is a robust property of MAST-U NT discharges. Interestingly, infiniten stability modeling suggests that access to this ELM-free scenario may be governed by different physics than in DIII-D or TCV, as is discussed in section .

Over the course of discharge #49463, the normalized confinement H 98y2 , as calculated with the TRANSP code [24], decreases steadily with decreasing κ (figure 2e) from H 98y2 ∼ 0 9 at 300 ms to . H 98y2 ∼ 0 66 during the ELM-. free phase. No significant changes in H 98y2 are observed near δ = 0 or δ = δ crit . The normalized pressure β N (figure 2b) drops significantly from β N ∼ 2 4 to . β N ∼ 2 when the inner strikeline moves across the divertor opening but remains constant at β N ∼ 2 for the rest of the discharge. Notably, little change in the impurity content is observed across this discharge, as indicated by the carbon emission measurements captured with the SPRED diagnostic shown in figure 2d [25]. Slight ( ∼ 10 -20%) increases in the line-averaged density and the C and He SPRED emissions are observed during the ELM-free phase, potentially suggestive of reduced impurity transport under ELM-free conditions.

Representative traces of the electron temperature ( T e ), density ( n e ) and pressure ( p e ) at the plasma core ( ψ N ∼ 0 25) and near the pedestal top ( . ψ N ∼ 0 92) are shown . in figure 3. Near the pedestal top, T e ped , drops significantly both when the inner strikeline transitions to the outboard side of the vessel around t ∼ 420 ms and then again when the ELM-free state is achieved at t ∼ 710 ms. The edge density n e ped , stays roughly constant as the strikeline shifts, but drops upon the transition to the ELM-free state before growing again towards the end of the discharge. Notably, the core traces display markedly different trends: while n e core , remains roughly constant after H-mode is achieved, T e core , actually rises steadily as

FIG. 3: The (a) electron temperature, (b) electron density and (c) electron pressure as a function of time for two radial locations: the core (blue) and the pedestal top (red). Note that the core profiles are scaled down by constant factors for visibility.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00180v1-figure-a1fc9f4b-3d10-4d8c-9396-90241854a504.jpg)

δ is decreased, reaching its highest values in the ELMfree portion of the discharge when δ is the most negative. This behavior is reproduced qualitatively in initial simulations with the STEP code [26], which recovers an increase in the core T e as δ is decreased at fixed density and injected power that results from a reduced TGLF/TGYRO SAT-2 [27] electron heat flux at more negative triangularity. The rise in core temperature at constant density translates to a steady rise in the core pressure as δ is decreased, as seen in figure 3c. Similar qualitative behavior demonstrating a decrease in edge pressure but a retention of high core pressures has been observed in NT plasmas on DIII-D [4, 6] and TCV [2, 28] and is often associated with a reduction in the electron thermal conductivity in NT scenarios.

Further characterization of the plasma performance (including transport mechanisms, β N limits and confinement scalings) as a function of δ &lt; 0 requires additional investigations over a dedicated set of NT discharges on MAST-U and are thus not considered here. Exploration of the heat-flux width, ELM behavior and shaping limits will also be included in future studies.

## DECORRELATION OF SECOND STABILITY CLOSURE AND ELM-FREE NT OPERATION

On the DIII-D and TCV tokamaks, which have conventional aspect ratios, access to ELM-free conditions via NT shaping occurs suddenly at a critical triangularity ( δ crit ). This phenomena has been strongly correlated with the closure of the second stability region for infiniten ballooning modes, which introduces an additional transport mechanism at the plasma edge, preventing the growth of steep pressure gradients in the NT edge [6, 8, 14, 17, 29]. However, as mentioned in [8, 29], it is in general not necessarily sufficient in tokamak plasmas that second stability access be closed to ensure L-mode operation. Broadly speaking, if the collisionality is high enough to suppress the bootstrap current, local shear can remain high enough to keep the edge in the first stable region even when weak peeling-ballooning modes are destabilized. This effect is particularly evident in devices with low aspect ratio, and it is not uncommon for MASTU H-mode discharges with positive triangularity to form pedestals that reside in the first-stable region [18]. On NSTX, the prevalence of first-stable H-modes can be observed in the bifurcation of pedestal widths into 'wide' and 'narrow' branches, which are in turn set by variations in the stability of kinetic ballooning modes [19, 30].

The ideal ballooning stability of tokamak discharges can be easily assessed using the BALOO code [11], which requires as input well-converged kinetic equilibrium reconstructions. For this work, kinetic equilibria were generated every 10 ms using the kineticEFITtime workflow provided in the OMFIT framework [31]. This workflow allows for iteration between profile fits to the electron and ion profiles and a free-boundary EFIT reconstruction [32]. The stability results from this study are presented in figure 4, which shows the normalized pressure gradient ( α ) corresponding to the stability boundaries (dashed lines) and experimental conditions (solid lines) for five key times throughout discharge #49463.

In figure 4, the first four times (yellow - 300 ms, orange - 400 ms, magenta - 480 ms and red - 590 ms) are all from periods of the discharge when the plasma is firmly in a Type-III ELMy regime. No stability boundary is found by the BALOO code during the earliest times ( t ≲ 400 ms), allowing the pressure gradient at the edge to grow to high values. However, the stability boundary rapidly closes around t ∼ 480 ms when the average triangularity measured at the last closed flux surface (LCFS) nears zero. Throughout the evolution of this discharge, the peak pressure gradient in the edge region decreases as δ is decreased. ELMy H-mode behavior is robustly sustained until a triangularity of δ ∼ -0 06 is . reached around t = 710ms (black lines in figure 4). An ELM-free time slice ( t = 800ms) is shown in blue for comparison. Notably, max( α ) remains similar for the entire ELM-free phase ( t = 720 -830 ms) and does not continue to drop

FIG. 4: (a) A zoom-in of the divertor region and (b) the corresponding infiniten ballooning stability for a series of times during MAST-U discharge #49463. ELMy timeslices are colored yellow, orange, red and magenta. The stability region (dashed lines) is open originally but closes well in advance of the ELM-free transition (black). The final ELM-free state is shown in blue. At no point during the evolution of this discharge is the experimental normalized pressure gradient (solid lines) directly limited by the first stability limit for infiniten modes.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00180v1-figure-2bed5c48-b890-4a8e-9871-cb492ec1358b.jpg)

as δ is further decreased during this time.

Strikingly, this analysis reveals a decorrelation between the closure of ballooning stability window and the prevention of ELMs by NT shaping. While previously suggested as a possible complication to the no-ELM NT scenario [8, 29, 33], this has not yet been observed on either the DIII-D or TCV tokamaks, where closure of second stability typically coincides with the loss of H-mode in NT [14, 17, 29, 34].

To compensate for any uncertainties in the equilibrium reconstructions that could impact this result, over 500 variations of the kinetic equilibrium representing the plasma at t = 500ms were generated by solving the fixedboundary using the CHEASE code [35] with arbitrary profile modifications, following the workflow established in [8]. Scaling the edge pressure gradient, collisionality and

FIG. 5: Calculations of the ballooning-critical pedestal (BCP; red, dashed) and gyrokinetic-critical pedestal (GCP; blacked, solid) for MAST-U discharge #49463 at 600 ms, which is ELMy at modest δ &lt; 0 with closed access to the second stability region.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00180v1-figure-42dd85e9-5b12-4f9b-8d93-5aef55dee684.jpg)

current separately by factors of five was not sufficient to open access to the second stability region in this shape, suggesting that the conclusion that this MAST-U discharge achieves a first-stable H-mode regime at low values of δ &lt; 0 is not an artifact of any workflow.

As seen in figure 2f, an ELM-free state is still achieved in this low aspect ratio discharge, suggesting that additional physics may be needed to fully explain the transition to ELM-free states via NT shaping over a wider variety of plasma conditions. As the plasma conditions ( n e , P inj and impurity content) were held constant over this scan, it is likely that the transition to an ELM-free state in discharge #49463 is still largely a result of the shape changes. Interestingly, and as a result of the unique shapes achieved in this discharge, while the triangularity measured at the last closed flux surface drops below zero around t = 590ms, the triangularity just inside of the LCFS (up to ψ N ∼ 0 9) remains positive until the ELM-. free transition at t = 710ms. This may suggest that the effects of NT must extend over a significant portion of the edge region to inhibit the transition to H-mode, which is likely true in conventional aspect ratio scenarios as well.

To begin to assess possible mechanisms for pedestal height reduction responsible for the ELM-free behavior observed at strong NT on MAST-U, the gk ped code [36] was utilized to determine the Gyrokinetic Critical Pedestal (GCP) in discharge #49463. This model combines linear gyrokinetics with self-consistent pedestal equilibrium variation to explicitly determine stability limits stemming from the Kinetic Ballooning Mode (KBM) in the plasma edge and has recently been utilized to reproduce experimental measurements of the pedestal width-height scalings in the NSTX tokamak, which also has low-aspect-ratio [19]. This GCP is calculated for an ELMing case with δ &gt; 0 (600 ms) in figure 5, showing excellent agreement with the experimentally measured pedestal height ( β θ, ped ) and width (∆ ped ), which are calculated from a modified-tanh fit to the edge pressure profile. Notably, the Ballooning Critical Pedestal (BCP), which is calculated directly from the infiniten ballooning stability boundary, lies well above this curve, in agreement with the BALOO modeling presented in figure 4.

This result carries several important implications. First, it offers insight into the Type-III ELMy nature of these discharges. Initial investigations with the ELITE code [37] suggest that these plasmas sit well within the nominal stability region for Type-I ELMs - a reduction in the diamagnetic stabilization factor by a factor of ∼ 1 / 10 from typical MAST-U values is needed to localize the experimental measurements along the predicted stability boundary. If the plasma is indeed limited by a KBM constraint that sits below the peeling-ballooning boundary, as is suggested by gk ped , it would offer an justification for the reduced pedestal height and Type-III ELMy operation. Second, the prevalence of this KBM limit lends some credibility to the theoretical arguments made in [12], which finds that the KBM will be a primary source of transport in ST NT scenarios with reactor-relevant parameters. Further study of the role of KBMs in setting the β limits for NT scenarios on STs should be the subject of future work.

Finally, the reduction of the modeled GCP limit from the BCP provides a potential explanation for the decorrelation of second stability closure and ELM-free NT operation that distinguishes this regime from counterparts on conventional aspect ratio tokamaks. While the mechanisms for ELM suppression are still unclear, reduction of the pedestal height below the first stability limit for infiniten ballooning modes occurs already in the ELMy phase of this discharge as a result of the KBM limit. While potentially more applicable in low-aspect ratio devices [19, 30], such a reduction due to kinetic effects could also be present in conventional aspect ratio devices and offer a potential avenue for gradient reduction below the ideal ballooning limit, as observed in [6, 29]. Future work, including similar KBM analysis on conventional aspect ratio machines and further study of the ELM-free NT regime on MAST-U across a broader parameter range, is needed to fully explore the validity of this model.

## DISCUSSION AND CONCLUSION

This Letter reports on the first NT plasma achieved in a spherical tokamak. While further experiments, including on other STs like NSTX-U and SMART, are needed to fully characterize the performance and stability of this regime, several interesting observations can already be made due to the fortuitously stable plasma conditions achieved in this discharge. First, the qualitative behavior

expected from NT operation on other machines concerning a sharp transition to an ELM-free regime without significant degradation of plasma performance is achieved. While H 98y2 decreases across the shape scan with decreasing κ , no obvious decrease in β N is associated with the transition to a shape with δ &lt; 0 or when achieving H-mode inhibition. This optimistically suggests that NT operation of low aspect ratio devices follow the intuition developed on larger aspect ratio devices, though gyrokinetic predictions do predict degradation of plasma performance at very low A [12, 13]. Appropriate characterization of the confinement properties of the plasmas demands additional experimental investigation.

Perhaps more interestingly, at less negative triangularities an extended period of Type-III ELMing H-mode persists even though second stability for ideal ballooning modes is closed. The decoupling of second stability closure and access to the NT ELM-free state has not yet been observed on either DIII-D or TCV, which typically transition directly from a Type-I ELMy H-mode into the NT ELM-free state when the ballooning stability bifurcates. Access to first-stable ELMy H-mode at small absolute values of δ &lt; 0 could have significant consequences on the feasibility of this approach as a pilot plant scenario for low aspect ratio machines. This supports arguments presented in reference [8] regarding consideration of the scalings presented therein only as bounding values on the required shaping parameters for access to the ELM-free state. If first-stable H-modes with Type-III ELMy can be achieved in a given configuration, more negative values of δ may be needed to achieve robust H-mode inhibition.

We note further that this key result - that additional physics mechanisms beyond the ideal ballooning stability likely impact the edge gradient limits in NT plasmas is not necessarily inconsistent with results obtained on other machines at larger aspect ratio [29]. On DIII-D, while access to ELM-free conditions is strongly correlated with the closure of ideal ballooning stability [14], a majority of ELM-free NT discharges sustain gradients well below the maximum gradient allowed by first stability limit for ballooning modes. This can be observed clearly in figure 3 of reference [6], which shows numerous ELMfree DIII-D equilibria sitting well under the first stability limit. Taken together, these results are indicative of additional transport mechanisms present in the NT edge that can both limit the edge gradient and robustly prevent ELMs. Future work on this subject should aim to identify this physics to enable reliable predictions of the NT edge under the conditions expected in future fusion power plants.

## ACKNOWLEDGMENTS

Some calculations performed in this study were completed through the OMFIT framework [31], and the generated data are available upon request. This material is based upon work supported by the U.S. Department of Energy, Office of Science, under Awards DESC0022272, DE-SC0018991, DE-AC05-00OR22725 and DE-FC02-04ER54698. This work has been funded in part by the EPSRC Energy Programme [grant number EP/W006839/1]. To obtain further information on the data and models underlying this paper please contact PublicationsManager@ukaea.uk.

- ∗ Corresponding author: a.o.nelson@columbia.edu
- [1] A. Marinoni, O. Sauter, and S. Coda. A brief history of negative triangularity tokamak plasmas. 5(1):6. ISSN 2367-3192. doi:10.1007/s41614-021-00054-0. URL https: //doi.org/10.1007/s41614-021-00054-0 .
- [2] Y. Camenen, A. Pochelon, R. Behn, A. Bottino, A. Bortolon, S. Coda, A. Karpushov, O. Sauter, and G. Zhuang. Impact of plasma triangularity and collisionality on electron heat transport in TCV l-mode plasmas. 47(7):510516. ISSN 00295515. doi:10.1088/0029-5515/47/7/002. URL https://doi.org/10.1088/0029-5515/47/7/002 .
- [3] M. E. Austin, A. Marinoni, M. L. Walker, M. W. Brookman, J. S. Degrassie, A. W. Hyatt, G. R. McKee, C. C. Petty, T. L. Rhodes, S. P. Smith, C. Sung, K. E. Thome, and A. D. Turnbull. Achievement of reactorrelevant performance in negative triangularity shape in the DIII-d tokamak. 122:115001. ISSN 10797114. doi: 10.1103/PhysRevLett.122.115001. URL https://doi. org/10.1103/PhysRevLett.122.115001 .
- [4] C. Paz-Soldan, C. Chrystal, P. Lunia, A. O. Nelson, M. E. Austin, T. B. Cote, A. W. Hyatt, A. Marinoni, T. H. Osborne, M. Pharr, O. Sauter, F. Scotti, K. E. Thome, T. M. Wilks, and H. S. Wilson. Simultaneous access to high normalized current, pressure, density, and confinement in strongly-shaped diverted negative triangularity plasmas. page in review.
- [5] S Coda, A Merle, O Sauter, L Porte, F Bagnato, J Boedo, T Bolzonella, O F´ evrier, B Labit, A Marinoni, A Pau, L Pigatto, U Sheikh, C Tsui, M Vallar, and T Vu. Enhanced confinement in diverted negative-triangularity lmode plasmas in TCV. 64(14004). ISSN 0741-3335, 1361-6587. doi:10.1088/1361-6587/ac3fec. URL https: //doi.org/10.1088/1361-6587/ac3fec .
- [6] A O Nelson, L Schmitz, C Paz-Soldan, K E Thome, T B Cote, N Leuthold, F Scotti, M E Austin, A Hyatt, and T Osborne. Robust avoidance of edge-localized modes alongside gradient formation in the negative triangularity tokamak edge. 131:195101, . URL https://doi.org/10. 1103/PhysRevLett.131.195101 .
- [7] S. Yu Medvedev, M. Kikuchi, L. Villard, T. Takizuka, P. Diamond, H. Zushi, K. Nagasaki, X. Duan, Y. Wu, A. A. Ivanov, A. A. Martynov, Yu Yu Poshekhonov, A. Fasoli, and O. Sauter. The negative triangularity tokamak: Stability limits and prospects as a fusion energy system. 55:063013, . ISSN 17414326. doi:10.1088/00295515/55/6/063013. URL https://doi.org/10.1088/ 0029-5515/55/6/063013 .
- [8] A. O. Nelson, Carlos Paz-Soldan, and Samuli Saarelma. Prospects for h-mode inhibition in negative triangularity

- tokamak reactor plasmas. 62:096020, . ISSN 0029-5515, 1741-4326. doi:10.1088/1741-4326/ac8064. URL https: //doi.org/10.1088/1741-4326/ac8064 .
- [9] G Rutherford and et al. MANTA: A negativetriangularity NASEM-compliant fusion pilot plant. MANTA: A Negative-Triangularity NASEM-Compliant Fusion Pilot Plant:in review.
- [10] M. Kikuchi, T. Takizuka, S. Medvedev, T. Ando, D. Chen, J. X. Li, M. Austin, O. Sauter, L. Villard, A. Merle, M. Fontana, Y. Kishimoto, and K. Imadera. L-mode-edge negative triangularity tokamak reactor. 59: 056017. ISSN 17414326. doi:10.1088/1741-4326/ab076d. URL https://doi.org/10.1088/1741-4326/ab076d .
- [11] MAMiller and et al. Power handling solutions for a negative triangularity radiative l-mode pilot plant. Power Handling Solutions for a Negative Triangularity Radiative L-mode Pilot Plant:in review.
- [12] R Davies, D Dickinson, and H Wilson. Kinetic ballooning modes as a constraint on plasma triangularity in commercial spherical tokamaks. 64(10):105001. ISSN 0741-3335, 1361-6587. doi:10.1088/1361-6587/ac8615. URL https://iopscience.iop.org/article/10.1088/ 1361-6587/ac8615 .
- [13] A Balestri, J Ball, S Coda, D J Cruz-Zabala, M GarciaMunoz, and E Viezzer. Physical insights from the aspect ratio dependence of turbulence in negative triangularity plasmas. 66(7):075012. ISSN 0741-3335, 1361-6587. doi: 10.1088/1361-6587/ad4d1d. URL https://iopscience. iop.org/article/10.1088/1361-6587/ad4d1d .
- [14] S Saarelma, M E Austin, M Knolker, A Marinoni, C PazSoldan, L Schmitz, and P B Snyder. Ballooning instability preventing the h-mode access in plasmas with negative triangularity shape on the DIII-d tokamak. 63 (105006). doi:10.1088/1361. URL https://doi.org/10. 1088/1361-6587/ac1ea4 .
- [15] C.M. Bishop. Stability of localized MHD modes in divertor tokamaks -a picture of the h-mode. 26 (8):1063-1072. ISSN 0029-5515, 1741-4326. doi: 10.1088/0029-5515/26/8/006. URL https://doi.org/ 10.1088/0029-5515/26/8/006 .
- [16] S.A. Sabbagh, M.H. Hughes, M.W. Phillips, A.M.M. Todd, and G.A. Navratil. Transition to the second region of ideal MHD stability. 29(3):423-435. ISSN 0029-5515, 1741-4326. doi:10.1088/0029-5515/29/3/006. URL https://iopscience.iop.org/article/10.1088/ 0029-5515/29/3/006 .
- [17] A. Merle, O. Sauter, and S. Yu Medvedev. Pedestal properties of h-modes with negative triangularity using the EPED-CH model. 59(10). ISSN 13616587. doi: 10.1088/1361-6587/aa7ac0. URL https://doi.org/10. 1088/1361-6587/aa7ac0 .
- [18] D Dickinson, S Saarelma, R Scannell, A Kirk, C M Roach, and H R Wilson. Towards the construction of a model to describe the inter-ELM evolution of the pedestal on MAST. 53(11):115010. ISSN 07413335, 1361-6587. doi:10.1088/0741-3335/53/11/115010. URL https://iopscience.iop.org/article/10.1088/ 0741-3335/53/11/115010 .
- [19] J.F. Parisi, W. Guttenfelder, A.O. Nelson, R. Gaur, A. Kleiner, M. Lampert, G. Avdeeva, J.W. Berkery, C. Clauser, M. Curie, A. Diallo, W. Dorland, S.M. Kaye, J. McClenaghan, and F.I. Parra. Kineticballooning-limited pedestals in spherical tokamak plasmas. 64(5):054002, . ISSN 0029-5515, 1741-4326. doi:
- 10.1088/1741-4326/ad39fb. URL https://iopscience. iop.org/article/10.1088/1741-4326/ad39fb .
- [20] C. Hansen, I G Stewart, D. Burgess, M. Pharr, S. Guizzo, F. Logak, A. O. Nelson, and C. Paz-Soldan. TokaMaker: An open-source time-dependent grad-shafranov tool for the design and modeling of axisymmetric fusion devices. 298:109111. URL https://doi.org/10.1016/ j.cpc.2024.109111 .
- [21] D. A. Humphreys, J. R. Ferron, M. Bakhtiari, J. A. Blair, Y. In, G. L. Jackson, H. Jhang, R. D. Johnson, J. S. Kim, R. J. Lahaye, J. A. Leuer, B. G. Penaflor, E. Schuster, M. L. Walker, H. Wang, A. S. Welander, and D. G. Whyte. Development of ITERrelevant plasma control solutions at DIII-d. 47(8):943951. ISSN 00295515. doi:10.1088/0029-5515/47/8/028. URL https://iopscience.iop.org/article/10.1088/ 0029-5515/47/8/028 .
- [22] Graham McArdle, Luigi Pangione, and Martin Kochan. The MAST upgrade plasma control system. 159:111764. ISSN 09203796. doi:10.1016/j.fusengdes.2020.111764. URL https://linkinghub.elsevier.com/retrieve/ pii/S0920379620303124 .
- [23] H. Anand, W. Wehner, D. Eldon, A. Welander, Z. Xing, A. Lvovskiy, J. Barr, E. Cho, B. Sammuli, D. Humphreys, N. Eidietis, A. Leonard, M. Kochan, C. Vincent, G. McArdle, G. Cunningham, A. Thornton, J. Harrison, V. Soukhanovskii, and J. Lovell. Realtime plasma equilibrium reconstruction and shape control for the MAST upgrade tokamak. 64(8):086051. ISSN 0029-5515, 1741-4326. doi:10.1088/1741-4326/ad5c80. URL https://iopscience.iop.org/article/10.1088/ 1741-4326/ad5c80 .
- [24] Joshua Breslau, Marina Gorelenkova, Francesca Poli, Jai Sachdev, and Xingqiu Yuan. TRANSP. URL https: //dx.doi.org/10.11578/dc.20180627.4 .
- [25] R. J. Fonck, A. T. Ramsey, and R. V. Yelle. Multichannel grazing-incidence spectrometer for plasma impurity diagnosis: SPRED. 21(12):2115-2123. doi: 10.1364/AO.21.002115. URL https://opg.optica.org/ ao/abstract.cfm?URI=ao-21-12-2115 .
- [26] B. C. Lyons, J. McClenaghan, T. Slendebroek, O. Meneghini, T. F. Neiser, S. P. Smith, D. B. Weisberg, E. A. Belli, J. Candy, J. M. Hanson, L. L. Lao, N. C. Logan, S. Saarelma, O. Sauter, P. B. Snyder, G. M. Staebler, K. E. Thome, and A. D. Turnbull. Flexible, integrated modeling of tokamak stability, transport, equilibrium, and pedestal physics. 30(9):092510. ISSN 1070664X, 1089-7674. doi:10.1063/5.0156877. URL https: //pubs.aip.org/pop/article/30/9/092510/2911814/ Flexible-integrated-modeling-of-tokamak-stability

.

- [27] G. M. Staebler, J. E. Kinsey, and R. E. Waltz. A theory-based transport model with comprehensive physics. 14(5):055909. ISSN 1070-664X, 1089-7674. doi:10.1063/1.2436852. URL https: //pubs.aip.org/pop/article/14/5/055909/929642/ A-theory-based-transport-model-with-comprehensive

.

- [28] O. Sauter, S. Brunner, D. Kim, G. Merlo, R. Behn, Y. Camenen, S. Coda, B. P. Duval, L. Federspiel, T. P. Goodman, A. Karpushov, A. Merle, and Tcv Team. On the non-stiffness of edge transport in l-mode tokamak plasmas. 21(55906). ISSN 10897674. doi: 10.1063/1.4876612. URL https://doi.org/10.1063/1. 4876612 .

- [29] A. O. Nelson, L. Schmitz, T. Cote, J. F. Parisi, S. Stewart, C. Paz-Soldan, K. E. Thome, M. E. Austin, F. Scotti, J. L. Barr, A. Hyatt, N. Leuthold, A. Marinoni, T. Neiser, T. Osborne, N. Richner, A. S. Welander, W. P. Wehner, R. Wilcox, T. M. Wilks, and J. Yang. Characterization of the ELM-free negative triangularity edge on DIII-d. page in review, . URL http://arxiv.org/abs/2405.11082 .
- [30] J F Parisi, A O Nelson, R Gaur, S M Kaye, F I Parra, J W Berkery, K Barada, C Clauser, A J Creely, A Diallo, W Guttenfelder, J W Hughes, L A Kogan, A Kleiner, A Q Kuang, M Lampert, T Macwan, J E Menard, and M A Miller. Kinetic-ballooning-bifurcation in tokamak pedestals across shaping and aspect-ratio. 31:030702, . URL https://doi.org/10.1063/5.0190818 .
- [31] O. Meneghini, S. P. Smith, L. L. Lao, O. Izacard, Q. Ren, J. M. Park, J. Candy, Z. Wang, C. J. Luna, V. A. Izzo, B. A. Grierson, P. B. Snyder, C. Holland, J. Penna, G. Lu, P. Raum, A. McCubbin, D. M. Orlov, E. A. Belli, N. M. Ferraro, R. Prater, T. H. Osborne, A. D. Turnbull, and G. M. Staebler. Integrated modeling applications for tokamak experiments with OMFIT. 55:083008. ISSN 17414326. doi:10.1088/0029-5515/55/8/083008. URL http://dx.doi.org/10.1088/0029-5515/55/8/083008 .
- [32] L. L. Lao, H. E. St John, R. D. Stambuch, A. G. Kellman, and W. W. Pfeiffer. Reconstruction of current profile parameters and plasma shapes in tokamaks. 25:1611. doi: 10.1088/0029-5515/25/11/007. URL https://doi.org/ 10.1088/0029-5515/25/11/007 .
- [33] Chenchen Qin, Maolin Mou, Shaoyong Chen, Yongjian Li, Yang Luo, Limin Zhu, Liangkang Dong, and Changjian Tang. Effect of negative triangularity
- on peeling-ballooning instability. ISSN 0031-8949, 1402-4896. doi:10.1088/1402-4896/acbdcb. URL https://iopscience.iop.org/article/10.1088/ 1402-4896/acbdcb .
- [34] S. Yu Medvedev, A. A. Ivanov, A. A. Martynov, Yu Yu Poshekhonov, R Behn, Y. R. Martin, A. Pochelon, O. Sauter, and L. Villard. Beta limits and edge stability for negative triangularity plasma in TCV tokamak. 32D:1.072, . URL https: //www.researchgate.net/profile/Antoine-Pochelon/ publication/37468729\_Beta\_limits\_and\_Edge\_
- Stability\_for\_Negative\_Triangularity\_Plasma\_ in\_TCV\_Tokamak/links/5fd9af1645851553a0bd74f0/ Beta-limits-and-Edge-Stability-for-Negative-Triangularitypdf .
- [35] H. L¨tjens, A. Bondeson, and O. Sauter. u The CHEASE code for toroidal MHD equilibria. 97(3):219-260. ISSN 00104655. doi:10.1016/0010-4655(96)00046-X. URL https://doi.org/10.1016/0010-4655(96)00046-X .
- [36] J.F. Parisi, A. O. Nelson, and R. Gaur. gk ped. In (available at: https://github.com/gafusion/OMFITsource/tree/gk ped) , .
- [37] P. B. Snyder, H. R. Wilson, J. R. Ferron, L. L. Lao, A. W. Leonard, T. H. Osborne, A. D. Turnbull, D. Mossessian, M. Murakami, and X. Q. Xu. Edge localized modes and the pedestal: A model based on coupled peeling-ballooning modes. 9:2037-2043. ISSN 1070664X. doi:10.1063/1.1449463. URL https://doi.org/ 10.1063/1.1449463 .
A. O. Nelson C. Vincent H. Anand J. Lovell J. F. Parisi H. S. Wilson K. Imada W. P. Wehner M. Kochan S. Blackmore G. McArdle S. Guizzo L. Rondini S. Freiberger C. Paz-Soldan
First Access to ELM-free Negative Triangularity at Low Aspect Ratio
A plasma scenario with negative triangularity (NT) shaping is achieved on
MAST-U for the first time. While edge localized modes (ELMs) are eventually
suppressed as the triangularity is decreased below $\delta$ < -0.06, an
extended period of H-mode operation with Type-III ELMs is sustained at less
negative $\delta$ even through access to the second stability region for ideal
ballooning modes is closed. This documents a qualitative difference from the
ELM-free access conditions documented in NT scenarios on conventional aspect
ratio machines. The electron temperature at the pedestal top drops across the
transition to ELM-free operation, but a steady rise in core temperature as
$\delta$ is decreased allows for similar normalized beta in the ELM-free NT and
H-mode positive triangularity shapes.