Showing 1–50 of 51 results for author: Sokolov, I

Solar Wind with Field Lines and Energetic Particles (SOFIE) Model: Application to Historical Solar Energetic Particle Events

Authors: Lulu Zhao, Igor Sokolov, Tamas Gombosi, David Lario, Kathryn Whitman, Zhenguang Huang, Gabor Toth, Ward Manchester, Bart van der Holst, Nishtha Sachdeva

Abstract: In this paper, we demonstrate the applicability of the data-driven and self-consistent solar energetic particle model, Solar-wind with FIeld-lines and Energetic-particles (SOFIE), to simulate acceleration and transport processes of solar energetic particles. SOFIE model is built upon the Space Weather Modeling Framework (SWMF) developed at the University of Michigan. In SOFIE, the background solar… ▽ More In this paper, we demonstrate the applicability of the data-driven and self-consistent solar energetic particle model, Solar-wind with FIeld-lines and Energetic-particles (SOFIE), to simulate acceleration and transport processes of solar energetic particles. SOFIE model is built upon the Space Weather Modeling Framework (SWMF) developed at the University of Michigan. In SOFIE, the background solar wind plasma in the solar corona and interplanetary space is calculated by the Aflvén Wave Solar-atmosphere Model(-Realtime) (AWSoM-R) driven by the near-real-time hourly updated Global Oscillation Network Group (GONG) solar magnetograms. In the background solar wind, coronal mass ejections (CMEs) are launched by placing an imbalanced magnetic flux rope on top of the parent active region, using the Eruptive Event Generator using Gibson-Low model (EEGGL). The acceleration and transport processes are modeled by the Multiple-Field-Line Advection Model for Particle Acceleration (M-FLAMPA). In this work, nine solar energetic particle events (Solar Heliospheric and INterplanetary Environment (SHINE) challenge/campaign events) are modeled. The three modules in SOFIE are validated and evaluated by comparing with observations, including the steady-state background solar wind properties, the white-light image of the CME, and the flux of solar energetic protons, at energies of > 10 MeV. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Three-dimensional, Time-dependent MHD Simulation of Disk-Magnetosphere-Stellar Wind Interaction in a T Tauri, Protoplanetary System

Authors: Ofer Cohen, Cecilia Garraffo, Jeremy Drake, Kristina Monsch, Igor Sokolov, Julian Alvarado-Gomez, Federico Fraschetti

Abstract: We present a three-dimensional, time-dependent, MHD simulation of the short-term interaction between a protoplanetary disk and the stellar corona in a T Tauri system. The simulation includes the stellar magnetic field, self-consistent coronal heating and stellar wind acceleration, and a disk rotating at sub-Keplerian velocity to induce accretion. We find that initially, as the system relaxes from… ▽ More We present a three-dimensional, time-dependent, MHD simulation of the short-term interaction between a protoplanetary disk and the stellar corona in a T Tauri system. The simulation includes the stellar magnetic field, self-consistent coronal heating and stellar wind acceleration, and a disk rotating at sub-Keplerian velocity to induce accretion. We find that initially, as the system relaxes from the assumed initial conditions, the inner part of the disk winds around and moves inward and close to the star as expected. However, the self-consistent coronal heating and stellar wind acceleration build up the original state after some time, significantly pushing the disk out beyond $10R_\star$. After this initial relaxation period, we do not find clear evidence of a strong, steady accretion flow funneled along coronal field lines, but only weak, sporadic accretion. We produce synthetic coronal X-ray line emission light curves which show flare-like increases that are not correlated with accretion events nor with heating events. These variations in the line emission flux are the result of compression and expansion due to disk-corona pressure variations. Vertical disk evaporation evolves above and below the disk. However, the disk - stellar wind boundary stays quite stable, and any disk material that reaches the stellar wind region is advected out by the stellar wind. △ Less

Submitted 24 March, 2023; v1 submitted 23 March, 2023; originally announced March 2023.

Comments: Accepted to ApJ, 12 pages, 11 figures

A Titov-Démoulin Type Eruptive Event Generator for $β>0$ Plasmas

Authors: Igor V. Sokolov, Tamas I Gombosi

Abstract: We provide exact analytical solutions for the magnetic field produced by prescribed current distributions located inside a toroidal filament of finite thickness. The solutions are expressed in terms of toroidal functions which are modifications of the Legendre functions. In application to the MHD equilibrium of a twisted toroidal current loop in the solar corona, the Grad-Shafranov equation is dec… ▽ More We provide exact analytical solutions for the magnetic field produced by prescribed current distributions located inside a toroidal filament of finite thickness. The solutions are expressed in terms of toroidal functions which are modifications of the Legendre functions. In application to the MHD equilibrium of a twisted toroidal current loop in the solar corona, the Grad-Shafranov equation is decomposed into an analytic solution describing an equilibrium configuration against the pinch-effect from its own current and an approximate solution for an external strapping field to balance the hoop force. Our solutions can be employed in numerical simulations of coronal mass ejections. When superimposed on the background solar coronal magnetic field, the excess magnetic energy of the twisted current loop configuration can be made unstable by applying flux cancellation to reduce the strapping field. Such loss of stability accompanied by the formation of an expanding flux rope is typical for the Titov & Démoulin (1999) eruptive event generator. The main new features of the proposed model are: (i) The filament is filled with finite $β$ plasma with finite mass and energy, (ii) The model describes an equilibrium solution that will spontaneously erupt due to magnetic reconnection of the strapping magnetic field arcade, and (iii) There are analytic expressions connecting the model parameters to the asymptotic velocity and total mass of the resulting CME, providing a way to connect the simulated CME properties to multipoint coronograph observations. △ Less

Submitted 13 April, 2023; v1 submitted 27 December, 2022; originally announced December 2022.

Comments: 20 pages, 7 figures

Solar wind modeling with the Alfven Wave Solar atmosphere Model driven by HMI-based Near-Real-Time maps by the National Solar Observatory

Authors: Nishtha Sachdeva, Ward B Manchester, Igor Sokolov, Zhenguang Huang, Alexander Pevtsov, Luca Bertello, Alexei A. Pevtsov, Gabor Toth, Bart van der Holst

Abstract: We explore model performance for the Alfven Wave Solar atmosphere Model (AWSoM) with near-real-time (NRT) synoptic maps of the photospheric vector magnetic field. These maps, produced by assimilating data from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), use a different method developed at the National Solar Observatory (NSO) to provide a near contemporaneou… ▽ More We explore model performance for the Alfven Wave Solar atmosphere Model (AWSoM) with near-real-time (NRT) synoptic maps of the photospheric vector magnetic field. These maps, produced by assimilating data from the Helioseismic Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO), use a different method developed at the National Solar Observatory (NSO) to provide a near contemporaneous source of data to drive numerical models. Here, we apply these NSO-HMI-NRT maps to simulate three Carrington rotations (CRs): 2107-2108 (centered on 2011/03/07 20:12 CME event), 2123 (integer CR) and 2218--2219 (centered on 2019/07/2 solar eclipse), which together cover a wide range of activity level for solar cycle 24. We show simulation results, which reproduce both extreme ultraviolet emission (EUV) from the low corona while simultaneously matching in situ observations at 1 au as well as quantify the total unsigned open magnetic flux from these maps. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Modeling the Solar Wind During Different Phases of the Last Solar Cycle

Authors: Zhenguang Huang, Gabor Toth, Nishtha Sachdeva, Lulu Zhao, Bartholomeus van der Holst, Igor Sokolov, Ward Manchester, Tamas Gombosi

Abstract: We describe our first attempt to systematically simulate the solar wind during different phases of the last solar cycle with the Alfvén Wave Solar atmosphere Model (AWSoM) developed at the University of Michigan. Key to this study is the determination of the optimal values of one of the most important input parameters of the model, the Poynting flux parameter, which prescribes the energy flux pass… ▽ More We describe our first attempt to systematically simulate the solar wind during different phases of the last solar cycle with the Alfvén Wave Solar atmosphere Model (AWSoM) developed at the University of Michigan. Key to this study is the determination of the optimal values of one of the most important input parameters of the model, the Poynting flux parameter, which prescribes the energy flux passing through the chromospheric boundary of the model in the form of Alfvén wave turbulence. It is found that the optimal value of the Poynting flux parameter is correlated with the area of the open magnetic field regions with the Spearman's correlation coefficient of 0.96 and anti-correlated with the average unsigned radial component of the magnetic field with the Spearman's correlation coefficient of -0.91. Moreover, the Poynting flux in the open field regions is approximately constant in the last solar cycle, which needs to be validated with observations and can shed light on how Alfvén wave turbulence accelerates the solar wind during different phases of the solar cycle. Our results can also be used to set the Poynting flux parameter for real-time solar wind simulations with AWSoM. △ Less

Submitted 24 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 15 pages, 1 table, 4 figures

Exoplanet Radio Transits as a Probe for Exoplanetary Magnetic Fields -- Time-dependent MHD Simulations

Authors: Soumitra Hazra, Ofer Cohen, Igor V. Sokolov

Abstract: We perform a series of time dependent Magnetohydrodynamic simulations of the HD 189733 star-planet system in order to predict radio transit modulations due to the interaction between the stellar wind and planetary magnetic field. The simulation combines a model for the stellar corona and wind with an exoplanet that is orbiting the star in a fully dynamic, time-dependent manner. Our simulations gen… ▽ More We perform a series of time dependent Magnetohydrodynamic simulations of the HD 189733 star-planet system in order to predict radio transit modulations due to the interaction between the stellar wind and planetary magnetic field. The simulation combines a model for the stellar corona and wind with an exoplanet that is orbiting the star in a fully dynamic, time-dependent manner. Our simulations generate synthetic radio images that enable us to obtain synthetic radio lightcurves in different frequencies. We find a clear evidence for the planetary motion in the radio light curves. Moreover, we find specific repeated features in the light curves that are attributed to the passage of the planetary magnetosphere in front of the star during transit. More importantly, we find a clear dependence in the magnitude and phase of these lightcurve features on the strength of the planetary magnetic field. Our work demonstrates that if radio transits could be observed, they could indeed provide information about the magnetic field strength of the transiting exoplanet. Future work to parameterize these lightcurve features and their dependence on the planetary field strength would provide tools to search for these features in radio observations datasets. As we only consider the thermal radio emission from the host star for our study, very sensitive radio interferometers are necessary to detect these kinds of planetary transit in radio. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 19 Pages, 10 Figures

