Source author record

Thomas Henning

Thomas Henning appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

259works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

259 published item(s)

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) III: The vertical structure of debris disks

Debris disks -- collisionally sustained belts of dust and sometimes gas around main sequence stars -- are remnants of planet formation processes and are found in systems ${\gtrsim}10$ Myr old. Millimeter-wavelength observations are particularly important, as the grains probed by these observations are not strongly affected by radiation pressure and stellar winds, allowing them to probe the dynamics of large bodies producing dust. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is analyzing high-resolution observations of 24 debris disks to enable the characterization of debris disk substructures across a large sample for the first time. For the most highly inclined disks, it is possible to recover the vertical structure of the disk. We aim to model and analyze the most highly inclined systems in the ARKS sample in order to uniformly extract the vertical dust distributions for a sample of well-resolved debris disks. We employed both parametric and nonparametric methods to constrain the vertical dust distributions for the most highly inclined ARKS targets. We find a broad range of aspect ratios, revealing a wide diversity in vertical structure, with a range of best-fit parametric values of $0.0026 \leq h_{\rm HWHM} \leq 0.193$ and a median best-fit value of $h_{\rm HWHM}=0.021$. The results obtained by nonparametric modeling are generally consistent with the parametric modeling results. We find that five of the 13 disks are consistent with having total disk masses less than that of Neptune (17 $M_{\oplus}$), assuming stirring by internal processes (self-stirring and collisional and frictional damping). Furthermore, most systems show a significant preference for a Lorentzian vertical profile rather than a Gaussian.

preprint2025arXiv

MIDIS: Quantifying the AGN component of X-ray-detected galaxies

