Algorithm 1018: FaVeST—Fast Vector Spherical Harmonic Transforms

Vector spherical harmonics on the unit sphere of ℝ 3 have broad applications in geophysics, quantum mechanics, and astrophysics. In the representation of a tangent vector field, one needs to evaluate the expansion and the Fourier coefficients of vector spherical harmonics. In this article, we develop fast algorithms (FaVeST) for vector spherical harmonic transforms on these evaluations. The forward FaVeST evaluates the Fourier coefficients and has a computational cost proportional to N log √ N for N number of evaluation points. The adjoint FaVeST, which evaluates a linear combination of vector spherical harmonics with a degree up to ⊡ M for M evaluation points, has cost proportional to M log √ M . Numerical examples of simulated tangent fields illustrate the accuracy, efficiency, and stability of FaVeST.

Base Point Freeness, Uniqueness of Decompositions and Double Points for Veronese and Segre Varieties

We prove a base point freeness result for linear systems of forms vanishing at general double points of the projective plane. For tensors we study the uniqueness problem for the representation of a tensor as a sum of terms corresponding to points and tangent vectors of the Segre variety associated with the format of the tensor. We give complete results for unions of one point and one tangent vector.

An application of symplectic integration for general relativistic planetary orbitography subject to non-gravitational forces

Spacecraft propagation tools describe the motion of near-Earth objects and interplanetary probes using Newton’s theory of gravity supplemented with the approximate general relativistic n-body Einstein–Infeld–Hoffmann equations of motion. With respect to the general theory of relativity and the long-standing recommendations of the International Astronomical Union for astrometry, celestial mechanics and metrology, we believe modern orbitography software is now reaching its limits in terms of complexity. In this paper, we present the first results of a prototype software titled General Relativistic Accelerometer-based Propagation Environment (GRAPE). We describe the motion of interplanetary probes and spacecraft using extended general relativistic equations of motion which account for non-gravitational forces using end-user supplied accelerometer data or approximate dynamical models. We exploit the unique general relativistic quadratic invariant associated with the orthogonality between four-velocity and acceleration and simulate the perturbed orbits for Molniya, Parker Solar Probe and Mercury Planetary Orbiter-like test particles subject to a radiation-like four-force. The accuracy of the numerical procedure is maintained using a 5-stage, $$10^\mathrm{th}$$ 10 th -order structure-preserving Gauss collocation symplectic integration scheme. GRAPE preserves the norm of the tangent vector to the test particle worldline at the order of $$10^{-32}$$ 10 - 32 .

The pseudo-null geometric phase along optical fiber

In this study, a pseudo-null space curve in Minkowski 3-space is used to describe an optical fiber that is injected into monochromatic linear polarized light. The direction of the electric field vector with respect to the Frenet frame of a pseudo-null curve determines the state polarization of a monochromatic linearly polarized light wave traveling along an optical fiber. For the Frenet frame of a pseudo-null curve in Minkowski 3-space, the polarization vector [Formula: see text] is assumed to be perpendicular to the tangent vector [Formula: see text] with respect to anholonomic coordinates. Anholonomic coordinates for the Frenet frame of a pseudo-null curve are used to describe pseudo-null electromagnetic curves in the normal and binormal directions along an optical fiber. For the Frenet frame of the pseudo-null curve, Lorentz force equations in the normal and binormal directions along the optical fiber are presented. Pseudo-normal and binormal Rytov parallel transport laws for electric fields in the normal and binormal directions along with the optical fiber for the Frenet frame of the pseudo-null curve via anholonomic coordinates are presented. For anholonomic coordinates in Minkowski 3-space, rotations of the polarization planes of a light wave traveling in the normal and binormal directions along with the optical fiber with respect to the Frenet frame of the pseudo-null curve are obtained. Finally, a pseudo-null curve’s Maxwellian evolution is determined.

Dimensionality Reduction of SPD Data Based on Riemannian Manifold Tangent Spaces and Isometry

