BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

James E. Owen

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAastro-ph.HEastro-ph.IMphysics.space-ph

Catalog footprint

What is connected

42works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

42 published item(s)

preprint2024arXiv

On the likely magnesium-iron silicate dusty tails of catastrophically evaporating rocky planets

Catastrophically evaporating rocky planets provide a unique opportunity to study the composition of small planets. The surface composition of these planets can be constrained via modelling their comet-like tails of dust. In this work, we present a new self-consistent model of the dusty tails: we physically model the trajectory of the dust grains after they have left the gaseous outflow, including an on-the-fly calculation of the dust cloud's optical depth. We model two catastrophically evaporating planets: KIC 1255b and K2-22b. For both planets, we find the dust is likely composed of magnesium-iron silicates (olivine and pyroxene), consistent with an Earth-like composition. We constrain the initial dust grain sizes to be $\sim$ 1.25-1.75 $μ$m and the average (dusty) planetary mass-loss rate to be $\sim$ 3$M_\oplus \mathrm{Gyr^{-1}}$. Our model shows the origin of the leading tail of dust of K2-22b is likely a combination of the geometry of the outflow and a low radiation pressure force to stellar gravitational force ratio. We find the optical depth of the dust cloud to be a factor of a few in the vicinity of the planet. Our composition constraint supports the recently suggested idea that the dusty outflows of these planets go through a greenhouse effect-nuclear winter cycle, which gives origin to the observed transit depth time variability. Magnesium-iron silicates have the necessary visible-to-infrared opacity ratio to give origin to this cycle in the high mass-loss state.

0 citations0 reviewsNo code linkedOpen access
preprint2024arXiv

Using helium 10830 Å transits to constrain planetary magnetic fields

Planetary magnetic fields can affect the predicted mass loss rate for close-in planets that experience large amounts of UV irradiation. In this work, we present a method to detect the magnetic fields of close-in exoplanets undergoing atmospheric escape using transit spectroscopy at the 10830 Angstrom line of helium. Motivated by previous work on hydrodynamic and magneto-hydrodynamic photoevaporation, we suggest that planets with magnetic fields that are too weak to control the outflow's topology lead to blue-shifted transits due to day-to-night-side flows. In contrast, strong magnetic fields prevent this day-to-night flow, as the gas is forced to follow the magnetic field's roughly dipolar topology. We post-process existing 2D photoevaporation simulations to test this concept, computing synthetic transit profiles in helium. As expected, we find that hydrodynamically dominated outflows lead to blue-shifted transits on the order of the sound speed of the gas. Strong surface magnetic fields lead to unshifted or slightly red-shifted transit profiles. High-resolution observations can distinguish between these profiles; however, eccentricity uncertainties generally mean that we cannot conclusively say velocity shifts are due to the outflow for individual planets. The majority of helium observations are blue-shifted, which could be a tentative indication that close-in planets generally have surface dipole magnetic field strengths $\lesssim 0.1$ gauss. More 3D hydrodynamic and magneto-hydrodynamic are needed to confirm this conclusion robustly.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An APEX search for carbon emission from NGC 1977 proplyds

We used the Atacama Pathfinder Experiment (APEX) telescope to search for CI 1-0 (492.16GHz) emission towards 8 proplyds in NGC 1977, which is an FUV radiation environment two orders of magnitude weaker than that irradiating the Orion Nebular Cluster (ONC) proplyds. CI is expected to enable us to probe the wind launching region of externally photoevaporating discs. Of the 8 targets observed, no 3$σ$ detections of the CI line were made despite reaching sensitivities deeper than the anticipated requirement for detection from prior APEX CI observations of nearby discs and models of external photoevaporation of quite massive discs. By comparing both the proplyd mass loss rates and CI flux constraints with a large grid of external photoevaporation simulations, we determine that the non-detections are in fact fully consistent with the models if the proplyd discs are very low mass. Deeper observations in CI and probes of the disc mass with other tracers (e.g. in the continuum and CO) can test this. If such a test finds higher masses, this would imply carbon depletion in the outer disc, as has been proposed for other discs with surprisingly low CI fluxes, though more massive discs would also be incompatible with models that can explain the observed mass loss rates and CI non-detections. The expected remaining lifetimes of the proplyds are estimated to be similar to those of proplyds in the ONC at 0.1Myr. Rapid destruction of discs is therefore also a feature of common, intermediate UV environments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extreme Pebble Accretion in Ringed Protoplanetary Discs

Axisymmetric dust rings containing tens to hundreds of Earth masses of solids have been observed in protoplanetary discs with (sub-)millimetre imaging. Here, we investigate the growth of a planetary embryo in a massive (150M$_\oplus$) axisymmetric dust trap through dust and gas hydrodynamics simulations. When accounting for the accretion luminosity of the planetary embryo from pebble accretion, the thermal feedback on the surrounding gas leads to the formation of an anticyclonic vortex. Since the vortex forms at the location of the planet, this has significant consequences for the planet's growth: as dust drifts towards the pressure maximum at the centre of the vortex, which is initially co-located with the planet, a rapid accretion rate is achieved, in a distinct phase of ``vortex-assisted'' pebble accretion. Once the vortex separates from the planet due to interactions with the disc, it accumulates dust, shutting off accretion onto the planet. We find that this rapid accretion, mediated by the vortex, results in a planet containing $\approx$ 100M$_\oplus$ of solids. We follow the evolution of the vortex, as well as the efficiency with which dust grains accumulate at its pressure maximum as a function of their size, and investigate the consequences this has for the growth of the planet as well as the morphology of the protoplanetary disc. We speculate that this extreme formation scenario may be the origin of giant planets which are identified to be significantly enhanced in heavy elements.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

One year of AU Mic with HARPS: I -- measuring the masses of the two transiting planets

The system of two transiting Neptune-sized planets around the bright, young M-dwarf AU Mic provides a unique opportunity to test models of planet formation, early evolution, and star-planet interaction. However, the intense magnetic activity of the host star makes measuring the masses of the planets via the radial velocity (RV) method very challenging. We report on a 1-year, intensive monitoring campaign of the system using 91 observations with the HARPS spectrograph, allowing for detailed modelling of the $\sim 600\, \rm{m\,s^{-1}}$ peak-to-peak activity-induced RV variations. We used a multidimensional Gaussian Process framework to model these and the planetary signals simultaneously. We detect the latter with semi-amplitudes of $\rm{K_b} = 5.8 \pm 2.5\, \rm{m\,s^{-1}}$ and $\rm{K_c} = 8.5 \pm 2.5\, \rm{m\, s^{-1}}$, respectively. The resulting mass estimates, $\rm{M_b} = 11.7 \pm 5.0\, \rm{M_{\rm \oplus}}$ and $\rm{M_c} = 22.2 \pm 6.7\, \rm{M_{\rm \oplus}}$, suggest that planet b might be less dense, and planet c considerably denser than previously thought. These results are in tension with the current standard models of core-accretion. They suggest that both planets accreted a H/He envelope that is smaller than expected, and the trend between the two planets' envelope fractions is the opposite of what is predicted by theory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

