33,975 research outputs found

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

    Get PDF
    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

    A Theoretical Model for the Mbh−σM_{\rm bh}-\sigma Relation for Supermassive Black Holes in Galaxies

    Full text link
    We construct a model for the formation of black holes within galactic bulges. The initial state is a slowly rotating isothermal sphere, characterized by effective transport speed \aeff and rotation rate Ω\Omega. The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole. This model reproduces the observed correlation between black hole masses and galactic velocity dispersions, \mbh \approx 10^8 M_\odot (\sigma/200 \kms)^4, where \sigma = \sqrt{2} \aeff. This model also predicts the ratio \mrat of black hole mass to host mass: \mrat ≈\approx 0.004 (\sigma/200 \kms).Comment: 9 pages, 2 figures, submitted to Astrophysical Journal Letter

    Geothermal studies - Yellowstone National Park /test site 11/, Wyoming

    Get PDF
    Summary report of diamond drilling in thermal areas of Yellowstone National Park, and method for determining heat flow in thermal area

    X-ray and Radio Interactions in the Cores of Cooling Flow Clusters

    Get PDF
    We present high resolution ROSAT x-ray and radio observations of three cooling flow clusters containing steep spectrum radio sources at their cores. All three systems exhibit strong signs of interaction between the radio plasma and the hot intracluster medium. Two clusters, A133 and A2626, show enhanced x-ray emission spatially coincident with the radio source whereas the third cluster, A2052, exhibits a large region of x-ray excess surrounding much of the radio source. Using 3-D numerical simulations, we show that a perturbed jet propagating through a cooling flow atmosphere can give rise to amorphous radio morphologies, particularly in the case where the jet was ``turned off'' and allowed to age passively. In addition, the simulated x-ray surface brightness produced both excesses and deficits as seen observationally.Comment: 25 pages, 10 figures, accepted for publication in A

    Theoretical uncertainty in baryon oscillations

    Full text link
    We discuss the systematic uncertainties in the recovery of dark energy properties from the use of baryon acoustic oscillations as a standard ruler. We demonstrate that while unknown relativistic components in the universe prior to recombination would alter the sound speed, the inferences for dark energy from low-redshift surveys are unchanged so long as the microwave background anisotropies can measure the redshift of matter-radiation equality, which they can do to sufficient accuracy. The mismeasurement of the radiation and matter densities themselves (as opposed to their ratio) would manifest as an incorrect prediction for the Hubble constant at low redshift. In addition, these anomalies do produce subtle but detectable features in the microwave anisotropies.Comment: 4 pages, REVTeX, 1 figure. Submitted to PR

    A Prediction of Observable Rotation in the ICM of Abell 3266

    Full text link
    We present a numerical Hydro+N-body model of A3266 whose X-ray surface brightness, temperature distribution, and galaxy spatial and velocity distribution data are consistent with the A3266 data. The model is an old (~3 Gyr), off-axis merger having a mass ratio of ~2.5:1. The less massive subcluster in the model is moving on a trajectory from southwest to northeast passing on the western side of the dominant cluster while moving into the plane of the sky at ~45 degrees. Off-axis mergers such as this one are an effective mechanism for transferring angular momentum to the intracluster medium (ICM), making possible a large scale rotation of the ICM. We demonstrate here that the ICM rotation predicted by our fully 3-dimensional model of A3266 is observable with current technology. As an example, we present simulated observations assuming the capabilities of the high resolution X-ray spectrometer (XRS) which was to have flown on Astro-E.Comment: 9 pages, 7 postscript figures, Fig. 3 and 6 are color postscript, Accepted for publication in the Astrophysical Journa

    Measured performance of the new University of California gamma ray telescope

    Get PDF
    The design of the new medium energy balloon-borne gamma ray telescope is discussed. This telescope is sensitive to 1-30 MeV gamma rays. The results of the initial calibration are described. The position and energy resolutions of 32 plastic and NaI(Tl) scintillator bars, each 100 cm long are discussed. The telescope's measured angular and energy resolutions as a function of incident angle are compared with detailed Monte Carlo calculations at 1.37, 2.75 and 6.13 MeV. The expected resolutions are 5 deg FHWM and 8% at 2.75 MeV. The expected area-efficiency is 250 cm

    Recent progress in the truncated Lanczos method : application to hole-doped spin ladders

    Full text link
    The truncated Lanczos method using a variational scheme based on Hilbert space reduction as well as a local basis change is re-examined. The energy is extrapolated as a power law function of the Hamiltonian variance. This systematic extrapolation procedure is tested quantitatively on the two-leg t-J ladder with two holes. For this purpose, we have carried out calculations of the spin gap and of the pair dispersion up to size 2x15.Comment: 5 pages, 4 included eps figures, submitted to Phys. Rev. B; revised versio
    • …
    corecore