Multi-particle azimuthal cumulants are measured as a function of centrality and transverse momentum using 470 mu b(-1) of Pb+Pb collisions at root s(NN) = 5.02TeV with the ATLAS detector at the LHC. These cumulants provide information on the event-by-event fluctuations of harmonic flow coefficients v(n) and correlated fluctuations between two harmonics v(n) and v(m). For the first time, a non-zero four-particle cumulant is observed for dipolar flow, v(1). The four-particle cumulants for elliptic flow, v(2), and triangular flow, v(3), exhibit a strong centrality dependence and change sign in ultra-central collisions. This sign change is consistent with significant non-Gaussian fluctuations in v(2) and v(3). The four-particle cumulant for quadrangular flow, v(4), is found to change sign in mid-central collisions. Correlations between two harmonics are studied with three- and four-particle mixed-harmonic cumulants, which indicate an anti-correlation between v(2) and v(3), and a positive correlation between v(2) and v(4). These correlations decrease in strength towards central collisions and either approach zero or change sign in ultra-central collisions. To investigate the possible flow fluctuations arising from intrinsic centrality or volume fluctuations, the results are compared between two different event classes used for centrality definitions. In peripheral and mid-central collisions where the cumulant signals are large, only small differences are observed. In ultra-central collisions, the differences are much larger and transverse momentum dependent. These results provide new information to disentangle flow fluctuations from the initial and final states, as well as new insights on the influence of centrality fluctuations.

The elliptic flow of muons from the decay of charm and bottom hadrons is measured in pp collisions at root s = 13 TeV using a data sample with an integrated luminosity of 150 pb(-1) recorded by the ATLAS detector at the LHC. The muons from heavy-flavor decay are separated from light-hadron decay muons using momentum imbalance between the tracking and muon spectrometers. The heavy-flavor decay muons are further separated into those from charm decay and those from bottom decay using the distance-of-closest-approach to the collision vertex. The measurement is performed for muons in the transverse momentum range 4-7 GeV and pseudorapidity range vertical bar eta vertical bar < 2.4. A significant nonzero elliptic anisotropy coefficient nu(2) is observed for muons from charm decays, while the nu(2) value for muons from bottom decays is consistent with zero within uncertainties.

A search for the electroweak production of charginos and sleptons decaying into final states with two electrons or muons is presented. The analysis is based on 139 fb(-1) of proton-proton collisions recorded by the ATLAS detector at the Large Hadron Collider at v s = 13 TeV. Three R-parity-conserving scenarios where the lightest neutralino is the lightest supersymmetric particle are considered: the production of chargino pairs with decays via eitherW bosons or sleptons, and the direct production of slepton pairs. The analysis is optimised for the first of these scenarios, but the results are also interpreted in the others. No significant deviations from the Standard Model expectations are observed and limits at 95% confidence level are set on the masses of relevant supersymmetric particles in each of the scenarios. For a massless lightest neutralino, masses up to 420 GeV are excluded for the production of the lightest-chargino pairs assuming W-boson-mediated decays and up to 1 TeV for slepton-mediated decays, whereas for slepton-pair production masses up to 700 GeV are excluded assuming three generations of mass-degenerate sleptons.

The production yield of Z bosons is measured in the electron and muon decay channels in Pb+Pb collisions at /S-NN = 5.02 TeV with the ATLAS detector. Data from the 2015 LHC run corresponding to an integrated luminosity of 0.49 nb(-1) are used for the analysis. The Z boson yield, normalised by the total number of minimum-bias events and the mean nuclear thickness function, is measured as a function of dilepton rapidity and event centrality. The measurements in Pb+Pb collisions are compared with similar measurements made in proton-proton collisions at the same centre-of-mass energy. The nuclear modification factor is found to be consistent with unity for all centrality intervals. The results are compared with theoretical predictions obtained at next-to-leading order using nucleon and nuclear parton distribution functions. The normalised Z boson yields in Pb+Pb collisions lie 1-3a above the predictions. The nuclear modification factor measured as a function of rapidity agrees with unity and is consistent with a next-to-leading-order QCD calculation including the isospin effect.

A search for the direct production of the supersymmetric partners of tau-leptons (staus) in final states with two hadronically decaying tau-leptons is presented. The analysis uses a dataset of pp collisions corresponding to an integrated luminosity of 139 fb(-1), recorded with the ATLAS detector at the Large Hadron Collider at a center-of-mass energy of 13 TeV. No significant deviation from the expected Standard Model background is observed. Limits are derived in scenarios of direct production of stau pairs with each stau decaying into the stable lightest neutralino and one tau-lepton in simplified models where the two stau mass eigenstates are degenerate. Stau masses from 120 GeV to 390 GeV are excluded at 95% confidence level for a massless lightest neutralino.

A search for the associated production of a heavy resonance with a top-quark or a top-antitop-quark pair, and decaying into a tt¯ pair is presented. The search uses the data recorded by the ATLAS detector in pp collisions at s=13 TeV at the Large Hadron Collider during the years 2015–2018, corresponding to an integrated luminosity of 139 fb-1. Events containing exactly one electron or muon are selected. The two hadronically decaying top quarks from the resonance decay are reconstructed using jets clustered with a large radius parameter of R=1. The invariant mass spectrum of the two top quark candidates is used to search for a resonance signal in the range of 1.0 TeV to 3.2 TeV. The presence of a signal is examined using an approach with minimal model dependence followed by a model-dependent interpretation. No significant excess is observed over the background expectation. Upper limits on the production cross section times branching ratio at 95% confidence level are provided for a heavy Z′ boson based on a simplified model, for Z′ mass between 1.0 TeV and 3.0 TeV. The observed (expected) limits range from 21 (14) fb to 119 (86) fb depending on the choice of model parameters.

