Harmonic analysis of cosmic microwave background data – II. From ring-sets to the sky

Despite the fact that the physics of the cosmic microwave background anisotropies is most naturally expressed in Fourier space, pixelized maps are almost always used in the analysis and simulation of microwave data. A complementary approach is investigated here, in which maps are used only in the visualization of the data, and the temperature anisotropies and polarization are only ever expressed in terms of their spherical multipoles. This approach has a number of advantages: there is no information loss (assuming a band-limited observation); deconvolution of asymmetric beam profiles and the temporal response of the instrument are naturally included; correlated noise can easily be taken into account, removing the need for additional 'destriping'; polarization is also analysed in the same framework; and reliable estimates of the spherical multipoles of the sky and their errors are obtained directly for subsequent component separation and power spectrum estimation. The formalism required to analyse experiments which survey the full sky by scanning on circles is derived here, with particular emphasis on the Planck mission. A number of analytical results are obtained in the limit of simple scanning strategies. Although there are non-trivial computational obstacles to be overcome before the techniques described here can be implemented at high resolution, if these can be overcome the method should allow for a more robust return from the next generation of full-sky microwave background experiments.     

Alignment and signed-intensity anomalies in Wilkinson Microwave Anisotropy Probe data

Significant alignment and signed-intensity anomalies of local features of the cosmic microwave background (CMB) are detected on the three-year Wilkinson Microwave Anisotropy Probe data, through a decomposition of the signal with steerable wavelets on the sphere. In addition to identifying local features of a signal at specific scales, steerable wavelets allow one to determine their local orientation and signed intensity. First, an alignment analysis identifies two mean preferred planes in the sky, both with normal axes close to the CMB dipole axis. The first plane is defined by the directions towards which local CMB features are anomalously aligned. A mean preferred axis is also identified in this plane, located very close to the ecliptic poles axis. The second plane is defined by the directions anomalously avoided by local CMB features. This alignment anomaly provides further insight on recent results. Secondly, a signed-intensity analysis identifies three mean preferred directions in the southern Galactic hemisphere with anomalously high or low temperature of local CMB features: a cold spot essentially identified with a known cold spot, a second cold spot lying very close to the southern end of the CMB dipole axis, and a hotspot lying close to the southern end of the ecliptic poles axis. In both analyses, the anomalies are observed at wavelet scales corresponding to angular sizes around 10° on the celestial sphere, with global significance levels around 1 per cent. Further investigation reveals that the alignment and signed-intensity anomalies are only very partially related. Instrumental noise, foreground emissions and some form of other systematics are strongly rejected as possible origins of the detections. An explanation might still be envisaged in terms of a global violation of the isotropy of the Universe, inducing an intrinsic statistical anisotropy of the CMB.     

Compressed sensing imaging techniques for radio interferometry

Radio interferometry probes astrophysical signals through incomplete and noisy Fourier measurements. The theory of compressed sensing demonstrates that such measurements may actually suffice for accurate reconstruction of sparse or compressible signals. We propose new generic imaging techniques based on convex optimization for global minimization problems defined in this context. The versatility of the framework notably allows introduction of specific prior information on the signals, which offers the possibility of significant improvements of reconstruction relative to the standard local matching pursuit algorithm CLEAN used in radio astronomy. We illustrate the potential of the approach by studying reconstruction performances on simulations of two different kinds of signals observed with very generic interferometric configurations. The first kind is an intensity field of compact astrophysical objects. The second kind is the imprint of cosmic strings in the temperature field of the cosmic microwave background radiation, of particular interest for cosmology.     

Extrapolation of the solar magnetic field within the potential-field approximation from full-disk magnetograms