TESS Hunt for Young and Maturing Exoplanets (THYME) VI: an 11 Myr giant planet transiting a very low-mass star in Lower Centaurus Crux

Mature super-Earths and sub-Neptunes are predicted to be $\simeq$Jovian radius when younger than 10 Myr. Thus, we expect to find 5-15$R_\oplus$ planets around young stars even if their older counterparts harbor none. We report the discovery and validation of TOI 1227 b, a $0.85\pm0.05R_J$ (9.5$R_\oplus$) planet transiting a very low-mass star ($0.170\pm0.015M_\odot$) every 27.4 days. TOI~1227's kinematics and strong lithium absorption confirm it is a member of a previously discovered sub-group in the Lower Centaurus Crux OB association, which we designate the Musca group. We derive an age of 11$\pm$2 Myr for Musca, based on lithium, rotation, and the color-magnitude diagram of Musca members. The TESS data and ground-based follow-up show a deep (2.5\%) transit. We use multiwavelength transit observations and radial velocities from the IGRINS spectrograph to validate the signal as planetary in nature, and we obtain an upper limit on the planet mass of $\simeq0.5 M_J$. Because such large planets are exceptionally rare around mature low-mass stars, we suggest that TOI 1227 b is still contracting and will eventually turn into one of the more common $<5R_\oplus$ planets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The California-Kepler Survey. X. The Radius Gap as a Function of Stellar Mass, Metallicity, and Age

In 2017, the California-Kepler Survey (CKS) published its first data release (DR1) of high-resolution optical spectra of 1305 planet hosts. Refined CKS planet radii revealed that small planets are bifurcated into two distinct populations: super-Earths (smaller than 1.5 $R_E$) and sub-Neptunes (between 2.0 and 4.0 $R_E$), with few planets in between (the "Radius Gap.") Several theoretical models of the Radius Gap predict variation with stellar mass, but testing these predictions are challenging with CKS DR1 due to its limited $M_\star$ range of 0.8-1.4 $M_\odot$. Here, we present CKS DR2 with 411 additional spectra and derived properties focusing on stars of 0.5-0.8 $M_\odot$. We found the Radius Gap follows $R_p \propto P^m$ with $m = -0.10 \pm 0.03$, consistent with predictions of XUV- and core-powered mass-loss mechanisms. We found no evidence that $m$ varies with $M_\star$. We observed a correlation between the average sub-Neptune size and $M_\star$. Over 0.5 to 1.4 $M_\odot$, the average sub-Neptune grows from 2.1 to 2.6 $R_E$, following $R_p \propto M_\star^α$ with $α= 0.25 \pm 0.03$. In contrast, there is no detectable change for super-Earths. These $M_\star$-$R_p$ trends suggests that protoplanetary disks can efficiently produce cores up to a threshold mass of $M_c$, which grows linearly with stellar mass according to $M_c \approx 10 M_E~(M_\star / M_\odot)$. There is no significant correlation between sub-Neptune size and stellar metallicity (over $-$0.5 to $+$0.5 dex) suggesting a weak relationship between planet envelope opacity and stellar metallicity. Finally, there is no significant variation in sub-Neptune size with stellar age (over 1 to 10 Gyr), which suggests that the majority of envelope contraction concludes after $\sim$1 Gyr.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The UV-SCOPE Mission: Ultraviolet Spectroscopic Characterization Of Planets and their Environments

UV-SCOPE is a mission concept to determine the causes of atmospheric mass loss in exoplanets, investigate the mechanisms driving aerosol formation in hot Jupiters, and study the influence of the stellar environment on atmospheric evolution and habitability. As part of these investigations, the mission will generate a broad-purpose legacy database of time-domain ultraviolet (UV) spectra for nearly 200 stars and planets. The observatory consists of a 60 cm, f/10 telescope paired to a long-slit spectrograph, yielding simultaneous, almost continuous coverage between 1203 Å and 4000 Å, with resolutions ranging from 6000 to 240. The efficient instrument provides throughputs > 4% (far-UV; FUV) and > 15% (near-UV; NUV), comparable to HST/COS and much better than HST/STIS, over the same spectral range. A key design feature is the LiF prism, which serves as a dispersive element and provides high throughput even after accounting for radiation degradation. The use of two delta-doped Electron-Multiplying CCD detectors with UV-optimized, single-layer anti-reflection coatings provides high quantum efficiency and low detector noise. From the Earth-Sun second Lagrangian point, UV-SCOPE will continuously observe planetary transits and stellar variability in the full FUV-to-NUV range, with negligible astrophysical background. All these features make UV-SCOPE the ideal instrument to study exoplanetary atmospheres and the impact of host stars on their planets. UV-SCOPE was proposed to NASA as a Medium Explorer (MidEx) mission for the 2021 Announcement of Opportunity. If approved, the observatory will be developed over a 5-year period. Its primary science mission takes 34 months to complete. The spacecraft carries enough fuel for 6 years of operations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Unveiling the Planet Population at Birth

The radius distribution of small, close-in exoplanets has recently been shown to be bimodal. The photoevaporation model predicted this bimodality. In the photoevaporation scenario, some planets are completely stripped of their primordial H/He atmospheres, whereas others retain them. Comparisons between the photoevaporation model and observed planetary populations have the power to unveil details of the planet population inaccessible by standard observations, such as the core mass distribution and core composition. In this work, we present a hierarchical inference analysis on the distribution of close-in exoplanets using forward-models of photoevaporation evolution. We use this model to constrain the planetary distributions for core composition, core mass and initial atmospheric mass fraction. We find that the core-mass distribution is peaked, with a peak-mass of $\sim 4$M$_\oplus$. The bulk core-composition is consistent with a rock/iron mixture that is ice-poor and ``Earth-like''; the spread in core-composition is found to be narrow ($\lesssim 16\%$ variation in iron-mass fraction at the 2$σ$ level) and consistent with zero. This result favours core formation in a water/ice poor environment. We find the majority of planets accreted a H/He envelope with a typical mass fraction of $\sim 4\%$; only a small fraction did not accrete large amounts of H/He and were ``born-rocky''. We find four-times as many super-Earths were formed through photoevaporation, as formed without a large H/He atmosphere. Finally, we find core-accretion theory over-predicts the amount of H/He cores would have accreted by a factor of $\sim 5$, pointing to additional mass-loss mechanisms (e.g. ``boil-off'') or modifications to core-accretion theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fingerprints of giant planets in the composition of solar twins