The Milky Way is expected to be embedded in a halo of dark matter particles, with the highest density in the central region, and decreasing density with the halo-centric radius. Dark matter might be indirectly detectable at Earth through a flux of stable particles generated in dark matter annihilations and peaked in the direction of the Galactic Center. We present a search for an excess flux of muon (anti-) neutrinos from dark matter annihilation in the Galactic Center using the cubic-kilometer-sized IceCube neutrino detector at the South Pole. There, the Galactic Center is always seen above the horizon. Thus, new and dedicated veto techniques against atmospheric muons are required to make the southern hemisphere accessible for IceCube. We used 319.7 live-days of data from IceCube operating in its 79-string configuration during 2010 and 2011. No neutrino excess was found and the final result is compatible with the background. We present upper limits on the self-annihilation cross-section, < sAv >, for WIMP masses ranging from 30GeV up to 10TeV, assuming cuspy (NFW) and flat-cored (Burkert) dark matter halo profiles, reaching down to similar or equal to 4 . 10(-24) cm(3) s(-1), and similar or equal to 2.6 . 10(-23) cm(3) s(-1) for the nu(nu) over bar channel, respectively.

Fully nonlinear propagation of ion-acoustic solitary waves in a collisionless dense/quantum electron-positron-ion plasma is investigated. The electrons and positrons are assumed to follow the Thomas-Fermi density distribution and the ions are described by the hydrodynamic equations. An energy balance-like equation involving a Sagdeev-type pseudo-potential is derived. Finite amplitude solutions are obtained numerically and their characteristics are discussed. The small-but finite-amplitude limit is also considered and an exact analytical solution is obtained. The present studies might be helpful to understand the excitation of nonlinear ion-acoustic solitary waves in a degenerate plasma such as in superdense white dwarfs.

By using the Thomas-Fermi electron density distribution for quantum degenerate electrons, the hydrodynamic equations for ions, and the Poisson equation, planar and nonplanar ion-acoustic solitary waves in an unmagnetized collisionless plasma are investigated. The reductive perturbation method is used to derive cylindrical and spherical Korteweg-de Vries equations. Numerical solutions of the latter are presented. The present results can be useful in understanding the features of small but finite amplitude localized ion-acoustic solitary pulses in a degenerate plasma.

We present the fourth Fermi Large Area Telescope catalog (4FGL) of gamma-ray sources. Based on the first eight years of science data from the Fermi Gamma-ray Space Telescope mission in the energy range from 50 MeV to 1 TeV, it is the deepest yet in this energy range. Relative to the 3FGL catalog, the 4FGL catalog has twice as much exposure as well as a number of analysis improvements, including an updated model for the Galactic diffuse gamma-ray emission, and two sets of light curves (one-year and two-month intervals). The 4FGL catalog includes 5064 sources above 4 sigma significance, for which we provide localization and spectral properties. Seventy-five sources are modeled explicitly as spatially extended, and overall, 358 sources are considered as identified based on angular extent, periodicity, or correlated variability observed at other wavelengths. For 1336 sources, we have not found plausible counterparts at other wavelengths. More than 3130 of the identified or associated sources are active galaxies of the blazar class, and 239 are pulsars.

This document presents a summary of the 2023 Terrestrial Very-Long-Baseline Atom Interferometry Workshop hosted by CERN. The workshop brought together experts from around the world to discuss the exciting developments in large-scale atom interferometer (AI) prototypes and their potential for detecting ultralight dark matter and gravitational waves. The primary objective of the workshop was to lay the groundwork for an international TVLBAI proto-collaboration. This collaboration aims to unite researchers from different institutions to strategize and secure funding for terrestrial large-scale AI projects. The ultimate goal is to create a roadmap detailing the design and technology choices for one or more kilometer–scale detectors, which will be operational in the mid-2030s. The key sections of this report present the physics case and technical challenges, together with a comprehensive overview of the discussions at the workshop together with the main conclusions.

First proton production rates from the d(d,p)t reaction in the Mega Amp Spherical Tokamak Upgrade (MAST-U) are measured. The data were taken during the MAST-U experimental campaign with an upgraded version of the proton detector (PD) previously used in MAST. The new detector array consists of three collimated silicon surface barrier detectors with a depletion depth of 300 mu m and a collimated 120 mu m thick diamond detector, mounted on the MAST-U reciprocating probe arm. This array measures the energies of unconfined energetic 3 MeV protons and 1 MeV tritons mainly produced by beam-thermal DD reactions during neutral beam injection heating. Diamond detectors have the potential to be uniquely suited to detect charged fusion products as they promise to be much more radiation resistant and much less sensitive to temperature variations compared to silicon-based detectors. Using silicon and diamond-based detectors simultaneously allowed us to directly compare the performance of these two detector types. PD particle rates measured during different plasma scenarios are presented and compared to neutron rates measured using the neutron camera upgrade and TRANSP predictions.

We apply the collective coordinate model framework to describe collisions of a kink and an antikink with nonzero total momentum, i.e., when the solitons possess different velocities. The minimal moduli space with only two coordinates (the mutual distance and the position of the center of mass) is of a wormhole type, whose throat shrinks to a point for symmetric kinks. In this case, a singularity is formed. For non-zero momentum, it prohibits solutions where the solitons pass through each other. We show that this unphysical feature can be cured by enlarging the dimension of the moduli space, e.g., by the inclusion of internal modes.