Stream-Aligned Magnetohydrodynamics for Solar Wind Simulations

Authors: Igor V Sokolov, Lulu Zhao, Tamas I Gombosi

Abstract: We present a reduced magnetohydrodynamic (MHD) mathematical model describing the dynamical behavior of highly conducting plasmas with frozen-in magnetic fields, constrained by the assumption that, there exists a frame of reference, where the magnetic field vector, $\mathbf{B}$, is aligned with the plasma velocity vector, $\mathbf{u}$, at each point. We call this solution "stream-aligned MHD" (SA-M… ▽ More We present a reduced magnetohydrodynamic (MHD) mathematical model describing the dynamical behavior of highly conducting plasmas with frozen-in magnetic fields, constrained by the assumption that, there exists a frame of reference, where the magnetic field vector, $\mathbf{B}$, is aligned with the plasma velocity vector, $\mathbf{u}$, at each point. We call this solution "stream-aligned MHD" (SA-MHD). Within the framework of this model, the electric field, $\mathbf{E} = -\mathbf{u} \times \mathbf{B} \equiv 0$, in the induction equation vanishes identically and so does the electromagnetic energy flux (Poynting flux), $\mathbf{E}\times\mathbf{B}\equiv0$, in the energy equation. At the same time, the force effect from the magnetic field on the plasma motion (the Ampere force) is fully taken into account in the momentum equation. Any steady-state solution of the proposed model is a legitimate solution of the full MHD system of equations. However, the converse statement is not true: in an arbitrary steady-state magnetic field the electric field does not have to vanish identically (its curl has to, though). Specifically, realistic tree-dimensional solutions for the steady-state ("ambient") solar atmosphere in the form of so-called Parker spirals, can be efficiently generated within the stream-aligned MHD (SA-MHD) with no loss in generality △ Less

Submitted 28 November, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 22 pages, 6 figures. Submitted to The Astrophysical Journal

Journal ref: ApJ 926 102 (2022)

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.

Comments: 105 pages, 36 figures, in press

Journal ref: J. Space Weather and Space Climate, 2021

Potential asteroid discoveries by the ESA Gaia mission: Results from follow-up observations

Authors: B. Carry, W. Thuillot, F. Spoto, P. David, J. Berthier, P. Tanga, F. Mignard, S. Bouquillon, R . A. Mendez, J. -P. Rivet, A. Le Van Suu, A. Dell'Oro, G. Fedorets, B. Frezouls, M. Granvik, J. Guiraud, K. Muinonen, C. Panem, T. Pauwels, W. Roux, G. Walmsley, J. -M. Petit, L. Abe, V. Ayv azian, K. Baillié , et al. (20 additional authors not shown)

Abstract: Since July 2014, the ESA Gaia mission has been surveying the entire sky down to magnitude 20.7 in the visible. In addition to the millions of stars, thousands of Solar System Objects (SSOs) are observed daily. By comparing their positions to those of known objects, a daily processing pipeline filters known objects from potential discoveries. However, owing to Gaia's specific scanning law designed… ▽ More Since July 2014, the ESA Gaia mission has been surveying the entire sky down to magnitude 20.7 in the visible. In addition to the millions of stars, thousands of Solar System Objects (SSOs) are observed daily. By comparing their positions to those of known objects, a daily processing pipeline filters known objects from potential discoveries. However, owing to Gaia's specific scanning law designed for stars, potential newly discovered moving objects are characterized by very few observations, acquired over a limited time. This aspect was recognized early in the design of the Gaia data processing. A daily processing pipeline dedicated to these candidate discoveries was set up to release calls for observations to a network of ground-based telescopes. Their aim is to acquire follow-up astrometry and to characterize these objects. From the astrometry measured by Gaia, preliminary orbital solutions are determined, allowing to predict the position of these potentially new discovered objects in the sky accounting for the large parallax between Gaia and the Earth (separated by 0.01 au). A specific task within the Gaia Consortium has been responsible for the distribution of requests for follow-up observations of potential Gaia SSO discoveries. Since late 2016, these calls for observations (called alerts) are published daily via a Web interface, freely available to anyone world-wide. Between November 2016 and July 2020, over 1700 alerts have been published, leading to the successful recovery of more than 200 objects. Among those, six have provisional designation assigned with the Gaia observations, the others being previously known objects with poorly characterized orbits, precluding identification at the time of Gaia observations. There is a clear trend for objects with a high inclination to be unidentified, revealing a clear bias in the current census of SSOs against high inclination populations. △ Less

Submitted 6 October, 2020; originally announced October 2020.

Comments: Submitted to Astronomy & Astrophysics

Journal ref: A&A 648, A96 (2021)

THESEUS$-$BTA cosmological tests using Multimessenger Gamma-Ray Bursts observations

Authors: S. I. Shirokov, I. V. Sokolov, V. V. Vlasyuk, L. Amati, V. V. Sokolov, Yu. V. Baryshev

Abstract: Modern Multimessenger Astronomy is a powerful instrument for performing cosmological crucial tests of the Standard Cosmological Model in the wide redshift interval up to $z \sim 10$. This is principally important for discussion related to discrepancies between local and global measurements of cosmological parameters. We present a review of multimessenger gamma-ray burst observations currently cond… ▽ More Modern Multimessenger Astronomy is a powerful instrument for performing cosmological crucial tests of the Standard Cosmological Model in the wide redshift interval up to $z \sim 10$. This is principally important for discussion related to discrepancies between local and global measurements of cosmological parameters. We present a review of multimessenger gamma-ray burst observations currently conducted and planed for THESEUS$-$BTA cooperative program. Such observations give a unique opportunity to test the fundamental foundations of cosmological models: gravitation theory; cosmological principle of homogeneity and isotropy of large-scale distribution of matter; and space expansion paradigm. Important role of various selection effects leading to systematic distortions of true cosmological relations is discussed. △ Less

Submitted 6 July, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 17 pages, 4 figures, 106 references, accepted for publication in Astrophys. Bulletin June 10 2020, v.70, iss.3

MSC Class: 85-02

Journal ref: 2020, Astrophysical Bulletin, 75(3), 207-218

High Redshift Long Gamma-Ray Bursts Hubble Diagram as a Test of Basic Cosmological Relations

Authors: S. I. Shirokov, I. V. Sokolov, N. Yu. Lovyagin, L. Amati, Yu. V. Baryshev, V. V. Sokolov, V. L. Gorokhov

Abstract: We examine the prospects of the high redshift Long Gamma Ray Bursts (LGRB) Hubble Diagram as a test of the basic cosmological principles. Analysis of the Hubble Diagram allows us to test several fundamental cosmological principles using the directly observed flux-distance-redshift relation. Modern LGRB data, together with the correlation between the spectral peak energy and the isotropic-equivalen… ▽ More We examine the prospects of the high redshift Long Gamma Ray Bursts (LGRB) Hubble Diagram as a test of the basic cosmological principles. Analysis of the Hubble Diagram allows us to test several fundamental cosmological principles using the directly observed flux-distance-redshift relation. Modern LGRB data, together with the correlation between the spectral peak energy and the isotropic-equivalent radiated energy (the so-called Amati relation) can be used for construction of the Hubble Diagram at the model-independent level. We emphasise observational selection effects, which inevitably exist and distort the theoretically predicted relations. An example is the weak and strong gravitational lensing bias effect for high redshift LGRB in the presence of limited observational sensitivity (Malmquist bias). After bias correction, there is a tendency to vacuum dominated models with $Ω_Λ\rightarrow 0.9$, $Ω_\text{m} \rightarrow 0.1$. Forthcoming gamma-ray observations by the space THESEUS mission together with ground and space based multimessenger facilities will allow us to improve essentially the restrictions on alternative basic principles of cosmological models. △ Less

Submitted 21 June, 2020; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: 16 pages, 8 figures, 3 tables, edited and references corrected, MNRAS accepted

MSC Class: 85-05

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 496, Issue 2, August 2020, Pages 1530-1544

Coupled MHD -- Hybrid Simulations of Space Plasmas

Authors: S. P. Moschou, I. V. Sokolov, O. Cohen, G. Toth, J. J. Drake, Z. Huang, C. Garraffo, J. D. Alvarado-Gómez, T. Gombosi

Abstract: Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of space plasmas. Hybrid codes, on the other hand, capture important physical processes, such as electric currents and effects of finite Larmor radius, but they can… ▽ More Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of space plasmas. Hybrid codes, on the other hand, capture important physical processes, such as electric currents and effects of finite Larmor radius, but they can be used locally only, since the limitations in available computational resources do not allow for their use throughout a global computational domain. In the present work, we present a new coupled scheme which allows to switch blocks in the block-adaptive grids from fluid MHD to hybrid simulations, without modifying the self-consistent computation of the electromagnetic fields acting on fluids (in MHD simulation) or charged ion macroparticles (in hybrid simulation). In this way, the hybrid scheme can refine the description in specified regions of interest without compromising the efficiency of the global MHD code. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: 13 pages, 3 figures, ASTRONUM 2019 refereed proceedings paper (in press)

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.

Comments: 50 pages, 9 figures

Photometric observations of the asteroid 3200 Phaethon using small and middle telescopes

Authors: A. Kartashova, M. Husarik, O. Ivanova, G. Kokhirova, E. Bakanas, I. Sokolov, U. Kh. Khamroev, A. A. Ibragimov

Abstract: The main aim of photometrical observations of the asteroid 3200 Phaethon was searching for its low-level cometary activity (possible coma and/or dust tail) in the pre-perihelion passage. We performed observational runs with telescopes ranging from 0.61-m to 2-m and BVR color imaging. Three longer photometric series were used for modeling of the 3D shape of Phaethon. The color indices and size of t… ▽ More The main aim of photometrical observations of the asteroid 3200 Phaethon was searching for its low-level cometary activity (possible coma and/or dust tail) in the pre-perihelion passage. We performed observational runs with telescopes ranging from 0.61-m to 2-m and BVR color imaging. Three longer photometric series were used for modeling of the 3D shape of Phaethon. The color indices and size of the asteroid were estimated. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: 6 pages, 2 figures, 1 table

The afterglow and kilonova of the short GRB 160821B

Authors: E. Troja, A. J. Castro-Tirado, J. Becerra Gonzalez, Y. Hu, G. S. Ryan, S. B. Cenko, R. Ricci, G. Novara, R. Sanchez-Ramirez, J. A. Acosta-Pulido, M. D. Caballero Garcia, S. Guziy, S. Jeong, A. Y. Lien, I. Marquez, S. B. Pandey, I. H. Park, J. C. Tello, T. Sakamoto, I. V. Sokolov, V. V. Sokolov, A. Tiengo, A. F. Valeev, B. B. Zhang, S. Veilleux

Abstract: GRB 160821B is a short duration gamma-ray burst (GRB) detected and localized by the Neil Gehrels Swift Observatory in the outskirts of a spiral galaxy at z=0.1613, at a projected physical offset of 16 kpc from the galaxy's center. We present X-ray, optical/nIR and radio observations of its counterpart and model them with two distinct components of emission: a standard afterglow, arising from the i… ▽ More GRB 160821B is a short duration gamma-ray burst (GRB) detected and localized by the Neil Gehrels Swift Observatory in the outskirts of a spiral galaxy at z=0.1613, at a projected physical offset of 16 kpc from the galaxy's center. We present X-ray, optical/nIR and radio observations of its counterpart and model them with two distinct components of emission: a standard afterglow, arising from the interaction of the relativistic jet with the surrounding medium, and a kilonova, powered by the radioactive decay of the sub-relativistic ejecta. Broadband modeling of the afterglow data reveals a weak reverse shock propagating backward into the jet, and a likely jet-break at 3.5 d. This is consistent with a structured jet seen slightly off-axis while expanding into a low-density medium. Analysis of the kilonova properties suggests a rapid evolution toward red colors, similar to AT2017gfo, and a low nIR luminosity, possibly due to the presence of a long-lived neutron star. The global properties of the environment, the inferred low mass (M_ej < 0.006 Msun) and velocities (v > 0.05 c) of lanthanide-rich ejecta are consistent with a binary neutron star merger progenitor. △ Less

Submitted 26 August, 2019; v1 submitted 3 May, 2019; originally announced May 2019.

Comments: 14 pages, 6 figures, MNRAS, in press. Moderate revision, added Figure 5 and X-ray data to Table 1

Synthetic Radio Imaging for Quiescent and CME-flare Scenarios

Authors: Sofia-Paraskevi Moschou, Igor Sokolov, Ofer Cohen, Jeremy J. Drake, Dmitry Borovikov, Justin C. Kasper, Julian D. Alvarado-Gomez, Cecilia Garraffo

Abstract: Radio observations grant access to a wide range of physical processes through different emission mechanisms. These processes range from thermal and quiescent to eruptive phenomena, such as shock waves and particle beams. We present a new synthetic radio imaging tool that calculates and visualizes the Bremsstrahlung radio emission. This tool works concurrently with state-of-the-art Magnetohydrodyna… ▽ More Radio observations grant access to a wide range of physical processes through different emission mechanisms. These processes range from thermal and quiescent to eruptive phenomena, such as shock waves and particle beams. We present a new synthetic radio imaging tool that calculates and visualizes the Bremsstrahlung radio emission. This tool works concurrently with state-of-the-art Magnetohydrodynamic (MHD) simulations of the solar corona using the code BATS-R-US. Our model produces results that are in good agreement with both high and low frequency observations of the solar disk. In this study, a ray-tracing algorithm is used and the radio intensity is computed along the actual curved ray trajectories. We illustrate the importance of refraction in locating the radio emitting source by comparison of the radio imaging illustrations when the line-of-sight instead of the refracted paths are considered. We are planning to incorporate non-thermal radio emission mechanisms in a future version of the radio imaging tool. △ Less

Submitted 25 September, 2018; originally announced September 2018.

Comments: 15 pages, 8 figures, Accepted by ApJ (NOTE: The current version includes extra material compared to the accepted to ApJ version.)

Exoplanet Modulation of Stellar Coronal Radio Emission

Authors: Ofer Cohen, Sofia-Paraskevi Moschou, Alex Glocer, Igor Sokolov, Tsevi Mazeh, Jeremy Drake, Cecilia Garraffo, Julian Alvarado-Gomez

Abstract: The search for exoplanets in the radio bands has been focused on detecting radio emissions produced by the interaction between magnetized planets and the stellar wind (auroral emission). Here we introduce a new tool, which is part of our MHD stellar corona model, to predict the ambient coronal radio emission and its modulations induced by a close planet. For simplicity, the present work assumes th… ▽ More The search for exoplanets in the radio bands has been focused on detecting radio emissions produced by the interaction between magnetized planets and the stellar wind (auroral emission). Here we introduce a new tool, which is part of our MHD stellar corona model, to predict the ambient coronal radio emission and its modulations induced by a close planet. For simplicity, the present work assumes that the exoplanet is stationary in the frame rotating with the stellar rotation. We explore the radio flux modulations using a limited parameter space of idealized cases by changing the magnitude of the planetary field, its polarity, the planetary orbital separation, and the strength of the stellar field. We find that the modulations induced by the planet could be significant and observable in the case of hot Jupiter planets --- above 100% modulation with respect to the ambient flux in the $10-100~MHz$ range in some cases, and 2-10% in the frequency bands above $250~MHz$ for some cases. Thus, our work indicates that radio signature of exoplanets might not be limited to low-frequency radio range. We find that the intensity modulations are sensitive to the planetary magnetic field polarity for short-orbit planets, and to the stellar magnetic field strength for all cases. The new radio tool, when applied to real systems, could provide predictions for the frequency range at which the modulations can be observed by current facilities. △ Less

Submitted 17 September, 2018; originally announced September 2018.

Comments: Accepted to AJ

Extended MHD modeling of the steady solar corona and the solar wind

Authors: Tamas I. Gombosi, Bart van der Holst, Ward B. Manchester, Igor V. Sokolov

Abstract: The history and present state of large-scale magnetohydrodynamic (MHD) modeling of the solar corona and the solar wind with steady or quasi-steady coronal physics is reviewed. We put the evolution of ideas leading to the recognition of the existence of an expanding solar atmosphere into historical context. The development and main features of the first generation of global corona and solar wind mo… ▽ More The history and present state of large-scale magnetohydrodynamic (MHD) modeling of the solar corona and the solar wind with steady or quasi-steady coronal physics is reviewed. We put the evolution of ideas leading to the recognition of the existence of an expanding solar atmosphere into historical context. The development and main features of the first generation of global corona and solar wind models are described in detail. This historical perspective is also applied to the present suite of global corona and solar wind models. We discuss the evolution of new ideas and their implementation into numerical simulation codes. We point out the scientific and computational challenges facing these models and discuss the ways various groups tried to overcome these challenges. Next, we discuss the latest, state-of-the art models and point to the expected next steps in modeling the corona and the interplanetary medium. △ Less

Submitted 1 July, 2018; originally announced July 2018.

Comments: Accepted for publication, Living Reviews in Solar Physics, 2018

Long-period SU UMa dwarf nova V1006 Cygni outburst activity and variability at different brightness states in 2015-2017

Authors: E. P. Pavlenko, S. Yu. Shugarov, A. O. Simon, A. A. Sosnovskij, K. A. Antonyuk, O. I. Antonyuk, A. V. Shchurova, A. V. Baklanov, Ju. V. Babina, A. S. Sklyanov, V. V. Vasylenko, V. G. Godunova, I. Sokolov, I. V. Rudakov

Abstract: CCD photometric observations of the dwarf nova V1006 Cyg were carried out in 2015-2017 with 11 telescopes located at 7 observatories. They covered the 2015 superoutburst with rebrightening, five normal outbursts of ~4-day duration and one wide outburst that lasted at least seven days. The interval between normal outbursts was 16 and 22 days, and between superoutbursts is expected to be longer than… ▽ More CCD photometric observations of the dwarf nova V1006 Cyg were carried out in 2015-2017 with 11 telescopes located at 7 observatories. They covered the 2015 superoutburst with rebrightening, five normal outbursts of ~4-day duration and one wide outburst that lasted at least seven days. The interval between normal outbursts was 16 and 22 days, and between superoutbursts is expected to be longer than 124 days. The positive superhumps with the mean period of 0^d.10544(10) and 0^d.10406(17) were detected during the 2015 superoutburst and during the short-term quiescence between rebrightening and the start of the first normal outburst, respectively. During a wide 2015 outburst the orbital period 0^d.09832(15) was found. The amplitude of this signal was ~2.5 times larger at the outburst decline than at its end. During the quiescence stage between the first and the second normal outbursts in 2017 we possibly detected the negative superhumps with the period of 0^d.09714(7). In all other cases of quiescence we found only the quasi-periodic brightness variations on a time scale of 20-30 minutes with a different degree of coherence and a variable amplitude reaching 0.5 mag in extremal cases. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: 17 pages, 10 figures, two tables, accepted for publication in COASP

The Excess Density of Field Galaxies near z=0.56 around the Gamma-Ray Burst GRB021004 Position

Authors: I. V. Sokolov, A. J. Castro-Tirado, O. P. Zhelenkova, I. A. Solovyev, O. V. Verkhodanov, V. V. Sokolov

Abstract: We test for reliability any signatures of field galaxies clustering in the GRB 021004 line of sight. The first signature is the GRB 021004 field photometric redshifts distribution based on the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences observations with a peak near z = 0.56 estimated from multicolor photometry in the GRB direction. The second signatur… ▽ More We test for reliability any signatures of field galaxies clustering in the GRB 021004 line of sight. The first signature is the GRB 021004 field photometric redshifts distribution based on the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences observations with a peak near z = 0.56 estimated from multicolor photometry in the GRB direction. The second signature is the Mg II 2796, 2803AA absorption doublet at z = 0.56 in VLT/UVES spectra obtained for the GRB 021004 afterglow. The third signature is the galaxy clustering in a larger (of about 3 sq.deg.) area around GRB 021004 with an effective peak near z = 0.56 for both the spectral and photometric redshifts from a few catalogs of clusters based on the Sloan Digital Sky Survey (SDSS) and Baryon Oscillation Spectroscopic Survey (BOSS) as a part of SDSS-III. From catalog data the size of the whole inhomogeneity in distribution of the galaxy cluster with the peak near z = 0.56 was also estimated as about 6--8 deg. or 140--190 Mpc. A possibility of inhomogeneity (a galaxy cluster) near the GRB 021004 direction can be also confirmed by an inhomogeneity in cosmic microwave background related with the Sunyaev-Zeldovich effect. △ Less

Submitted 3 December, 2018; v1 submitted 18 May, 2018; originally announced May 2018.

Comments: 23 pages, 13 figures

Journal ref: Astrophysical Bulletin, Volume 73, Issue 2, pp.111-123 (2018)

Chromosphere to 1 AU Simulation of the 2011 March 7th Event: A Comprehensive Study of Coronal Mass Ejection Propagation

Authors: M. Jin, W. B. Manchester, B. van der Holst, I. Sokolov, G. Toth, A. Vourlidas, C. A. de Koning, T. I. Gombosi

Abstract: We perform and analyze results of a global magnetohydrodyanmic (MHD) simulation of the fast coronal mass ejection (CME) that occurred on 2011 March 7. The simulation is made using the newly developed Alfvén Wave Solar Model (AWSoM), which describes the background solar wind starting from the upper chromosphere and extends to 24 R$_{\odot}$. Coupling AWSoM to an inner heliosphere (IH) model with th… ▽ More We perform and analyze results of a global magnetohydrodyanmic (MHD) simulation of the fast coronal mass ejection (CME) that occurred on 2011 March 7. The simulation is made using the newly developed Alfvén Wave Solar Model (AWSoM), which describes the background solar wind starting from the upper chromosphere and extends to 24 R$_{\odot}$. Coupling AWSoM to an inner heliosphere (IH) model with the Space Weather Modeling Framework (SWMF) extends the total domain beyond the orbit of Earth. Physical processes included in the model are multi-species thermodynamics, electron heat conduction (both collisional and collisionless formulations), optically thin radiative cooling, and Alfvén-wave turbulence that accelerates and heats the solar wind. The Alfvén-wave description is physically self-consistent, including non-Wentzel-Kramers-Brillouin (WKB) reflection and physics-based apportioning of turbulent dissipative heating to both electrons and protons. Within this model, we initiate the CME by using the Gibson-Low (GL) analytical flux rope model and follow its evolution for days, in which time it propagates beyond STEREO A. A detailed comparison study is performed using remote as well as \textit{in situ} observations. Although the flux rope structure is not compared directly due to lack of relevant ejecta observation at 1 AU in this event, our results show that the new model can reproduce many of the observed features near the Sun (e.g., CME-driven extreme ultraviolet (EUV) waves, deflection of the flux rope from the coronal hole, "double-front" in the white light images) and in the heliosphere (e.g., shock propagation direction, shock properties at STEREO A). △ Less

Submitted 27 November, 2016; originally announced November 2016.

Comments: 55 pages, 17 figures, 1 table, accepted by ApJ

Threaded-Field-Lines Model for the Low Solar Corona Powered by the Alfven Wave Turbulence

Authors: Igor V. Sokolov, Bart van der Holst, Ward B. Manchester, Doga Can Su Ozturk, Judit Szente, Aleksandre Taktakishvili, Gabor Tóth, Meng Jin, Tamas I. Gombosi

Abstract: We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfvén wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In clos… ▽ More We present an updated global model of the solar corona, including the transition region. We simulate the realistic tree-dimensional (3D) magnetic field using the data from the photospheric magnetic field measurements and assume the magnetohydrodynamic (MHD) Alfvén wave turbulence and its non-linear dissipation to be the only source for heating the coronal plasma and driving the solar wind. In closed field regions the dissipation efficiency in a balanced turbulence is enhanced. In the coronal holes we account for a reflection of the outward propagating waves, which is accompanied by generation of weaker counter-propagating waves. The non-linear cascade rate degrades in strongly imbalanced turbulence, thus resulting in colder coronal holes. The distinctive feature of the presented model is the description of the low corona as almost-steady-state low-beta plasma motion and heat flux transfer along the magnetic field lines. We trace the magnetic field lines through each grid point of the lower boundary of the global corona model, chosen at some heliocentric distance, $R=R_{b}\sim1.1\ R_\odot$ well above the transition region. One can readily solve the plasma parameters along the magnetic field line from 1D equations for the plasma motion and heat transport together with the Alfvén wave propagation, which adequately describe physics within the heliocentric distances range, $R_{\odot}<R<R_{b}$, in the low solar corona. By interfacing this threaded-field-lines model with the full MHD global corona model at $r=R_{b}$, we find the global solution and achieve a faster-than-real-time performance of the model on $\sim200$ cores. △ Less

Submitted 18 November, 2020; v1 submitted 8 September, 2016; originally announced September 2016.

Comments: 23 pages, 6 figures

Data Constrained Coronal Mass Ejections in A Global Magnetohydrodynamics Model

Authors: M. Jin, W. B. Manchester, B. van der Holst, I. Sokolov, G. Toth, R. E. Mullinix, A. Taktakishvili, A. Chulaki, T. I. Gombosi

Abstract: We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GO… ▽ More We present a first-principles-based coronal mass ejection (CME) model suitable for both scientific and operational purposes by combining a global magnetohydrodynamics (MHD) solar wind model with a flux rope-driven CME model. Realistic CME events are simulated self-consistently with high fidelity and forecasting capability by constraining initial flux rope parameters with observational data from GONG, SOHO/LASCO, and STEREO/COR. We automate this process so that minimum manual intervention is required in specifying the CME initial state. With the newly developed data-driven Eruptive Event Generator Gibson-Low (EEGGL), we present a method to derive Gibson-Low (GL) flux rope parameters through a handful of observational quantities so that the modeled CMEs can propagate with the desired CME speeds near the Sun. A test result with CMEs launched with different Carrington rotation magnetograms are shown. Our study shows a promising result for using the first-principles-based MHD global model as a forecasting tool, which is capable of predicting the CME direction of propagation, arrival time, and ICME magnetic field at 1 AU (see companion paper by Jin et al. 2016b). △ Less

Submitted 27 November, 2016; v1 submitted 17 May, 2016; originally announced May 2016.

Comments: 30 pages, 7 figures, 3 tables, accepted by ApJ

The optical identifcation of events with poorly defined locations: The case of the Fermi GBM GRB140801A

Authors: V. M. Lipunov, J. Gorosabel, M. V. Pruzhinskaya, A. de Ugarte Postigo, V. Pelassa, A. E. Tsvetkova, I. V. Sokolov, D. A. Kann, Dong Xu, E. S. Gorbovskoy, V. V. Krushinski, V. G. Kornilov, P. V. Balanutsa, S. V. Boronina, N. M. Budnev, Z. Cano, A. J. Castro-Tirado, V. V. Chazov, V. Connaughton, C. Delvaux, D. D. Frederiks, J. F. U. Fynbo, A. V. Gabovich, A. Goldstein, J. Greiner , et al. (28 additional authors not shown)

Abstract: We report the early discovery of the optical afterglow of gamma-ray burst (GRB) 140801A in the 137 deg$^2$ 3-$σ$ error-box of the Fermi Gamma-ray Burst Monitor (GBM). MASTER is the only observatory that automatically react to all Fermi alerts. GRB 140801A is one of the few GRBs whose optical counterpart was discovered solely from its GBM localization. The optical afterglow of GRB 140801A was found… ▽ More We report the early discovery of the optical afterglow of gamma-ray burst (GRB) 140801A in the 137 deg$^2$ 3-$σ$ error-box of the Fermi Gamma-ray Burst Monitor (GBM). MASTER is the only observatory that automatically react to all Fermi alerts. GRB 140801A is one of the few GRBs whose optical counterpart was discovered solely from its GBM localization. The optical afterglow of GRB 140801A was found by MASTER Global Robotic Net 53 sec after receiving the alert, making it the fastest optical detection of a GRB from a GBM error-box. Spectroscopy obtained with the 10.4-m Gran Telescopio Canarias and the 6-m BTA of SAO RAS reveals a redshift of $z=1.32$. We performed optical and near-infrared photometry of GRB 140801A using different telescopes with apertures ranging from 0.4-m to 10.4-m. GRB 140801A is a typical burst in many ways. The rest-frame bolometric isotropic energy release and peak energy of the burst is $E_\mathrm{iso} = 5.54_{-0.24}^{+0.26} \times 10^{52}$ erg and $E_\mathrm{p, rest}\simeq280$ keV, respectively, which is consistent with the Amati relation. The absence of a jet break in the optical light curve provides a lower limit on the half-opening angle of the jet $θ=6.1$ deg. The observed $E_\mathrm{peak}$ is consistent with the limit derived from the Ghirlanda relation. The joint Fermi GBM and Konus-Wind analysis shows that GRB 140801A could belong to the class of intermediate duration. The rapid detection of the optical counterpart of GRB 140801A is especially important regarding the upcoming experiments with large coordinate error-box areas. △ Less

Submitted 27 October, 2015; originally announced October 2015.

Comments: in press MNRAS, 2015

Journal ref: MNRAS, 455, 712 (2016)

A Steady-State Picture of Solar Wind Acceleration and Charge State Composition Derived from a Global Wave-Driven MHD Model

Authors: Rona Oran, Enrico Landi, Bart van der Holst, Susan T. Lepri, Alberto M. Vásquez, Federico. A. Nuevo, Richard Frazin, Ward B. Manchester IV, Igor V. Sokolov, Tamas I. Gombosi

Abstract: The higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSo… ▽ More The higher charge states found in slow ($<$400km s$^{-1}$) solar wind streams compared to fast streams have supported the hypothesis that the slow wind originates in closed coronal loops, and released intermittently through reconnection. Here we examine whether a highly ionized slow wind can also form along steady and open magnetic field lines. We model the steady-state solar atmosphere using AWSoM, a global magnetohydrodynamic model driven by Alfv{é}n waves, and apply an ionization code to calculate the charge state evolution along modeled open field lines. This constitutes the first charge states calculation covering all latitudes in a realistic magnetic field. The ratios $O^{+7}/O^{+6}$ and $C^{+6}/C^{+5}$ are compared to in-situ Ulysses observations, and are found to be higher in the slow wind, as observed; however, they are under-predicted in both wind types. The modeled ion fractions of S, Si, and Fe are used to calculate line-of-sight intensities, which are compared to EIS observations above a coronal hole. The agreement is partial, and suggests that all ionization rates are under-predicted. Assuming the presence of suprathermal electrons improved the agreement with both EIS and Ulysses observations; importantly, the trend of higher ionization in the slow wind was maintained. The results suggest there can be a sub-class of slow wind that is steady and highly ionized. Further analysis shows it originates from coronal hole boundaries (CHB), where the modeled electron density and temperature are higher than inside the hole, leading to faster ionization. This property of CHBs is global, and observationally supported by EUV tomography. △ Less

Submitted 29 December, 2014; originally announced December 2014.

Comments: Submitted to the Astrophysical Journal

Journal ref: ApJ, 806, 55, 2015

The dark nature of GRB 130528A and its host galaxy

Authors: S. Jeong, A. J. Castro-Tirado, M. Bremer, J. M. Winters, J. Gorosabel, S. Guziy, S. B. Pandey, M. Jelínek, R. Sánchez-Ramírez, Ilya V. Sokolov, N. V. Orekhova, A. S. Moskvitin, J. C. Tello, R. Cunniffe, O. Lara-Gil, S. R. Oates, D. Pérez-Ramírez, J. Bai, Y. Fan, C. Wang, I. H. Park

Abstract: We study the dark nature of GRB 130528A through multi-wavelength observations and conclude that the main reason for the optical darkness is local extinction inside of the host galaxy. Automatic observations were performed at BOOTES-4/MET robotic telescope. We also triggered target of opportunity (ToO) observation at the OSN, IRAM PdBI and the GTC+OSIRIS. The host galaxy photometric observations in… ▽ More We study the dark nature of GRB 130528A through multi-wavelength observations and conclude that the main reason for the optical darkness is local extinction inside of the host galaxy. Automatic observations were performed at BOOTES-4/MET robotic telescope. We also triggered target of opportunity (ToO) observation at the OSN, IRAM PdBI and the GTC+OSIRIS. The host galaxy photometric observations in optical to near-infrared (nIR) wavelengths were achieved through large ground-based aperture telescopes, such as the 10.4m GTC, the 4.2m WHT, 6m BTA, and the 2m LT. Based on these observations, spectral energy distributions (SED) for the host galaxy and afterglow were constructed. Thanks to mm observations at PdBI, we confirm the presence of a mm source within the XRT error circle that faded over the course of our observations and identify the host galaxy. However, we do not find any credible optical source within early observations with BOOTES-4/MET and 1.5m OSN telescopes. Spectroscopic observation of this galaxy by GTC showed a single faint emission line that likely corresponds to [OII] 3727\{AA} at a redshift of 1.250+/-0.001 implying a SFR(M_sun/yr) > 6.18 M_sun/yr without correcting for dust extinction. The probable extinction was revealed through analysis of the afterglow SED, resulting in a value of AV >= ~ 0.9 at the rest frame, this is comparable to extinction levels found among other dark GRBs. The SED of the host galaxy is explained well (chi2/d.o.f.=0.564) by a luminous (MB=-21.16), low-extinction (AV =0, rest frame), and aged (2.6 Gyr) stellar population. We can explain this apparent contradiction in global and line-of-sight extinction if the GRB birth place happened to lie in a local dense environment. In light of having relatively small specific SFR (SSFR) ~ 5.3 M_sun/yr (L/L_star)-1, this also could explain the age of the old stellar population of host galaxy. △ Less

Submitted 2 September, 2014; v1 submitted 3 April, 2014; originally announced April 2014.

Comments: 6 pages, 4 figures

Journal ref: A&A 569, A93 (2014)

Alfvén Wave Turbulence as a Coronal Heating Mechanism: Simultaneously Predicting the Heating Rate and the Wave-Induced Emission Line Broadening

Authors: R. Oran, E. Landi, B. van der Holst, I. V. Sokolov, T. I. Gombosi

Abstract: In the present work, we test the predictions of the AWSoM model, a global extended-MHD model capable of calculating the propagation and turbulent dissipation of Alfvén waves in any magnetic topology, against high resolution spectra of the quiescent off-disk solar corona. Wave dissipation is the only heating mechanism assumed in this model. Combining 3D model results with the CHIANTI atomic databas… ▽ More In the present work, we test the predictions of the AWSoM model, a global extended-MHD model capable of calculating the propagation and turbulent dissipation of Alfvén waves in any magnetic topology, against high resolution spectra of the quiescent off-disk solar corona. Wave dissipation is the only heating mechanism assumed in this model. Combining 3D model results with the CHIANTI atomic database, we were able to create synthetic line-of-sight spectra which include the effects of emission line broadening due to both thermal and wave-related non-thermal motions. To the best of our knowledge this is the first time a global model is used to obtain synthetic non-thermal line broadening. We obtained a steady-state solution driven by a synoptic magnetogram and compared the synthetic spectra with SUMER observations of a quiescent area above the solar west limb extending between 1.04 and 1.34 solar radii at the equator. Both the predicted line widths and the total line fluxes were consistent with the observations for 5 different ions. Using the 3D solution, we were able to locate the region that contributes the most to the emission used for measuring electron properties; we found that region to be a pseudo-streamer, whose modeled electron temperature and density are consistent with the measured ones. We conclude that the turbulent dissipation assumed in the AWSoM model can simultaneously account for the observed heating rate and the non-dissipated wave energy observed in this region. △ Less

Submitted 2 January, 2014; originally announced January 2014.

Comments: Submitted to the Astrophysical Journal

Alfven Wave Solar Model (AWSoM): Coronal Heating

Authors: Bart van der Holst, Igor V. Sokolov, Xing Meng, Meng Jin, Ward B. Manchester IV, Gabor Toth, Tamas I. Gombosi

Abstract: We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimen… ▽ More We present a new version of the Alfven Wave Solar Model (AWSoM), a global model from the upper chromosphere to the corona and the heliosphere. The coronal heating and solar wind acceleration are addressed with low-frequency Alfven wave turbulence. The injection of Alfven wave energy at the inner boundary is such that the Poynting flux is proportional to the magnetic field strength. The three-dimensional magnetic field topology is simulated using data from photospheric magnetic field measurements. This model does not impose open-closed magnetic field boundaries; those develop self-consistently. The physics includes: (1) The model employs three different temperatures, namely the isotropic electron temperature and the parallel and perpendicular ion temperatures. The firehose, mirror, and ion-cyclotron instabilities due to the developing ion temperature anisotropy are accounted for. (2) The Alfven waves are partially reflected by the Alfven speed gradient and the vorticity along the field lines. The resulting counter-propagating waves are responsible for the nonlinear turbulent cascade. The balanced turbulence due to uncorrelated waves near the apex of the closed field lines and the resulting elevated temperatures are addressed. (3) To apportion the wave dissipation to the three temperatures, we employ the results of the theories of linear wave damping and nonlinear stochastic heating. (4) We have incorporated the collisional and collisionless electron heat conduction. We compare the simulated multi-wavelength EUV images of CR2107 with the observations from STEREO/EUVI and SDO/AIA instruments. We demonstrate that the reflection due to strong magnetic fields in proximity of active regions intensifies the dissipation and observable emission sufficiently. △ Less

Submitted 2 December, 2013; v1 submitted 16 November, 2013; originally announced November 2013.

Comments: 26 pages, 5 figures; submitted to the Astrophysical Journal

A Global Wave-Driven MHD Solar Model with a Unified Treatment of Open and Closed Magnetic Field Topologies

Authors: R. Oran, B. van der Holst, E. Landi, M. Jin, I. V. Sokolov, T. I. Gombosi

Abstract: We describe, analyze and validate the recently developed Alfvén Wave Solar Model (AWSoM), a 3D global model starting from the top of the chromosphere and extending into interplanetary space (up to 1-2 AU). This model solves the extended two temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low frequency Alfvén waves. In this picture, heating and acceleration of the… ▽ More We describe, analyze and validate the recently developed Alfvén Wave Solar Model (AWSoM), a 3D global model starting from the top of the chromosphere and extending into interplanetary space (up to 1-2 AU). This model solves the extended two temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low frequency Alfvén waves. In this picture, heating and acceleration of the plasma are due to wave dissipation and wave pressure gradients, respectively. The dissipation process is described by a fully developed turbulent cascade of counter-propagating waves. We adopt a unified approach for calculating the wave dissipation in both open and closed magnetic field lines, allowing for a self-consistent treatment of any magnetic topology. Wave dissipation is the only heating mechanism assumed in the model, and no geometric heating functions are invoked. Electron heat conduction and radiative cooling are also included. We demonstrate that the large-scale, steady-state (in the co-rotating frame) properties of the solar environment are reproduced, using three adjustable parameters: the Poynting flux of chromospheric Alfvén waves, the perpendicular correlation length of the turbulence, and a pseudo-reflection coefficient. We compare model results for Carrington Rotation 2063 (November-December 2007) to remote observations in the EUV and X-ray ranges from STEREO, SOHO and Hinode spacecraft, as well as to in-situ measurements performed by Ulysses. The model results are in good agreement with observations. This is the first global model capable of simultaneously reproducing the multi-wavelength observations of the lower corona and the wind structure beyond Earth's orbit. △ Less

Submitted 17 July, 2013; originally announced July 2013.

Journal ref: Astrophysical Journal, 2013, 778, 176

MHD Waves and Coronal Heating: Unifying Empirical and MHD Turbulence Models

Authors: Igor V. Sokolov, Bart van der Holst, Rona Oran, Cooper Downs, Ilia I. Roussev, Meng Jin, Ward B. Manchester IV, Rebekah M. Evans, Tamas I. Gombosi

Abstract: We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only… ▽ More We present a new global model of the solar corona, including the low corona, the transition region and the top of chromosphere. The realistic 3D magnetic field is simulated using the data from the photospheric magnetic field measurements. The distinctive feature of the new model is incorporating the MHD Alfven wave turbulence. We assume this turbulence and its non-linear dissipation to be the only momentum and energy source for heating the coronal plasma and driving the solar wind. The difference between the turbulence dissipation efficiency in coronal holes and that in closed field regions is because the non-linear cascade rate degrades in strongly anisotropic turbulence in coronal holes (no inward propagating wave), thus resulting in colder coronal holes with the bi-modal solar wind originating from them. The detailed presentation of the theoretical model is illustrated with the synthetic images for multi-wavelength EUV emission compared with the observations from SDO AIA and Stereo EUVI instruments for the Carrington rotation 2107. △ Less

Submitted 5 December, 2012; v1 submitted 15 August, 2012; originally announced August 2012.

Comments: Submitted to the Astrophysical Journal

Simulating the long-term evolution of radiative shocks in shock tubes

Authors: B. van der Holst, G. Toth, I. V. Sokolov, B. R. Torralva, K. G. Powell, R. P. Drake

Abstract: We present the latest improvements in the Center for Radiative Shock Hydrodynamics (CRASH) code, a parallel block-adaptive-mesh Eulerian code for simulating high-energy-density plasmas. The implementation can solve for radiation models with either a gray or a multigroup method in the flux-limited-diffusion approximation. The electrons and ions are allowed to be out of temperature equilibrium and f… ▽ More We present the latest improvements in the Center for Radiative Shock Hydrodynamics (CRASH) code, a parallel block-adaptive-mesh Eulerian code for simulating high-energy-density plasmas. The implementation can solve for radiation models with either a gray or a multigroup method in the flux-limited-diffusion approximation. The electrons and ions are allowed to be out of temperature equilibrium and flux-limited electron thermal heat conduction is included. We have recently implemented a CRASH laser package with 3-D ray tracing, resulting in improved energy deposition evaluation. New, more accurate opacity models are available which significantly improve radiation transport in materials like xenon. In addition, the HYPRE preconditioner has been added to improve the radiation implicit solver. With this updated version of the CRASH code we study radiative shock tube problems. In our set-up, a 1 ns, 3.8 kJ laser pulse irradiates a 20 micron beryllium disk, driving a shock into a xenon-filled plastic tube. The electrons emit radiation behind the shock. This radiation from the shocked xenon preheats the unshocked xenon. Photons traveling ahead of the shock will also interact with the plastic tube, heat it, and in turn this can drive another shock off the wall into the xenon. We are now able to simulate the long term evolution of radiative shocks. △ Less

Submitted 6 June, 2012; originally announced June 2012.

Comments: submitted to High Energy Density Physics

Simulating radiative shocks in nozzle shock tubes

Authors: B. van der Holst, G. Toth, I. V. Sokolov, L. K. S. Daldorff, K. G. Powell, R. P. Drake

Abstract: We use the recently developed Center for Radiative Shock Hydrodynamics (CRASH) code to numerically simulate laser-driven radiative shock experiments. These shocks are launched by an ablated beryllium disk and are driven down xenon-filled plastic tubes. The simulations are initialized by the two-dimensional version of the Lagrangian Hyades code which is used to evaluate the laser energy deposition… ▽ More We use the recently developed Center for Radiative Shock Hydrodynamics (CRASH) code to numerically simulate laser-driven radiative shock experiments. These shocks are launched by an ablated beryllium disk and are driven down xenon-filled plastic tubes. The simulations are initialized by the two-dimensional version of the Lagrangian Hyades code which is used to evaluate the laser energy deposition during the first 1.1ns. The later times are calculated with the CRASH code. This code solves for the multi-material hydrodynamics with separate electron and ion temperatures on an Eulerian block-adaptive-mesh and includes a multi-group flux-limited radiation diffusion and electron thermal heat conduction. The goal of the present paper is to demonstrate the capability to simulate radiative shocks of essentially three-dimensional experimental configurations, such as circular and elliptical nozzles. We show that the compound shock structure of the primary and wall shock is captured and verify that the shock properties are consistent with order-of-magnitude estimates. The produced synthetic radiographs can be used for comparison with future nozzle experiments at high-energy-density laser facilities. △ Less

Submitted 18 September, 2011; originally announced September 2011.

Comments: submitted to High Energy Density Physics

The Dynamics of Stellar Coronae Harboring Hot-jupiters II. A Space Weather Event on A Hot-jupiter

Authors: O. Cohen, V. L. Kashyap, J. J. Drake, I. V. Sokolov, T. I. Gombosi

Abstract: We carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations of space weather events on Earth. The simulation includes the planetary orbital motion, which leads to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to t… ▽ More We carry out a numerical simulation depicting the effects of a Coronal Mass Ejection (CME) event on a close-in giant planet in an extrasolar system. We drive the CME in a similar manner as in simulations of space weather events on Earth. The simulation includes the planetary orbital motion, which leads to the forming of a comet-like planetary magnetotail which is oriented almost perpendicular to the direction of propagation of the CME. The combination of this feature and the fact that the CME does not expand much by the time it reaches the planet leads to a unique CME-magnetosphere interaction, where the CME itself is highly affected by the presence of the planetary magnetosphere. We find that the planet is well-shielded from CME penetration, even for a relatively weak internal magnetic field. The planetary angular momentum loss associated with such an event is negligible compared to the total planetary angular momentum. We also find that the energy which is deposited in the magnetosphere is much higher than in the case of the Earth, and our simulation suggests there is a large-scale change in the orientation of the magnetosphere-ionosphere current system during the CME event. △ Less

Submitted 7 July, 2011; v1 submitted 20 February, 2011; originally announced February 2011.

Comments: 16 pages, 10 figures, accepted to ApJ

The Dynamics of Stellar Coronae Harboring Hot-jupiters I. A Time-dependent MHD Simulation of the Interplanetary Environment in the HD 189733 Planetary System

Authors: O. Cohen, V. L. Kashyap, J. J. Drake, I. V. Sokolov, C. Garraffo, T. I. Gombosi

Abstract: We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sect… ▽ More We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar planetary system to study the interaction of the stellar magnetic field and wind with the planetary magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized by sectors with different plasma properties. The star-planet interaction varies in magnitude and complexity, depending on the planetary phase, planetary magnetic field strength, and the relative orientation of the stellar and planetary fields. It also reveals a long, comet-like tail which is a result of the wrapping of the planetary magnetospheric tail by its fast orbital motion. A reconnection event occurs at a specific orbital phase, causing mass loss from the planetary magnetosphere that can generate a hot spot on the stellar surface. The simulation also shows that the system has sufficient energy to produce hot-spots observed in Ca II lines in giant planet hosting stars. However, the short duration of the reconnection event suggests that such SPI cannot be observed persistently. △ Less

Submitted 25 January, 2011; originally announced January 2011.

Comments: Submitted to ApJ, 14 pages, 9 figures

Crash: A Block-Adaptive-Mesh Code for Radiative Shock Hydrodynamics - Implementation and Verification

Authors: B. van der Holst, G. Toth, I. V. Sokolov, K. G. Powell, J. P. Holloway, E. S. Myra, Q. Stout, M. L. Adams, J. E. Morel, R. P. Drake

Abstract: We describe the CRASH (Center for Radiative Shock Hydrodynamics) code, a block adaptive mesh code for multi-material radiation hydrodynamics. The implementation solves the radiation diffusion model with the gray or multigroup method and uses a flux limited diffusion approximation to recover the free-streaming limit. The electrons and ions are allowed to have different temperatures and we include a… ▽ More We describe the CRASH (Center for Radiative Shock Hydrodynamics) code, a block adaptive mesh code for multi-material radiation hydrodynamics. The implementation solves the radiation diffusion model with the gray or multigroup method and uses a flux limited diffusion approximation to recover the free-streaming limit. The electrons and ions are allowed to have different temperatures and we include a flux limited electron heat conduction. The radiation hydrodynamic equations are solved in the Eulerian frame by means of a conservative finite volume discretization in either one, two, or three-dimensional slab geometry or in two-dimensional cylindrical symmetry. An operator split method is used to solve these equations in three substeps: (1) solve the hydrodynamic equations with shock-capturing schemes, (2) a linear advection of the radiation in frequency-logarithm space, and (3) an implicit solve of the stiff radiation diffusion, heat conduction, and energy exchange. We present a suite of verification test problems to demonstrate the accuracy and performance of the algorithms. The CRASH code is an extension of the Block-Adaptive Tree Solarwind Roe Upwind Scheme (BATS-R-US) code with this new radiation transfer and heat conduction library and equation-of-state and multigroup opacity solvers. Both CRASH and BATS-R-US are part of the publicly available Space Weather Modeling Framework (SWMF). △ Less

Submitted 19 January, 2011; originally announced January 2011.

Comments: 51 pages, 19 figures; submitted to Astrophysical Journal

Study of Faint Galaxies in the Field of GRB 021004

Authors: Yu. V. Baryshev, I. V. Sokolov, A. S. Moskvitin, T. A. Fatkhullin, N. V. Nabokov, Brajesh Kumar

Abstract: We present an analysis of BVRcIc observations of the field sized around 4' x 4' centered at the host galaxy of the gamma-ray burst GRB 021004 with the 6-m BTA telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. We measured the magnitudes and constructed the color diagrams for 311 galaxies detected in the field (S/N > 3). The differential and integral counts of ga… ▽ More We present an analysis of BVRcIc observations of the field sized around 4' x 4' centered at the host galaxy of the gamma-ray burst GRB 021004 with the 6-m BTA telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences. We measured the magnitudes and constructed the color diagrams for 311 galaxies detected in the field (S/N > 3). The differential and integral counts of galaxies up to the limit, corresponding to 28.5 (B), 28.0 (V), 27.0 (Rc), 26.5 (Ic) were computed. We compiled the galaxy catalog, consisting of 183 objects, for which the photometric redshifts up to the limiting magnitudes 26.0 (B), 25.5 (V), 25.0 (Rc), 24.5 (Ic) were determined using the HyperZ code. We then examined the radial distribution of galaxies based on the z estimates. We have built the curves expected in the case of a uniform distribution of galaxies in space, and obtained the estimates for the size and contrast of the possible super-large-scale structures, which are accessible with the observations of this type. △ Less

Submitted 26 October, 2010; v1 submitted 19 October, 2010; originally announced October 2010.

Comments: 15 pages, 7 figures, 2 tables, to be published in Astrophysical Bulletin 65, 311-325 (2010)

The Impact of Hot Jupiters on the Spin-down of their Host Stars

Authors: O. Cohen, J. J Drake, V. L. Kashyap, I. V. Sokolov, T. I. Gombosi

Abstract: We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass loss rate drops by a factor of $\approx$1.5-2, while the angular momentum loss rate drops by a factor of $\approx$4 as the distance decreases past the Alfvén surface. This reduction in angular momentum los… ▽ More We present a numerical Magnetohydrodynamic (MHD) study of the dependence of stellar mass and angular momentum- loss rates on the orbital distance to close-in giant planets. We find that the mass loss rate drops by a factor of $\approx$1.5-2, while the angular momentum loss rate drops by a factor of $\approx$4 as the distance decreases past the Alfvén surface. This reduction in angular momentum loss is due to the interaction between the stellar and planetary Alfvén surfaces, which modifies the global structure of the stellar corona and stellar wind on the hemisphere facing the planet, as well as the opposite hemisphere. The simulation also shows that the magnitude of change in angular momentum loss rate depends mostly on the strength of the planetary magnetic field and not on its polarity. The interaction however, begins at greater separation if the overall field topology of the star and the planet are of anti-aligned. Our results are consistent with evidence for excess angular momentum in stars harboring close-in giant planets, and show that the reduction in wind-driven angular momentum loss can compete with, and perhaps dominate, spin-up due to tidal interaction. △ Less

Submitted 29 September, 2010; originally announced September 2010.

Comments: 6 pages, 3 figures

Algorithm for Tracing Radio Rays in Solar Corona and Chromosphere

Authors: Leonid Benkevitch, Igor Sokolov, Divya Oberoi, Thomas Zurbuchen

Abstract: In this paper a new efficient algorithm for computation of radio wave ray trajectories is described. The algorithm is based on an original second-order difference scheme with a specific "length-conservation" property, which allows to resolve the ray shape even in the regions where its curvature is high. Besides the scheme, the algorithm includes a number of mechanisms ensuring its correctness and… ▽ More In this paper a new efficient algorithm for computation of radio wave ray trajectories is described. The algorithm is based on an original second-order difference scheme with a specific "length-conservation" property, which allows to resolve the ray shape even in the regions where its curvature is high. Besides the scheme, the algorithm includes a number of mechanisms ensuring its correctness and stability. The algorithm is intended for obtaining multi-pixel solar images in wide range of radio frequencies, and it is designed to be used in studies of the solar environment with modern high-resolution radio interferometers and radio telescopes such as the Murchison Widefield Array △ Less

Submitted 18 October, 2010; v1 submitted 29 June, 2010; originally announced June 2010.

Comments: 14 pages, 9 figures; author Thomas Zurbuchen added; equivalence of the ray equation to that of Walker-Haselgrove shown; derivation of the ray equation from Fermat's principle shown in detail

Towards a Realistic, Data-Driven Thermodynamic MHD Model of the Global Solar Corona

Authors: Cooper Downs, Ilia I. Roussev, Bart van der Holst, Noé Lugaz, Igor V. Sokolov, Tamas I. Gombosi

Abstract: In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnet… ▽ More In this work we describe our implementation of a thermodynamic energy equation into the global corona model of the Space Weather Modeling Framework (SWMF), and its development into the new Lower Corona (LC) model. This work includes the integration of the additional energy transport terms of coronal heating, electron heat conduction, and optically thin radiative cooling into the governing magnetohydrodynamic (MHD) energy equation. We examine two different boundary conditions using this model; one set in the upper transition region (the Radiative Energy Balance model), as well as a uniform chromospheric condition where the transition region can be modeled in its entirety. Via observation synthesis from model results and the subsequent comparison to full sun extreme ultraviolet (EUV) and soft X-Ray observations of Carrington Rotation (CR) 1913 centered on Aug 27, 1996, we demonstrate the need for these additional considerations when using global MHD models to describe the unique conditions in the low corona. Through multiple simulations we examine ability of the LC model to asses and discriminate between coronal heating models, and find that a relative simple empirical heating model is adequate in reproducing structures observed in the low corona. We show that the interplay between coronal heating and electron heat conduction provides significant feedback onto the 3D magnetic topology in the low corona as compared to a potential field extrapolation, and that this feedback is largely dependent on the amount of mechanical energy introduced into the corona. △ Less

Submitted 11 December, 2009; originally announced December 2009.

Comments: 17 pages, 11 figures, Submitted to ApJ on 12/08/2009

Journal ref: Astrophys.J.712:1219-1231,2010

Solar-Terrestrial Simulations of CMEs with a Realistic Initiation Mechanism: Case Study for Active Region 10069

Authors: N. Lugaz, I. I. Roussev, I. V. Sokolov, C. Jacobs

Abstract: Most simulations of coronal mass ejections (CMEs) to date either focus on the interplanetary propagation of a giant plasma "blob" without paying too much attention to its origin and to the formation process or they focus on the complex evolution of the coronal magnetic field due to (sub-)photospheric motions which result in an eruption. Here, we present global simulations of CMEs where coronal m… ▽ More Most simulations of coronal mass ejections (CMEs) to date either focus on the interplanetary propagation of a giant plasma "blob" without paying too much attention to its origin and to the formation process or they focus on the complex evolution of the coronal magnetic field due to (sub-)photospheric motions which result in an eruption. Here, we present global simulations of CMEs where coronal motions are used to produce a realistic evolution of the coronal magnetic field and cause an eruption. We focus on active region 10069, which produced a number of eruptions in late August 2002, including the August 24, 2002 CME - a fast (~2000 km/s) eruption originating from W81-, as well as a slower eruption on August 22, 2002 (originating from W62). Using a three-dimensional magneto-hydrodynamic (MHD) simulation of these ejections with the Space Weather Modeling Framework (SWMF), we show how a realistic initiation mechanism enables us to study the deflection of the CME in the corona and in the heliosphere. Reconnection of the erupting magnetic field with that of neighboring streamers and active regions modify the solar connectivity of the field lines connecting to Earth and change the expected solar energetic particle fluxes. Comparing the results at 1 AU of our simulations with in situ observations by the ACE spacecraft, we propose an alternate solar origin for the shock wave observed at L1 on August 26. △ Less

Submitted 10 November, 2009; originally announced November 2009.

Comments: 4 pages, 2 figures, refereed proceedings for Solar Wind 12

Journal ref: AIP Conf.Proc.1216:440-443,2010

Interactions of the magnetospheres of stars and close-in giant planets

Authors: O. Cohen, J. J. Drake, V. L. Kashyap, S. H. Saar, I. V. Sokolov, W. B. Manchester IV, K. C. Hansen, T. I. Gombosi

Abstract: Since the first discovery of an extrasolar planetary system more than a decade ago, hundreds more have been discovered. Surprisingly, many of these systems harbor Jupiter-class gas giants located close to the central star, at distances of 0.1 AU or less. Observations of chromospheric 'hot spots' that rotate in phase with the planetary orbit, and elevated stellar X-ray luminosities,suggest that t… ▽ More Since the first discovery of an extrasolar planetary system more than a decade ago, hundreds more have been discovered. Surprisingly, many of these systems harbor Jupiter-class gas giants located close to the central star, at distances of 0.1 AU or less. Observations of chromospheric 'hot spots' that rotate in phase with the planetary orbit, and elevated stellar X-ray luminosities,suggest that these close-in planets significantly affect the structure of the outer atmosphere of the star through interactions between the stellar magnetic field and the planetary magnetosphere. Here we carry out the first detailed three-dimensional MagnetoHydroHynamics (MHD) simulation containing the two magnetic bodies and explore the consequences of such interactions on the steady-state coronal structure. The simulations reproduce the observable features of 1) increase in the total X-ray luminosity, 2) appearance of coronal hot spots, and 3) phase shift of these spots with respect to the direction of the planet. The proximate cause of these is an increase in the density of coronal plasma in the direction of the planet, which prevents the corona from expanding and leaking away this plasma via a stellar wind. The simulations produce significant low temperature heating. By including dynamical effects, such as the planetary orbital motion, the simulation should better reproduce the observed coronal heating. △ Less

