Search

We report the results of a $J$-band survey for photometric variability in a sample of young, low-gravity objects using the New Technology Telescope (NTT) and the United Kingdom InfraRed Telescope (UKIRT). Surface gravity is a key parameter in the atmospheric properties of brown dwarfs and this is the first large survey that aims to test the gravity dependence of variability properties. We do a full analysis of the spectral signatures of youth and assess the group membership probability of each target using membership tools from the literature. This results in a 30 object sample of young low-gravity brown dwarfs. Since we are lacking in objects with spectral types later than L9, we focus our statistical analysis on the L0-L8.5 objects. We find that the variability occurrence rate of L0-L8.5 low-gravity brown dwarfs in this survey is $30^{+16}_{-8}\%$. We reanalyse the results of Radigan 2014 and find that the field dwarfs with spectral types L0-L8.5 have a variability occurrence rate of $11^{+13}_{-4}\%$. We determine a probability of $98\%$ that the samples are drawn from different distributions. This is the first quantitative indication that the low-gravity objects are more likely to be variable than the field dwarf population. Furthermore, we present follow-up $J_S$ and $K_S$ observations of the young, planetary-mass variable object PSO 318.5-22 over three consecutive nights. We find no evidence of phase shifts between the $J_S$ and $K_S$ bands and find higher $J_S$ amplitudes. We use the $J_S$ lightcurves to measure a rotational period of $8.45\pm0.05~$hr for PSO 318.5-22., Comment: accepted for publication in MNRAS

We report direct observations of photon-mediated spin-exchange interactions in an atomic ensemble. Interactions extending over a distance of 500 microns are generated within a cloud of cold rubidium atoms coupled to a single mode of light in an optical resonator. We characterize the system via quench dynamics and imaging of the local magnetization, verifying the coherence of the interactions and demonstrating optical control of their strength and sign. Furthermore, by initializing the spin-1 system in the mF = 0 Zeeman state, we observe correlated pair creation in the mF = +/- 1 states, a process analogous to spontaneous parametric down-conversion and to spin mixing in Bose-Einstein condensates. Our work opens new opportunities in quantum simulation with long-range interactions and in entanglement-enhanced metrology.

We present 27 new L subdwarfs and classify five of them as esdL and 22 as sdL. Our L subdwarf candidates were selected with the UKIRT Infrared Deep Sky Survey and Sloan Digital Sky Survey. Spectroscopic follow-up was carried out primarily with the OSIRIS spectrograph on the Gran Telescopio Canarias. Some of these new objects were followed up with the X-shooter instrument on the Very Large Telescope. We studied the photometric properties of the population of known L subdwarfs using colour-spectral type diagrams and colour-colour diagrams, by comparison with L dwarfs and main-sequence stars, and identified new colour spaces for L subdwarf selection/study in current and future surveys. We further discussed the brown dwarf transition-zone and the observational stellar/substellar boundary. We found that about one-third of 66 known L subdwarfs are substellar objects, with two-thirds being very low-mass stars. We also present the Hertzsprung-Russell diagrams, spectral type-absolute magnitude corrections, and tangential velocities of 20 known L subdwarfs observed by the Gaia astrometry satellite. One of our L subdwarf candidates, ULAS J233227.03+123452.0, is a mildly metal-poor spectroscopic binary brown dwarf: a ~L6p dwarf and a ~T4p dwarf. This binary is likely a thick disc member according to its kinematics., Comment: 30 pages, 23 figures

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model., Comment: 33 pages, 25 figures. Accepted in Astronomy and Astrophysics

Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars., Comment: Accepted for publication in A&A. 26 pages, 22 figures, 12 tables