We combine the deepest X-ray survey from the Chandra Deep Field-South (CDF-S) `7-Ms' survey with the deepest mid-infrared (5.6$ μm$) image from the JWST/MIRI Deep Imaging Survey (MIDIS) in the Hubble Ultra-Deep Field (HUDF) to study the infrared counterparts and point-source emission of 31 X-ray sources with a median, intrinsic, rest-frame X-ray luminosity of $\log_{10}(L_{\rm Xc}^{\rm 0.5-7keV})$=42.04$\pm$0.22 erg $\rm s^{-1}$. The sample includes 24 AGN with a redshift range, as set by the X-ray detectability, of $z \simeq 0.5-3$. Through a multi-wavelength morphological decomposition, employing three separate classifications (visual, parametric and non-parametric) we separate (where present) the luminosity of the point-like AGN component from the remainder of the host-galaxy emission. The unprecedented mid-infrared sensitivity and imaging resolution of MIRI allows, in many cases, the direct characterisation of point-like (i.e. unresolved) components in the galaxies' emission. We establish a broad agreement between the three morphological classifications. At least 70% of the X-ray sources, including some classified as galaxies, show unresolved emission in the MIRI images, with the unresolved-to-total flux fraction at rest-frame 2$μm$ ranging from $\sim$0.2 to $\sim$0.9. At high X-ray luminosities ($\log_{10}(L_{\rm Xc}$)>43 erg $\rm s^{-1}$) we derive a consistent rest-frame near-infrared 2$ μm$ point-source luminosity to that derived for local AGN, whilst at lower X-ray luminosity we identify an excess in the 2$ μm$ emission compared to pre-JWST studies. We speculate this offset may be driven by a combination of Compton-thick AGN components and nuclear starburst, merger driven activity. Our observations highlight the complex nature of X-ray sources in the distant Universe and demonstrate the power of JWST/MIRI in quantifying their nuclear infrared emission. (Abridged)

preprint2025arXiv

MINDS. Young binary systems with JWST/MIRI: Variable water-rich primaries and extended emission

As part of the JWST GTO program MINDS, we analyze the mid-infrared emission of three Class II binary systems: VW Cha, WX Cha, and RW Aur, to investigate the impact of stellar multiplicity on the chemistry and physics of their inner disk. We analyze the 1D spectrum from JWST/MIRI-MRS for primary and secondary disks separately, extracted by combining forward modeling with a theoretical PSF and aperture photometry. We modeled the molecular lines with 0D slab models. We interpret the results by comparing our JWST spectra to VLT/CRIRES+, Spitzer/IRS, and ALMA. Primary and secondary disks are dramatically different in their mid-infrared emission, with primary disks showing H2O-rich spectra, and secondary disks being mostly line poor to the sensitivity of our spectra. When comparing MIRI-MRS to Spitzer/IRS, we observe large variability in the line emission of VW Cha A, as well as in the continuum of RW Aur A. The disks around VW Cha BC and RW Aur B show evidence of ionizing radiation, and a further comparison with ALMA at high angular resolution dust continuum suggest that the spectrum of RW Aur B is well explained by its ~4 au cavity. All the systems show [Ne II] jet emission, and three of them also show spatially resolved emission structures in H2, likely originated by outflows and dynamical interactions. Many of the observed features in the primary disks, such as enhanced water emission, could be linked to the increased accretion and radial drift produced by dynamical disk truncation. However, additional mechanisms are needed to explain the large differences between primary and secondary disks, potentially inner disk substructures. This work is an example of the need for combining multiple facilities to fully understand the observations from JWST.

preprint2025arXiv

Revisiting the atmosphere of HAT-P-70b with CARMENES high-resolution transmission spectroscopy

Owing to hot and inflated envelopes that facilitate atmospheric studies, ultra-hot Jupiters (UHJs) have attracted much attention. Significant progress has been achieved, from enlarging the sample size to broadening the studies to encompass diverse stellar types and ages. Here, we present a transmission spectroscopy study of HAT-P-70b, an UHJ orbiting a young A-type star, through high-resolution observations with CARMENES at the 3.5m Calar Alto telescope. By using the line-by-line technique, we confirm the previous detections of Ha, Na I, and Ca II, report a new tentative detection of K I, and impose an upper limit on the He triplet absorption. Through cross-correlation analysis, we identify the Ca II and Fe I absorptions, both blue-shifted by approximately 5 km/s, indicating a day-to-night side atmospheric wind. Additionally, we find a new tentative detection of K I. We do not see any significant atmospheric molecular signal in the near-infrared data. Putting HAT-P-70b in the context of UHJs from the literature, it turns out that (1) Ha absorption is more common on gas giants orbiting stars younger than 1 Gyr, with a relative detection probability of $P_{\rm Age<1\,Gyr}({\rm Ha})/P_{\rm Age\geq1\,Gyr}({\rm Ha})\sim 3$; (2) any UHJ is likely to exhibit Fe I absorption if it has Ca II.

preprint2023arXiv

Forbidden emission lines in protostellar outflows and jets with MUSE

Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral field spectrograph of the Multi Unit Spectroscopic Explorer (MUSE), at the Very Large Telescope (VLT). With a resolution of 0.025 X 0.025 arcsec$^{2}$, we aim to derive the position angle of the outflow/jet (PA$_{outflow/jet}$) that is connected with the inner disk. The forbidden emission lines analyzed here have their origin at the inner parts of the protoplanetary disk. From the maximum intensity emission along the outflow/jet in DL Tau, CI Tau, DS Tau, IP Tau, and IM Lup, we were able to reliably measure the PA$_{outflow/jet}$ for most of the identified lines. We found that our estimates agree with PA$_{dust}$ for most of the disks. These estimates depend on the signal-to-noise level and the collimation of the outflow (jet). The outflows/jets in CIDA 9, GO Tau, and GW Lup are too compact for a PA$_{outflow/jet}$ to be estimated. Based on our kinematics analysis, we confirm that DL Tau and CI Tau host a strong outflow/jet with line-of-sight velocities much greater than 100 km s$^{-1}$, whereas DS Tau, IP Tau, and IM Lup velocities are lower and their structures encompass low-velocity components to be more associated with winds. Our estimates for the mass-loss rate, $\dot{M}_{loss}$, range between (1.1-6.5)X10$^{-7}$-10$^{-8}$ $M_{\odot}$ yr$^{-1}$ for the disk-outflow/jet systems analyzed here. The outflow/jet systems analyzed here are aligned within around 1 degree between the inner and outer disk. Further observations are needed to confirm a potential misalignment in IM Lup.

preprint2022arXiv

A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508

We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass $m \sin i = 4.00^{+0.53}_{-0.55}\ M_{\oplus}$. We find no evidence of significant signals at the detected period in spectroscopic stellar activity indicators or MEarth photometry. The planet, Ross 508 b, has a semimajor-axis of $0.05366^{+0.00056}_{-0.00049}$ au. This gives an orbit-averaged insolation of $\approx$1.4 times the Earth's value, placing Ross 508 b near the inner edge of its star's habitable zone. We have explored the possibility that the planet has a high eccentricity and its host is accompanied by an additional unconfirmed companion on a wide orbit. Our discovery demonstrates that the near-infrared RV search can play a crucial role to find a low-mass planet around cool M dwarfs like Ross 508.

preprint2022arXiv

A transiting, temperate mini-Neptune orbiting the M dwarf TOI-1759 unveiled by TESS

We report the discovery and characterization of TOI-1759~b, a temperate (400 K) sub-Neptune-sized exoplanet orbiting the M~dwarf TOI-1759 (TIC 408636441). TOI-1759 b was observed by TESS to transit on sectors 16, 17 and 24, with only one transit observed per sector, creating an ambiguity on the orbital period of the planet candidate. Ground-based photometric observations, combined with radial-velocity measurements obtained with the CARMENES spectrograph, confirm an actual period of $18.85019 \pm 0.00014$ d. A joint analysis of all available photometry and radial velocities reveal a radius of $3.17 \pm 0.10\,R_\oplus$ and a mass of $10.8 \pm 1.5\,M_\oplus$. Combining this with the stellar properties derived for TOI-1759 ($R_\star = 0.597 \pm 0.015\,R_\odot$; $M_\star = 0.606 \pm 0.020\,M_\odot$; $T_{\textrm{eff}} = 4065 \pm 51$ K), we compute a transmission spectroscopic metric (TSM) value of over 80 for the planet, making it a good target for transmission spectroscopy studies. TOI-1759 b is among the top five temperate, small exoplanets ($T_\textrm{eq} < 500$ K, $R_p < 4 \,R_\oplus$) with the highest TSM discovered to date. Two additional signals with periods of 80 d and $>$ 200 d seem to be present in our radial velocities. While our data suggest both could arise from stellar activity, the later signal's source and periodicity are hard to pinpoint given the $\sim 200$ d baseline of our radial-velocity campaign with CARMENES. Longer baseline radial-velocity campaigns should be performed in order to unveil the true nature of this long period signal.

preprint2022arXiv

Dynamical architecture of the HD 107148 system

We present an independent Doppler validation and dynamical orbital analysis of the two-planet system HD 107148, which was recently announced in Rosenthal et al. (2021). Our detailed analyses are based on literature HIRES data and newly obtained HARPS and CARMENES radial velocity (RV) measurements as part of our survey in search for additional planets around single planet systems. We perform a periodogram analysis of the available HIRES and HARPS precise RVs and stellar activity indicators. We do not find any apparent correlation between the RV measurements and the stellar activity indicators, thus linking the two strong periodicities to a moderately compact multiple-planet system. We carry out orbital fitting analysis by testing various one- and two-planet orbital configurations and studying the posterior probability distribution of the fitted parameters. Our results solidify the existence of a Saturn-mass planet (HD 107148 b, discovered first) with a period $P_b\sim77.2$ d, and a second, eccentric ($e_c \sim$ 0.4), Neptune-mass exoplanet (HD 107148 c), with an orbital period of $P_c\sim18.3$ d. Finally, we investigate the two-planet system's long-term stability and overall orbital dynamics with the posterior distribution of our preferred orbital configuration. Our N-body stability simulations show that the system is long-term stable and exhibits large secular osculations in eccentricity but in no particular mean-motion resonance configuration. The HD 107148 system, consisting of a Solar-type main sequence star with two giant planets in a rare configuration, features a common proper motion white dwarf companion and is, therefore, a valuable target for understanding the formation and evolution of planetary systems.

preprint2022arXiv

Exploring refractory organics in extraterrestrial particles

The origin of organic compounds detected in meteorites and comets, some of which could serve as precursors of life on Earth, still remains an open question. The aim of the present study is to make one more step in revealing the nature and composition of organic materials of extraterrestrial particles by comparing infrared spectra of laboratory-made refractory organic residues to the spectra of cometary particles returned by the Stardust mission, interplanetary dust particles, and meteorites. Our results reinforce the idea of a pathway for the formation of refractory organics through energetic and thermal processing of molecular ices in the solar nebula. There is also the possibility that some of the organic material formed already in the parental molecular cloud before it entered the solar nebula. The majority of the IR "organic" bands of the studied extraterrestrial particles can be reproduced in the spectra of the laboratory organic residues. We confirm the detection of water, nitriles, hydrocarbons, and carbonates in extraterrestrial particles and link it to the formation location of the particles in the outer regions of the solar nebula. To clarify the genesis of the species, high-sensitivity observations in combination with laboratory measurements like those presented in this paper are needed. Thus, this study presents one more piece of the puzzle of the origin of water and organic compounds on Earth and motivation for future collaborative laboratory and observational projects.

preprint2022arXiv

Interpreting the atmospheric composition of exoplanets: sensitivity to planet formation assumptions

Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations.

preprint2022arXiv

LBT search for companions and sub-structures in the (pre)transitional disk of AB Aurigae

Multi-wavelengths high-resolution imaging of protoplanetary disks has revealed the presence of multiple, varied substructures in their dust and gas components which might be signposts of young, forming planetary systems. AB Aurigae bears an emblematic (pre)transitional disk showing spiral structures observed in the inner cavity of the disk in both the sub-millimeter (ALMA; 1.3mm, $^{12}$CO) and near-infrared (SPHERE; 1.5-2.5$μ$m) wavelengths which have been claimed to arise from dynamical interactions with a massive companion. In this work, we present new deep $K_s$ (2.16$μ$m) and L' (3.7$μ$m) band images of AB Aurigae obtained with LMIRCam on the Large Binocular Telescope, aimed for the detection of both planetary companions and extended disk structures. No point source is recovered, in particular at the outer regions of the disk, where a putative candidate ($ρ= 0.681", PA = 7.6^{\circ}$) had been previously claimed. The nature of a second innermost planet candidate ($ρ= 0.16'', PA = 203.9^{\circ}$) can not be investigated by the new data. We are able to derive 5$σ$ detection limits in both magnitude and mass for the system, going from 14 \Mjup at 0.3'' (49 au) down to 3-4 \Mjup at 0.6'' (98 au) and beyond, based on the ATMO 2020 evolutionary models. We detect the inner spiral structures (< 0.5'') resolved in both CO and polarimetric H-band observations. We also recover the ring structure of the system at larger separation (0.5-0.7") showing a clear south-east/north-west asymmetry. This structure, observed for the first time at L'-band, remains interior to the dust cavity seen at ALMA, suggesting an efficient dust trapping mechanism at play in the disk.

preprint2022arXiv

Locating dust and molecules in the inner circumstellar environment of R~Sculptoris with MATISSE

AGB stars are one of the main sources of dust production in the Galaxy. However, it is not clear what this process looks like and where the dust is condensing in the circumstellar environment. By characterizing the location of the dust and the molecules in the close environment of an AGB star, we aim to achieve a better understanding the history of the dust formation process. We observed the carbon star R Scl with the VLTI-MATISSE instrument in L- and N-bands. The high angular resolution of the VLTI observations, combined with a large uv-plane coverage allowed us to use image reconstruction methods. To constrain the dust and molecules' location, we used two different methods: MIRA image reconstruction and the 1D code RHAPSODY. We found evidence of C2H2 and HCN molecules between 1 and 3.4 Rstar which is much closer to the star than the location of the dust (between 3.8 and 17.0 Rstar). We also estimated a mass-loss rate of 1.2+-0.4x10-6 Msun per yr. In the meantime, we confirmed the previously published characteristics of a thin dust shell, composed of amorphous carbon (amC) and silicon carbide (SiC). However, no clear SiC feature has been detected in the MATISSE visibilities. This might be caused by molecular absorption that can affect the shape of the SiC band at 11.3 micron. The appearance of the molecular shells is in good agreement with predictions from dynamical atmosphere models. For the first time, we co-located dust and molecules in the environment of an AGB star. We confirm that the molecules are located closer to the star than the dust. The MIRA images unveil the presence of a clumpy environment in the fuzzy emission region beyond 4.0 Rstar. Furthermore, with the available dynamic range and angular resolution, we did not detect the presence of a binary companion. Additional observations combining MATISSE and SAM-VISIR instrument should enable this detection in future studies.

preprint2022arXiv

RV-detected planets around M dwarfs: Challenges for core accretion models

Planet formation is sensitive to the conditions in protoplanetary disks, for which scaling laws as a function of stellar mass are known. We aim to test whether the observed population of planets around low-mass stars can be explained by these trends, or if separate formation channels are needed. We address this question by confronting a state-of-the-art planet population synthesis model with a sample of planets around M dwarfs observed by the HARPS and CARMENES radial velocity (RV) surveys. To account for detection biases, we performed injection and retrieval experiments on the actual RV data to produce synthetic observations of planets that we simulated following the core accretion paradigm. These simulations robustly yield the previously reported high occurrence of rocky planets around M dwarfs and generally agree with their planetary mass function. In contrast, our simulations cannot reproduce a population of giant planets around stars less massive than 0.5 solar masses. This potentially indicates an alternative formation channel for giant planets around the least massive stars that cannot be explained with current core accretion theories. We further find a stellar mass dependency in the detection rate of short-period planets. A lack of close-in planets around the earlier-type stars ($M_\star \gtrsim 0.4\, M_\odot$) in our sample remains unexplained by our model and indicates dissimilar planet migration barriers in disks of different spectral subtypes. Both discrepancies can be attributed to gaps in our understanding of planet migration in nascent M dwarf systems. They underline the different conditions around young stars of different spectral subtypes, and the importance of taking these differences into account when studying planet formation.

preprint2022arXiv

The dusty heart of Circinus: I. Imaging the circumnuclear dust in N-band

Active galactic nuclei play a key role in the evolution of galaxies, but their inner workings and physical connection to the host are poorly understood due to a lack of angular resolution. Infrared interferometry makes it possible to resolve the circumnuclear dust in the nearby Seyfert 2 galaxy, Circinus. Previous observations have revealed complex structures and polar dust emission but interpretation was limited to simple models. MATISSE makes it possible to image these structures for the first time. We observed the Circinus Galaxy with VLTI/MATISSE, producing 150 correlated flux spectra and 100 closure phase spectra. We reconstructed images in the N-band at ~10 mas resolution. We fit blackbody functions with dust extinction to several aperture-extracted fluxes from the images to produce a temperature distribution of central dusty structures. We find significant substructure in the circumnuclear dust: central unresolved flux of ~0.5 Jy, a thin disk 1.9 pc in diameter oriented along ~45 deg,and a ~4x1.5 pc polar emission extending orthogonal to the disk. The polar emission exhibits patchiness, which we attribute to clumpy dust. Flux enhancements to the east and west of the disk are seen for the first time. We distinguish the temperature profiles of the disk and of the polar emission: the disk shows a steep temperature gradient indicative of denser material; the polar profile is flatter, indicating clumpiness and/or lower dust density. The unresolved flux is fitted with a high temperature, ~370 K. The polar dust remains warm (~200 K) out to 1.5 pc from the disk. The recovered morphology and temperature distribution resembles modeling of accretion disks with radiation-driven winds at large scales, but we placed new constraints on the subparsec dust. The subparsec features imaged here place new constraints on the physical modeling of circumnuclear dust in active galaxies.

preprint2022arXiv

The PEPSI Exoplanet Transit Survey (PETS) I: Investigating the presence of a silicate atmosphere on the super-Earth 55 Cnc e

The study of exoplanets and especially their atmospheres can reveal key insights on their evolution by identifying specific atmospheric species. For such atmospheric investigations, high-resolution transmission spectroscopy has shown great success, especially for Jupiter-type planets. Towards the atmospheric characterization of smaller planets, the super-Earth exoplanet 55 Cnc e is one of the most promising terrestrial exoplanets studied to date. Here, we present a high-resolution spectroscopic transit observation of this planet, acquired with the PEPSI instrument at the Large Binocular Telescope. Assuming the presence of Earth-like crust species on the surface of 55 Cnc e, from which a possible silicate-vapor atmosphere could have originated, we search in its transmission spectrum for absorption of various atomic and ionized species such as Fe , Fe+, Ca , Ca+, Mg and K , among others. Not finding absorption for any of the investigated species, we are able to set absorption limits with a median value of 1.9 x RP. In conclusion, we do not find evidence of a widely extended silicate envelope on this super-Earth reaching several planetary radii.

preprint2022arXiv

The TESS Grand Unified Hot Jupiter Survey. I. Ten TESS Planets

We report the discovery of ten short-period giant planets (TOI-2193A b, TOI-2207 b, TOI-2236 b, TOI-2421 b, TOI-2567 b, TOI-2570 b, TOI-3331 b, TOI-3540A b, TOI-3693 b, TOI-4137 b). All of the planets were identified as planet candidates based on periodic flux dips observed by NASA's Transiting Exoplanet Survey Satellite (TESS). The signals were confirmed to be from transiting planets using ground-based time-series photometry, high angular resolution imaging, and high-resolution spectroscopy coordinated with the TESS Follow-up Observing Program. The ten newly discovered planets orbit relatively bright F and G stars ($G < 12.5$,~$T_\mathrm{eff}$ between 4800 and 6200 K). The planets' orbital periods range from 2 to 10~days, and their masses range from 0.2 to 2.2 Jupiter masses. TOI-2421 b is notable for being a Saturn-mass planet and TOI-2567 b for being a ``sub-Saturn'', with masses of $0.322\pm 0.073$ and $0.195\pm 0.030$ Jupiter masses, respectively. In most cases, we have little information about the orbital eccentricities. Two exceptions are TOI-2207 b, which has an 8-day period and a detectably eccentric orbit ($e = 0.17\pm0.05$), and TOI-3693 b, a 9-day planet for which we can set an upper limit of $e < 0.052$. The ten planets described here are the first new planets resulting from an effort to use TESS data to unify and expand on the work of previous ground-based transit surveys in order to create a large and statistically useful sample of hot Jupiters.

preprint2022arXiv

The vertical structure of debris disks and the impact of gas

The vertical structure of debris disks provides clues about their dynamical evolution and the collision rate of the unseen planetesimals. Thanks to the ever-increasing angular resolution of contemporary instruments and facilities, we are beginning to constrain the scale height of a handful of debris disks, either at near-infrared or millimeter wavelengths. Nonetheless, this is often done for individual targets only. We present here the geometric modeling of eight disks close to edge-on, all observed with the same instrument (SPHERE) and using the same mode (dual-beam polarimetric imaging). Motivated by the presence of CO gas in two out of the eight disks, we then investigate the impact that gas can have on the scale height by performing N-body simulations including gas drag and collisions. We show that gas can quickly alter the dynamics of particles (both in the radial and vertical directions), otherwise governed by gravity and radiation pressure. We find that, in the presence of gas, particles smaller than a few tens of microns can efficiently settle toward the midplane at the same time as they migrate outward beyond the birth ring. For second generation gas ($M_\mathrm{gas} \leq 0.1$ $M_\oplus$), the vertical settling should be best observed in scattered light images compared to observations at millimeter wavelengths. But if the gas has a primordial origin ($M_\mathrm{gas} \geq 1$ $M_\oplus$), the disk will appear very flat both at near-infrared and sub-mm wavelengths. Finally, far beyond the birth ring, our results suggest that the surface brightness profile can be as shallow as $\sim -2.25$.

preprint2022arXiv

Towards RNA life on Early Earth: From atmospheric HCN to biomolecule production in warm little ponds

The origin of life on Earth involves the early appearance of an information-containing molecule such as RNA. The basic building blocks of RNA could have been delivered by carbon-rich meteorites, or produced in situ by processes beginning with the synthesis of hydrogen cyanide (HCN) in the early Earth's atmosphere. Here, we construct a robust physical and non-equilibrium chemical model of the early Earth atmosphere. The atmosphere is supplied with hydrogen from impact degassing of meteorites, sourced with water evaporated from the oceans, carbon dioxide from volcanoes, and methane from undersea hydrothermal vents, and in which lightning and external UV-driven chemistry produce HCN. This allows us to calculate the rain-out of HCN into warm little ponds (WLPs). We then use a comprehensive sources and sinks numerical model to compute the resulting abundances of nucleobases, ribose, and nucleotide precursors such as 2-aminooxazole resulting from aqueous and UV-driven chemistry within them. We find that at 4.4 bya (billion years ago) the limits of adenine concentrations in ponds for habitable surfaces is 0.05$μ$M in the absence of seepage. These concentrations can be maintained for over 100 Myr. Meteorite delivery of adenine to WLPs can provide boosts in concentration by 2-3 orders of magnitude, but these boosts deplete within months by UV photodissociation, seepage, and hydrolysis. The early evolution of the atmosphere is dominated by the decrease of hydrogen due to falling impact rates and atmospheric escape, and the rise of oxygenated species such as OH from H2O photolysis. Our work points to an early origin of RNA on Earth within ~200 Myr of the Moon-forming impact.

preprint2021arXiv

$\textit{TESS}$ Giants Transiting Giants I: A Non-inflated Hot Jupiter Orbiting a Massive Subgiant

While the population of confirmed exoplanets continues to grow, the sample of confirmed transiting planets around evolved stars is still limited. We present the discovery and confirmation of a hot Jupiter orbiting TOI-2184 (TIC 176956893), a massive evolved subgiant ($M_\star= 1.53 \pm 0.12 M_\odot$, $R_\star= 2.90 \pm 0.14 R_\odot$) in the $\textit{TESS}$ Southern Continuous Viewing Zone. The planet was flagged as a false positive by the $\textit{TESS}$ Quick-Look Pipeline due to periodic systematics introducing a spurious depth difference between even and odd transits. Using a new pipeline to remove background scattered light in $\textit{TESS}$ Full Frame Image (FFI) data, we combine space-based $\textit{TESS}$ photometry, ground-based photometry, and ground-based radial velocity measurements to report a planet radius of $R_p= 1.017 \pm 0.051 R_J$ and mass of $M_p= 0.65 \pm 0.16 M_J$. For a planet so close to its star, the mass and radius of TOI-2184b are unusually well matched to those of Jupiter. We find that the radius of TOI-2184b is smaller than theoretically predicted based on its mass and incident flux, providing a valuable new constraint on the timescale of post-main-sequence planet inflation. The discovery of TOI-2184b demonstrates the feasibility of detecting planets around faint ($\textit{TESS}$ magnitude $>12$) post-main sequence stars and suggests that many more similar systems are waiting to be detected in the $\textit{TESS}$ FFIs, whose confirmation may elucidate the final stages of planetary system evolution.

preprint2021arXiv

A wide-orbit giant planet in the high-mass b Centauri binary system

Planet formation occurs around a wide range of stellar masses and stellar system architectures. An improved understanding of the formation process can be achieved by studying it across the full parameter space, particularly toward the extremes. Earlier studies of planets in close-in orbits around high-mass stars have revealed an increase in giant planet frequency with increasing stellar mass until a turnover point at 1.9 solar masses, above which the frequency rapidly decreases. This could potentially imply that planet formation is impeded around more massive stars, and that giant planets around stars exceeding 3 solar masses may be rare or non-existent. However, the methods used to detect planets in small orbits are insensitive to planets in wide orbits. Here we demonstrate the existence of a planet at 560 times the Sun-Earth distance from the 6-10 solar mass binary b Centauri through direct imaging. The planet-to-star mass ratio of 0.10-0.17% is similar to the Jupiter-Sun ratio, but the separation of the detected planet is ~100 times wider than that of Jupiter. Our results show that planets can reside in much more massive stellar systems than what would be expected from extrapolation of previous results. The planet is unlikely to have formed in-situ through the conventional core accretion mechanism, but might have formed elsewhere and arrived to its present location through dynamical interactions, or might have formed via gravitational instability.

preprint2021arXiv

BEAST begins: Sample characteristics and survey performance of the B-star Exoplanet Abundance Study

While the occurrence rate of wide giant planets appears to increase with stellar mass at least up through the A-type regime, B-type stars have not been systematically studied in large-scale surveys so far. It therefore remains unclear up to what stellar mass this occurrence trend continues. The B-star Exoplanet Abundance Study (BEAST) is a direct imaging survey with the extreme adaptive optics instrument SPHERE, targeting 85 B-type stars in the young Scorpius-Centaurus (Sco-Cen) region with the aim to detect giant planets at wide separations and constrain their occurrence rate and physical properties. The statistical outcome of the survey will help determine if and where an upper stellar mass limit for planet formation occurs. In this work, we describe the selection and characterization of the BEAST target sample. Particular emphasis is placed on the age of each system, which is a central parameter in interpreting direct imaging observations. We implement a novel scheme for age dating based on kinematic sub-structures within Sco-Cen, which complements and expands upon previous age determinations in the literature. We also present initial results from the first epoch observations, including the detections of ten stellar companions, of which six were previously unknown. All planetary candidates in the survey will need follow up in second epoch observations, which are part of the allocated observational programme and will be executed in the near future.

preprint2021arXiv

Dynamical masses for two M1 + mid-M dwarf binaries monitored during the SPHERE-SHINE survey

We present orbital fits and dynamical masses for HIP 113201AB and HIP 36985AB, two M1 + mid-M dwarf binary systems monitored as part of the SPHERE SHINE survey. To robustly determine ages via gyrochronology, we undertook a photometric monitoring campaign for HIP 113201 and for GJ 282AB, the two wide K star companions to HIP 36985, using the 40 cm Remote Observatory Atacama Desert (ROAD) telescope. We adopt ages of 1.2$\pm$0.1 Gyr for HIP 113201AB and 750$\pm$100 Myr for HIP 36985AB. To derive dynamical masses for all components of these systems, we used parallel-tempering Markov Chain Monte Carlo sampling to fit a combination of radial velocity, direct imaging, and Gaia and Hipparcos astrometry. Fitting the direct imaging and radial velocity data for HIP 113201 yields a primary mass of 0.54$\pm$0.03 M$_{\odot}$, fully consistent with its M1 spectral type, and a secondary mass of 0.145$\pm$ M$_{\odot}$. The secondary masses derived with and without including Hipparcos/Gaia data are more massive than the 0.1 M$_{\odot}$ estimated mass from the photometry of the companion. An undetected brown dwarf companion to HIP 113201B could be a natural explanation for this apparent discrepancy. At an age $>$1 Gyr, a 30 M$_{Jup}$ companion to HIP 113201B would make a negligible ($<$1$\%$) contribution to the system luminosity, but could have strong dynamical impacts. Fitting the direct imaging, radial velocity, and Hipparcos/Gaia proper motion anomaly for HIP 36985AB, we find a primary mass of 0.54$\pm$0.01 M$_{\odot}$ and a secondary mass of 0.185$\pm$0.001 M$_{\odot}$ which agree well with photometric estimates of component masses, the masses estimated from $M_{K}$-- mass relationships for M dwarf stars, and previous dynamical masses in the literature.

preprint2021arXiv

Elemental abundances of nearby M dwarfs based on high-resolution near-infrared spectra obtained by the Subaru/IRD survey: Proof of concept

Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolution (~70,000) near-infrared (970-1750 nm) spectra to measure the abundances of Na, Mg, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Sr by the line-by-line analysis based on model atmospheres, with typical errors ranging from 0.2 dex for [Fe/H] to 0.3-0.4 dex for other [X/H]. We measure radial velocities from the spectra and combine them with Gaia astrometry to calculate the Galactocentric space velocities UVW. The resulting [Fe/H] values agree with previous estimates based on medium-resolution K-band spectroscopy, showing a wide distribution of metallicity (-0.6 < [Fe/H] < +0.4). The abundance ratios of individual elements [X/Fe] are generally aligned with the solar values in all targets. While the [X/Fe] distributions are comparable to those of nearby FGK stars, most of which belong to the thin disk population, the most metal-poor object, GJ 699, could be a thick disk star. The UVW velocities also support this. The results raise the prospect that near-infrared spectra of M dwarfs obtained in the planet search projects can be used to grasp the trend of elemental abundances and Galactic stellar population of nearby M dwarfs.

preprint2021arXiv

Mass determination of protoplanetary disks from dust evolution

The mass of protoplanetary disks is arguably one of their most important quantities shaping their evolution toward planetary systems, but it remains a challenge to determine this quantity. Using the high spatial resolution now available on telescopes such as the Atacama Large Millimeter/submillimeter Array (ALMA), recent studies derived a relation between the disk surface density and the location of the "dust lines". This is a new concept in the field, linking the disk size at different continuum wavelengths with the radial distribution of grain populations of different sizes. We aim to use a dust evolution model to test the dependence of the dust line location on disk gas mass. In particular, we are interested in the reliability of the method for disks showing radial substructures, as recent high-resolution observations revealed. Our models show that the determination of the dust line location is a promising approach to the mass estimate of protoplanetay disks, but the exact relation between the dust line location and disk mass depends on the structure of the particular disk. We calibrated the relation for disks without evidence of radial structures, while for more complex structures we ran a simple dust evolution model. However, this method fails when there is evidence of strong dust traps. It is possible to reveal when dust evolution is dominated by traps, providing the necessary information for when the method should be applied with caution.

preprint2021arXiv

Precise transit and radial-velocity characterization of a resonant pair: a warm Jupiter TOI-216c and eccentric warm Neptune TOI-216b

TOI-216 hosts a pair of warm, large exoplanets discovered by the TESS Mission. These planets were found to be in or near the 2:1 resonance, and both of them exhibit transit timing variations (TTVs). Precise characterization of the planets' masses and radii, orbital properties, and resonant behavior can test theories for the origins of planets orbiting close to their stars. Previous characterization of the system using the first six sectors of TESS data suffered from a degeneracy between planet mass and orbital eccentricity. Radial velocity measurements using HARPS, FEROS, and PFS break that degeneracy, and an expanded TTV baseline from TESS and an ongoing ground-based transit observing campaign increase the precision of the mass and eccentricity measurements. We determine that TOI-216c is a warm Jupiter, TOI-216b is an eccentric warm Neptune, and that they librate in the 2:1 resonance with a moderate libration amplitude of 60 +/- 2 degrees; small but significant free eccentricity of 0.0222 +0.0005/-0.0003 for TOI-216b; and small but significant mutual inclination of 1.2-3.9 degrees (95% confidence interval). The libration amplitude, free eccentricity, and mutual inclination imply a disturbance of TOI-216b before or after resonance capture, perhaps by an undetected third planet.

preprint2021arXiv

Thermal imaging of dust hiding the black hole in the Active Galaxy NGC 1068

In the widely accepted 'Unified Model' solution of the classification puzzle of Active Galactic Nuclei, the orientation of a dusty accretion torus around the central black hole dominates their appearance. In 'type-1' systems, the bright nucleus is visible at the centre of a face-on torus. In 'type-2' systems the thick, nearly edge-on torus hides the central engine. Later studies suggested evolutionary effects and added dusty clumps and polar winds but left the basic picture intact. However, recent high-resolution images of the archetypal type-2 galaxy NGC 1068 suggested a more radical revision. They displayed a ring-like emission feature which the authors advocated to be hot dust surrounding the black hole at the radius where the radiation from the central engine evaporates the dust. That ring is too thin and too far tilted from edge-on to hide the central engine, and ad hoc foreground extinction is needed to explain the type-2 classification. These images quickly generated reinterpretations of the type 1-2 dichotomy. Here we present new multi-band mid-infrared images of NGC1068 that detail the dust temperature distribution and reaffirm the original model. Combined with radio data, our maps locate the central engine below the previously reported ring and obscured by a thick, nearly edge-on disk, as predicted by the Unified Model. We also identify emission from polar flows and absorbing dust that is mineralogically distinct from that towards the Milky Way centre.

preprint2021arXiv

TOI-954 b and K2-329 b: Short-Period Saturn-Mass Planets that Test whether Irradiation Leads to Inflation

We report the discovery of two short-period Saturn-mass planets, one transiting the G subgiant TOI-954 (TIC 44792534, $ V = 10.343 $, $ T = 9.78 $) observed in TESS sectors 4 and 5, and one transiting the G dwarf K2-329 (EPIC 246193072, $ V = 12.70 $, $ K = 10.67 $) observed in K2 campaigns 12 and 19. We confirm and characterize these two planets with a variety of ground-based archival and follow-up observations, including photometry, reconnaissance spectroscopy, precise radial velocity, and high-resolution imaging. Combining all available data, we find that TOI-954 b has a radius of $0.852_{-0.062}^{+0.053} \, R_{\mathrm{J}}$ and a mass of $0.174_{-0.017}^{+0.018} \, M_{\mathrm{J}}$ and is in a 3.68 day orbit, while K2-329 b has a radius of $0.774_{-0.024}^{+0.026} \, R_{\mathrm{J}}$ and a mass of $0.260_{-0.022}^{+0.020} \, M_{\mathrm{J}}$ and is in a 12.46 day orbit. As TOI-954 b is 30 times more irradiated than K2-329 b but more or less the same size, these two planets provide an opportunity to test whether irradiation leads to inflation of Saturn-mass planets and contribute to future comparative studies that explore Saturn-mass planets at contrasting points in their lifetimes.

preprint2021arXiv

Unveiling wide-orbit companions to K-type stars in Sco-Cen with Gaia EDR3

Abbreviated. We aim to identify new low-mass companions to young stars using the astrometric measurements provided by the Gaia space mission and complementary VLT/SPHERE data. We identify companion candidates from a sample of K-type, pre-main sequence stars in the Scorpius Centaurus association using the early version of the third data release of the Gaia space mission. Based on the provided positions, proper motions, and magnitudes, we identify all objects within a predefined radius whose differential proper motions are consistent with a gravitationally bound system. We derive companion masses through comparison with evolutionary tracks. For seven identified companion candidates we use additional data collected with VLT/SPHERE and VLT/NACO to assess the accuracy of the properties of the companions based on Gaia photometry alone. We identify 110 comoving companions that have a companionship likelihood of more than $95\,\%$. We identify ten especially intriguing companions that have masses in the brown dwarf regime down to $20\,M_\mathrm{Jup}$. Our high-contrast imaging data confirm both astrometry and photometric masses derived from Gaia alone. We discover a new brown dwarf companion, TYC 8252-533-1 B, with a projected separation of approximately $570\,\mathrm{au}$ from its Sun-like primary. SED modeling provides a companion mass of $52^{+17}_{-11}\,M_\mathrm{Jup}$. We show that the Gaia database can identify low-mass companions at wide separations from their host stars. For K-type Sco-Cen members Gaia can detect sub-stellar objects at projected separations larger than $300\,\mathrm{au}$ and is sensitivity limited beyond $1,000\,\mathrm{au}$ with a lower mass limit down to $20\,M_\mathrm{Jup}$. A similar analysis of other star-forming regions could significantly enlarge the sample size of such objects and test formation and evolution theories of planetary systems.

preprint2020arXiv

A low-mass stellar companion to the young variable star RZ Psc

RZ Psc is a young Sun-like star with a bright and warm infrared excess that is occasionally dimmed significantly by circumstellar dust structures. Optical depth arguments suggest that the dimming events do not probe a typical sight line through the circumstellar dust, and are instead caused by structures that appear above an optically thick mid-plane. This system may therefore be similar to systems where an outer disk is shadowed by material closer to the star. Here we report the discovery that RZ Psc hosts a 0.12$M_\odot$ companion at a projected separation of 23au. We conclude that the disk must orbit the primary star. While we do not detect orbital motion, comparison of the angle of linear polarization of the primary with the companion's on-sky position angle provides circumstantial evidence that the companion and disc may not share the same orbital plane. Whether the companion severely disrupts the disc, truncates it, or has little effect at all, will require further observations of both the companion and disc.

preprint2020arXiv

Birth of convective low-mass to high-mass second Larson cores

Stars form as an end product of the gravitational collapse of cold, dense gas in magnetized molecular clouds. This multi-scale scenario occurs via the formation of two quasi-hydrostatic cores and involves complex physical processes, which require a robust, self-consistent numerical treatment. The aim of this study is to understand the formation and evolution of the second Larson core and the dependence of its properties on the initial cloud core mass. We used the PLUTO code to perform high resolution, 1D and 2D RHD collapse simulations. We include self-gravity and use a grey FLD approximation for the radiative transfer. Additionally, we use for the gas EOS density- and temperature-dependent thermodynamic quantities to account for the effects such as dissociation, ionisation, and molecular vibrations and rotations. Properties of the second core are investigated using 1D studies spanning a wide range of initial cloud core masses from 0.5 to 100 $M_{\odot}$. Furthermore, we expand to 2D collapse simulations for a few cases of 1, 5, 10, and 20 $M_{\odot}$. We follow the evolution of the second core for $\geq$ 100 years after its formation, for each of these non-rotating cases. Our results indicate a dependence of several second core properties on the initial cloud core mass. For the first time, due to an unprecedented resolution, our 2D non-rotating collapse studies indicate that convection is generated in the outer layers of the second core, which is formed due to the gravitational collapse of a 1 $M_{\odot}$ cloud core. Additionally, we find large-scale oscillations of the second accretion shock front triggered by the standing accretion shock instability, which has not been seen before in early evolutionary stages of stars. We predict that the physics within the second core would not be significantly influenced by the effects of magnetic fields or an initial cloud rotation.

preprint2020arXiv

Chemical modelling of FU Ori protoplanetary disks

Luminosity outbursts of the FUOri type stars, which have a magnitude of $\sim100 L_{\odot}$ and last for decades, may affect chemical composition of the surrounding protoplanetary disk. Using astrochemical modeling we analyze the changes induced by the outburst and search for species sensitive to the luminosity rise. Some changes in the disk molecular composition appear not only during the outburst itself but can also retain for decades after the end of the outburst. We analyze main chemical processes responsible for these effects and assess timescales at which chemically inert species return to the pre-outburst abundances.

preprint2020arXiv

Dust/ice mixing in cold regions and solid-state water in the diffuse interstellar medium

Whether ice in cold cosmic environments is physically separated from the silicate dust or mixed with individual silicate moieties is not known. However, different grain models give very different compositions and temperatures of grains. The aim of the present study is a comparison of the mid-IR spectra of laboratory silicate-grains/water-ice mixtures with astronomical observations to evaluate the presence of dust/ice mixtures in interstellar and circumstellar environments. The laboratory data can explain the observations assuming reasonable mass-averaged temperatures for the protostellar envelopes and protoplanetary disks demonstrating that a substantial fraction of water ice may be mixed with silicate grains. Based on the combination of laboratory data and infrared observations, we provide evidence of the presence of solid-state water in the diffuse interstellar medium. Our results have implications for future laboratory studies trying to investigate cosmic dust grain analogues and for future observations trying to identify the structure, composition, and temperature of grains in different astrophysical environments.

preprint2020arXiv

EDEN: Sensitivity Analysis and Transiting Planet Detection Limits for Nearby Late Red Dwarfs

Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches.

preprint2020arXiv

Equatorial anti-rotating day side wind flow in WASP-43b elicited by deep wind jets?

We present WASP-43b climate simulations with deep wind jets (down to 700~bar) that are linked to retrograde (westward) flow at the equatorial day side for $p<0.1$~bar. Retrograde flow inhibits efficient eastward heat transport and naturally explains the small hotspot shift and large day-night-side gradient of WASP-43b ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days) observed with Spitzer. We find that deep wind jets are mainly associated with very fast rotations ($P_{\text{rot}}=P_{\text{orb}}\leq 1.5$~days) which correspond to the Rhines length smaller than $2$ planetary radii. We also diagnose wave activity that likely gives rise to deviations from superrotation. Further, we show that we can achieve full steady state in our climate simulations by imposing a deep forcing regime for $p>10$~bar: convergence time scale $τ_{\text{conv}}=10^6-10^8$~s to a common adiabat, as well as linear drag at depth ($p\geq 200$~bar), which mimics to first order magnetic drag. Lower boundary stability and the deep forcing assumptions were also tested with climate simulations for HD~209458b ($P_{\text{orb}}=P_{\text{rot}}=3.5$~days). HD~209458b simulations always show shallow wind jets (never deeper than 100~bar) and unperturbed superrotation. If we impose a fast rotation ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days), also the HD~209458b-like simulation shows equatorial retrograde flow at the day side. We conclude that the placement of the lower boundary at $p=200$~bar is justified for slow rotators like HD~209458b, but we suggest that it has to be placed deeper for fast-rotating, dense hot Jupiters ($P_{\text{orb}}\leq 1.5$~days) like WASP-43b. Our study highlights that the deep atmosphere may have a strong influence on the observable atmospheric flow in some hot Jupiters.

preprint2020arXiv

Evidence of Three Mechanisms Explaining the Radius Anomaly of Hot Jupiters

The radii of hot Jupiters are still not fully understood and all of the proposed explanations are based on the idea that these close-in giant planets possess hot interiors. We approach the radius anomaly problem by adopting a statistical approach. We infer the internal luminosity for the sample of hot Jupiters, study its effect on the interior structure, and put constraints on which mechanism is the dominant one. We develop a flexible and robust hierarchical Bayesian model that couples the interior structure of exoplanets to their observed properties. We apply the model to 314 hot Jupiters and infer the internal luminosity distribution for each planet and study at the population level ({\it i}) the mass-luminosity-radius distribution and as a function of equilibrium temperature the distributions of the ({\it ii}) heating efficiency, ({\it iii}) internal temperature, and the ({\it iv}) pressure of the radiative-convective-boundary (RCB). We find that hot Jupiters tend to have high internal luminosity leading to hot interiors. This has important consequences on the cooling rate and we find that the RCB is located at low pressures. Assuming that the ultimate source of the extra heating is the irradiation from the host star, we illustrate that the heating efficiency follows a Gaussian distribution, in agreement with previous results. We discuss our findings in the context of the proposed heating mechanisms and illustrate that ohmic dissipation, advection of potential temperature, and thermal tides are in agreement with certain trends inferred from our analysis and thus all three models can explain aspects of the observations. We provide new insights on the interior structure of hot Jupiters and show that with our current knowledge it is still challenging to firmly identify the universal mechanism driving the inflated radii.

preprint2020arXiv

Following the TraCS of exoplanets with Pan-Planets: Wendelstein-1b and Wendelstein-2b

Hot Jupiters seem to get rarer with decreasing stellar mass. The goal of the Pan-Planets transit survey was the detection of such planets and a statistical characterization of their frequency. Here, we announce the discovery and validation of two planets found in that survey, Wendelstein-1b and Wendelstein-2b, which are two short-period hot Jupiters that orbit late K host stars. We validated them both by the traditional method of radial velocity measurements with the HIgh Resolution Echelle Spectrometer (HIRES) and the Habitable-zone Planet Finder (HPF) instruments and then by their Transit Color Signature (TraCS). We observed the targets in the wavelength range of $4000 - 24000$ Angström and performed a simultaneous multiband transit fit and additionally determined their thermal emission via secondary eclipse observations. Wendelstein-1b is a hot Jupiter with a radius of $1.0314_{-0.0061}^{+0.0061}$ $R_J$ and mass of $0.592_{-0.129}^{+0.165}$ $M_J$, orbiting a K7V dwarf star at a period of $2.66$ d, and has an estimated surface temperature of about $1727_{-90}^{+78}$ K. Wendelstein-2b is a hot Jupiter with a radius of $1.1592_{-0.0210}^{+0.0204}$ $R_J$ and a mass of $0.731_{-0.311}^{+0.541}$ $M_J$, orbiting a K6V dwarf star at a period of $1.75$ d, and has an estimated surface temperature of about $1852_{-140}^{+120}$ K. With this, we demonstrate that multiband photometry is an effective way of validating transiting exoplanets, in particular for fainter targets since radial velocity (RV) follow-up becomes more and more costly for those targets.

preprint2020arXiv

HCN production in Titan's Atmosphere: Coupling quantum chemistry and disequilibrium atmospheric modeling

Hydrogen cyanide (HCN) is a critical reactive source of nitrogen for building key biomolecules relevant for the origin of life. Still, many HCN reactions remain uncharacterized by experiments and theory, and the complete picture of HCN production in planetary atmospheres is not fully understood. To improve this situation, we develop a novel technique making use of computational quantum chemistry, experimental data, and atmospheric numerical simulations. First, we use quantum chemistry simulations to explore the entire field of possible reactions for a list of primary species in N2-, CH4-, and H2-dominated atmospheres. In this process, we discover 33 new reactions with no previously known rate coefficients. From here, we develop a consistent reduced atmospheric hybrid chemical network (CRAHCN) containing experimental values when available, and our calculated rate coefficients otherwise. Next, we couple CRAHCN to a 1D chemical kinetic model (ChemKM) to compute the HCN abundance as a function of atmospheric depth on Titan. Our simulated atmospheric HCN profile agrees very well with the Cassini observations. CRAHCN contains 104 reactions however nearly all of the simulated atmospheric HCN profile can be obtained using a scaled down network of only 19 dominant reactions. From here, we form a complete picture of HCN chemistry in Titan's atmosphere, from the dissociation of the main atmospheric species, down to the direct production of HCN along 4 major channels. One of these channels was first discovered and characterized in Pearce et al. (2019) and this work.

preprint2020arXiv

Ice coverage of dust grains in cold astrophysical environments

Surface processes on cosmic solids in cold astrophysical environments lead to gas phase depletion and molecular complexity. Most astrophysical models assume that the molecular ice forms a thick multilayer substrate, not interacting with the dust surface. In contrast, we present experimental results demonstrating the importance of the surface for porous grains. We show that cosmic dust grains may be covered by a few monolayers of ice only. This implies that the role of dust surface structure, composition, and reactivity in models describing surface processes in cold interstellar, protostellar, and protoplanetary environments has to be reevaluated.

preprint2020arXiv

Indications for very high metallicity and absence of methane for the eccentric exo-Saturn WASP-117b

We investigate the atmospheric composition of the long period ($P_{\rm orb}=$ 10 days), eccentric exo-Saturn WASP-117b. WASP-117b could be in atmospheric temperature and chemistry similar to WASP-107b. In mass and radius WASP-117b is similar to WASP-39b, which allows a comparative study of these planets. We analyze a near-infrared transmission spectrum of WASP-117b taken with Hubble Space Telescope/WFC3 G141, which was reduced with two independent pipelines. High resolution measurements were taken with VLT/ESPRESSO in the optical. We report the robust ($3σ$) detection of a water spectral feature. Using a 1D atmosphere model with isothermal temperature, uniform cloud deck and equilibrium chemistry, the Bayesian evidence of a retrieval analysis of the transmission spectrum indicates a preference for a high atmospheric metallicity ${\rm [Fe/H]}=2.58^{+0.26}_{-0.37}$ and clear skies. The data are also consistent with a lower-metallicity composition ${\rm [Fe/H]}<1.75$ and a cloud deck between $10^{-2.2} - 10^{-5.1}$ bar, but with weaker Bayesian preference. We retrieve a low CH$_4$ abundance of $<10^{-4}$ volume fraction within $1 σ$ and $<2\cdot 10^{-1}$ volume fraction within $3 σ$. We cannot constrain the equilibrium temperature between theoretically imposed limits of 700 and 1000~K. Further observations are needed to confirm quenching of CH$_4$ with $K_{zz}\geq 10^8$~cm$^2$/s. We report indications of Na and K in the VLT/ESPRESSO high resolution spectrum with substantial Bayesian evidence in combination with HST data.

preprint2020arXiv

NGTS-11 b / TOI-1847 b: A transiting warm Saturn recovered from a TESS single-transit event

We report the discovery of NGTS-11 b (=TOI-1847 b), a transiting Saturn in a 35.46-day orbit around a mid K-type star (Teff=5050 K). We initially identified the system from a single-transit event in a TESS full-frame image light-curve. Following seventy-nine nights of photometric monitoring with an NGTS telescope, we observed a second full transit of NGTS-11 b approximately one year after the TESS single-transit event. The NGTS transit confirmed the parameters of the transit signal and restricted the orbital period to a set of 13 discrete periods. We combined our transit detections with precise radial velocity measurements to determine the true orbital period and measure the mass of the planet. We find NGTS-11 b has a radius of 0.817+0.028-0.032 $R_J$, a mass of 0.344+0.092-0.073 $M_J$, and an equilibrium temperature of just 435+34-32 K, making it one of the coolest known transiting gas giants. NGTS-11 b is the first exoplanet to be discovered after being initially identified as a TESS single-transit event, and its discovery highlights the power of intense photometric monitoring in recovering longer-period transiting exoplanets from single-transit events.

preprint2020arXiv

Pebble-driven Planet Formation around Very Low-mass Stars and Brown Dwarfs

We conduct a pebble-driven planet population synthesis study to investigate the formation of planets around very low-mass stars and brown dwarfs, in the (sub)stellar mass range between $0.01 \ M_{\odot}$ and $0.1 \ M_{\odot}$. Based on the extrapolation of numerical simulations of planetesimal formation by the streaming instability, we obtain the characteristic mass of the planetesimals and the initial masses of the protoplanets (largest bodies from the planetesimal size distributions), in either the early self-gravitating phase or the later non-self-gravitating phase of the protoplanetary disk evolution. We find that the initial protoplanets form with masses that increase with host mass, orbital distance and decrease with disk age. Around late M-dwarfs of $0.1 \ M_{\odot}$, these protoplanets can grow up to Earth-mass planets by pebble accretion. However, around brown dwarfs of $0.01 \ M_{\odot}$, planets do not grow larger than Mars mass when the initial protoplanets are born early in self-gravitating disks, and their growth stalls at around $0.01$ Earth-mass when they are born late in non-self-gravitating disks. Around these low mass stars and brown dwarfs, we find no channel for gas giant planet formation because the solid cores remain too small. When the initial protoplanets form only at the water-ice line, the final planets typically have ${\gtrsim} 15\%$ water mass fraction. Alternatively, when the initial protoplanets form log-uniformly distributed over the entire protoplanetary disk, the final planets are either very water-rich (water mass fraction ${\gtrsim}15\%$) or entirely rocky (water mass fraction ${\lesssim}5\%$).

preprint2020arXiv

Pebbles in an Embedded Protostellar Disk: The Case of CB26

Planetary cores are thought to form in proto-planetary disks via the growth of dusty solid material. However, it is unclear how early this process begins. We study the physical structure and grain growth in the edge-on disk that surrounds the ~1 Myr old low-mass (~0.55 Msun) protostar embedded in the Bok Globule CB26 to examine how much grain growth has already occurred in the protostellar phase. We combine the SED between 0.9 $μ$m and 6.4 cm with high angular resolution continuum maps at 1.3, 2.9, and 8.1 mm, and use the radiative transfer code RADMC-3D to conduct a detailed modelling of the dust emission from the disk and envelope of CB 26. We infer inner and outer disk radii of around 16 au and 172$\pm$22 au, respectively. The total gas mass in the disk is ~0.076 Msun, which amounts to ~14% of the mass of the central star. The inner disk contains a compact free-free emission region, which could be related to either a jet or a photoevaporation region. The thermal dust emission from the outer disk is optically thin at mm wavelengths, while the emission from the inner disk midplane is moderately optically thick. Our best-fit radiative transfer models indicate that the dust grains in the disk have already grown to pebbles with diameters of the order of 10 cm in size. Residual 8.1 mm emission suggests the presence of even larger particles in the inner disk. For the optically thin mm dust emission from the outer disk, we derive a mean opacity slope of 0.6$\pm$0.4, which is consistent with the presence of large dust grains. The presence of cm-sized bodies in the CB 26 disk indicates that solids grow rapidly already during the first million years in a protostellar disk. It is thus possible that Class II disks are already seeded with large particles and may contain even planetesimals.

preprint2020arXiv

SCExAO/CHARIS Near-IR Integral Field Spectroscopy of the HD 15115 Debris Disk

We present new, near-infrared ($1.1 - 2.4$ $μm$) high-contrast imaging of the debris disk around HD 15115 with the Subaru Coronagraphic Extreme Adaptive Optics system (SCExAO) coupled with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS). SCExAO/CHARIS resolves the disk down to $ρ\sim 0.2''$ ($\rm{r_{proj}} \sim 10$ $\rm{au}$), a factor of $\sim 3-5$ smaller than previous recent studies. We derive a disk position angle of $\rm{PA}$ $\sim 279.4^\circ - 280.5^\circ$ and an inclination of $\rm{i}$ $\sim 85.3^\circ - 86.2^\circ$. While recent SPHERE/IRDIS imagery of the system could suggest a significantly misaligned two ring disk geometry, CHARIS imagery does not reveal conclusive evidence for this hypothesis. Moreover, optimizing models of both one and two ring geometries using differential evolution, we find that a single ring having a Hong-like scattering phase function matches the data equally well within the CHARIS field of view ($ρ\lesssim 1''$). The disk's asymmetry, well-evidenced at larger separations, is also recovered; the west side of the disk appears on average around 0.4 magnitudes brighter across the CHARIS bandpass between $0.25''$ and $1''$. Comparing STIS/50CCD optical photometry ($2000-10500$ $Å$) with CHARIS NIR photometry, we find a red (STIS/50CCD$-$CHARIS broadband) color for both sides of the disk throughout the $0.4'' - 1''$ region of overlap, in contrast to the blue color reported at similar wavelengths for regions exterior to $\sim 2''$. Further, this color may suggest a smaller minimum grain size than previously estimated at larger separations. Finally, we provide constraints on planetary companions, and discuss possible mechanisms for the observed inner disk flux asymmetry and color.

preprint2020arXiv

Sensitivity of gas-grain chemical models to surface reaction barriers: Effect from a key carbon-insertion reaction, C + H$_2$ $\rightarrow$ CH$_2$

The feasibility of contemporary gas-grain astrochemical models depends on the availability of accurate kinetics data, in particular, for surface processes. We study the sensitivity of gas-grain chemical models to the energy barrier Ea of the important surface reaction between some of the most abundant species: C and H2 (surface C + surface H2 = surface CH2). We used the gas-grain code ALCHEMIC to model the time-dependent chemical evolution over a 2D grid of densities (nH: 10^{3} - 10^{12}cm^{-3}) and temperatures (T: 10 - 300 K), assuming UV-dark (Av = 20 mag) and partly UV-irradiated (Av = 3 mag) conditions that are typical of the dense interstellar medium. We considered two values for the energy barrier of the surface reaction, Ea = 2500 K (as originally implemented in the networks) and Ea = 0 K (as measured in the laboratory and computed by quantum chemistry simulations). We find that if the C + H2 = CH2 surface reaction is barrierless, a more rapid conversion of the surface carbon atoms into methane ice occurs. Overproduction of the CHn hydrocarbon ices affects the surface formation of more complex hydrocarbons, cyanides and nitriles, and CS-bearing species at low temperatures < 10-15 K. The surface hydrogenation of CO and hence the synthesis of complex (organic) molecules become affected as well. As a result, important species whose abundances may change by more than a factor of two at 1 Myr include atomic carbon, small mono-carbonic (C1) and di-carbonic (C2) hydrocarbons, CO2, CN, HCN, HNC, HNCO, CS, H2CO, H2CS, CH2CO, and CH3OH (in either gas and/or ice). The abundances of key species, CO, H2O, and N2 as well as O, HCO+, N2H+, NH3, NO, and most of the S-bearing molecules, remain almost unaffected. Further accurate laboratory measurements and quantum chemical calculations of the surface reaction barriers will be crucial to improve the accuracy of astrochemical models.

preprint2020arXiv

Separate silicate and carbonaceous solids formed from mixed atomic and molecular species diffusing in neon ice

The formation and growth of refractory matter on pre-existing interstellar dust grain surfaces was studied experimentally by annealing neon-ice matrices in which potential precursors of silicate grains (Mg and Fe atoms, SiO and SiO$_2$ molecules) and of solid carbon (C$_n$ molecules, $n$ = 2-10) were initially isolated. Other molecules, mainly O$_3$, CO, CO$_2$, C$_3$O, and H$_2$O, were embedded at the same time in the matrices. The annealing procedure caused the cold dopants to diffuse and interact in the neon ice. Monitoring the procedure in situ with infrared spectroscopy revealed the disappearance of the silicon oxide and carbon molecules at temperatures lower than 13 K, and the rise of the Si-O stretching band of silicates. Ex situ electron microscopy confirmed the formation of silicate grains and showed that their structure was amorphous. It also showed that amorphous carbon matter was formed simultaneously next to the silicate grains, the two materials being chemically separated. The results of the experiments support the hypothesis that grains of complex silicates and of carbonaceous materials are re-formed in the cold ISM, as suggested by astronomical observations and evolution models of cosmic dust masses. Moreover, they show that the potential precursors of one material do not combine with those of the other at cryogenic temperatures, providing us with a clue as to the separation of silicates and carbon in interstellar grains.

preprint2020arXiv

Solving grain size inconsistency between ALMA polarization and VLA continuum in the Ophiuchus IRS 48 protoplanetary disk

The protoplanetary disk around Ophiuchus IRS 48 shows an azimuthally asymmetric dust distribution in (sub-)millimeter observations, which is interpreted as a vortex, where millimeter/centimeter-sized particles are trapped at the location of the continuum peak. In this paper, we present 860 $μ$m ALMA observations of polarized dust emission of this disk. The polarized emission was detected toward a part of the disk. The polarization vectors are parallel to the disk minor axis, and the polarization fraction was derived to be $1-2$\%. These characteristics are consistent with models of self-scattering of submillimeter-wave emission, which indicate a maximum grain size of $\sim100$ $μ$m. However, this is inconsistent with the previous interpretation of millimeter/centimeter dust particles being trapped by a vortex. To explain both, ALMA polarization and previous ALMA and VLA observations, we suggest that the thermal emission at 860 $μ$m wavelength is optically thick ($τ_{\rm abs}\sim7.3$) at the dust trap with the maximum observable grain size of $\sim100$ $μ$m rather than an optically thin case with $\sim$ cm dust grains. We note that we cannot rule out that larger dust grains are accumulated near the midplane if the 860 $μ$m thermal emission is optically thick.

preprint2020arXiv

Spitzer Variability Properties of Low-Gravity L Dwarfs

We present \textit{Spitzer Space Telescope} variability monitoring observations of three low-gravity L dwarfs with previous detections of variability in the near-IR, 2MASS J0045+16, 2MASS J0501-00 and 2MASS J1425-36. We detect significant, periodic variability in two of our targets, 2MASS J0045+16 and 2MASS J0501-00. We do not detect variability in 2MASS J1425-36. Combining our new rotation periods with rotational velocities, we calculate inclination angles of $22\pm1^{\circ}$, ${60^{+13 }_{-8}} ^{\circ}$ and $52^{+19}_{-13}~^{\circ}$ for 2MASS J0045+16, 2MASS J0501-00 and 2MASS J1425-36 respectively. Our three new objects are consistent with the tentative relations between inclination, amplitude and color anomaly previously reported. Objects with the highest variability amplitudes are inclined equator-on, while the maximum observed amplitude decreases as the inclination angle decreases. We also find a correlation between the inclination angle and $(J-K)_{\mathrm{2MASS}}$ color anomaly for the sample of objects with measured inclinations. Compiling the entire sample of brown dwarfs with \textit{Spitzer} variability detections, we find no enhancement in amplitude for young, early-L dwarfs compared to the field dwarf population. We find a possible enhancement in amplitude of low-gravity late-L dwarfs at $4.5~μ$m. We do not find a correlation between amplitude ratio and spectral type for field dwarfs or for the young population. Finally, we compile the rotation periods of a large sample of brown dwarfs with ages 1 Myr to 1 Gyr and compare the rotation rates predicted by evolutionary models assuming angular momentum conservation. We find that the rotation rates of the current sample of brown dwarfs fall within the expected range set by evolutionary models and breakup limits.

preprint2020arXiv

Star-Gas Surface Density Correlations in Twelve Nearby Molecular Clouds I: Data Collection and Star-Sampled Analysis

We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects from the most recent $Spitzer$ Extended Solar Neighborhood Archive (SESNA) catalog to probe the stellar surface density. Using a star-sampled nearest neighbor technique to probe the star-gas surface density correlations at the scale of a few parsecs, we find that the stellar mass surface density varies as a power-law of the gas mass surface density, with a power-law index of $\sim$2 in all the clouds. The consistent power-law index implies that star formation efficiency is directly correlated with gas column density, and no gas column density threshold for star formation is observed. We compare the observed correlations with the predictions from an analytical model of thermal fragmentation, and with the synthetic observations of a recent hydrodynamic simulation of a turbulent star-forming molecular cloud. We find that the observed correlations are consistent for some clouds with the thermal fragmentation model and can be reproduced using the hydrodynamic simulations.

preprint2020arXiv

The Herschel Orion Protostar Survey: Far-Infrared Photometry and Colors of Protostars and Their Variations across Orion A and B

The degree to which the properties of protostars are affected by environment remains an open question. To investigate this, we look at the Orion A and B molecular clouds, home to most of the protostars within 500 pc. At ~400 pc, Orion is close enough to distinguish individual protostars across a range of environments in terms of both the stellar and gas projected densities. As part of the Herschel Orion Protostar Survey (HOPS), we used the Photodetector Array Camera and Spectrometer (PACS) to map 108 partially overlapping square fields with edge lengths of 5 arcmin or 8 arcmin and measure the 70 micron and 160 micron flux densities of 338 protostars within them. In this paper we examine how these flux densities and their ratio depend on evolutionary state and environment within the Orion complex. We show that Class 0 protostars occupy a region of the 70 micron flux density versus 160 micron to 70 micron flux density ratio diagram that is distinct from their more evolved counterparts. We then present evidence that the Integral-Shaped Filament (ISF) and Orion B contain protostars with more massive envelopes than those in the more sparsely populated LDN 1641 region. This can be interpreted as evidence for increasing star formation rates in the ISF and Orion B or as a tendency for more massive envelopes to be inherited from denser birth environments. We also provide technical details about the map-making and photometric procedures used in the HOPS program.

preprint2020arXiv

The multi-planet system TOI-421 -- A warm Neptune and a super puffy mini-Neptune transiting a G9 V star in a visual binary

We report the discovery of a warm Neptune and a hot sub-Neptune transiting TOI-421 (BD-14 1137, TIC 94986319), a bright (V=9.9) G9 dwarf star in a visual binary system observed by the TESS space mission in Sectors 5 and 6. We performed ground-based follow-up observations -- comprised of LCOGT transit photometry, NIRC2 adaptive optics imaging, and FIES, CORALIE, HARPS, HIRES, and PFS high-precision Doppler measurements -- and confirmed the planetary nature of the 16-day transiting candidate announced by the TESS team. We discovered an additional radial velocity signal with a period of 5 days induced by the presence of a second planet in the system, which we also found to transit its host star. We found that the inner mini-Neptune, TOI-421b, has an orbital period of Pb =5.19672 +- 0.00049 days, a mass of Mb = 7.17 +- 0.66 Mearth and a radius of Rb = 2.68+0.19-0.18 Rearth, whereas the outer warm Neptune, TOI-421 c, has a period of Pc =16.06819 +- 0.00035 days, a mass of Mc = 16.42+1.06-1.04 Mearth, a radius of Rc = 5.09+0.16-0.15 Rearth and a density of rho_c =0.685+0.080-0.072 g cm-3 . With its characteristics the inner planet (rho_b=2.05+0.52-0.41 g cm-3) is placed in the intriguing class of the super-puffy mini-Neptunes. TOI-421b and TOI-421c are found to be well suitable for atmospheric characterization. Our atmospheric simulations predict significant Ly-alpha transit absorption, due to strong hydrogen escape in both planets, and the presence of detectable CH_4 in the atmosphere of TOI-421c if equilibrium chemistry is assumed.

preprint2020arXiv

The role of clouds on the depletion of methane and water dominance in the transmission spectra of irradiated exoplanets

Observations suggest an abundance of water and paucity of methane in the majority of observed exoplanetary atmospheres. We isolate the effect of atmospheric processes to investigate possible causes. Previously, we studied the effect of effective temperature, surface gravity, metallicity, carbon-to-oxygen ratio, and stellar type assuming cloud-free thermochemical equilibrium and disequilibrium chemistry. However, under these assumptions, methane remains a persisting spectral feature in the transmission spectra of exoplanets over a certain parameter space, the Methane Valley. In this work we investigate the role of clouds on this domain and we find that clouds change the spectral appearance of methane in two direct ways: 1) by heating-up the photosphere of colder planets, and 2) by obscuring molecular features. The presence of clouds also affects methane features indirectly: 1) cloud heating results in more evaporation of condensates and hence releases additional oxygen, causing water dominated spectra of colder carbon-poor exoplanets, and 2) HCN/CO production results in a suppression of depleted methane features by these molecules. The presence of HCN/CO and a lack of methane could be an indication of cloud formation on hot exoplanets. Cloud heating can also deplete ammonia. Therefore, a simultaneous depletion of methane and ammonia is not unique to photochemical processes. We propose that the best targets for methane detection are likely to be massive but smaller planets with a temperature around 1450 K orbiting colder stars. We also construct Spitzer synthetic color-maps and find that clouds can explain some of the high contrast observations by IRAC's channel 1 and 2.

preprint2020arXiv

The surprisingly low carbon mass in the debris disk around HD 32297

Gas has been detected in a number of debris disks. It is likely secondary, i.e. produced by colliding solids. Here, we report ALMA Band 8 observations of neutral carbon in the CO-rich debris disk around the 15--30 Myr old A-type star HD 32297. We find that C$^0$ is located in a ring at $\sim$110 au with a FWHM of $\sim$80 au, and has a mass of $(3.5\pm0.2)\times10^{-3}$ M$_\oplus$. Naively, such a surprisingly small mass can be accumulated from CO photo-dissociation in a time as short as $\sim$10$^4$ yr. We develop a simple model for gas production and destruction in this system, properly accounting for CO self-shielding and shielding by neutral carbon, and introducing a removal mechanism for carbon gas. We find that the most likely scenario to explain both C$^0$ and CO observations, is one where the carbon gas is rapidly removed on a timescale of order a thousand years and the system maintains a very high CO production rate of $\sim$15 M$_\oplus$ Myr$^{-1}$, much higher than the rate of dust grind-down. We propose a possible scenario to meet these peculiar conditions: the capture of carbon onto dust grains, followed by rapid CO re-formation and re-release. In steady state, CO would continuously be recycled, producing a CO-rich gas ring that shows no appreciable spreading over time. This picture might be extended to explain other gas-rich debris disks.

preprint2020arXiv

Thermal formation of ammonium carbamate on the surface of laboratory analogues of carbonaceous grains in protostellar envelopes and planet-forming disks

The catalytic role of dust grain surfaces in the thermal reaction CO2 + 2NH3 $\rightarrow$ NH4+NH2COO was recently demonstrated by our group. The rate coefficients for the reaction at 80 K on the surface of nanometre-sized carbon and silicate grains were measured to be up to three times higher compared to the reaction rate coefficients measured on KBr. In this study, the reaction was performed on carbon grains and on KBr in the extended temperature range of 50 - 80 K and with the addition of water ice. The reaction activation energy was found to be about 3 times lower on grains compared to the corresponding ice layer on KBr. Thus, the catalytic role of the dust grain surface in the studied reaction can be related to a reduction of the reaction barrier. Addition of water to NH3:CO2 ice on grains slowed the reaction down. At the H2O:CO2 ratio of 5:1, the reaction was not detected on the experimental timescale. This result calls into question the thermal formation of ammonium carbamate in dense molecular clouds and outer regions of protostellar and protoplanetary environments with dominating water ice mantle chemistry. However, it can still happen in inner regions of protostellar and protoplanetary environments in crystalline ices.

preprint2020arXiv

TOI 694 b and TIC 220568520 b: Two Low-Mass Companions Near the Hydrogen Burning Mass Limit Orbiting Sun-like Stars

We report the discovery of TOI 694 b and TIC 220568520 b, two low-mass stellar companions in eccentric orbits around metal-rich Sun-like stars, first detected by the Transiting Exoplanet Survey Satellite (TESS). TOI 694 b has an orbital period of 48.05131$\pm$0.00019 days and eccentricity of 0.51946$\pm$0.00081, and we derive a mass of 89.0$\pm$5.3 $M_J$ (0.0849$\pm$0.0051 $M_\odot$) and radius of 1.111$\pm$0.017 $R_J$ (0.1142$\pm$0.0017 $R_\odot$). TIC 220568520 b has an orbital period of 18.55769$\pm$0.00039 days and eccentricity of 0.0964$\pm$0.0032, and we derive a mass of 107.2$\pm$5.2 $M_J$ (0.1023$\pm$0.0050 $M_\odot$) and radius of 1.248$\pm$0.018 $R_J$ (0.1282$\pm$0.0019 $R_\odot$). Both binary companions lie close to and above the Hydrogen burning mass threshold that separates brown dwarfs and the lowest mass stars, with TOI 694 b being 2-$σ$ above the canonical mass threshold of 0.075 $M_\odot$. The relatively long periods of the systems mean that the magnetic fields of the low-mass companions are not expected to inhibit convection and inflate the radius, which according to one leading theory is common in similar objects residing in short-period tidally-synchronized binary systems. Indeed we do not find radius inflation for these two objects when compared to theoretical isochrones. These two new objects add to the short but growing list of low-mass stars with well-measured masses and radii, and highlight the potential of the TESS mission for detecting such rare objects orbiting bright stars.

preprint2020arXiv

Two intermediate-mass transiting brown dwarfs from the TESS mission

We report the discovery of two intermediate-mass brown dwarfs (BDs), TOI-569b and TOI-1406b, from NASA's Transiting Exoplanet Survey Satellite mission. TOI-569b has an orbital period of $P = 6.55604 \pm 0.00016$ days, a mass of $M_b = 64.1 \pm 1.9 M_J$, and a radius of $R_b = 0.75 \pm 0.02 R_J$. Its host star, TOI-569, has a mass of $M_\star = 1.21 \pm 0.03 M_\odot$, a radius of $R_\star = 1.47 \pm 0.03 R_\odot$, $\rm [Fe/H] = +0.29 \pm 0.09$ dex, and an effective temperature of $T_{\rm eff} = 5768 \pm 110K$. TOI-1406b has an orbital period of $P = 10.57415 \pm 0.00063$ days, a mass of $M_b =46.0 \pm 2.7 M_J$, and a radius of $R_b = 0.86 \pm 0.03 R_J$. The host star for this BD has a mass of $M_\star =1 .18 \pm 0.09 M_\odot$, a radius of $R_\star = 1.35 \pm 0.03 R_\odot$, $ \rm [Fe/H] = -0.08 \pm 0.09$ dex and an effective temperature of $T_{\rm eff} = 6290 \pm 100K$. Both BDs are in circular orbits around their host stars and are older than 3 Gyr based on stellar isochrone models of the stars. TOI-569 is one of two slightly evolved stars known to host a transiting BD (the other being KOI-415). TOI-1406b is one of three known transiting BDs to occupy the mass range of $40-50 M_J$ and one of two to have a circular orbit at a period near 10 days (with the first being KOI-205b).Both BDs have reliable ages from stellar isochrones in addition to their well-constrained masses and radii, making them particularly valuable as tests for substellar isochrones in the BD mass-radius diagram.

preprint2020arXiv

Two transiting hot Jupiters from the WASP survey: WASP-150b and WASP-176b

We report the discovery of two transiting exoplanets from the WASP survey, WASP-150b and WASP-176b. WASP-150b is an eccentric ($e$ = 0.38) hot Jupiter on a 5.6 day orbit around a $V$ = 12.03, F8 main-sequence host. The host star has a mass and radius of 1.4 $\rm M_{\odot}$ and 1.7 $\rm R_{\odot}$ respectively. WASP-150b has a mass and radius of 8.5 $\rm M_J$ and 1.1 $\rm R_J$, leading to a large planetary bulk density of 6.4 $\rm ρ_J$. WASP-150b is found to be $\sim3$ Gyr old, well below its circularisation timescale, supporting the eccentric nature of the planet. WASP-176b is a hot Jupiter planet on a 3.9 day orbit around a $V$ = 12.01, F9 sub-giant host. The host star has a mass and radius of 1.3 $\rm M_{\odot}$ and 1.9 $\rm R_{\odot}$. WASP-176b has a mass and radius of 0.86 $\rm M_J$ and 1.5 $\rm R_J$ respectively, leading to a planetary bulk density of 0.23 $\rm ρ_J$.

preprint2019arXiv

HD 213885b: A transiting 1-day-period super-Earth with an Earth-like composition around a bright ($V=7.9$) star unveiled by TESS

We report the discovery of the 1.008-day, ultra-short period (USP) super-Earth HD 213885b (TOI-141b) orbiting the bright ($V=7.9$) star HD 213885 (TOI-141, TIC 403224672), detected using photometry from the recently launched TESS mission. Using FEROS, HARPS and CORALIE radial-velocities, we measure a precise mass of $8.8\pm0.6$ $M_\oplus$ for this $1.74 \pm 0.05$ $R_\oplus$ exoplanet, which provides enough information to constrain its bulk composition, which is similar to Earth's but enriched in iron. The radius, mass and stellar irradiation of HD 213885b are, given our data, very similar to 55 Cancri e, making this exoplanet a good target to perform comparative exoplanetology of short period, highly irradiated super-Earths. Our precise radial-velocities reveal an additional $4.78$-day signal which we interpret as arising from a second, non-transiting planet in the system, HD 213885c (TOI-141c), whose minimum mass of $19.95\pm 1.4$ $M_\oplus$ makes it consistent with being a Neptune-mass exoplanet. The HD 213885 system is very interesting from the perspective of future atmospheric characterization, being the second brightest star to host an ultra-short period transiting super-Earth (with the brightest star being, in fact, 55 Cancri). Prospects for characterization with present and future observatories are discussed.

preprint2019arXiv

High-Resolution Near-Infrared Polarimetry and Sub-Millimeter Imaging of FS Tau A: Possible Streamers in Misaligned Circumbinary Disk System

We analyzed the young (2.8-Myr-old) binary system FS Tau A using near-infrared (H-band) high-contrast polarimetry data from Subaru/HiCIAO and sub-millimeter CO (J=2-1) line emission data from ALMA. Both the near-infrared and sub-millimeter observations reveal several clear structures extending to $\sim$240 AU from the stars. Based on these observations at different wavelengths, we report the following discoveries. One arm-like structure detected in the near-infrared band initially extends from the south of the binary with a subsequent turn to the northeast, corresponding to two bar-like structures detected in ALMA observations with an LSRK velocity of 1.19-5.64 km/s. Another feature detected in the near-infrared band extends initially from the north of the binary, relating to an arm-like structure detected in ALMA observations with an LSRK velocity of 8.17-16.43 km/s. From their shapes and velocities, we suggest that these structures can mostly be explained by two streamers that connect the outer circumbinary disk and the central binary components. These discoveries will be helpful for understanding the evolution of streamers and circumstellar disks in young binary systems.

preprint2019arXiv

TOI-222: a single-transit TESS candidate revealed to be a 34-day eclipsing binary with CORALIE, EulerCam and NGTS

We report the period, eccentricity, and mass determination for the TESS single-transit event candidate TOI-222, which displayed a single 3000 ppm transit in the TESS two-minute cadence data from Sector 2. We determine the orbital period via radial velocity measurements (P=33.9,days), which allowed for ground-based photometric detection of two subsequent transits. Our data show that the companion to TOI-222 is a low mass star, with a radius of $0.18_{-0.10}^{+0.39}$ Rsun and a mass of $0.23\pm0.01$ Msun. This discovery showcases the ability to efficiently discover long-period systems from TESS single transit events using a combination of radial velocity monitoring coupled with high precision ground-based photometry.

preprint2019arXiv

TOI-677 b: A Warm Jupiter (P=11.2d) on an eccentric orbit transiting a late F-type star

We report the discovery of TOI-677 b, first identified as a candidate in light curves obtained within Sectors 9 and 10 of the Transiting Exoplanet Survey Satellite (TESS) mission and confirmed with radial velocities. TOI-677 b has a mass of M_p = 1.236$^{+0.069}_{-0.067}$ M_J, a radius of R_p = 1.170 +- 0.03 R_J,and orbits its bright host star (V=9.8 mag) with an orbital period of 11.23660 +- 0.00011 d, on an eccentric orbit with e = 0.435 +- 0.024. The host star has a mass of M_* = 1.181 +- 0.058 M_sun, a radius of R_* = 1.28 +- 0.03 R_sun, an age of 2.92$^{+0.80}_{-0.73}$ Gyr and solar metallicity, properties consistent with a main sequence late F star with T_eff = 6295 +- 77 K. We find evidence in the radial velocity measurements of a secondary long term signal which could be due to an outer companion. The TOI-677 b system is a well suited target for Rossiter-Mclaughlin observations that can constrain migration mechanisms of close-in giant planets.

preprint2016arXiv

197 Candidates and 104 Validated Planets in K2's First Five Fields

We present 197 planet candidates discovered using data from the first year of the NASA K2 mission (Campaigns 0-4), along with the results of an intensive program of photometric analyses, stellar spectroscopy, high-resolution imaging, and statistical validation. We distill these candidates into sets of 104 validated planets (57 in multi-planet systems), 30 false positives, and 63 remaining candidates. Our validated systems span a range of properties, with median values of R_P = 2.3 R_E, P=8.6 d, Tef = 5300 K, and Kp=12.7 mag. Stellar spectroscopy provides precise stellar and planetary parameters for most of these systems. We show that K2 has increased by 30% the number of small planets known to orbit moderately bright stars (1-4 R_E, Kp=9-13 mag). Of particular interest are 37 planets smaller than 2 R_E, 15 orbiting stars brighter than Kp=11.5, five receiving Earth-like irradiation levels, and several multi-planet systems -- including four planets orbiting the M dwarf K2-72 near mean-motion resonances. By quantifying the likelihood that each candidate is a planet we demonstrate that our candidate sample has an overall false positive rate of 15-30%, with rates substantially lower for small candidates (< 2 R_E) and larger for candidates with radii > 8 R_E and/or with P < 3 d. Extrapolation of the current planetary yield suggests that K2 will discover between 500-1000 planets in its planned four-year mission -- assuming sufficient follow-up resources are available. Efficient observing and analysis, together with an organized and coherent follow-up strategy, is essential to maximize the efficacy of planet-validation efforts for K2, TESS, and future large-scale surveys.

preprint2016arXiv

A Census of Large-Scale ($\ge$ 10 pc), Velocity-Coherent, Dense Filaments in the Northern Galactic Plane: Automated Identification Using Minimum Spanning Tree

Large-scale gaseous filaments with length up to the order of 100 pc are on the upper end of the filamentary hierarchy of the Galactic interstellar medium. Their association with respect to the Galactic structure and their role in Galactic star formation are of great interest from both observational and theoretical point of view. Previous "by-eye" searches, combined together, have started to uncover the Galactic distribution of large filaments, yet inherent bias and small sample size limit conclusive statistical results to be drawn. Here, we present (1) a new, automated method to identify large-scale velocity-coherent dense filaments, and (2) the first statistics and the Galactic distribution of these filaments. We use a customized minimum spanning tree algorithm to identify filaments by connecting voxels in the position-position-velocity space, using the Bolocam Galactic Plane Survey spectroscopic catalog. In the range of $7.^{\circ}5 \le l \le 194^{\circ}$, we have identified 54 large-scale filaments and derived mass ($\sim 10^3 - 10^5 \, M_\odot$), length (10-276 pc), linear mass density (54-8625 $M_\odot \, \rm{pc}^{-1}$), aspect ratio, linearity, velocity gradient, temperature, fragmentation, Galactic location and orientation angle. The filaments concentrate along major spiral arms. They are widely distributed across the Galactic disk, with 50% located within $\pm$20 pc from the Galactic mid-plane and 27% run in the center of spiral arms (aka "bones"). An order of 1% of the molecular ISM is confined in large filaments. Massive star formation is more favorable in large filaments compared to elsewhere. This is the first comprehensive catalog of large filaments useful for a quantitative comparison with spiral structures and numerical simulations.

preprint2016arXiv

A Resolved Near-Infrared Image of The Inner Cavity in The GM Aur Transitional Disk

We present high-contrast H-band polarized intensity (PI) images of the transitional disk around the young solar-like star GM Aur. The near-infrared direct imaging of the disk was derived by polarimetric differential imaging using the Subaru 8.2-m Telescope and HiCIAO. An angular resolution and an inner working angle of 0."07 and r~0."05, respectively, were obtained. We clearly resolved a large inner cavity, with a measured radius of 18+/-2 au, which is smaller than that of a submillimeter interferometric image (28 au). This discrepancy in the cavity radii at near-infrared and submillimeter wavelengths may be caused by a 3-4M_Jup planet about 20 au away from the star, near the edge of the cavity. The presence of a near-infrared inner is a strong constraint on hypotheses for inner cavity formation in a transitional disk. A dust filtration mechanism has been proposed to explain the large cavity in the submillimeter image, but our results suggest that this mechanism must be combined with an additional process. We found that the PI slope of the outer disk is significantly different from the intensity slope obtained from HST/NICMOS, and this difference may indicate the grain growth process in the disk.

preprint2016arXiv

A Substellar Companion to Pleiades HII 3441

We find a new substellar companion to the Pleiades member star, Pleiades HII 3441, using the Subaru telescope with adaptive optics. The discovery is made as part of the high-contrast imaging survey to search for planetary-mass and substellar companions in the Pleiades and young moving groups. The companion has a projected separation of 0".49 +/- 0".02 (66 +/- 2 AU) and a mass of 68 +/- 5 M_J based on three observations in the J-, H-, and K_S-band. The spectral type is estimated to be M7 (~2700 K), and thus no methane absorption is detected in the H band. Our Pleiades observations result in the detection of two substellar companions including one previously reported among 20 observed Pleiades stars, and indicate that the fraction of substellar companions in the Pleiades is about 10.0 +26.1/-8.8 %. This is consistent with multiplicity studies of both the Pleiades stars and other open clusters.

preprint2016arXiv

C60 as a probe for astrophysical environments

The C60 molecule has been recently detected in a wide range of astrophysical environments through its four active intramolecular vibrational modes (T1u) near 18.9, 17.4, 8.5, and 7.0 microns. The strengths of the mid-infrared emission bands have been used to infer astrophysical conditions in the fullerene-rich regions. Widely varying values of the relative intrinsic strengths (RIS) of these four bands are reported in laboratory and theoretical papers, which impedes the derivation of the excitation mechanism of C60 in the astrophysical sources. The spectroscopic analysis of the C60 samples produced with our method delivers highly reproducible RIS values of 100, 25 +- 1, 26 +- 1, and 40 +- 4. A comparison of the inferred C60 emission band strengths with the astrophysical data shows that the observed strengths cannot be explained in terms of fluorescent or thermal emission alone. The large range in the observed 17.4/18.9 emission ratios indicates that either the emission bands contain significant contributions from emitters other than C60, or that the population distribution among the C60 vibrational modes is affected by physical processes other than thermal or UV excitation, such as chemo-luminescence from nascent C60 or possibly, Poincare fluorescence resulting from an inverse internal energy conversion. We have carefully analyzed the effect of the weakly-active fundamental modes and second order modes in the mid-infrared spectrum of C60 and propose that neutral C60 is the carrier of the unidentified emission band at 6.49 microns which has been observed in fullerene-rich environments.

preprint2016arXiv

Circumstellar Disks of the Most Vigorously Accreting Young Stars

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

preprint2016arXiv

Cold CO gas in the disk of the young eruptive star EX Lup

EX Lupi-type objects (EXors) form a sub-class of T Tauri stars, defined by sudden sporadic flare-ups of 1-5 magnitudes at optical wavelengths. These eruptions are attributed to enhanced mass accretion from the circumstellar disk to the star, and may constitute important events in shaping the structure of the inner disk and the forming planetary system. Although disk properties must play a fundamental role in driving the outbursts, they are surprisingly poorly known. In order to characterize the dust and gas components of EXor disks, here we report on observations of the $^{12}$CO J=3-2 and 4-3 lines, and the $^{13}$CO 3-2 line in EX Lup, the prototype of the EXor class. We reproduce the observed line fluxes and profiles with a line radiative transfer model, and compare the obtained parameters with corresponding ones of other T Tauri disks.

preprint2016arXiv

Constraining the Movement of the Spiral Features and the Locations of Planetary Bodies within the AB Aur System

We present new analysis of multi-epoch, H-band, scattered light images of the AB Aur system. We used a Monte Carlo, radiative transfer code to simultaneously model the system's SED and H-band polarized intensity imagery. We find that a disk-dominated model, as opposed to one that is envelope dominated, can plausibly reproduce AB Aur's SED and near-IR imagery. This is consistent with previous modeling attempts presented in the literature and supports the idea that at least a subset of AB Aur's spirals originate within the disk. In light of this, we also analyzed the movement of spiral structures in multi-epoch H-band total light and polarized intensity imagery of the disk. We detect no significant rotation or change in spatial location of the spiral structures in these data, which span a 5.8 year baseline. If such structures are caused by disk-planet interactions, the lack of observed rotation constrains the location of the orbit of planetary perturbers to be >47 AU.

preprint2016arXiv

Deep HST/STIS Visible-Light Imaging of Debris Systems around Solar Analog Hosts

We present new Hubble Space Telescope observations of three a priori known starlight-scattering circumstellar debris systems (CDSs) viewed at intermediate inclinations around nearby close-solar analog stars: HD 207129, HD 202628, and HD 202917. Each of these CDSs possesses ring-like components that are more-massive analogs of our solar system's Edgeworth- Kuiper belt. These systems were chosen for follow-up observations to provide higher-fidelity and better sensitivity imaging for the sparse sample of solar-analog CDSs that range over two decades in systemic ages with HD 202628 and HD 202917 (both ~ 2.3 Gyr) currently the oldest CDSs imaged in visible or near-IR light. These deep (10 - 14 ksec) observations, with six-roll point-spread-function template subtracted visible-light coronagraphy using the Space Telescope Imaging Spectrograph, were designed to better reveal their angularly large, diffuse/low surface brightness, debris rings, and for all targets probe their exo-ring environments for starlight-scattering materials that present observational challenges for current ground-based facilities and instruments. Contemporaneously also observing with a narrower occulter position, these observations additionally probe the CDS endo-ring environments seen to be relatively devoid of scatterers. We discuss the morphological, geometrical, and photometric properties of these CDSs also in the context of other FGK-star hosted CDSs we have previously imaged as a homogeneously observed ensemble. From this combined sample we report a general decay in quiescent disk F_disk/F_star optical brightness ~ t^-0.8, similar to what is seen in at thermal IR wavelengths, and CDSs with a significant diversity in scattering phase asymmetries, and spatial distributions of their starlight-scattering grains.

preprint2016arXiv

Extreme Asymmetry in the Disk of V1247 Ori

We present the first near-infrared scattered-light detection of the transitional disk around V1247 Ori, which was obtained using high-resolution polarimetric differential imaging observations with Subaru/HiCIAO. Our imaging in the H band reveals the disk morphology at separations of ~0.14"-0.86" (54-330 au) from the central star. The polarized intensity (PI) image shows a remarkable arc-like structure toward the southeast of the star, whereas the fainter northwest region does not exhibit any notable features. The shape of the arm is consistent with an arc of 0.28" $\pm$ 0.09" in radius (108 au from the star), although the possibility of a spiral arm with a small pitch angle cannot be excluded. V1247 Ori features an exceptionally large azimuthal contrast in scattered, polarized light; the radial peak of the southeastern arc is about three times brighter than the northwestern disk measured at the same distance from the star. Combined with the previous indication of an inhomogeneous density distribution in the gap at $\lesssim$46 au, the notable asymmetry in the outer disk suggests the presence of unseen companions and/or planet-forming processes ongoing in the arc.

preprint2016arXiv

HATS-18 b: An Extreme Short--Period Massive Transiting Planet Spinning Up Its Star

We report the discovery by the HATSouth network of HATS-18 b: a 1.980 +/- 0.077 Mj, 1.337 +0.102 -0.049 Rj planet in a 0.8378 day orbit, around a solar analog star (mass 1.037 +/- 0.047 Msun, and radius 1.020 +0.057 -0.031 Rsun) with V=14.067 +/- 0.040 mag. The high planet mass, combined with its short orbital period, implies strong tidal coupling between the planetary orbit and the star. In fact, given its inferred age, HATS-18 shows evidence of significant tidal spin up, which together with WASP-19 (a very similar system) allows us to constrain the tidal quality factor for Sun-like stars to be in the range 6.5 <= lg(Q*/k_2) <= 7 even after allowing for extremely pessimistic model uncertainties. In addition, the HATS-18 system is among the best systems (and often the best system) for testing a multitude of star--planet interactions, be they gravitational, magnetic or radiative, as well as planet formation and migration theories.

preprint2016arXiv

High-Contrast Imaging of Intermediate-Mass Giants with Long-Term Radial Velocity Trends

A radial velocity (RV) survey for intermediate-mass giants has been operated for over a decade at Okayama Astrophysical Observatory (OAO). The OAO survey has revealed that some giants show long-term linear RV accelerations (RV trends), indicating the presence of outer companions. Direct imaging observations can help clarify what objects generate these RV trends. We present the results of high-contrast imaging observations or six intermediate-mass giants with long-term RV trends using the Subaru Telescope and HiCIAO camera. We detected co-moving companions to $γ$ Hya B ($0.61^{+0.12}_{-0.14} M_\odot$), HD 5608 B ($0.10 \pm 0.01 M_\odot$), and HD 109272 B ($0.28 \pm 0.06 M_\odot$). For the remaining targets($ι$ Dra, 18 Del, and HD 14067) we exclude companions more massive than 30-60 $M_\mathrm{Jup}$ at projected separations of 1arcsec-7arcsec. We examine whether these directly imaged companions or unidentified long-period companions can account for the RV trends observed around the six giants. We find that the Kozai mechanism can explain the high eccentricity of the inner planets $ι$ Dra b, HD 5608 b, and HD 14067 b.

preprint2016arXiv

K2 Discovers a Busy Bee: An Unusual Transiting Neptune Found in the Beehive Cluster

Open clusters have been the focus of several exoplanet surveys but only a few planets have so far been discovered. The \emph{Kepler} spacecraft revealed an abundance of small planets around small, cool stars, therefore, such cluster members are prime targets for exoplanet transit searches. Kepler's new mission, K2, is targeting several open clusters and star-forming regions around the ecliptic to search for transiting planets around their low-mass constituents. Here, we report the discovery of the first transiting planet in the intermediate-age (800 Myr) Beehive cluster (Praesepe). K2-95 is a faint ($\mathrm{Kp = 15.5\,mag}$) $\mathrm{M3.0\pm0.5}$ dwarf from K2's Campaign 5 with an effective temperature of $\mathrm{3471 \pm 124\,K}$, approximately solar metallicity and a radius of $\mathrm{0.402 \pm 0.050 \,R_\odot}$. We detected a transiting planet with a radius of $\mathrm{3.47^{+0.78}_{-0.53} \, R_\oplus}$ and an orbital period of 10.134 days. We combined photometry, medium/high-resolution spectroscopy, adaptive optics/speckle imaging and archival survey images to rule out any false positive detection scenarios, validate the planet, and further characterize the system. The planet's radius is very unusual as M-dwarf field stars rarely have Neptune-sized transiting planets. The comparatively large radius of K2-95b is consistent with the other recently discovered cluster planets K2-25b (Hyades) and K2-33b (Upper Scorpius), indicating systematic differences in their evolutionary states or formation. These discoveries from K2 provide a snapshot of planet formation and evolution in cluster environments and thus make excellent laboratories to test differences between field-star and cluster planet populations.

preprint2016arXiv

Mass estimates for very cold (< 8 K) gas in molecular cloud cores

The mass of prestellar cores is an essential ingredient to understand the onset of star formation in the core. The low level of emission from cold dust may keep parts of it hidden from observation. We aim to determine the fraction of core mass in the temperature range < 8 K that can be expected for typical low- and high-mass star formation regions. We calculate the dust temperature within standard spherically symmetric prestellar cores for a grid of density powerlaws in the outer core regions, core masses, and variation in the external multi-component radiation field. The dust is assumed to be composed of amorphous silicate and carbon, and variations of its optical properties are discussed. As measure for the distribution of cores and clumps, we use core mass functions derived for various environments. In view of the high densities in very cold central regions, dust and gas temperatures are assumed to be equal. We find that the fraction of mass with temperatures < 8 K in typical low- and high-mass cores is < 20\%. It is possible to obtain higher fractions of very cold gas by placing intermediate- or high-mass cores in a typical low-mass star formation environment. We show that the mass uncertainty arising from FIR to mm modeling of very cold dust emission is smaller than the mass uncertainty due to the unknown dust opacities.

preprint2016arXiv

Mass measurements in protoplanetary disks from hydrogen deuteride

The total gas mass of a protoplanetary disk is a fundamental, but poorly determined, quantity. A new technique \citep{bergin+13} has been demonstrated to assess directly the bulk molecular gas reservoir of molecular hydrogen using the HD J=1-0 line at 112 $μ$m. In this work we present a {\it Herschel} Space Observatory\footnote{Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA.} survey of six additional T Tauri disks in the HD line. Line emission is detected at $>$3$σ$ significance in two cases: DM Tau and GM Aur. For the other four disks, we establish upper limits to the line flux. Using detailed disk structure and ray tracing models, we calculate the temperature structure and dust mass from modeling the observed spectral energy distributions, and include the effect of UV gas heating to determine the amount of gas required to fit the HD line. The range of gas masses are 1.0-4.7$\times10^{-2}$ for DM Tau and 2.5-20.4$\times10^{-2}$ for GM Aur. These values are larger than those found using CO for GM Aur, while the CO-derived gas mass for DM Tau is consistent with the lower end of our mass range. This suggests a CO chemical depletion from the gas phase of up to a factor of five for DM Tau and up to two orders of magnitude for GM Aur. We discuss how future analysis can narrow the mass ranges further.

preprint2016arXiv

Measuring Turbulence in TW Hya with ALMA: Methods and Limitations

We obtain high spatial and spectral resolution images of the CO J=2-1, CN N=2-1 and CS J=5-4 emission with ALMA in Cycle~2. The radial distribution of the turbulent broadening is derived with three approaches: two `direct' and one modelling. The first requires a single transition and derives \Tex{} directly from the line profile, yielding a \vturb{}. The second assumes two different molecules are co-spatial thus their relative linewidths allow for a calculation of \Tkin{} and \vturb{}. Finally we fit a parametric disk model where physical properties of the disk are described by power laws, to compare our `direct' methods with previous values. The two direct methods were limited to the outer $r > 40$~au disk due to beam smear. The direct method found \vturb{} ranging from $\approx$~\vel{130} at 40~au, dropping to $\approx$~\vel{50} in the outer disk, qualitatively recovered with the parametric model fitting. This corresponds to roughly $0.2 - 0.4~c_s$. CN was found to exhibit strong non-LTE effects outside $r \approx 140$~au, so \vturb{} was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed linewidths only within $r \lesssim 100$~au, within which \vturb{} was found to drop from \vel{100} ($\approx~0.4~c_s$) to nothing at 100~au. The parametric model yielded a near constant \vel{50} for CS ($0.2 - 0.4~c_s$). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require to reach an absolute calibration precision of order 3\%, or to find a suitable combination of light and heavy molecules which are co-located in the disk.

