Search

Solar modelling has long been split into ''internal'' and ''surface'' modelling, because of the lack of tools to connect the very different scales in space and time, as well as the widely different environments and dominating physical effects involved. Significant efforts have recently been put into resolving this disconnect. We address the outstanding bottlenecks in connecting internal convection zone and dynamo simulations to the surface of the Sun, and conduct a proof-of-concept high resolution global simulation of the convection zone of the Sun, using the task-based DISPATCH code framework. We present a new `volleyball' mesh decomposition, which has Cartesian patches tessellated on a sphere with no singularities. We use our new entropy based HLLS approximate Riemann solver to model magneto-hydrodynamics in a global simulation, ranging between 0.655 -- 0.995 R$_\odot$, with an initial ambient magnetic field set to 0.1 Gauss. The simulations develop convective motions with complex, turbulent structures. Small-scale dynamo action twists the ambient magnetic field and locally amplifies magnetic field magnitudes by more than two orders of magnitude within the initial run-time., Comment: 12 pages, 9 figures, submitted to A&A. Movies available online

With advance of supercomputers we can now afford simulations with very large range of scales. In astrophysical applications, e.g. simulating Solar, stellar and planetary atmospheres, physical quantities, like gas pressure, density, temperature and plasma $\beta$ can vary by orders of magnitude. This requires a robust solver, which can deal with a very wide range of conditions and be able to maintain hydrostatic equilibrium. We reformulate a Godunov-type HLLD Riemann solver so it would be suitable to maintain hydrostatic equilibrium in atmospheric applications and would be able to handle low and high Mach numbers, regimes where kinetic and magnetic energies dominate over thermal energy without any ad-hoc corrections. We change the solver to use entropy instead of total energy as the 'energy' variable in the system of MHD equations. The entropy is *not conserved*, it increases when kinetic and magnetic energy is converted to heat, as it should. We conduct a series of standard tests with varying conditions and show that the new formulation for the Godunot type Riemann solver works and is very promising., Comment: 12 pages, 14 figures, submitted to A&A

Solar coronal loops are commonly subject to oscillations. Observations of coronal oscillations are used to infer physical properties of the coronal plasma using coronal seismology. Excitation and evolution of oscillations in coronal loops is typically studied using highly idealised models of magnetic flux-tubes. In order to improve our understanding of coronal oscillations, it is necessary to consider the effect of realistic magnetic field topology and evolution. We study excitation and evolution of coronal oscillations in three-dimensional self-consistent simulations of solar atmosphere spanning from convection zone to solar corona using radiation-MHD code Bifrost. We use forward-modelled EUV emission and three-dimensional tracing of magnetic field to analyse oscillatory behaviour of individual magnetic loops. We further analyse the evolution of individual plasma velocity components along the loops using wavelet power spectra to capture changes in the oscillation periods. Various types of oscillations commonly observed in the corona are present in the simulation. We detect standing oscillations in both transverse and longitudinal velocity components, including higher order oscillation harmonics. We also show that self-consistent simulations reproduce existence of two distinct regimes of transverse coronal oscillations: rapidly decaying oscillations triggered by impulsive events and sustained small-scale oscillations showing no observable damping. No harmonic drivers are detected at the footpoints of oscillating loops. We show that coronal loop oscillations are abundant in self-consistent 3D MHD simulations of the solar atmosphere. The dynamic evolution and variability of individual magnetic loops suggest we need to reevaluate our models of monolithic and static coronal loops with constant lengths in favour of more realistic models., Comment: 13 pages, 11 figures, accepted for publication in A&A

NASA's Interface Region Imaging Spectrograph (IRIS) provides high resolution observations of the solar atmosphere through UV spectroscopy and imaging. Since the launch of IRIS in June 2013, we have conducted systematic observation campaigns in coordination with the Swedish 1-m Solar Telescope (SST) on La Palma. The SST provides complementary high-resolution observations of the photosphere and chromosphere. The SST observations include spectro-polarimetric imaging in photospheric Fe I lines and spectrally-resolved imaging in the chromospheric Ca II 8542 A, H-alpha, and Ca II K lines. We present a database of co-aligned IRIS and SST datasets that is open for analysis to the scientific community. The database covers a variety of targets including active regions, sunspots, plage, quiet Sun, and coronal holes., Comment: Accepted for publication in Astronomy & Astrophysics (July 21, 2020)

Coronal rain consists of cool and dense plasma condensations formed in coronal loops as a result of thermal instability. Previous numerical simulations of thermal instability and coronal rain formation have relied on artificially adding a coronal heating term to the energy equation. To reproduce large-scale characteristics of the corona, using more realistic coronal heating prescription is necessary. We analyse coronal rain formation and evolution in a 3-dimensional radiative magnetohydrodynamic simulation spanning from convection zone to corona which is self-consistently heated by magnetic field braiding as a result of convective motions. We investigate the spatial and temporal evolution of energy dissipation along coronal loops which become thermally unstable. Ohmic dissipation in the model leads to the heating events capable of inducing sufficient chromospheric evaporation into the loop to trigger thermal instability and condensation formation. The cooling of the thermally unstable plasma occurs on timescales comparable to the duration of the individual impulsive heating events. The impulsive heating has sufficient duration to trigger thermal instability in the loop but does not last long enough to lead to coronal rain limit cycles. We show that condensations can either survive and fall into the chromosphere or be destroyed by strong bursts of Joule heating associated with a magnetic reconnection events. In addition, we find that condensations can also form along open magnetic field lines., Comment: 10 pages, 8 figures, accepted by A&A

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q~0.03 which differs from that reported in Shvartzvald+2016. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s~0.15 or s~7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M-dwarf orbited by a massive Jupiter companion at very close (M_{host}=0.30^{+0.08}_{-0.06} M_{Sun}, M_{comp}=10.1^{+2.9}_{-2.2} M_{Jup}, a_{exp}=0.40^{+0.05}_{-0.04} au) or wide (M_{host}=0.28^{+0.10}_{-0.08} M_{Sun}, M_{comp}=9.9^{+3.8}_{-3.5}M_{Jup}, a_{exp}=18.0^{+3.2}_{-3.2} au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q=0.03 which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M_{Jup}. It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M_{Jup} in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M-dwarfs and even BDs which is challenging for other exoplanet detection methods., Comment: 22 pages, 10 figures, Accepted for publication in The Astronomical Journal

We present the analysis of stellar binary microlensing event OGLE-2015-BLG-0060 based on observations obtained from 13 different telescopes. Intensive coverage of the anomalous parts of the light curve was achieved by automated follow-up observations from the robotic telescopes of the Las Cumbres Observatory. We show that, for the first time, all main features of an anomalous microlensing event are well covered by follow-up data, allowing us to estimate the physical parameters of the lens. The strong detection of second-order effects in the event light curve necessitates the inclusion of longer-baseline survey data in order to constrain the parallax vector. We find that the event was most likely caused by a stellar binary-lens with masses $M_{\star1} = 0.87 \pm 0.12 M_{\odot}$ and $M_{\star2} = 0.77 \pm 0.11 M_{\odot}$. The distance to the lensing system is 6.41 $\pm 0.14$ kpc and the projected separation between the two components is 13.85 $\pm 0.16$ AU. Alternative interpretations are also considered., Comment: 13 pages, 5 figures, Published in MNRAS

We present broad-band photometry of eleven planetary transits of the hot Jupiter WASP-74b, using three medium-class telescopes and employing the telescope-defocussing technique. Most of the transits were monitored through I filters and one was simultaneously observed in five optical (U, g', r', i', z') and three near infrared (J, H, K) passbands, for a total of 18 light curves. We also obtained new high-resolution spectra of the host star. We used these new data to review the orbital and physical properties of the WASP-74 planetary system. We were able to better constrain the main system characteristics, measuring smaller radius and mass for both the hot Jupiter and its host star than previously reported in the literature. Joining our optical data with those taken with the HST in the near infrared, we built up an observational transmission spectrum of the planet, which suggests the presence of strong optical absorbers, as TiO and VO gases, in its atmosphere., Comment: 12 pages, 9 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

We investigate the effects of radiative energy transfer on a series of nested-grid, high-resolution hydrodynamic simulations of gas and particle dynamics in the vicinity of an Earth-mass planetary embryo. We include heating due to the accretion of solids and the subsequent convective motions. Using a constant embryo surface temperature, we show that radiative energy transport results in a tendency to reduce the entropy in the primordial atmosphere, but this tendency is alleviated by an increase in the strength of convective energy transport, triggered by a correspondingly increased super-adiabatic temperature gradient. As a consequence, the amplitude of the convective motions increase by roughly an order of magnitude in the vicinity of the embryo. In the cases investigated here, where the optical depth towards the disk surface is larger than unity, the reduction of the temperature in the outer parts of the Hill sphere relative to cases without radiative energy transport is only $\sim$100K, while the mass density increase is on the order of a factor of two in the inner parts of the Hill sphere. Our results demonstrate that, unless unrealistically low dust opacities are assumed, radiative cooling in the context of primordial rocky planet atmospheres can only become important after the disk surface density has dropped significantly below minimum-mass-solar-nebula values., Comment: 5 pages, 4 figures, 2 appendices; MNRAS Letters, accepted

(abridged) Using the particularly long gravitational microlensing event OGLE-2014-BLG-1186 with a time-scale $t_\mathrm{E}$ ~ 300 d, we present a methodology for identifying the nature of localised deviations from single-lens point-source light curves, which ensures that 1) the claimed signal is substantially above the noise floor, 2) the inferred properties are robustly determined and their estimation not subject to confusion with systematic noise in the photometry, 3) there are no alternative viable solutions within the model framework that might have been missed. Annual parallax and binarity could be separated and robustly measured from the wing and the peak data, respectively. We find matching model light curves that involve either a binary lens or a binary source. Our binary-lens models indicate a planet of mass $M_2$ = (45 $\pm$ 9) $M_\oplus$, orbiting a star of mass $M_1$ = (0.35 $\pm$ 0.06) $M_\odot$, located at a distance $D_\mathrm{L}$ = (1.7 $\pm$ 0.3) kpc from Earth, whereas our binary-source models suggest a brown-dwarf lens of $M$ = (0.046 $\pm$ 0.007) $M_\odot$, located at a distance $D_\mathrm{L}$ = (5.7 $\pm$ 0.9) kpc, with the source potentially being a (partially) eclipsing binary involving stars predicted to be of similar colour given the ratios between the luminosities and radii. The ambiguity in the interpretation would be resolved in favour of a lens binary by observing the luminous lens star separating from the source at the predicted proper motion of $\mu$ = (1.6 $\pm$ 0.3) mas yr$^{-1}$, whereas it would be resolved in favour of a source binary if the source could be shown to be a (partially) eclipsing binary matching the obtained model parameters. We experienced that close binary source stars pose a challenge for claiming the detection of planets by microlensing in events where the source passes very close to the lens star hosting the planet., Comment: 26 pages, 14 figures

We present nested-grid, high-resolution hydrodynamic simulations of gas and particle dynamics in the vicinity of Mars- to Earth-mass planetary embryos. The simulations extend from the surface of the embryos to a few vertical disk scale heights, with \rev{a spatial} dynamic range \rev{of} $\sim\! 1.4\times 10^5$. Our results confirm that "pebble"-sized particles are readily accreted, with accretion rates continuing to increase up to metre-size "boulders" for a 10\% MMSN surface density model. The gas mass flux in and out of the Hill sphere is consistent with the Hill rate, $\Sigma\Omega R_\mathrm{H}^2 = 4\, 10^{-3}$ M$_\oplus$ yr$^{-1}$. While smaller size particles mainly track the gas, a net accretion rate of $\approx 2\,10^{-5}$ M$_\oplus$ yr$^{-1}$ is reached for 0.3--1 cm particles, even though a significant fraction leaves the Hill sphere again. Effectively all pebble-sized particles that cross the Bondi sphere are accreted. The resolution of these simulations is sufficient to resolve accretion-driven convection. Convection driven by a nominal accretion rate of $10^{-6}$ M$_\oplus$ yr$^{-1}$ does not significantly alter the pebble accretion rate. We find that, due to cancellation effects, accretion rates of pebble-sized particles are nearly independent of disk surface density. As a result, we can estimate accurate growth times for specified particle sizes. For 0.3--1 cm size particles, the growth time from a small seed is $\sim$0.15 million years for an Earth mass planet at 1 AU and $\sim$0.1 million years for a Mars mass planet at 1.5 AU., Comment: 19 pages, 21 figures, update includes accepted version, MNRAS 479, 5136--5156 (2018)

We report the discovery of OGLE-2016-BLG-1190Lb, which is likely to be the first Spitzer microlensing planet in the Galactic bulge/bar, an assignation that can be confirmed by two epochs of high-resolution imaging of the combined source-lens baseline object. The planet's mass M_p= 13.4+-0.9 M_J places it right at the deuterium burning limit, i.e., the conventional boundary between "planets" and "brown dwarfs". Its existence raises the question of whether such objects are really "planets" (formed within the disks of their hosts) or "failed stars" (low mass objects formed by gas fragmentation). This question may ultimately be addressed by comparing disk and bulge/bar planets, which is a goal of the Spitzer microlens program. The host is a G dwarf M_host = 0.89+-0.07 M_sun and the planet has a semi-major axis a~2.0 AU. We use Kepler K2 Campaign 9 microlensing data to break the lens-mass degeneracy that generically impacts parallax solutions from Earth-Spitzer observations alone, which is the first successful application of this approach. The microlensing data, derived primarily from near-continuous, ultra-dense survey observations from OGLE, MOA, and three KMTNet telescopes, contain more orbital information than for any previous microlensing planet, but not quite enough to accurately specify the full orbit. However, these data do permit the first rigorous test of microlensing orbital-motion measurements, which are typically derived from data taken over <1% of an orbital period., Comment: 63 pages, 13 figures, 7 tables, AJ, in press

The formation and dynamical history of hot Jupiters is currently debated, with wide stellar binaries having been suggested as a potential formation pathway. Additionally, contaminating light from both binary companions and unassociated stars can significantly bias the results of planet characterisation studies, but can be corrected for if the properties of the contaminating star are known. We search for binary companions to known transiting exoplanet host stars, in order to determine the multiplicity properties of hot Jupiter host stars. We also characterise unassociated stars along the line of sight, allowing photometric and spectroscopic observations of the planetary system to be corrected for contaminating light. We analyse lucky imaging observations of 97 Southern hemisphere exoplanet host stars, using the Two Colour Instrument on the Danish 1.54m telescope. For each detected companion star, we determine flux ratios relative to the planet host star in two passbands, and measure the relative position of the companion. The probability of each companion being physically associated was determined using our two-colour photometry. A catalogue of close companion stars is presented, including flux ratios, position measurements, and estimated companion star temperature. For companions that are potential binary companions, we review archival and catalogue data for further evidence. For WASP-77AB and WASP-85AB, we combine our data with historical measurements to determine the binary orbits, showing them to be moderately eccentric and inclined to the line of sight and planetary orbital axis. Combining our survey with the similar Friends of Hot Jupiters survey, we conclude that known hot Jupiter host stars show a deficit of high mass stellar companions compared to the field star population; however, this may be a result of the biases in detection and target selection by ground-based surveys., Comment: 32 pages, 14 figures, 12 tables. The contents of tables 2, 3, 4, 9 and 10 are included in the arXiv source. Updated to correct typos and include tables in final format

We introduce a high-performance simulation framework that permits the semi-independent, task-based solution of sets of partial differential equations, typically manifesting as updates to a collection of `patches' in space-time. A hybrid MPI/OpenMP execution model is adopted, where work tasks are controlled by a rank-local `dispatcher' which selects, from a set of tasks generally much larger than the number of physical cores (or hardware threads), tasks that are ready for updating. The definition of a task can vary, for example, with some solving the equations of ideal magnetohydrodynamics (MHD), others non-ideal MHD, radiative transfer, or particle motion, and yet others applying particle-in-cell (PIC) methods. Tasks do not have to be grid-based, while tasks that are, may use either Cartesian or orthogonal curvilinear meshes. Patches may be stationary or moving. Mesh refinement can be static or dynamic. A feature of decisive importance for the overall performance of the framework is that time steps are determined and applied locally; this allows potentially large reductions in the total number of updates required in cases when the signal speed varies greatly across the computational domain, and therefore a corresponding reduction in computing time. Another feature is a load balancing algorithm that operates `locally' and aims to simultaneously minimise load and communication imbalance. The framework generally relies on already existing solvers, whose performance is augmented when run under the framework, due to more efficient cache usage, vectorisation, local time-stepping, plus near-linear and, in principle, unlimited OpenMP and MPI scaling., Comment: 17 pages, 8 figures. Accepted by MNRAS

We report 13 high-precision light curves of eight transits of the exoplanet WASP-52b, obtained by using four medium-class telescopes, through different filters, and adopting the defocussing technique. One transit was recorded simultaneously from two different observatories and another one from the same site but with two different instruments, including a multi-band camera. Anomalies were clearly detected in five light curves and modelled as starspots occulted by the planet during the transit events. We fitted the clean light curves with the jktebop code, and those with the anomalies with the prism+gemc codes in order to simultaneously model the photometric parameters of the transits and the position, size and contrast of each starspot. We used these new light curves and some from the literature to revise the physical properties of the WASP-52 system. Starspots with similar characteristics were detected in four transits over a period of 43 days. In the hypothesis that we are dealing with the same starspot, periodically occulted by the transiting planet, we estimated the projected orbital obliquity of WASP-52b to be lambda = 3.8 \pm 8.4 degree. We also determined the true orbital obliquity, psi = 20 \pm 50 degree, which is, although very uncertain, the first measurement of psi purely from starspot crossings. We finally assembled an optical transmission spectrum of the planet and searched for variations of its radius as a function of wavelength. Our analysis suggests a flat transmission spectrum within the experimental uncertainties., Comment: 16 pages, 12 figures, to appear in Monthly Notices of the Royal Astronomical Society

Aims. In this work we rigorously show the shortcomings of various simplifications that are used to calculate the total internal partition function. These shortcomings can lead to errors of up to 40 percent or more in the estimated partition function. These errors carry on to calculations of thermodynamic quantities. Therefore a more complicated approach has to be taken. Methods. Seven possible simplifications of various complexity are described, together with advantages and disadvantages of direct summation of experimental values. These were compared to what we consider the most accurate and most complete treatment (case 8). Dunham coefficients were determined from experimental and theoretical energy levels of a number of electronically excited states of H$_2$ . Both equilibrium and normal hydrogen was taken into consideration. Results. Various shortcomings in existing calculations are demonstrated, and the reasons for them are explained. New partition functions for equilibrium, normal, and ortho and para hydrogen are calculated and thermodynamic quantities are reported for the temperature range 1 - 20000 K. Our results are compared to previous estimates in the literature. The calculations are not limited to the ground electronic state, but include all bound and quasi-bound levels of excited electronic states. Dunham coefficients of these states of H$_2$ are also reported. Conclusions. For most of the relevant astrophysical cases it is strongly advised to avoid using simplifications, such as a harmonic oscillator and rigid rotor or ad hoc summation limits of the eigenstates to estimate accurate partition functions and to be particularly careful when using polynomial fits to the computed values. Reported internal partition functions and thermodynamic quantities in the present work are shown to be more accurate than previously available data., Comment: Accepted for publication in Astronomy & Astrophysics, 22 pages, 10 figures, 28 tables

We present a combined analysis of the observations of the gravitational microlensing event OGLE-2015-BLG-0479 taken both from the ground and by the {\it Spitzer Space Telescope}. The light curves seen from the ground and from space exhibit a time offset of $\sim 13$ days between the caustic spikes, indicating that the relative lens-source positions seen from the two places are displaced by parallax effects. From modeling the light curves, we measure the space-based microlens parallax. Combined with the angular Einstein radius measured by analyzing the caustic crossings, we determine the mass and distance of the lens. We find that the lens is a binary composed of two G-type stars with masses $\sim 1.0\ M_\odot$ and $\sim 0.9\ M_\odot$ located at a distance $\sim 3$ kpc. In addition, we are able to constrain the complete orbital parameters of the lens thanks to the precise measurement of the microlens parallax derived from the joint analysis. In contrast to the binary event OGLE-2014-BLG-1050, which was also observed by {\it Spitzer}, we find that the interpretation of OGLE-2015-BLG-0479 does not suffer from the degeneracy between $(\pm,\pm)$ and $(\pm,\mp)$ solutions, confirming that the four-fold parallax degeneracy in single-lens events collapses into the two-fold degeneracy for the general case of binary-lens events. The location of the blend in the color-magnitude diagram is consistent with the lens properties, suggesting that the blend is the lens itself. The blend is bright enough for spectroscopy and thus this possibility can be checked from future follow-up observations., Comment: 9 pages, 4 figures

Simultaneous observations of microlensing events from multiple locations allow for the breaking of degeneracies between the physical properties of the lensing system, specifically by exploring different regions of the lens plane and by directly measuring the "microlens parallax". We report the discovery of a 30-55$M_J$ brown dwarf orbiting a K dwarf in microlensing event OGLE-2015-BLG-1319. The system is located at a distance of $\sim$5 kpc toward the Galactic bulge. The event was observed by several ground-based groups as well as by $Spitzer$ and $Swift$, allowing the measurement of the physical properties. However, the event is still subject to an 8-fold degeneracy, in particular the well-known close-wide degeneracy, and thus the projected separation between the two lens components is either $\sim$0.25 AU or $\sim$45 AU. This is the first microlensing event observed by $Swift$, with the UVOT camera. We study the region of microlensing parameter space to which $Swift$ is sensitive, finding that while for this event $Swift$ could not measure the microlens parallax with respect to ground-based observations, it can be important for other events. Specifically, for detecting nearby brown dwarfs and free-floating planets in high magnification events., Comment: 28 pages, 6 figures, 2 tables. Submitted to ApJ

We show the benefits of using Electron-Multiplying CCDs and the shift-and-add technique as a tool to minimise the effects of the atmospheric turbulence such as blending between stars in crowded fields and to avoid saturated stars in the fields observed. We intend to complete, or improve, the census of the variable star population in globular cluster NGC~6715. Our aim is to obtain high-precision time-series photometry of the very crowded central region of this stellar system via the collection of better angular resolution images than has been previously achieved with conventional CCDs on ground-based telescopes. Observations were carried out using the Danish 1.54-m Telescope at the ESO La Silla observatory in Chile. The telescope is equipped with an Electron-Multiplying CCD that allowed to obtain short-exposure-time images (ten images per second) that were stacked using the shift-and-add technique to produce the normal-exposure-time images (minutes). The high precision photometry was performed via difference image analysis employing the DanDIA pipeline. We attempted automatic detection of variable stars in the field. We statistically analysed the light curves of 1405 stars in the crowded central region of NGC~6715 to automatically identify the variable stars present in this cluster. We found light curves for 17 previously known variable stars near the edges of our reference image (16 RR Lyrae and 1 semi-regular) and we discovered 67 new variables (30 RR Lyrae, 21 long-period irregular, 3 semi-regular, 1 W Virginis, 1 eclipsing binary, and 11 unclassified). Photometric measurements for these stars are available in electronic form through the Strasbourg Astronomical Data Centre., Comment: Accepted for publication in Astronomy & Astrophysics, 18 pages, 7 figures, 4 tables

(abridged) Context. Wide binaries are a potential pathway for the formation of hot Jupiters. The binary fraction among host stars is an important discriminator between competing formation theories, but has not been well characterised. Additionally, contaminating light from unresolved stars can significantly affect the accuracy of photometric and spectroscopic measurements in studies of transiting exoplanets. Aims. We observed 101 transiting exoplanet host systems in the Southern hemisphere in order to create a homogeneous catalogue of both bound companion stars and contaminating background stars. We investigate the binary fraction among the host stars in order to test theories for the formation of hot Jupiters, in an area of the sky where transiting exoplanetary systems have not been systematically searched for stellar companions. Methods. Lucky imaging observations from the Two Colour Instrument on the Danish 1.54m telescope at La Silla were used to search for previously unresolved stars at small angular separations. The separations and relative magnitudes of all detected stars were measured. For 12 candidate companions to 10 host stars, previous astrometric measurements were used to evaluate how likely the companions are to be physically associated. Results. We provide measurements of 499 candidate companions within 20 arcseconds of our sample of 101 planet host stars. 51 candidates are located within 5 arcseconds of a host star, and we provide the first published measurements for 27 of these. Calibrations for the plate scale and colour performance of the Two Colour Instrument are presented. Conclusions. We find that the overall multiplicity rate of the host stars is 38 +17 -13%, consistent with the rate among solar-type stars in our sensitivity range, suggesting that planet formation does not preferentially occur in long period binaries compared to a random sample of field stars., Comment: 19 pages, 9 figures. Accepted in A&A. Minor correction

Spitzer microlensing parallax observations of OGLE-2015-BLG-1212 decisively breaks a degeneracy between planetary and binary solutions that is somewhat ambiguous when only ground-based data are considered. Only eight viable models survive out of an initial set of 32 local minima in the parameter space. These models clearly indicate that the lens is a stellar binary system possibly located within the bulge of our Galaxy, ruling out the planetary alternative. We argue that several types of discrete degeneracies can be broken via such space-based parallax observations., Comment: 27 pages, 6 figures, accepted by ApJ

$K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $\pi_{\rm E}$ for $\gtrsim$127 microlensing events. These will include several that are planetary in nature as well as many short-timescale microlensing events, which are potentially indicative of free-floating planets (FFPs). These satellite parallax measurements will in turn allow for the direct measurement of the masses of and distances to the lensing systems. In this white paper we provide an overview of the $K2$C9 space- and ground-based microlensing survey. Specifically, we detail the demographic questions that can be addressed by this program, including the frequency of FFPs and the Galactic distribution of exoplanets, the observational parameters of $K2$C9, and the array of resources dedicated to concurrent observations. Finally, we outline the avenues through which the larger community can become involved, and generally encourage participation in $K2$C9, which constitutes an important pathfinding mission and community exercise in anticipation of $WFIRST$., Comment: 19 pages, 11 figures, 3 tables; submitted to PASP

The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 ($\mu_{\rm cl}$(N,E) = (+0.36+-0.10, +1.42+-0.10) mas/yr) as well as the source and show that the lens is not a cluster member. Even though this particular event does not probe the distribution of planets in globular clusters, other potential cluster lens events can be verified using our methodology. Additionally, we find that microlens parallax measured using OGLE photometry is consistent with the value found based on the light curve displacement between Earth and Spitzer., Comment: ApJ submitted, 12 pages, 8 figures, 2 tabels

Obtain time-series photometry of the very crowded central regions of Galactic globular clusters with better angular resolution than previously achieved with conventional CCDs on ground-based telescopes to complete, or improve, the census of the variable star population in those stellar systems. Images were taken using the Danish 1.54-m Telescope at the ESO observatory at La Silla in Chile. The telescope was equipped with an electron-multiplying CCD and the short-exposure-time images obtained (10 images per second) were stacked using the shift-and-add technique to produce the normal-exposure-time images (minutes). Photometry was performed via difference image analysis. Automatic detection of variable stars in the field was attempted. The light curves of 12541 stars in the cores of 10 globular clusters were statistically analysed in order to automatically extract the variable stars. We obtained light curves for 31 previously known variable stars (3 L, 2 SR, 20 RR Lyrae, 1 SX Phe, 3 cataclysmic variables, 1 EW and 1 NC) and we discovered 30 new variables (16 L, 7 SR, 4 RR Lyrae, 1 SX Phe and 2 NC)., Comment: 23 pages, 39 figures and 13 tables. Accepted for publication in the Astronomy & Astrophysics Journal (A&A)

We present 13 high-precision and four additional light curves of four bright southern-hemisphere transiting planetary systems: WASP-22, WASP-41, WASP-42 and WASP-55. In the cases of WASP-42 and WASP-55, these are the first follow-up observations since their discovery papers. We present refined measurements of the physical properties and orbital ephemerides of all four systems. No indications of transit timing variations were seen. All four planets have radii inflated above those expected from theoretical models of gas-giant planets; WASP-55b is the most discrepant with a mass of 0.63 Mjup and a radius of 1.34 Rjup. WASP-41 shows brightness anomalies during transit due to the planet occulting spots on the stellar surface. Two anomalies observed 3.1 d apart are very likely due to the same spot. We measure its change in position and determine a rotation period for the host star of 18.6 +/- 1.5 d, in good agreement with a published measurement from spot-induced brightness modulation, and a sky-projected orbital obliquity of lambda = 6 +/- 11 degrees. We conclude with a compilation of obliquity measurements from spot-tracking analyses and a discussion of this technique in the study of the orbital configurations of hot Jupiters., Comment: 13 pages, 8 tables, 11 figures. Version 2 is the final accepted version of the paper

Accurate measurements of the physical characteristics of a large number of exoplanets are useful to strongly constrain theoretical models of planet formation and evolution, which lead to the large variety of exoplanets and planetary-system configurations that have been observed. We present a study of the planetary systems WASP-45 and WASP-46, both composed of a main-sequence star and a close-in hot Jupiter, based on 29 new high-quality light curves of transits events. In particular, one transit of WASP-45 b and four of WASP-46 b were simultaneously observed in four optical filters, while one transit of WASP-46 b was observed with the NTT obtaining precision of 0.30 mmag with a cadence of roughly three minutes. We also obtained five new spectra of WASP-45 with the FEROS spectrograph. We improved by a factor of four the measurement of the radius of the planet WASP-45 b, and found that WASP-46 b is slightly less massive and smaller than previously reported. Both planets now have a more accurate measurement of the density (0.959 +\- 0.077 \rho Jup instead of 0.64 +\- 0.30 \rho Jup for WASP-45 b, and 1.103 +\- 0.052 \rho Jup instead of 0.94 +\- 0.11 \rho Jup for WASP-46 b). We tentatively detected radius variations with wavelength for both planets, in particular in the case of WASP-45 b we found a slightly larger absorption in the redder bands than in the bluer ones. No hints for the presence of an additional planetary companion in the two systems were found either from the photometric or radial velocity measurements., Comment: 14 pages, 10 figures. Accepted for publication in MNRAS

Context. Photometric monitoring of the variability of brown dwarfs can provide useful information about the structure of clouds in their cold atmospheres. The brown-dwarf binary system Luhman 16AB is an interesting target for such a study, as its components stand at the L/T transition and show high levels of variability. Luhman 16AB is also the third closest system to the Solar system, allowing precise astrometric investigations with ground-based facilities. Aims. The aim of the work is to estimate the rotation period and study the astrometric motion of both components. Methods. We have monitored Luhman 16AB over a period of two years with the lucky-imaging camera mounted on the Danish 1.54m telescope at La Silla, through a special i+z long-pass filter, which allowed us to clearly resolve the two brown dwarfs into single objects. An intense monitoring of the target was also performed over 16 nights, in which we observed a peak-to-peak variability of 0.20 \pm 0.02 mag and 0.34 \pm 0.02 mag for Luhman 16A and 16B, respectively. Results. We used the 16-night time-series data to estimate the rotation period of the two components. We found that Luhman 16B rotates with a period of 5.1 \pm 0.1 hr, in very good agreement with previous measurements. For Luhman 16A, we report that it rotates slower than its companion and, even though we were not able to get a robust determination, our data indicate a rotation period of roughly 8 hr. This implies that the rotation axes of the two components are well aligned and suggests a scenario in which the two objects underwent the same accretion process. The 2-year complete dataset was used to study the astrometric motion of Luhman 16AB. We predict a motion of the system that is not consistent with a previous estimate based on two months of monitoring, but cannot confirm or refute the presence of additional planetary-mass bodies in the system., Comment: 10 pages, 7 figures, accepted for publication by Astronomy & Astrophysics

For all exoplanet candidates, the reliability of a claimed detection needs to be assessed through a careful study of systematic errors in the data to minimize the false positives rate. We present a method to investigate such systematics in microlensing datasets using the microlensing event OGLE-2013-BLG-0446 as a case study. The event was observed from multiple sites around the world and its high magnification (A_{max} \sim 3000) allowed us to investigate the effects of terrestrial and annual parallax. Real-time modeling of the event while it was still ongoing suggested the presence of an extremely low-mass companion (\sim 3M_\oplus ) to the lensing star, leading to substantial follow-up coverage of the light curve. We test and compare different models for the light curve and conclude that the data do not favour the planetary interpretation when systematic errors are taken into account., Comment: accepted ApJ 2015

We report on the mass and distance measurements of two single-lens events from the 2015 \emph{Spitzer} microlensing campaign. With both finite-source effect and microlens parallax measurements, we find that the lens of OGLE-2015-BLG-1268 is very likely a brown dwarf. Assuming that the source star lies behind the same amount of dust as the Bulge red clump, we find the lens is a $45\pm7$ $M_{\rm J}$ brown dwarf at $5.9\pm1.0$ kpc. The lens of of the second event, OGLE-2015-BLG-0763, is a $0.50\pm0.04$ $M_\odot$ star at $6.9\pm1.0$ kpc. We show that the probability to definitively measure the mass of isolated microlenses is dramatically increased once simultaneous ground- and space-based observations are conducted., Comment: 10 papers, 4 figures, 2 tables; ApJ in press

Transits in the WASP-57 planetary system have been found to occur half an hour earlier than expected. We present ten transit light curves from amateur telescopes, on which this discovery was based, thirteen transit light curves from professional facilities which confirm and refine this finding, and high-resolution imaging which show no evidence for nearby companions. We use these data to determine a new and precise orbital ephemeris, and measure the physical properties of the system. Our revised orbital period is 4.5s shorter than found from the discovery data alone, which explains the early occurrence of the transits. We also find both the star and planet to be larger and less massive than previously thought. The measured mass and radius of the planet are now consistent with theoretical models of gas giants containing no heavy-element core, as expected for the sub-solar metallicity of the host star. Two transits were observed simultaneously in four passbands. We use the resulting light curves to measure the planet's radius as a function of wavelength, finding that our data are sufficient in principle but not in practise to constrain its atmospheric properties. We conclude with a discussion of the current and future status of transmission photometry studies for probing the atmospheres of gas-giant transiting planets., Comment: Accepted for publication in MNRAS. 13 pages, 8 figures, 7 tables. The data will be available at the first author's website and the system will be added to the TEPCat catalogue at http://www.astro.keele.ac.uk/jkt/tepcat/

We report the detection of a Cold Neptune m_planet=21+/-2MEarth orbiting a 0.38MSol M dwarf lying 2.5-3.3 kpc toward the Galactic center as part of a campaign combining ground-based and Spitzer observations to measure the Galactic distribution of planets. This is the first time that the complex real-time protocols described by Yee et al. (2015), which aim to maximize planet sensitivity while maintaining sample integrity, have been carried out in practice. Multiple survey and follow-up teams successfully combined their efforts within the framework of these protocols to detect this planet. This is the second planet in the Spitzer Galactic distribution sample. Both are in the near-to-mid disk and clearly not in the Galactic bulge., Comment: 28 pages, 3 figures, 2 tables, submitted to ApJ

We analyse 20 nights of CCD observations in the V and I bands of the globular cluster M68 (NGC 4590), using these to detect variable objects. We also obtained electron-multiplying CCD (EMCCD) observations for this cluster in order to explore its core with unprecedented spatial resolution from the ground. We reduced our data using difference image analysis, in order to achieve the best possible photometry in the crowded field of the cluster. In doing so, we showed that when dealing with identical networked telescopes, a reference image from any telescope may be used to reduce data from any other telescope, which facilitates the analysis significantly. We then used our light curves to estimate the properties of the RR Lyrae (RRL) stars in M68 through Fourier decomposition and empirical relations. The variable star properties then allowed us to derive the cluster's metallicity and distance. We determine new periods for the variable stars, and search for new variables, especially in the core of the cluster where our method performs particularly well. We detect an additional 4 SX Phe stars, and confirm the variability of another star, bringing the total number of confirmed variable stars in this cluster to 50. We also used archival data stretching back to 1951 in order to derive period changes for some of the single-mode RRL stars, and analyse the significant number of double-mode RRL stars in M68. Furthermore, we find evidence for double-mode pulsation in one of the SX Phe stars in this cluster. Using the different types of variables, we derived an estimate of the metallicity, [Fe/H]=$-2.07 \pm 0.06$ on the ZW scale, and 4 independent estimates of the distance modulus ($\mu_0 \sim$ 15.00 mag) for this cluster. Thanks to the first use of difference image analysis on time-series observations of M68, we are now confident that we have a complete census of the RRL stars in this cluster., Comment: 25 pages, 18 figures, A&A in press

We present microlens parallax measurements for 21 (apparently) isolated lenses observed toward the Galactic bulge that were imaged simultaneously from Earth and Spitzer, which was ~1 AU West of Earth in projection. We combine these measurements with a kinematic model of the Galaxy to derive distance estimates for each lens, with error bars that are small compared to the Sun's Galactocentric distance. The ensemble therefore yields a well-defined cumulative distribution of lens distances. In principle it is possible to compare this distribution against a set of planets detected in the same experiment in order to measure the Galactic distribution of planets. Since these Spitzer observations yielded only one planet, this is not yet possible in practice. However, it will become possible as larger samples are accumulated., Comment: accepted for publication in ApJ

We present 17 transit light curves of the ultra-short period planetary system WASP-103, a strong candidate for the detection of tidally-induced orbital decay. We use these to establish a high-precision reference epoch for transit timing studies. The time of the reference transit midpoint is now measured to an accuracy of 4.8s, versus 67.4s in the discovery paper, aiding future searches for orbital decay. With the help of published spectroscopic measurements and theoretical stellar models, we determine the physical properties of the system to high precision and present a detailed error budget for these calculations. The planet has a Roche lobe filling factor of 0.58, leading to a significant asphericity; we correct its measured mass and mean density for this phenomenon. A high-resolution Lucky Imaging observation shows no evidence for faint stars close enough to contaminate the point spread function of WASP-103. Our data were obtained in the Bessell $RI$ and the SDSS $griz$ passbands and yield a larger planet radius at bluer optical wavelengths, to a confidence level of 7.3 sigma. Interpreting this as an effect of Rayleigh scattering in the planetary atmosphere leads to a measurement of the planetary mass which is too small by a factor of five, implying that Rayleigh scattering is not the main cause of the variation of radius with wavelength., Comment: Accepted for publication in MNRAS. 12 pages, 8 tables, 7 figures. The reduced data will be available at http://www.astro.keele.ac.uk/~jkt/data-teps.html

In this paper, we present the analysis of time-series observations from 2013 and 2014 of five metal rich ([Fe/H] $>$ -1) globular clusters: NGC~6388, NGC~6441, NGC~6528, NGC~6638, and NGC~6652. The data have been used to perform a census of the variable stars in the central parts of these clusters. The observations were made with the electron multiplying CCD (EMCCD) camera at the Danish 1.54m Telescope at La Silla, Chile, and they were analysed using difference image analysis (DIA) to obtain high-precision light curves of the variable stars. It was possible to identify and classify all of the previously known or suspected variable stars in the central regions of the five clusters. Furthermore, we were able to identify, and in most cases classify 48, 49, 7, 8, and 2 previously unknown variables in NGC~6388, NGC~6441, NGC~6528, NGC~6638, and NGC~6652, respectively. Especially interesting is the case of NGC~6441, for which the variable star population of about 150 stars has been thoroughly examined by previous studies, including a Hubble Space Telescope study. In this paper we are able to present 49 new variable stars for this cluster, of which one (possibly two) are RR Lyrae stars, two are W Virginis stars, and the rest are long period semi-regular/irregular variables on the red giant branch. We have also detected the first double mode RR Lyrae in the cluster., Comment: Accepted to A&A, 24 pages, 23 figures, 11 tables. Updated to match version published in A&A

We present time-series photometric observations of thirteen transits in the planetary systems WASP-24, WASP-25 and WASP-26. All three systems have orbital obliquity measurements, WASP-24 and WASP-26 have been observed with Spitzer, and WASP-25 was previously comparatively neglected. Our light curves were obtained using the telescope-defocussing method and have scatters of 0.5 to 1.2 mmag relative to their best-fitting geometric models. We used these data to measure the physical properties and orbital ephemerides of the systems to high precision, finding that our improved measurements are in good agreement with previous studies. High-resolution Lucky Imaging observations of all three targets show no evidence for faint stars close enough to contaminate our photometry. We confirm the eclipsing nature of the star closest to WASP-24 and present the detection of a detached eclipsing binary within 4.25 arcmin of WASP-26., Comment: Accepted for publication in MNRAS. 14 pages, 10 figures, 8 tables. Data and supplementary information are available on request

The extrasolar planet WASP-67 b is the first hot Jupiter definitively known to undergo only partial eclipses. The lack of the second and third contact point in this planetary system makes it difficult to obtain accurate measurements of its physical parameters. Aims. By using new high-precision photometric data, we confirm that WASP-67 b shows grazing eclipses and compute accurate estimates of the physical properties of the planet and its parent star. Methods. We present high-quality, multi-colour, broad-band photometric observations comprising five light curves covering two transit events, obtained using two medium-class telescopes and the telescope-defocussing technique. One transit was observed through a Bessel-R filter and the other simultaneously through filters similar to Sloan griz. We modelled these data using jktebop. The physical parameters of the system were obtained from the analysis of these light curves and from published spectroscopic measurements. Results. All five of our light curves satisfy the criterion for being grazing eclipses. We revise the physical parameters of the whole WASP-67 system and, in particular, significantly improve the measurements of the planet's radius and density as compared to the values in the discovery paper. The transit ephemeris was also substantially refined. We investigated the variation of the planet's radius as a function of the wavelength, using the simultaneous multi-band data, finding that our measurements are consistent with a flat spectrum to within the experimental uncertainties., Comment: 10 pages, 8 figures, 5 tables, to appear in Astronomy & Astrophysics

WASP-80 is one of only two systems known to contain a hot Jupiter which transits its M-dwarf host star. We present eight light curves of one transit event, obtained simultaneously using two defocussed telescopes. These data were taken through the Bessell I, Sloan griz and near-infrared JHK passbands. We use our data to search for opacity-induced changes in the planetary radius, but find that all values agree with each other. Our data are therefore consistent with a flat transmission spectrum to within the observational uncertainties. We also measure an activity index of the host star of log R'_HK=-4.495, meaning that WASP-80A shows strong chromospheric activity. The non-detection of starspots implies that, if they exist, they must be small and symmetrically distributed on the stellar surface. We model all available optical transit light curves to obtain improved physical properties and orbital ephemerides for the system., Comment: 9 pages, 6 figures

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q 0.03, which differs from that reported in Shvartzvald et al. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s ~ 0.15 or s ~ 7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M dwarf orbited by a massive Jupiter companion at very close (M(host) = 0.30(sup +0.08, sub -0.06)Mʘ, M(comp) = 10.1 (sup +2.9, sub -2.2 M(jup), a(exp) = 0.40 (sup +0.05, sub -0.40)au) or wide (M(host) = 0.28 (sup +0.10, sub -0.80)Mʘ, M(comp) = 9.9 (sup +3.8, sub -3.5)M(jup), a(exp) = 18.0 (sup +3.2, sub -3.2)au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q = 0.03, which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M(Jup). It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M(Jup) in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M dwarfs and even BDs, which is challenging for other exoplanet detection methods.