We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

We present improved results of the measurement of the correlation between galaxies and the intergalactic medium (IGM) transmission at the end of reionisation. We have gathered a sample of $13$ spectroscopically confirmed Lyman-break galaxies (LBGs) and $21$ Lyman-$\alpha$ emitters (LAEs) at angular separations $20'' \lesssim \theta \lesssim 10'$ ($\sim 0.1-4$ pMpc at $z\sim 6$) from the sightlines to $8$ background $z\gtrsim 6$ quasars. We report for the first time the detection of an excess of Lyman-$\alpha$ transmission spikes at $\sim 10-60$ cMpc from LAEs ($3.2\sigma$) and LBGs ($1.9\sigma$). We interpret the data with an improved model of the galaxy-Lyman-$\alpha$ transmission and two-point cross-correlations which includes the enhanced photoionisation due to clustered faint sources, enhanced gas densities around the central bright objects and spatial variations of the mean free path. The observed LAE(LBG)-Lyman-$\alpha$ transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionisation to $\langle f_{\rm esc} \rangle_{M_{\rm UV}<-12} = 0.14_{-0.05}^{+0.28}\,(0.23_{-0.12}^{+0.46})$. We investigate if the 2PCCF measurement can determine whether bright or faint galaxies are the dominant contributors to reionisation. Our results show that a contribution from faint galaxies ($M_{\rm UV} > -20 \, (2\sigma)$) is necessary to reproduce the observed 2PCCF and that reionisation might be driven by different sub-populations around LBGs and LAEs at $z\sim 6$.

We anticipate noise from the Laser Interferometer Space Antenna (LISA) will exhibit nonstationarities throughout the duration of its mission due to factors such as antenna repointing, cyclostationarities from spacecraft motion, and glitches as highlighted by LISA Pathfinder. In this paper, we use a surrogate data approach to test the stationarity of a time series, with the goal of identifying noise nonstationarities in the future LISA mission. This will be necessary for determining how often the LISA noise power spectral density (PSD) will need to be updated for parameter estimation routines. We conduct a thorough simulation study illustrating the power/size of various versions of the hypothesis tests, and then apply these approaches to differential acceleration measurements from LISA Pathfinder. We also develop a data analysis strategy for addressing nonstationarities in the LISA PSD, where we update the noise PSD over time, while simultaneously conducting parameter estimation, with a focus on planned data gaps. We show that assuming stationarity when noise is nonstationary leads to systematic biases and large posterior variances in parameter estimates for galactic white dwarf binary gravitational wave signals.

Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$\alpha$ emitters' (LAEs) within the proximity zone of quasar J0836 at $z=5.802$ (\textit{Aerith A, B} and \textit{C}). \textit{Aerith A}, located behind the quasar with an impact parameter $D_\perp = 278$ pkpc, provides the first detection of a Lyman-$\alpha$ transverse proximity effect. We model the transmission and show it constrains the onset of J0836's quasar phase to $0.2 \text{Myr}<t<20\text{Myr}$ in the past. The second object, \textit{Aerith B} at a distance $D=750$ pkpc from the quasar, displays a bright, broad double-peaked \lal emission line. Based on relations calibrated at $z\leq3$, the peak separation implies a low ionizing $f_{\text{esc}} \lesssim 1\%$, the most direct such constraint on a reionization-era galaxy. We fit the Ly-$\alpha$ line with an outflowing shell model, finding a completely typical central density $\text{log N}_{\text{HI}}/\text{cm}^{-2} = 19.3_{-0.2}^{+0.8}$, outflow velocity $v=16_{-11}^{+4}$ km s$^{-1}$, and gas temperature $\text{log} T/\text{K} = 3.8_{-0.7}^{+0.8}$ compared to $2<z<3$ analogue LAEs. Finally, we detect an emission line at $\lambda=8177$ \AA\ in object \textit{Aerith C} which, if it is \lal at $z=5.726$, would correspond closely with the end of the quasar's proximity zone ($\Delta z<0.02$ from the boundary) and suggests the quasar influences the IGM up to $\sim85$ cMpc away, making it the largest quasar proximity zone. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during the epoch of reionization.

