Nonlinear model to study filamentation instability of circularly polarized dispersive Alfvén waves in solar wind plasmas

This paper examines the small-scale solar wind turbulence driven in view of the Alfvén waves subjected to ponderomotive nonlinearity. Filamentation instability is known to take place for the case of dispersive Alfvén wave (DAW) propagating parallel to the ambient magnetic field. The ponderomotive force associated with DAW is responsible for wave localization and these webs of filaments become more intense and irregular as one proceeds along the spatial domain. The ponderomotive force associated with pump changes with pump parameters giving rise to different evolution patterns. This paper studies in detail the nonlinear evolution of filamentation instability introduced by dispersive Alfven waves (DAWs) which becomes dispersive on account of the finite frequency of DAW i.e., pump frequency is comparable to the ion cyclotron frequency. We have explicitly obtained the perturbation dynamics and then examined the impact of pump magnitude on the driven magnetic turbulence using numerical simulation. The results show steepening at small scales with increasing pump amplitude. The compressibility associated with acoustic fluctuations may explain the variation in spectral scaling of solar wind turbulence as observed by Alexandrova et al. (Astrophys. J. 674:1157, 2008).     

Radio variability of the blazar S5 0716+714: a ∼6.1 year quasi-periodicity

We report a ～6.1 yr quasi-periodicity for the blazar S5 0716+714 using the radio light curves at 4.8, 8 and 14.5 GHz observed by the University of Michigan Radio Astronomical Observatory (UMRAO) from 1981 to 2012, by means of the Jurkevich, discrete correlation function (DCF) and power spectral analysis techniques. There are a general correlation among light curves at different frequencies and a time lag of \(170\pm 10\) days between 4.8 and 14.5 GHz light curves can be confirmed. We also estimate the orbit parameters assuming a binary black hole system, and the magnetic field strength under the jet comoving frame.     

A satellite relative motion model including $$J_2$$ and $$J_3$$J3 via Vinti’s intermediary

Vinti’s potential is revisited for analytical propagation of the main satellite problem, this time in the context of relative motion. A particular version of Vinti’s spheroidal method is chosen that is valid for arbitrary elliptical orbits, encapsulating \(J_2\), \(J_3\), and generally a partial \(J_4\) in an orbit propagation theory without recourse to perturbation methods. As a child of Vinti’s solution, the proposed relative motion model inherits these properties. Furthermore, the problem is solved in oblate spheroidal elements, leading to large regions of validity for the linearization approximation. After offering several enhancements to Vinti’s solution, including boosts in accuracy and removal of some singularities, the proposed model is derived and subsequently reformulated so that Vinti’s solution is piecewise differentiable. While the model is valid for the critical inclination and nonsingular in the element space, singularities remain in the linear transformation from Earth-centered inertial coordinates to spheroidal elements when the eccentricity is zero or for nearly equatorial orbits. The new state transition matrix is evaluated against numerical solutions including the \(J_2\) through \(J_5\) terms for a wide range of chief orbits and separation distances. The solution is also compared with side-by-side simulations of the original Gim–Alfriend state transition matrix, which considers the \(J_2\) perturbation. Code for computing the resulting state transition matrix and associated reference frame and coordinate transformations is provided online as supplementary material.     

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.     

Gravitational instability in isotropic MHD plasma waves

The effect of compressive viscosity, thermal conductivity and radiative heat-loss functions on the gravitational instability of infinitely extended homogeneous MHD plasma has been investigated. By taking in account these parameters we developed the six-order dispersion relation for magnetohydrodynamic (MHD) waves propagating in a homogeneous and isotropic plasma. The general dispersion relation has been developed from set of linearized basic equations and solved analytically to analyse the conditions of instability and instability of self-gravitating plasma embedded in a constant magnetic field. Our result shows that the presence of viscosity and thermal conductivity in a strong magnetic field substantially modifies the fundamental Jeans criterion of gravitational instability.