Search

Interstellar neutral (ISN) hydrogen is the most abundant species in the outer heliosheath and the very local interstellar medium (VLISM). Charge exchange collisions in the outer heliosheath result in filtration, reducing the ISN hydrogen density inside the heliosphere. Additionally, these atoms are intensively ionized close to the Sun, resulting in a substantial reduction of their density within a few au from the Sun. The products of this ionization - pickup ions (PUIs) - are detected by charged particle detectors. The Solar Wind Around Pluto (SWAP) instrument on New Horizons provides, for the first time, PUI observations from the distant heliosphere. We analyze the observations collected between 22 and 52 au from the Sun to find the ISN hydrogen density profile and compare the results with predictions from global heliosphere models. We conclude that the density profile derived from the observations is inconsistent with steady-state model predictions. This discrepancy is not explained by time variations close to the Sun and thus may be related to the temporal evolution of the outer boundaries or VLISM conditions. Furthermore, we show that the cold and hot models of ISN hydrogen distribution are not a good approximation closer to the termination shock. Therefore, we recommend a new fiduciary point based on the available New Horizons observations at 40 au from the Sun, at ecliptic direction (285.62{\deg}, 1.94{\deg}), where the ISN hydrogen density is 0.11 cm$^{-3}$. The continued operation of New Horizons should give better insight into the source of the discussed discrepancy., Comment: Author Accepted Manuscript. Accepted for publication in The Astrophysical Journal Letters. 11 pages, 4 figures, 1 table

We present a study of the influence of solar UV anisotropy on the heliospheric backscatter helioglow generated by resonant scattering of solar Lyman-alpha photons on interstellar hydrogen atoms around the Sun. Simulations based on the WawHelioGlow model suggest that the response of the helioglow pole-to-ecliptic ratio to the anisotropy is linear, but 15% of the anisotropy (polar darkening) generates 30-40% change in the ratio in the solar minimum and 15-20% in the solar maximum. We attribute this difference to an interplay between the solar UV anisotropy and the latitudinal structure of the solar wind in solar minima. The solar UV anisotropy also increases the helioglow intensity from the downwind direction by 5-10%, due to the influence of the anisotropy on the ionization losses and trajectories of atoms passing by the Sun in polar regions. Consequently, mid-latitude regions (in the heliographic and ecliptic coordinates) are least affected by the UV anisotropy. By comparison of the simulation results with observations of the SOHO/SWAN satellite instrument, we derive the day-by-day time evolution of the solar Lyman-alpha anisotropy for the north and south poles over two solar cycles from 1996 to 2022. The inferred anisotropy is ~5-10% in solar minima and ~15-25% in solar maxima, the northern anisotropy being stronger than the southern. Our study suggests that in solar minima a highly structured solar wind is associated with relatively small solar UV anisotropy, while in solar maxima the solar wind is more isotropic but a substantial solar UV anisotropy appears., Comment: 25 pages, 14 figures. Accepted for publication in Astrophysical Journal

Direct-sampling observations of interstellar neutral (ISN) species and their secondary populations give information about the physical state of the local interstellar medium and processes occuring in the outer heliosheath. Such observations are performed from Earth's orbit by the IBEX-Lo experiment on board the Interstellar Boundary Explorer (IBEX) mission. IBEX ISN viewing is restricted to directions close to perpendicular to the Earth-Sun line, which limits the observations of interstellar species to several months during the year. A greatly improved data set will be possible for the upcoming Interstellar Mapping and Acceleration Probe (IMAP) mission due to a novel concept of putting the IMAP ISN detector, IMAP-Lo, on a pivot platform that varies the angle of observation relative to the Sun-Earth line and the detector boresight. Here we suggest a 2 yr scenario for varying the viewing angle in such a way that all the necessary atom components can be observed sufficiently well to achieve the science goals of the nominal IMAP mission. This scenario facilitates, among others, removal of the correlation of the inflow parameters of interstellar gas, unambiguous analysis of the primary and secondary populations of interstellar helium (He), neon (Ne) and oxygen (O), and determination of the ionization rates of He and Ne free of possible calibration bias. The scheme is operationally simple, provides a good counting statistics, and synergizes observations of interstellar species and heliospheric energetic neutral atoms., Comment: Approved for ApjS

In situ observations of interstellar neutral (ISN) helium atoms by the IBEX-Lo instrument onboard the Interstellar Boundary Explorer (IBEX) mission are used to determine the velocity and temperature of the pristine very local interstellar medium (VLISM). Most ISN helium atoms penetrating the heliosphere, known as the primary population, originate in the pristine VLISM. As the primary atoms travel through the outer heliosheath, they charge exchange with He$^+$ ions in slowed and compressed plasma creating the secondary population. With more than 2.4 million ISN helium atoms sampled by IBEX during ISN seasons 2009-2020, we compare the observations with predictions of a parametrized model of ISN helium transport in the heliosphere. We account for the filtration of ISN helium atoms at the heliospheric boundaries by charge exchange and elastic collisions. We examine the sensitivity of the ISN helium fluxes to the interstellar conditions described by the pristine VLISM velocity, temperature, magnetic field, and composition. We show that comprehensive modeling of the filtration processes is critical for interpreting ISN helium observations, as the change in the derived VLISM conditions exceeds the statistical uncertainties when accounting for these effects. The pristine VLISM parameters found by this analysis are the flow speed (26.6 km s$^{-1}$), inflow direction in ecliptic coordinates (255.7$^\circ$, 5.04$^\circ$), temperature (7350 K), and B-V plane inclination to the ecliptic plane (53.7$^\circ$). The derived pristine VLISM He$^+$ density is $9.7\times10^3$ cm$^{-3}$. Additionally, we show a strong correlation between the interstellar plasma density and magnetic field strength deduced from these observations., Comment: 13 pages, 3 figures, 2 tables, accepted for publication in ApJ

Interstellar neutral hydrogen flows into the heliosphere as a mixture of the primary and secondary populations from two somewhat different directions due to splitting occurring in the magnetized outer heliosheath. The direction of inflow of interstellar neutral H observed in the inner heliosphere, confronted with that of the unperturbed flow of interstellar neutral helium, is important for understanding the geometry of the distortion of the heliosphere from axial symmetry. It is also needed for facilitating remote-sensing studies of the solar wind structure based on observations of the helioglow, such as those presently performed by SOHO/SWAN, and in a near future by IMAP/GLOWS. In the past, the only means to measure the flow direction of interstellar hydrogen were spectroscopic observations of the helioglow. Here, we propose a new method to determine this parameter based on a long series of photometric observations of the helioglow. The method is based on purely geometric considerations and does not depend on any model and absolute calibration of the measurements. We apply this method to sky maps of the helioglow available from the SOHO/SWAN experiment and derive the mean flow longitude of interstellar hydrogen. We obtain $253.1\degr \pm 2.8\degr$, which is in perfect agreement with the previously obtained results based on spectroscopic observations., Comment: Published in ApJ. Title adjusted to the published version