The nonlinear propagation of small but finite amplitude dust ion-acoustic waves (DIAWs) in an ion-beam driven plasma consisting of Boltzmannian electrons, positive ions, and stationary negatively charged dust grains is studied by using the standard reductive perturbation technique. It is shown that there exist two critical values (γ_c1) and (γ_c2) of ion beam to ion phase velocity ratio (γ), above and below which the beam generated solitons are not possible. The effects of the parameters, namely, γ, the ratio of the ion beam to plasma ion density (μ_i), the dust to ion density ratio (μ_d), and the ion beam to plasma ion mass ratio (μ) on both the amplitude and width of the stationary DIAWs, are analyzed numerically, and applications of the results to laboratory ion beam as well as space plasmas (e.g., auroral plasmas) are explained.

The quadrature error associated with a regular quadrature rule for evaluation of a layer potential increases rapidly when the evaluation point approaches the surface and the integral becomes nearly singular. Error estimates are needed to determine when the accuracy is insufficient and a more costly special quadrature method should be utilized.& nbsp;The final result of this paper are such quadrature error estimates for the composite Gauss-Legendre rule and the global trapezoidal rule, when applied to evaluate layer potentials defined over smooth curved surfaces in R-3. The estimates have no unknown coefficients and can be efficiently evaluated given the discretization of the surface, invoking a local one-dimensional root-finding procedure. They are derived starting with integrals over curves, using complex analysis involving contour integrals, residue calculus and branch cuts. By complexifying the parameter plane, the theory can be used to derive estimates also for curves in R3. These results are then used in the derivation of the estimates for integrals over surfaces. In this procedure, we also obtain error estimates for layer potentials evaluated over curves in R2. Such estimates combined with a local root-finding procedure for their evaluation were earlier derived for the composite Gauss-Legendre rule for layer potentials written in complex form [4]. This is here extended to provide quadrature error estimates for both complex and real formulations of layer potentials, both for the Gauss-Legendre and the trapezoidal rule.& nbsp;Numerical examples are given to illustrate the performance of the quadrature error estimates. The estimates for integration over curves are in many cases remarkably precise, and the estimates for curved surfaces in R-3 are also sufficiently precise, with sufficiently low computational cost, to be practically useful.

The quadrature error associated with a regular quadrature rule for evaluation of a layer potential increases rapidly when the evaluation point approaches the surface and the integral becomes nearly singular. Error estimates are needed to determine when the accuracy is insufficient and a more costly special quadrature method should be utilized.& nbsp;The final result of this paper are such quadrature error estimates for the composite Gauss-Legendre rule and the global trapezoidal rule, when applied to evaluate layer potentials defined over smooth curved surfaces in R-3. The estimates have no unknown coefficients and can be efficiently evaluated given the discretization of the surface, invoking a local one-dimensional root-finding procedure. They are derived starting with integrals over curves, using complex analysis involving contour integrals, residue calculus and branch cuts. By complexifying the parameter plane, the theory can be used to derive estimates also for curves in R3. These results are then used in the derivation of the estimates for integrals over surfaces. In this procedure, we also obtain error estimates for layer potentials evaluated over curves in R2. Such estimates combined with a local root-finding procedure for their evaluation were earlier derived for the composite Gauss-Legendre rule for layer potentials written in complex form [4]. This is here extended to provide quadrature error estimates for both complex and real formulations of layer potentials, both for the Gauss-Legendre and the trapezoidal rule.& nbsp;Numerical examples are given to illustrate the performance of the quadrature error estimates. The estimates for integration over curves are in many cases remarkably precise, and the estimates for curved surfaces in R-3 are also sufficiently precise, with sufficiently low computational cost, to be practically useful.

We design and develop a new Particle-in-Cell (PIC) method for plasma simulations using Deep-Learning (DL) to calculate the electric field from the electron phase space. We train a Multilayer Perceptron (MLP) and a Convolutional Neural Network (CNN) to solve the two-stream instability test. We verify that the DL-based MLP PIC method produces the correct results using the two-stream instability: the DL-based PIC provides the expected growth rate of the two-stream instability. The DL-based PIC does not conserve the total energy and momentum. However, the DL-based PIC method is stable against the cold-beam instability, affecting traditional PIC methods. This work shows that integrating DL technologies into traditional computational methods is a viable approach for developing next-generation PIC algorithms.

We report on the temporally and spatially resolved detection of the precursory stages that lead to the formation of an unmagnetized, supercritical collisionless shock in a laser-driven laboratory experiment. The measured evolution of the electrostatic potential associated with the shock unveils the transition from a current free double layer into a symmetric shock structure, stabilized by ion reflection at the shock front. Supported by a matching particle-in-cell simulation and theoretical considerations, we suggest that this process is analogous to ion reflection at supercritical collisionless shocks in supernova remnants.

We report on the experimental observation of the instability of a plasma shell, which formed during the expansion of a laser-ablated plasma into a rarefied ambient medium. By means of a proton radiography technique, the evolution of the instability is temporally and spatially resolved on a timescale much shorter than the hydrodynamic one. The density of the thin shell exceeds that of the surrounding plasma, which lets electrons diffuse outward. An ambipolar electric field grows on both sides of the thin shell that is antiparallel to the density gradient. Ripples in the thin shell result in a spatially varying balance between the thermal pressure force mediated by this field and the ram pressure force that is exerted on it by the inflowing plasma. This mismatch amplifies the ripples by the same mechanism that drives the hydrodynamic nonlinear thin-shell instability (NTSI). Our results thus constitute the first experimental verification that the NTSI can develop in colliding flows.

