Source author record

M. W. Kunz

M. W. Kunz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

10works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

10 published item(s)

preprint2022arXiv

The structure of cluster merger shocks: turbulent width and the electron heating timescale

We present a new 2 Ms Chandra observation of the cluster merger Abell 2146, which hosts two huge M~2 shock fronts each ~500 kpc across. For the first time, we resolve and measure the width of cluster merger shocks. The best-fit width for the bow shock is 17+/-1 kpc and for the upstream shock is 10.7+/-0.3 kpc. A narrow collisionless shock will appear broader in projection if its smooth shape is warped by local gas motions. We show that both shock widths are consistent with collisionless shocks blurred by local gas motions of 290+/-30 km/s. The upstream shock forms later on in the merger than the bow shock and is therefore expected to be significantly narrower. From the electron temperature profile behind the bow shock, we measure the timescale for the electrons and ions to come back into thermal equilibrium. We rule out rapid thermal equilibration of the electrons with the shock-heated ions at the 6 sigma level. The observed temperature profile instead favours collisional equilibration. For these cluster merger shocks, which have low sonic Mach numbers and propagate through a high $β$ plasma, we find no evidence for electron heating over that produced by adiabatic compression. Our findings are expected to be valid for collisionless shocks with similar parameters in other environments and support the existing picture from the solar wind and supernova remnants. The upstream shock is consistent with this result but has a more complex structure, including a ~2 keV increase in temperature ~50 kpc ahead of the shock.

preprint2020arXiv

Electron-only magnetic reconnection in plasma turbulence

Hybrid-Vlasov-Maxwell simulations of magnetized plasma turbulence including non-linear electron-inertia effects in a generalized Ohm's law are presented. When fluctuation energy is injected on scales sufficiently close to ion-kinetic scales, the ions efficiently become de-magnetized and electron-scale current sheets largely dominate the distribution of the emerging current structures, in contrast to the usual picture, where a full hierarchy of structure sizes is generally observed. These current sheets are shown to be the sites of electron-only reconnection (e-rec), in which the usual electron exhausts are unaccompanied by ion outflows and which are in qualitative agreement with those recently observed by MMS in the Earth's turbulent magnetosheath, downstream of the bow shock. Some features of the e-rec phenomenology are shown to be consistent with an electron magnetohydrodynamic description. Simulations suggest that this regime of collisionless reconnection may be found in turbulent systems where plasma processes, such as micro-instabilities and/or shocks, overpower the more customary turbulent cascade by directly injecting energy close to the ion-kinetic scales.

preprint2020arXiv

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle--particle collisions allows departures from local thermodynamic equilibrium. These departures exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present a numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart. If instead the parallel viscous stress is left unabated -- a situation relevant to recent kinetic simulations of the fluctuation dynamo and to the early stages of the dynamo in a magnetized ICM -- the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics of the saturated state of the large-Prandtl-number (${\rm Pm}\gtrsim{1}$) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches magnetic-energy spectra. A prediction of this model, confirmed by our simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. In all cases that do exhibit a dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

preprint2016arXiv

Multi-Species Measurements of the Firehose and Mirror Instability Thresholds in the Solar Wind

The firehose and mirror instabilities are thought to arise in a variety of space and astrophysical plasmas, constraining the pressure anisotropies and drifts between particle species. The plasma stability depends on all species simultaneously, meaning that a combined analysis is required. Here, we present the first such analysis in the solar wind, using the long-wavelength stability parameters to combine the anisotropies and drifts of all major species (core and beam protons, alphas, and electrons). At the threshold, the firehose parameter was found to be dominated by protons (67%), but also to have significant contributions from electrons (18%) and alphas (15%). Drifts were also found to be important, contributing 57% in the presence of a proton beam. A similar situation was found for the mirror, with contributions of 61%, 28%, and 11% for protons, electrons, and alphas, respectively. The parallel electric field contribution, however, was found to be small at 9%. Overall, the long-wavelength thresholds constrain the data well (<1% unstable), and the implications of this are discussed.

preprint2016arXiv

Suppression of thermal conduction in a mirror-unstable plasma

