Source author record

Willy Benz

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

23 published item(s)

preprint2022arXiv

A pair of Sub-Neptunes transiting the bright K-dwarf TOI-1064 characterised with CHEOPS

We report the discovery and characterisation of a pair of sub-Neptunes transiting the bright K-dwarf TOI-1064 (TIC 79748331), initially detected in TESS photometry. To characterise the system, we performed and retrieved CHEOPS, TESS, and ground-based photometry, HARPS high-resolution spectroscopy, and Gemini speckle imaging. We characterise the host star and determine $T_{\rm eff, \star}=4734\pm67$ K, $R_{\star}=0.726\pm0.007$ $R_{\odot}$, and $M_{\star}=0.748\pm0.032$ $M_{\odot}$. We present a novel detrending method based on PSF shape-change modelling and demonstrate its suitability to correct flux variations in CHEOPS data. We confirm the planetary nature of both bodies and find that TOI-1064 b has an orbital period of $P_{\rm b}=6.44387\pm0.00003$ d, a radius of $R_{\rm b}=2.59\pm0.04$ $R_{\oplus}$, and a mass of $M_{\rm b}=13.5_{-1.8}^{+1.7}$ $M_{\oplus}$, whilst TOI-1064 c has an orbital period of $P_{\rm c}=12.22657^{+0.00005}_{-0.00004}$ d, a radius of $R_{\rm c}=2.65\pm0.04$ $R_{\oplus}$, and a 3$σ$ upper mass limit of 8.5 ${\rm M_{\oplus}}$. From the high-precision photometry we obtain radius uncertainties of $\sim$1.6%, allowing us to conduct internal structure and atmospheric escape modelling. TOI-1064 b is one of the densest, well-characterised sub-Neptunes, with a tenuous atmosphere that can be explained by the loss of a primordial envelope following migration through the protoplanetary disc. It is likely that TOI-1064 c has an extended atmosphere due to the tentative low density, however further RVs are needed to confirm this scenario and the similar radii, different masses nature of this system. The high-precision data and modelling of TOI-1064 b are important for planets in this region of mass-radius space, and it allows us to identify a trend in bulk density-stellar metallicity for massive sub-Neptunes that may hint at the formation of this population of planets.

preprint2022arXiv

Uncovering the true periods of the young sub-Neptunes orbiting TOI-2076

Context: TOI-2076 is a transiting three-planet system of sub-Neptunes orbiting a bright (G = 8.9 mag), young ($340\pm80$ Myr) K-type star. Although a validated planetary system, the orbits of the two outer planets were unconstrained as only two non-consecutive transits were seen in TESS photometry. This left 11 and 7 possible period aliases for each. Aims: To reveal the true orbits of these two long-period planets, precise photometry targeted on the highest-probability period aliases is required. Long-term monitoring of transits in multi-planet systems can also help constrain planetary masses through TTV measurements. Methods: We used the MonoTools package to determine which aliases to follow, and then performed space-based and ground-based photometric follow-up of TOI-2076 c and d with CHEOPS, SAINT-EX, and LCO telescopes. Results: CHEOPS observations revealed a clear detection for TOI-2076 c at $P=21.01538^{+0.00084}_{-0.00074}$ d, and allowed us to rule out three of the most likely period aliases for TOI-2076 d. Ground-based photometry further enabled us to rule out remaining aliases and confirm the $P=35.12537\pm0.00067$ d alias. These observations also improved the radius precision of all three sub-Neptunes to $2.518\pm0.036$, $3.497\pm0.043$, and $3.232\pm0.063$ $R_\oplus$. Our observations also revealed a clear anti-correlated TTV signal between planets b and c likely caused by their proximity to the 2:1 resonance, while planets c and d appear close to a 5:3 period commensurability, although model degeneracy meant we were unable to retrieve robust TTV masses. Their inflated radii, likely due to extended H-He atmospheres, combined with low insolation makes all three planets excellent candidates for future comparative transmission spectroscopy with JWST.

