Charged particle collisionless transport near the X-point of the two-wire model

Collisionless charged particle motion and its transport in the two-wire model (TWM) with no axial magnetic fields is investigated numerically. The TWM configuration contains a magnetic X-point, and single particle motions in such a field have two conserved quantities: the total kinetic energy and the base field line value which is a quantity derived from the axial canonical momentum. As gyrating particles travel along the field lines, they may reach near the X-point region where the magnetic moment, the first adiabatic invariant, can be occasionally shifted due to a large gradient of the field. When the magnetic moment becomes large, resulting in a large Larmor radius, particles probabilistically cross the X-point to migrate to the opposite side of the TWM configuration. These phenomena are investigated with single particle simulations. We find that the statistical behaviour of the seemingly chaotic magnetic moment shifts are completely determined by the two aforementioned conserved quantities, and also that there exists a threshold energy, determined by the base field line value, allowing only particles with a higher energy to cross the separatrix and migrate. It is found that the crossing time is distributed exponentially, and that the migration confinement time, which is the average crossing time, is shorter for particles with a base field line closer to the separatrix and a higher energy. We provide an empirical expression, derived with the simulations, for estimating the collisionless migration confinement time.


Introduction
Cross-field (across the background magnetic field) transport of charged particles is gaining attention in low temperature plasmas (Curreli & Chen 2014) in addition to its strong relevance in magnetically confined high temperature plasmas.This is due to, for instance, its relevance to linear plasma generators (Owen et al. 2017), industry applications using magnetrons or helicons (Curreli & Chen 2014) and Langmuir probe diagnostics of magnetized plasmas (Popov et al. 2016).We investigate such cross-field transport in a curved magnetic field geometry including a magnetic X-point.The two-wire model (TWM) is the magnetic configuration generated by the axially flowing two parallel currents which form a magnetic X-point at the centre.The TWM is of research interest as it resembles the magnetic configurations of single-null diverted tokamaks and those of the magnetic reconnections (Auerbach & Boozer 1980).Understanding plasma transport in the TWM configuration, especially in the vicinity of the X-point, can help us to comprehend how plasmas are lost to the divertors as well as their spatial profiles in the scrape-off-layer region in tokamaks, which are active research topics for controlling the heat and particle fluxes associated with fusion-grade hot plasmas (Pitcher & Stangeby 1997; ITER Divertor Physics Expert Group on Divertor and ITER Physics Expert Group on Divertor Modelling and Database and ITER Physics Basis 1999; Loarte et al. 2007;Soukhanovskii 2017;Innocente et al. 2022).Furthermore, a collisionless transport phenomenon with the finite Larmor-radius (FLR) effect is deemed significant in a magnetic island with the gyrokinetic approach (Choi & Hahm 2022;Choi 2023).
Numerous studies about the TWM have been performed in the past.Reiman (1996) identified a set of singular surfaces in the TWM open field line region, which plays a role of rational surfaces and separatrices in toroidal plasmas.Ali, Punjabi & Boozer (2009) studied the magnetic field structure in the vicinity of the separatrix of the TWM.Auerbach & Boozer (1980) studied classical diffusion in the TWM.They superposed axial (toroidal) magnetic fields, making the Larmor radii small enough even at the X-point such that the FLR effects could be neglected.
We are, on the other hand, interested in properties of plasmas in the basic TWM, i.e. no superposed axial magnetic fields, as it intensifies the effects of spatially varying magnetic fields and the FLR on plasma behaviours near the X-point region where the total field strength approaches zero.Only two circular fields produced by the two current carrying wires are combined to make a figure-eight shaped two-dimensional TWM field in this work.Recently, a low temperature plasma device was developed to experimentally investigate plasmas in such fields (Lim et al. 2020), where interesting plasma profiles were reported around the X-point.
Owing to its large spatial gradient of the magnetic field and the FLR near the X-point, the study of plasma properties in the TWM configuration with kinetic theory or fluid models can be expensive; thus, we numerically investigate a single particle motion which allows us to understand plasma transport in the collisionless limit and to gain a physical intuition of such plasmas.In fact, there have been several related studies with single particle motions.Sonnerup (1971) studied adiabatic particle orbits in a magnetic null sheet and conserved quantities of interest.Kim & Cary (1983) analysed charged particle motions near a linear magnetic null.Numata & Yoshida (2002, 2003) studied collisionless resistivity in the magnetic null region by analysing chaotic orbits of particles.Chirikov (1960) reported resonant effects between the gyro-frequency and the bounce frequency in a magnetic trap that can cause variations of the magnetic moment and subsequent particle loss from a magnetic trap.Büchner & Zelenyi (1989) investigated non-adiabatic trapped particle motions in magnetic field reversals and showed how the curvature of the fields determines different types of particle trajectories.
The purpose of this paper is to numerically investigate collisionless transport of charged particles near the X-point in the TWM configuration based on single particle motion.In this study, the term collisionless specifically means single particle motion without considering the influence of electric and magnetic fields induced by other charged https://doi.org/10.1017/S0022377824000321Published online by Cambridge University Press particles.Gyrating particles are guided along the corresponding base field line (invariant), which is determined by the axial canonical momentum, one of the two conserved quantities in the TWM.The other conserved quantity is the total kinetic energy.As the particles travel along the base field line, they can enter a region where a spatial gradient of the field is large, i.e. near the X-point, and it causes the first adiabatic invariant, the magnetic moment μ = W ⊥ /B, to change (Stephens, Brzozowski & Jenko 2017), which, then, can lead to collisionless transport.Here, W ⊥ and B are the perpendicular kinetic energy and the magnetic field magnitude, respectively.Specifically, when the magnetic moment becomes large, resulting in a large Larmor radius, particles probabilistically cross the separatrix to migrate to the opposite side of the TWM configuration.
We find that a seemingly random shift of the magnetic moment near the X-point is, in fact, statistically well behaved and with a property that is completely determined by the two conserved quantities, i.e. the total kinetic energy and the axial canonical momentum or the base field line value.We also find that there exists a threshold energy allowing only particles with a higher energy to cross the separatrix and migrate, where the threshold energy is determined by the base field line value.Note that a migrating particle experiences the reversal of the magnetic field direction, resulting in a betatron orbit, similar to orbits found near null regions of field-reversed-configuration plasmas (Welch et al. 2010).The crossing time is found to be exponentially distributed and its average, which we denote as the migration confinement time, is shorter for particles with a base field line closer to the separatrix and a higher energy.The migration confinement time is heuristically formulated with the simulation results so that time scales of the collisionless transport can be readily estimated.
The outline of this paper is as follows.We first introduce the formula, geometry and characteristics of the TWM in detail in § 2.Then, with the equation of motion, invariants that govern the particle motion are explained in § 3 as well as the simulation method and its conditions.With an example, two phenomena of interests, i.e. magnetic moment shift and the migration, are also discussed in this section.In § 4, the statistical property of magnetic moment shift is analysed.In § 5, the migration phenomenon is discussed.The necessary condition for the migration to occur is derived, and the migration confinement times are empirically formulated.Then, the summary and conclusion are given in § 6.