The Sun shows a $\sim 10$% depletion in refractory elements relative to nearby solar twins. It has been suggested that this depletion is a signpost of planet formation. The exoplanet statistics are now good enough to show that the origin of this depletion does not arise from the sequestration of refractory material inside the planets themselves. This conclusion arises because most sun-like stars host close-in planetary systems that are on average more massive than the Sun's. Using evolutionary models for the protoplanetary discs that surrounded the young Sun and solar twins we demonstrate that the origin of the depletion likely arises due to the trapping of dust exterior to the orbit of a forming giant planet. In this scenario a forming giant planet opens a gap in the gas disc, creating a pressure trap. If the planet forms early enough, while the disc is still massive, the planet can trap $\gtrsim 100 $M$_\oplus$ of dust exterior to its orbit, preventing the dust from accreting onto the star in contrast to the gas. Forming giant planets can create refractory depletions of $\sim 5-15$%, with the larger values occurring for initial conditions that favour giant planet formation (e.g. more massive discs, that live longer). The incidence of solar-twins that show refractory depletion matches both the occurrence of giant planets discovered in exoplanet surveys and "transition" discs that show similar depletion patterns in the material that is accreting onto the star.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation of Orion Fingers

"Orion fingers" are a system of dozens of bowshocks, with the wings of shocks pointing to a common system of origin, which is centered on a dynamically disintegrating system of several massive stars. The shock heads propagate with velocities of up to 300-400 km/s, but the formation and physical properties of the "bullets" leading the shocks are not known. Here we summarize two possible scenarios for the formation of the "bullets" and the resulting bowshocks ("fingers"). In the first scenario, bullets are self-gravitating, Jupiter-mass objects which were formed rapidly and then ejected during the strong dynamical interactions of massive stars and their disks. This scenario naturally explains the similar timescales for the outflow of bullets and for the dynamical interaction of the massive stars, but has some difficulty explaining the observed high velocities of the bullets. In the second scenario, bullets are formed via hydrodynamic instabilities in a massive, infrared-driven wind, naturally explaining the high velocities and the morphology of outflow, but the bullets are not required to be self-gravitating. The processes that created the Orion fingers are likely not unique to this particular star-forming region and may result in free-floating, high-velocity, core-less planets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Massive discs around low-mass stars

We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Snow-lines can be thermally unstable

Volatile species in protoplanetary discs can undergo a phase change from vapour to solid. These "snow-lines" can play vital roles in planet formation at all scales, from dust coagulation to planetary migration. In the outer regions of protoplanetary discs, the temperature profile is set by the absorption of reprocessed stellar light by the solids. Further, the temperature profile sets the distribution of solids through sublimation and condensation at various snow-lines. Hence the snow-line position depends on the temperature profile and vice-versa. We show that this coupling can be thermally unstable, such that a patch of the disc at a snow-line will produce either run-away sublimation or condensation. This thermal instability arises at moderate optical depths, where heating by absorption of reprocessed stellar light from the disc's atmosphere is optically thick, yet cooling is optically thin. Since volatiles in the solid phase drift much faster than volatiles in the vapour phase, this thermal instability results in a limit-cycle. The snow-line progressively moves in, condensing volatiles, before receding, as the volatiles sublimate. Using numerical simulations, we study the evolution of the CO snow-line. We find the CO snow-line is thermally unstable under typical disc conditions and evolves inwards from $\sim50$ to $\sim30$~AU on timescales from 1,000-10,000 years. The CO snow-line spends between $\sim 10\%-50\%$ of its time at smaller separations, where the exact value is sensitive to the total opacity and turbulent viscosity. The evolving snow-line also creates ring-like structures in the solid distribution interior to the snow-line. Multiple ring-like structures created by moving snow-lines could potentially explain the sub-structures seen in many {\it ALMA} images.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The nature and origins of sub-Neptune size planets

Planets intermediate in size between the Earth and Neptune, and orbiting closer to their host stars than Mercury does the Sun, are the most common type of planet revealed by exoplanet surveys over the last quarter century. Results from NASA's Kepler mission have revealed a bimodality in the radius distribution of these objects, with a relative underabundance of planets between 1.5 and 2.0 $R_{\oplus}$. This bimodality suggests that sub-Neptunes are mostly rocky planets that were born with primary atmospheres a few percent by mass accreted from the protoplanetary nebula. Planets above the radius gap were able to retain their atmospheres ("gas-rich super-Earths"), while planets below the radius gap lost their atmospheres and are stripped cores ("true super-Earths"). The mechanism that drives atmospheric loss for these planets remains an outstanding question, with photoevaporation and core-powered mass loss being the prime candidates. As with the mass-loss mechanism, there are two contenders for the origins of the solids in sub-Neptune planets: the migration model involves the growth and migration of embryos from beyond the ice line, while the drift model involves inward-drifting pebbles that coagulate to form planets close-in. Atmospheric studies have the potential to break degeneracies in interior structure models and place additional constraints on the origins of these planets. However, most atmospheric characterization efforts have been confounded by aerosols. Observations with upcoming facilities are expected to finally reveal the atmospheric compositions of these worlds, which are arguably the first fundamentally new type of planetary object identified from the study of exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Observational Anatomy of Externally Photoevaporating Planet-Forming Discs I: Atomic Carbon

We demonstrate the utility of CI as a tracer of photoevaporative winds that are being driven from discs by their ambient UV environment. Commonly observed CO lines only trace these winds in relatively weak UV environments and are otherwise dissociated in the wind at the intermediate to high UV fields that most young stars experience. However, CI traces unsubtle kinematic signatures of a wind in intermediate UV environments ($\sim1000$G$_0$) and can be used to place constraints on the kinematics and temperature of the wind. In CI position-velocity diagrams external photoevaporation results in velocities that are faster than those from Keplerian rotation alone, as well as emission from quadrants of position-velocity space in which there would be no Keplerian emission. This is independent of viewing angle because the wind has components that are perpendicular to the azimuthal rotation of the disc. At intermediate viewing angles ($\sim 30-60^\circ$) moment 1 maps also exhibit a twisted morphology over large scales (unlike other processes that result in twists, which are typically towards the inner disc). CI is readily observable with ALMA, which means that it is now possible to identify and characterise the effect of external photoevaporation on planet-forming discs in intermediate UV environments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chondrule Transport in Protoplanetary Disks