This paper is concerned with the Laplace boundary-value problem with the directional derivative, corresponding to the specific nature of measurements of the longitudinal component of the photospheric magnetic field. The boundary conditions are specified by a distribution on the sphere of the projection of the magnetic field vector into a given direction, i.e., they exactly correspond to the data of daily magnetograms distributed across the full solar disk. It is shown that the solution of this problem exists in the form of a spherical harmonic expansion, and uniqueness of this solution is proved. A conceptual sketch of numerical determination of the harmonic series coefficients is given. The field of application of the method is analyzed with regard to the peculiarities of actual data. Results derived from calculating magnetic fields from real magnetograms are presented. Finally, we present differences in results derived from extrapolating the magnetic field from a synoptic map and a full-disk magnetogram.     

From planes to spheres: about gravitational lens magnifications

We discuss the classic theorem according to which a gravitational lens always produces at least one image with a magnification greater than unity. This theorem seems to contradict the conservation of total flux from a lensed source. The standard solution to this paradox is based on the exact definition of the reference 'unlensed' situation, in which the lens mass can be either removed or smoothly redistributed.
We calculate magnifications and amplifications (in photon number and energy flux density) for general lensing scenarios not limited to regions close to the optical axis. In this way the formalism is naturally extended from tangential planes for the source and lensed images to complete spheres. We derive the lensing potential theory on the sphere and find that the Poisson equation is modified by an additional source term that is related to the mean density and to the Newtonian potential at the positions of observer and source. This new term generally reduces the magnification, to below unity far from the optical axis, and ensures conservation of the total photon number received on a sphere around the source.
This discussion does not affect the validity of the focusing theorem , in which the unlensed situation is defined to have an unchanged affine distance between source and observer. The focusing theorem does not contradict flux conservation, because the mean total magnification (or amplification) directly corresponds to different areas of the source (or observer) sphere in the lensed and unlensed situation. We argue that a constant affine distance does not define an astronomically meaningful reference.
By exchanging source and observer, we confirm that magnification and amplification differ according to Etherington's reciprocity law, so that surface brightness is no longer strictly conserved. At this level we also have to distinguish between different surface brightness definitions that are based on photon number, photon flux and energy flux.     

The exact transformation from spherical harmonic to ellipsoidal harmonic coefficients for gravitational field modeling

The spherical and ellipsoidal harmonic series of the external gravitational potential for a given mass distribution are equivalent in their mutual region of uniform convergence. In an instructive case, the equality of the two series on the common coordinate surface of an infinitely large sphere reveals the exact correspondence between the spherical and ellipsoidal harmonic coefficients. The transformation between the two sets of coefficients can be accomplished via the numerical methods by Walter (Celest Mech 2:389–397, 1970) and Dechambre and Scheeres (Astron Astrophys 387:1114–1122, 2002), respectively. On the other hand, the harmonic coefficients are defined by the integrals of mass density moments in terms of the respective solid harmonics. This paper presents general algebraic formulas for expressing the solid ellipsoidal harmonics as a linear combination of the corresponding solid spherical harmonics. An exact transformation from spherical to ellipsoidal harmonic coefficients is found by incorporating these connecting expressions into the density integral. A computational procedure is proposed for the transformation. Numerical results based on the nearly ellipsoidal Martian moon, Phobos, are presented for validation of the method.     

Measuring statistical isotropy of CMB anisotropy

The statistical expectation values of the temperature fluctuations and polarization of cosmic microwave background (CMB) are assumed to be preserved under rotations of the sky. We investigate the statistical isotropy (SI) of the CMB maps recently measured by the Wilkinson microwave anisotropy probe (WMAP) using the bipolar spherical harmonic formalism proposed in Hajian and Souradeep [Hajian, A., Souradeep, T. (2003) Astrophys. J. Lett. 597, L5] for CMB temperature anisotropy and extended to CMB polarization in Basak, Hajian and Souradeep [Basak, S., Hajian, A., Souradeep, T. (2006) Phys. Rev. D74, 02130(R)]. The Bipolar Power Spectrum (BiPS) had been measured for the full sky CMB anisotropy maps of the first year WMAP data and now for the recently released three years of WMAP data. We also introduce and measure directional sensitive reduced Bipolar coefficients on the three year WMAP ILC map. Consistent with our published results from first year WMAP data we have no evidence for violation of statistical isotropy on large angular scales. Preliminary analysis of the recently released first WMAP polarization maps, however, indicate significant violation of SI even when the foreground contaminated regions are masked out. Further work is required to confirm a possible cosmic origin and rule out the (more likely) origin in observational artifact such as foreground residuals at high galactic latitude.     

