Solving Kepler’s equation via Smale’s \alpha -theory

We obtain an approximate solution $\tilde{E}=\tilde{E}(e,M)$ of Kepler’s equation $E-e\sin (E)=M$ for any $e\in [0,1)$ and $M\in [0,\pi ]$ . Our solution is guaranteed, via Smale’s $\alpha $ -theory, to converge to the actual solution $E$ through Newton’s method at quadratic speed, i.e. the $n$ -th iteration produces a value $E_n$ such that $|E_n-E|\le (\frac{1}{2})^{2^n-1}|\tilde{E}-E|$ . The formula provided for $\tilde{E}$ is a piecewise rational function with conditions defined by polynomial inequalities, except for a small region near $e=1$ and $M=0$ , where a single cubic root is used. We also show that the root operation is unavoidable, by proving that no approximate solution can be computed in the entire region $[0,1)\times [0,\pi ]$ if only rational functions are allowed in each branch.     

The broad band spectral index of TeV blazars detected by Fermi LAT

Using γ-ray data detected by Fermi Large Area Telescope (LAT) and multi-wave band data for 35 TeV blazars sample, we have studied the possible correlations between different broad band spectral indices ( $\alpha_{\rm r.ir}$ , $\alpha_{\rm{r.o}}$ , $\alpha_{\rm r.x}$ , $\alpha_{\rm r.\gamma}$ , $\alpha_{\rm{ir.o}}$ , $\alpha_{\rm ir.x}$ , $\alpha_{\rm ir.\gamma}$ , $\alpha_{\rm o.x}$ , $\alpha_{\rm o.\gamma}$ , $\alpha_{\rm r.x}$ , $\alpha_{\rm x.\gamma}$ ) in all states (average/high/low). Our results are as follows: (1) For our TeV blazars sample, the strong positive correlations were found between $\alpha_{\rm r.ir}$ and $\alpha_{\rm{r.o}}$ , between $\alpha_{\rm r.ir}$ and $\alpha_{\rm r.x}$ , between $\alpha_{\rm r.ir}$ and $\alpha_{\rm r.\gamma}$ in all states (average/high/low); (2) For our TeV blazars sample, the strong anti-correlations were found between $\alpha_{\rm r.ir}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm ir.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm o.\gamma}$ , between $\alpha_{\rm{r.o}}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\mathrm{ir.o}}$ and $\alpha_{\rm o.\gamma}$ , between $\alpha_{\rm r.x}$ and $\alpha_{\rm x.\gamma}$ , between $\alpha_{\rm ir.x}$ and $\alpha_{\rm x.\gamma}$ in all states (average/high/low). The results suggest that the synchrotron self-Compton radiation (SSC) is the main mechanism of high energy γ-ray emission and the inverse Compton scattering of circum-nuclear dust is likely to be a important complementary mechanism for TeV blazars. Our results also show that the possible correlations vary from state to state in the same pair of indices, Which suggest that there may exist differences in the emitting process and in the location of the emitting region for different states.     

A Statistical Study on Photospheric Magnetic Nonpotentiality of Active Regions and Its Relationship with Flares During Solar Cycles 22?–?23

A statistical study is carried out on the photospheric magnetic nonpotentiality in solar active regions and its relationship with associated flares. We select 2173 photospheric vector magnetograms from 1106 active regions observed by the Solar Magnetic Field Telescope at Huairou Solar Observing Station, National Astronomical Observatories of China, in the period of 1988??C?2008, which covers most of the 22nd and 23rd solar cycles. We have computed the mean planar magnetic shear angle ( $\overline{\Delta\phi}$ ), mean shear angle of the vector magnetic field ( $\overline{\Delta\psi}$ ), mean absolute vertical current density ( $\overline{|J_{z}|}$ ), mean absolute current helicity density ( $\overline{|h_{\mathrm{c}}|}$ ), absolute twist parameter (|?? _av|), mean free magnetic energy density ( $\overline{\rho_{\mathrm{free}}}$ ), effective distance of the longitudinal magnetic field (d _E), and modified effective distance (d _Em) of each photospheric vector magnetogram. Parameters $\overline{|h_{\mathrm{c}}|}$ , $\overline{\rho_{\mathrm{free}}}$ , and d _Em show higher correlations with the evolution of the solar cycle. The Pearson linear correlation coefficients between these three parameters and the yearly mean sunspot number are all larger than 0.59. Parameters $\overline {\Delta\phi}$ , $\overline{\Delta\psi}$ , $\overline{|J_{z}|}$ , |?? _av|, and d _E show only weak correlations with the solar cycle, though the nonpotentiality and the complexity of active regions are greater in the activity maximum periods than in the minimum periods. All of the eight parameters show positive correlations with the flare productivity of active regions, and the combination of different nonpotentiality parameters may be effective in predicting the flaring probability of active regions.     