The divertor solutions obtained with the plasma edge modelling tool SOLPS5.0 are discussed. The code results are benchmarked against carefully analysed L-mode discharges at various density levels with and without impurity seeding in the full-metal tokamaks ASDEX Upgrade and JET. The role of the cross-field drifts and currents in the solutions is analysed in detail, and the improvements achieved by fully activating the drift and current terms in view of matching the experimental signals are addressed. The persisting discrepancies are also discussed.

The divertor solutions obtained with the plasma edge modelling tool SOLPS5.0 are discussed. The code results are benchmarked against carefully analysed L-mode discharges at various density levels with and without impurity seeding in the full-metal tokamaks ASDEX Upgrade and JET. The role of the cross-field drifts and currents in the solutions is analysed in detail, and the improvements achieved by fully activating the drift and current terms in view of matching the experimental signals are addressed. The persisting discrepancies are also discussed.

The stability with respect to a peeling-ballooning mode (PBM) was investigated numerically with extended MHD simulation codes in JET, JT-60U and future JT-60SA plasmas. The MINERVA-DI code was used to analyze the linear stability, including the effects of rotation and ion diamagnetic drift (omega(*i)), in JET-ILW and JT-60SA plasmas, and the JOREK code was used to simulate nonlinear dynamics with rotation, viscosity and resistivity in JT-60U plasmas. It was validated quantitatively that the ELM trigger condition in JET-ILW plasmas can be reasonably explained by taking into account both the rotation and omega(*i) effects in the numerical analysis. When deuterium poloidal rotation is evaluated based on neoclassical theory, an increase in the effective charge of plasma destabilizes the PBM because of an acceleration of rotation and a decrease in omega(*i). The difference in the amount of ELM energy loss in JT-60U plasmas rotating in opposite directions was reproduced qualitatively with JOREK. By comparing the ELM affected areas with linear eigenfunctions, it was confirmed that the difference in the linear stability property, due not to the rotation direction but to the plasma density profile, is thought to be responsible for changing the ELM energy loss just after the ELM crash. A predictive study to determine the pedestal profiles in JT-60SA was performed by updating the EPED1 model to include the rotation and w*i effects in the PBM stability analysis. It was shown that the plasma rotation predicted with the neoclassical toroidal viscosity degrades the pedestal performance by about 10% by destabilizing the PBM, but the pressure pedestal height will be high enough to achieve the target parameters required for the ITER-like shape inductive scenario in JT-60SA.

Stability to the type-I edge localized mode (ELM) in JET plasmas was investigated numerically by analyzing the stability to a peeling-ballooning mode with the effects of plasma rotation and ion diamagnetic drift. The numerical analysis was performed by solving the extended Frieman-Rotenberg equation with the MINERVA-DI code. To take into account these effects in the stability analysis self-consistently, the procedure of JET equilibrium reconstruction was updated to include the profiles of ion temperature and toroidal rotation, which are determined based on the measurement data in experiments. With the new procedure and MINERVA-DI, it was identified that the stability analysis including the rotation effect can explain the ELM trigger condition in JET with ITER like wall (JET-ILW), though the stability in JET with carbon wall (JET-C) is hardly affected by rotation. The key difference is that the rotation shear in JET-ILW plasmas analyzed in this study is larger than that in JET-C ones, the shear which enhances the dynamic pressure destabilizing a peeling-ballooning mode. In addition, the increase of the toroidal mode number of the unstable MHD mode determining the ELM trigger condition is also important when the plasma density is high in JET-ILW. Though such modes with high toroidal mode number are strongly stabilized by the ion diamagnetic drift effect, it was found that plasma rotation can sometimes overcome this stabilizing effect and destabilizes the peeling-ballooning modes in JET-ILW.

Stability to the type-I edge localized mode (ELM) in JET plasmas was investigated numerically by analyzing the stability to a peeling-ballooning mode with the effects of plasma rotation and ion diamagnetic drift. The numerical analysis was performed by solving the extended Frieman-Rotenberg equation with the MINERVA-DI code. To take into account these effects in the stability analysis self-consistently, the procedure of JET equilibrium reconstruction was updated to include the profiles of ion temperature and toroidal rotation, which are determined based on the measurement data in experiments. With the new procedure and MINERVA-DI, it was identified that the stability analysis including the rotation effect can explain the ELM trigger condition in JET with ITER like wall (JET-ILW), though the stability in JET with carbon wall (JET-C) is hardly affected by rotation. The key difference is that the rotation shear in JET-ILW plasmas analyzed in this study is larger than that in JET-C ones, the shear which enhances the dynamic pressure destabilizing a peeling-ballooning mode. In addition, the increase of the toroidal mode number of the unstable MHD mode determining the ELM trigger condition is also important when the plasma density is high in JET-ILW. Though such modes with high toroidal mode number are strongly stabilized by the ion diamagnetic drift effect, it was found that plasma rotation can sometimes overcome this stabilizing effect and destabilizes the peeling-ballooning modes in JET-ILW.

