Omnidirectional induced compton scattering by relativistic electrons

Quasars, pulsars and other cosmic sources of intense radiation are known to have large brightness temperature (kT _b?mc ²) and relativistic electron density values. In this case the induced Compton scattering by relativistic electrons should be considered. The probability of scattering with decreasing radiation frequency is derived for isotropic radiation scattering. When induced scattering takes place, the relativistic electron obtains its energy by transforming high-frequency quanta into the low-frequency ones. In the most intensive sources electrons would receive energiesE?mc ² ××(kT _b/mc ²)^1/7 due to the heating rate proportional toE ^?5 with the cooling rate proportional toE ². Considerable distortion of the quasar spectrum is possible for reasonably large values of relativistic electron density (N?10⁶cm^?3) notwithstanding that the heating is negligible. In pulsars relativistic electron heating and spectrum distortion appear to depend more on the induced Compton scattering.     

Resonant Compton scattering associated with pair creation

Studies of Compton scattering by relativistic electrons in a strong magnetic field have been restricted to either incident photon angles θ′ aligned along the magnetic field B or incident photon energies ω′ below the first pair creation threshold $\omega'_{PC}$ . When these restrictions are relaxed there is a resonance in Compton scattering associated with pair creation (PC), that is analogous to but independent of known resonances associated with gyromagnetic absorption (GA). As with the GA resonances, that may be labeled by the Landau quantum numbers of the relevant states, there is a sequence of PC resonances where the scattering cross section diverges. In this paper, the lowest divergence is studied for incident photon energies satisfying ω′²sin² θ′/(2eB)?1, assuming that the scattering electron is in its ground (Landau) state. This lowest resonance affects only parallel-polarized photons.     

On the radiative and thermodynamic properties of the cosmic radiations using COBE FIRAS instrument data: IV. Sunyaev-Zel’dovich ($\mu$-type) distortion effect

The Sunyaev-Zel’dovich (SZ) effect represents a small spectral distortion to the cosmic microwave background (CMB) radiation, caused by the Compton scattering of CMB photons by the hot gas of galaxy clusters. In an early stage of universe, the SZ effect generates \(\mu\)-type of distortions for the CMB spectrum. A \(\mu\)-type distortion is created between the double Compton scattering decoupling (\(z \sim 10^{6}\)) and the thermalization decoupling by the Compton scattering (\(z \sim 10^{5}\)). In this case, to describe the small spectral distortion of the CMB spectrum, we use the Bose-Einstein (\(\mu\)-type) distribution with a non-zero chemical potential. At present, it is interesting to investigate the effect of this spectral distortion on the integral characteristics of the Bose-Einstein (\(\mu\)-type) spectrum. The thermal radiative and thermodynamic functions are such integral characteristics. These functions are as follows: a) the total radiation power per unit area; b) total energy density; c) number density of photons; d) grand potential density; e) Helmholtz free energy density; f) entropy density; g) heat capacity at constant volume; h) enthalpy density; and i) pressure. Precise analytical expressions are obtained for the temperature dependences of these functions. Using the observational data obtained by the COBE FIRAS, PIXIE, PRISM, and Planck missions, the thermal radiative and thermodynamic functions are calculated. A comparative analysis of the results obtained with the results for the same functions of the CMB spectrum at \(T = 2.72548~\mbox{K}\) is carried out. Very small distortions are observed for the thermal radiative and thermodynamic functions. In the redshift range \(10^{5} < z < 3 \times10^{6}\), these functions are calculated. The expressions are obtained for new astrophysical parameters, such as the entropy density/Boltzmann constant and number density, created by the Bose-Einstein (\(\mu\)-type) spectrum.     

Is the GeV-TeV emission of PKS 0447-439 from the proton synchrotron radiation?

We study the multi-wavelength emission features of PKS 0447-439 in the frame of the one-zone homogeneous lepto-hadronic model. In this model, we assumed that the steady power-laws with exponential cut-offs distributions of protons and electrons are injected into the source. The non-linear time-dependent kinematic equations, describing the evolution of protons, electrons and photons, are defined; these equations self-consistently involve synchrotron radiation of protons, photon-photon interaction, synchrotron radiation of electron/positron pairs, inverse Compton scattering and synchrotron self-absorption. The model is applied to reproduce the multi-wavelength spectrum of PKS 0447-439. Our results indicate that the spectral energy distribution (SED) of PKS 0447-439 can be reproduced well by the model. In particular, the GeV-TeV emission is produced by the synchrotron radiation of relativistic protons. The physically plausible solutions require the magnetic strength \(10~\text{G}\lesssim B \lesssim 100~\text{G}\). We found that the observed spectrum of PKS 0447-439 can be reproduced well by the model whether \(z = 0.16\) or \(z = 0.2\), and the acceptable upper limit of redshift is \(z=0.343\).     

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.     

