Showing 1–50 of 59 results for author: Zank, G

Characterization of Turbulent Fluctuations in the Sub-Alfvenic Solar Wind

Authors: Gary P. Zank, Lingling Zhao, Laxman Adhikari, Daniele Telloni, Prashant Baruwal, Prashrit Baruwal, Xingyu Zhu, Masaru Nakanotani, Alexander Pitna, Justin C. Kasper, Stuart D. Bale

Abstract: Parker Solar Probe (PSP) observed sub-Alfvenic solar wind intervals during encounters 8 - 14, and low-frequency magnetohydrodynamic turbulence in these regions may differ from that in super-Alfvenic wind. We apply a new mode-decomposition analysis (Zank et al 2023) to the sub-Alfvénic flow observed by PSP on 2021 April 28, identifying and characterizing entropy, magnetic islands, forward and backw… ▽ More Parker Solar Probe (PSP) observed sub-Alfvenic solar wind intervals during encounters 8 - 14, and low-frequency magnetohydrodynamic turbulence in these regions may differ from that in super-Alfvenic wind. We apply a new mode-decomposition analysis (Zank et al 2023) to the sub-Alfvénic flow observed by PSP on 2021 April 28, identifying and characterizing entropy, magnetic islands, forward and backward Alfvén waves, including weakly/non-propagating Alfvén vortices, forward and backward fast and slow magnetosonic modes. Density fluctuations are primarily and almost equally entropy and backward propagating slow magnetosonic modes. The mode-decomposition provides phase information (frequency and wavenumber k) for each mode. Entropy-density fluctuations have a wavenumber anisotropy k_{||} >> k_{perp} whereas slow mode density fluctuations have k_{perp} > k_{||}. Magnetic field fluctuations are primarily magnetic island modes (delta B^i) with an O(1) smaller contribution from uni-directionally propagating Alfven waves (delta B^{A+}) giving a variance anisotropy of <{δB^i}^2> / <delta B^A}^2> = 4.1. Incompressible magnetic fluctuations dominate compressible contributions from fast and slow magnetosonic modes. The magnetic island spectrum is Kolmogorov-like k_{perp}^{-1.6} in perpendicular wavenumber and the uni-directional Alfven wave spectra are k_{||}^{-1.6} and k_{perp}^{-1.5}. Fast magnetosonic modes propagate at essentially the Alfvén speed with anti-correlated transverse velocity and magnetic field fluctuations and are almost exclusively magnetic due to beta_p<<1. Transverse velocity fluctuations are the dominant velocity component in fast magnetosonic modes and longitudinal fluctuations dominate in slow modes. Mode-decomposition is an effective tool in identifying the basic building blocks of MHD turbulence and provides detailed phase information about each of the modes. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Properties of an interplanetary shock observed at 0.07 and 0.7 Astronomical Units by Parker Solar Probe and Solar Orbiter

Authors: D. Trotta, A. Larosa, G. Nicolaou, T. S. Horbury, L. Matteini, H. Hietala, X. Blanco-Cano, L. Franci, C. H. K. Chen, L. Zhao, G. P. Zank, C. M. S. Cohen, S. D. Bale, R. Laker, N. Fargette, F. Valentini, Y. Khotyaintsev, R. Kieokaew, N. Raouafi, E. Davies, R. Vainio, N. Dresing, E. Kilpua, T. Karlsson, C. J. Owen , et al. (1 additional authors not shown)

Abstract: The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with… ▽ More The Parker Solar Probe (PSP) and Solar Orbiter (SolO) missions opened a new observational window in the inner heliosphere, which is finally accessible to direct measurements. On September 05, 2022, a coronal mass ejection (CME)-driven interplanetary (IP) shock has been observed as close as 0.07 au by PSP. The CME then reached SolO, which was well radially-aligned at 0.7 au, thus providing us with the opportunity to study the shock properties at so different heliocentric distances. We characterize the shock, investigate its typical parameters and compare its small-scale features at both locations. Using the PSP observations, we investigate how magnetic switchbacks and ion cyclotron waves are processed upon shock crossing. We find that switchbacks preserve their V--B correlation while compressed upon the shock passage, and that the signature of ion cyclotron waves disappears downstream of the shock. By contrast, the SolO observations reveal a very structured shock transition, with a population of shock-accelerated protons of up to about 2 MeV, showing irregularities in the shock downstream, which we correlate with solar wind structures propagating across the shock. At SolO, we also report the presence of low-energy ($\sim$ 100 eV) electrons scattering due to upstream shocklets. This study elucidates how the local features of IP shocks and their environments can be very different as they propagate through the heliosphere. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: In review in ApJ

Coronal Heating Rate in the Slow Solar Wind

Authors: Daniele Telloni, Marco Romoli, Marco Velli, Gary P. Zank, Laxman Adhikari, Cooper Downs, Aleksandr Burtovoi, Roberto Susino, Daniele Spadaro, Lingling Zhao, Alessandro Liberatore, Chen Shi, Yara De Leo, Lucia Abbo, Federica Frassati, Giovanna Jerse, Federico Landini, Gianalfredo Nicolini, Maurizio Pancrazzi, Giuliana Russano, Clementina Sasso, Vincenzo Andretta, Vania Da Deppo, Silvano Fineschi, Catia Grimani , et al. (37 additional authors not shown)

Abstract: This Letter reports the first observational estimate of the heating rate in the slowly expanding solar corona. The analysis exploits the simultaneous remote and local observations of the same coronal plasma volume with the Solar Orbiter/Metis and the Parker Solar Probe instruments, respectively, and relies on the basic solar wind magnetohydrodynamic equations. As expected, energy losses are a mino… ▽ More This Letter reports the first observational estimate of the heating rate in the slowly expanding solar corona. The analysis exploits the simultaneous remote and local observations of the same coronal plasma volume with the Solar Orbiter/Metis and the Parker Solar Probe instruments, respectively, and relies on the basic solar wind magnetohydrodynamic equations. As expected, energy losses are a minor fraction of the solar wind energy flux, since most of the energy dissipation that feeds the heating and acceleration of the coronal flow occurs much closer to the Sun than the heights probed in the present study, which range from 6.3 to 13.3 solar radii. The energy deposited to the supersonic wind is then used to explain the observed slight residual wind acceleration and to maintain the plasma in a non-adiabatic state. As derived in the Wentzel-Kramers-Brillouin limit, the present energy transfer rate estimates provide a lower limit, which can be very useful in refining the turbulence-based modeling of coronal heating and subsequent solar wind acceleration. △ Less

Submitted 19 June, 2023; originally announced June 2023.

HelioSwarm: A Multipoint, Multiscale Mission to Characterize Turbulence

Authors: Kristopher G. Klein, Harlan Spence, Olga Alexandrova, Matthew Argall, Lev Arzamasskiy, Jay Bookbinder, Theodore Broeren, Damiano Caprioli, Anthony Case, Benjamin Chandran, Li-Jen Chen, Ivan Dors, Jonathan Eastwood, Colin Forsyth, Antoinette Galvin, Vincent Genot, Jasper Halekas, Michael Hesse, Butler Hine, Tim Horbury, Lan Jian, Justin Kasper, Matthieu Kretzschmar, Matthew Kunz, Benoit Lavraud , et al. (25 additional authors not shown)

Abstract: HelioSwarm (HS) is a NASA Medium-Class Explorer mission of the Heliophysics Division designed to explore the dynamic three-dimensional mechanisms controlling the physics of plasma turbulence, a ubiquitous process occurring in the heliosphere and in plasmas throughout the universe. This will be accomplished by making simultaneous measurements at nine spacecraft with separations spanning magnetohydr… ▽ More HelioSwarm (HS) is a NASA Medium-Class Explorer mission of the Heliophysics Division designed to explore the dynamic three-dimensional mechanisms controlling the physics of plasma turbulence, a ubiquitous process occurring in the heliosphere and in plasmas throughout the universe. This will be accomplished by making simultaneous measurements at nine spacecraft with separations spanning magnetohydrodynamic and sub-ion spatial scales in a variety of near-Earth plasmas. In this paper, we describe the scientific background for the HS investigation, the mission goals and objectives, the observatory reference trajectory and instrumentation implementation before the start of Phase B. Through multipoint, multiscale measurements, HS promises to reveal how energy is transferred across scales and boundaries in plasmas throughout the universe. △ Less

Submitted 10 June, 2023; originally announced June 2023.

Comments: Submitted to Space Science Reviews, 60 pages, 15 figures, 2 videos

Does Turbulence along the Coronal Current Sheet Drive Ion Cyclotron Waves?

Authors: Daniele Telloni, Gary P. Zank, Laxman Adhikari, Lingling Zhao, Roberto Susino, Ester Antonucci, Silvano Fineschi, Marco Stangalini, Catia Grimani, Luca Sorriso-Valvo, Daniel Verscharen, Raffaele Marino, Silvio Giordano, Raffaella D'Amicis, Denise Perrone, Francesco Carbone, Alessandro Liberatore, Roberto Bruno, Gaetano Zimbardo, Marco Romoli, Vincenzo Andretta, Vania Da Deppo, Petr Heinzel, John D. Moses, Giampiero Naletto , et al. (10 additional authors not shown)

