BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

T. D. Arber

T. D. Arber contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRphysics.plasm-phastro-phphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2016arXiv

Simulations of Alfven wave driving of the solar chromosphere - efficient heating and spicule launching

Two of the central problems in our understanding of the solar chromosphere are how the upper chromosphere is heated and what drives spicules. Estimates of the required chromospheric heating, based on radiative and conductive losses suggest a rate of ${\sim} 0.1 \mathrm{\:erg\:cm^{-3}\:s^{-1}}$ in the lower chromosphere dropping to ${\sim} 10^{-3} \mathrm{\:erg\:cm^{-3}\:s^{-1}}$ in the upper chromosphere (\citet{Avrett1981}). The chromosphere is also permeated by spicules, higher density plasma from the lower atmosphere propelled upwards at speeds of ${\sim} 10-20 \mathrm{\:km\:s^{-1}}$, for so called Type-I spicules (\citet{Pereira2012,Zhang2012}), reaching heights of ${\sim} 3000-5000 \mathrm{\:km}$ above the photosphere. A clearer understanding of chromospheric dynamics, its heating and the formation of spicules, is thus of central importance to solar atmospheric science. For over thirty years it has been proposed that photospheric driving of MHD waves may be responsible for both heating and spicule formation. This paper presents results from a high-resolution MHD treatment of photospheric driven Alfvén and kink waves propagating upwards into an expanding flux tube embedded in a model chromospheric atmosphere. We show that the ponderomotive coupling from Alfvén and kink waves into slow modes generates shocks which both heat the upper chromosphere and drive spicules. These simulations show that wave driving of the solar chromosphere can give a local heating rate which matches observations and drive spicules consistent with Type-I observations all within a single coherent model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Alfvén Wave Heating of the Solar Chromosphere: 1.5D models

Physical processes which may lead to solar chromospheric heating are analyzed using high-resolution 1.5D non-ideal MHD modelling. We demonstrate that it is possible to heat the chromospheric plasma by direct resistive dissipation of high-frequency Alfvén waves through Pedersen resistivity. However this is unlikely to be sufficient to balance radiative and conductive losses unless unrealistic field strengths or photospheric velocities are used. The precise heating profile is determined by the input driving spectrum since in 1.5D there is no possibility of Alfvén wave turbulence. The inclusion of the Hall term does not affect the heating rates. If plasma compressibility is taken into account, shocks are produced through the ponderomotive coupling of Alfvén waves to slow modes and shock heating dominates the resistive dissipation. In 1.5D shock coalescence amplifies the effects of shocks and for compressible simulations with realistic driver spectra the heating rate exceeds that required to match radiative and conductive losses. Thus while the heating rates for these 1.5D simulations are an overestimate they do show that ponderomotive coupling of Alfvén waves to sound waves is more important in chromospheric heating than Pedersen dissipation through ion-neutral collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dense electron-positron plasmas and bursts of gamma-rays from laser-generated QED plasmas

In simulations of a 12.5PW laser (focused intensity I = 4x10^23W/cm^2) striking a solid aluminium target 10% of the laser energy is converted to gamma-rays. A dense electron-positron plasma is generated with a maximum density of 10^26/m^3; seven orders of magnitude denser than pure e-e+ plasmas generated with 1PW lasers. When the laser power is increased to 320PW (I = 10^25W/cm^2) 40% of the laser energy is converted to gamma-ray photons and 10% to electron-positron pairs. In both cases there is strong feedback between the QED emission processes and the plasma physics; the defining feature of the new `QED-plasma' regime reached in these interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Modelling Gamma Ray Emission and Pair Production in High-Intensity Laser-Matter Interactions

In high-intensity (> 10^21W/cm^2) laser-matter interactions gamma-ray photon emission by the electrons can strongly affect the electron's dynamics and copious numbers of electron-positron pairs can be produced by the emitted photons. We show how these processes can be included in simulations by coupling a Monte-Carlo algorithm describing the emission to a particle-in-cell code. The Monte-Carlo algorithm includes quantum corrections to the photon emission, which we show must be included if the pair production rate is to be correctly determined. The accuracy, convergence and energy conservation properties of the Monte-Carlo algorithm are analysed in simple test problems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simulation of the Observed Coronal Kink Instability and Its Implications for the SDO/AIA

Srivastava et al. (2010) have observed a highly twisted coronal loop, which was anchored in AR10960 during the period 04:43 UT-04:52 UT on 4 June 2007. The loop length and radius are approximately 80 Mm and 4 Mm, with a twist of 11.5 $π$. These observations are used as initial conditions in a three dimensional nonlinear magnetohydrodynamic simulation with parallel thermal conduction included. The initial unstable equilibrium evolves into the kink instability, from which synthetic observables are generated for various high-temperature filters of SDO/AIA. These observables include temporal and spatial averaging to account for the resolution and exposure times of SDO/AIA images. Using the simulation results, we describe the implications of coronal kink instability as observables in SDO/AIA filters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Synchrotron radiation, pair production and longitudinal electron motion during 10-100PW laser solid interactions

At laser intensities above 1023W/cm2 the interaction of a laser with a plasma is qualitatively different to the interactions at lower intensities. In this intensity regime solid targets start to become relativistically underdense, gamma-ray production by synchrotron emission starts to become an important feature of the dynamics and, at even higher intensities, electron-positron pair production by the non-linear Breit-Wheeler process starts to occur. Previous work in this intensity regime has considered ion acceleration1,2, identified different mechanisms for the underlying plasma physics of laser generation of gamma-rays3,4,5 considered the effect of target parameters on gamma-ray generation6 and considered the creation of solid density positronium plasma3. However a complete linked understanding of the important new physics of this regime is still lacking. In this paper, an analysis is presented of the effects of target density, laser intensity, target preplasma properties and other parameters on the conversion efficiency, spectrum and angular distribution of gamma-rays by synchrotron emission. An analysis of the importance of Breit-Wheeler pair production is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dense Electron-Positron Plasmas and Ultra-Intense Bursts of Gamma-Rays from Laser-Irradiated Solids

In simulations of a 10PW laser striking a solid we demonstrate the possibility of producing a pure electron-positron plasma by the same processes as those thought to operate in high-energy astrophysical environments. A maximum positron density of 10^26/m^3 is achieved, seven orders of magnitude greater than achieved in previous experiments. Additionally, 35% of the laser energy is converted to a burst of gamma-rays of intensity 10^22W/cm^2, potentially the most intense gamma-ray source available in the laboratory. This absorption results in a strong feedback between both pair and gamma-ray production and classical plasma physics in the new `QED-plasma' regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear wave propagation and reconnection at magnetic X-points in the Hall MHD regime