On X-ray variability in narrow-line and broad-line active galactic nuclei

We propose a novel mathematical method to construct an exact polytropic sphere in self-gravitating hydrostatic equilibrium, improving the non-linear Poisson equation. The central boundary condition for the present equation requires a ratio of gas pressure to total one at the centre, which is uniquely identified by the whole mass and molecular weight of the system. The special solution derived from the Lane–Emden equation can be reproduced. This scheme is now available for modelling the molecular cloud cores in interstellar media. The mass–radius relation of the first core is found to be consistent with the recent results of radiation hydrodynamic simulations.     

Mapping on the HEALPix grid

The natural spherical projection associated with the Hierarchical Equal-Area and Isolatitude Pixelization (HEALPix) is described and shown to be one of a hybrid class that combines the cylindrical equal-area and Collignon projections, not previously documented in the cartographic literature. Projection equations are derived for the class in general and are used to investigate its properties. It is shown that the HEALPix projection suggests a simple method of (i) storing and (ii) visualizing data sampled on the grid of the HEALPix pixelization, and also suggests an extension of the pixelization that is better suited for these purposes. Potentially useful properties of other members of the class are described, and new triangular and hexagonal pixelizations are constructed from them. Finally, the standard formalism is defined for representing the celestial coordinate system for any member of the class in the FITS data format.     

Frame transformation of spherical-harmonic coefficients of differential particle intensity

We present a general second-order-correct frame transformation on spherical-harmonic coefficients of differential particle intensity. The transformation, valid for relativistic particles as well, provides a clear view of the Compton-Getting effect. It shows explicitly how each transformed harmonic coefficient depends on a subset of the original harmonic coefficients. The general expression for the first-order Compton-Getting vector anisotropy is derived and interpreted. In addition, we show how the new transformation allows one to simplify a current procedure for determining the directional intensity in a comoving frame. This involves the directional particle data measured on a spacecraft.     

A secondary resonance in Mercury’s rotation

The resonant rotation of Mercury can be modelised by a kernel model on which we can add perturbations. Our kernel model is a two-degree of freedom one written in Hamiltonian formalism. For this kernel, we consider that Mercury is solid and rotates on a Keplerian orbit. By introducing the perturbations due to the other planets of the Solar System, it appears that, in a particular case, our slow degree of freedom may enter into a 1:1 resonance with the Great Inequality of Jupiter and Saturn. Actually, as the moments of inertia of Mercury are still poorly known, this phenomenon cannot be excluded.     

Small body surface gravity fields via spherical harmonic expansions

Conventional gravity field expressions are derived from Laplace’s equation, the result being the spherical harmonic gravity field. This gravity field is said to be the exterior spherical harmonic gravity field, as its convergence region is outside the Brillouin (i.e., circumscribing) sphere of the body. In contrast, there exists its counterpart called the interior spherical harmonic gravity field for which the convergence region lies within the interior Brillouin sphere that is not the same as the exterior Brillouin sphere. Thus, the exterior spherical harmonic gravity field cannot model the gravitation within the exterior Brillouin sphere except in some special cases, and the interior spherical harmonic gravity field cannot model the gravitation outside the interior Brillouin sphere. In this paper, we will discuss two types of other spherical harmonic gravity fields that bridge the null space of the exterior/interior gravity field expressions by solving Poisson’s equation. These two gravity fields are obtained by assuming the form of Helmholtz’s equation to Poisson’s equation. This method renders the gravitational potentials as functions of spherical Bessel functions and spherical harmonic coefficients. We refer to these gravity fields as the interior/exterior spherical Bessel gravity fields and study their characteristics. The interior spherical Bessel gravity field is investigated in detail for proximity operation purposes around small primitive bodies. Particularly, we apply the theory to asteroids Bennu (formerly 1999 RQ36) and Castalia to quantify its performance around both nearly spheroidal and contact-binary asteroids, respectively. Furthermore, comparisons between the exterior gravity field, interior gravity field, interior spherical Bessel gravity field, and polyhedral gravity field are made and recommendations are given in order to aid planning of proximity operations for future small body missions.     

