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.     

Ultra-high energy neutrino fluxes from supermassive AGN black holes

We compute the ultra-high energy (UHE) neutrino fluxes from plausible accreting supermassive black holes closely linking to the 377 active galactic nuclei (AGNs). They have well-determined black hole masses collected from the literature. The neutrinos are produced via simple or modified URCA processes, even after the neutrino trapping, in superdense proto-matter medium. The resulting fluxes are ranging from: (1) (quark reactions)— $J^{q}_{\nu\varepsilon}/(\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1})\simeq8.29\times 10^{-16}$ to 3.18×10^?4, with the average $\overline{J}^{q}_{\nu\varepsilon}\simeq5.53\times 10^{-10}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ , where ε _d ～10^?12 is the opening parameter; (2) (pionic reactions)— $J^{\pi}_{\nu\varepsilon} \simeq0.112J^{q}_{\nu\varepsilon}$ , with the average $J^{\pi}_{\nu\varepsilon} \simeq3.66\times 10^{-11}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ ; and (3) (modified URCA processes)— $J^{URCA}_{\nu\varepsilon}\simeq7.39\times10^{-11} J^{q}_{\nu\varepsilon}$ , with the average $\overline{J}^{URCA}_{\nu\varepsilon} \simeq2.41\times10^{-20} \varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ . We conclude that the AGNs are favored as promising pure neutrino sources, because the computed neutrino fluxes are highly beamed along the plane of accretion disk, peaked at high energies and collimated in smaller opening angle θ～ε _d .     

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.     

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.     

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.     

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.     

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.     

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.     

The motion around the libration points in the restricted three-body problem with the effect of radiation and oblateness

This paper deals with the existence of libration points and their linear stability when the more massive primary is radiating and the smaller is an oblate spheroid. Our study includes the effects of oblateness of $\bar{J}_{2i}$ (i=1,2) with respect to the smaller primary in the restricted three-body problem. Under combining the perturbed forces that were mentioned before, the collinear points remain unstable and the triangular points are stable for 0<μ<μ _c , and unstable in the range $\mu_{c} \le\mu\le\frac{1}{2}$ , where $\mu_{c} \in(0,\frac{1}{2})$ , it is also observed that for these points the range of stability will decrease. The relations for periodic orbits around five libration points with their semimajor, semiminor axes, eccentricities, the frequencies of orbits and periods are found, furthermore for the orbits around the triangular points the orientation and the coefficients of long and short periodic terms also are found in the range 0<μ<μ _c .     

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.     

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 cosmogonic shadow effect

We consider the Alfvén-Arrhenius fall-down mechanism and describe an approximate model for the infall, capture and distribution of dust particles on a given magnetic field line and their possible neutralization at the ‘2’/3 points, the points at which the field aligned compnents of the gravitational and centrifugal forces are equal and opposite. We find that a small fraction (<10%) of an incoming particle distribution will actually contribute to the above ‘2’/3 fall-down process. We also show that if at the 2/3 points, the ratio of dust to plasma density is $$\frac{{n_D \left( {\tfrac{2}{3}} \right)}}{{n_p \left( {\tfrac{2}{3}} \right)}} > \frac{{10^{ - 3} }}{{r_{g_\mu } T_{eV} }}$$ . (r _gμ=radius of a grain in microns,T=plasma temperature in eV), then the dust particles will lose their charge, decouple from the field line and follow Keplerian orbits in accordance with the Alfvén-Arrhenius mechanism. We then determine the limits on the plasma parameters in order that rotation of a quasi-neutral plasma in thermal equilibrium be possible in the gravitational and dipole field of a rotating central body. The constraints imposed by the above conditions are rather weak, and the plasma parameters can have a wide range of values. For a plasma corotating with an angular velocity Ω～10^?4s^?1, we show that the plasma temperature and density must satisfy $$10^{ - 1}<< T_{(eV)}<< 10^2 ,10T_{eV}^2<< n^p \left( {cm^3 } \right)<< 10^6 $$ .     

The Wave–Driver System of the Off-Disk Coronal Wave of 17 January 2010