The stability with respect to a peeling-ballooning mode (PBM) was investigated numerically with extended MHD simulation codes in JET, JT-60U and future JT-60SA plasmas. The MINERVA-DI code was used to analyze the linear stability, including the effects of rotation and ion diamagnetic drift (omega(*i)), in JET-ILW and JT-60SA plasmas, and the JOREK code was used to simulate nonlinear dynamics with rotation, viscosity and resistivity in JT-60U plasmas. It was validated quantitatively that the ELM trigger condition in JET-ILW plasmas can be reasonably explained by taking into account both the rotation and omega(*i) effects in the numerical analysis. When deuterium poloidal rotation is evaluated based on neoclassical theory, an increase in the effective charge of plasma destabilizes the PBM because of an acceleration of rotation and a decrease in omega(*i). The difference in the amount of ELM energy loss in JT-60U plasmas rotating in opposite directions was reproduced qualitatively with JOREK. By comparing the ELM affected areas with linear eigenfunctions, it was confirmed that the difference in the linear stability property, due not to the rotation direction but to the plasma density profile, is thought to be responsible for changing the ELM energy loss just after the ELM crash. A predictive study to determine the pedestal profiles in JT-60SA was performed by updating the EPED1 model to include the rotation and w*i effects in the PBM stability analysis. It was shown that the plasma rotation predicted with the neoclassical toroidal viscosity degrades the pedestal performance by about 10% by destabilizing the PBM, but the pressure pedestal height will be high enough to achieve the target parameters required for the ITER-like shape inductive scenario in JT-60SA.

In order to get high deposition rate and good film properties, the stabilization of the transition zone between the metallic and compound modes is beneficial. We have shown earlier that at least in some cases, HiPIMS can reduce hysteresis effect in reactive sputtering. In our previous work, mechanisms for the suppression/elimination of the hysteresis effect have been suggested. Reactive HiPIMS can suppress/eliminate the hysteresis effect in the range of optimum frequency [1] lead to the process stability during the deposition with high deposition rate. The mechanisms behind this optimum frequency may relate with high erosion rate during the pulse [2,3] and gas rarefaction effect in front of the target [4]. 

In this contribution, reactive sputtering process using high power impulse magnetron sputtering (HiPIMS) has been studied with focus on the gas rarefaction. Through variations in the sputtering conditions such as pulse frequencies, peak powers, and target area, their effect on the shape of current waveforms have been analyzed. The current waveforms in compound mode are strongly affected. Our experiments show that the shape and amplitude of peak current cannot be explained by the change of the secondary electron yield due to target oxidation only. Reduced rarefaction in compound mode contributes to the observed very high peak current values.

A strategy that facilitates a substantial increase of carbon ionization in magnetron sputtering discharges is presented in this work. The strategy is based on increasing the electron temperature in a high power impulse magnetron sputtering discharge by using Ne as the sputtering gas. This allows for the generation of an energetic C+ ion population and a substantial increase in the C+ ion flux as compared to a conventional Ar-HiPIMS process. A direct consequence of the ionization enhancement is demonstrated by an increase in the mass density of the grown films up to 2.8 g/cm(3); the density values achieved are substantially higher than those obtained from conventional magnetron sputtering methods.

Intense fluctuations in the electric field at high altitudes in the auroral zone are frequently measured by the Viking satellite. We have made an analysis of the origin of electric and magnetic fluctuations in the frequency range of 0.1 - 1 Hz by assuming four different sources for the signals: (I) spatial structures, (2) spatial structures with a parallel potential drop below the satellite, (3) traveling; shear Alfven waves, and (4) interfering shear Alfven waves. We will shaw that these different sources of the signals may produce similar amplitude ratios and phase differences between the perpendicular electric and magnetic fields. Since the different sources have different frequency dependencies, this can be used as an additional test if the signals are broadband. In other cases, additional information is needed, for example, satellite particle measurements or ground; magnetic measurements. The ideas presented in the theory were tested for one Viking eveningside pass over Scandinavia, where ground-based magnetometer and EISCAT radar measurements were available. The magnetic conditions were active during this pass and several interfering shear Alfven waves were found. Also, a spatial structure with a parallel potential drop below the satellite was identified. The magnitude of the 10-km-wide potential drop was at least 2 kV and the upward field-aligned current 26 mu A m(-2) (value mapped to the ionospheric level). The held-aligned conductance was estimated as 1.3 - 2.2x10(-8) S m(-2).

We present Swarm satellite and EISCAT radar observations of electrodynamical parameters in the midnight sector at high latitudes. The most striking feature is a plasma flow channel located equatorward of the polar cap boundary within the dawn convection cell. The flow channel is 1.5 degrees wide in latitude and contains southward electric field of 150 mV/m, corresponding to eastward plasma velocities of 3,300 m/s in the F-region ionosphere. The theoretically computed ion temperature enhancement produced by the observed ion velocity is in accordance with the measured one by the EISCAT radar. The total width of the auroral oval is about 10 degrees in latitude. While the poleward part is electric field dominant with low conductivity and the flow channel, the equatorward part is conductivity dominant with at least five auroral arcs. The main part of the westward electrojet flows in the conductivity dominant part, but it extends to the electric field dominant part. According to Kamide and Kokubun (1996), the whole midnight sector westward electrojet is expected to be conductivity dominant, so the studied event challenges the traditional view. The flow channel is observed after substorm onset. We suggest that the observed flow channel, which is associated with a 13-kV horizontal potential difference, accommodates increased nightside plasma flows during the substorm expansion phase as a result of reconnection in the near-Earth magnetotail.