preprint2021arXiv

Fundamental physics with ESPRESSO: Precise limit on variations in the fine-structure constant towards the bright quasar HE 0515$-$4414

The strong intervening absorption system at redshift 1.15 towards the very bright quasar HE 0515$-$4414 is the most studied absorber for measuring possible cosmological variations in the fine-structure constant, $α$. We observed HE 0515$-$4414 for 16.1$\,$h with the Very Large Telescope and present here the first constraint on relative variations in $α$ with parts-per-million (ppm) precision from the new ESPRESSO spectrograph: $Δα/α= 1.3 \pm 1.3_{\rm stat} \pm 0.4_{\rm sys}\,{\rm ppm}$. The statistical uncertainty (1$σ$) is similar to the ensemble precision of previous large samples of absorbers, and derives from the high S/N achieved ($\approx$105 per 0.4$\,$km$\,$s$^{-1}$ pixel). ESPRESSO's design, and calibration of our observations with its laser frequency comb, effectively removed wavelength calibration errors from our measurement. The high resolving power of our ESPRESSO spectrum ($R=145000$) enabled the identification of very narrow components within the absorption profile, allowing a more robust analysis of $Δα/α$. The evidence for the narrow components is corroborated by their correspondence with previously detected molecular hydrogen and neutral carbon. The main remaining systematic errors arise from ambiguities in the absorption profile modelling, effects from redispersing the individual quasar exposures, and convergence of the parameter estimation algorithm. All analyses of the spectrum, including systematic error estimates, were initially blinded to avoid human biases. We make our reduced ESPRESSO spectrum of HE 0515$-$4414 publicly available for further analysis. Combining our ESPRESSO result with 28 measurements, from other spectrographs, in which wavelength calibration errors have been mitigated, yields a weighted mean $Δα/α= -0.5 \pm 0.5_{\rm stat} \pm 0.4_{\rm sys}\,$ppm at redshifts 0.6-2.4.

preprint2020arXiv

Expected performances of the Characterising Exoplanet Satellite (CHEOPS) II. The CHEOPS simulator

The CHaracterising ExOPlanet Satellite (CHEOPS) is a mission dedicated to the search for exoplanetary transits through high precision photometry of bright stars already known to host planets. The telescope will provide the unique capability of determining accurate radii for planets whose masses have already been measured from ground-based spectroscopic surveys. This will allow a first-order characterisation of the planets' internal structure through the determination of the bulk density, providing direct insight into their composition. The CHEOPS simulator has been developed to perform detailed simulations of the data which is to be received from the CHEOPS satellite. It generates accurately simulated images that can be used to explore design options and to test the on-ground data processing, in particular, the pipeline producing the photometric time series. It is, thus, a critical tool for estimating the photometric performance expected in flight and to guide photometric analysis. It can be used to prepare observations, consolidate the noise budget, and asses the performance of CHEOPS in realistic astrophysical fields that are difficult to reproduce in the laboratory. Images generated by CHEOPSim take account of many detailed effects, including variations of the incident signal flux and backgrounds, and detailed modelling of the satellite orbit, pointing jitter and telescope optics, as well as the CCD response, noise and readout. The simulator results presented in this paper have been used in the context of validating the data reduction processing chain, in which image time series generated by CHEOPSim were used to generate light curves for simulated planetary transits across real and simulated targets. Independent analysts were successfully able to detect the planets and measure their radii to an accuracy within the science requirements of the mission.

preprint2020arXiv

Nightside condensation of iron in an ultra-hot giant exoplanet

