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.     

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.     

Perspectives on effectively constraining the location of a massive trans-Plutonian object with the New Horizons spacecraft: a sensitivity analysis

The radio tracking apparatus of the New Horizons spacecraft, currently traveling to the Pluto system where its arrival is scheduled for July 2015, should be able to reach an accuracy of 10 m (range) and 0.1  $\text{ mm } \text{ s }^{-1}$ mm s ? 1 (range-rate) over distances up to 50 au. This should allow to effectively constrain the location of a putative trans-Plutonian massive object, dubbed Planet X (PX) hereafter, whose existence has recently been postulated for a variety of reasons connected with, e.g., the architecture of the Kuiper belt and the cometary flux from the Oort cloud. Traditional scenarios involve a rock-ice planetoid with $m_\mathrm{X}\approx 0.7\,m_{\oplus }$ m X ≈ 0.7 m ⊕ at some 100–200 au, or a Jovian body with $m_\mathrm{X}\lesssim 5\,m_\mathrm{J}$ m X ? 5 m J at about 10,000–20,000 au; as a result of our preliminary sensitivity analysis, they should be detectable by New Horizons since they would impact its range at a km level or so over a time span 6 years long. Conversely, range residuals statistically compatible with zero having an amplitude of 10 m would imply that PX, if it exists, could not be located at less than about 4,500 au ( $m_\mathrm{X}=0.7\,m_{\oplus }$ m X = 0.7 m ⊕ ) or 60,000 au ( $m_\mathrm{X}=5\,m_\mathrm{J}$ m X = 5 m J ), thus making a direct detection quite demanding with the present-day technologies. As a consequence, it would be appropriate to rename such a remote body as Thelisto. Also fundamental physics would benefit from this analysis since certain subtle effects predicted by MOND for the deep Newtonian regions of our Solar System are just equivalent to those of a distant pointlike mass.     

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.     

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.     

Electron capture of nuclei in the region surrounding magnetar

Effects of ultra-strong magnetic field on electron capture rates for ⁵⁷Fe, ⁵⁸Co and ⁵⁹Ni have been analyzed in the nuclear shell model and under the Landau energy levels quantized approximation in the ultra-strong magnetic field, the result increase about 3 orders magnitude. The rate of change of electron abundance, $\dot{Y}_{e}$ , for every nuclide and total $\dot{Y}_{e}$ in the condition without magnetic field and B=4.414×10¹⁵ G have been calculated, and exceed about 6 orders of magnitude generally. These conclusions play an important role in future studying the evolution of magnetar.     

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.     

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

Optical spectral decomposition of NGC 4051

Considering the host galaxy contribution, a spectral decomposition method is used to reanalyzed the archive data of optical spectra for a narrow line Seyfert 1 galaxy, NGC 4051. The light curves of the continuum f _λ (5100 Å), and Hβ, He ii, Fe ii emission lines are given. We find strong flux correlations between line emissions of Hβ, He ii, Fe ii and the continuum f _λ (5100 Å). These low-ionization lines (Hβ, Fe ii, He ii) have “inverse” intrinsic Baldwin effects. Using the methods of the cross-correlation function and the Monte Carlo simulation, we find the time delays, with respect to the continuum, are $3.45^{+12.0}_{-0.5}~\mbox{days}$ with the probability of 34 % for the intermediate component of Hβ, $6.45^{+13.0}_{-1.0}~\mbox{days}$ with the probability of 65 % for the intermediate component of He ii. From these intermediate components of Hβ and He ii, the calculated central black hole masses are $0.86^{+4.35}_{-0.33}\times 10^{6}$ and $0.82^{+3.12}_{-0.45}\times 10^{6}~M_{\odot }$ . We also find that the time delays for Fe ii are $9.7^{+3.0}_{-5.0}~\mbox{days}$ with the probability of 36 %, $8.45^{+1.0}_{-2.0}~\mbox{days}$ with the probability of 18 % for the total epochs and “subset 1” data, respectively. It seems that the Fe ii emission region is outside of the Hβ emission region.     

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.     

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.     

On the rotational energy distributions of reactive, non-polar species in the interstellar medium

A basic model for the formation of non-equilibrium rotational energy distributions is described for reactive, homo-polar diatomic molecules and ions in the interstellar medium. Kinetic models were constructed to calculate the rotational populations of $\mathrm{C}_{2}^{+}$ under the conditions it would experience in the diffuse interstellar medium. As the non-polar ion reacts with molecular hydrogen, but not atomic hydrogen, the thermalization of a hot nascent rotational population will be arrested by chemical reaction when the H₂ density begins to be significant. Populations that deviate strongly from the local thermodynamic equilibrium are predicted for $\mathrm{C}_{2}^{+}$ in environments where it may be detectable. Consequences of this are discussed and a new optical spectrum is calculated.     

