On the energy integral for first post-Newtonian approximation

The post-Newtonian approximation for general relativity is widely adopted by the geodesy and astronomy communities. It has been successfully exploited for the inclusion of relativistic effects in practically all geodetic applications and techniques such as satellite/lunar laser ranging and very long baseline interferometry. Presently, the levels of accuracy required in geodetic techniques require that reference frames, planetary and satellite orbits and signal propagation be treated within the post-Newtonian regime. For arbitrary scalar W and vector gravitational potentials \(W^j (j=1,2,3)\), we present a novel derivation of the energy associated with a test particle in the post-Newtonian regime. The integral so obtained appears not to have been given previously in the literature and is deduced through algebraic manipulation on seeking a Jacobi-like integral associated with the standard post-Newtonian equations of motion. The new integral is independently verified through a variational formulation using the post-Newtonian metric components and is subsequently verified by numerical integration of the post-Newtonian equations of motion.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

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.