Broadband gravitational-wave pulses from binary neutron stars in eccentric orbits

The gravitational-wave radiation from binary stars in elliptical orbits peaks at times close to the periastron passage. For a stationary distribution of binary neutron stars in the Galaxy, there are several systems with large orbital eccentricities and periods in the range from several tens of minutes to several days from which gravitational-wave radiation at periastron will be observed as a broad pulse in the frequency range 1–100 mHz. The LISA space interferometer will be able to record pulsed signals from these systems at a signal-to-noise ratio $S/N > 5\sqrt 5$ in the frequency range ～10^?3–10^?1 Hz. Algorithms for detecting such signals are discussed.     

Discrete frequency scattering in spherical layers surrounding a point source

The scattered radiation field in homogeneously absorbing and isotropically scattering spherical layers is studied, when the isotropic point source is at the centre. A complete frequency redistribution is assumed. It is shown, that on the inner boundaryr=R ₀ of the cavity, whenR ₀?1 (all radii are expressed in the path lengths), the source functionB～R ₀ ^?1 ln ^?1/2 R ₀ for the Doppler profile andB～R ₀ ^?3/2 for the Voigt and Lorentz profiles. The asymptotical behaviour of the source functionB(r) significantly differs from the analogous behaviour of solution for an infinite medium.     

Small-scale fluctuations of relic radiation

Perturbations of the matter density in a homogeneous and isotropic cosmological model which leads to the formation of galaxies should, at later stages of evolution, cause spatial fluctuations of relic radiation. Silk assumed that an adiabatic connection existed between the density perturbations at the moment of recombination of the initial plasma and fluctuations of the observed temperature of radiation δT/T=δ _?m /3 _?m . It is shown in this article that such a simple connection is not applicable due to:

1. The long time of recombination;
2. The fact that when regions withM<10¹⁵ M _⊙ become transparent for radiation, the optical depth to the observer is still large due to Thompson scattering;
3. The spasmodic increase of δ _?m/?m in recombination.

As a result the expected temperature fluctuations of relic radiation should be smaller than adiabatic fluctuations. In this article the value of δT/T arising from scattering of radiation on moving electrons is calculated; the velocity field is generated by adiabatic or entropy density perturbations. Fluctuations of the relic radiation due to secondary heating of the intergalactic gas are also estimated. A detailed investigation of the spectrum of fluctuations may, in principle, lead to an understanding of the nature of initial density perturbations since a distinct periodic dependence of the spectral density of perturbations on wavelength (mass) is peculiar to adiabatic perturbations. Practical observations are quite difficult due to the smallness of the effects and the presence of fluctuations connected with discrete sources of radio emission.     

Stability of triangular equilibrium points in the elliptic restricted problem of three bodies with radiating and triaxial primaries

This paper studies the stability of infinitesimal motions about the triangular equilibrium points in the elliptic restricted three body problem assuming bigger primary as a source of radiation and the smaller one a triaxial rigid body. The perturbation technique developed by Bennet (Icarus 4:177, 1965b) has been used for determination of characteristic exponents. This technique is based on Floquet’s Theory for determination of characteristic exponents in the system with periodic coefficients. The results of the study are analytical and numerical expressions are simulated for the transition curves bounding the region of stability in the μ–e plane, accurate to O(e ²). The unstable region is found to be divided into three parts. The effect of radiation parameter is significant. For small values of e, the results are in favor with the numerical analysis of Danby (Astron. J. 69:166, 1964), Bennet (Icarus 4:177, 1965b), Alfriend and Rand (AIAA J. 6:1024, 1969). The effect of radiation pressure is significant than the oblateness and triaxiality of the primaries.     

Symmetric and regularized coordinates on the plane triple collision manifold