Ultra-hot giant exoplanets receive thousands of times Earth's insolation. Their high-temperature atmospheres (>2,000 K) are ideal laboratories for studying extreme planetary climates and chemistry. Daysides are predicted to be cloud-free, dominated by atomic species and substantially hotter than nightsides. Atoms are expected to recombine into molecules over the nightside, resulting in different day-night chemistry. While metallic elements and a large temperature contrast have been observed, no chemical gradient has been measured across the surface of such an exoplanet. Different atmospheric chemistry between the day-to-night ("evening") and night-to-day ("morning") terminators could, however, be revealed as an asymmetric absorption signature during transit. Here, we report the detection of an asymmetric atmospheric signature in the ultra-hot exoplanet WASP-76b. We spectrally and temporally resolve this signature thanks to the combination of high-dispersion spectroscopy with a large photon-collecting area. The absorption signal, attributed to neutral iron, is blueshifted by -11+/-0.7 km s-1 on the trailing limb, which can be explained by a combination of planetary rotation and wind blowing from the hot dayside. In contrast, no signal arises from the nightside close to the morning terminator, showing that atomic iron is not absorbing starlight there. Iron must thus condense during its journey across the nightside.

preprint2020arXiv

Pebbles versus planetesimals: the outcomes of population synthesis models

In the core accretion scenario, a massive core forms first and then accretes an envelope. When discussing how this core forms some divergences appear. First scenarios of planet formation predict the accretion of km-sized bodies, called planetesimals, while more recent works suggest growth by accretion of pebbles, which are cm-sized objects. These two accretion models are often discussed separately and we aim here at comparing the outcomes of the two models with identical initial conditions. We use two distinct codes: one computing planetesimal accretion, the other pebble accretion. Using a population synthesis approach, we compare planet simulations and study the impact of the two solid accretion models, focussing on the formation of single planets. We find that the planetesimal model predicts the formation of more giant planets, while the pebble accretion model forms more super-Earth mass planets. This is due to the pebble isolation mass concept, which prevents planets formed by pebble accretion to accrete gas efficiently before reaching Miso. This translates into a population of planets that are not heavy enough to accrete a consequent envelope but that are in a mass range where type I migration is very efficient. We also find higher gas mass fractions for a given core mass for the pebble model compared to the planetesimal one caused by luminosity differences. This also implies planets with lower densities which could be confirmed observationally. Focusing on giant planets, we conclude that the sensitivity of their formation differs: for the pebble accretion model, the time at which the embryos are formed, as well as the period over which solids are accreted strongly impact the results, while for the planetesimal model it depends on the planetesimal size and on the splitting in the amount of solids available to form planetesimals.

preprint2019arXiv

Expected performances of the Characterising Exoplanet Satellite (CHEOPS). I. Photometric performances from ground-based calibration

The Characterising Exoplanet Satellite (CHEOPS) is a space mission designed to perform photometric observations of bright stars to obtain precise radii measurements of transiting planets. The high-precision photometry of CHEOPS relies on careful on-ground calibration of its payload. For that purpose, intensive pre-launch campaigns of measurements were carried out to calibrate the instrument and characterise its photometric performances. We report on main results of these campaigns, provide a complete analysis of data sets and estimate in-flight photometric performance by mean of end-to-end simulation. The on-ground photometric stability of the instrument is found to be of the order of 15 parts per million over 5 hours. Our end-to-end simulation shows that measurements of planet-to-star radii ratio with CHEOPS can be determined with a precision of 2% for a Neptune-size planet transiting a K-dwarf star and 5% for an Earth-size planet orbiting a Sun-like star. It corresponds to signal-to-noise ratios on the transit depths of 25 and 10 respectively, allowing the characterisation and detection of these planets. The pre-launch CHEOPS performances are shown to be compliant with the mission requirements.

preprint2016arXiv

Planet formation with envelope enrichment: new insights on planetary diversity