Coupling between corotation and Lindblad resonances in the presence of secular precession rates

We investigate the dynamics of two satellites with masses $\mu _s$ and $\mu '_s$ orbiting a massive central planet in a common plane, near a first order mean motion resonance $m+1{:}m$ (m integer). We consider only the resonant terms of first order in eccentricity in the disturbing potential of the satellites, plus the secular terms causing the orbital apsidal precessions. We obtain a two-degrees-of-freedom system, associated with the two critical resonant angles $\phi = (m+1)\lambda ' -m\lambda - \varpi $ and $\phi '= (m+1)\lambda ' -m\lambda - \varpi '$ , where $\lambda $ and $\varpi $ are the mean longitude and longitude of periapsis of $\mu _s$ , respectively, and where the primed quantities apply to $\mu '_s$ . We consider the special case where $\mu _s \rightarrow 0$ (restricted problem). The symmetry between the two angles $\phi $ and $\phi '$ is then broken, leading to two different kinds of resonances, classically referred to as corotation eccentric resonance (CER) and Lindblad eccentric Resonance (LER), respectively. We write the four reduced equations of motion near the CER and LER, that form what we call the CoraLin model. This model depends upon only two dimensionless parameters that control the dynamics of the system: the distance $D$ between the CER and LER, and a forcing parameter $\epsilon _L$ that includes both the mass and the orbital eccentricity of the disturbing satellite. Three regimes are found: for $D=0$ the system is integrable, for $D$ of order unity, it exhibits prominent chaotic regions, while for $D$ large compared to 2, the behavior of the system is regular and can be qualitatively described using simple adiabatic invariant arguments. We apply this model to three recently discovered small Saturnian satellites dynamically linked to Mimas through first order mean motion resonances: Aegaeon, Methone and Anthe. Poincaré surfaces of section reveal the dynamical structure of each orbit, and their proximity to chaotic regions. This work may be useful to explore various scenarii of resonant capture for those satellites.     

High order normal form stability estimates for co-orbital motion

We present estimates of the size of the analytic domain of stability for co-orbital motions obtained by a high order normal form in the framework of the elliptic restricted three body problem. As a demonstration example, we consider the motion of a Trojan body in an extrasolar planetary system with a giant planet of mass parameter $\mu =0.005$ μ = 0.005 and eccentricity $e^{\prime }=0.1$ e ′ = 0.1 . The analysis contains three basic steps: (i) derivation of an accurate expansion of the Hamiltonian, (ii) computation of the normal form up to an optimal order (in the Nekhoroshev sense), and (iii) computation of the optimal size of the remainder at various values of the action integrals (proper elements) of motion. We explain our choice of variables as well as the method used to expand the Hamiltonian so as to ensure a precise model. We then compute the normal form up to the normalisation order $r=50$ r = 50 by use of a computer-algebraic program. We finally estimate the size $||R||$ | | R | | of the remainder as a function of the normalization order, and compute the optimal normalization order at which the remainder becomes minimum. It is found that the optimal value $\log (||R_{opt}||)$ log ( | | R o p t | | ) can serve in order to construct a stability map for the domain of co-orbital motion using only series. This is compared to the stability map found by a purely numerical approach based on chaotic indicators.     

On the Maximum Separation of Visual Binaries

