MCMC methods for functions modifying old algorithms to make\ud them faster

Many problems arising in applications result in the need\ud to probe a probability distribution for functions. Examples include Bayesian nonparametric statistics and conditioned diffusion processes. Standard MCMC algorithms typically become arbitrarily slow under the mesh refinement dictated by nonparametric description of the unknown function. We describe an approach to modifying a whole range of MCMC methods which ensures that their speed of convergence is robust under mesh refinement. In the applications of interest the data is often sparse and the prior specification is an essential part of the overall modeling strategy. The algorithmic approach that we describe is applicable whenever the desired probability measure has density with respect to a Gaussian process or Gaussian random field prior, and to some useful non-Gaussian priors constructed through random truncation. Applications are shown in density estimation, data assimilation in fluid mechanics, subsurface geophysics and image registration. The key design principle is to formulate the MCMC method for functions. This leads to algorithms which can be implemented via minor modification of existing algorithms, yet which show enormous speed-up on a wide range of applied problems

Collapse of Uniformly Rotating Stars to Black Holes and the Formation of Disks

Simulations in general relativity show that the outcome of collapse of a marginally unstable, uniformly rotating star spinning at the mass-shedding limit depends critically on the equation of state. For a very stiff equation of state, which is likely to characterize a neutron star, essentially all of the mass and angular momentum of the progenitor are swallowed by the Kerr black hole formed during the collapse, leaving nearly no residual gas to form a disk. For a soft equation of state with an adiabatic index \Gamma - 4/3 << 1, which characterizes a very massive or supermassive star supported predominantly by thermal radiation pressure, as much as 10% of the mass of the progenitor avoids capture and goes into a disk about the central hole. We present a semi-analytic calculation that corroborates these numerical findings and shows how the final outcome of such a collapse may be determined from simple physical considerations. In particular, we employ a simple energy variational principle with an approximate, post-Newtonian energy functional to determine the structure of a uniformly rotating, polytropic star at the onset of collapse as a function of polytropic index n, where \Gamma = 1+1/n. We then use this data to calculate the mass and spin of the final black hole and ambient disk. We show that the fraction of the total mass that remains in the disk falls off sharply as 3-n (equivalently, \Gamma - 4/3) increases.Comment: 11 pages, 2 figures, 2 tables, AASTeX; accepted to appear in The Astrophysical Journa

A first direct measurement of the intergalactic medium temperature around a quasar at z=6

The thermal state of the intergalactic medium (IGM) provides an indirect probe of both the HI and HeII reionisation epochs. Current constraints on the IGM temperature from the Lya forest are restricted to the redshift range 2<z<4.5, limiting the ability to probe the thermal memory of HI reionisation toward higher redshift. In this work, we present the first direct measurement of the IGM temperature around a z=6 quasar by analysing the Doppler widths of Lya absorption lines in the proximity zone of SDSS J0818+1722. We use a high resolution (R= 40000) Keck/HIRES spectrum in combination with detailed numerical modelling to obtain the temperature at mean density, T_0=23600\pm^5000_6900K (\pm^9200_9300K) at 68 (95) per cent confidence assuming a prior probability 13500K<T_0<38500 K following HI and HeII reionisation. This enables us to place an upper limit on the redshift of HI reionisation, z_H, within 33 comoving Mpc of SDSS J0818+1722. If the quasar reionises the HeII in its vicinity, then in the limit of instantaneous reionisation we infer z_H<9.0 (11.0) at 68 (95) per cent confidence assuming photoheating is the dominant heat source and that HI reionisation is driven by ionising sources with soft spectra, typical of population II stars. If the HI and HeII in the IGM around SDSS J0818+1722 are instead reionised simultaneously by a population of massive metal-free stars, characterised by very hard ionising spectra, we obtain a tighter upper limit of z_H<8.4 (9.4). Initiating reionisation at higher redshifts produces temperatures which are too low with respect to our constraint unless the HI ionising sources or the quasar itself have spectra significantly harder than typically assumed.Comment: 15 pages, 9 figures, accepted to MNRA

Black Hole-Neutron Star Binaries in General Relativity: Quasiequilibrium Formulation

We present a new numerical method for the construction of quasiequilibrium models of black hole-neutron star binaries. We solve the constraint equations of general relativity, decomposed in the conformal thin-sandwich formalism, together with the Euler equation for the neutron star matter. We take the system to be stationary in a corotating frame and thereby assume the presence of a helical Killing vector. We solve these coupled equations in the background metric of a Kerr-Schild black hole, which accounts for the neutron star's black hole companion. In this paper we adopt a polytropic equation of state for the neutron star matter and assume large black hole--to--neutron star mass ratios. These simplifications allow us to focus on the construction of quasiequilibrium neutron star models in the presence of strong-field, black hole companions. We summarize the results of several code tests, compare with Newtonian models, and locate the onset of tidal disruption in a fully relativistic framework.Comment: 17 pages, 7 figures; added discussion, tables; PRD in pres

Use of Social Media by Sales Students for First Post-Graduation Job-Search: An Exploratory Study