Submitted 16 September, 2009; originally announced September 2009.

Journal ref: Astrophys.J.704:L85-L88,2009

The August 24, 2002 Coronal Mass Ejection: When a Western Limb Event Connects to Earth

Authors: N. Lugaz, I. I. Roussev, I. V. Sokolov

Abstract: We discuss how some coronal mass ejections (CMEs) originating from the western limb of the Sun are associated with space weather effects such as solar energetic particles (SEPs), shock or geo-effective ejecta at Earth. We focus on the August 24, 2002 coronal mass ejection, a fast (~ 2000 km/s) eruption originating from W81. Using a three-dimensional magneto-hydrodynamic simulation of this ejecti… ▽ More We discuss how some coronal mass ejections (CMEs) originating from the western limb of the Sun are associated with space weather effects such as solar energetic particles (SEPs), shock or geo-effective ejecta at Earth. We focus on the August 24, 2002 coronal mass ejection, a fast (~ 2000 km/s) eruption originating from W81. Using a three-dimensional magneto-hydrodynamic simulation of this ejection with the Space Weather Modeling Framework (SWMF), we show how a realistic initiation mechanism enables us to study the deflection of the CME in the corona and the heliosphere. Reconnection of the erupting magnetic field with that of neighboring streamers and active regions modify the solar connectivity of the field lines connecting to Earth and can also partly explain the deflection of the eruption during the first tens of minutes. Comparing the results at 1 AU of our simulation with observations by the ACE spacecraft, we find that the simulated shock does not reach Earth, but has a maximum angular span of about 120$^\circ$, and reaches 35$^\circ$ West of Earth in 58 hours. We find no significant deflection of the CME and its associated shock wave in the heliosphere, and we discuss the consequences for the shock angular span. △ Less