The ICM plasma is subject to firehose and mirror instabilities at scales of order the ion Larmor radius. The mirror instability generates fluctuations of magnetic-field strength $δB / B \sim 1$. These fluctuations act as magnetic traps for the heat-conducting electrons, suppressing their transport. We calculate the effective parallel thermal conductivity in the ICM in the presence of the mirror fluctuations for different stages of the evolution of the instability. The mirror fluctuations are limited in amplitude by the maximum and minimum values of the field strength, with no large deviations from the mean value. This key property leads to a finite suppression of thermal conduction at large scales. We find suppression down to $\approx 0.2$ of the Spitzer value for the secular phase of the perturbations' growth, and $\approx 0.3$ for their saturated phase. The effect operates in addition to other suppression mechanisms and independently of them. Globally, fluctuations $δB / B \sim 1$ can be present on much larger scales, of the order of the scale of turbulent motions. However, we do not expect large suppression of thermal conduction by these, because their scale is considerably larger than the collisional mean free path of the ICM electrons. The obtained suppression of thermal conduction by a factor of $\sim 5$ appears to be characteristic and potentially universal for a weakly collisional mirror-unstable plasma.

preprint2015arXiv

Inertial-Range Kinetic Turbulence in Pressure-Anisotropic Astrophysical Plasmas

A theoretical framework for low-frequency electromagnetic (drift-)kinetic turbulence in a collisionless, multi-species plasma is presented. The result generalises reduced magnetohydrodynamics (RMHD) and kinetic RMHD (Schekochihin et al. 2009) for pressure-anisotropic plasmas, allowing for species drifts---a situation routinely encountered in the solar wind and presumably ubiquitous in hot dilute astrophysical plasmas (e.g. intracluster medium). Two main objectives are achieved. First, in a non-Maxwellian plasma, the relationships between fluctuating fields (e.g., the Alfven ratio) are order-unity modified compared to the more commonly considered Maxwellian case, and so a quantitative theory is developed to support quantitative measurements now possible in the solar wind. The main physical feature of low-frequency plasma turbulence survives the generalisation to non-Maxwellian distributions: Alfvenic and compressive fluctuations are energetically decoupled, with the latter passively advected by the former; the Alfvenic cascade is fluid, satisfying RMHD equations (with the Alfven speed modified by pressure anisotropy and species drifts), whereas the compressive cascade is kinetic and subject to collisionless damping. Secondly, the organising principle of this turbulence is elucidated in the form of a generalised kinetic free-energy invariant. It is shown that non-Maxwellian features in the distribution function reduce the rate of phase mixing and the efficacy of magnetic stresses; these changes influence the partitioning of free energy amongst the various cascade channels. As the firehose or mirror instability thresholds are approached, the dynamics of the plasma are modified so as to reduce the energetic cost of bending magnetic-field lines or of compressing/rarefying them. Finally, it is shown that this theory can be derived as a long-wavelength limit of non-Maxwellian slab gyrokinetics.

preprint2015arXiv

Pressure-anisotropy-driven microturbulence and magnetic-field evolution in shearing, collisionless plasma

The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulting microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flow), so a new macroscale configuration is superimposed on the microscale state left behind by the previous one. There is a threshold value of plasma beta: when $β\llΩ/S$ (ion cyclotron frequency/shear rate), the emergence of firehose/mirror fluctuations driven unstable by shear and their disappearance when the shear is removed/reversed are quasi-instantaneous compared to the shear time, viz., the decay time of these fluctuations is $\simβ/Ω\ll 1/S$ (this result follows from the free decay of the fluctuations being constrained by the same marginal-stability thresholds as their growth). In contrast, when $β\gtrsimΩ/S$ ("ultra-high" beta), the old microscale state can only be removed on the shear timescale. In this regime, driven firehose fluctuations grow secularly to order-unity amplitudes, compensating for the decay of the mean field and so pinning pressure anisotropy at marginal stability with no appreciable scattering of particles---which is unlike what happens at moderate $β$. When the shear reverses, the shearing away of this firehose turbulence compensates for the increase in the mean field and thus prevents growth of the pressure anisotropy, stopping the system from going mirror-unstable. Therefore, at ultra-high beta, the system stays close to the firehose threshold, the mirror instability is largely suppressed, while the mean magnetic energy barely changes at all. Implications for plasma dynamo and thus the origin of cosmic magnetism are discussed.

preprint2015arXiv

The Effect of Anisotropic Viscosity on Cold Fronts in Galaxy Clusters

