If neutrinos have self-interactions, these will induce scatterings between astrophysical and cosmic neutrinos. Prior work proposed to look for possible resulting resonance features in astrophysical neutrino spectra in order to seek a neutrino self-interaction which can be either diagonal in the neutrino flavor space or couple different neutrino flavors. The calculation of the astrophysical spectra involves either a Monte Carlo simulation or a computationally intensive numerical integration of an integro-partial-differential equation. As a result only limited regions of the neutrino self-interaction parameter space have been explored, and only flavor-diagonal self-interactions have been considered. Here, we present a fully analytic form for the astrophysical neutrino spectra for arbitrary neutrino number and arbitrary self-coupling matrix that accurately obtains the resonance features in the observable neutrino spectra. The results can be applied to calculations of the diffuse supernova neutrino background and of the spectrum from high-energy astrophysical neutrino sources. We illustrate with a few examples.

Popular accounts of exciting discoveries often draw students to physics and astronomy, but at the introductory level it is challenging to connect with these in a meaningful way. The use of real astronomical data in the classroom can help bridge this gap and build valuable quantitative and scientific reasoning skills. This paper presents a strategy for studying Hubble's Law and the accelerating expansion of the universe using actual data. Along with understanding the physical concepts, an explicit goal is to develop skills for analyzing data in terms of a model.

After inflation, a period of preheating may have produced a stochastic background of high frequency gravitational waves (GWs) that would persist until today. The nature of the inflaton's coupling to Standard Model or other fields is unknown, so it is useful to ask what features such fields may typically have, and how these affect predictions for the GW's produced. Here we consider the inflaton to be coupled to a light scalar field, and show that even a very small quartic self-interaction term will reduce the amplitude of the GW spectrum. For self-coupling $\lambda_{\chi} \gtrsim g^2$, where $g^2$ is the inflaton-scalar coupling, the peak energy density goes as $\Omega_{\rm gw}^{(\lambda_{\chi})} / \Omega_{\rm gw}^{(\lambda_{\chi}=0)} \sim (g^2/\lambda_{\chi})^{2}$. A consequence is that if the universe reheats through an inflaton-Higgs coupling then the spectrum would be suppressed but the dynamics would be sensitive to the Higgs potential near the energy scale of inflation.

In hidden sector models with an extra U(1) gauge group, new fields can interact with the Standard Model only through gauge kinetic mixing and the Higgs portal. After the U(1) is spontaneously broken, these interactions couple the resultant cosmic strings to Standard Model particles. We calculate the spectrum of radiation emitted by these "dark strings" in the form of Higgs bosons, Z bosons, and Standard Model fermions assuming that string tension is above the TeV scale. We also calculate the scattering cross sections of Standard Model fermions on dark strings due to the Aharonov-Bohm interaction. These radiation and scattering calculations will be applied in a subsequent paper to study the cosmological evolution and observational signatures of dark strings.

We consider the Standard Model extended by a hidden sector U(1)_X symmetry that is spontaneously broken at the TeV scale by the vacuum expectation value of an additional scalar field. We study "dark string" solutions in this model and their properties due to the Higgs portal and gauge kinetic mixing operators. We find that dark strings effectively interact with Higgs and Z bosons by linear couplings, and with leptons and baryons via Aharonov-Bohm couplings, thus possibly leading to new cosmological constraints on dark matter models with spontaneously broken extra U(1) symmetry factors.

We examine the colour-magnitude relation of approximately 5000 Brightest Cluster Galaxies (BCGs) in the Sloan Digital Sky Survey, and compare with non-BCG E/S0 galaxies. The colour-magnitude and colour-sigma (velocity dispersion) relations are flatter in slope (by a factor of about 2) for BCGs than for non-BCG E/S0s, and the BCGs also tend to be redder by 0.01 magnitudes in g-r. We investigate radial colour gradients in both samples, using the ratio of the de Vaucouleurs radii in the g and r bands. We find BCGs have significantly flatter (by 23%) mean colour gradients than other high luminosity E/S0s. In early-type galaxies, the colour gradients are strongest at intermediate luminosities of Mr=-22. Colour gradients in E/S0s increase with radius (up to 10kpc) and are negatively correlated with 10sigma + Mr (velocity dispersion relative to luminosity) and with mass density. The gradients also tend to decrease with increasing stellar age. These trends are weak or not seen in BCGs, in which the mean colour gradient is low whatever the other properties. We discuss possible explanations, which involve a greater amount of dry merging in the formation history of the BCGs.