preprint2016arXiv

Near-Infrared Imaging Polarimetry of LkCa 15: A Possible Warped Inner Disk

We present high-contrast H-band polarized intensity images of the transitional disk around the young solar-like star LkCa 15. By utilizing Subaru/HiCIAO for polarimetric differential imaging, both the angular resolution and the inner working angle reach 0.07" and r=0.1", respectively. We obtained a clearly resolved gap (width <~ 27 AU) at ~ 48 AU from the central star. This gap is consistent with images reported in previous studies. We also confirmed the existence of a bright inner disk with a misaligned position angle of 13+/-4 degree with respect to that of the outer disk, i.e., the inner disk is possibly warped. The large gap and the warped inner disk both point to the existence of a multiple planetary system with a mass of <~1Mjup.

preprint2016arXiv

Prestellar core modeling in the presence of a filament - The dense heart of L1689B

Short version: We apply a new synergetic radiative transfer method: the derived 1D density profiles are both consistent with a cut through the Herschel PACS/SPIRE and JCMT SCUBA-2 continuum maps of L1689B and with a derived local interstellar radiation field. Choosing an appropriate cut along the filament major axis, we minimize the impact of the filament emission on the modeling. For the bulk of the core (5000-20000 au) an isothermal sphere model with a temperature of around 10 K provides the best fits. We show that the power law index of the density profile, as well as the constant temperature can be derived directly from the radial surface brightness profiles. For the inner region (< 5000 au), we find a range of densities and temperatures that are consistent with the surface brightness profiles and the local interstellar radiation field. Based on our core models, we find that pixel-by-pixel single temperature spectral energy distribution fits are incapable of determining dense core properties. We conclude that, to derive physical core properties, it is important to avoid an azimuthal average of core and filament. Correspondingly, derived core masses are too high since they include some mass of the filament, and might introduce errors when determining core mass functions. The forward radiative transfer methods also avoids the loss of information owing to smearing of all maps to the coarsest spatial resolution. We find the central core region to be colder and denser than estimated in recent inverse radiative transfer modeling, possibly indicating the start of star formation in L1689B.

preprint2016arXiv

Rotating Bullets from A Variable Protostar

We present SMA CO(2-1) observations toward the protostellar jet driven by SVS13A, a variable protostar in the NGC1333 star-forming region. The SMA CO(2-1) images show an extremely high-velocity jet composed of a series of molecular 'bullets'. Based on the SMA CO observations, we discover clear and large systematic velocity gradients, perpendicular to the jet axis, in the blueshifted and redshifted bullets. After discussing several alternative interpretations, such as twin-jets, jet precession, warped disk, and internal helical shock, we suggest that the systematic velocity gradients observed in the bullets result from the rotation of the SVS13A jet. From the SMA CO images, the measured rotation velocities are 11.7-13.7 km/s for the blueshifted bullet and 4.7+/-0.5 km/s for the redshifted bullet. The estimated specific angular momenta of the two bullets are comparable to those of dense cores, about 10 times larger than those of protostellar envelopes, and about 20 times larger than those of circumstellar disks. If the velocity gradients are due to the rotation of the SVS13A jet, the significant amount of specific angular momenta of the bullets indicates that the rotation of jets/outflows is a key mechanism to resolve the so-called 'angular momentum problem' in the field of star formation. The kinematics of the bullets suggests that the jet launching footprint on the disk has a radius of about 7.2-7.7 au, which appears to support the extended disk-wind model. We note that further observations are needed to comprehensively understand the kinematics of the SVS13A jet, in order to confirm the rotation nature of the bullets.

preprint2016arXiv

Spiral Density Waves in a Young Protoplanetary Disk

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

preprint2016arXiv

The imprint of exoplanet formation history on observable present-day spectra of hot Jupiters

The composition of a planet's atmosphere is determined by its formation, evolution, and present-day insolation. A planet's spectrum therefore may hold clues on its origins. We present a "chain" of models, linking the formation of a planet to its observable present-day spectrum. The chain links include (1) the planet's formation and migration, (2) its long-term thermodynamic evolution, (3) a variety of disk chemistry models, (4) a non-gray atmospheric model, and (5) a radiometric model to obtain simulated spectroscopic observations with JWST and ARIEL. In our standard chemistry model the inner disk is depleted in refractory carbon as in the Solar System and in white dwarfs polluted by extrasolar planetesimals. Our main findings are: (1) Envelope enrichment by planetesimal impacts during formation dominates the final planetary atmospheric composition of hot Jupiters. We investigate two, under this finding, prototypical formation pathways: a formation inside or outside the water iceline, called "dry" and "wet" planets, respectively. (2) Both the "dry" and "wet" planets are oxygen-rich (C/O<1) due to the oxygen-rich nature of the solid building blocks. The "dry" planet's C/O ratio is <0.2 for standard carbon depletion, while the "wet" planet has typical C/O values between 0.1 and 0.5 depending mainly on the clathrate formation efficiency. Only non-standard disk chemistries without carbon depletion lead to carbon-rich C/O ratios >1 for the "dry" planet. (3) While we consistently find C/O ratios <1, they still vary significantly. To link a formation history to a specific C/O, a better understanding of the disk chemistry is thus needed.

preprint2016arXiv

The Massive Stellar Population of W49: A Spectroscopic Survey

Massive stars form on different scales ranging from large, dispersed OB associations to compact, dense starburst clusters. The complex structure of regions of massive star formation, and the involved short timescales provide a challenge for our understanding of their birth and early evolution. As one of the most massive and luminous star-forming region in our Galaxy, W49 is the ideal place to study the formation of the most massive stars. By classifying the massive young stars deeply embedded into the molecular cloud of W49, we aim to investigate and trace the star formation history of this region. We analyse near-infrared $K$-band spectroscopic observations of W49 from LBT/LUCI combined with $JHK$ images obtained with NTT/SOFI and LBT/LUCI. Based on $JHK$-band photometry and K-band spectroscopy the massive stars are placed in a Hertzsprung Russell diagram. By comparison with evolutionary models, their age and hence the star formation history of W49 can be investigated. Fourteen O type stars as well as two young stellar objects (YSOs) are identified by our spectroscopic survey. Eleven O-stars are main sequence stars with subtypes ranging from O3 to O9.5, with masses ranging from $\sim 20$ $M_{\odot}$ to $\sim 120$ $M_{\odot}$. Three of the O-stars show strong wind features, and are considered to be Of-type supergiants with masses beyond 100 $M_{\odot}$ . The two YSOs show CO emission, indicative for the presence of circumstellar disks in the central region of the massive cluster. The age of the cluster is estimated as $\sim1.5$ Myr, with star formation still ongoing in different parts of the region. The ionising photons from the central massive stars have not yet cleared the molecular cocoon surrounding the cluster. W49 is comparable to extragalactic star-forming regions and provides us with an unique possibility to study a starburst in detail.

preprint2016arXiv

The SPHERE view of the planet-forming disk around HD100546

We image with unprecedented spatial resolution and sensitivity disk features that could be potential signs of planet-disk interaction. Two companion candidates have been claimed in the disk around the young Herbig Ae/Be star HD100546. Thus, this object serves as an excellent target for our investigation of the natal environment of giant planets. We exploit the power of extreme adaptive optics operating in conjunction with the new high-contrast imager SPHERE to image HD100546 in scattered light. We obtain the first polarized light observations of this source in the visible (with resolution as fine as 2 AU) and new H and K band total intensity images that we analyze with the Pynpoint package. The disk shows a complex azimuthal morphology, where multiple scattering of photons most likely plays an important role. High brightness contrasts and arm-like structures are ubiquitous in the disk. A double-wing structure (partly due to ADI processing) resembles a morphology newly observed in inclined disks. Given the cavity size in the visible (11 AU), the CO emission associated to the planet candidate 'c' might arise from within the circumstellar disk. We find an extended emission in the K band at the expected location of 'b'. The surrounding large-scale region is the brightest in scattered light. There is no sign of any disk gap associated to 'b'.

preprint2016arXiv

The VLA view of the HL Tau Disk - Disk Mass, Grain Evolution, and Early Planet Formation

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings.

preprint2016arXiv

Two Small Temperate Planets Transiting Nearby M Dwarfs in K2 Campaigns 0 and 1

The prime Kepler mission revealed that small planets (<4 R_earth) are common, especially around low-mass M dwarfs. K2, the re-purposed Kepler mission, continues this exploration of small planets around small stars. Here we combine K2 photometry with spectroscopy, adaptive optics imaging, and archival survey images to analyze two small planets orbiting the nearby, field age, M dwarfs K2-26 (EPIC 202083828) and K2-9. K2-26 is an M1.0 +/- 0.5 dwarf at 93 +/- 7 pc from K2 Campaign 0. We validate its 14.5665 d period planet and estimate a radius of 2.67^+0.46_-0.42 R_earth. K2-9 is an M2.5 +/- 0.5 dwarf at 110 +/- 12 pc from K2 Campaign 1. K2-9b was first identified by Montet et al. 2015; here we present spectra and adaptive optics imaging of the host star and independently validate and characterize the planet. Our analyses indicate K2-9b is a 2.25^+0.53_-0.96 R_earth planet with a 18.4498 d period. K2-26b exhibits a transit duration that is too long to be consistent with a circular orbit given the measured stellar radius. Thus, the long transits are likely due to the photoeccentric effect and our transit fits hint at an eccentric orbit. Both planets receive low incident flux from their host stars and have estimated equilibrium temperatures <500 K. K2-9b may receive approximately Earth-like insolation. However, its host star exhibits strong GALEX UV emission which could affect any atmosphere it harbors. K2-26b and K2-9b are representatives of a poorly studied class of small planets with cool temperatures that have radii intermediate to Earth and Neptune. Future study of these systems can provide key insight into trends in bulk composition and atmospheric properties at the transition from silicate dominated to volatile rich bodies.

preprint2015arXiv

A Lucky Imaging search for stellar companions to transiting planet host stars

The presence of stellar companions around planet hosting stars influences the architecture of their planetary systems. To find and characterise these companions and determine their orbits is thus an important consideration to understand planet formation and evolution. For transiting systems even unbound field stars are of interest if they are within the photometric aperture of the light curve measurement. Then they contribute a constant flux offset to the transit light curve and bias the derivation of the stellar and planetary parameters if their existence is unknown. Close stellar sources are, however, easily overlooked by common planet surveys due to their limited spatial resolution. We therefore performed high angular resolution imaging of 49 transiting exoplanet hosts to identify unresolved binaries, characterize their spectral type, and determine their separation. The observations were carried out with the Calar Alto 2.2m telescope using the Lucky Imaging camera AstraLux Norte. All targets were imaged in i' and z' passbands. We found new companion candidates to WASP-14 and WASP-58, and we re-observed the stellar companion candidates to CoRoT-2, CoRoT-3, CoRoT-11, HAT-P-7, HAT-P-8, HAT-P-41, KIC 10905746, TrES-2, TrES-4, and WASP-2. We deduce from the stellar density around all sources that two companion candidates out of the targets with the first position measurement (CoRoT-11, HAT-P-41, KIC 10905746, WASP-14 and WASP-58) are probably unbound. In addition, we re-analyse the influence of the sources close to WASP-14 and WASP-58 on the planetary parameters given in the literature and find no significant changes.

preprint2015arXiv

A multi-wavelength interferometric study of the massive young stellar object IRAS 13481-6124

We present new mid-infrared interferometric observations of the massive young stellar object IRAS 13481-6124, using VLTI/MIDI for spectrally-resolved, long-baseline measurements (projected baselines up to $\sim120$ m) and GSO/T-ReCS for aperture-masking interferometry in five narrow-band filters (projected baselines of $\sim1.8-6.4$ m) in the wavelength range of $7.5-13$ $μ$m. We combine these measurements with previously-published interferometric observations in the $K$ and $N$ bands in order to assemble the largest collection of infrared interferometric observations for a massive YSO to date. Using a combination of geometric and radiative-transfer models, we confirm the detection at mid-infrared wavelengths of the disk previously inferred from near-infrared observations. We show that the outflow cavity is also detected at both near- and mid-infrared wavelengths, and in fact dominates the mid-infrared emission in terms of total flux. For the disk, we derive the inner radius ($\sim1.8$ mas or $\sim6.5$ AU at 3.6 kpc), temperature at the inner rim ($\sim1760$ K), inclination ($\sim48$ deg) and position angle ($\sim107$ deg). We determine that the mass of the disk cannot be constrained without high-resolution observations in the (sub-)millimeter regime or observations of the disk kinematics, and could be anywhere from $\sim10^{-3}$ to $20$ M$_\odot$. Finally, we discuss the prospects of interpreting the spectral energy distributions of deeply-embedded massive YSOs, and warn against attempting to infer disk properties from the SED.

preprint2015arXiv

A nearby M star with three transiting super-Earths discovered by K2

Small, cool planets represent the typical end-products of planetary formation. Studying the archi- tectures of these systems, measuring planet masses and radii, and observing these planets' atmospheres during transit directly informs theories of planet assembly, migration, and evolution. Here we report the discovery of three small planets orbiting a bright (Ks = 8.6 mag) M0 dwarf using data collected as part of K2, the new transit survey using the re-purposed Kepler spacecraft. Stellar spectroscopy and K2 photometry indicate that the system hosts three transiting planets with radii 1.5-2.1 R_Earth, straddling the transition region between rocky and increasingly volatile-dominated compositions. With orbital periods of 10-45 days the planets receive just 1.5-10x the flux incident on Earth, making these some of the coolest small planets known orbiting a nearby star; planet d is located near the inner edge of the system's habitable zone. The bright, low-mass star makes this system an excellent laboratory to determine the planets' masses via Doppler spectroscopy and to constrain their atmospheric compositions via transit spectroscopy. This discovery demonstrates the ability of K2 and future space-based transit searches to find many fascinating objects of interest.

