
Mean field electrodynamics (MFE) facilitates practical modeling of secular,
large scale properties of astrophysical or laboratory systems with
fluctuations.Practitioners commonly assume wide scale separation between mean
and fluctuating quantities, to justify equality of ensemble and spatial or
temporal averages.Often however, real systems do not exhibit such scale
separation. This raises two questions: (I) what are the appropriate generalized
equations of MFE in the presence of mesoscale fluctuations? (II) how precise
are theoretical predictions from MFE? We address both by first deriving the
equations of MFE for different types of averaging, along with mesoscale
correction terms that depend on the ratio of averaging scale to variation scale
of the mean. We then show that even if these terms are small, predictions of
MFE can still have a significant precision error. This error has an intrinsic
contribution from the dynamo input parameters and a filtering contribution from
differences in the way observations and theory are projected through the
measurement kernel.Minimizing the sum of these contributions can produce an
optimal scale of averaging that makes the theory maximally precise.The
precision error is important to quantify when comparing to observations because
it quantifies the resolution of predictive power. We exemplify these principles
for galactic dynamos, comment on broader implications, and identify
possibilities for further work.

We present a new magnetohydrodynamic (MHD) simulation code with the aim of
providing accurate numerical solutions to astrophysical phenomena where
discontinuities, shock waves, and turbulence are inherently important. The code
implements the HLLD approximate Riemann solver, the
fifthordermonotonicitypreserving interpolation (MP5) scheme, and the
hyperbolic divergence cleaning method for a magnetic field. This choice of
schemes significantly improved numerical accuracy and stability, and saved
computational costs in multidimensional problems. Numerical tests of one and
twodimensional problems showed the advantages of using the highorder scheme
by comparing with results from a standard secondorder TVD MUSCL scheme. The
present code enabled us to explore longterm evolution of a threedimensional
accretion disk around a black hole, in which compressible MHD turbulence caused
continuous mass accretion via nonlinear growth of the magnetorotational
instability (MRI). Numerical tests with various computational cell sizes
exhibited a convergent picture of the early nonlinear growth of the MRI in a
global model, and indicated that the MP5 scheme has more than twice the
resolution of the MUSCL scheme in practical applications.

Along the development of freeelectron laser operating at the wavelength of
Xray, the importance of investigation on the radiation has increased. A
theoretical simulation is an essential tool for studying existing and proposed
experiments. The available simulations preserve net energy conservation over
the timestep, upholding causality between the electrons' power and the radiated
energy flux. However, according to WheelerFeynman timesymmetric theory, the
timestep is too short to ensure this. Therefore, the time evolution of
electrons and radiation field predicted by the simulations should contradict
the correct physics, regardless of whether the stimulated emission dominates
over the spontaneous emission or not.

We study the slowing down of particle beams passing through the dusty plasma
with powerlaw kappadistributions. Three plasma components, electrons, ions
and dust particles, can have different kappaparameter. We derive the
deceleration factor (the velocity moment equation) and the slowing down time of
a test particle, and numerically study the slowing down properties of an
electron beam, a proton beam and a dust particle beam, respectively, in the
kappadistributed dusty plasma. We show that the slowing down properties of
particle beams depend strongly on the kappaparameters of the plasma
components, and the dust component plays a dominant role in the slowing down.
And the slowing down also depends on mass and charge of the dust particles in
the dusty plasma. More detailed results are shown in 17 numerical graphs.

Timedependent nonlinear media, such as rapidly generated plasmas produced
via laser ionization of gases, can increase the energy of individual laser
photons and generate tunable highorder harmonic pulses. This phenomenon, known
as photon acceleration, has traditionally required extremeintensity laser
pulses and macroscopic propagation lengths. Here, we report on a novel
nonlinear material$$an ultrathin semiconductor metasurface$$that exhibits
efficient photon acceleration at low intensities. We observe a signature
nonlinear manifestation of photon acceleration: thirdharmonic generation of
nearinfrared photons with tunable frequencies reaching up to
$\approx3.1\omega$. A simple timedependent coupledmode theory, found to be in
good agreement with experimental results, is utilized to predict a new path
towards nonlinear radiation sources that combine resonant upconversion with
broadband operation.