The planar problem of three bodies is described by means of Murnaghan's symmetric variables (the sidesa _j of the triangle and an ignorable angle), which directly allow for the elimination of the nodes. Then Lemaitre's regularized variables \(\alpha _j = \sqrt {(\alpha ^2 - \alpha _j )}\) , where \(\alpha ^2 = \tfrac{1}{2}(a_1 + a_2 + a_3 )\) , as well as their canonically conjugated momenta are introduced. By finally applying McGehee's scaling transformation \(\alpha _j = r^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}} \tilde \alpha _j\) , wherer ² is the moment of inertia a system of 7 differential equations (with 2 first integrals) for the 5-dimensional triple collision manifold \(T\) is obtained. Moreover, the zero angular momentum solutions form a 4-dimensional invariant submanifold \(N \subset T\) represented by 6 differential equations with polynomial right-hand sides. The manifold \(N\) is of the topological typeS ²×S ² with 12 points removed, and it contains all 5 restpoint (each one in 8 copies). The flow on \(T\) is gradient-like with a Lyapounov function stationary in the 40 restpoints. These variables are well suited for numerical studies of planar triple collision.     

Spectrum of the cosmic X- and gamma ray background in the energy range 1 keV-1 MeV

Available satellite, rocket and balloon observations on cosmic X- and gamma ray background are critically examined to understand the spectral characteristics of the radiation. Appropriate corrections have been applied to the balloon observations to account for the multiple Compton scattering of X-rays in the atmosphere. It is shown that within the experimental uncertainties, all the available observations of cosmic X- and gamma ray background in the energy range 1 keV-1 MeV are consistent with a single spectrum of type $${\text{d}}N/{\text{d}}E = 30 E^{ - 2.0 \pm 0.2} {\text{photons cm}}^{{\text{ - 2}}} {\text{s}}^{{\text{ - 1}}} {\text{sr}}^{{\text{ - 1}}} {\text{keV}}^{{\text{ - 1}}} $$ .     

Non-thermal processes in large solar flares

We analyze particle acceleration processes in large solar flares, using observations of the August, 1972, series of large events. The energetic particle populations are estimated from the hard X-ray and γ-ray emission, and from direct interplanetary particle observations. The collisional energy losses of these particles are computed as a function of height, assuming that the particles are accelerated high in the solar atmosphere and then precipitate down into denser layers. We compare the computed energy input with the flare energy output in radiation, heating, and mass ejection, and find for large proton event flares that:

1. The ～10–10² keV electrons accelerated during the flash phase constitute the bulk of the total flare energy.
2. The flare can be divided into two regions depending on whether the electron energy input goes into radiation or explosive heating. The computed energy input to the radiative quasi-equilibrium region agrees with the observed flare energy output in optical, UV, and EUV radiation.
3. The electron energy input to the explosive heating region can produce evaporation of the upper chromosphere needed to form the soft X-ray flare plasma.
4. Very intense energetic electron fluxes can provide the energy and mass for interplanetary shock wave by heating the atmospheric gas to energies sufficient to escape the solar gravitational and magnetic fields. The threshold for shock formation appears to be ～10³¹ ergs total energy in >20 keV electrons, and all of the shock energy can be supplied by electrons if their spectrum extends down to 5–10 keV.
5. High energy protons are accelerated later than the 10–10² keV electrons and most of them escape to the interplanetary medium. The energetic protons are not a significant contributor to the energization of flare phenomena. The observations are consistent with shock-wave acceleration of the protons and other nuclei, and also of electrons to relativistic energies.
6. The flare white-light continuum emission is consistent with a model of free-bound transitions in a plasma with strong non-thermal ionization produced in the lower solar chromosphere by energetic electrons. The white-light continuum is inconsistent with models of photospheric heating by the energetic particles. A threshold energy of ～5×10³⁰ ergs in >20 keV electrons is required for detectable white-light emission.

The highly efficient electron energization required in these flares suggests that the flare mechanism consists of rapid dissipation of chromospheric and coronal field-aligned or sheet currents, due to the onset of current-driven Buneman anomalous resistivity. Large proton flares then result when the energy input from accelerated electrons is sufficient to form a shock wave.     

Cosmic-ray antimatter: A primary origin hypothesis

We examine the possibility that the observed cosmic-ray protons are of primary extragalactic origin. The present \(\bar p\) data are consistent with a primary extragalactic component having \(\bar p\) /p?3.2±0.7 x 10^-4 independent of energy. Following the suggestion that most extragalactic cosmic rays are from active galaxies, we propose that most of the observed \(\bar p\) 's are alos from the same sites. This would imply the possibility of destroying the corresponding \(\bar \alpha \) 'sat the source, thus leading to a flux ratio \(\bar \alpha \) /α< \(\bar p\) /p. We further predict an estimate for \(\bar \alpha \) α～10^-5, within the range of future cosmic-ray detectors. the cosmological implications of this proposal are discussed.     

15-digit accuracy calculations of Ambartsumian-Chandrasekhar’s $$-functions for four-term phase functions with the double-exponential formula