The IBEX-Lo instrument on the Interstellar Boundary Explorer (IBEX) mission measures interstellar neutral (ISN) helium atoms. The detection of helium atoms is made through negative hydrogen (H$^-$) ions sputtered by the helium atoms from the IBEX-Lo conversion surface. The energy spectrum of ions sputtered by ISN helium atoms is broad and overlaps the four lowest IBEX-Lo electrostatic analyzer (ESA) steps. Consequently, the energy response function for helium atoms does not correspond to the nominal energy step transmission. Moreover, laboratory calibration is incomplete because it is difficult to produce narrow-energy neutral atom beams that are expected for ISN helium atoms. Here, we analyze the ISN helium observations in ESA steps 1-4 to derive the relative in-flight response of IBEX-Lo to helium atoms. We compare the ratios of the observed count rates as a function of the mean ISN helium atom energy estimated using the Warsaw Test Particle Model (WTPM). The WTPM uses a global heliosphere model to calculate charge exchange gains and losses to estimate the secondary ISN helium population. We find that the modeled mean energies of ISN helium atoms, unlike their modeled fluxes, are not very sensitive to the very local interstellar medium parameters. The obtained relative responses supplement the laboratory calibration and enable more detailed quantitative studies of the ISN helium signal. A similar procedure that we applied to the IBEX-Lo observations may be used to complement laboratory calibration of the next-generation IMAP-Lo instrument on the Interstellar Mapping and Acceleration Probe (IMAP) mission., Comment: 17 pages, 8 figures, 1 table, accepted for publication in ApJS

A good understanding of the ionization rates of neutral species in the heliosphere is important for studies of the heliosphere and planetary atmospheres. So far, the intensities of the ionization reactions have been studied based on observations of the contributing phenomena, such as the solar spectral flux in the EUV band and the flux of the solar wind protons, alpha particles, and electrons. The results strongly depend on absolute calibration of these measurements, which, especially for the EUV measurements, is challenging. Here, we propose a novel method of determining the ionization rate of neutral species based on direct sampling of interstellar neutral gas from two locations in space distant to each other. In particular, we suggest performing observations from the vicinity of Earth's orbit and using ratios of fluxes of ISN He for the direct and indirect orbits of interstellar atoms. We identify the most favorable conditions and observations geometries, suitable for implementation on the forthcoming NASA mission Interstellar Mapping and Acceleration Probe., Comment: Accepted for ApJS

We analyze the reasons for the correlation between the temperature, direction, and speed of the interstellar neutral gas inflow into the heliosphere, obtained in analyzes of observations performed by the IBEX-Lo instrument onboard Interstellar Boundary Explorer (IBEX). We point out that this correlation is the combined result of the inability to measure the speed of the atoms that enter the instrument and the restriction of the observations to a short orbital arc around the Sun performed by the instrument during observation. We demonstrate that without the capability to measure the speed, but with the ability to perform observations along longer orbital arcs, or from at least two distant locations on the orbit around the Sun, it is possible to break the parameter correlation. This, however, requires a capability to adjust the boresight of the instrument relative to the spacecraft rotation axis, such as that of the planned IMAP-Lo camera onboard the Interstellar Mapping and Acceleration Probe (IMAP)., Comment: Accepted for ApJ

The IBEX-Lo instrument on board the Interstellar Boundary Explorer (IBEX) mission samples interstellar neutral (ISN) helium atoms penetrating the heliosphere from the very local interstellar medium (VLISM). In this study, we analyze the IBEX-Lo ISN helium observations covering a complete solar cycle, from 2009 through 2020 using a comprehensive uncertainty analysis including statistical and systematic sources.W e employ the Warsaw Test Particle Model to simulate ISN helium fluxes at IBEX, which are subsequently compared with the observed count rate in the three lowest energy steps of IBEX-Lo. The $\chi^2$ analysis shows that the ISN helium flows from ecliptic $(\lambda,\beta)=(255.59^{\circ}\pm0.23^{\circ}, 5.14^{\circ}\pm0.08^{\circ})$, with speed $v_\text{HP}=25.86\pm0.21$ km s$^{-1}$ and temperature $T_\text{HP}=7450\pm140$ K at the heliopause. Accounting for gravitational attraction and elastic collisions, the ISN helium speed and temperature in the pristine VLISM far from the heliopause are $v_\text{VLISM}=25.9$ km s$^{-1}$ and $T_\text{VLISM}=6150$ K, respectively. The time evolution of the ISN helium fluxes at 1 au over 12 years suggests significant changes in the IBEX-Lo detection efficiency, higher ionization rates of ISN helium atoms in the heliosphere than assumed in the model, or an additional unaccounted signal source in the analyzed observations. Nevertheless, we do not find any indication of the evolution of the derived parameters of ISN helium over the period analyzed. Finally, we argue that the continued operation of IBEX-Lo to overlap with the Interstellar Mapping and Acceleration Probe (IMAP) will be pivotal in tracking possible physical changes in the VLISM., Comment: 31 pages, 6 figures, 8 tables, accepted for publication in ApJS

The helioglow is a fluorescence of interstellar atoms inside the heliosphere, where they are excited by the solar EUV. Because the mean free path between collisions for the interstellar gas is comparable to the size of the heliosphere, the distribution function of this gas inside the heliosphere strongly varies in space and with time and is non-Maxwellian. Coupling between realistically modeled solar factors and the distribution function of interstellar neutral gas is accounted for in a helioglow model that we have developed. WawHelioGlow is presented in the accompanying Paper I. Here, we present the evolution of the gas density, solar illumination, helioglow source function, and other relevant parameters building up the helioglow signal for selected lines of sight observed at 1 au. We compare these elements for various phases of the solar cycle and we present the sensitivity of the results to heliolatitudinal anisotropy of the solar EUV output. We assume a realistic latitudinal anisotropy of the solar wind flux using results from analysis of interplanetary scintillations. We compare the simulated helioglow with with selected maps observed by the SOHO/SWAN instrument. We demonstrate that WawHelioGlow is able to reproduce fundamental features of the sky distribution of the helioglow. For some phases of the solar cycle, the model with an anisotropy of the solar EUV output better reproduces the observations, while for other phases no EUV anisotropy is needed. In all simulated cases, the solar wind anisotropy following insight from interplanetary scintillation measurements is present., Comment: Accepted for ApJ

The helioglow is a fluorescence of interstellar atoms inside the heliosphere, where they are excited by the solar EUV emission. So far, the helioglow of interstellar H and He have been detected. The helioglow features a characteristic distribution in the sky, which can be used to derive both the properties of interstellar neutral gas and those of the solar wind. This requires a simulation model capable of catching with a sufficient realism the essential coupling relations between the solar factors and interstellar. The solar factors include the solar wind flux and its variation with time and heliolatitude, as well as the heliolatitude and time variation of the solar EUV output. The ISN gas inside the heliosphere features a complex distribution function, which varies with time and location. The paper presents the first version of a WawHelioGlow simulation model for the helioglow flux using an optically thin, single scattering approximation. The helioglow computations are based on a sophisticated kinetic treatment of the distribution functions of interstellar H and He provided by the (n)WTPM model. The model takes into account heliolatitudinal and spectral variations of the solar EUV output from observations. We present a formulation of the model and the treatment of the solar spectral flux. The accompanying Paper II illustrates details of the line of sight evolution of the elements of the model and a brief comparison of results of the WawHelioGlow code with selected sky maps of the hydrogen helioglow, obtained by the SWAN instrument onboard the SOHO mission., Comment: Accepted for ApJ

