Showing 1–46 of 46 results for author: Prasad, A

From eruption to post-flare rain: a 2.5D MHD model

Authors: Samrat Sen, Avijeet Prasad, Valeriia Liakh, Rony Keppens

Abstract: The formation of the MFRs in the pre-flare stage, and how this leads to coronal rain in a post-eruption magnetic loop is not fully understood. We explore the formation, and eruption of MFRs, followed by the appearance of coronal rain in the post-flare loops, to understand the magnetic and thermodynamic properties of eruptive events and their multi-thermal aspects in the solar atmosphere. We perfor… ▽ More The formation of the MFRs in the pre-flare stage, and how this leads to coronal rain in a post-eruption magnetic loop is not fully understood. We explore the formation, and eruption of MFRs, followed by the appearance of coronal rain in the post-flare loops, to understand the magnetic and thermodynamic properties of eruptive events and their multi-thermal aspects in the solar atmosphere. We perform a resistive-magnetohydrodynamic (MHD) simulation with the open-source code \texttt{MPI-AMRVAC} to explore the evolution of sheared magnetic arcades that can lead to flux rope eruptions. The system is in mechanical imbalance at the initial state, and evolves self-consistently in a non-adiabatic atmosphere under the influence of radiative losses, thermal conduction, and background heating. We use an additional level of adaptive mesh refinement to achieve the smallest cell size of $\approx 32.7$ km in each direction to reveal the fine structures in the system. The system achieves a semi-equilibrium state after a short transient evolution from its initial mechanically imbalanced condition. A series of erupting MFRs is formed due to spontaneous magnetic reconnection, across current sheets created underneath the erupting flux ropes. Gradual development of thermal imbalance is noticed at a loop top in the post-eruption phase, which leads to catastrophic cooling and formation of condensations. We obtain plasma blobs which fall down along the magnetic loop in the form of coronal rain. The dynamical and thermodynamic properties of these cool-condensations are in good agreement with observations of post-flare coronal rain. The presented simulation supports the development and eruption of multiple MFRs, and the formation of coronal rain in post-flare loops, which is one of the key aspects to reveal the coronal heating mystery in the solar atmosphere. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Accepted for publication in Astronomy and Astrophysics journal. The animation associated with Figure 3 is available on request from the corresponding author. Comments are welcome

Characterizing the solar cycle variability using nonlinear time series analysis at different amounts of dynamo supercriticality: Solar dynamo is not highly supercritical

Authors: Aparup Ghosh, Pawan Kumar, Amrita Prasad, Bidya Binay Karak

Abstract: The solar dynamo is essentially a cyclic process in which the toroidal component of the magnetic field is converted into the poloidal one and vice versa. This cyclic loop is disturbed by some nonlinear and stochastic processes mainly operating in the toroidal to poloidal part. Hence, the memory of the polar field decreases in every cycle. On the other hand, the dynamo efficiency and, thus, the sup… ▽ More The solar dynamo is essentially a cyclic process in which the toroidal component of the magnetic field is converted into the poloidal one and vice versa. This cyclic loop is disturbed by some nonlinear and stochastic processes mainly operating in the toroidal to poloidal part. Hence, the memory of the polar field decreases in every cycle. On the other hand, the dynamo efficiency and, thus, the supercriticality of the dynamo decreases with the Sun's age. Previous studies have shown that the memory of the polar magnetic field decreases with the increase of supercriticality of the dynamo. In this study, we employ popular techniques of time series analysis, namely, compute Higuchi's fractal dimension, Hurst exponent, and Multi-Fractal Detrended Fluctuation Analysis, to the amplitude of the solar magnetic cycle obtained from dynamo models operating at near-critical and supercritical regimes. We show that the magnetic field in the near-critical regime is governed by strong memory, less stochasticity, intermittency, and breakdown of self-similarity. On the contrary, the magnetic field in the supercritical region has less memory, strong stochasticity, and shows a good amount of self-similarity. Finally, applying the same time series analysis techniques in the reconstructed sunspot data of 85 cycles and comparing their results with that from models, we conclude that the solar dynamo is possibly operating near the critical regime and not too much supercritical regime. Thus Sun may not be too far from the critical dynamo transition. △ Less

Submitted 10 March, 2024; originally announced March 2024.

Comments: Accepted for publication in the Astronomical Journal

VaTEST III: Validation of 8 Potential Super-Earths from TESS Data

Authors: Priyashkumar Mistry, Aniket Prasad, Mousam Maity, Kamlesh Pathak, Sarvesh Gharat, Georgios Lekkas, Surendra Bhattarai, Dhruv Kumar, Jack J. Lissauer, Joseph D. Twicken, Abderahmane Soubkiou, Francisco J. Pozuelos, Jon Jenkins, Keith Horne, Steven Giacalone, Khalid Barkaoui, Mathilde Timmermans, Cristilyn N. Watkins, Ramotholo Sefako, Karen A. Collins, Avi Shporer, Zouhair Benkhaldoun, Chris Stockdale, Emmanuël Jehin, Felipe Murgas , et al. (7 additional authors not shown)

Abstract: NASA's all-sky survey mission, the Transiting Exoplanet Survey Satellite (TESS), is specifically engineered to detect exoplanets that transit bright stars. Thus far, TESS has successfully identified approximately 400 transiting exoplanets, in addition to roughly 6000 candidate exoplanets pending confirmation. In this study, we present the results of our ongoing project, the Validation of Transitin… ▽ More NASA's all-sky survey mission, the Transiting Exoplanet Survey Satellite (TESS), is specifically engineered to detect exoplanets that transit bright stars. Thus far, TESS has successfully identified approximately 400 transiting exoplanets, in addition to roughly 6000 candidate exoplanets pending confirmation. In this study, we present the results of our ongoing project, the Validation of Transiting Exoplanets using Statistical Tools (VaTEST). Our dedicated effort is focused on the confirmation and characterization of new exoplanets through the application of statistical validation tools. Through a combination of ground-based telescope data, high-resolution imaging, and the utilization of the statistical validation tool known as \texttt{TRICERATOPS}, we have successfully discovered eight potential super-Earths. These planets bear the designations: TOI-238b (1.61$^{+0.09} _{-0.10}$ R$_\oplus$), TOI-771b (1.42$^{+0.11} _{-0.09}$ R$_\oplus$), TOI-871b (1.66$^{+0.11} _{-0.11}$ R$_\oplus$), TOI-1467b (1.83$^{+0.16} _{-0.15}$ R$_\oplus$), TOI-1739b (1.69$^{+0.10} _{-0.08}$ R$_\oplus$), TOI-2068b (1.82$^{+0.16} _{-0.15}$ R$_\oplus$), TOI-4559b (1.42$^{+0.13} _{-0.11}$ R$_\oplus$), and TOI-5799b (1.62$^{+0.19} _{-0.13}$ R$_\oplus$). Among all these planets, six of them fall within the region known as 'keystone planets,' which makes them particularly interesting for study. Based on the location of TOI-771b and TOI-4559b below the radius valley we characterized them as likely super-Earths, though radial velocity mass measurements for these planets will provide more details about their characterization. It is noteworthy that planets within the size range investigated herein are absent from our own solar system, making their study crucial for gaining insights into the evolutionary stages between Earth and Neptune. △ Less

Submitted 2 April, 2024; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: Accepted: Publications of the Astronomical Society of Australia

Formation of an observed eruptive flux rope above the torus instability threshold through tether-cutting magnetic reconnection

Authors: Avijeet Prasad, Sanjay Kumar, Alphonse C. Sterling, Ronald L. Moore, Guillaume Aulanier, R. Bhattacharyya, Qiang Hu

Abstract: Erupting magnetic flux ropes (MFRs) play a crucial role in producing solar flares. However, the formation of erupting MFRs in complex coronal magnetic configurations and their subsequent evolution in the flaring events are not fully understood. We performed an MHD simulation of active region NOAA 12241 to understand the formation of a rising MFR during the onset of an M6.9 flare on 2014 December 1… ▽ More Erupting magnetic flux ropes (MFRs) play a crucial role in producing solar flares. However, the formation of erupting MFRs in complex coronal magnetic configurations and their subsequent evolution in the flaring events are not fully understood. We performed an MHD simulation of active region NOAA 12241 to understand the formation of a rising MFR during the onset of an M6.9 flare on 2014 December 18, around 21:41 UT. The MHD simulation was initialised with an extrapolated non-force-free magnetic field generated from the photospheric vector magnetogram of the active region taken a few minutes before the flare. The initial magnetic field topology displays a pre-existing sheared arcade enveloping the polarity inversion line. The simulated dynamics exhibit the movement of the oppositely directed legs of the sheared arcade field lines towards each other due to the converging Lorentz force, resulting in the onset of tether-cutting magnetic reconnection that produces an underlying flare arcade and flare ribbons. Concurrently, an MFR above the flare arcade develops inside the sheared arcade and shows a rising motion. The MFR is found to be formed in a torus-unstable region, thereby explaining its eruptive nature. Interestingly, the location and rise of the rope are in good agreement with the corresponding observations seen in EUV channels. Furthermore, the foot points of the simulation's flare arcade match well with the location of the observed parallel ribbons of the flare. The presented simulation supports the development of the MFR by the tether-cutting magnetic reconnection inside the sheared coronal arcade during flare onset. The MFR is then found to extend along the polarity inversion line (PIL) through slip-running reconnection. The MFR's eruptive nature is ascribed both to its formation in the torus-unstable region and also to the runaway tether-cutting reconnection. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: 12 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A43 (2023)

Magnetohydrodynamics simulation of magnetic flux rope formation in a quadrupolar magnetic field configuration

Authors: Sanjay Kumar, Avijeet Prasad, Sushree S. Nayak, Satyam Agarwal, R. Bhattacharyya

Abstract: Magnetic flux ropes (MFRs) play an important role in high-energetic events like solar flares and coronal mass ejections in the solar atmosphere. Importantly, solar observations suggest an association of some flaring events with quadrupolar magnetic configurations. However, the formation and subsequent evolution of MFRs in such magnetic configurations still need to be fully understood. In this pape… ▽ More Magnetic flux ropes (MFRs) play an important role in high-energetic events like solar flares and coronal mass ejections in the solar atmosphere. Importantly, solar observations suggest an association of some flaring events with quadrupolar magnetic configurations. However, the formation and subsequent evolution of MFRs in such magnetic configurations still need to be fully understood. In this paper, we present idealized magnetohydrodynamics (MHD) simulations of MFR formation in a quadrupolar magnetic configuration. A suitable initial magnetic field having a quadrupolar configuration is constructed by modifying a three-dimensional (3D) linear force-free magnetic field. The initial magnetic field contains neutral lines, which consist of X-type null points. The simulated dynamics initially demonstrate the oppositely directed magnetic field lines located across the polarity inversion lines (PILs) moving towards each other, resulting in magnetic reconnections. Due to these reconnections, four highly twisted MFRs form over the PILs. With time, the foot points of the MFRs move towards the X-type neutral lines and reconnect, generating complex magnetic structures around the neutral lines, thus making the MFR topology more complex in the quadrupolar configuration than those formed in bipolar loop systems. Further evolution reveals the non-uniform rise of the MFRs. Importantly, the simulations indicate that the pre-existing X-type null points in magnetic configurations can be crucial to the evolution of the MFRs and may lead to the observed brightenings during the onset of some flaring events in the quadrupolar configurations. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 18 pages, 9 figures

