LASCO White-Light Observations of Eruptive Current Sheets Trailing CMEs

Many models of eruptive flares or coronal mass ejections (CMEs) involve formation of a current sheet connecting the ejecting CME flux rope with a magnetic loop arcade. However, there is very limited observational information on the properties and evolution of these structures, hindering progress in understanding eruptive activity from the Sun. In white-light images, narrow coaxial rays trailing the outward-moving CME have been interpreted as current sheets. Here, we undertake the most comprehensive statistical study of CME-rays to date. We use SOHO/LASCO data, which have a higher cadence, larger field of view, and better sensitivity than any previous coronagraph. We compare our results to a previous study of Solar Maximum Mission (SMM) CMEs, in 1984?–?1989, having candidate magnetic disconnection features at the CME base, about half of which were followed by coaxial bright rays. We examine all LASCO CMEs during two periods of minimum and maximum activity in Solar Cycle 23, resulting in many more events, \(\sim130\) CME-rays, than during SMM. Important results include: The occurrence rate of the rays is \(\sim11~\%\) of all CMEs during solar minimum, but decreases to \(\sim7~\%\) at solar maximum; this is most likely related to the more complex coronal background. The rays appear on average 3?–?4 hours after the CME core, and are typically visible for three-fourths of a day. The mean observed current sheet length over the ray lifetime is \(\sim12~R_{\odot}\), with the longest current sheet of \(18.5~R_{\odot}\). The mean CS growth rates are \(188~\mbox{km}\,\mathrm{s}^{-1}\) at minimum and \(324~\mbox{km}\,\mathrm{s}^{-1}\) at maximum. Outward-moving blobs within several rays, which are indicative of reconnection outflows, have average velocities of \(\sim350~\mbox{km}\,\mathrm{s}^{-1}\) with small positive accelerations. A pre-existing streamer is blown out in most of the CME-ray events, but half of these are observed to reform within \(\sim1\) day. The long lifetime and long lengths of the CME-rays challenge our current understanding of the evolution of the magnetic field in the aftermath of CMEs.     

Temporal and Periodic Variations of Sunspot Counts in Flaring and Non-Flaring Active Regions

We analyzed temporal and periodic variations of sunspot counts (SSCs) in flaring (C-, M-, or X-class flares), and non-flaring active regions (ARs) for nearly three solar cycles (1986 through 2016). Our main findings are as follows: i) temporal variations of monthly means of the daily total SSCs in flaring and non-flaring ARs behave differently during a solar cycle and the behavior varies from one cycle to another; during Solar Cycle 23 temporal SSC profiles of non-flaring ARs are wider than those of flaring ARs, while they are almost the same during Solar Cycle 22 and the current Cycle 24. The SSC profiles show a multi-peak structure and the second peak of flaring ARs dominates the current Cycle 24, while the difference between peaks is less pronounced during Solar Cycles 22 and 23. The first and second SSC peaks of non-flaring ARs have comparable magnitude in the current solar cycle, while the first peak is nearly absent in the case of the flaring ARs of the same cycle. ii) Periodic variations observed in the SSCs profiles of flaring and non-flaring ARs derived from the multi-taper method (MTM) spectrum and wavelet scalograms are quite different as well, and they vary from one solar cycle to another. The largest detected period in flaring ARs is \(113\pm 1.6~\mbox{days}\) while we detected much longer periodicities (\(327\pm 13\), \(312 \pm 11\), and \(256\pm 8~\mbox{days}\)) in the non-flaring AR profiles. No meaningful periodicities were detected in the MTM spectrum of flaring ARs exceeding \(55\pm 0.7~\mbox{days}\) during Solar Cycles 22 and 24, while a \(113\pm 1.3~\mbox{days}\) period was detected in flaring ARs of Solar Cycle 23. For the non-flaring ARs the largest detected period was only \(31\pm 0.2~\mbox{days}\) for Cycle 22 and \(72\pm 1.3~\mbox{days}\) for the current Cycle 24, while the largest measured period was \(327\pm 13~\mbox{days}\) during Solar Cycle 23.     

Turbulence and Heating in the Flank and Wake Regions of a Coronal Mass Ejection