In situ measurements of the heliospheric particle populations by the Voyager spacecraft can only be put in an appropriate context with remote-sensing observations of energetic and interstellar neutral atoms (ENAs and ISN, respectively) at 1 au when the time delay between the production and the observation times is taken into account. ENA times of flight from the production regions in the heliosheath are relatively easy to estimate because these atoms follow almost constant speed, force-free trajectories. For the ISN populations, dynamical and ballistic selection effects are important, and times of flight are much longer. We estimate these times for ISN He and H atoms observed by IBEX and in the future by IMAP using the WTPM model with synthesis method. We show that for the primary population atoms, the times of flight are on the order of three solar cycle periods, with a spread equivalent to one solar cycle. For the secondary populations, the times of flight are on the order of ten solar cycle periods, and during the past ten years of observations, IBEX has been collecting secondary He atoms produced in the OHS during almost entire 19th century. ISN atoms penetrating the heliopause at the time of Voyager crossing will become gradually visible about 2027, during the planned IMAP observations. Hypothetical variations in the ISN flow in the Local Interstellar Medium are currently not detectable. Nevertheless, we expect steady-state heliosphere models used with appropriately averaged solar wind parameters to be suitable for understanding the ISN observations., Comment: Accepted for ApJ

Energetic neutral atoms (ENA) are an important tool for investigating the structure of the heliosphere. Recently, it was observed that fluxes of ENAs (with energy $\le$ 55 keV) coming from the upwind and downwind regions of the heliosphere are similar in strength. This led the authors of these observations to hypothesize that the heliosphere is bubble-like rather than comet-like, meaning that it has no extended tail. We investigate the directional distribution of the ENA flux for a wide energy range (3--88 keV) including the observations from IBEX (Interstellar Boundary Explorer), INCA (Ion and Neutral Camera, on board Cassini), and HSTOF (High energy Suprathermal Time Of Flight sensor, on board SOHO, Solar and Heliospheric Observatory). An essential element is the model of pickup ion acceleration at the termination shock (TS) proposed by Zank. We use state of the art models of the global heliosphere, interstellar neutral gas density, and pickup ion distributions. The results, based on the "comet-like" model of the heliosphere, are close in flux magnitude to ENA observations by IBEX, HSTOF and partly by INCA (except for the 5.2-13.5 keV energy channel). We find that the ENA flux from the tail dominates at high energy (in agreement with HSTOF, but not INCA). At low energy, our comet-like model produces the similar strengths of the ENA fluxes from the upwind and downwind directions, which, therefore, removes this as a compelling argument for a bubble-like heliosphere., Comment: Accepted for publication in The Astrophysical Journal

We present candidates for non-pulsating stars lying in the classical Cepheid instability strip based on OGLE photometric maps combined with Str\"omgren photometry obtained with the 4.1-m SOAR telescope, and Gaia DR2 data in four fields in the Large Magellanic Cloud. We selected 19 candidates in total. After analysis of their light curves from OGLE surveys we found that all these stars appear to be photometrically stable at the level of a few mmag. Our results show that non-pulsating stars might constitute to about 21%-30% of the whole sample of giant stars located in the classical instability strip. Furthermore, we identified potential candidates for classical Cepheids with hot companions based on their Str\"omgren colours., Comment: 10 pages, 8 figures, 4 tables, published in MNRAS. Grant number in Acknowledgements updated. Added missing reference to a table in a text

Interstellar neutral gas atoms penetrate the heliopause and reach 1~au, where they are detected by IBEX. The flow of neutral interstellar helium through the perturbed interstellar plasma in the outer heliosheath (OHS) results in creation of the secondary population of interstellar He atoms, the so-called Warm Breeze, due to charge exchange with perturbed ions. The secondary population brings the imprint of the OHS conditions to the IBEX-Lo instrument. Based on a global simulation of the heliosphere with measurement-based parameters and detailed kinetic simulation of the filtration of He in the OHS, we find the number density of interstellar He$^+$ population at $(8.98\pm 0.12)\times 10^{-3}$~cm$^{-3}$. With this, we obtain the absolute density of interstellar H$^+$ $5.4\times 10^{-2}$~cm$^{-3}$ and electrons $6.3\times 10^{-2}$~cm$^{-3}$, and ionization degrees of H 0.26 and He 0.37. The results agree with estimates of the Very Local Interstellar Matter parameters obtained from fitting the observed spectra of diffuse interstellar EUV and soft X-Ray background., Comment: Accepted for ApJ

We present a comprehensive study of energetic neutral atoms (ENAs) of 10 eV to 2.5 keV from the downwind hemisphere of the heliosphere. These ENAs are believed to originate mostly from pickup protons and solar wind protons in the inner heliosheath. This study includes all low-energy observations made with the Interstellar Boundary Explorer over the first 8 years. Since the protons around 0.1 keV dominate the plasma pressure in the inner heliosheath in downwind direction, these ENA observations offer the unique opportunity to constrain the plasma properties and dimensions of the heliosheath where no in-situ observations are available. We first derive energy spectra of ENA intensities averaged over time for 49 macropixels covering the entire downwind hemisphere. The results confirm previous studies regarding integral intensities and the roll-over around 0.1 keV energy. With the expanded dataset we now find that ENA intensities at 0.2 and 0.1 keV seem to anti-correlate with solar activity. We then derive the product of total plasma pressure and emission thickness of protons in the heliosheath to estimate lower limits on the thickness of the inner heliosheath. The temporally averaged ENA intensities support a rather spherical shape of the termination shock and a heliosheath thickness between 150 and 210 au for most regions of the downwind hemisphere. Around the nominal downwind direction of 76{\deg} ecliptic longitude, the heliosheath is at least 280 au thick. There, the neutral hydrogen density seems to be depleted compared to upwind directions by roughly a factor of 2., Comment: Preprint of article in The Astrophysical Journal

The Interstellar Boundary Explorer (IBEX) has now operated in space for 7 years and returned nearly continuous observations that have led to scientific discoveries and reshaped our entire understanding of the outer heliosphere and its interaction with the local interstellar medium. Here we extend prior work, adding the 2014-2015 data for the first time, and examine, validate, initially analyze, and provide a complete 7-year set of Energetic Neutral Atom (ENA) observations from ~0.1 to 6 keV. The data, maps, and documentation provided here represent the 10th major release of IBEX data and include improvements to various prior corrections to provide the citable reference for the current version of IBEX data. We are now able to study time variations in the outer heliosphere and interstellar interaction over more than half a solar cycle. We find that the Ribbon has evolved differently than the globally distributed flux (GDF), with a leveling off and partial recovery of ENAs from the GDF, owing to solar wind output flattening and recovery. The Ribbon has now also lost its latitudinal ordering, which reflects the breakdown of solar minimum solar wind conditions and exhibits a greater time delay than for the surrounding GDF. Together, the IBEX observations strongly support a secondary ENA source for the Ribbon, and we suggest that this be adopted as the nominal explanation of the Ribbon going forward., Comment: Pre-print, 59 pages, 33 figures, published in ApJS