Abstract: Evidence for the presence of ion cyclotron waves, driven by turbulence, at the boundaries of the current sheet is reported in this paper. By exploiting the full potential of the joint observations performed by Parker Solar Probe and the Metis coronagraph on board Solar Orbiter, local measurements of the solar wind can be linked with the large-scale structures of the solar corona. The results sugge… ▽ More Evidence for the presence of ion cyclotron waves, driven by turbulence, at the boundaries of the current sheet is reported in this paper. By exploiting the full potential of the joint observations performed by Parker Solar Probe and the Metis coronagraph on board Solar Orbiter, local measurements of the solar wind can be linked with the large-scale structures of the solar corona. The results suggest that the dynamics of the current sheet layers generates turbulence, which in turn creates a sufficiently strong temperature anisotropy to make the solar-wind plasma unstable to anisotropy-driven instabilities such as the Alfvén ion-cyclotron, mirror-mode, and firehose instabilities. The study of the polarization state of high-frequency magnetic fluctuations reveals that ion cyclotron waves are indeed present along the current sheet, thus linking the magnetic topology of the remotely imaged coronal source regions with the wave bursts observed in situ. The present results may allow improvement of state-of-the-art models based on the ion cyclotron mechanism, providing new insights into the processes involved in coronal heating. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 12 pages

Journal ref: The Astrophysical Journal, 2023

Coronal Loop Heating by Nearly Incompressible Magnetohydrodynamic and Reduced Magnetohydrodynamic Turbulence Models

Authors: Mehmet Sarp Yalim, Gary P. Zank, Mahboubeh Asgari-Targhi

Abstract: The transport of waves and turbulence beyond the photosphere is central to the coronal heating problem. Turbulence in the quiet solar corona has been modeled on the basis of the nearly incompressible magnetohydrodynamic (NI MHD) theory to describe the transport of low-frequency turbulence in open magnetic field regions. It describes the evolution of the coupled majority quasi-2D and minority slab… ▽ More The transport of waves and turbulence beyond the photosphere is central to the coronal heating problem. Turbulence in the quiet solar corona has been modeled on the basis of the nearly incompressible magnetohydrodynamic (NI MHD) theory to describe the transport of low-frequency turbulence in open magnetic field regions. It describes the evolution of the coupled majority quasi-2D and minority slab component, driven by the magnetic carpet and advected by a subsonic, sub-Alfvenic flow from the lower corona. In this paper, we couple the NI MHD turbulence transport model with an MHD model of the solar corona to study the heating problem in a coronal loop. In a realistic benchmark coronal loop problem, we find that a loop can be heated to ~1.5 million K by transport and dissipation of MHD turbulence described by the NI MHD model. We also find that the majority 2D component is as important as the minority slab component in the heating of the coronal loop. We compare our coupled MHD/NI MHD model results with a reduced MHD (RMHD) model. An important distinction between these models is that RMHD solves for small-scale velocity and magnetic field fluctuations and obtains the actual viscous/resistive dissipation associated with their evolution whereas NI MHD evolves scalar moments of the fluctuating velocity and magnetic fields and approximates dissipation using an MHD turbulence phenomenology. Despite the basic differences between the models, their simulation results match remarkably well, yielding almost identical heating rates inside the corona. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: 23 pages, 6 figures, 1 table, The Astrophysical Journal, in press

Space Plasma Physics: A Review

Authors: Bruce T. Tsurutani, Gary P. Zank, Veerle J. Sterken, Kazunari Shibata, Tsugunobu Nagai, Anthony J. Mannucci, David M. Malaspina, Gurbax S. Lakhina, Shrikanth G. Kanekal, Keisuke Hosokawa, Richard B. Horne, Rajkumar Hajra, Karl-Heinz Glassmeier, C. Trevor Gaunt, Peng-Fei Chen, Syun-Ichi Akasofu

Abstract: Owing to the ever-present solar wind, our vast solar system is full of plasmas. The turbulent solar wind, together with sporadic solar eruptions, introduces various space plasma processes and phenomena in the solar atmosphere all the way to the Earth's ionosphere and atmosphere and outward to interact with the interstellar media to form the heliopause and termination shock. Remarkable progress has… ▽ More Owing to the ever-present solar wind, our vast solar system is full of plasmas. The turbulent solar wind, together with sporadic solar eruptions, introduces various space plasma processes and phenomena in the solar atmosphere all the way to the Earth's ionosphere and atmosphere and outward to interact with the interstellar media to form the heliopause and termination shock. Remarkable progress has been made in space plasma physics in the last 65 years, mainly due to sophisticated in-situ measurements of plasmas, plasma waves, neutral particles, energetic particles, and dust via space-borne satellite instrumentation. Additionally high technology ground-based instrumentation has led to new and greater knowledge of solar and auroral features. As a result, a new branch of space physics, i.e., space weather, has emerged since many of the space physics processes have a direct or indirect influence on humankind. After briefly reviewing the major space physics discoveries before rockets and satellites, we aim to review all our updated understanding on coronal holes, solar flares and coronal mass ejections, which are central to space weather events at Earth, solar wind, storms and substorms, magnetotail and substorms, emphasizing the role of the magnetotail in substorm dynamics, radiation belts/energetic magnetospheric particles, structures and space weather dynamics in the ionosphere, plasma waves, instabilities, and wave-particle interactions, long-period geomagnetic pulsations, auroras, geomagnetically induced currents (GICs), planetary magnetospheres and solar/stellar wind interactions with comets, moons and asteroids, interplanetary discontinuities, shocks and waves, interplanetary dust, space dusty plasmas and solar energetic particles and shocks, including the heliospheric termination shock. This paper is aimed to provide a panoramic view of space physics and space weather. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: Accepted for publication in IEEE Transactions on Plasma Science (2022)

Journal ref: IEEE Transactions on Plasma Science (2022)

Mixing Interstellar Clouds Surrounding the Sun

Authors: Paweł Swaczyna, Nathan A. Schwadron, Eberhard Möbius, Maciej Bzowski, Priscilla C. Frisch, Jeffrey L. Linsky, David J. McComas, Fatemeh Rahmanifard, Seth Redfield, Réka M. Winslow, Brian E. Wood, Gary P. Zank

Abstract: On its journey through the Galaxy, the Sun passes through diverse regions of the interstellar medium. High-resolution spectroscopic measurements of interstellar absorption lines in spectra of nearby stars show absorption components from more than a dozen warm partially ionized clouds within 15 pc of the Sun. The two nearest clouds - the Local Interstellar Cloud (LIC) and Galactic (G) cloud - move… ▽ More On its journey through the Galaxy, the Sun passes through diverse regions of the interstellar medium. High-resolution spectroscopic measurements of interstellar absorption lines in spectra of nearby stars show absorption components from more than a dozen warm partially ionized clouds within 15 pc of the Sun. The two nearest clouds - the Local Interstellar Cloud (LIC) and Galactic (G) cloud - move toward each other. Their bulk heliocentric velocities can be compared with the interstellar neutral helium flow velocity obtained from space-based experiments. We combine recent results from Ulysses, IBEX, and STEREO observations to find a more accurate estimate of the velocity and temperature of the very local interstellar medium. We find that, contrary to the widespread viewpoint that the Sun resides inside the LIC, the locally observed velocity of the interstellar neutral helium is consistent with a linear combination of the velocities of the LIC and G cloud, but not with either of these two velocities. This finding shows that the Sun travels through a mixed-cloud interstellar medium composed of material from both these clouds. Interactions between these clouds explain the substantially higher density of the interstellar hydrogen near the Sun and toward stars located within the interaction region of these two clouds. The observed asymmetry of the interstellar helium distribution function also supports this interaction. The structure and equilibrium in this region require further studies using in situ and telescopic observations. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 13 pages, 5 figures, 1 table, accepted for publication in ApJL

Journal ref: ApJL, 937 (2022), L32

Turbulence in the outer heliosphere

Authors: Federico Fraternale, Laxman Adhikari, Horst Fichtner, Tae K. Kim, Jens Kleimann, Sean Oughton, Nikolai V. Pogorelov, Vadim Roytershteyn, Charles W. Smith, Arcadi V. Usmanov, G P. Zank, Lingling Zhao

Abstract: The solar wind (SW) and local interstellar medium (LISM) are turbulent media. Their interaction is governed by complex physical processes and creates heliospheric regions with significantly different properties in terms of particle populations, bulk flow and turbulence. Our knowledge of the solar wind turbulence \nature and dynamics mostly relies on near-Earth and near-Sun observations, and has be… ▽ More The solar wind (SW) and local interstellar medium (LISM) are turbulent media. Their interaction is governed by complex physical processes and creates heliospheric regions with significantly different properties in terms of particle populations, bulk flow and turbulence. Our knowledge of the solar wind turbulence \nature and dynamics mostly relies on near-Earth and near-Sun observations, and has been increasingly improving in recent years due to the availability of a wealth of space missions, including multi-spacecraft missions. In contrast, the properties of turbulence in the outer heliosphere are still not completely understood. In situ observations by Voyager and New Horizons, and remote neutral atom measurements by IBEX strongly suggest that turbulence is one of the critical processes acting at the heliospheric interface. It is intimately connected to charge exchange processes responsible for the production of suprathermal ions and energetic neutral atoms. This paper reviews the observational evidence of turbulence in the distant SW and in the LISM, advances in modeling efforts, and open challenges. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: 82 pages, 31 figures. Paper accepted by Space Science Reviews (collection: The Heliosphere in the Local Interstellar Medium: Into the Unknown), and presented at the ISSI Workshop on November 8-12, 2021