Chondrule formation remains one of the most elusive early Solar System events. Here, we take the novel approach of employing numerical simulations to investigate chondrule origin beyond purely cosmochemical methods. We model the transport of generically-produced chondrules and dust in a 1D viscous protoplanetary disk model, in order to constrain the chondrule formation events. For a single formation event we are able to match analytical predictions of the memory chondrule and dust populations retain of each other (complementarity), finding that a large mass accretion rate ($\gtrsim 10^{-7}$~M$_\odot$~yr$^{-1}$) allows for delays on the order of the disk's viscous timescale between chondrule formation and chondrite accretion. Further, we find older disks to be severely diminished of chondrules, with accretion rates $\lesssim 10^{-9}$~M$_\odot$~yr$^{-1}$ for nominal parameters. We then characterize the distribution of chondrule origins in both space and time, as functions of disk parameters and chondrule formation rates, in runs with continuous chondrule formation and both static and evolving disks. Our data suggest that these can account for the observed diversity between distinct chondrite classes, if some diversity in accretion time is allowed for.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Disk-fed giant planet formation

Massive giant planets, such as the ones being discovered by direct imaging surveys, likely experience the majority of their growth through a circumplanetary disc. We argue that the entropy of accreted material is determined by boundary layer processes, unlike the "cold-" or "hot-start" hypotheses usually invoked in the core accretion and direct collapse scenarios. A simple planetary evolution model illustrates how a wide range of radius and luminosity tracks become possible, depending on details of the accretion process. Specifically, the proto-planet evolves towards "hot-start" tracks if the scale-height of the boundary layer is $\gtrsim0.24$, a value not much larger than the scale-height of the circumplanetary disc. Understanding the luminosity and radii of young giant planets will thus require detailed models of circumplanetary accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Habitability of Terrestrial-Mass Planets in the HZ of M Dwarfs. I. H/He-Dominated Atmospheres

The ubiquity of M dwarfs, combined with the relative ease of detecting terrestrial-mass planets around them, has made them prime targets for finding and characterising planets in the "Habitable Zone" (HZ). However, Kepler finds that terrestrial-mass exoplanets are often born with voluminous H/He envelopes, comprising mass-fractions ($M_{env}/M_{core}$) $\gtrsim 1$%. If these planets retain such envelopes over Gyr timescales, they will not be "habitable" even within the HZ. Given the strong X-ray/UV fluxes of M dwarfs, we study whether sufficient envelope mass can be photoevaporated away for these planets to become habitable. We improve upon previous work by using hydrodynamic models that account for radiative cooling as well as the transition from hydrodynamic to ballistic escape. Adopting a template active M dwarf XUV spectrum, including stellar evolution, and considering both evaporation and thermal evolution, we show that: (1) the mass-loss is (considerably) lower than previous estimates that use an "energy-limited" formalism and ignore the transition to Jeans escape, (2) at the inner edge of the HZ, planets with core mass $\lesssim 0.9$ M$_\oplus$, can lose enough H/He to become habitable if their initial envelope mass-fraction is $\sim$1%, (3) at the outer edge of the HZ, evaporation cannot remove a $\sim$1% H/He envelope even from cores down to 0.8 M$_\oplus$. Thus, if planets form with bulky H/He envelopes, only those with low-mass cores may eventually be habitable. Cores $\gtrsim$ 1 M$_\oplus$, with $\gtrsim$1% natal H/He envelopes, will not be habitable in the HZ of M dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rapid radiative clearing of protoplanetary discs

The lack of observed transition discs with inner gas holes of radii greater than ~50AU implies that protoplanetary discs dispersed from the inside out must remove gas from the outer regions rapidly. We investigate the role of photoevaporation in the final clearing of gas from low mass discs with inner holes. In particular, we study the so-called "thermal sweeping" mechanism which results in rapid clearing of the disc. Thermal sweeping was originally thought to arise when the radial and vertical pressure scale lengths at the X-ray heated inner edge of the disc match. We demonstrate that this criterion is not fundamental. Rather, thermal sweeping occurs when the pressure maximum at the inner edge of the dust heated disc falls below the maximum possible pressure of X-ray heated gas (which depends on the local X-ray flux). We derive new critical peak volume and surface density estimates for rapid radiative clearing which, in general, result in rapid dispersal happening less readily than in previous estimates. This less efficient clearing of discs by X-ray driven thermal sweeping leaves open the issue of what mechanism can clear gas from the outer disc sufficiently quickly to explain the non-detection of cold gas around weak line T Tauri stars.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The initial physical conditions of Kepler-36 b & c

The Kepler planetary system consists of two exoplanets at similar separations (0.115 & 0.128 AU), which have dramatically different densities. The inner planet has a density consistent with an Earth-like composition, while the outer planet is extremely low-density, such that it must contain a voluminous H/He envelope. Such a density difference would pose a problem for any formation mechanism if their current densities were representative of their composition at formation. However, both planets are at close enough separations to have undergone significant evaporation in the past. We constrain the core-mass, core composition, initial envelope-mass, and initial cooling-time of each planet using evaporation models conditioned on their present-day masses and radii, as inferred from Kepler photometry and transit timing analysis. The inner planet is consistent with being an evaporatively stripped core, while the outer planet has retained some of its initial envelope due to its higher core-mass. Therefore, both planets could have had a similar formation pathway, with the inner planet having an initial envelope-mass fraction of $\lesssim 10\%$ and core-mass of $\sim4.4$ M$_\oplus$, while the outer had an initial envelope-mass fraction of order 15-30\% and core-mass $\sim7.3$ M$_\oplus$. Finally, our results indicate that the outer planet had a long ($\gtrsim30$ Myr) initial cooling-time, much longer than would naively be predicted from simple timescale arguments. The long initial cooling-time could be evidence for a dramatic early cooling episode such as the recently proposed "boil-off" process.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Water loss from Earth-sized planets in the habitable zones of ultracool dwarfs: Implications for the planets of TRAPPIST-1

Ultracool dwarfs (UCD; $T_{\rm eff}<\sim3000~$K) cool to settle on the main sequence after $\sim$1 Gyr. For brown dwarfs, this cooling never stops. Their habitable zone (HZ) thus sweeps inward at least during the first Gyr of their lives. Assuming they possess water, planets found in the HZ of UCDs have experienced a runaway greenhouse phase too hot for liquid water prior to entering the HZ. It has been proposed that such planets are desiccated by this hot early phase and enter the HZ as dry worlds. Here we model the water loss during this pre-HZ hot phase taking into account recent upper limits on the XUV emission of UCDs and using 1D radiation-hydrodynamic simulations. We address the whole range of UCDs but also focus on the planets recently found around the $0.08~M_\odot$ dwarf TRAPPIST-1. Despite assumptions maximizing the FUV-photolysis of water and the XUV-driven escape of hydrogen, we find that planets can retain significant amounts of water in the HZ of UCDs, with a sweet spot in the $0.04$-$0.06~M_\odot$ range. We also studied the TRAPPIST-1 system using observed constraints on the XUV-flux. We find that TRAPPIST-1b and c may have lost as much as 15 Earth Oceans and planet d -- which might be inside the HZ -- may have lost less than 1 Earth Ocean. Depending on their initial water contents, they could have enough water to remain habitable. TRAPPIST-1 planets are key targets for atmospheric characterization and could provide strong constraints on the water erosion around UCDs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Atmospheres of low-mass planets: the "boil-off"