Two-wire model
To investigate the single charged particle motion, we first introduce the TWM configuration.The TWM is the one of the simplest analytic geometries that contains the magnetic X-point.In the TWM, we place two axial current carrying wires at the positions of x = 0 and y = ± .The current I is applied to both wires.Then, the field of the TWM can be expressed as where A 0 = Iμ 0 /2π is the system magnetic flux, μ 0 is the vacuum permeability, y ± = y ± , r 2 = x 2 + y 2 and r 2 ± = x 2 + y 2 ± .Normalization is performed to make all quantities dimensionless, so that the results of this investigation can be applied to systems with any scale.The normalized TWM field can be written as  where B 0 = A 0 / is the system field magnitude, x = x/ , ỹ = y/ , ỹ± = y/ ± 1, r2 = x2 + ỹ2 and r2 ± = x2 + ỹ2 ± .In this work, the subscript '0' is used to describe system characteristic quantities, whereas the tilde hat, i.e. 'x', is used to denote normalized dimensionless quantities.
Realizing that we can express the magnetic field as B = ẑ × ∇ψ due to there being no axial magnetic fields, a field line (flux surface) of the TWM can be formulated as where a constant s or ψ draws a field line.Notice that (2.4) is the equation for the Cassini oval (see figure 1).The separatrix is at s = 1, or equivalently ψ = 0. Finally, the normalized magnitude and gradient scale length of the magnetic field in the TWM configuration are, respectively, Figure 1(a) shows normalized magnitudes of the magnetic field, which becomes infinitely strong at the wires and becomes null at the X-point.
System quantities used to normalize various physical parameters in the TWM.
X-point, the magnetic moment μ cannot be conserved properly (Stephens et al. 2017).This phenomenon and its consequences are discussed in detail with numerical results in § 4 and § 5.