Modeling of Joint Parker Solar Probe - Metis/Solar Orbiter Observations

Authors: Laxman Adhikari, Gary P. Zank, Daniele Telloni, Lingling Zhao

Abstract: We present a first theoretical modeling of joint Parker Solar Probe (PSP) - Metis/Solar Orbiter (SolO) quadrature observations Telloni et al 2022c. The combined observations describe the evolution of a slow solar wind plasma parcel from the extended solar corona ($3.5-6.3$ R$_\odot$) to the very inner heliosphere (23.2 R$_\odot$). The Metis/SolO instrument remotely measures the solar wind speed fi… ▽ More We present a first theoretical modeling of joint Parker Solar Probe (PSP) - Metis/Solar Orbiter (SolO) quadrature observations Telloni et al 2022c. The combined observations describe the evolution of a slow solar wind plasma parcel from the extended solar corona ($3.5-6.3$ R$_\odot$) to the very inner heliosphere (23.2 R$_\odot$). The Metis/SolO instrument remotely measures the solar wind speed finding a range from $96-201$ kms$^{-1}$, and PSP measures the solar wind plasma in situ, observing a radial speed of 219.34 kms$^{-1}$. We find theoretically and observationally that the solar wind speed accelerates rapidly within 3.3 -- 4 R$_\odot$, and then increases more gradually with distance. Similarly, we find that the theoretical solar wind density is consistent with the remotely and in situ observed solar wind density. The normalized cross-helicity and normalized residual energy observed by PSP are 0.96 and -0.07, respectively, indicating that the slow solar wind is very Alfvénic. The theoretical NI/slab results are very similar to PSP measurements, which is a consequence of the highly magnetic field-aligned radial flow ensuring that PSP can measure slab fluctuations and not 2D. Finally, we calculate the theoretical 2D and slab turbulence pressure, finding that the theoretical slab pressure is very similar to that observed by PSP. △ Less

Submitted 27 June, 2022; originally announced June 2022.

Comments: 12 pages, 5 figures

Observation of Magnetic Switchback in the Solar Corona

Authors: Daniele Telloni, Gary P. Zank, Marco Stangalini, Cooper Downs, Haoming Liang, Masaru Nakanotani, Vincenzo Andretta, Ester Antonucci, Luca Sorriso-Valvo, Laxman Adhikari, Lingling Zhao, Raffaele Marino, Roberto Susino, Catia Grimani, Michele Fabi, Raffaella D'Amicis, Denise Perrone, Roberto Bruno, Francesco Carbone, Salvatore Mancuso, Marco Romoli, Vania Da Deppo, Silvano Fineschi, Petr Heinzel, John D. Moses , et al. (27 additional authors not shown)

Abstract: Switchbacks are sudden, large radial deflections of the solar wind magnetic field, widely revealed in interplanetary space by the Parker Solar Probe. The switchbacks' formation mechanism and sources are still unresolved, although candidate mechanisms include Alfvénic turbulence, shear-driven Kelvin-Helmholtz instabilities, interchange reconnection, and geometrical effects related to the Parker spi… ▽ More Switchbacks are sudden, large radial deflections of the solar wind magnetic field, widely revealed in interplanetary space by the Parker Solar Probe. The switchbacks' formation mechanism and sources are still unresolved, although candidate mechanisms include Alfvénic turbulence, shear-driven Kelvin-Helmholtz instabilities, interchange reconnection, and geometrical effects related to the Parker spiral. This Letter presents observations from the Metis coronagraph onboard Solar Orbiter of a single large propagating S-shaped vortex, interpreted as first evidence of a switchback in the solar corona. It originated above an active region with the related loop system bounded by open-field regions to the East and West. Observations, modeling, and theory provide strong arguments in favor of the interchange reconnection origin of switchbacks. Metis measurements suggest that the initiation of the switchback may also be an indicator of the origin of slow solar wind. △ Less

Submitted 9 June, 2022; v1 submitted 7 June, 2022; originally announced June 2022.

Turbulence in the Sub-Alfvénic Solar Wind

Authors: G. P. Zank, L. -L. Zhao, L. Adhikri, D. Telloni, J. C. Kasper, M. Stevens, A. Rahmati, S. D. Bale

Abstract: 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 P… ▽ More 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. △ Less

Submitted 5 February, 2022; originally announced February 2022.

Comments: Accepted in ApJL

Generalized Fluid Models of the Braginskii Type

Authors: P. Hunana, T. Passot, E. Khomenko, D. Martinez-Gomez, M. Collados, A. Tenerani, G. P. Zank, Y. Maneva, M. L. Goldstein, G. M. Webb

Abstract: 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, w… ▽ More 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. △ Less

Submitted 2 June, 2022; v1 submitted 27 January, 2022; originally announced January 2022.

Comments: Updated version. Title change from "Braginskii-type" to "Braginskii Type" to match the published version (accepted to ApJ Supplements, 29 Jan 2022). Bookmarks included, TOC is now at the beginning, numbering in Appendix I is fixed, few typos corrected

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

Authors: T. Getachew, D. J. McComas, C. J. Joyce, E. Palmerio, E. R. Christian, C. M. S. Cohen, M. I. Desai, J. Giacalone, M. E. Hill, W. H. Matthaeus, R. L. McNutt, D. G. Mitchell, J. G. Mitchell, J. S. Rankin, E. C. Roelof, N. A. Schwadron, J. R. Szalay, G. P. Zank, L. -L. Zhao, B. J. Lynch, T. D. Phan, S. D. Bale, P. L. Whittlesey, J. C. Kasper

Abstract: 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… ▽ More 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. △ Less

Submitted 8 December, 2021; originally announced December 2021.

Comments: 25 pages, 8 figures

Inertial range magnetic fluctuation anisotropy observed from PSP first seven orbits

Authors: L. -L. Zhao, G. P. Zank, L. Adhikari, M. Nakanotani

Abstract: 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… ▽ More 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. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: This paper is to be published in ApJL

On the Conservation of Turbulence Energy in Turbulence Transport Models

Authors: Bingbing Wang, Gary P. Zank, Laxman Adhikari, Lingling Zhao

Abstract: Zank et al. developed models describing the transport of low frequency incompressible and nearly incompressible turbulence in inhomogeneous flows. The formalism was based on expressing the fluctuating variables in terms of the Elsässar variables and then taking "moments" subject to various closure hypotheses. The turbulence transport models are different according to whether the plasma beta regime… ▽ More Zank et al. developed models describing the transport of low frequency incompressible and nearly incompressible turbulence in inhomogeneous flows. The formalism was based on expressing the fluctuating variables in terms of the Elsässar variables and then taking "moments" subject to various closure hypotheses. The turbulence transport models are different according to whether the plasma beta regime is large or of order 1 or smaller. Here, we show explicitly that the two sets of turbulence transport models admit a conservation representation that resembles the well-known WKB transport equation for Alfvén wave energy density after introducing appropriate definitions of the "pressure" associated with the turbulent fluctuations. This includes introducing a distinct turbulent pressure tensor for 3D incompressible turbulence (the large plasma beta limit) and pressure tensors for quasi-2D and slab turbulence (the plasma beta order 1 or small regimes) that generalize the form of the WKB pressure tensor. Various limits of the different turbulent pressure tensors are discussed. However, the analogy between the conservation form of the turbulence transport models and the WKB model is not close for multiple reasons, including that the turbulence models express fully nonlinear physical processes unlike the strictly linear WKB description. The analysis presented here serves both as a check on the validity and correctness of the turbulence transport models and provides greater transparency of the energy dissipation term and the "turbulent pressure" in our models, which is important for many practical applications. △ Less

Submitted 22 November, 2021; originally announced November 2021.

Comments: 19 pages, 0 figure

Signature of a heliotail organized by the solar magnetic field and the role of non-ideal processes in modeled IBEX ENA maps: a comparison of the BU and Moscow MHD models

Authors: M. Kornbleuth, M. Opher, I. Baliukin, M. A. Dayeh, E. Zirnstein, M. Gkioulidou, K. Dialynas, A. Galli, J. D. Richardson, V. Izmodenov, G. P. Zank, S. Fuselier