Based on the particleincell (PIC) plasma simulation method, the
speedlimited PIC (SLPIC) method delivers faster kinetic plasma simulation in
cases where the particle distributions evolve slowly compared with the maximum
stable PIC timestep. SLPIC thus offers more feasible, fully kinetic simulation
in regimes that historically have required fluid approaches, such as
magnetohydrodynamic (MHD), twofluid, or Boltzmann electron treatments. In
particular, SLPIC allows an explicit time advance with steps much larger than
the inverse plasma frequency, avoiding the instability explicit PIC faces with
large timesteps. SLPIC avoids this instability by slowing down fast particles
(e.g., electrons) in a way that is rigorously underpinned by an approximate
Vlasov equation; unlike MHD, twofluid, and Boltzmann electron approaches,
SLPIC does not fundamentally neglect any firstprinciples plasma physics,
although the choices of grid cell size, timestep, and number of macroparticles
per cell naturally limit the physical phenomena that can be accurately
represented. SLPIC can be implemented with minor modifications of a standard
PIC code and does not require an implicit time advance. It enables large
timesteps in firstprinciples kinetic plasma simulation of appropriately slow
phenomena, and it can handle many of the same complications as PIC, such as
boundary conditions and collisions. In an argon plasma sheath test problem, a
SLPIC simulation achieved a speedup of a factor of 160 over the corresponding
PIC simulation, without loss of accuracy.

Although turbulence has been conjectured to be important for magnetic
reconnection, still very little is known about its role in collisionless
plasmas. Previous attempts to quantify the effect of turbulence on reconnection
usually prescribed Alfv\'enic or other lowfrequency fluctuations or
investigated collisionless kinetic effects in just twodimensional
configurations and antiparallel magnetic fields. In view of this, we analyzed
the kinetic turbulence selfgenerated by threedimensional guidefield
reconnection through forcefree current sheets in frequency and wavenumber
spaces, utilizing 3D particleincell code numerical simulations. Our
investigations reveal reconnection rates and kinetic turbulence with features
similar to those obtained by current insitu spacecraft observations of MMS as
well as in the laboratory reconnection experiments MRX, VTF and
\textsc{Vineta}II. In particular we found that the kinetic turbulence
developing in the course of 3D guidefield reconnection exhibits a broadband
powerlaw spectrum extending beyond the lowerhybrid frequency and up to the
electron frequencies. In the frequency space the spectral index of the
turbulence appeared to be close to 2.8 at the reconnection Xline. In the
wavenumber space it also becomes 2.8 as soon as the normalized reconnection
rate reaches 0.1. The broadband kinetic turbulence is mainly due to
currentstreaming and electronflowshear instabilities excited in the
sufficiently thin current sheets of kinetic reconnection. The growth of the
kinetic turbulence corresponds to high reconnection rates which exceed those of
fast laminar, nonturbulent reconnection.

The Lagrangian mechanical consideration of the dynamics of ideal
incompressible hydrodynamic, magnetohydrodynamic, and Hall magnetohydrodynamic
media, which are formulated as dynamical systems in appropriate Lie groups
equipped with Riemannian metrics, leads to the notion of generalized
vorticities, as well as generalized coordinates, velocities, and momenta. The
action of each system is conserved against the integral path variation in the
direction of the generalized vorticity, and this invariance is associated with
the particle relabeling symmetry. The generalized vorticities are formulated by
the operation of integrodifferential operators on the generalized velocities.
The eigenfunctions of the operators provide sets of orthogonal functions, and
we obtain common mathematical expressions concerning these dynamical systems
using the orthogonoal functions. In particular, we find that the product of the
Riemannian metric, $g_{im}$, and the structure constants of the Lie group,
$C^{m}_{jk}$, is given by the product of the eigenvalue of the operator,
$\Lambda(i)$, and a certain totally antisymmetric tensor, $T_{ijk}$:
$g_{i\alpha}C^{\alpha}_{jk}=\Lambda(i)T_{ijk}.$ Its physical implications,
including the weak interaction conjecture of MHD turbulence, are also
discussed.

Modifications of the equations of ideal fluid dynamics with advected
quantities are introduced that allow selective decay of either the energy $h$
or the Casimir quantities $C$ in the LiePoisson formulation. The dissipated
quantity (energy or Casimir, respectively) is shown to decrease in time until
the modified system reaches an equilibrium state consistent with ideal
energyCasimir equilibria, namely $\delta(h+C)=0$. The result holds for
LiePoisson equations in general, independently of the Lie algebra and the
choice of Casimir. This selective decay process is illustrated with a number of
examples in 2D and 3D magnetohydrodynamics (MHD).

