Showing 1–50 of 91 results for author: Fülöp, T

Evolution of the autoresonant plasma wave excitation in two-dimensional particle-in-cell simulations

Authors: Mufei. Luo, Caterina. Riconda, Anna. Grassi, Ning. Wang, Jonathan S. Wurtele, Tünde Fülöp, István Pusztai

Abstract: The generation of an autoresonantly phase-locked high amplitude plasma waves to the chirped beat frequency of two driving lasers is studied in two dimensions using particle-in-cell simulations. The two-dimensional plasma and laser parameters correspond to those that optimized the plasma wave amplitude in one-dimensional simulations. Near the start of autoresonant locking, the two-dimensional simul… ▽ More The generation of an autoresonantly phase-locked high amplitude plasma waves to the chirped beat frequency of two driving lasers is studied in two dimensions using particle-in-cell simulations. The two-dimensional plasma and laser parameters correspond to those that optimized the plasma wave amplitude in one-dimensional simulations. Near the start of autoresonant locking, the two-dimensional simulations appear similar to one-dimensional particle-in-cell results [Luo et al., Phys. Rev. Res. 6, 013338 (2024)] with plasma wave amplitudes above the Rosenbluth-Liu limit. Later, just below wave-breaking, the two-dimensional simulation exhibits a Weibel-like instability and eventually laser beam filamentation. These limit the coherence of the plasma oscillation after the peak plasma wave field is obtained. In spite of the reduction of spatial coherence of the accelerating density structure, the acceleration of self-injected electrons in the case studied remains at $70\%$ to $80\%$ of that observed in one dimension. Other effects such as plasma wave bowing are discussed. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 23 pages, 13 figures

Control of autoresonant plasma beat-wave wakefield excitation

Authors: M. Luo, C. Riconda, I. Pusztai, A. Grassi, J. S. Wurtele, T. Fülöp

Abstract: Autoresonant phase-locking of the plasma wakefield to the beat frequency of two driving lasers offers advantages over conventional wakefield acceleration methods, since it requires less demanding laser parameters and is robust to variations in the target plasma density. Here, we investigate the kinetic and nonlinear processes that come into play during autoresonant plasma beat-wave acceleration of… ▽ More Autoresonant phase-locking of the plasma wakefield to the beat frequency of two driving lasers offers advantages over conventional wakefield acceleration methods, since it requires less demanding laser parameters and is robust to variations in the target plasma density. Here, we investigate the kinetic and nonlinear processes that come into play during autoresonant plasma beat-wave acceleration of electrons, their impact on the field amplitude of the accelerating structure, and on acceleration efficiency. Particle-in-Cell simulations show that the process depends on the plasma density in a non-trivial way but can be reliably modeled under specific conditions. Beside recovering previous fluid results in the deeply underdense plasma limit, we demonstrate that robust field excitation can be achieved within a fully kinetic self-consistent modeling. By adjusting the laser properties, we can amplify the electric field to the desired level, up to wave-breaking, and efficiently accelerate particles; we provide suggestions for optimized laser and plasma parameters. This versatile and efficient acceleration scheme, producing electrons from tens to hundreds of MeV energies, holds promise for a wide range of applications in research industry and medicine. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 14 pages, 13 figures

Fluid and kinetic studies of tokamak disruptions using Bayesian optimization

Authors: Ida Ekmark, Mathias Hoppe, Tünde Fülöp, Patrik Jansson, Liam Antonsson, Oskar Vallhagen, Istvan Pusztai

Abstract: When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter we derive and implement kinetic sources for… ▽ More When simulating runaway electron dynamics in tokamak disruptions, fluid models with lower numerical cost are often preferred to more accurate kinetic models. The aim of this work is to compare fluid and kinetic simulations of a large variety of different disruption scenarios in ITER. We consider both non-activated and activated scenarios; for the latter we derive and implement kinetic sources for the Compton scattering and tritium beta decay runaway electron generation mechanisms in our simulation tool DREAM [M. Hoppe et al 2021 Comp. Phys. Commun. 268, 108098]. To achieve a diverse set of disruption scenarios, Bayesian optimization is used to explore a range of massive material injection densities for deuterium and neon. The cost function is designed to distinguish between successful and unsuccessful disruption mitigation based on the runaway current, current quench time and transported fraction of the heat loss. In the non-activated scenarios, we find that fluid and kinetic disruption simulations can have significantly different runaway electron dynamics, due to an overestimation of the runaway seed by the fluid model. The primary cause of this is that the fluid hot-tail generation model neglects superthermal electron transport losses during the thermal quench. In the activated scenarios, the fluid and kinetic models give similar predictions, which can be explained by the significant influence of the activated sources on the RE dynamics and the seed. △ Less

Submitted 5 July, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 25 pages, 11 figures

Journal ref: Journal of Plasma Physics, 90(3), 905900306 (2024)

Runaway Electron Dynamics in ITER Disruptions with Shattered Pellet Injections

Authors: Oskar Vallhagen, Lise Hanebring, Javier Artola, Michael Lehnen, Eric Nardon, Tünde Fülöp, Mathias Hoppe, Sarah Newton, Istvan Pusztai

Abstract: This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well a… ▽ More This study systematically explores the parameter space of disruption mitigation through shattered pellet injection in ITER with a focus on runaway electron dynamics, using the disruption modelling tool DREAM. The physics fidelity is considerably increased compared to previous studies, by e.g., using realistic magnetic geometry, resistive wall configuration, thermal quench onset criteria, as well as including additional effects, such as ion transport and enhanced runaway electron transport during the thermal quench. The work aims to provide a fairly comprehensive coverage of experimentally feasible scenarios, considering plasmas representative of both non-activated and high-performance DT operation, different thermal quench onset criteria and transport levels, a wide range of hydrogen and neon quantities injected in one or two stages, and pellets with various characteristic shard sizes. Using a staggered injection scheme, with a pure hydrogen injection preceding a mixed hydrogen-neon injection, we find injection parameters leading to acceptable runaway electron currents in all investigated discharges without activated runaway sources. Dividing the injection into two stages is found to significantly enhance the assimilation and minimize runaway electron generation due to the hot-tail mechanism. However, while a staggered injection outperforms a single stage injection also in cases with radioactive runaway electron sources, no cases with acceptable runaway electron currents are found for a DT-plasma with a 15 MA plasma current. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 16 pages, 7 figures, submitted to Nuclear Fusion

Bayesian optimization of massive material injection for disruption mitigation in tokamaks

Authors: Istvan Pusztai, Ida Ekmark, Hannes Bergström, Peter Halldestam, Patrik Jansson, Mathias Hoppe, Oskar Vallhagen, Tünde Fülöp

Abstract: A Bayesian optimization framework is used to investigate scenarios for disruptions mitigated with combined deuterium and neon injection in ITER. The optimization cost function takes into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The aim is to explore the dependence of the cost function on injected densities, and provide in… ▽ More A Bayesian optimization framework is used to investigate scenarios for disruptions mitigated with combined deuterium and neon injection in ITER. The optimization cost function takes into account limits on the maximum runaway current, the transported fraction of the heat loss and the current quench time. The aim is to explore the dependence of the cost function on injected densities, and provide insights into the behaviour of the disruption dynamics for representative scenarios. The simulations are conducted using the numerical framework DREAM (Disruption Runaway Electron Analysis Model). We show that irrespective of the quantities of the material deposition, multi-megaampere runaway currents will be produced in the deuterium-tritium phase of operations, even in the optimal scenarios. However, the severity of the outcome can be influenced by tailoring the radial profile of the injected material; in particular if the injected neon is deposited at the edge region it leads to a significant reduction of both the final runaway current and the transported heat losses. The Bayesian approach allows us to map the parameter space efficiently, with more accuracy in favorable parameter regions, thereby providing us information about the robustness of the optima. △ Less

Submitted 24 February, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Drift of ablated material after pellet injection in a tokamak

Authors: O. Vallhagen, I. Pusztai, P. Helander, S. L. Newton, T. Fülöp

Abstract: Pellet injection is used for fuelling and controlling discharges in tokamaks, and it is foreseen in ITER. During pellet injection, a movement of the ablated material towards the low-field side (or outward major radius direction) occurs because of the inhomogeneity of the magnetic field. Due to the complexity of the theoretical models, computer codes developed to simulate the cross-field drift are… ▽ More Pellet injection is used for fuelling and controlling discharges in tokamaks, and it is foreseen in ITER. During pellet injection, a movement of the ablated material towards the low-field side (or outward major radius direction) occurs because of the inhomogeneity of the magnetic field. Due to the complexity of the theoretical models, computer codes developed to simulate the cross-field drift are computationally expensive. Here, we present a one-dimensional semi-analytical model for the radial displacement of ablated material after pellet injection, taking into account both the Alfvén and ohmic currents which short-circuit the charge separation creating the drift. The model is suitable for rapid calculation of the radial drift displacement, and can be useful for e.g. modelling of disruption mitigation via pellet injection. △ Less

Submitted 24 April, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: 22 pages, 4 figures. Submitted to Journal of Plasma Physics

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