Submitted 7 October, 2008; originally announced October 2008.

Comments: 7 pages, 4 figures, IAU 257 Symposium Proceedings

Journal ref: IAU Symposium (2009), 257, 391-398

Three-Dimensional MHD Simulation of the 2003 October 28 Coronal Mass Ejection: Comparison with LASCO Coronagraph Observations

Authors: Ward Manchester IV, Angelos Vourlidas, Gabor Toth, Noe Lugaz, Ilia Roussev, Tamas Gombosi, Igor Sokolov, Darren De Zeeuw, Merav Opher

Abstract: We numerically model the coronal mass ejection (CME) event of October 28, 2003 that erupted from active region 10486 and propagated to Earth in less than 20 hours causing severe geomagnetic storms. The magnetohydrodynamic (MHD) model is formulated by first arriving at a steady state corona and solar wind employing synoptic magnetograms. We initiate two CMEs from the same active region, one appro… ▽ More We numerically model the coronal mass ejection (CME) event of October 28, 2003 that erupted from active region 10486 and propagated to Earth in less than 20 hours causing severe geomagnetic storms. The magnetohydrodynamic (MHD) model is formulated by first arriving at a steady state corona and solar wind employing synoptic magnetograms. We initiate two CMEs from the same active region, one approximately a day earlier that preconditions the solar wind for the much faster CME on the 28th. This second CME travels through the corona at a rate of over 2500 km s$^{-1}$ driving a strong forward shock. We clearly identify this shock in an image produced by the Large Angle Spectrometric Coronagraph (LASCO) C3, and reproduce the shock and its appearance in synthetic white light images from the simulation. We find excellent agreement with both the general morphology and the quantitative brightness of the model CME with LASCO observations. These results demonstrate that the CME shape is largely determined by its interaction with the ambient solar wind and may not be sensitive to the initiation process. We then show how the CME would appear as observed by wide-angle coronagraphs onboard the Solar Terrestrial Relations Observatory (STEREO) spacecraft. We find complex time evolution of the white-light images as a result of the way in which the density structures pass through the Thomson sphere. The simulation is performed with the Space Weather Modeling Framework (SWMF). △ Less

