-
Interaction between a Coronal Mass Ejection and Comet 67P/Churyumov-Gerasimenko
Authors:
Zhenguang Huang,
Gabor Toth,
Tamas I. Gombosi,
Michael R. Combi,
Xianzhe Jia,
Yinsi Shou,
Valeriy Tenishev,
Kathrin Altwegg,
Martin Rubin
Abstract:
The interaction between a Coronal Mass Ejection (CME) and a comet has been observed several times by in-situ observations from the Rosetta Plasma Consortium (RPC), which is designed to investigate the cometary magnetosphere of comet 67P/Churyumov-Gerasimenko (CG). Goetz et al. (2019) reported a magnetic field of up to 300 nT measured in the inner coma, which is among the largest interplanetary mag…
▽ More
The interaction between a Coronal Mass Ejection (CME) and a comet has been observed several times by in-situ observations from the Rosetta Plasma Consortium (RPC), which is designed to investigate the cometary magnetosphere of comet 67P/Churyumov-Gerasimenko (CG). Goetz et al. (2019) reported a magnetic field of up to 300 nT measured in the inner coma, which is among the largest interplanetary magnetic fields observed in the solar system. They suggested the large magnetic field observations in the inner coma come from magnetic field pile-up regions, which are generated by the interaction between a CME and/or corotating interaction region and the cometary magnetosphere. However, the detailed interaction between a CME and the cometary magnetosphere of comet CG in the inner coma has not been investigated by numerical simulations yet. In this manuscript, we will use a numerical model to simulate the interaction between comet CG and a Halloween class CME and investigate its magnetospheric response to the CME. We find that the plasma structures change significantly during the CME event, and the maximum value of the magnetic field strength is more than 500nT close to the nucleus. Virtual satellites at similar distances as Rosetta show that the magnetic field strength can be as large as 250nT, which is slightly less than what Goetz et al. (2019) reported.
△ Less
Submitted 8 March, 2024;
originally announced March 2024.
-
Neutral gas coma dynamics: modeling of flows and attempts to link inner coma structures to properties of the nucleus
Authors:
Raphael Marschall,
Björn J. R. Davidsson,
Martin Rubin,
Valeriy Tenishev
Abstract:
Deriving properties of cometary nuclei from coma data is of significant importance for our understanding of cometary activity and has implications beyond. Ground-based data represent the bulk of measurements available for comets. Yet, to date these observations only access a comet's gas and dust coma at rather large distances from the surface and do not directly observe its surface. In contrast, s…
▽ More
Deriving properties of cometary nuclei from coma data is of significant importance for our understanding of cometary activity and has implications beyond. Ground-based data represent the bulk of measurements available for comets. Yet, to date these observations only access a comet's gas and dust coma at rather large distances from the surface and do not directly observe its surface. In contrast, spacecraft missions are one of the only tools that gain direct access to surface measurements. However, such missions are limited to roughly one per decade. We can overcome these challenges by recognising that the coma contains information about the nucleus's properties. In particular, the near-surface gas environment is most representative of the nucleus. It can inform us about the composition, regionality of activity, and sources of coma features and how they link to the topography, morphology, or other surface properties. The inner coma data is a particularly good proxy because it has not yet, or only marginally, been contaminated by coma chemistry or secondary gas sources. Additionally, the simultaneous observation of the innermost coma with the surface provides the potential to make a direct link between coma measurements and the nucleus. If we hope to link outer coma measurements obtained by Earth-based telescopes to the surface, we must first understand how the inner coma measurements are linked to the surface. Numerical models that describe the flow from the surface into the immediate surroundings are needed to make this connection. This chapter focuses on the advances made to understand the flow of the neutral gas coma from the surface to distances up to a few tens of nuclei radii. We describe both simple/heuristic models and state-of-the-art physically consistent models. Model limitations and what they each are best suited for are discussed.
△ Less
Submitted 10 August, 2023;
originally announced August 2023.
-
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.
-
Exosphere Modeling of Proxima b: A Case Study of Photochemical Escape with a Venus-like Atmosphere
Authors:
Yuni Lee,
Chuanfei Dong,
Valeriy Tenishev
Abstract:
Exoplanets orbiting M-dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our Solar System, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar…
▽ More
Exoplanets orbiting M-dwarfs within habitable zones are exposed to stellar environments more extreme than that terrestrial planets experience in our Solar System, which can significantly impact the atmospheres of the exoplanets and affect their habitability and sustainability. This study provides the first prediction of hot oxygen corona structure and the associated photochemical loss from a 1 bar CO2-dominated atmosphere of a Venus-like rocky exoplanet, where dissociative recombination of O2+ ions is assumed to be the major source reaction for the escape of neutral O atoms and formation of the hot O corona (or exospheres) as on Mars and Venus. We employ a 3D Monte Carlo code to simulate the exosphere of Proxima Centauri b (PCb) based on the ionosphere simulated by a 3D magnetohydrodynamic model. Our simulation results show that variability of the stellar wind dynamic pressure over one orbital period of PCb does not affect the overall spatial structure of the hot O corona but contributes to the change in the global hot O escape rate that varies by an order of magnitude. The escape increases dramatically when the planet possesses its intrinsic magnetic fields as the ionosphere becomes more extended with the presence of a global magnetic field. The extended hot O corona may lead to a more extended H exosphere through collisions between thermal H and hot O, which exemplifies the importance of considering nonthermal populations in exospheres to interpret future observations.
△ Less
Submitted 14 September, 2021;
originally announced September 2021.
-
Magnetohydrodynamic with Adaptively Embedded Particle-in-Cell model: MHD-AEPIC
Authors:
Yinsi Shou,
Valeriy Tenishev,
Yuxi Chen,
Gabor Toth,
Natalia Ganushkina
Abstract:
Space plasma simulations have seen an increase in the use of magnetohydrodynamic (MHD) with embedded Particle-in-Cell (PIC) models. This combined MHD-EPIC algorithm simulates some regions of interest using the kinetic PIC method while employing the MHD description in the rest of the domain. The MHD models are highly efficient and their fluid descriptions are valid for most part of the computationa…
▽ More
Space plasma simulations have seen an increase in the use of magnetohydrodynamic (MHD) with embedded Particle-in-Cell (PIC) models. This combined MHD-EPIC algorithm simulates some regions of interest using the kinetic PIC method while employing the MHD description in the rest of the domain. The MHD models are highly efficient and their fluid descriptions are valid for most part of the computational domain, thus making large-scale global simulations feasible. However, in practical applications, the regions where the kinetic effects are critical can be changing, appearing, disappearing and moving in the computational domain. If a static PIC region is used, this requires a much larger PIC domain than actually needed, which can increase the computational cost dramatically.
To address the problem, we have developed a new method that is able to dynamically change the region of the computational domain where a PIC model is applied. We have implemented this new MHD with Adaptively Embedded PIC (MHD-AEPIC) algorithm using the BATS-R-US Hall MHD and the Adaptive Mesh Particle Simulator (AMPS) as the semi-implicit PIC models. We describe the algorithm and present a test case of two merging flux ropes to demonstrate its accuracy. The implementation uses dynamic allocation/deallocation of memory and load balancing for efficient parallel execution. We evaluate the performance of MHD-AEPIC compared to MHD-EPIC and the scaling properties of the model to large number of computational cores.
△ Less
Submitted 11 August, 2021;
originally announced August 2021.
-
What Sustained Multi-Disciplinary Research Can Achieve: The Space Weather Modeling Framework
Authors:
Tamas I. Gombosi,
Yuxi Chen,
Alex Glocer,
Zhenguang Huang,
Xianzhe Jia,
Michael W. Liemohn,
Ward B. Manchester,
Tuija Pulkkinen,
Nishtha Sachdeva,
Qusai Al Shidi,
Igor V. Sokolov,
Judit Szente,
Valeriy Tenishev,
Gabor Toth,
Bart van der Holst,
Daniel T. Welling,
Lulu Zhao,
Shasha Zou
Abstract:
MHD-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environ…
▽ More
MHD-based global space weather models have mostly been developed and maintained at academic institutions. While the "free spirit" approach of academia enables the rapid emergence and testing of new ideas and methods, the lack of long-term stability and support makes this arrangement very challenging. This paper describes a successful example of a university-based group, the Center of Space Environment Modeling (CSEM) at the University of Michigan, that developed and maintained the Space Weather Modeling Framework (SWMF) and its core element, the BATS-R-US extended MHD code. It took a quarter of a century to develop this capability and reach its present level of maturity that makes it suitable for research use by the space physics community through the Community Coordinated Modeling Center (CCMC) as well as operational use by the NOAA Space Weather Prediction Center (SWPC).
△ Less
Submitted 27 May, 2021;
originally announced May 2021.
-
The Confinement of the Heliosheath Plasma by the Solar Magnetic Field as Revealed by Energetic Neutral Atom Simulations
Authors:
M. Kornbleuth,
M. Opher,
A. T. Michael,
J. M. Sokol,
G. Toth,
V. Tenishev,
J. F. Drake
Abstract:
Traditionally, the solar magnetic field has been considered to have a negligible effect in the outer regions of the heliosphere. Recent works have shown that the solar magnetic field may play a crucial role in collimating the plasma in the heliosheath. Interstellar Boundary Explorer (IBEX) observations of the heliotail indicated a latitudinal structure varying with energy in the energetic neutral…
▽ More
Traditionally, the solar magnetic field has been considered to have a negligible effect in the outer regions of the heliosphere. Recent works have shown that the solar magnetic field may play a crucial role in collimating the plasma in the heliosheath. Interstellar Boundary Explorer (IBEX) observations of the heliotail indicated a latitudinal structure varying with energy in the energetic neutral atom (ENA) fluxes. At energies ~1 keV, the ENA fluxes show an enhancement at low latitudes and a deficit of ENAs near the poles. At energies >2.7 keV, ENA fluxes had a deficit within low latitudes, and lobes of higher ENA flux near the poles. This ENA structure was initially interpreted to be a result of the latitudinal profile of the solar wind during solar minimum. We extend the work of Kornbleuth et al. (2018) by using solar minimum-like conditions and the recently developed SHIELD model. The SHIELD model couples the magnetohydrodynamic (MHD) plasma solution with a kinetic description of neutral hydrogen. We show that while the latitudinal profile of the solar wind during solar minimum contributes to the lobes in ENA maps, the collimation by the solar magnetic field is important in creating and shaping the two high latitude lobes of enhanced ENA flux observed by IBEX. This is the first work to explore the effect of the changing solar magnetic field strength on ENA maps. Our findings suggest that IBEX is providing the first observational evidence of the collimation of the heliosheath plasma by the solar magnetic field.
△ Less
Submitted 13 May, 2020;
originally announced May 2020.
-
The Solar-wind with Hydrogen Ion Exchange and Large-scale Dynamics (SHIELD) model: A Self-Consistent Kinetic-MHD Model of the Outer Heliosphere
Authors:
Adam T. Michael,
Merav Opher,
Gabor Toth,
Valeriy Tenishev,
Dmitry Borovikov
Abstract:
Neutral hydrogen has been shown to greatly impact the plasma flow in the heliopshere and the location of the heliospheric boundaries. We present the results of the Solar-wind with Hydrogen Ion Exchange and Large-scale Dynamics (SHIELD) model, a new, self-consistent, kinetic-MHD model of the outer heliosphere within the Space Weather Modeling Framework. The charge-exchange mean free path is on orde…
▽ More
Neutral hydrogen has been shown to greatly impact the plasma flow in the heliopshere and the location of the heliospheric boundaries. We present the results of the Solar-wind with Hydrogen Ion Exchange and Large-scale Dynamics (SHIELD) model, a new, self-consistent, kinetic-MHD model of the outer heliosphere within the Space Weather Modeling Framework. The charge-exchange mean free path is on order of the size of the heliosphere; therefore, the neutral atoms cannot be described as a fluid. The SHIELD model couples the MHD solution for a single plasma fluid to the kinetic solution from for neutral hydrogen atoms streaming through the system. The kinetic code is based on the Adaptive Mesh Particle Simulator (AMPS), a Monte Carlo method for solving the Boltzmann equation. The SHIELD model accurately predicts the increased filtration of interstellar neutrals into the heliosphere. In order to verify the correct implementation within the model, we compare the results of the SHIELD model to other, well-established kinetic-MHD models. The SHIELD model matches the neutral hydrogen solution of these studies as well as the shift in all heliospheric boundaries closer to the Sun in comparison the the multi-fluid treatment of the neutral hydrogen atoms. Overall the SHIELD model shows excellent agreement to these models and is a significant improvement to the fluid treatment of interstellar hydrogen.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
High Resolution Finite Volume Method for Kinetic Equations with Poisson Brackets
Authors:
Igor V. Sokolov,
Haomin Sun,
Gabor Toth,
Zhenguang Huang,
Valeriy Tenishev,
Lulu Zhao,
Jozsef Kota,
Ofer Cohen,
Tamas Gombosi
Abstract:
Simulation of plasmas in electromagnetic fields requires numerical solution of a kinetic equation that describes the time evolution of the particle distribution function. In this paper we propose a finite volume scheme based on integral relation for Poisson brackets to solve the Liouville equation, the most fundamental kinetic equation. The proposed scheme conserves the number of particles, mainta…
▽ More
Simulation of plasmas in electromagnetic fields requires numerical solution of a kinetic equation that describes the time evolution of the particle distribution function. In this paper we propose a finite volume scheme based on integral relation for Poisson brackets to solve the Liouville equation, the most fundamental kinetic equation. The proposed scheme conserves the number of particles, maintains the total-variation-diminishing (TVD) property, and provides high-quality numerical results. Other types of kinetic equations may be also formulated in terms of Poisson brackets and solved with the proposed method including the transport equations describing the acceleration and propagation of Solar Energetic Particles (SEPs), which is of practical importance, since the high energy SEPs produce radiation hazards. The proposed scheme is demonstrated to be accurate and efficient, which makes it applicable to global simulation systems analyzing space weather.
△ Less
Submitted 26 November, 2022; v1 submitted 25 October, 2019;
originally announced October 2019.
-
The Surface Distributions of the Production of the Major Volatile Species, H2O, CO2, CO and O2, from the Nucleus of Comet 67P/Churyumov-Gerasimenko throughout the Rosetta Mission as Measured by the ROSINA Double Focusing Mass Spectrometer
Authors:
Michael Combi,
Yinsi Shou,
Nicolas Fougere,
Valeriy Tenishev,
Kathrin Altwegg,
Martin Rubin,
Dominique Bockelée-Morvan,
Fabrizio Capaccioni,
Yu-Chi Cheng,
Uwe Fink,
Tamas Gombosi,
Kenneth C. Hansen,
Zhenguang Huang,
David Marshall,
Gabor Toth
Abstract:
The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) suite of instruments operated throughout the over two years of the Rosetta mission operations in the vicinity of comet 67P/Churyumov-Gerasimenko. It measured gas densities and composition throughout the comet's atmosphere, or coma. Here we present two-years' worth of measurements of the relative densities of the four major vola…
▽ More
The Rosetta Orbiter Spectrometer for Ion and Neutral Analysis (ROSINA) suite of instruments operated throughout the over two years of the Rosetta mission operations in the vicinity of comet 67P/Churyumov-Gerasimenko. It measured gas densities and composition throughout the comet's atmosphere, or coma. Here we present two-years' worth of measurements of the relative densities of the four major volatile species in the coma of the comet, H2O. CO2, CO and O2, by one of the ROSINA sub-systems called the Double Focusing Mass Spectrometer (DFMS). The absolute total gas densities were provided by the Comet Pressure Sensor (COPS), another ROSINA sub-system. DFMS is a very high mass resolution and high sensitivity mass spectrometer able to resolve at a tiny fraction of an atomic mass unit. We have analyzed the combined DFMS and COPS measurements using an inversion scheme based on spherical harmonics that solves for the distribution of potential surface activity of each species as the comet rotates, changing solar illumination, over short intervals and as the comet changes distance from the sun and orientation of its spin axis over long time intervals. We also use the surface boundary conditions derived from the inversion scheme to simulate the whole coma with our fully kinetic Direct Simulation Monte Carlo model and calculate the production rates of the four major species throughout the mission. We compare the derived production rates with revised remote sensing observations by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) as well as with published observations from the Microwave Instrument for the Rosetta Orbiter (MIRO). Finally we use the variation of the surface production of the major species to calculate the total mass loss over the mission and, for different estimates of the dust/gas ratio, calculate the variation of surface loss over the nucleus.
△ Less
Submitted 4 September, 2019;
originally announced September 2019.
-
A comparison between the two lobes of comet 67P/Churyumov-Gerasimenko based on D/H ratios in H2O measured with the Rosetta/ROSINA DFMS
Authors:
Isaac R. H. G. Schroeder I,
Kathrin Altwegg,
Hans Balsiger,
Jean-Jacques Berthelier,
Michael R. Combi,
Johan De Keyser,
Björn Fiethe,
Stephen A. Fuselier,
Tamas I. Gombosi,
Kenneth C. Hansen,
Martin Rubin,
Yinsi Shou,
Valeriy M. Tenishev,
Thierry Sémon,
Susanne F. Wampfler,
Peter Wurz
Abstract:
The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its…
▽ More
The nucleus of the Jupiter-family comet 67P/Churyumov-Gerasimenko was discovered to be bi-lobate in shape when the European Space Agency spacecraft Rosetta first approached it in July 2014. The bi-lobate structure of the cometary nucleus has led to much discussion regarding the possible manner of its formation and on how the composition of each lobe might compare with that of the other. During its two-year-long mission from 2014 to 2016, Rosetta remained in close proximity to 67P/Churyumov-Gerasimenko, studying its coma and nucleus in situ. Based on lobe-specific measurements of HDO and H2O performed with the ROSINA DFMS mass spectrometer on board Rosetta, the Deuterium-to-Hydrogen ratios in water from the two lobes could be compared. No appreciable difference was observed, suggesting that both lobes formed in the same region and are homogeneous in their Deuterium-to-Hydrogen ratios.
△ Less
Submitted 16 July, 2019;
originally announced July 2019.
-
Modeling Martian Atmospheric Losses over Time: Implications for Exoplanetary Climate Evolution and Habitability
Authors:
Chuanfei Dong,
Yuni Lee,
Yingjuan Ma,
Manasvi Lingam,
Stephen Bougher,
Janet Luhmann,
Shannon Curry,
Gabor Toth,
Andrew Nagy,
Valeriy Tenishev,
Xiaohua Fang,
David Mitchell,
David Brain,
Bruce Jakosky
Abstract:
In this Letter, we make use of sophisticated 3D numerical simulations to assess the extent of atmospheric ion and photochemical losses from Mars over time. We demonstrate that the atmospheric ion escape rates were significantly higher (by more than two orders of magnitude) in the past at $\sim 4$ Ga compared to the present-day value owing to the stronger solar wind and higher ultraviolet fluxes fr…
▽ More
In this Letter, we make use of sophisticated 3D numerical simulations to assess the extent of atmospheric ion and photochemical losses from Mars over time. We demonstrate that the atmospheric ion escape rates were significantly higher (by more than two orders of magnitude) in the past at $\sim 4$ Ga compared to the present-day value owing to the stronger solar wind and higher ultraviolet fluxes from the young Sun. We found that the photochemical loss of atomic hot oxygen dominates over the total ion loss at the current epoch whilst the atmospheric ion loss is likely much more important at ancient times. We briefly discuss the ensuing implications of high atmospheric ion escape rates in the context of ancient Mars, and exoplanets with similar atmospheric compositions around young solar-type stars and M-dwarfs.
△ Less
Submitted 14 May, 2018;
originally announced May 2018.
-
Solar wind interaction with the Martian upper atmosphere: Roles of the cold thermosphere and hot oxygen corona
Authors:
Chuanfei Dong,
Stephen W. Bougher,
Yingjuan Ma,
Yuni Lee,
Gabor Toth,
Andrew F. Nagy,
Xiaohua Fang,
Janet Luhmann,
Michael W. Liemohn,
Jasper S. Halekas,
Valeriy Tenishev,
David J. Pawlowski,
Michael R. Combi
Abstract:
We study roles of the thermosphere and exosphere on the Martian ionospheric structure and ion escape rates in the process of the solar wind-Mars interaction. We employ a four-species multifluid MHD (MF-MHD) model to simulate the Martian ionosphere and magnetosphere. The $cold$ thermosphere background is taken from the Mars Global Ionosphere Thermosphere Model (M-GITM) and the $hot$ oxygen exospher…
▽ More
We study roles of the thermosphere and exosphere on the Martian ionospheric structure and ion escape rates in the process of the solar wind-Mars interaction. We employ a four-species multifluid MHD (MF-MHD) model to simulate the Martian ionosphere and magnetosphere. The $cold$ thermosphere background is taken from the Mars Global Ionosphere Thermosphere Model (M-GITM) and the $hot$ oxygen exosphere is adopted from the Mars exosphere Monte Carlo model - Adaptive Mesh Particle Simulator (AMPS). A total of four cases with the combination of 1D (globally averaged) and 3D thermospheres and exospheres are studied.
The ion escape rates calculated by adopting 1D and 3D atmospheres are similar; however, the latter are required to adequately reproduce MAVEN ionospheric observations. In addition, our simulations show that the 3D hot oxygen corona plays an important role in preventing planetary molecular ions (O$_2^+$ and CO$_2^+$) escaping from Mars, mainly resulting from the mass loading of the high-altitude exospheric O$^+$ ions. The $cold$ thermospheric oxygen atom, however, is demonstrated to be the primary neutral source for O$^+$ ion escape during the relatively weak solar cycle 24.
△ Less
Submitted 22 June, 2018; v1 submitted 3 April, 2018;
originally announced April 2018.
-
Hall Effect in the coma of 67P/Churyumov-Gerasimenko
Authors:
Z. Huang,
G. Toth,
T. I. Gombosi,
X. Jia,
M. R. Combi,
K. C. Hansen,
N. Fougere,
Y. Shou,
V. Tenishev,
K. Altwegg,
M. Rubin
Abstract:
Magnetohydrodynamics simulations have been carried out in studying the solar wind and cometary plasma interactions for decades. Various plasma boundaries have been simulated and compared well with observations for comet 1P/Halley. The Rosetta mission, which studies comet 67P/Churyumov-Gerasimenko, challenges our understanding of the solar wind and comet interactions. The Rosetta Plasma Consortium…
▽ More
Magnetohydrodynamics simulations have been carried out in studying the solar wind and cometary plasma interactions for decades. Various plasma boundaries have been simulated and compared well with observations for comet 1P/Halley. The Rosetta mission, which studies comet 67P/Churyumov-Gerasimenko, challenges our understanding of the solar wind and comet interactions. The Rosetta Plasma Consortium observed regions of very weak magnetic field outside the predicted diamagnetic cavity. In this paper, we simulate the inner coma with the Hall magnetohydrodynamics equations and show that the Hall effect is important in the inner coma environment. The magnetic field topology becomes complex and magnetic reconnection occurs on the dayside when the Hall effect is taken into account. The magnetic reconnection on the dayside can generate weak magnetic filed regions outside the global diamagnetic cavity, which may explain the Rosetta Plasma Consortium observations. We conclude that the substantial change in the inner coma environment is due to the fact that the ion inertial length (or gyro radius) is not much smaller than the size of the diamagnetic cavity.
△ Less
Submitted 11 January, 2018;
originally announced January 2018.