Showing 1–36 of 36 results for author: Cridland, A

ALMA high-resolution observations unveil planet formation shaping molecular emission in the PDS 70 disk

Authors: L. Rampinelli, S. Facchini, M. Leemker, J. Bae, M. Benisty, R. Teague, C. J. Law, K. I. Öberg, B. Portilla-Revelo, A. J. Cridland

Abstract: With two directly detected protoplanets, the PDS 70 system is a unique source in which to study the complex interplay between forming planets and their natal environment. The large dust cavity carved by the two giant planets can affect the disk chemistry, and therefore the molecular emission morphology. On the other hand, chemical properties of the gas component of the disk are expected to leave a… ▽ More With two directly detected protoplanets, the PDS 70 system is a unique source in which to study the complex interplay between forming planets and their natal environment. The large dust cavity carved by the two giant planets can affect the disk chemistry, and therefore the molecular emission morphology. On the other hand, chemical properties of the gas component of the disk are expected to leave an imprint on the planetary atmospheres. In this work, we reconstruct the emission morphology of a rich inventory of molecular tracers in the PDS 70 disk, and we look for possible chemical signatures of the two actively accreting protoplanets, PDS b and c. We leverage Atacama Large Millimeter/submillimeter Array (ALMA) band 6 high-angular-resolution and deep-sensitivity line emission observations, together with image and $uv$-plane techniques, to boost the detection of faint lines. We robustly detect ring-shaped emission from $^{12}$CO, $^{13}$CO, C$^{18}$O, H$^{13}$CN, HC$^{15}$N, DCN, H$_2$CO, CS, C$_2$H, and H$^{13}$CO$^{+}$ lines in unprecedented detail. Most of the molecular tracers show a peak of the emission inside the millimeter dust peak. We interpret this as the direct impact of the effective irradiation of the cavity wall, as a result of the planet formation process. Moreover, we have found evidence of an O-poor gas reservoir in the outer disk, which is supported by the observations of bright C-rich molecules, the non-detection of SO, and a lower limit on the $\mathrm{CS/SO}$ ratio of $\sim1$. Eventually, we provide the first detection of the c-C$_3$H$_2$ transitions at 218.73 GHz, and the marginal detection of an azimuthal asymmetry in the higher-energy H$_2$CO (3$_{2,1}$-2$_{2,0}$) line, which could be due to accretion heating near PDS 70b. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 25 pages, 10 figures, 3 tables, 8 figures and one table in appendix. Accepted for publication in A&A

High contrast at short separation with VLTI/GRAVITY: Bringing Gaia companions to light

Authors: N. Pourré, T. O. Winterhalder, J. -B. Le Bouquin, S. Lacour, A. Bidot, M. Nowak, A. -L. Maire, D. Mouillet, C. Babusiaux, J. Woillez, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube , et al. (151 additional authors not shown)