Journal ref: Plasma Phys. Control. Fusion 2023 65 085008

Coronal Magnetic Field Extrapolation and Topological Analysis of Fine-Scale Structures during Solar Flare Precursors

Authors: Wen He, Qiang Hu, Ju Jing, Haimin Wang, Chaowei Jiang, Sushree S. Nayak, Avijeet Prasad

Abstract: Magnetic field plays an important role in various solar eruptions like flares, coronal mass ejections, etc. The formation and evolution of characteristic magnetic field topology in solar eruptions are critical problems that will ultimately help us understand the origination of these eruptions in the solar source regions. With the development of advanced techniques and instruments, observations wit… ▽ More Magnetic field plays an important role in various solar eruptions like flares, coronal mass ejections, etc. The formation and evolution of characteristic magnetic field topology in solar eruptions are critical problems that will ultimately help us understand the origination of these eruptions in the solar source regions. With the development of advanced techniques and instruments, observations with higher resolutions in different wavelengths and fields of view have provided more quantitative information for finer structures. So it is essential to improve our method to study the magnetic field topology in the solar source regions by taking advantage of high-resolution observations. In this study, we employ a nonlinear force-free field (NLFFF) extrapolation method based on a nonuniform grid setting for an M-class flare eruption event (SOL2015-06-22T17:39) with embedded magnetograms from the Solar Dynamics Observatory (SDO) and the Goode Solar Telescope (GST). The extrapolation results employing the embedded magnetogram for the bottom boundary are obtained by maintaining the native resolutions of the corresponding GST and SDO magnetograms. We compare the field line connectivity with the simultaneous GST/H$α$ and SDO/AIA observations for fine-scale structures associated with precursor brightenings. Then we perform a topological analysis of the field line connectivity corresponding to fine-scale magnetic field structures based on the extrapolation results. The results indicate that by combining the high-resolution GST magnetogram with a larger HMI magnetogram, the derived magnetic field topology is consistent with a scenario of magnetic reconnection among sheared field lines across the main polarity inversion line during solar flare precursors. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 17 pages, 9 figures, 2 tables, submitted to ApJ

VaTEST II: Statistical Validation of 11 TESS-Detected Exoplanets Orbiting K-type Stars

Authors: Priyashkumar Mistry, Kamlesh Pathak, Aniket Prasad, Georgios Lekkas, Surendra Bhattarai, Sarvesh Gharat, Mousam Maity, Dhruv Kumar, Karen A. Collins, Richard P. Schwarz, Christopher R. Mann, Elise Furlan, Steve B. Howell, David Ciardi, Allyson Bieryla, Elisabeth C. Matthews, Erica Gonzales, Carl Ziegler, Ian Crossfield, Steven Giacalone, Thiam-Guan Tan, Phil Evans, Krzysztof G. Helminiak, Kevin I. Collins, Norio Narita , et al. (26 additional authors not shown)

Abstract: NASA's Transiting Exoplanet Survey Satellite (TESS) is an all-sky survey mission designed to find transiting exoplanets orbiting nearby bright stars. It has identified more than 329 transiting exoplanets, and almost 6,000 candidates remain unvalidated. In this manuscript, we discuss the findings from the ongoing VaTEST (Validation of Transiting Exoplanets using Statistical Tools) project, which ai… ▽ More NASA's Transiting Exoplanet Survey Satellite (TESS) is an all-sky survey mission designed to find transiting exoplanets orbiting nearby bright stars. It has identified more than 329 transiting exoplanets, and almost 6,000 candidates remain unvalidated. In this manuscript, we discuss the findings from the ongoing VaTEST (Validation of Transiting Exoplanets using Statistical Tools) project, which aims to validate new exoplanets for further characterization. We validated 11 new exoplanets by examining the light curves of 24 candidates using the LATTE and TESS-Plot tools and computing the False Positive Probabilities using the statistical validation tool TRICERATOPS. These include planets suitable for atmospheric characterization using transmission spectroscopy (TOI-2194b), emission spectroscopy (TOI-3082b and TOI-5704b) and for both transmission and emission spectroscopy (TOI-672b, TOI- 1694b, and TOI-2443b); One super-Earth (TOI-2194b) orbiting a bright (V = 8.42 mag), metal-poor ([Fe/H] = -0.3720 $\pm$ 0.1) star; one short-period Neptune-like planet (TOI-5704) in the Hot Neptune Desert. In total, we validated 1 super-Earth, 7 sub-Neptunes, 1 Neptune-like, and 2 sub-Saturn or super-Neptune-like exoplanets. Additionally, we identify five likely planet candidates (TOI-323, TOI- 1180, TOI-2200, TOI-2408 and TOI-3913) which can be further studied to establish their planetary nature. △ Less

Submitted 13 May, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: Accepted for Publication in Astronomical Journal, 28 Pages, 7 Figures

Evolution of magnetic fields and energy release processes during homologous eruptive flares

Authors: Suraj Sahu, Bhuwan Joshi, Avijeet Prasad, Kyung-Suk Cho

Abstract: We explore the processes of repetitive build-up and explosive release of magnetic energy together with the formation of magnetic flux ropes that eventually resulted into three homologous eruptive flares of successively increasing intensities (i.e., M2.0, M2.6, and X1.0). The flares originated from NOAA active region 12017 during 2014 March 28-29. EUV observations and magnetogram measurements toget… ▽ More We explore the processes of repetitive build-up and explosive release of magnetic energy together with the formation of magnetic flux ropes that eventually resulted into three homologous eruptive flares of successively increasing intensities (i.e., M2.0, M2.6, and X1.0). The flares originated from NOAA active region 12017 during 2014 March 28-29. EUV observations and magnetogram measurements together with coronal magnetic field modeling suggest that the flares were triggered by the eruption of flux ropes embedded by a densely packed system of loops within a small part of the active region. In X-rays, the first and second events show similar evolution with single, compact sources, while the third event exhibits multiple emission centroids with a set of strong non-thermal conjugate sources at 50-100 keV during the HXR peak. The photospheric magnetic field over an interval of approximately 44 hr encompassing the three flares undergoes important phases of emergence and cancellation processes together with significant changes near the polarity inversion lines within the flaring region. Our observations point toward the tether-cutting mechanism as the plausible triggering process of the eruptions. Between the second and third event, we observe a prominent phase of flux emergence which temporally correlates with the build-up phase of free magnetic energy in the active region corona. In conclusion, our analysis reveals an efficient coupling between the rapidly evolving photospheric and coronal magnetic fields in the active region that led to a continued phase of the build-up of free energy, resulting into the homologous flares of successively increasing intensities. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: Accepted for publication in 'The Astrophysical Journal'

VaTEST I: Validation of Sub-Saturn Exoplanet TOI-181b in Narrow Orbit from its Host Star

Authors: Priyashkumar Mistry, Kamlesh Pathak, Georgios Lekkas, Aniket Prasad, Surendra Bhattarai, Mousam Maity, Charles A. Beichman, David R. Ciardi, Phil Evans, Allyson Bieryla, Jason D. Eastman, David W. Latham, Gilbert A. Esquerdo, Jennifer P. Lucero

Abstract: We present here a validation of sub-Saturn exoplanet TOI-181b orbiting a K spectral type star TOI-181 (Mass: 0.822 $\pm$ 0.04 M$_{\odot}$, Radius: 0.745 $\pm$ 0.02 R$_{\odot}$, Temperature: 4994 $\pm$ 50 K) as a part of Validation of Transiting Exoplanets using Statistical Tools (VaTEST) project. TOI-181b is a planet with radius 6.95 $\pm$ 0.08 R$_{\oplus}$, mass 46.16 $\pm$ 2.71 M$_{\oplus}$, orb… ▽ More We present here a validation of sub-Saturn exoplanet TOI-181b orbiting a K spectral type star TOI-181 (Mass: 0.822 $\pm$ 0.04 M$_{\odot}$, Radius: 0.745 $\pm$ 0.02 R$_{\odot}$, Temperature: 4994 $\pm$ 50 K) as a part of Validation of Transiting Exoplanets using Statistical Tools (VaTEST) project. TOI-181b is a planet with radius 6.95 $\pm$ 0.08 R$_{\oplus}$, mass 46.16 $\pm$ 2.71 M$_{\oplus}$, orbiting in a slightly eccentric orbit with eccentricity 0.15 $\pm$ 0.06 and semi-major axis of 0.054 $\pm$ 0.004 AU, with an orbital period of 4.5320 $\pm$ 0.000002 days. The transit photometry data was collected using Transiting Exoplanet Survey Satellite (TESS) and spectroscopic data for radial velocity analysis was collected using The European Southern Observatory's (ESO) High Accuracy Radial Velocity Planet Searcher (HARPS) telescope. Based on the radial velocity best-fit model we measured RV semi-amplitude to be 20.56 $\pm$ 2.37 m s$^{-1}$. Additionally, we used \texttt{VESPA} and \texttt{TRICERATOPS} to compute the False Positive Probability (FPP), and the findings were FPP values of $1.68\times10^{-14}$ and $3.81\times10^{-04}$, respectively, which are significantly lower than the 1% threshold. The finding of TOI-181b is significant in the perspective of future work on the formation and migration history of analogous planetary systems since warm sub-Saturns are uncommon in the known sample of exoplanets. △ Less

Submitted 10 January, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: 13 pages, 10 figures, submitted to MNRAS, Structural changes are made, New results

Magnetohydrodynamics evolution of three-dimensional magnetic null in NOAA active region 11515 initiated using non-force-free field extrapolation

Authors: Sanjay Kumar, Avijeet Prasad, Ranadeep Sarkar, Ramit Bhattacharyya

