Testing Binary Population Synthesis Models with Hot Subdwarfs

Models of binary star interactions have been successful in explaining the origin of field hot subdwarf (sdB) stars in short period systems, but longer-period systems that formed via Roche-lobe overflow (RLOF) mass transfer from the present sdB to its companion have received less attention. We map sets of initial binaries into present-day binaries that include sdBs and main-sequence stars, distinguishing "observable" sdBs from "hidden" ones. We aim to find out whether (1) the existing catalogues of sdBs are sufficiently fair samples of all the kinds of sdB binaries that theory predicts; or instead whether (2) large predicted hidden populations mandate the construction of new catalogues, perhaps using wide-field imaging surveys such as 2MASS, SDSS, and Galex. We also report on a pilot study to identify hidden subdwarfs, using 2MASS and GALEX data.Comment: 3 pages with 2 figures. Uses AIP style files. To appear in Future Directions in Ultraviolet Astronomy, ed. Michael E. VanSteenberg (AIP Conf Proc

The Ultrasensitivity of Living Polymers

Synthetic and biological living polymers are self-assembling chains whose chain length distributions (CLDs) are dynamic. We show these dynamics are ultrasensitive: even a small perturbation (e.g. temperature jump) non-linearly distorts the CLD, eliminating or massively augmenting short chains. The origin is fast relaxation of mass variables (mean chain length, monomer concentration) which perturbs CLD shape variables before these can relax via slow chain growth rate fluctuations. Viscosity relaxation predictions agree with experiments on the best-studied synthetic system, alpha-methylstyrene.Comment: 4 pages, submitted to Phys. Rev. Let

Nonlinear equation for anomalous diffusion: unified power-law and stretched exponential exact solution

The nonlinear diffusion equation $\frac{\partial \rho}{\partial t}=D \tilde{\Delta} \rho^\nu$ is analyzed here, where $\tilde{\Delta}\equiv \frac{1}{r^{d-1}}\frac{\partial}{\partial r} r^{d-1-\theta} \frac{\partial}{\partial r}$, and $d$, $\theta$ and $\nu$ are real parameters. This equation unifies the anomalous diffusion equation on fractals ($\nu =1$) and the spherical anomalous diffusion for porous media ($\theta=0$). Exact point-source solution is obtained, enabling us to describe a large class of subdiffusion ($\theta > (1-\nu)d$), normal diffusion ($\theta= (1-\nu)d$) and superdiffusion ($\theta < (1-\nu)d$). Furthermore, a thermostatistical basis for this solution is given from the maximum entropic principle applied to the Tsallis entropy.Comment: 3 pages, 2 eps figure

Non-Equilibrium in Adsorbed Polymer Layers

High molecular weight polymer solutions have a powerful tendency to deposit adsorbed layers when exposed to even mildly attractive surfaces. The equilibrium properties of these dense interfacial layers have been extensively studied theoretically. A large body of experimental evidence, however, indicates that non-equilibrium effects are dominant whenever monomer-surface sticking energies are somewhat larger than kT, a common case. Polymer relaxation kinetics within the layer are then severely retarded, leading to non-equilibrium layers whose structure and dynamics depend on adsorption kinetics and layer ageing. Here we review experimental and theoretical work exploring these non-equilibrium effects, with emphasis on recent developments. The discussion addresses the structure and dynamics in non-equilibrium polymer layers adsorbed from dilute polymer solutions and from polymer melts and more concentrated solutions. Two distinct classes of behaviour arise, depending on whether physisorption or chemisorption is involved. A given adsorbed chain belonging to the layer has a certain fraction of its monomers bound to the surface, f, and the remainder belonging to loops making bulk excursions. A natural classification scheme for layers adsorbed from solution is the distribution of single chain f values, P(f), which may hold the key to quantifying the degree of irreversibility in adsorbed polymer layers. Here we calculate P(f) for equilibrium layers; we find its form is very different to the theoretical P(f) for non-equilibrium layers which are predicted to have infinitely many statistical classes of chain. Experimental measurements of P(f) are compared to these theoretical predictions.Comment: 29 pages, Submitted to J. Phys.: Condens. Matte

Geometrical optics analysis of the short-time stability properties of the Einstein evolution equations

Many alternative formulations of Einstein's evolution have lately been examined, in an effort to discover one which yields slow growth of constraint-violating errors. In this paper, rather than directly search for well-behaved formulations, we instead develop analytic tools to discover which formulations are particularly ill-behaved. Specifically, we examine the growth of approximate (geometric-optics) solutions, studied only in the future domain of dependence of the initial data slice (e.g. we study transients). By evaluating the amplification of transients a given formulation will produce, we may therefore eliminate from consideration the most pathological formulations (e.g. those with numerically-unacceptable amplification). This technique has the potential to provide surprisingly tight constraints on the set of formulations one can safely apply. To illustrate the application of these techniques to practical examples, we apply our technique to the 2-parameter family of evolution equations proposed by Kidder, Scheel, and Teukolsky, focusing in particular on flat space (in Rindler coordinates) and Schwarzchild (in Painleve-Gullstrand coordinates).Comment: Submitted to Phys. Rev.

Mapping population synthesis event rates on model parameters II: Convergence and accuracy of multidimensional fits

Binary population synthesis calculations and associated predictions, especially event rates, are known to depend on a significant number of input model parameters with different degrees of sensitivity. At the same time, for systems with relatively low formation rates, such simulations are heavily computationally demanding and therefore the needed explorations of the high-dimensional parameter space require major -- often prohibitive -- computational resources. In the present study, to better understand several key event rates involving binary evolution and binaries with two compact objects in Milky Way-like galaxies and to provide ways of reducing the computational costs of complete parameter space explorations: (i) we perform a methodical parameter study of the \emph{StarTrack} population synthesis code ; and (ii) we develop a formalism and methodology for the derivation of {\em multi-dimensional fits} for event rates. We significantly generalize our earlier study, and we focus on ways of thoroughly assessing the accuracy of the fits. We anticipate that the efficient tools developed here can be applied in lieu of large-scale population calculations and will facilitate the exploration of the dependence of rate predictions on a wide range binary evolution parameters. Such explorations can then allow the derivation of constraints on these parameters, given empirical rate constraints and accounting for fitting errors. Here we describe in detail the principles and practice behind constructing these fits, estimating their accuracy, and comparing them with observations in a manner that accounts for their errors.Comment: 19 pages, 9 figures, formatted with emulateapj. Submitted to ApJ. v2 has been completely rewritten to improve clarit

Nonlinear anomalous diffusion equation and fractal dimension: Exact generalized gaussian solution

In this work we incorporate, in a unified way, two anomalous behaviors, the power law and stretched exponential ones, by considering the radial dependence of the $N$-dimensional nonlinear diffusion equation $\partial\rho /\partial{t}={\bf \nabla} \cdot (K{\bf \nabla} \rho^{\nu})-{\bf \nabla}\cdot(\mu{\bf F} \rho)-\alpha \rho ,$ where $K=D r^{-\theta}$, $\nu$, $\theta$, $\mu$ and $D$ are real parameters and $\alpha$ is a time-dependent source. This equation unifies the O'Shaugnessy-Procaccia anomalous diffusion equation on fractals ($\nu =1$) and the spherical anomalous diffusion for porous media ($\theta=0$). An exact spherical symmetric solution of this nonlinear Fokker-Planck equation is obtained, leading to a large class of anomalous behaviors. Stationary solutions for this Fokker-Planck-like equation are also discussed by introducing an effective potential.Comment: Latex, 6 pages. To appear in Phys. Rev.

Random walk through fractal environments

We analyze random walk through fractal environments, embedded in 3-dimensional, permeable space. Particles travel freely and are scattered off into random directions when they hit the fractal. The statistical distribution of the flight increments (i.e. of the displacements between two consecutive hittings) is analytically derived from a common, practical definition of fractal dimension, and it turns out to approximate quite well a power-law in the case where the dimension D of the fractal is less than 2, there is though always a finite rate of unaffected escape. Random walks through fractal sets with D less or equal 2 can thus be considered as defective Levy walks. The distribution of jump increments for D > 2 is decaying exponentially. The diffusive behavior of the random walk is analyzed in the frame of continuous time random walk, which we generalize to include the case of defective distributions of walk-increments. It is shown that the particles undergo anomalous, enhanced diffusion for D_F < 2, the diffusion is dominated by the finite escape rate. Diffusion for D_F > 2 is normal for large times, enhanced though for small and intermediate times. In particular, it follows that fractals generated by a particular class of self-organized criticality (SOC) models give rise to enhanced diffusion. The analytical results are illustrated by Monte-Carlo simulations.Comment: 22 pages, 16 figures; in press at Phys. Rev. E, 200

Impact of subdominant modes on the interpretation of gravitational-wave signals from heavy binary black hole systems

Over the past year, a handful of new gravitational wave models have been developed to include multiple harmonic modes thereby enabling for the first time fully Bayesian inference studies including higher modes to be performed. Using one recently developed numerical relativity surrogate model, NRHybSur3dq8, we investigate the importance of higher modes on parameter inference of coalescing massive binary black holes. We focus on examples relevant to the current three-detector network of observatories, with a detector-frame mass set to 120 M⊙ and with signal amplitude values that are consistent with plausible candidates for the next few observing runs. We show that for such systems the higher mode content will be important for interpreting coalescing binary black holes, reducing systematic bias, and computing properties of the remnant object. Even for comparable-mass binaries and at low signal amplitude, the omission of higher modes can influence posterior probability distributions. We discuss the impact of our results on source population inference and self-consistency tests of general relativity. Our work can be used to better understand asymmetric binary black hole merger events, such as GW190412. Higher modes are critical for such systems, and their omission usually produces substantial parameter biases