Single particle motion
Single charged particle motion in the TWM configuration is studied numerically at the collisionless limit.At this limit, interactions among particles through induced electric and magnetic fields are neglected.Therefore, only the externally applied TWM field governs the motion of the particles, and the equation of motion is given as where ω c0 = qB 0 /m is the system cyclotron frequency, q and m are, respectively, the charge and mass of the particle of interest.In addition, τ 0 = 2π/|ω c0 | is the system gyro-period, that normalizes temporal quantities in this study; v 0 = qA 0 /m = ω c0 and W 0 = mv 2 0 /2 are the system speed and energy, which normalize speed and energy quantities, respectively.As there are various normalizing quantities in this work, they are summarized in table 1.

Constants of motion: the total kinetic energy and the base field line value
There are two constants of motion in this system: the total kinetic energy W total , and the base field line value ψ base .Since the magnetic field alone does no work, the total energy, W total = mv 2 /2, is conserved.The Lagrangian L and the axial canonical momentum p z,can of the system are By the Euler-Lagrange equation, i.e. dp z,can /dt = ∂L/∂z = 0, we see that p z,can is a constant of motion in this system (Kim & Cary 1983).We define a base field line value ψ base using p z,can such that ψ base = p z,can /(−q) = ψ(x, y) − mv z /q; thus, ψ base is a conserved quantity in this system as well.The two invariant quantities, W total and ψ base , are normalized as The base field line can be treated as the guiding field line that a particle follows and gyrates about.So, ψbase determines the field line that the particle follows, and Wtotal governs the travel speed and the Larmor radius of the particle.Hence, these two invariants are used to distinguish particles in the TWM system.