We have established an iterative scheme to calculate with 15-digit accuracy the numerical values of Ambartsumian-Chandrasekhar’s \(H\)-functions for anisotropic scattering characterized by the four-term phase function: the method incorporates some advantageous features of the iterative procedure of Kawabata (Astrophys. Space Sci. 358:32, 2015) and the double-exponential integration formula (DE-formula) of Takahashi and Mori (Publ. Res. Inst. Math. Sci. Kyoto Univ. 9:721, 1974), which proved highly effective in Kawabata (Astrophys. Space Sci. 361:373, 2016). Actual calculations of the \(H\)-functions have been carried out employing 27 selected cases of the phase function, 56 values of the single scattering albedo \(\varpi_{0}\), and 36 values of an angular variable \(\mu(= \cos\theta)\), with \(\theta\) being the zenith angle specifying the direction of incidence and/or emergence of radiation. Partial results obtained for conservative isotropic scattering, Rayleigh scattering, and anisotropic scattering due to a full four-term phase function are presented. They indicate that it is important to simultaneously verify accuracy of the numerical values of the \(H\)-functions for \(\mu<0.05\), the domain often neglected in tabulation. As a sample application of the isotropic scattering \(H\)-function, an attempt is made in Appendix to simulate by iteratively solving the Ambartsumian equation the values of the plane and spherical albedos of a semi-infinite, homogeneous atmosphere calculated by Rogovtsov and Borovik (J. Quant. Spectrosc. Radiat. Transf. 183:128, 2016), who employed their analytical representations for these quantities and the single-term and two-term Henyey-Greenstein phase functions of appreciably high degrees of anisotropy. While our results are in satisfactory agreement with theirs, our procedure is in need of a faster algorithm to routinely deal with problems involving highly anisotropic phase functions giving rise to near-conservative scattering.     

Rational approximation formula for Chandrasekhar’s <Emphasis Type="Italic">H</Emphasis>-function for isotropic scattering

In this work, we first establish a simple procedure to obtain with 11-figure accuracy the values of Chandrasekhar’s H-function for isotropic scattering using a closed-form integral representation and the Gauss-Legendre quadrature. Based on the numerical values of the function produced by this method for various combinations of ? ₀, the single scattering albedo, and μ, the cosine of the zenith angle θ of the direction of radiation emergent from or incident upon a semi-infinite scattering-absorbing medium, we propose a rational approximation formula with μ ^1/4 and \(\sqrt{1-\varpi_{0}}\) as the independent variables. This allows us to reproduce the correct values of H(? ₀,μ) within a relative error of 2.1×10^?5 without recourse to any iterative procedure or root-finding process.     

Microwave emission above steady and moving sunspots

Two-dimensional maps of radio brightness temperature and polarization, computed assuming thermal emission with free-free and gyroresonance absorption, are compared with observations of active region 2502, performed at Westerbork at λ = 6.16 cm during a period of 3 days in June 1980. The computation is done assuming a homogeneous model in the whole field of view (5′ × 5′) and a force-free extrapolation of the photospheric magnetic field observed at MSFC with a resolution of 2″.34. The mean results are the following:

1. A very good agreement is found above the large leading sunspot of the group, assuming a potential extrapolation of the magnetic field and a constant conductive flux in the transition region ranging from 2 × 10⁶ to 10⁷ erg cm^?2s^?1.
2. A strong radio source, associated with a new-born moving sunspot, cannot be ascribed to thermal emission. It is suggested that this source may be due to synchrotron radiation by mildly relativistic electrons accelerated by resistive instabilities occurring in the evolving magnetic configuration. An order-of-magnitude computation of the expected number of accelerated particles seems to confirm this hypothesis.

     

Evolution of the \mathcal {L}_1 halo family in the radial solar sail circular restricted three-body problem

