Search

We present analytical expressions for the fluxes of energy and angular momentum from a point mass on an inclined spherical orbit about a Kerr black hole. The expressions are obtained using the method of Mano, Suzuki and Takasugi to construct analytical solutions of the Teukolsky equation, and are given as post-Newtonian expansions valid through 12PN, with arbitrary values for the inclination parameter $x$ and black hole spin $a$. We characterize the structure of the PN expansions in terms of their dependence on $x$ and $a$, and we validate our results against numerical calculations.

Wurtzite (Al,Sc)N ferroelectrics are attractive for microelectronics applications due to their chemical and epitaxial structural compatibility with wurtzite semiconductors such as GaN and (Al,Ga)N. However, the leakage current in epitaxial stacks reported to date should be reduced for reliable device operation. Following the tradition of other semiconductor heterostructures, crystalline structural quality, as measured by breadth of diffraction peaks and correlating with dislocation density, is commonly used as a proxy for leakage current, but we demonstrate here that the crystalline mosaicity that dominates the broadening of diffraction peaks in epitaxial Al$_{0.7}$Sc$_{0.3}$N stacks does not dominate leakage current. We report here well-saturated ferroelectric hysteresis loops and orders of magnitude lower leakage current (0.07 A cm$^{-2}$) compared to values reported in literature (1 ~ 19 A cm$^{-2}$) for sputter-deposited epitaxial Al$_{0.7}$Sc$_{0.3}$N/GaN of comparable crystalline quality to prior reports. Further, we show Al$_{0.7}$Sc$_{0.3}$N on lattice-matched InGaN buffers with improved structural characteristics exhibits increased leakage characteristics. This demonstration and understanding can help to guide further efforts towards reliable wurtzite ferroelectric devices and prioritize approaches targeting further leakage current reduction., Comment: 14 pages, 4 figures, 1 table

We present a numerical framework for modeling extended hyperelastic bodies based on a Lagrangian formulation of general relativistic elasticity theory. We use finite element methods to discretize the body, then use the semi--discrete action to derive ordinary differential equations of motion for the discrete nodes. The nodes are evolved in time using fourth--order Runge--Kutta. We validate our code against the normal modes of oscillation of a hyperelastic sphere, which are known analytically in the limit of small (linear), slow (Newtonian) oscillations. The algorithm displays second order convergence. This numerical framework can be used to obtain the orbital motion and internal dynamics of a hyperelastic body of any shape, for any spacetime metric, and for varying hyperelastic energy models., Comment: Corrected typos

We calculate the high-order post-Newtonian (PN) expansion of the energy and angular momentum fluxes onto the horizon of a nonspinning black hole primary in eccentric-orbit extreme-mass-ratio inspirals. The first-order black hole perturbation theory calculation uses \textsc{Mathematica} and makes an analytic expansion of the Regge-Wheeler-Zerilli functions using the Mano-Suzuki-Takasugi formalism. The horizon absorption, or tidal heating and torquing, is calculated to 18PN relative to the leading horizon flux (i.e., 22PN order relative to the leading quadrupole flux at infinity). Each PN term is a function of eccentricity $e$ and is calculated as a series to $e^{10}$. A second expansion, to 10PN horizon-relative order (or 14PN relative to the flux at infinity), is computed deeper in eccentricity to $e^{20}$. A number of resummed closed-form functions are found for the low PN terms in the series. Using a separate Teukolsky perturbation code, numerical comparisons are made to test how accurate the PN expansion is when extended to a close $p=10$ orbit. We find that the horizon absorption expansion is not as convergent as a previously computed infinity-side flux expansion. However, given that the horizon absorption is suppressed by 4PN, useful results can be obtained even with an orbit as tight as this for $e \le 1/2$. Combining the present results with our earlier expansion of the fluxes to infinity makes the knowledge of the total dissipation known to 19PN for eccentric-orbit nonspinning EMRIs., Comment: 19 pages, 1 figure

With the precision now afforded by modern space-based photometric observations from the retired K2 and current TESS missions, the effects of general relativity (GR) may be detectable in the light curves of pulsating white dwarfs (WDs). Almost all WD models are calculated using a Newtonian description of gravity and hydrodynamics. To determine if inclusion of GR leads to observable effects, we used idealized models of compact stars and made side-by-side comparison of mode periods computed using a (i) Newtonian and (ii) GR description of the equilibrium structure and nonradial pulsations. For application to white dwarfs, it is only necessary to include the first post-Newtonian (1PN) approximation to GR. The mathematical nature of the linear nonradial pulsation problem is then qualitatively unchanged and the GR corrections can be written as extensions of the classic Dziembowski equations. As such, GR effects might easily be included in existing asteroseismology codes. The idealized stellar models are (i) \pn1 relativistic polytropes and (ii) stars with cold degenerate-electron equation of state featuring a near-surface chemical transition from $\mu_e = 2$ to $\mu_e = 1$, simulating a surface hydrogen layer. Comparison of Newtonian and 1PN normal mode periods reveals fractional differences on the order of the surface gravitational redshift $z$. For a typical WD, this fractional difference is $\sim 10^{-4}$ and is greater than the period uncertainty $\sigma_{\Pi}/\Pi$ of many white dwarf pulsation modes observed by TESS. A consistent theoretical modeling of periods observed in these stars should in principle include GR effects to 1PN order.