We show that, for a low-mass planet that orbits its host star within a few tenths of an AU (like the majority of the {\it Kepler} planets), the atmosphere it was able to accumulate while embedded in the proto-planetary disk may not survive unscathed after the disk disperses. This gas envelope, if more massive than a few percent of the core (with a mass below $10 M_\oplus$), has a cooling time that is much longer than the time-scale on which the planet exits the disk. As such, it could not have contracted significantly from its original size, of order the Bondi radius. So a newly exposed proto-planet would be losing mass via a Parker wind that is catalyzed by the stellar continuum radiation. This represents an intermediate stage of mass-loss, occurring soon after the disc has dispersed, but before the EUV/X-ray driven photoevaporation becomes relevant. The surface mass-loss induces a mass movement within the envelope that advects internal heat outward. As a result, the planet atmosphere rapidly cools down and contracts, until it has reached a radius of order $0.1$ Bondi radius, at which time the mass-loss effectively shuts down. Within a million years after the disk disperses, we find a planet that has only about ten percent of its original envelope, and a Kelvin-Helmholtz time that is much longer than its actual age. We suggest that this "boil-off" process may be partially responsible for the lack of planets above a radius of $2.5 R_\oplus$ in the {\it Kepler} data, provided planet formation results in initial envelope masses of tens of percent.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hot Jupiter Breezes: Time-dependent Outflows from Extrasolar Planets

We explore the dynamics of magnetically controlled outflows from Hot Jupiters, where these flows are driven by UV heating from the central star. In these systems, some of the open field lines do not allow the flow to pass smoothly through the sonic point, so that steady-state solutions do not exist in general. This paper focuses on this type of magnetic field configuration, where the resulting flow becomes manifestly time-dependent. We consider the case of both steady heating and time-variable heating, and find the time scales for the corresponding time variations of the outflow. Because the flow cannot pass through the sonic transition, it remains subsonic and leads to so-called breeze solutions. One manifestation of the time variability is that the flow samples a collection of different breeze solutions over time, and the mass outflow rate varies in quasi-periodic fashion. Because the flow is subsonic, information can propagate inward from the outer boundary, which determines, in part, the time scale of the flow variability. This work finds the relationship between the outer boundary scale and the time scale of flow variations. In practice, the location of the outer boundary is set by the extent of the sphere of influence of the planet. The measured time variability can be used, in principle, to constrain the parameters of the system (e.g., the strengths of the surface magnetic fields).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The long-term evolution of photoevaporating transition discs with giant planets

Photo-evaporation and planet formation have both been proposed as mechanisms responsible for the creation of a transition disc. We have studied their combined effect through a suite of 2d simulations of protoplanetary discs undergoing X-ray photoevaporation with an embedded giant planet. In a previous work we explored how the formation of a giant planet triggers the dispersal of the inner disc by photo-evaporation at earlier times than what would have happened otherwise. This is particularly relevant for the observed transition discs with large holes and high mass accretion rates that cannot be explained by photo-evaporation alone. In this work we significantly expand the parameter space investigated by previous simulations. In addition, the updated model includes thermal sweeping, needed for studying the complete dispersal of the disc. After the removal of the inner disc the disc is a non accreting transition disc, an object that is rarely seen in observations. We assess the relative length of this phase, to understand if it is long lived enough to be found observationally. Depending on the parameters, especially on the X-ray luminosity of the star, we find that the fraction of time spent as a non-accretor greatly varies. We build a population synthesis model to compare with observations and find that in general thermal sweeping is not effective enough to destroy the outer disc, leaving many transition discs in a relatively long lived phase with a gas free hole, at odds with observations. We discuss the implications for transition disc evolution. In particular, we highlight the current lack of explanation for the missing non-accreting transition discs with large holes, which is a serious issue in the planet hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The origin and evolution of transition discs: successes, problems and open questions

Transition discs are protoplanetary discs that show evidence for large holes or wide gaps (with widths comparable to their radii) in their dust component. These discs could be giving us clues about the disc destruction mechanism or hints about the location and time-scales for the formation of planets. However, at the moment there remain key gaps in our theoretical understanding. The vast majority of transition discs are accreting onto their central stars, indicating that - at least close to the star - dust has been depleted from the gas by a very large amount. In this review, we discuss evidence for two distinct populations of transition discs: mm-faint - those with low mm-fluxes, small holes ($\lesssim$10AU) and low accretion rates (~1e-10-1e-9 Msun/yr) and mm-bright - discs with large mm-fluxes, large holes ($\gtrsim$20 AU) and high accretion rates ~1e-8 Msun/yr. MM-faint transition discs are consistent with what would naively be expected from a disc undergoing dispersal; however, mm-bright discs are not, and are likely to be rare and long-lived objects. We discuss the two commonly proposed mechanisms for creating transition discs: photoevaporation and planet-disc interactions, with a particular emphasis on how they would evolve in these models, comparing these predictions to the observed population. More theoretical work on explaining the lack of optically thick, non-accreting transition discs is required in both the photoevaporation and planetary hypothesis, before we can start to use transition discs to constrain models of planet formation. Finally, we suggest that the few discs with primordial looking SEDs, but serendipitously imaged, showing large cavities in the mm (e.g. MWC758 & WSB 60) may represent a hidden population of associated objects. Characterising and understanding how these objects fit into the overall paradigm may allow us to unravel the mystery of transition discs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Torque on an exoplanet from an anisotropic evaporative wind

Winds from short-period Earth and Neptune mass exoplanets, driven by high energy radiation from a young star, may evaporate a significant fraction of a planet's mass. If the momentum flux from the evaporative wind is not aligned with the planet/star axis, then it can exert a torque on the planet's orbit. Using steady-state one-dimensional evaporative wind models we estimate this torque using a lag angle that depends on the product of the speed of the planet's upper atmosphere and a flow timescale for the wind to reach its sonic radius. We also estimate the momentum flux from time-dependent one-dimensional hydrodynamical simulations. We find that only in a very narrow regime in planet radius, mass and stellar radiation flux is a wind capable of exerting a significant torque on the planet's orbit. Similar to the Yarkovsky effect, the wind causes the planet to drift outward if atmospheric circulation is prograde (super-rotating) and in the opposite direction if the circulation is retrograde. A close-in super Earth mass planet that loses a large fraction of its mass in a wind could drift a few percent of its semi-major axis. While this change is small, it places constraints on the evolution of resonant pairs such as Kepler 36 b and c.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

