A NEW METHOD OF STUDYING GALAXY VELOCITY FIELD: FIRST RESULTS

It is shown that the average magnitude - redshift diagram for a magnitude-limited sample of galaxies can be effectively used for mapping galaxy velocity field. This new approach is quantitatively introduced and applied to the CfA, ESO/LV and ZCAT samples of galaxies. Strong arguments are presented for infall velocities in the Virgo and Fornax superclusters and in the Great Wall structure.Comment: MNRAS submitted, 11 pages, uuencoded compressed Postscript file, including 7 figures. Postscript also available at anonymous ftp://sirius.astrouw.edu.pl/pub/submit/jch/p1.ps.

Extinction in the Galaxy from surface brightnesses of ESO-LV galaxies: testing 'standard' extinction maps

The relative extinction in the Galaxy computed with our new method (Choloniewski and Valentijn 2003, CV) is compared with three patterns: Schlegel, Finkbeiner and Davis (1998, SFD), Burstein and Heiles (1978, BH) and the cosecans law. It is shown that extinction of SFD is more reliable then that of BH since it stronger correlates with our new extinction. The smallest correlation coeffcient have been obtained for the cosecans law. Linear regression analysis show that SFD overestimate the extinction by a factor of 1.4. Our results clearly indicate that there is non-zero extinction at the Galactic South pole and that the extinction near the Galactic equator ($|b|<40^o$) is significantly larger in the Southern hemisphere than in the Northern.Comment: 15 pages, 10 figures, submitted to Acta Astronomic

Disentangling Morphology, Star Formation, Stellar Mass, and Environment in Galaxy Evolution

We present a study of the spectroscopic and photometric properties of galaxies in six nearby clusters. We perform a partial correlation analysis on our dataset to investigate whether the correlation between star formation rates in galaxies and their environment is merely another aspect of correlations of morphology, stellar mass, or mean stellar age with environment, or whether star formation rates vary independently of these other correlations. We find a residual correlation of ongoing star formation with environment, indicating that even galaxies with similar morphologies, stellar masses, and mean stellar ages have lower star formation rates in denser environments. Thus, the current star formation gradient in clusters is not just another aspect of the morphology-density, stellar mass-density, or mean stellar age-density relations. Furthermore, the star formation gradient cannot be solely the result of initial conditions, but must partly be due to subsequent evolution through a mechanism (or mechanisms) sensitive to environment. Our results constitute a true ``smoking gun'' pointing to the effect of environment on the later evolution of galaxies.Comment: 31 pages, including 5 figures; accepted for publication in Ap

The U-band Galaxy Luminosity Function of Nearby Clusters

Despite the great potential of the U-band galaxy luminosity function (GLF) to constrain the history of star formation in clusters, to clarify the question of variations of the GLF across filter bands, to provide a baseline for comparisons to high-redshift studies of the cluster GLF, and to estimate the contribution of bound systems of galaxies to the extragalactic near-UV background, determinations have so far been hampered by the generally low efficiency of detectors in the U-band and by the difficulty of constructing both deep and wide surveys. In this paper, we present U-band GLFs of three nearby, rich clusters to a limit of M_U=-17.5 (M*_U+2). Our analysis is based on a combination of separate spectroscopic and R-band and U-band photometric surveys. For this purpose, we have developed a new maximum-likelihood algorithm for calculating the luminosity function that is particularly useful for reconstructing the galaxy distribution function in multi-dimensional spaces (e.g., the number of galaxies as a simultaneous function of luminosity in different filter bands, surface brightness, star formation rate, morphology, etc.), because it requires no prior assumptions as to the shape of the distribution function. The composite luminosity function can be described by a Schechter function with characteristic magnitude M*_U=-19.82+/-0.27 and faint end slope alpha_U=-1.09+/-0.18. The total U-band GLF is slightly steeper than the R-band GLF, indicating that cluster galaxies are bluer at fainter magnitudes. Quiescent galaxies dominate the cumulative U-band flux for M_U<-14. The contribution of galaxies in nearby clusters to the U-band extragalactic background is <1% Gyr^-1 for clusters of masses ~3*10^14 to 2*10^15 M_solar.Comment: 44 pages, 11 figures, accepted for publication in Ap

