We present early results from a study addressing the question of how one treats the propagation of incertitude, that is, epistemic uncertainty, in input parameters in astrophysical simulations. As an example, we look at the propagation of incertitude in control parameters for stellar winds in MESA stellar evolution simulations. We apply two methods of incertitude propagation, the Cauchy Deviates method and the Quadratic Response Surface method, to quantify the output uncertainty in the final white dwarf mass given a range of values for wind parameters. The methodology we apply is applicable to the problem of propagating input incertitudes through any simulation code treated as a "black box," i.e. a code for which the algorithmic details are either inaccessible or prohibitively complicated. We have made the tools developed for this study freely available to the community.

We present a new, fully generative model of optical telescope image sets, along with a variational procedure for inference. Each pixel intensity is treated as a Poisson random variable, with a rate parameter dependent on latent properties of stars and galaxies. Key latent properties are themselves random, with scientific prior distributions constructed from large ancillary data sets. We check our approach on synthetic images. We also run it on images from a major sky survey, where it exceeds the performance of the current state-of-the-art method for locating celestial bodies and measuring their colors.

We reanalyze microlensing events in the published list of anomalous events that were observed from the OGLE lensing survey conducted during 2004-2008 period. In order to check the existence of possible degenerate solutions and extract extra information, we conduct analyses based on combined data from other survey and follow-up observation and consider higher-order effects. Among the analyzed events, we present analyses of 8 events for which either new solutions are identified or additional information is obtained. We find that the previous binary-source interpretations of 5 events are better interpreted by binary-lens models. These events include OGLE-2006-BLG-238, OGLE-2007-BLG-159, OGLE-2007-BLG-491, OGLE-2008-BLG-143, and OGLE-2008-BLG-210. With additional data covering caustic crossings, we detect finite-source effects for 6 events including OGLE-2006-BLG-215, OGLE-2006-BLG-238, OGLE-2006-BLG-450, OGLE-2008-BLG-143, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. Among them, we are able to measure the Einstein radii of 3 events for which multi-band data are available. These events are OGLE-2006-BLG-238, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. For OGLE-2008-BLG-143, we detect higher-order effect induced by the changes of the observer's position caused by the orbital motion of the Earth around the Sun. In addition, we present degenerate solutions resulting from the known close/wide or ecliptic degeneracy. Finally, we note that the masses of the binary companions of the lenses of OGLE-2006-BLG-450 and OGLE-2008-BLG-210 are in the brown-dwarf regime.

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.

The mass function and statistics of binaries provide important diagnostics of the star formation process. Despite this importance, the mass function at low masses remains poorly known due to observational difficulties caused by the faintness of the objects. Here we report the microlensing discovery and characterization of a binary lens composed of very low-mass stars just above the hydrogen-burning limit. From the combined measurements of the Einstein radius and microlens parallax, we measure the masses of the binary components of $0.10\pm 0.01\ M_\odot$ and $0.09\pm 0.01\ M_\odot$. This discovery demonstrates that microlensing will provide a method to measure the mass function of all Galactic populations of very low mass binaries that is independent of the biases caused by the luminosity of the population.

Microlensing can provide a useful tool to probe binary distributions down to low-mass limits of binary companions. In this paper, we analyze the light curves of 8 binary lensing events detected through the channel of high-magnification events during the seasons from 2007 to 2010. The perturbations, which are confined near the peak of the light curves, can be easily distinguished from the central perturbations caused by planets. However, the degeneracy between close and wide binary solutions cannot be resolved with a $3\sigma$ confidence level for 3 events, implying that the degeneracy would be an important obstacle in studying binary distributions. The dependence of the degeneracy on the lensing parameters is consistent with a theoretic prediction that the degeneracy becomes severe as the binary separation and the mass ratio deviate from the values of resonant caustics. The measured mass ratio of the event OGLE-2008-BLG-510/MOA-2008-BLG-369 is $q\sim 0.1$, making the companion of the lens a strong brown-dwarf candidate.