We calculate the eccentricity dependence of the high-order post-Newtonian (PN) expansion of the spin-precession invariant $\psi$ for eccentric-orbit extreme-mass-ratio inspirals with a Schwarzschild primary. The series is calculated in first-order black hole perturbation theory through direct analytic expansion of solutions in the Regge-Wheeler-Zerilli formalism, using a code written in \textsc{Mathematica}. Modes with small values of $l$ are found via the Mano-Suzuki-Takasugi (MST) analytic function expansion formalism for solutions to the Regge-Wheeler equation. Large-$l$ solutions are found by applying a PN expansion ansatz to the Regge-Wheeler equation. Previous work has given $\psi$ to 9.5PN order and to order $e^2$ (i.e., the near circular orbit limit). We calculate the expansion to 9PN but to $e^{16}$ in eccentricity. It proves possible to find a few terms that have closed-form expressions, all of which are associated with logarithmic terms in the PN expansion. We also compare the numerical evaluation of our PN expansion to prior numerical calculations of $\psi$ in close orbits to assess its radius of convergence. We find that the series is not as rapidly convergent as the one for the redshift invariant at $r \simeq 10M$ but still yielding $\sim 1\%$ accuracy for eccentricities $e \lesssim 0.25$., Comment: 23 pages, 2 figures

We calculate the eccentricity dependence of the high-order post-Newtonian (PN) series for the generalized redshift invariant $\langle u^t \rangle_\tau$ for eccentric-orbit extreme-mass-ratio inspirals on a Schwarzschild background. These results are calculated within first-order black hole perturbation theory (BHPT) using Regge-Wheeler-Zerilli (RWZ) gauge. Our \textsc{Mathematica} code is based on a familiar procedure, using PN expansion of the Mano-Suzuki-Takasugi (MST) analytic function formalism for $l$ modes up to a certain maximum and then using a direct general-$l$ PN expansion of the RWZ equation for arbitrarily high $l$. We calculate dual expansions in PN order and in powers of eccentricity, reaching 10PN relative order and $e^{20}$. Detailed knowledge of the eccentricity expansion at each PN order allows us to find within the eccentricity dependence numerous closed-form expressions and multiple infinite series with known coefficients. We find leading logarithm sequences in the PN expansion of the redshift invariant that reflect a similar behavior in the PN expansion of the energy flux to infinity. A set of flux terms and special functions that appear in the energy flux, like the Peters-Mathews flux itself, are shown to reappear in the redshift PN expansion., Comment: 26 pages, 2 figures

In this paper, we introduce new formal methods and provide empirical evidence to highlight a unique safety concern prevalent in reinforcement learning (RL)-based recommendation algorithms -- 'user tampering.' User tampering is a situation where an RL-based recommender system may manipulate a media user's opinions through its suggestions as part of a policy to maximize long-term user engagement. We use formal techniques from causal modeling to critically analyze prevailing solutions proposed in the literature for implementing scalable RL-based recommendation systems, and we observe that these methods do not adequately prevent user tampering. Moreover, we evaluate existing mitigation strategies for reward tampering issues, and show that these methods are insufficient in addressing the distinct phenomenon of user tampering within the context of recommendations. We further reinforce our findings with a simulation study of an RL-based recommendation system focused on the dissemination of political content. Our study shows that a Q-learning algorithm consistently learns to exploit its opportunities to polarize simulated users with its early recommendations in order to have more consistent success with subsequent recommendations that align with this induced polarization. Our findings emphasize the necessity for developing safer RL-based recommendation systems and suggest that achieving such safety would require a fundamental shift in the design away from the approaches we have seen in the recent literature., Comment: In proceedings of the 6th AAAI/ACM Conference on Artificial Intelligence, Ethics and Society (AIES '23)

