On the use of logistic regression for stellar classification

We are totally immersed in the Big Data era and reliable algorithms and methods for data classification are instrumental for astronomical research. Random Forest and Support Vector Machines algorithms have become popular over the last few years and they are widely used for different stellar classification problems. In this article, we explore an alternative supervised classification method scarcely exploited in astronomy, Logistic Regression, that has been applied successfully in other scientific areas, particularly biostatistics. We have applied this method in order to derive membership probabilities for potential T Tauri star candidates from ultraviolet-infrared colour-colour diagrams.     

The First Solar Seeing Profile Measurement with Two Apertures and Multiple Guide Regions

The sky brightness is a critical parameter for estimating the coronal observation conditions for a solar observatory. As part of a site-survey project in Western China, we measured the sky brightness continuously at the Lijiang Observatory in Yunnan province in 2011. A sky brightness monitor (SBM) was adopted to measure the sky brightness in a region extending from 4.5 to 7.0 apparent solar radii based on the experience of the Daniel K. Inouye Solar Telescope (DKIST) site survey. Every month, the data were collected manually for at least one week. We collected statistics of the sky brightness at four bandpasses located at 450, 530, 890, and 940 nm. The results indicate that aerosol scattering is of great importance for the diurnal variation of the sky brightness. For most of the year, the sky brightness remains under 20 millionths per airmass before local Noon. On average, the sky brightness is less than 20 millionths, which accounts for 40.41% of the total observing time on a clear day. The best observation time is from 9:00 to 13:00 (Beijing time). The Lijiang Observatory is therefore suitable for coronagraphs investigating the structures and dynamics of the corona.     

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.     

An Anisotropic Cosmic-Ray Enhancement Event on 07-June-2015: A Possible Origin

A usual event, called anisotropic cosmic-ray enhancement (ACRE), was observed as a small increase (\({\leq}\,5\%\)) in the count rates of polar neutron monitors during 12?–?19 UT on 07 June 2015. The enhancement was highly anisotropic, as detected only by neutron monitors with asymptotic directions in the southwest quadrant in geocentric solar ecliptic (GSE) coordinates. The estimated rigidity of the corresponding particles is \({\leq}\,1\) GV. No associated detectable increase was found in the space-borne data from the Geostationary Operational Environmental Satellite (GOES), the Energetic and Relativistic Nuclei and Electron (ERNE) on board the Solar and Heliospheric Observatory (SOHO), or the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA) instruments, whose sensitivity was not sufficient to detect the event. No solar energetic particles were present during that time interval. The heliospheric conditions were slightly disturbed, so that the interplanetary magnetic field strength gradually increased during the event, followed by an increase of the solar wind speed after the event. It is proposed that the event was related to a crossing of the boundary layer between two regions with different heliospheric parameters, with a strong gradient of low-rigidity (\({<}\,1\) GV) particles. It was apparently similar to another cosmic-ray enhancement (e.g., on 22 June 2015) that is thought to have been caused by the local anisotropy of Forbush decreases, with the difference that in our case, the interplanetary disturbance was not observed at Earth, but passed by southward for this event.     

Correlations among color,morphology, and luminosity for the luminous red galaxy (LRG) sample of the SDSS data release 5

We have investigated the correlations among color, morphology and luminosity for all LRGs, cut I LRGs, cut II LRGs, and Main galaxies that are also classified as LRGs. It is found that the morphology of LRGs is tightly correlated with luminosity. The rest-frame u-g color of cut I LRGs and cut II LRGs is nearly independent of luminosity, but the color of Main galaxies is correlated with luminosity. For cut I LRGs and Main galaxies, the early type proportion apparently changes with color: at the rest-frame u-g < 1.3 it increases strongly with increasing color, while at the rest-frame u-g > 1.3 it decreases with increasing color. We also notice that the morphology of cut II LRGs is only a weak function of color. Published in Astrofizika, Vol. 50, No. 3, pp. 335–345 (August 2007).     

Lie group variational integrators for the full body problem in orbital mechanics

Equations of motion, referred to as full body models, are developed to describe the dynamics of rigid bodies acting under their mutual gravitational potential. Continuous equations of motion and discrete equations of motion are derived using Hamilton’s principle. These equations are expressed in an inertial frame and in relative coordinates. The discrete equations of motion, referred to as a Lie group variational integrator, provide a geometrically exact and numerically efficient computational method for simulating full body dynamics in orbital mechanics; they are symplectic and momentum preserving, and they exhibit good energy behavior for exponentially long time periods. They are also efficient in only requiring a single evaluation of the gravity forces and moments per time step. The Lie group variational integrator also preserves the group structure without the use of local charts, reprojection, or constraints. Computational results are given for the dynamics of two rigid dumbbell bodies acting under their mutual gravity; these computational results demonstrate the superiority of the Lie group variational integrator compared with integrators that are not symplectic or do not preserve the Lie group structure.     