We compute, for the first time, self-consistent models of planet growth including the effect of envelope enrichment. The change of envelope metallicity is assumed to be the result of planetesimal disruption or icy pebble sublimation. We solve internal structure equations taking into account global energy conservation for the envelope to compute in-situ planetary growth. We consider different opacities and equations of state suited for a wide range of metallicities. We find that envelope enrichment speeds up the formation of gas giants. It also explains naturally the formation of low and intermediate mass objects with large fractions of H-He (~ 20 - 30 % in mass). High opacity models explain well the metallicity of the giant planets of the solar system, whereas low opacity models are suited for forming small mass objects with thick H-He envelopes and gas giants with sub-solar envelope metallicities. We find good agreement between our models and the estimated water abundance for WASP-43b. For HD 189733b, HD 209458b and WASP-12b we predict fractions of water larger than what is estimated from observations, by at least a factor ~ 2. Envelope enrichment by icy planetesimals is the natural scenario to explain the formation of a large variety of objects, ranging from mini-Neptunes, to gas giants. We predict that the total and envelope metallicity decrease with planetary mass.

preprint2015arXiv

Can we constrain interior structure of rocky exoplanets from mass and radius measurements?

We present an inversion method based on Bayesian analysis to constrain the interior structure of terrestrial exoplanets, in the form of chemical composition of the mantle and core size. Specifically, we identify what parts of the interior structure of terrestrial exoplanets can be determined from observations of mass, radius, and stellar elemental abundances. We perform a full probabilistic inverse analysis to formally account for observational and model uncertainties and obtain confidence regions of interior structure models. This enables us to characterize how model variability depends on data and associated uncertainties. We test our method on terrestrial solar system planets and find that our model predictions are consistent with independent estimates. Furthermore, we apply our method to synthetic exoplanets up to 10 Earth masses and up to 1.7 Earth radii as well as to exoplanet Kepler-36b. Importantly, the inversion strategy proposed here provides a framework for understanding the level of precision required to characterize the interior of exoplanets. Our main conclusions are: (1) observations of mass and radius are sufficient to constrain core size; (2) stellar elemental abundances (Fe, Si, Mg) are key constraints to reduce degeneracy in interior structure models and to constrain mantle composition; (3) the inherent degeneracy in determining interior structure from mass and radius observations does not only depend on measurement accuracies but also on the actual size and density of the exoplanet. We argue that precise observations of stellar elemental abundances are central in order to place constraints on planetary bulk composition and to reduce model degeneracy. [...]

preprint2014arXiv

From planetesimals to planets: volatile molecules

Solar and extrasolar planets are the subject of numerous studies aiming to determine their chemical composition and internal structure. In the case of extrasolar planets, the composition is important as it partly governs their potential habitability. Moreover, observational determination of chemical composition of planetary atmospheres are becoming available, especially for transiting planets. The present works aims at determining the chemical composition of planets formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data for models of planet formation and evolution, and future interpretation of chemical composition of solar and extrasolar planets. We have developed a model that computes the composition of ices in planets in different stellar systems with the use of models of ice and planetary formation. We provide the chemical composition, ice/rock mass ratio and C:O molar ratio for planets in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, we produce a wide variety of planetary chemical compositions as a function of the mass of the disk and distance to the star. The volatile species incorporated in planets are mainly composed of H2O, CO, CO2, CH3OH, and NH3. Icy or ocean planets have systematically higher values of molecular abundances compared to giant and rocky planets. Gas giant planets are depleted in highly volatile molecules such as CH4, CO, and N2 compared to icy or ocean planets. The ice/rock mass ratio in icy or ocean and gas giant planets is, respectively, equal at maximum to 1.01+-0.33 and 0.8+-0.5, and is different from the usual assumptions made in planet formation models, which suggested this ratio to be 2-3. The C:O molar ratio in the atmosphere of gas giant planets is depleted by at least 30% compared to solar value.

preprint2014arXiv

From stellar nebula to planetesimals