preprint2015arXiv

Absence of Significant Cool Disks in Young Stellar Objects Exhibiting Repetitive Optical Outbursts

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

preprint2015arXiv

Can star cluster environment affect dust input from massive AGB stars?

We examine the fraction of massive asymptotic giant branch (AGB) stars remaining bound in their parent star clusters and the effect of irradiation of these stars by intracluster ultraviolet (UV) field. We employ a set of N-body models of dynamical evolution of star clusters rotating in a galactic potential at the solar galactocentric radius. The cluster models are combined with stellar evolution formulae, a library of stellar spectra, and simple models for SiO photodissociation in circumstellar environment (CSE). The initial stellar masses of clusters are varied from $50\rm M_\odot$ to $10^{5}\rm M_\odot$. Results derived for individual clusters are combined using a mass distribution function for young star clusters. We find that about 30% of massive AGB stars initially born in clusters become members of the field population, while the rest evolves in star clusters. They are irradiated by strong intracluster UV radiation resulting in the decrease of the photodissociation radius of SiO molecules, in many stars down to the dust formation zone. In absence of dust shielding, the UV photons penetrate in the CSE deeper than $10R_*$ in 64% and deeper than $2 R_*$ in 42% of all massive AGB stars. If this suppresses following dust formation, the current injection rate of silicate dust from AGB stars in the local Galaxy decreases from $2.2 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ to $1.8 \times 10^{-4}\rm M_\odot\,kpc^{-2}\,Gyr^{-1}$ at most. A lower revised value of 40% for the expected fraction of presolar silicate grains from massive AGB stars is still high to explain the non-detection of these grains in meteorites.

preprint2015arXiv

Chemistry in Disks. IX. Observations and modeling of HCO+ and DCO+ in DM Tau

We present resolved Plateau de Bure Array observations of DM Tau in lines of HCO+ (3-2), (1-0) and DCO+ (3-2). A power-law fitting approach allowed a derivation of column densities of these two molecules. A chemical inner hole of ~50 AU was found in both HCO+ and DCO+ with DCO+ emission extending to only 450 AU. An isotopic ratio of R_D = N(DCO+) / N(HCO+) was found to range from 0.1 at 50 AU and 0.2 at 450 AU. Chemical modeling allowed an exploration of the sensitivity of these molecular abundances to physical parameters out with temperature, finding that X-rays were the domination ionization source in the HCO+ molecular region and that R_D also is sensitive to the CO depletion. The ionization fraction, assuming a steady state system, was found to be x(e-) ~ 10$^{-7}$. Modeling suggests that HCO+ is the dominant charged molecule in the disk but its contribution to ionization fraction is dwarfed by atmoic ions such as C+, S+ and H+.

preprint2015arXiv

Cold condensation of dust in the ISM

The condensation of complex silicates with pyroxene and olivine composition at conditions prevailing in molecular clouds has been experimentally studied. For this purpose, molecular species comprising refractory elements were forced to accrete on cold substrates representing the cold surfaces of surviving dust grains in the interstellar medium. The efficient formation of amorphous and homogeneous magnesium iron silicates at temperatures of about 12 K has been monitored by IR spectroscopy. The gaseous precursors of such condensation processes in the interstellar medium are formed by erosion of dust grains in supernova shock waves. In the laboratory, we have evaporated glassy silicate dust analogs and embedded the released species in neon ice matrices that have been studied spectroscopically to identify the molecular precursors of the condensing solid silicates. A sound coincidence between the 10 micron band of the interstellar silicates and the 10 micron band of the low-temperature siliceous condensates can be noted.

preprint2015arXiv

Detailed structure of the outer disk around HD 169142 with polarized light in H-band

Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at λ= 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.

preprint2015arXiv

Detection of Sharp Symmetric Features in the Circumbinary Disk Around AK Sco

The Search for Planets Orbiting Two Stars (SPOTS) survey aims to study the formation and distribution of planets in binary systems by detecting and characterizing circumbinary planets and their formation environments through direct imaging. With the SPHERE Extreme Adaptive Optics instrument, a good contrast can be achieved even at small (<300 mas) separations from bright stars, which enables studies of planets and disks in a separation range that was previously inaccessible. Here, we report the discovery of resolved scattered light emission from the circumbinary disk around the well-studied young double star AK Sco, at projected separations in the ~13--40 AU range. The sharp morphology of the imaged feature is surprising, given the smooth appearance of the disk in its spectral energy distribution. We show that the observed morphology can be represented either as a highly eccentric ring around AK Sco, or as two separate spiral arms in the disk, wound in opposite directions. The relative merits of these interpretations are discussed, as well as whether these features may have been caused by one or several circumbinary planets interacting with the disk.

preprint2015arXiv

Dissociative photoionization of polycyclic aromatic hydrocarbon molecules carrying an ethynyl group

The life cycle of the population of interstellar polycyclic aromatic hydrocarbon (PAH) molecules depends partly on the photostability of the individual species. We have studied the dissociative photoionization of two ethynyl-substituted PAH species, namely, 9-ethynylphenanthrene and 1-ethynylpyrene. Their adiabatic ionization energy and the appearance energy of fragment ions have been measured with the photoelectron photoion coincidence (PEPICO) spectroscopy technique. The adiabatic ionization energy has been found at 7.84 +/- 0.02 eV for 9-ethynylphenanthrene and at 7.41 +/- 0.02 eV for 1-ethynylpyrene. These values are similar to those determined for the corresponding non-substituted PAH molecules phenanthrene and pyrene. The appearance energy of the fragment ion indicative of the loss of a H atom following photoionization is also similar for either ethynyl-substituted PAH molecule and its non-substituted counterpart. The measurements are used to estimate the critical energy for the loss of a H atom by the PAH cations and the stability of ethynyl-substituted PAH molecules upon photoionization. We conclude that these PAH derivatives are as photostable as the non-substituted species in HI regions. If present in the interstellar medium, they may play an important role in the growth of interstellar PAH molecules.

preprint2015arXiv

Grain Growth in the Circumstellar Disks of the Young Stars CY Tau and DoAr 25

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at λ = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at λ = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, β, is consistent with the observed dust emission in both disks, with low-β in the inner disk and high-β in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

preprint2015arXiv

Model atmospheres of irradiated exoplanets: The influence of stellar parameters, metallicity, and the C/O ratio

Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen.

preprint2015arXiv

Near-Infrared Polarimetry of the GG Tauri A Binary System

A high angular resolution near-infrared polarized-intensity image of the GG Tau A binary system was obtained with the Subaru Telescope. The image shows the circumbinary disk scattering the light from the central binary. The azimuthal profile of the polarized intensity of the circumbinary disk is roughly reproduced by a simple disk model with the Henyey-Greenstein function and the Rayleigh function, indicating small dust grains at the surface of the disk. Combined with a previous observation of the circumbinary disk, our image indicates that the gap structure in the circumbinary disk orbits anti-clockwise, while material in the disk orbit clockwise. We propose a shadow of material located between the central binary and the circumbinary disk. The separations and position angles of the stellar components of the binary in the past 20 years are consistent with the binary orbit with a = 33.4 AU and e = 0.34.

preprint2015arXiv

Near-IR High-Resolution Imaging Polarimetry of the SU Aur Disk: Clues for Tidal Tails?

We present new high-resolution ($\sim$0\farcs09) $H$-band imaging observations of the circumstellar disk around the T Tauri star SU Aur. Our observations with Subaru-HiCIAO have revealed the presence of scattered light as close as 0\farcs15 ($\sim$20 AU) to the star. Within our image, we identify bright emission associated with a disk with a minimum radius of $\sim$90 AU, an inclination of $\sim$35\degr from the plane of the sky, and an approximate P.A. of 15\degr for the major axis. We find a brightness asymmetry between the northern and southern sides of the disk due to a non-axisymmetric disk structure. We also identify a pair of asymmetric tail structures extending east and west from the disk. The western tail extends at least 2\farcs5 (350 AU) from the star, and is probably associated with a reflection nebula previously observed at optical and near-IR wavelengths. The eastern tail extends at least 1\arcsec (140 AU) at the present signal-to-noise. These tails are likely due to an encounter with an unseen brown dwarf, but our results do not exclude the explanation that these tails are outflow cavities or jets.

preprint2015arXiv

Near-IR Polarized Scattered Light Imagery of the DoAr 28 Transitional Disk

We present the first spatially resolved polarized scattered light H-band detection of the DoAr 28 transitional disk. Our two epochs of imagery detect the scattered light disk from our effective inner working angle of 0.10" (13 AU) out to 0.50" (65 AU). This inner working angle is interior to the location of the system's gap inferred by previous studies using SED modeling (15 AU). We detected a candidate point source companion 1.08" northwest of the system; however, our second epoch of imagery strongly suggests that this object is a background star. We constructed a grid of Monte Carlo Radiative Transfer models of the system, and our best fit models utilize a modestly inclined (50 deg), 0.01 Msun disk that has a partially depleted inner gap from the dust sublimation radius out to ~8 AU. Subtracting this best fit, axi-symmetric model from our polarized intensity data reveals evidence for two small asymmetries in the disk, which could be attributable to variety of mechanisms.

preprint2015arXiv

Optical imaging polarimetry of the LkCa 15 protoplanetary disk with SPHERE ZIMPOL

We present the first optical (590--890 nm) imaging polarimetry observations of the pre-transitional protoplanetary disk around the young solar analog LkCa 15, addressing a number of open questions raised by previous studies. We detect the previously unseen far side of the disk gap, confirm the highly eccentric scattered-light gap shape that was postulated from near-infrared imaging, at odds with the symmetric gap inferred from millimeter interferometry. Furthermore, we resolve the inner disk for the first time and trace it out to 30 AU. This new source of scattered light may contribute to the near-infrared interferometric signal attributed to the protoplanet candidate LkCa 15 b, which lies embedded in the outer regions of the inner disk. Finally, we present a new model for the system architecture of LkCa 15 that ties these new findings together. These observations were taken during science verification of SPHERE ZIMPOL and demonstrate this facility's performance for faint guide stars under adverse observing conditions.

preprint2015arXiv

The Kinematic and Chemical Properties of a Potential Core-Forming Clump: Perseus B1-E

We present 13CO and C18O (1-0), (2-1), and (3-2) maps towards the core-forming Perseus B1-E clump using observations from the James Clerk Maxwell Telescope (JCMT), Submillimeter Telescope (SMT) of the Arizona Radio Observatory, and IRAM 30 m telescope. We find that the 13CO and C18O line emission both have very complex velocity structures, indicative of multiple velocity components within the ambient gas. The (1-0) transitions reveal a radial velocity gradient across B1-E of 1 km/s/pc that increases from north-west to south-east, whereas the majority of the Perseus cloud has a radial velocity gradient increasing from south-west to north-east. In contrast, we see no evidence of a velocity gradient associated with the denser Herschel-identified substructures in B1-E. Additionally, the denser substructures have much lower systemic motions than the ambient clump material, which indicates that they are likely decoupled from the large-scale gas. Nevertheless, these substructures themselves have broad line widths (0.4 km/s) similar to that of the C18O gas in the clump, which suggests they inherited their kinematic properties from the larger-scale, moderately dense gas. Finally, we find evidence of C18O depletion only toward one substructure, B1-E2, which is also the only object with narrow (transonic) line widths. We suggest that as prestellar cores form, their chemical and kinematic properties are linked in evolution, such that these objects must first dissipate their turbulence before they deplete in CO.

preprint2015arXiv

The LEECH Exoplanet Imaging Survey: Characterization of the Coldest Directly Imaged Exoplanet, GJ 504 b, and Evidence for Super-Stellar Metallicity

As gas giant planets and brown dwarfs radiate away the residual heat from their formation, they cool through a spectral type transition from L to T, which encompasses the dissipation of cloud opacity and the appearance of strong methane absorption. While there are hundreds of known T-type brown dwarfs, the first generation of directly-imaged exoplanets were all L-type. Recently, Kuzuhara et al. (2013) announced the discovery of GJ 504 b, the first T dwarf exoplanet. GJ 504 b provides a unique opportunity to study the atmosphere of a new type of exoplanet with a ~500 K temperature that bridges the gap between the first directly imaged planets (~1000 K) and our own Solar System's Jupiter (~130 K). We observed GJ 504 b in three narrow L-band filters (3.71, 3.88, and 4.00 microns), spanning the red end of the broad methane fundamental absorption feature (3.3 microns) as part of the LEECH exoplanet imaging survey. By comparing our new photometry and literature photometry to a grid of custom model atmospheres, we were able to fit GJ 504 b's unusual spectral energy distribution for the first time. We find that GJ 504 b is well-fit by models with the following parameters: T_eff=544+/-10 K, g<600 m/s^2, [M/H]=0.60+/-0.12, cloud opacity parameter of f_sed=2-5, R=0.96+/-0.07 R_Jup, and log(L)=-6.13+/-0.03 L_Sun, implying a hot start mass of 3-30 M_jup for a conservative age range of 0.1-6.5 Gyr. Of particular interest, our model fits suggest that GJ 504 b has a super-stellar metallicity. Since planet formation can create objects with non-stellar metallicities, while binary star formation cannot, this result suggests that GJ 504 b formed like a planet, not like a binary companion.

preprint2015arXiv

The LEECH Exoplanet Imaging Survey: Orbit and Component Masses of the Intermediate Age, Late-Type Binary NO UMa

We present high-resolution Large Binocular Telescope LBTI/LMIRcam images of the spectroscopic and astrometric binary NO UMa obtained as part of the LBTI Exozodi Exoplanet Common Hunt (LEECH) exoplanet imaging survey. Our H, K$_s$, and L'-band observations resolve the system at angular separations <0.09". The components exhibit significant orbital motion over a span of ~7 months. We combine our imaging data with archival images, published speckle interferometry measurements, and existing spectroscopic velocity data to solve the full orbital solution and estimate component masses. The masses of the K2.0$\pm$0.5 primary and K6.5$\pm$0.5 secondary are 0.83$\pm$0.02 M$_{\odot}$ and 0.64$\pm$0.02 M$_{\odot}$, respectively. We also derive a system distance of d = 25.87$\pm$0.02 pc and revise the Galactic kinematics of NO UMa. Our revised Galactic kinematics confirm NO UMa as a nuclear member of the ~500 Myr old Ursa Major moving group and it is thus a mass and age benchmark. We compare the masses of the NO UMa binary components to those predicted by five sets of stellar evolution models at the age of the Ursa Major group. We find excellent agreement between our measured masses and model predictions with little systematic scatter between the models. NO UMa joins the short list of nearby, bright, late-type binaries having known ages and fully characterized orbits.

preprint2015arXiv

Variability in a Young, L/T Transition Planetary-Mass Object

As part of our ongoing NTT SoFI survey for variability in young free-floating planets and low mass brown dwarfs, we detect significant variability in the young, free-floating planetary mass object PSO J318.5-22, likely due to rotational modulation of inhomogeneous cloud cover. A member of the 23$\pm$3 Myr $β$ Pic moving group, PSO J318.5-22 has T$_\mathrm{eff}$ = 1160$^{+30}_{-40}$ K and a mass estimate of 8.3$\pm$0.5 M$_{Jup}$ for a 23$\pm$3 Myr age. PSO J318.5-22 is intermediate in mass between 51 Eri b and $β$ Pic b, the two known exoplanet companions in the $β$ Pic moving group. With variability amplitudes from 7-10$\%$ in J$_{S}$ at two separate epochs over 3-5 hour observations, we constrain the rotational period of this object to $>$5 hours. In K$_{S}$, we marginally detect a variability trend of up to 3$\%$ over a 3 hour observation. This is the first detection of weather on an extrasolar planetary mass object. Among L dwarfs surveyed at high-photometric precision ($<$3$\%$) this is the highest amplitude variability detection. Given the low surface gravity of this object, the high amplitude preliminarily suggests that such objects may be more variable than their high mass counterparts, although observations of a larger sample is necessary to confirm this. Measuring similar variability for directly imaged planetary companions is possible with instruments such as SPHERE and GPI and will provide important constraints on formation. Measuring variability at multiple wavelengths can help constrain cloud structure.

preprint2014arXiv

An Enigmatic Pointlike Feature within the HD 169142 Transitional Disk

We report the detection of a faint pointlike feature possibly related to ongoing planet-formation in the disk of the transition disk star HD 169142. The pointlike feature has a $Δ$mag(L)$\sim$6.4, at a separation of $\sim$0.11" and PA$\sim$0$^{\circ}$. Given its lack of an H or K$_{S}$ counterpart despite its relative brightness, this candidate cannot be explained by purely photospheric emission and must be a disk feature heated by an as yet unknown source. Its extremely red colors make it highly unlikely to be a background object, but future multi-wavelength followup is necessary for confirmation and characterization of this feature.

preprint2014arXiv

Broad-band transmission spectrum and K-band thermal emission of WASP-43b as observed from the ground

(Abridged) We observed one transit and one occultation of the hot Jupiter WASP-43b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2-meter telescope. From the transit event, we have independently derived WASP-43's system parameters with high precision, and improved the period to be 0.81347437(13) days. No significant variation in transit depths is detected, with the largest deviations coming from the i', H, and K bands. Given the observational uncertainties, the broad-band transmission spectrum can be explained by either a flat featureless straight line that indicates thick clouds, synthetic spectra with absorption signatures of atomic Na/K or molecular TiO/VO that indicate cloud-free atmosphere, or a Rayleigh scattering profile that indicates high-altitude hazes. From the occultation event, we have detected planetary dayside thermal emission in the K-band with a flux ratio of 0.197 +/- 0.042%, which confirms previous detections obtained in the 2.09 micron narrow band and Ks-band. The K-band brightness temperature 1878 +108/-116 K favors an atmosphere with poor day- to night-side heat redistribution. We also have a marginal detection in the i'-band (0.037 +0.023/-0.021%), which is either a false positive, a signature of non-blackbody radiation at this wavelength, or an indication of reflective hazes at high altitude.

preprint2014arXiv

Characterizing the Youngest Herschel-detected Protostars I. Envelope Structure Revealed by CARMA Dust Continuum Observations

We present CARMA 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample (Stutz et al.). These objects are characterized by very red 24 \micron\ to 70 \micron\ colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from 4 protostars and 1 starless core in the fields toward the PBRS; we also report 1 new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of $\sim$5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv-distance, with more than 50\% of the source emission arising from radii $<$ 1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with $ρ$~$\propto$~R$^{-2.5}$. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power-law. The large amount of mass on scales $<$ 1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

preprint2014arXiv

Fragmentation and Kinematics of dense molecular cores in the filamentary infrared-dark cloud G011.11-0.12

We present new Plateau de Bure Interferometer observations of a region in the filamentary infrared-dark cloud (IRDC) G011.11-0.12 containing young, star-forming cores. In addition to the 3.2mm continuum emission from cold dust, we map this region in the N$_2$H$^+$(1-0) line to trace the core kinematics with an angular resolution of 2" and velocity resolution of 0.2km s$^{-1}$. These data are presented in concert with recent {\em Herschel} results, single-dish N$_2$H$^+$(1-0) data, SABOCA 350$μ$m continuum data, and maps of the C$^{18}$O (2-1) transition obtained with the IRAM 30m telescope. We recover the star-forming cores at 3.2mm continuum, while in N$_2$H$^+$ they appear at the peaks of extended structures. The mean projected spacing between N$_2$H$^+$ emission peaks is 0.18pc, consistent with simple isothermal Jeans fragmentation. The 0.1pc-sized cores have low virial parameters on the criticality borderline, while on the scale of the whole region, we infer that it is undergoing large-scale collapse. The N$_2$H$^+$ linewidth increases with evolutionary stage, while CO isotopologues show no linewidth variation with core evolution. Centroid velocities of all tracers are in excellent agreement, except in the starless region where two N$_2$H$^+$ velocity components are detected, one of which has no counterpart in C$^{18}$O. We suggest that gas along this line of sight may be falling into the quiescent core, giving rise to the second velocity component, possibly connected to the global collapse of the region.

preprint2014arXiv

Giant Molecular Filaments in the Milky Way

Throughout the Milky Way, molecular clouds typically appear filamentary, and mounting evidence indicates that this morphology plays an important role in star formation. What is not known is to what extent the dense filaments most closely associated with star formation are connected to the surrounding diffuse clouds up to arbitrarily large scales. How are these cradles of star formation linked to the Milky Way's spiral structure? Using archival Galactic plane survey data, we have used multiple datasets in search of large-scale, velocity-coherent filaments in the Galactic plane. In this paper, we present our methods employed to identify coherent filamentary structures first in extinction and confirmed using Galactic Ring Survey data. We present a sample of seven Giant Molecular Filaments (GMFs) that have lengths of order $\sim$100 pc, total masses of 10$^4$ - 10$^5$ M$_{\odot}$, and exhibit velocity coherence over their full length. The GMFs we study appear to be inter-arm clouds and may be the Milky Way analogues to spurs observed in nearby spiral galaxies. We find that between 2 and 12% of the total mass (above $\sim$10$^{20}$ cm$^{-2}$) is "dense" (above 10$^{22}$ cm$^{-2}$), where filaments near spiral arms in the Galactic midplane tend to have higher dense gas mass fractions than those further from the arms.

preprint2014arXiv

Ground-based detection of the near-infrared emission from the dayside of WASP-5b

(Abridged) WASP-5b is a highly irradiated dense hot Jupiter orbiting a G4V star every 1.6 days. We observed two secondary eclipses of WASP-5b in the J, H and K bands simultaneously. Thermal emission of WASP-5b is detected in the J and K bands. The retrieved planet-to-star flux ratios in the J and K bands are 0.168 +0.050/-0.052% and 0.269+/-0.062%, corresponding to brightness temperatures of 2996 +212/-261K and 2890 +246/-269K, respectively. No thermal emission is detected in the H band, with a 3-sigma upper limit of 0.166%, corresponding to a maximum temperature of 2779K. On the whole, our J, H, K results can be explained by a roughly isothermal temperature profile of ~2700K in the deep layers of the planetary dayside atmosphere that are probed at these wavelengths. Together with Spitzer observations, which probe higher layers that are found to be at ~1900K, a temperature inversion is ruled out in the range of pressures probed by the combined data set. While an oxygen-rich model is unable to explain all the data, a carbon-rich model provides a reasonable fit but violates energy balance.

preprint2014arXiv

Herbig-Haro objects and mid-infrared outflows in the Vela C molecular cloud

We have performed a deep [SII]6717/6731 wide field Herbig-Haro (HH) object survey toward the Vela C molecular cloud with a sky coverage of about 2 deg2. In total, 18 new HH objects, HH 1090-1107, are discovered and the two previously known HH objects, HH 73-74, are also detected in our [SII] images. We also present an investigation of mid-infrared outflows in the Vela C molecular cloud using the Wide-field Infrared Survey Explorer images taken from AllWISE data release. Using the method suggested by Zhang & Wang, eleven extended green objects (EGOs) are identified to be the mid-infrared outflows, including 6 new mid-infrared outflows that have not been detected previously at other wavelengths and 5 mid-infrared counterparts of the HH objects detected in this work. Using the AllWISE Source Catalog and the source classification scheme suggested by Koenig et al., we have identified 56 young stellar object (YSO) candidates in the Vela C molecular cloud. The possible driving sources of the HH objects and EGOs are discussed based on the morphology of HH objects and EGOs and the locations of HH objects, EGOs and YSO candidates. Finally we associate 12 HH objects and five EGOs with 10 YSOs and YSO candidates. The median length of the outflows in Vela C is 0.35pc and the outflows seem to be oriented randomly.

preprint2014arXiv

Hierarchical fragmentation and differential star formation in the Galactic "Snake": infrared dark cloud G11.11-0.12

We present Submillimeter Array (SMA) $λ=$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

preprint2014arXiv

High contrast imaging at the LBT: the LEECH exoplanet imaging survey

In Spring 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its $\sim$130-night campaign from the Large Binocular Telescope (LBT) atop Mt Graham, Arizona. This survey benefits from the many technological achievements of the LBT, including two 8.4-meter mirrors on a single fixed mount, dual adaptive secondary mirrors for high Strehl performance, and a cold beam combiner to dramatically reduce the telescope's overall background emissivity. LEECH neatly complements other high-contrast planet imaging efforts by observing stars at L' (3.8 $μ$m), as opposed to the shorter wavelength near-infrared bands (1-2.4 $μ$m) of other surveys. This portion of the spectrum offers deep mass sensitivity, especially around nearby adolescent ($\sim$0.1-1 Gyr) stars. LEECH's contrast is competitive with other extreme adaptive optics systems, while providing an alternative survey strategy. Additionally, LEECH is characterizing known exoplanetary systems with observations from 3-5$μ$m in preparation for JWST.

preprint2014arXiv

