Source author record

G. P. Zank

G. P. Zank 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

33works
11topics
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

33 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.

preprint2016arXiv

Modeling the Solar Wind at the Ulysses, Voyager, and New Horizons Spacecraft

The outer heliosphere is a dynamic region shaped largely by the interaction between the solar wind and the interstellar medium. While interplanetary magnetic field and plasma observations by the Voyager spacecraft have significantly improved our understanding of this vast region, modeling the outer heliosphere still remains a challenge. We simulate the three-dimensional, time-dependent solar wind flow from 1 to 80 astronomical units (AU), where the solar wind is assumed to be supersonic, using a two-fluid model in which protons and interstellar neutral hydrogen atoms are treated as separate fluids. We use 1-day averages of the solar wind parameters from the OMNI data set as inner boundary conditions to reproduce time-dependent effects in a simplified manner which involves interpolation in both space and time. Our model generally agrees with Ulysses data in the inner heliosphere and Voyager data in the outer heliosphere. Ultimately, we present the model solar wind parameters extracted along the trajectory of New Horizons spacecraft. We compare our results with in situ plasma data taken between 11 and 33 AU and at the closest approach to Pluto on July 14, 2015.

preprint2015arXiv

Multi-Symplectic Magnetohydrodynamics: II, Addendum and Erratum

A recent paper arXiv:1312.4890 on multi-symplectic magnetohydrodynamics (MHD) using Clebsch variables in an Eulerian action principle with constraints is further extended. We relate a class of symplecticity conservation laws to a vorticity conservation law, and provide a corrected form of the Poincaré-Cartan differential form formulation of the system. We also correct some typographical errors (omissions) in arXiv:1312.4890. We show that the vorticity-symplecticity conservation law, that arises as a compatibility condition on the system, expressed in terms of the Clebsch variables is equivalent to taking the curl of the conservation form of the MHD momentum equation. We use the Cartan-Poincaré form to obtain a class of differential forms that represent the system using Cartan's geometric theory of partial differential equations.

preprint2015arXiv

Small-scale magnetic islands in the solar wind and their role in particle acceleration. Part 1: Dynamics of magnetic islands near the heliospheric current sheet

Increases of ion fluxes in the keV-MeV range are sometimes observed near the heliospheric current sheet (HCS) during periods when other sources are absent. These resemble solar energetic particle (SEP) events, but the events are weaker and apparently local. Conventional explanations based on either shock acceleration of charged particles or particle acceleration due to magnetic reconnection at interplanetary current sheets are not persuasive. We suggest instead that recurrent magnetic reconnection occurs at the HCS and smaller current sheets in the solar wind (Zharkova & Khabarova 2012), of which a consequence is particle energization by the dynamically evolving secondary current sheets and magnetic islands (Zank et al. 2014; Drake et al. 2006a). The effectiveness of the trapping and acceleration process associated with magnetic islands depends in part on the topology of the HCS. We show that the HCS possesses ripples superimposed on the large-scale flat or wavy structure. We conjecture that the ripples can efficiently confine plasma and provide tokamak-like conditions that are favorable for the appearance of small-scale magnetic islands that merge and/or contract. Particles trapped in the vicinity of merging islands and experiencing multiple small-scale reconnection events are accelerated by the induced electric field, and experience first-order Fermi acceleration in contracting magnetic islands (Zank et al. 2014). We present multi-spacecraft observations of magnetic island merging and particle energization in the absence of other sources, providing support for theory and simulations that show particle energization by reconnection related processes of magnetic island merging and contraction.

preprint2014arXiv

Local and Nonlocal Advected Invariants and Helicities in Magnetohydrodynamics and Gas Dynamics II: Noether's Theorems and Casimirs

Conservation laws in ideal gas dynamics and magnetohydrodynamics (MHD) associated with fluid relabelling symmetries are derived using Noether's first and second theorems. Lie dragged invariants are discussed in terms of the MHD Casimirs. A nonlocal conservation law for fluid helicity applicable for a non-barotropic fluid involving Clebsch variables is derived using Noether's theorem, in conjunction with a fluid relabelling symmetry and a gauge transformation. A nonlocal cross helicity conservation law involving Clebsch potentials, and the MHD energy conservation law are derived by the same method. An Euler Poincaré variational approach is also used to derive conservation laws associated with fluid relabelling symmetries using Noether's second theorem.