Submitted 23 May, 2008; originally announced May 2008.

Comments: 14 pages, 10 figures

Whether can the decision of the orthopositronium problem to stimulate studying the problem of a dark matter in the Universe?

Authors: B. M. Levin, V. I. Sokolov

Abstract: The nature of a dark matter in the Universe can be connected with realization of the (plus/minus) Planck mass in final state of the beta-plus-decay of the nuclei Na-22, Ga-68, etc. and in a gravitational field of sufficient force. This hypothesis (the topological quantum transition - a generalized "displacement current") follows out of the phenomenological model that has explained the anomalies… ▽ More The nature of a dark matter in the Universe can be connected with realization of the (plus/minus) Planck mass in final state of the beta-plus-decay of the nuclei Na-22, Ga-68, etc. and in a gravitational field of sufficient force. This hypothesis (the topological quantum transition - a generalized "displacement current") follows out of the phenomenological model that has explained the anomalies of the orthopositronium, forming in substance by the beta-decay positrons. Vacuumlike state of substance of the positive mass (Gliner, 1965) is structured in this model (the "microstructure" of vacuum) as the long-range "atom" with the Planck mass in the limited macroscopic "volume" of space-time. Accordingly, the negative mass C-field (Hoyle & Narlikar, 1964) plays the role of a compensating field. △ Less