Abstract: Energetic neutral atom (ENA) models typically require post-processing routines to convert the distributions of plasma and H atoms into ENA maps. Here we investigate how two different kinetic-MHD models of the heliosphere (the BU and Moscow models) manifest in modeled ENA maps using the same prescription and how they compare with Interstellar Boundary Explorer (IBEX) observations. Both MHD models t… ▽ More Energetic neutral atom (ENA) models typically require post-processing routines to convert the distributions of plasma and H atoms into ENA maps. Here we investigate how two different kinetic-MHD models of the heliosphere (the BU and Moscow models) manifest in modeled ENA maps using the same prescription and how they compare with Interstellar Boundary Explorer (IBEX) observations. Both MHD models treat the solar wind as a single-ion plasma for protons, which include thermal solar wind ions, pick-up ions (PUIs), and electrons. Our ENA prescription partitions the plasma into three distinct ion populations (thermal solar wind, PUIs transmitted and ones energized at the termination shock) and models the populations with Maxwellian distributions. Both kinetic-MHD heliospheric models produce a heliotail with heliosheath plasma organized by the solar magnetic field into two distinct north and south columns that become lobes of high mass flux flowing down the heliotail, though in the BU model the ISM flows between the two lobes at distances in the heliotail larger than 300 AU. While our prescription produces similar ENA maps for the two different plasma and H atom solutions at the IBEX-Hi energy range (0.5 - 6 keV), the modeled ENA maps require a scaling factor of ~2 to be in agreement with the data. This problem is present in other ENA models with the Maxwellian approximation of multiple ion species and indicates that a higher neutral density or some acceleration of PUIs in the heliosheath is required. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 29 pages, 3 tables, 9 figures, accepted to ApJ

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

Authors: M. Kornbleuth, M. Opher, I. Baliukin, M. Gkioulidou, J. D. Richardson, G. P. Zank, A. T. Michael, G. Toth, V. Tenishev, V. Izmodenov, D. Alexashov, S. Fuselier, J. F. Drake, K. Dialynas

Abstract: 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 B… ▽ More 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. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 2 tables, 12 figures, accepted to ApJ

Exploring the Solar Wind from its Source on the Corona into the Inner Heliosphere during the First Solar Orbiter - Parker Solar Probe Quadrature

Authors: Daniele Telloni, Vincenzo Andretta, Ester Antonucci, Alessandro Bemporad, Giuseppe E. Capuano, Silvano Fineschi, Silvio Giordano, Shadia Habbal, Denise Perrone, Rui F. Pinto, Luca Sorriso-Valvo, Daniele Spadaro, Roberto Susino, Lloyd D. Woodham, Gary P. Zank, Marco Romoli, Stuart D. Bale, Justin C. Kasper, Frédéric Auchère, Roberto Bruno, Gerardo Capobianco, Anthony W. Case, Chiara Casini, Marta Casti, Paolo Chioetto , et al. (46 additional authors not shown)

Abstract: This Letter addresses the first Solar Orbiter (SO) -- Parker Solar Probe (PSP) quadrature, occurring on January 18, 2021, to investigate the evolution of solar wind from the extended corona to the inner heliosphere. Assuming ballistic propagation, the same plasma volume observed remotely in corona at altitudes between 3.5 and 6.3 solar radii above the solar limb with the Metis coronagraph on SO ca… ▽ More This Letter addresses the first Solar Orbiter (SO) -- Parker Solar Probe (PSP) quadrature, occurring on January 18, 2021, to investigate the evolution of solar wind from the extended corona to the inner heliosphere. Assuming ballistic propagation, the same plasma volume observed remotely in corona at altitudes between 3.5 and 6.3 solar radii above the solar limb with the Metis coronagraph on SO can be tracked to PSP, orbiting at 0.1 au, thus allowing the local properties of the solar wind to be linked to the coronal source region from where it originated. Thanks to the close approach of PSP to the Sun and the simultaneous Metis observation of the solar corona, the flow-aligned magnetic field and the bulk kinetic energy flux density can be empirically inferred along the coronal current sheet with an unprecedented accuracy, allowing in particular estimation of the Alfvén radius at 8.7 solar radii during the time of this event. This is thus the very first study of the same solar wind plasma as it expands from the sub-Alfvénic solar corona to just above the Alfvén surface. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 10 pages, 4 figures

Journal ref: Telloni, D., Andretta, V., Antonucci, E., et al. 2021, ApJL, 920, L14

Turbulence/wave transmission at an ICME-driven shock observed by Solar Orbiter and Wind

Authors: L. L. Zhao, G. P. Zank, J. S. He, D. Telloni, Q. Hu, G. Li, M. Nakanotani, L. Adhikari, E. K. J. Kilpua, T. S. Horbury, H. O'Brien, V. Evans, V. Angelini

Abstract: Solar Orbiter observed an interplanetary coronal mass ejection (ICME) event at 0.8 AU on 2020 April 19. The ICME was also observed by Wind at 1 AU on 2020 April 20. An interplanetary shock wave was driven in front of the ICME. We focus on the transmission of the magnetic fluctuations across the shock and analyze the characteristic wave modes of solar wind turbulence near the shock observed by both… ▽ More Solar Orbiter observed an interplanetary coronal mass ejection (ICME) event at 0.8 AU on 2020 April 19. The ICME was also observed by Wind at 1 AU on 2020 April 20. An interplanetary shock wave was driven in front of the ICME. We focus on the transmission of the magnetic fluctuations across the shock and analyze the characteristic wave modes of solar wind turbulence near the shock observed by both spacecraft. The ICME event is characterized by a magnetic helicity based technique. The shock normal is determined by magnetic coplanarity method for Solar Orbiter and using a mixed coplanarity approach for Wind. The power spectra of magnetic field fluctuations are generated by applying both a fast Fourier transform and Morlet wavelet analysis. To understand the nature of waves observed near the shock, we use the normalized magnetic helicity as a diagnostic parameter. The wavelet reconstructed magnetic field fluctuation hodograms are used to further study the polarization properties of waves. We find that the ICME-driven shock observed by Solar Orbiter and Wind is a fast forward oblique shock with a more perpendicular shock angle at 1 AU. After the shock crossing, the magnetic field fluctuation power increases. Most of the magnetic field fluctuation power resides in the transverse fluctuations. In the vicinity of the shock, both spacecraft observe right-hand polarized waves in the spacecraft frame. The upstream wave signatures fall in a relatively broad and low-frequency band, which might be attributed to low-frequency MHD waves excited by the streaming particles. For the downstream magnetic wave activity, we find oblique kinetic Alfven waves with frequencies near the proton cyclotron frequency in the spacecraft frame. The frequency of the downstream waves increases by a factor of 7-10 due to the shock compression and the Doppler effect. △ Less

Submitted 5 February, 2021; originally announced February 2021.

Journal ref: A&A 656, A3 (2021)

Detection of small magnetic flux ropes from the third and fourth Parker Solar Probe encounters

Authors: L. -L. Zhao, G. P. Zank, Q. Hu, D. Telloni, Y. Chen, L. Adhikari, M. Nakanotani, J. C. Kasper, J. Huang, S. D. Bale, K. E. Korreck, A. W. Case, M. Stevens, J. W. Bonnell, T. Dudok de Wit, K. Goetz, P. R. Harvey, R. J. MacDowall, D. M. Malaspina, M. Pulupa, D. E. Larson, R. Livi, P. Whittlesey, K. G. Klein, N. E. Raouafi

Abstract: We systematically search for magnetic flux rope structures in the solar wind to within the closest distance to the Sun of 0.13 AU, using data from the third and fourth orbits of the Parker Solar Probe. We extend our previous magnetic helicity based technique of identifying magnetic flux rope structures. The method is improved upon to incorporate the azimuthal flow, which becomes larger as the spac… ▽ More We systematically search for magnetic flux rope structures in the solar wind to within the closest distance to the Sun of 0.13 AU, using data from the third and fourth orbits of the Parker Solar Probe. We extend our previous magnetic helicity based technique of identifying magnetic flux rope structures. The method is improved upon to incorporate the azimuthal flow, which becomes larger as the spacecraft approaches the Sun. A total of 21 and 34 magnetic flux ropes are identified during the third (21 days period) and fourth (17 days period) orbits of the Parker Solar Probe, respectively. We provide a statistical analysis of the identified structures, including their relation to the streamer belt and heliospheric current sheet crossing. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Journal ref: A&A 650, A12 (2021)

A data-driven MHD model of the weakly-ionized chromosphere

Authors: Mehmet Sarp Yalim, Avijeet Prasad, Nikolai Pogorelov, Gary Zank, Qiang Hu

Abstract: The physics of the solar chromosphere is complex from both theoretical and modeling perspectives. The plasma temperature from the photosphere to corona increases from ~5,000 K to ~1 million K over a distance of only ~10,000 km from the chromosphere and the transition region. Certain regions of the solar atmosphere have sufficiently low temperature and ionization rates to be considered as weakly-io… ▽ More The physics of the solar chromosphere is complex from both theoretical and modeling perspectives. The plasma temperature from the photosphere to corona increases from ~5,000 K to ~1 million K over a distance of only ~10,000 km from the chromosphere and the transition region. Certain regions of the solar atmosphere have sufficiently low temperature and ionization rates to be considered as weakly-ionized. In particular, this is true at the lower chromosphere. In this paper, we present an overview of our data-driven magnetohydrodynamics model for the weakly-ionized chromosphere and show a benchmark result. It utilizes the Cowling resistivity which is orders of magnitude greater than the Coulomb resistivity. Ohm's law therefore includes anisotropic dissipation. We investigate the effects of the Cowling resistivity on heating and magnetic reconnection in the chromosphere as the flare-producing active region (AR) 11166 evolves. In particular, we analyze a C2.0 flare emerging from AR11166 and find a normalized reconnection rate of 0.12. △ Less