High-Resolution Submillimeter and Near-Infrared Studies of the Transition Disk around Sz 91

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(3--2) observations with the Submillimeter Array ($\sim1\arcsec$--3$\arcsec$ resolution), and high-resolution imaging of polarized intensity at the $K_s$-band by using the HiCIAO instrument on the Subaru Telescope ($0\farcs25$ resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 AU and 170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H$_2$ mass of $2.4\times10^{-3}$ $M_\sun$ in the cold ($T<$30 K) outer part at $65<r<170$ AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount ($>3\times10^{-9}$ $M_\sun$) of hot ($T\sim$180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3--2) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.

preprint2014arXiv

Indications of M-dwarf Deficits in the Halo and Thick Disk of the Galaxy

We compared the number of faint stars detected in deep survey fields with the current stellar distribution model of the Galaxy and found that the detected number in the H band is significantly smaller than the predicted number. This indicates that M-dwarfs, the major component, are fewer in the halo and the thick disk. We used archived data of several surveys in both the north and south field of GOODS (Great Observatories Origins Deep Survey), MODS in GOODS-N, and ERS and CANDELS in GOODS-S. The number density of M-dwarfs in the halo has to be 20+/-13% relative to that in the solar vicinity, in order for the detected number of stars fainter than 20.5 mag in the H band to match with the predicted value from the model. In the thick disk, the number density of M-dwarfs must be reduced (52+/-13%) or the scale height must be decreased (~600 pc). Alternatively, overall fractions of the halo and thick disks can be significantly reduced to achieve the same effect, because our sample mainly consists of faint M-dwarfs. Our results imply that the M-dwarf population in regions distant from the Galactic plane is significantly smaller than previously thought. We then discussed the implications this has on the suitability of the model predictions for the prediction of non-companion faint stars in direct imaging extrasolar planet surveys by using the best-fit number densities.

preprint2014arXiv

Infrared H3+ and CO Studies of the Galactic Core: GCIRS 3 and GCIRS 1W

We have obtained improved spectra of key fundamental band lines of H3+, R(1,1)l, R(3,3)l, and R(2,2)l, and ro-vibrational transitions of CO on sightlines toward the luminous infrared sources GCIRS 3 and GCIRS 1W, each located in the Central Cluster of the Galactic center within several arcseconds of Sgr A*. The spectra reveal absorption occurring in three kinds of gaseous environments: (1) cold dense and diffuse gas associated with foreground spiral/lateral arms; (2) warm and diffuse gas absorbing over a wide and mostly negative velocity range, which appears to fill a significant fraction of the Galaxy's Central Molecular Zone (CMZ); and (3) warm, dense and compact clouds with velocities near +50 km s^-1 probably within 1-2 pc of the center. The absorptions by the first two cloud types are nearly identical for all the sources in the Central Cluster, and are similar to those previously observed on sightlines from Sgr A* to 30 pc east of it. Cloud type (3), which has only been observed toward the Central Cluster, shows distinct differences between the sightlines to GCIRS 3 and GCIRS 1W, which are separated on the sky by only 0.33 pc in projection. We identify this material as part of an inward extension of the Circumnuclear Disk previously known from HCN mapping. Lower limits on the products of the hydrogen ionization rate zeta and the path length L are 2.3 x 10^5 cm s^-1 and 1.5 x 10^3 cm s^-1 for the warm and diffuse CMZ gas and for the warm and dense clouds in the core, respectively. The limits indicate that the ionization rates in these regions are well above 10^-15 s^-1.

preprint2014arXiv

Kepler-432 b: a massive planet in a highly eccentric orbit transiting a red giant

We report the first disclosure of the planetary nature of Kepler-432 b (aka Kepler object of interest KOI-1299.01). We accurately constrained its mass and eccentricity by high-precision radial velocity measurements obtained with the CAFE spectrograph at the CAHA 2.2-m telescope. By simultaneously fitting these new data and Kepler photometry, we found that Kepler-432 b is a dense transiting exoplanet with a mass of Mp = 4.87 +/- 0.48 MJup and radius of Rp = 1.120 +/- 0.036 RJup. The planet revolves every 52.5 d around a K giant star that ascends the red giant branch, and it moves on a highly eccentric orbit with e = 0.535 +/- 0.030. By analysing two NIR high-resolution images, we found that a star is located at 1.1 from Kepler-432, but it is too faint to cause significant effects on the transit depth. Together with Kepler-56 and Kepler-91, Kepler-432 occupies an almost-desert region of parameter space, which is important for constraining the evolutionary processes of planetary systems.

preprint2014arXiv

LEECH: A 100 Night Exoplanet Imaging Survey at the LBT

In February 2013, the LEECH (LBTI Exozodi Exoplanet Common Hunt) survey began its 100-night campaign from the Large Binocular Telescope atop Mount Graham in Arizona. LEECH nearly complements other high-contrast planet imaging efforts by observing stars in L' band (3.8 microns) as opposed to the shorter wavelength near-infrared bands (1-2.3 microns). This part of the spectrum offers deeper mass sensitivity for intermediate age (several hundred Myr-old) systems, since their Jovian-mass planets radiate predominantly in the mid-infrared. In this proceedings, we present the science goals for LEECH and a preliminary contrast curve from some early data.

preprint2014arXiv

Life Cycle of Dust in the Magellanic Clouds and the Milky Way

To a great extent, our understanding of the life cycle of dust is based on the observational and theoretical studies of the Milky Way and the Magellanic Clouds, which will be the topic of this contribution. Over past years, a large volume of observations with unprecedented spatial resolution has been accumulated for the Milky Way. It permits investigations of different stages of the life cycle of dust, from its formation in stellar sources to destruction in star-forming regions and supernovae shocks. Observations of dust emission, extinction, polarisation of light, and interstellar element depletions in the solar neighbourhood provide the most accurate constraints for the reference dust models applied to study extragalactic systems. However, global spatial studies of the circumstellar and interstellar dust are complicated in the Milky Way disk because of high extinction, confusion along the line of sight and large uncertainties in distances. In contrast, the favourable location in the sky and the proximity of the Magellanic Clouds allow detailed multi-wavelength studies of the dust-forming stellar populations and the investigation of variations of the interstellar grain properties for the entire galaxies. They enable the first comparison between the global stardust production rates from theoretical calculations and those from observations, which confirm discrepancy between accumulated stardust mass and observed interstellar dust mass - "the missing dust-source" problem. Modelling of the life cycle of dust in the Large Magellanic Cloud showed that dust growth by mantle accretion in the ISM, a major dust source in the Milky Way, can be responsible for the existing dust mass in the LMC. We will present comparison of the dust input from different sources to the dust budgets of the Milky Way and the Magellanic Clouds, which reveals how the role of these dust sources depends on metallicity.

preprint2014arXiv

Observed spectral energy distribution of the thermal emission from the dayside of WASP-46b

We aim to construct a spectral energy distribution (SED) for the emission from the dayside atmosphere of the hot Jupiter WASP-46b and to investigate its energy budget. We observed a secondary eclipse of WASP-46b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2m telescope. Eclipse depths of the acquired light curves were derived to infer the brightness temperatures at multibands that cover the SED peak. We report the first detection of the thermal emission from the dayside of WASP-46b in the K band at 4.2-sigma level and tentative detections in the H (2.5-sigma) and J (2.3-sigma) bands, with flux ratios of 0.253 +0.063/-0.060%, 0.194 +/- 0.078%, and 0.129 +/- 0.055%, respectively. The derived brightness temperatures (2306 +177/-187K, 2462 +245/-302K, and 2453 +198/-258K, respectively) are consistent with an isothermal temperature profile of 2386K, which is significantly higher than the dayside-averaged equilibrium temperature, indicative of very poor heat redistribution efficiency. We also investigate the tentative detections in the g'r'i' bands and the 3-sigma upper limit in the z' band, which might indicate the existence of reflective clouds if these tentative detections do not arise from systematics.

preprint2014arXiv

Orbital Monitoring of the AstraLux Large M-dwarf Multiplicity Sample

Orbital monitoring of M-type binaries is essential for constraining their fundamental properties. This is particularly useful in young systems, where the extended pre-main sequence evolution can allow for precise isochronal dating. Here, we present the continued astrometric monitoring of the more than 200 binaries of the AstraLux Large Multiplicity Survey, building both on our previous work, archival data, and new astrometric data spanning the range of 2010-2012. The sample is very young overall -- all included stars have known X-ray emission, and a significant fraction (18%) of them have recently also been identified as members of young moving groups in the Solar neighborhood. We identify ~30 targets that both have indications of being young and for which an orbit either has been closed or appears possible to close in a reasonable timeframe (a few years to a few decades). One of these cases, GJ 4326, is however identified as probably being substantially older than has been implied from its apparent moving group membership, based on astrometric and isochronal arguments. With further astrometric monitoring, these targets will provide a set of empirical isochrones, against which theoretical isochrones can be calibrated, and which can be used to evaluate the precise ages of nearby young moving groups.

preprint2014arXiv

Planetary population synthesis coupled with atmospheric escape: a statistical view of evaporation

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates.

preprint2014arXiv

Probing for Exoplanets Hiding in Dusty Debris Disks: Disk Imaging, Characterization, and Exploration with HST/STIS Multi-Roll Coronagraphy

Spatially resolved scattered-light images of circumstellar (CS) debris in exoplanetary systems constrain the physical properties and orbits of the dust particles in these systems. They also inform on co-orbiting (but unseen) planets, systemic architectures, and forces perturbing starlight-scattering CS material. Using HST/STIS optical coronagraphy, we have completed the observational phase of a program to study the spatial distribution of dust in ten CS debris systems, and one "mature" protoplanetrary disk all with HST pedigree, using PSF-subtracted multi-roll coronagraphy. These observations probe stellocentric distances > 5 AU for the nearest stars, and simultaneously resolve disk substructures well beyond, corresponding to the giant planet and Kuiper belt regions in our Solar System. They also disclose diffuse very low-surface brightness dust at larger stellocentric distances. We present new results inclusive of fainter disks such as HD92945 confirming, and better revealing, the existence of a narrow inner debris ring within a larger diffuse dust disk. Other disks with ring-like sub-structures, significant asymmetries and complex morphologies include: HD181327 with a posited spray of ejecta from a recent massive collision in an exo-Kuiper belt; HD61005 suggested interacting with the local ISM; HD15115 & HD32297, discussed also in the context of environmental interactions. These disks, and HD15745, suggest debris system evolution cannot be treated in isolation. For AU Mic's edge-on disk, out-of-plane surface brightness asymmetries at > 5 AU may implicate one or more planetary perturbers. Time resolved images of the MP Mus proto-planetary disk provide spatially resolved temporal variability in the disk illumination. These and other new images from our program enable direct inter-comparison of the architectures of these exoplanetary debris systems in the context of our own Solar System.

preprint2014arXiv

Recovery of the Candidate Protoplanet HD 100546 b with Gemini/NICI and Detection of Additional (Planet-Induced?) Disk Structure at Small Separations

We report the first independent, second-epoch (re-)detection of a directly-imaged protoplanet candidate. Using $L^\prime$ high-contrast imaging of HD 100546 taken with the Near-Infrared Coronagraph and Imager (NICI) on Gemini South, we recover `HD 100546 b' with a position and brightness consistent with the original VLT/NaCo detection from Quanz et al, although data obtained after 2013 will be required to decisively demonstrate common proper motion. HD 100546 b may be spatially resolved, up to $\approx$ 12-13 AU in diameter, and is embedded in a finger of thermal IR bright, polarized emission extending inwards to at least 0.3". Standard hot-start models imply a mass of $\approx$ 15 $M_{J}$. But if HD 100546 b is newly formed or made visible by a circumplanetary disk, both of which are plausible, its mass is significantly lower (e.g. 1--7 $M_{J}$). Additionally, we discover a thermal IR-bright disk feature, possibly a spiral density wave, at roughly the same angular separation as HD 100546 b but 90 degrees away. Our interpretation of this feature as a spiral arm is not decisive, but modeling analyses using spiral density wave theory implies a wave launching point exterior to $\approx$ 0.45" embedded within the visible disk structure: plausibly evidence for a second, hitherto unseen wide-separation planet. With one confirmed protoplanet candidate and evidence for 1--2 others, HD 100546 is an important evolutionary precursor to intermediate-mass stars with multiple super-jovian planets at moderate/wide separations like HR 8799.

preprint2014arXiv

Resolved Multifrequency Radio Observations of GG Tau

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

preprint2014arXiv

SEEDS Adaptive Optics Imaging of the Asymmetric Transition Disk Oph IRS 48 in Scattered Light

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$μ$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

preprint2014arXiv

SPOTS: The Search for Planets Orbiting Two Stars. I. Survey description and first observations

Direct imaging surveys for exoplanets commonly exclude binary stars from their target lists, leaving a large part of the overall planet demography unexplored. To address this gap in our understanding of planet formation and evolution, we have launched the first direct imaging survey dedicated to circumbinary planets: SPOTS, the Search for Planets Orbiting Two Stars. In this paper, we discuss the theoretical context, scientific merit, and technical feasibility of such observations, describe the target sample and observational strategy of our survey, and report on the first results from our pilot survey of 26 targets with the VLT NaCo facility. While we have not found any confirmed substellar companions to date, a number of promising candidate companions remain to be tested for common proper motion in upcoming follow-up observations. We also report on the astrometry of the three resolved binaries in our target sample. This pilot survey constitutes a successful proof of concept for our survey strategy and paves the way for a second stage of exploratory observations with VLT SPHERE.

preprint2014arXiv

Surface Geometry of Protoplanetary Disks Inferred From Near-Infrared Imaging Polarimetry

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

preprint2014arXiv

The AstraLux Multiplicity Survey: Extension to Late M-dwarfs

The distribution of multiplicity among low-mass stars is a key issue to understanding the formation of stars and brown dwarfs, and recent surveys have yielded large enough samples of nearby low-mass stars to study this issue statistically to good accuracy. Previously, we have presented a multiplicity study of ~700 early/mid M-type stars observed with the AstraLux high-resolution Lucky Imaging cameras. Here, we extend the study of multiplicity in M-type stars through studying 286 nearby mid/late M-type stars, bridging the gap between our previous study and multiplicity studies of brown dwarfs. Most of the targets have been observed more than once, allowing us to assess common proper motion to confirm companionship. We detect 68 confirmed or probable companions in 66 systems, of which 41 were previously undiscovered. Detections are made down to the resolution limit of ~100 mas of the instrument. The raw multiplicity in the AstraLux sensitivity range is 17.9%, leading to a total multiplicity fraction of 21-27% depending on the mass ratio distribution, which is consistent with being flat down to mass ratios of ~0.4, but cannot be stringently constrained below this value. The semi-major axis distribution is well represented by a log-normal function with mu_a = 0.78 and sigma_a = 0.47, which is narrower and peaked at smaller separations than for a Sun-like sample. This is consistent with a steady decrease in average semi-major axis from the highest-mass binary stars to the brown dwarf binaries.

preprint2014arXiv

The Discovery of a Very Massive Star in W49

Very massive stars (M>100 M$_{\odot}$) are very rare objects, but have a strong influence on their environment. The formation of this kind of objects is of prime importance in star formation, but observationally still poorly constrained. We report on the identification of a very massive star in the central cluster of the star-forming region W49. We investigate near-infrared K-band spectroscopic observations of W49 from VLT/ISAAC together with JHK images obtained with NTT/SOFI and LBT/LUCI. We derive a spectral type of W49nr1, the brightest star in the dense core of the central cluster of W49. On the basis of its K-band spectrum, W49nr1 is classified as an O2-3.5If* star with a K-band absolute magnitude of -6.27$\pm$0.10 mag. The effective temperature and bolometric correction are estimated from stars of similar spectral type. After comparison to the Geneva evolutionary models, we find an initial mass between 100 M$_{\odot}$ and 180 M$_{\odot}$. Varying the extinction law results in a larger initial mass range of 90 - 250 M$_{\odot}$.

preprint2014arXiv

The extraordinary far-infrared variation of a protostar: Herschel/PACS observations of LRLL54361

We report Herschel/PACS photometric observations at 70 μm and 160 μm of LRLL54361 - a suspected binary protostar that exhibits periodic (P=25.34 days) flux variations at shorter wavelengths (3.6 μm and 4.5 μm) thought to be due to pulsed accretion caused by binary motion. The PACS observations show unprecedented flux variation at these far-infrared wavelengths that are well cor- related with the variations at shorter wavelengths. At 70 μm the object increases its flux by a factor of six while at 160μm the change is about a factor of two, consistent with the wavelength dependence seen in the far-infrared spectra. The source is marginally resolved at 70 μm with varying FWHM. Deconvolved images of the sources show elongations exactly matching the outflow cavities traced by the scattered light observations. The spatial variations are anti-correlated with the flux variation indicating that a light echo is responsible for the changes in FWHM. The observed far-infrared flux variability indicates that the disk and en- velope of this source is periodically heated by the accretion pulses of the central source, and suggests that such long wavelength variability in general may provide a reasonable proxy for accretion variations in protostars.

preprint2014arXiv

The Herschel/PACS view of the Cep OB2 region: Global protoplanetary disk evolution and clumpy star formation

(Abridged) We use Herschel PACS observations at 70 and 160$μ$m to probe the protoplanetary disks around young stars in the CepOB2 clusters Tr37 and NGC7160 and to trace the small-scale cloud structure. We detect 95 protoplanetary disks at 70$μ$m, 41 at 160$μ$m, and obtain upper limits for over 130 objects. The detection fraction at 70$μ$m depends on the spectral type (88% for K4 or earlier, 17% for M3 or later stars) and on the disk type ($\sim$50% for full and pre-transitional disks, $\sim$35% for transitional disks, no low-excess/depleted disks detected). Non-accreting disks are consistent with significantly lower masses. Accreting transition and pre-transition disks have higher 70$μ$m excesses than full disks, suggestive of more massive, flared and/or thicker disks. Herschel data also reveal several mini-clusters in Tr37, small, compact structures containing a few young stars surrounded by nebulosity. Far-IR data are an excellent probe of the evolution of disks that are too faint for submillimetre observations. We find a strong link between far-IR emission and accretion, and between the inner and outer disk structure. Herschel confirms the dichotomy between accreting and non-accreting transition disks. Substantial mass depletion and global evolution need to occur to shut down accretion in a protoplanetary disk, even if the disk has inner holes. Disks likely follow different evolutionary paths: Low disk masses do not imply opening inner holes, and having inner holes does not require low disk masses. The mini-clusters reveal multi-episodic star formation in Tr37. The long survival of mini-clusters suggest that they formed from the fragmentation of the same core. Their various morphologies favour different formation/triggering mechanisms acting within the same cluster. Herschel also unveils what could be the first heavy mass loss episode of the O6.5 star HD206267.

preprint2014arXiv

The Structure of Pre-transitional Protoplanetary Disks. II. Azimuthal Asymmetries, Different Radial Distributions of Large and Small Dust Grains in PDS~70

The formation scenario of a gapped disk, i.e., transitional disk, and its asymmetry is still under debate. Proposed scenarios such as disk-planet interaction, photoevaporation, grain growth, anticyclonic vortex, eccentricity, and their combinations would result in different radial distributions of the gas and the small (sub-$μ$m size) and large (millimeter size) dust grains as well as asymmetric structures in a disk. Optical/near-infrared (NIR) imaging observations and (sub-)millimeter interferometry can trace small and large dust grains, respectively; therefore multi-wavelength observations could help elucidate the origin of complicated structures of a disk. Here we report SMA observations of the dust continuum at 1.3~mm and $^{12}$CO~$J=2\rightarrow1$ line emission of the pre-transitional protoplanetary disk around the solar-mass star PDS~70. PDS~70, a weak-lined T Tauri star, exhibits a gap in the scattered light from its disk with a radius of $\sim$65~AU at NIR wavelengths. However, we found a larger gap in the disk with a radius of $\sim$80~AU at 1.3~mm. Emission from all three disk components (the gas and the small and large dust grains) in images exhibits a deficit in brightness in the central region of the disk, in particular, the dust-disk in small and large dust grains has asymmetric brightness. The contrast ratio of the flux density in the dust continuum between the peak position to the opposite side of the disk reaches 1.4. We suggest the asymmetries and different gap-radii of the disk around PDS~70 are potentially formed by several (unseen) accreting planets inducing dust filtration.

preprint2014arXiv

The Transiting Exoplanet Survey Satellite

The Transiting Exoplanet Survey Satellite (TESS) will search for planets transiting bright and nearby stars. TESS has been selected by NASA for launch in 2017 as an Astrophysics Explorer mission. The spacecraft will be placed into a highly elliptical 13.7-day orbit around the Earth. During its two-year mission, TESS will employ four wide-field optical CCD cameras to monitor at least 200,000 main-sequence dwarf stars with I = 4-13 for temporary drops in brightness caused by planetary transits. Each star will be observed for an interval ranging from one month to one year, depending mainly on the star's ecliptic latitude. The longest observing intervals will be for stars near the ecliptic poles, which are the optimal locations for follow-up observations with the James Webb Space Telescope. Brightness measurements of preselected target stars will be recorded every 2 min, and full frame images will be recorded every 30 min. TESS stars will be 10-100 times brighter than those surveyed by the pioneering Kepler mission. This will make TESS planets easier to characterize with follow-up observations. TESS is expected to find more than a thousand planets smaller than Neptune, including dozens that are comparable in size to the Earth. Public data releases will occur every four months, inviting immediate community-wide efforts to study the new planets. The TESS legacy will be a catalog of the nearest and brightest stars hosting transiting planets, which will endure as highly favorable targets for detailed investigations.

preprint2014arXiv

The young binary HD 102077: Orbit, spectral type, kinematics, and moving group membership

The K-type binary star HD 102077 was proposed as a candidate member of the TW Hydrae Association (TWA) which is a young (5-15 Myr) moving group in close proximity (~50 pc) to the solar system. The aim of this work is to verify this hypothesis by different means. We first combine diffraction-limited observations from the ESO NTT 3.5m telescope in SDSS-i' and -z' passbands and ESO 3.6m telescope in H-band with literature data to obtain a new, amended orbit fit, estimate the spectral types of both components, and reanalyse the Hipparcos parallax and proper motion taking the orbital motion into account. Moreover, we use two high-resolution spectra of HD 102077 obtained with the fibre-fed optical echelle spectrograph FEROS at the MPG/ESO 2.2m telescope to determine the radial velocity and the lithium equivalent width of the system. The trajectory of HD 102077 is well constrained and we derive a total system mass of $2.6 \pm 0.8\,$ M$_{\odot}$ and a semi-major axis of $14.9 \pm 1.6\,$AU. From the i'-z' colours we infer an integrated spectral type of K2V, and individual spectral types of K0 +/- 1 and K5 +/- 1. The radial velocity corrected for the orbital motion of the system is $17.6 \pm 2\,$km/s. Even though the parallax determination from the Hipparcos data is not influenced by the orbital motion, the proper motion changes to $μ_α*\cos(δ) = -137.84 \pm 1.26\,$ mas/yr and $μ_δ= -33.53 \pm 1.45 \,$mas/yr. With the resultant space motion, the probability of HD 102077 being a member of TWA is less than 1%. Furthermore, the lithium equivalent width of $200 \pm 4\,$mÅ$\,$ is consistent with an age between 30 Myr and 120 Myr and thus older than the predicted age of TWA. In conclusion, HD 102077's age, galactic space motion, and position do not fit TWA or any other young moving group.

preprint2014arXiv

Unfolding the Laws of Star Formation: The Density Distribution of Molecular Clouds

The formation of stars shapes the structure and evolution of entire galaxies. The rate and efficiency of this process are affected substantially by the density structure of the individual molecular clouds in which stars form. The most fundamental measure of this structure is the probability density function of volume densities (rho-PDF), which determines the star formation rates predicted with analytical models. This function has remained unconstrained by observations. We have developed an approach to quantify rho-PDFs and establish their relation to star formation. The rho-PDFs instigate a density threshold of star formation and allow us to quantify the star formation efficiency above it. The rho-PDFs provide new constraints for star formation theories and correctly predict several key properties of the star-forming interstellar medium.

preprint2013arXiv

A Discovery of a Candidate Companion to a Transiting System KOI-94: A Direct Imaging Study for a Possibility of a False Positive

We report a discovery of a companion candidate around one of {\it Kepler} Objects of Interest (KOIs), KOI-94, and results of our quantitative investigation of the possibility that planetary candidates around KOI-94 are false positives. KOI-94 has a planetary system in which four planetary detections have been reported by {\it Kepler}, suggesting that this system is intriguing to study the dynamical evolutions of planets. However, while two of those detections (KOI-94.01 and 03) have been made robust by previous observations, the others (KOI-94.02 and 04) are marginal detections, for which future confirmations with various techniques are required. We have conducted high-contrast direct imaging observations with Subaru/HiCIAO in $H$ band and detected a faint object located at a separation of $\sim0.6''$ from KOI-94. The object has a contrast of $\sim 1\times 10^{-3}$ in $H$ band, and corresponds to an M type star on the assumption that the object is at the same distance of KOI-94. Based on our analysis, KOI-94.02 is likely to be a real planet because of its transit depth, while KOI-94.04 can be a false positive due to the companion candidate. The success in detecting the companion candidate suggests that high-contrast direct imaging observations are important keys to examine false positives of KOIs. On the other hand, our transit light curve reanalyses lead to a better period estimate of KOI-94.04 than that on the KOI catalogue and show that the planetary candidate has the same limb darkening parameter value as the other planetary candidates in the KOI-94 system, suggesting that KOI-94.04 is also a real planet in the system.

preprint2013arXiv

A Herschel and APEX Census of the Reddest Sources in Orion: Searching for the Youngest Protostars

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 um and 160 um that are either too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 um band. We find that the 11 reddest protostar candidates with log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 um fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our analysis reveals that the PBRs sample is composed of Class 0 like sources characterized by very red SEDs (T_bol < 45 K) and large values of sub-millimeter fluxes (L_smm/L_bol > 0.6%). Modified black-body fits to the SEDs provide lower limits to the envelope masses of 0.2 M_sun to 2 M_sun and luminosities of 0.7 L_sun to 10 L_sun. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates.

preprint2013arXiv

A Herschel view of IC 1396 A: Unveiling the different sequences of star formation

The IC1396A globule in the young cluster Tr37, hosting many young stars and protostars, is assumed to be a site of triggered star formation. We mapped IC1396A with Herschel/PACS at 70 and 160 micron. The Herschel maps trace in great detail the very embedded protostellar objects and the structure of the cloud. PACS data reveal a previously unknown Class 0 object (IC1396A-PACS-1) located behind the ionization front. IC1396A-PACS-1 is not detectable with Spitzer, but shows marginal X-ray emission. The data also allowed to study three of the Class I intermediate-mass objects within the cloud. We derived approximate cloud temperatures to study the effect and potential interactions between the protostars and the cloud. The Class 0 object is associated with the densest and colder part of IC1396A. Heating in the cloud is dominated by the winds and radiation of the O6.5 star HD 206267 and, to a lesser extent, by the effects of the Herbig Ae star V 390 Cep. The surroundings of the Class I and Class II objects embedded in the cloud also appear warmer than the sourceless areas, although most of the low-mass objects cannot be individually extracted due to distance and beam dilution. The observations suggest that at least two episodes of star formation have occurred in IC1396A. One would have originated the known, ~1 Myr-old Class I and II objects in the cloud, and a new wave of star formation would have produced the Class 0 source at the tip of the brigth-rimmed cloud. From its location and properties, IC1396A-PACS-1 is consistent with triggering via radiative driven implosion (RDI) induced by HD 206267. The mechanisms behind the formation of the more evolved population of Class I/II/III objects in the cloud are uncertain. Heating of most of the remaining cloud by Class I/Class II objects and by HD 206267 itself may preclude further star formation in the region.

preprint2013arXiv

Accretion of gas onto gap-opening planets and circumplanetary flow structure in magnetized turbulent disks

We have performed three-dimensional magneto-hydrodynamical simulations of stellar accretion disks, using the PLUTO code, and studied the accretion of gas onto a Jupiter-mass planet and the structure of the circumplanetary gas flow after opening a gap in the disk. We compare our results with simulations of laminar, yet viscous disks with different levels of an $α$-type viscosity. In all cases, we find that the accretion flow across the surface of the Hill sphere of the planet is not spherically or azimuthally symmetric, and is predominantly restricted to the mid-plane region of the disk. Even in the turbulent case, we find no significant vertical flow of mass into the Hill sphere. The outer parts of the circumplanetary disk are shown to rotate significantly below Keplerian speed, independent of viscosity, while the circumplanetary disk density (therefore the angular momentum) increases with viscosity. For a simulation of a magnetized turbulent disk, where the global averaged alpha stress is $α_{MHD}=10^{-3}$, we find the accretion rate onto the planet to be $\dot{M}\approx 2\times10^{-6}M_{J}yr^{-1}$ for a gap surface density of $12 g cm^{-2}$. This is about a third of the accretion rate obtained in a laminar viscous simulation with equivalent $α$ parameter.

preprint2013arXiv

An Analysis of the Environments of FU Orionis Objects with Herschel

We present Herschel-HIFI, SPIRE, and PACS 50-670 μm imaging and spectroscopy of six FU Orionis-type objects and candidates (FU Orionis, V1735 Cyg, V1515 Cyg, V1057 Cyg, V1331 Cyg, and HBC 722), ranging in outburst date from 1936-2010, from the "FOOSH" (FU Orionis Objects Surveyed with Herschel) program, as well as ancillary results from Spitzer-IRS and the Caltech Submillimeter Observatory. In their system properties (Lbol, Tbol, line emission), we find that FUors are in a variety of evolutionary states. Additionally, some FUors have features of both Class I and II sources: warm continuum consistent with Class II sources, but rotational line emission typical of Class I, far higher than Class II sources of similar mass/luminosity. Combining several classification techniques, we find an evolutionary sequence consistent with previous mid-IR indicators. We detect [O I] in every source at luminosities consistent with Class 0/I protostars, much greater than in Class II disks. We detect transitions of 13CO (J_up of 5 to 8) around two sources (V1735 Cyg and HBC 722) but attribute them to nearby protostars. Of the remaining sources, three (FU Ori, V1515 Cyg, and V1331 Cyg) exhibit only low-lying CO, but one (V1057 Cyg) shows CO up to J = 23 - 22 and evidence for H2O and OH emission, at strengths typical of protostars rather than T Tauri stars. Rotational temperatures for "cool" CO components range from 20-81 K, for ~ 10^50 total CO molecules. We detect [C I] and [N II] primarily as diffuse emission.

preprint2013arXiv

APEX/SABOCA observations of small-scale structure of infrared-dark clouds I. Early evolutionary stages of star-forming cores

Infrared-dark clouds (IRDCs) harbor the early phases of cluster and high-mass star formation and are comprised of cold (~20 K), dense (n > 10$^4$ cm$^{-3}$) gas. The spectral energy distribution (SED) of IRDCs is dominated by the far-infrared and millimeter wavelength regime, and our initial Herschel study examined IRDCs at the peak of the SED with high angular resolution. Here we present a follow-up study using the SABOCA instrument on APEX which delivers 7.8" angular resolution at 350 micron, matching the resolution we achieved with Herschel/PACS, and allowing us to characterize substructure on ~0.1pc scales. Our sample of 11 nearby IRDCs are a mix of filamentary and clumpy morphologies, and the filamentary clouds show significant hierarchical structure, while the clumpy IRDCs exhibit little hierarchical structure. All IRDCs, regardless of morphology, have about 14% of their total mass in small scale core-like structures which roughly follow a trend of constant volume density over all size scales. Out of the 89 protostellar cores we identified in this sample with Herschel, we recover 40 of the brightest and re-fit their SEDs and find their properties agree fairly well with our previous estimates (<T> ~ 19K). We detect a new population of "cold cores" which have no 70 micron counterpart, but are 100 and 160 micron-bright, with colder temperatures (<T> ~ 16K). This latter population, along with SABOCA-only detections, are predominantly low-mass objects, but their evolutionary diagnostics are consistent with the earliest starless or prestellar phase of cores in IRDCs.

preprint2013arXiv

Chemistry in Protoplanetary Disks

This comprehensive review summarizes our current understanding of the evolution of gas, solids and molecular ices in protoplanetary disks. Key findings related to disk physics and chemistry, both observationally and theoretically, are highlighted. We discuss which molecular probes are used to derive gas temperature, density, ionization state, kinematics, deuterium fractionation, and study organic matter in protoplanetary disks.

preprint2013arXiv

Connection between dense gas mass fraction, turbulence driving, and star formation efficiency of molecular clouds

We examine the physical parameters that affect the accumulation of gas in molecular clouds to high column densities where the formation of stars takes place. In particular, we analyze the dense gas mass fraction (DGMF) in a set of self-gravitating, isothermal, magnetohydrodynamic turbulence simulations including sink particles to model star formation. We find that the simulations predict close to exponential DGMFs over the column density range N(H2) = 3-25 x 10^{21} cm^{-2} that can be easily probed via, e.g., dust extinction measurements. The exponential slopes correlate with the type of turbulence driving and also with the star formation efficiency. They are almost uncorrelated with the sonic Mach number and magnetic-field strength. The slopes at early stages of cloud evolution are steeper than at the later stages. A comparison of these predictions with observations shows that only simulations with relatively non-compressive driving (b ~< 0.4) agree with the DGMFs of nearby molecular clouds. Massive infrared dark clouds can show DGMFs that are in agreement with more compressive driving. The DGMFs of molecular clouds can be significantly affected by how compressive the turbulence is on average. Variations in the level of compression can cause scatter to the DGMF slopes, and some variation is indeed necessary to explain the spread of the observed DGMF slopes. The observed DGMF slopes can also be affected by the clouds' star formation activities and statistical cloud-to-cloud variations.

preprint2013arXiv

DIGIT survey of far-infrared lines from protoplanetary discs: II. CO

CO is an important component of a protoplanetary disc as it is one of the most abundant gas phase species. Furthermore, observations of the CO transitions can be used as a diagnostic of the gas, tracing conditions in both the inner and outer disc. We present Herschel/PACS spectroscopy of a sample of 22 Herbig Ae/Be (HAEBEs) and 8 T Tauri stars (TTS), covering the pure rotational CO transitions from J = 14-13 up to J = 49-48. CO is detected in only 5 HAEBEs, AB Aur, HD 36112, HD 97048, HD 100546 and IRS 48, and in 4 TTS, AS 205, S CrA, RU Lup, and DG Tau. The highest transition seen is J = 36-35, with Eup of 3669 K, detected in HD 100546 and DG Tau. We construct rotational diagrams for the discs with at least 3 CO detections to derive Trot, and find average temperatures of 270 K for the HAEBEs and 485 K for the TTS. HD 100546 requires an extra temperature component with Trot ~ 900-1000 K, suggesting a range of temperatures in its disc atmosphere, consistent with thermo-chemical disc models. In HAEBEs, the objects with CO detections all have flared discs in which the gas and dust are thermally decoupled. We use a small model grid to analyse our observations and find that an increased amount of flaring means higher line flux, as it increases the mass in warm gas. CO is not detected in our flat discs as the emission is below the detection limit. We find that HAEBE sources with CO detections have high Luv and strong PAH emission, again connected to the heating of the gas. In TTS, the objects with CO detections are all sources with evidence for a disc wind or outflow. For both groups of objects, sources with CO detections generally have high UV luminosity (either stellar in HAEBEs or due to accretion in TTS), but this is not a sufficient condition for the detection of the far-IR CO lines.

preprint2013arXiv

Direct Imaging Search for Extrasolar Planets in the Pleiades

We carried out an imaging survey for extrasolar planets around stars in the Pleiades (125 Myr, 135 pc) in the $H$ and $K_{S}$ bands using HiCIAO combined with the adaptive optics, AO188, on the Subaru telescope. We found 13 companion candidates fainter than 14.5 mag in the $H$ band around 9 stars. Five of these 13 were confirmed to be background stars by measurement of their proper motion. One was not found in the second epoch observation, and thus was not a background or companion object. One had multi-epoch image, but the precision of its proper motion was not sufficient to conclude whether it was background object. Four other candidates are waiting for second epoch observations to determine their proper motion. Finally, the remaining 2 were confirmed to be 60 $M_{J}$ brown dwarf companions orbiting around HD 23514 (G0) and HII 1348 (K5) respectively, as had been reported in previous studies. In our observations, the average detection limit for a point source was 20.3 mag in the $H$ band beyond 1''.5 from the central star. On the basis of this detection limit, we calculated the detection efficiency to be 90% for a planet with 6 to 12 Jovian masses and a semi-major axis of 50--1000 AU. For this we extrapolated the distribution of planet mass and semi-major axis derived from RV observations and adopted the planet evolution model of Baraffe et al. (2003). As there was no detection of a planet, we estimated the frequency of such planets to be less than 17.9% ($2σ$) around one star of the Pleiades cluster.

preprint2013arXiv

Dust input from AGB stars in the Large Magellanic Cloud

The dust-forming population of AGB stars and their input to the interstellar dust budget of the Large Magellanic Cloud (LMC) are studied with evolutionary dust models with the main goals (1) to investigate how the amount and composition of dust from AGB stars vary over galactic history; (2) to characterise the mass and metallicity distribution of the present population of AGB stars; (3) to quantify the contribution of AGB stars of different mass and metallicity to the present stardust population in the interstellar medium (ISM). We use models of the stardust lifecycle in the ISM developed and tested for the Solar neighbourhood. The first global spatially resolved reconstruction of the star formation history of the LMC from the Magellanic Clouds Photometric Survey is employed to calculate the stellar populations in the LMC. The dust input from AGB stars is dominated by carbon grains from stars with masses < 4 Msun almost over the entire history of the LMC. The production of silicate, silicon carbide and iron dust is delayed until the ISM is enriched to about half the present metallicity in the LMC. For the first time, theoretically calculated dust production rates of AGB stars are compared to those derived from IR observations of AGB stars for the entire galaxy. We find good agreement within scatter of various observational estimates. We show that the majority of silicate and iron grains in the present stardust population originate from a small population of intermediate-mass stars consisting of only about 4% of the total number of stars, whereas in the Solar neighbourhood they originate from low-mass stars. With models of the lifecycle of stardust grains in the ISM we confirm a large discrepancy between dust input from stars and the existing interstellar dust mass in the LMC reported in Matsuura et al. 2009.

preprint2013arXiv

G048.66-0.29: Physical State of an Isolated Site of Massive Star Formation

We present continuum observations of the infrared dark cloud (IRDC) G48.66-0.22 (G48) obtained with Herschel, Spitzer, and APEX, in addition to several molecular line observations. The Herschel maps are used to derive temperature and column density maps of G48 using a model based on a modified blackbody. We find that G48 has a relatively simple structure and is relatively isolated; thus this IRDC provides an excellent target to study the collapse and fragmentation of a filamentary structure in the absence of complicating factors such as strong external feedback. The derived temperature structure of G48 is clearly non-isothermal from cloud to core scale. The column density peaks are spatially coincident with the lowest temperatures (~ 17.5 K) in G48. A total cloud mass of ~390Msun is derived from the column density maps. By comparing the luminosity-to-mass ratio of 13 point sources detected in the Herschel/PACS bands to evolutionary models, we find that two cores are likely to evolve into high-mass stars (M>8 Msun). The derived mean projected separation of point sources is smaller than in other IRDCs but in good agreement with theoretical predications for cylindrical collapse. We detect several molecular species such as CO, HCO+, HCN, HNC and N2H+. CO is depleted by a factor of ~3.5 compared to the expected interstellar abundance, from which we conclude that CO freezes out in the central region. Furthermore, the molecular clumps, associated with the sub-millimeter peaks in G48, appear to be gravitationally unbound or just pressure confined. The analysis of critical line masses in G48 show that the entire filament is collapsing, overcoming any internal support.

preprint2013arXiv

Ground-based detections of thermal emission from the dense hot Jupiter WASP-43b in H and Ks-bands

We report new detections of thermal emission from the transiting hot Jupiter WASP-43b in the H and Ks-bands as observed at secondary eclipses. The observations were made with the WIRCam instrument on the CFHT. We obtained a secondary eclipse depth of 0.103$_{-0.017}^{+0.017}%$ and 0.194$_{-0.029}^{+0.029}%$ in the H and Ks-bands, respectively. The Ks band depth is consistent with previous measurement in the narrow band centered at 2.09um by Gillon et al. (2012). Our eclipse depths in both bands are consistent with a blackbody spectrum with a temperature of ~1850 K, slightly higher than the dayside equilibrium temperature without day-night energy redistribution. Based on theoretical models of the dayside atmosphere of WASP-43b, our data constrain the day-night energy redistribution in the planet to be $\lesssim 15-25$%, depending on the metal content in the atmosphere. Combined with energy balance arguments our data suggest that a strong temperature inversion is unlikely in the dayside atmosphere of WASP-43b. However, a weak inversion cannot be strictly ruled out at the current time. Future observations are required to place detailed constraints on the chemical composition of the atmosphere.

preprint2013arXiv

High-Contrast Near-Infrared Imaging Polarimetry of the Protoplanetary Disk around RY Tau

We present near-infrared coronagraphic imaging polarimetry of RY Tau. The scattered light in the circumstellar environment was imaged at H-band at a high resolution (~0".05) for the first time, using Subaru-HiCIAO. The observed polarized intensity (PI) distribution shows a butterfly-like distribution of bright emission with an angular scale similar to the disk observed at millimeter wavelengths. This distribution is offset toward the blueshifted jet, indicating the presence of a geometrically thick disk or a remnant envelope, and therefore the earliest stage of the Class II evolutionary phase. We perform comparisons between the observed PI distribution and disk models with: (1) full radiative transfer code, using the spectral energy distribution (SED) to constrain the disk parameters; and (2) monochromatic simulations of scattered light which explore a wide range of parameters space to constrain the disk and dust parameters. We show that these models cannot consistently explain the observed PI distribution, SED, and the viewing angle inferred by millimeter interferometry. We suggest that the scattered light in the near-infrared is associated with an optically thin and geometrically thick layer above the disk surface, with the surface responsible for the infrared SED. Half of the scattered light and thermal radiation in this layer illuminates the disk surface, and this process may significantly affect the thermal structure of the disk.

preprint2013arXiv

Magnetized Accretion and Dead Zones in Protostellar Disks

The edges of magnetically-dead zones in protostellar disks have been proposed as locations where density bumps may arise, trapping planetesimals and helping form planets. Magneto-rotational turbulence in magnetically-active zones provides both accretion of gas on the star and transport of mass to the dead zone. We investigate the location of the magnetically-active regions in a protostellar disk around a solar-type star, varying the disk temperature, surface density profile, and dust-to-gas ratio. We also consider stellar masses between 0.4 and 2 $M_\odot$, with corresponding adjustments in the disk mass and temperature. The dead zone's size and shape are found using the Elsasser number criterion with conductivities including the contributions from ions, electrons, and charged fractal dust aggregates. The charged species' abundances are found using the approach proposed by S. Okuzumi. The dead zone is in most cases defined by the ambipolar diffusion. In our maps, the dead zone takes a variety of shapes, including a fish-tail pointing away from the star and islands located on and off the midplane. The corresponding accretion rates vary with radius, indicating locations where the surface density will increase over time, and others where it will decrease. We show that density bumps do not readily grow near the dead zone's outer edge, independently of the disk parameters and the dust properties. Instead, the accretion rate peaks at the radius where the gas-phase metals freeze out. This could lead to clearing a valley in the surface density, and to a trap for pebbles located just outside the metal freeze-out line.

preprint2013arXiv

Mapping H-band Scattered Light Emission in the Mysterious SR21 Transitional Disk

We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21. These images were obtained with the Subaru HiCIAO camera, adaptive optics and the polarized differential imaging (PDI) technique. We resolve the disk in scattered light at H-band for stellocentric 0.1"<r<0.6" (12<r<75AU). We compare our results with previously published spatially-resolved 880 micron continuum Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor invoked to explain SR21's sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more moderate depletions (delta=0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (delta~10^-6) are inconsistent with our H-band images when they are assumed to carry over to small grains, suggesting that surface grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. In fact, the radial polarized intensity profile of our H-band observations is smooth and steeply inwardly-increasing (r^-3), with no evidence of a break at the 36AU sub-mm cavity wall. We hypothesize that this profile is dominated by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, and find that we can neither exclude nor verify this scenario. This study demonstrates the power of multiwavelength imaging of transitional disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for clearing these disks of their large midplane grains.

preprint2013arXiv

Polyynyl-substituted PAH molecules and DIB carriers

Polycyclic aromatic hydrocarbon (PAH) molecules have been long considered promising candidates for the carriers of the diffuse interstellar bands (DIBs). The PAH-DIB hypothesis, however, raises two major issues. First, the number of interstellar PAH species is potentially orders of magnitude larger than the number of DIBs. Second, the absorption spectrum of a PAH is in general dominated by bands found at UV wavelengths while, conversely, DIBs are absent from the UV wavelength domain and arise at visible and near IR wavelengths. These issues do not necessarily weaken the PAH-DIB hypothesis and can actually allow us to refine it. In that context, we analyze the UV/vis absorption spectra of PAH molecules isolated in Ne matrices and propose that polyynyl-substituted PAHs, or similar species, are valid candidates for the carriers of the DIBs. Finally, a possible lifecycle for DIB-carrying PAHs is presented.

preprint2013arXiv

Sculpting the disk around T Cha: an interferometric view

(Abridged) Circumstellar disks are believed to be the birthplace of planets and are expected to dissipate on a timescale of a few Myr. The processes responsible for the removal of the dust and gas will strongly modify the radial distribution of the dust and consequently the SED. In particular, a young planet will open a gap, resulting in an inner disk dominating the near-IR emission and an outer disk emitting mostly in the far-IR. We analyze a full set of data (including VLTI/Pionier, VLTI/Midi, and VLT/NaCo/Sam) to constrain the structure of the transition disk around TCha. We used the Mcfost radiative transfer code to simultaneously model the SED and the interferometric observations. We find that the dust responsible for the emission in excess in the near-IR must have a narrow temperature distribution with a maximum close to the silicate sublimation temperature. This translates into a narrow inner dusty disk (0.07-0.11 AU). We find that the outer disk starts at about 12 AU and is partially resolved by the Pionier, Sam, and Midi instruments. We show that the Sam closure phases, interpreted as the signature of a candidate companion, may actually trace the asymmetry generated by forward scattering by dust grains in the upper layers of the outer disk. These observations help constrain the inclination and position angle of the outer disk. The presence of matter inside the gap is difficult to assess with present-day observations. Our model suggests the outer disk contaminates the interferometric signature of any potential companion that could be responsible for the gap opening, and such a companion still has to be unambiguously detected. We stress the difficulty to observe point sources in bright massive disks, and the consequent need to account for disk asymmetries (e.g. anisotropic scattering) in model-dependent search for companions.

preprint2013arXiv

SMA Observations of Class 0 Protostars: A High-Angular Resolution Survey of Protostellar Binary Systems

We present high angular resolution 1.3 mm and 850 um dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.5 arcsec, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64+/-0.08 and 0.91+/-0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I YSOs, and approximately three (for MF) and four (for CSF) times larger than the values found among MS stars, with a similar range of separations. Furthermore, the observed fraction of high order multiple systems to binary systems in Class 0 protostars (0.50+/-0.09) is also larger than the fractions found in Class I YSOs (0.31+/-0.07) and MS stars (< 0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary systems shows a general trend in which companion star fraction increases with decreasing companion separation. We find that 67%+/-8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation.

preprint2013arXiv

The abundances of hydrocarbon functional groups in the interstellar medium inferred from laboratory spectra of hydrogenated and methylated polycyclic aromatic hydrocarbons

Infrared (IR) absorption spectra of individual polycyclic aromatic hydrocarbons (PAHs) containing methyl (-CH3), methylene (>CH2), or diamond-like *CH groups and IR spectra of mixtures of methylated and hydrogenated PAHs prepared by gas phase condensation were measured at room temperature (as grains in pellets) and at low temperature (isolated in Ne matrices). In addition, the PAH blends were subjected to an in-depth molecular structure analysis by means of high-performance liquid chromatography, nuclear magnetic resonance spectroscopy, and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Supported by calculations at the density functional theory level, the laboratory results were applied to analyze in detail the aliphatic absorption complex of the diffuse interstellar medium at 3.4 mu-m and to determine the abundances of hydrocarbon functional groups. Assuming that the PAHs are mainly locked in grains, aliphatic CHx groups (x = 1,2,3) would contribute approximately in equal quantities to the 3.4 mu-m feature (N_{CHx} / N_{H} approx 10^{-5} - 2 * 10^{-5}). The abundances, however, may be two to four times lower if a major contribution to the 3.4 mu-m feature comes from molecules in the gas phase. Aromatic =CH groups seem to be almost absent from some lines of sight, but can be nearly as abundant as each of the aliphatic components in other directions (N_{=CH} / N_{H} < 2 * 10^{-5}; upper value for grains). Due to comparatively low binding energies, astronomical IR emission sources do not display such heavy excess hydrogenation. At best, especially in proto-planetary nebulae, >CH2 groups bound to aromatic molecules, i.e., excess hydrogens on the molecular periphery only, can survive the presence of a nearby star.

preprint2013arXiv

The Link between Magnetic Fields and Filamentary Clouds: Bimodal Cloud Orientations in the Gould Belt

The orientations of filamentary molecular clouds in the Gould Belt and their local ICM (inter-cloud media) magnetic fields are studied using near-infrared dust extinction maps and optical stellar polarimetry data. These filamentary clouds are a few-to-ten parsecs in length, and we find that their orientations tend to be either parallel or perpendicular to the mean field directions of the local ICM. This bimodal distribution is not found in cloud simulations with super-Alfvenic turbulence, in which the cloud orientations should be random. ICM magnetic fields that are dynamically important compared to inertial-range turbulence and self-gravity can readily explain both field-filament configurations. Previous studies commonly recognize that strong magnetic fields can guide gravitational contraction and result in filaments perpendicular to them, but few discuss the fact that magnetic fields can also channel sub-Alfvenic turbulence to form filaments aligned with them. This strong-field scenario of cloud formation is also consistent with the constant field strength observed from ICM to clouds (Crutcher et al. 2010) and is possible to explain the "hub-filament" cloud structure (Myers 2009) and the density threshold of cloud gravitational contraction (Kainulainen et al. 2009).

preprint2013arXiv

The low-mass stellar population in the young cluster Tr37: Disk evolution, accretion, and environment

We present a study of accretion and protoplanetary disks around M-type stars in the 4 Myr-old cluster Tr37. With a well-studied solar-type population, Tr37 is a benchmark for disk evolution. We used low-resolution spectroscopy to identify 141 members (78 new) and 64 probable members, mostly M-type stars. Hα emission provides information about accretion. Optical, 2MASS, Spitzer, and WISE data are used to trace the SEDs. We construct radiative transfer models to explore the structures of full-disks, pre-transition, transition, and dust-depleted disks. Including the new and previously known members, we confirm that a substantial fraction (~2/5) of disks show signs of evolution, either as radial dust evolution (transition/pre-transition disks) or as a more global evolution (low small-dust masses, dust settling, and weak/absent accretion signatures). Accretion is strongly dependent on the SED type. About half of the transition objects are consistent with no accretion, and dust-depleted disks have weak (or undetectable) accretion signatures, especially among M-type stars. The analysis of accretion and disk structure suggests a parallel evolution of dust and gas. We find several distinct classes of evolved disks, based on SED type and accretion, pointing to different disk dispersal mechanisms and probably different evolutionary paths. Dust depletion and opening of inner holes appear to be independent processes: most transition disks are not dust-depleted, and most dust-depleted disks do not require inner holes. The differences in disk structure between M-type and solar-type stars in Tr37 (4 Myr) are not as remarkable as in the young, sparse, Coronet cluster (1-2 Myr), suggesting that other factors, like the environment/interactions, are likely to play a role in the disk evolution and dispersal. Finally, we also find some evidence of clumpy star formation or mini-clusters within Tr37.

preprint2013arXiv

The SEEDS Direct Imaging Survey for Planets and Scattered Dust Emission in Debris Disk Systems

Debris disks around young main-sequence stars often have gaps and cavities which for a long time have been interpreted as possibly being caused by planets. In recent years, several giant planet discoveries have been made in systems hosting disks of precisely this nature, further implying that interactions with planets could be a common cause of such disk structures. As part of the SEEDS high-contrast imaging survey, we are surveying a population of debris disk-hosting stars with gaps and cavities implied by their spectral energy distributions, in order to attempt to spatially resolve the disk as well as to detect any planets that may be responsible for the disk structure. Here we report on intermediate results from this survey. Five debris disks have been spatially resolved, and a number of faint point sources have been discovered, most of which have been tested for common proper motion, which in each case has excluded physical companionship with the target stars. From the detection limits of the 50 targets that have been observed, we find that beta Pic b-like planets (~10 Mjup planets around G--A-type stars) near the gap edges are less frequent than 15--30%, implying that if giant planets are the dominant cause of these wide (27 AU on average) gaps, they are generally less massive than beta Pic b.

preprint2013arXiv

The twofold debris disk around HD 113766 A - Warm and cold dust as seen with VLTI/Midi and Herschel/Pacs

Warm debris disks are a sub-sample of the large population of debris disks, and display excess emission in the mid-IR. Around solar-type stars, very few objects show emission features in mid-IR spectroscopic observations, that are attributed to small, warm silicate dust grains. The origin of this warm dust can possibly be explained either by a collision between several bodies or by transport from an outer belt. We present and analyse new far-IR Herschel/Pacs observations, supplemented by ground-based data in the mid-IR (VLTI/Midi and VLT/Visir), for one of these rare systems: the 10-16 Myr old debris disk around HD 113766 A. We improve an existing model to account for these new observations, and better constrain the spatial distribution of the dust and its composition. We underline the limitations of SED modelling and the need for spatially resolved observations. We find that the system is best described by an inner disk located within the first AU, well constrained by the Midi data, and an outer disk located between 9-13 AU. In the inner dust belt, our previous finding of Fe-rich crystalline olivine grains still holds. We do not observe time variability of the emission features over at least a 8 years time span, in a environment subjected to strong radiation pressure. The time stability of the emission features indicates that μm-sized dust grains are constantly replenished from the same reservoir, with a possible depletion of sub-μm-sized grains. We suggest that the emission features may arise from multi-composition aggregates. We discuss possible scenarios concerning the origin of the warm dust. The compactness of the innermost regions as probed by Midi, as well as the dust composition, suggest that we are witnessing the outcomes of (at least) one collision between partially differentiated bodies, in an environment possibly rendered unstable by terrestrial planetary formation.

preprint2013arXiv

The VLTI/MIDI survey of massive young stellar objects - Sounding the inner regions around intermediate- and high-mass young stars using mid-infrared interferometry

We aim to characterize the distribution and composition of circumstellar material around young massive stars, and to investigate exactly which physical structures in these objects are probed by long-baseline mid-infrared interferometric observations. We used the two-telescope interferometric instrument MIDI of the Very Large Telescope Interferometer of the European Southern Observatory to observe a sample of 24 intermediate- and high-mass young stellar objects in the N band (8-13 micron). We had successful fringe detections for 20 objects, and present spectrally-resolved correlated fluxes and visibility levels for projected baselines of up to 128 m. We fit the visibilities with geometric models to derive the sizes of the emitting regions, as well as the orientation and elongation of the circumstellar material. Fourteen objects in the sample show the 10 micron silicate feature in absorption in the total and correlated flux spectra. For 13 of these objects, we were able to fit the correlated flux spectra with a simple absorption model, allowing us to constrain the composition and absorptive properties of the circumstellar material. Nearly all of the massive young stellar objects observed show significant deviations from spherical symmetry at mid-infrared wavelengths. In general, the mid-infrared emission can trace both disks and outflows, and in many cases it may be difficult to disentangle these components on the basis of interferometric data alone, because of the sparse spatial frequency coverage normally provided by current long-baseline interferometers. For the majority of the objects in this sample, the absorption occurs on spatial scales larger than those probed by MIDI. Finally, the physical extent of the mid-infrared emission around these sources is correlated with the total luminosity, albeit with significant scatter.

preprint2013arXiv

Weather on the Nearest Brown Dwarfs: Resolved Simultaneous Multi-Wavelength Variability Monitoring of WISE J104915.57-531906.1AB

We present two epochs of MPG/ESO 2.2m GROND simultaneous 6-band ($r'i'z'JHK$) photometric monitoring of the closest known L/T transition brown dwarf binary WISE J104915.57-531906.1AB. We report here the first resolved variability monitoring of both the T0.5 and L7.5 components. We obtained 4 hours of focused observations on the night of UT 2013-04-22, as well as 4 hours of defocused (unresolved) observations on the night of UT 2013-04-16. We note a number of robust trends in our light curves. The $r'$ and $i'$ light curves appear to be anticorrelated with $z'$ and $H$ for the T0.5 component and in the unresolved lightcurve. In the defocused dataset, $J$ appears correlated with $z'$ and $H$ and anticorrelated with $r'$ and $i'$, while in the focused dataset we measure no variability for $J$ at the level of our photometric precision, likely due to evolving weather phenomena. In our focused T0.5 component lightcurve, the $K$ band lightcurve displays a significant phase offset relative to both $H$ and $z'$. We argue that the measured phase offsets are correlated with atmospheric pressure probed at each band, as estimated from 1D atmospheric models. We also report low-amplitude variability in $i'$ and $z'$ intrinsic to the L7.5 component.

preprint2013arXiv

Young Stellar Objects in Lynds 1641: Disks, Accretion, and Star Formation History

We investigate the young stellar objects (YSOs) in the Lynds~1641 (L1641) cloud using multi-wavelength data including Spitzer, WISE, 2MASS, and XMM covering ~1390 YSOs across a range of evolutionary stages. In addition, we targeted a sub-sample of YSOs for optical spectroscopy. We use this data, along with archival photometric data, to derive spectral types, extinction values, masses, ages, as well as accretion rates. We obtain a disk fraction of ~50% in L1641. The disk frequency is almost constant as a function of stellar mass with a slight peak at log(M_*/M_sun)\approx-0.25. The analysis of multi-epoch spectroscopic data indicates that the accretion variability of YSOs cannot explain the two orders of magnitude of scatter for YSOs with similar masses. Forty-six new transition disk (TD) objects are confirmed in this work, and we find that the fraction of accreting TDs is lower than for optically thick disks (40--45% vs. 77--79% respectively). We confirm our previous result that the accreting TDs have a similar median accretion rate to normal optically thick disks. We confirm that two star formation modes (isolated vs. clustered) exist in L1641. We find that the diskless YSOs are statistically older than the YSOs with optically-thick disks and the transition disk objects have a median age which is intermediate between the two populations. We tentatively study the star formation history in L1641 based on the age distribution and find that star formation started to be active 2--3 Myr ago.

preprint2012arXiv

A Common Proper Motion Stellar Companion to HAT-P-7

We report that HAT-P-7 has a common proper motion stellar companion. The companion is located at $\sim3.9$ arcsec to the east and estimated as an M5.5V dwarf based on its colors. We also confirm the presence of the third companion, which was first reported by Winn et al. (2009), based on long-term radial velocity measurements. We revisit the migration mechanism of HAT-P-7b given the presence of those companions, and propose sequential Kozai migration as a likely scenario in this system. This scenario may explain the reason for an outlier in the discussion of the spin-orbit alignment timescale for HAT-P-7b by Albrecht et al. (2012).

preprint2012arXiv

A Herschel Survey of Cold Dust in Disks Around Brown Dwarfs and Low-Mass Stars

We report the complete photometric results from our Herschel study which is the first comprehensive program to search for far-infrared emission from cold dust around young brown dwarfs. We surveyed 50 fields containing 51 known or suspected brown dwarfs and very low mass stars that have evidence of circumstellar disks based on Spitzer photometry and/or spectroscopy. The objects with known spectral types range from M3 to M9.5. Four of the candidates were subsequently identified as extragalactic objects. Of the remaining 47 we have successfully detected 36 at 70micron and 14 at 160micron with S/N greater than 3, as well as several additional possible detections with low S/N. The objects exhibit a range of [24]--[70] micron colors suggesting a range in mass and/or structure of the outer disk. We present modeling of the spectral energy distributions of the sample and discuss trends visible in the data. Using two Monte Carlo radiative transfer codes we investigate disk masses and geometry. We find a very wide range in modeled total disk masses from less than 1e-6 solar masses up to 1e-3 solar masses with a median disk mass of order 3e-5 solar masses, suggesting that the median ratio of disk mass to central object mass may be lower than for T Tauri stars. The disk scale heights and flaring angles, however, cover a range consistent with those seen around T Tauri stars. The host clouds in which the young brown dwarfs and low-mass stars are located span a range in estimated age from ~1-3 Myr to ~10 Myr and represent a variety of star-forming environments. No obvious dependence on cloud location or age is seen in the disk properties, though the statistical significance of this conclusion is not strong.

preprint2012arXiv

A Likely Close-In Low-Mass Stellar Companion to the Transitional Disk Star HD 142527

With the uniquely high contrast within 0.1" (Δmag(L') = 5-6.5 magnitudes) available using Sparse Aperture Masking (SAM) with NACO at VLT, we detected asymmetry in the flux from the Herbig Fe star HD 142527 with a barycenter emission situated at a projected separation of 88+/-5 mas (12.8+/-1.5 AU at 145 pc) and flux ratios in H, K, and L' of 0.016+/-0.007, 0.012+/-0.008, 0.0086+/-0.0011 respectively (3-σerrors) relative to the primary star and disk. After extensive closure-phase modeling, we interpret this detection as a close-in, low-mass stellar companion with an estimated mass of ~0.1-0.4 M_Sun. HD 142527 has a complex disk structure, with an inner gap imaged in both the near and mid-IR as well as a spiral feature in the outer disk in the near-IR. This newly detected low-mass stellar companion may provide a critical explanation of the observed disk structure.

preprint2012arXiv

A solution to the radiation pressure problem in the formation of massive stars

We review our recent studies demonstrating that the radiation pressure problem in the formation of massive stars can be circumvented via an anisotropy of the thermal radiation field. Such an anisotropy naturally establishes with the formation of a circumstellar disk. The required angular momentum transport within the disk can be provided by developing gravitational torques. Radiative Rayleigh-Taylor instabilities in the cavity regions - as previously suggested in the literature - are not required and are shown to be not occurring in the context of massive star formation.

preprint2012arXiv

Constraints on the Radial Variation of Grain Growth in the AS 209 Circumstellar Disk

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity κ_ν. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (β), we find that β(R) increases from β<0.5 at \sim20 AU to β>1.5 for R>80 AU, inconsistent with a constant value of β across the disk (at the 10σ level). Furthermore, if radial variations of κ_ν are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

preprint2012arXiv

Experimental and theoretical study on the infrared spectroscopy of astrophysically relevant PAH derivatives 2- and 9-vinylanthracene

We propose to evaluate the contribution of polycyclic aromatic hydrocarbon molecules that carry side groups to the mid-infrared emission spectra. Within this framework, the IR absorption spectra of 2- and 9-vinylanthracene were measured in Ar matrices at 12 K and in CsI and polyethylene pellets at room temperature. The laboratory spectra were analyzed with the support of simulations based on the density functional theory. For each PAH molecule, eight IR spectra were computed by combining the B3LYP functional with as many different basis sets, namely 4-31G, 4-31G(d), 6-31G, 6-311G, 6-31G(d), 6-31G(d,p), 6-31+G(d,p), and 6-31++G(d,p). The comparison of the theoretical spectra with the laboratory data allowed us to determine the most suitable combinations for modeling the IR spectra of neutral PAH molecules that carry a vinyl side group. It was concluded from the examples of 2- and 9-VA that the optimum basis set is 6-31G unless a steric interaction has to be taken into account, in which case the optimum basis set is 6-31G(d). Thus, in the presence of such an interaction, the use of d-type polarization functions is recommended. We discuss the possibility for neutral vinyl-substituted PAHs to contribute to the mid-infrared emission spectra and find that their specific features do not match with the mid-infrared aromatic emission bands.

preprint2012arXiv

Instantaneous starburst of the massive clusters Westerlund 1 and NGC 3603 YC

We present a new method to determine the age spread of resolved stellar populations in a starburst cluster. The method relies on a two-step process. In the first step, kinematic members of the cluster are identified based on multi-epoch astrometric monitoring. In the second step, a Bayesian analysis is carried out, comparing the observed photometric sequence of cluster members with sets of theoretical isochrones. When applying this methodology to optical and near-infrared high angular resolution Hubble Space Telescope (HST) and adaptive optics observations of the ~5 Myr old starburst cluster Westerlund 1 and ~2 Myr old starburst cluster NGC 3603 YC, we derive upper limits for the age spreads of 0.4 and 0.1 Myr, respectively. The results strongly suggest that star formation in these starburst clusters happened almost instantaneously.

preprint2012arXiv

Multiwavelength Observations of V2775 Ori, an Outbursting Protostar in L 1641: Exploring the Edge of the FU Orionis Regime

Individual outbursting young stars are important laboratories for studying the physics of episodic accretion and the extent to which this phenomenon can explain the luminosity distribution of protostars. We present new and archival data for V2775 Ori (HOPS 223), a protostar in the L 1641 region of the Orion molecular clouds that was discovered by Caratti o Garatti et al. (2011) to have recently undergone an order-of-magnitude increase in luminosity. Our near-infrared spectra of the source have strong blueshifted He I 10830 absorption, strong H2O and CO absorption, and no H I emission, all typical of FU Orionis sources. With data from IRTF, 2MASS, HST, Spitzer, WISE, Herschel, and APEX that span from 1 to 70 microns pre-outburst and from 1 to 870 microns post-outburst, we estimate that the outburst began between 2005 April and 2007 March. We also model the pre- and post-outburst spectral energy distributions of the source, finding it to be in the late stages of accreting its envelope with a disk-to-star accretion rate that increased from about 2x10^-6 M_sun/yr to about 10^-5 M_sun/yr during the outburst. The post-outburst luminosity at the epoch of the FU Orionis-like near-IR spectra is 28 L_sun, making V2775 Ori the least luminous documented FU Orionis outburster with a protostellar envelope. The existence of low-luminosity outbursts supports the notion that a range of episiodic accretion phenomena can partially explain the observed spread in protostellar luminosities.

preprint2012arXiv

On the massive young stellar object AFGL4176: High-spatial-resolution multi-wavelength observations and modeling

Deeply embedded and at distances of several kiloparsecs, massive young stellar objects (MYSOs) present numerous challenges for observation and study. In this work, we present spatially-resolved observations of one MYSO, AFGL 4176, together with survey and literature data, ranging from interferometric observations with VLTI/MIDI in the mid-infrared, to single-dish Herschel measurements in the far-infrared, and sub-millimeter data from APEX. We consider this spatially-resolved, multi-wavelength data set in terms of both radiative transfer and geometric models. We find that the observations are well described by one-dimensional models overall, but there are also substantial deviations from spherical symmetry at scales of tens to hundreds of astronomical units, which are revealed by the mid-infrared interferometric measurements. We use a multiple-component, geometric modeling approach to explain the mid-infrared emission on scales of tens to hundreds of astronomical units, and find the MIDI measurements are well described by a model consisting of a one-dimensional Gaussian halo and an inclined (θ=60 deg) circumstellar disk extending out to several hundred astronomical units along a position angle of 160 deg. Finally, we compare our results both with previous models of this source, and with those of other MYSOs, and discuss the present situation with mid-infrared interferometric observations of massive stars.

preprint2012arXiv

Optical spectroscopy of EX Lupi during quiescence and outburst: Infall, wind, and dynamics in the accretion flow

We explore the accretion mechanisms in EX Lupi, prototype of EXor variables, during its quiescence and outburst phases. We analyse high-resolution optical spectra taken before, during, and after its 2008 outburst. In quiescence and outburst, the star presents many permitted emission lines, including typical CTTS lines and numerous neutral and ionized metallic lines. During the outburst, the number of emission lines increases to over a thousand, with narrow plus broad component structure (NC+BC). The BC profile is highly variable on short timescales (24-72h). An active chromosphere can explain the metallic lines in quiescence and the outburst NC. The dynamics of the BC line profiles suggest an origin in a hot, dense, non-axisymmetric, and non-uniform accretion column that suffers velocity variations along the line-of-sight on timescales of days. Assuming Keplerian rotation, the emitting region would be located at ~0.1-0.2 AU, consistent with the inner disk rim, but the velocity profiles of the lines reveal a combination of rotation and infall. Line ratios of ions and neutrals can be reproduced with a temperature of T~6500 K for electron densities of a few times 10$^{12}$cm$^{-3}$ in the line-emitting region. The data confirm that the 2008 outburst was an episode of increased accretion, albeit much stronger than previous EX Lupi and typical EXors outbursts. The line profiles are consistent with the infall/rotation of a non-axisymmetric structure that could be produced by clumpy accretion during the outburst phase. A strong inner disk wind appears in the epochs of higher accretion. The rapid recovery of the system after the outburst and the similarity between the pre-outburst and post-outburst states suggest that the accretion channels are similar during the whole period, and only the accretion rate varies, providing a superb environment for studying the accretion processes.

preprint2012arXiv

Planetesimal formation by sweep-up: How the bouncing barrier can be beneficial to growth

The formation of planetesimals is often accredited to collisional sticking of dust grains. The exact process is unknown, as collisions between larger aggregates tend to lead to fragmentation or bouncing rather than sticking. Recent laboratory experiments have however made great progress in the understanding and mapping of the complex physics involved in dust collisions. We want to study the possibility of planetesimal formation using the results from the latest laboratory experiments, particularly by including the fragmentation with mass transfer effect, which might lead to growth even at high impact velocities. We present a new experimentally and physically motivated dust collision model capable of predicting the outcome of a collision between two particles of arbitrary masses and velocities. It is used together with a continuum dust-size evolution code that is both fast in terms of execution time and able to resolve the dust well at all sizes, allowing for all types of interactions to be studied without biases. We find that for the general dust population, bouncing collisions prevent the growth above millimeter-sizes. However, if a small number of cm-sized particles are introduced, for example due to vertical mixing or radial drift, they can act as a catalyst and start to sweep up the smaller particles. At a distance of 3 AU, 100-meter-sized bodies are formed on a timescale of 1 Myr. We conclude that direct growth of planetesimals might be a possibility thanks to a combination of the existence of a bouncing barrier and the fragmentation with mass transfer effect. The bouncing barrier is here even beneficial, as it prevents the growth of too many large particles that would otherwise only fragment among each other, and creates a reservoir of small particles that can be swept up by larger bodies. However, for this process to work, a few seeds of cm in size or larger have to be introduced.

preprint2012arXiv

Polarimetric Imaging of Large Cavity Structures in the Pre-transitional Protoplanetary Disk around PDS 70: Observations of the disk

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is ~70 AU. Our data show that the geometric center of the disk shifts by ~6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at ~30 to ~50MJ within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap.

preprint2012arXiv

Protostars, multiplicity, and disk evolution in the Corona Australis region: A Herschel Gould Belt Study

The CrA region and the Coronet cluster form a nearby (138 pc), young (1-2 Myr) star-forming region hosting a moderate population of YSO. We present Herschel PACS photometry at 100 and 160 micron, obtained as part of the Herschel Gould Belt Survey. The Herschel maps reveal the cluster members with high sensitivity and high dynamic range. Many of the cluster members are detected, including some embedded, very low-mass objects, several protostars, and substantial emission from the surrounding cloud. The Herschel data provide sufficient spatial resolution to detect small-scale details, such as bright filaments around the IRS5 protostar complex and a bubble-shaped rim associated with the Class I object IRS2. The disks around the Class II objects display a wide range of mid- and far-IR excesses consistent with different disk structures. We have modeled the disks using the RADMC radiative transfer code, finding an interesting mixture of objects for a young and presumably coeval cluster. Some of them are consistent with flared, massive, relatively primordial disks (SCrA, TCrA). Others display significant evidence for inside-out evolution, consistent with the presence of inner holes/gaps (G-85, G-87). Finally, we find disks with a dramatic small dust depletion (G-1, HBC677) that, in some cases, could be related to truncation or to the presence of large gaps in a flared disk (CrA-159). The derived masses for the disks around the low-mass stars are found to be below the typical values in Taurus, in agreement with previous Spitzer observations. Given the high degree of multiplicity and interactions observed among the protostars in the region, the diversity of disks may be a consequence of the early star formation history, which should also be taken into account when studying the disk properties in similar sparsely populated clusters.

preprint2012arXiv

SMA and Spitzer Observations of Bok Glouble CB17: A Candidate First Hydrostatic Core?

We present high angular resolution SMA and Spitzer observations toward the Bok globule CB17. SMA 1.3mm dust continuum images reveal within CB17 two sources with an angular separation of about 21" (about 5250 AU at a distance of 250 pc). The northwestern continuum source, referred to as CB17 IRS, dominates the infrared emission in the Spitzer images, drives a bipolar outflow extending in the northwest-southeast direction, and is classified as a low luminosity Class0/I transition object (L_bol ~ 0.5 L_sun). The southeastern continuum source, referred to as CB17 MMS, has faint dust continuum emission in the SMA 1.3mm observations (about 6 sigma detection; ~3.8 mJy), but is not detected in the deep Spitzer infrared images at wavelengths from 3.6 to 70 micron. Its bolometric luminosity and temperature, estimated from its spectral energy distribution, are less than 0.04 L_sun and 16 K, respectively. The SMA CO(2-1) observations suggest that CB17 MMS may drive a low-velocity molecular outflow (about 2.5 km/s), extending in the east-west direction. Comparisons with prestellar cores and Class0 protostars suggest that CB17 MMS is more evolved than prestellar cores but less evolved than Class0 protostars. The observed characteristics of CB17 MMS are consistent with the theoretical predictions from radiative/magneto hydrodynamical simulations of a first hydrostatic core, but there is also the possibility that CB17 MMS is an extremely low luminosity protostar deeply embedded in an edge-on circumstellar disk. Further observations are needed to study the properties of CB17 MMS and to address more precisely its evolutionary stage.

preprint2012arXiv

Stellar companions to exoplanet host stars: Lucky Imaging of transiting planet hosts

Observed properties of stars and planets in binary/multiple star systems provide clues to planet formation and evolution. We extended our survey for visual stellar companions to the hosts of transiting exoplanets by 21 stars, using the Lucky Imaging technique with the two AstraLux instruments: AstraLux Norte at the Calar Alto 2.2-m telescope, and AstraLux Sur at the ESO 3.5-m New Technology Telescope at La Silla. We present observations of two previously unknown binary candidate companions, to the transiting planet host stars HAT-P-8 and WASP-12, and derive photometric and astrometric properties of the companion candidates. The common proper motions of the previously discovered candidate companions with the exoplanet host stars TrES-4 and WASP-2 are confirmed from follow-up observations. A Bayesian statistical analysis of 31 transiting exoplanet host stars observed with AstraLux suggests that the companion star fraction of planet hosts is not significantly different from that of solar-type field stars, but that the binary separation is on average larger for planet host stars.

preprint2012arXiv

Synthetic observations of first hydrostatic cores in collapsing low-mass dense cores II. Simulated ALMA dust emission maps

First hydrostatic cores are predicted by theories of star formation, but their existence has never been demonstrated convincingly by (sub)millimeter observations. Furthermore, the multiplicity at the early phases of the star formation process is poorly constrained. The purpose of this paper is twofold. First, we seek to provide predictions of ALMA dust continuum emission maps from early Class 0 objects. Second, we show to what extent ALMA will be able to probe the fragmentation scale in these objects. Following our previous paper (Commerçon et al. 2012, hereafter paper I), we post-process three state-of-the-art radiation-magneto-hydrodynamic 3D adaptive mesh refinement calculations to compute the emanating dust emission maps. We then produce synthetic ALMA observations of the dust thermal continuum from first hydrostatic cores. We present the first synthetic ALMA observations of dust continuum emission from first hydrostatic cores. We analyze the results given by the different bands and configurations and we discuss for which combinations of the two the first hydrostatic cores would most likely be observed. We also show that observing dust continuum emission with ALMA will help in identifying the physical processes occurring within collapsing dense cores. If the magnetic field is playing a role, the emission pattern will show evidence of a pseudo-disk and even of a magnetically driven outflow, which pure hydrodynamical calculations cannot reproduce. The capabilities of ALMA will enable us to make significant progress towards understanding fragmentation at the early Class 0 stage and discovering first hydrostatic cores.

preprint2012arXiv

Synthetic observations of first hydrostatic cores in collapsing low-mass dense cores. I. Spectral energy distributions and evolutionary sequence

The low-mass star formation evolutionary sequence is relatively well-defined both from observations and theoretical considerations. The first hydrostatic core is the first protostellar equilibrium object that is formed during the star formation process. Using state-of-the-art radiation-magneto-hydrodynamic 3D adaptive mesh refinement calculations, we aim to provide predictions for the dust continuum emission from first hydrostatic cores. We investigate the collapse and the fragmentation of magnetized one solar mass prestellar dense cores and the formation and evolution of first hydrostatic cores using the RAMSES code. We use three different magnetization levels for the initial conditions, which cover a large variety of early evolutionary morphology, e.g., the formation of a disk or a pseudo-disk, outflow launching, and fragmentation. We post-process the dynamical calculations using the 3D radiative transfer code RADMC-3D. We compute spectral energy distributions and usual evolutionary stage indicators such as bolometric luminosity and temperature. We find that the first hydrostatic core lifetimes depend strongly on the initial magnetization level of the parent dense core. We derive, for the first time, spectral energy distribution evolutionary sequences from high-resolution radiation-magneto-hydrodynamic calculations. We show that under certain conditions, first hydrostatic cores can be identified from dust continuum emission at 24 microns and 70 microns. We also show that single spectral energy distributions cannot help to distinguish between the formation scenarios of the first hydrostatic core, i.e., between the magnetized and non-magnetized models. Spectral energy distributions are a first useful and direct way to target first hydrostatic core candidates but high-resolution interferometry is definitively needed to determine the evolutionary stage of the observed sources.

preprint2012arXiv

The AstraLux Large M-dwarf Multiplicity Survey

We present the results of an extensive high-resolution imaging survey of M-dwarf multiplicity using the Lucky Imaging technique. The survey made use of the AstraLux Norte camera at the Calar Alto 2.2m telescope and the AstraLux Sur camera at the ESO New Technology Telescope in order to cover nearly the full sky. In total, 761 stars were observed (701 M-type and 60 late K-type), among which 182 new and 37 previously known companions were detected in 205 systems. Most of the targets have been observed during two or more epochs, and could be confirmed as physical companions through common proper motion, often with orbital motion being confirmed in addition. After accounting for various bias effects, we find a total M-dwarf multiplicity fraction of 27+/-3% within the AstraLux detection range of 0.08-6" (semi-major axes of ~3-227 AU at a median distance of 30 pc). We examine various statistical multiplicity properties within the sample, such as the trend of multiplicity fraction with stellar mass and the semi-major axis distribution. The results indicate that M-dwarfs are largely consistent with constituting an intermediate step in a continuous distribution from higher-mass stars down to brown dwarfs. Along with other observational results in the literature, this provides further indications that stars and brown dwarfs may share a common formation mechanism, rather than being distinct populations.

preprint2012arXiv

The Earliest Phases of Star Formation (EPoS): A Herschel Key Program - The precursors to high-mass stars and clusters

(Abridged) We present an overview of the sample of high-mass star and cluster forming regions observed as part of the Earliest Phases of Star Formation (EPoS) Herschel Guaranteed Time Key Program. A sample of 45 infrared-dark clouds (IRDCs) were mapped at PACS 70, 100, and 160 micron and SPIRE 250, 350, and 500 micron. In this paper, we characterize a population of cores which appear in the PACS bands and place them into context with their host cloud and investigate their evolutionary stage. We construct spectral energy distributions (SEDs) of 496 cores which appear in all PACS bands, 34% of which lack counterparts at 24 micron. From single-temperature modified blackbody fits of the SEDs, we derive the temperature, luminosity, and mass of each core. These properties predominantly reflect the conditions in the cold, outer regions. Taking into account optical depth effects and performing simple radiative transfer models, we explore the origin of emission at PACS wavelengths. The core population has a median temperature of 20K and has masses and luminosities that span four to five orders of magnitude. Cores with a counterpart at 24 micron are warmer and bluer on average than cores without a 24 micron counterpart. We conclude that cores bright at 24 micron are on average more advanced in their evolution, where a central protostar(s) have heated the outer bulk of the core, than 24 micron-dark cores. The 24 micron emission itself can arise in instances where our line of sight aligns with an exposed part of the warm inner core. About 10% of the total cloud mass is found in a given cloud's core population. We uncover over 300 further candidate cores which are dark until 100 micron. These are candidate starless objects, and further observations will help us determine the nature of these very cold cores.

preprint2012arXiv

The Herschel DIGIT Survey of Weak-line T Tauri Stars: implications for disk evolution and dissipation

As part of the "Dust, Ice, and Gas In Time (DIGIT)" Herschel Open Time Key Program, we present Herschel photometry (at 70, 160, 250, 350 and 500 micron) of 31 Weak-Line T Tauri star (WTTS) candidates in order to investigate the evolutionary status of their circumstellar disks. Thirteen stars in our sample had circumstellar disks previously known from infrared observations at shorter wavelengths, while eighteen of them had no previous evidence for a disk. We detect a total of 15 disks as all previously known disks are detected at one or more Herschel wavelengths and two additional disks are identified for the first time. The spectral energy distributions (SEDs) of our targets seem to trace the dissipation of the primordial disk and the transition to the debris disk regime. Seven of the 15 disks appear to be optically thick primordial disks, including two objects with SEDs indistinguishable from those of typical Classical T Tauri stars, four objects that have significant deficit of excess emission at all IR wavelengths, and one "pre-transitional" object with a known gap in the disk. Despite their previous WTTS classification, we find that the seven targets in our sample with optically thick disks show evidence for accretion. The remaining eight disks have weaker IR excesses similar to those of optically thin debris disks. Six of them are warm and show significant 24 micron Spitzer excesses, while the last two are newly identified cold debris-like disks with photospheric 24 micron fluxes, but significant excess emission at longer wavelengths. The Herschel photometry also places strong constraints on the non-detections, where systems with F70/F70,star > 5 - 15 and L,disk/L,star > 1xE-3 to 1xE-4 can be ruled out. We present preliminary models for both the optically thick and optically thin disks and discuss our results in the context of the evolution and dissipation of circumstellar disks.

preprint2012arXiv

The high-mass disk candidates NGC7538IRS1 and NGC7538S

Context: The nature of embedded accretion disks around forming high-mass stars is one of the missing puzzle pieces for a general understanding of the formation of the most massive and luminous stars. Methods: Using the Plateau de Bure Interferometer at 1.36mm wavelengths in its most extended configuration we probe the dust and gas emission at ~0.3",corresponding to linear resolution elements of ~800AU. Results: NGC7538IRS1 remains a single compact and massive gas core with extraordinarily high column densities, corresponding to visual extinctions on the order of 10^5mag, and average densities within the central 2000AU of ~2.1x10^9cm^-3 that have not been measured before. We identify a velocity gradient across in northeast-southwest direction that is consistent with the mid-infrared emission, but we do not find a gradient that corresponds to the proposed CH3OH maser disk. The spectral line data toward NGC7538IRS1 reveal strong blue- and red-shifted absorption toward the mm continuum peak position. The red-shifted absorption allows us to estimate high infall rates on the order of 10^-2 Msun/yr. Although we cannot prove that the gas will be accreted in the end, the data are consistent with ongoing star formation activity in a scaled-up low-mass star formation scenario. Compared to that, NGC7538S fragments in a hierarchical fashion into several sub-sources. While the kinematics of the main mm peak are dominated by the accompanying jet, we find rotational signatures from a secondary peak. Furthermore, strong spectral line differences exist between the sub-sources which is indicative of different evolutionary stages within the same large-scale gas clump.

preprint2012arXiv

The smoothness of the interstellar extinction curve in the UV. Comparison with recent laboratory measurements of PAH mixtures

Context. As revealed by high-resolution spectral investigations in the wavelength range between 300 and 400 nm, the interstellar extinction curve does not display any of the sharp electronic absorption bands that are characteristic for large polyatomic molecules, such as polycyclic aromatic hydrocarbons (PAHs), which belong to the most abundant interstellar molecules. Aims. We aim to verify whether the absorption curves of mixtures of medium-sized PAHs produced in the laboratory can explain the astronomical observations. Methods. The PAH mixtures were synthesized by infrared laser pyrolysis and subsequent chemical extraction and size separation. The matrix isolation technique was used to study the absorption spectra of isolated molecules at low temperature. Results. Our experimental results demonstrate that the UV-visible absorption curves of PAH mixtures can be very smooth, displaying no sharp bands, if the molecular diversity is sufficiently high. Conclusions. In view of the absence of sharp electronic features on the interstellar extinction curve for 300 < λ < 400 nm, we conclude from our experimental findings that the interstellar PAH population must be very diverse. The low fractional abundances of individual species prevent their detection on the basis of spectral fingerprints in the UV.

preprint2012arXiv

Transient dust in warm debris disks - Detection of Fe-rich olivine grains

(Abridged) Debris disks trace remnant reservoirs of leftover planetesimals in planetary systems. A handful of "warm" debris disks have been discovered in the last years, where emission in excess starts in the mid-infrared. An interesting subset within these warm debris disks are those where emission features are detected in mid-IR spectra, which points towards the presence of warm micron-sized dust grains. Given the ages of the host stars, the presence of these grains is puzzling, and questions their origin and survival in time. This study focuses on determining the mineralogy of the dust around 7 debris disks with evidence for warm dust, based on Spitzer/IRS spectroscopic data, in order to provide new insights into the origin of the dust grains. We present a new radiative transfer code dedicated to SED modeling of optically thin disks. We make use of this code on the SEDs of seven warm debris disks, in combination with recent laboratory experiments on dust optical properties. We find that most, if not all, debris disks in our sample are experiencing a transient phase, suggesting a production of small dust grains on relatively short timescales. From a mineralogical point of view, we find that enstatite grains have small abundances compared to crystalline olivine grains. The main result of our study is that we find evidences for Fe-rich crystalline olivine grains (Fe / [Mg + Fe] ~ 0.2) for several debris disks. This finding contrasts with studies of gas-rich protoplanetary disks. The presence of Fe-rich olivine grains, and the overall differences between the mineralogy of dust in Class II disks compared to debris disks suggest that the transient crystalline dust is of a new generation. We discuss possible crystallization routes to explain our results, and comment on the mechanisms that may be responsible for the production of small dust grains.

preprint2012arXiv

WASP-4b Transit Observations With GROND

Ground-based simultaneous multiband transit observations allow an accurate system parameters determination and may lead to the detection and characterization of additional bodies via the transit timing variations (TTVs) method. We aimed to (i) characterize the heavily bloated WASP-4b hot Jupiter and its star by measuring system parameters and the dependence of the planetary radius as a function of four (Sloan g', r', i', z') wavelengths and (ii) search for TTVs. We recorded 987 images during three complete transits with the GROND instrument, mounted on the MPG/ESO-2.2m telescope at La Silla Observatory. Assuming a quadratic law for the stellar limb darkening we derive system parameters by fitting a composite transit light curve over all bandpasses simultaneously. To compute uncertainties of the fitted parameters, we employ the Bootstrap Monte Carlo Method. The three central transit times are measured with precision down to 6 s. We find a planetary radius Rp = 1.413+/-0.020 RJup, an orbital inclination i = 88deg57'+/-0.45deg and calculate a new ephemeris, a period P = 1.33823144+/-0.00000032days and a reference transit epoch T0 = 2454697.798311+/-0.000046 (BJD). Analysis of the new transit mid-times in combination with previous measurements shows no sign of a TTV signal greater than 20s. We perform simplified numerical simulations to place upper-mass limits of a hypothetical perturber in the WASP-4b system.

preprint2011arXiv

3D MHD Simulations of Planet Migration in Turbulent Stratified Disks

We performed 3D MHD simulations of planet migration in stratified disks using the Godunov code PLUTO, where the disk is turbulent due to the magnetorotational instability. We study the migration for planets with different planet-star mass ratios $q=M_{p}/M_{s}$. In agreement with previous studies, for the low-mass planet cases ($q=5\times10^{-6}$ and $10^{-5}$), migration is dominated by random fluctuations in the torque. For a Jupiter-mass planet $(q=M_{p}/M_{s}=10^{-3}$ for $M_{s}=1M_{\odot})$, we find a reduction of the magnetic stress inside the orbit of the planet and around the gap region. After an initial stage where the torque on the planet is positive, it reverses and we recover migration rates similar to those found in disks where the turbulent viscosity is modelled by an $α$ viscosity. For the intermediate-mass planets ($q=5\times10^{-5}, 10^{-4}$ and $2\times10^{-4}$) we find a new and so far unexpected behavior. In some cases they experience sustained and systematic outwards migration for the entire duration of the simulation. For this case, the horseshoe region is resolved and torques coming from the corotation region can remain unsaturated due to the stresses in the disk. These stresses are generated directly by the magnetic field. The magnitude of the horseshoe drag can overcome the negative Lindblad contribution when the local surface density profile is flat or increasing outwards, which we see in certain locations in our simulations due to the presence of a zonal flow. The intermediate-mass planet is migrating radially outwards in locations where there is a positive gradient of a pressure bump (zonal flow).

preprint2011arXiv

A Herschel Search For Cold Dust in Brown Dwarf Disks: First Results

We report initial results from a {\it Herschel} program to search for far-infrared emission from cold dust around a statistically significant sample of young brown dwarfs. The first three objects in our survey are all detected at 70\micron, and we report the first detection of a brown dwarf at 160\micron. The flux densities are consistent with the presence of substantial amounts of cold dust in the outer disks around these objects. We modeled the SED's with two different radiative transfer codes. We find that a broad range of model parameters provides a reasonable fit to the SED's, but that the addition of our 70\micron, and especially the 160\micron\ detection enables strong lower limits to be placed on the disk masses since most of the mass is in the outer disk. We find likely disk masses in the range of a few $\times 10^{-6}$ to $10^{-4}$ \msun. Our models provide a good fit to the SED's and do not require dust settling.

preprint2011arXiv

Can neutral and ionized PAHs be carriers of the UV extinction bump and the diffuse interstellar bands?

Up to now, no laboratory-based study has investigated polycyclic aromatic hydrocarbon (PAH) species as potential carriers of both the diffuse interstellar bands (DIBs) and the 2175 A UV bump. We examined the proposed correlation between these two features by applying experimental and theoretical techniques on two specific medium-sized/large PAHs (dibenzorubicene C30H14 and hexabenzocoronene C42H18) in their neutral and cationic states. It was already shown that mixtures of sufficiently large, neutral PAHs can partly or even completely account for the UV bump. We investigated how the absorption bands are altered upon ionization of these molecules by interstellar UV photons. The experimental studies presented here were realized by performing matrix isolation spectroscopy with subsequent far-UV irradiation. The main effects were found to be a broadening of the absorption bands in the UV combined with slight red shifts. The position of the complete pi - pi* absorption structure around 217.5 nm, however, remains more or less unchanged which could explain the observed position invariance of the interstellar bump for different lines of sight. This favors the assignment of this feature to the interstellar PAH population. As far as the DIBs are concerned, neither our investigations nor the laboratory studies carried out by other research groups support a possible connection with this class of molecules. Instead, there are reasonable arguments that neutral and singly ionized cationic PAHs cannot be made responsible for the DIBs.

preprint2011arXiv

CID: Chemistry in disks VI.sulfur-bearing molecules in the protoplanetary disks surrounding LkCa15, MWC480, DM Tau, and GO Tau

We study the content in S-bearing molecules of protoplanetary disks around low-mass stars. We used the new IRAM 30-m receiver EMIR to perform simultaneous observations of the $1_{10}-1_{01}$ line of H$_2$S at 168.8 GHz and $2_{23}-1_{12}$ line of SO at 99.3 GHz. We compared the observational results with predictions coming from the astrochemical code NAUTILUS, which has been adapted to protoplanetary disks. The data were analyzed together with existing CS J=3-2 observations. We fail to detect the SO and H$_2$S lines, although CS is detected in LkCa15, DM\,Tau, and GO\,Tau but not in MWC\,480. However, our new upper limits are significantly better than previous ones and allow us to put some interesting constraints on the sulfur chemistry. Our best modeling of disks is obtained for a C/O ratio of 1.2, starting from initial cloud conditions of H density of $2\times 10^5$ cm$^{-3}$ and age of $10^6$ yr. The results agree with the CS data and are compatible with the SO upper limits, but fail to reproduce the H$_2$S upper limits. The predicted H$_2$S column densities are too high by at least one order of magnitude. H$_2$S may remain locked onto grain surfaces and react with other species, thereby preventing the desorption of H$_2$S.

preprint2011arXiv

Collapse of massive magnetized dense cores using radiation-magneto-hydrodynamics: early fragmentation inhibition

We report the results of radiation-magneto-hydrodynamics calculations in the context of high mass star formation, using for the first time a self-consistent model for photon emission (i.e. via thermal emission and in radiative shocks) and with the high resolution necessary to resolve properly magnetic braking effects and radiative shocks on scales <100 AU. We investigate the combined effects of magnetic field, turbulence, and radiative transfer on the early phases of the collapse and the fragmentation of massive dense cores. We identify a new mechanism that inhibits initial fragmentation of massive dense cores, where magnetic field and radiative transfer interplay. We show that this interplay becomes stronger as the magnetic field strength increases. Magnetic braking is transporting angular momentum outwards and is lowering the rotational support and is thus increasing the infall velocity. This enhances the radiative feedback owing to the accretion shock on the first core. We speculate that highly magnetized massive dense cores are good candidates for isolated massive star formation, while moderately magnetized massive dense cores are more appropriate to form OB associations or small star clusters.

preprint2011arXiv

Disentangling the Environment of the FU Orionis Candidate HBC 722 with Herschel

We analyze the submillimeter emission surrounding the new FU Orionis-type object, HBC 722. We present the first epoch of observations of the active environs of HBC 722, with imaging and spectroscopy from PACS, SPIRE, and HIFI aboard the Herschel Space Observatory, as well as CO J= 2-1 and 350 um imaging (SHARC-II) with the Caltech Submillimeter Observatory. The primary source of submillimeter continuum emission in the region -- 2MASS 20581767+4353310 -- is located 16$\arcsec$ south-southeast of the optical flaring source while the optical and near-IR emission is dominated by HBC 722. A bipolar outflow extends over HBC 722; the most likely driver is the submillimeter source. We detect warm (100 K) and hot (246 K) CO emission in the surrounding region, evidence of outflow-driven heating in the vicinity. The region around HBC 722 itself shows little evidence of heating driven by the new outbursting source itself.

preprint2011arXiv

Disk Evolution in OB Associations - Deep Spitzer/IRAC Observations of IC 1795

We present a deep Spitzer/IRAC survey of the OB association IC 1795 carried out to investigate the evolution of protoplanetary disks in regions of massive star formation. Combining Spitzer/IRAC data with Chandra/ACIS observations, we find 289 cluster members. An additional 340 sources with an infrared excess, but without X-ray counterpart, are classified as cluster member candidates. Both surveys are complete down to stellar masses of about 1 Msun. We present pre-main sequence isochrones computed for the first time in the Spitzer/IRAC colors. The age of the cluster, determined via the location of the Class III sources in the [3.6]-[4.5]/[3.6] color-magnitude diagram, is in the range of 3 - 5 Myr. As theoretically expected, we do not find any systematic variation in the spatial distribution of disks within 0.6 pc of either O-type star in the association. However, the disk fraction in IC 1795 does depend on the stellar mass: sources with masses >2 Msun have a disk fraction of ~20%, while lower mass objects (2-0.8 Msun) have a disk fraction of ~50%. This implies that disks around massive stars have a shorter dissipation timescale.

preprint2011arXiv

Dissecting high-mass star-forming regions; tracing back their complex formation history

We present near-infrared JHKs imaging as well as K-band multi-object spectroscopy of the massive stellar content of W3 Main using LUCI at the LBT. We confirm 13 OB stars by their absorption line spectra in W3 Main and spectral types between O5V and B4V have been found. Three massive Young Stellar Objects are identified by their emission line spectra and near-infrared excess. From our spectrophotometric analysis of the massive stars and the nature of their surrounding HII regions we derive the evolutionary sequence of W3 Main and we find evidence of an age spread of at least 2-3 Myr. While the most massive star (IRS2) is already evolved, indications for high-mass pre--main-sequence evolution is found for another star (IRS N1), deeply embedded in an ultra compact HII region, in line with the different evolutionary phases observed in the corresponding HII regions. We have detected the photospheres of OB stars from the more evolved diffuse HII region to the much younger UCHII regions, suggesting that the OB stars have finished their formation and cleared away their possible circumstellar disks very fast. Only in the hyper-compact HII region (IRS5), the early type stars are still surrounded by circumstellar material.

preprint2011arXiv

Dust properties and disk structure of evolved protoplanetary disks in Cep OB2: Grain growth, settling, gas and dust mass, and inside-out evolution

We present Spitzer/IRS spectra of 31 TTS and IRAM/1.3mm observations for 34 low- and intermediate-mass stars in the Cep OB2 region. Including our previously published data, we analyze 56 TTS and the 3 intermediate-mass stars with silicate features in Tr 37 (~4 Myr) and NGC 7160 (~12 Myr). The silicate emission features are well reproduced with a mixture of amorphous (with olivine, forsterite, and silica stoichiometry) and crystalline grains (forsterite, enstatite). We explore grain size and disk structure using radiative transfer disk models, finding that most objects have suffered substantial evolution (grain growth, settling). About half of the disks show inside-out evolution, with either dust-cleared inner holes or a radially-dependent dust distribution, typically with larger grains and more settling in the innermost disk. The typical strong silicate features require nevertheless the presence of small dust grains, and could be explained by differential settling according to grain size, anomalous dust distributions, and/or optically thin dust populations within disk gaps. M-type stars tend to have weaker silicate emission and steeper SEDs than K-type objects. The inferred low dust masses are in a strong contrast with the relatively high gas accretion rates, suggesting global grain growth and/or an anomalous gas to dust ratio. Transition disks (TD) in the Cep OB2 region display strongly processed grains, suggesting that they are dominated by dust evolution and settling. Finally, the presence of rare but remarkable disks with strong accretion at old ages reveals that some very massive disks may still survive to grain growth, gravitational instabilities, and planet formation.

preprint2011arXiv

HD 135344B: a young star has reached its rotational limit

We search for periodic variations in the radial velocity of the young Herbig star HD 135344B with the aim to determine a rotation period. We analyzed 44 high-resolution optical spectra taken over a time range of 151 days. The spectra were acquired with FEROS at the 2.2m MPG/ESO telescope in La Silla. The stellar parameters of HD 135344B are determined by fitting synthetic spectra to the stellar spectrum. In order to obtain radial velocity measurements, the stellar spectra have been cross-correlated with a theoretical template computed from determined stellar parameters. We report the first direct measurement of the rotation period of a Herbig star from radial-velocity measurements. The rotation period is found to be 0.16 d (3.9 hr), which makes HD 135344B a rapid rotator at or close to its break-up velocity. The rapid rotation could explain some of the properties of the circumstellar environment of HD 135344B such as the presence of an inner disk with properties (composition, inclination), that are significantly different from the outer disk.

preprint2011arXiv

Herschel Observations of the T Cha Transition Disk: Constraining the Outer Disk Properties

T Cha is a nearby (d = 100 pc) transition disk known to have an optically thin gap separating optically thick inner and outer disk components. Huelamo et al. (2011) recently reported the presence of a low-mass object candidate within the gap of the T Cha disk, giving credence to the suspected planetary origin of this gap. Here we present the Herschel photometry (70, 160, 250, 350, and 500 micron) of T Cha from the "Dust, Ice, and Gas in Time" (DIGIT) Key Program, which bridges the wavelength range between existing Spitzer and millimeter data and provide important constraints on the outer disk properties of this extraordinary system. We model the entire optical to millimeter wavelength spectral energy distribution (SED) of T Cha (19 data points between 0.36 and 3300 micron without any major gaps in wavelength coverage). T Cha shows a steep spectral slope in the far-IR, which we find clearly favors models with outer disks containing little or no dust beyond 40 AU. The full SED can be modeled equally well with either an outer disk that is very compact (only a few AU wide) or a much larger one that has a very steep surface density profile. That is, T Cha's outer disk seems to be either very small or very tenuous. Both scenarios suggest a highly unusual outer disk and have important but different implications for the nature of T Cha. Spatially resolved images are needed to distinguish between the two scenarios.

preprint2011arXiv

High-resolution simulations of planetesimal formation in turbulent protoplanetary discs

We present high-resolution computer simulations of dust dynamics and planetesimal formation in turbulence generated by the magnetorotational instability. We show that the turbulent viscosity associated with magnetorotational turbulence in a non-stratified shearing box increases when going from 256^3 to 512^3 grid points in the presence of a weak imposed magnetic field, yielding a turbulent viscosity of $α\approx0.003$ at high resolution. Particles representing approximately meter-sized boulders concentrate in large-scale high-pressure regions in the simulation box. The appearance of zonal flows and particle concentration in pressure bumps is relatively similar at moderate (256^3) and high (512^3) resolution. In the moderate-resolution simulation we activate particle self-gravity at a time when there is little particle concentration, in contrast with previous simulations where particle self-gravity was activated during a concentration event. We observe that bound clumps form over the next ten orbits, with initial birth masses of a few times the dwarf planet Ceres. At high resolution we activate self-gravity during a particle concentration event, leading to a burst of planetesimal formation, with clump masses ranging from a significant fraction of to several times the mass of Ceres. We present a new domain decomposition algorithm for particle-mesh schemes. Particles are spread evenly among the processors and the local gas velocity field and assigned drag forces are exchanged between a domain-decomposed mesh and discrete blocks of particles. We obtain good load balancing on up to 4096 cores even in simulations where particles sediment to the mid-plane and concentrate in pressure bumps.

preprint2011arXiv

Infrared spectroscopy of intermediate mass young stellar objects

In this paper we present Spitzer Infrared Spectrograph spectroscopy for 14 intermediate-mass young stellar objects. We use Spitzer spectroscopy to investigate the physical properties of these sources and their environments. Our sample can be divided into two types of objects: young isolated, embedded objects with spectra that are dominated by ice and silicate absorption bands, and more evolved objects that are dominated by extended emission from polycyclic aromatic hydrocarbons (PAHs) and pure H2 rotational lines. We are able to constrain the illuminating FUV fields by classifying the PAH bands below 9micron. For most of the sources we are able to detect several atomic fine structure lines. In particular, the [NeII] line appearing in two regions could originate from unresolved photodissociation regions (PDRs) or J-shocks. We relate the identified spectral features to observations obtained from NIR through submillimeter imaging. The spatial extent of several H2 and PAH bands is matched with morphologies identified in previous Spitzer/IRAC observations. This also allows us to distinguish between the different H2 excitation mechanisms. In addition, we calculate the optical extinction from the silicate bands and use this to constrain the spectral energy distribution fit, allowing us to estimate the masses of these YSOs.

preprint2011arXiv

Measuring Organic Molecular Emission in Disks with Low Resolution Spitzer Spectroscopy

We explore the extent to which Spitzer IRS spectra taken at low spectral resolution can be used in quantitative studies of organic molecular emission from disks surrounding low mass young stars. We use Spitzer IRS spectra taken in both the high and low resolution modules for the same sources to investigate whether it is possible to define line indices that can measure trends in the strength of the molecular features in low resolution data. We find that trends in HCN emission strength seen in the high resolution data can be recovered in low resolution data. In examining the factors that influence the HCN emission strength, we find that the low resolution HCN flux is modestly correlated with stellar accretion rate and X-ray luminosity. Correlations of this kind are perhaps expected based on recent observational and theoretical studies of inner disk atmospheres. Our results demonstrate the potential of using the large number of low resolution disk spectra that reside in the Spitzer archive to study the factors that influence the strength of molecular emission from disks. Such studies would complement results for the much smaller number of circumstellar disks that have been observed at high resolution with IRS.

preprint2011arXiv

MIDI interferometry of massive YSOs: Updates on the MPIA programme

Information about the inner structure of young stellar objects is crucial for understanding how the central forming stars gain their mass. However, especially for observations of (usually further away) high-mass young stellar objects, conventional imaging has limitations in spatial resolution. By means of mid-infrared interferometry, we can peek deeply into the strongly extinguished central 100 AU of such objects. Here, we report on new data we obtained within our programme using MIDI at the VLTI. Shown are preliminary results on the known outflow sources AFGL 2136 IRS 1 and Mon R2 IRS 3A. In particular, we describe how quantities like differential phases give additional geometrical structure information. We demonstrate how the combined interpretation of single-telescope and interferometric data at different wavelength regimes can lead to a more complete picture of the nature of such MYSOs.

preprint2011arXiv

Observations and modeling of the massive young star AFGL 4176: From large scales to small

We present spatially-resolved mid-infrared interferometric observations of the massive young stellar object AFGL 4176, together with literature and survey data. We model these observations using a simple, spherically-symmetric radiative transfer model, and find that the observational data are consistent with a highly luminous star surrounded by a thick envelope.

preprint2011arXiv

On the stability of radiation-pressure-dominated cavities

Context: When massive stars exert a radiation pressure onto their environment that is higher than their gravitational attraction, they launch a radiation-pressure-driven outflow. It has been claimed that a radiative Rayleigh-Taylor instability should lead to the collapse of the outflow cavity and foster the growth of massive stars. Aims: We investigate the stability of idealized radiation-pressure-dominated cavities, focusing on its dependence on the radiation transport approach for the stellar radiation feedback. Methods: We compare two different methods for stellar radiation feedback: gray flux-limited diffusion (FLD) and ray-tracing (RT). We also derive simple analytical models to support our findings. Results: Only the FLD cases lead to prominent instability in the cavity shell. The RT cases do not show such instability. The gray FLD method underestimates the opacity at the location of the cavity shell and leads to extended epochs of marginal Eddington equilibrium in the cavity shell, making them prone to the radiative Rayleigh-Taylor instability. In the RT cases, the radiation pressure exceeds gravity by 1-2 orders of magnitude. The radiative Rayleigh-Taylor instability is then consequently suppressed. Conclusions: Treating the stellar irradiation in the gray FLD approximation underestimates the radiative forces acting on the cavity shell. This can lead artificially to situations that are affected by the radiative Rayleigh-Taylor instability. The proper treatment of direct stellar irradiation by massive stars is crucial for the stability of radiation-pressure-dominated cavities.

preprint2011arXiv

Piercing the Glare: Direct Imaging Search for Planets in the Sirius System

Astrometric monitoring of the Sirius binary system over the past century has yielded several predictions for an unseen third system component, the most recent one suggesting a \leq50 MJup object in a ~6.3-year orbit around Sirius A. Here we present two epochs of high-contrast imaging observations performed with Subaru IRCS and AO188 in the 4.05 \mum narrow-band Br alpha filter. These data surpass previous observations by an order of magnitude in detectable companion mass, allowing us to probe the relevant separation range down to the planetary mass regime (6-12 M_Jup at 1", 2-4 M_Jup at 2", and 1.6 M_Jup beyond 4"). We complement these data with one epoch of M-band observations from MMT/AO Clio, which reach comparable performance. No dataset reveals any companion candidates above the 5-sigma level, allowing us to refute the existence of Sirius C as suggested by the previous astrometric analysis. Furthermore, our Br alpha photometry of Sirius B confirms the lack of an infrared excess beyond the white dwarf's blackbody spectrum.

preprint2011arXiv

Pre--Main-Sequence stellar populations across Shapley Constellation III. I. Photometric Analysis and Identification

We present our investigation of pre--main-sequence (PMS) stellar populations in the Large Magellanic Cloud (LMC) from imaging with Hubble Space Telescope WFPC2 camera. Our targets of interest are four star-forming regions located at the periphery of the super-giant shell LMC 4 (Shapley Constellation III). The PMS stellar content of the regions is revealed through the differential Hess diagrams and the observed color-magnitude diagrams (CMDs). Further statistical analysis of stellar distributions along cross-sections of the faint part of the CMDs allowed the quantitative assessment of the PMS stars census, and the isolation of faint PMS stars as the true low-mass stellar members of the regions. These distributions are found to be well represented by a double Gaussian function, the first component of which represents the main-sequence field stars and the second the native PMS stars of each region. Based on this result, a cluster membership probability was assigned to each PMS star according to its CMD position. The higher extinction in the region LH 88 did not allow the unambiguous identification of its native stellar population. The CMD distributions of the PMS stars with the highest membership probability in the regions LH 60, LH 63 and LH 72 exhibit an extraordinary similarity among the regions, suggesting that these stars share common characteristics, as well as common recent star formation history. Considering that the regions are located at different areas of the edge of LMC 4, this finding suggests that star formation along the super-giant shell may have occurred almost simultaneously.

preprint2011arXiv

Protostars and stars in the Coronet cluster: Age, evolution, and cluster structure

We present new optical spectroscopy with FLAMES/VLT, near-IR imaging with HAWK-I/VLT, and 870 micron mapping with APEX/LABOCA of the Coronet cluster. The optical data allow to estimate spectral types, extinction and the presence of accretion in 6 more M-type members, in addition to the 12 that we had previously studied. The submillimeter maps and near-IR data reveal the presence of nebular structures and high extinction regions, which are in some cases associated to known IR, optical, and X-ray sources. Most star formation is associated to two elongated structures crossing in the central part of the cluster. Placing all the 18 objects with known spectral types and extinction in the HR diagram suggests that the cluster is younger than previously thought (<2 Myr, and probably ~0.5-1 Myr). The new age estimate is in agreement with the evolutionary status of the various protostars in the region and with its compactness (<1.3 pc across), but results in a conflict with the low disk and accretion fraction (only 50-65% of low-mass stars appear to have protoplanetary disks, and most transitional and homologously depleted disks are consistent with no accretion) and with the evolutionary features observed in the mid-IR spectra and spectral energy distributions of the disks.

preprint2011arXiv

Radiation pressure feedback in the formation of massive stars

We investigate the radiation pressure feedback in the formation of massive stars in 1, 2, and 3D radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. In contrast to previous research, we consider frequency dependent stellar radiation feedback, resolve the dust sublimation front in the vicinity of the forming star down to 1.27 AU, compute the evolution for several 10^5 yrs covering the whole accretion phase of the forming star, and perform a comprehensive survey of the parameter space. The most fundamental result is that the formation of a massive accretion disk in slowly rotating cores preserves a high anisotropy in the radiation field. The thermal radiation escapes through the optically thin atmosphere, effectively diminishing the radiation pressure feedback onto the accretion flow. Gravitational torques in the self-gravitating disk drive a sufficiently high accretion rate to overcome the residual radiation pressure. Simultaneously, the radiation pressure launches an outflow in the bipolar direction, which grows in angle with time and releases a substantial fraction of the initial core mass from the star-disk system. Summarized, for an initial core mass of 60, 120, 240, and 480 Msol these mechanisms allow the star to grow up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

preprint2011arXiv

Resolving the inner regions of the HD97048 circumstellar disk with VLT/NACO polarimetric differential imaging

Circumstellar disks are the cradles of planetary systems and their physical and chemical properties directly influence the planet formation process. As most planets supposedly form in the inner disk regions, i.e., within a few tens of AU, it is crucial to study circumstellar disk on these scales to constrain the conditions for planet formation. Our aims are to characterize the inner regions of the circumstellar disk around the young Herbig Ae/Be star HD97048 in polarized light. We use VLT/NACO to observe HD97048 in polarimetric differential imaging (PDI) mode in the H and Ks band. We spatially resolve the disk around HD97048 in polarized flux in both filters on scales between ~0.1"-1.0" corresponding to the inner ~16-160 AU. Fitting isophots to the flux calibrated H-band image between 13 - 14 mag/arcsec^2 and 14 - 15 mag/arcsec^2 we derive a apparent disk inclination angle of 34+-5 deg and 47+-2 deg, respectively. The disk position angle in both brightness regimes is almost identical and roughly 80 deg. Along the disk major axis the surface brightness of the polarized flux drops from ~11 mag/arcsec^2 at ~0.1" (~16 AU) to ~15.3 mag/arcsec^2 at ~1.0" (~160 AU). The brightness profiles along the major axis are fitted with power-laws falling off as ~r^(-1.78+-0.02) in H and ~r^(-2.34+-0.04) in Ks. As the surface brightness drops off more rapidly in Ks compared to H the disks becomes relatively bluer at larger separations possibly indicating changing dust grain properties as a function of radius. For the first time the inner ~0.1"-1.0" (~16-160 AU) of the surface layer of the HD97048 circumstellar disk have been imaged in scattered light demonstrating the power of ground-based imaging polarimetry. Our data fill an important gap in a large collection of existing data including resolved thermal dust and PAH emission images as well as resolved gas emission lines.

preprint2011arXiv

The 2008 outburst of EX Lup - silicate crystals in motion

EX Lup is the prototype of the EXor class of eruptive young stars. These objects show optical outbursts which are thought to be related to runaway accretion onto the star. In a previous study we observed in-situ crystal formation in the disk of EX Lup during its latest outburst in 2008, making the object an ideal laboratory to investigate circumstellar crystal formation and transport. This outburst was monitored by a campaign of ground-based and Spitzer Space Telescope observations. Here we modeled the spectral energy distribution of EX Lup in the outburst from optical to millimeter wavelengths with a 2D radiative transfer code. Our results showed that the shape of the SED at optical wavelengths was more consistent with a single temperature blackbody than a temperature distribution. We also found that this single temperature component emitted 80-100 % of the total accretion luminosity. We concluded that a thermal instability, the most widely accepted model of EXor outbursts, was likely not the triggering mechanism of the 2008 outburst of EX Lup. Our mid-infrared Spitzer spectra revealed that the strength of all crystalline bands between 8 and 30 um increased right after the end of the outburst. Six months later, however, the crystallinity in the 10 um silicate feature complex decreased. Our modeling of the mid-infrared spectral evolution of EXLup showed that, although vertical mixing should be stronger during the outburst than in the quiescent phase, fast radial transport of crystals (e.g., by stellar/disk wind) was required to reproduce the observed mid-infrared spectra.

preprint2011arXiv

The alignment of molecular cloud magnetic fields with the spiral arms in M33

The formation of molecular clouds, which serve as stellar nurseries in galaxies, is poorly understood. A class of cloud formation models suggests that a large-scale galactic magnetic field is irrelevant at the scale of individual clouds, because the turbulence and rotation of a cloud may randomize the orientation of its magnetic field. Alternatively, galactic fields could be strong enough to impose their direction upon individual clouds, thereby regulating cloud accumulation and fragmentation, and affecting the rate and efficiency of star formation. Our location in the disk of the Galaxy makes an assessment of the situation difficult. Here we report observations of the magnetic field orientation of six giant molecular cloud complexes in the nearby, almost face-on, galaxy M33. The fields are aligned with the spiral arms, suggesting that the large-scale field in M33 anchors the clouds.

preprint2011arXiv

The Initial Mass Function of the Orion Nebula Cluster across the H-burning limit

We present a new census of the Orion Nebula Cluster (ONC) over a large field of view (>30'x30'), significantly increasing the known population of stellar and substellar cluster members with precisely determined properties. We develop and exploit a technique to determine stellar effective temperatures from optical colors, nearly doubling the previously available number of objects with effective temperature determinations in this benchmark cluster. Our technique utilizes colors from deep photometry in the I-band and in two medium-band filters at lambda~753 and 770nm, which accurately measure the depth of a molecular feature present in the spectra of cool stars. From these colors we can derive effective temperatures with a precision corresponding to better than one-half spectral subtype, and importantly this precision is independent of the extinction to the individual stars. Also, because this technique utilizes only photometry redward of 750nm, the results are only mildly sensitive to optical veiling produced by accretion. Completing our census with previously available data, we place some 1750 sources in the Hertzsprung-Russel diagram and assign masses and ages down to 0.02 solar masses. At faint luminosities, we detect a large population of background sources which is easily separated in our photometry from the bona fide cluster members. The resulting initial mass function of the cluster has good completeness well into the substellar mass range, and we find that it declines steeply with decreasing mass. This suggests a deficiency of newly formed brown dwarfs in the cluster compared to the Galactic disk population.

preprint2011arXiv

The substellar mass function in the central region of the open cluster Praesepe from deep LBT observations

Studies of the mass function (MF) of open clusters of different ages allow us to probe the efficiency with which brown dwarfs evaporate from clusters to populate the field. Surveys of older clusters (age >100 Myr) are not affected so severely by several problems encountered in young clusters, such as intra-cluster extinction and large uncertainties in brown dwarf models. We present the results of a deep photometric survey to study the MF of the central region of the old open cluster Praesepe (age~590$^{+150}_{-120}$ Myr, distance ~190$^{+6.0}_{-5.8}$ pc), down to the substellar regime. We performed an optical ($riz$ and $Y$-band) photometric survey of Praesepe using the Large Binocular Telescope Camera covering an area of ~0.59 deg$^2$ in the cluster centre from $i$~19.0 mag (~100 M_Jup) down to a 5$σ$ detection limit at $i$~25.6mag (~40 M_Jup). The survey is approximately 95% complete at $i=23.8$ mag and $z=22.0$ mag (~55 M_Jup). We identify 59 cluster member candidates, of which 37 are substellar, by comparing with the predictions of a dusty atmosphere model. The MF of those candidates rises from the substellar boundary until ~67 M_Jup and then declines. This is quite different from the form inferred for other open clusters older than 50 Myr, but seems to be similar to those found in very young open clusters, the MFs of which peak at ~10 M_Jup. Either Praesepe really does have a different MF from other clusters or they had similar initial MFs but a different dynamical evolution. Since most of the candidates are faint, we lack astrometric or spectroscopic follow-ups to test their memberships. However, the contaminations by field dwarfs, galaxies, or giants are found to have little effect on the shape of MF and therefore the MF of `real' cluster members should have similar characteristics.

preprint2011arXiv

The Ubiquity of Micrometer-Sized Dust Grains in the Dense Interstellar Medium

Cold molecular clouds are the birthplaces of stars and planets, where dense cores of gas collapse to form protostars. The dust mixed in these clouds is thought to be made of grains of an average size of 0.1 micrometer. We report the widespread detection of the coreshine effect as a direct sign of the existence of grown, micrometer-sized dust grains. This effect is seen in half of the cores we have analyzed in our survey, spanning all Galactic longitudes, and is dominated by changes in the internal properties and local environment of the cores, implying that the coreshine effect can be used to constrain fundamental core properties such as the three-dimensional density structure and ages and also the grain characteristics themselves.

preprint2011arXiv

Three-dimensional simulation of massive star formation in the disk accretion scenario

The most massive stars can form via standard disk accretion - despite of the radiation pressure generated - due to the fact that the massive accretion disk yields a strong anisotropy in the radiation field, releasing most of the radiation pressure perpendicular to the disk accretion flow. Here, we analyze the self-gravity of the forming circumstellar disk as the potential major driver of the angular momentum transport in such massive disks responsible for the high accretion rates needed for the formation of massive stars. For this purpose, we perform self-gravity radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. The formation and evolution of the resulting circumstellar disk is investigated in 1.) axially symmetric simulations using an alpha-shear-viscosity prescription and 2.) a three-dimensional simulation, in which the angular momentum transport is provided self-consistently by developing gravitational torques in the self-gravitating accretion disk. The simulation series of different strength of the alpha-viscosity shows that the accretion history of the forming star is mostly independent of the alpha-viscosity-parameter. The accretion history of the three-dimensional run driven by self-gravity is more time-dependent than the viscous disk evolution in axial symmetry. The mean accretion rate, i.e. the stellar mass growth, is nearly identical to the alpha-viscosity models. We conclude that the development of gravitational torques in self-gravitating disks around forming massive stars provides a self-consistent mechanism to efficiently transport the angular momentum to outer disk radii. Also the formation of the most massive stars can therefore be understood in the standard accretion disk scenario.