An example of numerical simulation
To simulate a particle following (3.1), the Boris method is implemented for its algorithmic efficiency (Qin et al. 2013;Wei et al. 2015).Unless stated otherwise, simulations in this work have the TWM parameters of I = 1 kA and = 0.1 m, resulting in A 0 = 200 μTm and B 0 = 2000 μT.Particles are electrons, which give q/m = −1.76× 10 11 C kg −1 , τ 0 = 17.86 ns and W 0 = 3517.64eV.Note that our results can be readily applied to any charged particles such as ions since we use normalized quantities.The time step of the simulation is set to t = 1 ps, which is much smaller than the system gyro-period τ 0 , in order to minimize numerical errors.
Simulation results are shown in figure 2 as an example.For this case, we have simulated two single particles that follow the base field line of ψbase = −0.05(the magenta field line) with a total energy of Wtotal = 0.0035, whose corresponding dimensional values are ψ base = −10 μTm and W total = 12.31 eV.The two particles start their motions from the top of ψ = −0.05 with tan −1 (v y /v x ) = 45 • (red) and 46 • (green), respectively.Figure 2(a,b) shows the traces of the particle motions in the time range of t = 0−250 ns and t = 250−500 ns, respectively.Figure 2(c) shows temporal evolutions of various parameters which are, from the top, ỹ, ψ, ψbase , Wtotal , W ⊥ /W total , B and μ, where the red and green lines correspond to the red and green particles from figure 2(a,b), respectively.Here, the normalized magnetic moment μ is defined as μ = (W ⊥ /W total )/ B.
There are two interesting phenomena to observe in this example, which are the main topics of this study, i.e. the magnetic moment shifts and the associated migrations.First, as can be seen from the last plot of figure 2(c), the particles' normalized magnetic moments μ make seemingly random shifts to various values.The μ shifts occur whenever the particle travels near the X-point, where the field gradient is large.This observation including its sensitivity to the initial injection angle, i.e. difference in green and red particles, is further detailed in § 4.
Second, a particle may jump to follow other closed field lines.For ψbase < 0, the base field line is two closed circles separated by the separatrix, as indicated by the magenta lines in figure 2(a,b).A particle typically gyrates about and travels along the closed field line on one side.However, if the Larmor radius becomes large enough, due to occasional magnetic moment shifts, the particle can transport to the other closed field line in the TWM configuration, as can be seen with the red particle in figure 2(a).During this transport, due to the reversal of the magnetic field direction, the particle follows one cycle of the betatron orbit.This transport phenomenon is termed 'migration' in this study.The first temporal plot of figure 2(c) shows that the red particle migrates at t ≈ 100 ns, that is, ỹ becomes negative from positive, and stays on the negative side.Note that the green particle does not migrate, at least up to t = 1000 ns in this example; however, it can eventually migrate to the other side as well for it satisfies the necessary condition for the migration, which is discussed in § 5.

Magnetic moment shift
We have observed that the behaviours of the μ shift are different for different initial injection angles with the same two invariants (see the μ traces in figure 2c).This can be explained with an inhomogeneous magnetic field causing chaotic particle orbits (Büchner & Zelenyi 1989;Numata & Yoshida 2002, 2003), which result in different realizations of μ.Due to its chaotic behaviour, the statistical properties of the μ shift are investigated.

Histogram of the magnetic moment
We first consider the overall consequences of the initial injection angles for the μ shift.This is done with seven particles whose invariants are the same, i.e. ψbase = −0.05 and Wtotal = 0.002.All of them start from the top of ψ = −0.05,but with different initial injection angles from 0 • to 30 • with an interval of 5 • .Figure 3(a) shows histograms of μ for the seven particles obtained by simulating them up to t = 10 ms.It clearly shows that the different initial injection angles all lead to the same histogram, while their individual realizations are different, as shown in figure 3(b).
Although not shown here, many more cases with random positions and injection angles (0 to 2π) for the same invariant pair are simulated, and the resultant histograms are found to be identical.The same holds true for cases with different invariant pairs, and the invariant pair is the only factor that determines the shape of the histogram.Therefore, we can infer that the path, which a particle covers in a long time, is independent of the initial condition and distinguished by the invariant pair of the particle.