With the velocity vector and temperature of the pristine interstellar neutral (ISN) He recently obtained with high precision from a coordinated analysis summarized by McComas et al.2015b, we analyzed the IBEX observations of neutral He left out from this analysis. These observations were collected during the ISN observation seasons 2010---2014 and cover the region in the Earth's orbit where the Warm Breeze persists. We used the same simulation model and a very similar parameter fitting method to that used for the analysis of ISN He. We approximated the parent population of the Warm Breeze in front of the heliosphere with a homogeneous Maxwell-Boltzmann distribution function and found a temperature of $\sim 9\,500$ K, an inflow speed of 11.3 km s$^{-1}$, and an inflow longitude and latitude in the J2000 ecliptic coordinates $251.6^\circ$, $12.0^\circ$. The abundance of the Warm Breeze relative to the interstellar neutral He is 5.7\% and the Mach number is 1.97. The newly found inflow direction of the Warm Breeze, the inflow directions of ISN H and ISN He, and the direction to the center of IBEX Ribbon are almost perfectly co-planar, and this plane coincides within relatively narrow statistical uncertainties with the plane fitted only to the inflow directions of ISN He, ISN H, and the Warm Breeze. This co-planarity lends support to the hypothesis that the Warm Breeze is the secondary population of ISN He and that the center of the Ribbon coincides with the direction of the local interstellar magnetic field. The common plane for the direction of inflow of ISN gas, ISN H, the Warm Breeze, and the local interstellar magnetic field %includes the Sun and is given by the normal direction: ecliptic longitude $349.7^\circ \pm 0.6^\circ$ and latitude $35.7^\circ \pm 0.6^\circ$ in the J2000 coordinates, with the correlation coefficient of 0.85., Comment: Accepted for publication in Astrophysical Journal Supplement Series

We present the discovery of the first Neptune analog exoplanet or super-Earth with Neptune-like orbit, MOA-2013-BLG-605Lb. This planet has a mass similar to that of Neptune or a super-Earth and it orbits at $9\sim 14$ times the expected position of the snow-line, $a_{\rm snow}$, which is similar to Neptune's separation of $ 11\,a_{\rm snow}$ from the Sun. The planet/host-star mass ratio is $q=(3.6\pm0.7)\times 10^{-4}$ and the projected separation normalized by the Einstein radius is $s=2.39\pm0.05$. There are three degenerate physical solutions and two of these are due to a new type of degeneracy in the microlensing parallax parameters, which we designate "the wide degeneracy". The three models have (i) a Neptune-mass planet with a mass of $M_{\rm p}=21_{-7}^{+6} M_{Earth}$ orbiting a low-mass M-dwarf with a mass of $M_{\rm h}=0.19_{-0.06}^{+0.05} M_\odot$, (ii) a mini-Neptune with $M_{\rm p}= 7.9_{-1.2}^{+1.8} M_{Earth}$ orbiting a brown dwarf host with $M_{\rm h}=0.068_{-0.011}^{+0.019} M_\odot$ and (iii) a super-Earth with $M_{\rm p}= 3.2_{-0.3}^{+0.5} M_{Earth}$ orbiting a low-mass brown dwarf host with $M_{\rm h}=0.025_{-0.004}^{+0.005} M_\odot$ which is slightly favored. The 3-D planet-host separations are 4.6$_{-1.2}^{+4.7}$ AU, 2.1$_{-0.2}^{+1.0}$ AU and 0.94$_{-0.02}^{+0.67}$ AU, which are $8.9_{-1.4}^{+10.5}$, $12_{-1}^{+7}$ or $14_{-1}^{+11}$ times larger than $a_{\rm snow}$ for these models, respectively. The Keck AO observation confirm that the lens is faint. This discovery suggests that low-mass planets with Neptune-like orbit are common. So processes similar to the one that formed Neptune in our own Solar System or cold super-Earth may be common in other solar systems., Comment: 54 pages, 10 figures, 13 tables, Accepted for publication in the ApJ

This paper is one of three companion papers presenting the results of our in-depth analysis of the interstellar neutral helium (ISN He) observations carried out using the IBEX-Lo during the first six Interstellar Boundary Explorer (IBEX) observation seasons. We derive corrections for losses due to the limited throughput of the interface buffer and determine the IBEX spin-axis pointing. We develop an uncertainty system for the data, taking into account the resulting correlations between the data points. This system includes uncertainties due to Poisson statistics, background, spin-axis determination, systematic deviation of the boresight from the prescribed position, correction for the interface buffer losses, and the expected Warm Breeze (WB) signal. Subsequently, we analyze the data from 2009 to examine the role of various components of the uncertainty system. We show that the ISN He flow parameters are in good agreement with the values obtained by the original analysis. We identify the WB as the principal contributor to the global $\chi^2$ values in previous analyses. Other uncertainties have a much milder role and their contributions are comparable to each other. The application of this uncertainty system reduced the minimum $\chi^2$ value 4-fold. The obtained $\chi^2$ value, still exceeding the expected value, suggests that either the uncertainty system may still be incomplete or the adopted physical model lacks a potentially important element, which is likely an imperfect determination of the WB parameters. The derived corrections and uncertainty system are used in the accompanying paper by Bzowski et al. in an analysis of the data from six seasons., Comment: 43 pages, 9 figures

We analyzed observations of interstellar neutral helium (ISN~He) obtained from the Interstellar Boundary Explorer (IBEX) satellite during its first six years of operation. We used a refined version of the ISN~He simulation model, presented in the companion paper by Sokol_et al. 2015, and a sophisticated data correlation and uncertainty system and parameter fitting method, described in the companion paper by Swaczyna et al 2015. We analyzed the entire data set together and the yearly subsets, and found the temperature and velocity vector of ISN~He in front of the heliosphere. As seen in the previous studies, the allowable parameters are highly correlated and form a four-dimensional tube in the parameter space. The inflow longitudes obtained from the yearly data subsets show a spread of ~6 degree, with the other parameters varying accordingly along the parameter tube, and the minimum chi-square value is larger than expected. We found, however, that the Mach number of the ISN~He flow shows very little scatter and is thus very tightly constrained. It is in excellent agreement with the original analysis of ISN~He observations from IBEX and recent reanalyses of observations from Ulysses. We identify a possible inaccuracy in the Warm Breeze parameters as the likely cause of the scatter in the ISN~He parameters obtained from the yearly subsets, and we suppose that another component may exist in the signal, or a process that is not accounted for in the current physical model of ISN~He in front of the heliosphere. From our analysis, the inflow velocity vector, temperature, and Mach number of the flow are equal to lambda_ISNHe = 255.8 +/- 0.5 degree, beta_ISNHe = 5.16 +/- 0.10 degree, T_ISNHe = 7440 +/- 260 K, v_ISNHe = 25.8 +/- 0.4$ km/s, and M_ISNHe = 5.079 +/- 0.028, with uncertainties strongly correlated along the parameter tube., Comment: Updated references