In this paper, an efficient algorithm is established for computing the maximum (minimum) angular separation ρ _max(ρ _min), the corresponding apparent position angles ( $\theta|_{\rho_{\rm max}}$ , $\theta|_{\rho_{\rm min}}$ ) and the individual masses of visual binary systems. The algorithm uses Reed’s formulae (1984) for the masses, and a technique of one-dimensional unconstrained minimization, together with the solution of Kepler’s equation for $(\rho_{\rm max}, \theta|_{\rho_{\rm max}})$ and $(\rho_{\rm min}, \theta|_{\rho_{\rm min}})$ . Iterative schemes of quadratic coverage up to any positive integer order are developed for the solution of Kepler’s equation. A sample of 110 systems is selected from the Sixth Catalog of Orbits (Hartkopf et al. 2001). Numerical studies are included and some important results are as follows: (1) there is no dependence between ρ _max and the spectral type and (2) a minor modification of Giannuzzi’s (1989) formula for the upper limits of ρ _max functions of spectral type of the primary.     

Continua of central configurations with a negative mass in the n-body problem

The number of equivalence classes of central configurations of $n \le 4$ bodies of positive mass is known to be finite, but it remains to be shown if this is true for $n \ge 5$ . By allowing one mass to be negative, Gareth Roberts constructed a continuum of inequivalent planar central configurations of $n = 5$ bodies. We reinterpret Roberts’ example and generalize the construction of his continuum to produce a family of continua of central configurations, each with a single negative mass. These new continua exist in even dimensional spaces $\mathbb R ^k$ for $k \ge 4$ .     

On the “blue sky catastrophe” termination in the restricted four-body problem

The restricted three-body problem (R3BP) possesses the property that some classes of doubly asymptotic (i.e., homoclinic or heteroclinic) orbits are limit members of families of periodic orbits, this phenomenon has been known as the “blue sky catastrophe” termination principle. A similar case occurs in the restricted four body problem for the collinear equilibrium point $L_{2}$ L 2 . In the restricted four body problem with primaries in a triangle relative equilibrium, we show that the same phenomenon observed in the R3BP occurs. We prove that there exists a critical value of the mass parameter $\mu _{b}$ μ b such that for $\mu =\mu _{b}$ μ = μ b a Hamiltonian Hopf bifurcation takes place. Moreover we show that for $\mu >\mu _{b}$ μ > μ b the stable and unstable manifolds of $L_{2}$ L 2 intersect transversally and the spectrum corresponds to a complex saddle. This proves that Henrard’s theorem applies at least for $\mu $ μ close to $\mu _{b}$ μ b . In particular there exists a family of periodic orbits having the homoclinic orbit as a limit.     

Fundamental Emission of Type III Bursts Produced in Non-Maxwellian Coronal Plasmas with Kappa-Distributed Background Particles

Detailed simulations based on quasi-linear theory are presented for fundamental ( $f_{\rm p}$ ) emission of type III bursts produced in non-Maxwellian, suprathermal, background coronal plasma by injection of energetic electrons during flares with a power-law or Maxwellian velocity distribution, where $f_{\rm p}$ is the electron plasma frequency. The background plasma is assumed to have a kappa (κ) distribution, as inferred from solar wind data and proposed by theories for the corona and solar wind. The predicted type III beam speeds, Langmuir wave levels, and the drift rate and flux of $f_{\rm p}$ emission are strongly sensitive to the presence of suprathermal background electrons in the corona. The simulations show the following results. i) Fast beams with speeds $v_{\rm b}>0.5c$ are produced for coronal background electrons with small κ (κ?5) by injected electrons with power-law spectra. ii) Moderately fast beams with $v_{\rm b} \approx0.3\,\mbox{--}\,0.5c$ are generated in coronal plasma with κ?8 by injections of power-law or Maxwellian electrons. iii) Slow beams with $v_{\rm b}<0.3c$ are produced for coronal background electrons with large κ (κ>8), including the asymptotic limit κ→∞ where the electrons are Maxwellian, for both power-law and Maxwellian injections. The observation of fast type III beams (with $v_{\rm b}>0.5c$ ) thus suggests that these beams are produced in coronal regions where the background electron distribution has small κ by injected electrons with power-law spectra, at least when such beams are observed. The simulations, from the viewpoint of type III bursts, thus support: i) the presence, at least sometimes, of suprathermal background electrons in the corona and the associated mechanisms for coronal heating and solar wind acceleration; ii) power-law spectra for injected energetic electrons, consistent with observations of such electrons in situ and of X-ray emission.     