Correlation between the current and next magnetic moments
Next, we statistically investigate how the current value of μ is correlated with the next one.Figure 4(a) shows the trajectory (red) of a single particle having the invariant values of ψbase = −0.05(magenta) and Wtotal = 0.00225 with the initial injection angle of tan −1 (v y /v x ) = 45 • .We find that μ shifts abruptly whenever the particle passes near the large gradient region around the X-point, indicated by the blue line in figure 4(a), which is the reference line connecting the positions of the two wires (cf.figure 1c).It can be seen clearly from figure 4(b) that the abrupt shifts of μ coincide with the blue vertical lines that indicate the times when the particle passes the reference line.Thus, we define an expectation value of the μ, denoted as μ k , to be the arithmetic average of μ values in between two consecutive (kth and k + 1th) crossings of the reference line.
This single particle is simulated up to 30 ms to obtain the correlation between μ k and μ k+1 .This is visualized as a two-dimensional normalized histogram with 400 × 400 square bins in figure 4(c).As attested by the histogram, the behaviour of the μ shifts is not completely random, indicating that the next μ can be probabilistically predicted given a current μ , albeit with stochasticity.In addition, multiple simulations with the same pair of invariants are performed to conclude that the resultant histograms are independent of the initial injection positions or angles.
We have also investigated how different sets of the invariants affect the correlation between the current and next μ .For this purpose, we introduce a new quantity, the 'base energy', defined as Wbase = ψ2 base .Basically, Wbase is the normalized threshold energy for a particle to cross the separatrix, i.e.only if Wtotal / Wbase is greater than unity can a particle cross the separatrix.Its role as the threshold and the physical meaning of Wbase are derived and discussed in § 5. From here on, we use ψbase and Wtotal / Wbase as a pair of the two invariants to distinguish a particle.
Figure 5 shows the two-dimensional histograms illustrating correlations between the current and next μ for different sets of the invariants.Figure 5(a-d) is of ψ = −0.05,while figure 5(e-h) is of ψ = −0.20 with different values of Wtotal / Wbase .The particles with a lower energy ( Wtotal / Wbase = 0.1) have the sharper histograms (figure 5a,e), meaning that the amounts of μ shifts are relatively small and the shifts are more predictable.This is because the particle travels more slowly, and the magnetic field change it feels is slower.The particles with Wtotal / Wbase = 1.0 move faster and can reach the X-point.Hence, they can feel more rapid changes of the magnetic field, and so the histograms are wider (figure 5c,g).Therefore, we conclude that, in the TWM configuration, the magnetic moment shift is more stochastic for a higher energy since such particles travel faster and go closer to the large gradient region.

Migration
As noted earlier, a normalized field line value ψ less than zero corresponds to two separate closed field lines located above and below the magnetic X-point inside the separatrix (see magenta lines in figure 2a,b).Migration is defined, in this work, as the phenomenon where a particle jumps to follow the other field line, i.e. the sign of the ỹ-position of a particle flips, with the same ψ in the TWM configuration.

Threshold energy for the migration
We find that probability of migration depends on the two invariants, and there is a minimum required energy.These findings can be well explained with a possible particle trajectory described by the effective potential (Kim & Cary 1983).Since there is no z−component of the magnetic field in the TWM configuration, we divide the equation of motion (3.1) into xy-and z-components to derive the effective potential as Then, taking the xy-component, we have where (5.3)Note that the 'base energy' was defined earlier (see § 4) as Wbase = ψ2 base = W base /W 0 .The total energy of a particle is the sum of the kinetic energy in the xy-plane and the effective potential, i.e.W total = W xy (v x , v y ) + U eff (x, y), which means that U eff cannot exceed W total .Therefore, a particle's trajectory is restricted to be within a certain range of ψbase , where its minimum and maximum bounds, denoted as ψmin and ψmax , are  This tells us that a particle with ψbase < 0 can have a trajectory outside of or on the separatrix only if Wtotal / Wbase ≥ 1.Only in this case may the particle cross the separatrix and the X-point.Therefore, this is precisely the necessary condition for migration; in other words, Wbase is the threshold energy.
Figure 6(a-d) shows contour plots of U eff /W total for different energies Wtotal / Wbase with the same ψbase = −0.10(magenta lines).This shows the two bounds, i.e. ψmin and ψmax , for the particle motion.It is impossible for a particle to travel into the region of U eff /W total > 1, depicted as the grey area.It clearly shows that the accessible area is divided into upper and lower sides, and the two sides are connected when Wtotal / Wbase ≥ 1.This means that migration is only possible when this condition is satisfied, which is the same conclusion discussed above.
A particle with Wtotal / Wbase ≥ 1 travels on one side of the TWM configuration, experiences the occasional magnetic moment shifts (discussed in § 4) and, eventually, it can have a Larmor radius large enough to cross the ỹ = 0 axis to migrate to the other side.This is how the migration occurs, and it is significant for it represents a possibility of collisionless transport loss of particles that reside inside the separatrix in addition to a turbulence driven collisionless transport (Choi & Hahm 2022;Choi 2023).