Abstract: Magnetohydrodynamics simulation of active region NOAA 11515 is performed to examine the initiation of the M5.6 flaring event that starts around 10:43 UT on 2012 July 2. The simulation is conducted using an extrapolated non-force-free magnetic field generated from the photospheric vector magnetogram of the active region as the initial magnetic field. The magnetic field shows the presence of a three… ▽ More Magnetohydrodynamics simulation of active region NOAA 11515 is performed to examine the initiation of the M5.6 flaring event that starts around 10:43 UT on 2012 July 2. The simulation is conducted using an extrapolated non-force-free magnetic field generated from the photospheric vector magnetogram of the active region as the initial magnetic field. The magnetic field shows the presence of a three-dimensional (3D) magnetic null with the corresponding dome overlying a filament and a low-lying magnetic flux rope, observed in 304~Å~ and 131~Å~ respectively. The simulated dynamics, triggered by the initial Lorentz force, lead to the bifurcations of the flux rope, which is similar to the observed bifurcation in the 131 Å~ brightenings. Additionally, the rope exhibits a rise and reconnects at the 3D null. These reconnections convert field lines of the rope into the anchored outer spine of the 3D null -- explaining the occurrence of a nearby confined C-class flare. Further, the results show that the field lines of the flux rope reach the vicinity of the filament and become non-parallel to the field lines of the filament. This initiates the reconnections between the rope and the field lines of the filament -- activating the filament for the eruption. This interesting interaction of the flux rope and filament seems to contribute to the onset of the M-class flare. △ Less

Submitted 8 October, 2022; originally announced October 2022.

Comments: 10 Figures

Spectroscopic study of solar transition region oscillations in the quiet-Sun observed by IRIS using Si IV spectral line

Authors: Kartika Sangal, A. K. Srivastava, P. Kayshap, T. J Wang, J. J. González-Avilés, Abhinav Prasad