As a coronal mass ejection (CME) passes, the flank and wake regions are typically strongly disturbed. Various instruments, including the Large Angle and Spectroscopic Coronagraph (LASCO), the Atmospheric Imaging Assembly (AIA), and the Coronal Multi-channel Polarimeter (CoMP), observed a CME close to the east limb on 26 October 2013. A hot (\({\approx}\,10~\mbox{MK}\)) rising blob was detected on the east limb, with an initial ejection flow speed of \({\approx}\, 330~\mbox{km}\,\mbox{s}^{-1}\). The magnetic structures on both sides and in the wake of the CME were strongly distorted, showing initiation of turbulent motions with Doppler-shift oscillations enhanced from \({\approx}\, \pm 3~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\, \pm 15~\mbox{km}\,\mbox{s}^{-1}\) and effective thermal velocities from \({\approx}\,30~\mbox{km}\,\mbox{s}^{-1}\) to \({\approx}\,60~\mbox{km}\,\mbox{s}^{-1}\), according to the CoMP observations at the Fe?xiii line. The CoMP Doppler-shift maps suggest that the turbulence behaved differently at various heights; it showed clear wave-like torsional oscillations at lower altitudes, which are interpreted as the antiphase oscillation of an alternating red/blue Doppler shift across the strands at the flank. The turbulence seems to appear differently in the channels of different temperatures. Its turnover time was \({\approx}\,1000\) seconds for the Fe 171 Å channel, while it was \({\approx}\,500\) seconds for the Fe 193 Å channel. Mainly horizontal swaying rotations were observed in the Fe 171 Å channel, while more vertical vortices were seen in the Fe 193 Å channel. The differential-emission-measure profiles in the flank and wake regions have two components that evolve differently: the cool component decreased over time, evidently indicating a drop-out of cool materials due to ejection, while the hot component increased dramatically, probably because of the heating process, which is suspected to be a result of magnetic reconnection and turbulence dissipation. These results suggest a new turbulence-heating scenario of the solar corona and solar wind.     

Density variation effect on multi-ions with kinetic Alfven wave around cusp region—a kinetic approach

The kinetic Alfven waves in the presence of homogeneous magnetic field plasma with multi-ions effect are investigated. The dispersion relation and normalised damping rate are derived for low-\(\beta\) plasma using kinetic theory. The effect of density variation of \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions is observed on frequency and damping rate of the wave. The variation of frequency (\(\omega\)) and normalised damping rate (\(\gamma / \varOmega_{H^{ +}} \)) of the wave are studied with respect to \(k_{ \bot} \rho_{j}\), where \(k_{ \bot} \) is the perpendicular wave number, \(\rho_{j}\) is the ion gyroradius and \(j \) denotes \(\text{H}^{+}\), \(\text{He}^{+}\) and \(\text{O}^{+}\) ions. The variation with \(k_{ \bot} \rho_{j}\) is considered over wide range. The parameters appropriate to cusp region are used for the explanation of results. It is found that with hydrogen and helium ions gyration, the frequency of wave is influenced by the density variation of \(\text{H}^{+}\) and \(\text{He}^{+}\) ions but remains insensitive to the change in density of \(\text{O}^{+}\) ions. For oxygen ion gyration, the frequency of wave varies over a short range only for \(\text{O}^{+}\) ion density variation. The wave shows damping at lower altitude due to variation in density of lighter \(\text{H}^{+}\) and \(\text{He}^{+}\) ions whereas at higher altitude only heavy \(\text{O}^{+}\) ions contribute in wave damping. The damping of wave may be due to landau damping or energy transfer from wave to particles. The present study signifies that the both lighter and heavier ions dominate differently to change the characteristics of kinetic Alfven wave and density variation is also an important parameter to understand wave phenomena in cusp region.     

An interacting new holographic dark energy in the framework of fractal cosmology