Submitted 24 July, 2020; originally announced July 2020.

Comments: 14 pages, 4 figures, Journal of Physics Conference Series - Proceedings of the 19th Annual International Astrophysics Conference, Santa Fe, NM, March 9-13, 2020 (accepted)

Effects of Cowling Resistivity in the Weakly-Ionized Chromosphere

Authors: Mehmet Sarp Yalim, Avijeet Prasad, Nikolai Pogorelov, Gary Zank, Qiang Hu

Abstract: The physics of the solar chromosphere is complex from both theoretical and modeling perspectives. The plasma temperature from the photosphere to corona increases from ~5,000 K to ~1 million K over a distance of only ~10,000 km from the chromosphere and the transition region. Certain regions of the solar atmosphere have sufficiently low temperature and ionization rates to be considered as weakly-io… ▽ More The physics of the solar chromosphere is complex from both theoretical and modeling perspectives. The plasma temperature from the photosphere to corona increases from ~5,000 K to ~1 million K over a distance of only ~10,000 km from the chromosphere and the transition region. Certain regions of the solar atmosphere have sufficiently low temperature and ionization rates to be considered as weakly-ionized. In particular, this is true at the lower chromosphere. As a result, the Cowling resistivity is orders of magnitude greater than the Coulomb resistivity. Ohm's law therefore includes anisotropic dissipation. To evaluate the Cowling resistivity, we need to know the external magnetic field strength and to estimate the neutral fraction as a function of the bulk plasma density and temperature. In this study, we determine the magnetic field topology using the non-force-free field (NFFF) extrapolation technique based on SDO/HMI SHARP vector magnetogram data, and the stratified density and temperature profiles from the Maltby-M umbral core model for sunspots. We investigate the variation and effects of Cowling resistivity on heating and magnetic reconnection in the chromosphere as the flare-producing active region (AR) 11166 evolves. In particular, we analyze a C2.0 flare emerging from AR11166 and find a normalized reconnection rate of 0.051. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: 15 pages, 5 figures, The Astrophysical Journal Letters (in press)

An ACE/CRIS-observation-based Galactic Cosmic Rays heavy nuclei spectra model II

Authors: Shuai Fu, LingLing Zhao, Gary P. Zank, Miao Wang, Yong Jiang

Abstract: An observation-based Galactic Cosmic Ray (GCR) spectral model for heavy nuclei is developed. Zhao and Qin (J. Geophys. Res. Space Phys.118, 1837(2013)) proposed an empirical elemental GCR spectra model for nuclear charge 5-28 over the energy range from 30 to 500 MeV/nuc, which is proved to be successful in predicting yearly averaged GCR heavy nuclei spectra.Based on the latest highly statistically… ▽ More An observation-based Galactic Cosmic Ray (GCR) spectral model for heavy nuclei is developed. Zhao and Qin (J. Geophys. Res. Space Phys.118, 1837(2013)) proposed an empirical elemental GCR spectra model for nuclear charge 5-28 over the energy range from 30 to 500 MeV/nuc, which is proved to be successful in predicting yearly averaged GCR heavy nuclei spectra.Based on the latest highly statistically precise measurements from ACE/CRIS,a further elemental GCR model with monthly averaged spectra is presented. The model can reproduce the past and predict the futureGCR intensity monthly by correlating model parameters with thecontinuous sunspot number (SSN) record. The effects of solar activity on GCR modulation are considered separately for odd and even solar cycles. Compared with other comprehensive GCR models, our modeling results are satisfyingly consistent with the GCR spectral measurements from ACE/SIS and IMP-8, and have comparable prediction accuracy as the Badhwar & O'Neill 2014 model.A detailed error analysis is also provided.Finally, the GCR carbon and iron nuclei fluxes for the subsequent two solar cycles (SC 25 and 26) are predicted and they show a potential trend in reduced flux amplitude, which is suspected to be relevant to possible weak solar cycles. △ Less

Submitted 1 May, 2020; originally announced May 2020.

Journal ref: Sci. China-Phys. Mech. Astron. 63, 219511 (2020)

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

Authors: S. Matsukiyo, T. Noumi, G. P. Zank, H. Washimi, T. Hada

Abstract: 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… ▽ More 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. △ Less

Submitted 18 December, 2019; originally announced January 2020.

Comments: 14 figures, accepted for publication in ApJ

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

Authors: T. K. Kim, N. V. Pogorelov, C. N. Arge, C. J. Henney, S. I. Jones-Mecholsky, W. P. Smith, S. D. Bale, J. W. Bonnell, T. Dudok de Wit, K. Goetz, P. R. Harvey, R. J. MacDowall, D. M. Malaspina, M. Pulupa, J. C. Kasper, K. E. Korreck, M. Stevens, A. W. Case, P. Whittlesey, R. Livi, D. E. Larson, K. G. Klein, G. P. Zank

Abstract: 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 magnet… ▽ More 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. △ Less

Submitted 5 December, 2019; originally announced December 2019.

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

Authors: L. -L. Zhao, G. P. Zank, L. Adhikari, Q. Hu, J. C. Kasper, S. D. Bale, K. E. Korreck, A. W. Case, M. Stevens, J. W. Bonnell, T. Dudok de Wit, K. Goetz, P. R. Harvey, R. J. MacDowall, D. M. Malaspina, M. Pulupa, D. E. Larson, R. Livi, P. Whittlesey, K. G. Klein

Abstract: 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… ▽ More 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. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 22 pages, 5 figures, accepted in ApJS

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

Authors: A. Czechowski, M. Bzowski, J. M. Sokół, M. A. Kubiak, J. Heerikhuisen, E. J. Zirnstein, N. V. Pogorelov, N. A. Schwadron, M. Hilchenbach, J. Grygorczuk, G. P. Zank

Abstract: 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 h… ▽ More 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. △ Less

Submitted 22 November, 2019; originally announced November 2019.

Comments: Accepted for publication in The Astrophysical Journal

Energy Dissipation and Entropy in Collisionless Plasma

Authors: Senbei Du, Gary P. Zank, Fan Guo, Xiaocan Li

Abstract: It is well known that collisionless systems are dissipation free from the perspective of particle collision and thus conserve entropy. On the other hand, processes such as magnetic reconnection and turbulence appear to convert large-scale magnetic energy into heat. In this paper, we investigate the energization and heating of collisionless plasma. The dissipation process is discussed in terms of f… ▽ More It is well known that collisionless systems are dissipation free from the perspective of particle collision and thus conserve entropy. On the other hand, processes such as magnetic reconnection and turbulence appear to convert large-scale magnetic energy into heat. In this paper, we investigate the energization and heating of collisionless plasma. The dissipation process is discussed in terms of fluid entropy in both isotropic and gyrotropic forms. Evolution equations for the entropy are derived and they reveal mechanisms that lead to changes in fluid entropy. These equations are verified by a collisionless particle-in-cell simulation of multiple reconnecting current sheets. In addition to previous findings regarding the pressure tensor, we emphasize the role of heat flux in the dissipation process. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: 19 pages, 4 figure

Journal ref: Phys. Rev. E 101, 033208 (2020)

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

Authors: H. Che, G. P. Zank

Abstract: 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… ▽ More 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. △ Less

Submitted 24 August, 2019; originally announced August 2019.

Comments: in publication

Effect of Star Rotation Rates on the Characteristics of Energetic Particle Events

Authors: Shuai Fu, Yong Jiang, Vladimir Airapetian, Junxiang Hu, Gang Li, Gary Zank

Abstract: Recent detection of superflares on solar-type stars by Kepler mission raised a possibility that they can be associated with energetic coronal mass ejections (CMEs) and energetic particle events (SEPs). These space weather events can impact habitability of exoplanets around these stars. Here we use the improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, to model the time… ▽ More Recent detection of superflares on solar-type stars by Kepler mission raised a possibility that they can be associated with energetic coronal mass ejections (CMEs) and energetic particle events (SEPs). These space weather events can impact habitability of exoplanets around these stars. Here we use the improved Particle Acceleration and Transport in the Heliosphere (iPATH) model, to model the time intensity profile and spectrum of SEPs accelerated at CME-driven shocks from stars of different ages traced by their rotation rates. We consider a solar-like (G-type) star with 6 different rotation rates varying from 0.5 to 3.0 times rotation speed of the sun. In all 6 cases, a fast CME is launched with the same speed of 1500 km/sec and the resulting time intensity profiles at 3 locations and and energy spectra at 5 locations at 1 AU are obtained. The maximum particle energy at the shock front as a function of r is also shown. Our results suggest that within 0.8 AU the maximum particle energy at the shock front increases with the rotation rate of the star. However, event integrated spectra for the five selected locations along the CME path show complicated patterns. This is because the Parker magnetic field for rapidly rotating stars is more tightly winded. Our results can be used in estimating the radiation environments of terrestrial-type exoplanets around solar-type stars. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: Manuscript has been accepted by The Astrophysical Journal Letters

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