The highly dynamical, complex nature of the solar atmosphere naturally implies the presence of waves in a topologically varied magnetic environment. Here, the interaction of waves with topological features such as null points is inevitable and potentially important for energetics. The low resistivity of the solar coronal plasma implies that non-MHD effects should be considered in studies of magnetic energy release in this environment. This paper investigates the role of the Hall term in the propagation and dissipation of waves, their interaction with 2D magnetic X-points and the nature of the resulting reconnection. A Lagrangian remap shock-capturing code (Lare2d) is used to study the evolution of an initial fast magnetoacoustic wave annulus for a range of values of the ion skin depth in resistive Hall MHD. A magnetic null-point finding algorithm is also used to locate and track the evolution of the multiple null-points that are formed in the system. Depending on the ratio of ion skin depth to system size, our model demonstrates that Hall effects can play a key role in the wave-null interaction. In particular, the initial fast-wave pulse now consists of whistler and ion-cyclotron components; the dispersive nature of the whistler wave leads to (i) earlier interaction with the null, (ii) the creation of multiple additional, transient nulls and, hence, an increased number of energy release sites. In the Hall regime, the relevant timescales (such as the onset of reconnection and the period of the oscillatory relaxation) of the system are reduced significantly, and the reconnection rate is enhanced.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observational signatures of the coronal kink instability with thermal conduction

It is known from numerical simulations that thermal conduction along magnetic field lines plays an important role in the evolution of the kink instability in coronal loops. This study presents the observational signatures of the kink instability in long coronal loops when parallel thermal conduction is included. The 3D nonlinear magnetohydrodynamic equations are solved numerically to simulate the evolution of a coronal loop that is initially in an unstable equilibrium. The loop has length 80 Mm, width 8 Mm and an initial maximum twist of Phi = 11.5 pi, where Phi is a function of the radius. The initial loop parameters are obtained from a highly twisted loop observed in the TRACE 171 A waveband. Synthetic observables are generated from the data. These observables include spatial and temporal averaging to account for the resolution and exposure times of TRACE images. Parallel thermal conduction reduces the maximum local temperature by up to an order of magnitude. This means that different spectral lines are formed and different internal loop structures are visible with or without the inclusion of thermal conduction. However, the response functions sample a broad range of temperatures. The result is that the inclusion of parallel thermal conductivity does not have as large an impact on observational signatures as the order of magnitude reduction in the maximum temperature would suggest; the net effect is a blurring of internal features of the loop structure.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Emergence of a flux tube through a partially ionised solar atmosphere

For a magnetic flux tube, or indeed any flux, to emerge into the Solar corona from the convection zone it must pass through the partially ionised layers of the lower atmosphere: the photosphere and the chromosphere. In such regions the ion-neutral collisions lead to an increased resistivity for currents flowing across magnetic field lines. This Cowling resistivity can exceed the Spitzer resistivity by orders of magnitude and in 2.5D simulations has been shown to be sufficient to remove all cross field current from emerging flux. Here we extend this modelling into 3D. Once again it is found that the Cowling resistivity removes perpendicular current. However the presence of 3D structure prevents the simple comparison possible in 2.5D simulations. With a fully ionised atmosphere the flux emergence leads to an unphysically low temperature region in the overlying corona, lifting of chromospheric material and the subsequent onset of the Rayleigh-Taylor instability. Including neutrals removes the low temperature region, lifts less chromospheric matter and shows no signs of the Rayleigh-Taylor instability. Simulations of flux emergence therefore should include such a neutral layer in order to obtain the correct perpendicular current, remove the Rayleigh-Taylor instability and get the correct temperature profile. In situations when the temperature is not important, i.e. when no simulated spectral emission is required, a simple model for the neutral layer is demonstrated to adequately reproduce the results of fully consistent simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Aspects of Electron Acoustic Wave Physics in Laser Backscatter from Plasmas

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.

0 citations0 reviewsNo code linkedOpen access