We present a detailed investigation of the dramatic changes that occur in the \(\mathcal {L}_1\) halo family when radiation pressure is introduced into the Sun–Earth circular restricted three-body problem (CRTBP). This photo-gravitational CRTBP can be used to model the motion of a solar sail orientated perpendicular to the Sun-line. The problem is then parameterized by the sail lightness number, the ratio of solar radiation pressure acceleration to solar gravitational acceleration. Using boundary-value problem numerical continuation methods and the AUTO software package (Doedel et al. in Int J Bifurc Chaos 1:493–520, 1991) the families can be fully mapped out as the parameter \(\beta \) is increased. Interestingly, the emergence of a branch point in the retrograde satellite family around the Earth at \(\beta \approx 0.0387\) acts to split the halo family into two new families. As radiation pressure is further increased one of these new families subsequently merges with another non-planar family at \(\beta \approx 0.289\) , resulting in a third new family. The linear stability of the families changes rapidly at low values of \(\beta \) , with several small regions of neutral stability appearing and disappearing. By using existing methods within AUTO to continue branch points and period-doubling bifurcations, and deriving a new boundary-value problem formulation to continue the folds and Krein collisions, we track bifurcations and changes in the linear stability of the families in the parameter \(\beta \) and provide a comprehensive overview of the halo family in the presence of radiation pressure. The results demonstrate that even at small values of \(\beta \) there is significant difference to the classical CRTBP, providing opportunity for novel solar sail trajectories. Further, we also find that the branch points between families in the solar sail CRTBP provide a simple means of generating certain families in the classical case.     

On the velocity of jets powering hot spots similar to those in Cygnus A

Hot spots similar to those in the radio galaxy Cygnus A can be explained by the strong shock produced by a supersonic but classical jet \(\left( {u_{jet}< c/\sqrt 3 } \right)\) . The high integrated radio luminosity (L?2×10⁴⁴ erg s^?1) and the strength of mean magnetic field (B?2×10^?4 G) suggest the hot spots are the downstream flow of a very strong shock which generates the ultrarelativistic electrons of energy ?≥20 MeV. The fully-developed subsonic turbulence amplifies the magnetic field of the jet up to 1.6×10^?4 G by the dynamo effect. If we assume that the post-shock pressure is dominated by relativistic particles, the ratio between the magnetic energy density to the energy density in relativistic particles is found to be ?2×10^?2, showing that the generally accepted hypothesis of equipartition is not valid for hot spots. The current analysis allows the determination of physical parameters inside hot spots. It is found that:

1. The velocity of the upstream flow in the frame of reference of the shock isu ₁?0.2c. Radio observations indicate that the velocity of separation of hot spots isu _sep?0.05c, so that the velocity of the jet isu _jet=u ₁+u _sep?0.25c.
2. The density of the thermal electrons inside the hot spot isn ₂?5×10^?3 e ^? cm^?3 and the mass ejected per year to power the hot spot is ?4M ₀yr^?1.
3. The relativistic electron density is less than 20% of the thermal electron density inside the hot spot and the spectrum is a power law which continues to energies as low as 30 MeV.
4. The energy density of relativistic protons is lower than the energy density of relativistic electrons unlike the situation for cosmic rays in the Galaxy.

     

Fourier analysis of the light curves of eclipsing variables,XXIII

The method of evaluating the photometric perturbationsB _2m of eclipsing variables in the frequency domain, developed by Kopal (1959, 1975e, 1978) for an interpretation of mutual eclipses in systems whose components are distorted by axial rotation and mutual tidal action. The aim of the present paper has been to establish explicit expressions for the photometric perturbationB _2m in such systems, regardless of the kind of eclipses and non-integral values ofm. Recently, Kopal (1978) introduced two different kinds of integrals with respect to associated α-functions andI-integrals which have been expressed in terms of certain general types of series that can be easily programmed for automatic computation within seconds of real time on highspeed computers. Following a brief introduction (Section 1) in which the need of this new approach will be expounded, in Section 3 we shall deduce the integral $$\int_0^{\theta \prime } {\tfrac{{\alpha _n^\prime }}{\delta }} d(sin^{2m} \theta )$$ in terms of a certain general type of series and also β-function, which should enable us to evaluate explicit expressions forf _* ^(h) ,f ₁ ^(h) ,f ₂ ^(h) as well asB _2m .     

A new approximate solution of parallel whistler-mode dispersion equation

A new approximation of the real part of the nonrelativistic plasma dispersion function ?Z of a real argument ξ₀ is proposed: namely, $$\Re Z_A (\xi _0 ) = \frac{{2\xi _0 }}{{1 - 2\xi _0^2 }}.$$ This approximation gives the exact value for ?Z when d²?Z/dξ ₀ ² = 0 and gives the correct expressions for the first two terms of its expansion for large ξ₀. On the basis of this approximation, a new approximate expression for whistler-mode refractive index is derived for the case of wave propagation parallel to the magnetic field in a hot anisotropic and dense plasma. Under certain conditions this expression is more convenient for practical applications in magnetospheric and/or astrophysical conditions than other approximate expressions. The approximation ?Z _A was also used in deriving the approximate expression for whistler-mode instability or damping (γ), although in this case it seems to have little merit when compared with the previously used expressions for γ.