Numerical integration for fractal measures

We find estimates for the error in replacing an integral $\int f d\mu$ with respect to a fractal measure $\mu$ with a discrete sum $\sum_{x \in E} w(x) f(x)$ over a given sample set $E$ with weights $w$. Our model is the classical Koksma-Hlawka theorem for integrals over rectangles, where the error is estimated by a product of a discrepancy that depends only on the geometry of the sample set and weights, and variance that depends only on the smoothness of $f$. We deal with p.c.f self-similar fractals, on which Kigami has constructed notions of energy and Laplacian. We develop generic results where we take the variance to be either the energy of $f$ or the $L^1$ norm of $\Delta f$, and we show how to find the corresponding discrepancies for each variance. We work out the details for a number of interesting examples of sample sets for the Sierpi\'{n}ski gasket, both for the standard self-similar measure and energy measures, and for other fractals.Comment: 50 pages, figures all contained in the main tex fil

Malmquist Bias and the Distance to the Virgo Cluster

This paper investigates the impact of Malmquist bias on the distance to the Virgo cluster determined by the H_0 Key Project using M100, and consequently on the derived value of H_0. Malmquist bias is a volume-induced statistical effect which causes the most probable distance to be different from the raw distance measured. Consideration of the bias in the distance to the Virgo cluster raises this distance and lowers the calculated value of H_0. Monte Carlo simulations of the cluster have been run for several possible distributions of spirals within the cluster and of clusters in the local universe. Simulations consistent with known information regarding the cluster and the errors of measurement result in a bias of about 6.5%-8.5%. This corresponds to an unbiased distance of 17.2-17.4 Mpc and a value of H_0 in the range 80-82 km/s/Mpc. The problem of determining the bias to Virgo illustrates several key points regarding Malmquist bias. Essentially all conventional astronomical distance measurements are subject to this bias. In addition, the bias accumulates when an attempt is made to construct "distance ladders" from measurements which are individually biased. As will be shown in the case of Virgo, the magnitude and direction of the bias are sensitive to the spatial distribution of the parent poputation from which the observed object is drawn - a distribution which is often poorly known. This leads to uncertainty in the magnitude of the bias, and adds to the importance of minimizing the number of steps in "distance ladders".Comment: 19 pages, 3 figures, Latex, To appear in Ap

A maximum-likelihood method for improving faint source flux and color estimates

Flux estimates for faint sources or transients are systematically biased high because there are far more truly faint sources than bright. Corrections which account for this effect are presented as a function of signal-to-noise ratio and the (true) slope of the faint-source number-flux relation. The corrections depend on the source being originally identified in the image in which it is being photometered. If a source has been identified in other data, the corrections are different; a prescription for calculating the corrections is presented. Implications of these corrections for analyses of surveys are discussed; the most important is that sources identified at signal-to-noise ratios of four or less are practically useless.Comment: 9 pp., accepted for publication in PAS

Self-Assessment and Planned Change of Placement and Career Services Center

In the 1990s, college and university career services and placement units face many challenges which may influence their success including (a) increased budgetary constraints, (b) changing student demographics, (c) increased availability of computer technologies, (d) new service delivery models, and (e) changing employer recruitment practices. In an effort to address these issues and examine its role within the university (Roth, 1994), the Placement and Career Information Center (PCIC) at Central Michigan University undertook a program of applied research, self-assessment and planned change. The purpose of this article is to briefly report our experiences in conducting this program of applied research. The assessment strategy is presented in the first part of this report. Next, the methods used to collect data and assessment, sampling procedures and response rates are described. Following this, highlights of the assessment results are presented including a summary of some changes already made and those planned for the future. The results of this study are presented in detail in Adams, et at. (1994)

The Correlation Between Galaxy HI Linewidths and K' Luminosities

The relationship between galaxy luminosities and rotation rates is studied with total luminosities in the K' band. Extinction problems are essentially eliminated at this band centered at 2.1 micron. A template luminosity-linewidth relation is derived based on 65 galaxies drawn from two magnitude-limited cluster samples. The zero-point is determined using 4 galaxies with accurately known distances. The calibration is applied to give the distance to the Pisces Cluster (60 Mpc) at a redshift in the CMB frame of 4771 km/s. The resultant value of the Hubble Constant is 81 km/s/Mpc. The largest sources of uncertainty arises from the small number of zero-point calibrators at this time at K' and present application to only one cluster.Comment: 13 pages including 5 figures and 2 tables. Accepted for publication in Astrophysical Journa

Biases in Virial Black Hole Masses: An SDSS Perspective

We compile black hole (BH) masses for $\sim 60,000$ quasars in the redshift range $0.1 \lesssim z \lesssim 4.5$ included in the Fifth Data Release of the Sloan Digital Sky Survey (SDSS), using virial BH mass estimators based on the \hbeta, \MgII, and \CIV emission lines. We find that: (1) within our sample, the widths of the three lines follow log-normal distributions, with means and dispersions that do not depend strongly on luminosity or redshift;(2) the \MgII- and \hbeta-estimated BH masses are consistent with one another; and (3) the \CIV BH mass estimator may be more severely affected by a disk wind component than the \MgII and \hbeta estimators, giving a positive bias in mass correlated with the \CIV-\MgII blueshift. Most SDSS quasars have virial BH masses in the range $10^8-10^9 M_\odot$. There is a clear upper mass limit of $\sim 10^{10} M_\odot$ for active BHs at $z \gtrsim 2$, decreasing at lower redshifts. Making the reasonable assumptions that the underlying BH mass distribution decreases with mass and that the Eddington ratio distribution at fixed BH mass has non-zero width, we show that the measured virial BH mass distribution and Eddington ratio distribution are subject to Malmquist bias. A radio quasar subsample (with $1.5\lesssim z\lesssim 2.3$) has mean virial BH mass larger by $\sim 0.12$ dex than the whole sample. A broad absorption line (BAL) quasar subsample (with $1.7\lesssim z\lesssim 2.2$) has identical virial mass distribution as the nonBAL sample, with no mean offset. (Abridged)Comment: Updated virial mass measurements; improved presentation of the MC simulation; added new discussion sections; conclusions unchanged. The full table1 is available at http://www.astro.princeton.edu/~yshen/BH_mass/datafile1.txt.tar.g