In this paper, we study an interacting holographic dark energy model in the framework of fractal cosmology. The features of fractal cosmology could pass ultraviolet divergencies and also make a better understanding of the universe in different dimensions. We discuss a fractal FRW universe filled with the dark energy and cold dark matter interacting with each other. It is observed that the Hubble parameter embraces the recent observational range while the deceleration parameter demonstrates an accelerating universe and a behavior similar to \(\Lambda \mbox{CDM}\). Plotting the equation of state shows that it lies in phantom region for interaction mode. We use \(\mathit{Om}\)-diagnostic tool and it shows a phantom behavior of dark energy which is a condition of avoiding the formation of black holes. Finally we execute the StateFinder diagnostic pair and all the trajectories for interacting and non-interacting state of the model meet the fixed point \(\Lambda \mbox{CDM}\) at the start of the evolution. A behavior similar to Chaplygin gas also can be observed in statefinder plane. We find that new holographic dark energy model (NHDE) in fractal cosmology expressed the consistent behavior with recent observational data and can be considered as a model to avoid the formation of black holes in comparison with the main model of NHDE in the simple FRW universe. It has also been observed that for the interaction term varying with matter density, the model generates asymptotic de-Sitter solution. However, if the interaction term varies with energy density, then the model shows Big-Rip singularity. Using our modified CAMB code, we observed that the interacting model suppresses the CMB spectrum at low multipoles \(l<50\) and enhances the acoustic peaks. Based on the observational data sets used in this paper and using Metropolis-Hastings method of MCMC numerical calculation, it seems that the best value with \(1\sigma \) and \(2\sigma \) confidence interval are \(\Omega _{m0}=0.278^{+0.008~+0.010} _{-0.007~-0.009}\), \(H_{0}=69.9^{+0.95~+1.57}_{-0.95~-1.57}\), \(r_{c}=0.08^{+0.02~+0.027}_{-0.002~-0.0027}\), \(\beta =0.496^{+0.005~+0.009} _{-0.005~-0.009}\), \(c= 0.691^{+0.024~+0.039}_{-0.025~-0.037}\) and \(b^{2}=0.035\) according to which we find that the proposed model in the presence of interaction is compatible with the recent observational data.     

Analysis of the exoplanet containing system Kepler-13

