Cosmological probes pose an inverse problem where the measurement result is obtained through observations, and the objective is to infer values of model parameters which characterize the underlying physical system -- our Universe. Modern cosmological probes increasingly rely on measurements of the small-scale structure, and the only way to accurately model physical behavior on those scales, roughly 65 Mpc/h or smaller, is via expensive numerical simulations. In this paper, we provide a detailed description of a novel statistical framework for obtaining accurate parameter constraints by combining observations with a very limited number of cosmological simulations. The proposed framework utilizes multi-output Gaussian process emulators that are adaptively constructed using Bayesian optimization methods. We compare several approaches for constructing multi-output emulators that enable us to take possible inter-output correlations into account while maintaining the efficiency needed for inference. Using Lyman alpha forest flux power spectrum, we demonstrate that our adaptive approach requires considerably fewer --- by a factor of a few in Lyman alpha P(k) case considered here --- simulations compared to the emulation based on Latin hypercube sampling, and that the method is more robust in reconstructing parameters and their Bayesian credible intervals.

We report the discovery of a new quasar: SDSS J022155.26-064916.6. This object was discovered while reducing spectra of a sample of stars being considered as spectrophotometric standards for the Dark Energy Survey. The flux and wavelength calibrated spectrum is presented with four spectral lines identified. From these lines, the redshift is determined to be z is approximately equal to 0.806. In addition, the rest-frame u-, g-, and r-band luminosity, determined using a k-correction obtained with synthetic photometry of a proxy QSO, are reported as 7.496 $\times 10^{13}$ solar luminosities, 2.049 $\times 10^{13}$ solar luminosities, and $1.896 \times 10^{13}$ solar luminosities, respectively.

We report the observation and physical characterization of the possible dwarf planet \UZ\ ("DeeDee"), a dynamically detached trans-Neptunian object discovered at 92 AU. This object is currently the second-most distant known trans-Neptunian object with reported orbital elements, surpassed in distance only by the dwarf planet Eris. The object was discovered with an $r$-band magnitude of 23.0 in data collected by the Dark Energy Survey between 2014 and 2016. Its 1140-year orbit has $(a,e,i) = (109~\mathrm{AU}, 0.65, 26.8^{\circ})$. It will reach its perihelion distance of 38 AU in the year 2142. Integrations of its orbit show it to be dynamically stable on Gyr timescales, with only weak interactions with Neptune. We have performed followup observations with ALMA, using 3 hours of on-source integration time to measure the object's thermal emission in the Rayleigh-Jeans tail. The signal is detected at 7$\sigma$ significance, from which we determine a $V$-band albedo of $13.1^{+3.3}_{-2.4}\mathrm{(stat)}^{+2.0}_{-1.4}\mathrm{(sys)}$ percent and a diameter of $635^{+57}_{-61}\mathrm{(stat)}^{+32}_{-39}\mathrm{(sys)}$~km, assuming a spherical body with uniform surface properties.

Finding the optimal ordering of k-subsets with respect to an objective function is known to be an extremely challenging problem. In this paper we introduce a new objective for this task, rooted in the problem of star identification on spacecrafts: subsets of detected spikes are to be generated in an ordering that minimizes time to detection of a valid star constellation. We carry out an extensive analysis of the combinatorial optimization problem, and propose multiple algorithmic solutions, offering different quality-complexity trade-offs. Three main approaches are investigated: exhaustive search (branch and prune), goal-driven (greedy scene elimination, minimally intersecting subsets), and stateless algorithms which implicitly seek to satisfy the problem's goals (pattern shifting, base unrank). In practical terms, these last algorithms are found to provide satisfactory approximations to the ideal performance levels, at small computational costs.

We analyze the impact of the installation of a lunar laser ranging device on the NTT 3.6m ESO telescope. With such an instrument, the scientific communities of fundamental physics and solar system formation will highly benefit of the only LLR station in the Southern Hemisphere. The quality of the NTT 3.6 meter telescope will also greatly complement the LLR 3.5 meter Apache Point telescope (3.5 m) instrument in the Northern Hemisphere (USA) which is the best instrument for tracking the Moon since 2006. Finally, we also consider the technical characteristics of such installation including the observational constraints.

The second 2+ state of 12C, predicted over fifty years ago as an excitation of the Hoyle state, has been unambiguously identified using the 12C(g,a_0)8Be reaction. The alpha particles produced by the photodisintegration of 12C were detected using an Optical Time Projection Chamber (O-TPC). Data were collected at beam energies between 9.1 and 10.7 MeV using the intense nearly mono-energetic gamma-ray beams at the HIgS facility. The measured angular distributions determine the cross section and the E1-E2 relative phases as a function of energy leading to an unambiguous identification of the second 2+ state in 12C at 10.03(11) MeV, with a total width of 800(130) keV and a ground state gamma-decay width of 60(10) meV; B(E2: 2+ ---> gs) = 0.73(13) e2fm4 [or 0.45(8) W.u.]. The Hoyle state and its rotational 2+ state that are more extended than the ground state of 12C presents a challenge and constraints for models attempting to reveal the nature of three alpha particle states in 12C. Specifically it challenges the ab-initio Lattice Effective Field Theory (L-EFT) calculations that predict similar r.m.s. radii for the ground state and the Hoyle state.

It has been argued that there is biological and modeling evidence that a non-linear diffusion coefficient of the type D(b) = D_0 b^{k} underlies the formation of a number of growth patterns of bacterial colonies. We study a reaction-diffusion system with a non-linear diffusion coefficient introduced by Ben-Jacob et al. Due to the fact that the bacterial diffusion coefficient vanishes when the bacterial density b -> 0, the standard linear stability analysis for fronts cannot be used. We introduce an extension of the stability analysis which can be applied to such singular fronts, map out the region of stability in the D-k-plane and derive an interfacial approximation in some limits. Our linear stability analysis and sharp interface formulation will also be applicable to other examples of interface formation due to nonlinear diffusion, like in porous media or in the problem of vortex motion in superconductors.