sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Güdel, Manuel

Normalized to: Güdel, M.

78 article(s) in total. 680 co-authors, from 1 to 20 common article(s). Median position in authors list is 5,0.

[1]  [pdf] - 2104969
Feedback of molecular outflows from protostars in NGC 1333, revealed by Herschel and Spitzer spectro-imaging observations
Comments: 19 pages, 14 figures, 5 tables, A&A accepted
Submitted: 2020-06-01
Large scale spectral maps of star forming regions enable the comparative study of the gas excitation around an ensemble of sources at a common frame of reference, providing direct insights in the multitude of processes involved. In this paper we employ spectral-line maps to decipher the excitation, the kinematical and dynamical processes in NGC 1333 as revealed by a number of different emission lines, aiming to set a reference for the applicability of tracers in constraining diverse physical processes. We reconstruct line maps for H$_2$ , CO, H$_2$O and C$^+$ using data obtained with the Spitzer-IRS and Herschel HIFI-SPIRE. We compare the morphological features of the maps and derive the gas excitation for regions of interest employing LTE and non-LTE methods. We also calculate the kinematical and dynamical properties for each outflow tracer consistently for all outflows in NGC 1333. We finally measure the water abundance in outflows with respect to carbon monoxide and molecular hydrogen. CO and H$_2$ are highly excited around B-stars and at lower levels trace protostellar outflows. H$_2$O emission is dominated by a moderately fast component associated with outflows. Intermediate J CO lines appear brightest at the locations traced by a narrow H$_2$O component, indicating that beyond the dominating collisional processes, a secondary, radiative excitation component can also be active. The morphology, kinematics, excitation and abundance variations of water are consistent with its excitation and partial dissociation in shocks. Water abundance ranges between 5 x 10$^{-7}$ and 10$^{-5}$, with the lower values being more representative. Water is brightest and most abundant around IRAS 4A which is consistent with the latter hosting a hot corino source. Finally, the outflow mass flux is found highest for CO and decreases by one and two orders of magnitude for H$_2$ and H$_2$O, respectively.
[2]  [pdf] - 2076955
The solar wind from a stellar perspective: how do low-resolution data impact the determination of wind properties?
Comments: accepted for publication in A&A
Submitted: 2020-02-24
Alfv\'en-wave-driven 3D magnetohydrodynamic (MHD) models, which are increasingly used to predict stellar wind properties, contain unconstrained parameters and rely on low-resolution stellar magnetograms. We explore the effects of the input Alfv\'en wave energy flux and the surface magnetogram on the wind properties predicted by the Alfv\'en Wave Solar Model (AWSoM). We lowered the resolution of two solar magnetograms during solar cycle maximum and minimum using spherical harmonic decomposition. The Alfv\'en wave energy was altered based on non-thermal velocities determined from a far ultraviolet (FUV) spectrum of the solar twin 18 Sco. Additionally, low-resolution magnetograms of three solar analogues were obtained using Zeeman Doppler imaging (ZDI). Finally, the simulated wind properties were compared to Advanced Composition Explorer (ACE) observations. AWSoM simulations using well constrained input parameters taken from solar observations can reproduce the observed solar wind mass and angular momentum loss rates. The resolution of the magnetogram has a small impact on the wind properties and only during cycle maximum. However, variation in Alfv\'en wave energy influences the wind properties irrespective of the solar cycle activity level. Furthermore, solar wind simulations carried out using the low-resolution magnetogram of the three stars instead of the solar magnetogram could lead to an order of a magnitude difference in the simulated wind properties. The choice in Alfv\'en energy has a stronger influence on the wind output compared to the magnetogram resolution. The influence could be even stronger for stars whose input boundary conditions are not as well constrained as those of the Sun. Unsurprisingly, replacing the solar magnetogram with a stellar magnetogram could lead to completely inaccurate solar wind properties, and should be avoided in solar and stellar wind simulations.
[3]  [pdf] - 1995544
Escape and evolution of Mars' CO2 atmosphere: Influence of suprathermal atoms
Submitted: 2019-11-07
With a Monte-Carlo model we investigate the escape of hot oxygen and carbon from the martian atmosphere for four points in time in its history corresponding to 1, 3, 10, and 20 times the present solar EUV flux. We study and discuss different sources of hot oxygen and carbon atoms in the thermosphere and their changing importance with the EUV flux. The increase of the production rates due to higher densities resulting from the higher EUV flux competes against the expansion of the thermosphere and corresponding increase in collisions. We find that the escape due to photodissociation increases with increasing EUV level. However, for the escape via some other reactions, e.g.~dissociative recombination of O$_2^+$, this is only true until the EUV level reaches 10 times the present EUV flux, and then the rates start to decrease. Furthermore, our results show that Mars could not have had a dense atmosphere at the end of the Noachian epoch, since such an atmosphere would not have been able to escape until today. In the pre-Noachian era, most of a magma ocean and volcanic activity related outgassed CO$_2$ atmosphere could have been lost thermally until the Noachian epoch, when non-thermal loss processes such as suprathermal atom escape became dominant. Thus, early Mars could have been hot and wet during the pre-Noachian era with surface CO$_2$ pressures larger than 1 bar during the first 300 Myr after the planet's origin.
[4]  [pdf] - 1994441
Solar XUV and ENA-driven water loss from early Venus' steam atmosphere
Comments: 6 figures and 2 tables
Submitted: 2019-11-06
The influence of the hydrogen hydrodynamic upper atmosphere escape, driven by the solar soft X-ray and extreme ultraviolet radiation (XUV) flux, on an expected magma ocean outgassed steam atmosphere of early Venus is studied. By assuming that the young Sun was either a weak or moderate active young G star, we estimated the water loss from a hydrogen dominated thermosphere due to the absorption of the solar XUV flux and the precipitation of solar wind produced energetic hydrogen atoms (ENAs). The production of ENAs and their interaction with the hydrodynamic extended upper atmosphere, including collision-related feedback processes, have been calculated by means of Monte Carlo models. ENAs that collide in the upper atmosphere deposit their energy and heat the surrounding gas mainly above the main XUV energy deposition layer. It is shown that precipitating ENAs modify the thermal structure of the upper atmosphere, but the enhancement of the thermal escape rates caused by these energetic hydrogen atoms is negligible. Our results also indicate that the majority of oxygen arising from dissociated H$_2$O molecules is left behind during the first 100 Myr. It is thus suggested that the main part of the remaining oxygen has been absorbed by crustal oxidation.
[5]  [pdf] - 1979713
Infrared and sub-mm observations of outbursting young stars with Herschel and Spitzer
Comments: 23 pages, 17 figures, A&A accepted
Submitted: 2019-09-02
Episodic accretion plays an important role in the evolution of young stars. Although it has been under investigation for a long time, the origin of such episodic accretion events is not yet understood. We investigate the dust and gas emission of a sample of young outbursting sources in the infrared to get a better understanding of their properties and circumstellar material, and we use the results in a further work to model the objects. We used Herschel data, from our PI program of 12 objects and complemented with archival observations to obtain the spectral energy distributions (SEDs) and spectra of our targets. We report here the main characteristics of our sample, focussing on the SED properties and on the gas emission lines detected in the PACS and SPIRE spectra. The SEDs of our sample show the diversity of the outbursting sources, with several targets showing strong emission in the far-infrared from the embedded objects. Most of our targets reside in a complex environment, which we discuss in detail. We detected several atomic and molecular lines, in particular rotational CO emission from several transitions from J=38-37 to J=4-3. We constructed rotational diagrams for the CO lines, and derived in three domains of assumed local thermodynamic equilibrium (LTE) temperatures and column densities, ranging mainly between 0-100 K and 400-500K. We confirm correlation in our sample between intense CO $J=16-15$ emission and the column density of the warm domain of CO, N(warm). We notice a strong increase in luminosity of HH 381 IRS and a weaker increase for PP 13 S, which shows the beginning of an outburst.
[6]  [pdf] - 1931032
Modelling atmospheric escape and MgII near-ultraviolet absorption of the highly irradiated hot Jupiter WASP-12b
Comments: 23 pages, 8 figures
Submitted: 2019-08-07
We present two-dimensional multi-fluid numerical modelling of the upper atmosphere of the hot Jupiter WASP-12b. The model includes hydrogen chemistry, and self-consistently describes the expansion of the planetary upper atmosphere and mass loss due to intensive stellar irradiation, assuming a weakly magnetized planet. We simulate the planetary upper atmosphere and its interaction with the stellar wind (SW) with and without the inclusion of tidal force and consider different XUV irradiation conditions and SW parameters. With the inclusion of tidal force, even for a fast SW, the escaping planetary material forms two streams, propagating towards and away from the star. The atmospheric escape and related mass loss rate reaching the value of 10^12 gs^-1 appear to be mostly controlled by the stellar gravitational pull. We computed the column density and dynamics of MgII ions considering three different sets of SW parameters and XUV fluxes. The simulations enable to compute the absorption at the position of the Mg h line and to reproduce the times of ingress and egress. In case of a slow SW and without accounting for tidal force, the high orbital velocity leads to the formation of a shock approximately in the direction of the planetary orbital motion. In this case, mass loss is proportional to the stellar XUV flux. At the same time, ignoring of tidal effects for WASP-12b is a strong simplification, so the scenario with a shock, altogether is an unrealistic one.
[7]  [pdf] - 1914753
Detecting volcanically produced tori along orbits of exoplanets using UV spectroscopy
Comments: White paper submitted on the use of HST for TESS follow-up observations, 5 pages, 2 figures
Submitted: 2019-07-11
We suggest to use the Hubble Space Telescople (HST) follow-up observations of the TESS targets for detecting possible plasma tori along the orbits of exoplanets orbiting M dwarfs. The source of the torus could be planetary volcanic activity due to tidal or electromagnetic induction heating. Fast losses to space for planets orbiting these active stars can lead to the lost material forming a torus along the planetary orbit, similar to the Io plasma torus. We show that such torus would be potentially detectable by the HST in the UV.
[8]  [pdf] - 1898019
Stellar activity and planetary atmosphere evolution in tight binary star systems
Comments: Accepted for publication by A&A
Submitted: 2019-04-17
Context. In tight binary star systems, tidal interactions can significantly influence the rotational and orbital evolution of both stars, and therefore their activity evolution. This can have strong effects on the atmospheric evolution of planets that are orbiting the two stars. Aims. In this paper, we aim to study the evolution of stellar rotation and of X-ray and ultraviolet (XUV) radiation in tight binary systems consisting of two solar mass stars and use our results to study planetary atmosphere evolution in the habitable zones of these systems. Methods. We have applied a rotation model developed for single stars to binary systems, taking into account the effects of tidal interactions on the rotational and orbital evolution of both stars. We used empirical rotation-activity relations to predict XUV evolution tracks for the stars, which we used to model hydrodynamic escape of hydrogen dominated atmospheres. Results. When significant, tidal interactions increase the total amount of XUV energy emitted, and in the most extreme cases by up to factor of $\sim$50. We find that in the systems that we study, habitable zone planets with masses of 1~M$_\oplus$ can lose huge hydrogen atmospheres due to the extended high levels of XUV emission, and the time that is needed to lose these atmospheres depends on the binary orbital separation.For some orbital separations, and when the stars are born as rapid rotators, it is also possible for tidal interactions to protect atmospheres from erosion by quickly spinning down the stars. For very small orbital separations, the loss of orbital angular momentum by stellar winds causes the two stars to merge. We suggest that the merging of the two stars could cause previously frozen planets to become habitable due to the habitable zone boundaries moving outwards.
[9]  [pdf] - 1875445
Extreme hydrodynamic losses of Earth-like atmospheres in the habitable zones of very active stars
Comments: Accepted for publication by A&A Letters
Submitted: 2019-04-01
Aims. In this letter, we calculate for the first time the full transonic hydrodynamic escape of an Earth-like atmosphere. We consider the case of an Earth-mass planet with an atmospheric composition identical to that of the current Earth orbiting at 1 AU around a young and very active solar mass star. Methods. To model the upper atmosphere, we used the Kompot Code, which is a first-principles model that calculates the physical structures of the upper atmospheres of planets, taking into account hydrodynamics and the main chemical and thermal processes taking place in the upper atmosphere of a planet. This model enabled us to calculate the 1D vertical structure of the atmosphere using as input the high-energy spectrum of a young and active Sun. Results. The atmosphere has the form of a transonic hydrodynamic Parker wind, which has an outflow velocity at the upper boundary of our computational domain that exceeds the escape velocity. The outflowing gas is dominated by atomic nitrogen and oxygen and their ion equivalents and has a maximum ionization fraction of 20%. The mass outflow rate is found to be 1.8x10^9 g s^-1, which would erode the modern Earth's atmosphere in less than 0.1 Myr. Conclusions. This extreme mass loss rate suggests that an Earth-like atmosphere cannot form when the planet is orbiting within the habitable zone of a very active star. Instead, such an atmosphere can only form after the activity of the star has decreased to a much lower level. This happened in the early atmosphere of the Earth, which was likely dominated by other gases such as CO2. Since the time it takes for the activity of a star to decay is highly dependent on its mass, this is important for understanding possible formation timescales for planets orbiting low-mass stars.
[10]  [pdf] - 1854239
The fastest components in stellar jets
Comments: Submitted to the 2020 Astronomy decadal survey as science white paper
Submitted: 2019-03-22
Young stars accrete mass from a circumstellar disk, but at the same time disk and star eject outflows and jets. These outflows have an onion-like structure where the innermost and fastest layers are surrounded by increasingly lower velocity components. The outer layers are probably photo-evaporative and magnetocentrifugally launched disk winds, but the nature of the inner winds is still uncertain. Since the fastest components carry only a small fraction of the mass, they are best observed at high-energies (X-ray and UV) as the slower, more massive components do not reach plasma temperatures sufficient for relevant X-ray or UV emission. Outflows are the most likely way in which a star or its disk can shed angular momentum and allow accretion to proceed; thus we cannot understand the accretion and the rotation rate of young stars if we cannot solve the origin of the inner jet components. Stellar jets share characteristics with their counterparts in more massive astrophysical objects, such as stellar mass black holes and AGN, with the added benefit that young stars are found at much closer distances and thus scales not accessible in other types of objects can be resolved. To understand the origin and impact of the inner jets, sub-arcsecond imaging and spectroscopy in the UV and X-rays is required, together with theory and modelling to interpret existing and future observations.
[11]  [pdf] - 1855786
Transit Ly-$\alpha$ signatures of terrestrial planets in the habitable zones of M dwarfs
Comments: 17 pages, 12 figures, accepted for publication in Astronomy & Astrophysics
Submitted: 2019-03-06
We modeled the transit signatures in the Lya line of a putative Earth-sized planet orbiting in the HZ of the M dwarf GJ436. We estimated the transit depth in the Lya line for an exo-Earth with three types of atmospheres: a hydrogen-dominated atmosphere, a nitrogen-dominated atmosphere, and a nitrogen-dominated atmosphere with an amount of hydrogen equal to that of the Earth. We calculated the in-transit absorption they would produce in the Lya line. We applied it to the out-of-transit Lya observations of GJ 436 obtained by the HST and compared the calculated in-transit absorption with observational uncertainties to determine if it would be detectable. To validate the model, we also used our method to simulate the deep absorption signature observed during the transit of GJ 436b and showed that our model is capable of reproducing the observations. We used a DSMC code to model the planetary exospheres. The code includes several species and traces neutral particles and ions. At the lower boundary of the DSMC model we assumed an atmosphere density, temperature, and velocity obtained with a hydrodynamic model for the lower atmosphere. We showed that for a small rocky Earth-like planet orbiting in the HZ of GJ436 only the hydrogen-dominated atmosphere is marginally detectable with the STIS/HST. Neither a pure nitrogen atmosphere nor a nitrogen-dominated atmosphere with an Earth-like hydrogen concentration in the upper atmosphere are detectable. We also showed that the Lya observations of GJ436b can be reproduced reasonably well assuming a hydrogen-dominated atmosphere, both in the blue and red wings of the Lya line, which indicates that warm Neptune-like planets are a suitable target for Lya observations. Terrestrial planets can be observed in the Lya line if they orbit very nearby stars, or if several observational visits are available.
[12]  [pdf] - 2034325
Constraining the stellar energetic particle flux in young solar-like stars
Comments: 2 pages, 1 figure, accepted for IAUS345 "Origins: From the Protosun to the First Steps of Life" proceedings
Submitted: 2019-02-03
Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a large number of energetic particles ($E\gtrsim10\,$MeV), often called solar cosmic rays. The particle flux of the contemporary Sun cannot explain these anomalies, but, similar to \mbox{T Tauri} stars, the young Sun was more active and probably produced enough high energy particles. However, the stellar particle (SP) flux of young stars is essentially unknown. We model the impact of high-energy ionization sources on the chemistry of the circumstellar environment (disks and envelopes). The model includes X-ray radiative transfer and makes use of particle transport models to calculate the individual molecular hydrogen ionization rates. We study the impact on the chemistry via the ionization tracers HCO$^+$ and N$_2$H$^+$. We argue that spatially resolved observations of those molecules combined with detailed models allow for disentangling the contribution of the individual high-energy ionization sources and to put constraints on the SP flux in young stars.
[13]  [pdf] - 1811290
Episodic excursions of low-mass protostars on the Hertzsprung-Russell diagram
Comments: 15 pages, 12 figures
Submitted: 2018-12-28
Following our recent work devoted to the effect of accretion on the pre-main-sequence evolution of low-mass stars, we perform a detailed analysis of episodic excursions of low-mass protostars in the Hertzsprung-Russell (H-R) diagram triggered by strong mass accretion bursts typical of FU Orionis-type objects (FUors). These excursions reveal themselves as sharp increases in the stellar total luminosity and/or effective temperature of the protostar and can last from hundreds to a few thousands of years, depending on the burst strength and characteristics of the protostar. During the excursions, low-mass protostars occupy the same part of the H-R diagram as young intermediate-mass protostars in the quiescent phase of accretion. Moreover, the time spent by low-mass protostars in these regions is on average a factor of several longer than that spent by the intermediate-mass stars in quiescence. During the excursions, low-mass protostars pass close to the position of most known FUors in the H-R diagram, but owing to intrinsic ambiguity the model stellar evolutionary tracks are unreliable in determining the FUor properties. We find that the photospheric luminosity in the outburst state may dominate the accretion luminosity already after a few years after the onset of the outburst, meaning that the mass accretion rates of known FUors inferred from the bolometric luminosity may be systematically overestimated, especially in the fading phase.
[14]  [pdf] - 1795998
Observatory science with eXTP
Zand, Jean J. M. in 't; Bozzo, Enrico; Qu, Jinlu; Li, Xiang-Dong; Amati, Lorenzo; Chen, Yang; Donnarumma, Immacolata; Doroshenko, Victor; Drake, Stephen A.; Hernanz, Margarita; Jenke, Peter A.; Maccarone, Thomas J.; Mahmoodifar, Simin; de Martino, Domitilla; De Rosa, Alessandra; Rossi, Elena M.; Rowlinson, Antonia; Sala, Gloria; Stratta, Giulia; Tauris, Thomas M.; Wilms, Joern; Wu, Xuefeng; Zhou, Ping; Agudo, Iván; Altamirano, Diego; Atteia, Jean-Luc; Andersson, Nils A.; Baglio, M. Cristina; Ballantyne, David R.; Baykal, Altan; Behar, Ehud; Belloni, Tomaso; Bhattacharyya, Sudip; Bianchi, Stefano; Bilous, Anna; Blay, Pere; Braga, João; Brandt, Søren; Brown, Edward F.; Bucciantini, Niccolò; Burderi, Luciano; Cackett, Edward M.; Campana, Riccardo; Campana, Sergio; Casella, Piergiorgio; Cavecchi, Yuri; Chambers, Frank; Chen, Liang; Chen, Yu-Peng; Chenevez, Jérôme; Chernyakova, Maria; Chichuan, Jin; Ciolfi, Riccardo; Costantini, Elisa; Cumming, Andrew; D'Aı, Antonino; Dai, Zi-Gao; D'Ammando, Filippo; De Pasquale, Massimiliano; Degenaar, Nathalie; Del Santo, Melania; D'Elia, Valerio; Di Salvo, Tiziana; Doyle, Gerry; Falanga, Maurizio; Fan, Xilong; Ferdman, Robert D.; Feroci, Marco; Fraschetti, Federico; Galloway, Duncan K.; Gambino, Angelo F.; Gandhi, Poshak; Ge, Mingyu; Gendre, Bruce; Gill, Ramandeep; Götz, Diego; Gouiffès, Christian; Grandi, Paola; Granot, Jonathan; Güdel, Manuel; Heger, Alexander; Heinke, Craig O.; Homan, Jeroen; Iaria, Rosario; Iwasawa, Kazushi; Izzo, Luca; Ji, Long; Jonker, Peter G.; José, Jordi; Kaastra, Jelle S.; Kalemci, Emrah; Kargaltsev, Oleg; Kawai, Nobuyuki; Keek, Laurens; Komossa, Stefanie; Kreykenbohm, Ingo; Kuiper, Lucien; Kunneriath, Devaky; Li, Gang; Liang, En-Wei; Linares, Manuel; Longo, Francesco; Lu, Fangjun; Lutovinov, Alexander A.; Malyshev, Denys; Malzac, Julien; Manousakis, Antonios; McHardy, Ian; Mehdipour, Missagh; Men, Yunpeng; Méndez, Mariano; Mignani, Roberto P.; Mikusincova, Romana; Miller, M. Coleman; Miniutti, Giovanni; Motch, Christian; Nättilä, Joonas; Nardini, Emanuele; Neubert, Torsten; O'Brien, Paul T.; Orlandini, Mauro; Osborne, Julian P.; Pacciani, Luigi; Paltani, Stéphane; Paolillo, Maurizio; Papadakis, Iossif E.; Paul, Biswajit; Pellizzoni, Alberto; Peretz, Uria; Torres, Miguel A. Pérez; Perinati, Emanuele; Prescod-Weinstein, Chanda; Reig, Pablo; Riggio, Alessandro; Rodriguez, Jerome; Rodrıguez-Gil, Pablo; Romano, Patrizia; Różańska, Agata; Sakamoto, Takanori; Salmi, Tuomo; Salvaterra, Ruben; Sanna, Andrea; Santangelo, Andrea; Savolainen, Tuomas; Schanne, Stéphane; Schatz, Hendrik; Shao, Lijing; Shearer, Andy; Shore, Steven N.; Stappers, Ben W.; Strohmayer, Tod E.; Suleimanov, Valery F.; Svoboda, Jirı; Thielemann, F. -K.; Tombesi, Francesco; Torres, Diego F.; Torresi, Eleonora; Turriziani, Sara; Vacchi, Andrea; Vercellone, Stefano; Vink, Jacco; Wang, Jian-Min; Wang, Junfeng; Watts, Anna L.; Weng, Shanshan; Weinberg, Nevin N.; Wheatley, Peter J.; Wijnands, Rudy; Woods, Tyrone E.; Woosley, Stan E.; Xiong, Shaolin; Xu, Yupeng; Yan, Zhen; Younes, George; Yu, Wenfei; Yuan, Feng; Zampieri, Luca; Zane, Silvia; Zdziarski, Andrzej A.; Zhang, Shuang-Nan; Zhang, Shu; Zhang, Shuo; Zhang, Xiao; Zingale, Michael
Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)
Submitted: 2018-12-10
In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.
[15]  [pdf] - 1882475
Consistent dust and gas models for protoplanetary disks III. Models for selected objects from the FP7 DIANA project
Comments: 70 pages, 35 figures, 22 tables, accepted by PASP
Submitted: 2018-12-06
The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (
[16]  [pdf] - 1842409
Meter- to Millimeter Emission from Cool Stellar Systems: Latest Results, Synergies Across the Spectrum, and Outlook for the Next Decade
Comments: Splinter session summary, to appear in the proceedings of the 20th Cambridge Workshop on Cool Stars, Stellar Systems, and the Sun (ed. S. J. Wolk)
Submitted: 2018-11-08
Radio observations of cool stellar systems provide unique information on their magnetic fields, high-energy processes, and chemistry. Buoyed by powerful new instruments (e.g. ALMA, JVLA, LOFAR), advances in related fields (e.g., the Gaia astrometric revolution), and above all a renewed interest in the relevant stellar astrophysics, stellar radio astronomy is experiencing a renaissance. In this splinter session, participants took stock of the present state of stellar radio astronomy to chart a course for the field's future.
[17]  [pdf] - 1767793
Science with an ngVLA: Resolving the Radio Complexity of EXor and FUor-type Systems with the ngVLA
Comments: 8 pages, 1 figure, To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA)
Submitted: 2018-10-15
Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, and unprecedented sensitivity will allow for the detailed analysis of outburst systems. The proposed frequency range of the ngVLA allows for observations of the gas, dust, and non-thermal emission from the star and disk.
[18]  [pdf] - 1763104
Interaction of infalling solid bodies with primordial atmospheres of disk-embedded planets
Submitted: 2018-09-14
Planets that form early enough to be embedded in the circumstellar gas disk accumulate thick atmospheres of nebular gas. Models of these atmospheres need to specify the surface luminosity (i.e. energy loss rate) of the planet. This luminosity is usually associated with a continuous inflow of solid bodies, where the gravitational energy released from these bodies is the source of energy. However, if these bodies release energy in the atmosphere instead of at the surface, this assumption might not be justified. Our aim is to explore the interactions of infalling planetesimals with primordial atmospheres at an embedded phase of evolution. We investigate effects of atmospheric interaction on the planetesimals (mass loss) and the atmosphere (heating/cooling). We used atmospheric parameters from a snapshot of time-dependent evolution simulations for embedded atmospheres and simulated purely radial, infall events of siliceous planetesimals in a 1D, explicit code. We implemented energy transfer between friction, radiation transfer by the atmosphere and the body and thermal ablation; this gives us the possibility to examine the effects on the planetesimals and the atmosphere. We find that a significant amount of gravitational energy is indeed dissipated into the atmosphere, especially for larger planetary cores, which consequently cannot contribute to the atmospheric planetary luminosity. Furthermore, we examine that planetesimal infall events for cores, $M_\mathrm{C} > 2$M$_{\oplus}$, which actually result in a local cooling of the atmosphere; this is totally in contradiction with the classical model.
[19]  [pdf] - 1767468
Multiepoch, multiwavelength study of accretion onto T Tauri: X-ray versus optical and UV accretion tracers
Comments: 16 pages, 11 figures, A&A accepted
Submitted: 2018-06-18
Classical T Tauri stars (CTTSs) accrete matter from the inner edge of their surrounding circumstellar disks. The impact of the accretion material on the stellar atmosphere results in a strong shock, which causes emission from the X-ray to the near-infrared (NIR) domain. Shock velocities of several 100 km s$^{-1}$ imply that the immediate post shock plasma emits mainly in X-rays. Indeed, two X-ray diagnostics, the so-called soft excess and the high densities observed in He-like triplets, differentiate CTTSs from their non-accreting siblings. However, accretion shock properties derived from X-ray diagnostics often contradict established ultraviolet (UV)-NIR accretion tracers and a physical model simultaneously explaining both, X-ray and UV-NIR accretion tracers, is not yet available. We present new XMM-Newton and Chandra grating observations of the CTTS T Tauri combined with UV and optical data. During all epochs, the soft excess is large and the densities derived from the O VII and Ne IX He-like triplets are compatible with coronal densities. This confirms that the soft X-ray emission cannot originate in accretion funnels that carry the bulk of the accretion rate despite T Tauri's large soft excess. Instead, we propose a model of radially density stratified accretion columns to explain the density diagnostics and the soft excess. In addition, accretion rate and X-ray luminosity are inversely correlated in T Tauri over several epochs. Such an anti-correlation has been observed in samples of stars. Hence the process causing it must be intrinsic to the accretion process, and we speculate that the stellar magnetic field configuration on the visibile hemisphere affects both the accretion rate and the coronal emission, eventually causing the observed anti-correlation.
[20]  [pdf] - 1755951
The Upper Atmospheres of Terrestrial Planets: Carbon Dioxide Cooling and the Earth's Thermospheric Evolution
Comments: 37 pages, 14 figures, to be published in A&A
Submitted: 2018-06-18
Context: The thermal and chemical structures of the upper atmospheres of planets crucially influence losses to space and must be understood to constrain the effects of losses on atmospheric evolution. Aims: We develop a 1D first-principles hydrodynamic atmosphere model that calculates atmospheric thermal and chemical structures for arbitrary planetary parameters, chemical compositions, and stellar inputs. We apply the model to study the reaction of the Earth's upper atmosphere to large changes in the CO$_2$ abundance and to changes in the input solar XUV field due to the Sun's activity evolution from 3~Gyr in the past to 2.5~Gyr in the future. Methods: For the thermal atmosphere structure, we consider heating from the absorption of stellar X-ray, UV, and IR radiation, heating from exothermic chemical reactions, electron heating from collisions with non-thermal photoelectrons, Joule heating, cooling from IR emission by several species, thermal conduction, and energy exchanges between the neutral, ion, and electron gases. For the chemical structure, we consider $\sim$500 chemical reactions, including 56 photoreactions, eddy and molecular diffusion, and advection. In addition, we calculate the atmospheric structure by solving the hydrodynamic equations. To solve the equations in our model, we develop the Kompot code and provide detailed descriptions of the numerical methods used in the appendices. Results: We verify our model by calculating the structures of the upper atmospheres of the modern Earth and Venus. By varying the CO$_2$ abundances at the lower boundary (65~km) of our Earth model, we show that the atmospheric thermal structure is significantly altered. [Abstract Truncated]
[21]  [pdf] - 1622341
X-ray radiative transfer in protoplanetary disks - The role of dust and X-ray background fields
Comments: 16 pages, 12 figures, accepted for publication in A&A
Submitted: 2017-11-20
The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code ProDiMo, which includes X-ray scattering and absorption by both the gas and dust component. For the X-ray radiative transfer, we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well-established disk ionization tracers N2H+ and HCO+. For evolved dust populations, X-ray opacities are mostly dominated by the gas; only for photon energies $E\gtrsim5-10\,$keV, dust opacities become relevant. Consequently, the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii $r\gtrsim20\,$au. However, the N2H+ and HCO+ column densities are only significantly affected in case of low cosmic-ray ionization rates, or if the background flux is at least a factor of ten higher than the flux level expected for clusters typical for the solar vicinity. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field. For the majority of the cluster members, the X-ray background field has only little impact on the disk chemical structure.
[22]  [pdf] - 1585643
Radio Exploration of Planetary Habitability: Conference Summary
Comments: Five pages; conference Web site:
Submitted: 2017-07-07
Radio Exploration of Planetary Habitability was the fifth in the series of American Astronomical Society's Topical Conference Series. Notable aspects of the conference included the interdisciplinary nature of both the topics and the intellectual breadth of the participants, the diversity of approaches to studying this topic presented by recent discoveries and of the participants themselves, the expanding meaning of the topic of "star-planet interactions," and the expectation of an increasingly statistical approach to the topic. Potential areas of future research include the actual extent to which planetary magnetic fields shield planetary atmospheres; the planetary dynamo process itself, particularly once multiple extrasolar planetary magnetic fields are confirmed; and "planet-star interactions." A major major topic of the conference concerned observational opportunities, highlighted by a number of new or upcoming, specialized observatories to observe exoplanets especially at radio wavelengths. This article summarizes these main points of the conference and expands briefly upon these potential avenues for future investigation. A future meeting on this topic, given the variety of data sets being generated over the next few years, is warranted.
[23]  [pdf] - 1583318
The chemistry of episodic accretion in embedded objects. 2D radiation thermo-chemical models of the post-burst phase
Comments: 17 pages, 6 figures, 2 tables. Accepted for publication in A&A
Submitted: 2017-05-10, last modified: 2017-06-01
Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it is possible to identify targets that experienced a recent luminosity burst. We present a new model for episodic accretion chemistry based on the 2D, radiation thermo-chemical disk code ProDiMo. We have extended ProDiMo with a proper treatment for envelope structures. For a representative Class I model, we calculated the chemical abundances in the post-burst phase and produced synthetic observables like intensity maps and radial profiles. During a burst many chemical species, like CO, sublimate from the dust surfaces. As the burst ends they freeze out again (post-burst phase). This freeze-out happens from inside-out due to the radial density gradient in the disk and envelope structure. This inside-out freeze-out produces clear observational signatures in spectral line emission, like rings and distinct features in the slope of radial intensity profiles. We fitted synthetic C18O J=2-1 observations with single and two component fits and find that post-burst images are much better matched by the latter. Comparing the quality of such fits allows identification of post-burst targets in a model-independent way. Our models confirm that it is possible to identify post-burst objects from spatially resolved CO observations. However, to derive proper statistics, like frequencies of bursts, from observations it is important to consider aspects like the inclination and structure of the target and also dust properties as those have a significant impact on the freeze-out timescale.
[24]  [pdf] - 1546405
Radio emission and mass loss rate limits of four young solar-type stars
Submitted: 2017-02-27
Observations of free-free continuum radio emission of four young main-sequence solar-type stars (EK Dra, Pi1 UMa, Chi1 Ori; and Kappa1 Cet) are studied to detect stellar winds or at least to place upper limits on their thermal radio emission, which is dominated by the ionized wind. These stars are excellent proxies for representing the young Sun. Upper limits on mass loss rates are calculated using their observational radio emission. Our aim is to re-examine the faint young Sun paradox by assuming that the young Sun was more massive in its past, and hence to find a possible solution for this famous problem. The observations of our sample are performed with the Karl G. Jansky VLA with excellent sensitivity, using the C-band and the Ku-band. ALMA observations are performed at 100 GHz. For the estimation of the mass loss limits, spherically symmetric winds and stationary, anisotropic, ionized winds are assumed. We compare our results to 1) mass loss rate estimates of theoretical rotational evolution models, and 2) to results of the indirect technique of determining mass loss rates: Lyman-alpha absorption. We are able to derive the most stringent direct upper limits on mass loss so far from radio observations. Two objects, EK Dra and Chi1 Ori, are detected at 6 and 14 GHz down to an excellent noise level. These stars are very active and additional radio emission identified as non-thermal emission was detected, but limits for the mass loss rates of these objects are still derived. The stars Pi1 UMa and Kappa1 Cet were not detected in either C-band or in Ku-band. For these objects we give upper limits to their radio free-free emission and calculate upper limits to their mass loss rates. Finally, we reproduce the evolution of the Sun and derive an estimate for the solar mass of the Sun at a younger age.
[25]  [pdf] - 1581617
Stellar energetic particle ionization in protoplanetary disks around T Tauri stars
Comments: 15 pages, 11 figures, 5 tables. Accepted for publication in A&A
Submitted: 2017-02-24
Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a high amount of energetic particles (E$\gtrsim$10 MeV). The particle flux of the contemporary Sun cannot explain these anomalies. However, similar to T Tauri stars the young Sun was more active and probably produced enough high energy particles to explain those anomalies. We want to study the interaction of stellar energetic particles with the gas component of the disk and identify possible observational tracers of this interaction. We use a 2D radiation thermo-chemical protoplanetary disk code to model a disk representative for T Tauri stars. We use a particle energy distribution derived from solar flare observations and an enhanced stellar particle flux proposed for T Tauri stars. For this particle spectrum we calculate the stellar particle ionization rate throughout the disk with an accurate particle transport model. We study the impact of stellar particles for models with varying X-ray and cosmic-ray ionization rates. We find that stellar particle ionization has a significant impact on the abundances of the common disk ionization tracers HCO$^+$ and N$_2$H$^+$, especially in models with low cosmic-ray ionization rates. In contrast to cosmic rays and X-rays, stellar particles cannot reach the midplane of the disk. Therefore molecular ions residing in the disk surface layers are more affected by stellar particle ionization than molecular ions tracing the cold layers/midplane of the disk. Spatially resolved observations of molecular ions tracing different vertical layers of the disk allow to disentangle the contribution of stellar particle ionization from other competing ionization sources. Modeling such observations with a model like the one presented here allows to constrain the stellar particle flux in disks around T Tauri stars.
[26]  [pdf] - 1531689
A gas density drop in the inner 6 AU of the transition disk around the Herbig Ae star HD 139614: Further evidence for a giant planet inside the disk?
Comments: accepted for publication in A&A; 30 pages (including Appendix), 24 figure; v2-v4: text corrections
Submitted: 2016-09-21, last modified: 2016-12-09
Context: Quantifying the gas content inside the dust gaps of transition disks is important to establish their origin. Aims: We seek to constrain the surface density of warm gas in the disk of HD 139614, a Herbig Ae star with a transition disk exhibiting a dust gap from 2.3 to 6 AU. Methods: We have obtained ESO/VLT CRIRES high-resolution spectra of CO ro-vibrational emission. We derived constraints on the disk's structure by modeling the line-profiles, the spectroastrometric signal, and the rotational diagrams using flat Keplerian disk models. Results: We detected v=1-0 12CO, 2-1 12CO, 1-0 13CO, 1-0 C18O, and 1-0 C17O ro-vibrational lines. 12CO v=1-0 lines have an average width of 14 km/s, Tgas of 450 K and an emitting region from 1 to 15 AU. 13CO and C18O lines are on average 70 and 100 K colder, 1 and 4 km/s narrower, and are dominated by emission at R>6 AU. The 12CO v=1-0 line-profile indicates that if there is a gap in the gas it must be narrower than 2 AU. We find that a drop in the gas surface density (delta_gas) at R<5-6 AU is required to simultaneously reproduce the line-profiles and rotational diagrams of the three CO isotopologs. Delta_gas can range from 10^-2 to 10^-4 depending on the gas-to-dust ratio of the outer disk. We find that at 1<R<6 AU the gas surface density profile is flat or increases with radius. We derive a gas column density at 1<R<6 AU of NH=3x10^19 - 10^21 cm^-2. We find a 5sigma upper limit on NCO at R<1 AU of 5x10^15 cm^-2 (NH<5x10^19 cm^-2). Conclusions: The dust gap in the disk of HD 139614 has gas. The gas surface density in the disk at R<6 AU is significantly lower than the surface density expected from HD 139614's accretion rate assuming a viscous alpha-disk model. The gas density drop, the non-negative density gradient of the gas inside 6 AU, and the absence of a wide (>2 AU) gas gap suggest the presence of an embedded <2 MJ planet at around 4 AU.
[27]  [pdf] - 1532876
The gas disk: Evolution and chemistry
Comments: 45 pages, 15 figures. Accepted for publication as a chapter in the International Space Science Institute (ISSI) Book on "The Disk in Relation to the Formation of Planets and their Proto-atmospheres", published in Space Science Reviews by Springer
Submitted: 2016-11-29
Protoplanetary disks are the birthplaces of planetary systems. The evolution of the star-disk system and the disk chemical composition determines the initial conditions for planet formation. Therefore a comprehensive understanding of the main physical and chemical processes in disks is crucial for our understanding of planet formation. We give an overview of the early evolution of disks, discuss the importance of the stellar high-energy radiation for disk evolution and describe the general thermal and chemical structure of disks. Finally we provide an overview of observational tracers of the gas component and disk winds.
[28]  [pdf] - 1531329
The feedback of atomic jets from embedded protostars in NGC 1333
Comments: 14 pages, 8 figures, 3 tables. Accepted for publication in A&A
Submitted: 2016-08-22
Star-formation feedback onto the parent cloud is conventionally examined through the study of molecular outflows. Little is however known on the effect that atomic ejecta, tracing fast shocks, can have on the global cloud properties. In this study we employ Herschel/PACS [OI] and [CII] spectral line maps of the NGC 1333 star-forming region to assess the relative influence of atomic jets onto the star-formation process. Atomic line maps are compared against molecular outflow tracers and atomic ejecta are associated to individual driving sources. We study the detailed morphology and velocity distribution of [OI] line using channel and line-centroid maps and derive the momentum, energy, and mass flux for all the bipolar jets traced by [OI] line emission. We find that the line-centroid maps can trace velocity structures down to 5 km s$^{-1}$ which is a factor of $\sim$20 beyond the nominal velocity resolution reached by Herschel/PACS. These maps reveal an unprecedented degree of details that assist significantly in the association and characterization of jets and outflows. Comparisons of the dynamical and kinematical properties shows that [OI] momentum accounts for only $\sim$1% of the momentum carried by the large scale CO outflows but the energy released through the jets corresponds to 50 - 100% of the energy released in outflows. The estimated ratios of the jet to the outflow momenta and energies are consistent with the results of two-component, nested jet/outflow simulations, where jets are associated to episodic accretion events. Under this scenario, the energy from atomic jets to the cloud is as important as the energy output from outflows in maintaining turbulence and dissipating the cloud gas.
[29]  [pdf] - 1521239
Athena Wide Field Imager Key Science Drivers
Comments: 11 pages, 10 figures, 1 table, SPIE Astronomical Telescopes and Instrumentation, Edinburgh (26 June - 1 July, 2016)
Submitted: 2016-07-04
The Wide Field Imager (WFI) is one of two instruments for the Advanced Telescope for High-ENergy Astrophysics (Athena). In this paper we summarise three of the many key science objectives for the WFI - the formation and growth of supermassive black holes, non-gravitational heating in clusters of galaxies, and spin measurements of stellar mass black holes - and describe their translation into the science requirements and ultimately instrument requirements. The WFI will be designed to provide excellent point source sensitivity and grasp for performing wide area surveys, surface brightness sensitivity, survey power, and absolute temperature and density calibration for in-depth studies of the outskirts of nearby clusters of galaxies and very good high-count rate capability, throughput, and low pile-up, paired with very good spectral resolution, for detailed explorations of bright Galactic compact objects.
[30]  [pdf] - 1396808
Thermal mass loss of protoplanetary cores with hydrogen-dominated atmospheres: The influences of ionization and orbital distance
Comments: 8 pages, 2 figures, 1 table, submitted to MNRAS
Submitted: 2016-01-04, last modified: 2016-02-23
We investigate the loss rates of the hydrogen atmospheres of terrestrial planets with a range of masses and orbital distances by assuming a stellar extreme ultraviolet (EUV) luminosity that is 100 times stronger than that of the current Sun. We apply a 1D upper atmosphere radiation absorption and hydrodynamic escape model that takes into account ionization, dissociation and recombination to calculate hydrogen mass loss rates. We study the effects of the ionization, dissociation and recombination on the thermal mass loss rates of hydrogen-dominated super-Earths and compare the results to those obtained by the energy-limited escape formula which is widely used for mass loss evolution studies. Our results indicate that the energy-limited formula can to a great extent over- or underestimate the hydrogen mass loss rates by amounts that depend on the stellar EUV flux and planetary parameters such as mass, size, effective temperature, and EUV absorption radius.
[31]  [pdf] - 1449880
The Evolution of Stellar Rotation and the hydrogen atmospheres of habitable-zone Terrestrial Planets
Comments: Accepted for publication in ApJ Letters
Submitted: 2015-11-11
Terrestrial planets formed within gaseous protoplanetary disks can accumulate significant hydrogen envelopes. The evolution of such an atmosphere due to XUV driven evaporation depends on the activity evolution of the host star, which itself depends sensitively on its rotational evolution, and therefore on its initial rotation rate. In this letter, we derive an easily applicable method for calculating planetary atmosphere evaporation that combines models for a hydrostatic lower atmosphere and a hydrodynamic upper atmosphere. We show that the initial rotation rate of the central star is of critical importance for the evolution of planetary atmospheres and can determine if a planet keeps or loses its primordial hydrogen envelope. Our results highlight the need for a detailed treatment of stellar activity evolution when studying the evolution of planetary atmospheres.
[32]  [pdf] - 1351565
Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs
Comments: accepted by A&A, 36 pages, 30 figures, 5 tables
Submitted: 2015-11-11
We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.
[33]  [pdf] - 1331062
Mid-IR spectra of Pre-Main Sequence Herbig stars: an explanation for the non-detections of water lines
Submitted: 2015-10-27
The mid-IR detection rate of water lines in disks around Herbig stars disks is about 5\%, while it is around 50\% for disks around TTauri stars. The reason for this is still unclear. In this study, we want to find an explanation for the different detection rates between low mass and high mass pre-main-sequence stars (PMSs) in the mid-IR regime. We run disk models with stellar parameters adjusted to spectral types B9 through M2, using the radiation thermo-chemical disk modeling code ProDiMo. We produce convolved spectra at the resolution of Spitzer IRS, JWST MIRI and VLT VISIR spectrographs. We apply random noise derived from typical Spitzer spectra for a direct comparison with observations. The strength of the mid-IR water lines correlates directly with the luminosity of the central star. We explored a small parameter space around a standard disk model, considering dust-to-gas mass ratio, disk gas mass, mixing coefficient for dust settling, flaring index, dust maximum size and size power law distribution index. The models show that it is possible to suppress the water emission, however, current observations are not sensitive enough to detect mid-IR lines in disks for most of the explored parameters. The presence of noise in the spectra, combined with the high continuum flux (noise level is proportional to the continuum flux), is the most likely explanation for the non detections towards Herbig stars. Mid-IR spectra with resolution higher than 20000 are needed to investigate water in protoplanetary disks. Intrinsic differences in disk structure, e.g. inner gaps, gas-to-dust ratio, dust size and distribution, and inner disk scale height, between Herbig and TTauri star disks are able to explain a lower water detection rate in disks around Herbig stars.
[34]  [pdf] - 1224999
Origin and Stability of Exomoon Atmospheres - Implications for Habitability
Comments: 27 pages, 3 figures
Submitted: 2015-06-22
We study the origin and escape of catastrophically outgassed volatiles (H$_2$O, CO$_2$) from exomoons with Earth-like densities and masses of $0.1M_{\oplus}$, $0.5M_{\oplus}$ and $1M_{\oplus}$ orbiting an extra-solar gas giant inside the habitable zone of a young active solar-like star. We apply a radiation absorption and hydrodynamic upper atmosphere model to the three studied exomoon cases. We model the escape of hydrogen and dragged dissociation products O and C during the activity saturation phase of the young host star. Because the soft X-ray and EUV radiation of the young host star may be up to $\sim$100 times higher compared to today's solar value during the first 100 Myr after the system's origin, an exomoon with a mass $ < 0.25M_{\oplus}$ located in the HZ may not be able to keep an atmosphere because of its low gravity. Depending on the spectral type and XUV activity evolution of the host star, exomoons with masses between $\sim0.25-0.5M_{\oplus}$ may evolve to Mars-like habitats. More massive bodies with masses $ > 0.5M_{\oplus}$, however, may evolve to habitats that are a mixture of Mars-like and Earth-analogue habitats, so that life may originate and evolve at the exomoon's surface.
[35]  [pdf] - 1219522
Stellar magnetic activity and their influence on the habitability of exoplanets
Comments: 6 pages, 1 figure
Submitted: 2015-06-19
Stellar magnetism, explorable via polarimetry, is a crucial driver of activity, ionization, photodissociation, chemistry and winds in stellar environments. Thus it has an important impact on the atmospheres and magnetospheres of surrounding planets. Modeling of stellar magnetic fields and their winds is extremely challenging, both from the observational and the theoretical points of view, and only recent ground breaking advances in observational instrumentation - as were discussed during this Symposium - and a deeper theoretical understanding of magnetohydrodynamic processes in stars enable us to model stellar magnetic fields and winds and the resulting influence on surrounding planets in more and more detail. We have initiated a national and international research network (NFN): 'Pathways to Habitability - From Disks to Active Stars, Planets to Life', to address questions on the formation and habitability of environments in young, active stellar/planetary systems. In this contribution we discuss the work we are carrying out within this project and focus on how stellar magnetic fields, their winds and the relation to stellar rotation can be assessed observationally with relevant techniques such as Zeeman Doppler Imaging (ZDI), field extrapolation and wind simulations.
[36]  [pdf] - 1245882
Short-period stellar activity cycles with Kepler photometry
Comments: 37 pages, 17 figures
Submitted: 2015-06-18
We study the short-periodic component of stellar activity with a cycle periods Pcyc up to 1000 days using the Kepler mission photometry of fast-rotating (rotational periods from 1 to 4 days) stars with spectra of M4V to F3V. Applying the originally developed two non-spectral methods, we measured the effective period of stellar cycles in 462 objects. The obtained results are in accordance with previous measurements by Vida et al. (2014), do not seem to result from a beating effect. The performed measurements of Pcyc cluster in a specific branch which covers the previously unstudied region in the Saar-Brandenburg (1999) diagram, and connects the branch of inactive stars with the area populated by super-active objects. It is shown that the formation of the discovered branch is due to the alpha-quenching effect, which saturates the magnetic dynamo and decreases the cycle periods with the increase of inverted Rossby number. This finding is important in the context of the discussion on catastrophic quenching and other heuristic approximations of the non-linear alpha-effect.
[37]  [pdf] - 1044562
The Coronal Temperatures of Low-Mass Main-Sequence Stars
Comments: 5 pages, 3 figures, 1 table; accepted for publication in A&A
Submitted: 2015-05-04, last modified: 2015-05-19
Aims. We study the X-ray emission of low-mass main-sequence stars to derive a reliable general scaling law between coronal temperature and the level of X-ray activity. Methods. We collect ROSAT measurements of hardness ratios and X-ray luminosities for a large sample of stars to derive which stellar X-ray emission parameter is most closely correlated with coronal temperature. We calculate average coronal temperatures for a sample of 24 low-mass main-sequence stars with measured emission measure distributions (EMDs) collected from the literature. These EMDs are based on high-resolution X-ray spectra measured by XMM-Newton and Chandra. Results. We confirm that there is one universal scaling relation between coronal average temperature and surface X-ray flux, Fx, that applies to all low-mass main-sequence stars. We find that coronal temperature is related to Fx by Tcor=0.11 Fx^0.26, where Tcor is in MK and Fx is in erg/s/cm^2.
[38]  [pdf] - 1044543
Stellar Winds on the Main-Sequence II: the Evolution of Rotation and Winds
Comments: 26 pages, 14 figures, 2 tables, to be published in A&A
Submitted: 2015-03-25, last modified: 2015-05-19
Aims: We study the evolution of stellar rotation and wind properties for low-mass main-sequence stars. Our aim is to use rotational evolution models to constrain the mass loss rates in stellar winds and to predict how their properties evolve with time on the main-sequence. Methods: We construct a rotational evolution model that is driven by observed rotational distributions of young stellar clusters. Fitting the free parameters in our model allows us to predict how wind mass loss rate depends on stellar mass, radius, and rotation. We couple the results to the wind model developed in Paper I of this series to predict how wind properties evolve on the main-sequence. Results: We estimate that wind mass loss rate scales with stellar parameters as $\dot{M}_\star \propto R_\star^2 \Omega_\star^{1.33} M_\star^{-3.36}$. We estimate that at young ages, the solar wind likely had a mass loss rate that is an order of magnitude higher than that of the current solar wind. This leads to the wind having a higher density at younger ages; however, the magnitude of this change depends strongly on how we scale wind temperature. Due to the spread in rotation rates, young stars show a large range of wind properties at a given age. This spread in wind properties disappears as the stars age. Conclusions: There is a large uncertainty in our knowledge of the evolution of stellar winds on the main-sequence, due both to our lack of knowledge of stellar winds and the large spread in rotation rates at young ages. Given the sensitivity of planetary atmospheres to stellar wind and radiation conditions, these uncertainties can be significant for our understanding of the evolution of planetary environments.
[39]  [pdf] - 1055964
The Extreme Ultraviolet and X-Ray Sun in Time: High-Energy Evolutionary Tracks of a Solar-Like Star
Comments: 4 pages, 4 figures, accepted for publication in A&A
Submitted: 2015-04-17
Aims. We aim to describe the pre-main sequence and main-sequence evolution of X-ray and extreme-ultaviolet radiation of a solar mass star based on its rotational evolution starting with a realistic range of initial rotation rates. Methods. We derive evolutionary tracks of X-ray radiation based on a rotational evolution model for solar mass stars and the rotation-activity relation. We compare these tracks to X-ray luminosity distributions of stars in clusters with different ages. Results. We find agreement between the evolutionary tracks derived from rotation and the X-ray luminosity distributions from observations. Depending on the initial rotation rate, a star might remain at the X-ray saturation level for very different time periods, approximately from 10 Myr to 300 Myr for slow and fast rotators, respectively. Conclusions. Rotational evolution with a spread of initial conditions leads to a particularly wide distribution of possible X-ray luminosities in the age range of 20 to 500 Myrs, before rotational convergence and therefore X-ray luminosity convergence sets in. This age range is crucial for the evolution of young planetary atmospheres and may thus lead to very different planetary evolution histories.
[40]  [pdf] - 985565
Stellar Winds on the Main-Sequence I: Wind Model
Comments: 24 pages, 13 figures, 2 tables, Accepted for publication in A&A
Submitted: 2015-03-23
Aims: We develop a method for estimating the properties of stellar winds for low-mass main-sequence stars between masses of 0.4 and 1.1 solar masses at a range of distances from the star. Methods: We use 1D thermal pressure driven hydrodynamic wind models run using the Versatile Advection Code. Using in situ measurements of the solar wind, we produce models for the slow and fast components of the solar wind. We consider two radically different methods for scaling the base temperature of the wind to other stars: in Model A, we assume that wind temperatures are fundamentally linked to coronal temperatures, and in Model B, we assume that the sound speed at the base of the wind is a fixed fraction of the escape velocity. In Paper II of this series, we use observationally constrained rotational evolution models to derive wind mass loss rates. Results: Our model for the solar wind provides an excellent description of the real solar wind far from the solar surface, but is unrealistic within the solar corona. We run a grid of 1200 wind models to derive relations for the wind properties as a function of stellar mass, radius, and wind temperature. Using these results, we explore how wind properties depend on stellar mass and rotation. Conclusions: Based on our two assumptions about the scaling of the wind temperature, we argue that there is still significant uncertainty in how these properties should be determined. Resolution of this uncertainty will probably require both the application of solar wind physics to other stars and detailed observational constraints on the properties of stellar winds. In the final section of this paper, we give step by step instructions for how to apply our results to calculate the stellar wind conditions far from the stellar surface.
[41]  [pdf] - 1005643
Colliding Winds in Low-Mass Binary Star Systems: wind interactions and implications for habitable planets
Comments: 15 pages, 11 figures, to be published in A&A
Submitted: 2015-02-11
Context. In binary star systems, the winds from the two components impact each other, leading to strong shocks and regions of enhanced density and temperature. Potentially habitable circumbinary planets must continually be exposed to these interactions regions. Aims. We study, for the first time, the interactions between winds from low-mass stars in a binary system, to show the wind conditions seen by potentially habitable circumbinary planets. Methods. We use the advanced 3D numerical hydrodynamic code Nurgush to model the wind interactions of two identical winds from two solar mass stars with circular orbits and a binary separation of 0.5 AU. As input into this model, we use a 1D hydrodynamic simulation of the solar wind, run using the Versatile Advection Code. We derive the locations of stable and habitable orbits in this system to explore what wind conditions potentially habitable planets will be exposed to during their orbits. Results. Our wind interaction simulations result in the formation of two strong shock waves separated by a region of enhanced density and temperature. The wind-wind interaction region has a spiral shape due to Coriolis forces generated by the orbital motions of the two stars. The stable and habitable zone in this system extends from approximately 1.4 AU to 2.4 AU. (TRUNCATED)
[42]  [pdf] - 920037
Stellar flares observed by LOFT: implications for the physics of coronae and for the "space weather" environment of extrasolar planets
Comments: White Paper in Support of the Mission Concept of the Large Observatory for X-ray Timing
Submitted: 2015-01-12
This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of stellar flares. For a summary, we refer to the paper.
[43]  [pdf] - 923407
Impact inducted surface heating by planetesimals on early Mars
Comments: submitted to A&A
Submitted: 2014-05-22
We investigate the influence of impacts of large planetesimals and small planetary embryos on the early Martian surface on the hydrodynamic escape of an early steam atmosphere that is exposed to the high soft X-ray and EUV flux of the young Sun. Impact statistics in terms of number, masses, velocities, and angles of asteroid impacts onto the early Mars are determined via n-body integrations. Based on these statistics, smoothed particle hydrodynamics (SPH) simulations result in estimates of energy transfer into the planetary surface material and according surface heating. For the estimation of the atmospheric escape rates we applied a soft X-ray and EUV absorption model and a 1-D upper atmosphere hydrodynamic model to a magma ocean-related catastrophically outgassed steam atmosphere with surface pressure values of 52 bar H2O and 11 bar CO2. The estimated impact rates and energy deposition onto an early Martian surface can account for substantial heating. The energy influx and conversion rate into internal energy is most likely sufficient to keep a shallow magma ocean liquid for an extended period of time. Higher surface temperatures keep the outgassed steam atmosphere longer in vapor form and therefore enhance its escape to space within about 0.6 Myr after its formation.
[44]  [pdf] - 840013
First Detection of Thermal Radio Emission from Solar-Type Stars with the Karl G. Jansky Very Large Array
Comments: 9 pages, 2 figures. Accepted for publication in ApJ
Submitted: 2014-05-09
We present the first detections of thermal radio emission from the atmospheres of solar-type stars {\tau} Cet, {\eta} Cas A, and 40 Eri A. These stars all resemble the Sun in age and level of magnetic activity, as indicated by X-ray luminosity and chromospheric emission in calcium-II H and K lines. We observed these stars with the Karl G. Jansky Very Large Array with sensitivities of a few {\mu}Jy at combinations of 10.0, 15.0, and 34.5 GHz. {\tau} Cet, {\eta} Cas A, and 40 Eri A are all detected at 34.5 GHz with signal-to-noise ratios of 6.5, 5.2, and 4.5, respectively. 15.0-GHz upper limits imply a rising spectral index greater than 1.0 for {\tau} Cet and 1.6 for {\eta} Cas A, at the 95% confidence level. The measured 34.5-GHz flux densities correspond to stellar disk-averaged brightness temperatures of roughly 10,000 K, similar to the solar brightness temperature at the same frequency. We explain this emission as optically- thick thermal free-free emission from the chromosphere, with possible contributions from coronal gyroresonance emission above active regions and coronal free-free emission. These and similar quality data on other nearby solar-type stars, when combined with ALMA observations, will enable the construction of temperature profiles of their chromospheres and lower transition regions.
[45]  [pdf] - 1179670
The PLATO 2.0 Mission
Rauer, H.; Catala, C.; Aerts, C.; Appourchaux, T.; Benz, W.; Brandeker, A.; Christensen-Dalsgaard, J.; Deleuil, M.; Gizon, L.; Goupil, M. -J.; Güdel, M.; Janot-Pacheco, E.; Mas-Hesse, M.; Pagano, I.; Piotto, G.; Pollacco, D.; Santos, N. C.; Smith, A.; -C., J.; Suárez; Szabó, R.; Udry, S.; Adibekyan, V.; Alibert, Y.; Almenara, J. -M.; Amaro-Seoane, P.; Eiff, M. Ammler-von; Asplund, M.; Antonello, E.; Ball, W.; Barnes, S.; Baudin, F.; Belkacem, K.; Bergemann, M.; Bihain, G.; Birch, A. C.; Bonfils, X.; Boisse, I.; Bonomo, A. S.; Borsa, F.; Brandão, I. M.; Brocato, E.; Brun, S.; Burleigh, M.; Burston, R.; Cabrera, J.; Cassisi, S.; Chaplin, W.; Charpinet, S.; Chiappini, C.; Church, R. P.; Csizmadia, Sz.; Cunha, M.; Damasso, M.; Davies, M. B.; Deeg, H. J.; DÍaz, R. F.; Dreizler, S.; Dreyer, C.; Eggenberger, P.; Ehrenreich, D.; Eigmüller, P.; Erikson, A.; Farmer, R.; Feltzing, S.; Fialho, F. de Oliveira; Figueira, P.; Forveille, T.; Fridlund, M.; García, R. A.; Giommi, P.; Giuffrida, G.; Godolt, M.; da Silva, J. Gomes; Granzer, T.; Grenfell, J. L.; Grotsch-Noels, A.; Günther, E.; Haswell, C. A.; Hatzes, A. P.; Hébrard, G.; Hekker, S.; Helled, R.; Heng, K.; Jenkins, J. M.; Johansen, A.; Khodachenko, M. L.; Kislyakova, K. G.; Kley, W.; Kolb, U.; Krivova, N.; Kupka, F.; Lammer, H.; Lanza, A. F.; Lebreton, Y.; Magrin, D.; Marcos-Arenal, P.; Marrese, P. M.; Marques, J. P.; Martins, J.; Mathis, S.; Mathur, S.; Messina, S.; Miglio, A.; Montalban, J.; Montalto, M.; Monteiro, M. J. P. F. G.; Moradi, H.; Moravveji, E.; Mordasini, C.; Morel, T.; Mortier, A.; Nascimbeni, V.; Nelson, R. P.; Nielsen, M. B.; Noack, L.; Norton, A. J.; Ofir, A.; Oshagh, M.; Ouazzani, R. -M.; Pápics, P.; Parro, V. C.; Petit, P.; Plez, B.; Poretti, E.; Quirrenbach, A.; Ragazzoni, R.; Raimondo, G.; Rainer, M.; Reese, D. R.; Redmer, R.; Reffert, S.; Rojas-Ayala, B.; Roxburgh, I. W.; Salmon, S.; Santerne, A.; Schneider, J.; Schou, J.; Schuh, S.; Schunker, H.; Silva-Valio, A.; Silvotti, R.; Skillen, I.; Snellen, I.; Sohl, F.; Sousa, S. G.; Sozzetti, A.; Stello, D.; Strassmeier, K. G.; Švanda, M.; Szabó, Gy. M.; Tkachenko, A.; Valencia, D.; van Grootel, V.; Vauclair, S. D.; Ventura, P.; Wagner, F. W.; Walton, N. A.; Weingrill, J.; Werner, S. C.; Wheatley, P. J.; Zwintz, K.
Comments: 63 pages, 17 figures, submitted to Experimental Astronomy (ExA)
Submitted: 2013-10-02, last modified: 2014-03-04
PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.
[46]  [pdf] - 1944175
Jets and Outflows From Star to Cloud: Observations Confront Theory
Comments: Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning
Submitted: 2014-02-14
In this review we focus on the role jets and outflows play in the star and planet formation process. Our essential question can be posed as follows: are jets/outflows merely an epiphenomenon associated with star formation or do they play an important role in mediating the physics of assembling stars both individually and globally? We address this question by reviewing the current state of observations and their key points of contact with theory. Our review of jet/outflow phenomena is organized into three length-scale domains: Source and Disk Scales ($0.1-10^2$ au) where the connection with protostellar and disk evolution theories is paramount; Envelope Scales ($10^2-10^5$ au) where the chemistry and propagation shed further light on the jet launching process, its variability and its impact on the infalling envelope; Parent Cloud Scales ($10^5-10^6$ au) where global momentum injection into cluster/cloud environments become relevant. Issues of feedback are of particular importance on the smallest scales where planet formation regions in a disk may be impacted by the presence of disk winds, irradiation by jet shocks or shielding by the winds. Feedback on envelope scales may determine the final stellar mass (core-to-star efficiency) and envelope dissipation. Feedback also plays an important role on the larger scales with outflows contributing to turbulent support within clusters including alteration of cluster star formation efficiencies (feedback on larger scales currently appears unlikely). A particularly novel dimension of our review is that we consider results on jet dynamics from the emerging field of High Energy Density Laboratory Astrophysics (HEDLA). HEDLA is now providing direct insights into the 3-D dynamics of fully magnetized, hypersonic, radiative outflows.
[47]  [pdf] - 1202299
Stellar wind interaction and pick-up ion escape of the Kepler-11 "super-Earths"
Comments: 8 pages, 3 figures, accepted to A&A
Submitted: 2013-12-17, last modified: 2014-01-22
We study the interactions between stellar wind and the extended hydrogen-dominated upper atmospheres of planets and the resulting escape of planetary pick-up ions from the 5 "super-Earths" in the compact Kepler-11 system and compare the escape rates with the efficiency of the thermal escape of neutral hydrogen atoms. Assuming the stellar wind of Kepler-11 is similar to the solar wind, we use a polytropic 1D hydrodynamic wind model to estimate the wind properties at the planetary orbits. We apply a Direct Simulation Monte Carlo Model to model the hydrogen coronae and the stellar wind plasma interaction around Kepler-11b-f within a realistic expected heating efficiency range of 15-40%. The same model is used to estimate the ion pick-up escape from the XUV heated and hydrodynamically extended upper atmospheres of Kepler-11b-f. From the interaction model we study the influence of possible magnetic moments, calculate the charge exchange and photoionization production rates of planetary ions and estimate the loss rates of pick-up H+ ions for all five planets. We compare the results between the five "super-Earths" and in a more general sense also with the thermal escape rates of the neutral planetary hydrogen atoms. Our results show that for all Kepler-11b-f exoplanets, a huge neutral hydrogen corona is formed around the planet. The non-symmetric form of the corona changes from planet to planet and is defined mostly by radiation pressure and gravitational effects. Non-thermal escape rates of pick-up ionized hydrogen atoms for Kepler-11 "super-Earths" vary between approximately 6.4e30 1/s and 4.1e31 1/s depending on the planet's orbital location and assumed heating efficiency. These values correspond to non-thermal mass loss rates of approximately 1.07e7 g/s and 6.8e7 g/s respectively, which is a few percent of the thermal escape rates.
[48]  [pdf] - 791889
Atomic jet from SMM1 (FIRS1) in Serpens uncovers non-coeval binary companion
Comments: 10 pages, 7 figures, 1 table. Accepted for publication in Astronomy \& Astrophysics
Submitted: 2014-01-14
We report on the detection of an atomic jet associated with the protostellar source SMM1 (FIRS1) in Serpens. The jet is revealed in [FeII] and [NeII] line maps observed with Spitzer/IRS, and further confirmed in HiRes IRAC and MIPS images. It is traced very close to SMM1 and peaks at ~5 arcsec" from the source at a position angle of $\sim 125 degrees. In contrast, molecular hydrogen emission becomes prominent at distances > 5" from the protostar and extends at a position angle of 160 degrees. The morphological differences suggest that the atomic emission arises from a companion source, lying in the foreground of the envelope surrounding the embedded protostar SMM1. In addition the molecular and atomic Spitzer maps disentangle the large scale CO (3-2) emission observed in the region into two distinct bipolar outflows, giving further support to a proto-binary source setup. Analysis at the peaks of the [FeII] jet show that emission arises from warm and dense gas (T ~1000 K, n(electron) 10^5 - 10^6 cm^-3). The mass flux of the jet derived independently for the [FeII] and [NeII] lines is 10^7 M(sun)/yr, pointing to a more evolved Class~I/II protostar as the driving source. All existing evidence converge to the conclusion that SMM1 is a non-coeval proto-binary source.
[49]  [pdf] - 791878
Origin and Loss of nebula-captured hydrogen envelopes from "sub"- to "super-Earths" in the habitable zone of Sun-like stars
Comments: 17 pages, 5 figures, accepted to MNRAS
Submitted: 2014-01-13
We investigate the origin and loss of captured hydrogen envelopes from protoplanets between `sub-Earth'-like bodies of 0.1$M_{\oplus}$ up to `super-Earths' with 5$M_{\oplus}$ in the HZ of a Sun like G star, assuming their rocky cores had formed before the nebula dissipated. We model the gravitational accumulation of nebula gas around a core as a function of protoplanetary luminosity during accretion and calculate the resulting surface temperature by solving the hydrostatic structure equations for the protoplanetary nebula. Depending on nebular properties and resulting luminosities, for planetary bodies of 0.1--1$M_{\oplus}$ we obtain hydrogen envelopes with masses between $\sim 2.5\times 10^{19}$--$1.5\times 10^{26}$ g. For `super-Earths' with masses between 2--5$M_{\oplus}$ hydrogen envelopes within the mass range of $\sim 7.5\times 10^{23}$--$1.5\times 10^{28}$ g can be captured. To study the escape of these hydrogen-dominated protoatmospheres, we apply a hydrodynamic upper atmosphere model and calculate the loss rates due to the heating by the high XUV flux of the young star. Our results indicate that under most nebula conditions `sub-Earth' and Earth-mass planets can lose their envelopes by thermal escape during the first $100$ Myr after the disk dissipated. However, if a nebula has a low dust depletion factor or low accretion rates resulting in low protoplanetary luminosities, it is possible that even protoplanets with Earth-mass cores may keep their hydrogen envelopes during their whole lifetime. In contrast to lower mass protoplanets, `super-Earths' accumulate a huge amount of nebula gas and lose only tiny fractions of their primordial envelopes. Our results agree with the fact that Venus, Earth, and Mars are not surrounded by dense hydrogen envelopes, as well as with the recent discoveries of low density `super-Earths' that most likely could not get rid of their protoatmospheres.
[50]  [pdf] - 1179571
Constraining Stellar Winds of Young Sun-like Stars
Comments: 2 pages, 1 figure, IAUS 302
Submitted: 2013-09-29
As part of the project Pathways to Habitability (, we study the properties of the stellar winds of low-mass and Sun-like stars, and their influences on the atmospheres of potentially habitable planets. For this purpose, we combine mapping of stellar magnetic fields with magnetohydrodynamic wind models.
[51]  [pdf] - 1179096
The magnetosphere of the close accreting PMS binary V4046 Sgr AB
Comments: 2 pages, proceedings of IAUS 302 "Magnetic Fields throughout Stellar Evolution", Biarritz, France, August 2013. CUP, in press
Submitted: 2013-09-09
We present a preliminary 3D potential field extrapolation model of the joint magnetosphere of the close accreting PMS binary V4046 Sgr. The model is derived from magnetic maps obtained as part of a coordinated optical and X-ray observing program.
[52]  [pdf] - 1178919
The magnetosphere of the close accreting PMS binary V4046 Sgr
Comments: Proceedings of the "Physics at the Magnetospheric Boundary" conference, Geneva, June 2013
Submitted: 2013-08-31
V4046 Sagittarii AB is a close short-period classical T Tauri binary. It is a circularised and synchronised system accreting from a circumbinary disk. In 2009 it was observed as part of a coordinated program involving near-simultaneous spectropolarimetric observations with ESPaDOnS at the Canada-France-Hawai'i Telescope and high-resolution X-ray observations with XMM-Newton. Magnetic maps of each star were derived from Zeeman-Doppler imaging. After briefly highlighting the most significant observational findings, we present a preliminary 3D model of the binary magnetosphere constructed from the magnetic maps using a newly developed binary magnetic field extrapolation code. The large-scale fields (the dipole components) of both stars are highly tilted with respect to their rotation axes, and their magnetic fields are linked.
[53]  [pdf] - 1173135
Escape of the martian protoatmosphere and initial water inventory
Comments: 47 pages, 10 figures, 3 tables, submitted to PSS
Submitted: 2013-08-01
Latest research in planet formation indicate that Mars formed within a few million years (Myr) and remained a planetary embryo that never grew to a more massive planet. It can also be expected from dynamical models, that most of Mars' building blocks consisted of material that formed in orbital locations just beyond the ice line which could have contained ~0.1-0.2 wt. % of H2O. By using these constraints, we estimate the nebula-captured and catastrophically outgassed volatile contents during the solidification of Mars' magma ocean and apply a hydrodynamic upper atmosphere model for the study of the soft X-ray and extreme ultraviolet (XUV) driven thermal escape of the martian protoatmosphere during the early active epoch of the young Sun. The amount of gas that has been captured from the protoplanetary disk into the planetary atmosphere is calculated by solving the hydrostatic structure equations in the protoplanetary nebula. Depending on nebular properties such as the dust grain depletion factor, planetesimal accretion rates and luminosities, hydrogen envelopes with masses >=3x10^{19} g to <=6.5x10^{22} g could have been captured from the nebula around early Mars. Depending of the before mentioned parameters, due to the planets low gravity and a solar XUV flux that was ~100 times stronger compared to the present value, our results indicate that early Mars would have lost its nebular captured hydrogen envelope after the nebula gas evaporated, during a fast period of ~0.1-7.5 Myr. After the solidification of early Mars' magma ocean, catastrophically outgassed volatiles with the amount of ~50-250 bar H2O and ~10-55 bar CO2 could have been lost during ~0.4-12 Myr, if the impact related energy flux of large planetesimals and small embryos to the planet's surface lasted long enough, that the steam atmosphere could have been prevented from condensing. If this was not the case... (continued)
[54]  [pdf] - 1158638
XUV exposed non-hydrostatic hydrogen-rich upper atmospheres of terrestrial planets. Part I: Atmospheric expansion and thermal escape
Comments: 50 pages, 9 figures, 4 tables, submitted to Astrobiology
Submitted: 2012-12-20, last modified: 2013-06-18
The recently discovered low-density "super-Earths" Kepler-11b, Kepler-11f, Kepler-11d, Kepler-11e, and planets such as GJ 1214b represent most likely planets which are surrounded by dense H/He envelopes or contain deep H2O oceans also surrounded by dense hydrogen envelopes. Although these "super-Earths" are orbiting relatively close to their host stars, they have not lost their captured nebula-based hydrogen-rich or degassed volatile-rich steam protoatmospheres. Thus it is interesting to estimate the maximum possible amount of atmospheric hydrogen loss from a terrestrial planet orbiting within the habitable zone of late main sequence host stars. For studying the thermosphere structure and escape we apply a 1-D hydrodynamic upper atmosphere model which solves the equations of mass, momentum and energy conservation for a planet with the mass and size of the Earth and for a "super-Earth" with a size of 2 R_Earth and a mass of 10 M_Earth. We calculate volume heating rates by the stellar soft X-ray and EUV radiation and expansion of the upper atmosphere, its temperature, density and velocity structure and related thermal escape rates during planet's life time. Moreover, we investigate under which conditions both planets enter the blow-off escape regime and may therefore experience loss rates which are close to the energy-limited escape. Finally we discuss the results in the context of atmospheric evolution and implications for habitability of terrestrial planets in general.
[55]  [pdf] - 1158612
XUV exposed, non-hydrostatic hydrogen-rich upper atmospheres of terrestrial planets II: Hydrogen coronae and ion escape
Comments: 53 pages, 6 figures, 10 tables, submitted to Astrobiology
Submitted: 2012-12-19, last modified: 2013-06-18
We study the interactions between the stellar wind plasma flow of a typical M star, such as GJ 436, and hydrogen-rich upper atmospheres of an Earth-like planet and a "super-Earth" with the radius of 2 R_Earth and a mass of 10 M_Earth, located within the habitable zone at ~0.24 AU. We investigate the formation of extended atomic hydrogen coronae under the influences of the stellar XUV flux (soft X-rays and EUV), stellar wind density and velocity, shape of a planetary obstacle (e.g., magnetosphere, ionopause), and the loss of planetary pick-up ions on the evolution of hydrogen-dominated upper atmospheres. Stellar XUV fluxes which are 1, 10, 50 and 100 times higher compared to that of the present-day Sun are considered and the formation of high-energy neutral hydrogen clouds around the planets due to the charge-exchange reaction under various stellar conditions have been modeled. Charge-exchange between stellar wind protons with planetary hydrogen atoms, and photoionization, leads to the production of initially cold ions of planetary origin. We found that the ion production rates for the studied planets can vary over a wide range, from ~1.0x10^{25} s^{-1} to ~5.3x10^{30} s^{-1}, depending on the stellar wind conditions and the assumed XUV exposure of the upper atmosphere. Our findings indicate that most likely the majority of these planetary ions are picked up by the stellar wind and lost from the planet. Finally, we estimate the long-time non-thermal ion pick-up escape for the studied planets and compare them with the thermal escape. According to our estimates, non-thermal escape of picked up ionized hydrogen atoms over a planet's lifetime varies between ~0.4 Earth ocean equivalent amounts of hydrogen (EO_H) to <3 EO_H and usually is several times smaller in comparison to the thermal atmospheric escape rates.
[56]  [pdf] - 678288
The Hot and Energetic Universe: Star formation and evolution
Comments: Supporting paper for the science theme "The Hot and Energetic Universe" to be implemented by the Athena+ X-ray observatory ( 10 pages, 8 figures
Submitted: 2013-06-10
Stars over a wide range of masses and evolutionary stages are nowadays known to emit X-rays. This X-ray emission is a unique probe of the most energetic phenomena occurring in the circumstellar environment of these stars, and provides precious insight on magnetic phenomena or hydrodynamic shocks. Owing to its large collecting area, Athena+ will open up an entirely new window on these phenomena. Indeed, Athena+ will not only allow us to study many more objects with an unprecedented spectral resolution, but will also pioneer the study of the dynamics of these objects via time-resolved high-resolution spectroscopy. In this way, Athena+ will be a unique tool to study accretion processes in TTauri stars, flaring activity in young stars, dynamos in ultra-cool dwarfs, small and large-scale structures in the winds of single massive stars, wind interactions in massive binary systems, hot bubbles in planetary nebula... All these studies will lead to a deeper understanding of yet poorly understood processes which have profound impact in star and planetary system formation as well as in feedback processes on Galactic scale.
[57]  [pdf] - 678287
The Hot and Energetic Universe: Solar system and exoplanets
Comments: Supporting paper for the science theme "The Hot and Energetic Universe" to be implemented by the Athena+ X-ray observatory ( 8 pages, 6 figures
Submitted: 2013-06-10
The high resolution non-dispersive spectroscopy and unprecedented sensitivity of Athena+ will revolutionize solar system observing: the origin of the ions producing Jupiter's X-ray aurorae via charge exchange will be conclusively established, as well as their dynamics, giving clues to their acceleration mechanisms. X-ray aurorae on Saturn will be searched for to a depth unattainable by current Earth-bound observatories. The X-ray Integral Field Unit of Athena+ will map Mars' expanding exosphere, which has a line-rich solar wind charge exchange spectrum, under differing solar wind conditions and through the seasons; relating Mars' X-ray emission to its atmospheric loss will have significant impact also on the study of exoplanet atmospheres. Spectral mapping of cometary comae, which are spectacular X-ray sources with extremely line-rich spectra, will probe solar wind composition and speed at varying distances from the Sun. Athena+ will provide unique contributions also to exoplanetary astrophysics. Athena+ will pioneer the study of ingress/eclipse/egress effects during planetary orbits of hot-Jupiters, and will confirm/improve the evidence of Star-Planet Interactions (SPI) in a wider sample of planetary systems. Finally Athena+ will drastically improve the knowledge of the X-ray incident radiation on exoplanets, a key element for understanding the effects of atmospheric mass loss and of the chemical and physical evolution of planet atmospheres, particularly relevant in the case of young systems.
[58]  [pdf] - 663130
Characterizing Exoplanets in the Visible and Infrared: A Spectrometer Concept for the EChO Space Mission
Comments: 20 pages, 8 figures, accepted for publication in the Journal of Astronomical Instrumentation
Submitted: 2013-05-14
Transit-spectroscopy of exoplanets is one of the key observational techniques to characterize the extrasolar planet and its atmosphere. The observational challenges of these measurements require dedicated instrumentation and only the space environment allows an undisturbed access to earth-like atmospheric features such as water or carbon-dioxide. Therefore, several exoplanet-specific space missions are currently being studied. One of them is EChO, the Exoplanet Characterization Observatory, which is part of ESA's Cosmic Vision 2015-2025 program, and which is one of four candidates for the M3 launch slot in 2024. In this paper we present the results of our assessment study of the EChO spectrometer, the only science instrument onboard this spacecraft. The instrument is a multi-channel all-reflective dispersive spectrometer, covering the wavelength range from 400 nm to 16 microns simultaneously with a moderately low spectral resolution. We illustrate how the key technical challenge of the EChO mission - the high photometric stability - influences the choice of spectrometer concept and drives fundamentally the instrument design. First performance evaluations underline the fitness of the elaborated design solution for the needs of the EChO mission.
[59]  [pdf] - 636050
VLA Observations of DG Tau's Radio Jet: A highly collimated thermal outflow
Comments: 19 pages, 7 figures
Submitted: 2013-02-04
The active young protostar DG Tau has an extended jet that has been well studied at radio, optical, and X-ray wavelengths. We report sensitive new VLA full-polarization observations of the core and jet between 5 GHz and 8 GHz. Our high angular resolution observation at 8 GHz clearly shows an unpolarized inner jet with a size 42 AU (0.35") extending along a position angle similar to the optical-X ray outer jet. Using our nearly coeval 2012 VLA observations, we find a spectral-index=+0.46+/-0.05, which combined with the lack of polarization, is consistent with bremsstrahlung (free-free) emission, with no evidence for a non-thermal coronal component. By identifying the end of the radio jet as the optical depth unity surface, and calculating the resulting emission measure, we find our radio results are in agreement with previous optical line studies of electron density and consequent mass-loss rate. We also detect a weak radio knot at 5 GHz located 7" from the base of the jet, coincident with the inner radio knot detected by Rodriguez et al. (2012) in 2009 but at lower surface brightness. We interpret this as due to expansion of post-shock ionized gas in the three years between observations.
[60]  [pdf] - 1159073
Discovery of X-ray emission from young suns in the Small Magellanic Cloud
Comments: Accepted for publication in ApJ
Submitted: 2013-01-15
We report the discovery of extended X-ray emission within the young star cluster NGC 602 in the Wing of the Small Magellanic Cloud (SMC) based on observations obtained with the Chandra X-ray Observatory. X-ray emission is detected from the cluster core area with the highest stellar density and from a dusty ridge surrounding the HII region. We use a census of massive stars in the cluster to demonstrate that a cluster wind or wind-blown bubble is unlikely to provide a significant contribution to the X-ray emission detected from the central area of the cluster. We therefore suggest that X-ray emission at the cluster core originates from an ensemble of low- and solar-mass pre-main-sequence (PMS) stars, each of which would be too weak in X-rays to be detected individually. We attribute the X-ray emission from the dusty ridge to the embedded tight cluster of the new-born stars known in this area from infrared studies. Assuming that the levels of X-ray activity in young stars in the low-metallicity environment of NGC 602a are comparable to their Galactic counterparts, then the detected spatial distribution, spectral properties, and level of X-ray emission are largely consistent with those expected from low- and solar-mass PMS stars and young stellar objects (YSOs). This is the first discovery of X-ray emission attributable to PMS stars and YSOs in the SMC, which suggests that the accretion and dynamo processes in young, low-mass objects in the SMC resemble those in the Galaxy.
[61]  [pdf] - 1117768
The close T Tauri binary system V4046 Sgr: Rotationally modulated X-ray emission from accretion shocks
Comments: paper accepted by ApJ
Submitted: 2012-04-04
We report initial results from a quasi-simultaneous X-ray/optical observing campaign targeting V4046 Sgr, a close, synchronous-rotating classical T Tauri star (CTTS) binary in which both components are actively accreting. V4046 Sgr is a strong X-ray source, with the X-rays mainly arising from high-density (n_e ~ 10^(11-12) cm^(-3)) plasma at temperatures of 3-4 MK. Our multiwavelength campaign aims to simultaneously constrain the properties of this X-ray emitting plasma, the large scale magnetic field, and the accretion geometry. In this paper, we present key results obtained via time-resolved X-ray grating spectra, gathered in a 360 ks XMM-Newton observation that covered 2.2 system rotations. We find that the emission lines produced by this high-density plasma display periodic flux variations with a measured period, 1.22+/-0.01 d, that is precisely half that of the binary star system (2.42 d). The observed rotational modulation can be explained assuming that the high-density plasma occupies small portions of the stellar surfaces, corotating with the stars, and that the high-density plasma is not azimuthally symmetrically distributed with respect to the rotational axis of each star. These results strongly support models in which high-density, X-ray-emitting CTTS plasma is material heated in accretion shocks, located at the base of accretion flows tied to the system by magnetic field lines.
[62]  [pdf] - 1084123
GSC 07396-00759 = V4046 Sgr C[D]: a Wide-separation Companion to the Close T Tauri Binary System V4046 Sgr AB
Comments: 14 pages, 3 figures; to appear in the Astrophysical Journal (Letters)
Submitted: 2011-09-15
We explore the possibility that GSC 07396-00759 (spectral type M1e) is a widely separated (~2.82', or projected separation ~12,350 AU) companion to the "old" (age ~12 Myr) classical T Tauri binary system V4046 Sgr AB, as suggested by the proximity and similar space motions of the two systems. If the two systems are equidistant and coeval, then GSC 07396--00759, like V4046 Sgr AB, must be a spectroscopic binary with nearly equal-mass components, and V4046 Sgr must be at least ~8 Myr old. Analysis of a serendipitous Chandra X-ray gratings spectrum and light curve as well as XMM-Newton light curves and CCD spectra of GSC 07396-00759 obtained during long exposures targeting V4046 Sgr AB reveals a relatively hard (T_X ~ 10^7 K) X-ray spectrum, strong flaring, and relatively low-density plasma. These X-ray characteristics of GCS 07396--00759 are indicative of a high level of coronal activity, consistent with its apparent weak-lined T Tauri star status. Interactions between V4046 Sgr AB and GCS 07396-00759 when the two systems were more closely bound may be responsible for (a) their dissolution ~10^6 yr ago, (b) the present tight, circular orbit of V4046 Sgr AB, and (c) the persistence of the gaseous circumbinary disk still orbiting V4046 Sgr AB.
[63]  [pdf] - 1076238
Radio Astrometry Of The Triple Systems Algol And UX Arietis
Submitted: 2011-04-26
We have used multi-epoch long-baseline radio interferometry to determine the proper motion and orbital elements of Algol and UX Arietis, two radio-bright, close binary stellar systems with distant tertiary components. For Algol, we refine the proper motion and outer orbit solutions, confirming the recent result of Zavala et al. (2010) that the inner orbit is retrograde. The radio centroid closely tracks the motion of the KIV secondary. In addition, the radio morphology varies from double-lobed at low flux level to crescent-shaped during active periods. These results are most easily interpreted as synchrotron emission from a large, co-rotating meridional loop centered on the K-star. If this is correct, it provides a radio-optical frame tie candidate with an uncertainty {\pm}0.5 mas. For UX Arietis, we find a outer orbit solution that accounts for previous VLBI observations of an acceleration term in the proper motion fit. The outer orbit solution is also consistent with previously published radial velocity curves and speckle observations of a third body. The derived tertiary mass, 0.75 solar masses, is consistent with the K1 main-sequence star detected spectroscopically. The inner orbit solution favors radio emission from the active K0IV primary only. The radio morphology, consisting of a single, partially resolved emission region, may be associated with the persistent polar spot observed using Doppler imaging.
[64]  [pdf] - 1052654
Disentangling the Environment of the FU Orionis Candidate HBC 722 with Herschel
Comments: Accepted by ApJL 10 March 2011; 14 pages, 4 figures, 1 table
Submitted: 2011-03-10
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.
[65]  [pdf] - 1051542
Searching for gas emission lines in Spitzer Infrared Spectrograph (IRS) spectra of young stars in Taurus
Comments: 24 pages, 10 figures Accepted for publication in Astronomy and Astrophysics
Submitted: 2011-01-19
Our knowledge of circumstellar disks has traditionally been based on studies of dust. However, gas dominates the disk mass and its study is key to understand the star and planet formation process. Spitzer can access gas emission lines in the mid-infrared, providing new diagnostics of the physical conditions in accretion disks and outflows. We have studied the spectra of 64 pre-main-sequence stars in Taurus using Spitzer/IRS observations. We have detected H2 (17.03, 28.22 $\mu$m) emission in 6 objects, [Ne II] (12.81 $\mu$m) in 18 objects, and [Fe II] (17.93, 25.99 $\mu$m) in 7 objects. [Ne II] detections are found primarily in Class II objects. The luminosity of the [Ne II] line, is in general higher for objects known to drive jets than for those without known jets, but the two groups are not statistically distinguishable. We have searched for correlations between the line luminosities and different parameters related to the star-disk system. The [Ne II] luminosity is correlated with X-ray luminosity for Class II objects. The [NeII] luminosity is correlated with disk mass and accretion rate when the sample is divided into high and low accretors. We also find correlations between [NeII] luminosity and mid-IR continuum luminosity and with luminosity of the [O I] (6300 \AA) line, the latter being an outflow tracer. [Fe II] luminosity correlates with mass accretion rate. No correlations were found between H2 luminosity and several tested parameters. Our study reveals a general trend toward accretion-related phenomena as the origin of the gas emission lines. Shocks in jets and outflowing material are more likely to play a significant role than shocks in infalling material. The role of X-ray irradiation is less prominent but still present for [Ne II], in particular for Class II sources, the lack of correlation between [Fe II] and [Ne II] points toward different emitting mechanisms.
[66]  [pdf] - 955755
The Radio-X-ray Relation in Cool Stars: Are we headed toward a divorce?
Comments: 12 pages, Cool Stars 16 splinter session summary, to appear in conference proceedings
Submitted: 2010-12-07
This Cool Stars 16 splinter session was devoted to reviewing our current knowledge of correlated X-ray and radio emission from cool stars in order to prepare for new large radio observatories such as the EVLA. A key interest was to discuss why the X-ray and radio luminosities of some cool stars are in clear breach of a correlation that holds for other active stars, the so-called G\"udel-Benz relation. This article summarizes the contributions whereas the actual presentations can be accessed on the splinter website (
[67]  [pdf] - 1032599
First results of the Herschel Key Program 'Dust, Ice and Gas in Time': Dust and Gas Spectroscopy of HD 100546
Submitted: 2010-05-19
We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.
[68]  [pdf] - 1032479
Dust, Ice and Gas in Time (DIGIT) Herschel program first results: A full PACS-SED scan of the gas line emission in protostar DK Cha
Comments: 4 Page letter, To appear in A&A special issue on Herschel
Submitted: 2010-05-14
DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.
[69]  [pdf] - 1026102
Chandra Detects the Rare Oxygen-type Wolf-Rayet Star WR 142 and OB Stars in Berkeley 87
Comments: To appear in the Astrophysical Journal
Submitted: 2010-04-03
We present first results of a Chandra X-ray observation of the rare oxygen-type Wolf-Rayet star WR 142 (= Sand 5 = St 3) harbored in the young, heavily-obscured cluster Berkeley 87. Oxygen type WO stars are thought to be the most evolved of the WRs and progenitors of supernovae or gamma ray bursts. As part of an X-ray survey of supposedly single Wolf-Rayet stars, we observed WR 142 and the surrounding Berkeley 87 region with Chandra ACIS-I. We detect WR 142 as a faint, yet extremely hard X-ray source. Due to weak emission, its nature as a thermal or nonthermal emitter is unclear and thus we discuss several emission mechanisms. Additionally, we report seven detections and eight non-detections by Chandra of massive OB stars in Berkeley 87, two of which are bright yet soft X-ray sources whose spectra provide a dramatic contrast to the hard emission from WR 142.
[70]  [pdf] - 1018816
A multi-wavelength study of the young star V1118 Orionis in outburst
Comments: Accepted in A&A, Tables will be published online
Submitted: 2009-12-16
Abriged version for astroph: The young late-type star V1118 Orionis was in outburst from 2005 to 2006. We followed the outburst with optical and near-infrared photometry; the X-ray emission was further probed with observations taken with XMM-Newton and Chandra during and after the outburst. In addition, we obtained mid-infrared photometry and spectroscopy with Spitzer at the peak of the outburst and in the post-outburst phase. The spectral energy distribution of V1118 Ori varied significantly over the course of the outburst. The optical flux showed the largest variations, most likely due to enhanced emission by a hot spot. The latter dominated the optical and near-infrared emission at the peak of the outburst, while the disk emission dominated in the mid-infrared. The X-ray flux correlated with the optical and infrared fluxes, indicating that accretion affected the magnetically active corona and the stellar magnetosphere. The thermal structure of the corona was variable with some indication of a cooling of the coronal temperature in the early phase of the outburst with a gradual return to normal values. Color-color diagrams in the optical and infrared showed variations during the outburst, with no obvious signature of reddening due to circumstellar matter. Using MC realizations of star+disk+hotspot models to fit the SED in ``quiescence'' and at the peak of the outburst, we determined that the mass accretion rate varied from about 2.5E-7 Msun/yr to 1E-6 Msun/yr; in addition the fractional area of the hotspot increased significantly as well. The multi-wavelength study of the V1118 Ori outburst helped us to understand the variations in spectral energy distributions and demonstrated the interplay between the disk and the stellar magnetosphere in a young, strongly accreting star.
[71]  [pdf] - 315638
HDE 245059: A Weak-Lined T Tauri Binary Revealed by Chandra and Keck
Comments: 12 pages, 10 figures Accepted for publication in ApJ
Submitted: 2009-02-20
We present the Chandra High Energy Transmission Grating Spectrometer (HETGS) and Keck observations of HDE 245059, a young weak-lined T Tauri star (WTTS), member of the pre-main sequence group in the Lambda Orionis Cluster. Our high spatial resolution, near-infrared observations with Keck reveal that HDE 245059 a binary separated by 0.87". Based on this new information we have obtained an estimate of the masses of the binary components; 3M_{sun} and 2.5M_{sun} for the north and south components, respectively. We have estimated the age of the system to be ~2-3 Myr. We detect both components of the binary in the zeroth order Chandra image and in the grating spectra. Our fits to the spectrum of the binary have shown that the emission is dominated by a plasma between 8 and 15 MK, a soft component at 4 MK and a hard component at 50 MK are also detected. The value of the hydrogen column density was low, 8 x 10^{19} cm^{-2}, likely due to the clearing of the inner region of the Lambda Orionis cloud. The abundance pattern shows an inverse First Ionization Potential (FIP) effect for all elements from O to Fe, the only exception being Ca. A 3-T model was fitted to the individual zeroth order spectra using the abundances derived for the binary. We have also obtained several lines fluxes from the grating spectra. The fits to the triplets show no evidence of high densities. We conclude that the X-ray properties of the weak-lined T Tau binary HDE 245059 are similar to those generally observed in other weak-lined T Tau stars. Although its accretion history may have been affected by the clearing of the interstellar material around Lambda Ori, its coronal properties appears not to have been strongly modified.
[72]  [pdf] - 11979
Multiwavelength studies of the gas and dust disc of IRAS 04158+2805
Comments: 8 figures, 11 pages, accepted by A&A
Submitted: 2008-04-22
We present a study of the circumstellar environment of IRAS 04158+2805 based on multi-wavelength observations and models. Images in the optical and near-infrared, a polarisation map in the optical, and mid-infrared spectra were obtained with VLT-FORS1, CFHT-IR, and Spitzer-IRS. Additionally we used an X-ray spectrum observed with Chandra. We interpret the observations in terms of a central star surrounded by an axisymmetric circumstellar disc, but without an envelope, to test the validity of this simple geometry. We estimate the structural properties of the disc and its gas and dust content. We modelled the dust disc with a 3D continuum radiative transfer code, MCFOST, based on a Monte-Carlo method that provides synthetic scattered light images and polarisation maps, as well as spectral energy distributions. We find that the disc images and spectral energy distribution narrowly constrain many of the disc model parameters, such as a total dust mass of 1.0-1.75x10^-4 sollar masses and an inclination of 62-63 degrees. The maximum grain size required to fit all available data is of the order of 1.6-2.8 microns although the upper end of this range is loosely constrained. The observed optical polarisation map is reproduced well by the same disc model, suggesting that the geometry we find is adequate and the optical properties are representative of the visible dust content. We compare the inferred dust column density to the gas column density derived from the X-ray spectrum and find a gas-to-dust ratio along the line of sight that is consistent with the ISM value. To our knowledge, this measurement is the first to directly compare dust and gas column densities in a protoplanetary disc.
[73]  [pdf] - 6001
Scaling laws of solar and stellar flares
Comments: 9 Figs., (paper in press, The Astrophsycial Journal)
Submitted: 2007-10-12
In this study we compile for the first time comprehensive data sets of solar and stellar flare parameters, including flare peak temperatures T_p, flare peak volume emission measures EM_p, and flare durations t_f from both solar and stellar data, as well as flare length scales L from solar data. Key results are that both the solar and stellar data are consistent with a common scaling law of EM_p ~ T_p^4.7, but the stellar flares exhibit ~250 times higher emission measures (at the same flare peak temperature). For solar flares we observe also systematic trends for the flare length scale L(T_p) ~ T_p^0.9 and the flare duration t_F(T_p) ~ T_p^0.9 as a function of the flare peak temperature. Using the theoretical RTV scaling law and the fractal volume scaling observed for solar flares, i.e., V(L) ~ L^2.4, we predict a scaling law of EM_p ~ T_p^4.3, which is consistent with observations, and a scaling law for electron densities in flare loops, n_p ~ T_p^2/L ~ T_p^1.1. The RTV-predicted electron densities were also found to be consistent with densities inferred from total emission measures, n_p=(EM_p/q_V*V)^1/2, using volume filling factors of q_V=0.03-0.08 constrained by fractal dimensions measured in solar flares. Our results affect also the determination of radiative and conductive cooling times, thermal energies, and frequency distributions of solar and stellar flare energies.
[74]  [pdf] - 84785
Unbinned maximum-likelihood estimators for low-count data: Applications to faint X-ray spectra in the Taurus Molecular Cloud
Comments: accepted for publication in A&A
Submitted: 2006-09-07
Traditional binned statistics such as $\chi^2$ suffer from information loss and arbitrariness of the binning procedure. We point out that the underlying statistical quantity (the log likelihood $L$) does not require any binning beyond the one implied by instrumental readout channels, and we propose to use it for low-count data. The performance of $L$ in the model classification and point estimation problems is explored by Monte-Carlo simulations of Chandra and XMM X-ray spectra, and is compared to the performances of the binned Poisson statistic ($C$), Pearson's $\chi^2$ and Neyman's $\chi^2_N$, the Kolmogorov- Smirnov, and Kuiper' statistics. It is found that the unbinned log likelihood $L$ performs best with regard to the expected chi-square distance between true and estimated spectra, the chance of a successful identification among discrete candidate models, the area under the receiver-operator curve of reduced (two-model) binary classification problems, and generally also with regard to the mean square errors of individual spectrum parameters. The $\chi^2$ ($\chi^2_{\rm N}$) statistics should only be used if more than 10 (15) predicted counts per bin are available. From the practical point of view, the computational cost of evaluating $L$ is smaller than for any of the alternative methods if the forward model is specified in terms of a Poisson intensity and normalization is a free parameter. The maximum-$L$ method is applied to 14 observations from the Taurus Molecular Cloud, and the unbinned results are compared to binned XSPEC results, and found to generally agree, with exceptions explained by instability under re-binning and by background fine structures. The maximum-$L$ method has no lower limit on the available counts, and allows to treat weak sources which are beyond the means of binned methods.
[75]  [pdf] - 84619
A U-band survey of brown dwarfs in the Taurus Molecular Cloud with the XMM-Newton Optical/UV Monitor
Comments: 6 pages and 4 Figures. Accepted by A&A, to appear in a special section/issue dedicated to the XMM-Newton Extended Survey of the Taurus Molecular Cloud (XEST)
Submitted: 2006-09-01
We aim to characterize the U-band variability of young brown dwarfs in the Taurus Molecular Cloud and discuss its origin. We used the XMM-Newton Extended Survey of the Taurus Molecular Cloud, where a sample of 11 young bona fide brown dwarfs (spectral type later than M6) were observed simultaneously in X-rays with XMM-Newton and in the U-band with the XMM-Newton Optical/UV Monitor (OM). We obtained upper limits to the U-band emission of 10 brown dwarfs (U>19.6-20.6 mag), whereas 2MASSJ04141188+2811535 was detected in the U-band. Remarkably, the magnitude of this brown dwarf increased regularly from U~19.5 mag at the beginning of the observation, peaked 6h later at U~18.4 mag, and then decreased to U~18.65 mag in the next 2h. The first OM U-band measurement is consistent with the quiescent level observed about one year later thanks to ground follow-up observations. This brown dwarf was not detected in X-rays by XMM-Newton during the OM observation. We discuss the possible sources of U-band variability for this young brown dwarf, namely a magnetic flare, non-steady accretion onto the substellar surface, and rotational modulation of a hot spot. We conclude that this event is related to accretion from a circumsubstellar disk, where the mass accretion rate was about a factor of 3 higher than during the quiescent level.
[76]  [pdf] - 66073
On the sizes of stellar X-ray coronae
Comments: 22 pages, 13 figuges, 11 tables (A&A, in press)
Submitted: 2004-07-12, last modified: 2004-08-13
Spatial information from stellar X-ray coronae can be inferred by use of scaling laws from the solar corona. This requires measurement of the physical parameters temperature and density. We concentrate on the newly available density measurements from line fluxes of X-ray lines measured for a large sample of 64 grating spectra of 42 stellar coronae obtained with Chandra and XMM-Newton. Line counts of strong H-like and He-like ions and FeXXI lines were measured with the CORA line fitting tool by Ness and Wichmann 2002. Densities are estimated from He-like f/i flux ratios of OVII and NeIX representing only cooler (1-6MK) plasma components. The densities scatter between log n_e = 9.5-11 from the OVII triplet and between log n_e = 10.5-12 from the NeIX triplet, but we caution that the latter triplet may be biased by contamination from FeXIX and FeXXI lines. We find that low-activity stars tend to show densities of no more than a few times 10^10cm^-3, whereas no definitive trend is found for the more active stars. Investigating the densities of the hotter plasma with various FeXXI line ratios, we found no consistent proof for the presence of very high densities. This result is in line with constant pressure in the emitting active regions. We focus on the commonly used RTV78 scaling law to derive loop lengths from temperatures and densities assuming loop-like structures as identical building blocks. We derive emitting volumes from direct measurements of ion-specific emission measures and densities. Available volumes are calculated from the loop-lengths and stellar radii. The ratios of these volumes depict volume filling factors, which we find similar up to 0.1 for all stages of activity for OVII and increasing for NeIX.
[77]  [pdf] - 65045
Gamma-Ray Burst Polarization: Limits from RHESSI Measurements
Comments: 39 pages, 11 figures, accepted for publication by the Astrophysical Journal
Submitted: 2004-05-26, last modified: 2004-06-09
Using the RHESSI satellite as a Compton polarimeter, a recent study claimed that the prompt emission of GRB021206 was almost fully linearly polarized. This was challenged by a subsequent reanalysis. We present an novel approach, applying our method to the same data. We identify Compton scattering candidates by carefully filtering events in energy, time, and scattering geometry. Our polarization search is based on time dependent scattering rates in perpendicular directions, thus optimally excluding systematic errors. We perform simulations to obtain the instrument's polarimetric sensitivity, and these simulations include photon polarization. For GRB021206, we formally find a linear polarization degree of 41% (+57% -44%), concluding that the data quality is insufficient to constrain the polarization degree in this case. We further applied our analysis to GRB030519B and found again a null result.
[78]  [pdf] - 59695
Coronae of Cool Stars
Comments: To appear in IAU Symp. 219, ``Stars as Suns: Activity, Evolution, Planets'', ed. A. K. Dupree & A. O. Benz (San Francisco: ASP)
Submitted: 2003-10-01
We present preliminary results of grating observations of YY Mensae and V824 Arae by Chandra and XMM-Newton. Spectral features are presented in the context of the emission measure distributions, the coronal abundances, and plasma electron densities. In particular, we observe a coronal N/C enhancement in YY Men believed to reflect the photospheric composition (CN cycle). Finally, we interpret line broadening in YY Men as Doppler thermal broadening in its very hot corona.