We have applied the close binary system analysis program WinFitter, with its physically detailed fitting function, to an intensive study of the complex multiple system Kepler-13 using photometry data from all 13 short cadence quarters downloaded from the NASA Exoplanet Archive (NEA) (http://exoplanetarchive.ipac.caltech.edu). The data-point error of our normalized, phase-sequenced and binned (380 points per bin: 0.00025 phase interval) flux values, at 14 ppm, allows the model’s specification for the mean reference flux level of the system to a precision better than 1 ppm. Our photometrically derived values for the mass and radius of KOI13.01 are \(6.8\pm0.6~\mbox{M}_{\mathrm{J}}\) and \(1.44\pm0.04~\mbox{R}_{\mathrm{J}}\). The star has a radius of \(1.67\pm0.05~\mbox{R}_{\odot}\). Our modelling sets the mean of the orbital inclination \(i\) at \(94.35\pm0.14^{\circ}\), with the star’s mean precession angle \(\phi_{p}\)—\(49.1\pm5.0^{\circ}\) and obliquity \(\theta_{o}\)\(67.9 \pm 3.0^{\circ}\), though there are known ambiguities about the sense in which such angles are measured.Our findings did not confirm secular variation in the transit modelling parameters greater than their full correlated errors, as argued by previous authors, when each quarter’s data was best-fitted with a determinable parameter set without prejudice. However, if we accept that most of the parameters remain the same for each transit, then we could confirm a small but steady diminution in the cosine of the orbital inclination over the 17 quarter timespan. This is accompanied by a slight increase of the star’s precession angle (less negative), but with no significant change in the obliquity of its spin axis. There are suggestions of a history of strong dynamical interaction with a highly distorted planet rotating in a 3:2 resonance with its revolution, together with a tidal lag of \(\sim30~\mbox{deg}\). The mean precessional period is derived to be about 1000 y, but at the present time the motion of the star’s rotation axis appears to be supporting the gravitational torque, rather than providing the balance against it that would be expected over long periods of time.The planet has a small but detectable backwarming effect on the star, which helps to explain the difference in brightness just after transit and just before occultation eclipses. In assessing these findings it is recognized that sources of uncertainty remain, notably with possible inherent micropulsational effects, variations from other components of the multiple star, stellar activity, differential rotation and the neglect of higher order terms (than \(r_{1}^{5}\)) in the fitting function, where \(r_{1}\) is the ratio of the radius of the star to the mean orbital separation of planet and host star.     

GONG p-Mode Parameters Through Two Solar Cycles

We investigate the parameters of global solar p-mode oscillations, namely damping width \(\Gamma\), amplitude \(A\), mean squared velocity \(\langle v^{2}\rangle\), energy \(E\), and energy supply rate \(\mathrm{d}E/\mathrm{d}t\), derived from two solar cycles’ worth (1996?–?2018) of Global Oscillation Network Group (GONG) time series for harmonic degrees \(l=0\,\mbox{--}\,150\). We correct for the effect of fill factor, apparent solar radius, and spurious jumps in the mode amplitudes. We find that the amplitude of the activity-related changes of \(\Gamma\) and \(A\) depends on both frequency and harmonic degree of the modes, with the largest variations of \(\Gamma\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le3300~\upmu\mbox{Hz}\) and \(31\le l \le60\) with a minimum-to-maximum variation of \(26.6\pm0.3\%\) and of \(A\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le 3300~\upmu\mbox{Hz}\) and \(61\le l \le100\) with a minimum-to-maximum variation of \(27.4\pm0.4\%\). The level of correlation between the solar radio flux \(F_{10.7}\) and mode parameters also depends on mode frequency and harmonic degree. As a function of mode frequency, the mode amplitudes are found to follow an asymmetric Voigt profile with \(\nu_{\text{max}}=3073.59\pm0.18~\upmu\mbox{Hz}\). From the mode parameters, we calculate physical mode quantities and average them over specific mode frequency ranges. In this way, we find that the mean squared velocities \(\langle v^{2}\rangle\) and energies \(E\) of p modes are anticorrelated with the level of activity, varying by \(14.7\pm0.3\%\) and \(18.4\pm0.3\%\), respectively, and that the mode energy supply rates show no significant correlation with activity. With this study we expand previously published results on the temporal variation of solar p-mode parameters. Our results will be helpful to future studies of the excitation and damping of p modes, i.e., the interplay between convection, magnetic field, and resonant acoustic oscillations.     

A study of the geomagnetic indices asymmetry based on the interplanetary magnetic field polarities

Data of geomagnetic indices (aa, Kp, Ap, and Dst) recorded near 1 AU over the period 1967–2016, have been studied based on the asymmetry between the interplanetary magnetic field (IMF) directions above and below of the heliospheric current sheet (HCS). Our results led to the following conclusions: (i) Throughout the considered period, 31 random years (62%) showed apparent asymmetries between Toward (\(\mathbf{T}\)) and Away (\(\mathbf{A}\)) polarity days and 19 years (38%) exhibited nearly a symmetrical behavior. The days of \(\mathbf{A}\) polarity predominated over the \(\mathbf{T}\) polarity days by 4.3% during the positive magnetic polarity epoch (1991–1999). While the days of \(\mathbf{T}\) polarity exceeded the days of \(\mathbf{A}\) polarity by 5.8% during the negative magnetic polarity epoch (2001–2012). (ii) Considerable yearly North–South (N–S) asymmetries of geomagnetic indices observed throughout the considered period. (iii) The largest toward dominant peaks for \(aa\) and \(Ap\) indices occurred in 1995 near to minimum of solar activity. Moreover, the most substantial away dominant peaks for \(aa\) and \(Ap\) indices occurred in 2003 (during the descending phase of the solar cycle 23) and in 1991 (near the maximum of solar activity cycle) respectively. (iv) The N–S asymmetry of \(Kp\) index indicated a most significant away dominant peak occurred in 2003. (v) Four of the away dominant peaks of Dst index occurred at the maxima of solar activity in the years 1980, 1990, 2000, and 2013. The largest toward dominant peak occurred in 1991 (at the reversal of IMF polarity). (vi) The geomagnetic indices (aa, Ap, and \(Kp\)) all have northern dominance during positive magnetic polarity epoch (1971–1979), while the asymmetries shifts to the southern solar hemisphere during negative magnetic polarity epoch (2001–2012).     

A Check on the Validity of Magnetic Field Reconstructions

We investigate a method to test whether a numerically computed model coronal magnetic field \({\boldsymbol {B}}\) departs from the divergence-free condition (also known as the solenoidality condition). The test requires a potential field \({\boldsymbol {B}}_{0}\) to be calculated, subject to Neumann boundary conditions, given by the normal components of the model field \({\boldsymbol {B}}\) at the boundaries. The free energy of the model field may be calculated using \(\frac{1}{2\mu _{0}}\int ({\boldsymbol {B}}-{\boldsymbol {B}}_{0})^{2}\mathrm{d}V\), where the integral is over the computational volume of the model field. A second estimate of the free energy is provided by calculating \(\frac{1}{2\mu _{0}}\int {\boldsymbol {B}}^{2}\,\mathrm{d}V-\frac{1}{2\mu _{0}}\int {\boldsymbol {B}}_{0}^{2}\,\mathrm{d}V\). If \({\boldsymbol {B}}\) is divergence free, the two estimates of the free energy should be the same. A difference between the two estimates indicates a departure from \(\nabla \cdot {\boldsymbol {B}}=0\) in the volume. The test is an implementation of a procedure proposed by Moraitis et al. (Solar Phys.289, 4453, 2014) and is a simpler version of the Helmholtz decomposition procedure presented by Valori et al. (Astron. Astrophys.553, A38, 2013). We demonstrate the test in application to previously published nonlinear force-free model fields, and also investigate the influence on the results of the test of a departure from flux balance over the boundaries of the model field. Our results underline the fact that, to make meaningful statements about magnetic free energy in the corona, it is necessary to have model magnetic fields that satisfy the divergence-free condition to a good approximation.     

Effects of Geometries and Substructures of ICMEs on Geomagnetic Storms

To better understand geomagnetic storm generations by ICMEs, we consider the effect of substructures (magnetic cloud, MC, and sheath) and geometries (impact location of flux-rope at the Earth) of the ICMEs. We apply the toroidal magnetic flux-rope model to 59 CDAW CME–ICME pairs to identify their substructures and geometries, and select 20 MC-associated and five sheath-associated storm events. We investigate the relationship between the storm strength indicated by minimum Dst index \((\mathrm{Dst}_{\mathrm{min}})\) and solar wind conditions related to a southward magnetic field. We find that all slopes of linear regression lines for sheath-storm events are steeper (\({\geq}\,1.4\)) than those of the MC-storm events in the relationship between \(\mathrm{Dst}_{\mathrm{min}}\) and solar wind conditions, implying that the efficiency of sheath for the process of geomagnetic storm generations is higher than that of MC. These results suggest that different general solar wind conditions (sheaths have a higher density, dynamic and thermal pressures with a higher fluctuation of the parameters and higher magnetic fields than MCs) have different impact on storm generation. Regarding the geometric encounter of ICMEs, 100% (2/2) of major storms (\(\mathrm{Dst}_{\mathrm{min}} \leq -100~\mbox{nT}\)) occur in the regions at negative \(P_{Y}\) (relative position of the Earth trajectory from the ICME axis in the \(Y\) component of the GSE coordinate) when the eastern flanks of ICMEs encounter the Earth. We find similar statistical trends in solar wind conditions, suggesting that the dependence of geomagnetic storms on 3D ICME–Earth impact geometries is caused by asymmetric distributions of the geoeffective solar wind conditions. For western flank events, 80% (4/5) of the major storms occur in positive \(P_{Y}\) regions, while intense geoeffective solar wind conditions are not located in the positive \(P_{Y}\). These results suggest that the strength of geomagnetic storms depends on ICME–Earth impact geometries as they determine the solar wind conditions at Earth.     

Electron acceleration behind a wavy dipolarization front

In this paper, with the in-situ observations from the Time History of Events and Macroscale Interactions during Substorms (THEMIS) probes we report a wavy dipolarization front (DF) event, where the DF has different magnetic structures and electron distributions at different \(y\) positions in the Geocentric Solar Magnetospheric (GSM) coordinates. At \(y \sim2.1R_{E}\) (\(R_{E}\) is the radius of Earth), the DF has a relatively simple structure, which is similar to that of a conventional DF. At \(y \sim3.0R_{E}\), the DF is revealed to have a multiple DF structure, where the plasma exhibits a vortex flow. Such a wavy DF could be the results of the interchange instability. The different structure of such a wavy DF at different sites has a great effect on electron acceleration. Fermi acceleration can occur at the site of the DF with a simple or multiple DF structure, while betatron acceleration as a local process has the contribution to energetic electrons only at the site of the DF with a simple structure.     

Secular tidal changes in lunar orbit and Earth rotation

Small tidal forces in the Earth–Moon system cause detectable changes in the orbit. Tidal energy dissipation causes secular rates in the lunar mean motion n, semimajor axis a, and eccentricity e. Terrestrial dissipation causes most of the tidal change in n and a, but lunar dissipation decreases eccentricity rate. Terrestrial tidal dissipation also slows the rotation of the Earth and increases obliquity. A tidal acceleration model is used for integration of the lunar orbit. Analysis of lunar laser ranging (LLR) data provides two or three terrestrial and two lunar dissipation parameters. Additional parameters come from geophysical knowledge of terrestrial tides. When those parameters are converted to secular rates for orbit elements, one obtains dn/dt = \(-25.97\pm 0.05 ''/\)cent\(^{2}\), da/dt = 38.30 ± 0.08 mm/year, and di/dt = ?0.5 ± 0.1 \(\upmu \)as/year. Solving for two terrestrial time delays and an extra de/dt from unspecified causes gives \(\sim \) \(3\times 10^{-12}\)/year for the latter; solving for three LLR tidal time delays without the extra de/dt gives a larger phase lag of the N2 tide so that total de/dt = \((1.50 \pm 0.10)\times 10^{-11}\)/year. For total dn/dt, there is \(\le \)1 % difference between geophysical models of average tidal dissipation in oceans and solid Earth and LLR results, and most of that difference comes from diurnal tides. The geophysical model predicts that tidal deceleration of Earth rotation is \(-1316 ''\)/cent\(^{2}\) or 87.5 s/cent\(^{2}\) for UT1-AT, a 2.395 ms/cent increase in the length of day, and an obliquity rate of 9 \(\upmu \)as/year. For evolution during past times of slow recession, the eccentricity rate can be negative.     

Fourier Power Spectra of Solar Noise Storms

We analyzed three noise storms recorded on 200?–?400 MHz Trieste Callisto radio spectra on 2 July 2012, 8 July 2012, and 16 July 2012 by the Fourier method. We divided intervals of the noise storms into five-minute intervals, and in these intervals we computed the mean Fourier spectra as a function of the wave numbers in the frequency and height-scale spaces. We found that these Fourier spectra, where the spectrum from the quiet-activity interval was subtracted, are power-law spectra. The mean power-law index of these spectra in the range \(\ln(k_{z}) = [1.8, 2.9]\) (where \(k_{z}\) is the wave number in the height-scale space) is \(-1.7\pm0.14\), \(-1.6\pm0.14\), and \(-1.5 \pm0.12\) for the 2 July 2012, the 8 July 2012, and the 16 July 2012 noise storms, respectively. It appears that as the number of Type-I bursts in the studied interval increases, the power-law index becomes closer to \(-5/3\); this is known as the Kolmogorov spectral index. The power-law index of the noise storms is very similar to that of the narrowband dm-spikes found in our previous studies. Furthermore, we found a break in the power spectra at \(\ln(k_{z}) \approx2.9\), and the mean power-law index values above this break are \(-2.9\pm0.46\), \(-3.1\pm0.65\), and \(-3.4\pm0.98\), respectively.     

Role of primordial black holes in the direct collapse scenario of supermassive black hole formation at high redshifts

In this paper, we explore the possibility of accreting primordial black holes as the source of heating for the collapsing gas in the context of the direct collapse black hole scenario for the formation of super-massive black holes (SMBHs) at high redshifts, \(z\sim \) 6–7. One of the essential requirements for the direct collapse model to work is to maintain the temperature of the in-falling gas at \(\approx \)10\(^4\) K. We show that even under the existing abundance limits, the primordial black holes of masses \(\gtrsim \)10\(^{-2}M_\odot \), can heat the collapsing gas to an extent that the \(\mathrm{H}_2\) formation is inhibited. The collapsing gas can maintain its temperature at \(10^4\) K till the gas reaches a critical density \(n_{{c}} \,{\approx }\, 10^3~\hbox {cm}^{-3}\), at which the roto-vibrational states of \(\mathrm{H}_2\) approaches local thermodynamic equilibrium and \(\mathrm{H}_2\) cooling becomes inefficient. In the absence of \(\mathrm{H}_2\) cooling, the temperature of the collapsing gas stays at \(\approx \)10\(^4\) K even as it collapses further. We discuss scenarios of subsequent angular momentum removal and the route to find collapse through either a supermassive star or a supermassive disk.     

Addendum to: Strength of the Solar Coronal Magnetic Field – A Comparison of Independent Estimates Using Contemporaneous Radio and White-Light Observations

This addendum uses an alternate fit for the electron density distribution \(N(r)\) (see Figure 1) and estimates the coronal magnetic field using the new model. We find that the estimates of the magnetic field are in close agreement using both the models.

We have fit the \(N(r)\) distribution obtained from STEREO-A/COR1 and SOHO/LASCO-C2 using a fifth-order polynomial (see Figure 1). The expression can be written as

$$\begin{aligned} N_{\text{cor}}(r) &= 1.43 \times 10^{9} r^{-5} - 1.91 \times 10^{9} r^{-4} + 1.07 \times 10^{9} r^{-3} - 2.87 \times 10^{8} r^{-2} \\ &\quad {} + 3.76 \times 10^{7} r^{-1} - 1.91 \times 10^{6} , \end{aligned}$$

(1)

where \(N_{\text{cor}}(r)\) is in units of cm^?3 and \(r\) is in units of \(\mathrm{R}_{\odot}\). The background coronal electron density is enhanced by a factor of 5.5 at 2.63 \(\mathrm{R}_{\odot}\) during the coronal mass ejection (CME). The estimated coronal magnetic field strength (\(B\)) using radio data indicates that \(B(r) \approx(0.51\text{\,--\,}0.48) \pm 0.02\ \mathrm{G}\) in the range \(r \approx2.65\text{\, --\,}2.82\ \mathrm{R}_{\odot}\). The field strengths for STEREO-A/COR1 and SOHO/LASCO-C2 are ≈?0.32 G at \(r \approx 3.11\ \mathrm{R}_{\odot}\) and ≈?0.12 G at \(r \approx 4.40\ \mathrm{R}_{\odot}\), respectively.

     

Fine-resolution analysis of exoplanetary distributions by wavelets: hints of an overshooting iceline accumulation

We investigate 1D exoplanetary distributions using a novel analysis algorithm based on the continuous wavelet transform. The analysis pipeline includes an estimation of the wavelet transform of the probability density function (p.d.f.) without pre-binning, use of optimized wavelets, a rigorous significance testing of the patterns revealed in the p.d.f., and an optimized minimum-noise reconstruction of the p.d.f. via matching pursuit iterations.In the distribution of orbital periods, \(P\), our analysis revealed a narrow subfamily of exoplanets within the broad family of “warm Jupiters”, or massive giants with \(P\gtrsim 300~\mbox{d}\), which are often deemed to be related with the iceline accumulation in a protoplanetary disk. We detected a p.d.f. pattern that represents an upturn followed by an overshooting peak spanning \(P\sim 300\mbox{--}600~\mbox{d}\), right beyond the “period valley”. It is separated from the other planets by p.d.f. concavities from both sides. It has at least 2-sigma significance.In the distribution of planet radii, \(R\), and using the California Kepler Survey sample properly cleaned, we confirm the hints of a bimodality with two peaks about \(R=1.3R_{\oplus }\) and \(R=2.4R_{ \oplus }\), and the “evaporation valley” between them. However, we obtain just a modest significance for this pattern, 2-sigma only at the best. Besides, our follow-up application of the Hartigan and Hartigan dip test for unimodality returns 3 per cent false alarm probability (merely 2.2-sigma significance), contrary to 0.14 per cent (or 3.2-sigma), as claimed by Fulton et al. (2017).     

Statistical approach to triple systems in three-dimensional motion

This paper considers disruption of triple close approaches with low initial velocities and equal masses in the framework of statistical escape theory in a three-dimensional space. The statistical escape theory is based on the assumption that the phase trajectory of a triple system is quasi-ergodic. This system is described by allowing for both energy and angular momentum conservation in the phase space. In this paper, “possibility of escape” is derived with the formation of a binary on the basis of relative distances of the participating bodies. The complete statistical solutions (i.e. the semi-major axis \(a\), the distributions of eccentricity \(e\) of the binary, binary energy \({E}_{{b}}\), escape energy \({E}_{{s}}\) of escaper, and its escape velocity \({v}_{{s}}\)) of the system are derived from the allowable phase space volumes and are in good agreement with the numerical results in the range of perturbing velocities \({v}_{{i}}\)(\(10^{ - 1} \le {v}_{{i}} \le 10^{ - 10}\)) and directions of \({v}_{{i}}(0 \le \alpha _{{i}},\beta _{{i}},\gamma _{{i}} \le \pi )\), \({i} = 1,2,3\). In this paper, the double limit process has been applied to approximate the escape probability. Through this process, it is observed that the perturbing velocity \({v}_{{i}} \to 0^{ +} \), as the product of the semi-major axis \(a\) of the final binary and the square of the escape velocity \({v}_{{s}}\) approach 2/3, i.e. \({a} {v}_{{s}}^{2} \to 2 / 3\), whatever direction of \(\mathbf{v}_{{i}}\) may be.     

A Major Geoeffective CME from NOAA 12371: Initiation,CME–CME Interactions,and Interplanetary Consequences

In this article, we present a multi-wavelength and multi-instrument investigation of a halo coronal mass ejection (CME) from active region NOAA 12371 on 21 June 2015 that led to a major geomagnetic storm of minimum \(\mathrm{Dst} = -204\) nT. The observations from the Atmospheric Imaging Assembly onboard the Solar Dynamics Observatory in the hot EUV channel of 94 Å confirm the CME to be associated with a coronal sigmoid that displayed an intense emission (\(T \sim6\) MK) from its core before the onset of the eruption. Multi-wavelength observations of the source active region suggest tether-cutting reconnection to be the primary triggering mechanism of the flux rope eruption. Interestingly, the flux rope eruption exhibited a two-phase evolution during which the “standard” large-scale flare reconnection process originated two composite M-class flares. The eruption of the flux rope is followed by the coronagraphic observation of a fast, halo CME with linear projected speed of 1366 km?s^?1. The dynamic radio spectrum in the decameter-hectometer frequency range reveals multiple continuum-like enhancements in type II radio emission which imply the interaction of the CME with other preceding slow speed CMEs in the corona within \(\approx10\)?–?\(90~\mbox{R} _{\odot}\). The scenario of CME–CME interaction in the corona and interplanetary medium is further confirmed by the height–time plots of the CMEs occurring during 19?–?21 June. In situ measurements of solar wind magnetic field and plasma parameters at 1 AU exhibit two distinct magnetic clouds, separated by a magnetic hole. Synthesis of near-Sun observations, interplanetary radio emissions, and in situ measurements at 1 AU reveal complex processes of CME–CME interactions right from the source active region to the corona and interplanetary medium that have played a crucial role towards the large enhancement of the geoeffectiveness of the halo CME on 21 June 2015.     

Determination of the Coronal and Interplanetary Magnetic Field Strength and Radial Profiles from Large-Scale Photospheric Magnetic Fields

We propose a new model for the magnetic field at different distances from the Sun during different phases of the solar cycle. The model depends on the observed large-scale non-polar (\({\pm}\, 55^{\circ }\)) photospheric magnetic field and on the magnetic field measured at polar regions from \(55^{\circ }\) N to \(90^{\circ }\) N and from \(55^{\circ }\) S to \(90^{\circ }\) S, which are the visible manifestations of cyclic changes in the toroidal and poloidal components of the global magnetic field of the Sun. The modeled magnetic field is determined as the superposition of the non-polar and polar photospheric magnetic field and considers cycle variations. The agreement between the model predictions and magnetic fields derived from direct in situ measurements at different distances from the Sun, obtained with different methods and at different solar activity phases, is quite satisfactory. From a comparison of the magnetic fields as observed and calculated from the model at 1 AU, we conclude that the model magnetic field variations adequately explain the main features of the interplanetary magnetic field (IMF) radial, \(B_{\mathrm{x}}\), component cycle evolution at Earth’s orbit. The modeled magnetic field averaged over a Carrington rotation (CR) correlates with the IMF \(B_{\mathrm{x}}\) component also averaged over a CR at Earth’s orbit with a coefficient of 0.691, while for seven CR-averaged data, the correlation reaches 0.81. The radial profiles of the modeled magnetic field are compared with those of already existing models. In contrast to existing models, ours provides realistic magnetic-field radial distributions over a wide range of heliospheric distances at different cycle phases, taking into account the cycle variations of the solar toroidal and poloidal magnetic fields. The model is a good approximation of the cycle behavior of the magnetic field in the heliosphere. In addition, the decrease in the non-polar and polar photospheric magnetic fields is shown. Furthermore, the magnetic field during solar cycle maxima and minima decreased from Cycle 21 to Cycle 24. This implies that both the toroidal and poloidal components, and therefore the solar global magnetic field, decreased from Cycle 21 to Cycle 24.