Solar and extrasolar comets and extrasolar planets are the subject of numerous studies in order to determine their chemical composition and internal structure. In the case of planetesimals, their compositions are important as they govern in part the composition of future planets. The present works aims at determining the chemical composition of icy planetesimals, believed to be similar to present day comets, formed in stellar systems of solar chemical composition. The main objective of this work is to provide valuable theoretical data on chemical composition for models of planetesimals and comets, and models of planet formation and evolution. We have developed a model that calculates the composition of ices formed during the cooling of the stellar nebula. Coupled with a model of refractory element formation, it allows us to determine the chemical composition and mass ratio of ices to rocks in icy planetesimals throughout in the protoplanetary disc. We provide relationships for ice line positions (for different volatile species) in the disc, and chemical compositions and mass ratios of ice relative to rock for icy planetesimals in stellar systems of solar chemical composition. From an initial homogeneous composition of the nebula, a wide variety of chemical compositions of planetesimals were produced as a function of the mass of the disc and distance to the star. Ices incorporated in planetesimals are mainly composed of H2O, CO, CO2, CH3OH, and NH3. The ice/rock mass ratio is equal to 1+-0.5 in icy planetesimals following assumptions. This last value is in good agreement with observations of solar system comets, but remains lower than usual assumptions made in planet formation models, taking this ratio to be of 2-3.

preprint2011arXiv

Extrasolar planet population synthesis III. Formation of planets around stars of different masses

We extend the models presented in Mordasini et al. (2009) to the formation of planets orbiting stars of different masses. We discuss the properties of the resulting synthetic planet population in terms of mass, orbit, and metallicity distributions. The population synthesis calculations we use are based on the planet formation model developed in Alibert et al. (2005a), which self-consistently takes into account planetary growth and migration in an evolving proto-planetary disk. Using this model, we generate synthetic populations of planets by following their growth in a large number of proto-planetary disks, whose properties (mass and lifetime) are selected in a Monte Carlo fashion using probability distributions derived from observations. We show that the scaling of the proto-planetary disk mass with the mass of the central star has a direct and large influence on the properties of the resulting planet population. In particular, the observed paucity of large mass planets orbiting 0.5 Msun stars can directly be explained as resulting from a only slightly steeper than linear scaling. The observed lack of short period planet orbiting 2.0 Msun stars can also be attributed to this scaling but only if associated with a decrease of the mean disk lifetime for stars more massive than 1.5 Msun. Finally, we show that the distribution of minimum mass and semi-major axis of our synthetic planets are statistically comparable with observations.

preprint2011arXiv

Giant planet formation: episodic impacts vs. gradual core growth

We describe the growth of gas giant planets in the core accretion scenario. The core growth is not modeled as a gradual accretion of planetesimals but as episodic impacts of large mass ratios, i.e. we study impacts of 0.02 - 1 Earth masses onto cores of 1-15 Earth masses. Such impacts could deliver the majority of solid matter in the giant impact regime. We focus on the thermal response of the envelope to the energy delivery. Previous studies have shown that sudden shut off of core accretion can dramatically speed up gas accretion. We therefore expect that giant impacts followed by periods of very low core accretion will result in a net increase in gas accretion rate. This study aims at modelling such a sequence of events and to understand the reaction of the envelope to giant impacts in more detail. To model this scenario, we spread the impact energy deposition over a time that is long compared to the sound crossing time, but very short compared to the Kelvin-Helmholtz time. The simulations are done in spherical symmetry and assume quasi-hydrostatic equilibrium. Results confirm what could be inferred from previous studies: gas can be accreted faster onto the core for the same net core growth speed while at the same time rapid gas accretion can occur for smaller cores -- significantly smaller than the usual critical core mass. Furthermore our simulations show, that significant mass fractions of the envelope can be ejected by such an impact.

preprint2010arXiv

The HARPS search for southern extrasolar planets. XXIII. 8 planetary companions to low-activity solar-type stars

In this paper, we present our HARPS radial-velocity data for eight low-activity solar-type stars belonging to the HARPS volume-limited sample: HD6718, HD8535, HD28254, HD290327, HD43197, HD44219, HD148156, and HD156411. Keplerian fits to these data reveal the presence of low-mass companions around these targets. With minimum masses ranging from 0.58 to 2.54 MJup, these companions are in the planetary mass domain. The orbital periods of these planets range from slightly less than one to almost seven years. The eight orbits presented in this paper exhibit a wide variety of eccentricities: from 0.08 to above 0.8.

preprint2010arXiv