On the connection between the Nekhoroshev theorem and Arnold diffusion

The analytical techniques of the Nekhoroshev theorem are used to provide estimates on the coefficient of Arnold diffusion along a particular resonance in the Hamiltonian model of Froeschlé et al. (Science 289:2108–2110, 2000). A resonant normal form is constructed by a computer program and the size of its remainder ||R _opt || at the optimal order of normalization is calculated as a function of the small parameter ${\epsilon}$ . We find that the diffusion coefficient scales as ${D \propto ||R_{opt}||^3}$ , while the size of the optimal remainder scales as ${||R_{opt}|| \propto {\rm exp}(1/\epsilon^{0.21})}$ in the range ${10^{-4} \leq \epsilon \leq 10^{-2}}$ . A comparison is made with the numerical results of Lega et al. (Physica D 182:179–187, 2003) in the same model.     

An analytical proof on certain determinants connected with the collinear central configurations in the n-body problem

In this paper we give a short analytical proof of the inequalities proved by Albouy–Moeckel through computer algebra, in the cases $n=5$ and $n=6$ . These inequalities guarantee that, in the $n$ -body problem, the family of mass vectors making a given collinear configuration a central configuration is 2-dimensional. The induction techniques here may be used to prove the inequalities for general $n$ with more subtle estimation but currently the inequalities still remains unproved for $n\ge 7$ .     

Integral wide-field spectroscopy in astronomy: the Imaging FTS solution

Long-slit grating spectrometers in scanning mode and Fabry–Perot interferometers as tunable filters are commonly used to perform integral wide-field spectroscopy on extended astrophysical objects as HII regions and nearby galaxies. The goal of this paper is to demonstrate, by comparison, through a thorough review of the imaging Fourier transform spectrometer (IFTS) properties, that this instrument represents another interesting solution. After a brief recall of the performances, regarding FOV and spectral resolution, of the grating spectrometer, without and with integral field units (IFU), and of the imaging Fabry–Perot, it is demonstrated that for an IFTS the product of the maximum resolution R by the entrance beam étendue U is equal to $2.6\,N\times S_I$ with $N\,\times \,N$ the number of pixels of the detector array and S $_I$ the area of the interferometer beamsplitter. As a consequence, the IFTS offers the most flexible choice of field size and spectral resolution, up to high values for both parameters. It also presents on a wide field an important multichannel advantage in comparison to integral field grating spectrometers, even with multiple IFUs. To complete, the few astronomical IFTSs, built behind ground-based telescopes and in space, for the visible range up to the sub-millimetric domain, are presented. Through two wide-field IFTS projects, one in the visible, the other one in the mid-infrared, the question is addressed of the practical FOV and resolution limits, set by the optical design of the instrument, which can be achieved. Within the 0.3 to $\sim $ 2.5 $\upmu$ m domain, a Michelson interferometer with wide-field diopric collimators provides the easiest solution. This design is illustrated by a $11^{\prime}\times 11^{\prime}$ -field IFTS in the 0.35–0.90 $\upmu$ m range around an off-axis interferometer, called SITELLE, proposed for the 3.6-m CFH Telescope. At longer wavelengths, an all-mirror optics is required, as studied for a spaceborne IFTS, H2EX, for the 8–29 $\upmu$ m range, a $20^{\prime} \times 20^{\prime}$ field, and a high resolution of $\simeq 3\times 10^4$ at 10 $\upmu$ m. To comply with these characteristics, the interferometer is designed with cat’s eye retroreflectors. In the same domain and up to the far infrared, if the instrument aims only at a low spectral resolution (few thousands) and a smaller field (few arcmins $^2$ ), roof-top or corner cube mirrors, as for the IFTS SPIRE on the Herschel space telescope, are usable. At last, perspectives are opened, behind an ELT in the visible and the near infrared with the SITELLE optical combination, in the 2–5 $\upmu$ m on the Antarctic plateau or in space up to longer wavelengths, with the H2EX design, to provide the missing capability of global high spectral resolution studies of extended sources, from comets to distant galaxy clusters.     

The Chandra X-ray galaxy clusters at z<1.4: constraints on the evolution of L X −T−M g relations