Migration confinement time (collisionless transport time scale)
We have recorded the times a particle spent in between the consecutive migrations, and histograms of such times are shown in figure 7 for [ ψbase , Wtotal / Wbase ] equal to (a) [−0.10, 1.20] and (b) [−0.10, 1.50].As attested by the figure, the times are distributed exponentially, indicating that migration is an almost memoryless phenomenon.There exists a sharp outlier in a very short time, which may be caused by the gyrating orbits of a particle very close to the X-point.The average time is termed as the 'migration confinement time', τ mig , and we have τ mig = 594 ns and 275 ns for figure 7(a,b), respectively.This quantity represents how probable migration is for a particle with a certain pair of the invariants.We see that the migration confinement time is shorter for a particle with a higher energy.
We have seen that collisionless transport is possible due to the magnetic moment shift as a particle with enough energy travels into the region of a large magnetic field gradient, such as around the magnetic X-point.Thus, we have estimated the normalized migration confinement time τmig = τ mig /τ 0 for 600 different pairs of the invariants such that the overall quantitative behaviour of τ mig can be empirically formulated.We have simulated a particle with six values of ψbase , which are −0.005,−0.010, −0.020, −0.050, −0.100 and −0.200.For each ψbase , we have 100 different values of Wtotal from Wtotal − Wbase equal to 0.01 Wbase to unity.
Figure 8 shows the estimated τmig as a function of (a) Wtotal − Wbase , i.e. how much excess energy a particle has, and (b) Wtotal / Wbase − 1, i.e. how much fractional excess energy with respect to the base energy Wbase (or the threshold energy) a particle has.We now construct an empirical expression for τmig based on the numerical results.As can be seen from figure 8, the trend follows two different power laws.First of all, there seems to be, from figure 8

Application to tokamak plasmas
This migration phenomenon causes particles to become completely lost from the system if there is a loss boundary on the other side of the TWM configuration, such as a divertor in a tokamak.To investigate how relevant this collisionless transport can be, we create a similar environment, although it is crude for there are no axial (or toroidal) magnetic fields, as in a tokamak with the independent system quantities of I = 100 kA and = 1.0 m with a deuteron species (q/m = 4.79 × 10 7 C kg −1 ).We then have τ 0 = 6555.33ns and W 0 = 9584.85eV and, therefore, t = 1 ns is used for this simulation.We randomly place 30 000 deuterons on the upper side of ψbase = −0.1 with isotropically distributed initial injection angles.These ions are initialized with a Maxwellian distribution having T = 2/3 W total = 100 eV (see figure 10d).Since all the ions have ψbase = −0.1,we have W base = 95.85 eV, which means only some of the deuterons with a high enough energy may migrate to the lower side and eventually be lost to the loss boundary (virtual divertor) located at ỹ = −0.5.
Figure 10(a-c) shows the spatial distribution of deuterons (green dots) at times t = 0.0, 0.2 and 0.5 ms, respectively.Figure 10(d-f ) shows the W total histograms at these times.It is clear that, as time progresses, the high tail of the distribution is quickly lost to the virtual divertor since the migration confinement times of the high energy particles are shorter.This can be also seen with the decreasing average energy of the particles, denoted as W total .By the time t = 0.5 ms, most of the deuterons with W total > W base have migrated to the lower side, and only the low energy deuterons remain confined on the upper side.Figure 10(g,h) shows temporal evolutions of the number of particles and W total , respectively, for the upper (red) and lower (blue) sides of the TWM configuration.
We note that, even at around t = 2.0 ms, there exist particles migrating to the lower side due to the random nature of the phenomenon.Through this example we can see that, even if particles are following the base field line inside the separatrix, they may migrate to the other side and become lost to a boundary.
As mentioned earlier, this result may not be directly applied to tokamak edge physics as we lack not only the axial (toroidal) magnetic fields but also collisions, which restore a Maxwellian distribution.Nevertheless, as a typical transport time and an ion collision time in a tokamak are of the order of a few to tens of ms (Wesson 2011), we see that the collisionless transport time scale associated with the migration is much faster.If there were axial magnetic fields, they would have slowed down the migration transport, while the collisions would have escalated it.Thus, the migration may play a non-negligible role for determining an edge transport time scale as well as a radial profile in the scrape-off-layer region of a tokamak.