Abstract: Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working… ▽ More Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($Δ\mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($Δ\mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 16 pages, 14 figures. Submitted to A&A

Four-of-a-kind? Comprehensive atmospheric characterisation of the HR 8799 planets with VLTI/GRAVITY

Authors: E. Nasedkin, P. Mollière, S. Lacour, M. Nowak, L. Kreidberg, T. Stolker, J. J. Wang, W. O. Balmer, J. Kammerer, J. Shangguan, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli , et al. (73 additional authors not shown)

Abstract: With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels… ▽ More With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels of model flexibility to understand the temperature structure, chemistry and clouds of each planet using both petitRADTRANS atmospheric retrievals and fits to self-consistent radiative-convective equilibrium models. Using Bayesian Model Averaging to combine multiple retrievals, we find that the HR 8799 planets are highly enriched in metals, with [M/H] $\gtrsim$1, and have stellar to super-stellar C/O ratios. The C/O ratio increases with increasing separation from $0.55^{+0.12}_{-0.10}$ for d to $0.78^{+0.03}_{-0.04}$ for b, with the exception of the innermost planet which has a C/O ratio of $0.87\pm0.03$. By retrieving a quench pressure and using a disequilibrium chemistry model we derive vertical mixing strengths compatible with predictions for high-metallicity, self-luminous atmospheres. Bayesian evidence comparisons strongly favour the presence of HCN in HR 8799 c and e, as well as CH$_{4}$ in HR 8799 c, with detections at $>5σ$ confidence. All of the planets are cloudy, with no evidence for patchiness. The clouds of c, d and e are best fit by silicate clouds lying above a deep iron cloud layer, while the clouds of the cooler HR 8799 b are more likely composed of Na$_{2}$S. With well defined atmospheric properties, future exploration of this system is well positioned to unveil further detail in these planets, extending our understanding of the composition, structure, and formation history of these siblings. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 45 pages, 25 figures, submitted to Astronomy & Astrophysics

Combining Gaia and GRAVITY: Characterising five new Directly Detected Substellar Companions

Authors: T. O. Winterhalder, S. Lacour, A. Mérand, A. -L. Maire, J. Kammerer, T. Stolker, N. Pourré, C. Babusiaux, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (74 additional authors not shown)

Abstract: Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observati… ▽ More Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at an angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. This translates to a semi-major axis of (0.938 $\pm$ 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 $\pm$ 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia's upcoming DR4 and the advent of GRAVITY+. △ Less

Submitted 24 June, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

A catalogue of dual-field interferometric binary calibrators

Authors: M. Nowak, S. Lacour, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, B. Charnay, G. Chauvin, A. Chavez, E. Choquet, V. Christiaens, Y. Clénet, V. Coudé du Foresto, A. Cridland , et al. (75 additional authors not shown)

Abstract: Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to with… ▽ More Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for proper pointing of the instrument. Up until now, no list of properly vetted calibrators was available for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list, and make it available to the community. We identify a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the Southern sky. We observe them as part of a dedicated calibration programme, and determine whether these objects are true binaries (excluding higher multiplicities resolved interferometrically but unseen by imaging), and extract measurements of the separation vectors. We combine these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. We compile a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always at least two binary calibrators at airmass < 2 in the sky over the Paranal observatory, at any point in time. Any Principal Investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of "whereistheplanet" to predict the astrometry of these calibrators, which seamlessly integrates with "p2Gravity" for VLTI/GRAVITY dual-field observing material preparation. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages, 6 figures. Submitted to A&A

VLTI/GRAVITY Provides Evidence the Young, Substellar Companion HD 136164 Ab formed like a "Failed Star"

Authors: William O. Balmer, L. Pueyo, S. Lacour, J. J. Wang, T. Stolker, J. Kammerer, N. Pourré, M. Nowak, E. Rickman, S. Blunt, A. Sivaramakrishnan, D. Sing, K. Wagner, G. -D. Marleau, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, M. S. Bordoni , et al. (71 additional authors not shown)

Abstract: Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or c… ▽ More Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or core fragmentation)? Additionally, their orbits can yield model-independent dynamical masses, and when paired with wide wavelength coverage and accurate system age estimates, can constrain evolutionary models in a regime where the models have a wide dispersion depending on initial conditions. We present new interferometric observations of the $16\,\mathrm{Myr}$ substellar companion HD~136164~Ab (HIP~75056~Ab) with VLTI/GRAVITY and an updated orbit fit including proper motion measurements from the Hipparcos-Gaia Catalogue of Accelerations. We estimate a dynamical mass of $35\pm10\,\mathrm{M_J}$ ($q\sim0.02$), making HD~136164~Ab the youngest substellar companion with a dynamical mass estimate. The new mass and newly constrained orbital eccentricity ($e=0.44\pm0.03$) and separation ($22.5\pm1\,\mathrm{au}$) could indicate that the companion formed via the low-mass tail of the Initial Mass Function. Our atmospheric fit to the \texttt{SPHINX} M-dwarf model grid suggests a sub-solar C/O ratio of $0.45$, and $3\times$ solar metallicity, which could indicate formation in the circumstellar disk via disk fragmentation. Either way, the revised mass estimate likely excludes ``bottom-up" formation via core accretion in the circumstellar disk. HD~136164~Ab joins a select group of young substellar objects with dynamical mass estimates; epoch astrometry from future \textit{Gaia} data releases will constrain the dynamical mass of this crucial object further. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astronomical Journal. 9 figures, 3 tables

First VLTI/GRAVITY Observations of HIP 65426 b: Evidence for a Low or Moderate Orbital Eccentricity

Authors: S. Blunt, W. O. Balmer, J. J. Wang, S. Lacour, S. Petrus, G. Bourdarot, J. Kammerer, N. Pourré, E. Rickman, J. Shangguan, T. Winterhalder, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (73 additional authors not shown)

Abstract: Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging… ▽ More Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior-dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best fit ExoREM model with solar metallicity and C/O=0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available at https://github.com/sblunt/hip65426. △ Less

Submitted 6 October, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 28 pages, 14 figures. Revised and resubmitted to AJ

VLTI/GRAVITY Observations and Characterization of the Brown Dwarf Companion HD 72946 B

Authors: W. O. Balmer, L. Pueyo, T. Stolker, H. Reggiani, S. Lacour, A. -L. Maire, P. Mollière, M. Nowak, D. Sing, N. Pourré, S. Blunt, J. J. Wang, E. Rickman, Th. Henning, K. Ward-Duong, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet , et al. (74 additional authors not shown)

Abstract: Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. W… ▽ More Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using orbitize! and derive a strong dynamical mass constraint $\mathrm{M_B}=69.5\pm0.5~\mathrm{M_{Jup}}$ assuming a strong prior on the host star mass $\mathrm{M_A}=0.97\pm0.01~\mathrm{M_\odot}$ from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent BT-Settl-CIFIST model atmospheres, and perform atmospheric retrievals using petitRADTRANS. A dynamical mass prior only marginally influences the sampled distribution on effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure-temperature or cloud parameters within our current retrieval setup. Independent of cloud prescription and prior assumptions, we find agreement within $\pm2\,σ$ between the C/O ratio of the host ($0.52\pm0.05)$ and brown dwarf ($0.43$ to $0.63$), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high ($0.6-0.8$) in light of an excellent agreement of the data with the solar abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface-gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis. △ Less

Submitted 15 September, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astrophysical Journal. 32 pages, 14 figures, 7 tables. v2&3 correct errors in co-author's affiliations, figure rendering, and some grant acknowledgements

Planet formation in the PDS 70 system: Constraining the atmospheric chemistry of PDS 70b and c

Authors: A. J. Cridland, S. Facchini, E. F. van Dishoeck, M. Benisty

Abstract: Understanding the chemical link between protoplanetary disks and planetary atmospheres is complicated by the fact that the popular targets in the study of disks and planets are widely separated both in space and time. The 5 Myr PDS 70 systems offers a unique opportunity to directly compare the chemistry of a giant planet's atmosphere to the chemistry of its natal disk. To that end, we derive our c… ▽ More Understanding the chemical link between protoplanetary disks and planetary atmospheres is complicated by the fact that the popular targets in the study of disks and planets are widely separated both in space and time. The 5 Myr PDS 70 systems offers a unique opportunity to directly compare the chemistry of a giant planet's atmosphere to the chemistry of its natal disk. To that end, we derive our current best physical and chemical model for the PDS 70 disk through forward modelling of the $^{12}$CO, C$^{18}$O, and C$_2$H emission radial profiles with the thermochemical code DALI and find a volatile C/O ratio above unity in the outer disk. Using what we know of the PDS 70 disk today, we analytically estimate the properties of the disk as it was 4 Myr in the past when we assume that the giant planets started their formation, and compute a chemical model of the disk at that time. We compute the formation of PDS 70b and PDS 70c using the standard core accretion paradigm and account for the accretion of volatile and refractory sources of carbon and oxygen to estimate the resulting atmospheric carbon-to-oxygen number ratio (C/O) for these planets. Our inferred C/O ratio of the gas in the PDS 70 disk indicates that it is marginally carbon rich relative to the stellar C/O = 0.44 which we derive from an empirical relation between stellar metallicity and C/O. Under the assumption that the disk has been carbon rich for most of its lifetime, we find that the planets acquire a super-stellar C/O in their atmospheres. If the carbon-rich disk is a relatively recent phenomenon (i.e. developed after the formation of the planets at $\sim 1$ Myr) then the planets should have close to the stellar C/O in their atmospheres. This work lays the groundwork to better understand the disk in the PDS 70 system as well as the planet formation scenario that produce its planets. △ Less

Submitted 31 March, 2023; originally announced March 2023.

Comments: 18 pages, 14 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 674, A211 (2023)

Direct discovery of the inner exoplanet in the HD206893 system. Evidence for deuterium burning in a planetary-mass companion

Authors: S. Hinkley, S. Lacour, G. -D. Marleau, A. M. Lagrange, J. J. Wang, J. Kammerer, A. Cumming, M. Nowak, L. Rodet, T. Stolker, W. -O. Balmer, S. Ray, M. Bonnefoy, P. Mollière, C. Lazzoni, G. Kennedy, C. Mordasini, R. Abuter, S. Aigrain, A. Amorim, R. Asensio-Torres, C. Babusiaux, M. Benisty, J. -P. Berger, H. Beust , et al. (89 additional authors not shown)

Abstract: Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoi… ▽ More Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoing precision RV measurements with the HARPS spectrograph, as well as Gaia host star astrometry, we report a high significance detection of the companion HD206893c over three epochs, with clear evidence for Keplerian orbital motion. Our astrometry with $\sim$50-100 $μ$arcsec precision afforded by GRAVITY allows us to derive a dynamical mass of 12.7$^{+1.2}_{-1.0}$ M$_{\rm Jup}$ and an orbital separation of 3.53$^{+0.08}_{-0.06}$ au for HD206893c. Our fits to the orbits of both companions in the system utilize both Gaia astrometry and RVs to also provide a precise dynamical estimate of the previously uncertain mass of the B component, and therefore derive an age of $155\pm15$ Myr. We find that theoretical atmospheric/evolutionary models incorporating deuterium burning for HD206893c, parameterized by cloudy atmospheres provide a good simultaneous fit to the luminosity of both HD206893B and c. In addition to utilizing long-term RV information, this effort is an early example of a direct imaging discovery of a bona fide exoplanet that was guided in part with Gaia astrometry. Utilizing Gaia astrometry is expected to be one of the primary techniques going forward to identify and characterize additional directly imaged planets. Lastly, this discovery is another example of the power of optical interferometry to directly detect and characterize extrasolar planets where they form at ice-line orbital separations of 2-4\,au. △ Less

Submitted 3 April, 2023; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted to A&A

Early planet formation in embedded protostellar disks: Setting the stage for the first generation of planetesimals

Authors: A. J. Cridland, G. P. Rosotti, B. Tabone, L. Tychoniec, M. McClure, E. F. van Dishoeck

Abstract: (Abridged) Recent surveys of young star formation regions have shown that the average Class II object does not have enough dust mass to make the cores of giant planets. Younger Class 0/I objects have enough dust in their embedded disk, which begs the questions: can the first steps of planet formation occur in these younger systems? The first step is building the first planetesimals, generally beli… ▽ More (Abridged) Recent surveys of young star formation regions have shown that the average Class II object does not have enough dust mass to make the cores of giant planets. Younger Class 0/I objects have enough dust in their embedded disk, which begs the questions: can the first steps of planet formation occur in these younger systems? The first step is building the first planetesimals, generally believed to be the product of the streaming instability. Hence the question can be restated: are the physical conditions of embedded disks conducive to the growth of the streaming instability? Here we model the collapse of a `dusty' proto-stellar cloud to show that if there is sufficient drift between the falling gas and dust, regions of the embedded disk can become sufficiently enhanced in dust to drive the streaming instability. We include four models, three with different dust grain sizes and one with a different initial cloud angular momentum to test a variety of collapse trajectories. We find a `sweet spot' for planetesimal formation for grain sizes of a few 10s of micron since they fall sufficiently fast relative to the gas to build a high dust-to-gas ratio along the disk midplane, but have slow enough radial drift speeds in the embedded disk to maintain the high dust-to-gas ratio. Unlike the gas, which is held in hydrostatic equilibrium for a time due to gas pressure, the dust can begin collapsing from all radii at a much earlier time. The streaming instability can produce at least between 7-35 M$_\oplus$ of planetesimals in the Class 0/I phase of our smooth embedded disks, depending on the size of the falling dust grains. This first generation of planetesimals could represent the first step in planet formation, and occurs earlier in the lifetime of the young star than is traditionally thought. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A, 22 pages, 19 figures

Journal ref: A&A 662, A90 (2022)

MHD disc winds can reproduce fast disc dispersal and the correlation between accretion rate and disc mass in Lupus

Authors: Benoît Tabone, Giovanni P. Rosotti, Giuseppe Lodato, Philip J. Armitage, Alexander J. Cridland, Ewine F. van Dishoeck

Abstract: The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solut… ▽ More The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solutions of disc evolution and a disc population synthesis approach. We show that MHD disc winds alone can account for both disc dispersal and accretion properties. The decline of disc fraction over time is reproduced by assuming that the initial accretion timescale (a generalization of the viscous timescale) varies from disc to disc and that the decline of the magnetic field strength is slower than that of the gas. The correlation between accretion rate and disc mass, and the dispersion of the data around the mean trend as observed in Lupus is then naturally reproduced. The model also accounts for the rapidity of the disc dispersal. This paves the way for planet formation models in the paradigm of wind-driven accretion. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS Letters

Secular evolution of MHD wind-driven discs: analytical solutions in the expanded α-framework

Authors: Benoît Tabone, Giovanni P. Rosotti, Alexander J. Cridland, Philip J. Armitage, Giuseppe Lodato

Abstract: The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscou… ▽ More The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscous evolution for the MHD wind case. The transport of angular momentum and mass induced by the wind is parameterized by an $α$-like parameter and by the magnetic lever arm parameter $λ$. Extensions of the paradigmatic Lynden-Bell and Pringle similarity solutions to the wind case are presented. We show that wind-driven accretion leads to a steeper decrease in the disc mass and accretion rate than in viscous models due to the absence of disc spreading. If the decline of the magnetic field strength is slower than that of the gas surface density, the disc is dispersed after a finite time. The evolution of the disc in the $\dot{M}_*-M_D$ plane is sensitive to the wind and turbulence parameters. A disc population evolving under the action of winds can exhibit a correlation between $\dot{M}_*$ and $M_D$ depending on the initial conditions. The simplified framework proposed in this work opens to a new avenue to test the effectiveness of wind-driven accretion from the observed disc demographics and constitutes an important step to include wind-driven accretion in planet population synthesis models. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS

Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories

Authors: Jessica Speedie, Ralph E. Pudritz, Alex J. Cridland, Farzana Meru, Richard A. Booth

Abstract: ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed thro… ▽ More ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed throughout the disk, assuming either a high or low turbulent $α$ viscosity. We have discovered that there is a bifurcation in the migration history of forming planets as a consequence of varying the disk viscosity. In our high viscosity disks, inward migration prevails and yields compact planetary systems, tempered only by planet trapping at the water iceline around 5 au. In our lower viscosity models however, low mass planets can migrate outward to twice their initial orbital radii, driven by a radially extended region of strong outward-directed corotation torques located near the heat transition (where radiative heating of the disk by the star is comparable to viscous heating) - before eventually migrating inwards. We derive analytic expressions for the planet mass at which the corotation torque dominates, and find that this "corotation mass" scales as $M_{\rm p, corot} \sim α^{2/3}$. If disk winds dominate the corotation torque, the corotation mass scales linearly with wind strength. We propose that the observed bifurcation in disk demographics into a majority of compact dust disks and a minority of extended ring/gap systems is a consequence of a distribution of viscosity across the disk population. △ Less

Submitted 4 February, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: 25 pages, 15 figures, 10 movies (links embedded). Published in MNRAS

Journal ref: MNRAS, Volume 510, Issue 4, March 2022, Pages 6059-6084

The mass of Beta Pictoris c from Beta Pictoris b orbital motion

Authors: S. Lacour, J. J. Wang, L. Rodet, M. Nowak, J. Shangguan, H. Beust, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, S. Blunt, A. Boccaletti, A. Bohn, M. -L. Bolzer, M. Bonnefoy, H. Bonnet, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay, G. Chauvin, E. Choquet , et al. (74 additional authors not shown)

Abstract: We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was… ▽ More We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was also reobserved at a separation of 96\,mas, confirming the previous orbital estimations. From the astrometry of planet b only, we can (i) detect the presence of $β$ Pictoris c and (ii) constrain its mass to $10.04^{+4.53}_{-3.10}\,M_{\rm Jup}$. If one adds the astrometry of $β$ Pictoris c, the mass is narrowed down to $9.15^{+1.08}_{-1.06}\,M_{\rm Jup}$. The inclusion of radial velocity measurements does not affect the orbital parameters significantly, but it does slightly decrease the mass estimate to $8.89^{+0.75}_{-0.75}\,M_{\rm Jup}$. With a semimajor axis of $2.68\pm0.02$\,au, a period of $1221\pm15$ days, and an eccentricity of $0.32\pm0.02$, the orbital parameters of $β$ Pictoris c are now constrained as precisely as those of $β$ Pictoris b. The orbital configuration is compatible with a high-order mean-motion resonance (7:1). The impact of the resonance on the planets' dynamics would then be negligible with respect to the secular perturbations, which might have played an important role in the eccentricity excitation of the outer planet. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: Accepted in A&A Letter

Journal ref: A&A 654, L2 (2021)

GRAVITY K-band spectroscopy of HD 206893 B: brown dwarf or exoplanet

Authors: J. Kammerer, S. Lacour, T. Stolker, P. Mollière, D. K. Sing, E. Nasedkin, P. Kervella, J. J. Wang, K. Ward-Duong, M. Nowak, R. Abuter, A. Amorim, R. Asensio-Torres, M. Bauböck, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, A. Bohn, M. -L. Bolzer, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube , et al. (72 additional authors not shown)

Abstract: We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectra… ▽ More We aim to reveal the nature of the reddest known substellar companion HD 206893 B by studying its near-infrared colors and spectral morphology and by investigating its orbital motion. We fit atmospheric models for giant planets and brown dwarfs and perform spectral retrievals with petitRADTRANS and ATMO on the observed GRAVITY, SPHERE, and GPI spectra of HD 206893 B. To recover its unusual spectral features, we include additional extinction by high-altitude dust clouds made of enstatite grains in the atmospheric model fits. We also infer the orbital parameters of HD 206893 B by combining the $\sim 100~μ\text{as}$ precision astrometry from GRAVITY with data from the literature and constrain the mass and position of HD 206893 C based on the Gaia proper motion anomaly of the system. The extremely red color and the very shallow $1.4~μ\text{m}$ water absorption feature of HD 206893 B can be fit well with the adapted atmospheric models and spectral retrievals. Altogether, our analysis suggests an age of $\sim 3$-$300~\text{Myr}$ and a mass of $\sim 5$-$30~\text{M}_\text{Jup}$ for HD 206893 B, which is consistent with previous estimates but extends the parameter space to younger and lower-mass objects. The GRAVITY astrometry points to an eccentric orbit ($e = 0.29^{+0.06}_{-0.11}$) with a mutual inclination of $< 34.4~\text{deg}$ with respect to the debris disk of the system. While HD 206893 B could in principle be a planetary-mass companion, this possibility hinges on the unknown influence of the inner companion on the mass estimate of $10^{+5}_{-4}~\text{M}_\text{Jup}$ from radial velocity and Gaia as well as a relatively small but significant Argus moving group membership probability of $\sim 61\%$. However, we find that if the mass of HD 206893 B is $< 30~\text{M}_\text{Jup}$, then the inner companion HD 206893 C should have a mass between $\sim 8$-$15~\text{M}_\text{Jup}$. △ Less

Submitted 15 June, 2021; originally announced June 2021.

Comments: 28 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 652, A57 (2021)

The ExoGRAVITY project: using single mode interferometry to characterize exoplanets

Authors: S. Lacour, J. J. Wang, M. Nowak, L. Pueyo, F. Eisenhauer, A. -M. Lagrange, P. Mollière, R. Abuter, A. Amorim, R. Asensio-Torres, M. Bauböck, M. Benisty, J. P. Berger, H. Beust, S. Blunt, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay, G. Chauvin, E. Choquet , et al. (67 additional authors not shown)

Abstract: Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orb… ▽ More Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances. △ Less

Submitted 19 January, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

Comments: SPIE 2020, invited talk

Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

Authors: J. J. Wang, A. Vigan, S. Lacour, M. Nowak, T. Stolker, R. J. De Rosa, S. Ginzburg, P. Gao, R. Abuter, A. Amorim, R. Asensio-Torres, M. Baubck, M. Benisty, J. P. Berger, H. Beust, J. -L. Beuzit, S. Blunt, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (79 additional authors not shown)

Abstract: We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbi… ▽ More We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $μ$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b. △ Less

Submitted 3 February, 2021; v1 submitted 11 January, 2021; originally announced January 2021.

Comments: 30 pages, 9 figures, Accepted for publication in AJ

Direct confirmation of the radial-velocity planet $β$ Pic c

Authors: M. Nowak, S. Lacour, A. -M. Lagrange, P. Rubini, J. Wang, T. Stolker, A. Amorim, R. Asensio-Torres, M. Bauböck, M. Benisty, J. P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube, B. Charnay, E. Choquet, V. Christiaens, Y. Clénet, V. Coudé du Foresto, A. Cridland, P. T. de Zeeuw , et al. (68 additional authors not shown)

Abstract: Methods used to detect giant exoplanets can be broadly divided into two categories: indirect and direct. Indirect methods are more sensitive to planets with a small orbital period, whereas direct detection is more sensitive to planets orbiting at a large distance from their host star. %, and thus on long orbital period. This dichotomy makes it difficult to combine the two techniques on a single ta… ▽ More Methods used to detect giant exoplanets can be broadly divided into two categories: indirect and direct. Indirect methods are more sensitive to planets with a small orbital period, whereas direct detection is more sensitive to planets orbiting at a large distance from their host star. %, and thus on long orbital period. This dichotomy makes it difficult to combine the two techniques on a single target at once. Simultaneous measurements made by direct and indirect techniques offer the possibility of determining the mass and luminosity of planets and a method of testing formation models. Here, we aim to show how long-baseline interferometric observations guided by radial-velocity can be used in such a way. We observed the recently-discovered giant planet $β$ Pictoris c with GRAVITY, mounted on the Very Large Telescope Interferometer (VLTI). This study constitutes the first direct confirmation of a planet discovered through radial velocity. We find that the planet has a temperature of $T = 1250\pm50$\,K and a dynamical mass of $M = 8.2\pm0.8\,M_{\rm Jup}$. At $18.5\pm2.5$\,Myr, this puts $β$ Pic c close to a 'hot start' track, which is usually associated with formation via disk instability. Conversely, the planet orbits at a distance of 2.7\,au, which is too close for disk instability to occur. The low apparent magnitude ($M_{\rm K} = 14.3 \pm 0.1$) favours a core accretion scenario. We suggest that this apparent contradiction is a sign of hot core accretion, for example, due to the mass of the planetary core or the existence of a high-temperature accretion shock during formation. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 10 pages

Connecting planet formation and astrochemistry: C/O and N/O of warm giant planets and Jupiter-analogs

Authors: Alex J. Cridland, Ewine F. van Dishoeck, Matthew Alessi, Ralph E. Pudritz

Abstract: (Abridged) The chemical composition of planetary atmospheres has long been thought to store information regarding where and when a planet accretes its material. Predicting this chemical composition theoretically is a crucial step in linking observational studies to the underlying physics that govern planet formation. As a follow-up to a study of hot Jupiters in our previous work, we present a popu… ▽ More (Abridged) The chemical composition of planetary atmospheres has long been thought to store information regarding where and when a planet accretes its material. Predicting this chemical composition theoretically is a crucial step in linking observational studies to the underlying physics that govern planet formation. As a follow-up to a study of hot Jupiters in our previous work, we present a population of warm Jupiters (semi-major axis between 0.5-4 AU) extracted from the same planetesimal formation population synthesis model as used in our previous work. We compute the astrochemical evolution of the protoplanetary disks included in this population to predict the carbon-to-oxygen (C/O) and nitrogen-to-oxygen (N/O) ratio evolution of the disk gas, ice, and refractory sources, the accretion of which greatly impacts the resulting C/O and N/O in the atmosphere of giant planets. We confirm that the main sequence (between accreted solid mass and atmospheric C/O) we found previously is largely reproduced by the presented population of synthetic warm Jupiters. And as a result, the majority of the population fall along the empirically derived mass-metallicity relation when the natal disk has solar or lower metallicity. Planets forming from disks with high metallicity ([Fe/H] $>$ 0.1) result in more scatter in chemical properties which could explain some of the scatter found in the mass-metallicity relation. Combining predicted C/O and N/O ratios shows that Jupiter does not fall among our population of synthetic planets, suggesting that it likely did not form in the inner 5 AU of the solar system before proceeding into a Grand Tack. This result is consistent with recent analysis of the chemical composition of Jupiter's atmosphere which suggests that it accreted most of its heavy element abundance farther than tens of AU away from the Sun. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A, 18 pages, 15 figures

Journal ref: A&A 642, A229 (2020)

Tracing bulk elemental ratios in exoplanetary atmospheres with TiO chemistry

Authors: Vanesa Ramírez, Alex J. Cridland, Paul Mollière

Abstract: Knowing the bulk elemental abundances of an exoplanetary atmosphere is not an easy task, but it is crucial to understand the formation history of planets. The purpose of this work is to show that the observability of TiO features at optical wavelengths in the transmission spectra of hot Jupiter atmospheres is sensitive to the bulk chemical properties of the atmosphere. For that, we run a grid of c… ▽ More Knowing the bulk elemental abundances of an exoplanetary atmosphere is not an easy task, but it is crucial to understand the formation history of planets. The purpose of this work is to show that the observability of TiO features at optical wavelengths in the transmission spectra of hot Jupiter atmospheres is sensitive to the bulk chemical properties of the atmosphere. For that, we run a grid of chemical models which include TiO formation and destruction, for the ultra-hot Jupiter WASP-19b and an ultra-hot version of HD~209458b. We take into account non-equilibrium chemistry and changes in the temperature and pressure structure of these atmospheres caused by different C/O ratios. We calculate synthetic transmission spectra for these models, and study the relative strengths of TiO and \ce{H2O} features quantitatively. To compare with observations, we use a model independent metric for molecular abundances, $ΔZ_{\rm TiO-H_2O}/H_{\rm eq}$ that has been previously used in observational studies of exoplanetary atmospheres. We find that with this metric we can differentiate between different chemical models and place constraints on the atmosphere's bulk carbon and oxygen abundance. From chemical considerations we expected that the TiO abundance would depend on the bulk nitrogen, however found that any change to N/H did not result in changes to the resulting TiO. We apply our method to a set of known exoplanets that have been observed in the relevant optical wavelengths and find good agreement between low resolution observations and our model for WASP-121b, marginally good agreement with WASP-79b, WASP-76b, and WASP-19b, and poorer agreement with HD 209458b. Our method can be particularly helpful for indirect studies of the bulk abundances of carbon and oxygen. △ Less

Submitted 2 July, 2020; originally announced July 2020.

Comments: 14 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 641, A87 (2020)

Dust masses of young disks: constraining the initial solid reservoir for planet formation

Authors: Łukasz Tychoniec, Carlo F. Manara, Giovanni P. Rosotti, Ewine F. van Dishoeck, Alexander J. Cridland, Tien-Hao Hsieh, Nadia M. Murillo, Dominique Segura-Cox, Sierk E. van Terwisga, John J. Tobin

Abstract: In recent years evidence has been building that planet formation starts early, in the first $\sim$ 0.5 Myr. Studying the dust masses available in young disks enables understanding the origin of planetary systems since mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. We aim to determine if disks in the embedded stage… ▽ More In recent years evidence has been building that planet formation starts early, in the first $\sim$ 0.5 Myr. Studying the dust masses available in young disks enables understanding the origin of planetary systems since mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka-band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities. Using the DIANA opacity model including large grains, with a dust opacity value of $κ_{\rm 9\ mm}$ = 0.28 cm$^{2}$ g$^{-1}$, the median dust masses of the embedded disks in Perseus are 158 M$_\oplus$ for Class 0 and 52 M$_\oplus$ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M$_\oplus$ and 12 M$_\oplus$ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value $κ_{\rm 1.3mm}$ = 2.3 cm$^{2}$ g$^{-1}$. The dust masses of young Class 0 and I disks are larger by at least a factor of 10 and 3, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of $\sim$ 15%. Higher efficiency of $\sim$ 30% is necessary if the planet formation is set to start in Class I disks. △ Less

Submitted 5 June, 2020; v1 submitted 4 June, 2020; originally announced June 2020.

Comments: 16 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 640, A19 (2020)

Formation of Planetary Populations II: Effects of Initial Disk Size & Radial Dust Drift

Authors: Matthew Alessi, Ralph E. Pudritz, Alex J. Cridland

Abstract: Recent ALMA observations indicate that while a range of disk sizes exist, typical disk radii are small, and that radial dust drift affects the distribution of solids in disks. Here we explore the consequences of these features in planet population synthesis models. A key feature of our model is planet traps - barriers to otherwise rapid type-I migration of forming planets - for which we include th… ▽ More Recent ALMA observations indicate that while a range of disk sizes exist, typical disk radii are small, and that radial dust drift affects the distribution of solids in disks. Here we explore the consequences of these features in planet population synthesis models. A key feature of our model is planet traps - barriers to otherwise rapid type-I migration of forming planets - for which we include the ice line, heat transition, and outer edge of the dead zone. We find that the ice line plays a fundamental role in the formation of warm Jupiters. In particular, the ratio of super Earths to warm Jupiters formed at the ice line depend sensitively on the initial disk radius. Initial gas disk radii of $\sim$50 AU results in the largest super Earth populations, while both larger and smaller disk sizes result in the ice line producing more gas giants near 1 AU. This transition between typical planet class formed at the ice line at various disk radii confirms that planet formation is fundamentally linked to disk properties (in this case, disk size), and is a result that is only seen when dust evolution effects are included in our models. Additionally, we find that including radial dust drift results in the formation of more super Earths between 0.1 - 1 AU, having shorter orbital radii than those produced in models where dust evolution effects are not included. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Comments: 24 pages, 13 figures. Submitted to MNRAS; revised in response to referee

Vertical gas accretion impacts the carbon-to-oxygen ratio of gas giant atmospheres

Authors: Alex J. Cridland, Arthur D. Bosman, Ewine F. van Dishoeck

Abstract: Recent theoretical, numerical, and observational work have suggested that when a growing planet opens a gap in its disk the flow of gas into the gap is dominated by gas falling vertically from a height of at least one gas scale height. Our primary objective is to include, for the first time, the chemical impact that accreting gas above the midplane will have on the resulting C/O. We compute the ac… ▽ More Recent theoretical, numerical, and observational work have suggested that when a growing planet opens a gap in its disk the flow of gas into the gap is dominated by gas falling vertically from a height of at least one gas scale height. Our primary objective is to include, for the first time, the chemical impact that accreting gas above the midplane will have on the resulting C/O. We compute the accretion of gas onto planetary cores beginning at different disk radii and track the chemical composition of the gas and small icy grains to predict the resulting carbon-to-oxygen ratio (C/O) in their atmospheres. In our model, all of the planets which began their evolution inward of 60 AU open a gap in the gas disk, and hence are chemically affected by the vertically accreting gas. Two important conclusions follow from this vertical flow: (1) more oxygen rich icy dust grains become available for accretion onto the planetary atmosphere. (2) The chemical composition of the gas dominates the final C/O of planets in the inner ($<$ 20 AU) part of the disk. This implies that with the launch of the James Webb Space Telescope we can trace the disk material that sets the chemical composition of exoplanetary atmospheres. △ Less

Submitted 16 January, 2020; originally announced January 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 635, A68 (2020)

Jupiter formed as a pebble pile around the N$_2$ ice line

Authors: Arthur D. Bosman, Alex J. Cridland, Yamila Miguel

Abstract: The region around the H$_2$O ice line, due to its higher surface density, seems to be the ideal location to form planets. The core of Jupiter, as well as the cores of close in gas giants are thus thought to form in this region of the disk. Actually constraining the formation location of individual planets has proven to be difficult, however. We aim to use the Nitrogen abundance in Jupiter, which i… ▽ More The region around the H$_2$O ice line, due to its higher surface density, seems to be the ideal location to form planets. The core of Jupiter, as well as the cores of close in gas giants are thus thought to form in this region of the disk. Actually constraining the formation location of individual planets has proven to be difficult, however. We aim to use the Nitrogen abundance in Jupiter, which is around 4 times solar, in combination with \textit{Juno} constraints on the total mass of heavy elements in Jupiter, to narrow down its formation scenario. Different pathways of enrichment of Jupiter's atmosphere, such as the accretion of enriched gas, pebbles or planetesimals are considered and their implications for the oxygen abundance of Jupiter is discussed. The super solar Nitrogen abundance in Jupiter necessitates the accretion of extra N$_2$ from the proto-solar nebula. The only location of the disk that this can happen is outside, or just inside the N$_2$ ice line. These constraints favor a pebble accretion origin of Jupiter, both from the composition as well as from a planet formation perspective. We predict that Jupiter's oxygen abundance is between 3.6 and 4.5 times solar. △ Less

Submitted 25 November, 2019; originally announced November 2019.

Comments: 5 pages, 2 figures

Journal ref: A&A 632, L11 (2019)

Connecting planet formation and astrochemistry: A main sequence for C/O in hot-exoplanetary atmospheres

Authors: Alex J. Cridland, Ewine F. van Dishoeck, Matthew Alessi, Ralph E. Pudritz

Abstract: To understand the role that planet formation history has on the observable atmospheric carbon-to-oxygen ratio (C/O) we have produced a population of astrochemically evolving protoplanetary disks. Based on the parameters used in a pre-computed population of growing planets their combination allows us to trace the molecular abundances of the gas that is being collected into planetary atmospheres. We… ▽ More To understand the role that planet formation history has on the observable atmospheric carbon-to-oxygen ratio (C/O) we have produced a population of astrochemically evolving protoplanetary disks. Based on the parameters used in a pre-computed population of growing planets their combination allows us to trace the molecular abundances of the gas that is being collected into planetary atmospheres. We include atmospheric pollution of incoming (icy) planetesimals as well as the effect of refractory carbon erosion noted to exist in our own solar system. We find that the carbon and oxygen content of Neptune-mass planets are determined primarily through solid accretion and result in more oxygen-rich (by roughly two orders of magnitude) atmospheres than Hot-Jupiters, whose C/O are primarily determined by gas accretion. Generally we find a `main-sequence' between the fraction of planetary mass accreted through solid accretion and the resulting atmospheric C/O - with planets of higher solid accretion fraction having lower C/O. Hot-Jupiters whose atmospheres have been chemically characterized agree well with our population of planets, and our results suggest that Hot Jupiter formation typically begins near the water ice line. Lower mass hot-Neptunes are observed to be much more carbon-rich (with 0.33 $\lesssim$ C/O $\lesssim$ 1) than is found in our models (C/O $\sim 10^{-2}$), and suggest that some form of chemical processing may affect their observed C/O over the few Gyrs between formation and observation. Our population reproduces the general mass-metallicity trend of the solar system and qualitatively reproduces the C/O-metallicity anti-correlation that has been inferred for the population of characterized exoplanetary atmospheres. △ Less

Submitted 29 October, 2019; originally announced October 2019.

Comments: Resubmitted to A&A after the second round of referee comments, 20 pages, 13 figures, 3 tables

Journal ref: A&A 632, A63 (2019)

Diverse outcomes of planet formation and composition around low-mass stars and brown dwarfs

Authors: Y. Miguel, A. Cridland, C. W. Ormel, J. J. Fortney, S. Ida

Abstract: The detection of Earth-size exoplanets around low-mass stars -- in stars such as Proxima Centauri and TRAPPIST-1 -- provide an exceptional chance to improve our understanding of the formation of planets around M stars and brown dwarfs. We explore the formation of such planets with a population synthesis code based on a planetesimal-driven model previously used to study the formation of the Jovian… ▽ More The detection of Earth-size exoplanets around low-mass stars -- in stars such as Proxima Centauri and TRAPPIST-1 -- provide an exceptional chance to improve our understanding of the formation of planets around M stars and brown dwarfs. We explore the formation of such planets with a population synthesis code based on a planetesimal-driven model previously used to study the formation of the Jovian satellites. Because the discs have low mass and the stars are cool, the formation is an inefficient process that happens at short periods, generating compact planetary systems. Planets can be trapped in resonances and we follow the evolution of the planets after the gas has dissipated and they undergo orbit crossings and possible mergers. We find that formation of planets above Mars mass and in the planetesimal accretion scenario, is only possible around stars with masses $M_{\star} \ge 0.07 M_{sun}$ and discs of $M_{disc} \ge 10^{-2}~M_{sun} $. We find that planets above Earth-mass form around stars with masses larger than $0.15 M_{\oplus}$ while planets larger than $5 M_{\oplus}$do not form in our model, even not under the most optimal conditions (massive disc), showing that planets such as GJ 3512b form with another, more efficient mechanism. Our results show that the majority of planets form with a significant water fraction; that most of our synthetic planetary systems have 1, 2, or 3 planets, but those with 4, 5, 6, and 7 planets are also common, confirming that compact planetary systems with many planets should be a relatively common outcome of planet formation around small stars. △ Less

Submitted 18 December, 2019; v1 submitted 26 September, 2019; originally announced September 2019.

Comments: published in MNRAS

Journal ref: MNRAS 491, 1998-2009

Connecting planet formation and astrochemistry: Refractory carbon depletion leading to super-stellar C/O in giant planetary atmospheres

Authors: Alex J. Cridland, Christian Eistrup, Ewine F. van Dishoeck

Abstract: [Abridged] Combining a time-dependent astrochemical model with a model of planet formation and migration, we compute the carbon-to-oxygen ratio (C/O) of a range of planetary embryos starting their formation in the inner solar system (1-3 AU). The volatile and ice abundance of relevant carbon and oxygen bearing molecular species are determined through a complex chemical kinetic code which includes… ▽ More [Abridged] Combining a time-dependent astrochemical model with a model of planet formation and migration, we compute the carbon-to-oxygen ratio (C/O) of a range of planetary embryos starting their formation in the inner solar system (1-3 AU). The volatile and ice abundance of relevant carbon and oxygen bearing molecular species are determined through a complex chemical kinetic code which includes both gas and grain surface chemistry. This is combined with a model for the abundance of the refractory dust grains to compute the total carbon and oxygen abundance in the protoplanetary disk available for incorporation into a planetary atmosphere. We include the effects of the refractory carbon depletion that has been observed in our solar system, and posit two models that would put this missing carbon back into the gas phase. This excess gaseous carbon then becomes important in determining the final planetary C/O because the gas disk now becomes more carbon rich relative to oxygen (high gaseous C/O). One model, where the carbon excess is maintained throughout the lifetime of the disk results in Hot Jupiters that have super-stellar C/O. The other model deposits the excess carbon early in the disk life and allows it to advect with the bulk gas. In this model the excess carbon disappears into the host star within 0.8 Myr, returning the gas disk to its original (sub-stellar) C/O, so the Hot Jupiters all exclusively have sub-stellar C/O. This shows that while the solids will tend to be oxygen rich, Hot Jupiters can have super-stellar C/O if a carbon excess can be maintained by some chemical processing of the dust grains. Whether the carbon and oxygen content of the atmosphere was accreted primarily by gas or solid accretion is heavily dependent on the mass of the atmosphere and where in the disk the growing planet accreted. △ Less

Submitted 25 January, 2019; originally announced January 2019.

Comments: 13 pages, 7 figures, resubmitted to A&A after referee's comments

Journal ref: A&A 627, A127 (2019)

Physics of Planet Trapping with Applications to HL Tau

Authors: Alex J. Cridland, Ralph E. Pudritz, Matthew Alessi

Abstract: We explore planet formation in the HL Tau disk and possible origins of the prominent gaps and rings observed by ALMA. We investigate whether dust gaps are caused by dynamically trapped planetary embryos at the ice lines of abundant volatiles. The global properties of the HL Tau disk (total mass, size) at its current age are used to constrain an evolving analytic disk model describing its temperatu… ▽ More We explore planet formation in the HL Tau disk and possible origins of the prominent gaps and rings observed by ALMA. We investigate whether dust gaps are caused by dynamically trapped planetary embryos at the ice lines of abundant volatiles. The global properties of the HL Tau disk (total mass, size) at its current age are used to constrain an evolving analytic disk model describing its temperature and density profiles. By performing a detailed analysis of the planet-disk interaction for a planet near the water ice line including a rigorous treatment of the dust opacity, we confirm that water is sufficiently abundant (1.5x10^-4 molecules per H) to trap planets at its ice line due to an opacity transition. When the abundance of water is reduced by 50% planet trapping disappears. We extend our analysis to other planet traps: the heat transition, dead zone edge, and the CO_2 ice line and find similar trapping. The formation of planets via planetesimal accretion is computed for dynamically trapped embryos at the water ice line, dead zone, and heat transition. The end products orbit in the inner disk (R < 3 AU), unresolved by ALMA, with masses that range between sub-Earth to 5 Jupiter masses. While we find that the dust gaps correspond well with the radial positions of the CO_2 , CH_4 , and CO ice lines, the planetesimal accretion rates at these radii are too small to build large embryos within 1 Myr. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

Magnetically induced termination of giant planet formation

Authors: A. J. Cridland

Abstract: Here a physical model for terminating giant planet formation is outlined and compared to other methods of late-stage giant planet formation. As has been pointed out before, gas accreting into a gap and onto the planet will encounter the planetary dynamo-generated magnetic field. The planetary magnetic field produces an effective cross section through which gas is accreted. Gas outside this cross s… ▽ More Here a physical model for terminating giant planet formation is outlined and compared to other methods of late-stage giant planet formation. As has been pointed out before, gas accreting into a gap and onto the planet will encounter the planetary dynamo-generated magnetic field. The planetary magnetic field produces an effective cross section through which gas is accreted. Gas outside this cross section is recycled into the protoplanetary disk, hence only a fraction of mass that is accreted into the gap remains bound to the planet. This cross section inversely scales with the planetary mass, which naturally leads to stalled planetary growth late in the formation process. We show that this method naturally leads to Jupiter-mass planets and does not invoke any artificial truncation of gas accretion, as has been done in some previous population synthesis models. The mass accretion rate depends on the radius of the growing planet after the gap has opened, and we show that so-called hot-start planets tend to become more massive than cold-start planets. When this result is combined with population synthesis models, it might show observable signatures of cold-start versus hot-start planets in the exoplanet population. △ Less

Submitted 20 September, 2018; v1 submitted 12 September, 2018; originally announced September 2018.

Comments: 8 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 619, A165 (2018)

Connecting Planetary Composition with Formation

Authors: Ralph E. Pudritz, Alex J. Cridland, Matthew Alessi

Abstract: The rapid advances in observations of the different populations of exoplanets, the characterization of their host stars and the links to the properties of their planetary systems, the detailed studies of protoplanetary disks, and the experimental study of the interiors and composition of the massive planets in our solar system provide a firm basis for the next big question in planet formation theo… ▽ More The rapid advances in observations of the different populations of exoplanets, the characterization of their host stars and the links to the properties of their planetary systems, the detailed studies of protoplanetary disks, and the experimental study of the interiors and composition of the massive planets in our solar system provide a firm basis for the next big question in planet formation theory. How do the elemental and chemical compositions of planets connect with their formation? The answer to this requires that the various pieces of planet formation theory be linked together in an end-to-end picture that is capable of addressing these large data sets. In this review, we discuss the critical elements of such a picture and how they affect the chemical and elemental make up of forming planets. Important issues here include the initial state of forming and evolving disks, chemical and dust processes within them, the migration of planets and the importance of planet traps, the nature of angular momentum transport processes involving turbulence and/or MHD disk winds, planet formation theory, and advanced treatments of disk astrochemistry. All of these issues affect, and are affected by the chemistry of disks which is driven by X-ray ionization of the host stars. We discuss how these processes lead to a coherent end-to-end model and how this may address the basic question. △ Less

Submitted 11 April, 2018; originally announced April 2018.

Comments: Invited review, accepted for publication in the 'Handbook of Exoplanets', eds. H.J. Deeg and J.A. Belmonte, Springer (2018). 46 pages, 10 figures

Composition of Early Planetary Atmospheres II: Coupled Dust and Chemical Evolution in Protoplanetary Disks

Authors: A. J. Cridland, Ralph E. Pudritz, Tilman Birnstiel, L. Ilsedore Cleeves, Edwin A. Bergin

Abstract: We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter… ▽ More We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter and a Super-Earth. We study the dependency of the final orbital radius, mass, and atmospheric chemistry of planets forming in disk models with initial disk masses that vary by 0.02 $M_\odot$ above and below our fiducial model ($M_{disk,0} = 0.1 ~M_\odot$). We compute C/O and C/N for the atmospheres formed in our 3 models and find that C/O$_{\rm planet}\sim$ C/O$_{\rm disk}$, which does not vary strongly between different planets formed in our model. The nitrogen content of atmospheres can vary in planets that grow in different disk models. These differences are related to the formation history of the planet, the time and location that the planet accretes its atmosphere, and are encoded in the bulk abundance of NH$_3$. These results suggest that future observations of atmospheric NH$_3$ and an estimation of the planetary C/O and C/N can inform the formation history of particular planetary systems. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: Accepted for publication in MNRAS

Radial Drift of Dust in Protoplanetary Disks: The Evolution of Ice lines and Dead zones

Authors: Alex J. Cridland, Ralph E. Pudritz, Tilman Birnstiel

Abstract: We have developed a new model for the astrochemical structure of a viscously evolving protoplanetary disk that couples an analytic description of the disk's temperature and density profile, chemical evolution, and an evolving dust distribution. We compute evolving radial distributions for a range of dust grain sizes, which depend on coagulation, fragmentation and radial drift processes. In particu… ▽ More We have developed a new model for the astrochemical structure of a viscously evolving protoplanetary disk that couples an analytic description of the disk's temperature and density profile, chemical evolution, and an evolving dust distribution. We compute evolving radial distributions for a range of dust grain sizes, which depend on coagulation, fragmentation and radial drift processes. In particular we find that the water ice line plays an important role in shaping the radial distribution of the maximum grain size because ice coated grains are significantly less susceptible to fragmentation than their dry counterparts. This in turn has important effects on disk ionization and therefore on the location of dead zones. In comparison to a simple constant gas-to-dust ratio model for the dust as an example, we find that the new model predicts an outer dead zone edge that moves in by a factor of about 3 at 1 Myr (to 5 AU) and by a factor of about 14 by 3 Myr (to 0.5 AU). We show that the changing position of the dead zone and heat transition traps have important implications for the formation and trapping of planets in protoplanetary disks. Finally, we consider our results in light of recent ALMA observations of HL Tau and TW Hya. △ Less

Submitted 8 November, 2016; originally announced November 2016.

Comments: 16 pages, 24 figures, submitted to MNRAS after one round of revisions

On the Formation and Chemical Composition of Super Earths

Authors: Matthew Alessi, Ralph E. Pudritz, Alex J. Cridland

Abstract: Super Earths are the largest population of exoplanets and are seen to exhibit a rich diversity of compositions as inferred through their mean densities. Here we present a model that combines equilibrium chemistry in evolving disks with core accretion that tracks materials accreted onto planets during their formation. In doing so, we aim to explain why super Earths form so frequently and how they a… ▽ More Super Earths are the largest population of exoplanets and are seen to exhibit a rich diversity of compositions as inferred through their mean densities. Here we present a model that combines equilibrium chemistry in evolving disks with core accretion that tracks materials accreted onto planets during their formation. In doing so, we aim to explain why super Earths form so frequently and how they acquire such a diverse range of compositions. A key feature of our model is disk inhomogeneities, or planet traps, that act as barriers to rapid type-I migration. The traps we include are the dead zone, which can be caused by either cosmic ray or X-ray ionization, the ice line, and the heat transition. We find that in disks with sufficiently long lifetimes ($\gtrsim$ 4 Myr), all traps produce Jovian planets. In these disks, planet formation in the heat transition and X-ray dead zone produces hot Jupiters while the ice line and cosmic ray dead zones produce Jupiters at roughly 1 AU. Super Earth formation takes place within short-lived disks ($\lesssim$ 2 Myr), whereby the disks are photoevaporated while planets are in a slow phase of gas accretion. We find that super Earth compositions range from dry and rocky ($<$ 6 % ice by mass) to those with substantial water contents ($>$ 30 % ice by mass). The traps play a crucial role in our results, as they dictate where in the disk particular planets can accrete from, and what compositions they are able to acquire. △ Less

Submitted 13 September, 2016; v1 submitted 29 June, 2016; originally announced June 2016.

Comments: 28 pages, 15 figures. Accepted for publication in MNRAS

Composition of Early Planetary Atmospheres I: Connecting Disk Astrochemistry to the Formation of Planetary Atmospheres

Authors: Alex J. Cridland, Ralph E. Pudritz, Matthew Alessi

Abstract: We present a model of the early chemical composition and elemental abundances of planetary atmospheres based on the cumulative gaseous chemical species that are accreted onto planets forming by core accretion from evolving protoplanetary disks. The astrochemistry of the host disk is computed using an ionization driven, non-equilibrium chemistry network within viscously evolving disk models. We acc… ▽ More We present a model of the early chemical composition and elemental abundances of planetary atmospheres based on the cumulative gaseous chemical species that are accreted onto planets forming by core accretion from evolving protoplanetary disks. The astrochemistry of the host disk is computed using an ionization driven, non-equilibrium chemistry network within viscously evolving disk models. We accrete gas giant planets whose orbital evolution is controlled by planet traps using the standard core accretion model and track the chemical composition of the material that is accreted onto the protoplanet. We choose a fiducial disk model and evolve planets in 3 traps - water ice line, dead zone and heat transition. For a disk with a lifetime of 4.1 Myr we produce two Hot Jupiters (M = 1.43, 2.67 M$_{\rm Jupiter}$, r = 0.15, 0.11 AU) in the heat transition and ice line trap and one failed core (M = 0.003 M$_{\rm Jupiter}$, r =3.7 AU) in the dead zone. These planets are found with mixing ratios for CO and H$_2$O of $1.99\times 10^{-4}$, $5.0\times 10^{-4}$ respectively for both Hot Jupiters. Additionally for these planets we find CO$_2$ and CH$_4$, with mixing ratios of $1.8\times 10^{-6}\rightarrow 9.8\times 10^{-10}$ and $1.1\times 10^{-8}\rightarrow 2.3\times 10^{-10}$ respectively. These ranges correspond well with the mixing ratio ranges that have been inferred through the detection of emission spectra from Hot Jupiters by multiple authors. We compute a carbon-to-oxygen ratio of 0.227 for the ice line planet and 0.279 for the heat transition planet. These planets accreted their gas inside the ice line, hence the sub-solar C/O. △ Less

Submitted 21 June, 2016; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: 25 pages, 16 figures, accepted for publication in MNRAS

Cold dust but warm gas in the unusual elliptical galaxy NGC 4125

Authors: C. D. Wilson, A. Cridland, K. Foyle, T. J. Parkin, E. Mentuch Cooper, H. Roussel, M. Sauvage, M. W. L. Smith, M. Baes, G. Bendo, M. Boquien, A. Boselli, L. Ciesla, D. L. Clements, A. Cooray, I. De Looze, M. Galametz, W. Gear, V. Lebouteiller, S. Madden, M. Pereira-Santaella, A. Remy-Ruyer

Abstract: Data from the Herschel Space Observatory have revealed an unusual elliptical galaxy, NGC 4125, which has strong and extended submillimeter emission from cold dust but only very strict upper limits to its CO and HI emission. Depending on the dust emissivity, the total dust mass is 2-5x10^6 Msun. While the neutral gas-to-dust mass ratio is extremely low (< 12-30), including the ionized gas traced by… ▽ More Data from the Herschel Space Observatory have revealed an unusual elliptical galaxy, NGC 4125, which has strong and extended submillimeter emission from cold dust but only very strict upper limits to its CO and HI emission. Depending on the dust emissivity, the total dust mass is 2-5x10^6 Msun. While the neutral gas-to-dust mass ratio is extremely low (< 12-30), including the ionized gas traced by [CII] emission raises this limit to < 39-100. The dust emission follows a similar r^{1/4} profile to the stellar light and the dust to stellar mass ratio is towards the high end of what is found in nearby elliptical galaxies. We suggest that NGC 4125 is currently in an unusual phase where evolved stars produced in a merger-triggered burst of star formation are pumping large amounts of gas and dust into the interstellar medium. In this scenario, the low neutral gas-to-dust mass ratio is explained by the gas being heated to temperatures >= 10^4 K faster than the dust is evaporated. If galaxies like NGC 4125, where the far-infrared emission does not trace neutral gas in the usual manner, are common at higher redshift, this could have significant implications for our understanding of high redshift galaxies and galaxy evolution. △ Less

Submitted 24 September, 2013; originally announced September 2013.

Comments: Accepted to ApJ Letters