Shock location and CME 3D reconstruction of a solar type II radio burst with LOFAR

Context. Type II radio bursts are evidence of shocks in the solar atmosphere and inner heliosphere that emit radio waves ranging from sub-meter to kilometer lengths. These shocks may be associated with coronal mass ejections (CMEs) and reach speeds higher than the local magnetosonic speed. Radio imaging of decameter wavelengths (20–90 MHz) is now possible with the Low Frequency Array (LOFAR), opening a new radio window in which to study coronal shocks that leave the inner solar corona and enter the interplanetary medium and to understand their association with CMEs. Aims. To this end, we study a coronal shock associated with a CME and type II radio burst to determine the locations at which the radio emission is generated, and we investigate the origin of the band-splitting phenomenon. Methods. Thetype II shock source-positions and spectra were obtained using 91 simultaneous tied-array beams of LOFAR, and the CME was observed by the Large Angle and Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory (SOHO) and by the COR2A coronagraph of the SECCHI instruments on board the Solar Terrestrial Relation Observatory(STEREO). The 3D structure was inferred using triangulation of the coronographic observations. Coronal magnetic fields were obtained from a 3D magnetohydrodynamics (MHD) polytropic model using the photospheric fields measured by the Heliospheric Imager (HMI) on board the Solar Dynamic Observatory (SDO) as lower boundary. Results. The type II radio source of the coronal shock observed between 50 and 70 MHz was found to be located at the expanding flank of the CME, where the shock geometry is quasi-perpendicular with θBn ~ 70°. The type II radio burst showed first and second harmonic emission; the second harmonic source was cospatial with the first harmonic source to within the observational uncertainty. This suggests that radio wave propagation does not alter the apparent location of the harmonic source. The sources of the two split bands were also found to be cospatial within the observational uncertainty, in agreement with the interpretation that split bands are simultaneous radio emission from upstream and downstream of the shock front. The fast magnetosonic Mach number derived from this interpretation was found to lie in the range 1.3–1.5. The fast magnetosonic Mach numbers derived from modelling the CME and the coronal magnetic field around the type II source were found to lie in the range 1.4–1.6

Cassiopeia A, Cygnus A, Taurus A, and Virgo A at ultra-low radio frequencies

Context. The four persistent radio sources in the northern sky with the highest flux density at metre wavelengths are Cassiopeia A, Cygnus A, Taurus A, and Virgo A; collectively they are called the A-team. Their flux densities at ultra-low frequencies (< 100 MHz) can reach several thousands of janskys, and they often contaminate observations of the low-frequency sky by interfering with image processing. Furthermore, these sources are foreground objects for all-sky observations hampering the study of faint signals, such as the cosmological 21 cm line from the epoch of reionisation. Aims. We aim to produce robust models for the surface brightness emission as a function of frequency for the A-team sources at ultra-low frequencies. These models are needed for the calibration and imaging of wide-area surveys of the sky with low-frequency interferometers. This requires obtaining images at an angular resolution better than 15\u2033 with a high dynamic range and good image fidelity. Methods. We observed the A-team with the Low Frequency Array (LOFAR) at frequencies between 30 MHz and 77 MHz using the Low Band Antenna system. We reduced the datasets and obtained an image for each A-team source. Results. The paper presents the best models to date for the sources Cassiopeia A, Cygnus A, Taurus A, and Virgo A between 30 MHz and 77 MHz. We were able to obtain the aimed resolution and dynamic range in all cases. Owing to its compactness and complexity, observations with the long baselines of the International LOFAR Telescope will be required to improve the source model for Cygnus A further

Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz. However, it is technically and logistically challenging to process LOFAR data at this resolution. To date only a handful of publications have exploited this capability. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. It is implemented in a pipeline using mostly dedicated LOFAR software tools and the same processing framework as the LOFAR Two-metre Sky Survey (LoTSS). We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LoTSS pointing with an 8 h observation and 13 international stations. We perform in-field delay calibration, solution referencing to other calibrators in the field, self-calibration of these calibrators, and imaging of example directions of interest in the field. We find that for this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5° away, while phase solution transferral works well over ~1°. We also demonstrate a check of the astrometry and flux density scale with the in-field delay calibrator source. Imaging in 17 directions, we find the restoring beam is typically ~0.3′′ ×0.2′′ although this varies slightly over the entire 5 deg2 field of view. We find we can achieve ~80–300 μJy bm−1 image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 μJy bm−1 for the 8 h observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image roughly 900 sources per LoTSS pointing. This equates to ~ 3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution makes this estimate a lower limit

The Role of Radioactivities in Astrophysics

I present both a history of radioactivity in astrophysics and an introduction to the major applications of radioactive abundances to astronomy

On the mechanisms governing gas penetration into a tokamak plasma during a massive gas injection

A new 1D radial fluid code, IMAGINE, is used to simulate the penetration of gas into a tokamak plasma during a massive gas injection (MGI). The main result is that the gas is in general strongly braked as it reaches the plasma, due to mechanisms related to charge exchange and (to a smaller extent) recombination. As a result, only a fraction of the gas penetrates into the plasma. Also, a shock wave is created in the gas which propagates away from the plasma, braking and compressing the incoming gas. Simulation results are quantitatively consistent, at least in terms of orders of magnitude, with experimental data for a D 2 MGI into a JET Ohmic plasma. Simulations of MGI into the background plasma surrounding a runaway electron beam show that if the background electron density is too high, the gas may not penetrate, suggesting a possible explanation for the recent results of Reux et al in JET (2015 Nucl. Fusion 55 093013)

Velocity-space sensitivity of the time-of-flight neutron spectrometer at JET

The velocity-space sensitivities of fast-ion diagnostics are often described by so-called weight functions. Recently, we formulated weight functions showing the velocity-space sensitivity of the often dominant beam-target part of neutron energy spectra. These weight functions for neutron emission spectrometry (NES) are independent of the particular NES diagnostic. Here we apply these NES weight functions to the time-of-flight spectrometer TOFOR at JET. By taking the instrumental response function of TOFOR into account, we calculate time-of-flight NES weight functions that enable us to directly determine the velocity-space sensitivity of a given part of a measured time-of-flight spectrum from TOFOR

Relationship of edge localized mode burst times with divertor flux loop signal phase in JET

A phase relationship is identified between sequential edge localized modes (ELMs) occurrence times in a set of H-mode tokamak plasmas to the voltage measured in full flux azimuthal loops in the divertor region. We focus on plasmas in the Joint European Torus where a steady H-mode is sustained over several seconds, during which ELMs are observed in the Be II emission at the divertor. The ELMs analysed arise from intrinsic ELMing, in that there is no deliberate intent to control the ELMing process by external means. We use ELM timings derived from the Be II signal to perform direct time domain analysis of the full flux loop VLD2 and VLD3 signals, which provide a high cadence global measurement proportional to the voltage induced by changes in poloidal magnetic flux. Specifically, we examine how the time interval between pairs of successive ELMs is linked to the time-evolving phase of the full flux loop signals. Each ELM produces a clear early pulse in the full flux loop signals, whose peak time is used to condition our analysis. The arrival time of the following ELM, relative to this pulse, is found to fall into one of two categories: (i) prompt ELMs, which are directly paced by the initial response seen in the flux loop signals; and (ii) all other ELMs, which occur after the initial response of the full flux loop signals has decayed in amplitude. The times at which ELMs in category (ii) occur, relative to the first ELM of the pair, are clustered at times when the instantaneous phase of the full flux loop signal is close to its value at the time of the first ELM