We provide fits to the distribution of galaxy luminosity, size, velocity dispersion and stellar mass as a function of concentration index C_r and morphological type in the SDSS. We also quantify how estimates of the fraction of `early' or `late' type galaxies depend on whether the samples were cut in color, concentration or light profile shape, and compare with similar estimates based on morphology. Our fits show that Es account for about 20% of the r-band luminosity density, rho_Lr, and 25% of the stellar mass density, rho_*; including S0s and Sas increases these numbers to 33% and 40%, and 50% and 60%, respectively. Summed over all galaxy types, we find rho_* ~ 3 * 10^8 M_Sun Mpc^{-3} at z ~ 0. This is in good agreement with expectations based on integrating the star formation history. However, compared to most previous work, we find an excess of objects at large masses, up to a factor of ~ 10 at M_* ~ 5*10^{11} M_Sun. The stellar mass density further increases at large masses if we assume different IMFs for Es and spiral galaxies, as suggested by some recent chemical evolution models, and results in a better agreement with the dynamical mass function. We also show that the trend for ellipticity to decrease with luminosity is primarily because the E/S0 ratio increases at large L. However, the most massive galaxies, M_* > 5 * 10^{11} M_Sun, are less concentrated and not as round as expected if one extrapolates from lower L, and they are not well-fit by pure deVaucouleur laws. This suggests formation histories with recent radial mergers. Finally, we show that the age-size relation is flat for Es of fixed dynamical mass, but, at fixed M_dyn, S0s and Sas with large sizes tend to be younger. Explaining this difference between E and S0 formation is a new challenge for models of early-type galaxy formation.

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

We select a sample of 70378 E/S0 (early-type) galaxies at 0<z<0.36 from the Sloan Digital Sky Survey, excluding disk and star-forming galaxies. We estimate g and r magnitudes in the observer- and rest-frames directly from the SDSS DR6 spectra; this provides an object-by-object estimate of the k-correction. We use the k-corrections from the spectra to study the evolution of the rest-frame colour-magnitude relation (CMR) and colour-sigma (velocity dispersion) relation (CsigmaR). The evolution is very sensitive to the k-correction. Both the CMR and CsigmaR evolve blueward with increasing redshift, approximately in agreement with passive evolution models with age approx. 12 Gyr. The slope and zero-point of the CMR depend on whether colours were defined in fixed physical or angular apertures, a consequence of the fact that the centers of these objects tend to be redder: the relation is steeper for fixed angular apertures. One the other hand, the CsigmaR slope does not show this dependence on the aperture in which the colour was defined, suggesting that colour gradients are correlated with residuals from the sigma-M_r relation. As these residuals are age indicators, our findings suggest that colour gradients correlate with the age of the stellar population.

From a sample of ~50000 early-type galaxies from the SDSS, we measured the traditional Fundamental Plane in four bands. We then replaced luminosity with stellar mass, and measured the "stellar mass" FP. The FP steepens slightly as one moves from shorter to longer wavelengths: the orthogonal fit has slope 1.40 in g and 1.47 in z. The FP is thinner at longer wavelengths: scatter is 0.062 dex in g, 0.054 dex in z. The scatter is larger at small galaxy sizes/masses; at large masses measurement errors account for essentially all of the observed scatter. The FP steepens further when luminosity is replaced with stellar mass, to slope ~ 1.6. The intrinsic scatter also reduces further, to 0.048 dex. Since color and stellar mass-to-light ratio are closely related, this explains why color can be thought of as the fourth FP parameter. However, the slope of the stellar mass FP remains shallower than the value of 2 associated with the virial theorem. This is because the ratio of dynamical to stellar mass increases at large masses as M_d^0.17. The face-on view of the stellar mass kappa-space suggests that there is an upper limit to the stellar density for a given dynamical mass, and this decreases at large masses: M_*/R_e^3 ~ M_d^-4/3. We also study how the estimated coefficients a and b of the FP are affected by other selection effects (e.g. excluding small sigma biases a high; excluding fainter L biases a low). These biases are seen in FPs which have no intrinsic curvature, so the observation that a and b scale with L and sigma is not, by itself, evidence that the Plane is warped. We show that the FP appears to curve sharply downwards at the small mass end, and more gradually downwards towards larger masses. Whereas the drop at small sizes is real, most of the latter effect is due to correlated errors.