preprint2011arXiv

VLT/NACO Polarimetric Differential Imaging of HD100546 - Disk Structure and Dust Grain Properties between 10-140 AU

We present polarimetric differential imaging (PDI) data of the circumstellar disk around the Herbig Ae/Be star HD100546 obtained with VLT/NACO. We resolve the disk in polarized light in the H and Ks filter between ~0.1-1.4" (i.e., ~10-140 AU). The innermost disk regions are directly imaged for the first time and the mean apparent disk inclination and position angle are derived. The surface brightness along the disk major axis drops off roughly with S(r) ~ r^-3 but has a maximum around 0.15" suggesting a marginal detection of the main disk inner rim at ~15 AU. We find a significant brightness asymmetry along the disk minor axis in both filters with the far side of the disk appearing brighter than the near side. This enhanced backward scattering and a low total polarization degree of the scattered disk flux of 14%(+19%/-8%) suggests that the dust grains on the disk surface are larger than typical ISM grains. Empirical scattering functions reveal the backward scattering peak at the largest scattering angles and a second maximum for the smallest scattering angles. This indicates a second dust grain population preferably forward scattering and smaller in size. It shows that, relatively, in the inner disk regions (40-50 AU) a higher fraction of larger grains is found compared to the outer disk regions (100-110 AU). Finally, our images reveal distinct substructures between 25-35 AU physical separation from the star and we discuss the possible origin for the two features in the context of ongoing planet formation.