Extreme-mass ratio inspirals (EMRIs), compact binaries with small mass ratios $\epsilon\ll 1$, will be important sources for low-frequency gravitational wave detectors. Almost all EMRIs will evolve through important transient orbital $r\theta$ resonances, which will enhance or diminish their gravitational wave flux, thereby affecting the phase evolution of the waveforms at $O(\epsilon^{1/2})$ relative to leading order. While modeling the local gravitational self-force (GSF) during resonances is essential for generating accurate EMRI waveforms, so far the full GSF has not been calculated for an $r\theta$-resonant orbit owing to computational demands of the problem. As a first step we employ a simpler model, calculating the scalar self-force (SSF) along $r\theta$-resonant geodesics in Kerr spacetime. We demonstrate two ways of calculating the $r\theta$-resonant SSF (and likely GSF), with one method leaving the radial and polar motions initially independent as if the geodesic is nonresonant. We illustrate results by calculating the SSF along geodesics defined by three $r\theta$-resonant ratios (1:3, 1:2, 2:3). We show how the SSF and averaged evolution of the orbital constants vary with the initial phase at which an EMRI enters resonance. We then use our SSF data to test a previously proposed integrability conjecture, which argues that conservative effects vanish at adiabatic order during resonances. We find prominent contributions from the conservative SSF to the secular evolution of the Carter constant, $\langle \dot{\mathcal{Q}}\rangle$, but these nonvanishing contributions are on the order of, or less than, the estimated uncertainties of our self-force results. The uncertainties come from residual, incomplete removal of the singular field in the regularization process. Higher order regularization parameters, once available, will allow definitive tests of the integrability conjecture., Comment: 32 pages, 15 figures; Updated to reflect published version

In a recent paper we showed that for eccentric-orbit extreme-mass-ratio inspirals the analytic forms of the leading-logarithm energy and angular momentum post-Newtonian (PN) flux terms (radiated to infinity) can, to arbitrary PN order, be determined by sums over the Fourier spectrum of the Newtonian quadrupole moment. We further showed that an essential part of the eccentricity dependence of the related subleading-logarithm PN sequences, at lowest order in the symmetric mass ratio $\nu$, stems as well from the Newtonian quadrupole moment. Once that part is factored out, the remaining eccentricity dependence is more easily determined by black hole perturbation theory. In this paper we show how the sequences that are the 1PN corrections to the entire leading-logarithm series, namely terms that appear at PN orders $x^{3k+1} \log^k(x)$ and $x^{3k+5/2} \log^k(x)$ (for PN compactness parameter $x$ and integers $k\ge 0$), at lowest order in $\nu$, are determined by the Fourier spectra of the Newtonian mass octupole, Newtonian current quadrupole, and 1PN part of the mass quadrupole moments. We also develop a conjectured (but plausible) form for 1PN correction to the leading logs at second order in $\nu$. Further, in analogy to the first paper, we show that these same source multipole moments also yield nontrivial parts of the 1PN correction to the subleading-logarithm series, and that the remaining eccentricity dependence (at lowest order in $\nu$) can then more easily be determined using black hole perturbation theory. We use this method to determine the entire analytic eccentricity dependence of the perturbative (i.e., lowest order in $\nu$) 4PN non-log terms, $\mathcal{R}_4(e_t)$ and $\mathcal{Z}_4(e_t)$, for energy and angular momentum respectively., Comment: 36 pages, 1 figure

We present an extension of work in an earlier paper showing high precision comparisons between black hole perturbation theory and post-Newtonian (PN) theory in their region of overlapping validity for bound, eccentric-orbit, Schwarzschild extreme-mass-ratio inspirals. As before we apply a numerical fitting scheme to extract eccentricity coefficients in the PN expansion of the gravitational wave fluxes, which are then converted to exact analytic form using an integer-relation algorithm. In this work, however, we fit to individual $lmn$ modes to exploit simplifying factorizations that lie therein. Since the previous paper focused solely on the energy flux, here we concentrate initially on analyzing the angular momentum flux to infinity. A first step involves finding convenient forms for hereditary contributions to the flux at low-PN order, analogous to similar terms worked out previously for the energy flux. We then apply the upgraded techniques to find new PN terms through 9PN order and (at many PN orders) to $e^{30}$ in the power series in eccentricity. With the new approach applied to angular momentum fluxes, we return to the energy fluxes at infinity to extend those previous results. Like before, the underlying method uses a \textsc{Mathematica} code based on use of the Mano-Suzuki-Takasugi (MST) function expansion formalism to represent gravitational perturbations and spectral source integration (SSI) to find numerical results at arbitrarily high precision., Comment: 36 pages, 1 figure