Symmetric positive definite (SPD) data have become a hot topic in machine learning. Instead of a linear Euclidean space, SPD data generally lie on a nonlinear Riemannian manifold. To get over the problems caused by the high data dimensionality, dimensionality reduction (DR) is a key subject for SPD data, where bilinear transformation plays a vital role. Because linear operations are not supported in nonlinear spaces such as Riemannian manifolds, directly performing Euclidean DR methods on SPD matrices is inadequate and difficult in complex models and optimization. An SPD data DR method based on Riemannian manifold tangent spaces and global isometry (RMTSISOM-SPDDR) is proposed in this research. The main contributions are listed: (1) Any Riemannian manifold tangent space is a Hilbert space isomorphic to a Euclidean space. Particularly for SPD manifolds, tangent spaces consist of symmetric matrices, which can greatly preserve the form and attributes of original SPD data. For this reason, RMTSISOM-SPDDR transfers the bilinear transformation from manifolds to tangent spaces. (2) By log transformation, original SPD data are mapped to the tangent space at the identity matrix under the affine invariant Riemannian metric (AIRM). In this way, the geodesic distance between original data and the identity matrix is equal to the Euclidean distance between corresponding tangent vector and the origin. (3) The bilinear transformation is further determined by the isometric criterion guaranteeing the geodesic distance on high-dimensional SPD manifold as close as possible to the Euclidean distance in the tangent space of low-dimensional SPD manifold. Then, we use it for the DR of original SPD data. Experiments on five commonly used datasets show that RMTSISOM-SPDDR is superior to five advanced SPD data DR algorithms.

Inverse potential problems in divergence form for measures in the plane

We study inverse potential problems with source term the divergence of some unknown (R 3 -valued) measure supported in a plane; e.g., inverse magnetization problems for thin plates. We investigate methods for recovering a magnetization μ by penalizing the measure-theoretic total variation norm kμk T V , and appealing to the decomposition of divergence-free measures in the plane as superpositions of unit tangent vector fields on rectifiable Jordan curves. In particular, we prove for magnetizations supported in a plane that T V -regularization schemes always have a unique minimizer, even in the presence of noise. It is further shown that T V -norm minimization (among magnetizations generating the same field) uniquely recovers planar magnetizations in the following two cases: (i) when the magnetization is carried by a collection of sufficiently separated line segments and a set that is purely 1-unrectifiable; (ii) when a superset of the support is tree- like. We note that such magnetizations can be recovered via T V -regularization schemes in the zero noise limit by taking the regularization parameter to zero. This suggests definitions of sparsity in the present infinite dimensional context, that generate results akin to compressed sensing.

Interior product, Lie derivative and Wilson line in the KBc subsector of open string field theory

The open string field theory of Witten (SFT) has a close formal similarity with Chern-Simons theory in three dimensions. This similarity is due to the fact that the former theory has concepts corresponding to forms, exterior derivative, wedge product and integration over the manifold. In this paper, we introduce the interior product and the Lie derivative in the KBc subsector of SFT. The interior product in SFT is specified by a two-component “tangent vector” and lowers the ghost number by one (like the ordinary interior product maps a p-form to (p − 1)-form). The Lie derivative in SFT is defined as the anti-commutator of the interior product and the BRST operator. The important property of these two operations is that they respect the KBc algebra.Deforming the original (K, B, c) by using the Lie derivative, we can consider an infinite copies of the KBc algebra, which we call the KBc manifold. As an application, we construct the Wilson line on the manifold, which could play a role in reproducing degenerate fluctuation modes around a multi-brane solution.

Twisted submanifolds of $${\mathbb {R}}^n$$

We propose a general procedure to construct noncommutative deformations of an embedded submanifold M of $${\mathbb {R}}^n$$ R n determined by a set of smooth equations $$f^a(x)=0$$ f a ( x ) = 0 . We use the framework of Drinfel’d twist deformation of differential geometry of Aschieri et al. (Class Quantum Gravity 23:1883, 2006); the commutative pointwise product is replaced by a (generally noncommutative) $$\star $$ ⋆ -product determined by a Drinfel’d twist. The twists we employ are based on the Lie algebra $$\Xi _t$$ Ξ t of vector fields that are tangent to all the submanifolds that are level sets of the $$f^a$$ f a (tangent infinitesimal diffeomorphisms); the twisted Cartan calculus is automatically equivariant under twisted $$\Xi _t$$ Ξ t . We can consistently project a connection from the twisted $${\mathbb {R}}^n$$ R n to the twisted M if the twist is based on a suitable Lie subalgebra $${\mathfrak {e}}\subset \Xi _t$$ e ⊂ Ξ t . If we endow $${\mathbb {R}}^n$$ R n with a metric, then twisting and projecting to the normal and tangent vector fields commute, and we can project the Levi–Civita connection consistently to the twisted M, provided the twist is based on the Lie subalgebra $${\mathfrak {k}}\subset {\mathfrak {e}}$$ k ⊂ e of the Killing vector fields of the metric; a twisted Gauss theorem follows, in particular. Twisted algebraic manifolds can be characterized in terms of generators and $$\star $$ ⋆ -polynomial relations. We present in some detail twisted cylinders embedded in twisted Euclidean $${\mathbb {R}}^3$$ R 3 and twisted hyperboloids embedded in twisted Minkowski $${\mathbb {R}}^3$$ R 3 [these are twisted (anti-)de Sitter spaces $$dS_2,AdS_2$$ d S 2 , A d S 2 ].