Abstract: In the present paper, we use Si IV 1393.755 Å spectral line observed by the Interface Region Imaging Spectrograph (IRIS) in the quiet-Sun to determine physical nature of the solar transition region (TR) oscillations. We analyze the properties of these oscillations using wavelet tools (e.g., power, cross-power, coherence, and phase difference) along with the stringent noise model (i.e., power-law +… ▽ More In the present paper, we use Si IV 1393.755 Å spectral line observed by the Interface Region Imaging Spectrograph (IRIS) in the quiet-Sun to determine physical nature of the solar transition region (TR) oscillations. We analyze the properties of these oscillations using wavelet tools (e.g., power, cross-power, coherence, and phase difference) along with the stringent noise model (i.e., power-law + constant). We estimate the period of the intensity and Doppler velocity oscillations at each chosen location in the quiet-Sun (QS) and quantify the distribution of the statistically significant power and associated periods in one bright and two dark regions. In the bright TR region, the mean periods in intensity and velocity are 7 min, and 8 min respectively. In the dark region, the mean periods in intensity and velocity are 7 min, and 5.4 min respectively. We also estimate the phase difference between the intensity and Doppler velocity oscillations at each location. The statistical distribution of phase difference is estimated, which peaks at -119\degree $\pm$ 13\degree, 33\degree $\pm$ 10\degree, 102\degree $\pm$ 10\degree\ in the bright region, while at -153\degree $\pm$ 13\degree, 6\degree $\pm$ 20\degree, 151\degree $\pm$ 10\degree\ in the dark region. The statistical distribution reveals that the oscillations are caused by propagating slow magnetoacoustic waves encountered with the TR. Some of these locations may also be associated with the standing slow waves. Even, in the given time domain, several locations exhibit presence of both propagating and standing oscillations at different frequencies. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 16 pages, 11 figures

Validation and interpretation of three-dimensional configuration of a magnetic cloud flux rope

Authors: Qiang Hu, Chunming Zhu, Wen He, Jiong Qiu, Lan K. Jian, Avijeet Prasad

Abstract: One "strong" magnetic cloud (MC) with the magnetic field magnitude reaching $\sim$ 40 nT at 1 au during 2012 June 16-17 is examined in association with a pre-existing magnetic flux rope (MFR) identified on the Sun. The MC is characterized by a quasi-three dimensional (3D) flux rope model based on in situ measurements from the Wind spacecraft. The magnetic flux contents and other parameters are qua… ▽ More One "strong" magnetic cloud (MC) with the magnetic field magnitude reaching $\sim$ 40 nT at 1 au during 2012 June 16-17 is examined in association with a pre-existing magnetic flux rope (MFR) identified on the Sun. The MC is characterized by a quasi-three dimensional (3D) flux rope model based on in situ measurements from the Wind spacecraft. The magnetic flux contents and other parameters are quantified. In addition, a correlative study with the corresponding measurements of the same structure crossed by the Venus Express (VEX) spacecraft at a heliocentric distance 0.7 au and with an angular separation $\sim 6^\circ$ in longitude is performed to validate the MC modeling results. The spatial variation between the Wind and VEX magnetic field measurements is attributed to the 3D configuration of the structure as featured by a knotted bundle of flux. The comparison of the magnetic flux contents between the MC and the source region on the Sun indicates that the 3D reconnection process accompanying an M1.9 flare may correspond to the magnetic reconnection between the field lines of the pre-existing MFR rooted in the opposite polarity footpoints. Such a process reduces the amount of the axial magnetic flux in the erupted flux rope, by approximately 50\%, in this case. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: Submitted to ApJ

The magnetic topology of the inverse Evershed flow

Authors: A. Prasad, M. Ranganathan, C. Beck, D. P. Choudhary, Q. Hu

Abstract: The inverse Evershed flow (IEF) is a mass motion towards sunspots at chromospheric heights. We combined high-resolution observations of NOAA 12418 from the Dunn Solar Telescope and vector magnetic field measurements from the Helioseismic and Magnetic Imager (HMI) to determine the driver of the IEF. We derived chromospheric line-of-sight (LOS) velocities from spectra of H$α$ and Ca II IR. The HMI d… ▽ More The inverse Evershed flow (IEF) is a mass motion towards sunspots at chromospheric heights. We combined high-resolution observations of NOAA 12418 from the Dunn Solar Telescope and vector magnetic field measurements from the Helioseismic and Magnetic Imager (HMI) to determine the driver of the IEF. We derived chromospheric line-of-sight (LOS) velocities from spectra of H$α$ and Ca II IR. The HMI data were used in a non-force-free magnetic field extrapolation to track closed field lines near the sunspot in the active region. We determined their length and height, located their inner and outer foot points, and derived flow velocities along them. The magnetic field lines related to the IEF reach on average a height of 3 Mm over a length of 13 Mm. The inner (outer) foot points are located at 1.2 (1.9) sunspot radii. The average field strength difference $ΔB$ between inner and outer foot points is +400 G. The temperature difference $ΔT$ is anti-correlated with $ΔB$ with an average value of -100 K. The pressure difference $Δp$ is dominated by $ΔB$ and is primarily positive with a driving force towards the inner foot points of 1.7 kPa on average. The velocities predicted from $Δp$ reproduce the LOS velocities of 2-10 km s$^{-1}$ with a square-root dependence. We find that the IEF is driven along magnetic field lines connecting network elements with the outer penumbra by a gas pressure difference that results from a difference in field strength as predicted by the classical siphon flow scenario. △ Less

Submitted 5 March, 2022; originally announced March 2022.

Comments: 23 pages, 27 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 662, A25 (2022)

Quantitative Characterization of Magnetic Flux Rope Properties for Two Solar Eruption Events

Authors: Wen He, Qiang Hu, Chaowei Jiang, Jiong Qiu, Avijeet Prasad

Abstract: In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR du… ▽ More In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR due to the lack of direct magnetic field measurement in the corona. In this study, we investigate the MFR structures, originating from two active regions (ARs), AR 11719 and AR 12158, and estimate their magnetic properties quantitatively. We perform the nonlinear force-free field extrapolations with preprocessed photospheric vector magnetograms. In addition, remote-sensing observations are employed to find indirect evidence of MFRs on the Sun and to analyze the time evolution of magnetic reconnection flux associated with the flare ribbons during the eruption. A coherent "pre-existing" MFR structure prior to the flare eruption is identified quantitatively for one event from the combined analysis of the extrapolation and observation. Then the characteristics of MFRs for two events on the Sun before and during the eruption, forming the CME-MFR, including the axial magnetic flux, field-line twist, and reconnection flux, are estimated and compared with the corresponding in situ modeling results. We find that the magnetic reconnection associated with the accompanying flares for both events injects significant amount of flux into the erupted CME-MFRs. △ Less

Submitted 9 January, 2022; originally announced January 2022.

Comments: 28 pages, 10 figures, 3 tables, submitted to ApJ

Role of Non-Ideal Dissipation with Heating-Cooling Misbalance on the Phase Shifts of Standing Slow Magnetohydrodynamic Waves

Authors: Abhinav Prasad, A. K. Srivastava, Tongjiang Wang, Kartika Sangal

Abstract: We analyse the phase shifts of standing, slow magnetohydrodynamic (MHD) waves in solar coronal loops using a linear MHD model taking into account the role of thermal conductivity, compressive viscosity, radiative losses, and heating-cooling misbalance. We estimate the phase shifts in time and space of density and temperature perturbations with respect to velocity perturbations and also calculate t… ▽ More We analyse the phase shifts of standing, slow magnetohydrodynamic (MHD) waves in solar coronal loops using a linear MHD model taking into account the role of thermal conductivity, compressive viscosity, radiative losses, and heating-cooling misbalance. We estimate the phase shifts in time and space of density and temperature perturbations with respect to velocity perturbations and also calculate the phase difference between density and temperature perturbations. The overall significance of compressive viscosity is found to be negligible for most of the loops considered in the study. For loops with high background density and/or low background temperature, the role of radiative losses (with heating-cooling misbalance) is found to be more significant. Also the effect of heating-cooling misbalance with a temperature- and density-dependent heating function is found to be more significant in the case of longer loop lengths ($L=500$\, Mm). We derived a general expression for the polytropic index [$γ_{\rm eff}$] and found that under linear MHD the effect of compressive viscosity on polytropic index is negligible. The radiative losses with constant heating lead to a monotonic increase of $γ_{\rm eff}$ with increasing density whereas the consideration of an assumed heating function [$H(ρ,T) \propto ρ^{a}T^{b}$, where $a=-0.5$ and $b=-3$] makes the $γ_{\rm eff}$ peak at a certain loop density. We also explored the role of different heating functions by varying the free parameters $a$ and $b$ for a fixed loop of $ρ_0 = 10^{-11}$\, kg $\text{m}^{-3}$, $T_0 = 6.3$\, MK and loop length $L= 180$\, Mm. We find that the consideration of different heating functions [$H(ρ,T)$] leads to a significant variation in the phase difference between density and temperature perturbations; however, the polytropic index remains close to a value of 1.66. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: Solar Physics (in press)

Comparison of the Hall Magnetohydrodynamics and Magnetohydrodynamics evolution of a flaring solar active region

Authors: Kamlesh Bora, Ramit Bhattacharyya, Avijeet Prasad, Bhuwan Joshi, Qiang Hu

Abstract: This work analyzes the Hall magnetohydrodynamics (HMHD) and magnetohydrodynamics (MHD) numerical simulations of a flaring solar active region as a testbed while idealizing the coronal Alfvén speed to be of two orders of magnitude lesser. HMHD supports faster magnetic reconnection and shows richer complexity in magnetic field line evolution compared to the MHD. The magnetic reconnections triggering… ▽ More This work analyzes the Hall magnetohydrodynamics (HMHD) and magnetohydrodynamics (MHD) numerical simulations of a flaring solar active region as a testbed while idealizing the coronal Alfvén speed to be of two orders of magnitude lesser. HMHD supports faster magnetic reconnection and shows richer complexity in magnetic field line evolution compared to the MHD. The magnetic reconnections triggering the flare are explored by numerical simulations augmented with relevant multi-wavelength observations. The initial coronal magnetic field is constructed by non-force-free extrapolation of photospheric vector magnetic field. Magnetic structure involved in the flare is identified to be a flux rope, with its overlying magnetic field lines constituting the quasi-separatrix layers (QSLs) along with a three-dimensional null point and a null line. Compared to the MHD simulation, the HMHD simulation shows a higher and faster ascend of the rope together with the overlying field lines, which further reconnect at the QSL located higher up in the corona. The foot points of the field lines match better with the observations for the HMHD case with the central part of the flare ribbon located at the chromosphere. Additionally, field lines are found to rotate in a circular pattern in the HMHD, whereas no such rotation is seen in the MHD results. Interestingly, plasma is also observed to be rotating in a co-spatial chromospheric region, which makes the HMHD simulation more credible. Based on the aforementioned agreements, HMHD simulation is found to agree better with observations and, thus, opens up a novel avenue to explore. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: 31 pages, 17 figures (accepted for publication in The Astrophysical Journal)

Heating of the solar chromosphere in a sunspot light bridge by electric currents

Authors: Rohan E. Louis, Avijeet Prasad, Christian Beck, Debi Prasad Choudhary, Mehmet S. Yalim

Abstract: Context: Resistive Ohmic dissipation has been suggested as a mechanism for heating the solar chromosphere, but few studies have established this association. Aim: We aim to determine how Ohmic dissipation by electric currents can heat the solar chromosphere. Methods: We combine high-resolution spectroscopic Ca II data from the Dunn Solar Telescope and vector magnetic field observations from the He… ▽ More Context: Resistive Ohmic dissipation has been suggested as a mechanism for heating the solar chromosphere, but few studies have established this association. Aim: We aim to determine how Ohmic dissipation by electric currents can heat the solar chromosphere. Methods: We combine high-resolution spectroscopic Ca II data from the Dunn Solar Telescope and vector magnetic field observations from the Helioseismic and Magnetic Imager (HMI) to investigate thermal enhancements in a sunspot light bridge. The photospheric magnetic field from HMI was extrapolated to the corona using a non-force-free field technique that provided the three-dimensional distribution of electric currents, while an inversion of the chromospheric Ca II line with a local thermodynamic equilibrium and a nonlocal thermodynamic equilibrium spectral archive delivered the temperature stratifications from the photosphere to the chromosphere. Results: We find that the light bridge is a site of strong electric currents, of about 0.3 A/m^2 at the bottom boundary, which extend to about 0.7 Mm while decreasing monotonically with height. These currents produce a chromospheric temperature excess of about 600-800 K relative to the umbra. Only the light bridge, where relatively weak and highly inclined magnetic fields emerge over a duration of 13 hr, shows a spatial coincidence of thermal enhancements and electric currents. The temperature enhancements and the Cowling heating are primarily confined to a height range of 0.4-0.7 Mm above the light bridge. The corresponding increase in internal energy of 200 J/m^3 can be supplied by the heating in about 10 min. Conclusions: Our results provide direct evidence for currents heating the lower solar chromosphere through Ohmic dissipation. △ Less

Submitted 26 July, 2021; originally announced July 2021.

Comments: Accepted for publication in Astronomy & Astrophysics Letters: 7 pages, 5 figures

Journal ref: A&A 652, L4 (2021)

Effect of Thermal Conductivity, Compressive Viscosity and Radiative Cooling on the Phase Shift of Propagating Slow Waves with and without Heating-Cooling Imbalance

Authors: Abhinav Prasad, A. K. Srivastava, Tongjiang Wang

Abstract: We study phase shifts of the propagating slow waves in coronal loops invoking the effects of thermal conductivity, compressive viscosity, radiative losses and heating-cooling imbalance. We derive a general dispersion relation and solve it to determine phase shifts of density and temperature perturbations relative to the velocity and their dependence on equilibrium parameters ($ρ_0$, $T_0$). We est… ▽ More We study phase shifts of the propagating slow waves in coronal loops invoking the effects of thermal conductivity, compressive viscosity, radiative losses and heating-cooling imbalance. We derive a general dispersion relation and solve it to determine phase shifts of density and temperature perturbations relative to the velocity and their dependence on equilibrium parameters ($ρ_0$, $T_0$). We estimate phase difference ($Δφ$) between density and temperature perturbations and its dependence on $ρ_0$ and $T_0$. The effect of viscosity on the phase shifts was found negligible. The role of radiative losses along with h/c imbalance for chosen specific heating function ($H(ρ, T) \propto ρ^{-0.5} T^{-3}$) in determining phase shifts, is found to be significant for the high density and low temperature loops. The h/c imbalance can increase the phase difference ($Δφ\approx 140^\circ$) for low temperature loops compared to the constant heating case ($Δφ\approx 30^\circ$). We derive a general expression for the polytropic index. We find that in the presence of thermal conduction alone, the polytropic index remains close to its classical value for all the considered $ρ_0$ and $T_0$. However, it reduces to a value $1.2$ when loop density is decreased by an order of magnitude compared to its normal coronal value. We find that the inclusion of radiative losses, with or without h/c imbalance, cannot explain the observed polytropic index. The thermal ratio ($d$) needs to be enhanced by an order of magnitude, in order to explain its observed value $1.1 \pm 0.02$ in the solar loops. We also explore the role of different heating functions for typical coronal parameters and found that although the polytropic index remains close to its classical value, the phase difference is highly dependent on the form of heating function (The abstract is restructured for arxiv). △ Less

Submitted 30 June, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 32 Pages; 8 Figures; 1 Table; The final paper (arxiv version) as appeared in "Solar Physics" Journal after the editorial and proof correction (Prasad et al., Sol Phys., 296, 105 (2021))

Journal ref: Solar Physics, 296, 105 (2021)

Are the brightest coronal loops always rooted in mixed-polarity magnetic flux?

Authors: Sanjiv K. Tiwari, Caroline L. Evans, Navdeep K. Panesar, Avijeet Prasad, Ronald L. Moore

Abstract: A recent study demonstrated that freedom of convection and strength of magnetic field in the photospheric feet of active-region (AR) coronal loops, together, can engender or quench heating in them. Other studies stress that magnetic flux cancellation at the loop-feet potentially drives heating in loops. We follow 24-hour movies of a bipolar AR, using EUV images from SDO/AIA and line-of-sight (LOS)… ▽ More A recent study demonstrated that freedom of convection and strength of magnetic field in the photospheric feet of active-region (AR) coronal loops, together, can engender or quench heating in them. Other studies stress that magnetic flux cancellation at the loop-feet potentially drives heating in loops. We follow 24-hour movies of a bipolar AR, using EUV images from SDO/AIA and line-of-sight (LOS) magnetograms from SDO/HMI, to examine magnetic polarities at the feet of 23 of the brightest coronal loops. We derived FeXVIII emission (hot-94) images (using the Warren et al. method) to select the hottest/brightest loops, and confirm their footpoint locations via non-force-free field extrapolations. From 6"$\times$6" boxes centered at each loop foot in LOS magnetograms we find that $\sim$40\% of the loops have both feet in unipolar flux, and $\sim$60\% of the loops have at least one foot in mixed-polarity flux. The loops with both feet unipolar are $\sim$15\% shorter lived on average than the loops having mixed-polarity foot-point flux, but their peak-intensity averages are equal. The presence of mixed-polarity magnetic flux in at least one foot of majority of the loops suggests that flux cancellation at the footpoints may drive most of the heating. But, the absence of mixed-polarity magnetic flux (to the detection limit of HMI) in $\sim$40\% of the loops suggests that flux cancellation may not be necessary to drive heating in coronal loops -- magnetoconvection and field strength at both loop feet possibly drive much of the heating, even in the cases where a loop foot presents mixed-polarity magnetic flux. △ Less

Submitted 19 February, 2021; originally announced February 2021.

Comments: 11 pages, 5 figures

Journal ref: ApJ, 2021, 908, 151

Magnetic reconnections in the presence of three-dimensional magnetic nulls and quasi-separatrix layers

Authors: Sanjay Kumar, Sushree S. Nayak, Avijeet Prasad, Ramit Bhattacharyya

Abstract: Three-dimensional (3D) magnetohydrodynamic simulations are carried out to explore magnetic reconnections in the presence of 3D magnetic nulls and quasi-separatrix layers (QSLs). The initial magnetic fields are created by superposing uniform vertical magnetic fields of two different magnitudes on a linear force-free field. The interior of the numerical box contains two 3D nulls with separatrix dome… ▽ More Three-dimensional (3D) magnetohydrodynamic simulations are carried out to explore magnetic reconnections in the presence of 3D magnetic nulls and quasi-separatrix layers (QSLs). The initial magnetic fields are created by superposing uniform vertical magnetic fields of two different magnitudes on a linear force-free field. The interior of the numerical box contains two 3D nulls with separatrix domes separated by a quasi-separator (or hyperbolic flux tube) with QSLs. In the first simulation, the uniform vertical field is so large that the nulls are located at low heights and the domes are separate. Initially unbalanced Lorentz forces drive rotational flows that form strong electric currents and strong torsional fan reconnection at the 3D nulls and weak QSL reconnection at the hyperbolic flux tube. Flipping or slipping of field lines is observed in both cases. In the second simulation, with a weaker vertical field and larger domes, the separatrix surfaces meet at the central quasi-separator and their rotation drives stronger QSL reconnection than before. △ Less

Submitted 30 December, 2020; originally announced January 2021.

Comments: 24 pages, 14 figures, Accepted for the publication in Solar Physics

Role of Compressive Viscosity and Thermal Conductivity on the Damping of Slow Waves in the Coronal Loops With and Without Heating-Cooling Imbalance

Authors: Abhinav Prasad, A. K. Srivastava, T. J. Wang

Abstract: In the present paper, we derive a new dispersion relation for slow magnetoacoustic waves invoking the effect of thermal conductivity, compressive viscosity, radiation and unknown heating term along with the consideration of heating cooling imbalance from linearized MHD equations. We solve the general dispersion relation to understand role of compressive viscosity and thermal conductivity in dampin… ▽ More In the present paper, we derive a new dispersion relation for slow magnetoacoustic waves invoking the effect of thermal conductivity, compressive viscosity, radiation and unknown heating term along with the consideration of heating cooling imbalance from linearized MHD equations. We solve the general dispersion relation to understand role of compressive viscosity and thermal conductivity in damping of the slow waves in coronal loops with and without heating cooling imbalance. We have analyzed wave damping for the range of loop length $L$=50-500 Mm, temperature $T$=5-30 MK, and density $ρ$=10$^{-11}$-10$^{-9}$ kg m$^{-3}$. It was found that inclusion of compressive viscosity along with thermal conductivity significantly enhances the damping of fundamental mode oscillations in shorter (e.g., $L$=50 Mm) and super-hot ($T>$10 MK) loops. However, role of the viscosity in damping is insignificant in longer (e.g., $L$=500 Mm) and hot loops (T$\leq$10 MK) where, instead, thermal conductivity along with the presence of heating cooling imbalance plays a dominant role. For the shorter loops at the super-hot regime of the temperature, increment in loop density substantially enhances damping of the fundamental modes due to thermal conductivity when the viscosity is absent, however, when the compressive viscosity is added the increase in density substantially weakens damping. Thermal conductivity alone is found to play a dominant role in longer loops at lower temperatures (T$\leq$10 MK), while compressive viscosity dominates in damping at super-hot temperatures ($T>$10 MK) in shorter loops. The predicted scaling law between damping time ($τ$) and wave period ($P$) is found to better match to observed SUMER oscillations when heating cooling imbalance is taken into account in addition to thermal conductivity and compressive viscosity for the damping of the fundamental slow mode oscillations. △ Less

Submitted 5 January, 2021; v1 submitted 29 November, 2020; originally announced November 2020.

Comments: 40 Pages; 13 Figures; 01 Table; Solar Physics (in press)

Magnetohydrodynamic Simulation of Magnetic Null-point Reconnections and Coronal dimmings during the X2.1 flare in NOAA AR 11283

Authors: Avijeet Prasad, Karin Dissauer, Qiang Hu, R. Bhattacharyya, Astrid M. Veronig, Sanjay Kumar, Bhuwan Joshi

Abstract: The magnetohydrodynamics of active region NOAA 11283 is simulated using an initial non-force-free magnetic field extrapolated from its photospheric vector magnetogram. We focus on the magnetic reconnections at a magnetic null point that participated in the X2.1 flare on 2011 September 6 around 22:21 UT (SOL2011-09-06T22:21X2.1) followed by the appearance of circular flare ribbons and coronal dimmi… ▽ More The magnetohydrodynamics of active region NOAA 11283 is simulated using an initial non-force-free magnetic field extrapolated from its photospheric vector magnetogram. We focus on the magnetic reconnections at a magnetic null point that participated in the X2.1 flare on 2011 September 6 around 22:21 UT (SOL2011-09-06T22:21X2.1) followed by the appearance of circular flare ribbons and coronal dimmings. The initial magnetic field from extrapolation displays a three-dimensional (3D) null topology overlying a sheared arcade. Prior to the flare, magnetic loops rise due to the initial Lorentz force, and reconnect at the 3D null, leading to expansion and loss of confined plasma that produce the observed pre-flare coronal dimmings. Further, the simulated dynamics documents the transfer of twist from the arcade to the overlying loops through reconnections, developing a flux rope. The non-parallel field lines comprising the rope and lower-lying arcades form an X-type geometry. Importantly, the simultaneous reconnections at the 3D null and the X-type geometry can explain the observed circular and parallel flare ribbons. Reconnections at the 3D null transform closed inner spine field lines into open field lines of the outer spine. The footpoints of these open field lines correspond to a ring-shaped coronal dimming region, tracing the dome. Further, the flux rope bifurcates because of these reconnections which also results in the generation of open magnetic field lines. The plasma loss along the open field lines can potentially explain the observed coronal dimming. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 29 pages, 13 figures, accepted for publication in The Astrophysical Journal

An Eruptive Circular-ribbon Flare with Extended Remote Brightenings

Authors: Chang Liu, Avijeet Prasad, Jeongwoo Lee, Haimin Wang

Abstract: We study an eruptive X1.1 circular-ribbon flare on 2013 November 10, combining multiwavelength observations with a coronal field reconstruction using a non-force-free field method. In the first stage, a filament forms via magnetic reconnection between two mildly twisted sheared arcades, which are embedded under the fan dome associated with a null point. This reconnection seems to be driven by phot… ▽ More We study an eruptive X1.1 circular-ribbon flare on 2013 November 10, combining multiwavelength observations with a coronal field reconstruction using a non-force-free field method. In the first stage, a filament forms via magnetic reconnection between two mildly twisted sheared arcades, which are embedded under the fan dome associated with a null point. This reconnection seems to be driven by photospheric shearing and converging flows around the inner two arcade footpoints, consistent with the flare-related changes of transverse field. The southern portion of the filament rises upward due to torus instability and pushes against the null point. The induced null point reconnection then generates the circular ribbon and the initial remote brightening in the west, as accelerated electrons precipitate along the fan and propagate outward along quasi-separatix surfaces with high values of the squashing factor (Q) in the envelope fields, which have a curtain-like shape here. In the second stage, the southern end of the flux rope breaks away from the surface, sequentially disrupts the dome and overlying fields, and erupts in a whipping-like fashion to become a partial halo coronal mass ejection. This leads to an enhanced flare emission and fast-moving remote brightenings at the footpoints of the magnetic curtain, which span a remarkably broad region and are also associated with coronal dimmings. This is a rare example of eruptive circular-ribbon flares, in which the evolution of a flux rope from its formation to successful eruption out of the dome and the resulting unusually extended remote brightenings are completely observed. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Comments: 12 pages, 7 figures, Accepted to ApJ

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

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

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

Submitted 24 July, 2020; originally announced July 2020.

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

Effects of Cowling Resistivity in the Weakly-Ionized Chromosphere

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

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

Submitted 23 July, 2020; originally announced July 2020.

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

Forecasting Solar Cycle 25 using Deep Neural Networks

Authors: B. Benson, W. D. Pan, A. Prasad, G. A. Gary, Q. Hu

Abstract: With recent advances in the field of machine learning, the use of deep neural networks for time series forecasting has become more prevalent. The quasi-periodic nature of the solar cycle makes it a good candidate for applying time series forecasting methods. We employ a combination of WaveNet and LSTM neural networks to forecast the sunspot number using the years 1749 to 2019 and total sunspot are… ▽ More With recent advances in the field of machine learning, the use of deep neural networks for time series forecasting has become more prevalent. The quasi-periodic nature of the solar cycle makes it a good candidate for applying time series forecasting methods. We employ a combination of WaveNet and LSTM neural networks to forecast the sunspot number using the years 1749 to 2019 and total sunspot area using the years 1874 to 2019 time series data for the upcoming Solar Cycle 25. Three other models involving the use of LSTMs and 1D ConvNets are also compared with our best model. Our analysis shows that the WaveNet and LSTM model is able to better capture the overall trend and learn the inherent long and short term dependencies in time series data. Using this method we forecast 11 years of monthly averaged data for Solar Cycle 25. Our forecasts show that the upcoming Solar Cycle 25 will have a maximum sunspot number around 106 $\pm$ 19.75 and maximum total sunspot area around 1771 $\pm$ 381.17. This indicates that the cycle would be slightly weaker than Solar Cycle 24. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 17 pages, 16 figures

Journal ref: Solar Physics volume 295, Article number: 65 (2020)

A Solar Magnetic-fan Flaring Arch Heated by Non-thermal Particles and Hot Plasma from an X-ray Jet Eruption

Authors: Kyoung-Sun Lee, Hirohisa Hara, Kyoko Watanabe, Anand D. Joshi, David H. Brooks, Shinsuke Imada, Avijeet Prasad, Phillip Dang, Toshifumi Shimizu, Sabrina L. Savage, Ronald Moore, Navdeep K. Panesar, Jeffrey W. Reep

Abstract: We have investigated an M1.3 limb flare, which develops as a magnetic loop/arch that fans out from an X-ray jet. Using Hinode/EIS, we found that the temperature increases with height to a value of over 10$^{7}$ K at the loop-top during the flare. The measured Doppler velocity (redshifts of 100$-$500 km s$^{-1}$) and the non-thermal velocity ($\geq$100 km s$^{-1}$) from Fe XXIV also increase with l… ▽ More We have investigated an M1.3 limb flare, which develops as a magnetic loop/arch that fans out from an X-ray jet. Using Hinode/EIS, we found that the temperature increases with height to a value of over 10$^{7}$ K at the loop-top during the flare. The measured Doppler velocity (redshifts of 100$-$500 km s$^{-1}$) and the non-thermal velocity ($\geq$100 km s$^{-1}$) from Fe XXIV also increase with loop height. The electron density increases from $0.3\times10^{9}$ cm$^{-3}$ early in the flare rise to $1.3\times10^{9}$ cm$^{-3}$ after the flare peak. The 3-D structure of the loop derived with STEREO/EUVI indicates that the strong redshift in the loop-top region is due to upflowing plasma originating from the jet. Both hard X-ray and soft X-ray emission from RHESSI were only seen as footpoint brightenings during the impulsive phase of the flare, then, soft X-ray emission moves to the loop-top in the decay phase. Based on the temperature and density measurements and theoretical cooling models, the temperature evolution of the flare arch is consistent with impulsive heating during the jet eruption followed by conductive cooling via evaporation and minor prolonged heating in the top of the fan loop. Investigating the magnetic field topology and squashing factor map from SDO/HMI, we conclude that the observed magnetic-fan flaring arch is mostly heated from low atmospheric reconnection accompanying the jet ejection, instead of from reconnection above the arch as expected in the standard flare model. △ Less

Submitted 20 May, 2020; originally announced May 2020.

Comments: 55 pages, 21 Figures, accepted for publication in ApJ

Identification of Pre-flare Processes and Their Possible Role in Driving a Large-scale Flux Rope Eruption with Complex M-class Flare in the Active Region NOAA 12371

Authors: Prabir K. Mitra, Bhuwan Joshi, Avijeet Prasad

Abstract: In this article, we study the origin of precursor flare activity and investigate its role towards triggering the eruption of a flux rope which resulted into a dual-peak M-class flare (SOL2015-06-21T02:36) in the active region NOAA 12371. The flare evolved in two distinct phases with peak flux levels of M2.1 and M2.6 at an interval of $\approx$54 min. The active region exhibited striking moving mag… ▽ More In this article, we study the origin of precursor flare activity and investigate its role towards triggering the eruption of a flux rope which resulted into a dual-peak M-class flare (SOL2015-06-21T02:36) in the active region NOAA 12371. The flare evolved in two distinct phases with peak flux levels of M2.1 and M2.6 at an interval of $\approx$54 min. The active region exhibited striking moving magnetic features (MMFs) along with sunspot rotation. Non-linear force free field (NLFFF) modelling of the active region corona reveals a magnetic flux rope along the polarity inversion line in the trailing sunspot group which is observationally manifested by the co-spatial structures of an active region filament and a hot channel identified in the 304 and 94 Å images, respectively, from the Atmospheric Imaging Assembly (AIA). The active region underwent a prolonged phase of flux enhancement followed by a relatively shorter period of flux cancellation prior to the onset of the flare which led to the build up and activation of the flux rope. Extreme ultra-violet (EUV) images reveal localised and structured pre-flare emission, from the region of MMFs, adjacent to the location of the main flare. Our analysis reveals strong, localised regions of photospheric currents of opposite polarities at the precursor location, thereby making the region susceptible to small-scale magnetic reconnection. Precursor reconnection activity from this location most likely induced a slipping reconnetion towards the northern leg of the hot channel which led to the destabilization of the flux rope. The application of magnetic virial theorem suggests that there was an overall growth of magnetic free energy in the active region during the prolonged pre-flare phase which decayed rapidly after the hot channel eruption and its successful transformation into a halo coronal mass ejection (CME). △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 31 pages, 14 figures, 1 table. This is a pre-print of an article published in Solar Physics. The final authenticated version is available online at: https://doi.org/10.1007/s11207-020-1596-2

Three-dimensional magnetic field structure of a flux emerging region in the solar atmosphere

Authors: Rahul Yadav, J. de la Cruz Rodríguez, C. J. Díaz Baso, Avijeet Prasad, Tine Libbrecht, Carolina Robustini, A. Asensio Ramos

Abstract: We analyze high-resolution spectropolarimetric observations of a flux emerging region (FER) in order to understand its magnetic and kinematic structure. Our spectropolarimetric observations in the He I 1083.0 nm spectral region of a FER are recorded with GRIS at the 1.5 m aperture GREGOR telescope. A Milne-Eddington based inversion code was employed to extract the photospheric information of the S… ▽ More We analyze high-resolution spectropolarimetric observations of a flux emerging region (FER) in order to understand its magnetic and kinematic structure. Our spectropolarimetric observations in the He I 1083.0 nm spectral region of a FER are recorded with GRIS at the 1.5 m aperture GREGOR telescope. A Milne-Eddington based inversion code was employed to extract the photospheric information of the Si I spectral line, whereas the He I triplet line was analyzed with the Hazel inversion code, which takes into account the joint action of the Hanle and the Zeeman effect. The spectropolarimetric analysis of Si I line displays a complex magnetic structure near the vicinity of FER. Moreover, we find supersonic downflows of 40 km/sec appears near the footpoints of loops connecting two pores of opposite polarity, whereas a strong upflows of 22 km/sec appears near the apex of the loops. Furthermore, non-force-free field extrapolations were performed separately at two layers in order to understand the magnetic field topology of the FER. We determine, using extrapolations from the photosphere and the observed chromospheric magnetic field, that the average formation height of the He triplet line is 2 Mm from the solar surface. The reconstructed loops using photospheric extrapolations along an arch filament system have a maximum height of 10.5 Mm from the solar surface with a foot-points separation of 19 Mm, whereas the loops reconstructed using chromospheric extrapolations are around 8.4 Mm high from the solar surface with a foot-point separation of 16 Mm at the chromospheric height. The magnetic topology in the FER suggests the presence of small-scale loops beneath the large loops. Under suitable conditions, due to magnetic reconnection, these loops can trigger various heating events in the vicinity of the FER. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: 14 pages, 13 figure, 1 table, accepted for publication in A&A

Journal ref: A&A 632, A112 (2019)

Godbillon-Vey Helicity and Magnetic Helicity in Magnetohydrodynamics

Authors: G. M. Webb, A. Prasad, S. C. Anco, Q. Hu

Abstract: The Godbillon-Vey invariant occurs in homology theory, and algebraic topology, when conditions for a co-dimension 1, foliation of a 3D manifold are satisfied. The magnetic Godbillon-Vey helicity invariant in magnetohydrodynamics (MHD) is a higher order helicity invariant that occurs for flows, in which the magnetic helicity density… ▽ More The Godbillon-Vey invariant occurs in homology theory, and algebraic topology, when conditions for a co-dimension 1, foliation of a 3D manifold are satisfied. The magnetic Godbillon-Vey helicity invariant in magnetohydrodynamics (MHD) is a higher order helicity invariant that occurs for flows, in which the magnetic helicity density $h_m={\bf A}{\bf\cdot}{\bf B}={\bf A}{\bf\cdot}(\nabla\times{\bf A})=0$, where ${\bf A}$ is the magnetic vector potential and ${\bf B}$ is the magnetic induction. This paper obtains evolution equations for the magnetic Godbillon-Vey field $\boldsymbolη={\bf A}\times{\bf B}/|{\bf A}|^2$ and the Godbillon-Vey helicity density $h_{gv}=\boldsymbolη{\bf\cdot}(\nabla\times{\boldsymbolη})$ in general MHD flows in which either $h_m=0$ or $h_m\neq 0$. A conservation law for $h_{gv}$ occurs in flows for which $h_m=0$. For $h_m\neq 0$ the evolution equation for $h_{gv}$ contains a source term in which $h_m$ is coupled to $h_{gv}$ via the shear tensor of the background flow. The transport equation for $h_{gv}$ also depends on the electric field potential $ψ$, which is related to the gauge for ${\bf A}$, which takes its simplest form for the advected ${\bf A}$ gauge in which $ψ={\bf A\cdot u}$ where ${\bf u}$ is the fluid velocity. An application of the Godbillon-Vey magnetic helicity to nonlinear force-free magnetic fields used in solar physics is investigated. The possible uses of the Godbillon-Vey helicity in zero helicity flows in ideal fluid mechanics, and in zero helicity Lagrangian kinematics of three-dimensional advection are discussed. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: 42 pages, 10 figures, accepted for publication in Journal of Plasma Physics

Effects of dark matter in star formation

Authors: Arun Kenath, S B Gudennavar, A Prasad, C Sivaram

Abstract: The standard model for the formation of structure assumes that there existed small fluctuations in the early universe that grew due to gravitational instability. The origins of these fluctuations are as yet unclear. In this work we propose the role of dark matter in providing the seed for star formation in the early universe. Very recent observations also support the role of dark matter in the for… ▽ More The standard model for the formation of structure assumes that there existed small fluctuations in the early universe that grew due to gravitational instability. The origins of these fluctuations are as yet unclear. In this work we propose the role of dark matter in providing the seed for star formation in the early universe. Very recent observations also support the role of dark matter in the formation of these first stars. With this we set observable constraints on luminosities, temperatures, and lifetimes of these early stars with an admixture of dark matter. △ Less

Submitted 5 February, 2019; originally announced February 2019.

Comments: 16 pages, 5 figures, To appear as Astrophysics and Space Science (2019) 364:24

Successive flux rope eruptions from $δ$-Sunspots region of NOAA 12673 and associated X-class eruptive flares on 2017 September 6

Authors: Prabir K. Mitra, Bhuwan Joshi, Avijeet Prasad, Astrid M. Veronig, R. Bhattacharyya

Abstract: In this paper, we present a multi-wavelength analysis of two X-class solar eruptive flares of classes X2.2 and X9.3 that occurred in the sigmoidal active region NOAA 12673 on 2017 September 6, by combining observations of Atmospheric Imaging Assembly and Helioseismic Magnetic Imager instruments on board the Solar Dynamics Observatory. On the day of the reported activity, the photospheric structure… ▽ More In this paper, we present a multi-wavelength analysis of two X-class solar eruptive flares of classes X2.2 and X9.3 that occurred in the sigmoidal active region NOAA 12673 on 2017 September 6, by combining observations of Atmospheric Imaging Assembly and Helioseismic Magnetic Imager instruments on board the Solar Dynamics Observatory. On the day of the reported activity, the photospheric structure of the active region displayed a very complex network of $δ$-sunspots that gave rise to the formation of a coronal sigmoid observed in the hot EUV channels. Both X-class flares initiated from the core of the sigmoid sequentially within an interval of $\sim$3 hours and progressed as a single "sigmoid--to--arcade" event. Differential emission measure analysis reveals strong heating of plasma at the core of the active region right from the pre-flare phase which further intensified and spatially expanded during each event. The identification of a pre-existing magnetic null by non-force-free-field modeling of the coronal magnetic fields at the location of early flare brightenings and remote faint ribbon-like structures during the pre-flare phase, which were magnetically connected with the core region, provide support for the breakout model of solar eruption. The magnetic extrapolations also reveal flux rope structures prior to both flares which are subsequently supported by the observations of the eruption of hot EUV channels. The second X-class flare diverged from the standard flare scenario in the evolution of two sets of flare ribbons, that are spatially well separated, providing firm evidence of magnetic reconnections at two coronal heights. △ Less

Submitted 31 October, 2018; originally announced October 2018.

Comments: 24 pages, 13 figures. Accepted for publication by The Astrophysical Journal (ApJ)

Magnetohydrodynamic simulation of magnetic null-point reconnections in NOAA AR12192 initiated with an extrapolated non-force-free-field

Authors: A. Prasad, R. Bhattacharyya, Qiang Hu, Sanjay Kumar, Sushree S. Nayak

Abstract: Magnetohydrodynamics of the solar corona is simulated numerically. The simulation is initialized with an extrapolated non-force-free magnetic field using the vector magnetogram of the active region (AR) NOAA 12192 obtained on the solar photosphere. Particularly, we focus on the magnetic reconnections occurring close to a magnetic null-point that resulted in appearance of circular chromospheric fla… ▽ More Magnetohydrodynamics of the solar corona is simulated numerically. The simulation is initialized with an extrapolated non-force-free magnetic field using the vector magnetogram of the active region (AR) NOAA 12192 obtained on the solar photosphere. Particularly, we focus on the magnetic reconnections occurring close to a magnetic null-point that resulted in appearance of circular chromospheric flare ribbons on October 24, 2014 around 21:21 UT, after peak of an X3.1 flare. The extrapolated field lines show the presence of the three-dimensional (3D) null near one of the polarity inversion lines---where the flare was observed. In the subsequent numerical simulation, we find magnetic reconnections occurring near the null point, where the magnetic field lines from the fan-plane of the 3D null form a X-type configuration with underlying arcade field lines. The footpoints of the dome-shaped field lines, inherent to the 3D null, show high gradients of the squashing factor. We find slipping reconnections at these quasi-separatrix layers, which are co-located with the post-flare circular brightening observed at the chromospheric heights. This demonstrates the viability of the initial non-force-free field along with the dynamics it initiates. Moreover, the initial field and its simulated evolution is found to be devoid of any flux rope, which is in congruence with the confined nature of the flare. △ Less

Submitted 2 May, 2018; originally announced May 2018.

Comments: 28 pages, 14 figures, accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, 860:96 (16pp), 2018 June 20

Magnetic field topology from non-force free extrapolation and magnetohydrodynamic simulation of its eventual dynamics

Authors: Sushree S. Nayak, R. Bhattacharyya, A. Prasad, Qiang Hu

Abstract: Magnetic reconnections (MRs) for various magnetic field line (MFL) topologies are believed to be the initiators of various solar eruptive events like flares and coronal mass ejections (CMEs). Consequently, important is a thorough understanding and quantification of the MFL topology and their evolution which leads to MRs. Contemporary standard is to extrapolate the coronal MFLs using equilibrium mo… ▽ More Magnetic reconnections (MRs) for various magnetic field line (MFL) topologies are believed to be the initiators of various solar eruptive events like flares and coronal mass ejections (CMEs). Consequently, important is a thorough understanding and quantification of the MFL topology and their evolution which leads to MRs. Contemporary standard is to extrapolate the coronal MFLs using equilibrium models where the Lorentz force on the coronal plasma is zero everywhere, because either there is no current or the current is parallel to the magnetic field. In tandem, a non-force-free-field (NFFF) extrapolation scheme has evolved and allows for a Lorentz force which is non-zero only at the photosphere but asymptotically vanishes with height. The paper reports magnetohydrodynamic (MHD)- simulations initiated by NFFF extrapolation of the coronal MFLs for the active region NOAA 11158. Interestingly, quasi-separatrix layers (QSLs) which facilitate MRs are detected in the extrapolated MFLs. The AR 11158 is flare producing and, the paper makes an attempt to asses the role of QSLs in the flare onsets. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: Proc. IAU Symp. 340. 2 pages, 3 figures, accepted for publication

Study of magnetic field topology of active region 12192 using an extrapolated non-force-free magnetic field

Authors: A. Prasad, R. Bhattacharyya, Q. Hu, S. S. Nayak, Sanjay Kumar

Abstract: The solar active region (AR) 12192 was one of the most flare productive region of solar cycle 24, which produced many X-class flares; the most energetic being an X3.1 flare on October 24, 2014 at 21:10 UT. Customarily, such events are believed to be triggered by magnetic reconnection in coronal magnetic fields. Here we use the vector magnetograms from solar photosphere, obtained from Heliospheric… ▽ More The solar active region (AR) 12192 was one of the most flare productive region of solar cycle 24, which produced many X-class flares; the most energetic being an X3.1 flare on October 24, 2014 at 21:10 UT. Customarily, such events are believed to be triggered by magnetic reconnection in coronal magnetic fields. Here we use the vector magnetograms from solar photosphere, obtained from Heliospheric Magnetic Imager (HMI) to investigate the magnetic field topology prior to the X3.1 event, and ascertain the conditions that might have caused the flare. To infer the coronal magnetic field, a novel non-force-free field (NFFF) extrapolation technique of the photospheric field is used, which suitably mimics the Lorentz forces present in the photospheric plasma. We also highlight the presence of magnetic null points and quasi-separatrix layers (QSLs) in the magnetic field topology, which are preferred sites for magnetic reconnections and discuss the probable reconnection scenarios. △ Less

Submitted 12 April, 2018; originally announced April 2018.

Comments: Proc. IAU Symp. 340. 2 pages, 4 figures, accepted for publication

Topological and statistical properties of nonlinear force-free fields

Authors: A. Mangalam, A. Prasad

Abstract: We use our semi-analytic solution of the nonlinear force-free field equation to construct three-dimensional magnetic fields that are applicable to the solar corona and study their statistical properties for estimating the degree of braiding exhibited by these fields. We present a new formula for calculating the winding number and compare it with the formula for the crossing number. The comparison… ▽ More We use our semi-analytic solution of the nonlinear force-free field equation to construct three-dimensional magnetic fields that are applicable to the solar corona and study their statistical properties for estimating the degree of braiding exhibited by these fields. We present a new formula for calculating the winding number and compare it with the formula for the crossing number. The comparison is shown for a toy model of two helices and for realistic cases of nonlinear force-free fields; conceptually the formulae are nearly the same but the resulting distributions calculated for a given topology can be different. We also calculate linkages, which are useful topological quantities that are independent measures of the contribution of magnetic braiding to the total free energy and relative helicity of the field. Finally, we derive new analytical bounds for the free energy and relative helicity for the field configurations in terms of the linking number. These bounds will be of utility in estimating the braided energy available for nano-flares or for eruptions. △ Less

Submitted 2 December, 2017; v1 submitted 8 November, 2017; originally announced November 2017.

Comments: 22 pages, 8 figures, 1 table, to appear in a special issue of Advances in Space Research on the "Dynamic Sun"

Magnetohydrodynamic modeling of solar coronal dynamics with an initial non-force-free magnetic field

Authors: A. Prasad, R. Bhattacharyya, Sanjay Kumar

Abstract: The magnetic fields in the solar corona are generally neither force-free nor axisymmetric and have complex dynamics that are difficult to characterize. Here we simulate the topological evolution of solar coronal magnetic field lines (MFLs) using a magnetohydrodynamic model. The simulation is initialized with a non-axisymmetric non-force-free magnetic field that best correlates with the observed ve… ▽ More The magnetic fields in the solar corona are generally neither force-free nor axisymmetric and have complex dynamics that are difficult to characterize. Here we simulate the topological evolution of solar coronal magnetic field lines (MFLs) using a magnetohydrodynamic model. The simulation is initialized with a non-axisymmetric non-force-free magnetic field that best correlates with the observed vector magnetograms of solar active regions (ARs). To focus on these ideas, simulations are performed for the flaring AR 11283 noted for its complexity and well-documented dynamics. The simulated dynamics develops as the initial Lorentz force pushes the plasma and facilitates successive magnetic reconnections at the two X-type null lines present in the initial field. Importantly, the simulation allows for the spontaneous development of mass flow, unique among contemporary works, that preferentially reconnects field lines at one of the X-type null lines. Consequently, a flux rope consisting of low-lying twisted MFLs, which approximately traces the major polarity inversion line, undergoes an asymmetric monotonic rise. The rise is attributed to a reduction in the magnetic tension force at the region overlying the rope, resulting from the reconnection. A monotonic rise of the rope is in conformity with the standard scenario of flares. Importantly, the simulated dynamics leads to bifurcations of the flux rope, which, being akin to the observed filament bifurcation in AR 11283, establishes the appropriateness of the initial field in describing ARs. △ Less

Submitted 7 May, 2017; originally announced May 2017.

Comments: 25 pages, 11 figures, published in The Astrophysical Journal

Journal ref: The Astrophysical Journal,840:37(12pp), 2017

A viable non-axisymmetric non-force-free field to represent solar active regions

Authors: A. Prasad, R. Bhattacharyya

Abstract: A combination of analytical calculations and vectormagnetogram data are utilized to develop a non-axisymmetric non-force-free magnetic field and asses its viability in describing solar active regions. For the purpose, we construct a local spherical shell where a planar surface, tangential to the inner sphere, represents a Cartesian cutout of an active region. The magnetic field defined on the surf… ▽ More A combination of analytical calculations and vectormagnetogram data are utilized to develop a non-axisymmetric non-force-free magnetic field and asses its viability in describing solar active regions. For the purpose, we construct a local spherical shell where a planar surface, tangential to the inner sphere, represents a Cartesian cutout of an active region. The magnetic field defined on the surface is then correlated with magnetograms. The analysis finds the non-axisymmetric non-force-free magnetic field, obtained by a superposition of two linear-force-free fields, correlates reasonably well with magnetograms. △ Less

Submitted 18 November, 2016; originally announced November 2016.

Comments: 5 pages, 3 figures, Published in Physics of Plasmas

Journal ref: Phys. Plasmas 23, 114504 (2016)

A global galactic dynamo with a corona constrained by relative helicity

Authors: A. Prasad, A. Mangalam

Abstract: We present a model for a global axisymmetric turbulent dynamo operating in a galaxy with a corona which treats the supernovae (SNe) and magneto-rotational instability (MRI) driven turbulence parameters under a common formalism. The nonlinear quenching of the dynamo is alleviated by inclusion of small-scale advective and diffusive magnetic helicity fluxes, which allow the gauge invariant magnetic h… ▽ More We present a model for a global axisymmetric turbulent dynamo operating in a galaxy with a corona which treats the supernovae (SNe) and magneto-rotational instability (MRI) driven turbulence parameters under a common formalism. The nonlinear quenching of the dynamo is alleviated by inclusion of small-scale advective and diffusive magnetic helicity fluxes, which allow the gauge invariant magnetic helicity to be transferred outside the disk and consequently build up a corona during the course of dynamo action. The time-dependent dynamo equations are expressed in a separable form and solved through an eigenvector expansion constructed using the steady-state solutions of the dynamo equation. The parametric evolution of the dynamo solution allows us to estimate the final structure of the global magnetic field and the saturated value of the turbulence parameter $α_m$, even before solving the dynamical equations for evolution of magnetic fields in the disk and the corona, along with $α$-quenching. We then solve these equations simultaneously to study the saturation of large-scale magnetic field, its dependence on the small-scale magnetic helicity fluxes and corresponding evolution of the force-free field in the corona. The quadrupolar large-scale magnetic field in the disk is found to reach equipartition strength within a timescale of 1 Gyr. The large-scale magnetic field in the corona obtained is much weaker in strength compared to the field inside the disk and has only a weak impact on the dynamo operation. △ Less

Submitted 9 December, 2015; v1 submitted 11 November, 2015; originally announced November 2015.

Comments: 59 pages, 14 figures, 37 sub-figures, 2 tables and 7 appendices. Accepted for publication in ApJ

Journal ref: 2016 ApJ 817 12

Separable solutions of force-free spheres and applications to solar active regions

Authors: A. Prasad, A. Mangalam, B. Ravindra

Abstract: In this paper, we present a systematic study of the force-free field equation for simple axisymmetric configurations in spherical geometry and apply it to the solar active regions. The condition of separability of solutions in the radial and angular variables leads to two classes of solutions: linear and nonlinear force-free fields. We have studied these linear solutions and extended the nonlinear… ▽ More In this paper, we present a systematic study of the force-free field equation for simple axisymmetric configurations in spherical geometry and apply it to the solar active regions. The condition of separability of solutions in the radial and angular variables leads to two classes of solutions: linear and nonlinear force-free fields. We have studied these linear solutions and extended the nonlinear solutions for the radial power law index to the irreducible rational form $n= p/q$, which is allowed for all cases of odd $p$ and cases of $q>p$ for even $p$, where the poloidal flux $ψ\propto1/r^n$ and field $\mathbf{B}\propto 1/r^{n+2}$). We apply these solutions to simulate photospheric vector magnetograms obtained using the spectropolarimeter on board Hinode. The effectiveness of our search strategy is first demonstrated on test inputs of dipolar, axisymmetric,and non axisymmetric linear force-free fields. Using the best-fit to these magnetograms, we build three-dimensional axisymmetric field configurations and calculate the energy and relative helicity with two independent methods, which are in agreement. We have analyzed five magnetograms for AR 10930 spanning a period of three days during which two X-class flares occurred which allowed us to find the free energy and relative helicity of the active region before and after the flare; our analysis indicates a peak in these quantities before the flare events which is consistent with the results mentioned in literature. We also analyzed single-polarity regions AR 10923 and 10933, which showed very good fits with potential fields. This method can provide useful reconstruction of the nonlinear force-free (NLFF) fields as well as reasonably good input fields for other numerical techniques. △ Less

