Researcher profile

G. P. Zank

G. P. Zank contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

18 published item(s)

preprint2022arXiv

Action Principles and Conservation Laws for Chew-Goldberger-Low Anisotropic Plasmas

The ideal CGL plasma equations, including the double adiabatic conservation laws for the parallel ($p_\parallel$) and perpendicular pressure ($p_\perp$), are investigated using a Lagrangian variational principle. An Euler-Poincaré variational principle is developed and the non-canonical Poisson bracket is obtained, in which the non-canonical variables consist of the mass flux ${\bf M}$, the density $ρ$, three entropy variables, $σ=ρS$, $σ_\parallel=ρS_\parallel$, $σ_\perp=ρS_\perp$ ($S_\parallel$ and $S_\perp$ are the two scalar entropy invariants), and the magnetic induction ${\bf B}$. Conservation laws of the CGL plasma equations are derived via Noether's theorem. The Galilean group leads to conservation of energy, momentum, center of mass, and angular momentum. Cross helicity conservation arises from a fluid relabeling symmetry, and is local or nonlocal depending on whether the entropy gradients of $S_\parallel$, $S_\perp$ and $S$ are perpendicular to ${\bf B}$ or otherwise. The point Lie symmetries of the CGL system are shown to comprise the Galilean transformations and scalings.

preprint2022arXiv

First comparison of composite 0.52-55 keV ENA spectra observed by IBEX and Cassini/INCA with simulated ENAs inferred by proton hybrid simulations downstream of the termination shock

We present a first comparison of Energetic Neutral Atom (ENA) heliosheath measurements, remotely sensed by the Interstellar Boundary Explorer (IBEX) mission and the Ion and Neutral Camera (INCA) on the Cassini mission, with modeled ENA inferred from interstellar pickup protons that have been accelerated at the termination shock, using hybrid simulations. The observed ENA intensities are an average value over the time period from 2009 to the end of 2012, along the Voyager 2 trajectory. The hybrid simulations parameters for the solar wind, interstellar pickup ions (PUIs), and magnetic field upstream of the termination shock, where Voyager 2 crossed, are based on observations. We report an energy dependent discrepancy between observed and simulated ENA fluxes, with the observed ENA fluxes, being consistently higher than the simulated ones, and discuss possible causes of this discrepancy.

preprint2022arXiv

Generalized Fluid Models of the Braginskii Type