We analyzed the luminosity-temperature-mass of gas (L _X ?T?M _g ) relations for a sample of 21 Chandra galaxy clusters. We used the standard approach (β?model) to evaluate these relations for our sample that differs from other catalogues since it considers galaxy clusters at higher redshifts (0.4<z<1.4). We assumed power-law relations in the form $L_{X} \sim(1 +z)^{A_{L_{X}T}} T^{\beta_{L_{X}T}}$ , $M_{g} \sim(1 + z)^{A_{M_{g}T}} T^{\beta_{M_{g}T}}$ , and $M_{g} \sim(1 + z)^{A_{M_{g}L_{X}}} L^{\beta_{M_{g}L_{X}}}$ . We obtained the following fitting parameters with 68 % confidence level: $A_{L_{X}T} = 1.50 \pm0.23$ , $\beta_{L_{X}T} = 2.55 \pm0.07$ ; $A_{M_{g}T} = -0.58 \pm0.13$ and $\beta_{M_{g}T} = 1.77 \pm0.16$ ; $A_{M_{g}L_{X}} \approx-1.86 \pm0.34$ and $\beta_{M_{g}L_{X}} = 0.73 \pm0.15$ , respectively. We found that the evolution of the M _g ?T relation is small, while the M _g ?L _X relation is strong for the cosmological parameters Ω _m =0.27 and Ω _Λ =0.73. In overall, the clusters at high-z have stronger dependencies between L _X ?T?M _g correlations, than those for clusters at low-z. For most of galaxy clusters (first of all, from MACS and RCS surveys) these results are obtained for the first time.     

The planetary spin and rotation period: a modern approach

Using a new approach, we have obtained a formula for calculating the rotation period and radius of planets. In the ordinary gravitomagnetism the gravitational spin (S) orbit (L) coupling, $\vec{L}\cdot\vec{S}\propto L^{2}$ , while our model predicts that $\vec{L}\cdot\vec{S}\propto\frac{m}{M}L^{2}$ , where M and m are the central and orbiting masses, respectively. Hence, planets during their evolution exchange L and S until they reach a final stability at which MS∝mL, or $S\propto\frac{m^{2}}{v}$ , where v is the orbital velocity of the planet. Rotational properties of our planetary system and exoplanets are in agreement with our predictions. The radius (R) and rotational period (D) of tidally locked planet at a distance a from its star, are related by, $D^{2}\propto\sqrt{\frac{M}{m^{3}}}R^{3}$ and that $R\propto\sqrt{\frac {m}{M}}a$ .     

Secular dynamics of a lunar orbiter: a global exploration using Prony’s frequency analysis

We study the secular dynamics of lunar orbiters, in the framework of high-degree gravity models. To achieve a global view of the dynamics, we apply a frequency analysis (FA) technique which is based on Prony’s method. This allows for an extensive exploration of the eccentricity ( $e$ )—inclination ( $i$ ) space, based on short-term integrations ( $\sim $ 8 months) over relatively high-resolution grids of initial conditions. Different gravity models are considered: 3rd, 7th and 10th degree in the spherical harmonics expansion, with the main perturbations from the Earth being added. Since the dynamics is mostly regular, each orbit is characterised by a few parameters, whose values are given by the spectral decomposition of the orbital elements time series. The resulting frequency and amplitude maps in ( $e_0,i_0$ ) are used to identify the dominant perturbations and deduce the “minimum complexity” model necessary to capture the essential features of the long-term dynamics. We find that the 7th degree zonal harmonic ( $J_7$ term) is of profound importance at low altitudes as, depending on the initial secular phases, it can lead to collision with the Moon’s surface within a few months. The 3rd-degree non-axisymmetric terms are enough to describe the deviations from the 1 degree-of-freedom zonal problem; their main effect is to modify the equilibrium value of the argument of periselenium, $\omega $ , with respect to the “frozen” solution ( $\omega =\pm 90^{\circ }, \forall \Omega $ , where $\Omega $ is the nodal longitude). Finally, we show that using FA on a fine grid of initial conditions, set around a suitably chosen ‘first guess’, one can compute an accurate approximation of the initial conditions of a periodic orbit.     

Observational constraints on G-corrected holographic dark energy using a Markov chain Monte Carlo method