Submitted 6 May, 2014; v1 submitted 2 April, 2014; originally announced April 2014.

Comments: 39 pages, 11 figures; published, Astrophysical Journal, 2014

Journal ref: Prasad, A., Mangalam, A. & Ravindra B., 2014, ApJ, 786, 81

Models of force-free spheres and applications to solar active regions

Authors: A. Prasad, A. Mangalam

Abstract: Here we present a systematic study of force-free field equation for simple axisymmetric configurations in spherical geometry. The condition of separability of solutions in radial and angular variables leads to two classes of solutions: linear and non-linear force-free fields. We have studied these linear solutions Chandrasekhar (1956) and extended the non-linear solutions given in Low \& Lou (1990… ▽ More Here we present a systematic study of force-free field equation for simple axisymmetric configurations in spherical geometry. The condition of separability of solutions in radial and angular variables leads to two classes of solutions: linear and non-linear force-free fields. We have studied these linear solutions Chandrasekhar (1956) and extended the non-linear solutions given in Low \& Lou (1990) to the irreducible rational form $n= p/q$, which is allowed for all cases of odd $p$ and to cases of $q>p$ for even $p$. We have further calculated their energies and relative helicities for magnetic field configurations in finite and infinite shell geometries. We demonstrate here a method here to be used to fit observed magnetograms as well as to provide good exact input fields for testing other numerical codes used in reconstruction on the non-linear force-free fields. △ Less

Submitted 2 April, 2014; originally announced April 2014.

Comments: 7 pages, 2 figures

Journal ref: Bulletin of Astronomical Society of India, ASI conference series, 2013, vol 10, p 53-59

Visibility graph analysis of solar wind velocity

Authors: Vinita Suyal, Awadhesh Prasad, Harinder P. Singh

Abstract: We analyze in situ measurements of solar wind velocity obtained by Advanced Composition Explorer (ACE) spacecraft and Helios spacecraft during the years 1998-2012 and 1975-1983 respectively. The data belong to mainly solar cycle 23 (1996-2008) and solar cycle 21 (1976-1986) respectively. We use Directed Horizontal Visibility graph (DHVg) algorithm and estimate a graph functional, namely, the degre… ▽ More We analyze in situ measurements of solar wind velocity obtained by Advanced Composition Explorer (ACE) spacecraft and Helios spacecraft during the years 1998-2012 and 1975-1983 respectively. The data belong to mainly solar cycle 23 (1996-2008) and solar cycle 21 (1976-1986) respectively. We use Directed Horizontal Visibility graph (DHVg) algorithm and estimate a graph functional, namely, the degree distance (D) as the Kullback-Leibler divergence (KLD) argument to understand time irreversibility of solar wind time series. We estimate this degree distance irreversibility parameter for these time series at different phases of solar activity cycle. Irreversibility parameter is first established for known dynamical data and then applied for solar wind velocity time series. It is observed that irreversibility in solar wind velocity fluctuations show similar behaviour at 0.3 AU (Helios data) and 1 AU (ACE data). Moreover it changes over the different phases of solar activity cycle. △ Less

Submitted 8 June, 2013; originally announced June 2013.

Comments: 14 pages, 9 figures, 3 tables, accepted for publication in Solar Physics

Symbolic analysis of slow solar wind data using rank order statistics

Authors: Vinita Suyal, Awadhesh Prasad, Harinder P. Singh

Abstract: We analyze time series data of the fluctuations of slow solar wind velocity using rank order statistics. We selected a total of 18 datasets measured by the Helios spacecraft at a distance of 0.32 AU from the sun in the inner heliosphere. The datasets correspond to the years 1975-1982 and cover the end of the solar activity cycle 20 to the middle of the activity cycle 21. We first apply rank order… ▽ More We analyze time series data of the fluctuations of slow solar wind velocity using rank order statistics. We selected a total of 18 datasets measured by the Helios spacecraft at a distance of 0.32 AU from the sun in the inner heliosphere. The datasets correspond to the years 1975-1982 and cover the end of the solar activity cycle 20 to the middle of the activity cycle 21. We first apply rank order statistics to time series from known nonlinear systems and then extend the analysis to the solar wind data. We find that the underlying dynamics governing the solar wind velocity remains almost unchanged during an activity cycle. However, during a solar activity cycle the fluctuations in the slow solar wind time series increase just before the maximum of the activity cycle △ Less

Submitted 22 December, 2011; originally announced December 2011.

Comments: 10 pages, 7 figures, 1 table

Journal ref: Planetary and Space Science, S.N0. 0032-0633, 2011

Hysteresis in a Solar Activity Cycle

Authors: Vinita Suyal, Awadhesh Prasad, Harinder P. Singh

Abstract: We analyze in situ measurements of solar wind velocity obtained by the Advanced Composition Explorer (ACE) spacecraft during the solar activity cycle 23. We calculated a robust complexity measure, the permutation entropy (S) of solar wind time series at different phases of a solar activity cycle. The permutation entropy measure is first tested on the known dynamical data before its application to… ▽ More We analyze in situ measurements of solar wind velocity obtained by the Advanced Composition Explorer (ACE) spacecraft during the solar activity cycle 23. We calculated a robust complexity measure, the permutation entropy (S) of solar wind time series at different phases of a solar activity cycle. The permutation entropy measure is first tested on the known dynamical data before its application to solar wind time series. It is observed that complexity of solar wind velocity fluctuations at 1 AU shows hysteresis phenomenon while following the ascending and descending phases of the activity cycle. This indicates the presence of multistability in the dynamics governing the solar wind velocity over a solar activity cycle. △ Less

Submitted 22 December, 2011; originally announced December 2011.

Comments: 10 pages, 5 figures; Solar Physics, 2011

Nonlinear Time Series Analysis of Sunspot Data

Authors: Vinita Suyal, Awadhesh Prasad, Harinder P. Singh

Abstract: This paper deals with the analysis of sunspot number time series using the Hurst exponent. We use the rescaled range (R/S) analysis to estimate the Hurst exponent for 259-year and 11360-year sunspot data. The results show a varying degree of persistence over shorter and longer time scales corresponding to distinct values of the Hurst exponent. We explain the presence of these multiple Hurst expo… ▽ More This paper deals with the analysis of sunspot number time series using the Hurst exponent. We use the rescaled range (R/S) analysis to estimate the Hurst exponent for 259-year and 11360-year sunspot data. The results show a varying degree of persistence over shorter and longer time scales corresponding to distinct values of the Hurst exponent. We explain the presence of these multiple Hurst exponents by their resemblance to the deterministic chaotic attractors having multiple centers of rotation. △ Less

Submitted 29 September, 2009; v1 submitted 23 September, 2009; originally announced September 2009.

Comments: 10 pages, 6 figures, accepted for publication in Solar Physics, journal style corrections done in this version

Spitzer IRS mapping of the central kpc of Centaurus A

Authors: Alice C. Quillen, Joss Bland-Hawthorn, Joel Green, J. D. Smith, D. Amelia Prasad, Almudena Alonso-Herrero, Mairi H. Brookes, Kieran Cleary, Charles R. Lawrence

Abstract: We report on the results of infrared spectroscopic mapping observations carried out in the nuclear region of Centaurus A (NGC5128). The 500 pc bipolar dust shell discovered by Quillen et al.(2006) is even more clearly seen in the 11.3 micron dust emission feature than previous broad band imaging. The pure rotational lines of molecular hydrogen other than the S(0) line are detected above the dust… ▽ More We report on the results of infrared spectroscopic mapping observations carried out in the nuclear region of Centaurus A (NGC5128). The 500 pc bipolar dust shell discovered by Quillen et al.(2006) is even more clearly seen in the 11.3 micron dust emission feature than previous broad band imaging. The pure rotational lines of molecular hydrogen other than the S(0) line are detected above the dusty disk and associated with the oval dust shell. The molecular hydrogen transitions indicate the presence of warm gas at temperatures 250--720K. The ratio of the dust emission features at 7.7 and 11.3 micron and the ratio of the [NeII](12.8) and 11.3 dust emission feature are lower in the 500 pc dust shell than in the star forming disk. The clearer shell morphology at 11.3 micron, warm molecular hydrogen emission in the shell, and variation in line ratios in the shell compared to those in the disk, confirm spectroscopically that this shell is a separate coherent entity and is unlikely to be a chance superposition of dust filaments. The physical conditions in the shell are most similar to Galactic supernova remnants where blast waves encounter molecular clouds. The lines requiring the highest level of ionization, [NeV] and [OIV], are detected 20--25 arcsec north-east and south-west of the nucleus and at position angles near the radio jet axis. Fine structure line ratios and limits from this region suggest that the medium is low density and illuminated by a hard radiation field at low ionization parameter. These higher S molecular hydrogen pure rotational transitions are also particularly bright in the same region as the [OIV] and [NeV] emission. This suggests that the gas associated with the dust shell has been excited near the jet axis and is part of an ionization cone. △ Less

Submitted 26 October, 2007; originally announced October 2007.

Comments: submitted to MNRAS

Journal ref: 2008, MNRAS, 384, 1469