Finding high-order analytic post-Newtonian parameters from a high-precision numerical self-force calculation

We present a novel analytic extraction of high-order post-Newtonian (pN) parameters that govern quasi-circular binary systems. Coefficients in the pN expansion of the energy of a binary system can be found from corresponding coefficients in an extreme-mass-ratio inspiral (EMRI) computation of the change $\Delta U$ in the redshift factor of a circular orbit at fixed angular velocity. Remarkably, by computing this essentially gauge-invariant quantity to accuracy greater than one part in $10^{225}$, and by assuming that a subset of pN coefficients are rational numbers or products of $\pi$ and a rational, we obtain the exact analytic coefficients. We find the previously unexpected result that the post-Newtonian expansion of $\Delta U$ (and of the change $\Delta\Omega$ in the angular velocity at fixed redshift factor) have conservative terms at half-integral pN order beginning with a 5.5 pN term. This implies the existence of a corresponding 5.5 pN term in the expansion of the energy of a binary system. Coefficients in the pN series that do not belong to the subset just described are obtained to accuracy better than 1 part in $10^{265-23n}$ at $n$th pN order. We work in a radiation gauge, finding the radiative part of the metric perturbation from the gauge-invariant Weyl scalar $\psi_0$ via a Hertz potential. We use mode-sum renormalization, and find high-order renormalization coefficients by matching a series in $L=\ell+1/2$ to the large-$L$ behavior of the expression for $\Delta U$. The non-radiative parts of the perturbed metric associated with changes in mass and angular momentum are calculated in the Schwarzschild gauge

Covariant Uniform Acceleration

We show that standard Relativistic Dynamics Equation F=dp/d\tau is only partially covariant. To achieve full Lorentz covariance, we replace the four-force F by a rank 2 antisymmetric tensor acting on the four-velocity. By taking this tensor to be constant, we obtain a covariant definition of uniformly accelerated motion. We compute explicit solutions for uniformly accelerated motion which are divided into four types: null, linear, rotational, and general. For null acceleration, the worldline is cubic in the time. Linear acceleration covariantly extends 1D hyperbolic motion, while rotational acceleration covariantly extends pure rotational motion. We use Generalized Fermi-Walker transport to construct a uniformly accelerated family of inertial frames which are instantaneously comoving to a uniformly accelerated observer. We explain the connection between our approach and that of Mashhoon. We show that our solutions of uniformly accelerated motion have constant acceleration in the comoving frame. Assuming the Weak Hypothesis of Locality, we obtain local spacetime transformations from a uniformly accelerated frame K' to an inertial frame K. The spacetime transformations between two uniformly accelerated frames with the same acceleration are Lorentz. We compute the metric at an arbitrary point of a uniformly accelerated frame. We obtain velocity and acceleration transformations from a uniformly accelerated system K' to an inertial frame K. We derive the general formula for the time dilation between accelerated clocks. We obtain a formula for the angular velocity of a uniformly accelerated object. Every rest point of K' is uniformly accelerated, and its acceleration is a function of the observer's acceleration and its position. We obtain an interpretation of the Lorentz-Abraham-Dirac equation as an acceleration transformation from K' to K.Comment: 36 page

Nonaxisymmetric Neutral Modes in Rotating Relativistic Stars

We study nonaxisymmetric perturbations of rotating relativistic stars. modeled as perfect-fluid equilibria. Instability to a mode with angular dependence $\exp(im\phi)$ sets in when the frequency of the mode vanishes. The locations of these zero-frequency modes along sequences of rotating stars are computed in the framework of general relativity. We consider models of uniformly rotating stars with polytropic equations of state, finding that the relativistic models are unstable to nonaxisymmetric modes at significantly smaller values of rotation than in the Newtonian limit. Most strikingly, the m=2 bar mode can become unstable even for soft polytropes of index $N \leq 1.3$, while in Newtonian theory it becomes unstable only for stiff polytropes of index $N \leq 0.808$. If rapidly rotating neutron stars are formed by the accretion-induced collapse of white dwarfs, instability associated with these nonaxisymmetric, gravitational-wave driven modes may set an upper limit on neutron-star rotation. Consideration is restricted to perturbations that correspond to polar perturbations of a spherical star. A study of axial perturbations is in progress.Comment: 57 pages, 9 figure

Quantum Lattice Fluctuations and Luminescence in C_60