preprint2010arXiv

A Giant Planet Around a Metal-poor Star of Extragalactic Origin

Stars in their late stage of evolution, such as Horizontal Branch stars, are still largely unexplored for planets. We report the detection of a planetary companion around HIP 13044, a very metal-poor star on the red Horizontal Branch, based on radial velocity observations with a high-resolution spectrograph at the 2.2-m MPG/ESO telescope. The star's periodic radial velocity variation of P=16.2 days caused by the planet can be distinguished from the periods of the stellar activity indicators. The minimum mass of the planet is 1.25 Jupiter masses and its orbital semi-major axis 0.116 AU. Because HIP 13044 belongs to a group of stars that have been accreted from a disrupted satellite galaxy of the Milky Way, the planet most likely has an extragalactic origin.

preprint2010arXiv

A ground-based near-infrared emission spectrum of the exoplanet HD 189733b

Detection of molecules using infrared spectroscopy probes the conditions and compositions of exoplanet atmospheres. Water (H2O), methane (CH4), carbon dioxide (CO2), and carbon monoxide (CO) have been detected in two hot Jupiters. These previous results relied on space-based telescopes that do not provide spectroscopic capability in the 2.4 - 5.2 micron spectral region. Here we report ground-based observations of the dayside emission spectrum for HD 189733b between 2.0-2.4 micron and 3.1-4.1 micron, where we find a bright emission feature. Where overlap with space-based instruments exists, our results are in excellent agreement with previous measurements. A feature at ~3.25 micron is unexpected and difficult to explain with models that assume local thermodynamic equilibrium (LTE) conditions at the 1 bar to 1 x 10-6 bar pressures typically sampled by infrared measurements. The most likely explanation for this feature is that it arises from non-LTE emission from CH4, similar to what is seen in the atmospheres of planets in our own Solar System. These results suggest that non-LTE effects may need to be considered when interpreting measurements of strongly irradiated exoplanets.

preprint2010arXiv

A large, massive, rotating disk around an isolated young stellar object

We present multi-wavelengths observations and a radiative transfer model of a newly discovered massive circumstellar disk of gas and dust which is one of the largest disks known today. Seen almost edge-on, the disk is resolved in high-resolution near-infrared (NIR) images and appears as a dark lane of high opacity intersecting a bipolar reflection nebula. Based on molecular line observations we estimate the distance to the object to be 3.5 kpc. This leads to a size for the dark lane of ~10500 AU but due to shadowing effects the true disk size could be smaller. In Spitzer/IRAC 3.6 micron images the elongated shape of the bipolar reflection nebula is still preserved and the bulk of the flux seems to come from disk regions that can be detected due to the slight inclination of the disk. At longer IRAC wavelengths, the flux is mainly coming from the central regions penetrating directly through the dust lane. Interferometric observations of the dust continuum emission at millimeter wavelengths with the SMA confirm this finding as the peak of the unresolved mm-emission coincides perfectly with the peak of the Spitzer/IRAC 5.8 micron flux and the center of the dark lane seen in the NIR images. Simultaneously acquired CO data reveal a molecular outflow along the northern part of the reflection nebula which seems to be the outflow cavity. An elongated gaseous disk component is also detected and shows signs of rotation. The emission is perpendicular to the molecular outflow and thus parallel to but even more extended than the dark lane in the NIR images. Based on the dust continuum and the CO observations we estimate a disk mass of up to a few solar masses depending on the underlying assumptions. Whether the disk-like structure is an actual accretion disk or rather a larger-scale flattened envelope or pseudodisk is difficult to discriminate with the current dataset (abridged).