We report the discovery of an esdL3 subdwarf, ULAS J020858.62+020657.0, and a usdL4.5 subdwarf, ULAS J230711.01+014447.1. They were identified as L subdwarfs by optical spectra obtained with the Gran Telescopio Canarias, and followed up by optical-to-near-infrared spectroscopy with the Very Large Telescope. We also obtained an optical-to-near-infrared spectrum of a previously known L subdwarf, ULAS J135058.85+081506.8, and reclassified it as a usdL3 subdwarf. These three objects all have typical halo kinematics. They have $T_{\rm eff}$ around 2050$-$2250 K, $-$1.8 $\leq$ [Fe/H] $\leq -$1.5, and mass around 0.0822$-$0.0833 M$_{\odot}$, according to model spectral fitting and evolutionary models. These sources are likely halo transitional brown dwarfs with unsteady hydrogen fusion, as their masses are just below the hydrogen-burning minimum mass, which is $\sim$ 0.0845 M$_{\odot}$ at [Fe/H] = $-$1.6 and $\sim$ 0.0855 M$_{\odot}$ at [Fe/H] = $-$1.8. Including these, there are now nine objects in the `halo brown dwarf transition zone', which is a `substellar subdwarf gap' that spans a wide temperature range within a narrow mass range of the substellar population., Comment: 9 pages, 5 figures

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project., Comment: Minor changes to author list and title. Comments and feedback welcome through the paper source repository: https://github.com/astropy/astropy-v2.0-paper For more information about Astropy, see http://www.astropy.org/

We present simultaneous HST WFC3 + Spitzer IRAC variability monitoring for the highly-variable young ($\sim$20 Myr) planetary-mass object PSO J318.5-22. Our simultaneous HST + Spitzer observations covered $\sim$2 rotation periods with Spitzer and most of a rotation period with HST. We derive a period of 8.6$\pm$0.1 hours from the Spitzer lightcurve. Combining this period with the measured $v sin i$ for this object, we find an inclination of 56.2$\pm 8.1^{\circ}$. We measure peak-to-trough variability amplitudes of 3.4$\pm$0.1$\%$ for Spitzer Channel 2 and 4.4 - 5.8$\%$ (typical 68$\%$ confidence errors of $\sim$0.3$\%$) in the near-IR bands (1.07-1.67 $\mu$m) covered by the WFC3 G141 prism -- the mid-IR variability amplitude for PSO J318.5-22 one of the highest variability amplitudes measured in the mid-IR for any brown dwarf or planetary mass object. Additionally, we detect phase offsets ranging from 200--210$^{\circ}$ (typical error of $\sim$4$^{\circ}$) between synthesized near-IR lightcurves and the Spitzer mid-IR lightcurve, likely indicating depth-dependent longitudinal atmospheric structure in this atmosphere. The detection of similar variability amplitudes in wide spectral bands relative to absorption features suggests that the driver of the variability may be inhomogeneous clouds (perhaps a patchy haze layer over thick clouds), as opposed to hot spots or compositional inhomogeneities at the top-of-atmosphere level., Comment: 48 pages, 22 figures, accepted to AJ

The atmospheres of cool stars are temporally and spatially inhomogeneous due to the effects of convection.The influence of this inhomogeneity, referred to as granulation, on colours has never been investigated over a large range of effective temperatures and gravities. We aim to study, in a quantitative way, the impact of granulation on colours. We use the CIFIST grid of CO5BOLD hydrodynamical models to compute emerging fluxes. These in turn are used to compute theoretical colours in the UBVRI,2MASS,Hipparcos,Gaia and SDSS systems. Every CO5BOLD model has a corresponding one dimensional (1D) plane-parallel LHD model computed for the same atmospheric parameters, which we used to define a "3D correction" that can be applied to colours computed from fluxes computed from any 1D model atmosphere code. The 3D corrections on colours are generally small, of the order of a few hundredths of a magnitude, yet they are far from negligible. We find that ignoring granulation effects can lead to underestimation of Teff by up to 200K and overestimation of gravity by up to 0.5dex, when using colours as diagnostics. We have identified a major shortcoming in how scattering is treated in the current version of the CIFIST grid, which could lead to offsets of the order 0.01mag, especially for colours involving blue and UV bands.We have investigated the Gaia and Hipparcos photometric systems and found that the (G-H_p),(BP-RP) diagram is immune to the effects of granulation. In addition, we point to the potential of the RVS photometry as a metallicity diagnostic. Our investigation shows that the effects of granulation should not be neglected if one wants to use colours as diagnostics of the stellar parameters of F,G,K stars. A limitation is that scattering is treated as true absorption in our current computations, thus our 3D corrections are likely an upper limit to the true effect. (Abridged), Comment: Accepted for publication by A&A, typo in the title corrected, cross reference to Paper II