preprint2013arXiv

Characteristic Length of Energy-Containing Structures at the Base of a Coronal Hole

An essential parameter for models of coronal heating and fast solar wind acceleration that rely on the dissipation of MHD turbulence is the characteristic energy-containing length $λ_{\bot}$ of the squared velocity and magnetic field fluctuations ($u^2$ and $b^2$) transverse to the mean magnetic field inside a coronal hole (CH) at the base of the corona. The characteristic length scale defines directly the heating rate. We use a time series analysis of solar granulation and magnetic field measurements inside two CHs obtained with the New Solar Telescope (NST) at Big Bear Solar Observatory. A data set for transverse magnetic fields obtained with the Solar Optical Telescope/Spectro-Polarimeter (SOT/SP) aboard {\it Hinode} spacecraft was utilized to analyze the squared transverse magnetic field fluctuations $b_t^2$. Local correlation tracking (LCT) was applied to derive the squared transverse velocity fluctuations $u^2$. We find that for $u^2$-structures, Batchelor integral scale $λ$ varies in a range of 1800 - 2100 km, whereas the correlation length $ς$ and the $e$-folding length $L$ vary between 660 and 1460 km. Structures for $b_t^2$ yield $λ\approx 1600$ km, $ς\approx 640$ km, and $L \approx 620$ km. An averaged (over $λ, ς$, and $L$) value of the characteristic length of $u^2$-fluctuations is 1260$\pm$500 km, and that of $b_t^2$ is 950$\pm$560 km. The characteristic length scale in the photosphere is approximately 1.5-50 times smaller than that adopted in previous models (3-30$\times10^3$ km). Our results provide a critical input parameter for current models of coronal heating and should yield an improved understanding of fast solar wind acceleration.

preprint2010arXiv

Can IBEX Identify Variations in the Galactic Environment of the Sun using Energetic Neutral Atom (ENAs)?

The Interstellar Boundary Explorer (IBEX) spacecraft is providing the first all-sky maps of the energetic neutral atoms (ENAs) produced by charge-exchange between interstellar neutral \HI\ atoms and heliospheric solar wind and pickup ions in the heliosphere boundary regions. The 'edge' of the interstellar cloud presently surrounding the heliosphere extends less than 0.1 pc in the upwind direction, terminating at an unknown distance, indicating that the outer boundary conditions of the heliosphere could change during the lifetime of the IBEX satellite. Using reasonable values for future outer heliosphere boundary conditions, ENA fluxes are predicted for one possible source of ENAs coming from outside of the heliopause. The ENA production simulations use three-dimensional MHD plasma models of the heliosphere that include a kinetic description of neutrals and a Lorentzian distribution for ions. Based on this ENA production model, it is then shown that the sensitivities of the IBEX 1.1 keV skymaps are sufficient to detect the variations in ENA fluxes that are expected to accompany the solar transition into the next upwind cloud. Approximately 20% of the IBEX 1.1 keV pixels appear capable of detecting the predicted model differences at the $ 3 σ$ level, with these pixels concentrated in the Ribbon region. Regardless of the detailed ENA production model, the success of the modeled \BdotR\ directions in reproducing the Ribbon locus, together with our results, indicate that the Ribbon phenomenon traces the variations in the heliosphere distortion caused by the relative pressures of the interstellar magnetic and gaseous components.

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.}

preprint2010arXiv

Modulation of waves due to charge-exchange collisions in magnetized partially ionized space plasma

A nonlinear time dependent fluid simulation model is developed that describes the evolution of magnetohydrodynamic waves in the presence of collisional and charge exchange interactions of a partially ionized plasma. The partially ionized plasma consists of electrons, ions and a significant number of neutral atoms. In our model, the electrons and ions are described by a single fluid compressible magnetohydrodynamic (MHD) model and are coupled self-consistently to the neutral gas, described by the compressible hydrodynamic equations. Both the plasma and neutral fluids are treated with different energy equations that describe thermal energy exchange processes between them. Based on our self-consistent model, we find that propagating Alfvénic and fast/slow modes grow and damp alternately through a nonlinear modulation process. The modulation appears to be robust and survives strong damping by the neutral component.