The gyrotron is the high power millimeter wave source in the electron
cyclotron resonance heating system for the tokamak. The cathode power supply is
one of the most important ancillary devices for gyrotron. Some interesting
transient phenomena about the cathode voltage and the cathode current was found
in the gyrotron operation in the electron cyclotron resonance heating system
for the Experimental Advanced Superconducting Tokamak. The cathode voltage
drops to about 10% of the original value by about 90 {\mu}s in the overcurrent
case, which is much longer than the 25 {\mu}s in the normal case. In order to
explain these phenomena, an equivalent circuit model of the magnetron injection
gun was proposed. The equivalent circuit is composed of parallel resistors and
capacitors, and it can explain the test results very well. Using the equivalent
circuit model to analyze gyrotrons may become an effective means of developing
and operating an electron cyclotron heating system.

Using spatiotemporal spectra we show direct evidence of excitation of
magnetosonic and Alfv\'en waves in threedimensional compressible
magnetohydrodynamic turbulence at small Mach numbers. For the plasma pressure
dominated regime, or high $\beta$ regime (with $\beta$ the ratio between fluid
and magnetic pressure), and for the magnetic pressure dominated regime, or low
$\beta$ regime, we study magnetic field fluctuations parallel and perpendicular
to a guide magnetic field $\textbf{B}_0$. In the low $\beta$ case we find
excitation of compressible and incompressible fluctuations, with a transfer of
energy towards Alfv\'enic modes and to a lesser extent towards magnetosonic
modes. In particular, we find signatures of the presence of fast magnetosonic
waves in a scenario compatible with that of weak turbulence. In the high
$\beta$ case, fast and slow magnetosonic waves are present, with no clear trace
of Alfv\'en waves, and a significant part of the energy is carried by
twodimensional turbulent eddies.

An example of the nonequilibrium phase transition is the formation of lanes
when one kind of particles is driven against the other. According to
experimental observation, lane formation in binary complex plasmas occurs when
the smaller particles are driven through the stationary dust cloud of the
larger particles. We calculate the driving force acting on a probe particle
that finds itself in a quiescent cloud of particles in complex plasma of the
lowpressure radio frequency discharge under microgravity conditions. It is
shown that the nonzero driving force is a result of the dependence of the ion
mean free path on the particle number density. If this effect is properly
included in the model of similar complex plasmas then one arrives at the
driving force that changes its sign at the point where the probe and the dust
particles have equal radii. If the probe is smaller than the dust particle then
the driving force is directed toward the discharge center and vice versa, in
accordance with experiment. Obtained results can serve as the ansatz for future
investigation of the lane formation in complex plasmas.

The effects of the stellar rotation and the consequent azimuthal stellar wind
flow in Parker's [9] hydrodynamic stellar wind model are discussed. Of special
interest is the emergence of a whole new hydrodynamic physics via a new
critical point in the stellar wind flow, which supersedes the critical point in
Parker's [9] hydrodynamic model. The effect of the stellar rotation is shown to
cause the new critical point to occur lower in the corona, so the stellar wind
experiences a stronger afterburner (as in an aircraft jet engine) action in the
corona. For strong rotators, the new critical point is shown to occur at a
fixed location for a given star, determined only by the basic stellar
parameters like the mass M and the angular velocity, the variations in the
stellar wind environment notwithstanding. The stellar rotation leads to
stronger density falloff and enhanced acceleration of the stellar wind at
large distances from the star  this effect materializes even close to the
star, for strong rotators. The stellar rotation causes the physical throat
section of the effective de Laval nozzle associated with the stellar wind flow
to become narrower, indicative of an enhanced flow acceleration. Thus, the
stellar rotation leads to tenuous and faster stellar wind flows without change
in the mass flux, and hence provides an efficient physical mechanism for
protostars and strong rotators to lose their angular momentum quickly.

Ab initio simulations of the propagation in a plasma of a soon to be
available relativistic electronpositron beam or fireball beam provide an
effective mean for the study of microphysics relevant to astrophysical
scenarios. We show that the current filamentation instability associated with
some of these scenarios reaches saturation after only 10 cm of propagation in a
typical laboratory plasma with a density 10^17/cc. The different regimes of the
instability, from the purely transverse to the mixed mode filamentation, can be
accessed by varying the background plasma density. The instability generates
large local plasma gradients, intense transverse magnetic fields, and enhanced
emission of radiation. We suggest that these effects may be observed
experimentally for the first time.