Submitted 3 October, 2006; originally announced October 2006.

Comments: 8 pages

Density Perturbation in the Universe with Nontrivial Topology of Space-Time

Authors: S. S. Moskaliuk, A. V. Nesteruk, I. I. Sokolov

Abstract: A space spectrum of density perturbation In the Universe with nontrivial topology of space-time is shown to become discrete. A space spectrum of density perturbation In the Universe with nontrivial topology of space-time is shown to become discrete. △ Less

Submitted 12 February, 2006; originally announced February 2006.

Comments: 9 pages, Latex, no figures

Journal ref: Hadronic J. 18 (1995) 392-401

Magnetic Effects Change Our View of the Heliosheath

Authors: M. Opher, P. C. Liewer, M. Velli, T. I. Gombosi, W. Manchester, D. L. DeZeeuw, G. Toth, I. Sokolov

Abstract: There is currently a controversy as to whether Voyager 1 has already crossed the Termination Shock, the first boundary of the Heliosphere. The region between the Termination Shock and the Heliopause, the Helisheath, is one of the most unknown regions theoretically. In the Heliosheath magnetic effects are crucial, as the solar magnetic field is compressed at the Termination Shock by the slowing f… ▽ More There is currently a controversy as to whether Voyager 1 has already crossed the Termination Shock, the first boundary of the Heliosphere. The region between the Termination Shock and the Heliopause, the Helisheath, is one of the most unknown regions theoretically. In the Heliosheath magnetic effects are crucial, as the solar magnetic field is compressed at the Termination Shock by the slowing flow. Recently, our simulations showed that the Heliosheath presents remarkable dynamics, with turbulent flows and the presence of a jet flow at the current sheet that is unstable due to magnetohydrodynamic instabilities \cite{opher,opher1}. In this paper we review these recent results, and present an additional simulation with constant neutral atom background. In this case the jet is still present but with reduced intensity. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how the Heliosheath behaves. Already we can say that this region presents remarkable dynamics, with turbulent flows, indicating that the Heliosheath might be very different from what we previously thought. △ Less