preprint2010arXiv

Accretion in Evolved and Transitional Disks in Cep OB2: Looking for the Origin of the Inner Holes

We present accretion rates for a large number of solar-type stars in the Cep OB2 region, based on U band observations. Our study comprises 95 members of the ~4 Myr-old cluster Tr 37 (including 20 "transition" objects; TO), as well as the only CTTS in the ~12 Myr-old cluster NGC 7160. The stars show different disk morphologies, with the majority of them having evolved and flattened disks. The typical accretion rates are about one order of magnitude lower than in regions aged 1-2 Myr, and we find no strong correlation between disk morphology and accretion rates. Although half of the TO are not accreting, the median accretion rates of normal CTTS and accreting "transition" disks are similar (~3 10^{-9} and 2 10^{-9} Msun/yr, respectively). Comparison with other regions suggests that the TO observed at different ages do not necessarily represent the same type of objects, which is consistent with the fact that the different processes that can lead to reduced IR excess/inner disk clearing (e.g., binarity, dust coagulation/settling, photoevaporation, giant planet formation) do not operate on the same timescales. Accreting TO in Tr 37 are probably suffering strong dust coagulation/settling. Regarding the equally large number of non-accreting TO in the region, other processes, like photoevaporation, the presence of stellar/substellar companions, and/or giant planet formation may account for their "transitional" SEDs and negligible accretion rates.

preprint2010arXiv

Circumventing the radiation pressure barrier in the formation of massive stars via disk accretion

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase (up to the end of the accretion disk epoch) for the forming massive star and to perform a broad scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. While confirming the upper mass limit of spherically symmetric accretion, our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field. This yields to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations, far beyond the upper mass limit of spherical accretion. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msun the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msun respectively.

preprint2010arXiv

Co-evolution of atmospheres, life, and climate

After Earth's origin, our host star, the Sun, was shining 20 to 25 percent less brightly than today. Without greenhouse-like conditions to warm the atmosphere, our early planet would have been an ice ball and life may never have evolved. But life did evolve, which indicates that greenhouse gases must have been present on early Earth to warm the planet. Evidence from the geologic record indicates an abundance of the greenhouse gas CO2. CH4 was probably present as well, and in this regard methanogenic bacteria, which belong to a diverse group of anaerobic procaryotes that ferment CO 2 plus H2 to CH4, may have contributed to modification of the early atmosphere. Molecular oxygen was not present, as is indicated by the study of rocks from that era, which contain iron carbonate rather than iron oxide. Multicellular organisms originated as cells within colonies that became increasingly specialized. The development of photosynthesis allowed the Sun's energy to be harvested directly by life forms. The resultant oxygen accumulated in the atmosphere and formed the ozone layer in the upper atmosphere. Aided by the absorption of harmful UV radiation in the ozone layer, life colonized Earth's surface. Our own planet is a very good example of how life forms modified the atmosphere over the planets' life time. We show that these facts have to be taken into account when we discover and characterize atmospheres of Earth-like exoplanets. If life has originated and evolved on a planet, then it should be expected that a strong co-evolution occurred between life and the atmosphere, the result of which is the planets' climate.

preprint2010arXiv

Discovery of a stellar companion to the nearby solar-analogue HD 104304

Sun-like stars are promising candidates to host exoplanets and are often included in exoplanet surveys by radial velocity (RV) and direct imaging. In this paper we report on the detection of a stellar companion to the nearby solar-analogue star HD 104304, which previously was considered to host a planetary mass or brown dwarf companion. We searched for close stellar and substellar companions around extrasolar planet host stars with high angular resolution imaging to characterize planet formation environments. The detection of the stellar companion was achieved by high angular resolution measurements, using the "Lucky Imaging" technique at the ESO NTT 3.5m with the AstraLux Sur instrument. We combined the results with VLT/NACO archive data, where the companion could also be detected. The results were compared to precise RV measurements of HD 104304, obtained at the Lick and Keck observatories from 2001-2010. We confirmed common proper motion of the binary system. A spectral type of M4V of the companion and a mass of 0.21 M_Sun was derived. Due to comparison of the data with RV measurements of the unconfirmed planet candidate listed in the Extrasolar Planets Encyclopaedia, we suggest that the discovered companion is the origin of the RV trend and that the inclination of the orbit of approximately 35 degrees explains the relatively small RV signal.

preprint2010arXiv

Dissecting the Moth: Discovery of an off-centered ring in the HD 61005 debris disk with high-resolution imaging

The debris disk known as "The Moth" is named after its unusually asymmetric surface brightness distribution. It is located around the ~90 Myr old G8V star HD 61005 at 34.5 pc and has previously been imaged by the HST at 1.1 and 0.6 microns. Polarimetric observations suggested that the circumstellar material consists of two distinct components, a nearly edge-on disk or ring, and a swept-back feature, the result of interaction with the interstellar medium. We resolve both components at unprecedented resolution with VLT/NACO H-band imaging. Using optimized angular differential imaging techniques to remove the light of the star, we reveal the disk component as a distinct narrow ring at inclination i=84.3 \pm 1.0°. We determine a semi-major axis of a=61.25 \pm 0.85 AU and an eccentricity of e=0.045 \pm 0.015, assuming that periastron is located along the apparent disk major axis. Therefore, the ring center is offset from the star by at least 2.75 \pm 0.85 AU. The offset, together with a relatively steep inner rim, could indicate a planetary companion that perturbs the remnant planetesimal belt. From our imaging data we set upper mass limits for companions that exclude any object above the deuterium-burning limit for separations down to 0.3". The ring shows a strong brightness asymmetry along both the major and minor axis. A brighter front side could indicate forward-scattering grains, while the brightness difference between the NE and SW components can be only partly explained by the ring center offset, suggesting additional density enhancements on one side of the ring. The swept-back component appears as two streamers originating near the NE and SW edges of the debris ring.

preprint2010arXiv

Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.

preprint2010arXiv

L1448 IRS2E: A candidate first hydrostatic core

Intermediate between the prestellar and Class 0 protostellar phases, the first core is a quasi-equilibrium hydrostatic object with a short lifetime and an extremely low luminosity. Recent MHD simulations suggest that the first core can even drive a molecular outflow before the formation of the second core (i.e., protostar). Using the Submillimeter Array and the Spitzer Space Telescope, we present high angular resolution observations towards the embedded dense core IRS2E in L1448. We find that source L1448 IRS2E is not visible in the sensitive Spitzer infrared images (at wavelengths from 3.6 to 70 um), and has weak (sub-)millimeter dust continuum emission. Consequently, this source has an extremely low bolometric luminosity (< 0.1 L_sun). Infrared and (sub-)millimeter observations clearly show an outflow emanating from this source; L1448 IRS2E represents thus far the lowest luminosity source known to be driving a molecular outflow. Comparisons with prestellar cores and Class 0 protostars suggest that L1448 IRS2E is more evolved than prestellar cores but less evolved than Class 0 protostars, i.e., at a stage intermediate between prestellar cores and Class 0 protostars. All these results are consistent with the theoretical predictions of the radiative/magneto hydrodynamical simulations, making L1448 IRS2E the most promising candidate of the first hydrostatic core revealed so far.

preprint2010arXiv

Mass segregation and elongation of the starburst cluster Westerlund 1

Massive stellar clusters are the best available laboratories to study the mass function of stars. Based on NTT/SofI near-infrared photometry, we have investigated the properties of the massive young cluster Westerlund 1. From comparison with stellar models, we derived an extinction A_{Ks} = 0.91 +/- 0.05 mag, an age τ= 4 +/- 0.5 Myr and a distance d = 4.0 +/- 0.2 kpc for Westerlund 1, as well as a total mass of M_{Wd1} = 4.91_{-0.49}^{+1.79} x 10^4 M_{sun}. Using spatially dependent completeness corrections we performed a 2D study of the cluster's IMF and, in addition, of the stellar density profiles of the cluster as a function of mass. From both IMF slope variations and stellar density, we find strong evidence of mass segregation. For a cluster with some 10^5 stars, this is not expected at such a young age as the result of two-body relaxation alone. We also confirm previous findings on the elongation of Westerlund 1; assuming an elliptical density profile, we found an axis ratio of a:b = 3:2. Rapid mass segregation and elongation could be well explained as the results of subclusters merging during the formation of Westerlund 1.

preprint2010arXiv

Optical properties of cosmic dust analogs: A review

Nanometer- and micrometer-sized solid particles play an important role in the evolutionary cycle of stars and interstellar matter. The optical properties of cosmic grains determine the interaction of the radiation field with the solids, thereby regulating the temperature structure and spectral appearance of dusty regions. Radiation pressure on dust grains and their collisions with the gas atoms and molecules can drive powerful winds. The analysis of observed spectral features, especially in the infrared wavelength range, provides important information on grain size, composition and structure as well as temperature and spatial distribution of the material. The relevant optical data for interstellar, circumstellar, and protoplanetary grains can be obtained by measurements on cosmic dust analogs in the laboratory or can be calculated from grain models based on optical constants. Both approaches have made progress in the last years, triggered by the need to interpret increasingly detailed high-quality astronomical observations. The statistical theoretical approach, spectroscopic experiments at variable temperature and absorption spectroscopy of aerosol particulates play an important role for the successful application of the data in dust astrophysics.

preprint2010arXiv

Recent Star Formation in Low Metallicities. The Star-Forming Region NGC 346/N66 in the Small Magellanic Cloud from Near-Infrared VLT/ISAAC Observations

We perform a photometric analysis on deep, seeing-limited near-IR VLT images of the region NGC 346/N66 and a nearby control field of the SMC in order to locate the centres of active high- and intermediate-mass star formation through the identification of near-IR bright objects as candidate stellar sources under formation.We use archived imaging data obtained with the high-resolution camera ISAAC at VLT of NGC 346/N66 and we construct the near-IR colour-magnitude (CMD) and colour-colour diagrams (C-CD) of all detected sources. We investigate the nature of all stellar populations in the observed CMDs, and we identify all stellar sources that show significant near-IR excess emission in the observed C-CD. We select, thus, the best candidates for being young stellar sources. Based on their near-IR colours we select 263 candidate young stellar sources. This sample comprises a variety of objects such as intermediate-mass PMS and Herbig Ae/Be stars and possibly massive YSOs, providing original near-IR colours for them. The spatial distribution of the selected candidate sources shows that they are located along the dusty filamentary structures of N66 seen in mid- and far-IR dust emission and agrees very well with that of previously detected candidate YSOs and PMS stars. Our study provides an original accurate set of near-IR colours for candidate young stellar sources. While this adds significant information on the star formation process in NGC 346/N66, it is still necessary to reveal the exact nature of these objects, a process that requires the combination of multi-wavelength data, so that complete Spectral Energy Distributions can be constructed for individual sources. This would be a quite important follow-up study to that presented here.

preprint2010arXiv

Search for Outer Massive Bodies around Transiting Planetary Systems: Candidates of Faint Stellar Companions around HAT-P-7

We present results of direct imaging observations for HAT-P-7 taken with the Subaru HiCIAO and the Calar Alto AstraLux. Since the close-in transiting planet HAT-P-7b was reported to have a highly tilted orbit, massive bodies such as giant planets, brown dwarfs, or a binary star are expected to exist in the outer region of this system. We show that there are indeed two candidates for distant faint stellar companions around HAT-P-7. We discuss possible roles played by such companions on the orbital evolution of HAT-P-7b. We conclude that as there is a third body in the system as reported by Winn et al. (2009, ApJL, 763, L99), the Kozai migration is less likely while planet-planet scattering is possible.

preprint2010arXiv

Star formation in the Taurus filament L1495: From Dense Cores to Stars

We present a study of dense structures in the L1495 filament in the Taurus Molecular Cloud and examine its star-forming properties. In particular we construct a dust extinction map of the filament using deep near-infrared observations, exposing its small-scale structure in unprecedented detail. The filament shows highly fragmented substructures and a high mass-per-length value of M_line=17 Msun/pc, reflecting star-forming potential in all parts of it. However, a part of the filament, namely B211, is remarkably devoid of young stellar objects. We argue that in this region the initial filament collapse and fragmentation is still taking place and star formation is yet to occur. In the star-forming part of the filament, we identify 39 cores with masses from 0.4...10 Msun and preferred separations in agreement with the local Jeans length. Most of these cores exceed the Bonnor-Ebert critical mass, and are therefore likely to collapse and form stars. The Dense Core Mass Function follows a power law with exponent Gamma=1.2, a form commonly observed in star-forming regions.

preprint2010arXiv

The role of accretion disks in the formation of massive stars

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase of several 10^5 yr for the forming massive star and to perform a comprehensive scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. Our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field, yielding to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msol the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

preprint2010arXiv

The seeds of star formation in the filamentary infrared-dark cloud G011.11-0.12

Infrared-dark clouds (IRDCs) are the precursors to massive stars and stellar clusters. G011.11-0.12 is a well-studied filamentary IRDC, though, to date, the absence of far-infrared data with sufficient spatial resolution has limited the understanding of the structure and star-formation activity. We use Herschel to study the embedded population of young pre- and protostellar cores in this IRDC. We examine the cloud structure, which appears in absorption at short wavelength and in emission at longer wavelength. We derive the properties of the massive cores from the spectral energy distributions of bright far-infrared point sources detected with the PACS instrument aboard Herschel. We report on the detection and characterization of pre- and protostellar cores in a massive filamentary infrared-dark cloud G011.11-0.12 using PACS. We characterize 18 cores directly associated with the filament, two of which have masses over 50 Msun, making them the best candidates to become massive stars in G011.11-0.12. These cores are likely at various stages of protostar formation, showing elevated temperature (<T> ~ 22 K) with respect to the ambient gas reservoir. The core masses (<M> ~ 24 Msun) are small compared to that in the cold filament. The mean core separation is 0.9 pc, well in excess of the Jeans length in the filament. We confirm that star formation in IRDCs is underway and diverse, and IRDCs have the capability of forming massive stars and clusters.

preprint2010arXiv

The visitor from an ancient galaxy: A planetary companion around an old, metal-poor red horizontal branch star

We report the detection of a planetary companion around HIP 13044, a metal-poor red horizontal branch star belonging to a stellar halo stream that results from the disruption of an ancient Milky Way satellite galaxy. The detection is based on radial velocity observations with FEROS at the 2.2-m MPG/ESO telescope. The periodic radial velocity variation of P=16.2 days can be distinguished from the periods of the stellar activity indicators. We computed a minimum planetary mass of 1.25 Jupiter masses and an orbital semimajor axis of 0.116 AU for the planet. This discovery is unique in three aspects: First, it is the first planet detection around a star with a metallicity much lower than few percent of the solar value; second, the planet host star resides in a stellar evolutionary stage that is still unexplored in the exoplanet surveys; third, the planetary system HIP 13044 most likely has an extragalactic origin in a disrupted former satellite of the Milky Way.

preprint2010arXiv

Trapping Solids at the Inner Edge of the Dead Zone: 3-D Global MHD Simulations

The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treating the turbulence driven by the magneto-rotational instability. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the `butterfly pattern' seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible.

preprint2010arXiv

Variable accretion as a mechanism for brightness variations in T Tau S

(Note: this is a shortened version of the original A&A-style structured abstract). The physical nature of the strong photometric variability of T Tau Sa, the more massive member of the Southern "infrared companion" to T Tau, has long been debated. Intrinsic luminosity variations due to variable accretion were originally proposed but later challenged in favor of apparent fluctuations due to time-variable foreground extinction. In this paper we use the timescale of the variability as a diagnostic for the underlying physical mechanism. Because the IR emission emerging from Sa is dominantly thermal emission from circumstellar dust at <=1500K, we can derive a minimum size of the region responsible for the time-variable emission. In the context of the variable foreground extinction scenario, this region must be (un-) covered within the variability timescale, which implies a minimum velocity for the obscuring foreground material. If this velocity supercedes the local Kepler velocity we can reject foreground extinction as a valid variability mechanism. The variable accretion scenario allows for shorter variability timescales since the variations in luminosity occur on much smaller scales, essentially at the surface of the star, and the disk surface can react almost instantly on the changing irradiation with a higher or lower dust temperature and according brightness. We have detected substantial variations at long wavelengths in T Tau S: +26% within four days at 12.8 micron. We show that this short-term variability cannot be due to variable extinction and instead must be due to variable accretion. Using a radiative transfer model of the Sa disk we show that variable accretion can in principle also account for the much larger (several magnitude) variations observed on timescales of several years. For the long-term variability, however, also variable foreground extinction is a viable mechanism.

preprint2009arXiv

Discovery of the Coldest Imaged Companion of a Sun-Like Star

We present the discovery of a brown dwarf or possible planet at a projected separation of 1.9" = 29 AU around the star GJ 758, placing it between the separations at which substellar companions are expected to form by core accretion (~5 AU) or direct gravitational collapse (typically >100 AU). The object was detected by direct imaging of its thermal glow with Subaru/HiCIAO. At 10-40 times the mass of Jupiter and a temperature of 550-640 K, GJ 758 B constitutes one of the few known T-type companions, and the coldest ever to be imaged in thermal light around a Sun-like star. Its orbit is likely eccentric and of a size comparable to Pluto's orbit, possibly as a result of gravitational scattering or outward migration. A candidate second companion is detected at 1.2" at one epoch.

preprint2009arXiv

Do we really know the dust? Systematics and uncertainties of the mid-infrared spectral analysis methods

The spectral region around 10 micrometer, showing prominent dust emission bands, is commonly used to derive the chemical composition of protoplanetary dust. Different analysis methods have been proposed for this purpose, but so far, no comparative test has been performed to test their validity. We calculated model spectra of disk models with different geometries and central sources, using a 2D radiative transfer code. These spectra were then fitted in a blind test using four different spectral decomposition methods. We studied the effect of disk structure (flared vs. flat), inclination angle, size of the inner disk hole and stellar luminosity on the fitted chemical composition. Our results show that the derived dust compositions by all methods deviate systematically from the real chemical composition. Out of the four tested spectral decomposition methods, our new two-layer temperature distribution method, differs the least from the input dust composition and the results show the weakest systematic effects. The reason for the deviations of the results given by other methods lies in their simplifying assumptions (e.g. single average grain temperature or one component continuum). We also tested the influence of different noise levels on the results of the spectral decomposition methods. We find that, for ground-based observations (8-13 micrometer) the expected uncertainty in the value of the crystallinity is about 11% for a signal-to-noise ratio of 100, while for space-based observations (7-17 micrometer) the uncertainty is about 5%. On the basis of our results, we propose a recipe for the analysis and interpretation of dust spectroscopy data in the mid-infrared which should be especially valuable for analysing Spitzer data and ground-based infrared spectroscopy data in the 10 micrometer window.

preprint2009arXiv

Dust amorphization in protoplanetary disks

High-energy irradiation of the circumstellar material might impact the structure and the composition of a protoplanetary disk and hence the process of planet formation. In this paper, we present a study on the possible influence of the stellar irradiation, indicated by X-ray emission, on the crystalline structure of the circumstellar dust. The dust crystallinity is measured for 42 class II T Tauri stars in the Taurus star-forming region using a decomposition fit of the 10 micron silicate feature, measured with the Spitzer IRS instrument. Since the sample includes objects with disks of various evolutionary stages, we further confine the target selection, using the age of the objects as a selection parameter. We correlate the X-ray luminosity and the X-ray hardness of the central object with the crystalline mass fraction of the circumstellar dust and find a significant anti-correlation for 20 objects within an age range of approx. 1 to 4.5 Myr. We postulate that X-rays represent the stellar activity and consequently the energetic ions of the stellar winds which interact with the circumstellar disk. We show that the fluxes around 1 AU and ion energies of the present solar wind are sufficient to amorphize the upper layer of dust grains very efficiently, leading to an observable reduction of the crystalline mass fraction of the circumstellar, sub-micron sized dust. This effect could also erase other relations between crystallinity and disk/star parameters such as age or spectral type.

preprint2009arXiv

Multiple low-turbulence starless cores associated with intermediate- to high-mass star formation

To characterize the initial conditions for intermediate- to high-mass star formation, we observed two Infrared Dark Clouds (IRDCs) that remain absorption features up to 70mum wavelength, with the PdBI in the 3.23mm dust continuum as well as the N2H+(1--0) and 13CS(2-1) line emission. While IRDC19175-4 is clearly detected in the 3.23mm continuum, the second source in the field, IRDC19175-5, is only barely observable above the 3sigma continuum detection threshold. However, the N2H+(1-0) observations reveal 17 separate sub-sources in the vicinity of the two IRDCs. Most of them exhibit low levels of turbulence (dv \leq 1km/s), indicating that the fragmentation process in these cores may be dominated by the interplay of thermal pressure and gravity, but not so much by turbulence. Combining the small line widths with the non-detection up to 70mum and the absence of other signs of star formation activity, most of these 17 cores with masses between sub-solar to ~10M_sun are likely still in a starless phase. Furthermore, we find a large CS depletion factor of the order 100. Although the strongest line and continuum peak is close to virial equilibrium, its slightly broader line width compared to the other cores is consistent with it being in a contraction phase potentially at the verge of star formation. The relative peak velocities between neighboring cores are usually below 1km/s, and we do not identify streaming motions along the filamentary structures. Average densities are between 10^5 and 10^6cm^{-3} (one to two orders of magnitude larger than for example in the Pipe nebula) implying relatively small Jeans-lengths that are consistent with the observed core separations of the order 5000AU. The quest for high-mass starless cores prior to any star formation activity remains open.

preprint2009arXiv

Search for very low-mass brown dwarfs and free-floating planetary-mass objects in Taurus

The number of low-mass brown dwarfs and even free floating planetary mass objects in young nearby star-forming regions and associations is continuously increasing, offering the possibility to study the low-mass end of the IMF in greater detail. In this paper, we present six new candidates for (very) low-mass objects in the Taurus star-forming region one of which was recently discovered in parallel by Luhman et al. (2009). The underlying data we use is part of a new database from a deep near-infrared survey at the Calar Alto observatory. The survey is more than four magnitudes deeper than the 2MASS survey and covers currently ~1.5 square degree. Complementary optical photometry from SDSS were available for roughly 1.0 square degree. After selection of the candidates using different color indices, additional photometry from Spitzer/IRAC was included in the analysis. In greater detail we focus on two very faint objects for which we obtained J-band spectra. Based on comparison with reference spectra we derive a spectral type of L2+/-0.5 for one object, making it the object with the latest spectral type in Taurus known today. From models we find the effective temperature to be 2080+/-140 K and the mass 5-15 Jupiter masses. For the second source the J-band spectrum does not provide a definite proof of the young, low-mass nature of the object as the expected steep water vapor absorption at 1.33 micron is not present in the data. We discuss the probability that this object might be a background giant or carbon star. If it were a young Taurus member, however, a comparison to theoretical models suggests that it lies close to or even below the deuterium burning limit (<13 Jupiter masses) as well. A first proper motion analysis for both objects shows that they are good candidates for being Taurus members.

preprint2009arXiv

The far future of exoplanet direct characterization

In this outlook we describe what could be the next steps of the direct characterization of habitable exoplanets after first the medium and large mission projects and investigate the benefits of the spectroscopic and direct imaging approaches. We show that after third and fourth generation missions foreseeable for the next 100 years, we will face a very long era before being able to see directly the morphology of extrasolar organisms.

preprint2008arXiv

AstraLux - the Calar Alto 2.2-m telescope Lucky Imaging camera

AstraLux is a Lucky Imaging camera for the Calar Alto 2.2-m telescope, based on an electron-multiplying high speed CCD. By selecting only the best 1-10% of several thousand short exposure frames, AstraLux provides nearly diffraction limited imaging capabilities in the SDSS i' and z' filters over a field of view of 24x24 arcseconds. By choosing commercially available components wherever possible, the instrument could be built in short time and at comparably low cost. We briefly present the instrument design, the data reduction pipeline, and summarise the performance and characteristics

preprint2008arXiv

IRAM-PdBI Observations of Binary Protostars I: The Hierarchical System SVS13 in NGC1333

We present millimeter interferometric observations of the young stellar object SVS13 in NCG1333 in the N2H+(1-0) line and at 1.4 and 3mm dust continuum, using the IRAM Plateau de Bure interferometer. The results are complemented by infrared data from the Spitzer Space Telescope. The millimeter dust continuum images resolve four sources (A, B, C, and VLA3) in SVS13. With the dust continuum images, we derive gas masses of 0.2-1.1 M_sun for the sources. N2H+(1-0) line emission is detected and spatially associated with the dust continuum sources B and VLA3. The observed mean line width is ~0.48 km/s and the estimated virial mass is ~0.7 M_sun. By simultaneously fitting the seven hyperfine line components of N2H+, we derive the velocity field and find a symmetric velocity gradient of about 28 km/s/pc across sources B and VLA3, which could be explained by core rotation. The velocity field suggests that sources B and VLA3 are forming a physically bound protobinary system embedded in a common N2H+ core. Spitzer images show mid-infrared emission from sources A and C, which is spatially associated with the mm dust continuum emission. No infrared emission is detected from source B, implying that the source is deeply embedded. Based on the morphologies and velocity structure, we propose a hierarchical fragmentation picture for SVS13 where the three sources (A, B, and C) were formed by initial fragmentation of a filamentary prestellar core, while the protobinary system (sources B and VLA3) was formed by rotational fragmentation of a single collapsing sub-core.

preprint2008arXiv

Star-forming cores embedded in a massive cold clump: Fragmentation, collapse and energetic outflows

The fate of massive cold clumps, their internal structure and collapse need to be characterised to understand the initial conditions for the formation of high-mass stars, stellar systems, and the origin of associations and clusters. We explore the onset of star formation in the 75 M_sun SMM1 clump in the region ISOSS J18364-0221 using infrared and (sub-)millimetre observations including interferometry. This contracting clump has fragmented into two compact cores SMM1 North and South of 0.05 pc radius, having masses of 15 and 10 M_sun, and luminosities of 20 and 180 L_sun. SMM1 South harbours a source traced at 24 and 70um, drives an energetic molecular outflow, and appears supersonically turbulent at the core centre. SMM1 North has no infrared counterparts and shows lower levels of turbulence, but also drives an outflow. Both outflows appear collimated and parsec-scale near-infrared features probably trace the outflow-powering jets. We derived mass outflow rates of at least 4E-5 M_sun/yr and outflow timescales of less than 1E4 yr. Our HCN(1-0) modelling for SMM1 South yielded an infall velocity of 0.14 km/s and an estimated mass infall rate of 3E-5 M_sun/yr. Both cores may harbour seeds of intermediate- or high-mass stars. We compare the derived core properties with recent simulations of massive core collapse. They are consistent with the very early stages dominated by accretion luminosity.

preprint2008arXiv

The Initial Mass Function of the Stellar Association NGC 602 in the Small Magellanic Cloud with Hubble Space Telescope ACS Observations

We present our photometric study of the stellar association NGC 602 in the wing of the Small Magellanic Cloud (SMC). The data were taken in the filters F555W and F814W using the Advanced Camera for Surveys (ACS) on-board the Hubble Space Telescope (HST). Photometry was performed using the ACS module of the stellar photometry package DOLPHOT. We detected more than 5,500 stars with a magnitude range of 14 \lsim m_{555} \lsim 28 mag. Three prominent stellar concentrations are identified with star counts in the observed field, the association NGC 602 itself, and two clusters, one of them not being currently in any known catalog. The Color-Magnitude Diagrams (CMDs) of both clusters show features typical for young open clusters, while that of the association reveals bright main sequence (MS) and faint pre-main sequence (PMS) stars as the members of the system. We construct the initial mass spectrum (IMS) of the association by applying an age-independent method of counting the PMS stars within evolutionary tracks, while for the bright MS stars we transform their magnitudes to masses with the use of mass-luminosity relations. The IMS of NGC 602 is found to be well represented by a single-power law, corresponding to an Initial Mass Function (IMF) of slope Γ\approx -1.2 for 1 \lsim M/M{\solar} \lsim 45. This indicates that the shape of the IMF of a star forming system in the SMC for stars with masses higher than 1 M{\solar} seems to be quite similar to the field IMF in the solar neighborhood.

preprint2007arXiv

Clustered Star Formation in the Small Magellanic Cloud. A Spitzer/IRAC View of the Star-Forming Region NGC 602/N 90

We present Spitzer/IRAC photometry on the star-forming HII region N 90, related to the young stellar association NGC 602 in the Small Magellanic Cloud. Our photometry revealed bright mid-infrared sources, which we classify with the use of a scheme based on templates and models of red sources in the Milky Way, and criteria recently developed from the Spitzer Survey of the SMC for the selection of candidate Young Stellar Objects (YSOs). We detected 57 sources in all four IRAC channels in a 6.2' x 4.8' field-of-view centered on N 90; 22 of these sources are classified as candidate YSOs. We compare the locations of these objects with the position of optical sources recently found in the same region with high-resolution HST/ACS imaging of NGC 602, and we find that 17 candidate YSOs have one or more optical counterparts. All of these optical sources are identified as pre-main sequence stars, indicating, thus, ongoing clustered star formation events in the region. The positions of the detected YSOs and their related PMS clusters give a clear picture of the current star formation in N 90, according to which the young stellar association photo-ionizes the surrounding interstellar medium, revealing the HII nebula, and triggering sequential star formation events mainly along the eastern and southern rims of the formed cavity of the parental molecular cloud.

preprint2007arXiv

Rapid planetesimal formation in turbulent circumstellar discs

The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.

preprint2007arXiv

Silicate Dust in Evolved Protoplanetary Disks: Growth, Sedimentation, and Accretion

We present the Spitzer IRS spectra for 33 young stars in Tr 37 and NGC 7160. The sample includes the high- and intermediate-mass stars with MIPS 24 microns excess, the only known active accretor in the 12 Myr-old cluster NGC 7160, and 19 low-mass stars with disks in the 4 Myr-old cluster Tr 37. We examine the 10 microns silicate feature, present in the whole sample of low-mass star and in 3 of the high- and intermediate-mass targets, and we find that PAH emission is detectable only in the Herbig Be star. We analyze the composition and size of the warm photospheric silicate grains by fitting the 10 microns silicate feature, and study the possible correlations between the silicate characteristics and the stellar and disk properties (age, SED slope, accretion rate, spectral type). We find indications of dust settling with age and of the effect of turbulent enrichment of the disk atmosphere with large grains. Crystalline grains are only small contributors to the total silicate mass in all disks, and do not seem to correlate with any other property, except maybe binarity. We also observe that spectra with very weak silicate emission are at least 3 times more frequent among M stars than among earlier spectral types, which may be an evidence of inner disk evolution. Finally, we find that 5 of the high- and intermediate-mass stars have SEDs and IRS spectra consistent with debris disk models involving planet formation, which could indicate debris disk formation at ages as early as 4 Myr.

preprint2006arXiv

Are Debris Disks and Massive Planets Correlated?

Using data from the Spitzer Space Telescope Legacy Science Program ``Formation and Evolution of Planetary Systems'' (FEPS), we have searched for debris disks around 9 FGK stars (2-10 Gyr), known from radial velocity (RV) studies to have one or more massive planets. Only one of the sources, HD 38529, has excess emission above the stellar photosphere; at 70 micron the signal-to-noise ratio in the excess is 4.7 while at wavelengths < 30 micron there is no evidence of excess. The remaining sources show no excesses at any Spitzer wavelengths. Applying survival tests to the FEPS sample and the results for the FGK survey published in Bryden et al. (2006), we do not find a significant correlation between the frequency and properties of debris disks and the presence of close-in planets. We discuss possible reasons for the lack of a correlation.

preprint2006arXiv

Chemistry in disks I - Deep search for N$_2$H$^+$ in the protoplanetary disks around LkCa 15, MWC 480, and DM Tau

Aims: To constrain the ionization fraction in protoplanetary disks, we present new high-sensitivity interferometric observations of N$_2$H$^+$ in three disks surrounding DM Tau, LkCa 15, and MWC 480. Methods: We used the IRAM PdBI array to observe the N$_2$H$^+$ J=1-0 line and applied a $χ^2$-minimization technique to estimate corresponding column densities. These values are compared, together with HCO$^+$ column densities, to results of a steady-state disk model with a vertical temperature gradient coupled to gas-grain chemistry. Results: We report two \ndhp detections for LkCa 15 and DM Tau at the $5 σ$ level and an upper limit for MWC 480. The column density derived from the data for LkCa 15 is much lower than previously reported. The [N$_2$H$^+$/HCO$^+$] ratio is on the order of 0.02--0.03. So far, HCO$^+$ remains the most abundant observed molecular ion in disks. Conclusions: All the observed values generally agree with the modelled column densities of disks at an evolutionary stage of a few million years (within the uncertainty limits), but the radial distribution of the molecules is not reproduced well. The low inferred concentration of N$_2$H$^+$ in three disks around low-mass and intermediate-mass young stars implies that this ion is not a sensitive tracer of the overall disk ionization fraction.

preprint2006arXiv

Direct imaging of the young spectroscopic binary HD 160934

We report on the direct detection of a close companion to HD 160934, a young active star, SB1 spectroscopic binary, and suggested member of the AB Doradus moving group. High angular resolution at the Calar Alto 2.2m telescope was achieved by means of the Lucky Imaging technique, allowing direct imaging close to the diffraction limit in the SDSS z' band. Our results are combined with pre-discovery HST archive data, own UBVRI broadband photometry, published JHK magnitudes, and available radial velocity measurements to constrain the physical properties of the HD 160934 close binary. We suggest that the direct detection may be identical to the spectroscopically discovered companion.

preprint2004arXiv

Direct Detection of Sub-stellar Companions with MIDI

Current detection methods of planetary companions do not allow retrieving their spectral information properties. The method of Differential Interferometry has the potential to complement such information by comparing interferometric observables in various spectral channels. We outline the basic aspects of this method and how it can soon be realised by the Mid-infrared Interferometric instrument (MIDI) at the Very Large Telescope Interferometer (VLTI) observatory. A set of possible candidates for such direct observation has been selected among currently identified planetary companions. Differential Interferometry with MIDI would be complementary to other ground-based programs. The method may also be an alternative for future space-based planetary spectroscopy.

preprint2003arXiv

Orbital Migration and Mass Accretion of Protoplanets in 3D Global Computations with Nested Grids

We investigate the evolution of protoplanets with different masses embedded in an accretion disk, via global fully three-dimensional hydrodynamical simulations. We consider a range of planetary masses extending from one and a half Earth's masses up to one Jupiter's mass, and we take into account physically realistic gravitational potentials of forming planets. In order to calculate accurately the gravitational torques exerted by disk material and to investigate the accretion process onto the planet, the flow dynamics has to be thoroughly resolved on long as well as short length scales. We achieve this strict resolution requirement by applying a nested-grid refinement technique which allows to greatly enhance the local resolution. Our results from altogether 51 simulations show that for large planetary masses, approximately above a tenth of the Jupiter's mass, migration rates are relatively constant, as expected in type II migration regime and in good agreement with previous two-dimensional calculations. In a range between seven and fifteen Earth's masses, we find a dependency of the migration speed on the planetary mass that yields time scales considerably longer than those predicted by linear analytical theories. This property may be important in determining the overall orbital evolution of protoplanets. The growth time scale is minimum around twenty Earth-masses, but it rapidly increases for both smaller and larger mass values. Significant differences between two- and three-dimensional calculations are found in particular for objects with masses smaller than ten Earth-masses. We also derive an analytical approximation for the numerically computed mass growth rates.

preprint2003arXiv

Thermo-Hydrodynamics of Circumstellar Disks with High-mass Planets

With a series of numerical simulations, we analyze the thermo-hydrodynamical evolution of circumstellar disks containing Jupiter-size protoplanets. In the framework of the two-dimensional approximation, we consider an energy equation that includes viscous heating and radiative effects in a simplified, yet consistent form. Multiple nested grids are used in order to study both global and local features around the planet. By means of different viscosity prescriptions, we investigate various temperature regimes. A planetary mass range from 0.1 to 1 Mj is examined. Computations show that gap formation is a general property which affects density, pressure, temperature, optical thickness, and radiated flux distributions. However, it remains a prominent feature only when the kinematic viscosity is on the order of 10^(15) cm^2/s or lower. Around accreting planets, a circumplanetary disk forms that has a surface density profile decaying exponentially with the distance and whose mass is 5-6 orders of magnitudes smaller than Jupiter's mass. Circumplanetary disk temperature profiles decline roughly as the inverse of the distance from the planet. Temperatures range from some 10 to ~1000 K. Planetary accretion and migration rates depend on the viscosity regime, with discrepancies within an order of magnitude. Estimates of growth and migration time scales inferred by these models are on the same orders of magnitude as those previously obtained with locally isothermal simulations.

preprint2001arXiv

Nested-grid calculations of disk-planet interaction

We study the evolution of embedded protoplanets in a protostellar disk using very high resolution nested-grid computations. This method allows us to perform global simulations of planets orbiting in disks and, at the same time, to resolve in detail the dynamics of the flow inside the Roche lobe of the planet. The primary interest of this work lies in the analysis of the gravitational torque balance acting on the planet. For this purpose we study planets of different masses, ranging from one Earth-mass up to one Jupiter-mass, assuming typical parameters of the protostellar disk. The high resolution of the method allows a precise determination of the mass flow onto the planet and the resulting torques. The obtained migration time scales are in the range from few times 10^4 years, for intermediate mass planets, to 10^6 years, for very low and high mass planets. Typical growth time scales depend strongly on the planetary mass, ranging from a few hundred years, in the case of Earth-type planets, to several ten thousand years, in the case of Jupiter-type planets.