A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data

We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l, followed by a particular deterministic rescaling operation of the sky signal. The acceptance probability for this joint move depends on the sky map only through the difference of chi-squared between the original and proposed sky sample, which is close to unity in the low signal-to-noise regime. The algorithm is completed by alternating this move with a standard Gibbs move. Together, these two proposals constitute a computationally efficient algorithm for mapping out the full joint CMB posterior, both in the high and low signal-to-noise regimes.Comment: Submitted to Ap

A Hybrid N-body--Coagulation Code for Planet Formation

We describe a hybrid algorithm to calculate the formation of planets from an initial ensemble of planetesimals. The algorithm uses a coagulation code to treat the growth of planetesimals into oligarchs and explicit N-body calculations to follow the evolution of oligarchs into planets. To validate the N-body portion of the algorithm, we use a battery of tests in planetary dynamics. Several complete calculations of terrestrial planet formation with the hybrid code yield good agreement with previously published calculations. These results demonstrate that the hybrid code provides an accurate treatment of the evolution of planetesimals into planets.Comment: Astronomical Journal, accepted; 33 pages + 11 figure

The upper critical field of filamentary Nb3Sn conductors

We have examined the upper critical field of a large and representative set of present multi-filamentary Nb3Sn wires and one bulk sample over a temperature range from 1.4 K up to the zero field critical temperature. Since all present wires use a solid-state diffusion reaction to form the A15 layers, inhomogeneities with respect to Sn content are inevitable, in contrast to some previously studied homogeneous samples. Our study emphasizes the effects that these inevitable inhomogeneities have on the field-temperature phase boundary. The property inhomogeneities are extracted from field-dependent resistive transitions which we find broaden with increasing inhomogeneity. The upper 90-99 % of the transitions clearly separates alloyed and binary wires but a pure, Cu-free binary bulk sample also exhibits a zero temperature critical field that is comparable to the ternary wires. The highest mu0Hc2 detected in the ternary wires are remarkably constant: The highest zero temperature upper critical fields and zero field critical temperatures fall within 29.5 +/- 0.3 T and 17.8 +/- 0.3 K respectively, independent of the wire layout. The complete field-temperature phase boundary can be described very well with the relatively simple Maki-DeGennes model using a two parameter fit, independent of composition, strain state, sample layout or applied critical state criterion.Comment: Accepted Journal of Applied Physics Few changes to shorten document, replaced eq. 7-

Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps

We apply a previously developed Gibbs sampling framework to the foreground corrected 3-yr WMAP polarization data and compute the power spectrum and residual foreground template amplitude posterior distributions. We first analyze the co-added Q- and V-band data, and compare our results to the likelihood code published by the WMAP team. We find good agreement, and thus verify the numerics and data processing steps of both approaches. However, we also analyze the Q- and V-bands separately, allowing for non-zero EB cross-correlations and including two individual foreground template amplitudes tracing synchrotron and dust emission. In these analyses, we find tentative evidence of systematics: The foreground tracers correlate with each of the Q- and V-band sky maps individually, although not with the co-added QV map; there is a noticeable negative EB cross-correlation at l <~ 16 in the V-band map; and finally, when relaxing the constraints on EB and BB, noticeable differences are observed between the marginalized band powers in the Q- and V-bands. Further studies of these features are imperative, given the importance of the low-l EE spectrum on the optical depth of reionization tau and the spectral index of scalar perturbations n_s.Comment: 5 pages, 4 figures, submitted to ApJ

Ground-State SiO Maser Emission Toward Evolved Stars