We select a sample of about 50,000 early-type galaxies from the Sloan Digital Sky Survey (SDSS), calibrate fitting formulae which correct for known problems with photometric reductions of extended objects, apply these corrections, and then measure a number of pairwise scaling relations in the corrected sample. We show that, because they are not seeing corrected, the use of Petrosian-based quantities in magnitude limited surveys leads to biases, and suggest that this is one reason why Petrosian-based analyses of BCGs have failed to find significant differences from the bulk of the early-type population. These biases are not present when seeing-corrected parameters derived from deVaucouleur fits are used. Most of the scaling relations we study show evidence for curvature: the most luminous galaxies have smaller velocity dispersions, larger sizes, and fainter surface brightnesses than expected if there were no curva-ture. These statements remain true if we replace luminosities with stellar masses; they suggest that dissipation is less important at the massive end. There is curvature in the dynamical to stellar mass relation as well: the ratio of dynamical to stellar mass increases as stellar mass increases, but it curves upwards from this scaling both at small and large stellar masses. In all cases, the curvature at low masses becomes apparent when the sample becomes dominated by objects with stellar masses smaller than 3 x 10^10 M_Sun. We quantify all these trends using second order polynomials; these generally provide significantly better description of the data than linear fits, except at the least luminous end.

We used the Advanced Camera for Surveys on board the Hubble Space Telescope to obtain high resolution i-band images of the centers of 23 single galaxies, which were selected because they have SDSS velocity dispersions larger than 350 km/s. The surface brightness profiles of the most luminous of these objects (M_i<-24) have well-resolved `cores' on scales of 150-1000 pc, and share similar properties to BCGs. The total luminosity of the galaxy is a better predictor of the core size than is the velocity dispersion. The correlations of luminosity and velocity dispersion with core size agree with those seen in previous studies of galaxy cores. Because of high velocity dispersions, our sample of galaxies can be expected to harbor the most massive black holes, and thus have large cores with large amounts of mass ejection. The mass-deficits inferred from core-Sersic fits to the surface-brightness profiles are approximately double the black-hole masses inferred from the M_bh-sigma relation and the same as those inferred from the M_bh-L relation. The less luminous galaxies (M_i>-23) tend to have steeper `power-law' inner profiles, higher-ellipticity, diskier isophotes, and bulge-to-total ratios of order 0.5 -- all of which suggest that they are `fast-rotators' and rotational motions could have contaminated the velocity dispersion estimate. There are obvious dust features within about 300 pc of the center in about 35% of the sample, predominantly in power-law rather than core galaxies.

We study a sample of 43 early-type galaxies, selected from the SDSS because they appeared to have velocity dispersion > 350 km/s. High-resolution photometry in the SDSS i passband using HRC-ACS on board the HST shows that just less than half of the sample is made up of superpositions of two or three galaxies, so the reported velocity dispersion is incorrect. The other half of the sample is made up of single objects with genuinely large velocity dispersions. None of these objects has sigma larger than 426 +- 30 km/s. These objects define rather different relations than the bulk of the early-type galaxy population: for their luminosities, they are the smallest, most massive and densest galaxies in the Universe. Although the slopes of the scaling relations they define are rather different from those of the bulk of the population, they lie approximately parallel to those of the bulk "at fixed sigma". These objects appear to be of two distinct types: the less luminous (M_r>-23) objects are rather flattened and extremely dense for their luminosities -- their properties suggest some amount of rotational support and merger histories with abnormally large amounts of gaseous dissipation. The more luminous objects (M_r<-23) tend to be round and to lie in or at the centers of clusters. Their properties are consistent with the hypothesis that they are BCGs. Models in which BCGs form from predominantly radial mergers having little angular momentum predict that they should be prolate. If viewed along the major axis, such objects would appear to have abnormally large sigma for their sizes, and to be abnormally round for their luminosities. This is true of the objects in our sample once we account for the fact that the most luminous galaxies (M_r<-23.5), and BCGs, become slightly less round with increasing luminosity.