The dynamics of the polar cap boundary and auroral oval in the nightside ionosphere are studied during late expansion and recovery of a substorm from the region between Tromsø (66.6 degree cgmLat) and Longyearbyen (75.2 degree cgmLat) on 27 February 2004 by using the coordinated EISCAT incoherent scatter radar, MIRACLE magnetometer and Cluster satellite measurements. During the late substorm expansion/early recovery phase, the polar cap boundary (PCB) made zig-zag-type motion with amplitude of 2.5 degree cgmLat and period of about 30 min near magnetic midnight. We suggest that the poleward motions of the PCB were produced by bursts of enhanced reconnection at the near-Earth neutral line (NENL). The subsequent equatorward motions of the PCB would then represent the recovery of the merging line towards the equilibrium state (Cowley and Lockwood, 1992). The observed bursts of enhanced westward electrojet just equatorward of the polar cap boundary during poleward expansions were produced plausibly by particles accelerated in the vicinity of the neutral line and thus lend evidence to the Cowley-Lockwood paradigm.

During the substorm recovery phase, the footpoints of the Cluster satellites at a geocentric distance of 4.4R_E mapped in the vicinity of EISCAT measurements. Cluster data indicate that outflow of H⁺ and O⁺ ions took place within the plasma sheet boundary layer (PSBL) as noted in some earlier studies as well. We show that in this case the PSBL corresponded to a region of enhanced electron temperature in the ionospheric F region. It is suggested that the ion outflow originates from the F region as a result of increased ambipolar diffusion. At higher altitudes, the ions could be further energized by waves, which at Cluster altitudes were observed as BBELF (broad band extra low frequency) fluctuations.

The four-satellite configuration of Cluster revealed a sudden poleward expansion of the PSBL by 2 degree during ˜ 5 min. The beginning of the poleward motion of the PCB was associated with an intensification of the downward FAC at the boundary. We suggest that the downward FAC sheet at the PCB is the high-altitude counterpart of the Earthward flowing FAC produced in the vicinity of the magnetotail neutral line by the Hall effect (Sonnerup, 1979) during a short-lived reconnection pulse.

The polar cap boundary (PCB) location and motion in the nightside ionosphere has been studied by using measurements from the EISCAT radars and the MIRACLE magnetometers during a period of four substorms on 18 February 2004. The OMNI database has been used for observations of the solar wind and the Geotail satellite for magnetospheric measurements. In addition, the event was modelled by the GUMICS-4 MHD simulation. The simulation of the PCB location was in a rather good agreement with the experimental estimates at the EISCAT longitude. During the first three substorm expansion phases, neither the local observations nor the global simulation showed any poleward motions of the PCB, even though the electrojets intensified. Rapid poleward motions of the PCB took place only in the early recovery phases of the substorms. Hence, in these cases the nightside reconnection rate was locally higher in the recovery phase than in the expansion phase.

In addition, we suggest that the IMF Bz component correlated with the nightside tail inclination angle and the PCB location with about a 17-min delay from the bow shock. By taking the delay into account, the IMF northward turnings were associated with dipolarizations of the magnetotail and poleward motions of the PCB in the recovery phase. The mechanism behind this effect should be studied further.

The 2004 tracer experiment of JET with the injection of (CH4)-C-13 into H-mode plasma at the outer divertor has been modelled with the Monte Carlo impurity transport code ERO. EDGE2D solutions for inter-ELM and ELM-peak phases were used as plasma backgrounds. Local two-dimensional (2D) deposition patterns at the vertical outer divertor target plate were obtained for comparison with post-mortem surface analyses. ERO also provides emission profiles for comparison with radially resolved spectroscopic measurements. Modelling indicates that enhanced re-erosion of deposited carbon layers is essential in explaining the amount of local deposition. Assuming negligible effective sticking of hydrocarbons, the measured local deposition of 20-34% is reproduced if re-erosion of deposits is enhanced by a factor of 2.5-7 compared to graphite erosion. If deposits are treated like the substrate, the modelled deposition is 55%. Deposition measurements at the shadowed area around injectors can be well explained by assuming negligible re-erosion but similar sticking behaviour there as on plasma-wetted surfaces.

Migration of beryllium into the divertor and deposition on tungsten in the final phase of the first ITER-like-wall campaign of JET are modelled with the 3D Monte Carlo impurity transport code ERO. The simulation covers the inner wall and the inner divertor. To generate the plasma background for Monte Carlo tracing of impurity particles, we use the EDGE2D/EIRENE code set. At the relevant regions of the wall, the estimated plasma conditions vary around T-e approximate to 5eV and n(e) 2 x 10(17) m(-3) (far-scrape-off layer; more than 10 cm away from the LCFS). We calculate impurity distributions in the plasma using the main chamber source as a free parameter in modelling and attempt to reproduce inter-ELM spectroscopic BeII line (527 nm) profiles at the divertor. The present model reproduces the level of emission close to the inner wall, but further work is needed to match also the measured emission peak values and ultimately link the modelled poloidal net deposition profiles of beryllium to post mortem data.

We present the first comparison of multiple global simulations of the solar wind interaction with Mercury's dayside magnetosphere, conducted in the framework of the international collaborative project SHOTS - Studies on Hermean magnetosphere Oriented Theories and Simulations. Two 3D magnetohydrodynamic and two 3D hybrid simulation codes are used to investigate the global response of the Hermean magnetosphere without its exosphere to a northward-oriented interplanetary magnetic field. We cross-compare the results of the four codes for a theoretical case and a MESSENGER orbit with similar upstream plasma conditions. The models agree on bowshock and magnetopause locations at 2.1 ​± ​0.11 and 1.4 ​± ​0.1 Mercury planetary radii, respectively. The latter locations may be influenced by subtle differences in the treatment of the plasma boundary at the planetary surface. The predicted magnetosheath thickness varies less between the codes. Finally, we also sample the plasma data along virtual trajectories of BepiColombo's Magnetospheric and Planetary Orbiter. Our ability to accurately predict the structure of the Hermean magnetosphere aids the analysis of the onboard plasma measurements of past and future magnetospheric missions.