preprint2009arXiv

3D Hall MHD Modeling of Solar Wind Plasma Spectra

We present fully self consistent 3D simulations of compressible Hall MHD plasma that describe spectral features relevant to the solar wind plasma. We find that a $k^{-7/3}$ spectrum sets in for the fluctuations that are smaller than ion gyro radius. We further investigate scale dependent anisotropy led by nonlinear processes relevant to the solar wind plasma. Our work is important particularly in understanding the role of wave and nonlinear cascades in the evolution of the solar wind, structure formation at the largest scales.

preprint2009arXiv

Exclusion of Tiny Interstellar Dust Grains from the Heliosphere

The distribution of interstellar dust grains (ISDG) observed in the Solar System depends on the nature of the interstellar medium-solar wind interaction. The charge of the grains couples them to the interstellar magnetic field (ISMF) resulting in some fraction of grains being excluded from the heliosphere while grains on the larger end of the size distribution, with gyroradii comparable to the size of the heliosphere, penetrate the termination shock. This results in a skewing the size distribution detected in the Solar System. We present new calculations of grain trajectories and the resultant grain density distribution for small ISDGs propagating through the heliosphere. We make use of detailed heliosphere model results, using three-dimensional (3-D) magnetohydrodynamic/kinetic models designed to match data on the shape of the termination shock and the relative deflection of interstellar neutral H and He flowing into the heliosphere. We find that the necessary inclination of the ISMF relative to the inflow direction results in an asymmetry in the distribution of the larger grains (0.1 micron) that penetrate the heliopause. Smaller grains (0.01 micron) are completely excluded from the Solar System at the heliopause.

preprint2009arXiv

Self-consistent Simulations of Plasma-Neutral in a Partially Ionized Astrophysical Turbulent Plasma

A local turbulence model is developed to study energy cascades in the heliosheath and outer heliosphere (OH) based on self-consistent two-dimensional fluid simulations. The model describes a partially ionized magnetofluid OH that couples a neutral hydrogen fluid with a plasma primarily through charge-exchange interactions. Charge-exchange interactions are ubiquitous in warm heliospheric plasma, and the strength of the interaction depends largely on the relative speed between the plasma and the neutral fluid. Unlike small-length scale linear collisional dissipation in a single fluid, charge-exchange processes introduce channels that can be effective on a variety of length scales that depend on the neutral and plasma densities, temperature, relative velocities, charge-exchange cross section, and the characteristic length scales. We find, from scaling arguments and nonlinear coupled fluid simulations, that charge-exchange interactions modify spectral transfer associated with large-scale energy-containing eddies. Consequently, the turbulent cascade rate prolongs spectral transfer among inertial range turbulent modes. Turbulent spectra associated with the neutral and plasma fluids are therefore steeper than those predicted by Kolmogorov's phenomenology. Our work is important in the context of the global heliospheric interaction, the energization and transport of cosmic rays, gamma-ray bursts, interstellar density spectra, etc. Furthermore, the plasma-neutral coupling is crucial in understanding the energy dissipation mechanism in molecular clouds and star formation processes.

preprint2009arXiv

Self-Organization in a Driven Dissipative Plasma System

We perform a fully self-consistent 3-D numerical simulation for a compressible, dissipative magneto-plasma driven by large-scale perturbations, that contain a fairly broader spectrum of characteristic modes, ranging from largest scales to intermediate scales and down to the smallest scales, where the energy of the system are dissipated by collisional (Ohmic) and viscous dissipations. Additionally, our simulation includes nonlinear interactions amongst a wide range of fluctuations that are initialized with random spectral amplitudes, leading to the cascade of spectral energy in the inertial range spectrum, and takes into account large scale as well as small scale perturbation that may have been induced by the background plasma fluctuations, also the non adiabatic exchange of energy leading to the migration of energy from the energy containing modes or randomly injected energy driven by perturbations and further dissipated by the smaller scales. Besides demonstrating the comparative decays of total energy and dissipation rate of energy, our results show the existence of a perpendicular component of current, thus clearly confirming that the self-organized state is non-force free.

