No Arabic abstract
The discrete prolate spheroidal sequences (DPSSs) are a set of orthonormal sequences in $ell_2(mathbb{Z})$ which are strictly bandlimited to a frequency band $[-W,W]$ and maximally concentrated in a time interval ${0,ldots,N-1}$. The timelimited DPSSs (sometimes referred to as the Slepian basis) are an orthonormal set of vectors in $mathbb{C}^N$ whose discrete time Fourier transform (DTFT) is maximally concentrated in a frequency band $[-W,W]$. Due to these properties, DPSSs have a wide variety of signal processing applications. The DPSSs are the eigensequences of a timelimit-then-bandlimit operator and the Slepian basis vectors are the eigenvectors of the so-called prolate matrix. The eigenvalues in both cases are the same, and they exhibit a particular clustering behavior -- slightly fewer than $2NW$ eigenvalues are very close to $1$, slightly fewer than $N-2NW$ eigenvalues are very close to $0$, and very few eigenvalues are not near $1$ or $0$. This eigenvalue behavior is critical in many of the applications in which DPSSs are used. There are many asymptotic characterizations of the number of eigenvalues not near $0$ or $1$. In contrast, there are very few non-asymptotic results, and these dont fully characterize the clustering behavior of the DPSS eigenvalues. In this work, we establish two novel non-asymptotic bounds on the number of DPSS eigenvalues between $epsilon$ and $1-epsilon$. Also, we obtain bounds detailing how close the first $approx 2NW$ eigenvalues are to $1$ and how close the last $approx N-2NW$ eigenvalues are to $0$. Furthermore, we extend these results to the eigenvalues of the prolate spheroidal wave functions (PSWFs), which are the continuous-time version of the DPSSs. Finally, we present numerical experiments demonstrating the quality of our non-asymptotic bounds on the number of DPSS eigenvalues between $epsilon$ and $1-epsilon$.
Prolate spheroidal wave functions provide a natural and effective tool for computing with bandlimited functions defined on an interval. As demonstrated by Slepian et al., the so called generalized prolate spheroidal functions (GPSFs) extend this apparatus to higher dimensions. While the analytical and numerical apparatus in one dimension is fairly complete, the situation in higher dimensions is less satisfactory. This report attempts to improve the situation by providing analytical and numerical tools for GPSFs, including the efficient evaluation of eigenvalues, the construction of quadratures, interpolation formulae, etc. Our results are illustrated with several numerical examples.
In order to produce high dynamic range images in radio interferometry, bright extended sources need to be removed with minimal error. However, this is not a trivial task because the Fourier plane is sampled only at a finite number of points. The ensuing deconvolution problem has been solved in many ways, mainly by algorithms based on CLEAN. However, such algorithms that use image pixels as basis functions have inherent limitations and by using an orthonormal basis that span the whole image, we can overcome them. The construction of such an orthonormal basis involves fine tuning of many free parameters that define the basis functions. The optimal basis for a given problem (or a given extended source) is not guaranteed. In this paper, we discuss the use of generalized prolate spheroidal wave functions as a basis. Given the geometry (or the region of interest) of an extended source and the sampling points on the visibility plane, we can construct the optimal basis to model the source. Not only does this gives us the minimum number of basis functions required but also the artifacts outside the region of interest are minimized.
Motivated by the discovery of prolate rotation of stars in Andromeda II, a dwarf spheroidal companion of M31, we study its origin via mergers of disky dwarf galaxies. We simulate merger events between two identical dwarfs changing the initial inclination of their disks with respect to the orbit and the amount of orbital angular momentum. On radial orbits the amount of prolate rotation in the merger remnants correlates strongly with the inclination of the disks and is well understood as due to the conservation of the angular momentum component of the disks along the merger axis. For non-radial orbits prolate rotation may still be produced if the orbital angular momentum is initially not much larger than the intrinsic angular momentum of the disks. The orbital structure of the remnants with significant rotation is dominated by box orbits in the center and long-axis tubes in the outer parts. The frequency analysis of stellar orbits in the plane perpendicular to the major axis reveals the presence of two families roughly corresponding to inner and outer long-axis tubes. The fraction of inner tubes is largest in the remnant forming from disks oriented most vertically initially and is responsible for the boxy shape of the galaxy. We conclude that prolate rotation results from mergers with a variety of initial conditions and no fine tuning is necessary to reproduce this feature. We compare the properties of our merger remnants to those of dwarfs resulting from the tidal stirring scenario and the data for Andromeda II.
We have carried out calculations of the triple-differential cross section for one-photon double ionization of molecular hydrogen for a central photon energy of $75$~eV, using a fully {it ab initio}, nonperturbative approach to solve the time-dependent Schro equation in prolate spheroidal coordinates. The spatial coordinates $xi$ and $eta$ are discretized in a finite-element discrete-variable representation. The wave packet of the laser-driven two-electron system is propagated in time through an effective short iterative Lanczos method to simulate the double ionization of the hydrogen molecule. For both symmetric and asymmetric energy sharing, the present results agree to a satisfactory level with most earlier predictions for the absolute magnitude and the shape of the angular distributions. A notable exception, however, concerns the predictions of the recent time-independent calculations based on the exterior complex scaling method in prolate spheroidal coordinates [Phys.~Rev.~A~{bf 82}, 023423 (2010)]. Extensive tests of the numerical implementation were performed, including the effect of truncating the Neumann expansion for the dielectronic interaction on the description of the initial bound state and the predicted cross sections. We observe that the dominant escape mode of the two photoelectrons dramatically depends upon the energy sharing. In the parallel geometry, when the ejected electrons are collected along the direction of the laser polarization axis, back-to-back escape is the dominant channel for strongly asymmetric energy sharing, while it is completely forbidden if the two electrons share the excess energy equally.
We present an approach to the design of distribution functions that depend on the phase-space coordinates through the action integrals. The approach makes it easy to construct a dynamical model of a given stellar component. We illustrate the approach by deriving distribution functions that self-consistently generate several popular stellar systems, including the Hernquist, Jaffe, and Navarro, Frenk and White models. We focus on non-rotating spherical systems, but extension to flattened and rotating systems is trivial. Our distribution functions are easily added to each other and to previously published distribution functions for discs to create self-consistent multi-component galaxies. The models this approach makes possible should prove valuable both for the interpretation of observational data and for exploring the non-equilibrium dynamics of galaxies via N-body simulation.