Theory of planet formation

We review the current theoretical understanding how growth from micro-meter sized dust to massive giant planets occurs in disks around young stars. After introducing a number of observational constraints from the solar system, from observed protoplanetary disks, and from the extrasolar planets, we simplify the problem by dividing it into a number of discrete stages which are assumed to occur in a sequential way. In the first stage - the growth from dust to kilometer sized planetesimals - the aerodynamics of the bodies are of central importance. We discuss both a purely coagulative growth mode, as well as a gravoturbulent mode involving a gravitational instability of the dust. In the next stage, planetesimals grow to protoplanets of roughly 1000 km in size. Gravity is now the dominant force. The mass accretion can be strongly non-linear, leading to the detachment of a few big bodies from the remaining planetesimals. In the outer planetary system (outside a few AU), some of these bodies can become so massive that they eventually accrete a large gaseous envelope. This is the stage of giant planet formation, as understood within the core accretion-gas capture paradigm. We also discuss the direct gravitational collapse model where giant planets are thought to form directly via a gravitational fragmentation of the gas disk. In the inner system, protoplanets collide in the last stage - probably after the dispersal of the gaseous disk - in giant impacts until the separations between the remaining terrestrial planets become large enough to allow long term stability. We finish the review with some selected questions.

preprint2009arXiv

Extrasolar planet population synthesis I: Method, formation tracks and mass-distance distribution

With the high number of extrasolar planets discovered by now, it becomes possible to constrain theoretical formation models in a statistical sense. This paper is the first in a series in which we carry out a large number of planet population synthesis calculations. We begin the series with a paper mainly dedicated to the presentation of our approach, but also the discussion of a representative synthetic planetary population of solar like stars. Based as tightly as possible on observational data, we have derived probability distributions for the most important initial conditions for the planetary formation process. We then draw sets of initial conditions from these distributions and obtain the corresponding synthetic planets with our formation model. Although the main purpose of this paper is the description of our methods, we present some key results: We find that the variation of the initial conditions in the limits occurring in nature leads to the formation of planets of large diversity. This formation process is best visualized in planetary formation tracks, where different phases of concurrent growth and migration can be identified. These phases lead to the emergence of sub-populations of planets distinguishable in a mass-semimajor axis diagram. The most important ones are the "failed cores", a vast group of core-dominated low mass planets, the "horizontal branch", a sub-population of Neptune mass planets extending out to 6 AU, and the "main clump", a concentration of giant gaseous giants planets at around 0.3-2 AU.

preprint2009arXiv

Extrasolar planet population synthesis II: Statistical comparison with observation

This is the second paper in a series of papers showing the results of extrasolar planet population synthesis calculations. In the companion paper (Paper I), we have presented in detail our methods. By applying an observational detection bias for radial velocity surveys, we identify the potentially detectable synthetic planets. The properties of these planets are compared in quantitative statistical tests with the properties of a carefully selected sub-population of actual exoplanets. We use a two dimensional Kolmogorov-Smirnov test to compare the mass-distance distributions of synthetic and observed planets, as well as 1D KS tests to compare the mass, the semimajor axis and the [Fe/H] distributions. We find that some models can account to a reasonable degree of significance for the observed properties. We concurrently account for many other observed features, e.g. the "metallicity effect". This gives us confidence that our model captures several essential features of giant planet formation. Our simulations allow us also to extract the properties of the underlying exoplanet population that are not yet detectable. For example, we have derived the planetary initial mass function (PIMF) and have been led to conclude that the planets detected so far represent only the tip of the iceberg. The PIMF can also be used to predict how the detectable extrasolar planet population will change as the precision of radial velocity surveys improves to an extreme precision of 0.1 m/s.

preprint2009arXiv

Fragment properties at the catastrophic disruption threshold: The effect of the parent body's internal structure

Numerical simulations of asteroid break-ups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by Q*D, which results in the escape of half of the target's mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize Q*D for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of Q*D, and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of Q*D on the impact speed of the projectile. In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. We then propose some power-law relationships between Q*D and both target's size and impact speed that can be used in collisional evolution models.