Neutral interstellar helium has been observed by the Interstellar Boundary Explorer (IBEX) since 2009 with a signal-to-noise ratio well above 1000. Because of the geometry of the observations, the signal observed from January to March each year is the easiest to identify. However, as we show via simulations, the portion of the signal in the range of intensities from 10^{-3} to 10^{-2} of the peak value, previously mostly left out from the analysis, may bring important information about the details of the distribution function of interstellar He gas in front of the heliosphere. In particular, these observations may inform us about possible departures of the parent interstellar He population from equilibrium. We compare the expected distribution of the signal for the canonical assumption of a single Maxwell-Boltzmann population with the distributions for a superposition of the Maxwell-Boltzmann primary population and the recently discovered Warm Breeze, and for a single primary population given by a kappa function. We identify the regions on the sky where the differences between those cases are expected to be the most visible against the background. We discuss the diagnostic potential of the fall peak of the interstellar signal, reduced by a factor of 50 due to the Compton-Getting effect but still above the detection limit of IBEX. We point out the strong energy dependence of the fall signal and suggest that searching for this signal in the data could bring an independent assessment of the low-energy measurement threshold of the IBEX-Lo sensor., Comment: This paper is a part of the special issue of Astrophysical Journal Supplement Series on interstellar neutrals measured by IBEX

We present optical observations of the peculiar Type Ibn supernova (SN Ibn) OGLE-2012-SN-006, discovered and monitored by the OGLE-IV survey, and spectroscopically followed by PESSTO at late phases. Stringent pre-discovery limits constrain the explosion epoch with fair precision to JD = 2456203.8 +- 4.0. The rise time to the I-band light curve maximum is about two weeks. The object reaches the peak absolute magnitude M(I) = -19.65 +- 0.19 on JD = 2456218.1 +- 1.8. After maximum, the light curve declines for about 25 days with a rate of 4 mag per 100d. The symmetric I-band peak resembles that of canonical Type Ib/c supernovae (SNe), whereas SNe Ibn usually exhibit asymmetric and narrower early-time light curves. Since 25 days past maximum, the light curve flattens with a decline rate slower than that of the 56Co to 56Fe decay, although at very late phases it steepens to approach that rate. An early-time spectrum is dominated by a blue continuum, with only a marginal evidence for the presence of He I lines marking this SN Type. This spectrum shows broad absorptions bluewards than 5000A, likely O II lines, which are similar to spectral features observed in super-luminous SNe at early epochs. The object has been spectroscopically monitored by PESSTO from 90 to 180 days after peak, and these spectra show the typical features observed in a number of SN 2006jc-like events, including a blue spectral energy distribution and prominent and narrow (v(FWHM) ~ 1900 km/s) He I emission lines. This suggests that the ejecta are interacting with He-rich circumstellar material. The detection of broad (10000 km/s) O I and Ca II features likely produced in the SN ejecta (including the [O I] 6300A,6364A doublet in the latest spectra) lends support to the interpretation of OGLE-2012-SN-006 as a core-collapse event., Comment: 16 pages, 9 figures, 4 tables. Accepted by MNRAS

We have analyzed a sample of 27,258 fundamental-mode RR Lyrae variable stars (type RRab) detected recently toward the Galactic bulge by the Optical Gravitational Lensing Experiment (OGLE) survey. The data support our earlier claim that these metal-poor stars trace closely the barred structure formed of intermediate-age red clump giants. The distance to the Galactic center (GC) inferred from the bulge RR Lyrae stars is R_0=8.27+/-0.01(stat)+/-0.40(sys) kpc. We show that their spatial distribution has the shape of a triaxial ellipsoid with an major axis located in the Galactic plane and inclined at an angle of i=20+/-3 deg to the Sun-GC line of sight. The obtained scale-length ratio of the major axis to the minor axis in the Galactic plane and to the axis vertical to the plane is 1:0.49(2):0.39(2). We do not see the evidence for the bulge RR Lyrae stars forming an X-shaped structure. Based on the light curve parameters, we derive metallicities of the RRab variables and show that there is a very mild but statistically significant radial metallicity gradient. About 60% of the bulge RRab stars form two very close sequences on the period-amplitude (or Bailey) diagram, which we interpret as two major old bulge populations: A and B. Their metallicities likely differ. Population A is about four times less abundant than the slightly more metal-poor population B. Most of the remaining stars seem to represent other, even more metal-poor populations of the bulge. The presence of multiple old populations indicates that the Milky Way bulge was initially formed through mergers., Comment: ApJ, in press

We report the discovery of two distinct types of ellipsoidal binary systems occupying, so called, sequence E on the period-luminosity (P-L) diagram. We propose that steeper P-L relation is composed of giant-dwarf binaries, while the other consists of giant-giant binary systems. Analysis is based on a sample of 5334 objects, which we select from the OGLE-III survey data toward the the Large Magellanic Cloud. We show that one of the components of ellipsoidal binaries is typically either a Red Clump or a Red Giant Branch star, which leads to clear separation split of the sequence E at P = 40 d. In its short-period part, we identify two subsequences corresponding to the two types of binary systems (E2 and E3), while in the longer-period part the two groups merge forming a single subsequence E1. We extract a group of 271 ellipsoidal systems with eccentric orbits, from our sample. We present the period-luminosity relation they follow., Comment: 16 pages, 15 figures, submitted to Acta Astronomica

We report the discovery of a Jupiter-mass planet orbiting an M-dwarf star that gave rise to the microlensing event OGLE-2011-BLG-0265. Such a system is very rare among known planetary systems and thus the discovery is important for theoretical studies of planetary formation and evolution. High-cadence temporal coverage of the planetary signal combined with extended observations throughout the event allows us to accurately model the observed light curve. The final microlensing solution remains, however, degenerate yielding two possible configurations of the planet and the host star. In the case of the preferred solution, the mass of the planet is $M_{\rm p} = 0.9\pm 0.3\ M_{\rm J}$, and the planet is orbiting a star with a mass $M = 0.22\pm 0.06\ M_\odot$. The second possible configuration (2$\sigma$ away) consists of a planet with $M_{\rm p}=0.6\pm 0.3\ M_{\rm J}$ and host star with $M=0.14\pm 0.06\ M_\odot$. The system is located in the Galactic disk 3 -- 4 kpc towards the Galactic bulge. In both cases, with an orbit size of 1.5 -- 2.0 AU, the planet is a "cold Jupiter" -- located well beyond the "snow line" of the host star. Currently available data make the secure selection of the correct solution difficult, but there are prospects for lifting the degeneracy with additional follow-up observations in the future, when the lens and source star separate., Comment: 10 pages, 2 tables, 5 figures. Accepted in ApJ