Results obtained in the Martian ionosphere by the Langmuir Probe and Waves instrument aboard the Mars Atmosphere and Volatile EvolutioN Mission spacecraft are presented. The results include ionospheric electron densities determined from the frequency of Langmuir waves. Since the amplitude of thermal Langmuir waves is often below the instrument's detection level, Langmuir Probe and Waves excites these waves by injecting into the plasma a 3.3-V white noise signal. Electric field spectral measurements obtained shortly after the excitation show a resonance line at frequencies slightly below the local plasma frequency. The observed resonance line is interpreted to originate from plasma waves excited in the wake behind the spacecraft. These results reveal an important phenomenon in electron density estimation from stimulated Langmuir waves. The observed phenomenon, not previously reported by earlier missions, may be a common process in active sounding that can affect in situ electron density measurements.

The purpose of this project is to investigate the spectral properties of turbulence in the magnetosheath of Saturn, using in-situ magnetic field measurements from the Cassini spacecraft. According to models of incompressible, turbulent fluids, the energy spectrum in the inertial range scales as the frequency to the power of -5/3, which has been observed in the near-Earth Solar wind but not in the Terrestrial magnetosheath unless close to the magnetopause. 120 time intervals for when Cassini is inside the magnetosheath are identified — 40 in each category of behind quasi-perpendicular bow shocks, behind quasi-parallel bow shocks, and inside the middle of the magnetosheath. The power spectral density is thereafter calculated for each interval, with logarithmic regressions performed at the MHD and sub-ion scales separated by the ion gyrofrequency. The results seem to indicate similar behaviour as in the magnetosheath of Earth, without significant difference between quasi-perpendicular and quasi-parallel cases except somewhat steeper exponents at the MHD scale for the former. These observations confirm the role of the bow shock in destroying the fully developed turbulence of the Solar wind, thus explaining the absence of the inertial range.

The transmission of a sheath region driven by an interplanetary coronal mass ejection into the Earth's magnetosheath is studied by investigating in situ magnetic field measurements upstream and downstream of the bow shock during an ICME sheath passage on 15 May 2005. We observe three distinct intervals in the immediate upstream region that included a southward magnetic field component and are traveling foreshocks. These traveling foreshocks were observed in the quasi-parallel bow shock that hosted backstreaming ions and magnetic fluctuations at ultralow frequencies. The intervals constituting traveling foreshocks in the upstream survive transmission to the Earth's magnetosheath, where their magnetic field, and particularly the southward component, was significantly amplified. Our results further suggest that the magnetic field fluctuations embedded in an ICME sheath may survive the transmission if their frequency is below ∼0.01 Hz. Although one of the identified intervals was coherent, extending across the ICME sheath and being long-lived, predicting ICME sheath magnetic fields that may transmit to the Earth's magnetosheath from the upstream at L1 observations has ambiguity. This can result from the strong spatial variability of the ICME sheath fields in the longitudinal direction, or alternatively from the ICME sheath fields developing substantially within the short time it takes the plasma to propagate from L1 to the bow shock. This study demonstrates the complex interplay ICME sheaths have with the Earth's magnetosphere when passing by the planet.

We report on a statistical analysis of the occurrence and properties of Alfven ion cyclotron (AIC) waves in sheath regions driven by interplanetary coronal mass ejections (ICMEs). We have developed an automated algorithm to identify AIC wave events from magnetic field data and apply it to investigate 91 ICME sheath regions recorded by the Wind spacecraft. Our analysis focuses on waves generated by the ion cyclotron instability. AIC waves are observed to be frequent structures in ICME-driven sheaths, and their occurrence is the highest in the vicinity of the shock. Together with previous studies, our results imply that the shock compression has a crucial role in generating wave activity in ICME sheaths. AIC waves tend to have their frequency below the ion cyclotron frequency, and, in general, occur in plasma that is stable with respect to the ion cyclotron instability and has lower ion beta(parallel to) than mirror modes. The results suggest that the ion beta anisotropy beta(perpendicular to)/beta(parallel to) > 1 appearing in ICME sheaths is regulated by both ion cyclotron and mirror instabilities.

We present a comprehensive statistical analysis of mirror mode waves and the properties of their plasma surroundings in sheath regions driven by interplanetary coronal mass ejection (ICME). We have constructed a semi-automated method to identify mirror modes from the magnetic field data. We analyze 91 ICME sheath regions from January 1997 to April 2015 using data from the Wind spacecraft. The results imply that similarly to planetary magnetosheaths, mirror modes are also common structures in ICME sheaths. However, they occur almost exclusively as dip-like structures and in mirror stable plasma. We observe mirror modes throughout the sheath, from the bow shock to the ICME leading edge, but their amplitudes are largest closest to the shock. We also find that the shock strength (measured by Alfven Mach number) is the most important parameter in controlling the occurrence of mirror modes. Our findings suggest that in ICME sheaths the dominant source of free energy for mirror mode generation is the shock compression. We also suggest that mirror modes that are found deeper in the sheath are remnants from earlier times of the sheath evolution, generated also in the vicinity of the shock.

