-
Observations and modelling of ion cyclotron emission observed in JET plasmas using a sub-harmonic arc detection system during ion cyclotron resonance heating
Authors:
K G McClements,
A Brisset,
B Chapman,
S C Chapman,
R O Dendy,
P Jacquet,
V G Kiptily,
M Mantsinen,
B C G Reman,
JET Contributors
Abstract:
Measurements are reported of electromagnetic emission close to the cyclotron frequency of energetic ions in JET plasmas heated by waves in the ion cyclotron range of frequencies (ICRF). Hydrogen was the majority ion species in all of these plasmas. The measurements were obtained using a sub-harmonic arc detection (SHAD) system in the transmission lines of one of the ICRF antennas. The measured ion…
▽ More
Measurements are reported of electromagnetic emission close to the cyclotron frequency of energetic ions in JET plasmas heated by waves in the ion cyclotron range of frequencies (ICRF). Hydrogen was the majority ion species in all of these plasmas. The measurements were obtained using a sub-harmonic arc detection (SHAD) system in the transmission lines of one of the ICRF antennas. The measured ion cyclotron emission (ICE) spectra were strongly filtered by the antenna system, and typically contained sub-structure, consisting of sets of peaks with a separation of a few kHz, suggesting the excitation of compressional Alfven eigenmodes (CAEs) closely spaced in frequency. In most cases the energetic ions can be clearly identified as ICRF-accelerated He-3 minority ions, although in two pulses the emission may have been produced by energetic He-4 ions, originating from third harmonic ICRF wave acceleration. It is proposed that the emission close to the He-3 cyclotron frequency was produced by energetic ions of this species undergoing drift orbit excursions to the outer midplane plasma edge. Particle-in-cell and hybrid (kinetic ion, fluid electron) simulations using plasma parameters corresponding to edge plasma conditions in these JET pulses, and energetic particle parameters inferred from the cyclotron resonance location, indicate strong excitation of waves at He-3 cyclotron harmonics, including the fundamental, which is identified with the observed emission. These results underline the potential importance of ICE measurements as a method of studying confined fast particles that are strongly suprathermal but have insufficient energies or are not present in sufficient numbers to excite detectable levels of gamma-ray emission or other collective instabilities.
△ Less
Submitted 13 June, 2018;
originally announced June 2018.
-
Sub-microsecond temporal evolution of edge density during edge localized modes in KSTAR tokamak plasmas inferred from ion cyclotron emission
Authors:
B. Chapman,
R. O. Dendy,
K. G. McClements,
S. C. Chapman,
G. S. Yun,
S. G. Thatipamula,
M. H. Kim
Abstract:
Ion cyclotron emission (ICE) is detected during edge localised modes (ELMs) in the KSTAR tokamak at harmonics of the proton cyclotron frequency in the outer plasma edge. The emission typically chirps downward (occasionally upward) during ELM crashes, and is driven by confined 3MeV fusion-born protons that have large drift excursions from the plasma core. We exploit fully kinetic simulations at mul…
▽ More
Ion cyclotron emission (ICE) is detected during edge localised modes (ELMs) in the KSTAR tokamak at harmonics of the proton cyclotron frequency in the outer plasma edge. The emission typically chirps downward (occasionally upward) during ELM crashes, and is driven by confined 3MeV fusion-born protons that have large drift excursions from the plasma core. We exploit fully kinetic simulations at multiple plasma densities to match the time-evolving features of the chirping ICE. This yields a unique, very high time resolution diagnostic of the collapsing edge pedestal density.
△ Less
Submitted 6 February, 2017;
originally announced February 2017.
-
Stimulated emission of fast Alfvén waves within magnetically confined fusion plasmas
Authors:
J W S Cook,
R O Dendy,
S C Chapman
Abstract:
A fast Alfvén wave with finite amplitude is shown to grow by a stimulated emission process that we propose for exploitation in toroidal magnetically confined fusion plasmas. Stimulated emission occurs while the wave propagates inward through the outer mid-plane plasma, where a population inversion of the energy distribution of fusion-born ions is observed to arise naturally. Fully nonlinear first…
▽ More
A fast Alfvén wave with finite amplitude is shown to grow by a stimulated emission process that we propose for exploitation in toroidal magnetically confined fusion plasmas. Stimulated emission occurs while the wave propagates inward through the outer mid-plane plasma, where a population inversion of the energy distribution of fusion-born ions is observed to arise naturally. Fully nonlinear first principles simulations, which self-consistently evolve particles and fields under the Maxwell-Lorentz system, demonstrate this novel "alpha-particle channelling" scenario for the first time.
△ Less
Submitted 1 June, 2016;
originally announced June 2016.
-
Quantifying fusion born ion populations in magnetically confined plasmas using ion cyclotron emission
Authors:
L. Carbajal,
R. O. Dendy,
S. C. Chapman,
J. W. S. Cook
Abstract:
Ion cyclotron emission (ICE) offers unique promise as a diagnostic of the fusion born alpha-particle population in magnetically confined plasmas. Pioneering observations from JET and TFTR found that ICE intensity $P_{ICE}$ scales approximately linearly with the measured neutron flux from fusion reactions, and with the inferred concentration, $n_α/n_i$, of fusion-born alpha-particles confined withi…
▽ More
Ion cyclotron emission (ICE) offers unique promise as a diagnostic of the fusion born alpha-particle population in magnetically confined plasmas. Pioneering observations from JET and TFTR found that ICE intensity $P_{ICE}$ scales approximately linearly with the measured neutron flux from fusion reactions, and with the inferred concentration, $n_α/n_i$, of fusion-born alpha-particles confined within the plasma. We present fully nonlinear self-consistent kinetic simulations that reproduce this scaling for the first time. This resolves a longstanding question in the physics of fusion alpha-particle confinement and stability in MCF plasmas. It confirms the magnetoacoustic cyclotron instability (MCI) as the likely emission mechanism and greatly strengthens the basis for diagnostic exploitation of ICE in future burning plasmas.
△ Less
Submitted 1 June, 2016;
originally announced June 2016.
-
The global build-up to intrinsic ELM bursts seen in divertor full flux loops in Jet
Authors:
S. C. Chapman,
R. O. Dendy,
T. N. Todd,
N. W. Watkins,
F. A. Calderon,
J. Morris,
JET Contributors
Abstract:
A global signature of the build-up to an intrinsic ELM is found in the phase of signals measured in full flux azimuthal loops in the divertor region of JET. Full flux loop signals provide a global measurement proportional to the voltage induced by changes in poloidal magnetic flux; they are electromagnetically induced by the dynamics of spatially integrated current density. We perform direct time-…
▽ More
A global signature of the build-up to an intrinsic ELM is found in the phase of signals measured in full flux azimuthal loops in the divertor region of JET. Full flux loop signals provide a global measurement proportional to the voltage induced by changes in poloidal magnetic flux; they are electromagnetically induced by the dynamics of spatially integrated current density. We perform direct time-domain analysis of the high time-resolution full flux loop signals VLD2 and VLD3. We analyze plasmas where a steady H-mode is sustained over several seconds, during which all the observed ELMs are intrinsic; there is no deliberate intent to pace the ELMing process by external means. ELM occurrence times are determined from the Be II emission at the divertor. We previously found that the occurrence times of intrinsic ELMs correlate with specific phases of the VLD2 and VLD3 signals. Here, we investigate how the VLD2 and VLD3 phases vary with time in advance of the ELM occurrence time. We identify a build-up to the ELM in which the VLD2 and VLD3 signals progressively align to the phase at which ELMs preferentially occur, on a ~ 2 -5ms timescale. At the same time, the VLD2 and VLD3 signals become phase synchronized with each other, consistent with the emergence of coherent global dynamics in the integrated current density. In a plasma that remains close to a global magnetic equilibrium, this can reflect bulk displacement or motion of the plasma. This build-up signature to an intrinsic ELM can be extracted from a time interval of data that does not extend beyond the ELM occurrence time, so that these full flux loop signals could assist in ELM prediction or mitigation.
△ Less
Submitted 10 March, 2015;
originally announced March 2015.
-
Fast particle-driven ion cyclotron emission (ICE) in tokamak plasmas and the case for an ICE diagnostic in ITER
Authors:
K. G. McClements,
R. D'Inca,
R. O. Dendy,
L. Carbajal,
S. C. Chapman,
J. W. S. Cook,
R. W. Harvey,
W. W. Heidbrink,
S. D. Pinches
Abstract:
Fast particle-driven waves in the ion cyclotron frequency range (ion cyclotron emission or ICE) have provided a valuable diagnostic of confined and escaping fast ions in many tokamaks. This is a passive, non-invasive diagnostic that would be compatible with the high radiation environment of deuterium-tritium plasmas in ITER, and could provide important information on fusion α-particles and beam io…
▽ More
Fast particle-driven waves in the ion cyclotron frequency range (ion cyclotron emission or ICE) have provided a valuable diagnostic of confined and escaping fast ions in many tokamaks. This is a passive, non-invasive diagnostic that would be compatible with the high radiation environment of deuterium-tritium plasmas in ITER, and could provide important information on fusion α-particles and beam ions in that device. In JET, ICE from confined fusion products scaled linearly with fusion reaction rate over six orders of magnitude and provided evidence that α-particle confinement was close to classical. In TFTR, ICE was observed from super-Alfvénic α-particles in the plasma edge. The intensity of beam-driven ICE in DIII-D is more strongly correlated with drops in neutron rate during fishbone excitation than signals from more direct beam ion loss diagnostics. In ASDEX Upgrade ICE is produced by both super-Alfvénic DD fusion products and sub-Alfvénic deuterium beam ions.
△ Less
Submitted 12 December, 2014;
originally announced December 2014.
-
Electron kinetics inferred from observations of microwave bursts during edge localised modes in the Mega-Amp Spherical Tokamak
Authors:
S. J. Freethy,
K. G. McClements,
S. C. Chapman,
R. O. Dendy,
W. N. Lai,
S. J. P. Pamela,
V. F. Shevchenko,
R. G. L. Vann
Abstract:
Recent measurements of microwave and X-ray emission during edge localised mode (ELM) activity in tokamak plasmas provide a fresh perspective on ELM physics. It is evident that electron kinetics, which are not incorporated in standard (fluid) models for the instability that drives ELMs, play a key role in the new observations. These effects should be included in future models for ELMs and the ELM c…
▽ More
Recent measurements of microwave and X-ray emission during edge localised mode (ELM) activity in tokamak plasmas provide a fresh perspective on ELM physics. It is evident that electron kinetics, which are not incorporated in standard (fluid) models for the instability that drives ELMs, play a key role in the new observations. These effects should be included in future models for ELMs and the ELM cycle. The observed radiative effects paradoxically imply acceleration of electrons parallel to the magnetic field combined with rapid acquisition of perpendicular momentum. It is shown that this paradox can be resolved by the action of the anomalous Doppler instability which enables fast collective radiative relaxation, in the perpendicular direction, of electrons accelerated in the parallel direction by inductive electric fields generated by the initial ELM instability.
△ Less
Submitted 4 November, 2014;
originally announced November 2014.
-
Ion cyclotron emission from fusion-born ions in large tokamak plasmas: a brief review from JET and TFTR to ITER
Authors:
R O Dendy,
K G McClements
Abstract:
Ion cyclotron emission (ICE) was the first collective radiative instability, driven by confined fusion-born ions, observed from deuterium-tritium plasmas in JET and TFTR. ICE comprises strongly suprathermal emission, which has spectral peaks at multiple ion cyclotron harmonic frequencies as evaluated at the outer mid-plane edge of tokamak plasmas. The measured intensity of ICE spectral peaks scale…
▽ More
Ion cyclotron emission (ICE) was the first collective radiative instability, driven by confined fusion-born ions, observed from deuterium-tritium plasmas in JET and TFTR. ICE comprises strongly suprathermal emission, which has spectral peaks at multiple ion cyclotron harmonic frequencies as evaluated at the outer mid-plane edge of tokamak plasmas. The measured intensity of ICE spectral peaks scaled linearly with measured fusion reactivity in JET. In other large tokamak plasmas, ICE is currently used as an indicator of fast ions physics. The excitation mechanism for ICE is the magnetoacoustic cyclotron instability (MCI); in the case of JET and TFTR, the MCI is driven by a set of centrally born fusion products, lying just inside the trapped-passing boundary in velocity space, whose drift orbits make large radial excursions to the outer mid-plane edge. Diagnostic exploitation of ICE in future experiments therefore rests in part on deep understanding of the MCI, and recent advances in computational plasma physics have led to substantial recent progress, reviewed here.
△ Less
Submitted 22 October, 2014;
originally announced October 2014.
-
Time-Resonant Tokamak Plasma Edge Instabilities?
Authors:
A. J. Webster,
R. O. Dendy,
F. A. Calderon,
S. C. Chapman,
E. Delabie,
D. Dodt,
R. Felton,
T. N. Todd,
F. Maviglia,
J. Morris,
V. Riccardo,
B. Alper,
S Brezinsek,
P. Coad,
J. Likonen,
M. Rubel,
JET EFDA Contributors
Abstract:
For a two week period during the Joint European Torus (JET) 2012 experimental campaign, the same high confinement plasma was repeated 151 times. The dataset was analysed to produce a probability density function (pdf) for the waiting times between edge-localised plasma instabilities ("ELMS"). The result was entirely unexpected. Instead of a smooth single peaked pdf, a succession of 4-5 sharp maxim…
▽ More
For a two week period during the Joint European Torus (JET) 2012 experimental campaign, the same high confinement plasma was repeated 151 times. The dataset was analysed to produce a probability density function (pdf) for the waiting times between edge-localised plasma instabilities ("ELMS"). The result was entirely unexpected. Instead of a smooth single peaked pdf, a succession of 4-5 sharp maxima and minima uniformly separated by 7-8 millisecond intervals was found. Here we explore the causes of this newly observed phenomenon, and conclude that it is either due to a self-organised plasma phenomenon or an interaction between the plasma and a real-time control system. If the maxima are a result of "resonant" frequencies at which ELMs can be triggered more easily, then future ELM control techniques can, and probably will, use them. Either way, these results demand a deeper understanding of the ELMing process.
△ Less
Submitted 1 October, 2013;
originally announced October 2013.
-
Transitions to improved confinement regimes induced by changes in heating in zero-dimensional models for tokamak plasmas
Authors:
H. Zhu,
S. C. Chapman,
R. O. Dendy,
K. Itoh
Abstract:
It is shown that rapid substantial changes in heating rate can induce transitions to improved energy confinement regimes in zero-dimensional models for tokamak plasma phenomenology. We examine for the first time the effect of step changes in heating rate in the models of E-J.Kim and P.H.Diamond, Phys.Rev.Lett. 90, 185006 (2003) and M.A.Malkov and P.H.Diamond, Phys.Plasmas 16, 012504 (2009) which n…
▽ More
It is shown that rapid substantial changes in heating rate can induce transitions to improved energy confinement regimes in zero-dimensional models for tokamak plasma phenomenology. We examine for the first time the effect of step changes in heating rate in the models of E-J.Kim and P.H.Diamond, Phys.Rev.Lett. 90, 185006 (2003) and M.A.Malkov and P.H.Diamond, Phys.Plasmas 16, 012504 (2009) which nonlinearly couple the evolving temperature gradient, micro-turbulence and a mesoscale flow; and in the extension of H.Zhu, S.C.Chapman and R.O.Dendy, Phys.Plasmas 20, 042302 (2013), which couples to a second mesoscale flow component. The temperature gradient rises, as does the confinement time defined by analogy with the fusion context, while micro-turbulence is suppressed. This outcome is robust against variation of heating rise time and against introduction of an additional variable into the model. It is also demonstrated that oscillating changes in heating rate can drive the level of micro-turbulence through a period-doubling path to chaos, where the amplitude of the oscillatory component of the heating rate is the control parameter.
△ Less
Submitted 18 June, 2014; v1 submitted 7 September, 2013;
originally announced September 2013.
-
Identifying low-dimensional dynamics in Type-I edge-localised-mode processes in JET plasmas
Authors:
F. A. Calderon,
R. O. Dendy,
S. C. Chapman,
A. J. Webster,
B. Alper,
R. M. Nicol,
JET EFDA Contributors
Abstract:
Edge localised mode (ELM) measurements from reproducibly similar plasmas in the Joint European Torus (JET) tokamak, which differ only in their gas puffing rate, are analysed in terms of the pattern in the sequence of inter-ELM time intervals. It is found that the category of ELM defined empirically as Type I - typically more regular, less frequent, and having larger amplitude than other ELM types…
▽ More
Edge localised mode (ELM) measurements from reproducibly similar plasmas in the Joint European Torus (JET) tokamak, which differ only in their gas puffing rate, are analysed in terms of the pattern in the sequence of inter-ELM time intervals. It is found that the category of ELM defined empirically as Type I - typically more regular, less frequent, and having larger amplitude than other ELM types - embraces substantially different ELMing processes. By quantifying the structure in the sequence of inter-ELM time intervals using delay time plots, we reveal transitions between distinct phase space dynamics, implying transitions between distinct underlying physical processes. The control parameter for these transitions between these different ELMing processes is the gas puffing rate.
△ Less
Submitted 8 April, 2013;
originally announced April 2013.
-
Statistical Characterisation and Classification of Edge Localised Plasma Instabilities
Authors:
A. J. Webster,
R. O. Dendy
Abstract:
The statistics of edge-localised plasma instabilities (ELMs) in toroidal magnetically confined fusion plasmas are considered. From first principles, standard experimentally motivated assumptions are shown to determine a specific probability distribution for the waiting times between ELMs: the Weibull distribution. This is confirmed empirically by a statistically rigorous comparison with a large da…
▽ More
The statistics of edge-localised plasma instabilities (ELMs) in toroidal magnetically confined fusion plasmas are considered. From first principles, standard experimentally motivated assumptions are shown to determine a specific probability distribution for the waiting times between ELMs: the Weibull distribution. This is confirmed empirically by a statistically rigorous comparison with a large data set from the Joint European Torus (JET). The successful characterisation of ELM waiting times enables future work to progress in various ways. Here we present a quantitative classification of ELM types, complementary to phenomenological approaches. It also informs us about the nature of ELMing processes, such as whether they are random or deterministic.
△ Less
Submitted 26 February, 2013;
originally announced February 2013.
-
Robustness of predator-prey models for confinement regime transitions in fusion plasmas
Authors:
H. Zhu,
S. C. Chapman,
R. O. Dendy
Abstract:
Energy transport and confinement in tokamak fusion plasmas is usually determined by the coupled nonlinear interactions of small-scale drift turbulence and larger scale coherent nonlinear structures, such as zonal flows, together with free energy sources such as temperature gradients. Zero-dimensional models, designed to embody plausible physical narratives for these interactions, can help identify…
▽ More
Energy transport and confinement in tokamak fusion plasmas is usually determined by the coupled nonlinear interactions of small-scale drift turbulence and larger scale coherent nonlinear structures, such as zonal flows, together with free energy sources such as temperature gradients. Zero-dimensional models, designed to embody plausible physical narratives for these interactions, can help identify the origin of enhanced energy confinement and of transitions between confinement regimes. A prime zero-dimensional paradigm is predator-prey or Lotka-Volterra. Here we extend a successful three-variable (temperature gradient; microturbulence level; one class of coherent structure) model in this genre [M A Malkov and P H Diamond, Phys. Plasmas 16, 012504 (2009)], by adding a fourth variable representing a second class of coherent structure. This requires a fourth coupled nonlinear ordinary differential equation. We investigate the degree of invariance of the phenomenology generated by the model of Malkov and Diamond, given this additional physics. We study and compare the long-time behaviour of the three-equation and four-equation systems, their evolution towards the final state, and their attractive fixed points and limit cycles. We explore the sensitivity of paths to attractors. It is found that, for example, an attractive fixed point of the three-equation system can become a limit cycle of the four-equation system. Addressing these questions - which we together refer to as "robustness" for convenience, is particularly important for models which, as here, generate sharp transitions in the values of system variables which may replicate some key features of confinement transitions. Our results help establish the robustness of the zero-dimensional model approach to capturing observed confinement phenomenology in tokamak fusion plasmas.
△ Less
Submitted 4 February, 2013;
originally announced February 2013.
-
Gyrobunching and wave-particle resonance in the lower hybrid drift instability
Authors:
J W S Cook,
R O Dendy,
S C Chapman
Abstract:
We report a first principles study of the coupled evolution of energetic ions, background majority ions, electrons and electromagnetic fields in magnetised plasma during the linear phase of the lower hybrid drift instability. A particle-in-cell code, with one spatial and three velocity space co-ordinates, is used to analyse the evolving distribution of a drifting ring-beam population of energetic…
▽ More
We report a first principles study of the coupled evolution of energetic ions, background majority ions, electrons and electromagnetic fields in magnetised plasma during the linear phase of the lower hybrid drift instability. A particle-in-cell code, with one spatial and three velocity space co-ordinates, is used to analyse the evolving distribution of a drifting ring-beam population of energetic protons in physical space and gyrophase angle. This analysis is carried out for plasma parameters that approximate to edge conditions in large tokamaks, in a scenario that is motivated by observations of ion cyclotron emission and may be relevant to alpha channelling. Resonant energy transfer occurs at the two gyrophase angles at which the instantaneous speed of an energetic proton on its cyclotron orbit precisely matches the phase velocity of the lower hybrid wave along the simulation domain. Electron space-charge oscillations determine the wavelength of the propagating lower hybrid wave, and thereby govern the spatial distribution of gyrobunching of the energetic protons that drive the instability.
△ Less
Submitted 11 January, 2011;
originally announced January 2011.
-
Self-consistent kinetic simulations of lower hybrid drift instability resulting in electron current driven by fusion products in tokamak plasmas
Authors:
J W S Cook,
S C Chapman,
R O Dendy
Abstract:
We present particle-in-cell (PIC) simulations of minority energetic protons in deuterium plasmas, which demonstrate a collective instability responsible for emission near the lower hybrid frequency and its harmonics. The simulations capture the lower hybrid drift instability in a regime relevant to tokamak fusion plasmas, and show further that the excited electromagnetic fields collectively and co…
▽ More
We present particle-in-cell (PIC) simulations of minority energetic protons in deuterium plasmas, which demonstrate a collective instability responsible for emission near the lower hybrid frequency and its harmonics. The simulations capture the lower hybrid drift instability in a regime relevant to tokamak fusion plasmas, and show further that the excited electromagnetic fields collectively and collisionlessly couple free energy from the protons to directed electron motion. This results in an asymmetric tail antiparallel to the magnetic field. We focus on obliquely propagating modes under conditions approximating the outer mid-plane edge in a large tokamak, through which there pass confined centrally born fusion products on banana orbits that have large radial excursions. A fully self-consistent electromagnetic relativistic PIC code representing all vector field quantities and particle velocities in three dimensions as functions of a single spatial dimension is used to model this situation, by evolving the initial antiparallel travelling ring-beam distribution of 3MeV protons in a background 10keV Maxwellian deuterium plasma with realistic ion-electron mass ratio. The simulations thus demonstrate a key building block of alpha channelling scenarios for burning fusion plasmas in tokamaks.
△ Less
Submitted 11 October, 2010;
originally announced October 2010.
-
Electron current drive by fusion-product-excited lower hybrid drift instability
Authors:
J. W. S. Cook,
S. C. Chapman,
R. O. Dendy
Abstract:
We present first principles simulations of the direct collisionless coupling of the free energy of fusion-born ions into electron current in a magnetically confined fusion plasma. These simulations demonstrate, for the first time, a key building block of some "alpha channelling" scenarios for tokamak experiments. A fully self-consistent electromagnetic 1D3V particle-in-cell code is used to evolve…
▽ More
We present first principles simulations of the direct collisionless coupling of the free energy of fusion-born ions into electron current in a magnetically confined fusion plasma. These simulations demonstrate, for the first time, a key building block of some "alpha channelling" scenarios for tokamak experiments. A fully self-consistent electromagnetic 1D3V particle-in-cell code is used to evolve a parallel drifting ring-beam distribution of 3MeV protons in a 10keV thermal deuterium-electron plasma with realistic mass ratio. Collective instability gives rise to electromagnetic field activity in the lower hybrid range of frequencies. These spontaneously excited obliquely propagating waves undergo Landau damping on resonant electrons, drawing out an asymmetric tail in the distribution of electron parallel velocities, which carries a current.
△ Less
Submitted 6 September, 2010;
originally announced September 2010.
-
Quantifying the Anisotropy and Solar Cycle Dependence of "$1/f$" Solar Wind Fluctuations Observed by Ace
Authors:
R. M. Nicol,
S. C. Chapman,
R. O. Dendy
Abstract:
The power spectrum of the evolving solar wind shows evidence of a spectral break between an inertial range of turbulent fluctuations at higher frequencies and a "$1/f$" like region at lower frequencies. In the ecliptic plane at $\sim 1$ AU, this break occurs approximately at timescales of a few hours, and is observed in the power spectra of components of velocity and magnetic field. The "$1/f$"…
▽ More
The power spectrum of the evolving solar wind shows evidence of a spectral break between an inertial range of turbulent fluctuations at higher frequencies and a "$1/f$" like region at lower frequencies. In the ecliptic plane at $\sim 1$ AU, this break occurs approximately at timescales of a few hours, and is observed in the power spectra of components of velocity and magnetic field. The "$1/f$" energy range is of more direct coronal origin than the inertial range, and carries signatures of the complex magnetic field structure of the solar corona, and of footpoint stirring in the solar photosphere. To quantify the scaling properties we use generic statistical methods such as generalised structure functions and PDFs, focusing on solar cycle dependence and on anisotropy with respect to the background magnetic field. We present structure function analysis of magnetic and velocity field fluctuations, using a novel technique to decompose the fluctuations into directions parallel and perpendicular to the mean local background magnetic field. Whilst the magnetic field is close to "$1/f$", we show that the velocity field is "$1/f^α$" with $α\neq1$. For the velocity, the value of $α$ varies between parallel and perpendicular fluctuations and with the solar cycle. There is also variation in $α$ with solar wind speed. We have examined the PDFs in the fast, quiet solar wind and intriguingly, whilst parallel and perpendicular are distinct, both the ${\boldmath$B$}$ field and velocity show the same PDF of their perpendicular flucutations, which is close to gamma or inverse Gumbel. These results point to distinct physical processes in the corona, and to their mapping out into the solar wind. The scaling exponents obtained constrain the models for these processes.
△ Less
Submitted 18 June, 2009;
originally announced June 2009.
-
The signature of evolving turbulence in quiet solar wind as seen by ULYSSES
Authors:
R. M. Nicol,
S. C. Chapman,
R. O. Dendy
Abstract:
Solar wind fluctuations, such as magnetic field or velocity, show power law power spectra suggestive both of an inertial range of intermittent turbulence (with $\sim -5/3$ exponent) and at lower frequencies, of fluctuations of coronal origin (with $\sim -1$ exponent). The ULYSSES spacecraft spent many months in the quiet fast solar wind above the Sun's polar coronal holes in a highly ordered mag…
▽ More
Solar wind fluctuations, such as magnetic field or velocity, show power law power spectra suggestive both of an inertial range of intermittent turbulence (with $\sim -5/3$ exponent) and at lower frequencies, of fluctuations of coronal origin (with $\sim -1$ exponent). The ULYSSES spacecraft spent many months in the quiet fast solar wind above the Sun's polar coronal holes in a highly ordered magnetic field. We use statistical analysis methods such as the generalized structure function (GSF) and extended self-similarity (ESS) to quantify the scaling of the moments of the probability density function of fluctuations in the magnetic field. The GSFs give power law scaling in the ``$f^{-1}$'' range of the form $<| y(t+τ)-y(t)|^{m}>\simτ^{ζ(m)}$, but ESS is required to reveal scaling in the inertial range, which is of the form $<| y(t+τ)-y(t)|^{m}>\sim [g(τ)]^{ζ(m)}$. We find that $g(τ)$ is independent of spacecraft position and $g(τ)\simτ^{-log_{10}(\tildeλτ)}$. The ``$f^{-1}$'' scaling fluctuates with radial spacecraft position. This confirms that, whereas the ``$f^{-1}$'' fluctuations are directly influenced by the corona, the inertial range fluctuations are consistent with locally evolving turbulence, but with an ``envelope'' $g(τ)$, which captures the formation of the quiet fast solar wind.
△ Less
Submitted 18 February, 2008; v1 submitted 5 December, 2007;
originally announced December 2007.
-
Solar cycle dependence of spatial correlation in the solar wind
Authors:
R. T. Wicks,
S. C. Chapman,
R. O. Dendy
Abstract:
We investigate the spatial correlation properties of the solar wind using simultaneous observations by the ACE and WIND spacecraft. We use mutual information as a nonlinear measure of correlation and compare this to linear correlation. We find that the correlation lengthscales of fluctuations in density and magnetic field magnitude vary strongly with the solar cycle, whereas correlation lengths…
▽ More
We investigate the spatial correlation properties of the solar wind using simultaneous observations by the ACE and WIND spacecraft. We use mutual information as a nonlinear measure of correlation and compare this to linear correlation. We find that the correlation lengthscales of fluctuations in density and magnetic field magnitude vary strongly with the solar cycle, whereas correlation lengths of fluctuations in B field components do not. We find the correlation length of |B| ~ 120 Re at solar minimum and ~ 270 Re at maximum and the correlation length of density ~ 75 Re at minimum and ~ 170 Re at minimum. The components of the B field have correlation lengths ~ correlation length |B| at minimum.
△ Less
Submitted 3 December, 2007; v1 submitted 29 November, 2007;
originally announced November 2007.
-
Mutual Information as a Tool for Identifying Phase Transitions in Dynamical Complex Systems With Limited Data
Authors:
R. T. Wicks,
S. C. Chapman,
R. O. Dendy
Abstract:
We use a well known model (T. Vicsek et al. Phys Rev Lett 15, 1226 (1995)) for flocking to test mutual information as a tool for detecting order-disorder transitions, in particular when observations of the system are limited. We show that mutual information is a sensitive indicator of the phase transition location, in terms of the natural dimensionless parameters of the system which we have iden…
▽ More
We use a well known model (T. Vicsek et al. Phys Rev Lett 15, 1226 (1995)) for flocking to test mutual information as a tool for detecting order-disorder transitions, in particular when observations of the system are limited. We show that mutual information is a sensitive indicator of the phase transition location, in terms of the natural dimensionless parameters of the system which we have identified. When only a few particles are tracked, and when only a subset of the positional and velocity components are available, mutual information provides a better measure of the phase transition location than the susceptibility of the data.
△ Less
Submitted 13 April, 2007; v1 submitted 20 December, 2006;
originally announced December 2006.
-
Aspects of Electron Acoustic Wave Physics in Laser Backscatter from Plasmas
Authors:
N. J. Sircombe,
T. D. Arber,
R. O. Dendy
Abstract:
Recent experimental results from the Trident laser confirm the importance of kinetic effects in determining laser reflectivities at high intensities. Examples observed include scattering from low frequency electron acoustic waves (EAWs), and the first few stages of a cascade towards turbulence through the Langmuir decay instability. Interpretive and predictive computational capability in this ar…
▽ More
Recent experimental results from the Trident laser confirm the importance of kinetic effects in determining laser reflectivities at high intensities. Examples observed include scattering from low frequency electron acoustic waves (EAWs), and the first few stages of a cascade towards turbulence through the Langmuir decay instability. Interpretive and predictive computational capability in this area is assisted by the development of Vlasov codes, which offer high velocity space resolution in high energy regions of particle phase space, and do not require analytical pre-processing of the fundamental equations. A direct Vlasov solver, capable of resolving these kinetic processes, is used here to address fundamental aspects of the existence and stability of the electron acoustic wave, together with its collective scattering properties. These simulations are extended to realistic laser and plasma parameters characteristic of single hot-spot experiments. Results are in qualitative agreement with experiments displaying both stimulated Raman and stimulated electron acoustic scattering. The amplitude of simulated EAWs is greater than that observed experimentally, and is accompanied by a higher phase velocity. These minor differences can be attributed to the limitations of a one-dimensional collisionless model.
△ Less
Submitted 24 March, 2006;
originally announced March 2006.
-
Phase speed of electrostatic waves: The critical parameter for efficient electron surfing acceleration
Authors:
M E Dieckmann,
N J Sircombe,
M Parviainen,
P K Shukla,
R O Dendy
Abstract:
Particle acceleration by means of non-linear plasma wave interactions is of great topical interest. Accordingly, in this paper we focus on the electron surfing process. Self-consistent kinetic simulations, using both relativistic Vlasov and PIC (Particle In Cell) approaches, show here that electrons can be accelerated to highly relativistic energies (up to 100 m_e c^2) if the phase speed of the…
▽ More
Particle acceleration by means of non-linear plasma wave interactions is of great topical interest. Accordingly, in this paper we focus on the electron surfing process. Self-consistent kinetic simulations, using both relativistic Vlasov and PIC (Particle In Cell) approaches, show here that electrons can be accelerated to highly relativistic energies (up to 100 m_e c^2) if the phase speed of the electrostatic wave is mildly relativistic (0.6c to 0.9c for the magnetic field strengths considered). The acceleration is strong because of relativistic stabilisation of the nonlinearly saturated electrostatic wave, seen in both relativistic Vlasov and PIC simulations. An inverse power law momentum distribution can arise for the most strongly accelerated electrons. These results are of relevance to observed rapid changes in the radio synchrotron emission intensities from microquasars, gamma ray bursts and other astrophysical objects that require rapid acceleration mechanisms for electrons.
△ Less
Submitted 15 February, 2006;
originally announced February 2006.
-
Recurrence plot statistics and the effect of embedding
Authors:
T. K. March,
S. C. Chapman,
R. O. Dendy
Abstract:
Recurrence plots provide a graphical representation of the recurrent patterns in a timeseries, the quantification of which is a relatively new field. Here we derive analytical expressions which relate the values of key statistics, notably determinism and entropy of line length distribution, to the correlation sum as a function of embedding dimension. These expressions are obtained by deriving th…
▽ More
Recurrence plots provide a graphical representation of the recurrent patterns in a timeseries, the quantification of which is a relatively new field. Here we derive analytical expressions which relate the values of key statistics, notably determinism and entropy of line length distribution, to the correlation sum as a function of embedding dimension. These expressions are obtained by deriving the transformation which generates an embedded recurrence plot from an unembedded plot. A single unembedded recurrence plot thus provides the statistics of all possible embedded recurrence plots. If the correlation sum scales exponentially with embedding dimension, we show that these statistics are determined entirely by the exponent of the exponential. This explains the results of Iwanski and Bradley (Chaos 8 [1998] 861-871) who found that certain recurrence plot statistics are apparently invariant to embedding dimension for certain low-dimensional systems. We also examine the relationship between the mutual information content of two timeseries and the common recurrent structure seen in their recurrence plots. This allows time-localized contributions to mutual information to be visualized. This technique is demonstrated using geomagnetic index data; we show that the AU and AL geomagnetic indices share half their information, and find the timescale on which mutual features appear.
△ Less
Submitted 9 February, 2005;
originally announced February 2005.
-
Mutual information between geomagnetic indices and the solar wind as seen by WIND - implications for propagation time estimates
Authors:
T. K. March,
S. C. Chapman,
R. O. Dendy
Abstract:
The determination of delay times of solar wind conditions at the sunward libration point to effects on Earth is investigated using mutual information. This measures the amount of information shared between two timeseries. We consider the mutual information content of solar wind observations, from WIND, and the geomagnetic indices. The success of five commonly used schemes for estimating interpla…
▽ More
The determination of delay times of solar wind conditions at the sunward libration point to effects on Earth is investigated using mutual information. This measures the amount of information shared between two timeseries. We consider the mutual information content of solar wind observations, from WIND, and the geomagnetic indices. The success of five commonly used schemes for estimating interplanetary propagation times is examined. Propagation assuming a fixed plane normal at 45 degrees to the GSE x-axis (i.e. the Parker Spiral estimate) is found to give optimal mutual information. The mutual information depends on the point in space chosen as the target for the propagation estimate, and we find that it is maximized by choosing a point in the nightside rather than dayside magnetosphere. In addition, we employ recurrence plot analysis to visualize contributions to the mutual information, this suggests that it appears on timescales of hours rather than minutes.
△ Less
Submitted 8 February, 2005;
originally announced February 2005.
-
The scaling properties of dissipation in incompressible isotropic three-dimensional magnetohydrodynamic turbulence
Authors:
J. A. Merrifield,
W. -C. Muller,
S. C. Chapman,
R. O. Dendy
Abstract:
The statistical properties of the dissipation process constrain the analysis of large scale numerical simulations of three dimensional incompressible magnetohydrodynamic (MHD) turbulence, such as those of Biskamp and Muller [Phys. Plasmas 7, 4889 (2000)]. The structure functions of the turbulent flow are expected to display statistical self-similarity, but the relatively low Reynolds numbers att…
▽ More
The statistical properties of the dissipation process constrain the analysis of large scale numerical simulations of three dimensional incompressible magnetohydrodynamic (MHD) turbulence, such as those of Biskamp and Muller [Phys. Plasmas 7, 4889 (2000)]. The structure functions of the turbulent flow are expected to display statistical self-similarity, but the relatively low Reynolds numbers attainable by direct numerical simulation, combined with the finite size of the system, make this difficult to measure directly. However, it is known that extended self-similarity, which constrains the ratio of scaling exponents of structure functions of different orders, is well satisfied. This implies the extension of physical scaling arguments beyond the inertial range into the dissipation range. The present work focuses on the scaling properties of the dissipation process itself. This provides an important consistency check in that we find that the ratio of dissipation structure function exponents is that predicted by the She and Leveque [Phys. Rev. Lett 72, 336 (1994)] theory proposed by Biskamp and Muller. This supplies further evidence that the cascade mechanism in three dimensional MHD turbulence is non-linear random eddy scrambling, with the level of intermittency determined by dissipation through the formation of current sheets.
△ Less
Submitted 12 January, 2005; v1 submitted 2 December, 2004;
originally announced December 2004.
-
Accelerated electron populations formed by Langmuir wave-caviton interactions
Authors:
N. J. Sircombe,
T. D. Arber,
R. O. Dendy
Abstract:
Direct numerical simulations of electron dynamics in externally driven electrostatic waves have been carried out using a relativistic two-fluid one-dimensional Vlasov-Poisson code. When the driver wave has sufficiently large amplitude, ion density holes (cavitons) form. The interaction between these cavitons and other incoming Langmuir waves gives rise to substantial local acceleration of groups…
▽ More
Direct numerical simulations of electron dynamics in externally driven electrostatic waves have been carried out using a relativistic two-fluid one-dimensional Vlasov-Poisson code. When the driver wave has sufficiently large amplitude, ion density holes (cavitons) form. The interaction between these cavitons and other incoming Langmuir waves gives rise to substantial local acceleration of groups of electrons, and fine jet-like structures arise in electron phase space. We show that these jets are caused by wave-breaking when finite amplitude Langmuir waves experience the ion density gradient at the leading edge of the holes, and are not caused by caviton burn-out. An analytical two-fluid model gives the critical density gradient and caviton depth for which this process can occur. In particular, the density gradient critically affects the rate at which a Langmuir wave, moving into the caviton, undergoes Landau damping. This treatment also enables us to derive analytical estimates for the maximum energy of accelerated electrons, and for the energy spectrum along a phase-space jet. These are confirmed by direct numerical simulations.
△ Less
Submitted 20 October, 2004;
originally announced October 2004.
-
Ion acceleration processes at reforming collisionless shocks
Authors:
R. E. Lee,
S. C. Chapman,
R. O. Dendy
Abstract:
The identification of pre-acceleration mechanisms for cosmic ray ions in supernova remnant shocks is an important problem in astrophysics. Recent particle-in-cell (PIC) shock simulations have shown that inclusion of the full electron kinetics yields non-time-stationary solutions, in contrast to previous hybrid (kinetic ions, fluid electrons) simulations. Here, by running a PIC code at high phase…
▽ More
The identification of pre-acceleration mechanisms for cosmic ray ions in supernova remnant shocks is an important problem in astrophysics. Recent particle-in-cell (PIC) shock simulations have shown that inclusion of the full electron kinetics yields non-time-stationary solutions, in contrast to previous hybrid (kinetic ions, fluid electrons) simulations. Here, by running a PIC code at high phase space resolution, ion acceleration mechanisms associated with the time dependence of a supercritical collisionless perpendicular shock are examined. In particular the components of $\int \mathbf{F} \cdot \mathbf{v} dt$ are analysed along trajectories for ions that reach both high and low energies. Selection mechanisms for the ions that reach high energies are also examined. In contrast to quasi-stationary shock solutions, the suprathermal protons are selected from the background population on the basis of the time at which they arrive at the shock, and thus are generated in bursts.
△ Less
Submitted 25 October, 2004;
originally announced October 2004.
-
Statistical characterisation of full-disk EUV/XUV solar irradiance and correlation with solar activity
Authors:
J. Greenhough,
S. C. Chapman,
R. O. Dendy,
V. M. Nakariakov,
G. Rowlands
Abstract:
We investigate the distribution of fluctuations in solar irradiance when integrated over the full disk, obtained using extreme ultraviolet/soft X-ray observations from the SOHO CELIAS/SEM instrument. This time series sums over both the contributions of single distinguishable flares, and of many other processes. By detrending we select events with timescales of less than a few hours such as waves…
▽ More
We investigate the distribution of fluctuations in solar irradiance when integrated over the full disk, obtained using extreme ultraviolet/soft X-ray observations from the SOHO CELIAS/SEM instrument. This time series sums over both the contributions of single distinguishable flares, and of many other processes. By detrending we select events with timescales of less than a few hours such as waves, slow flows, and CMEs. The statistics generated by this range of phenomena can be characterised by power-law-tailed distributions. We show that (i) during the high-activity period 2000 Jan-June the tail exponent a(T)=1.5 +/- 0.1; (ii) during the low-activity period 1996 Jan-June a(T)=3.0 +/- 0.2; and (iii) in general a(T) decreases with increasing activity.
△ Less
Submitted 28 July, 2003; v1 submitted 14 February, 2003;
originally announced February 2003.
-
Solar Flares as Cascades of Reconnecting Magnetic Loops
Authors:
D. Hughes,
M. Paczuski,
R. O. Dendy,
P. Helander,
K. G. McClements
Abstract:
A model for the solar coronal magnetic field is proposed where multiple directed loops evolve in space and time. Loops injected at small scales are anchored by footpoints of opposite polarity moving randomly on a surface. Nearby footpoints of the same polarity aggregate, and loops can reconnect when they collide. This may trigger a cascade of further reconnection, representing a solar flare. Num…
▽ More
A model for the solar coronal magnetic field is proposed where multiple directed loops evolve in space and time. Loops injected at small scales are anchored by footpoints of opposite polarity moving randomly on a surface. Nearby footpoints of the same polarity aggregate, and loops can reconnect when they collide. This may trigger a cascade of further reconnection, representing a solar flare. Numerical simulations show that a power law distribution of flare energies emerges, associated with a scale free network of loops, indicating self-organized criticality.
△ Less
Submitted 21 February, 2003; v1 submitted 9 October, 2002;
originally announced October 2002.
-
Identification of a 12-17 day timescale in X-ray observations of GRS 1915+105
Authors:
J. Greenhough,
S. C. Chapman,
S. Chaty,
R. O. Dendy,
G. Rowlands
Abstract:
Measurement of the energy emitted from accreting astrophysical systems provides an observational constraint on the plasma processes that may be operating within the disk. Here we examine the continual time variation over the past six years of the total X-ray flux from the microquasar GRS 1915+105. The application of differencing and rescaling techniques to RXTE/ASM data shows that the small ampl…
▽ More
Measurement of the energy emitted from accreting astrophysical systems provides an observational constraint on the plasma processes that may be operating within the disk. Here we examine the continual time variation over the past six years of the total X-ray flux from the microquasar GRS 1915+105. The application of differencing and rescaling techniques to RXTE/ASM data shows that the small amplitude fluctuations scale up to 12-17 days. A 17-day timescale in the X-ray fluctuations corresponds to half the measured binary orbital period of this system (33.5 +/- 1.5 days). While this may be coincidental, it is possible that these two timescales may be linked by, for example, a turbulent cascade in the accretion disk driven by a tidally-induced two-armed spiral shock corotating with the binary system. Temporal scaling is found only in the ever-present small fluctuations, and not in the intermittent larger-amplitude fluctuations. This is consistent with the basic model for this source which consists of a steady, cold outer disk and an unstable inner disk.
△ Less
Submitted 9 December, 2002; v1 submitted 29 June, 2002;
originally announced July 2002.
-
Self Organization of Edge and Internal Pedestals in a Sandpile
Authors:
S. C. Chapman,
R. O. Dendy,
B. Hnat
Abstract:
The temperature profiles of magnetically confined plasmas can display distinctive longlived pedestals at the edge and internally. Here we show that such structures can arise naturally through avalanching transport in a sandpile model. A single control parameter that is constant across the sandpile determines the occurrence and regularity of these effects, as well as the entrainment of global con…
▽ More
The temperature profiles of magnetically confined plasmas can display distinctive longlived pedestals at the edge and internally. Here we show that such structures can arise naturally through avalanching transport in a sandpile model. A single control parameter that is constant across the sandpile determines the occurrence and regularity of these effects, as well as the entrainment of global confinement to edge pedestal dynamics. The system dynamics indicate that the pedestals are a consequence of an inverse cascade in real space, and that self organization is necessary for their occurrence.
△ Less
Submitted 17 October, 2001;
originally announced October 2001.
-
Characterising anomalous transport in accretion disks from X-ray observations
Authors:
J. Greenhough,
S. C. Chapman,
S. Chaty,
R. O. Dendy,
G. Rowlands
Abstract:
Whilst direct observations of internal transport in accretion disks are not yet possible, measurement of the energy emitted from accreting astrophysical systems can provide useful information on the physical mechanisms at work. Here we examine the unbroken multi-year time variation of the total X-ray flux from three sources: Cygnus X-1, the microquasar GRS1915+105, and for comparison the nonaccr…
▽ More
Whilst direct observations of internal transport in accretion disks are not yet possible, measurement of the energy emitted from accreting astrophysical systems can provide useful information on the physical mechanisms at work. Here we examine the unbroken multi-year time variation of the total X-ray flux from three sources: Cygnus X-1, the microquasar GRS1915+105, and for comparison the nonaccreting Crab nebula. To complement previous analyses, we demonstrate that the application of advanced statistical methods to these observational time-series reveals important contrasts in the nature and scaling properties of the transport processes operating within these sources. We find the Crab signal resembles Gaussian noise; the Cygnus X-1 signal is a leptokurtic random walk whose self-similar properties persist on timescales up to three years; and the GRS1915+105 signal is similar to that from Cygnus X-1, but with self-similarity extending possibly to only a few days. This evidence of self-similarity provides a robust quantitative characterisation of anomalous transport occuring within the systems.
△ Less
Submitted 9 October, 2001; v1 submitted 4 July, 2001;
originally announced July 2001.
-
Shock Acceleration of Cosmic Rays - a critical review
Authors:
J. G. Kirk,
R. O. Dendy
Abstract:
Motivated by recent unsuccessful efforts to detect the predicted flux of TeV gamma-rays from supernova remnants, we present a critical examination of the theory on which these predictions are based. Three crucial problems are identified: injection, maximum achievable particle energy and spectral index. In each case significant new advances in understanding have been achieved, which cast doubt on…
▽ More
Motivated by recent unsuccessful efforts to detect the predicted flux of TeV gamma-rays from supernova remnants, we present a critical examination of the theory on which these predictions are based. Three crucial problems are identified: injection, maximum achievable particle energy and spectral index. In each case significant new advances in understanding have been achieved, which cast doubt on prevailing paradigms such as Bohm diffusion and single-fluid MHD. This indicates that more realistic analytical models, backed by more sophisticated numerical techniques should be employed to obtain reliable predictions. Preliminary work on incorporating the effects of anomalous transport suggest that the resulting spectrum should be significantly softer than that predicted by conventional theory.
△ Less
Submitted 11 January, 2001;
originally announced January 2001.
-
A sandpile model with tokamak-like enhanced confinement phenomenology
Authors:
S. C. Chapman,
R. O. Dendy,
B. Hnat
Abstract:
Confinement phenomenology characteristic of magnetically confined plasmas emerges naturally from a simple sandpile algorithm when the parameter controlling redistribution scalelength is varied. Close analogues are found for enhanced confinement, edge pedestals, and edge localised modes (ELMs), and for the qualitative correlations between them. These results suggest that tokamak observations of a…
▽ More
Confinement phenomenology characteristic of magnetically confined plasmas emerges naturally from a simple sandpile algorithm when the parameter controlling redistribution scalelength is varied. Close analogues are found for enhanced confinement, edge pedestals, and edge localised modes (ELMs), and for the qualitative correlations between them. These results suggest that tokamak observations of avalanching transport are deeply linked to the existence of enhanced confinement and ELMs.
△ Less
Submitted 2 February, 2001; v1 submitted 12 October, 2000;
originally announced October 2000.
-
Testing the SOC hypothesis for the magnetosphere
Authors:
N. W. Watkins,
M. P. Freeman,
S. C. Chapman,
R. O. Dendy
Abstract:
As noted by Chang, the hypothesis of Self-Organised Criticality provides a theoretical framework in which the low dimensionality seen in magnetospheric indices can be combined with the scaling seen in their power spectra and the recently-observed plasma bursty bulk flows. As such, it has considerable appeal, describing the aspects of the magnetospheric fuelling:storage:release cycle which are ge…
▽ More
As noted by Chang, the hypothesis of Self-Organised Criticality provides a theoretical framework in which the low dimensionality seen in magnetospheric indices can be combined with the scaling seen in their power spectra and the recently-observed plasma bursty bulk flows. As such, it has considerable appeal, describing the aspects of the magnetospheric fuelling:storage:release cycle which are generic to slowly-driven, interaction-dominated, thresholded systems rather than unique to the magnetosphere. In consequence, several recent numerical "sandpile" algorithms have been used with a view to comparison with magnetospheric observables. However, demonstration of SOC in the magnetosphere will require further work in the definition of a set of observable properties which are the unique "fingerprint" of SOC. This is because, for example, a scale-free power spectrum admits several possible explanations other than SOC.
A more subtle problem is important for both simulations and data analysis when dealing with multiscale and hence broadband phenomena such as SOC. This is that finite length systems such as the magnetosphere or magnetotail will by definition give information over a small range of orders of magnitude, and so scaling will tend to be narrowband. Here we develop a simple framework in which previous descriptions of magnetospheric dynamics can be described and contrasted. We then review existing observations which are indicative of SOC, and ask if they are sufficient to demonstrate it unambiguously, and if not, what new observations need to be made?
△ Less
Submitted 28 March, 2000; v1 submitted 17 February, 2000;
originally announced February 2000.
-
Electron acceleration due to high frequency instabilities at supernova remnant shocks
Authors:
M. E. Dieckmann,
K. G. McClements,
S. C. Chapman,
R. O. Dendy,
L. O'C. Drury
Abstract:
Observations of synchrotron radiation across a wide range of wavelengths provide clear evidence that electrons are accelerated to relativistic energies in supernova remnants (SNRs). However, a viable mechanism for the pre-acceleration of such electrons to mildly relativistic energies has not yet been established. In this paper an electromagnetic particle-in-cell (PIC) code is used to simulate ac…
▽ More
Observations of synchrotron radiation across a wide range of wavelengths provide clear evidence that electrons are accelerated to relativistic energies in supernova remnants (SNRs). However, a viable mechanism for the pre-acceleration of such electrons to mildly relativistic energies has not yet been established. In this paper an electromagnetic particle-in-cell (PIC) code is used to simulate acceleration of electrons from background energies to tens of keV at perpendicular collisionless shocks associated with SNRs. Free energy for electron energization is provided by ions reflected from the shock front, with speeds greater than the upstream electron thermal speed. The PIC simulation results contain several new features, including: the acceleration, rather than heating, of electrons via the Buneman instability; the acceleration of electrons to speeds exceeding those of the shock-reflected ions producing the instability; and strong acceleration of electrons perpendicular to the magnetic field. Electron energization takes place through a variety of resonant and non-resonant processes, of which the strongest involves stochastic wave-particle interactions. In SNRs the diffusive shock process could then supply the final step required for the production of fully relativistic electrons. The mechanisms identified in this paper thus provide a possible solution to the electron pre-acceleration problem.
△ Less
Submitted 17 February, 2000;
originally announced February 2000.
-
On the role of self-organised criticality in accretion systems
Authors:
R. O. Dendy,
P. Helander,
M. Tagger
Abstract:
Self-organised criticality (SOC) has been suggested as a potentially powerful unifying paradigm for interpreting the structure of, and signals from, accretion systems. After reviewing the most promising sites where SOC might be observable, we consider the theoretical arguments for supposing that SOC can occur in accretion discs. Perhaps the most rigorous evidence is provided by numerical modelli…
▽ More
Self-organised criticality (SOC) has been suggested as a potentially powerful unifying paradigm for interpreting the structure of, and signals from, accretion systems. After reviewing the most promising sites where SOC might be observable, we consider the theoretical arguments for supposing that SOC can occur in accretion discs. Perhaps the most rigorous evidence is provided by numerical modelling of energy dissipation due to magnetohydrodynamic turbulence in accretion discs by G Geertsema & A Achterberg (A&A {\bf 255}, 427 (1992)); we investigate how "sandpile"-type dynamics arise in this model. It is concluded that the potential sites for SOC in accretion systems are numerous and observationally accessible, and that theoretical support for the possible occurrence of SOC can be derived from first principles.
△ Less
Submitted 5 July, 1999;
originally announced July 1999.
-
Anomalous Transport and Particle Acceleration at Shocks
Authors:
P. Duffy,
J. G. Kirk,
Y. A. Gallant,
R. O. Dendy
Abstract:
The theory of first order Fermi acceleration at shocks assumes that particles diffuse due to scattering off slow-moving magnetic irregularities. However, cosmic rays are closely tied to magnetic field lines, and the transport process, particularly across the direction of the field, is likely to be more complicated. To describe cross field transport we employ recent extensions of the Rechester-Ro…
▽ More
The theory of first order Fermi acceleration at shocks assumes that particles diffuse due to scattering off slow-moving magnetic irregularities. However, cosmic rays are closely tied to magnetic field lines, and the transport process, particularly across the direction of the field, is likely to be more complicated. To describe cross field transport we employ recent extensions of the Rechester-Rosenbluth theory in localised stochastic regions of magnetic field. During acceleration at a shock and when the motion along the field is diffusive, there is a transition at a critical energy from \lq\lq sub-diffusive\rq\rq\ motion, where the mean square displacement of a particle increases with time as $t^{1/2}$, to compound diffusion, a combined process involving diffusion along a magnetic field which is itself wandering. Requiring this critical energy to be less than the system cut-off in a SNR of radius $R$ places an upper limit on the coherence length, $λ$, of the magnetic field for which diffusive shock acceleration can occur at quasi-perpendicular shocks when field line wandering dominates the effective transport of particles: $λ<R(U/v)^{1/2}$ where $U$ is the shock speed and $v$ the particle speed. The acceleration rate for particles in the sub-diffusive regime is derived and its implication for the origin of high energy cosmic rays is discussed.
△ Less
Submitted 11 September, 1995;
originally announced September 1995.