Summary and conclusion
We have numerically investigated a collisionless single charged particle motion in the TWM configuration whose magnetic geometry follows a Cassini oval with a magnetic X-point.We have focused on the magnetic moment shift and its subsequent transport as a particle travels near the X-point, where the conservation of the magnetic moment is no longer valid, i.e. a region with large spatial gradients of the magnetic field magnitude.
From the Lagrangian analysis, two system invariants which can be used to distinguish particles are identified, i.e. the total kinetic energy and the base field line value.We find that, although the particles' initial conditions are different, causing various realizations of chaotic magnetic moment shifts, their long time trajectories eventually cover the same area as long as their pairs of the invariants are the same.Detailed numerical analyses show that seemingly random magnetic moment shifts are, in fact, statistically well behaved, that is, knowledge of the current magnetic moment allows us to predict the next one with a certain probability distribution, which is completely determined by the two invariants.
We have observed that, when the magnetic moment shifts to a larger value, leading to a larger Larmor radius, the particle probabilistically crosses the separatrix to follow the base field line on the other side of the TWM configuration.This phenomenon is denoted as migration in this work.When a migration occurs, the particle can no longer be considered confined inside the separatrix, resulting in a collisionless loss of the confined particles.Using the two-dimensional effective potential, we have found a necessary condition for migration that the total energy of a particle must be larger than a threshold energy determined by one of the invariants, i.e. the base field line value.
As evidenced by the exponentially distributed separatrix-crossing time, migration is a memoryless random process.Therefore, we have estimated the average crossing time, the migration confinement time, and it is found to be only a function of the two invariants.Since the migration confinement time is relevant to the collisionless transport, we have formulated its empirical expression.This has been done by simulating particles with 600 different pairs of the invariants from which we can estimate the migration confinement time.
To examine whether or not the migration phenomenon investigated in this work is relevant for the edge transport of tokamak plasmas, we have simulated particles with similar environments to that in a typical tokamak.Albeit our results may not be directly used to explain the edge transport due to rather crude simplifications, we have found that the migration has a certain degree of potential for determining an edge transport time scale and a radial profile in the scrape-off-layer region of a tokamak.

FIGURE 1 .
FIGURE 1. Generated magnetic field with the TWM where the single particle motion is investigated in this work.(a) Normalized magnitudes of the magnetic field B with overlaid ψ drawn in black lines from −1 to 1 (−1, −0.5, −0.1, −0.05, 0, 0.05, 0.1, 0.5, 1).The field lines inside the figure-eight shaped separatrix have negative ψ.(b) Normalized field line values ψ.It is zero, i.e. white colour, at the separatrix.(c) Normalized gradient scale lengths of the magnetic field |∇ ln B| with the overlaid ψ drawn in black lines whose values are the same as in (a).
Figure1(a) shows normalized magnitudes of the magnetic field, which becomes infinitely strong at the wires and becomes null at the X-point.Figure1(b)shows the corresponding normalized field line values ψ, which are negative (positive) inside (outside) the separatrix where the field line value is zero.Figure1(c)shows normalized gradient scale lengths of the magnetic field, and the gradient becomes infinitely large at the X-point.Since the magnetic field magnitude changes fast in space around the