The ability to perform detailed chemical analysis of Sun-like F-, G-, and K-type stars is a powerful tool with many applications including studying the chemical evolution of the Galaxy and constraining planet formation theories. Unfortunately, complications in modeling cooler stellar atmospheres hinders similar analysis of M-dwarf stars. Empirically-calibrated methods to measure M dwarf metallicity from moderate-resolution spectra are currently limited to measuring overall metallicity and rely on astrophysical abundance correlations in stellar populations. We present a new, empirical calibration of synthetic M dwarf spectra that can be used to infer effective temperature, Fe abundance, and Ti abundance. We obtained high-resolution (R~25,000), Y-band (~1 micron) spectra of 29 M dwarfs with NIRSPEC on Keck II. Using the PHOENIX stellar atmosphere modeling code (version 15.5), we generated a grid of synthetic spectra covering a range of temperatures, metallicities, and alpha-enhancements. From our observed and synthetic spectra, we measured the equivalent widths of multiple Fe I and Ti I lines and a temperature-sensitive index based on the FeH bandhead. We used abundances measured from widely-separated solar-type companions to empirically calibrate transformations to the observed indices and equivalent widths that force agreement with the models. Our calibration achieves precisions in Teff, [Fe/H], and [Ti/Fe] of 60 K, 0.1 dex, and 0.05 dex, respectively and is calibrated for 3200 K < Teff < 4100 K, -0.7 < [Fe/H] < +0.3, and -0.05 < [Ti/Fe] < +0.3. This work is a step toward detailed chemical analysis of M dwarfs at a similar precision achieved for FGK stars., Comment: accepted for publication in ApJ, all synthetic spectra available at http://people.bu.edu/mveyette/phoenix/

M dwarfs are an important source of information when studying and probing the lower end of the Hertzsprung-Russell (HR) diagram, down to the hydrogen-burning limit. Being the most numerous and oldest stars in the galaxy, they carry fundamental information on its chemical history. The presence of molecules in their atmospheres, along with various condensed species, complicates our understanding of their physical properties and thus makes the determination of their fundamental stellar parameters more challenging and difficult. The aim of this study is to perform a detailed spectroscopic analysis of the high-resolution H-band spectra of M dwarfs in order to determine their fundamental stellar parameters and to validate atmospheric models. The present study will also help us to understand various processes, including dust formation and depletion of metals onto dust grains in M dwarf atmospheres. The high spectral resolution also provides a unique opportunity to constrain other chemical and physical processes that occur in a cool atmosphere The high-resolution APOGEE spectra of M dwarfs, covering the entire H-band, provide a unique opportunity to measure their fundamental parameters. We have performed a detailed spectral synthesis by comparing these high-resolution H-band spectra to that of the most recent BT-settl model and have obtained fundamental parameters such as effective temperature, surface gravity, and metallicity (Teff, log g and [Fe/H]) respectively., Comment: 15 pages, 10 figures, accepted for publication in A&A