Spherical needlets for cosmic microwave background data analysis

We discuss spherical needlets and their properties. Needlets are a form of spherical wavelets which do not rely on any kind of tangent plane approximation and enjoy good localization properties in both pixel and harmonic space; moreover needlet coefficients are asymptotically uncorrelated at any fixed angular distance, which makes their use in statistical procedures very promising. In view of these properties, we believe needlets may turn out to be especially useful in the analysis of cosmic microwave background (CMB) data on the incomplete sky, as well as of other cosmological observations. As a final advantage, we stress that the implementation of needlets is computationally very convenient and may rely completely on standard data analysis packages such as healp ix.     

An energy-conserving formalism for adaptive gravitational force softening in smoothed particle hydrodynamics and N-body codes

In this paper, we describe an adaptive softening length formalism for collisionless N -body and self-gravitating smoothed particle hydrodynamics (SPH) calculations which conserves momentum and energy exactly. This means that spatially variable softening lengths can be used in N -body calculations without secular increases in energy. The formalism requires the calculation of a small additional term to the gravitational force related to the gradient of the softening length. The extra term is similar in form to the usual SPH pressure force (although opposite in direction) and is therefore straightforward to implement in any SPH code at almost no extra cost. For N -body codes, some additional cost is involved as the formalism requires the computation of the density through a summation over neighbouring particles using the smoothing kernel. The results of numerical tests demonstrate that, for homogeneous mass distributions, the use of adaptive softening lengths gives a softening which is always close to the 'optimal' choice of fixed softening parameter, removing the need for fine-tuning. For a heterogeneous mass distribution (as may be found in any large-scale N -body simulation), we find that the errors on the least-dense component are lowered by an order of magnitude compared to the use of a fixed softening length tuned to the densest component. For SPH codes, our method presents a natural and an elegant choice of softening formalism which makes a small improvement to both the force resolution and the total energy conservation at almost zero additional cost.     

Resistance of motion to a small, hypervelocity sphere, sputtering through a gas

Earlier work on the resistance acting on a small sphere moving through a gas is reviewed. A model for the resistance encountered by a sphere, the surface molecules of which are sputtered off during collisions with the gas molecules, is derived and compared with the case of specular reflection. The sputtering model is applied to the case of small 10-μm radius meteoroids entering the Earth's atmosphere. A possible link between the results obtained and the recent discovery of unheated, organic grains at an altitude of 40 km in the Earth's atmosphere is considered.     

Spherical Magnetic Vortex in an External Potential Field: A Dissipative Contraction

