skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Riesz energy, L2$L^2$ discrepancy, and optimal transport of determinantal point processes on the sphere and the flat torus
Abstract Determinantal point processes exhibit an inherent repulsive behavior, thus providing examples of very evenly distributed point sets on manifolds. In this paper, we study the so‐called harmonic ensemble, defined in terms of Laplace eigenfunctions on the sphere and the flat torus , and the so‐called spherical ensemble on , which originates in random matrix theory. We extend results of Beltrán, Marzo, and Ortega‐Cerdà on the Riesz ‐energy of the harmonic ensemble to the nonsingular regime , and as a corollary find the expected value of the spherical cap discrepancy via the Stolarsky invariance principle. We find the expected value of the discrepancy with respect to axis‐parallel boxes and Euclidean balls of the harmonic ensemble on . We also show that the spherical ensemble and the harmonic ensemble on and with points attain the optimal rate in expectation in the Wasserstein metric , in contrast to independent and identically distributed random points, which are known to lose a factor of .  more » « less
Award ID(s):
2202877
PAR ID:
10500211
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Oxford University Press (OUP)
Date Published:
Journal Name:
Mathematika
Volume:
70
Issue:
2
ISSN:
0025-5793
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Consider two half-spaces$$H_1^+$$ H 1 + and$$H_2^+$$ H 2 + in$${\mathbb {R}}^{d+1}$$ R d + 1 whose bounding hyperplanes$$H_1$$ H 1 and$$H_2$$ H 2 are orthogonal and pass through the origin. The intersection$${\mathbb {S}}_{2,+}^d:={\mathbb {S}}^d\cap H_1^+\cap H_2^+$$ S 2 , + d : = S d H 1 + H 2 + is a spherical convex subset of thed-dimensional unit sphere$${\mathbb {S}}^d$$ S d , which contains a great subsphere of dimension$$d-2$$ d - 2 and is called a spherical wedge. Choosenindependent random points uniformly at random on$${\mathbb {S}}_{2,+}^d$$ S 2 , + d and consider the expected facet number of the spherical convex hull of these points. It is shown that, up to terms of lower order, this expectation grows like a constant multiple of$$\log n$$ log n . A similar behaviour is obtained for the expected facet number of a homogeneous Poisson point process on$${\mathbb {S}}_{2,+}^d$$ S 2 , + d . The result is compared to the corresponding behaviour of classical Euclidean random polytopes and of spherical random polytopes on a half-sphere. 
    more » « less
  2. Abstract Filamentary structures in neutral hydrogen (Hi) emission are well aligned with the interstellar magnetic field, so Hiemission morphology can be used to construct templates that strongly correlate with measurements of polarized thermal dust emission. We explore how the quantification of filament morphology affects this correlation. We introduce a new implementation of the Rolling Hough Transform (RHT) using spherical harmonic convolutions, which enables efficient quantification of filamentary structure on the sphere. We use this Spherical RHT algorithm along with a Hessian-based method to construct Hi-based polarization templates. We discuss improvements to each algorithm relative to similar implementations in the literature and compare their outputs. By exploring the parameter space of filament morphologies with the Spherical RHT, we find that the most informative Histructures for modeling the magnetic field structure are the thinnest resolved filaments. For this reason, we find a ∼10% enhancement in theB-mode correlation with polarized dust emission with higher-resolution Hiobservations. We demonstrate that certain interstellar morphologies can produce parity-violating signatures, i.e., nonzeroTBandEB, even under the assumption that filaments are locally aligned with the magnetic field. Finally, we demonstrate thatBmodes from interstellar dust filaments are mostly affected by the topology of the filaments with respect to one another and their relative polarized intensities, whereasEmodes are mostly sensitive to the shapes of individual filaments. 
    more » « less
  3. Abstract The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational-wave background (GWB) in its 15 yr data set. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy. By modeling the angular power distribution as a sum over spherical harmonics (where the coefficients are not bound to always generate positive power everywhere), we find that the Bayesian 95% upper limit on the level of dipole anisotropy is (Cl=1/Cl=0) < 27%. This is similar to the upper limit derived under the constraint of positive power everywhere, indicating that the dipole may be close to the data-informed regime. By contrast, the constraints on anisotropy at higher spherical-harmonic multipoles are strongly prior dominated. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15 yr data set and show that this data set has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. 
    more » « less
  4. High harmonic spectroscopy utilizes the extremely nonlinear optical process of high-order harmonic generation (HHG) to measure complex attosecond-scale dynamics within the emitting atom or molecule subject to a strong laser field. However, it can be difficult to compare theory and experiment, since the dynamics under investigation are often very sensitive to the laser intensity, which inevitably varies over the Gaussian profile of a typical laser beam. This discrepancy would usually be resolved by so-called macroscopic HHG simulations, but such methods almost always use a simplified model of the internal dynamics of the molecule, which is not necessarily applicable for high harmonic spectroscopy. In this Letter, we extend the existing framework of macroscopic HHG so that high-accuracyab initiocalculations can be used as the microscopic input. This new (to the best of our knowledge) approach is applied to a recent theoretical prediction involving the HHG spectra of open-shell molecules undergoing nonadiabatic dynamics. We demonstrate that the predicted features in the HHG spectrum unambiguously survive macroscopic response calculations, and furthermore they exhibit a nontrivial angular pattern in the far field. 
    more » « less
  5. Recent research in the theory of overparametrized learning has sought to establish generalization guarantees in the interpolating regime. Such results have been established for a few common classes of methods, but so far not for ensemble methods. We devise an ensemble classification method that simultaneously interpolates the training data, and is consistent for a broad class of data distributions. To this end, we define the manifold-Hilbert kernel for data distributed on a Riemannian manifold. We prove that kernel smoothing regression and classification using the manifold-Hilbert kernel are weakly consistent in the setting of Devroye et al. [19]. For the sphere, we show that the manifold-Hilbert kernel can be realized as a weighted random partition kernel, which arises as an infinite ensemble of partition-basedclassifiers. 
    more » « less