We present the results of a model-independent investigation of the rest-frame UV spectra from a comprehensive sample of $394$ quasars in the redshift range $1.5\leq z \leq 7.5$. We fit the main Broad Emission Lines (BELs) in the rest-frame range $1280 \text{\AA} \leq \lambda \leq 3000 \text{\AA}$ (O I, C II, Si IV, C III, C IV and Mg II) with a lightly-supervised spline fitting technique. Redshifts are derived from the peaks of each fitted BEL and used to compute relative velocity shifts. We show that our method gives unbiased velocity shifts and is insensitive to spectral resolution and instrumental parameters. It is found that the average blueshift of the \cfour\, line with respect to several low-ionisation lines in luminosity-matched samples does not significantly evolve over $1.5\gtrsim z\gtrsim6$. However, the average blueshift increases significantly by a factor $\sim 2.5$ at $z\gtrsim 6$. We propose that this redshift evolution can be explained by \cfour\, winds launched perpendicularly to an accretion disk with increased torus opacity at high-redshift, coupled with a potential orientation-driven selection bias. Our results open new exciting avenues of investigation into young quasars in the reionisation epoch.

Using integral field spectroscopy, we demonstrate that gravity-sensitive absorption features in the zJ-band (0.9--1.35 \micron) can constrain the low-mass stellar initial mass function (IMF) in the cores of two elliptical galaxies, M85 and M87. Compared to the visible bands, the near-infrared (NIR) is more sensitive to light from low-mass dwarf stars, whose relative importance is the primary subject of the debate over IMF variations in nearby galaxies. Our analysis compares the observed spectra to the latest stellar population synthesis models by employing two different methods: equivalent widths and spectral fitting. We find that the IMF slopes in M85 are similar to the canonical Milky Way IMF with a median IMF-mismatch parameter $\alpha_{K} = 1.26$. In contrast, we find that the IMF in M87 is steeper than a Salpeter IMF with $\alpha_{K} = 2.77$. The derived stellar population parameters, including the IMF slopes, are consistent with those from recent results in the visible bands based on spectroscopic and kinematic techniques. Certain elemental abundances, e.g. Na and Fe, have dramatic effects on the IMF-sensitive features and therefore the derived IMF slopes. We show evidence for a high [Na/H] $\sim$ 0.65 dex in the core of M85 from two independent \ion{Na}{1} absorption features. The high Na abundance may be the result of a recent galactic merger involving M85. This suggests that including [Na/H] in the stellar population model parameters is critical for constraining the IMF slopes in M85. These results confirm the viability of using NIR absorption features to investigate IMF variation in nearby galaxies.