We have made the first unambiguous detection of vibrational ground-state maser emission from SiO toward six evolved stars. Using the Very Large Array, we simultaneously observed the v=0, J=1-0, 43.4-GHz, ground-state and the v=1, J=1-0, 43.1-GHz, first excited-state transitions of SiO toward the oxygen-rich evolved stars IRC+10011, o Ceti, W Hya, RX Boo, NML Cyg, and R Cas and the S-type star chi Cyg. We detected at least one v=0 SiO maser feature from six of the seven stars observed, with peak maser brightness temperatures ranging from 10,000 K to 108,800 K. In fact, four of the seven v=0 spectra show multiple maser peaks, a phenomenon which has not been previously observed. Ground-state thermal emission was detected for one of the stars, RX Boo, with a peak brightness temperature of 200 K. Comparing the v=0 and the v=1 transitions, we find that the ground-state masers are much weaker with spectral characteristics different from those of the first excited-state masers. For four of the seven stars the velocity dispersion is smaller for the v=0 emission than for the v=1 emission, for one star the dispersions are roughly equivalent, and for two stars (one of which is RX Boo) the velocity spread of the v=0 emission is larger. In most cases, the peak flux density in the v=0 emission spectrum does not coincide with the v=1 maser peak. Although the angular resolution of these VLA observations were insufficient to completely resolve the spatial structure of the SiO emission, the SiO spot maps produced from the interferometric image cubes suggest that the v=0 masers are more extended than their v=1 counterparts

The effect of interleukin-10 and transforming growth factor beta-1 on HLA-DR expression in colonic epithelial cells.

The aim of this study was to assess whether interleukin-10 (IL-10) and/or transforming growth factor beta-1 (TGFbeta1) downregulate HLA-DR expression using the HT29 cell line as a model of colonic epithelial cells. HLA-DR expression was induced in HT29 cells with gamma-interferon. The effects of IL-10 alone, TGFbeta1 alone, and IL-10 and TGFbeta1 in combination were studied. HLA-DR expression was assessed using flow cytometric analysis. Gamma-interferon induced HLA-DR expression in a dose-dependent fashion. In the absence of gamma-interferon, neither IL-10 nor TGFbeta1 induced HLA-DR expression. In isolation, neither IL-10 nor TGFbeta1 downregulated HLA-DR expression. When IL-10 and TGFbeta1 were added in combination, small (6-30%) statistically significant reductions in HLA-DR expression were seen. The biological significance is unclear

Bayesian Power Spectrum Analysis of the First-Year WMAP data

We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotropies. We present the complete probability distributions for each C_l including any non-Gaussianities of the power spectrum likelihood. While we find good overall agreement with the previously published WMAP spectrum, our analysis uncovers discrepancies in the power spectrum estimates at low l multipoles. For example we claim the best-fit Lambda-CDM model is consistent with the C_2 inferred from our combined Q+V+W analysis with a 10% probability of an even larger theoretical C_2. Based on our exact analysis we can therefore attribute the "low quadrupole issue" to a statistical fluctuation.Comment: 5 pages. 4 figures. For additional information and data see http://www.astro.uiuc.edu/~iodwyer/research#wma

Incidence and treatment costs of severe bacterial infections among people who inject heroin: A cohort study in South London, England.

BACKGROUND: People who inject drugs often get bacterial infections. Few longitudinal studies have reported the incidence and treatment costs of these infections. METHODS: For a cohort of 2335 people who inject heroin entering treatment for drug dependence between 2006 and 2017 in London, England, we reported the rates of hospitalisation or death with primary causes of cutaneous abscess, cellulitis, phlebitis, septicaemia, osteomyelitis, septic arthritis, endocarditis, or necrotising fasciitis. We compared these rates to the general population. We also used NHS reference costs to calculate the cost of admissions. RESULTS: During a median of 8.0 years of follow-up, 24 % of patients (570/2335) had a severe bacterial infection, most commonly presenting with cutaneous abscesses or cellulitis. Bacterial infections accounted for 13 % of all hospital admissions. The rate was 73 per 1000 person-years (95 % CI 69-77); 50 times the general population, and the rate remained high throughout follow-up. The rate of severe bacterial infections for women was 1.50 (95 % CI 1.32-1.69) times the rate for men. The mean cost per admission was £4980, and we estimate that the annual cost of hospital treatment for people who inject heroin in London is £4.5 million. CONCLUSIONS: People who inject heroin have extreme and long-term risk of severe bacterial infections