AcDc - A new code for the NLTE spectral analysis of accretion discs: application to the helium CV AM CVn

We present a recently developed code for detailed NLTE calculations of accretion disc spectra of cataclysmic variables and compact X-ray binaries. Assuming a radial structure of a standard alpha-disc, the disc is divided into concentric rings. For each disc ring the solution of the radiation transfer equation and the structure equations, comprising the hydrostatic and radiative equilibrium, the population of the atomic levels as well as charge and particle conservation, is done self-consistently. Metal-line blanketing and irradiation by the central object are taken into account. As a first application, we show the influence of different disc parameters on the disc spectrum for the helium cataclysmic variable AM CVn.Comment: 7 pages, 11 figures to be published in A&

Random Matrix Theory and the Sixth Painlev\'e Equation

A feature of certain ensembles of random matrices is that the corresponding measure is invariant under conjugation by unitary matrices. Study of such ensembles realised by matrices with Gaussian entries leads to statistical quantities related to the eigenspectrum, such as the distribution of the largest eigenvalue, which can be expressed as multidimensional integrals or equivalently as determinants. These distributions are well known to be $\tau$-functions for Painlev\'e systems, allowing for the former to be characterised as the solution of certain nonlinear equations. We consider the random matrix ensembles for which the nonlinear equation is the $\sigma$ form of \PVI. Known results are reviewed, as is their implication by way of series expansions for the distributions. New results are given for the boundary conditions in the neighbourhood of the fixed singularities at $t=0,1,\infty$ of $\sigma$\PVI displayed by a generalisation of the generating function for the distributions. The structure of these expansions is related to Jimbo's general expansions for the $\tau$-function of $\sigma$\PVI in the neighbourhood of its fixed singularities, and this theory is itself put in its context of the linear isomonodromy problem relating to \PVI.Comment: Dedicated to the centenary of the publication of the Painlev\'e VI equation in the Comptes Rendus de l'Academie des Sciences de Paris by Richard Fuchs in 190

A Complete Catalog of Swift GRB Spectra and Durations: Demise of a Physical Origin for Pre-Swift High-Energy Correlations

We calculate durations and spectral paramaters for 218 Swift bursts detected by the BAT instrument between and including GRBs 041220 and 070509, including 77 events with measured redshifts. Incorporating prior knowledge into the spectral fits, we are able to measure the characteristic $\nu F_{\nu}$ spectral peak energy $E_{\rm pk,obs}$ and the isotropic equivalent energy $E_{\rm iso}$ (1--$10^4$ keV) for all events. This complete and rather extensive catalog, analyzed with a unified methodology, allows us to address the persistence and origin of high-energy correlations suggested in pre-Swift observations. We find that the $E_{\rm pk,obs}$-$E_{\rm iso}$ correlation is present in the Swift sample; however, the best-fit powerlaw relation is inconsistent with the best-fit pre-Swift relation at >5 sigma significance. Moreover, it has a factor >~ 2 larger intrinsic scatter, after accounting for large errors on $E_{\rm pk,obs}$. A large fraction of the Swift events are hard and subluminous relative to (and inconsistent with) the pre-Swift relation, in agreement with indications from BATSE GRBs without redshift. Moreover, we determine an experimental threshold for the BAT detector and show how the $E_{\rm pk,obs}$--$E_{\rm iso}$ correlation arises artificially due to partial correlation with the threshold. We show that pre-Swift correlations found by Amati et al.(2002), Yonetoku et al. (2004), Firmani et al.(2006) (and independently by others) are likely unrelated to the physical properties of GRBs and are likely useless for tests of cosmology. Also, an explanation of these correlations in terms of a detector threshold provides a natural and quantitative explanation for why short-duration GRBs and events at low redshift tend to be outliers to the correlations.Comment: 25 pages, 9 figures, 2 tables, Accepted to Ap

Volume limited dependent Galactic model parameters

We estimated 34 sets of Galactic model parameters for three intermediate latitude fields with Galactic longitudes l=60, l=90, and l=180, and we discussed their dependence on the volume. Also, we confirmed the variation of these parameters with absolute magnitude and Galactic longitude. The star samples in two fields are restricted with bright and unit absolute magnitude intervals, (4,5], and (5,6], whereas for the third field a larger absolute magnitude interval is adopted, (4,10]. The limiting apparent magnitudes of star samples are g=15 and g=22.5 mag which provide space densities within distances in the line of sight 0.9 and 25 kpc. The Galactic model parameters for the thin disc are not volume dependent. However, the ones for thick disc and halo do show spectacular trends in their variations with volume, except for the scalelength of the thick disc. The local space density of the thick disc increases, whereas the scaleheight of the same Galactic component decreases monotonically. However, both model parameters approach asymptotic values at large distances. The axial ratio of the halo increases abruptly for the volumes where thick disc is dominant, whereas it approaches an asymptotic value gradually for larger volumes, indicating a continuous transition from disclike structure to a spherical one at the outermost region of the Galaxy. The variation of the Galactic model parameters with absolute magnitude can be explained by their dependence on the stellar luminosity, whereas the variation with volume and Galactic longitude at short distances is a bias in analysis.Comment: 12 pages, including 8 figures and 5 tables, accepted for publication in PAS