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.     

Fundamental Plane of FSRQs from SDSS DR5 Quasar Catalogue

In this work, we explore a sample of 362 flat-spectrum radio quasars (FSRQs) to investigate the jet formation. We find that the fundamental plane for our FSRQs can be expressed as $L_{\rm 5~GHz}\propto M_{\rm bh}^{-0.19}L_{\rm 2~keV}^{1.08}$ . We also find that the 5?GHz luminosities are tightly related to both black hole mass and Eddington ratio, which is established as $L_{\rm 5~GHz}\propto M_{\rm bh}^{0.67}(L_{\rm bol}/ L_{\rm EDD})^{1.32}$ .     

The effect of wind on magnetized advection dominated accretion flows

We study the effects of winds on advection dominated accretion flows in the presence of a global magnetic field under a self-similar treatment. The disk gas is assumed to be isothermal. For a steady state structure of such accretion flows a set of self similar solutions are presented. We consider the wind in a general magnetic field with three components (r,φ,z) in advection-dominated accretion flows. The mass-accretion rate $\dot{M}$ decreases with radius r as $\dot{M}\propto r^{s+1/2}$ , where s is an arbitrary constant. We will see, by increasing the wind parameter s, radial and rotational velocities increase.     

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$ .     

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.     

Binary pulsars in magnetic field versus spin period diagram

We analyzed 186 binary pulsars (BPSRs) in the magnetic field versus spin period (B-P) diagram, where their relations to the millisecond pulsars (MSPs) can be clearly shown. Generally, both BPSRs and MSPs are believed to be recycled and spun-up in binary accreting phases, and evolved below the spin-up line setting by the Eddington accretion rate ( $\dot{M}{\simeq}10^{18}~\mbox{g/s}$ ). It is noticed that most BPSRs are distributed around the spin-up line with mass accretion rate $\dot{M}=10^{16}~\mbox{g/s}$ and almost all MSP samples lie above the spin-up line with $\dot{M}\sim10^{15}~\mbox{g/s}$ . Thus, we calculate that a minimum accretion rate ( $\dot{M}\sim10^{15}~\mbox{g/s}$ ) is required for the MSP formation, and physical reasons for this are proposed. In the B-P diagram, the positions of BPSRs and their relations to the binary parameters, such as the companion mass, orbital period and eccentricity, are illustrated and discussed. In addition, for the seven BPSRs located above the limit spin-up line, possible causes are suggested.     

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.     

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$ .     

Virial mass in warped DGP-inspired \mathcal{L}(R) gravity

A version of the virial theorem is derived in a brane-world scenario in the framework of a warped DGP model where the action on the brane is an arbitrary function of the Ricci scalar, $\mathcal{L}(R)$ . The extra terms in the modified Einstein equations generate an equivalent mass term (geometrical mass), which give an effective contribution to the gravitational energy and offer viable explanation to account for the virial mass discrepancy in clusters of galaxies. We also obtain the radial velocity dispersion of galaxy clusters and show that it is compatible with the radial velocity dispersion profile of such clusters. Finally, we compare the result of the model with $\mathcal{L}(R)$ gravity theories.     

Alfvén Reflection and Reverberation in the Solar Atmosphere

Magneto-atmospheres with Alfvén speed [a] that increases monotonically with height are often used to model the solar atmosphere, at least out to several solar radii. A?common example involves a uniform vertical or inclined magnetic field in an isothermal atmosphere, for which the Alfvén speed is exponential. We address the issue of internal reflection in such atmospheres, both for time-harmonic and for transient waves. It is found that a mathematical boundary condition may be devised that corresponds to perfect absorption at infinity, and, using this, that many atmospheres where a(x) is analytic and unbounded present no internal reflection of harmonic Alfvén waves. However, except for certain special cases, such solutions are accompanied by a wake, which may be thought of as a kind of reflection. For the initial-value problem where a harmonic source is suddenly switched on (and optionally off), there is also an associated transient that normally decays with time as $\mathcal{O}(t^{-1})$ or $\mathcal{O}(t^{-1}\ln t)$ , depending on the phase of the driver. Unlike the steady-state harmonic solutions, the transient does reflect weakly. Alfvén waves in the solar corona driven by a finite-duration train of p-modes are expected to leave such transients.     

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.     

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.     

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}^{-}}$ .     

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$ .     

Estimating the size of a radio quiet zone for the radio astronomy service

The size of a radio quiet zone (RQZ) is largely determined by transmission losses of interfering signals, which can be divided into free space loss and diffraction loss. The free space loss is dominant. The diffraction loss presented in this paper is described as unified smooth spherical and knife edge diffractions, which is a function of minimum path clearance. We present a complete method to calculate the minimum path clearance. The cumulative distribution of the lapse rate of refractivity (g _n ), between the earth surface and 1 km above, is studied by using Chinese radio climate data. Because the size of an RQZ is proportional to g _n , the cumulative distribution of g _n can be used as an approximation for the size of the RQZ. When interference originates from mobile communication or television transmissions at a frequency of 408 MHz, and $\overline {g_n } $ is 40 N/km, where the refractivity $N=\left( {n-1} \right) \times 10^6$ , the size of the RQZ would be 180 km for a mobile source or 210 km for a television source, with a probability in the range of 15–100% in different months and for different stations. When speaking of the size of an RQZ, the radius in the case of a circular zone is implied. It results that a size of an RQZ is mainly influenced by transmission loss rather than effective radiated power. In the case where the distance between an interfering source and a radio astronomical observatory is about 100 km, at a frequency of 408 MHz, the allowable effective radiated power of the interfering source should be less than ?30 dBW with a probability of about 85% for $\overline {g_n } $ equals 40 N/km, or ?42 dBW with a probability less than 1 % for $\overline {g_n } $ equals 80 N/km.     

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.     

Oblique shock dynamics in nonextensive magnetized plasma

The detailed properties of the classical electromagnetic Weibel instability in a semi-relativistic anisotropic plasma are investigated for Maxwellian distribution. In this article, the effects of one particular factor affecting the growth rate of Weibel instability, Coulomb collision effect of electrons and ions, is studied and discussed based on the equilibrium semi-relativistic Maxwellian distribution function, in a dense and unmagnetized anisotropic plasma. An analytical expression is derived for the growth rate of the Weibel instability. The two limiting cases ( $\left| \xi\right| \ll1$ and $\left| \xi\right| \gg1$ ) are considered. It is shown that in the limit $\left| \xi\right| \ll 1$ , the quantity η, which is due to the collision term, will appear in the growth and in the conditions of the rate of the Weibel instability. The quantity χ symbolizes the contribution from relativistic terms which becomes unity as we approach the non-relativistic Maxwellian case, leading to the standard Weibel instability scenario. The increasing of η leads to decreasing of the growth rate, and with the decreasing of η the growth rate will increase.     

Revisiting Noether gauge symmetry for F(R) theory of gravity

Noether gauge symmetry for F(R) theory of gravity has been explored recently. The fallacy is that, even after setting gauge to vanish, the form of F(R)∝R ⁿ (where n≠1 is arbitrary) obtained in the process, has been claimed to be an outcome of gauge Noether symmetry. On the contrary, earlier works proved that any nonlinear form other than $F(R) \propto R^{\frac{3}{2}}$ is obscure. Here, we show that, setting gauge term zero, Noether equations are satisfied only for n=2, which again does not satisfy the field equations. Thus, as noticed earlier, the only form that Noether symmetry admits is $F(R) \propto R^{\frac{3}{2}}$ . Noether symmetry with non-zero gauge has also been studied explicitly here, to show that it does not produce anything new.     

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.