We consider luminescence in photo-excited neutral C_60 using the Su-Schrieffer-Heeger model applied to a single C_60 molecule. To calculate the luminescence we use a collective coordinate method where our collective coordinate resembles the displacement of the carbon atoms of the Hg(8) phonon mode and extrapolates between the ground state "dimerisation" and the exciton polaron. There is good agreement for the existing luminescence peak spacing and fair agreement for the relative intensity. We predict the existence of further peaks not yet resolved in experiment. PACS Numbers : 78.65.Hc, 74.70.Kn, 36.90+

TreeGrad: Transferring Tree Ensembles to Neural Networks

Gradient Boosting Decision Tree (GBDT) are popular machine learning algorithms with implementations such as LightGBM and in popular machine learning toolkits like Scikit-Learn. Many implementations can only produce trees in an offline manner and in a greedy manner. We explore ways to convert existing GBDT implementations to known neural network architectures with minimal performance loss in order to allow decision splits to be updated in an online manner and provide extensions to allow splits points to be altered as a neural architecture search problem. We provide learning bounds for our neural network.Comment: Technical Report on Implementation of Deep Neural Decision Forests Algorithm. To accompany implementation here: https://github.com/chappers/TreeGrad. Update: Please cite as: Siu, C. (2019). "Transferring Tree Ensembles to Neural Networks". International Conference on Neural Information Processing. Springer, 2019. arXiv admin note: text overlap with arXiv:1909.1179

The Quantum Propagator for a Nonrelativistic Particle in the Vicinity of a Time Machine

We study the propagator of a non-relativistic, non-interacting particle in any non-relativistic ``time-machine'' spacetime of the type shown in Fig.~1: an external, flat spacetime in which two spatial regions, $V_-$ at time $t_-$ and $V_+$ at time $t_+$, are connected by two temporal wormholes, one leading from the past side of $V_-$ to t the future side of $V_+$ and the other from the past side of $V_+$ to the future side of $V_-$. We express the propagator explicitly in terms of those for ordinary, flat spacetime and for the two wormholes; and from that expression we show that the propagator satisfies completeness and unitarity in the initial and final ``chronal regions'' (regions without closed timelike curves) and its propagation from the initial region to the final region is unitary. However, within the time machine it satisfies neither completeness nor unitarity. We also give an alternative proof of initial-region-to-final-region unitarity based on a conserved current and Gauss's theorem. This proof can be carried over without change to most any non-relativistic time-machine spacetime; it is the non-relativistic version of a theorem by Friedman, Papastamatiou and Simon, which says that for a free scalar field, quantum mechanical unitarity follows from the fact that the classical evolution preserves the Klein-Gordon inner product

Path Integrals, Density Matrices, and Information Flow with Closed Timelike Curves

Two formulations of quantum mechanics, inequivalent in the presence of closed timelike curves, are studied in the context of a soluable system. It illustrates how quantum field nonlinearities lead to a breakdown of unitarity, causality, and superposition using a path integral. Deutsch's density matrix approach is causal but typically destroys coherence. For each of these formulations I demonstrate that there are yet further alternatives in prescribing the handling of information flow (inequivalent to previous analyses) that have implications for any system in which unitarity or coherence are not preserved.Comment: 25 pages, phyzzx, CALT-68-188

Models of helically symmetric binary systems

Results from helically symmetric scalar field models and first results from a convergent helically symmetric binary neutron star code are reported here; these are models stationary in the rotating frame of a source with constant angular velocity omega. In the scalar field models and the neutron star code, helical symmetry leads to a system of mixed elliptic-hyperbolic character. The scalar field models involve nonlinear terms that mimic nonlinear terms of the Einstein equation. Convergence is strikingly different for different signs of each nonlinear term; it is typically insensitive to the iterative method used; and it improves with an outer boundary in the near zone. In the neutron star code, one has no control on the sign of the source, and convergence has been achieved only for an outer boundary less than approximately 1 wavelength from the source or for a code that imposes helical symmetry only inside a near zone of that size. The inaccuracy of helically symmetric solutions with appropriate boundary conditions should be comparable to the inaccuracy of a waveless formalism that neglects gravitational waves; and the (near zone) solutions we obtain for waveless and helically symmetric BNS codes with the same boundary conditions nearly coincide.Comment: 19 pages, 7 figures. Expanded version of article to be published in Class. Quantum Grav. special issue on Numerical Relativit

Singularities of Nonlinear Elliptic Systems

Through Morrey's spaces (plus Zorko's spaces) and their potentials/capacities as well as Hausdorff contents/dimensions, this paper estimates the singular sets of nonlinear elliptic systems of the even-ordered Meyers-Elcrat type and a class of quadratic functionals inducing harmonic maps.Comment: 18 pages Communications in Partial Differential Equation