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.

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.

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 \AA\ 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).

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.

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{\eta}={\bf A}\times{\bf B}/|{\bf A}|^2$ and the Godbillon-Vey helicity density $h_{gv}=\boldsymbol{\eta}{\bf\cdot}(\nabla\times{\boldsymbol\eta})$ 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 $\psi$, which is related to the gauge for ${\bf A}$, which takes its simplest form for the advected ${\bf A}$ gauge in which $\psi={\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.

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.

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 $\delta$-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.

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.

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.

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.

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.

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.

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.

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 $\alpha_m$, even before solving the dynamical equations for evolution of magnetic fields in the disk and the corona, along with $\alpha$-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.

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 $\psi\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.

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.

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.

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

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.

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.