Morphologies of AGN host galaxies using HST/ACS in the CDFS-GOODS field

Using HST/ACS images in four bands F435W, F606W, F775W and F850LP, we identify optical counterparts to the X-ray sources in the Chandra Deep Field South in the GOODS South field. A detailed study has been made of these sources to study their morphological types. We use methods like decomposition of galaxy luminosity profiles, color maps and visual inspection of 192 galaxies which are identified as possible optical counterparts of Chandra X-ray sources in the CDFS-GOODS field. We find that most moderate luminosity AGN hosts are bulge dominated in the redshift range (z≈0.4–1.3), but not merging/interacting galaxies. This implies probable fueling of the moderate luminosity AGN by mechanisms other than those merger driven.     

Forecasting the Peak Amplitude of the 24th and 25th Sunspot Cycles and Accompanying Geomagnetic Activity

Forecasting solar and geomagnetic levels of activity is essential to help plan missions and to design satellites that will survive for their useful lifetimes. Therefore, amplitudes of the upcoming solar cycles and the geomagnetic activity were forecasted using the neuro-fuzzy approach. Results of this work allow us to draw the following conclusions: Two moderate cycles are estimated to approach their maximum sunspot numbers, 110 and 116 in 2011 and 2021, respectively. However, the predicted geomagnetic activity shown to be in phase with the peak of the 24th sunspot cycle will reach its minimum three years earlier, then it will rise sharply to reach the 25th maximum a year earlier (i.e., 2020). Our analysis of the three-century long sunspot number data-set suggests that the quasi-periodic variation of the long-term evolution of solar activity could explain the irregularity of the short-term cycles seen during the past decades.     

Evidence of Traveling Waves in Filament Threads

High-resolution Hα filtergrams (0.2″) obtained with the Swedish 1-m Solar Telescope resolve numerous very thin, thread-like structures in solar filaments. The threads are believed to represent thin magnetic flux tubes that must be longer than the observable threads. We report on evidence for small-amplitude (1 – 2 km s⁻¹) waves propagating along a number of threads with an average phase velocity of 12 km s⁻¹ and a wavelength of 4″. The oscillatory period of individual threads vary from 3 to 9 minutes. Temporal variation of the Doppler velocities averaged over a small area containing a number of individual threads shows a short-period (3.6 minutes) wave pattern. These short-period oscillations could possibly represent fast modes in accordance with numerical fibril models proposed by Díaz et al. (Astron. Astrophys. 379, 1083, 2001). In some cases, it is clear that the propagating waves are moving in the same direction as the mass flows.     

Bodily tides near spin–orbit resonances

Spin–orbit coupling can be described in two approaches. The first method, known as the “MacDonald torque”, is often combined with a convenient assumption that the quality factor Q is frequency-independent. This makes the method inconsistent, because derivation of the expression for the MacDonald torque tacitly fixes the rheology of the mantle by making Q scale as the inverse tidal frequency. Spin–orbit coupling can be treated also in an approach called “the Darwin torque”. While this theory is general enough to accommodate an arbitrary frequency-dependence of Q, this advantage has not yet been fully exploited in the literature, where Q is often assumed constant or is set to scale as inverse tidal frequency, the latter assertion making the Darwin torque equivalent to a corrected version of the MacDonald torque. However neither a constant nor an inverse-frequency Q reflect the properties of realistic mantles and crusts, because the actual frequency-dependence is more complex. Hence it is necessary to enrich the theory of spin–orbit interaction with the right frequency-dependence. We accomplish this programme for the Darwin-torque-based model near resonances. We derive the frequency-dependence of the tidal torque from the first principles of solid-state mechanics, i.e., from the expression for the mantle’s compliance in the time domain. We also explain that the tidal torque includes not only the customary, secular part, but also an oscillating part. We demonstrate that the lmpq term of the Darwin–Kaula expansion for the tidal torque smoothly passes zero, when the secondary traverses the lmpq resonance (e.g., the principal tidal torque smoothly goes through nil as the secondary crosses the synchronous orbit). Thus, we prepare a foundation for modeling entrapment of a despinning primary into a resonance with its secondary. The roles of the primary and secondary may be played, e.g., by Mercury and the Sun, correspondingly, or by an icy moon and a Jovian planet. We also offer a possible explanation for the “improper” frequency-dependence of the tidal dissipation rate in the Moon, discovered by LLR.