UV driven evaporation of close-in planets: energy-limited; recombination-limited and photon-limited flows

We have investigated the evaporation of close-in exoplanets irradiated by ionizing photons. We find that the properties of the flow are controlled by the ratio of the recombination time to the flow time-scale. When the recombination time-scale is short compared to the flow time-scale the the flow is in approximate local ionization equilibrium with a thin ionization front, where the photon mean free path is short compared to flow scale. In this "recombination limited" flow the mass-loss scales roughly with the square root of the incident flux. When the recombination time is long compared to the flow time-scale the ionization front becomes thick and encompasses the entire flow, with the mass-loss rate scaling linearly with flux. If the planet's potential is deep the flow is approximately "energy-limited"; however, if the planet's potential is shallow we identify a new limiting mass-loss regime, which we term "photon-limited". In this scenario the mass-loss rate is purely limited by the incoming flux of ionizing photons. We have developed a new numerical approach that takes into account the frequency dependence of the incoming ionizing spectrum and performed a large suite of 1D simulations to characterise UV driven mass-loss around low mass planets. We find the flow is "recombination-limited" at high fluxes but becomes "energy-limited" at low fluxes; however, the transition is broad occurring over several order of magnitude in flux. Finally, we point out the transitions between the different flow types does not occur at a single flux value, but depends on the planet's properties, with higher mass planets becoming "energy-limited" at lower fluxes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accreting planets as dust dams in `transition' discs

We investigate under what circumstances an embedded planet in a protoplanetary disc may sculpt the dust distribution such that it observationally presents as a `transition' disc. We concern ourselves with `transition' discs that have large holes ($\gtrsim 10$ AU) and high accretion rates ($\sim 10^{-9}-10^{-8}$ M$_\odot$ yr$^{-1}$). Particularly, those discs which photoevaporative models struggle to explain. Assuming the standard picture for how massive planets sculpt their parent discs, along with the observed accretion rates in `transition' discs, we find that the accretion luminosity from the forming planet is significant, and can dominate over the stellar luminosity at the gap edge. This planetary accretion luminosity can apply a significant radiation pressure to small ($s\lesssim 1μ$m) dust particles provided they are suitably decoupled from the gas. Secular evolution calculations that account for the evolution of the gas and dust components in a disc with an embedded, accreting planet, show that only with the addition of the radiation pressure can we explain the full observed characteristics of a `transition' disc (NIR dip in the SED, mm cavity and high accretion rate). At suitably high planet masses ($\gtrsim 3-4$ M$_J$), radiation pressure from the accreting planet is able to hold back the small dust particles, producing a heavily dust-depleted inner disc that is optically thin (vertically and radially) to Infra-Red radiation. We use our models to calculate synthetic observations and present a observational evolutionary scenario for a forming planet, sculpting its parent disc. The planet-disc system will present as a `transition' disc with a dip in the SED, only when the planet mass and planetary accretion rate is high enough. At other times it will present as a disc with a primordial SED, but with a cavity in the mm, as observed in a handful of protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Astro & cosmo-chemical consequences of accretion bursts I: the D/H ratio of water

The D/H ratio of water in protostellar systems is a result of both inheritance from the parent molecular cloud and isotopic exchange in the disc. A possibly widespread feature of disc evolution, ignored in previous studies, is accretion bursts (or FU Orionis outbursts), which may thermally process a large fraction of the water. One proposed underlying mechanism for FU Orionis outbursts relies on the presence of a magnetically dead zone. Here we examine the evolution of (D/H)$_{\rm water}$ in 1D simulations of a disc's evolution that include dead zones and infall from an envelope with given D/H ratio in the infalling water ($\sim 10^{-3}$), and compare the results with similar calculations without dead zones. We find that the accretion bursts result in a significantly lower (D/H)$_{\rm water}$ ratio and a more extended region (radius up to $\sim 1-3$ AU) where water is equilibrated with hydrogen gas (D/H=$2\times 10^{-5}$), when compared to burst-free models. Solar system constraints suggest that our solar nebula either experienced no accretion bursts and had a Schmidt number $\lesssim 0.2$ or had a Schmidt number closer to "nominal" values ($\sim 1$) and experienced several accretion bursts. Finally, future observations of (D/H)$_{\rm water}$ in protoplanetary discs will allow inferences about angular momentum properties of the disc during disc building and the role of accretion bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Importance of thermal diffusion in the gravo-magnetic limit cycle

We consider the role of thermal diffusion due to turbulence and radiation on accretion bursts that occur in protoplanetary discs which contain dead zones. Using 1D viscous disc models we show that diffusive radial transport of heat is important during the gravo-magnetic limit cycle, and can strongly modify the duration and frequency of accretion outbursts. When the Prandtl number is large - such that turbulent diffusion of heat is unimportant - radial radiative diffusion reduces the burst duration compared to models with no diffusive transport of heat. When the Prandtl number is small ($\lesssim 25$) we find that turbulent diffusion dominates the radial transport of heat, reducing the burst duration to $\lesssim 10^3$ years as well as increasing the burst frequency. Furthermore, inclusion of radial transport of heat extends the range of infall rates under which the disc undergoes accretion bursts from $10^{-8}$ to $10^{-6}$ M$_\odot$ yr$^{-1}$ with no diffusion, to $10^{-8}$ to $\gtrsim10^{-4}$ M$_\odot$ yr$^{-1}$ with radiative and strong turbulent diffusion. The relative roles of radiative and turbulent thermal diffusion are likely to vary during an accretion burst, but simple estimates suggest that the expected Prandtl numbers are of the order of 10 in protoplanetary discs, and hence that turbulent diffusion is likely to be an important process during accretion outbursts due to the gravo-magnetic limit cycle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetically controlled mass loss from extrasolar planets in close orbits

We consider the role magnetic fields play in guiding and controlling mass-loss via evaporative outflows from exoplanets that experience UV irradiation. First we present analytic results that account for planetary and stellar magnetic fields, along with mass-loss from both the star and planet. We then conduct series of numerical simulations for gas giant planets, and vary the planetary field strength, background stellar field strength, UV heating flux, and planet mass. These simulations show that the flow is magnetically controlled for moderate field strengths and even the highest UV fluxes, i.e., planetary surface fields $B_P\gtrsim 0.3$ gauss and fluxes $F_{UV}\sim10^{6}$ erg s$^{-1}$. We thus conclude that outflows from all hot Jupiters with moderate surface fields are magnetically controlled. The inclusion of magnetic fields highly suppresses outflow from the night-side of the planet. Only the magnetic field lines near the pole are open and allow outflow to occur. The fraction of open field lines depends sensitively on the strength (and geometry) of the background magnetic field from the star, along with the UV heating rate. The net effect of the magnetic field is to suppress the mass loss rate by (approximately) an order of magnitude. Finally, some open field lines do not allow the flow to pass smoothly through the sonic point; flow along these streamlines does not reach steady-state, resulting in time-variable mass-loss.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetically Controlled Outflows from Planets