Submitted 7 June, 2004; originally announced June 2004.

Comments: 6 pages, 5 figures, to appear in IGPP 3rd Annual International Astrophysics Conference, "PHYSICS OF THE OUTER HELIOSPHERE"

Magnetic Effects at the Edge of the Solar System: MHD Instabilities, the de Laval nozzle Effect and an Extended Jet

Authors: M. Opher, P. C. Liewer, M. Velli, L. Bettarini, T. I. Gombosi, W. Manchester, D. L. DeZeeuw, G. Toth, I. Sokolov

Abstract: To model the interaction between the solar wind and the interstellar wind, magnetic fields must be included. Recently Opher et al. 2003 found that, by including the solar magnetic field in a 3D high resolution simulation using the University of Michigan BATS-R-US code, a jet-sheet structure forms beyond the solar wind Termination Shock. Here we present an even higher resolution three-dimensional… ▽ More To model the interaction between the solar wind and the interstellar wind, magnetic fields must be included. Recently Opher et al. 2003 found that, by including the solar magnetic field in a 3D high resolution simulation using the University of Michigan BATS-R-US code, a jet-sheet structure forms beyond the solar wind Termination Shock. Here we present an even higher resolution three-dimensional case where the jet extends for $150AU$ beyond the Termination Shock. We discuss the formation of the jet due to a de Laval nozzle effect and it's su bsequent large period oscillation due to magnetohydrodynamic instabilities. To verify the source of the instability, we also perform a simplified two dimensional-geometry magnetohydrodynamic calculation of a plane fluid jet embedded in a neutral sheet with the profiles taken from our 3D simulation. We find remarkable agreement with the full three-dimensional evolution. We compare both simulations and the temporal evolution of the jet showing that the sinuous mode is the dominant mode that develops into a velocity-shear-instability with a growth rate of $5 \times 10^{-9} sec^{-1}=0.027 years^{-1}$. As a result, the outer edge of the heliosphere presents remarkable dynamics, such as turbulent flows caused by the motion of the jet. Further study, e.g., including neutrals and the tilt of the solar rotation from the magnetic axis, is required before we can definitively address how this outer boundary behaves. Already, however, we can say that the magnetic field effects are a major player in this region changing our previous notion of how the solar system ends. △ Less

