Highly eccentric inspirals into a black hole

We model the inspiral of a compact stellar-mass object into a massive nonrotating black hole including all dissipative and conservative first-order-in-the-mass-ratio effects on the orbital motion. The techniques we develop allow inspirals with initial eccentricities as high as $e\sim0.8$ and initial separations as large as $p\sim 50$ to be evolved through many thousands of orbits up to the onset of the plunge into the black hole. The inspiral is computed using an osculating elements scheme driven by a hybridized self-force model, which combines Lorenz-gauge self-force results with highly accurate flux data from a Regge-Wheeler-Zerilli code. The high accuracy of our hybrid self-force model allows the orbital phase of the inspirals to be tracked to within $\sim0.1$ radians or better. The difference between self-force models and inspirals computed in the radiative approximation is quantified.Comment: Updated to reflect published versio

Fast spectral source integration in black hole perturbation calculations

This paper presents a new technique for achieving spectral accuracy and fast computational performance in a class of black hole perturbation and gravitational self-force calculations involving extreme mass ratios and generic orbits. Called \emph{spectral source integration} (SSI), this method should see widespread future use in problems that entail (i) point-particle description of the small compact object, (ii) frequency domain decomposition, and (iii) use of the background eccentric geodesic motion. Frequency domain approaches are widely used in both perturbation theory flux-balance calculations and in local gravitational self-force calculations. Recent self-force calculations in Lorenz gauge, using the frequency domain and method of extended homogeneous solutions, have been able to accurately reach eccentricities as high as $e \simeq 0.7$. We show here SSI successfully applied to Lorenz gauge. In a double precision Lorenz gauge code, SSI enhances the accuracy of results and makes a factor of three improvement in the overall speed. The primary initial application of SSI--for us its \emph{raison d'\^{e}tre}--is in an arbitrary precision \emph{Mathematica} code that computes perturbations of eccentric orbits in the Regge-Wheeler gauge to extraordinarily high accuracy (e.g., 200 decimal places). These high accuracy eccentric orbit calculations would not be possible without the exponential convergence of SSI. We believe the method will extend to work for inspirals on Kerr, and will be the subject of a later publication. SSI borrows concepts from discrete-time signal processing and is used to calculate the mode normalization coefficients in perturbation theory via sums over modest numbers of points around an orbit. A variant of the idea is used to obtain spectral accuracy in solution of the geodesic orbital motion.Comment: 15 pages, 7 figure

Evolution of small-mass-ratio binaries with a spinning secondary

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 figure

Spontaneous Breaking of Translational Invariance in One-Dimensional Stationary States on a Ring

We consider a model in which positive and negative particles diffuse in an asymmetric, CP-invariant way on a ring. The positive particles hop clockwise, the negative counterclockwise and oppositely-charged adjacent particles may swap positions. Monte-Carlo simulations and analytic calculations suggest that the model has three phases; a "pure" phase in which one has three pinned blocks of only positive, negative particles and vacancies, and in which translational invariance is spontaneously broken, a "mixed" phase with a non-vanishing current in which the three blocks are positive, negative and neutral, and a disordered phase without blocks.Comment: 7 pages, LaTeX, needs epsf.st

Direct conversion of rheological compliance measurements into storage and loss moduli

We remove the need for Laplace/inverse-Laplace transformations of experimental data, by presenting a direct and straightforward mathematical procedure for obtaining frequency-dependent storage and loss moduli ($G'(\omega)$ and $G"(\omega)$ respectively), from time-dependent experimental measurements. The procedure is applicable to ordinary rheological creep (stress-step) measurements, as well as all microrheological techniques, whether they access a Brownian mean-square displacement, or a forced compliance. Data can be substituted directly into our simple formula, thus eliminating traditional fitting and smoothing procedures that disguise relevant experimental noise.Comment: 4 page

Managing Risk After Intracerebral Hemorrhage in Concomitant Atrial Fibrillation and Cerebral Amyloid Angiopathy.

The Dunhill Medical Trust (Grant ID: RTF44/0114)This is the author accepted manuscript. The final version is available from the American Heart Association via http://dx.doi.org/10.1161/STROKEAHA.116.01332

Steady-state conduction in self-similar billiards

The self-similar Lorentz billiard channel is a spatially extended deterministic dynamical system which consists of an infinite one-dimensional sequence of cells whose sizes increase monotonically according to their indices. This special geometry induces a nonequilibrium stationary state with particles flowing steadily from the small to the large scales. The corresponding invariant measure has fractal properties reflected by the phase-space contraction rate of the dynamics restricted to a single cell with appropriate boundary conditions. In the near-equilibrium limit, we find numerical agreement between this quantity and the entropy production rate as specified by thermodynamics

The Chemical Compositions of the Type II Cepheids -- The BL Her and W Vir Variables

Abundance analyses from high-resolution optical spectra are presented for 19 Type II Cepheids in the Galactic field. The sample includes both short-period (BL Her) and long-period (W Vir) stars. This is the first extensive abundance analysis of these variables. The C, N, and O abundances with similar spreads for the BL Her and W Vir show evidence for an atmosphere contaminated with $3\alpha$-process and CN-cycling products. A notable anomaly of the BL Her stars is an overabundance of Na by a factor of about five relative to their presumed initial abundances. This overabundance is not seen in the W Vir stars. The abundance anomalies running from mild to extreme in W Vir stars but not seen in the BL Her stars are attributed to dust-gas separation that provides an atmosphere deficient in elements of high condensation temperature, notably Al, Ca, Sc, Ti, and $s$-process elements. Such anomalies have previously been seen among RV Tau stars which represent a long-period extension of the variability enjoyed by the Type II Cepheids. Comments are offered on how the contrasting abundance anomalies of BL Her and W Vir stars may be explained in terms of the stars' evolution from the blue horizontal branch.Comment: 41 pages including 11 figures and 4 tables; Accepted for publication in Ap

Superconductivity from perturbative one-gluon exchange in high density quark matter

We study color superconductivity in QCD at asymptotically large chemical potential. In this limit, pairing is dominated by perturbative one-gluon exchange. We derive the Eliashberg equation for the pairing gap and solve this equation numerically. Taking into account both magnetic and electric gluon exchanges, we find $\Delta\sim g^{-5}\exp(-c/g)$ with $c=3\pi^2/\sqrt{2}$, verifying a recent result by Son. For chemical potentials that are of physical interest, $\mu< 1$ GeV, the calculation ceases to be reliable quantitatively, but our results suggest that the gap can be as large as 100 MeV.Comment: 19 pages, 6 figures. I accidentally replaced the paper with an outdated version. This version has typos corrected and will appear in PR