preprint2009arXiv

The Rossiter-McLaughlin effect of CoRoT-3b & HD189733b

We present radial-velocity sequences acquired during three transits of the exoplanet HD 189733b and one transit of the CoRoT-3b. We applied a combined Markov-Chain Monte Carlo analysis of spectroscopic and photometric data on these stars, to determine a full set of system parameters including the project spin-orbit misalignement angle of HD 189733b to an unprecedented precision via the Rossiter-McLaughlin effect: beta = 0.85 degrees (+0.32 -0.28) . This small but non-zero inclination of the planetary orbit is important to understand the origin of the system. On CoRoT-3b, results seem to point towards a non-zero inclination as well with beta = 37.6 degrees (+10.0 -22.3), but this remains marginal. Systematic effects due to non-gaussian cross-correlation functions appear to be the main cause of significant residuals that prevent an accurate determination of the projected stellar rotation velocity V sin(I) for both stars.

preprint2004arXiv

Models of Giant Planet formation with migration and disc evolution

We present a new model of giant planet formation that extends the core-accretion model of Pollack etal (1996) to include migration, disc evolution and gap formation. We show that taking into account these effects can lead to a much more rapid formation of giant planets, making it compatible with the typical disc lifetimes inferred from observations of young circumstellar discs. This speed up is due to the fact that migration prevents the severe depletion of the feeding zone as observed in in situ calculations. Hence, the growing planet is never isolated and it can reach cross-over mass on a much shorter timescale. To illustrate the range of planets that can form in our model, we describe a set of simulations in which we have varied some of the initial parameters and compare the final masses and semi-major axes with those inferred from observed extra-solar planets.

preprint1999arXiv

Towards Understanding Jovian Planet Migration

We present 2D hydrodynamic simulations of circumstellar disks around protostars using a `Piecewise Parabolic Method' (PPM) code. We include a point mass embedded within the disk and follow the migration of that point mass through the disk. Companions with masses $M_c\ga 0.5M_J$ can open a gap in the disk sufficient to halt rapid migration through the disk. Lower mass companions open gaps, but migration continues because sufficient disk mass remains close to the disk to exert large tidal torques. We find that the torques which dominate the migration of low mass planets originate within a radial region within 1-2 Hill radii of the planet's orbit radius, a distance smaller than the thickness of the disk. We conclude that a very high resolution 3D treatment will be required to adequately describe the planet's migration.

preprint1995arXiv

The Dynamics and Outcomes of Rapid Infall onto Neutron Stars

We present an extensive study of accretion onto neutron stars in which the velocity of the neutron star and structure of the surrounding medium is such that the Bondi-Hoyle accretion exceeds .001 Msun/y. For most cases, hypercritical accretion due to rapid neutrino cooling allows the neutron star to accrete above the Bondi-Hoyle rate as previously pointed out by Chevalier. However, for a subset of simulations which corresponds to evolutionarily common events, convection driven by neutrino heating can lead to explosions by a mechanism similar to that found in core-collapse supernovae. Armed with the results from our calculations, we are in a position to predict the fate of a range of rapid-infall neutron star accretors present in certain low-mass X-ray binaries, common envelope systems, supernova fallbacks and Thorne-Zytkow objects (TZOs). A majority of the common envelope systems that we considered led to explosions expelling the envelope, halting the neutron star's inward spiral, and allowing the formation of close binary systems. Smothered neutron stars produced in collisions also lead to explosions, preventing them from forming millisecond pulsars. For supernovae in which the fallback of material towards the neutron star is large, we find that a black hole is formed within a few seconds. Finally, we argue that the current set of TZO formation scenarios is inadequate and leads instead to hypercritical accretion and black hole formation. Moreover, it appears that many of the current TZ models have structures ill-suited for modeling by mixing length convection. This has prompted us to develop a simple test to determine the viability of this approximation for a variety of convective systems.