Submitted 7 June, 2004; originally announced June 2004.

Comments: 24 pages, 13 figures, accepted for publication in Astrophysical Journal (2004)

Journal ref: Astrophys.J. 611 (2004) 575-586

Probing the Edge of the Solar System: Formation of an Unstable Jet-Sheet

Authors: M. Opher, P. C. Liewer, T. I. Gombosi, W. Manchester, D. L. DeZeeuw, I. Sokolov, G. Toth

Abstract: The Voyager spacecraft is now approaching the edge of the solar system. Near the boundary between the solar system and the interstellar medium we find that an unstable ``jet-sheet'' forms. The jet-sheet oscillates up and down due to a velocity shear instability. This result is due to a novel application of a state-of-art 3D Magnetohydrodynamic (MHD) code with a highly refined grid. We assume as… ▽ More The Voyager spacecraft is now approaching the edge of the solar system. Near the boundary between the solar system and the interstellar medium we find that an unstable ``jet-sheet'' forms. The jet-sheet oscillates up and down due to a velocity shear instability. This result is due to a novel application of a state-of-art 3D Magnetohydrodynamic (MHD) code with a highly refined grid. We assume as a first approximation that the solar magnetic and rotation axes are aligned. The effect of a tilt of the magnetic axis with respect to the rotation axis remains to be seen. We include in the model self-consistently magnetic field effects in the interaction between the solar and interstellar winds. Previous studies of this interaction had poorer spatial resolution and did not include the solar magnetic field. This instability can affect the entry of energetic particles into the solar system and the intermixing of solar and interstellar material. The same effect found here is predicted for the interaction of rotating magnetized stars possessing supersonic winds and moving with respect to the interstellar medium, such as O stars. △ Less

Submitted 21 May, 2003; originally announced May 2003.

Comments: 9 pages, 4 figures, accepted for publication in ApJL

Journal ref: Astrophys.J. 591 (2003) L61-L66

On Gravitational Collapse in the Nonsymmetric Gravitational Theory

Authors: J. W. Moffat, I. Yu. Sokolov

Abstract: The analytical structure of the difference between the static vacuum solution in the nonsymmetric gravitational theory (NGT) and the Schwarzschild solution of Einstein's gravitational theory (EGT) is studied. It is proved that a smooth matching of the solutions does not exist in the range $0 < r \leq 2M$, for any non-zero values of the parameters $M$ and $s$ of the NGT solution. This means that… ▽ More The analytical structure of the difference between the static vacuum solution in the nonsymmetric gravitational theory (NGT) and the Schwarzschild solution of Einstein's gravitational theory (EGT) is studied. It is proved that a smooth matching of the solutions does not exist in the range $0 < r \leq 2M$, for any non-zero values of the parameters $M$ and $s$ of the NGT solution. This means that one cannot consider the difference between the two solutions using perturbation theory in this range of $r$. Assuming that the exterior solution in gravitational collapse is a small, time dependent perturbation of the static solution for a non-zero, generic NGT source ($s\not=0$) and mass density, it is shown that the matching of the interior and exterior solutions will not lead to black hole event horizons. △ Less

Submitted 12 October, 1995; originally announced October 1995.

Comments: 16 pages, REVTeX 3.0

Report number: UTPT-95-21

Galaxy Dynamics Predictions in the Nonsymmetric Gravitational Theory

Authors: J. W. Moffat, I. Yu. Sokolov

Abstract: In the weak field approximation, the nonsymmetric gravitational theory (NGT) has, in addition to the Newtonian gravitational potential, a Yukawa-like potential produced by the exchange of a spin $1^+$ boson between fermions. If the range $r_0=μ^{-1}$ is $25$ kpc, then this additional potential due to the interaction with matter in the halos of galaxies can explain the flat rotation curves of gal… ▽ More In the weak field approximation, the nonsymmetric gravitational theory (NGT) has, in addition to the Newtonian gravitational potential, a Yukawa-like potential produced by the exchange of a spin $1^+$ boson between fermions. If the range $r_0=μ^{-1}$ is $25$ kpc, then this additional potential due to the interaction with matter in the halos of galaxies can explain the flat rotation curves of galaxies and the Tully-Fisher law ($L\sim v^4$) without the dark matter hypothesis. Possible fits to clusters of galaxies and gravitational lensing observations are discussed. The results are based on a linear approximation to a new perturbatively consistent version of the NGT field equations, which does not violate the weak equivalence principle. △ Less

Submitted 3 October, 1995; v1 submitted 27 September, 1995; originally announced September 1995.

Comments: 16 pages, REVTeX 3.0; correction of a misprinting

Report number: UTPT-95-17

Journal ref: Phys.Lett.B378:59-67,1996