Diffusion as mixing mechanism in granular materials

We present several numerical results on granular mixtures. In particular, we examine the efficiency of diffusion as a mixing mechanism in these systems. The collisions are inelastic and to compensate the energy loss, we thermalize the grains by adding a random force. Starting with a segregated system, we show that uniform agitation (heating) leads to a uniform mixture of grains of different sizes. We define a characteristic mixing time, $\tau_{mix}$, and study theoretically and numerically its dependence on other parameters like the density. We examine a model for bidisperse systems for which we can calculate some physical quantities. We also examine the effect of a temperature gradient and demonstrate the appearance of an expected segregation.Comment: 15 eps figures, include

Phenomenological glass model for vibratory granular compaction

A model for weakly excited granular media is derived by combining the free volume argument of Nowak et al. [Phys. Rev. E 57, 1971 (1998)] and the phenomenological model for supercooled liquids of Adam and Gibbs [J. Chem. Phys. 43, 139 (1965)]. This is made possible by relating the granular excitation parameter \Gamma, defined as the peak acceleration of the driving pulse scaled by gravity, to a temperature-like parameter \eta(\Gamma). The resulting master equation is formally identical to that of Bouchaud's trap model for glasses [J. Phys. I 2, 1705 (1992)]. Analytic and simulation results are shown to compare favourably with a range of known experimental behaviour. This includes the logarithmic densification and power spectrum of fluctuations under constant \eta, the annealing curve when \eta is varied cyclically in time, and memory effects observed for a discontinuous shift in \eta. Finally, we discuss the physical interpretation of the model parameters and suggest further experiments for this class of systems.Comment: 2 references added; some figure labels tweaked. To appear in PR

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO's second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=3.47×10-25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering. © 2019 American Physical Society

Search for Tensor, Vector, and Scalar Polarizations in the Stochastic Gravitational-Wave Background

The detection of gravitational waves with Advanced LIGO and Advanced Virgo has enabled novel tests of general relativity, including direct study of the polarization of gravitational waves. While general relativity allows for only two tensor gravitational-wave polarizations, general metric theories can additionally predict two vector and two scalar polarizations. The polarization of gravitational waves is encoded in the spectral shape of the stochastic gravitational-wave background, formed by the superposition of cosmological and individually unresolved astrophysical sources. Using data recorded by Advanced LIGO during its first observing run, we search for a stochastic background of generically polarized gravitational waves. We find no evidence for a background of any polarization, and place the first direct bounds on the contributions of vector and scalar polarizations to the stochastic background. Under log-uniform priors for the energy in each polarization, we limit the energy densities of tensor, vector, and scalar modes at 95% credibility to Ω0T<5.58×10-8, Ω0V<6.35×10-8, and Ω0S<1.08×10-7 at a reference frequency f0=25 Hz. © 2018 American Physical Society

Erratum: "A Gravitational-wave Measurement of the Hubble Constant Following the Second Observing Run of Advanced LIGO and Virgo" (2021, ApJ, 909, 218)

[no abstract available

Basic Methods for Computing Special Functions

This paper gives an overview of methods for the numerical evaluation of special functions, that is, the functions that arise in many problems from mathematical physics, engineering, probability theory, and other applied sciences. We consider in detail a selection of basic methods which are frequently used in the numerical evaluation of special functions: converging and asymptotic series, including Chebyshev expansions, linear recurrence relations, and numerical quadrature. Several other methods are available and some of these will be discussed in less detail. We give examples of recent software for special functions where these methods are used. We mention a list of new publications on computational aspects of special functions available on our website

On the progenitor of binary neutron star merger GW170817

On 2017 August 17 the merger of two compact objects with masses consistent with two neutron stars was discovered through gravitational-wave (GW170817), gamma-ray (GRB 170817A), and optical (SSS17a/AT 2017gfo) observations. The optical source was associated with the early-type galaxy NGC 4993 at a distance of just ∼40 Mpc, consistent with the gravitational-wave measurement, and the merger was localized to be at a projected distance of ∼2 kpc away from the galaxy's center. We use this minimal set of facts and the mass posteriors of the two neutron stars to derive the first constraints on the progenitor of GW170817 at the time of the second supernova (SN). We generate simulated progenitor populations and follow the three-dimensional kinematic evolution from binary neutron star (BNS) birth to the merger time, accounting for pre-SN galactic motion, for considerably different input distributions of the progenitor mass, pre-SN semimajor axis, and SN-kick velocity. Though not considerably tight, we find these constraints to be comparable to those for Galactic BNS progenitors. The derived constraints are very strongly influenced by the requirement of keeping the binary bound after the second SN and having the merger occur relatively close to the center of the galaxy. These constraints are insensitive to the galaxy's star formation history, provided the stellar populations are older than 1 Gyr