We present the most comprehensive picture ever obtained of the central parts of the Milky Way probed with RR Lyrae variable stars. This is a collection of 38257 RR Lyr stars detected over 182 square degrees monitored photometrically by the Optical Gravitational Lensing Experiment (OGLE) in the most central regions of the Galactic bulge. The sample consists of 16804 variables found and published by the OGLE collaboration in 2011 and 21453 RR Lyr stars newly detected in the photometric databases of the fourth phase of the OGLE survey (OGLE-IV). 93% of the OGLE-IV variables were previously unknown. The total sample consists of 27258 RRab, 10825 RRc, and 174 RRd stars. We provide OGLE-IV I- and V-band light curves of the variables along with their basic parameters. About 300 RR Lyr stars in our collection are plausible members of 15 globular clusters. Among others, we found the first pulsating variables that may belong to the globular cluster Terzan 1 and the first RRd star in the globular cluster M54. Our survey also covers the center and outskirts of the Sagittarius Dwarf Spheroidal Galaxy enabling studies of the spatial distribution of the old stellar population from this galaxy. A group of double-mode RR Lyr stars with period ratios around 0.740 form a stream in the sky that may be a relic of a cluster or a dwarf galaxy tidally disrupted by the Milky Way. Three of our RR Lyr stars experienced a pulsation mode switching from double-mode to single fundamental mode or vice versa. We also present the first known RRd stars with large-amplitude Blazhko effect., Comment: 21 pages, 9 figures, submitted to Acta Astronomica

We present the design and first results of a real-time search for transients within the 650 sq. deg. area around the Magellanic Clouds, conducted as part of the OGLE-IV project and aimed at detecting supernovae, novae and other events. The average sampling of about 4 days from September to May, yielded a detection of 238 transients in 2012/2013 and 2013/2014 seasons. The superb photometric and astrometric quality of the OGLE data allows for numerous applications of the discovered transients. We use this sample to prepare and train a Machine Learning-based automated classifier for early light curves, which distinguishes major classes of transients with more than 80% of correct answers. Spectroscopically classified 49 supernovae Type Ia are used to construct a Hubble Diagram with statistical scatter of about 0.3 mag and fill the least populated region of the redshifts range in the Union sample. We investigate the influence of host galaxy environments on supernovae statistics and find the mean host extinction of A_I=0.19+-0.10 mag and A_V=0.39+-0.21 mag based on a subsample of supernovae Type Ia. We show that the positional accuracy of the survey is of the order of 0.5 pixels (0.13 arcsec) and that the OGLE-IV Transient Detection System is capable of detecting transients within the nuclei of galaxies. We present a few interesting cases of nuclear transients of unknown type. All data on the OGLE transients are made publicly available to the astronomical community via the OGLE website., Comment: Submitted to Acta Astronomica. Webpage: http://ogle.astrouw.edu.pl

We present the analysis of the gravitational microlensing event OGLE-2013-BLG-0102. The light curve of the event is characterized by a strong short-term anomaly superposed on a smoothly varying lensing curve with a moderate magnification $A_{\rm max}\sim 1.5$. It is found that the event was produced by a binary lens with a mass ratio between the components of $q = 0.13$ and the anomaly was caused by the passage of the source trajectory over a caustic located away from the barycenter of the binary. From the analysis of the effects on the light curve due to the finite size of the source and the parallactic motion of the Earth, the physical parameters of the lens system are determined. The measured masses of the lens components are $M_{1} = 0.096 \pm 0.013~M_{\odot}$ and $M_{2} = 0.012 \pm 0.002~M_{\odot}$, which correspond to near the hydrogen-burning and deuterium-burning mass limits, respectively. The distance to the lens is $3.04 \pm 0.31~{\rm kpc}$ and the projected separation between the lens components is $0.80 \pm 0.08~{\rm AU}$., Comment: 6 figures, 2 tables, ApJ submitted

We detect a cold, terrestrial planet in a binary-star system using gravitational microlensing. The planet has low mass (2 Earth masses) and lies projected at $a_{\perp,ph}$ ~ 0.8 astronomical units (AU) from its host star, similar to the Earth-Sun distance. However, the planet temperature is much lower, T<60 Kelvin, because the host star is only 0.10--0.15 solar masses and therefore more than 400 times less luminous than the Sun. The host is itself orbiting a slightly more massive companion with projected separation $a_{\perp,ch}=$10--15 AU. Straightforward modification of current microlensing search strategies could increase their sensitivity to planets in binary systems. With more detections, such binary-star/planetary systems could place constraints on models of planet formation and evolution. This detection is consistent with such systems being very common., Comment: Published in Science, Main and supplementary material combined

We investigate the signals from neutral He atoms observed from Earth orbit in 2010 by IBEX. The full He signal observed during the 2010 observation season can be explained as a superposition of pristine neutral interstellar He gas and an additional population of neutral He that we call the Warm Breeze. The Warm Breeze is approximately two-fold slower and 2.5 times warmer than the primary interstellar He population, and its density in front of the heliosphere is ~7% that of the neutral interstellar helium. The inflow direction of the Warm Breeze differs by ~19deg from the inflow direction of interstellar gas. The Warm Breeze seems a long-term feature of the heliospheric environment. It has not been detected earlier because it is strongly ionized inside the heliosphere, which brings it below the threshold of detection via pickup ion and heliospheric backscatter glow observations, as well as by the direct sampling of GAS/Ulysses. Possible sources for the Warm Breeze include (1) the secondary population of interstellar helium, created via charge exchange and perhaps elastic scattering of neutral interstellar He atoms on interstellar He+ ions in the outer heliosheath, or (2) a gust of interstellar He originating from a hypothetic wave train in the Local Interstellar Cloud. A secondary population is expected from models, but the characteristics of the Warm Breeze do not fully conform to modeling results. If, nevertheless, this is the explanation, IBEX-Lo observations of the Warm Breeze provide key insights into the physical state of plasma in the outer heliosheath. If the second hypothesis is true, the source is likely to be located within a few thousand of AU from the Sun, which is the propagation range of possible gusts of interstellar neutral helium with the Warm Breeze characteristics against dissipation via elastic scattering in the Local Cloud., Comment: submitted to ApJS