Authors: R. A. Tinguely, I. Pusztai, V. A. Izzo, K. S'"arkimäki, T. Fülöp, D. T. Garnier, R. S. Granetz, M. Hoppe, C. Paz-Soldan, A. Sundström, R. Sweeney

Abstract: In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2… ▽ More In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1 MA compared to ~8.7 MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.3 and with kinetic energies 0.2-15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 $\mathrm{m^2/s}$, is needed in the core to fully mitigate these REs. Properly accounting for (i) the CQ electric field's effect on RE transport in islands and (ii) the contribution of significant RE currents to disruption MHD may help achieve this. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Thermodynamically extended symplectic numerical simulation of viscoelastic, thermal expansion and heat conduction phenomena in solids

Authors: Donát M. Takács, Áron Pozsár, Tamás Fülöp

Abstract: Symplectic numerical schemes for reversible dynamical systems predict the solution reliably over large times as well, and are a good starting point for extension to schemes for simulating irreversible situations like viscoelastic wave propagation and heat conduction coupled via thermal expansion occuring in rocks, plastics, biological samples etc. Dissipation error (artificial nonpreservation of e… ▽ More Symplectic numerical schemes for reversible dynamical systems predict the solution reliably over large times as well, and are a good starting point for extension to schemes for simulating irreversible situations like viscoelastic wave propagation and heat conduction coupled via thermal expansion occuring in rocks, plastics, biological samples etc. Dissipation error (artificial nonpreservation of energies and amplitudes) of the numerical solution should be as small as possible since it should not be confused with the real dissipation occuring in the irreversible system. In addition, the other well-known numerical artefact, dispersion error (artificial oscillations emerging at sharp changes), should also be minimal to avoid confusion with the true wavy behaviour. The continuum thermodynamical aspects (respect for balances with fluxes, systematic constitutive relationships between intensive quantities and fluxes, the second law of thermodynamics with positive definite entropy production, and the spacetime-based kinematic viewpoint) prove valuable for obtaining such extended schemes and for monitoring the solutions. Generalizing earlier works in this direction, here, we establish and investigate such a numerical scheme for one-dimensional viscoelastic wave propagation in the presence of heat conduction coupled via thermal expansion, demonstrating long-term reliability and the applicability of thermodynamics-based quantities in supervising the quality of the solution. △ Less

Submitted 26 October, 2023; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: 13 pages, 7 figures, 3 tables

Journal ref: Continuum Mechanics and Thermodynamics, 2024

Runaway dynamics in reactor-scale spherical tokamak disruptions

Authors: E. Berger, I. Pusztai, S. L. Newton, M. Hoppe, O. Vallhagen, A. Fil, T. Fülöp

Abstract: Understanding generation and mitigation of runaway electrons in disruptions is important for the safe operation of future tokamaks. In this paper we investigate runaway dynamics in reactor-scale spherical tokamaks. We study both the severity of runaway generation during unmitigated disruptions, as well as the effect that typical mitigation schemes based on massive material injection have on runawa… ▽ More Understanding generation and mitigation of runaway electrons in disruptions is important for the safe operation of future tokamaks. In this paper we investigate runaway dynamics in reactor-scale spherical tokamaks. We study both the severity of runaway generation during unmitigated disruptions, as well as the effect that typical mitigation schemes based on massive material injection have on runaway production. The study is conducted using the numerical framework DREAM (Disruption Runaway Electron Analysis Model). We find that, in many cases, mitigation strategies are necessary to prevent the runaway current from reaching multi-megaampere levels. Our results indicate that with a suitably chosen deuterium-neon mixture for mitigation, it is possible to achieve a tolerable runaway current and ohmic current evolution. With such parameters, however, the majority of the thermal energy loss happens through radial transport rather than radiation, which poses a risk of unacceptable localised heat loads. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Comments: 19 pages, 7 figures

Bayesian approach for validation of runaway electron simulations

Authors: Aaro Järvinen, Tünde Fülöp, Eero Hirvijoki, Mathias Hoppe, Adam Kit, Jan Åström

Abstract: Plasma-terminating disruptions in future fusion reactors may result in conversion of the initial current to a relativistic runaway electron beam. Validated predictive tools are required to optimize the scenarios and mitigation actuators to avoid the excessive damage that can be caused by such events. Many of the simulation tools applied in fusion energy research require the user to specify several… ▽ More Plasma-terminating disruptions in future fusion reactors may result in conversion of the initial current to a relativistic runaway electron beam. Validated predictive tools are required to optimize the scenarios and mitigation actuators to avoid the excessive damage that can be caused by such events. Many of the simulation tools applied in fusion energy research require the user to specify several input parameters that are not constrained by the available experimental information. Hence, a typical validation exercise requires multiparameter optimization to calibrate the uncertain input parameters for the best possible representation of the investigated physical system. The conventional approach, where an expert modeler conducts the parameter calibration based on domain knowledge, is prone to lead to an intractable validation challenge. For a typical simulation, conducting exhaustive multiparameter investigations manually to ensure a globally optimal solution and to rigorously quantify the uncertainties is an unattainable task, typically covered only partially and unsystematically. Bayesian inference algorithms offer a promising alternative approach that naturally includes uncertainty quantification and is less subjective to user bias in choosing the input parameters. The main challenge in using these methods is the computational cost of simulating enough samples to construct the posterior distributions for the uncertain input parameters. This challenge can be overcome by combining probabilistic surrogate modelling, such as Gaussian Process regression, with Bayesian optimization, which can reduce the number of required simulations by several orders of magnitude. Here, we implement this type of Bayesian optimization framework for a model for analysis of disruption runaway electrons, and explore for simulations of current quench in a JET plasma discharge with an argon induced disruption. △ Less

Submitted 3 August, 2022; originally announced August 2022.

Runaway electron generation during tokamak start-up

Authors: M. Hoppe, I. Ekmark, E. Berger, T. Fülöp

Abstract: Tokamak start-up is characterized by low electron densities and strong electric fields, in order to quickly raise the plasma current and temperature, allowing the plasma to fully ionize and magnetic flux surfaces to form. Such conditions are ideal for the formation of superthermal electrons, which may reduce the efficiency of ohmic heating and prevent the formation of a healthy thermal fusion plas… ▽ More Tokamak start-up is characterized by low electron densities and strong electric fields, in order to quickly raise the plasma current and temperature, allowing the plasma to fully ionize and magnetic flux surfaces to form. Such conditions are ideal for the formation of superthermal electrons, which may reduce the efficiency of ohmic heating and prevent the formation of a healthy thermal fusion plasma. This is of particular concern in ITER where engineering limitations put restrictions on the allowable electric fields and limit the prefill densities during start-up. In this study, we present a new 0D burn-through simulation tool called STREAM (STart-up Runaway Electron Analysis Model), which self-consistently evolves the plasma density, temperature and electric field, while accounting for the generation and loss of relativistic runaway electrons. After verifying the burn-through model, we investigate conditions under which runaway electrons can form during tokamak start-up as well as their effects on the plasma initiation. We find that Dreicer generation plays a crucial role in determining whether a discharge becomes runaway-dominated or not, and that a large number of runaway electrons could limit the ohmic heating of the plasma, thus preventing successful burn-through or further ramp-up of the plasma current. The runaway generation can be suppressed by raising the density via gas fuelling, but only if done sufficiently early. Otherwise a large runaway seed may have already been built up, which can avalanche even at relatively low electric fields and high densities. △ Less

Submitted 29 June, 2022; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 23 pages, 8 figures

Journal ref: Journal of Plasma Physics, 88(3), 905880317 (2022)

Attosecond dispersion as a diagnostics tool for solid-density laser-generated plasmas

Authors: Andréas Sundström, István Pusztai, Per Eng-Johnsson, Tünde Fülöp

Abstract: Extreme-ultraviolet pulses can propagate through ionised solid-density targets, unlike optical pulses and, thus, have the potential to probe the interior of such plasmas on sub-femtosecond timescales. We present a synthetic diagnostic method for solid-density laser-generated plasmas based on the dispersion of an extreme-ultraviolet attosecond probe pulse, in a pump--probe scheme. We demonstrate th… ▽ More Extreme-ultraviolet pulses can propagate through ionised solid-density targets, unlike optical pulses and, thus, have the potential to probe the interior of such plasmas on sub-femtosecond timescales. We present a synthetic diagnostic method for solid-density laser-generated plasmas based on the dispersion of an extreme-ultraviolet attosecond probe pulse, in a pump--probe scheme. We demonstrate the theoretical feasibility of this approach through calculating the dispersion of an extreme-ultraviolet probe pulse propagating through a laser-generated plasma. The plasma dynamics is calculated using a particle-in-cell simulation, whereas the dispersion of the probe is calculated with an external pseudo-spectral wave solver, allowing for high accuracy when calculating the dispersion. The application of this method is illustrated on thin-film plastic and aluminium targets irradiated by a high-intensity pump pulse. By comparing the dispersion of the probe pulse at different delays relative to the pump pulse, it is possible to follow the evolution of the plasma as it disintegrates. The high-frequency end of the dispersion provides information on the line-integrated electron density, whereas lower frequencies are more affected by the highest density encountered along the path of the probe. In addition, the presence of thin-film interference could be used to study the evolution of the plasma surface. △ Less

Submitted 26 April, 2022; v1 submitted 1 February, 2022; originally announced February 2022.

Comments: Published in Journal of Plasma Physics

Journal ref: Journal of Plasma Physics, vol. 88, 905880211 (2022)

Effect of Two-Stage Shattered Pellet Injection on Tokamak Disruptions

Authors: O. Vallhagen, I. Pusztai, M. Hoppe, S. L. Newton, T. Fülöp

Abstract: An effective disruption mitigation system in a tokamak reactor should limit the exposure of the wall to localized heat losses and to the impact of high current runaway electron beams, and avoid excessive forces on the structure. We evaluate with respect to these aspects a two-stage deuterium-neon shattered pellet injection in an ITER-like plasma, using simulations with the DREAM framework [M. Hopp… ▽ More An effective disruption mitigation system in a tokamak reactor should limit the exposure of the wall to localized heat losses and to the impact of high current runaway electron beams, and avoid excessive forces on the structure. We evaluate with respect to these aspects a two-stage deuterium-neon shattered pellet injection in an ITER-like plasma, using simulations with the DREAM framework [M. Hoppe et al (2021) Comp. Phys. Commun. 268, 108098]. To minimize the obtained runaway currents an optimal range of injected deuterium quantities is found. This range is sensitive to the opacity of the plasma to Lyman radiation, which affects the ionization degree of deuterium, and thus avalanche runaway generation. The two-stage injection scheme, where dilution cooling is produced by deuterium before a radiative thermal quench caused by neon, reduces both the hot-tail seed and the localized transported heat load on the wall. However, during nuclear operation, additional runaway seed sources from the activated wall and tritium make it difficult to reach tolerably low runaway currents. △ Less

Submitted 26 January, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 17 pages, 7 figures, to be submitted to Nuclear Fusion

Modeling the complete prevention of disruption-generated runaway electron beam formation with a passive 3D coil in SPARC

Authors: RA Tinguely, VA Izzo, DT Garnier, A Sundström, K Särkimäki, O Embréus, T Fülöp, RS Granetz, M Hoppe, I Pusztai, R Sweeney

Abstract: The potential formation of multi-mega-ampere beams of relativistic "runaway" electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetr… ▽ More The potential formation of multi-mega-ampere beams of relativistic "runaway" electrons (REs) during sudden terminations of tokamak plasmas poses a significant challenge to the tokamak's development as a fusion energy source. Here, we use state-of-the-art modeling of disruption magnetohydrodynamics coupled with a self-consistent evolution of RE generation and transport to show that a non-axisymmetric in-vessel coil will passively prevent RE beam formation during disruptions in the SPARC tokamak, a compact, high-field, high-current device capable of achieving a fusion gain Q > 2 in deuterium-tritium plasmas. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Compatibility condition for the Eulerian left Cauchy-Green deformation tensor field

Authors: Tamás Fülöp

Abstract: An explicit compatibility condition formula is presented for the Eulerian left Cauchy-Green deformation tensor field. It is shown to be the appropriate finite-strain counterpart of Saint-Venant's compatibility condition. The difference between the Eulerian problem and the Lagrangian one is highlighted. An explicit compatibility condition formula is presented for the Eulerian left Cauchy-Green deformation tensor field. It is shown to be the appropriate finite-strain counterpart of Saint-Venant's compatibility condition. The difference between the Eulerian problem and the Lagrangian one is highlighted. △ Less

Submitted 11 October, 2022; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: 9 pages

MSC Class: 74A05

The hot-tail runaway seed landscape during the thermal quench in tokamaks

Authors: Ida Svenningsson, Ola Embreus, Mathias Hoppe, Sarah L Newton, Tünde Fülöp

Abstract: Runaway electron populations seeded from the hot-tail generated by the rapid cooling in plasma-terminating disruptions are a serious concern for next-step tokamak devices such as ITER. Here, we present a comprehensive treatment of the thermal quench, including the superthermal electron dynamics, heat and particle transport, atomic physics, and radial losses due to magnetic perturbations: processes… ▽ More Runaway electron populations seeded from the hot-tail generated by the rapid cooling in plasma-terminating disruptions are a serious concern for next-step tokamak devices such as ITER. Here, we present a comprehensive treatment of the thermal quench, including the superthermal electron dynamics, heat and particle transport, atomic physics, and radial losses due to magnetic perturbations: processes that are strongly linked and essential for the evaluation of the runaway seed in disruptions mitigated by material injection. We identify limits on the injected impurity density and magnetic perturbation level for which the runaway seed current is acceptable without excessive thermal energy being lost to the wall via particle impact. The consistent modelling of generation and losses shows that runaway beams tend to form near the edge of the plasma, where they could be deconfined via external perturbations. △ Less

Submitted 3 June, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Lett. 127, 035001 (2021)

DREAM: a fluid-kinetic framework for tokamak disruption runaway electron simulations

Authors: M. Hoppe, O. Embreus, T. Fülöp

Abstract: Avoidance of the harmful effects of runaway electrons (REs) in plasma-terminating disruptions is pivotal in the design of safety systems for magnetic fusion devices. Here, we describe a computationally efficient numerical tool, that allows for self-consistent simulations of plasma cooling and associated RE dynamics during disruptions. It solves flux-surface averaged transport equations for the pla… ▽ More Avoidance of the harmful effects of runaway electrons (REs) in plasma-terminating disruptions is pivotal in the design of safety systems for magnetic fusion devices. Here, we describe a computationally efficient numerical tool, that allows for self-consistent simulations of plasma cooling and associated RE dynamics during disruptions. It solves flux-surface averaged transport equations for the plasma density, temperature and poloidal flux, using a bounce-averaged kinetic equation to self-consistently provide the electron current, heat, density and RE evolution, as well as the electron distribution function. As an example, we consider disruption scenarios with material injection and compare the electron dynamics resolved with different levels of complexity, from fully kinetic to fluid modes. △ Less

Submitted 10 August, 2021; v1 submitted 30 March, 2021; originally announced March 2021.

Comments: 33 pages, 11 figures

Journal ref: Computer Physics Communications 268, November, 108098 (2021)

A case study of non-Fourier heat conduction using Internal Variables and GENERIC

Authors: Mátyás Szücs, Michal Pavelka, Róbert Kovács, Tamás Fülöp, Péter Ván, Miroslav Grmela

Abstract: Applying simultaneously the methodology of Non-Equilibrium Thermodynamics with Internal Variables (NET-IV) and the framework of General Equation for the Non-Equilibrium Reversible-Irreversible Coupling (GENERIC), we demonstrate that, in heat conduction theories, entropy current multipliers can be interpreted as relaxed state variables. Fourier's law and its various extensions -- the Maxwell-Cattan… ▽ More Applying simultaneously the methodology of Non-Equilibrium Thermodynamics with Internal Variables (NET-IV) and the framework of General Equation for the Non-Equilibrium Reversible-Irreversible Coupling (GENERIC), we demonstrate that, in heat conduction theories, entropy current multipliers can be interpreted as relaxed state variables. Fourier's law and its various extensions -- the Maxwell-Cattaneo-Vernotte, Guyer-Krumhansl, Jeffreys type, Ginzburg-Landau (Allen-Cahn) type and ballistic-diffusive -- heat conduction equations are derived in both formulations. Along these lines, a comparison of NET-IV and GENERIC is also performed. Our results may pave the way for microscopic/multiscale understanding of beyond-Fourier heat conduction, and open new ways for numerical simulations of heat-conduction problems. △ Less

Submitted 19 July, 2021; v1 submitted 4 March, 2021; originally announced March 2021.

Comments: 28 pages

Journal ref: Journal of Non-Equilibrium Thermodynamics, 47/1, 31-60, 2022

Size effects and beyond-Fourier heat conduction in room-temperature experiments

Authors: Anna Fehér, Norbert Lukács, László Somlai, Tamás Fodor, Mátyás Szücs, Tamás Fülöp, Péter Ván, Róbert Kovács

Abstract: It is a long-lasting task to understand heat conduction phenomena beyond Fourier. Besides the low-temperature experiments on extremely pure crystals, it has turned out recently that heterogeneous materials with macro-scale size can also show thermal effects that cannot be modelled by the Fourier equation. This is called over-diffusive propagation, different from low-temperature observations, and i… ▽ More It is a long-lasting task to understand heat conduction phenomena beyond Fourier. Besides the low-temperature experiments on extremely pure crystals, it has turned out recently that heterogeneous materials with macro-scale size can also show thermal effects that cannot be modelled by the Fourier equation. This is called over-diffusive propagation, different from low-temperature observations, and is found in numerous samples made from metal foam, rocks, and composites. The measured temperature history is indeed similar to what Fourier's law predicts but the usual evaluation cannot provide reliable thermal parameters. This paper is a report on our experiments on several rock types, each type having multiple samples with different thicknesses. We show that size-dependent thermal behaviour can occur for both Fourier and non-Fourier situations. Moreover, based on the present experimental data, we find an empirical relation between the Fourier and non-Fourier parameters, which may be helpful in later experiments to develop a more robust and reliable evaluation procedure. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Journal ref: Journal of Non-Equilibrium Thermodynamics, 46/4, 403--411, 2021

Modelling of runaway electron dynamics during argon-induced disruptions in ASDEX Upgrade and JET

Authors: K. Insulander Björk, O. Vallhagen, G. Papp, C. Reux, O. Embreus, E. Rachlew, T. Fülöp, the ASDEX Upgrade Team, JET contributors, the EUROfusion MST1 Team

Abstract: Disruptions in tokamak plasmas may lead to the generation of runaway electrons that have the potential to damage plasma-facing components. Improved understanding of the runaway generation process requires interpretative modelling of experiments. In this work we simulate eight discharges in the ASDEX Upgrade and JET tokamaks, where argon gas was injected to trigger the disruption. We use a fluid mo… ▽ More Disruptions in tokamak plasmas may lead to the generation of runaway electrons that have the potential to damage plasma-facing components. Improved understanding of the runaway generation process requires interpretative modelling of experiments. In this work we simulate eight discharges in the ASDEX Upgrade and JET tokamaks, where argon gas was injected to trigger the disruption. We use a fluid modelling framework with the capability to model the generation of runaway electrons through the hot-tail, Dreicer and avalanche mechanisms, as well as runaway electron losses. Using experimentally based initial values of plasma current and electron temperature and density, we can reproduce the plasma current evolution using realistic assumptions about temperature evolution and assimilation of the injected argon in the plasma. The assumptions and results are similar for the modelled discharges in ASDEX Upgrade and JET, indicating that the implemented models are applicable to machines of varying size, which is important for the modelling of future, larger machines. For the modelled discharges in ASDEX Upgrade, where the initial temperature was comparatively high, we had to assume that a large fraction of the hot-tail runaway electrons were lost in order to reproduce the measured current evolution. △ Less

Submitted 30 June, 2021; v1 submitted 6 January, 2021; originally announced January 2021.

Comments: 17 pages, 7 figures

Generation of attosecond electron bunches and X-ray pulses from few-cycle femtosecond laser pulses

Authors: Julien Ferri, Vojtech Horný, Tünde Fülöp

Abstract: Laser-plasma electron accelerators can be used to produce high-intensity X-rays, as electrons accelerated in wakefields emit radiation due to betatron oscillations.Such X-ray sources inherit the features of the electron beam; sub-femtosecond electron bunches produce betatron sources of the same duration, which in turn allow probing matter on ultrashort time scales. In this paper we show, via Parti… ▽ More Laser-plasma electron accelerators can be used to produce high-intensity X-rays, as electrons accelerated in wakefields emit radiation due to betatron oscillations.Such X-ray sources inherit the features of the electron beam; sub-femtosecond electron bunches produce betatron sources of the same duration, which in turn allow probing matter on ultrashort time scales. In this paper we show, via Particle-in-Cell simulations, that attosecond electron bunches can be obtained using low-energy, ultra-short laser beams both in the self-injection and the controlled injection regimes at low plasma densities. However, only in the controlled regime does the electron injection lead to a stable, isolated attosecond electron bunch. Such ultrashort electron bunches are shown to emit attosecond X-ray bursts with high brilliance △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 12 pages, 5 figures

Four spacetime dimensional simulation of rheological waves in solids and the merits of thermodynamics

Authors: Áron Pozsár, Mátyás Szücs, Róbert Kovács, Tamás Fülöp

Abstract: The recent results attained from a thermodynamically conceived numerical scheme applied on wave propagation in viscoelastic/rheological solids are generalized here, both in the sense that the scheme is extended to four spacetime dimensions and in the aspect of the virtues of a thermodynamical approach. Regarding the scheme, the arrangement of which quantity is represented where in discretized spac… ▽ More The recent results attained from a thermodynamically conceived numerical scheme applied on wave propagation in viscoelastic/rheological solids are generalized here, both in the sense that the scheme is extended to four spacetime dimensions and in the aspect of the virtues of a thermodynamical approach. Regarding the scheme, the arrangement of which quantity is represented where in discretized spacetime, including the question of appropriately realizing the boundary conditions, is nontrivial. In parallel, placing the problem in the thermodynamical framework proves to be beneficial in regards to monitoring and controlling numerical artefacts - instability, dissipation error, and dispersion error. This, in addition to the observed preciseness, speed, and resource-friendliness, makes the thermodynamically extended symplectic approach that is presented here advantageous above commercial finite element software solutions. △ Less

Submitted 7 March, 2024; v1 submitted 1 November, 2020; originally announced December 2020.

Comments: 24 pages, 101 figures; v3: 4 typos corrected

Journal ref: Entropy 2020, 22, 1376

Effects of magnetic perturbations and radiation on the runaway avalanche

Authors: P. Svensson, O. Embreus, S. L. Newton, K. Särkimäki, O. Vallhagen, T. Fülöp

Abstract: The electron runaway phenomenon in plasmas depends sensitively on the momentum-space dynamics. However, efficient simulation of the global evolution of systems involving runaway electrons typically requires a reduced fluid description. This is needed for example in the design of essential runaway mitigation methods for tokamaks. In this paper, we present a method to include the effect of momentum-… ▽ More The electron runaway phenomenon in plasmas depends sensitively on the momentum-space dynamics. However, efficient simulation of the global evolution of systems involving runaway electrons typically requires a reduced fluid description. This is needed for example in the design of essential runaway mitigation methods for tokamaks. In this paper, we present a method to include the effect of momentum-dependent spatial transport in the runaway avalanche growth rate. We quantify the reduction of the growth rate in the presence of electron diffusion in stochastic magnetic fields and show that the spatial transport can raise the effective critical electric field. Using a perturbative approach we derive a set of equations that allows treatment of the effect of spatial transport on runaway dynamics in the presence of radial variation in plasma parameters. This is then used to demonstrate the effect of spatial transport in current quench simulations for ITER-like plasmas with massive material injection. We find that in scenarios with sufficiently slow current quench, due to moderate impurity and deuterium injection, the presence of magnetic perturbations reduces the final runaway current considerably. Perturbations localized at the edge are not effective in suppressing the runaways, unless the runaway generation is off-axis, in which case they may lead to formation of strong current sheets at the interface of the confined and perturbed regions. △ Less

Submitted 29 December, 2020; v1 submitted 14 October, 2020; originally announced October 2020.

Comments: 26 pages, 9 figures

Journal ref: J. Plasma Phys. 87 (2021) 905870207

Multimillijoule terahertz radiation from laser interactions with microplasma-waveguides

Authors: Ke Hu, Longqing Yi, Tünde Fülöp

Abstract: When a relativistic, femtosecond laser pulse enters a waveguide, the pulse energy is coupled into waveguide optical modes. The longitudinal laser field effectively accelerates electrons along the axis of the channel, while the asymmetric transverse electromagnetic fields tend to expel fast electrons radially outwards. At the exit of the waveguide, the $\sim$${\rm nC}$, $\sim$$10\ {\rm MeV}$ electr… ▽ More When a relativistic, femtosecond laser pulse enters a waveguide, the pulse energy is coupled into waveguide optical modes. The longitudinal laser field effectively accelerates electrons along the axis of the channel, while the asymmetric transverse electromagnetic fields tend to expel fast electrons radially outwards. At the exit of the waveguide, the $\sim$${\rm nC}$, $\sim$$10\ {\rm MeV}$ electron bunch converts its energy to a $\sim$$10\ {\rm mJ}$ terahertz (THz) laser pulse through coherent diffraction radiation. In this paper, we present 3D particle-in-cell simulations and theoretical analyses of the aforementioned interaction process. We investigate the process of longitudinal acceleration and radial expulsion of fast electrons, as well as the dependence of the properties of the resulting THz radiation on laser and plasma parameters and the effects of the preplasma. The simulation results indicate that the conversion efficiency of energy can be over $5\%$ if the waveguide length is optimal and a high contrast pump laser is used. These results guide the design of more intense and powerful THz sources. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Attosecond betatron radiation pulse train

Authors: Vojtěch Horný, Miroslav Krůs, Tünde Fülop

Abstract: High-intensity X-ray sources are essential diagnostic tools for science, technology and medicine. Such X-ray sources can be produced in laser-plasma accelerators, where electrons emit short-wavelength radiation due to their betatron oscillations in the plasma wake of a laser pulse. Contemporary available betatron radiation X-ray sources can deliver a collimated X-ray pulse of duration on the order… ▽ More High-intensity X-ray sources are essential diagnostic tools for science, technology and medicine. Such X-ray sources can be produced in laser-plasma accelerators, where electrons emit short-wavelength radiation due to their betatron oscillations in the plasma wake of a laser pulse. Contemporary available betatron radiation X-ray sources can deliver a collimated X-ray pulse of duration on the order of several femtoseconds from a source size of the order of several micrometres. In this paper we demonstrate, through particle-in-cell simulations, that the temporal resolution of such a source can be enhanced by an order of magnitude by a spatial modulation of the emitting relativistic electron bunch. The modulation is achieved by the interaction of the that electron bunch with a co-propagating laser beam which results in the generation of a train of equidistant sub-femtosecond X-ray pulses. The distance between the single pulses of a train is tuned by the wavelength of the modulation laser pulse. The modelled experimental setup is achievable with current technologies. Potential applications include stroboscopic sampling of ultrafast fundamental processes. △ Less

Submitted 12 August, 2020; v1 submitted 10 June, 2020; originally announced June 2020.

Comments: 19 pages, 5 figures, revised manuscript submitted to Scientific Reports

Assessing energy dependence of the transport of relativistic electrons in perturbed magnetic fields with orbit-following simulations

Authors: Konsta Särkimäki, Ola Embreus, Eric Nardon, Tünde Fülöp, JET Contributors

Abstract: Experimental observations, as well as theoretical predictions, indicate that the transport of energetic electrons decreases with energy. This reduction in transport is attributed to finite orbit width (FOW) effects. Using orbit-following simulations in perturbed tokamak magnetic fields that have an ideal homogeneous stochastic layer at the edge, we quantify the energy dependence of energetic elect… ▽ More Experimental observations, as well as theoretical predictions, indicate that the transport of energetic electrons decreases with energy. This reduction in transport is attributed to finite orbit width (FOW) effects. Using orbit-following simulations in perturbed tokamak magnetic fields that have an ideal homogeneous stochastic layer at the edge, we quantify the energy dependence of energetic electrons transport and confirm previous theoretical estimates. However, using magnetic configurations characteristic of JET disruptions, we find no reduction in RE transport at higher energies, which we attribute to the mode widths being comparable to the minor radius, making the FOW effects negligible. Instead, the presence of islands and nonuniform magnetic perturbations are found to be more important. The diffusive-advective transport coefficients calculated in this work, based on simulations for electron energies 10 keV -- 100 MeV, can be used in reduced kinetic models to account for the transport due to the magnetic field perturbations. △ Less

Submitted 5 June, 2020; originally announced June 2020.

Spatiotemporal analysis of the runaway distribution function from synchrotron images in an ASDEX Upgrade disruption

Authors: M. Hoppe, L. Hesslow, O. Embreus, L. Unnerfelt, G. Papp, I. Pusztai, T. Fülöp, O. Lexell, T. Lunt, E. Macusova, P. J. McCarthy, G. Pautasso, G. I. Pokol, G. Por, P. Svensson, the ASDEX Upgrade team, the EUROfusion MST1 team

Abstract: Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analyzed using the synchrotron diagnostic tool SOFT and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quenc… ▽ More Synchrotron radiation images from runaway electrons (REs) in an ASDEX Upgrade discharge disrupted by argon injection are analyzed using the synchrotron diagnostic tool SOFT and coupled fluid-kinetic simulations. We show that the evolution of the runaway distribution is well described by an initial hot-tail seed population, which is accelerated to energies between 25-50 MeV during the current quench, together with an avalanche runaway tail which has an exponentially decreasing energy spectrum. We find that, although the avalanche component carries the vast majority of the current, it is the high-energy seed remnant that dominates synchrotron emission. With insights from the fluid-kinetic simulations, an analytic model for the evolution of the runaway seed component is developed and used to reconstruct the radial density profile of the RE beam. The analysis shows that the observed change of the synchrotron pattern from circular to crescent shape is caused by a rapid redistribution of the radial profile of the runaway density. △ Less

Submitted 5 February, 2021; v1 submitted 29 May, 2020; originally announced May 2020.

Comments: 25 pages, 12 figures

Journal ref: Journal of Plasma Physics, Volume 87, Issue 1, February 2021, 855870102

Kinetic modelling of runaway electron generation in argon-induced disruptions in ASDEX Upgrade

Authors: K. Insulander Björk, G. Papp, O. Embreus, L. Hesslow, T. Fülöp, O. Vallhagen, A. Lier, G. Pautasso, A. Bock, the ASDEX Upgrade Team, the EUROfusion MST1 Team

Abstract: Massive material injection has been proposed as a way to mitigate the formation of a beam of relativistic runaway electrons that may result from a disruption in tokamak plasmas. In this paper we analyse runaway generation observed in eleven ASDEX Upgrade discharges where disruption was triggered using massive gas injection. We present numerical simulations in scenarios characteristic of on-axis pl… ▽ More Massive material injection has been proposed as a way to mitigate the formation of a beam of relativistic runaway electrons that may result from a disruption in tokamak plasmas. In this paper we analyse runaway generation observed in eleven ASDEX Upgrade discharges where disruption was triggered using massive gas injection. We present numerical simulations in scenarios characteristic of on-axis plasma conditions, constrained by experimental observations, using a description of the runaway dynamics with self-consistent electric field and temperature evolution in two-dimensional momentum space and zero-dimensional real space. We describe the evolution of the electron distribution function during the disruption, and show that the runaway seed generation is dominated by hot-tail in all of the simulated discharges. We reproduce the observed dependence of the current dissipation rate on the amount of injected argon during the runaway plateau phase. Our simulations also indicate that above a threshold amount of injected argon, the current density after the current quench depends strongly on the argon densities. This trend is not observed in the experiments, which suggests that effects not captured by 0D kinetic modeling -- such as runaway seed transport -- are also important. △ Less

Submitted 27 August, 2020; v1 submitted 20 April, 2020; originally announced May 2020.

Comments: 17 pages, 15 figures, published in Journal of Plasma Physics (Invited Contributions from the 18th European Fusion Theory Conference)

Journal ref: Journal of Plasma Physics (2020) 86(4), 855860401

Runaway dynamics in the DT phase of ITER operations in the presence of massive material injection

Authors: O. Vallhagen, O Embreus, I Pusztai, L Hesslow, T Fülöp

Abstract: A runaway avalanche can result in a conversion of the initial plasma current into a relativistic electron beam in high current tokamak disruptions. We investigate the effect of massive material injection of deuterium-noble gas mixtures on the coupled dynamics of runaway generation, resistive diffusion of the electric field, and temperature evolution during disruptions in the DT phase of ITER opera… ▽ More A runaway avalanche can result in a conversion of the initial plasma current into a relativistic electron beam in high current tokamak disruptions. We investigate the effect of massive material injection of deuterium-noble gas mixtures on the coupled dynamics of runaway generation, resistive diffusion of the electric field, and temperature evolution during disruptions in the DT phase of ITER operations. We explore the dynamics over a wide range of injected concentrations and find substantial runaway currents, unless the current quench time is intolerably long. The reason is that the cooling associated with the injected material leads to high induced electric fields that, in combination with a significant recombination of hydrogen isotopes, leads to a large avalanche generation. Balancing Ohmic heating and radiation losses provides qualitative insights into the dynamics, however, an accurate modeling of the temperature evolution based on energy balance appears crucial for quantitative predictions. △ Less

Submitted 13 July, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: 24 pages, 8 figures

Effects of oblique incidence and colliding pulses on laser-driven proton acceleration from relativistically transparent ultrathin targets

Authors: Julien Ferri, Evangelos Siminos, Laurent Gremillet, Tünde Fülöp

Abstract: The use of ultrathin solid foils offers optimal conditions for accelerating protons from laser-matter interactions. When the target is thin enough that relativistic self-induced transparency (RSIT) sets in, all of the target electrons get heated to high energies by the laser, which maximizes the accelerating electric field and therefore the final ion energy. In this work, we first investigate how… ▽ More The use of ultrathin solid foils offers optimal conditions for accelerating protons from laser-matter interactions. When the target is thin enough that relativistic self-induced transparency (RSIT) sets in, all of the target electrons get heated to high energies by the laser, which maximizes the accelerating electric field and therefore the final ion energy. In this work, we first investigate how ion acceleration by ultraintense femtosecond laser pulses in transparent CH$_2$ solid foils is modified when turning from normal to oblique ($45^\circ$) incidence. Due to stronger electron heating, we find that higher proton energies can be obtained at oblique incidence but in thinner optimum targets. We then show that proton acceleration can be further improved by splitting the laser pulse into two half-pulses focused at opposite incidence angles. An increase by $\sim 30\,\%$ in the maximum proton energy and by a factor of $\sim 4$ in the high-energy proton charge is reported compared to the reference case of a single normally incident pulse. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 11 pages, 7 figures

Runaway electron synchrotron radiation in a vertically translated plasma

Authors: M. Hoppe, G. Papp, T. Wikjamp, A. Perek, J. Decker, B. Duval, O. Embreus, T. Fülöp, U. A. Sheikh

Abstract: Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron di… ▽ More Synchrotron radiation observed from runaway electrons (REs) in tokamaks depends upon the position and size of the RE beam, the RE energy and pitch distributions, as well as the location of the observer. We show that experimental synchrotron images of a vertically moving runaway electron beam sweeping past the detector in the TCV tokamak agree well with predictions from the synthetic synchrotron diagnostic Soft. This experimental validation lends confidence to the theory underlying the synthetic diagnostics which are used for benchmarking theoretical models of and probing runaway dynamics. We present a comparison of synchrotron measurements in TCV with predictions of kinetic theory for runaway dynamics in uniform magnetic fields. We find that to explain the detected synchrotron emission, significant non-collisional pitch angle scattering as well as radial transport of REs would be needed. Such effects could be caused by the presence of magnetic perturbations, which should be further investigated in future TCV experiments. △ Less

Submitted 18 August, 2020; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: 7 pages, 4 figures. Accepted for publication in Nuclear Fusion

Journal ref: Nucl. Fusion 60 094002 (2020)

Evaluation of the Dreicer runaway generation rate in the presence of high-Z impurities using a neural network

Authors: L Hesslow, L Unnerfelt, O Vallhagen, O Embreus, M Hoppe, G Papp, T Fülöp

Abstract: Integrated modelling of electron runaway requires computationally expensive kinetic models that are self-consistently coupled to the evolution of the background plasma parameters. The computational expense can be reduced by using parameterized runaway generation rates rather than solving the full kinetic problem. However, currently available generation rates neglect several important effects; in p… ▽ More Integrated modelling of electron runaway requires computationally expensive kinetic models that are self-consistently coupled to the evolution of the background plasma parameters. The computational expense can be reduced by using parameterized runaway generation rates rather than solving the full kinetic problem. However, currently available generation rates neglect several important effects; in particular, they are not valid in the presence of partially ionized impurities. In this work, we construct a multilayer neural network for the Dreicer runaway generation rate which is trained on data obtained from kinetic simulations performed for a wide range of plasma parameters and impurities. The neural network accurately reproduces the Dreicer runaway generation rate obtained by the kinetic solver. By implementing it in a fluid runaway electron modelling tool, we show that the improved generation rates lead to significant differences in the self-consistent runaway dynamics as compared to the results using the previously available formulas for the runaway generation rate. △ Less

Submitted 24 January, 2020; v1 submitted 1 October, 2019; originally announced October 2019.

Comments: 16 pages, 5 figures

Journal ref: Journal of Plasma Physics 85, 475850601 (2019)

Effect of plasma elongation on current dynamics during tokamak disruptions

Authors: T. Fülöp, P. Helander, O. Vallhagen, O. Embréus, L. Hesslow, P. Svensson, A. J. Creely, N. T. Howard, P. Rodriguez-Fernandez

Abstract: Plasma terminating disruptions in tokamaks may result in relativistic runaway electron beams with potentially serious consequences for future devices with large plasma currents. In this paper we investigate the effect of plasma elongation on the coupled dynamics of runaway generation and resistive diffusion of the electric field. We find that elongated plasmas are less likely to produce large runa… ▽ More Plasma terminating disruptions in tokamaks may result in relativistic runaway electron beams with potentially serious consequences for future devices with large plasma currents. In this paper we investigate the effect of plasma elongation on the coupled dynamics of runaway generation and resistive diffusion of the electric field. We find that elongated plasmas are less likely to produce large runaway currents, partly due to the lower induced electric fields associated with larger plasmas, and partly due to direct shaping effects, which mainly lead to a reduction in the runaway avalanche gain. △ Less

Submitted 6 January, 2020; v1 submitted 30 September, 2019; originally announced September 2019.

Comments: 11 pages, 3 figures

Spin-orbit interaction in a high-power laser irradiated micro-scale plasma waveguide

Authors: Longqing Yi, Ke Hu, Tünde Fülöp

Abstract: Light carries angular momentum as spin and orbital components. The spin-orbit interaction (SOI) of light refers to phenomena in which the spin (left or right circular polarisation) affects the spatial degrees of freedom. Recently, interest in SOI has surged, as it provides physical insight into the behaviour of polarised light at subwavelength scales, and allows for spin-controlled manipulation of… ▽ More Light carries angular momentum as spin and orbital components. The spin-orbit interaction (SOI) of light refers to phenomena in which the spin (left or right circular polarisation) affects the spatial degrees of freedom. Recently, interest in SOI has surged, as it provides physical insight into the behaviour of polarised light at subwavelength scales, and allows for spin-controlled manipulation of light. Most studies are performed with low intensity, leaving the role SOI plays in the relativistic laser-plasma interaction, characterised by nonlinearity, less understood. Here, using 3D particle-in-cell simulations, we report SOI effects in this unprecedented regime. Specifically, a circularly polarised Gaussian laser irradiating a micro-scale plasma waveguide drives a spin-controlled chiral surface wave, allowing the reflected harmonic photons to gain orbital angular momentum (OAM). These effects produce intense optical vortices in the extreme ultraviolet regime - an area of significant fundamental and applied physics potential. △ Less

Submitted 6 September, 2019; originally announced September 2019.

Comments: 7 pages, 3 figures

Thermodynamical extension of a symplectic numerical scheme with half space and time shift demonstrated on rheological waves in solids

Authors: Tamás Fülöp, Róbert Kovács, Mátyás Szücs, Mohammad Fawaier

Abstract: On the example of the Poynting-Thomson-Zener rheological model for solids, which exhibits both dissipation and wave propagation - with nonlinear dispersion relation -, we introduce and investigate a finite difference numerical scheme. Our goal is to demonstrate its properties and to ease the computations in later applications for continuum thermodynamical problems. The key element is the positioni… ▽ More On the example of the Poynting-Thomson-Zener rheological model for solids, which exhibits both dissipation and wave propagation - with nonlinear dispersion relation -, we introduce and investigate a finite difference numerical scheme. Our goal is to demonstrate its properties and to ease the computations in later applications for continuum thermodynamical problems. The key element is the positioning of the discretized quantities with shifts by half space and time steps with respect to each other. The arrangement is chosen according to the spacetime properties of the quantities and of the equations governing them. Numerical stability, dissipative error and dispersive error are analysed in detail. With the best settings found, the scheme is capable of making precise and fast predictions. Finally, the proposed scheme is compared to a commercial finite element software, COMSOL, which demonstrates essential differences even on the simplest - elastic - level of modelling. △ Less

Submitted 7 January, 2020; v1 submitted 21 August, 2019; originally announced August 2019.

Comments: 24 pages, 47 figures

Enhancement of laser-driven ion acceleration in non-periodic nanostructured targets

Authors: I. Thiele, J. Ferri, E. Siminos, L. Gremillet, E. Smetanina, A. Dmitriev, G. Cantono, C. -G. Wahlström, T. Fülöp

Abstract: Using particle-in-cell simulations, we demonstrate an improvement of the target normal sheath acceleration (TNSA) of protons in non-periodically nanostructured targets with micron-scale thickness. Compared to standard flat foils, an increase in the proton cutoff energy by up to a factor of two is observed in foils coated with nanocones or perforated with nanoholes. The latter nano-perforated foils… ▽ More Using particle-in-cell simulations, we demonstrate an improvement of the target normal sheath acceleration (TNSA) of protons in non-periodically nanostructured targets with micron-scale thickness. Compared to standard flat foils, an increase in the proton cutoff energy by up to a factor of two is observed in foils coated with nanocones or perforated with nanoholes. The latter nano-perforated foils yield the highest enhancement, which we show to be robust over a broad range of foil thicknesses and hole diameters. The improvement of TNSA performance results from more efficient hot-electron generation, caused by a more complex laser-electron interaction geometry and increased effective interaction area and duration. We show that TNSA is optimized for a nanohole distribution of relatively low areal density and that is not required to be periodic, thus relaxing the manufacturing constraints. △ Less

Submitted 16 July, 2019; v1 submitted 27 May, 2019; originally announced May 2019.

Comments: 11 pages, 8 figures

Coherent diffraction radiation of relativistic terahertz pulses from a laser-driven micro-plasma-waveguide

Authors: Longqing Yi, Tünde Fülöp

Abstract: We propose a method to generate isolated relativistic terahertz (THz) pulses using a high-power laser irradiating a mirco-plasma-waveguide (MPW). When the laser pulse enters the MPW, high-charge electron bunches are produced and accelerated to ~ 100 MeV by the transverse magnetic modes. A substantial part of the electron energy is transferred to THz emission through coherent diffraction radiation… ▽ More We propose a method to generate isolated relativistic terahertz (THz) pulses using a high-power laser irradiating a mirco-plasma-waveguide (MPW). When the laser pulse enters the MPW, high-charge electron bunches are produced and accelerated to ~ 100 MeV by the transverse magnetic modes. A substantial part of the electron energy is transferred to THz emission through coherent diffraction radiation as the electron bunches exit the MPW. We demonstrate this process with three-dimensional particle-in-cell simulations. The frequency of the radiation is determined by the incident laser duration, and the radiated energy is found to be strongly correlated to the charge of the electron bunches, which can be controlled by the laser intensity and micro-engineering of the MPW target. Our simulations indicate that 100-mJ level relativistic-intense THz pulses with tunable frequency can be generated at existing laser facilities, and the overall efficiency reaches 1%. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 123, 094801 (2019)

Influence of massive material injection on avalanche runaway generation during tokamak disruptions

Authors: L Hesslow, O Embréus, O Vallhagen, T Fülöp

Abstract: In high-current tokamak devices such as ITER, a runaway avalanche can cause a large amplification of a seed electron population. We show that disruption mitigation by impurity injection may significantly increase the runaway avalanche growth rate in such devices. This effect originates from the increased number of target electrons available for the avalanche process in weakly ionized plasmas, whic… ▽ More In high-current tokamak devices such as ITER, a runaway avalanche can cause a large amplification of a seed electron population. We show that disruption mitigation by impurity injection may significantly increase the runaway avalanche growth rate in such devices. This effect originates from the increased number of target electrons available for the avalanche process in weakly ionized plasmas, which is only partially compensated by the increased friction force on fast electrons. We derive an expression for the avalanche growth rate in partially ionized plasmas and investigate the effects of impurity injection on the avalanche multiplication factor and on the final runaway current for ITER-like parameters. For impurity densities relevant for disruption mitigation, the maximum amplification of a runaway seed can be increased by tens of orders of magnitude compared to previous predictions. This motivates careful studies to determine the required densities and impurity species to obtain tolerable current quench parameters, as well as more detailed modeling of the runaway dynamics including transport effects. △ Less

Submitted 10 June, 2019; v1 submitted 1 April, 2019; originally announced April 2019.

Comments: 6 pages, 2 figures

Journal ref: Nuclear Fusion 59, 084004 (2019)

Kluitenberg-Verhás rheology of solids in the GENERIC framework

Authors: Mátyás Szücs, Tamás Fülöp

Abstract: The internal variable methodology of nonequilibrium thermodynamics, with a symmetric tensorial internal variable, provides an important rheological model family for solids, the so-called Kluitenberg-Verhás model family [1]. This model family is distinguished not only from theoretical aspects but also on experimental grounds (see [2] for plastics and [3, 4, 5] for rocks). In this article, we presen… ▽ More The internal variable methodology of nonequilibrium thermodynamics, with a symmetric tensorial internal variable, provides an important rheological model family for solids, the so-called Kluitenberg-Verhás model family [1]. This model family is distinguished not only from theoretical aspects but also on experimental grounds (see [2] for plastics and [3, 4, 5] for rocks). In this article, we present and discuss how the internal variable formulation of the Kluitenberg-Verhás model family can be presented in the nonequilibrium thermodynamical framework GENERIC (General Equation for the Non-Equilibrium Reversible-Irreversible Coupling) [6, 7, 8, 9], for the benefit of both thermodynamical methodologies as well as for promising practical applications. △ Less

Submitted 14 June, 2019; v1 submitted 21 October, 2018; originally announced December 2018.

Comments: 12 pages

Radiation emission in laser-wakefields driven by structured laser pulses with orbital angular momentum

Authors: J. L. Martins, J. Vieira, J. Ferri, T. Fülöp

Abstract: High-intensity X-ray sources are invaluable tools, enabling experiments at the forefront of our understanding of materials science, chemistry, biology, and physics. Laser-plasma electron accelerators are sources of high-intensity X-rays, as electrons accelerated in wakefields emit short-wavelength radiation due to betatron oscillations. While applications such as phase-contrast imaging with these… ▽ More High-intensity X-ray sources are invaluable tools, enabling experiments at the forefront of our understanding of materials science, chemistry, biology, and physics. Laser-plasma electron accelerators are sources of high-intensity X-rays, as electrons accelerated in wakefields emit short-wavelength radiation due to betatron oscillations. While applications such as phase-contrast imaging with these betatron sources have already been demonstrated, others would require higher photon number and would benefit from increased tunability. In this paper we demonstrate, through detailed 3D simulations, a novel configuration for a laser-wakefield betatron source that increases the energy of the X-ray emission and also provides increased flexibility in the tuning of the X-ray photon energy. This is made by combining two Laguerre-Gaussian pulses with non-zero net orbital angular momentum, leading to a rotation of the intensity pattern, and hence, of the driven wakefields. The helical motion driven by the laser rotation is found to dominate the radiation emission, rather than the betatron oscillations. Moreover, the radius of this helical motion can be controlled through the laser spot size and orbital angular momentum indexes, meaning that the radiation can be tuned fully independently of the plasma parameters. △ Less

Submitted 16 October, 2018; originally announced October 2018.

Comments: 12 pages, 5 figures

Analytical solution method for rheological problems of solids

Authors: Tamás Fülöp, Mátyás Szücs

Abstract: In classical continuum theory, Volterra's principle [1, 2] is a long-known method to solve linear rheological (viscoelastic) problems derived from the corresponding elastic ones. Here, we introduce and present another approach that is simpler to apply (no operator inverse is required to compute but only linear ordinary differential equations to solve). Our method starts with the known elastic solu… ▽ More In classical continuum theory, Volterra's principle [1, 2] is a long-known method to solve linear rheological (viscoelastic) problems derived from the corresponding elastic ones. Here, we introduce and present another approach that is simpler to apply (no operator inverse is required to compute but only linear ordinary differential equations to solve). Our method starts with the known elastic solution, replaces the elasticity coefficients with time dependent functions, derives differential equations on them, and determines the solution corresponding to the initial conditions. We present several examples solved via this new method, like tunnels and spherical hollows opened in various initial stress states, and pressurizing of thick-walled tubes and spherical tanks. These examples are useful for applications and, in parallel, are suitable for testing and validating numerical methods of various kinds. △ Less

Submitted 17 September, 2021; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: 49 pages, 23 figures

Proton Acceleration in a Laser-induced Relativistic Electron Vortex

Authors: Longqing Yi, István Pusztai, Alexander Pukhov, Baifei Shen, Tünde Fülöp

Abstract: We show that when a solid plasma foil with a density gradient on the front surface is irradiated by an intense laser pulse at a grazing angle, around 80 degrees, a relativistic electron vortex is excited in the near-critical-density layer after the laser pulse depletion. The vortex structure and dynamics are studied using particle-in-cell simulations. Due to the asymmetry introduced by nonuniform… ▽ More We show that when a solid plasma foil with a density gradient on the front surface is irradiated by an intense laser pulse at a grazing angle, around 80 degrees, a relativistic electron vortex is excited in the near-critical-density layer after the laser pulse depletion. The vortex structure and dynamics are studied using particle-in-cell simulations. Due to the asymmetry introduced by nonuniform background density, the vortex drifts at a constant velocity, typically 0.2 to 0.3 times the speed of light. The strong magnetic field inside the vortex leads to significant charge separation; in the corresponding electric field initially stationary protons can be captured and accelerated to twice the velocity of the vortex (100-200 MeV). A representative scenario - with laser intensity of 10^21 W/cm^2 -is discussed: two dimensional simulations suggest that a quasi-monoenergetic proton beam can be obtained with a mean energy 140 MeV and an energy spread of about 10%. We derive an analytical estimate for the vortex velocity in terms of laser and plasma parameters, demonstrating that the maximum proton energy can be controlled by the incidence angle of the laser and the plasma density gradient. △ Less

Submitted 17 July, 2019; v1 submitted 14 August, 2018; originally announced August 2018.

Comments: 15 pages, 8 figures

Journal ref: Journal of Plasma Physics 85 (2019)

Generalized collision operator for fast electrons interacting with partially ionized impurities

Authors: L Hesslow, O Embréus, M Hoppe, T C DuBois, G Papp, M Rahm, T Fülöp

Abstract: Accurate modelling of the interaction between fast electrons and partially ionized atoms is important for evaluating tokamak disruption mitigation schemes based on material injection. This requires accounting for the effect of screening of the impurity nuclei by the cloud of bound electrons. In this paper, we detail the derivation of a generalized collision operator including the effect of partial… ▽ More Accurate modelling of the interaction between fast electrons and partially ionized atoms is important for evaluating tokamak disruption mitigation schemes based on material injection. This requires accounting for the effect of screening of the impurity nuclei by the cloud of bound electrons. In this paper, we detail the derivation of a generalized collision operator including the effect of partial screening. We calculate the effective ion length-scales, which are needed in the components of the collision operator, for a number of ion species commonly appearing in fusion experiments. We show that for high electric fields, the secondary runaway growth rate can be substantially larger than in a fully ionized plasma with the same effective charge, although the growth rate is significantly reduced at near-critical electric fields. Furthermore, by comparison with the Boltzmann collision operator, we show that the Fokker--Planck formalism is valid even for large impurity content. △ Less

Submitted 27 November, 2018; v1 submitted 13 July, 2018; originally announced July 2018.

Comments: Peer-reviewed version

Journal ref: J. Plasma Phys. (2018), vol . 84, 905840605

Dynamics of positrons during relativistic electron runaway

Authors: Ola Embréus, Linnea Hesslow, Mathias Hoppe, Gergely Papp, Katya Richards, Tünde Fülöp

Abstract: Sufficiently strong electric fields in plasmas can accelerate charged particles to relativistic energies. In this paper we describe the dynamics of positrons accelerated in such electric fields, and calculate the fraction of created positrons that become runaway accelerated, along with the amount of radiation that they emit. We derive an analytical formula that shows the relative importance of the… ▽ More Sufficiently strong electric fields in plasmas can accelerate charged particles to relativistic energies. In this paper we describe the dynamics of positrons accelerated in such electric fields, and calculate the fraction of created positrons that become runaway accelerated, along with the amount of radiation that they emit. We derive an analytical formula that shows the relative importance of the different positron production processes, and show that above a certain threshold electric field the pair production by photons is lower than that by collisions. We furthermore present analytical and numerical solutions to the positron kinetic equation; these are applied to calculate the fraction of positrons that become accelerated or thermalized, which enters into rate equations that describe the evolution of the density of the slow and fast positron populations. Finally, to indicate operational parameters required for positron detection during runaway in tokamak discharges, we give expressions for the parameter dependencies of detected annihilation radiation compared to bremsstrahlung detected at an angle perpendicular to the direction of runaway acceleration. Using the full leading order pair production cross section, we demonstrate that previous related work has overestimated the collisional pair production by at least a factor of four. △ Less

Submitted 12 July, 2018; originally announced July 2018.

Electron beam driven generation of frequency-tunable isolated relativistic sub-cycle pulses

Authors: Illia Thiele, Evangelos Siminos, Tünde Fülöp

Abstract: We propose a novel scheme for frequency-tunable sub-cycle electromagnetic pulse generation. To this end a pump electron beam is injected into an electromagnetic seed pulse as the latter is reflected by a mirror. The electron beam is shown to be able to amplify the field of the seed pulse while upshifting its central frequency and reducing its number of cycles. We demonstrate the amplification by m… ▽ More We propose a novel scheme for frequency-tunable sub-cycle electromagnetic pulse generation. To this end a pump electron beam is injected into an electromagnetic seed pulse as the latter is reflected by a mirror. The electron beam is shown to be able to amplify the field of the seed pulse while upshifting its central frequency and reducing its number of cycles. We demonstrate the amplification by means of 1D and 2D particle-in-cell simulations. In order to explain and optimize the process, a model based on fluid theory is proposed. We estimate that using currently available electron beams and terahertz pulse sources, our scheme is able to produce mJ-strong mid-infrared sub-cycle pulses. △ Less

Submitted 14 February, 2019; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: 8 figures

Journal ref: Phys. Rev. Lett. 122, 104803 (2019)

Enhanced target normal sheath acceleration using colliding laser pulses

Authors: J. Ferri, E. Siminos, T. Fülöp

Abstract: Laser-solid interaction can lead to the acceleration of protons to tens of MeV. Here, we show that a strong enhancement of this acceleration can be achieved by splitting the laser pulse to two parts of equal energy and opposite incidence angles. Through the use of two- and three-dimensional Particle-In-Cell simulations, we find that the multi-pulse interaction leads to a standing wave pattern at t… ▽ More Laser-solid interaction can lead to the acceleration of protons to tens of MeV. Here, we show that a strong enhancement of this acceleration can be achieved by splitting the laser pulse to two parts of equal energy and opposite incidence angles. Through the use of two- and three-dimensional Particle-In-Cell simulations, we find that the multi-pulse interaction leads to a standing wave pattern at the front side of the target, with an enhanced electric field and a substantial modification of the hot electron generation process. This in turn leads to significant improvement of the proton spectra, with an almost doubling of the accelerated proton energy and five-fold enhancement of the number of protons. The proposed scheme is robust with respect to incidence angles for the laser pulses, providing flexibility to the scheme, which should facilitate its experimental implementation. △ Less

Submitted 8 March, 2019; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: 10 pages, 7 figures

Proton acceleration by a pair of successive ultraintense femtosecond laser pulses

Authors: Julien Ferri, Lovisa Senje, Malay Dalui, Kristoffer Svensson, Bastian Aurand, Martin Hansson, Anders Persson, Olle Lundh, Claes-Göran Wahlström, Laurent Gremillet, Evangelos Siminos, Timothy DuBois, Longqing Yi, Joana Martins, Tünde Fülöp

Abstract: We investigate the target normal sheath acceleration of protons in thin aluminum targets irradiated at relativistic intensity by two time-separated ultrashort (35 fs) laser pulses. For identical laser pulses and target thicknesses of 3 and 6 $μ$m, we observe experimentally that the second pulse boosts the maximum energy and charge of the proton beam produced by the first pulse for time delays belo… ▽ More We investigate the target normal sheath acceleration of protons in thin aluminum targets irradiated at relativistic intensity by two time-separated ultrashort (35 fs) laser pulses. For identical laser pulses and target thicknesses of 3 and 6 $μ$m, we observe experimentally that the second pulse boosts the maximum energy and charge of the proton beam produced by the first pulse for time delays below $\sim0.6-1$ ps. By using two-dimensional particle-in-cell simulations we examine the variation of the proton energy spectra with respect to the time-delay between the two pulses. We demonstrate that the expansion of the target front surface caused by the first pulse significantly enhances the hot-electron generation by the second pulse arriving after a few hundreds of fs time delay. This enhancement, however, does not suffice to further accelerate the fastest protons driven by the first pulse once three-dimensional quenching effects have set in. This implies a limit to the maximum time delay that leads to proton energy enhancement, which we theoretically determine. △ Less

Submitted 20 February, 2018; originally announced February 2018.

Comments: 14 pages, 11 figures

Effect of partially ionized impurities and radiation on the effective critical electric field for runaway generation

Authors: Linnea Hesslow, Ola Embréus, George J Wilkie, Gergely Papp, Tünde Fülöp

Abstract: We derive a formula for the effective critical electric field for runaway generation and decay that accounts for the presence of partially ionized impurities in combination with synchrotron and bremsstrahlung radiation losses. We show that the effective critical field is drastically larger than the classical Connor-Hastie field, and even exceeds the value obtained by replacing the free electron de… ▽ More We derive a formula for the effective critical electric field for runaway generation and decay that accounts for the presence of partially ionized impurities in combination with synchrotron and bremsstrahlung radiation losses. We show that the effective critical field is drastically larger than the classical Connor-Hastie field, and even exceeds the value obtained by replacing the free electron density by the total electron density (including both free and bound electrons). Using a kinetic equation solver with an inductive electric field, we show that the runaway current decay after an impurity injection is expected to be linear in time and proportional to the effective critical electric field in highly inductive tokamak devices. This is relevant for the efficacy of mitigation strategies for runaway electrons since it reduces the required amount of injected impurities to achieve a certain current decay rate. △ Less

Submitted 2 May, 2018; v1 submitted 2 February, 2018; originally announced February 2018.

Comments: 17 pages, 6 figures. Minor revisions of original manuscript

Journal ref: Plasma Physics and Controlled Fusion 60, 074010 (2018)

Interpretation of runaway electron synchrotron and bremsstrahlung images

Authors: M. Hoppe, O. Embréus, C. Paz-Soldan, R. A. Moyer, T. Fülöp

Abstract: The crescent spot shape observed in DIII-D runaway electron synchrotron radiation images is shown to result from the high degree of anisotropy in the emitted radiation, the finite spectral range of the camera and the distribution of runaways. The finite spectral camera range is found to be particularly important, as the radiation from the high-field side can be stronger by a factor $10^6$ than the… ▽ More The crescent spot shape observed in DIII-D runaway electron synchrotron radiation images is shown to result from the high degree of anisotropy in the emitted radiation, the finite spectral range of the camera and the distribution of runaways. The finite spectral camera range is found to be particularly important, as the radiation from the high-field side can be stronger by a factor $10^6$ than the radiation from the low-field side in DIII-D. By combining a kinetic model of the runaway dynamics with a synthetic synchrotron diagnostic we see that physical processes not described by the kinetic model (such as radial transport) are likely to be limiting the energy of the runaways. We show that a population of runaways with lower dominant energies and larger pitch-angles than those predicted by the kinetic model provide a better match to the synchrotron measurements. Using a new synthetic bremsstrahlung diagnostic we also simulate the view of the Gamma Ray Imager (GRI) diagnostic used at DIII-D to resolve the spatial distribution of runaway-generated bremsstrahlung. △ Less

Submitted 12 December, 2017; originally announced December 2017.

Comments: 21 pages, 11 figures

Origins of plateau formation in ion energy spectra under target normal sheath acceleration

Authors: Timothy C. DuBois, Evangelos Siminos, Julien Ferri, Laurent Gremillet, Tünde Fülöp

Abstract: Target normal sheath acceleration (TNSA) is a method employed in laser--matter interaction experiments to accelerate light ions (usually protons). Laser setups with durations of a few 10 fs and relatively low intensity contrasts observe plateau regions in their ion energy spectra when shooting on thin foil targets with thicknesses of order 10 $\mathrmμ$m. In this paper we identify a mechanism whic… ▽ More Target normal sheath acceleration (TNSA) is a method employed in laser--matter interaction experiments to accelerate light ions (usually protons). Laser setups with durations of a few 10 fs and relatively low intensity contrasts observe plateau regions in their ion energy spectra when shooting on thin foil targets with thicknesses of order 10 $\mathrmμ$m. In this paper we identify a mechanism which explains this phenomenon using one dimensional particle-in-cell simulations. Fast electrons generated from the laser interaction recirculate back and forth through the target, giving rise to time-oscillating charge and current densities at the target backside. Periodic decreases in the electron density lead to transient disruptions of the TNSA sheath field: peaks in the ion spectra form as a result, which are then spread in energy from a modified potential driven by further electron recirculation. The ratio between the laser pulse duration and the recirculation period (dependent on the target thickness, including the portion of the pre-plasma which is denser than the critical density) determines if a plateau forms in the energy spectra. △ Less

Submitted 28 November, 2017; v1 submitted 10 October, 2017; originally announced October 2017.

Comments: 11 pages, 12 figures