A Statistical Method for Estimating Luminosity Functions using Truncated Data

The observational limitations of astronomical surveys lead to significant statistical inference challenges. One such challenge is the estimation of luminosity functions given redshift $z$ and absolute magnitude $M$ measurements from an irregularly truncated sample of objects. This is a bivariate density estimation problem; we develop here a statistically rigorous method which (1) does not assume a strict parametric form for the bivariate density; (2) does not assume independence between redshift and absolute magnitude (and hence allows evolution of the luminosity function with redshift); (3) does not require dividing the data into arbitrary bins; and (4) naturally incorporates a varying selection function. We accomplish this by decomposing the bivariate density into nonparametric and parametric portions. There is a simple way of estimating the integrated mean squared error of the estimator; smoothing parameters are selected to minimize this quantity. Results are presented from the analysis of a sample of quasars.Comment: 30 pages, 9 figures, Accepted for publication in Ap

Tests of Statistical Methods for Estimating Galaxy Luminosity Function and Applications to the Hubble Deep Field

We studied the statistical methods for the estimation of the luminosity function (LF) of galaxies. We focused on four nonparametric estimators: $1/V_{\rm max}$ estimator, maximum-likelihood estimator of Efstathiou et al. (1988), Cho{\l}oniewski's estimator, and improved Lynden-Bell's estimator. The performance of the $1/V_{\rm max}$ estimator has been recently questioned, especially for the faint-end estimation of the LF. We improved these estimators for the studies of the distant Universe, and examined their performances for various classes of functional forms by Monte Carlo simulations. We also applied these estimation methods to the mock 2dF redshift survey catalog prepared by Cole et al. (1998). We found that $1/V_{\rm max}$ estimator yields a completely unbiased result if there is no inhomogeneity, but is not robust against clusters or voids. This is consistent with the well-known results, and we did not confirm the bias trend of $1/V_{\rm max}$ estimator claimed by Willmer (1997) in the case of homogeneous sample. We also found that the other three maximum-likelihood type estimators are quite robust and give consistent results with each other. In practice we recommend Cho{\l}oniewski's estimator for two reasons: 1. it simultaneously provides the shape and normalization of the LF; 2. it is the fastest among these four estimators, because of the algorithmic simplicity. Then, we analyzed the photometric redshift data of the Hubble Deep Field prepared by Fern\'{a}ndez-Soto et al. (1999) using the above four methods. We also derived luminosity density $\rho_{\rm L}$ at $B$- and $I$-band. Our $B$-band estimation is roughly consistent with that of Sawicki, Lin, & Yee (1997), but a few times lower at $2.0 < z < 3.0$. The evolution of $\rho_{\rm L}(I)$ is found to be less prominent.Comment: To appear in ApJS July 2000 issue. 36 page

Lifting the Dusty Veil With Near- and Mid-Infrared Photometry: III. Two-Dimensional Extinction Maps of the Galactic Midplane Using the Rayleigh-Jeans Color Excess Method

We provide new, high-resolution A(Ks) extinction maps of the heavily reddened Galactic midplane based on the Rayleigh-Jeans Color Excess ("RJCE") method. RJCE determines star-by-star reddening based on a combination of near- and mid-infrared photometry. The new RJCE-generated maps have 2 x 2 arcmin pixels and span some of the most severely extinguished regions of the Galaxy -- those covered with Spitzer+IRAC imaging by the GLIMPSE-I, -II, -3D, and Vela-Carina surveys, from 256<l<65 deg and, in general, for |b| <= 1-1.5 deg (extending up to |b|<=4 deg in the bulge). Using RJCE extinction measurements, we generate dereddened color-magnitude diagrams and, in turn, create maps based on main sequence, red clump, and red giant star tracers, each probing different distances and thereby providing coarse three-dimensional information on the relative placement of dust cloud structures. The maps generated from red giant stars, which reach to ~18-20 kpc, probe beyond most of the Milky Way extinction in most directions and provide close to a "total Galactic extinction" map -- at minimum they provide high angular resolution maps of lower limits on A(Ks). Because these maps are generated directly from measurements of reddening by the very dust being mapped, rather than inferred on the basis of some less direct means, they are likely the most accurate to date for charting in detail the highly patchy differential extinction in the Galactic midplane. We provide downloadable FITS files and an IDL tool for retrieving extinction values for any line of sight within our mapped regions.Comment: 23 pages, 5 figures, accepted for publication in ApJ