We report the discovery of a young stellar bridge, that forms a continuous connection between the Magellanic Clouds. This finding is based on number density maps for stellar populations found in data gathered by OGLE-IV, that fully cover over 270 deg^2 of the sky in the Magellanic Bridge area. This is the most extensive optical survey of this region up to date. We find that the young population is present mainly in the western half of the MBR, which, together with the newly discovered young population in the eastern Bridge, form a continuous stream of stars connecting both galaxies along Dec ~ -73.5 deg. The young population distribution is clumped, with one of the major densities close to the SMC, and the other, fairly isolated and located approximately mid-way between the Clouds, which we call the OGLE island. These overdensities are well matched by HI surface density contours, although the newly found young population in the eastern Bridge is offset by ~2 deg north from the highest HI density contour. We observe a continuity of red clump stars between the Magellanic Clouds, which represent an intermediate-age population. Red clump stars are present mainly in the southern and central parts of the Magellanic Bridge, below its gaseous part, and their presence is reflected by a strong deviation from the radial density profiles of the two galaxies. This may indicate either a tidal stream of stars, or that the stellar halos of the two galaxies overlap. On the other hand, we do not observe such an overlap within an intermediate-age population represented by the top of the red giant branch and the asymptotic giant branch stars. We also see only minor mixing of the old populations of the Clouds in the southern part of the Bridge, represented by the lowest part of the red giant branch., Comment: 18 pages, 14 Figures. Published in ApJ. On-line data available from http://ogle.astrouw.edu.pl

We analyse long-term optical photometry for two Galactic recurrent novae (V745 Sco and V3890 Sgr) and one Large Magellanic Cloud object (Nova LMC 1990b) observed over several years by the Optical Gravitational Lensing Experiment (OGLE) sky survey. We do not find variability with the previously claimed orbital period of V745 Sco. This voids previous findings based on this value (e.g. the distance determination). The quiescence variability of this object is dominated by semiregular pulsations of the red giant secondary (with periods of 136.5 and 77.4 d). The photometry of Nova LMC 1990b reveals an unnoticed eruption in 2010 and eclipse-like variability in quiescence with a period of 1.26432(8) d. The photometric properties make this object very similar to U Sco. Finally, we describe the eruptions of two likely symbiotic novae, V5590 Sgr and OGLE-2011-BLG-1444. The secondary of V5590 Sgr is a Mira star with pulsation period of 236 d., Comment: 8 pages, 7 figures, published in MNRAS

Analysis of IBEX measurements of neutral interstellar He flux brought the inflow velocity vector different from the results of earlier analysis of observations from GAS/Ulysses. Recapitulation of results on the helium inflow direction from the past ~40 years suggested that the inflow direction may be changing with time. We reanalyze the old Ulysses data and reprocess them to increase the accuracy of the instrument pointing to investigate if the GAS observations support the hypothesis that the interstellar helium inflow direction is changing. We employ a similar analysis method as in the analysis of the IBEX data. We seek a parameter set that minimizes reduced chi-squared, using the Warsaw Test Particle Model for the interstellar He flux at Ulysses with a state of the art model of neutral He ionization in the heliosphere, and precisely reproducing the observation conditions. We also propose a supplementary method of constraining the parameters based on cross-correlations of parameters obtained from analysis of carefully selected subsets of data. We find that the ecliptic longitude and speed of interstellar He are in a very good agreement with the values reported in the original GAS analysis. We find, however, that the temperature is markedly higher. The 3-seasons optimum parameter set is lambda = 255.3, beta = 6, v = 26.0 km/s, T = 7500 K. We find no evidence that it is varying with time, but the uncertainty range is larger than originally reported. The originally-derived parameters of interstellar He from GAS are in good agreement with presently derived, except for the temperature, which seems to be appreciably higher, in good agreement with interstellar absorption line results. While the results of the present analysis are in marginal agreement with the earlier reported results from IBEX, the most likely values from the two analyses differ for reasons that are still not understood., Comment: submitted for publication in Astronomy & Astrophysics

We report the discovery of a massive planet OGLE-2008-BLG-355Lb. The light curve analysis indicates a planet:host mass ratio of q = 0.0118 +/- 0.0006 at a separation of 0.877 +/- 0.010 Einstein radii. We do not measure a significant microlensing parallax signal and do not have high angular resolution images that could detect the planetary host star. Therefore, we do not have a direct measurement of the host star mass. A Bayesian analysis, assuming that all host stars have equal probability to host a planet with the measured mass ratio implies a host star mass of M_h = 0.37_{-0.17}^{+0.30} M_Sun and a companion of mass M_P = 4.6^{+3.7}_{-2.2} M_Jup, at a projected separation of r_proj = 1.70^{+0.29}_{-0.30} AU. The implied distance to the planetary system is D_L = 6.8 +/- 1.1 kpc. A planetary system with the properties preferred by the Bayesian analysis would be a challenge to the core-accretion model of planet formation, as the core-accretion model predicts that massive planets are far more likely to form around more massive host stars. This core accretion model prediction is not consistent with our Bayesian prior of an equal probability of host stars of all masses to host a planet with the measured mass ratio. So, if the core accretion model prediction is right, we should expect that follow-up high angular resolution observations will detect a host star with a mass in the upper part of the range allowed by the Bayesian analysis. That is, the host would probably be a K or G dwarf., Comment: 42 pages, 12 figures, submitted for publication in ApJ

We report the discovery of an RR Lyrae star that experienced a switching of its pulsation mode. OGLE-BLG-RRLYR-12245 was discovered as a double-mode RRd star from the observations conducted in years 2001-2006 during the third phase of the Optical Gravitational Lensing Experiment (OGLE-III). The OGLE-IV observations carried out since 2010 reveal that this object is now a fundamental-mode RRab star, with no sign of the first-overtone pulsation. The analysis of the OGLE photometry shows that the final stage of the mode switching occurred on a relatively short timescale of a few months in 2005. We study the behavior of the star during this process, showing changes of the pulsational amplitudes and periods. We also discuss possible causes for the mode switching in RR Lyr stars., Comment: 9 pages, 4 figures, Acta Astronomica accepted

The mass of the lenses giving rise to Galactic microlensing events can be constrained by measuring the relative lens-source proper motion and lens flux. The flux of the lens can be separated from that of the source, companions to the source, and unrelated nearby stars with high-resolution images taken when the lens and source are spatially resolved. For typical ground-based adaptive optics (AO) or space-based observations, this requires either inordinately long time baselines or high relative proper motions. We provide a list of microlensing events toward the Galactic Bulge with high relative lens-source proper motion that are therefore good candidates for constraining the lens mass with future high-resolution imaging. We investigate all events from 2004 -- 2013 that display detectable finite-source effects, a feature that allows us to measure the proper motion. In total, we present 20 events with mu >~ 8 mas/yr. Of these, 14 were culled from previous analyses while 6 are new, including OGLE-2004-BLG-368, MOA-2005-BLG-36, OGLE-2012-BLG-0211, OGLE-2012-BLG-0456, MOA-2012-BLG-532, and MOA-2013-BLG-029. In <~12 years the lens and source of each event will be sufficiently separated for ground-based telescopes with AO systems or space telescopes to resolve each component and further characterize the lens system. Furthermore, for the most recent events, comparison of the lens flux estimates from images taken immediately to those estimated from images taken when the lens and source are resolved can be used to empirically check the robustness of the single-epoch method currently being used to estimate lens masses for many events., Comment: 10 pages, 8 figures, 6 tables, submitted to ApJ. For a brief video explaining the key results of this paper, please visit: https://www.youtube.com/watch?v=i_dzT8NydJI