Authors: V. S. Airapetian, R. Barnes, O. Cohen, G. A. Collinson, W. C. Danchi, C. F. Dong, A. D. Del Genio, K. France, K. Garcia-Sage, A. Glocer, N. Gopalswamy, J. L. Grenfell, G. Gronoff, M. G"udel, K. Herbst, W. G. Henning, C. H. Jackman, M. Jin, C. P. Johnstone, L. Kaltenegger, C. D. Kay, K. Kobayashi, W. Kuang, G. Li, B. J. Lynch , et al. (21 additional authors not shown)

Abstract: 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 astro… ▽ More 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. △ Less

Submitted 19 May, 2019; v1 submitted 9 May, 2019; originally announced May 2019.

Comments: 206 pages, 24 figures, 1 table; Review paper. International Journal of Astrobiology (2019)

Reconstructing Extreme Space Weather from Planet Hosting Stars

Authors: V. S. Airapetian, V. Adibekyan, M. Ansdell, D. Alexander, T. Bastian, S. Boro Saikia, A. S. Brun, O. Cohen, M. Cuntz, W. Danchi, J. Davenport, J. DeNolfo, R. DeVore, C. F. Dong, J. J. Drake, K. France, F. Fraschetti, K. Herbst, K. Garcia-Sage, M. Gillon, A. Glocer, J. L. Grenfell, G. Gronoff, N. Gopalswamy, M. Guedel , et al. (58 additional authors not shown)

Abstract: The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditi… ▽ More The field of exoplanetary science is making rapid progress both in statistical studies of exoplanet properties as well as in individual characterization. As space missions provide an emerging picture of formation and evolution of exoplanetary systems, the search for habitable worlds becomes one of the fundamental issues to address. To tackle such a complex challenge, we need to specify the conditions favorable for the origin, development and sustainment of life as we know it. This requires the understanding of global (astrospheric) and local (atmospheric, surface and internal) environments of exoplanets in the framework of the physical processes of the interaction between evolving planet-hosting stars along with exoplanetary evolution over geological timescales, and the resulting impact on climate and habitability of exoplanets. Feedbacks between astrophysical, physico-chemical atmospheric and geological processes can only be understood through interdisciplinary studies with the incorporation of progress in heliophysics, astrophysics, planetary, Earth sciences, astrobiology, and the origin of life communities. The assessment of the impacts of host stars on the climate and habitability of terrestrial (exo)planets and potential exomoons around them may significantly modify the extent and the location of the habitable zone and provide new directions for searching for signatures of life. Thus, characterization of stellar ionizing outputs becomes an important task for further understanding the extent of habitability in the universe. The goal of this white paper is to identify and describe promising key research goals to aid the theoretical characterization and observational detection of ionizing radiation from quiescent and flaring upper atmospheres of planet hosts as well as properties of stellar coronal mass ejections and stellar energetic particle events. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: White Paper submitted to the Astronomy and Astrophysics Decadal Survey (Astro2020), 8 pages, 1 figure

Plasma Energization in Colliding Magnetic Flux Ropes

Authors: Senbei Du, Fan Guo, Gary P. Zank, Xiaocan Li, Adam Stanier

Abstract: Magnetic flux ropes are commonly observed throughout the heliosphere, and recent studies suggest that interacting flux ropes are associated with some energetic particle events. In this work, we carry out 2D particle-in-cell (PIC) simulations to study the coalescence of two magnetic flux ropes (or magnetic islands), and the subsequent plasma energization processes. The simulations are initialized w… ▽ More Magnetic flux ropes are commonly observed throughout the heliosphere, and recent studies suggest that interacting flux ropes are associated with some energetic particle events. In this work, we carry out 2D particle-in-cell (PIC) simulations to study the coalescence of two magnetic flux ropes (or magnetic islands), and the subsequent plasma energization processes. The simulations are initialized with two magnetic islands embedded in a reconnecting current sheet. The two islands collide and eventually merge into a single island. Particles are accelerated during this process as the magnetic energy is released and converted to the plasma energy, including bulk kinetic energy increase by the ideal electric field, and thermal energy increase by the fluid compression and the non-ideal electric field. We find that contributions from these different energization mechanisms are all important and comparable with each other. Fluid shear and a non-gyrotropic pressure tensor also contribute to the energy conversion process. For simulations with different box sizes ranging from $L_x \sim 25$--$100 d_i$ and ion-to-electron mass ratios $m_i / m_e = 25$, 100 and 400, we find that the general evolution is qualitatively the same for all runs, and the energization depends only weakly on either the system size or the mass ratio. The results may help us understand plasma energization in solar and heliospheric environments. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Exploring Extreme Space Weather Factors of Exoplanetary Habitability

Authors: V. S. Airapetian, V. Adibekyan, M. Ansdell, O. Cohen, M. Cuntz, W. Danchi, C. F. Dong, J. J. Drake, A. Fahrenbach, K. France, K. Garcia-Sage, A. Glocer, J. L. Grenfell, G. Gronoff, H. Hartnett, W. Henning, N. R. Hinkel, A. G. Jensen, M. Jin, P. Kalas, S. R. Kane, K. Kobayashi, R. Kopparapu, J. Leake, M. López-Puertas , et al. (24 additional authors not shown)

Abstract: It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems.… ▽ More It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years. △ Less

Submitted 9 March, 2018; originally announced March 2018.

Comments: 6 pages, the white paper submitted to the US National Academy of Sciences call on Exoplanet Science Strategy

Life Beyond the Solar System: Space Weather and Its Impact on Habitable Worlds

Authors: V. S. Airapetian, W. C. Danchi, C. F. Dong, S. Rugheimer, M. Mlynczak, K. B. Stevenson, W. G. Henning, J. L. Grenfell, M. Jin, A. Glocer, G. Gronoff, B. Lynch, C. Johnstone, T. Lueftinger, M. Guedel, K. Kobayashi, A. Fahrenbach, G. Hallinan, V. Stamenkovic, O. Cohen, W. Kuang, B. van der Holst, C. Manchester, G. Zank, O. Verkhoglyadova , et al. (11 additional authors not shown)

Abstract: The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the sear… ▽ More The search of life in the Universe is a fundamental problem of astrobiology and a major priority for NASA. A key area of major progress since the NASA Astrobiology Strategy 2015 (NAS15) has been a shift from the exoplanet discovery phase to a phase of characterization and modeling of the physics and chemistry of exoplanetary atmospheres, and the development of observational strategies for the search for life in the Universe by combining expertise from four NASA science disciplines including heliophysics, astrophysics, planetary science and Earth science. The NASA Nexus for Exoplanetary System Science (NExSS) has provided an efficient environment for such interdisciplinary studies. Solar flares, coronal mass ejections and solar energetic particles produce disturbances in interplanetary space collectively referred to as space weather, which interacts with the Earth upper atmosphere and causes dramatic impact on space and ground-based technological systems. Exoplanets within close in habitable zones around M dwarfs and other active stars are exposed to extreme ionizing radiation fluxes, thus making exoplanetary space weather (ESW) effects a crucial factor of habitability. In this paper, we describe the recent developments and provide recommendations in this interdisciplinary effort with the focus on the impacts of ESW on habitability, and the prospects for future progress in searching for signs of life in the Universe as the outcome of the NExSS workshop held in Nov 29 - Dec 2, 2016, New Orleans, LA. This is one of five Life Beyond the Solar System white papers submitted by NExSS to the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe. △ Less

Submitted 16 January, 2018; originally announced January 2018.

Comments: 5 pages, the white paper was submitted to the National Academy of Sciences in support of the Astrobiology Science Strategy for the Search for Life in the Universe

Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields

Authors: N. V. Pogorelov, H. Fichtner, A. Czechowski, A. Lazarian, B. Lembege, J. A. le Roux, M. S. Potgieter, K. Scherer, E. C. Stone, R. D. Strauss, T. Wiengarten, P. Wurz, G. P. Zank, M. Zhang

Abstract: This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. T… ▽ More This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. The importance of magnetic field, charge exchange between atoms and ions, and solar cycle on the heliopause topology and observed heliocentric distances to different heliospheric discontinuities are discussed. It is shown that time-dependent boundary conditions are necessary to describe the heliospheric asymmetries detected by the Voyager spacecraft. We also discuss the structure of the HP, especially due to its instability and magnetic reconnection. It is demonstrated that the Rayleigh-Taylor instability of the nose of the HP creates consecutive layers of the interstellar and heliospheric plasma which are magnetically connected to different sources. This may be a possible explanation of abrupt changes in the galactic and anomalous cosmic ray fluxes observed by Voyager 1 when it was crossing the HP structure for a period of about one month in the summer of 2012. This paper also discusses the plausibility of fitting simulation results to a number of observational data sets obtained by in situ and remote measurements. The distribution of magnetic field in the vicinity of the HP is discussed in the context of Voyager measurements. We discuss the transport of energetic particles in the inner and outer heliosheath, concentrating on the anisotropic spatial diffusion diffusion tensor and the pitch-angle dependence of perpendicular diffusion and demonstrate that the latter can explain the observed pitch-angle anisotropies of both the anomalous and galactic cosmic rays in the outer heliosheath. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: Publication resulting from an International ISSI Team of the same name. Submitted to Space Science Reviews