We present new results on the analytic eccentricity dependence of several sequences of gravitational wave flux terms at high post-Newtonian (PN) order for extreme-mass-ratio inspirals. These sequences are the leading logarithms, which appear at PN orders $x^{3k} \log^k(x)$ and $x^{3k+3/2} \log^k(x)$ for integers $k\ge 0$ ($x$ a PN compactness parameter), and the subleading logarithms, which appear at orders $x^{3k} \log^{k-1}(x)$ and $x^{3k+3/2} \log^{k-1}(x)$ ($k\ge 1$), in both the energy and angular momentum radiated to infinity. For the energy flux leading logarithms, we show that to arbitrarily high PN order their eccentricity dependence is determined by particular sums over the function $g(n,e_t)$, derived from the Newtonian mass quadrupole moment, that normally gives the spectral content of the Peters-Mathews flux as a function of radial harmonic $n$. An analogous power spectrum $\tilde{g}(n,e_t)$ determines the leading logarithms of the angular momentum flux. For subleading logs, the quadrupole power spectra are again shown to play a role, providing a distinguishable part of the eccentricity dependence of these flux terms to high PN order. With the quadrupole contribution understood, the remaining analytic eccentricity dependence of the subleading logs can in principle be determined more easily using black hole perturbation theory. We show this procedure in action, deriving the complete analytic structure of the $x^6 \log(x)$ subleading-log term and an analytic expansion of the $x^{9/2}$ subleading log to high order in a power series in eccentricity. We discuss how these methods might be extended to other sequences of terms in the PN expansion involving logarithms., Comment: 24 pages. Submitted to Phys. Rev. D

We report development of a code to calculate the scalar self-force on a scalar-charged particle moving on generic bound orbits in the Kerr spacetime. The scalar self-force model allows rapid development of computational techniques relevant to generic gravitational extreme-mass-ratio inspirals (EMRIs). Our frequency-domain calculations are made with arbitrary numerical precision code written in \textsc{Mathematica}. We extend spectral source integration techniques to the Kerr spacetime, increasing computational efficiency. We model orbits with nearly arbitrary inclinations $0\leq\iota<\pi/2$ and eccentricities up to $e \lesssim 0.8$. This effort extends earlier work by Warburton and Barack where motion was restricted to the equatorial plane or to inclined spherical orbits. Consistent with a recent discovery by Thornburg and Wardell \cite{ThorWard17} in time-domain calculations, we observe self-force oscillations during the radially-outbound portion of highly eccentric orbits around a rapidly rotating black hole. As noted previously, these oscillations reflect coupling into the self-force by quasinormal modes excited during pericenter passage. Our results confirm the effect with a frequency-domain code. \emph{More importantly, we find that quasinormal bursts (QNBs) appear directly in the waveform following each periastron passage.} These faint bursts are shown to be a superposition of the least-damped overtone (i.e., fundamental) of at least four ($l=m \le 4$) quasinormal modes. Our results suggest that QNBs should appear in gravitational waveforms, and thus provide a gauge-invariant signal. Potentially observable in high signal-to-noise ratio EMRIs, QNBs would provide high-frequency components to the parameter estimation problem that would complement low-frequency elements of the waveform., Comment: 28 pages, 11 figures, 5 tables; Updated to reflect published version

We describe a procedure for constructing initial data for boosted black holes in the moving-punctures approach to numerical relativity that endows the initial time slice from the outset with trumpet geometry within the black hole interiors. We then demonstrate the procedure in numerical simulations using an evolution code from the Einstein Toolkit that employs 1+log slicing. The Lorentz boost of a single black hole can be precisely specified and multiple, widely separated black holes can be treated approximately by superposition of single hole data. There is room within the scheme for later improvement to re-solve (iterate) the constraint equations in the multiple black hole case. The approach is shown to yield an initial trumpet slice for one black hole that is close to, and rapidly settles to, a stationary trumpet geometry. Initial data in this new approach is shown to contain initial transient (or "junk") radiation that is suppressed by as much as two orders of magnitude relative to that in comparable Bowen-York initial data., Comment: 18 pages, 18 figures

We calculate the evolution and gravitational-wave emission of a spinning compact object inspiraling into a substantially more massive (non-rotating) black hole. We extend our previous model for a non-spinning binary [Phys. Rev. D 93, 064024] to include the Mathisson-Papapetrou-Dixon spin-curvature force. For spin-aligned binaries we calculate the dephasing of the inspiral and associated waveforms relative to models that do not include spin-curvature effects. We find this dephasing can be either positive or negative depending on the initial separation of the binary. For binaries in which the spin and orbital angular momentum are not parallel, the orbital plane precesses and we use a more general osculating element prescription to compute inspirals., Comment: 17 pages, 6 figures

An exact correspondence between a black hole and an accelerating mirror is demonstrated. It is shown that for a massless minimally coupled scalar field the same Bogolubov coefficients connecting the "in" and "out" states occur for a (1+1)D flat spacetime with a particular perfectly reflecting accelerating boundary trajectory and a (1+1)D curved spacetime in which a null shell collapses to form a black hole. Generalization of the latter to the (3+1)D case is discussed. The spectral dynamics is computed in both (1+1)-dimensional spacetimes along with the energy flux in the spacetime with a mirror. It is shown that the approach to equilibrium is monotonic, asymmetric in terms of the rate, and there is a specific time which characterizes the system when it is the most out-of-equilibrium., Comment: 25 pages, 7 figures