We study the 17 January 2010 flare–CME–wave event by using STEREO/SECCHI-EUVI and -COR1 data. The observational study is combined with an analytic model that simulates the evolution of the coronal wave phenomenon associated with the event. From EUV observations, the wave signature appears to be dome shaped having a component propagating on the solar surface ( $\overline{v}\approx280~\mathrm{km}\,\mathrm{s}^{-1}$ ) as well as one off-disk ( $\overline{v}\approx 600~\mathrm{km}\,\mathrm{s}^{-1}$ ) away from the Sun. The off-disk dome of the wave consists of two enhancements in intensity, which conjointly develop and can be followed up to white-light coronagraph images. Applying an analytic model, we derive that these intensity variations belong to a wave–driver system with a weakly shocked wave, initially driven by expanding loops, which are indicative of the early evolution phase of the accompanying CME. We obtain the shock standoff distance between wave and driver from observations as well as from model results. The shock standoff distance close to the Sun (<?0.3 R _⊙ above the solar surface) is found to rapidly increase with values of ≈?0.03?–?0.09 R _⊙, which gives evidence of an initial lateral (over)expansion of the CME. The kinematical evolution of the on-disk wave could be modeled using input parameters that require a more impulsive driver (duration t=90 s, acceleration a=1.7 km?s^?2) compared to the off-disk component (duration t=340 s, acceleration a=1.5 km?s^?2).     

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.     

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.     

Investigation of X-ray cavities in the cooling flow system Abell 1991

We present results based on the systematic analysis of Chandra archive data on the X-ray bright Abell Richness class-I type cluster Abell 1991 with an objective to investigate properties of the X-ray cavities hosted by this system. The unsharp masked image as well as 2-d β model subtracted residual image of Abell 1991 reveals a pair of X-ray cavities and a region of excess emission in the central ～12 kpc region. Both the cavities are of ellipsoidal shape and exhibit an order of magnitude deficiency in the X-ray surface brightness compared to that in the undisturbed regions. Spectral analysis of X-ray photons extracted from the cavities lead to the temperature values equal to $1.77_{-0.12}^{+0.19}~\mathrm{keV}$ for N-cavity and $1.53_{-0.06}^{+0.05}~\mathrm{keV}$ for S-cavity, while that for the excess X-ray emission region is found to be equal to $2.06_{-0.07}^{+0.12}~\mathrm{keV}$ . Radial temperature profile derived for Abell 1991 reveals a positive temperature gradient, reaching to a maximum of 2.63 keV at ～76 kpc and then declines in outward direction. 0.5–2.0 keV soft band image of the central 15′′ region of Abell 1991 reveals relatively cooler three different knot like features that are about 10′′ off the X-ray peak of the cluster. Total power of the cavities is found to be equal to ${\sim}8.64\times10^{43}~\mathrm{erg\,s}^{-1}$ , while the X-ray luminosity within the cooling radius is found to be $6.04 \times10^{43}~\mathrm{erg\,s}^{-1}$ , comparison of which imply that the mechanical energy released by the central AGN outburst is sufficient to balance the radative loss.     

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.     

Effects of Field-Aligned Flows on Standing Kink and Sausage Modes Supported by Coronal Loops

Fundamental standing modes and their overtones play an important role in coronal seismology. We examine the effects of a significant field-aligned flow on standing modes that are supported by coronal loops, which are modeled here as cold magnetic slabs. Of particular interest are the period ratios of the fundamental to its (n?1)th overtone [P ₁/nP _n ] for kink and sausage modes, and the threshold half-width-to-length ratio for sausage modes. For standing kink modes, the flow significantly reduces P ₁/nP _n in general, the effect being particularly strong for higher n and weaker density contrast [ $\rho_{0}/\rho_{\rm e}$ ] between loops and their surroundings. That said, even when $\rho_{0}/\rho_{\rm e}$ approaches infinity, this effect is still substantial, reducing the minimal P ₁/nP _n by up to 13.7?% (24.5?%) for n=2 (n=4) relative to the static case, when the Alfvén Mach number [M _A] reaches 0.8, where M _A measures the loop flow speed in units of the internal Alfvén speed. Although it is not negligible for standing sausage modes, the flow effect in reducing P ₁/nP _n is not as strong. However, the threshold half-width-to-length ratio is considerably higher in the flowing case than in its static counterpart. For $\rho_{0}/\rho_{\rm e}$ in the range [9,1024] and M _A in the range [0,0.5], an exhaustive parameter study yields that this threshold is well fitted by $(d/L)_{\rm cutoff, fit} = \frac{1}{2}\sqrt{\frac{1}{\rho_{0}/\rho_{\rm e}-1}} \exp (3.7 M_{\mathrm{A}}^{2} )$ , which involves the two parameters in a simple way. This allows one to analytically constrain the combination $(\rho_{0}/\rho_{\rm e}, M_{\mathrm {A}})$ for a loop with a known width-to-length ratio when a standing sausage oscillation is identified. It also allows one to examine the idea of partial sausage modes in more detail, and the flow is found to significantly reduce the spatial extent where partial modes are allowed.