Hot Jupiters can experience mass loss driven by heating from UV radiation from their host stars, and this flow is often controlled by magnetic fields. More specifically, near the planetry surface, the magnetic pressure dominates the ram pressure of the outflow by several orders of magnitude. After leaving the vicinity of the planet, the flow must connect onto the background environment provided by the stellar wind and the stellar magnetic field. This contribution considers magnetically controlled planetary outflows and extends previous work by comparing two different geometries for the background magnetic field provided by the star. In the first case, stellar field is assumed to retain the form of a dipole, which is anti-aligned with the dipole field of the planet. In the second case, the stellar outflow opens up the stellar magnetic field structure so that the background field at the location of the planet is perpendicular to the planetary dipole. Using numerical simulations, we consider the launch of the planetary wind with these field configurations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Snow-lines as probes of turbulent diffusion in protoplanetary discs

Sharp chemical discontinuities can occur in protoplanetary discs, particularly at `snow-lines' where a gas-phase species freezes out to form ice grains. Such sharp discontinuities will diffuse out due to the turbulence suspected to drive angular momentum transport in accretion discs. We demonstrate that the concentration gradient - in the vicinity of the snow-line - of a species present outside a snow-line but destroyed inside is strongly sensitive to the level of turbulent diffusion (provided the chemical and transport time-scales are decoupled) and provides a direct measurement of the radial `Schmidt number' (the ratio of the angular momentum transport to radial turbulent diffusion). Taking as an example the tracer species N$_2$H$^+$, which is expected to be destroyed inside the CO snow-line (as recently observed in TW Hya) we show that ALMA observations possess significant angular resolution to constrain the Schmidt number. Since different turbulent driving mechanisms predict different Schmidt numbers, a direct measurement of the Schmidt number in accretion discs would allow inferences about the nature of the turbulence to be made.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Characterising thermal sweeping: a rapid disc dispersal mechanism

(Abridged) We consider the properties of protoplanetary discs that are undergoing inside-out clearing by photoevaporation. In particular, we aim to characterise the conditions under which a protoplanetary disc may undergo `thermal sweeping', a rapid (< 1e4 years) disc destruction mechanism proposed to occur when a clearing disc reaches sufficiently low surface density at its inner edge and where the disc is unstable to runaway penetration by the X-rays. We use a large suite of 1D radiation-hydrodynamic simulations to probe the observable parameter space, which is unfeasible in higher dimensions. These models allow us to determine the surface density at which thermal sweeping will take over the disc's evolution and to evaluate this critical surface density as a function of X-ray luminosity, stellar mass and inner hole radius. We find that this critical surface density scales linearly with X-ray luminosity, increases with inner hole radius and decreases with stellar mass and we develop an analytic model that reproduces these results. This surface density criterion is then used to determine the evolutionary state of protoplanetary discs at the point that they become unstable to destruction by thermal sweeping. We find that transition discs created by photoevaporation will undergo thermal sweeping when their inner holes reach 20-40 AU, implying that transition discs with large holes and no accretion (which were previously a predicted outcome of the later stages of all flavours of photoevaporation model) will not form. We emphasise that the surface density criteria that we have developed apply to all situations where the disc develops an inner hole that is optically thin to X-rays. It thus applies not only to the case of holes originally created by photoevaporation but also to holes formed, for example, by the tidal influence of planets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kepler planets: a tale of evaporation

(Abridged) Inspired by the Kepler planet discoveries, we consider the thermal contraction of planets close to their parent star, under the influence of evaporation. The mass-loss rates are based on hydrodynamic models of evaporation that include both X-ray and EUV irradiation. We find that only low-mass planets with hydrogen envelopes are significantly affected by evaporation, with evaporation being able to remove massive hydrogen envelopes inward of 0.1 AU for Neptune-mass objects. We construct a theoretical population of planets with varying core masses, envelope masses, orbital separations, and stellar spectral types, and compare these against the sizes and densities measured for low-mass planets, both in the Kepler mission and from radial velocity surveys. This exercise leads us to conclude that evaporation is the driving force of evolution for close-in Kepler planets. In fact, some 50% of the Kepler planet candidates may have been significantly eroded. Evaporation explains two striking correlations observed in these objects: a lack of large radius/low density planets close to the stars, and a bimodal distribution in planet sizes with a deficit of planets around 2R_E. Planets that have experienced high X-ray exposures are generally smaller than this size, and those with lower X-ray exposures are typically larger. A bimodal distribution is naturally explained by the evaporation model, where, depending on their X-ray exposure, close-in planets can either hold on to hydrogen envelopes 1% in mass, or be stripped entirely. To quantitatively reproduce the observed features, we argue that not only do low-mass Kepler planets need to be made of rocky cores overlaid with hydrogen envelopes, but few of them should have initial masses above 20 M_E, and the majority of them should have core masses of a few Earth masses.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Testing protoplanetary disc dispersal with radio emission

We consider continuum free-free radio emission from the upper atmosphere of protoplanetary discs as a probe of the ionized luminosity impinging upon the disc. Making use of previously computed hydrodynamic models of disc photoevaporation within the framework of EUV and X-ray irradiation, we use radiative transfer post-processing techniques to predict the expected free-free emission from protoplanetary discs. In general, the free-free luminosity scales roughly linearly with ionizing luminosity in both EUV and X-ray driven scenarios, where the emission dominates over the dust tail of the disc and is partial optically thin at cm wavelengths. We perform a test observation of GM Aur at 14-18 Ghz and detect an excess of radio emission above the dust tail to a very high level of confidence. The observed flux density and spectral index are consistent with free-free emission from the ionized disc in either the EUV or X-ray driven scenario. Finally, we suggest a possible route to testing the EUV and X-ray driven dispersal model of protoplanetary discs, by combining observed free-free flux densities with measurements of mass-accretion rates. On the point of disc dispersal one would expect to find a M_dot^2 scaling with free-free flux in the case of EUV driven disc dispersal or a M_dot scaling in the case of X-ray driven disc dispersal.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The interplay between X-ray photoevaporation and planet formation

We assess the potential of planet formation instigating the early formation of a photoevaporation driven gap, up to radii larger than typical for photoevaporation alone. For our investigation we make use of hydrodynamics models of photoevaporating discs with a giant planet embedded. We find that, by reducing the mass accretion flow onto the star, discs that form giant planets will be dispersed at earlier times than discs without planets by X-ray photoevaporation. By clearing the portion of the disc inner of the planet orbital radius, planet formation induced photoevaporation (PIPE) is able to produce transition disc that for a given mass accretion rate have larger holes when compared to standard X-ray photoevaporation. This constitutes a possible route for the formation of the observed class of accreting transition discs with large holes, which are otherwise difficult to explain by planet formation or photoevaporation alone. Moreover, assuming that a planet is able to filter dust completely, PIPE produces a transition disc with a large hole and may provide a mechanism to quickly shut down accretion. This process appears to be too slow however to explain the observed desert in the population of transition disc with large holes and low mass accretion rates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planetary evaporation by UV & X-ray radiation: basic hydrodynamics

We consider the evaporation of close in planets by the star's intrinsic EUV and X-ray radiation. We calculate evaporation rates by solving the hydrodynamical problem for planetary evaporation including heating from both X-ray and EUV radiation. We show that most close-in planets (a<0.1 AU) are evaporating hydrodynamically, with the evaporation occurring in two distinct regimes: X-ray driven, in which the X-ray heated flow contains a sonic point, and EUV driven, in which the X-ray region is entirely sub-sonic. The mass-loss rates scale as L_X/a^2 for X-ray driven evaporation, and as Phi_*^{1/2}/a for EUV driven evaporation at early times, with mass-loss rates of order 10e10-10e14 g/s. No exact scaling exists for the mass-loss rate with planet mass and planet radius, however, in general evaporation proceeds more rapidly for planets with lower densities and higher masses. Furthermore, we find that in general the transition from X-ray driven to EUV driven evaporation occurs at lower X-ray luminosities for planets closer to their parent stars and for planets with lower densities. Coupling our evaporation models to the evolution of the high energy radiation - which falls with time - we are able to follow the evolution of evaporating planets. We find that most planets start off evaporating in the X-ray driven regime, but switch to EUV driven once the X-ray luminosity falls below a critical value. The evolution models suggest that while `hot Jupiters' are evaporating, they are not evaporating at a rate sufficient to remove the entire gaseous envelope on Gyr time-scales. However, we do find that close in Neptune mass planets are more susceptible to complete evaporation of their envelopes. Thus we conclude that planetary evaporation is more important for lower mass planets, particularly those in the `hot Neptune'/`super Earth' regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiative Transfer in Star Formation: Testing FLD and Hybrid Methods

We perform a comparison between two radiative transfer algorithms commonly employed in hydrodynamical calculations of star formation: grey flux limited diffusion and the hybrid scheme, in addition we compare these algorithms to results from the Monte-Carlo radiative transfer code MOCASSIN. In disc like density structures the hybrid scheme performs significantly better than the FLD method in the optically thin regions, with comparable results in optically thick regions. In the case of a forming high mass star we find the FLD method significantly underestimates the radiation pressure by a factor of ~100.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the theory of disc photoevaporation

We discuss a hydrodynamical model for the dispersal of protoplanetary discs around young, low mass (<1.5 M_sun) stars by photoevaporation from the central object's energetic radiation, which considers the far-ultraviolet as well as the X-ray component of the radiation field. We present analytical scaling relations and derive estimates for the total mass-loss rates, as well as discussing the existence of similarity solutions for flows from primordial discs and discs with inner holes. Furthermore, we perform numerical calculations, which span a wide range of parameter space and allow us to provide accurate scalings of the mass-loss rates with the physical parameters of the systems (X-ray and FUV luminosity, stellar mass, disc mass, disc temperature and inner hole radius). The model suggest that the X-ray component dominates the photoevaporative mass-loss rates from the inner disc. The mass-loss rates have values in the range from 10e-11 to 10e-7 M_sun/yr and scale linearly with X-ray luminosity, with only a weak dependence on the other parameters explored. However, in the case of high FUV to X-ray (L_FUV/L_X>100) luminosity ratios, the FUV constricts the X-ray flow and may dominate the mass-loss. Simulations of low mass discs with inner holes demonstrate a further stage of disc clearing, which we call `thermal sweeping'. This process occurs when the mid-plane pressure drops to sufficiently low values. At this stage a bound, warm, X-ray heated region becomes sufficiently large and unstable, such that the remaining disc material is cleared on approximately dynamical time-scales. This process significantly reduces the time taken to clear the outer regions of the disc, resulting in an expected transition disc population that will be dominated by accreting objects, as indicated by recent observations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The imprint of photoevaporation on edge-on discs