FIGURE 2 .
FIGURE 2. Two particles with the same invariants ( ψbase = −0.05,Wtotal = 0.0035) are simulated in the TWM configuration.These particles start from the top of ψ = −0.05 with tan −1 (v y /v x ) = 45 • (red) and 46 • (green).The magenta field line indicates the location of ψbase = −0.05.Traces of the particles are shown in (a) for t = 0-250 ns and in (b) for t = 250-500 ns.Temporal traces of various physical quantities of interest for the red and the green particles are shown in (c) from t = 0 to 1000 ns as red and green lines, respectively.
FIGURE 3. (a) Histograms of μ for seven particles obtained by simulating them up to t = 10 ms.Invariant pairs of the seven particles are identical, i.e. ψbase = −0.05 and Wtotal = 0.002.They start from the top of ψ = −0.05having different initial injection angles from 0 • to 30 • with a 5 • interval.(b) Each realization of μ for the seven particles from t = 0 to 1000 ns, showing that they are all different while the histograms are the same.
FIGURE 4. (a) Trajectory (red) of a single particle starting from the top with ψbase = −0.05(magenta), Wtotal = 0.00225 and tan −1 (v y /v x ) = 45 • from t = 0 to 400 ns.The vertical reference line (blue) connects the positions of the two wires where the spatial gradient of the magnetic field is large (cf.figure 1c).(b) Temporal evolution of μ from t = 0 to 2000 ns, and the blue vertical lines indicate when the particles crosses the reference line, and (c) normalized two-dimensional histogram for μ k (abscissa) and μ k+1 (ordinate) indicating the correlation between the two consecutive μ values.

FIGURE 5 .
FIGURE 5. Normalized two-dimensional histograms for μ k (abscissa) and μ k+1 (ordinate) indicating the correlation between the two consecutive μ values with different pairs of the two invariants denoted as [ ψbase , Wtotal / Wbase ] on each plot.
doi.org/10.1017/S0022377824000321Published online by Cambridge University Press FIGURE 6. (a-d) Contour plots of U eff /W total with various pairs of invariants denoted as [ ψbase , Wtotal / Wbase ] on each plot, showing that a particle can migrate to the opposite side of the closed field lines only if Wtotal / Wbase ≥ 1. Grey regions indicate inaccessible areas for the particle.Magenta lines correspond to ψbase = −0.10.

FIGURE 7 .
FIGURE 7. Histograms of time a particle spent in between the consecutive migrations for the invariant pair [ ψbase , Wtotal / Wbase ] of (a) [−0.10, 1.20] and (b) [−0.10, 1.50].Migration confinement time, τ mig , which is the average value from the histogram, is depicted in (a) and (b) with blue vertical dotted lines.
(a), a converging power law at large Wtotal − Wbase with different values of ψbase ; thus, we introduce a α( Wtotal − Wbase ) β factor in the formula.Second, from figure 8(b) we observe a similar slope for different values of ψbase for small Wtotal / Wbase − 1 with different offsets.Therefore, we include a [( Wtotal / Wbase − 1) γ − Wδ base + 1] factor in the formula as well.Together, we have τmig ( ψbase , Wtotal ) = α( Wtotal − Wbase ) β parameters α, β, γ and δ are obtained using the numerically obtained τmig .The black lines in figure 8(a,b) are the fitted lines, which show good agreements.Equipped with the empirical expression (5.5), we show the quantitative trend of τmig as a function of ψbase and Wtotal in figure9.Note that we use the absolute value of ψbase since the confined particles have ψbase < 0. It is clear that τmig is finite only when Wtotal > Wbase = ψ2base as a particle cannot migrate otherwise.The value of τmig decreases with

FIGURE 9 .
FIGURE 9. Value of τmig as a function of | ψbase | and Wtotal obtained by the empirical expression (5.5), showing the quantitative trend of τmig with the two invariants of a particle.

FIGURE 10 .
FIGURE 10.An example of particle losses due to the collisionless migration phenomenon with the loss boundary, acting as a divertor in a tokamak, located at ỹ = −0.5 (thick grey horizontal line).Initially, 30 000 Maxwellian deuterons (green dots) are randomly placed on ψbase = −0.1 (magenta line) on the upper side of the TWM configuration.(a-c) Show the spatial distribution of deuterons at times t = 0.0, 0.2 and 0.5 ms, respectively.(d-f ) Are the histograms of W total at these times.(g,h) Show temporal evolutions of the number of particles and W total for the upper and lower sides, respectively.