Pickup Ion Effect of the Solar Wind Interaction with the Local Interstellar Medium

Authors: Nikolai V. Pogorelov, Matthew C. Bedford, Igor A. Kryukov, Gary P. Zank

Abstract: Pickup ions are created when interstellar neutral atoms resonantly exchange charge with the solar wind (SW) ions, especially in the supersonic part of the wind, where they carry most of the plasma pressure. Here we present numerical simulation results of the 3D heliospheric interface treating pickup ions as a separate proton fluid. To satisfy the fundamental conservation laws, we solve the system… ▽ More Pickup ions are created when interstellar neutral atoms resonantly exchange charge with the solar wind (SW) ions, especially in the supersonic part of the wind, where they carry most of the plasma pressure. Here we present numerical simulation results of the 3D heliospheric interface treating pickup ions as a separate proton fluid. To satisfy the fundamental conservation laws, we solve the system of equations describing the flow of the mixture of electrons, thermal protons, and pickup ions. To find the density and pressure of pickup ions behind the termination shock, we employ simple boundary conditions that take into account the \emph{Voyager} observations that showed that the decrease in the kinetic energy of the mixture at the termination shock predominantly contributed to the increase in the pressure of pickup ions. We show that this model adequately describes the flow of the plasma mixture and results in a noticeable decrease in the heliosheath width. △ Less

Submitted 10 September, 2016; originally announced September 2016.

Comments: Accepted for publication in Journal of Physics, Conference Series

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

Authors: V. I. Abramenko, G. P. Zank, A. Dosch, V. B. Yurchyshyn, P. R. Goode, K. Ahn, W. Cao

Abstract: 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 dir… ▽ More 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. △ Less

Submitted 16 July, 2013; originally announced July 2013.

Comments: 21 pages, 9 figures

On the spectral hardening at ~> 300 keV in solar flares

Authors: Gang. Li, XiangLiang Kong, Gary Zank, Yao Chen

Abstract: It has been noted for a long time that the spectra of observed continuum emissions in many solar flares are consistent with double power laws with a hardening at energies $\sim > $ 300 keV. It is now largely believed that at least in electron-dominated events the hardening in photon spectrum reflects an intrinsic hardening in the source electron spectrum. In this paper, we point out that a power l… ▽ More It has been noted for a long time that the spectra of observed continuum emissions in many solar flares are consistent with double power laws with a hardening at energies $\sim > $ 300 keV. It is now largely believed that at least in electron-dominated events the hardening in photon spectrum reflects an intrinsic hardening in the source electron spectrum. In this paper, we point out that a power law spectrum of electron with a hardening at high energies can be explained by diffusive shock acceleration of electrons at a termination shock with a finite width. Our suggestion is based on an early analytical work by Drury et al., where the steady state transport equation at a shock with a tanh profile was solved for a $p$-independent diffusion coefficient. Numerical simulations with a $p$-dependent diffusion coefficient show hardenings in the accelerated electron spectrum which are comparable with observations. One necessary condition for our proposed scenario to work is that high energy electrons resonate with the inertial range of the MHD turbulence and low energy electrons resonate with the dissipation range of the MHD turbulence at the acceleration site, and the spectrum of the dissipation range $\sim k^{-2.7}$. A $\sim k^{-2.7}$ dissipation range spectrum is consistent with recent solar wind observations. △ Less

Submitted 24 March, 2013; originally announced March 2013.

Comments: 20 pages, 2 figures, accepted by ApJ

Trajectories and Distribution of Interstellar Dust Grains in the Heliosphere

Authors: Jonathan D. Slavin, Priscilla C. Frisch, Hans-Reinhard Mueller, Jacob Heerikhuisen, Nikolai V. Pogorelov, William T. Reach, Gary Zank

Abstract: The solar wind carves a bubble in the surrounding interstellar medium (ISM), known as the heliosphere. Charged interstellar dust grains (ISDG) encountering the heliosphere may be diverted around the heliopause or penetrate it depending on their charge-to-mass ratio. We present new calculations of trajectories of ISDG in the heliosphere, and the dust density distributions that result. We include up… ▽ More The solar wind carves a bubble in the surrounding interstellar medium (ISM), known as the heliosphere. Charged interstellar dust grains (ISDG) encountering the heliosphere may be diverted around the heliopause or penetrate it depending on their charge-to-mass ratio. We present new calculations of trajectories of ISDG in the heliosphere, and the dust density distributions that result. We include up-to-date grain charging calculations using a realistic UV radiation field and full 3-D magnetohydrodynamic fluid + kinetic models for the heliosphere. Models with two different (constant) polarities for the solar wind magnetic field (SWMF) are used, with the grain trajectory calculations done separately for each polarity. Small grains a_gr ~ 0.01 micron are completely excluded from the inner heliosphere. Large grains, a_gr ~ 1.0 micron pass into the inner solar system and are concentrated near the Sun by its gravity. Trajectories of intermediate size grains depend strongly on the SWMF polarity. When the field has magnetic north pointing to ecliptic north, the field de-focuses the grains resulting in low densities in the inner heliosphere, while for the opposite polarity the dust is focused near the Sun. The ISDG density outside the heliosphere inferred from applying the model results to in situ dust measurements is inconsistent with local ISM depletion data for both SWMF polarities, but is bracketed by them. This result points to the need to include the time variation in the SWMF polarity during grain propagation. Our results provide valuable insights for interpretation of the in situ dust observations from Ulysses. △ Less

Submitted 3 October, 2012; originally announced October 2012.

Comments: 34 pages, 15 figures, Accepted for publication in the Astrophysical Journal

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

Authors: P. C. Frisch, J. Heerikhuisen, N. V. Pogorelov, B. DeMajistre, G. B. Crew, H. O. Funsten, P. Janzen, D. J. McComas, E. Moebius, H. -R. Mueller, D. B. Reisenfeld, N. A. Schwadron, J. D. Slavin, G. P. Zank

Abstract: 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,… ▽ More 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. △ Less

Submitted 22 November, 2010; originally announced November 2010.

Journal ref: Published in ApJ 2010 volume 719 pages 1984-1992

Exclusion of Tiny Interstellar Dust Grains from the Heliosphere

Authors: J. D. Slavin, P. C. Frisch, J. Heerikhuisen, N. V. Pogorelov, H. -R. Mueller, W. T. Reach, G. P. Zank, B. Dasgupta, K. Avinash

Abstract: 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 s… ▽ More 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. △ Less

Submitted 9 November, 2009; originally announced November 2009.

Comments: 5 pages, 5 figures, accepted for publication in the Solar Wind 12 conference proceedings

Journal ref: AIP Conf.Proc.1216:497-501,2010

The Heliosphere in Time

Authors: H. -R. Müller, P. C. Frisch, B. D. Fields, G. P. Zank

Abstract: 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 s… ▽ More 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. △ Less

Submitted 2 October, 2008; originally announced October 2008.

Comments: 12 pages; to appear in Space Science Reviews, "From the Outer Heliosphere to the Local Bubble", J. Linsky et al. (eds); ISSI workshop held Oct. 2007

Journal ref: Space Science Reviews 143, 415-425 (2009)

The effects of a kappa-distribution in the heliosheath on the global heliosphere and ENA flux at 1 AU

Authors: J. Heerikhuisen, N. V. Pogorelov, V. Florinski, G. P. Zank, J. A. le Roux