We have performed hydrodynamic and radiative transfer calculations of a photoevaporating disc around a Herbig Ae/Be star to determine the evolution and observational impact of dust entrained in the wind. We find that the wind selectively entrains grains of different sizes at different radii resulting in a dust population that varies spatially and increases with height above the disc at radii > 10 AU. This variable grain population results in a 'wingnut' morphology to the dust density distribution. We calculate images of this dust distribution at NIR wavelengths that also show a wingnut morphology at all wavelengths considered. We have also considered the contribution that crystalline dust grains will have in the wind and show that a photoevaporative wind can result in a significant crystallinity fraction at all radii, when the disc is edge-on. However, when the disc's photosphere is unobscured, a photoevaporative wind makes no contribution to the observable crystallinity fraction in the disc. Finally, we conclude that the analysis of extended emission around edge-on discs could provide a new and independent method of testing photoevaporation models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Theoretical spectra of photoevaporating protoplanetary discs: An atlas of atomic and low-ionisation emission lines

We present a calculation of the atomic and low-ionisation emission line spectra of photoevaporating protoplanetary discs. Line luminosities and profiles are obtained from detailed photoionisation calculations of the disc and wind structures surrounding young active solar-type stars. The disc and wind density and velocity fields were obtained from the recently developed radiation-hydrodynamic models of Owen et al., that include stellar X-ray and EUV irradiation of protoplanetary discs at various stages of clearing, from primordial sources to inner hole sources of various hole sizes. Our models compare favourably with currently available observations, lending support to an X-ray driven photoevaporation model for disc dispersal. In particular, we find that X-rays drive a warm, predominantly neutral flow where the OI 6300A line can be produced by neutral hydrogen collisional excitation. Our models can, for the first time, provide a very good match to both luminosities and profiles of the low-velocity component of the OI 6300A line and other forbidden lines observed by Hartigan et al., which covered a large sample of T-Tauri stars. We find that the OI 6300A and the NeII 12.8um lines are predominantly produced in the X-ray-driven wind and thus appear blue-shifted by a few km/s for some of the systems when observed at non-edge-on inclinations. We note however that blue-shifts are only produced under certain conditions: X-ray luminosity, spectral shape and inner hole size all affect the location of the emitting region and the physical conditions in the wind. We caution therefore that while a blueshifted line is a tell-tale sign of an outflow, the lack of a blueshift should not be necessarily interpreted as a lack of outflow.

0 citations0 reviewsNo code linkedOpen access