We present a sample of 12 Quasi-Stellar Objects (QSOs) potentially acting as strong gravitational lenses on background Emission Line Galaxies (ELG) or Lyman-$\alpha$ Emitters (LAEs) selected selected through a systematic search of the 297301 QSOs in the Sloan Digital Sky Survey (SDSS)-III Data Release 12. Candidates are identified by looking for compound spectra, where emission lines at a redshift larger than that of the quasar can be identified in the residuals after a QSO spectral template is subtracted from the observed spectra. The narrow diameter of BOSS fibers (2") then ensures that the object responsible for the additional emission lines must lie close to the line of sight of the QSO and hence provides a high probability of lensing. Among the 12 candidates identified, 9 have definite evidence for the presence of a background ELG identified by at least 4 higher-redshift nebular emission lines. The remaining $3$ probable candidates present a strong asymmetrical emission line attributed to a background Lyman-$\alpha$ emitter (LAE). The QSO-ELG (QSO-LAE) lens candidates have QSO lens redshifts in the range $0.24\lesssim z_{QSO} \lesssim 0.66$ ($0.75 \lesssim z_{QSO} \lesssim 1.23$ ) and background galaxy redshifts in the range $0.48 \lesssim z_{S,ELG} \lesssim 0.94$ ($2.17 \lesssim z_{S,LAE} \lesssim 4.48$). We show that the algorithmic search is complete at >90% for QSO-ELG systems, whereas it falls at 40-60% for QSO-LAE, depending on the redshift of the source. Upon confirmation of the lensing nature of the systems, this sample may quadruple the number of known QSOs acting as strong lenses. We have determined the completeness of our search, which allows future studies to compute lensing probabilities of galaxies by QSOs and differentiate between different QSO models. (Abridged)

We present a new method to study the contribution of faint sources to the UV background using the 1D correlation of metal absorbers with the intergalactic medium (IGM) transmission in a Quasi Stellar Object (QSO) sightline. We take advantage of a sample of $25$ high signal-to-noise ratio QSO spectra to retrieve $150$ triply-ionised carbon (\cfour) absorbers at $4.5\lesssim z\lesssim 6.2$, of which $37$ systems whose expected H{~\small I} absorption lie in the Lyman-$\alpha$ forest. We derive improved constraints on the cosmic density of \cfour \,at $4.3< z < 6.2$ and infer from abundance-matching that \cfour \,absorbers trace $M_{\text{UV}}\lesssim -16$ galaxies. Correlation with the Lyman-$\alpha$ forest of the QSOs indicates that these objects are surrounded by a highly opaque region at $r\lesssim 5 $ cMpc/h followed by an excess of transmission at $r\gtrsim 10$ cMpc/h detected at $2.7\sigma$. This is in contrast to equivalent measurements at lower redshifts where only the opaque trough is detected. We interpret this excess as a statistical enhancement of the local photoionisation rate due to clustered faint galaxies around the \cfour \,absorbers. Using the analytical framework described in Paper I of this series, we derive a constraint on the average product of the escape fraction and the Lyman continuum photon production efficiency of the galaxy population clustered around the \cfour \,absorbers, $\log \langle f_{\text{esc}}\xi_{\text{ion}}\rangle /[{\rm erg^{-1}~Hz}] = 25.01^{+0.30}_{-0.19}$. This implies that faint galaxies beyond the reach of current facilities may have harder radiation fields and/or larger escape fractions than currently detected objects at the end of the reionisation epoch.

Bayesian statistical inference has become increasingly important for the analysis of observations from the Advanced LIGO and Advanced Virgo gravitational-wave detectors. To this end, iterative simulation techniques, in particular nested sampling and parallel tempering, have been implemented in the software library LALInference to sample from the posterior distribution of waveform parameters of compact binary coalescence events. Nested sampling was mainly developed to calculate the marginal likelihood of a model but can produce posterior samples as a by-product. Thermodynamic integration is employed to calculate the evidence using samples generated by parallel tempering but has been found to be computationally demanding. Here we propose the stepping-stone sampling algorithm, originally proposed by Xie et al. (2011) in phylogenetics and a special case of path sampling, as an alternative to thermodynamic integration. The stepping-stone sampling algorithm is also based on samples from the power posteriors of parallel tempering but has superior performance as fewer temperature steps and thus computational resources are needed to achieve the same accuracy. We demonstrate its performance and computational costs in comparison to thermodynamic integration and nested sampling in a simulation study and a case study of computing the marginal likelihood of a binary black hole signal model applied to simulated data from the Advanced LIGO and Advanced Virgo gravitational wave detectors. To deal with the inadequate methods currently employed to estimate the standard errors of evidence estimates based on power posterior techniques, we propose a novel block bootstrap approach and show its potential in our simulation study and LIGO application.

We have recently commissioned a novel infrared ($0.9-1.7$ $\mu$m) integral field spectrograph (IFS) called the Wide Integral Field Infrared Spectrograph (WIFIS). WIFIS is a unique instrument that offers a very large field-of-view (50$^{\prime\prime}$ x 20$^{\prime\prime}$) on the 2.3-meter Bok telescope at Kitt Peak, USA for seeing-limited observations at moderate spectral resolving power. The measured spatial sampling scale is $\sim1\times1^{\prime\prime}$ and its spectral resolving power is $R\sim2,500$ and $3,000$ in the $zJ$ ($0.9-1.35$ $\mu$m) and $H_{short}$ ($1.5-1.7$ $\mu$m) modes, respectively. WIFIS's corresponding etendue is larger than existing near-infrared (NIR) IFSes, which are mostly designed to work with adaptive optics systems and therefore have very narrow fields. For this reason, this instrument is specifically suited for studying very extended objects in the near-infrared such as supernovae remnants, galactic star forming regions, and nearby galaxies, which are not easily accessible by other NIR IFSes. This enables scientific programs that were not originally possible, such as detailed surveys of a large number of nearby galaxies or a full accounting of nucleosynthetic yields of Milky Way supernova remnants. WIFIS is also designed to be easily adaptable to be used with larger telescopes. In this paper, we report on the overall performance characteristics of the instrument, which were measured during our commissioning runs in the second half of 2017. We present measurements of spectral resolving power, image quality, instrumental background, and overall efficiency and sensitivity of WIFIS and compare them with our design expectations. Finally, we present a few example observations that demonstrate WIFIS's full capability to carry out infrared imaging spectroscopy of extended objects, which is enabled by our custom data reduction pipeline.

We introduce a new method for determining the influence of galaxies and active galactic nuclei (AGN) on the physical state of the intergalactic medium (IGM) at high redshift and illustrate its potential via a first application to the field of the $z=6.42$ QSO J1148+5251. By correlating the spatial positions of spectroscopically-confirmed Lyman break galaxies (LBGs) with fluctuations in the Lyman alpha forest seen in the high signal-to-noise spectrum of a background QSO, we provide a statistical measure of the typical escape fraction of Lyman continuum photons close to the end of cosmic reionisation. Here we use Keck DEIMOS spectroscopy to locate 7 colour-selected LBGs in the redshift range $5.3\lesssim z\lesssim 6.4$ and confirm a faint $z=5.701$ AGN. We then examine the spatial correlation between this sample and Ly$\alpha$/Ly$\beta$ transmission fluctuations in a Keck ESI spectrum of the QSO. Interpreting the statistical HI proximity effect as arising from faint galaxies clustered around the detected LBGs, we translate the observed mean Ly$\alpha$ transmitted flux around an average detected LBG into a constraint on the mean escape fraction $\langle f_{\rm esc}\rangle\geq0.08$ at $z\simeq6$. We also report evidence of the individual transverse HI proximity effect of a $z=6.177$ luminous LBG via a Ly$\beta$ transmission spike and two broad Ly$\alpha$ transmission spikes around the $z=5.701$ AGN. We discuss the possible origin of such associations which suggest that while faint galaxies are primarily driving reionisation, luminous galaxies and AGN may provide important contributions to the UV background or thermal fluctuations of the IGM at $z\simeq6$. Although a limited sample, our results demonstrate the potential of making progress using this method in resolving one of the most challenging aspects of the contribution of galaxies and AGN to cosmic reionisation.

We describe an update to the Herschel-SPIRE Fourier-Transform Spectrometer (FTS) calibration for extended sources, which incorporates a correction for the frequency-dependent far-field feedhorn efficiency, $\eta_\mathrm{FF}$. This significant correction affects all FTS extended-source calibrated spectra in sparse or mapping mode, regardless of the spectral resolution. Line fluxes and continuum levels are underestimated by factors of 1.3-2 in the Spectrometer Long-Wavelength band (SLW, 447-1018 GHz; 671-294 $\mu$m) and 1.4-1.5 in the Spectrometer Short-Wavelength band (SSW, 944-1568 GHz; 318-191 $\mu$m). The correction was implemented in the FTS pipeline version 14.1 and has also been described in the SPIRE Handbook since Feb 2017. Studies based on extended-source calibrated spectra produced prior to this pipeline version should be critically reconsidered using the current products available in the Herschel Science Archive. Once the extended-source calibrated spectra are corrected for $\eta_\mathrm{FF}$, the synthetic photometry and the broadband intensities from SPIRE photometer maps agree within 2-4% -- similar levels to the comparison of point-source calibrated spectra and photometry from point-source calibrated maps. The two calibration schemes for the FTS are now self-consistent: the conversion between the corrected extended-source and point-source calibrated spectra can be achieved with the beam solid angle and a gain correction that accounts for the diffraction loss.

It is well known that stellar companions can influence the evolution of a protoplanetary disk. Nevertheless, previous disk surveys did not - and could not - consistently exclude binaries from their samples. We present a study dedicated to investigating the frequency of ongoing disk accretion around single stars in a star-forming region. We obtained near-infrared spectroscopy of 54 low-mass stars selected from a high-angular resolution survey in the 2-3 Myr-old Chamaeleon I region to determine the presence of Brackett-$\gamma$ emission, taking the residual chance of undetected multiplicity into account, which we estimate to be on the order of 30%. The result is compared with previous surveys of the same feature in binary stars of the same region to provide a robust estimate of the difference between the accretor fractions of single stars and individual components of binary systems. We find Br$\gamma$ emission among $39.5^{+14.0}_{-9.9}$% of single stars, which is a significantly higher fraction than for binary stars in Chamaeleon I. In particular, close binary systems with separations <100 AU show emission in only $6.5^{+16.5}_{-3.0}$% of the cases according to the same analysis. The emitter frequency of wider binaries appears consistent with the single star value. Interpreting Br$\gamma$ emission as a sign of ongoing accretion and correcting for sensitivity bias, we infer an accretor fraction of single stars of F_acc=$47.8^{+14.0}_{-9.9}$%. This is slightly higher but consistent with previous estimates that do not clearly exclude binaries from their samples. Through our robust and consistent analysis, we confirm that the fraction of young single stars harboring accretion disks is much larger than that of close binaries at the same age. Our findings have important implications for the timescales of disk evolution and planet formation.

Using the latest numerical simulations of rotating stellar core collapse, we present a Bayesian framework to extract the physical information encoded in noisy gravitational wave signals. We fit Bayesian principal component regression models with known and unknown signal arrival times to reconstruct gravitational wave signals, and subsequently fit known astrophysical parameters on the posterior means of the principal component coefficients using a linear model. We predict the ratio of rotational kinetic energy to gravitational energy of the inner core at bounce by sampling from the posterior predictive distribution, and find that these predictions are generally very close to the true parameter values, with $90\%$ credible intervals $\sim 0.04$ and $\sim 0.06$ wide for the known and unknown arrival time models respectively. Two supervised machine learning methods are implemented to classify precollapse differential rotation, and we find that these methods discriminate rapidly rotating progenitors particularly well. We also introduce a constrained optimization approach to model selection to find an optimal number of principal components in the signal reconstruction step. Using this approach, we select 14 principal components as the most parsimonious model.

Extreme mass ratio inspirals (EMRIs) are thought to be one of the most exciting gravitational wave sources to be detected with LISA. Due to their complicated nature and weak amplitudes the detection and parameter estimation of such sources is a challenging task. In this paper we present a statistical methodology based on Bayesian inference in which the estimation of parameters is carried out by advanced Markov chain Monte Carlo (MCMC) algorithms such as parallel tempering MCMC. We analysed high and medium mass EMRI systems that fall well inside the low frequency range of LISA. In the context of the Mock LISA Data Challenges, our investigation and results are also the first instance in which a fully Markovian algorithm is applied for EMRI searches. Results show that our algorithm worked well in recovering EMRI signals from different (simulated) LISA data sets having single and multiple EMRI sources and holds great promise for posterior computation under more realistic conditions. The search and estimation methods presented in this paper are general in their nature, and can be applied in any other scenario such as AdLIGO, AdVIRGO and Einstein Telescope with their respective response functions.

We review a construction of holographic geometries dual to N=4 SYM theory on a Friedmann-Robertson-Walker background and in the presence or absence of a gluon condensate and instanton density. We find the most general solution with arbitrary scale factor and show that it is diffeomorphic to topological black holes. We introduce a time-dependent boundary cosmological constant \lambda(t) and show energy-momentum conservation in this background. For constant \lambda, the deconfinement properties of the temporal Wilson loop are analysed. In most cases the Wilson loop confines throughout cosmological evolution. However, there is an exceptional case which shows a transition from deconfinement at early times to confinement at late times. We classify the presence or absence of horizons, with important implications for the Wilson loop.

For type IIB supergravity with a running axio-dilaton, we construct bulk solutions which admit a cosmological background metric of Friedmann-Robertson-Walker type. These solutions include both a dark radiation term in the bulk as well as a four-dimensional (boundary) cosmological constant, while gravity at the boundary remains non-dynamical. We holographically calculate the stress-energy tensor, showing that it consists of two contributions: The first one, generated by the dark radiation term, leads to the thermal fluid of N = 4 SYM theory, while the second, the conformal anomaly, originates from the boundary cosmological constant. Conservation of the boundary stress tensor implies that the boundary cosmological constant is time-independent, such that there is no exchange between the two stress-tensor contributions. We then study (de)confinement by evaluating the Wilson loop in these backgrounds. While the dark radiation term favours deconfinement, a negative cosmological constant drives the system into a confined phase. When both contributions are present, we find an oscillating universe with negative cosmological constant which undergoes periodic (de)confinement transitions as the scale of three space expands and re-contracts.

Presented in this paper is a technique that we propose for extracting the physical parameters of a rotating stellar core collapse from the observation of the associated gravitational wave signal from the collapse and core bounce. Data from interferometric gravitational wave detectors can be used to provide information on the mass of the progenitor model, precollapse rotation and the nuclear equation of state. We use waveform libraries provided by the latest numerical simulations of rotating stellar core collapse models in general relativity, and from them create an orthogonal set of eigenvectors using principal component analysis. Bayesian inference techniques are then used to reconstruct the associated gravitational wave signal that is assumed to be detected by an interferometric detector. Posterior probability distribution functions are derived for the amplitudes of the principal component analysis eigenvectors, and the pulse arrival time. We show how the reconstructed signal and the principal component analysis eigenvector amplitude estimates may provide information on the physical parameters associated with the core collapse event.

Inspiral signals from binary compact objects (black holes and neutron stars) are primary targets of the ongoing searches by ground-based gravitational-wave interferometers (LIGO, Virgo, GEO-600 and TAMA-300). We present parameter-estimation simulations for inspirals of black-hole--neutron-star binaries using Markov-chain Monte-Carlo methods. For the first time, we have both estimated the parameters of a binary inspiral source with a spinning component and determined the accuracy of the parameter estimation, for simulated observations with ground-based gravitational-wave detectors. We demonstrate that we can obtain the distance, sky position, and binary orientation at a higher accuracy than previously suggested in the literature. For an observation of an inspiral with sufficient spin and two or three detectors we find an accuracy in the determination of the sky position of typically a few tens of square degrees.

We present a Markov-chain Monte-Carlo (MCMC) technique to study the source parameters of gravitational-wave signals from the inspirals of stellar-mass compact binaries detected with ground-based gravitational-wave detectors such as LIGO and Virgo, for the case where spin is present in the more massive compact object in the binary. We discuss aspects of the MCMC algorithm that allow us to sample the parameter space in an efficient way. We show sample runs that illustrate the possibilities of our MCMC code and the difficulties that we encounter.

We report on the analysis of selected single source data sets from the first round of the Mock LISA Data Challenges (MLDC) for white dwarf binaries. We implemented an end-to-end pipeline consisting of a grid-based coherent pre-processing unit for signal detection, and an automatic Markov Chain Monte Carlo post-processing unit for signal evaluation. We demonstrate that signal detection with our coherent approach is secure and accurate, and is increased in accuracy and supplemented with additional information on the signal parameters by our Markov Chain Monte Carlo approach. We also demonstrate that the Markov Chain Monte Carlo routine is additionally able to determine accurately the noise level in the frequency window of interest.

We present Keck near-infrared and WIYN optical photometry of a sample of galaxies detected by near-infrared narrowband imaging in the fields of quasar metal absorption line systems at z~2.4. Wide separations (0.6-1.6/h Mpc) from the quasars indicates that they are not directly responsible for the absorption systems. From the color excess of the galaxies we derived line fluxes, star formation rates, and equivalent widths. The data are consistent with one source having an active nucleus and two sources containing regions of star formation. The blue (R-K) colors for the sources suggest relatively lesser dust content. We discuss possible projects using current wide-field infrared instruments, which can cover an order of magnitude greater area with modest allocations of telescope time.

We present a strategy for a statistically rigorous Bayesian approach to the problem of determining cosmological parameters from the results of observations of anisotropies in the cosmic microwave background. Our strategy relies on Markov chain Monte Carlo methods, specifically the Metropolis-Hastings algorithm, to perform the necessary high-dimensional integrals. We describe the Metropolis-Hastings algorithm in detail and discuss the results of our test on simulated data.

We present a strategy for a statistically rigorous Bayesian approach to the problem of determining cosmological parameters from the results of observations of anisotropies in the cosmic microwave radiation background. We propose the application of Markov chain Monte Carlo methods, specifically the Metropolis-Hastings algorithm, to estimate the parameters. A complete statistical analysis is presented, with the Metropolis-Hastings algorithm described in detail.