We report the discovery of MOA-2013-BLG-220Lb, which has a super-Jupiter mass ratio $q=3.01\pm 0.02\times 10^{-3}$ relative to its host. The proper motion, $\mu=12.5\pm 1\, {\rm mas}\,{\rm yr}^{-1}$, is one of the highest for microlensing planets yet discovered, implying that it will be possible to separately resolve the host within $\sim 7$ years. Two separate lines of evidence imply that the planet and host are in the Galactic disk. The planet could have been detected and characterized purely with follow-up data, which has important implications for microlensing surveys, both current and into the LSST era., Comment: 2 tables, 2 figures

Characterizing a microlensing planet is done from modeling an observed lensing light curve. In this process, it is often confronted that solutions of different lensing parameters result in similar light curves, causing difficulties in uniquely interpreting the lens system, and thus understanding the causes of different types of degeneracy is important. In this work, we show that incomplete coverage of a planetary perturbation can result in degenerate solutions even for events where the planetary signal is detected with a high level of statistical significance. We demonstrate the degeneracy for an actually observed event OGLE-2012-BLG-0455/MOA-2012-BLG-206. The peak of this high-magnification event $(A_{\rm max}\sim400)$ exhibits very strong deviation from a point-lens model with $\Delta\chi^{2}\gtrsim4000$ for data sets with a total number of measurement 6963. From detailed modeling of the light curve, we find that the deviation can be explained by four distinct solutions, i.e., two very different sets of solutions, each with a two-fold degeneracy. While the two-fold (so-called "close/wide") degeneracy is well-understood, the degeneracy between the radically different solutions is not previously known. The model light curves of this degeneracy differ substantially in the parts that were not covered by observation, indicating that the degeneracy is caused by the incomplete coverage of the perturbation. It is expected that the frequency of the degeneracy introduced in this work will be greatly reduced with the improvement of the current lensing survey and follow-up experiments and the advent of new surveys., Comment: 5 pages, 3 figures, ApJ accepted

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons., Comment: 32 pages with 9 included figures

We present the results of a search for pulsating stars in the 7.12 deg^2 OGLE-III Galactic disk area in the direction tangent to the Centaurus Arm. We report the identification of 20 Classical Cepheids, 45 RR Lyr type stars, 31 Long-Period Variables, such as Miras and Semi-Regular Variables, one pulsating white dwarf, and 58 very likely delta Sct type stars. Based on asteroseismic models constructed for one quadruple-mode and six triple-mode delta Sct type pulsators, we estimated masses, metallicities, ages, and distance moduli to these objects. The modeled stars have masses in the range 0.9-2.5 M_sun and are located at distances between 2.5 kpc and 6.2 kpc. Two triple-mode and one double-mode pulsators seem to be Population II stars of the SX Phe type, probably from the Galactic halo. Our sample also includes candidates for Type II Cepheids and unclassified short-period (P<0.23 d) multi-mode stars which could be either delta Sct or beta Cep type stars. One of the detected variables is a very likely delta Sct star with an exceptionally high peak-to-peak I-band amplitude of 0.35 mag at the very short period of 0.0196 d. All reported pulsating variables but one object are new discoveries. They are included in the OGLE-III Catalog of Variable Stars. Finally, we introduce the on-going OGLE-IV Galactic Disk Survey, which covers more than half of the Galactic plane. For the purposes of future works on the spiral structure and star formation history of the Milky Way, we have already compiled a list of known Galactic Classical Cepheids., Comment: published in Acta Astronomica

We report analysis of high microlensing event MOA-2008-BLG-379, which has a strong microlensing anomaly at its peak, due to a massive planet with a mass ratio of q = 6.9 x 10^{-3}. Because the faint source star crosses the large resonant caustic, the planetary signal dominates the light curve. This is unusual for planetary microlensing events, and as a result, the planetary nature of this light curve was not immediately noticed. The planetary nature of the event was found when the MOA Collaboration conducted a systematic study of binary microlensing events previously identified by the MOA alert system. We have conducted a Bayesian analysis based on a standard Galactic model to estimate the physical parameters of the lens system. This yields a host star mass of M_L = 0.66_{-0.33}^{+0.29} M_Sun orbited by a planet of mass m_P = 4.8_{-2.4}^{+2.1} M_Jup at an orbital separation of a = 4.1_{-1.5}^{+1.9} AU at a distance of D_L = 3.6 +/- 1.3 kpc. The faint source magnitude of I_S = 21.30 and relatively high lens-source relative proper motion of mu_rel = 7.6 +/- 1.6 mas/yr implies that high angular resolution adaptive optics or Hubble Space Telescope observations are likely to be able to detect the source star, which would determine the masses and distance of the planet and its host star., Comment: 25 pages, 5 figures, accepted for publication in ApJ

We present a variability study of 4646 massive stars in the Small Magellanic Cloud (SMC) with known spectral types from the catalog of Bonanos et al. (2010) using the light curves from the OGLE-III database. The goal is to exploit the time domain information available through OGLE-III to gain insight into the processes that govern the evolution of massive stars. This variability survey of massive stars with known spectral types is larger than any previous survey by a factor of 7. We find that 60% of our sample (2766 stars) show no significant variability and 40% (1880 stars) exhibit variability distributed as follows: 807 stars display low-amplitude stochastic variability with fluctuations in I-band of up to 0.05 mag, 443 stars present irregular variability of higher amplitude (76% of these are reported as variables for the first time), 205 are eclipsing binaries (including 101 newly discovered systems), 50 are candidate rotating variables, 126 are classical Cepheids, 188 stars exhibit short-term sinusoidal periodicity (P < 3 days) making them candidate "slowly pulsating B stars" and non-radial Be pulsators, and 61 periodic stars exhibit longer periods. We demonstrate the wealth of information provided in the time domain, by doubling the number of known massive eclipsing binary systems and identifying 189 new candidate early-type Be and 20 Oe stars in the SMC. In addition, we find that ~80% of Be stars are photometrically variable in the OGLE-III time domain and provide evidence that short-term pulsating stars with additional photometric variability are rotating close to their break-up velocity., Comment: 46 pages, 18 figures, 11 tables. A&A in press. See http://media.wix.com/ugd/d2ba94_1596d7db762b496c89f21d03891f46c3.pdf for a version with full resolution figures

We present a large sample of eclipsing binary stars detected in the Small Magellanic Cloud fields covering about 14 square degrees that have been monitored for eight years during the third phase of the OGLE survey. This is the largest set of such variables containing 6138 objects, of which 777 are contact and 5361 non-contact binaries. The estimated completeness of this sample is around 82%. We analyze the statistical properties of the sample and present selected interesting objects: 32 systems having eccentric orbit with visible apsidal motion, one Transient Eclipsing Binary, ten RS CVn type stars, 22 still unexplained Double-Periodic Variable stars, and 15 candidates for doubly eclipsing quadruple systems. Based on the OGLE-III proper motions, we classified 47 binaries from our sample as foreground Galactic stars. We also list candidates suitable for the SMC distance determination., Comment: 16 pages, 12 figures, Acta Astronomica in-print