We present the analysis result of a gravitational binary-lensing event OGLE-2005-BLG-018. The light curve of the event is characterized by 2 adjacent strong features and a single weak feature separated from the strong features. The light curve exhibits noticeable deviations from the best-fit model based on standard binary parameters. To explain the deviation, we test models including various higher-order effects of the motions of the observer, source, and lens. From this, we find that it is necessary to account for the orbital motion of the lens in describing the light curve. From modeling of the light curve considering the parallax effect and Keplerian orbital motion, we are able to measure not only the physical parameters but also a complete orbital solution of the lens system. It is found that the event was produced by a binary lens located in the Galactic bulge with a distance $6.7\pm 0.3$ kpc from the Earth. The individual lens components with masses $0.9\pm 0.3\ M_\odot$ and $0.5\pm 0.1\ M_\odot$ are separated with a semi-major axis of $a=2.5 \pm 1.0$ AU and orbiting each other with a period $P=3.1 \pm 1.3$ yr. The event demonstrates that it is possible to extract detailed information about binary lens systems from well-resolved lensing light curves.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

To improve the planet detection efficiency, current planetary microlensing experiments are focused on high-magnification events searching for planetary signals near the peak of lensing light curves. However, it is known that central perturbations can also be produced by binary companions and thus it is important to distinguish planetary signals from those induced by binary companions. In this paper, we analyze the light curves of microlensing events OGLE-2007-BLG-137/MOA-2007-BLG-091, OGLE-2007-BLG-355/MOA-2007-BLG-278, and MOA-2007-BLG-199/OGLE-2007-BLG-419, for all of which exhibit short-term perturbations near the peaks of the light curves. From detailed modeling of the light curves, we find that the perturbations of the events are caused by binary companions rather than planets. From close examination of the light curves combined with the underlying physical geometry of the lens system obtained from modeling, we find that the short time-scale caustic-crossing feature occurring at a low or a moderate base magnification with an additional secondary perturbation is a typical feature of binary-lens events and thus can be used for the discrimination between the binary and planetary interpretations.

We analyze OGLE-2007-BLG-050, a high magnification microlensing event (A ~ 432) whose peak occurred on 2 May, 2007, with pronounced finite-source and parallax effects. We compute planet detection efficiencies for this event in order to determine its sensitivity to the presence of planets around the lens star. Both finite-source and parallax effects permit a measurement of the angular Einstein radius \theta_E = 0.48 +/- 0.01 mas and the parallax \pi_E = 0.12 +/- 0.03, leading to an estimate of the lens mass M = 0.50 +/- 0.14 M_Sun and its distance to the observer D_L = 5.5 +/- 0.4 kpc. This is only the second determination of a reasonably precise (<30%) mass estimate for an isolated unseen object, using any method. This allows us to calculate the planetary detection efficiency in physical units (r_\perp, m_p), where r_\perp is the projected planet-star separation and m_p is the planet mass. When computing planet detection efficiency, we did not find any planetary signature and our detection efficiency results reveal significant sensitivity to Neptune-mass planets, and to a lesser extent Earth-mass planets in some configurations. Indeed, Jupiter and Neptune-mass planets are excluded with a high confidence for a large projected separation range between the planet and the lens star, respectively [0.6 - 10] and [1.4 - 4] AU, and Earth-mass planets are excluded with a 10% confidence in the lensing zone, i.e. [1.8 - 3.1] AU.

We combine all available information to constrain the nature of OGLE-2005-BLG-071Lb, the second planet discovered by microlensing and the first in a high-magnification event. These include photometric and astrometric measurements from Hubble Space Telescope, as well as constraints from higher order effects extracted from the ground-based light curve, such as microlens parallax, planetary orbital motion and finite-source effects. Our primary analysis leads to the conclusion that the host of Jovian planet OGLE-2005-BLG-071Lb is an M dwarf in the foreground disk with mass M= 0.46 +/- 0.04 Msun, distance D_l = 3.3 +/- 0.4 kpc, and thick-disk kinematics v_LSR ~ 103 km/s. From the best-fit model, the planet has mass M_p = 3.8 +/- 0.4 M_Jup, lies at a projected separation r_perp = 3.6 +/- 0.2 AU from its host and so has an equilibrium temperature of T ~ 55 K, i.e., similar to Neptune. A degenerate model less favored by \Delta\chi^2 = 2.1 (or 2.2, depending on the sign of the impact parameter) gives similar planetary mass M_p = 3.4 +/- 0.4 M_Jup with a smaller projected separation, r_\perp = 2.1 +/- 0.1 AU, and higher equilibrium temperature T ~ 71 K. These results from the primary analysis suggest that OGLE-2005-BLG-071Lb is likely to be the most massive planet yet discovered that is hosted by an M dwarf. However, the formation of such high-mass planetary companions in the outer regions of M-dwarf planetary systems is predicted to be unlikely within the core-accretion scenario. There are a number of caveats to this primary analysis, which assumes (based on real but limited evidence) that the unlensed light coincident with the source is actually due to the lens, that is, the planetary host. However, these caveats could mostly be resolved by a single astrometric measurement a few years after the event.