We present Hobby-Eberly Telescope (HET) observations for galaxies at redshift z < 0.3 from the Sloan Digital Sky Survey (SDSS) showing large velocity dispersions while appearing to be single galaxies in HST images. The high signal-to-noise HET spectra provide more definitive velocity dispersions. The maximum velocity dispersion we find is 444 km/s. Emission-line widths in QSOs indicate that black holes can exist with masses exceeding 5 billion solar masses, implying velocity dispersions greater than 500 km/s by the local black hole mass - velocity dispersion relationship. This suggests either that QSO black hole masses are overestimated or that the black hole - bulge relationship changes at high black hole mass. The latter option is consistent with evidence that the increase in velocity dispersion with luminosity levels off for the brightest elliptical galaxies.

Black hole masses are tightly correlated with the stellar velocity dispersions of the bulges which surround them, and slightly less-well correlated with the bulge luminosity. It is common to use these correlations to estimate the expected abundance of massive black holes. This is usually done by starting from an observed distribution of velocity dispersions or luminosities and then changing variables. This procedure neglects the fact that there is intrinsic scatter in these black hole mass--observable correlations. Accounting for this scatter results in estimates of black hole abundances which are larger by almost an order of magnitude at masses >10^9 M_sun. Including this scatter is particularly important for models which seek to infer quasar lifetimes and duty cycles from the local black hole mass function. However, even when scatter has been accounted for, the M_bh-sigma relation predicts fewer massive black holes than does the M_bh-L relation. This is because the sigma-L relation in the black hole samples currently available is inconsistent with that in the SDSS sample from which the distributions of L or sigma are based: the black hole samples have smaller L for a given sigma, or larger sigma for a given L. The sigma-L relation in the black hole samples is similarly discrepant with that in other samples of nearby early-type galaxies. This suggests that current black hole samples are biased: if this is a selection rather than physical effect, then the M_bh-sigma and M_bh-L relations currently in the literature are also biased from their true values.

It is common to estimate black hole abundances by using a measured correlation between black hole mass and another more easily measured observable such as the velocity dispersion or luminosity of the surrounding bulge. The correlation is used to transform the distribution of the observable into an estimate of the distribution of black hole masses. However, different observables provide different estimates: the Mbh-sigma relation predicts fewer massive black holes than does the Mbh-L relation. This is because the sigma-L relation in black hole samples currently available is inconsistent with that in the SDSS sample, from which the distributions of L or sigma are based: the black hole samples have smaller L for a given sigma or have larger sigma for a given L. This is true whether L is estimated in the optical or in the NIR. If this is a selection rather than physical effect, then the Mbh-sigma and Mbh-L relations currently in the literature are also biased from their true values. We provide a framework for describing the effect of this bias. We then combine it with a model of the bias to make an estimate of the true intrinsic relations. While we do not claim to have understood the source of the bias, our simple model is able to reproduce the observed trends. If we have correctly modeled the selection effect, then our analysis suggests that the bias in the <Mbh|sigma> relation is likely to be small, whereas the <Mbh|L> relation is biased towards predicting more massive black holes for a given luminosity. In addition, it is likely that the Mbh-L relation is entirely a consequence of more fundamental relations between Mbh and sigma, and between sigma and L. The intrinsic relation we find suggests that at fixed luminosity, older galaxies tend to host more massive black holes.

The size-luminosity relation of early-type Brightest Cluster Galaxies (BCGs), R_e ~ L^0.88, is steeper than that for the bulk of the early-type galaxy population, for which R_e ~ L^0.68. In addition, although BCGs are hardly offset from the Fundamental Plane defined by the bulk of the early-type population, they show considerably smaller scatter. The larger than expected sizes of BCGs, and the increased homogeneity, are qualitatively consistent with models which seek to explain the colors of the most massive galaxies by invoking dry dissipationless mergers, since dissipation tends to reduce the sizes of galaxies, and wet mergers which result in star formation would tend to increase the scatter in luminosity at fixed size and velocity dispersion. Furthermore, BCGs define the same g-r color-magnitude relation as the bulk of the early-type population. If BCGs formed from dry mergers, then BCG progenitors must have been red for their magnitudes, suggesting that they hosted older stellar populations than typical for their luminosities. Our findings have two other consequences. First, the R_e-L relation of the early-type galaxy population as a whole (i.e., normal plus BCG) exhibits some curvature. Some of this curvature must be a consequence of the fact that an increasing fraction of the most luminous galaxies are BCGs. The second consequence is suggested by the fact that, despite following a steeper size-luminosity relation, BCGs tend to define a tight relation between dynamical mass R_e sigma^2/G and luminosity. As consequence, we find that BCGs define a shallower sigma-L relation than the bulk of the early-type galaxy population.