We consider the dissipative evolution of a spherical magnetic vortex with a force-free internal structure, located in a resistive medium and held in equilibrium by the potential external field. The magnetic field inside the sphere is force-free (the model of Chandrasekhar in Proc. Natl. Acad. Sci. 42, 1, 1956). Topologically, it is a set of magnetic toroids enclosed in spherical layers. A new exact MHD solution has been derived, describing a slow, uniform, radial compression of a magnetic spheroid under the pressure of an ambient field, when the plasma density and pressure are growing inside it. There is no dissipation in the potential field outside the sphere, but inside the sphere, where the current density can be high enough, the magnetic energy is continuously converted into heat. Joule dissipation lowers the magnetic pressure inside the sphere, which balances the pressure of the ambient field. This results in radial contraction of the magnetic sphere with a speed defined by the conductivity of the plasma and the characteristic spatial scale of the magnetic field inside the sphere. Formally, the sphere shrinks to zero within a finite time interval (magnetic collapse). The time of compression can be relatively small, within a day, even for a sphere with a radius of about 1 Mm, if the magnetic helicity trapped initially in the sphere (which is proportional to the number of magnetic toroids in the sphere) is quite large. The magnetic system is open along its axis of symmetry. On this axis, the magnetic and electric fields are strictly radial and sign-variable along the radius, so the plasma will be ejected along the axis of magnetic sphere outwards in both directions (as jets) at a rate much higher than the diffusive one, and the charged particles will be accelerated unevenly, in spurts, creating quasi-regular X-ray spikes. The applications of the solution to solar flares are discussed.     

黑洞候选体MAXI J1535-571爆发期间的时变现象

利用“慧眼”(Hard X-ray Modulation Telescope, Insight-HXMT)卫星在2017年9月对黑洞候选体MAXI J1535-571的观测数据,研究了该源在爆发期内的时变现象.当源处于不同的爆发谱态时,功率密度谱的谱型存在明显差异.在硬中间态,有明显的限带噪声(band-limited noise)成分和QPO (Quasi-Periodic Oscillation)成分.分析结果表明:低频限带噪声的特征频率随能量的变化呈现正相关,即软能段光子的特征频率小于硬能段光子的特征频率. 0.1–0.5 Hz频率区间的限带噪声RMS (Root Mean Square)谱在硬中间态和软中间态均出现峰值,且在高能端存在差异,可能是主导噪声RMS的能谱成分占比不同.当谱态由硬中间态过渡到软中间态时, C型QPO的RMS谱保持相似趋势,但限带噪声RMS谱存在谱态依赖现象,暗示着噪声和QPO有不同的起源机制.     

Saturated magnetic field amplification at supernova shocks

Cosmic ray streaming instabilities at supernova shocks are discussed in the quasi-linear diffusion formalism which takes into account the feedback effect of wave growth on the cosmic ray streaming motion. In particular, the non-resonant instability that leads to magnetic field amplification in the short wavelength regime is considered. The linear growth rate is calculated using kinetic theory for a streaming distribution. We show that the non-resonant instability is actually driven by a compensating current in the background plasma. The non-resonant instability can develop into a non-linear regime generating turbulence. The saturation of the amplified magnetic fields due to particle diffusion in the turbulence is derived analytically. It is shown that the evolution of parallel and perpendicular cosmic ray pressures is predominantly determined by non-resonant diffusion. However, the saturation is determined by resonant diffusion which tends to reduce the streaming motion through pitch angle scattering. The saturated level can exceed the mean background magnetic field.     

A method for aligning the plastic scintillator detector on DAMPE

The Plastic Scintillator Detector(PSD) onboard the DArk Matter Particle Explorer(DAMPE)is designed to measure cosmic ray charge(Z) and to act as a veto detector for gamma ray identification.To fully exploit the charge identification potential of PSD and to enhance its capability to identify gamma ray events, we develop an alignment method for the PSD. The path length of a given track in the volume of a PSD bar is derived taking into account the shift and rotation alignment corrections. By examining energy spectra of corner-passing events and fully contained events, position shifts and rotations of all PSD bars are obtained, and are found to be on average about 1 mm and 0.0015 radian respectively. To validate the alignment method, we introduce artificial shifts and rotations of PSD bars into the detector simulation.These shift and rotation parameters can be recovered successfully by the alignment procedure. As a result of the PSD alignment procedure, the charge resolution of the PSD is improved from 4% to 8%, depending on the nuclei.