As students who studied/majored in professional sales approach graduation and start looking for their first post-graduation/entry-level positions, what avenues are used to pursue these opportunities? One might suspect that with the explosion of social media and the millennial generation’s prolific use of said media this would be their preferred choice for job hunting. Are sales students engaging social media and other online activities to land entry-level career positions or are they still using traditional methods (i.e. College/University Career Fairs, College/University Career Services, Internships, etc.) to find these positons? This paper looks at the avenues millennial students who studied/majored in professional sales use to land their first post-college graduation sales positions

Merger of white dwarf-neutron star binaries: Prelude to hydrodynamic simulations in general relativity

White dwarf-neutron star binaries generate detectable gravitational radiation. We construct Newtonian equilibrium models of corotational white dwarf-neutron star (WDNS) binaries in circular orbit and find that these models terminate at the Roche limit. At this point the binary will undergo either stable mass transfer (SMT) and evolve on a secular time scale, or unstable mass transfer (UMT), which results in the tidal disruption of the WD. The path a given binary will follow depends primarily on its mass ratio. We analyze the fate of known WDNS binaries and use population synthesis results to estimate the number of LISA-resolved galactic binaries that will undergo either SMT or UMT. We model the quasistationary SMT epoch by solving a set of simple ordinary differential equations and compute the corresponding gravitational waveforms. Finally, we discuss in general terms the possible fate of binaries that undergo UMT and construct approximate Newtonian equilibrium configurations of merged WDNS remnants. We use these configurations to assess plausible outcomes of our future, fully relativistic simulations of these systems. If sufficient WD debris lands on the NS, the remnant may collapse, whereby the gravitational waves from the inspiral, merger, and collapse phases will sweep from LISA through LIGO frequency bands. If the debris forms a disk about the NS, it may fragment and form planets.Comment: 28 pages, 25 figures, 6 table

Precision Predictions for (Un)Stable W+W- Pair Production At and Beyond LEP2 Energies Beyond LEP2 Energies

We present precision calculations of the processes e+e- -> 4-fermions in which the double resonant W+W- intermediate state occurs. Referring to this latter intermediate state as the 'signal process', we show that, by using the YFS Monte Carlo event generators YFSWW3-1.14 and KORALW1.42 in an appropriate combination, we achieve a physical precision on the signal process, as isolated with LEP2 MC Workshop cuts, below 0.5 per cent. We stress the full gauge invariance of our calculations and we compare our results with those of other authors where appropriate. In particular, sample Monte Carlo data are explicitly illustrated and compared with the results of the program RacoonWW of Dittmaier {\it et al.}. In this way, we show that the total (physical plus technical) precision tag for the WW signal process cross section is 0.4 per cent for 200 GeV, for example. Results are also given for 500 GeV with an eye toward the LC.Comment: 19 pages, 2 figs;corrected Tab. 3;improved refs.,figs.,text;improved refs.,text;improved tex

Quasiequilibrium black hole-neutron star binaries in general relativity

We construct quasiequilibrium sequences of black hole-neutron star binaries in general relativity. We solve Einstein's constraint equations in the conformal thin-sandwich formalism, subject to black hole boundary conditions imposed on the surface of an excised sphere, together with the relativistic equations of hydrostatic equilibrium. In contrast to our previous calculations we adopt a flat spatial background geometry and do not assume extreme mass ratios. We adopt a Gamma=2 polytropic equation of state and focus on irrotational neutron star configurations as well as approximately nonspinning black holes. We present numerical results for ratios of the black hole's irreducible mass to the neutron star's ADM mass in isolation of M_{irr}^{BH}/M_{ADM,0}^{NS} = 1, 2, 3, 5, and 10. We consider neutron stars of baryon rest mass M_B^{NS}/M_B^{max} = 83% and 56%, where M_B^{max} is the maximum allowed rest mass of a spherical star in isolation for our equation of state. For these sequences, we locate the onset of tidal disruption and, in cases with sufficiently large mass ratios and neutron star compactions, the innermost stable circular orbit. We compare with previous results for black hole-neutron star binaries and find excellent agreement with third-order post-Newtonian results, especially for large binary separations. We also use our results to estimate the energy spectrum of the outgoing gravitational radiation emitted during the inspiral phase for these binaries.Comment: 17 pages, 15 figures, published in Phys. Rev.

Gravitational Collapse with a Cosmological Constant

We consider the effect of a positive cosmological constant on spherical gravitational collapse to a black hole for a few simple, analytic cases. We construct the complete Oppenheimer-Snyder-deSitter (OSdS) spacetime, the generalization of the Oppenheimer-Snyder solution for collapse from rest of a homogeneous dust ball in an exterior vacuum. In OSdS collapse, the cosmological constant may affect the onset of collapse and decelerate the implosion initially, but it plays a diminishing role as the collapse proceeds. We also construct spacetimes in which a collapsing dust ball can bounce, or hover in unstable equilibrium, due to the repulsive force of the cosmological constant. We explore the causal structure of the different spacetimes and identify any cosmological and black hole event horizons which may be present.Comment: 7 pages, 10 figures; To appear in Phys. Rev.