Abstract: We investigate heliosheath energetic neutral atom (ENA) fluxes at keV energies, by assuming that the heliosheath proton distribution can be approximated by a kappa-distribution. The choice of the kappa parameter derives from observational data of the solar wind (SW). This has direct applications to the upcoming IBEX mission. We will look at all-sky ENA maps within the IBEX energy range (10 eV to… ▽ More We investigate heliosheath energetic neutral atom (ENA) fluxes at keV energies, by assuming that the heliosheath proton distribution can be approximated by a kappa-distribution. The choice of the kappa parameter derives from observational data of the solar wind (SW). This has direct applications to the upcoming IBEX mission. We will look at all-sky ENA maps within the IBEX energy range (10 eV to 6 keV), as well as ENA energy spectra in several directions. We find that the use of kappa, as opposed to a Maxwellian, gives rise to greatly increased ENA fluxes above 1 keV, while medium energy fluxes are somewhat reduced. We show how IBEX data can be used to estimate the spectral slope in the heliosheath, and that the use of kappa reduces the differences between ENA maps at different energies. We also investigate the effect introducing a kappa-distribution has on the global interaction between the SW and the local interstellar medium (LISM), and find that there is generally an increase in energy transport from the heliosphere into the LISM, due to the modified profile of ENA's energies. This results in a termination shock that moves out by 4 AU, a heliopause that moves in by 9 AU and a bow shock 25 AU farther out, in the nose direction. △ Less

Submitted 17 March, 2008; originally announced March 2008.

Probing Heliospheric Asymmetries with an MHD-kinetic Model

Authors: Nikolai V. Pogorelov, Jacob Heerikhuisen, Gary P. Zank

Abstract: New solar wind data from the Voyager 1 and Voyager 2 spacecraft, together with the SOHO SWAN measurements of the direction that neutral hydrogen enters into the inner heliosheath and neutral helium measurements provided by multiple observations are expected to provide more reliable constraints on the ionization ratio of the local interstellar medium (LISM) and the direction and magnitude of the… ▽ More New solar wind data from the Voyager 1 and Voyager 2 spacecraft, together with the SOHO SWAN measurements of the direction that neutral hydrogen enters into the inner heliosheath and neutral helium measurements provided by multiple observations are expected to provide more reliable constraints on the ionization ratio of the local interstellar medium (LISM) and the direction and magnitude of the interstellar magnetic field (ISMF). In this paper we use currently the most sophisticated numerical model of the heliospheric interface, which is based on an MHD treatment of the ion flow and kinetic modeling of neutral particles, to analyze an ISMF-induced asymmetry of the heliosphere in the presence of the interplanetary magnetic field and neutral particles. It is shown that secondary hydrogen atoms modify the LISM properties leading to its shock-free deceleration at the heliopause. We determine the deflection of hydrogen atoms from their original trajectory in the unperturbed LISM and show that it occurs not only in the plane defined by the ISMF and LISM velocity vectors, but also, to a lesser extent, perpendicular to this plane. We also consider the possibility of using 2-3 kHz radio emission data to further constrain the ISMF direction. △ Less

Submitted 27 January, 2008; originally announced January 2008.

Comments: Paper accepted for publication in the Astrophysical Journal Letters

Heliospheric Response to Different Possible Interstellar Environments

Authors: Hans-R. Mueller, Priscilla C. Frisch, Vladimir Florinski, Gary P. Zank

Abstract: At present, the heliosphere is embedded in a warm low density interstellar cloud that belongs to a cloud system flowing through the local standard of rest with a velocity near ~18 km/s. The velocity structure of the nearest interstellar material (ISM), combined with theoretical models of the local interstellar cloud (LIC), suggest that the Sun passes through cloudlets on timescales of < 10^3 - 1… ▽ More At present, the heliosphere is embedded in a warm low density interstellar cloud that belongs to a cloud system flowing through the local standard of rest with a velocity near ~18 km/s. The velocity structure of the nearest interstellar material (ISM), combined with theoretical models of the local interstellar cloud (LIC), suggest that the Sun passes through cloudlets on timescales of < 10^3 - 10^4 yr, so the heliosphere has been, and will be, exposed to different interstellar environments over time. By means of a multi-fluid model that treats plasma and neutral hydrogen self-consistently, the interaction of the solar wind with a variety of partially ionized ISM is investigated, with the focus on low density cloudlets such as are currently near the Sun. Under the assumption that the basic solar wind parameters remain/were as they are today, a range of ISM parameters (from cold neutral to hot ionized, with various densities and velocities) is considered. In response to different interstellar boundary conditions, the heliospheric size and structure change, as does the abundance of interstellar and secondary neutrals in the inner heliosphere, and the cosmic ray level in the vicinity of Earth. Some empirical relations between interstellar parameters and heliospheric boundary locations, as well as neutral densities, are extracted from the models. △ Less

Submitted 26 July, 2006; originally announced July 2006.

Comments: 24 pages, 9 figures, 2 tables

Journal ref: Astrophys.J.647:1491-1505,2006

New Mass Loss Measurements from Astrospheric Lyman-alpha Absorption

Authors: Brian E. Wood, Hans-Reinhard Mueller, Gary P. Zank, Jeffrey L. Linsky, Seth Redfield

Abstract: Measurements of stellar mass loss rates are used to assess how wind strength varies with coronal activity and age for solar-like stars. Mass loss generally increases with activity, but we find evidence that winds suddenly weaken at a certain activity threshold. Very active stars are often observed to have polar starspots, and we speculate that the magnetic field geometry associated with these sp… ▽ More Measurements of stellar mass loss rates are used to assess how wind strength varies with coronal activity and age for solar-like stars. Mass loss generally increases with activity, but we find evidence that winds suddenly weaken at a certain activity threshold. Very active stars are often observed to have polar starspots, and we speculate that the magnetic field geometry associated with these spots may be inhibiting the winds. Our inferred mass-loss/age relation represents an empirical estimate of the history of the solar wind. This result is important for planetary studies as well as solar/stellar astronomy, since solar wind erosion may have played an important role in the evolution of planetary atmospheres. △ Less

Submitted 16 June, 2005; originally announced June 2005.

Comments: 12 pages, 3 figures, AASTEX v5.0 plus EPSF extensions in mkfig.sty; accepted by ApJ Letters

Journal ref: Astrophys.J.628:L143-L146,2005

Stellar Lyman-alpha Emission Lines in the Hubble Space Telescope Archive: Intrinsic Line Fluxes and Absorption from the Heliosphere and Astrospheres

Authors: Brian E. Wood, Seth Redfield, Jeffrey L. Linsky, Hans-Reinhard Mueller, Gary P. Zank

Abstract: We search the Hubble Space Telescope (HST) archive for previously unanalyzed observations of stellar H I Lyman-alpha emission lines, our primary purpose being to look for new detections of Lyman-alpha absorption from the outer heliosphere, and to also search for analogous absorption from the astrospheres surrounding the observed stars. The astrospheric absorption is of particular interest becaus… ▽ More We search the Hubble Space Telescope (HST) archive for previously unanalyzed observations of stellar H I Lyman-alpha emission lines, our primary purpose being to look for new detections of Lyman-alpha absorption from the outer heliosphere, and to also search for analogous absorption from the astrospheres surrounding the observed stars. The astrospheric absorption is of particular interest because it can be used to study solar-like stellar winds that are otherwise undetectable. We find and analyze 33 HST Lyman-alpha spectra in the archive. All the spectra were taken with the E140M grating of the Space Telescope Imaging Spectrograph (STIS) instrument on board HST. The HST/STIS spectra yield 4 new detections of heliospheric absorption (70 Oph, Xi Boo, 61 Vir, and HD 165185) and 7 new detections of astrospheric absorption (EV Lac, 70 Oph, Xi Boo, 61 Vir, Delta Eri, HD 128987, and DK UMa), doubling the previous number of heliospheric and astrospheric detections. When combined with previous results, 10 of 17 lines of sight within 10 pc yield detections of astrospheric absorption. This high detection fraction implies that most of the ISM within 10 pc must be at least partially neutral, since the presence of H I within the ISM surrounding the observed star is necessary for an astrospheric detection. In contrast, the detection percentage is only 9.7% (3 out of 31) for stars beyond 10 pc. Our Lyman-alpha analyses provide measurements of ISM H I and D I column densities for all 33 lines of sight, and we discuss some implications of these results. Finally, we measure chromospheric Lyman-alpha fluxes from the observed stars. We use these fluxes to determine how Lyman-alpha flux correlates with coronal X-ray and chromospheric Mg II emission, and we also study how Lyman-alpha emission depends on stellar rotation. △ Less

Submitted 16 March, 2005; originally announced March 2005.

Comments: 56 pages, 15 figures; AASTEX v5.0 plus EPSF extensions in mkfig.sty; accepted by ApJS

Investigating a Possible Spectral Signature of the Wind-ISM Interaction Region of Alpha Tau

Authors: Brian E. Wood, Graham M. Harper, Hans-Reinhard Mueller, Gary P. Zank

Abstract: Ultraviolet spectra from the GHRS instrument on board the Hubble Space Telescope reveal the presence of a mysterious absorption feature in the Mg II h & k lines of the nearby (d=20.0 pc) K5 III star Alpha Tau. The narrow absorption looks like an interstellar absorption feature but it is in the wrong location based on our knowledge of the local ISM flow vector. Since the absorption is close to th… ▽ More Ultraviolet spectra from the GHRS instrument on board the Hubble Space Telescope reveal the presence of a mysterious absorption feature in the Mg II h & k lines of the nearby (d=20.0 pc) K5 III star Alpha Tau. The narrow absorption looks like an interstellar absorption feature but it is in the wrong location based on our knowledge of the local ISM flow vector. Since the absorption is close to the rest frame of the star, it has been interpreted as being from the interaction region between Alpha Tau's massive, cool wind and the interstellar medium, i.e., Alpha Tau's "astrosphere". We compute hydrodynamic models of the Alpha Tau astrosphere in order to see if the models can reproduce the Mg II absorption feature. These models do predict that stellar wind material heated, decelerated, and compressed after passing through a termination shock a few thousand AU from the star should produce a Mg II absorption feature with about the right width at roughly the right velocity. However, our first models underestimate the Mg II column density by an order of magnitude. A much larger parameter search is necessary to see whether the observed Mg II absorption can be reproduced by acceptable changes to the adopted stellar wind and ISM properties. △ Less

Submitted 4 September, 2002; originally announced September 2002.

Comments: 8 pages, 4 figures; LaTex plus PASP style file newpasp.sty; to appear in conference proceedings for IAU symposium 210 on Modelling of Stellar Atmospheres