We report a newly found twostage mechanism of electron acceleration near
Xlines of 3D collisionless guidefield magnetic reconnection in the
nonrelativistic regime typical, e.g., for stellar coronae. We found that after
electrons are first preaccelerated during the linear growth of reconnection,
they become additionally accelerated in the course of the nonlinear stage of 3D
guidefield magnetic reconnection. This additional acceleration is due to the
filamentation of electric and magnetic fields caused by streaming
instabilities. In addition to enhanced parallel electric fields, the
filamentation leads to additional curvaturedriven electron acceleration in the
guidefield direction. As a result, part of the the accelerated electron
spectra becomes a power law with a spectral index of $\sim1.6$ near the
Xline. This second stage of acceleration due to nonlinear reconnection is
relevant for the production of energetic electrons in, e.g., thin current
sheets of stellar coronae.

The dispersion relation of plasma and plasmadielectric photonic multilayer
structures is approached in terms of a onedimensional ParticleinCell
simulation. For several plasmadielectric configurations, the system response
is obtained using a pulsed excitation and a subsequent twodimensional
frequency analysis. It is first shown that the dispersion relation of a single,
homogeneous plasma slab is well described by the coldplasma model even at low
pressures of 1 Pa. The study is extended to the simulation of plasma photonic
crystals with a variety of configurations, based on the work of Hojo and Mase
[J. Plasma Fusion Res. 80, 89 (2004)]. Considering a onedimensional plasma
photonic crystal made from alternating layers of dielectric and homogeneous
plasma slabs, it is shown that the assumption of a coldplasma description is
well justified also in this case. Moreover, in this work the results are
reformatted and analyzed in a band diagram representation, in particular based
on the lattice constant $a$. Based on these considerations a scaling invariant
representation is presented, utilizing a generalized set of parameters. The
study is completed with an exemplary comparison of three plasmadielectric
photonic crystal configurations and their corresponding band diagrams.

Using a helium ion hitting various metal surfaces as a model system, we
describe a general quantumkinetic approach for calculating ioninduced
secondary electron emission spectra at impact energies where the emission is
driven by the internal potential energy of the ion. It is based on an effective
model of the AndersonNewnstype for the subset of electronic states of the
ionsurface system most strongly affected by the collision. Central to our
approach is a pseudoparticle representation for the electronic configurations
of the projectile which enables us, by combining it with two additional
auxiliary bosons, to describe in a single Hamiltonian emission channels
involving electronic configurations with different internal potential energies.
It is thus possible to treat Auger neutralization of the ion on an equal
footing with Auger deexcitation of temporarily formed radicals and/or negative
ions. From the Dyson equations for the projectile propagators and an
approximate evaluation of the selfenergies, rate equations are obtained for
the probabilities with which the projectile configurations occur and an
electron is emitted in the course of the collision. Encouraging numerical
results, especially for the heliumtungsten system, indicate the potential of
the approach.

We present an approach for carrying out nonadiabatic molecular dynamics
simulations of systems in which nonadiabatic transitions arise from the
coupling between the classical atomic motions and a quasicontinuum of
electronic quantum states. Such conditions occur in many research areas,
including chemistry at metal surfaces, radiation damage of materials, and warm
dense matter physics. The classical atomic motions are governed by stochastic
Langevinlike equations, while the quantum electron dynamics is described by a
master equation for the populations of the electronic states. These working
equations are obtained from a firstprinciple derivation. Remarkably, unlike
the widely used Ehrenfest and surfacehopping methods, the approach naturally
satisfies the principle of detailed balance at equilibrium and, therefore, can
describe the evolution to thermal equilibrium from an arbitrary initial state.
In addition, unlike other schemes, there is no need to explicitly propagate
wave functions in time.

The MaxwellProca electrodynamics corresponding to a finite photon mass
causes a substantial change of the Maxwell stress tensor and, under certain
circumstances, may cause the electromagnetic stresses to act effectively as
"negative pressure." The paper describes a model where this negative pressure
imitates gravitational pull and may produce forces comparable to gravity and
even become dominant. The effect is associated with the random magnetic fields
in the galactic disk with a scale exceeding the photon Compton wavelength. The
presence of a weaker regular field does not affect the forces under
consideration. The stresses act predominantly on the interstellar gas and cause
an additional force pulling the gas towards the center and towards the galactic
plane. The stars do not experience any significant direct force but get
involved in this process via a "recycling loop" where rapidly evolving massive
stars are formed from the gas undergoing galactic rotation and then lose their
masses back to the gas within a time shorter than roughly 1/6 of the rotation
period. This makes their dynamics inseparable from that of the rotating gas.
The lighter, slowly evolving stars, as soon as they are formed, lose connection
to the gas and are confined within the galaxy only gravitationally. Numerical
examples based on the parameters of our galaxy reveal both opportunities and
challenges of this model and motivate further analysis. The critical issue is
the plausibility of formation of the irregular magnetic field that would be
force free. Another challenge is developing a predictive model of the evolution
of the gaseous and stellar population of the galaxy under the aforementioned
scenario. It may be interesting to also explore possible broader cosmological
implications of the negativepressure model.