Several generalizations of the well-known fluid model of Braginskii (Rev. of Plasma Phys., 1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multi-species plasmas with arbitrary masses and temperatures, where all the fluid moments are described by their evolution equations. The 21-moment model contains two "heat flux vectors" (3rd and 5th-order moments) and two "viscosity-tensors" (2nd and 4th-order moments). The Braginskii model is then obtained as a particular case of a one ion-electron plasma with similar temperatures, with de-coupled heat fluxes and viscosity-tensors expressed in a quasi-static approximation. We provide all the numerical values of the Braginskii model in a fully analytic form (together with the 4th and 5th-order moments). For multi-species plasmas, the model makes calculation of transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasi-static approximation which makes some Braginskii coefficients divergent in a weakly-collisional regime. Importantly, we show that the heat fluxes and viscosity-tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the 4th-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several Appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad.

preprint2022arXiv

Particle acceleration in an MHD-scale system of multiple current sheets

We investigate particle acceleration in an MHD-scale system of multiple current sheets by performing 2D and 3D MHD simulations combined with a test particle simulation. The system is unstable for the tearing-mode instability, and magnetic islands are produced by magnetic reconnection. Due to the interaction of magnetic islands, the system turns into a turbulent state. The 2D (3D) case yields both $-5/3$ ($-11/3$ and $-7/3$) power-law spacetra for magnetic and velocity fluctuations. Particles are efficiently energized by the generated turbulence, and it forms a power-law tail with an index of $-2.2$ and $-4.2$ in the energy distribution function for the 2D and 3D case, respectively. We find more energetic particles outside magnetic islands rather than inside. We observe super-diffusion in the 2D ($\sim t^{2.27}$) and 3D ($\sim t^{1.2}$) case in the energy space of energetic particles.

preprint2022arXiv

Turbulence in the Sub-Alfvénic Solar Wind

Parker Solar Probe (PSP) entered a region of the sub-Alfvenic solar wind during encounter 8 and we present the first detailed analysis of low-frequency turbulence properties in this novel region. The magnetic field and flow velocity vectors were highly aligned during this interval. By constructing spectrograms of the normalized magnetic helicity, cross helicity, and residual energy, we find that PSP observed primarily Alfvenic fluctuations, a consequence of the highly field-aligned flow that renders quasi-2D fluctuations unobservable to PSP. We extend Taylor hypothesis to sub and super Alfvenic flows. Spectra for the fluctuating forward and backward Elsasser variables are presented. The observed spectra are well fitted using a spectral theory for nearly incompressible magnetohydrodynamics assuming a wave number anisotropy. The density spectrum is a power law that resembles neither the Elsasser spectra nor the compressible magnetic field spectrum, suggesting that these are advected entropic rather than magnetosonic modes and not due to the parametric decay instability. Spectra in the neighboring modestly super Alfvenic intervals are similar.

preprint2021arXiv

Inertial range magnetic fluctuation anisotropy observed from PSP first seven orbits

Solar wind turbulence is anisotropic with respect to the mean magnetic field. Anisotropy leads to ambiguity when interpreting in-situ turbulence observations in the solar wind because an apparent change in the measurements could be due either to the change of intrinsic turbulence properties or to a simple change of the spacecraft sampling direction. We demonstrate the ambiguity using the spectral index and magnetic compressibility in the inertial range observed by the Parker Solar Probe during its first seven orbits ranging from 0.1 to 0.6 AU. To unravel the effects of the sampling direction, we assess whether the wavevector anisotropy is consistent with a two-dimensional (2D) plus slab turbulence transport model and determine the fraction of power in the 2D versus slab component. Our results confirm that the 2D plus slab model is consistent with the data and the power ratio between 2D and slab components depends on radial distance, with the relative power in 2D fluctuations becoming smaller closer to the Sun.

preprint2021arXiv

PSP/IS$\odot$IS Observation of a Solar Energetic Particle Event Associated With a Streamer Blowout Coronal Mass Ejection During Encounter 6

In this paper we examine a low-energy SEP event observed by IS$\odot$IS's Energetic Particle Instrument-Low (EPI-Lo) inside 0.18 AU on September 30, 2020. This small SEP event has a very interesting time profile and ion composition. Our results show that the maximum energy and peak in intensity is observed mainly along the open radial magnetic field. The event shows velocity dispersion, and strong particle anisotropies are observed throughout the event showing that more particles are streaming outward from the Sun. We do not see a shock in the in-situ plasma or magnetic field data throughout the event. Heavy ions, such as O and Fe were detected in addition to protons and 4He, but without significant enhancements in 3He or energetic electrons. Our analysis shows that this event is associated with a slow streamer-blowout coronal mass ejection (SBO-CME) and the signatures of this small CME event are consistent with those typical of larger CME events. The time-intensity profile of this event shows that PSP encountered the western flank of the SBO-CME. The anisotropic and dispersive nature of this event in a shockless local plasma give indications that these particles are most likely accelerated remotely near the Sun by a weak shock or compression wave ahead of the SBO-CME. This event may represent direct observations of the source of low-energy SEP seed particle population.

preprint2021arXiv

The development of a split-tail heliosphere and the role of non-ideal processes: a comparison of the BU and Moscow models

Global models of the heliosphere are critical tools used in the interpretation of heliospheric observations. There are several three-dimensional magnetohydrodynamic (MHD) heliospheric models that rely on different strategies and assumptions. Until now only one paper has compared global heliosphere models, but without magnetic field effects. We compare the results of two different MHD models, the BU and Moscow models. Both models use identical boundary conditions to compare how different numerical approaches and physical assumptions contribute to the heliospheric solution. Based on the different numerical treatments of discontinuities, the BU model allows for the presence of magnetic reconnection, while the Moscow model does not. Both models predict collimation of the solar outflow in the heliosheath by the solar magnetic field and produce a split-tail where the solar magnetic field confines the charged solar particles into distinct north and south columns that become lobes. In the BU model, the ISM flows between the two lobes at large distances due to MHD instabilities and reconnection. Reconnection in the BU model at the port flank affects the draping of the interstellar magnetic field in the immediate vicinity of the heliopause. Different draping in the models cause different ISM pressures, yielding different heliosheath thicknesses and boundary locations, with the largest effects at high latitudes. The BU model heliosheath is 15% thinner and the heliopause is 7% more inwards at the north pole relative to the Moscow model. These differences in the two plasma solutions may manifest themselves in energetic neutral atom measurements of the heliosphere.

preprint2019arXiv

A Brief Review on Particle Acceleration in Multi-island Magnetic Reconnection

The basic physics and recent progresses in theoretical and particle-in-cell (PIC) simulation studies of particle acceleration in multi-island magnetic reconnection are briefly reviewed. Particle acceleration in multi-island magnetic reconnection is considered a plausible mechanism for the acceleration of energetic particles in solar flares and the solar wind. Theoretical studies have demonstrated that such a mechanism can produce the observed power-law energy distribution of energetic particles if the particle motion is sufficiently randomized in the reconnection event. However, PIC simulations seem to suggest that the first-order Fermi acceleration mechanism is unable to produce a power-law particle energy distribution function in mildly relativistic multi-island magnetic reconnections. On the other hand, while simulations of highly relativistic reconnections appear to be able to produce a power-law energy spectrum, the spectral indices obtained are generally harder than the soft power-law spectra with indices $\sim -5$ commonly observed in the solar wind and solar flare events. In addition, the plasma heating due to kinetic instabilities in 3D magnetic reconnection may "thermalize" the power-law particles, making it even more difficult for multi-island reconnections to generate a power-law spectrum. We discuss the possible reasons that may lead to these problems.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 2 -- Kinetic theory, Padé approximants and Landau fluid closures

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(ζ)$, and the associated plasma response function $R(ζ)=1+ζZ(ζ)=-Z'(ζ)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Padé approximants of $R(ζ)$ up to the 8th-pole order, with many given in an analytic form.

preprint2019arXiv

Heliospheric structure as revealed by the 3 -- 88 keV H ENA spectra

Energetic neutral atoms (ENA) are an important tool for investigating the structure of the heliosphere. Recently, it was observed that fluxes of ENAs (with energy $\le$ 55 keV) coming from the upwind and downwind regions of the heliosphere are similar in strength. This led the authors of these observations to hypothesize that the heliosphere is bubble-like rather than comet-like, meaning that it has no extended tail. We investigate the directional distribution of the ENA flux for a wide energy range (3--88 keV) including the observations from IBEX (Interstellar Boundary Explorer), INCA (Ion and Neutral Camera, on board Cassini), and HSTOF (High energy Suprathermal Time Of Flight sensor, on board SOHO, Solar and Heliospheric Observatory). An essential element is the model of pickup ion acceleration at the termination shock (TS) proposed by Zank. We use state of the art models of the global heliosphere, interstellar neutral gas density, and pickup ion distributions. The results, based on the "comet-like" model of the heliosphere, are close in flux magnitude to ENA observations by IBEX, HSTOF and partly by INCA (except for the 5.2-13.5 keV energy channel). We find that the ENA flux from the tail dominates at high energy (in agreement with HSTOF, but not INCA). At low energy, our comet-like model produces the similar strengths of the ENA fluxes from the upwind and downwind directions, which, therefore, removes this as a compelling argument for a bubble-like heliosphere.

preprint2019arXiv

Identification of Magnetic Flux Ropes from Parker Solar Probe Observations during the First Encounter

The Parker Solar Probe (PSP) observed an interplanetary coronal mass ejection (ICME) event during its first orbit around the sun, among many other events. This event is analyzed by applying a wavelet analysis technique to obtain the reduced magnetic helicity, cross helicity, and residual energy, the first two of which are magnetohydrodynamics (MHD) invariants. Our results show that the ICME, as a large scale magnetic flux rope, possesses high magnetic helicity, very low cross helicity, and highly negative residual energy, thus pointing to a magnetic fluctuation dominated structure. Using the same technique, we also search for small-scale coherent magnetic flux rope structures during the period from 2018/10/22--2018/11/21, which are intrinsic to quasi-2D MHD turbulence in the solar wind. Multiple structures with duration between 8 and 300 minutes are identified from PSP in-situ spacecraft measurements. The location and scales of these structures are characterized by wavelet spectrograms of the normalized reduced magnetic helicity, normalized cross helicity and normalized residual energy. Transport theory suggests that these small-scale magnetic flux ropes may contribute to the acceleration of charged particles through magnetic reconnection processes, and the dissipation of these structures may be important for understanding the coronal heating processes.

preprint2019arXiv

Impact of Space Weather on Climate and Habitability of Terrestrial Type Exoplanets

The current progress in the detection of terrestrial type exoplanets has opened a new avenue in the characterization of exoplanetary atmospheres and in the search for biosignatures of life with the upcoming ground-based and space missions. To specify the conditions favorable for the origin, development and sustainment of life as we know it in other worlds, we need to understand the nature of astrospheric, atmospheric and surface environments of exoplanets in habitable zones around G-K-M dwarfs including our young Sun. Global environment is formed by propagated disturbances from the planet-hosting stars in the form of stellar flares, coronal mass ejections, energetic particles, and winds collectively known as astrospheric space weather. Its characterization will help in understanding how an exoplanetary ecosystem interacts with its host star, as well as in the specification of the physical, chemical and biochemical conditions that can create favorable and/or detrimental conditions for planetary climate and habitability along with evolution of planetary internal dynamics over geological timescales. A key linkage of (astro) physical, chemical, and geological processes can only be understood in the framework of interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary and Earth sciences. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets will significantly expand the current definition of the habitable zone to the biogenic zone and provide new observational strategies for searching for signatures of life. The major goal of this paper is to describe and discuss the current status and recent progress in this interdisciplinary field and to provide a new roadmap for the future development of the emerging field of exoplanetary science and astrobiology.

preprint2019arXiv

PIC simulation of a shock tube: Implications for wave transmission in the heliospheric boundary region

A shock tube problem is solved numerically by using one-dimensional full particle-in-cell simulations under the condition that a relatively tenuous and weakly magnetized plasma is continuously pushed by a relatively dense and strongly magnetized plasma having supersonic relative velocity. A forward and a reverse shock and a contact discontinuity are self-consistently reproduced. The spatial width of the contact discontinuity increases as the angle between the discontinuity normal and ambient magnetic field decreases. The inner structure of the discontinuity shows different profiles between magnetic field and plasma density, or pressure, which is caused by a non-MHD effect of the local plasma. The region between the two shocks is turbulent. The fluctuations in the relatively dense plasma are compressible and propagating away from the contact discontinuity, although the fluctuations in the relatively tenuous plasma contain both compressible and incompressible components. The source of the compressible fluctuations in the relatively dense plasma is in the relatively tenuous plasma. Only compressible fast mode fluctuations generated in the relatively tenuous plasma are transmitted through the contact discontinuity and propagate in the relatively dense plasma. These fast mode fluctuations are steepened when passing the contact discontinuity. This wave steepening and probably other effects may cause the broadening of the wave spectrum in the very local interstellar medium plasma. The results are discussed in the context of the heliospheric boundary region or heliopause.

preprint2019arXiv

Predicting the Solar Wind at Parker Solar Probe Using an Empirically Driven MHD Model

Since the launch on 2018/08/12, Parker Solar Probe (PSP) has completed its first and second orbits around the Sun, having reached down to 35.7 solar radii at each perihelion. In anticipation of the exciting new data at such unprecedented distances, we have simulated the global 3D heliosphere using an MHD model coupled with a semi-empirical coronal model using the best available photospheric magnetograms as input. We compare our heliospheric MHD simulation results with in situ measurements along the PSP trajectory from its launch to the completion of the second orbit, with particular emphasis on the solar wind structure around the first two solar encounters. Furthermore, we show our model prediction for the third perihelion, which occurred on 2019/09/01. Comparison of the MHD results with PSP observations provides a new insight on the solar wind acceleration. Moreover, PSP observations reveal how accurately the ADAPT-WSA predictions work throughout the inner heliosphere.

preprint2010arXiv

From Micro- to Macro-scales in the Heliosphere and Magnetospheres

From a broader perspective, the heliosphere and planetary magnetospheres provide a test bed to explore the plasma physics of the Universe. In particular, the underlying nonlinear coupling of different spatial and temporal scales plays a key role in determining the structure and dynamics of space plasmas and electromagnetic fields. Plasmas and fields exhibit both laminar and turbulent properties, corresponding to either well organized or disordered states, and the development of quantitative theoretical and analytical descriptions from physics based first principles is a profound challenge. Limited observations and complications introduced by geometry and physical parameters conspire to complicate the problem. Dimensionless scaling analysis and statistical methods are universally applied common approaches that allow for the application of related ideas to multiple physical problems. We discuss several examples of the interplay between the scales in a variety of space plasma environments, as exemplified in the presentations of the session \textit{From Micro- to Macro-scales in the Heliosphere and Magnetospheres.}

preprint2008arXiv

The Heliosphere in Time

Because of the dynamic nature of the interstellar medium, the Sun should have encountered a variety of different interstellar environments in its lifetime. As the solar wind interacts with the surrounding interstellar medium to form a heliosphere, different heliosphere shapes, sizes, and particle contents result from the different environments. Some of the large possible interstellar parameter space (density, velocity, temperature) is explored here with the help of global heliosphere models, and the features in the resulting heliospheres are compared and discussed. The heliospheric size, expressed as distance of the nose of the heliopause to the Sun, is set by the solar wind - interstellar pressure balance, even for extreme cases. Other heliospheric boundary locations and neutral particle results correlate with the interstellar parameters as well. If the H0 clouds identified in the Millennium Arecibo survey are typical of clouds encountered by the Sun, then the Sun spends ~99.4% of the time in warm low density ISM, where the typical upwind heliosphere radii are up to two orders of magnitude larger than at present.