We present an analysis of OGLE 2004-BLG-254, a high-magnification and relatively short duration microlensing event in which the source star, a Bulge K3-giant, has been spatially resolved by a point-like lens. We have obtained dense photometric coverage of the event light curve with OGLE and PLANET telescopes, as well as a high signal-to-noise ratio spectrum taken while the source was still magnified by 20, using the UVES/VLT spectrograph. Our dense coverage of this event allows us to measure limb darkening of the source star in the I and R bands. We also compare previous measurements of linear limb-darkening coefficients involving GK-giant stars with predictions from ATLAS atmosphere models. We discuss the case of K-giants and find a disagreement between limb-darkening measurements and model predictions, which may be caused by the inadequacy of the linear limb-darkening law.

In the favoured core-accretion model of formation of planetary systems, solid planetesimals accumulate to build up planetary cores, which then accrete nebular gas if they are sufficiently massive. Around M-dwarf stars (the most common stars in our Galaxy), this model favours the formation of Earth-mass to Neptune-mass planets with orbital radii of 1 to 10 astronomical units (AU), which is consistent with the small number of gas giant planets known to orbit M-dwarf host stars. More than 170 extrasolar planets have been discovered with a wide range of masses and orbital periods, but planets of Neptune's mass or less have not hitherto been detected at separations of more than 0.15 AU from normal stars. Here we report the discovery of a 5.5 (+5.5/-2.7) M_earth planetary companion at a separation of 2.6 (+1.5/-0.6) AU from a 0.22 (+0.21/-0.11) M_solar M-dwarf star. (We propose to name it OGLE-2005-BLG-390Lb, indicating a planetary mass companion to the lens star of the microlensing event.) The mass is lower than that of GJ876d, although the error bars overlap. Our detection suggests that such cool, sub-Neptune-mass planets may be more common than gas giant planets, as predicted by the core accretion theory.

It is generally assumed that the two dark components of the energy density of the universe, a smooth component called dark energy and a fluid of nonrelativistic weakly interacting particles called dark matter, are independent of each other and interact only through gravity. In this paper, we consider a class of models in which the dark matter and dark energy interact directly. The dark matter particle mass is proportional to the value of a scalar field, and the energy density of this scalar field comprises the dark energy. We study the phenomenology of these models and calculate the luminosity distance as a function of redshift and the CMB anisotropy spectrum for several cases. We find that the phenomenology of these models can differ significantly from the standard case, and current observations can already rule out the simplest models.

Models of dark energy are conveniently characterized by the equation-of-state parameter w=p/\rho, where \rho is the energy density and p is the pressure. Imposing the Dominant Energy Condition, which guarantees stability of the theory, implies that w \geq -1. Nevertheless, it is conceivable that a well-defined model could (perhaps temporarily) have w<-1, and indeed such models have been proposed. We study the stability of dynamical models exhibiting w<-1 by virtue of a negative kinetic term. Although naively unstable, we explore the possibility that these models might be phenomenologically viable if thought of as effective field theories valid only up to a certain momentum cutoff. Under our most optimistic assumptions, we argue that the instability timescale can be greater than the age of the universe, but only if the cutoff is at or below 100 MeV. We conclude that it is difficult, although not necessarily impossible, to construct viable models of dark energy with w<-1; observers should keep an open mind, but the burden is on theorists to demonstrate that any proposed new models are not ruled out by rapid vacuum decay.

The observables $T/S$ and $n-1$ are key to testing and understanding inflation. ($T$, $S$, and $n-1$ respectively quantify the gravity-wave and density-perturbation contributions to CMB anisotropy and the deviation of the density perturbations from the scale-invariant form.) Absent a standard model, there is no definite prediction for, or relation between, $T/S$ and $n-1$. By reformulating the equations governing inflation we show that models generally predict $T/S \approx -5(n-1)$ or 0, and in particular, if $n>0.85$, $T/S$ is expected to be $>10^{-3}$.