preprint2009arXiv

Spectral features of solar wind turbulent plasma

Spectral properties of a fully compressible solar wind Hall Magnetohydrodynamic plasma are investigated by means of time dependent three dimensional Hall MHD simulations. Our simulations, in agreement with spacecraft data, identify a spectral break in turbulence spectra at characteristic length-scales associated with electromagnetic fluctuations that are smaller than the ion gyroradius. In this regime, our 3D simulations show that turbulent spectral cascades in the presence of a mean magnetic field follow an omnidirectional anisotropic inertial range spectrum close to $k^{-7/3}$. The onset of the spectral break in our simulations can be ascribed to the presence of nonlinear Hall interactions that modify the spectral cascades. Our simulations further show that the underlying charachteristic turbulent fluctuations are spectrally anisotropic, the extent of which depends critically on the local wavenumber. The fluctuations associated with length scales smaller than the ion gyroradius are highly compressible and tend to exhibit a near equipartition in the velocity and magnetic fields. Finally, we find that the orientation of velocity and magnetic field fluctuations critically determine the character of nonlinear interactions that predominantly govern a Hall MHD plasma, like the solar wind.

preprint2009arXiv

Turbulent Density Spectrum in Solar Wind Plasma

The density fluctuation spectrum in the solar wind reveals a Kolmogorov-like scaling with a spectral slope of -5/3 in wavenumber space. The energy transfer process in the magnetized solar wind, characterized typically by MHD turbulence, over extended length-scales remains an unresolved paradox of modern turbulence theories, raising the question of how a compressible magnetofluid exhibits a turbulent spectrum that is characteristic of an incompressible hydrodynamic fluid. To address these questions, we have undertaken three-dimensional time dependent numerical simulations of a compressible magnetohydrodynamic fluid describing super-Alfvénic, supersonic and strongly magnetized plasma fluid. It is shown that a Kolmogorov-like density spectrum can develop by plasma motions that are dominated by Alfvénic cascades whereas compressive modes are dissipated.

preprint2009arXiv

Turbulent Spectra in the Solar Wind Plasma

Observations of interstellar scintillations at radio wavelengths reveal a Kolmogorov-like scaling of the electron density spectrum with a spectral slope of -5/3 over six decades in wavenumber space. A similar turbulent density spectrum in the solar wind plasma has been reported. The energy transfer process in the magnetized solar wind plasma over such extended length-scales remains an unresolved paradox of modern turbulence theories raising the especially intriguing question of how a compressible magnetized solar wind exhibits a turbulent spectrum that is a characteristic of an incompressible hydrodynamic fluid. To address these questions, we have undertaken three-dimensional time dependent numerical simulations of a compressible magnetohydrodynamic fluid describing super-Alfvénic, supersonic and strongly magnetized plasma. It is shown that the observed Kolmogorov-like (-5/3) spectrum can develop in the solar wind plasma by supersonic plasma motions that dissipate into highly subsonic motion that passively convect density fluctuations.

preprint2009arXiv

Whistler Wave Turbulence in Solar Wind Plasma

Whistler waves are present in solar wind plasma. These waves possess characteristic turbulent fluctuations that are characterized typically by the frequency and length scales that are respectively bigger than ion gyro frequency and smaller than ion gyro radius. The electron inertial length is an intrinsic length scale in whistler wave turbulence that distinguishably divides the high frequency solar wind turbulent spectra into scales smaller and bigger than the electron inertial length. We present nonlinear three dimensional, time dependent, fluid simulations of whistler wave turbulence to investigate their role in solar wind plasma. Our simulations find that the dispersive whistler modes evolve entirely differently in the two regimes. While the dispersive whistler wave effects are stronger in the large scale regime, they do not influence the spectral cascades which are describable by a Kolmogorov-like $k^{-7/3}$ spectrum. By contrast, the small scale turbulent fluctuations exhibit a Navier-Stokes like evolution where characteristic turbulent eddies exhibit a typical $k^{-5/3}$ hydrodynamic turbulent spectrum. By virtue of equipartition between the wave velocity and magnetic fields, we quantify the role of whistler waves in the solar wind plasma fluctuations.

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.