Characterization of the atmospheres of transiting exoplanets relies on accurate measurements of the extent of the optically thick area of the planet at multiple wavelengths with a precision $\lesssim$100 parts per million (ppm). Next-generation instruments onboard the James Webb Space Telescope (JWST) are expected to achieve $\sim$10 ppm precision for several tens of targets. A similar precision can be obtained in modelling only if other astrophysical effects, including the stellar limb-darkening, are accounted for properly. In this paper, we explore the limits on precision due to the mathematical formulas currently adopted to approximate the stellar limb-darkening, and to the use of limb-darkening coefficients obtained either from stellar-atmosphere models or empirically. We propose a new limb-darkening law with two coefficients, `power-2', which outperforms other two-coefficient laws adopted in the literature in most cases, and particularly for cool stars. Empirical limb-darkening based on two-coefficient formulas can be significantly biased, even if the light-curve residuals are nearly photon-noise limited. We demonstrate an optimal strategy to fitting for the four-coefficients limb-darkening in the visible, using prior information on the exoplanet orbital parameters to break some of the degeneracies that otherwise would prevent the convergence of the fit. Infrared observations taken with the James Webb Space Telescope (JWST) will provide accurate measurements of the exoplanet orbital parameters with unprecedented precision, which can be used as priors to improve the stellar limb-darkening characterization, and therefore the inferred exoplanet parameters, from observations in the visible, such as those taken with Kepler/K2, JWST, other past and future instruments.

51 Eridani b is an exoplanet around a young (20 Myr) nearby (29.4 pc) F0-type star, recently discovered by direct imaging. Being only 0.5" away from its host star it is well suited for spectroscopic analysis using integral field spectrographs. We aim to refine the atmospheric properties of this and to further constrain the architecture of the system by searching for additional companions. Using the SPHERE instrument at the VLT we extend the spectral coverage of the planet to the complete Y- to H-band range and provide photometry in the K12-bands (2.11, 2.25 micron). The object is compared to other cool and peculiar dwarfs. Furthermore, the posterior probability distributions of cloudy and clear atmospheric models are explored using MCMC. We verified our methods by determining atmospheric parameters for the two benchmark brown dwarfs Gl 570D and HD 3651B. For probing the innermost region for additional companions, archival VLT-NACO (L') SAM data is used. We present the first spectrophotometric measurements in the Y- and K-bands for the planet and revise its J-band flux to values 40% fainter than previous measurements. Cloudy models with uniform cloud coverage provide a good match to the data. We derive the temperature, radius, surface gravity, metallicity and cloud sedimentation parameter f_sed. We find that the atmosphere is highly super-solar (Fe/H~1.0) with an extended, thick cloud cover of small particles. The model radius and surface gravity suggest planetary masses of about 9 M_jup. The evolutionary model only provides a lower mass limit of >2 M_jup (for pure hot-start). The cold-start model cannot explain the planet's luminosity. The SPHERE and NACO/SAM detection limits probe the 51 Eri system at Solar System scales and exclude brown-dwarf companions more massive than 20 M_jup beyond separations of ~2.5 au and giant planets more massive than 2 M_jup beyond 9 au., Comment: 29 pages, 31 figures, accepted for publication in A&A

We conducted a multi-wavelength, multi-instrument observational characterisation of the candidate free-floating planet CFBDSIR~J214947.2-040308.9, a late T-dwarf with possible low-gravity features, in order to constrain its physical properties. We analyzed 9 hours of X-Shooter spectroscopy with signal detectable from 0.8--2.3$\mu$m, as well as additional photometry in the mid-infrared using the Spitzer Space Telescope. Combined with a VLT/HAWK-I astrometric parallax, this enabled a full characterisation of the absolute flux from the visible to 5$\mu$m, encompassing more than 90\% of the expected energy emitted by such a cool late T-type object. Our analysis of the spectrum also provided the radial velocity and therefore the determination of its full 3-D kinematics. While our new spectrum confirms the low gravity and/or high metallicity of CFBDSIR2149, the parallax and kinematics safely rule out membership to any known young moving group, including AB~Doradus. We use the equivalent width of the KI doublet at 1.25$\mu$m as a promising tool to discriminate the effects of low-gravity from the effects of high-metallicity on the emission spectra of cool atmospheres. In the case of CFBDSIR2149, the observed KI doublet clearly favours the low-gravity solution. CFBDSIR2149 is therefore a peculiar late-T dwarf that is probably a young, planetary-mass object (2--13Mjup, $<$500Myr) possibly similar to the exoplanet 51Erib, or perhaps a 2--40Mjup brown dwarf with super-solar metallicity., Comment: Accepted in A&A

SDSS J010448.46+153501.8 has previously been classified as an sdM9.5 subdwarf. However, its very blue $J-K$ colour ($-0.15 \pm 0.17$) suggests a much lower metallicity compared to normal sdM9.5 subdwarfs. Here, we re-classify this object as a usdL1.5 subdwarf based on a new optical and near-infrared spectrum obtained with X-shooter on the Very Large Telescope. Spectral fitting with BT-Settl models leads to $T_{\rm eff}$ = 2450 $\pm$ 150 K, [Fe/H] = $-$2.4 $\pm$ 0.2 and log $g$ = 5.5 $\pm$ 0.25. We estimate a mass for SDSS J010448.46+153501.8 of 0.086 $\pm$ 0.0015 M$_{\odot}$ which is just below the hydrogen-burning minimum mass at [Fe/H] = $-$2.4 ($\sim$0.088 M$_{\odot}$) according to evolutionary models. Our analysis thus shows SDSS J0104+15 to be the most metal-poor and highest mass substellar object known to-date. We found that SDSS J010448.46+153501.8 is joined by another five known L subdwarfs (2MASS J05325346+8246465, 2MASS J06164006$-$6407194, SDSS J125637.16$-$022452.2, ULAS J151913.03$-$000030.0 and 2MASS J16262034+3925190) in a 'halo brown dwarf transition zone' in the $T_{\rm eff}-$[Fe/H] plane, which represents a narrow mass range in which unsteady nuclear fusion occurs. This halo brown dwarf transition zone forms a 'substellar subdwarf gap' for mid L to early T types., Comment: 11 pages, 9 figures

Context. Tighter constraints on metal-poor stars we observe are needed to better understand the chemical processes of the early Universe. Computing a stellar spectrum in 3D allows one to model complex stellar behaviours, which cannot be replicated in 1D. Aims. We examine the effect that the intrinsic CNO abundances have on a 3D model structure and the resulting 3D spectrum synthesis. Methods. Model atmospheres were computed in 3D for three distinct CNO chemical compositions using the CO 5 BOLD model atmosphere code, and their internal structures were examined. Synthetic spectra were computed from these models using Linfor3D and they were compared. New 3D abundance corrections for the G-band and a selection of UV OH lines were also computed. Results. The varying CNO abundances change the metal content of the 3D models. This had an effect on the model structure and the resulting synthesis. However, it was found that the C/O ratio had a larger effect than the overall metal content of a model. Conclusions. Our results suggest that varying the C/O ratio has a substantial impact on the internal structure of the 3D model, even in the hot turn-off star models explored here. This suggests that bespoke 3D models, for specific CNO abundances should be sought. Such effects are not seen in 1D at these temperature regimes., Comment: Accepted for publication as a letter in A&A. 4 figures, 2 tables, 6 pages

The recent discovery of a diffuse cosmic neutrino flux extending up to PeV energies raises the question of which astrophysical sources generate this signal. One class of extragalactic sources which may produce such high-energy neutrinos are blazars. We present a likelihood analysis searching for cumulative neutrino emission from blazars in the 2nd Fermi-LAT AGN catalogue (2LAC) using an IceCube neutrino dataset 2009-12 which was optimised for the detection of individual sources. In contrast to previous searches with IceCube, the populations investigated contain up to hundreds of sources, the largest one being the entire blazar sample in the 2LAC catalogue. No significant excess is observed and upper limits for the cumulative flux from these populations are obtained. These constrain the maximum contribution of the 2LAC blazars to the observed astrophysical neutrino flux to be $27 \%$ or less between around 10 TeV and 2 PeV, assuming equipartition of flavours at Earth and a single power-law spectrum with a spectral index of $-2.5$. We can still exclude that the 2LAC blazars (and sub-populations) emit more than $50 \%$ of the observed neutrinos up to a spectral index as hard as $-2.2$ in the same energy range. Our result takes into account that the neutrino source count distribution is unknown, and it does not assume strict proportionality of the neutrino flux to the measured 2LAC $\gamma$-ray signal for each source. Additionally, we constrain recent models for neutrino emission by blazars., Comment: 18 pages, 22 figures