In the present work, we have studied the collective behaviour of a chiral
plasma with first and second order dissipative corrections allowed by conformal
symmetry. We have derived dispersion relations for ideal, first order and
second order MHD for the chiral plasma. We have also used our results to
explain the observed pulsar kicks.

Using the freefree continuum selfemission spectrum at photon energies above
15 keV is one of the most promising concepts for assessing the electron
temperature in ICF experiments. However, these photons are due to suprathermal
electrons whose meanfreepath is much larger than thermal, making their
distribution deviate from Maxwellian in a finitesize hotspot. The first study
of the freefree Xray emission from an ICF implosion is conducted with the
kinetic modifications to the electron distribution accounted for. These
modifications are found to result in qualitatively new features in the hard
Xray spectral continuum. Inference of the electron temperature as if the
emitting electrons are Maxwellian is shown to give a lower value than the
actual one.

Asymmetric Fano line profiles are frequently encountered, e.g., in the
photoionization spectra of atoms and ions. For the fitting of spectral line
profiles to experimental spectra the line profiles have to be convolved with
the experimental window function. The latter is often taken to be a Gaussian.
It is shown that the convolution can be represented by a rather simple analytic
expression involving the Faddeeva function for the evaluation of which
efficient and accurate numerical algorithms are available.

In contrast to hydrodynamic vortices, vortices in plasma contain an electric
current circulating around the center of the vortex, which generates a magnetic
field localized inside. Using computer simulations, we demonstrate that the
magnetic field associated with the vortex gives rise to a mechanism of
dissipation of the vortex pair in a collisionless plasma, leading to fast
annihilation of the magnetic field with its energy transforming into the energy
of fast electrons, secondary vortices, and plasma waves. Two major contributors
to the energy damping of double vortex system, namely, magnetic field
annihilation and secondary vortex formation, are regulated by the size of the
vortex with respect to the electron skin depth, which scales with the electron
gammafactor, $\gamma_e$, as $R/d_e \propto \gamma_e^{1/2}$. Magnetic field
annihilation appears to be dominant in mildly relativistic vortices, while for
the ultrarelativistic case, secondary vortex formation is the main channel for
damping of the initial double vortex system.

The description of the dynamics of an electron in an external electromagnetic
field of arbitrary intensity is one of the most fundamental outstanding
problems in electrodynamics. Remarkably, to date there is no unanimously
accepted theoretical solution for ultrahigh intensities and little or no
experimental data. The basic challenge is the inclusion of the selfinteraction
of the electron with the field emitted by the electron itself  the socalled
radiation reaction force. We report here on the experimental evidence of strong
radiation reaction, in an alloptical experiment, during the propagation of
highly relativistic electrons (maximum energy exceeding 2 GeV) through the
field of an ultraintense laser (peak intensity of $4\times10^{20}$ W/cm$^2$).
In their own rest frame, the highest energy electrons experience an electric
field as high as one quarter of the critical field of quantum electrodynamics
and are seen to lose up to 30% of their kinetic energy during the propagation
through the laser field. The experimental data show signatures of quantum
effects in the electron dynamics in the external laser field, potentially
showing departures from the constant cross field approximation.

Pulses of extreme ultraviolet (XUV) light, with wavelengths between 10 and
100$\,$nm, can be used to image and excite ultrafast phenomena such as the
motion of atomic electrons. Here we show that the illumination of plasma with
nearcritical electron density may be used as a source of relativistically
intense XUV radiation, providing the means for novel XUVpumpXUVprobe
experiments in the nonlinear regime. We describe how the optimal regime may be
reached by tailoring the lasertarget interaction parameters and by the
presence of preplasma. Our results indicate that currently available laser
facilities are capable of producing XUV pulses with duration $\sim
10~\text{fs}$, brilliance in excess of $10^{23}$ photons/s/mm$^2$/mrad$^2$
(0.1% bandwidth) and intensity $I\lambda^2 \gtrsim
10^{19}~\text{W}\text{cm}^{2}\mu\text{m}^2$.