Context. Clouds are ubiquitous in exoplanet atmospheres and they represent a challenge for the model interpretation of their spectra. When generating a large number of model spectra, complex cloud ...models often prove too costly numerically, whereas more efficient models may be overly simplified. Aims. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. Methods. We used our radiative transfer code petitRADTRANS for generating the spectra, which we coupled to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. Results. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE, and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that C/O = 0.60 −0.08 +0.07 . Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: either cloudy or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, without constraining C/O. Conclusions. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO 2 or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.
Context. The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced ...stage of disk clearing and may harbor direct signatures of disk evolution. Aims. We aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. Methods. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with VLT/SPHERE at 1.04 μm to resolve scattered light at a smaller inner working angle (0.093′′) and a higher angular resolution (0.027′′) than previously achieved. Results. We observe polarized scattered light within 0.53′′ (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r ~ 0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Conclusions. Our observations reveal the complex morphology of the disk MWC 758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.
In this paper, we review the impact of small sample statistics on detection thresholds and corresponding confidence levels (CLs) in high-contrast imaging at small angles. When looking close to the ...star, the number of resolution elements decreases rapidly toward small angles. This reduction of the number of degrees of freedom dramatically affects CLs and false alarm probabilities. Naively using the same ideal hypothesis and methods as for larger separations, which are well understood and commonly assume Gaussian noise, can yield up to one order of magnitude error in contrast estimations at fixed CL. The statistical penalty exponentially increases toward very small inner working angles. Even at 5-10 resolution elements from the star, false alarm probabilities can be significantly higher than expected. Here we present a rigorous statistical analysis that ensures robustness of the CL, but also imposes a substantial limitation on corresponding achievable detection limits (thus contrast) at small angles. This unavoidable fundamental statistical effect has a significant impact on current coronagraphic and future high-contrast imagers. Finally, the paper concludes with practical recommendations to account for small number statistics when computing the sensitivity to companions at small angles and when exploiting the results of direct imaging planet surveys.
High-contrast imaging and spectroscopy provide unique constraints for exoplanet formation models as well as for planetary atmosphere models. Instrumentation techniques in this field have greatly ...improved over the last two decades, with the development of stellar coronagraphy, in parallel with specific methods of wavefront sensing and control. Next generation space- and ground-based telescopes will enable the characterization of cold solar-system-like planets for the first time and maybe even in situ detection of bio-markers. However, the growth of primary mirror diameters, necessary for these detections, comes with an increase of their complexity (segmentation, secondary mirror features). These discontinuities in the aperture can greatly limit the performance of coronagraphic instruments. In this context, we introduced a new technique, Active Correction of Aperture Discontinuities-Optimized Stroke Minimization (ACAD-OSM), to correct for the diffractive effects of aperture discontinuities in the final image plane of a coronagraph, using deformable mirrors. In this paper, we present several tools that can be used to optimize the performance of this technique for its application to future large missions. In particular, we analyzed the influence of the deformable setup (size and separating distance) and found that there is an optimal point for this setup, optimizing the performance of the instrument in contrast and throughput while minimizing the strokes applied to the deformable mirrors. These results will help us design future coronagraphic instruments to obtain the best performance.
We present a new method for numerical propagation through Lyot-style coronagraphs using finite occulting masks. Standard methods for coronagraphic simulations involve Fast Fourier Transforms (FFT) of ...very large arrays, and computing power is an issue for the design and tolerancing of coronagraphs on segmented Extremely Large Telescopes (ELT) in order to handle both the speed and memory requirements. Our method combines a semi-analytical approach with non-FFT based Fourier transform algorithms. It enables both fast and memory-efficient computations without introducing any additional approximations. Typical speed improvements based on computation costs are of about ten to fifty for propagations from pupil to Lyot plane, with thirty to sixty times less memory needed. Our method makes it possible to perform numerical coronagraphic studies even in the case of ELTs using a contemporary commercial laptop computer, or any standard commercial workstation computer.
Future searches for bio-markers on habitable exoplanets will rely on telescope instruments that achieve extremely high contrast at small planet-to-star angular separations. Coronagraphy is a ...promising starlight suppression technique, providing excellent contrast and throughput for off-axis sources on clear apertures. However, the complexity of space- and ground-based telescope apertures goes on increasing over time, owing to the combination of primary mirror segmentation, the secondary mirror, and its support structures. These discontinuities in the telescope aperture limit the coronagraph performance. In this paper, we present ACAD-OSM, a novel active method to correct for the diffractive effects of aperture discontinuities in the final image plane of a coronagraph. Active methods use one or several deformable mirrors that are controlled with an interaction matrix to correct for the aberrations in the pupil. However, they are often limited by the amount of aberrations introduced by aperture discontinuities. This algorithm relies on the recalibration of the interaction matrix during the correction process to overcome this limitation. We first describe the ACAD-OSM technique and compare it to the previous active methods for the correction of aperture discontinuities. We then show its performance in terms of contrast and off-axis throughput for static aperture discontinuities (segmentation, struts) and for some aberrations evolving over the life of the instrument (residual phase aberrations, artifacts in the aperture, misalignments in the coronagraph design). This technique can now obtain the Earth-like planet detection threshold of contrast on any given aperture over at least a 10% spectral bandwidth, with several coronagraph designs.
We present an analysis of the orbital motion of the four substellar objects orbiting HR 8799. Our study relies on the published astrometric history of this system augmented with an epoch obtained ...with the Project 1640 coronagraph with an integral field spectrograph (IFS) installed at the Palomar Hale telescope. We first focus on the intricacies associated with astrometric estimation using the combination of an extreme adaptive optics system (PALM-3000), a coronagraph, and an IFS. We introduce two new algorithms. The first one retrieves the stellar focal plane position when the star is occulted by a coronagraphic stop. The second one yields precise astrometric and spectrophotometric estimates of faint point sources even when they are initially buried in the speckle noise. The second part of our paper is devoted to studying orbital motion in this system. In order to complement the orbital architectures discussed in the literature, we determine an ensemble of likely Keplerian orbits for HR 8799bcde, using a Bayesian analysis with maximally vague priors regarding the overall configuration of the system. Although the astrometric history is currently too scarce to formally rule out coplanarity, HR 8799d appears to be misaligned with respect to the most likely planes of HR 8799bce orbits. This misalignment is sufficient to question the strictly coplanar assumption made by various authors when identifying a Laplace resonance as a potential architecture. Finally, we establish a high likelihood that HR 8799de have dynamical masses below 13 M sub(jup), using a loose dynamical survival argument based on geometric close encounters. We illustrate how future dynamical analyses will further constrain dynamical masses in the entire system.
SUMMARY
As 3-D geological models become more numerous and widely available, the opportunity arises to combine them into large regional compilations. One of the biggest challenges facing these ...compilations is the connection and alignment of individual models, especially in less explored areas or across political borders. In this regard, gravity modelling is suitable for revealing additional subsurface information that can support a harmonization of structural models. Here, we present an integrated geological and gravity modelling approach to support the harmonization process of two geological 3-D models of the North German Basin in the cross-border region between the federal states of Saxony-Anhalt and Brandenburg. Gravity gradient calculation, filtering and Euler deconvolution are utilized to reveal new insights into the local fault system and gravity anomaly sources. The independent models are merged and harmonized during 3-D forward and inverse gravity modelling. Herein, density gradients for individual layers are incorporated in the framework of model parametrization. The resulting geological 3-D model consists of harmonized interfaces and is consistent with the observed gravity field. To demonstrate the plausibility of the derived model, we discuss the new geophysical findings on the sedimentary and crustal structures of the cross-border region in the context of the regional geological setting. The cross-border region is dominated by an NW–SE oriented fault system that coincides with the Elbe Fault System. We interpret a low-density zone within the basement of the Mid-German Crystalline Rise as a northward continuation of the Pretzsch–Prettin Crystalline Complex into the basement of the North German Basin. Additionally, we observe two types of anticlines within the basin, which we link to provinces of contrasting basement rigidity. Our gravity modelling implies that the Zechstein salt has mostly migrated into the deeper parts of the basin west of the Seyda Fault. Finally, we identify a pronounced syncline that accommodates a narrow and up to 800 m deep Cenozoic basin.
Context.
β
Pictoris is arguably one of the most studied stellar systems outside of our own. Some 30 yr of observations have revealed a highly-structured circumstellar disk, with rings, belts, and a ...giant planet:
β
Pictoris b. However very little is known about how this system came into being.
Aims.
Our objective is to estimate the C/O ratio in the atmosphere of
β
Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry.
Methods.
We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain
K
-band spectro-interferometric data on
β
Pic b. We extracted a medium resolution (
R
= 500)
K
-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio.
Results.
Our new astrometry disfavors a circular orbit for
β
Pic b (
e
= 0.15
−0.04
+0.05
). Combined with previous results and with H
IPPARCOS
/
Gaia
measurements, this astrometry points to a planet mass of
M
= 12.7 ± 2.2
M
Jup
. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of
β
Pic b. In particular, the C/O ratios derived with the two codes are identical (0.43 ± 0.05 vs. 0.43
−0.03
+0.04
). We argue that if the stellar C/O in
β
Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.
Despite the large number of magnetostratigraphic studies in the South Pyrenean Basin aiming to calibrate the basin chronostratigraphy and the biostratigraphic scales, the South Eastern Jaca Basin ...remains unexplored from this perspective, and its relation with the Ainsa Basin is not fully understood. In this work we contribute with new magnetostratigraphic data from the 950m thick Las Bellostas section, located in the northern hinge of the Balzes anticline. Well-proven primary signal (positive fold test and two pseudo-antiparallel polarities) supported by numerous primary data in the surroundings allow us building a reliable local polarity sequence of eight magnetozones (from R1 to N4). Additionally, seven new biostratigraphic samples (Nummulites and Assilina) in the lower part of the section (marine environment) allows tightening the section to the Geomagnetic Polarity Time Scale (GPTS) and proposing a refined age model for the southeastern Jaca Basin. The section starts in the Boltaña Fm., of Cuisian age (Shallow Benthic Zone, SBZ11), is followed by a sedimentary gap from C22n to C20n as witnessed by biostratigraphic data (SBZ11 underneath the hiatus and SBZ16 just atop). The deltaic Sobrarbe Fm. can be tracked until the C19n (Late Lutetian). From this point (200m) until the top of the section (950m), at least, the entire C18n chron can be recognized within the molassic Campodarbe Fm. (C18n.2n-C18n.1r and C18n.1n) (Bartonian) equivalent to the West to the prodeltaic Arguis marls Fm. The Middle Cusian (SBZ11) to Middle Lutetian (SBZ15) stratigraphic hiatus is, in part, enhanced by the structural position at the hinge of the Balzes anticline. These new chronostratigraphic constraints help refining the W-E and N-S stratigraphic relationships in the eastern Jaca Basin and in the Ainsa Basin. This section also allows us to accurately refine the kinematics of the rotational activity in the eastern External Sierras. The important difference in magnetic declination along the section and neighboring paleomagnetic data from the Balzes anticline (from ≈70º clockwise at the base of the stratigraphic section to non-significant at the top) together with the new age model for the Eastern Jaca Basin help characterizing the rotational activity of the Balzes thrust sheet. The rotation took place between chrons C20r (Middle Lutetian; 45Ma) and C17 (Lower Priabonian 37-38Ma) in agreement to nearby structures (Boltaña, Picodel Aguila anticlines) but clearly diachronic to western ones (Santo Domingo anticline). Besides, the rotational activity seems to follow a linear and continuous pattern (velocity 9-11º/Ma, R: 0.83-0.96) in contrast to closer structures that show two distinct rotational velocities (i.e. Boltaña). These new data still let open the debate on the rotational kinematics along the South Pyrenean basal thrust.