Photometric AGN reverberation mapping - an efficient tool for BLR sizes, black hole masses and host-subtracted AGN luminosities

Photometric reverberation mapping employs a wide bandpass to measure the AGN continuum variations and a suitable band, usually a narrow band (NB), to trace the echo of an emission line in the broad line region (BLR). The narrow band catches both the emission line and the underlying continuum, and one needs to extract the pure emission line light curve. We performed a test on two local AGNs, PG0003+199 (=Mrk335) and Ark120, observing well-sampled broad- (B, V) and narrow-band light curves with the robotic 15cm telescope VYSOS-6 on Cerro Armazones, Chile. In PG0003+199, H_alpha dominates the flux in the NB by 85%, allowing us to measure the time lag of H_alpha against B without the need to correct for the continuum contribution. In Ark120, H_beta contributes only 50% to the flux in the NB. The cross correlation of the B and NB light curves shows two distinct peaks of similar strength, one at lag zero from the autocorrelated continuum and one from the emission line at tau_cent = 47.5 +/- 3.4 days. We constructed a synthetic H_beta light curve, by subtracting a scaled V light curve, which traces the continuum, from the NB light curve. The cross correlation of this synthetic H_beta light curve with the B light curve shows only one major peak at tau_cent = 48.0 +/- 3.3 days, while the peak from the autocorrelated continuum at lag zero is absent. We conclude that, as long as the emission line contributes at least 50% to the bandpass, the pure emission line light curve can be reconstructed from photometric monitoring data so that the time lag can be measured. For both objects the lags we find are consistent with spectroscopic reverberation results. While the dense sampling (median 2 days) enables us to determine tau_cent with small (10%) formal errors, we caution that gaps in the light curves may lead to much larger systematic uncertainties. (Abstract shortened, see the manuscript.)Comment: 12 pages, 15 figures, accepted for publication in Astronomy and Astrophysic

A robust method for fitting peculiar velocity field models

We present a new method for fitting peculiar velocity models to complete flux limited magnitude-redshifts catalogues, using the luminosity function of the sources as a distance indicator.The method is characterised by its robustness. In particular, no assumptions are made concerning the spatial distribution of sources and their luminosity function. Moreover, selection effects in redshift are allowed. Furthermore the inclusion of additional observables correlated with the absolute magnitude -- such as for example rotation velocity information as described by the Tully-Fisher relation -- is straightforward. As an illustration of the method, the predicted IRAS peculiar velocity model characterised by the density parameter beta is tested on two samples. The application of our method to the Tully-Fisher MarkIII MAT sample leads to a value of beta=0.6 \pm 0.125, fully consistent with the results obtained previously by the VELMOD and ITF methods on similar datasets. Unlike these methods however, we make a very conservative use of the Tully-Fisher information. Specifically, we require to assume neither the linearity of the Tully-Fisher relation nor a gaussian distribution of its residuals. Moreover, the robustness of the method implies that no Malmquist corrections are required. A second application is carried out, using the fluxes of the IRAS 1.2 Jy sample as the distance indicator. In this case the effective depth of the volume in which the velocity model is compared to the data is almost twice the effective depth of the MarkIII MAT sample. The results suggest that the predicted IRAS velocity model, while successful in reproducing locally the cosmic flow, fails to describe the kinematics on larger scales.Comment: 10 pages, 14 figures, MNRAS in pres