Photoionization Cross-Section of Chlorine-like Iron

Semi-relativistic calculations are performed for the photoionization of Fe X (an important coronal ion) from its ground state 3s²3p $^{5}(^{2}P^{0}_{3/2})$ and the first two excited states 3s²3p $^{5}(^{2}P^{0}_{1/2})$ and 3s3p $^{6}(^{2}S_{1/2})$ using the Breit?CPauli R-matrix method. A lowest 41 state eigenfunction expansion for Fe XI is employed to ensure an extensive treatment of auto ionizing resonances that affect the effective cross-sections. In the present calculations, we have considered all the important physical effects like channel coupling, exchange and short range correlation. The present calculations using the lowest 41 target levels of Fe XI in the LSJ coupling scheme are reported and we expect that the present results should enable more accurate modelling of the emission spectrum of highly excited plasma from the optical to the far UV region.     

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

A new look at the origin of the 6.67 hr period X-ray pulsar 1E 161348-5055

The point X-ray source 1E 161348-5055 is observed to display pulsations with the period 6.67?hr and $|\dot{P}| \leq1.6 \times10^{-9}\,{\rm s\,s^{-1}}$ . It is associated with the supernova remnant RCW?103 and is widely believed to be a ～2000?yr old neutron star. Observations give no evidence for the star to be a member of a binary system. Nevertheless, it resembles an accretion-powered pulsar with the magnetospheric radius ～3000?km and the mass-accretion rate $\sim 10^{14}\,{\rm g\,s^{-1}}$ . This situation could be described in terms of accretion from a (residual) fossil disk established from the material falling back towards the star after its birth. However, current fall-back accretion scenarios encounter major difficulties explaining an extremely long spin period of the young neutron star. We show that the problems can be avoided if the accreting material is magnetized. The star in this case is surrounded by a fossil magnetic slab in which the material is confined by the magnetic field of the accretion flow itself. We find that the surface magnetic field of the neutron star within this scenario is ～10¹²?G and that a presence of $\gtrsim10^{-7}\,{\rm M_{\odot}}$ magnetic slab would be sufficient to explain the origin and current state of the pulsar.     

Tidal synchronization of close-in satellites and exoplanets. A rheophysical approach

This paper presents a new theory of the dynamical tides of celestial bodies. It is founded on a Newtonian creep instead of the classical delaying approach of the standard viscoelastic theories and the results of the theory derive mainly from the solution of a non-homogeneous ordinary differential equation. Lags appear in the solution but as quantities determined from the solution of the equation and are not arbitrary external quantities plugged in an elastic model. The resulting lags of the tide components are increasing functions of their frequencies (as in Darwin’s theory), but not small quantities. The amplitudes of the tide components depend on the viscosity of the body and on their frequencies; they are not constants. The resulting stationary rotations (often called pseudo-synchronous) have an excess velocity roughly proportional to $6ne^2/(\chi ^2+\chi ^{-2})$ ( $\chi $ is the mean-motion in units of one critical frequency—the relaxation factor—inversely proportional to the viscosity) instead of the exact $6ne^2$ of standard theories. The dissipation in the pseudo-synchronous solution is inversely proportional to $(\chi +\chi ^{-1})$ ; thus, in the inviscid limit, it is roughly proportional to the frequency (as in standard theories), but that behavior is inverted when the viscosity is high and the tide frequency larger than the critical frequency. For free rotating bodies, the dissipation is given by the same law, but now $\chi $ is the frequency of the semi-diurnal tide in units of the critical frequency. This approach fails, however, to reproduce the actual tidal lags on Earth. In this case, to reconcile theory and observations, we need to assume the existence of an elastic tide superposed to the creeping tide. The theory is applied to several Solar System and extrasolar bodies and currently available data are used to estimate the relaxation factor $\gamma $ (i.e. the critical frequency) of these bodies.     

A method for judging decay or growth of the magnetic field of pulsar

This paper provides a method for judging growth or decay of the magnetic field of pulsar by using pulse period P, or frequency ν, and its first and second derivatives $ \dot P,\ddot P $ or $ \dot v,\ddot v $ . The author uses this method to judge the growth or decay of the magnetic field of Crab pulsar. The judged result for Crab pulsar is that the magnetic field of Crab pulsar is growing now, but it is not decaying. The result corresponds with the actual case of Crab pulsar.     

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.     

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.