Cold fronts -- contact discontinuities in the intracluster medium (ICM) of galaxy clusters -- should be disrupted by Kelvin-Helmholtz (K-H) instabilities due to the associated shear velocity. However, many observed cold fronts appear stable. This opens the possibility to place constraints on microphysical mechanisms that stabilize them, such as the ICM viscosity and/or magnetic fields. We performed exploratory high-resolution simulations of cold fronts arising from subsonic gas sloshing in cluster cores using the grid-based Athena MHD code, comparing the effects of isotropic Spitzer and anisotropic Braginskii viscosity (expected in a magnetized plasma). Magnetized simulations with full Braginskii viscosity or isotropic Spitzer viscosity reduced by a factor f ~ 0.1 are both in qualitative agreement with observations in terms of suppressing K-H instabilities. The RMS velocity of turbulence within the sloshing region is only modestly reduced by Braginskii viscosity. We also performed unmagnetized simulations with and without viscosity and find that magnetic fields have a substantial effect on the appearance of the cold fronts, even if the initial field is weak and the viscosity is the same. This suggests that determining the dominant suppression mechanism of a given cold front from X-ray observations (e.g. viscosity or magnetic fields) by comparison with simulations is not straightforward. Finally, we performed simulations including anisotropic thermal conduction, and find that including Braginskii viscosity in these simulations does not significant affect the evolution of cold fronts; they are rapidly smeared out by thermal conduction, as in the inviscid case.

preprint2010arXiv

A thermally stable heating mechanism for the intracluster medium: turbulence, magnetic fields and plasma instabilities

We consider the problem of self-regulated heating and cooling in galaxy clusters and the implications for cluster magnetic fields and turbulence. Viscous heating of a weakly collisional magnetised plasma is regulated by the pressure anisotropy with respect to the local direction of the magnetic field. The intracluster medium is a high-beta plasma, where pressure anisotropies caused by the turbulent stresses and the consequent local changes in the magnetic field will trigger very fast microscale instabilities. We argue that the net effect of these instabilities will be to pin the pressure anisotropies at a marginal level, controlled by the plasma beta parameter. This gives rise to local heating rates that turn out to be comparable to the radiative cooling rates. Furthermore, we show that a balance between this heating and Bremsstrahlung cooling is thermally stable, unlike the often conjectured balance between cooling and thermal conduction. Given a sufficient (and probably self-regulating) supply of turbulent power, this provides a physical mechanism for mitigating cooling flows and preventing cluster core collapse. For observed density and temperature profiles, the assumed balance of viscous heating and radiative cooling allows us to predict magnetic-field strengths, turbulent velocities and turbulence scales as functions of distance from the centre. Specific predictions and comparisons with observations are given for several different clusters. Our predictions can be further tested by future observations of cluster magnetic fields and turbulent velocities.

preprint2010arXiv

The nonisothermal stage of magnetic star formation. II. Results

In a previous paper we formulated the problem of the formation and evolution of fragments (or cores) in magnetically-supported, self-gravitating molecular clouds in axisymmetric geometry, accounting for the effects of ambipolar diffusion and Ohmic dissipation, grain chemistry and dynamics, and radiative transfer. Here we present results of star formation simulations that accurately track the evolution of a protostellar fragment over eleven orders of magnitude in density (from 300 cm^-3 to \approx 10^14 cm^-3), i.e., from the early ambipolar-diffusion--initiated fragmentation phase, through the magnetically supercritical, dynamical-contraction phase and the subsequent magnetic decoupling stage, to the formation of a protostellar core in near hydrostatic equilibrium. As found by Fiedler & Mouschovias (1993), gravitationally-driven ambipolar diffusion leads to the formation and subsequent dynamic contraction of a magnetically supercritical core. Moreover, we find that ambipolar diffusion, not Ohmic dissipation, is responsible for decoupling all the species except the electrons from the magnetic field, by a density \approx 3 x 10^12 cm^-3. Magnetic decoupling precedes the formation of a central stellar object and ultimately gives rise to a concentration of magnetic flux (a `magnetic wall') outside the hydrostatic core --- as also found by Tassis & Mouschovias (2005a,b) through a different approach. At approximately the same density at which Ohmic dissipation becomes more important than ambipolar diffusion (\gtrsim 7 x 10^12 cm^-3), the grains carry most of the electric charge as well as the electric current. The prestellar core remains disclike down to radii ~ 10 AU, inside which thermal pressure becomes important. The magnetic flux problem of star formation is resolved for at least strongly magnetic newborn stars by this stage of the evolution, i.e., by a central density \approx 10^14 cm^-3. The hydrostatic core has radius \approx 2 AU, density \approx 10^14 cm^-3, temperature \approx 300 K, magnetic field strength \approx 0.2 G, magnetic flux \approx 5 x 10^18 Wb, luminosity ~ 10^-3 L_\odot, and mass ~ 10^-2 M_\odot.