We constrain holographic dark energy (HDE) with time varying gravitational coupling constant in the framework of the modified Friedmann equations using cosmological data from type Ia supernovae, baryon acoustic oscillations, cosmic microwave background radiation and X-ray gas mass fraction. Applying a Markov Chain Monte Carlo (MCMC) simulation, we obtain the best fit values of the model and cosmological parameters within 1σ confidence level (CL) in a flat universe as: $\varOmega_{b}h^{2}=0.0222^{+0.0018}_{-0.0013}$ , $\varOmega_{c}h^{2}=0.1121^{+0.0110}_{-0.0079}$ , $\alpha_{G}\equiv \dot{G}/(HG) =0.1647^{+0.3547}_{-0.2971}$ and the HDE constant $c=0.9322^{+0.4569}_{-0.5447}$ . Using the best fit values, the equation of state of the dark component at the present time w _d0 at 1σ CL can cross the phantom boundary w=?1.     

Stability of the triangular Lagrange points beyond Gascheau��s value

We examine the stability of the triangular Lagrange points L ₄ and L ₅ for secondary masses larger than the Gascheau??s value ${\mu_{\rm G}= (1-\sqrt{23/27}/2)= 0.0385208\ldots}$ (also known as the Routh value) in the restricted, planar circular three-body problem. Above that limit the triangular Lagrange points are linearly unstable. Here we show that between??? _G and ${\mu \approx 0.039}$ , the L ₄ and L ₅ points are globally stable in the sense that a particle released at those points at zero velocity (in the corotating frame) remains in the vicinity of those points for an indefinite time. We also show that there exists a family of stable periodic orbits surrounding L ₄ or L ₅ for ${\mu \ge \mu_G}$ . We show that??? _G is actually the first value of a series ${\mu_0 (=\mu_G), \mu_1,\ldots, \mu_i,\ldots}$ corresponding to successive period doublings of the orbits, which exhibit ${1, 2, \ldots, 2^i,\ldots}$ cycles around L ₄ or L ₅. Those orbits follow a Feigenbaum cascade leading to disappearance into chaos at a value ${\mu_\infty = 0.0463004\ldots}$ which generalizes Gascheau??s work.     

Stability of ion-acoustic waves in a pair-ion plasma with a third species of ions: application to cometary plasmas

A popular model of a cometary plasma is hydrogen (H⁺) with positively charged oxygen (O⁺) as a heavier ion component. However, the discovery of negatively charged oxygen (O^?) ions enables one to model a cometary plasma as a pair-ion plasma (of O⁺ and O^?) with hydrogen as a third ion constituent. We have, therefore, studied the stability of the ion-acoustic wave in such a pair-ion plasma with hydrogen and electrons streaming with velocities $V_{d\mathrm{H}^{+}}$ and V _de , respectively, relative to the oxygen ions. We find the calculated frequency of the ion-acoustic wave with this model to be in good agreement with the observed frequencies. The ion-acoustic wave can also be driven unstable by the streaming velocity of the hydrogen ions. The growth rate increases with increasing hydrogen density $n_{\mathrm{H}^{+}}$ , and streaming velocities $V_{d\mathrm{H}^{+}}$ and V _de . It, however, decreases with increasing oxygen ion densities $n_{\mathrm{O}^{+}}$ and $n_{\mathrm{O}^{-}}$ .     

Non-minimal derivatively coupled quintessence in the Palatini formalism

The quintessence dark energy model with a kinetic coupling to gravity within the Palatini formalism is studied in this paper. Two different coupling forms: $\hat{R}\partial^{\mu}\phi\partial_{\mu}\phi$ and $\hat {R}_{\mu\nu}\partial^{\mu}\phi\partial^{\nu}\phi$ are analyzed, respectively. We find that both the model with the $\hat{R}\partial^{\mu}\phi\partial_{\mu}\phi$ coupling and the one with the $\hat{R}_{\mu\nu}\partial^{\mu}\phi\partial^{\nu}\phi$ coupling can realize the phantom divide line crossing from phantom to quintessence at late time for its effective equation-of-state. Furthermore, the former can behave like phantom. These features are different from those found in the $\hat {R}\phi^{2}$ coupling case.