The objective of the SQUID project is to develop and in flight verify a miniature version of a wire boom deployment mechanism to be used for electric field measurements in the ionosphere. In February 2011 a small ejectable payload, built by a team of students from The Royal Institute of Technology (KTH), was launched from Esrange on-board the REXUS-10 sounding rocket. The payload separated from the rocket, deployed and retracted the wire booms, landed with a parachute and was subsequently recovered. Here the design of the experiment and post fight analysis are presented.

Induced magnetospheres form around planetary bodies with atmospheres through the interaction of the solar wind with their ionosphere. Induced magnetospheres are highly dependent on the solar wind conditions and have only been studied with single spacecraft missions in the past. Without simultaneous measurements of solar wind variations and phenomena in the magnetosphere, establishing a link between both can only be done indirectly, using statistics over a large set of measurements. This gap in knowledge could be addressed by a multi-spacecraft plasma mission, optimized for studying global spatial and temporal variations in the magnetospheric system around Venus, which hosts the most prominent example of an induced magnetosphere in our solar system. The MVSE mission comprises four satellites, of which three are identical scientific spacecraft, carrying the same suite of instruments probing different regions of the induced magnetosphere and the solar wind simultaneously. The fourth spacecraft is the transfer vehicle which acts as a relay satellite for communications at Venus. In this way, changes in the solar wind conditions and extreme solar events can be observed, and their effects can be quantified as they propagate through the Venusian induced magnetosphere. Additionally, energy transfer in the Venusian induced magnetosphere can be investigated. The scientific payload includes instrumentation to measure the magnetic field, electric field, and ion–electron velocity distributions. This study presents the scientific motivation for the mission as well as requirements and the resulting mission design. Concretely, a mission timeline along with a complete spacecraft design, including mass, power, communication, propulsion and thermal budgets are given. This mission was initially conceived at the Alpbach Summer School 2022 and refined during a week-long study at ESA's Concurrent Design Facility in Redu, Belgium.

Induced magnetospheres form around planetary bodies with atmospheres through the interaction of the solar wind with their ionosphere. Induced magnetospheres are highly dependent on the solar wind conditions and have only been studied with single spacecraft missions in the past. Without simultaneous measurements of solar wind variations and phenomena in the magnetosphere, establishing a link between both can only be done indirectly, using statistics over a large set of measurements. This gap in knowledge could be addressed by a multi-spacecraft plasma mission, optimized for studying global spatial and temporal variations in the magnetospheric system around Venus, which hosts the most prominent example of an induced magnetosphere in our solar system. The MVSE mission comprises four satellites, of which three are identical scientific spacecraft, carrying the same suite of instruments probing different regions of the induced magnetosphere and the solar wind simultaneously. The fourth spacecraft is the transfer vehicle which acts as a relay satellite for communications at Venus. In this way, changes in the solar wind conditions and extreme solar events can be observed, and their effects can be quantified as they propagate through the Venusian induced magnetosphere. Additionally, energy transfer in the Venusian induced magnetosphere can be investigated. The scientific payload includes instrumentation to measure the magnetic field, electric field, and ion–electron velocity distributions. This study presents the scientific motivation for the mission as well as requirements and the resulting mission design. Concretely, a mission timeline along with a complete spacecraft design, including mass, power, communication, propulsion and thermal budgets are given. This mission was initially conceived at the Alpbach Summer School 2022 and refined during a week-long study at ESA&apos;s Concurrent Design Facility in Redu, Belgium.

Turbulence, intermittency, and self-organized structures in space plasmas can be investigated by using a multifractal formalism mostly based on the canonical structure function analysis with fixed constraints about stationarity, linearity, and scales. Here, the Empirical Mode Decomposition (EMD) method is firstly used to investigate timescale fluctuations of the solar wind magnetic field components; then, by exploiting the local properties of fluctuations, the structure function analysis is used to gain insights into the scaling properties of both inertial and kinetic/dissipative ranges. Results show that while the inertial range dynamics can be described in a multifractal framework, characterizing an unstable fixed point of the system, the kinetic/dissipative range dynamics is well described by using a monofractal approach, because it is a stable fixed point of the system, unless it has a higher degree of complexity and chaos.

We study ghost-free multimetric theories for (N + 1) tensor fields with a coupling to matter and maximal global symmetry group S-N x (Z(2))(N). Their mass spectra contain a massless mode, the graviton, and N massive spin-2 modes. One of the massive modes is distinct by being the heaviest, the remaining (N - 1) massive modes are simply identical copies of each other. All relevant physics can therefore be understood from the case N = 2. Focussing on this case, we compute the full perturbative action up to cubic order and derive several features that hold to all orders in perturbation theory. The lighter massive mode does not couple to matter and neither of the massive modes decay into massless gravitons. We propose the lighter massive particle as a candidate for dark matter and investigate its phenomenology in the parameter region where the matter coupling is dominated by the massless graviton. The relic density of massive spin-2 can originate from a freeze-in mechanism or from gravitational particle production, giving rise to two different dark matter scenarios. The allowed parameter regions are very different from those in scenarios with only one massive spin-2 field and more accessible to experiments.

Starting from the leading soft term of the 5-point amplitude, involving a graviton and two Kerr black holes, that factorises into the product of the elastic amplitude without the graviton and the leading soft factor, we compute the infrared divergent contribution to the imaginary part of the two-loop eikonal. Then, using analyticity and crossing symmetry, we determine the radiative contribution to the real part of the two-loop eikonal and from it the radiative part of the deflection angle for spins aligned to the orbital angular momentum, the loss of angular momentum and the zero frequency limit of the energy spectrum for any spin and for any spin orientation. For spin one we find perfect agreement with recent results obtained with the supersymmetric worldline formalism.