| 1. |
- Chauvin, G., et al.
(författare)
-
Discovery of a warm, dusty giant planet around HIP 65426
- 2017
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 605
-
Tidskriftsartikel (refereegranskat)abstract
- Aims. The SHINE program is a high-contrast near-infrared survey of 600 young, nearby stars aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular-resolution imaging capabilities. It is also intended to place statistical constraints on the rate, mass and orbital distributions of the giant planet population at large orbits as a function of the stellar host mass and age to test planet-formation theories.Methods. We used the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-contrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP 65426. It is a member of the similar to 17 Myr old Lower Centaurus-Crux association. Results. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP 65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2 mu m indicate a warm, dusty atmosphere characteristic of young low-surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 M-Jup, T-eff = 1300-1600K and R = 1.5 +/- 0.1 R-Jup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g) = 4.0-5.0 with smaller radii (1.0-1.3 R-Jup).Conclusions. Given its physical and spectral properties, HIP 65426 b occupies a rather unique placement in terms of age, mass, and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.
|
|
| 2. |
- Gratton, R., et al.
(författare)
-
Investigating three Sirius-like systems with SPHERE
- 2021
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 646
-
Tidskriftsartikel (refereegranskat)abstract
- Context. Sirius-like systems are relatively wide binaries with a separation from a few to hundreds of au; they are composed of a white dwarf (WD) and a companion of a spectral type earlier than M0. Here we consider main sequence (MS) companions, where the WD progenitor evolves in isolation, but its wind during the former asymptotic giant branch (AGB) phase pollutes the companion surface and transfers some angular momentum. They are rich laboratories to constrain stellar models and binary evolution.Aims. Within the SpHere INfrared survey for Exoplanet survey that uses the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument at the Very Large Telescope, our goal is to acquire high contrast multi-epoch observations of three Sirius-like systems, HD 2133, HD 114174, and CD-56 7708 and to combine this data with archive high resolution spectra of the primaries, TESS archive, and literature data.Methods. These WDs are easy targets for SPHERE and were used as spectrophotometric standards. We performed very accurate abundance analyses for the MS stars using methods considered for solar analogs. Whenever possible, WD parameters and orbits were obtained using Monte Carlo Markov chain methods.Results. We found brighter J and K magnitudes for HD 114174B than obtained previously and extended the photometry down to 0.95 μm. Our new data indicate a higher temperature and then shorter cooling age (5.57 ± 0.02 Gyr) and larger mass (0.75 ± 0.03 M⊙) for this WD than previously assumed. Together with the oldest age for the MS star connected to the use of the Gaia DR2 distance, this solved the discrepancy previously found with the age of the MS star. The two other WDs are less massive, indicating progenitors of ∼1.3 M⊙ and 1.5 − 1.8 M⊙ for HD 2133B and CD-56 7708B, respectively. In spite of the rather long periods, we were able to derive useful constraints on the orbit for HD 114174 and CD-56 7708. They are both seen close to edge-on, which is in agreement with the inclination of the MS stars that are obtained coupling the rotational periods, stellar radii, and the projected rotational velocity from spectroscopy. The composition of the MS stars agrees fairly well with expectations from pollution by the AGB progenitors of the WDs: HD 2133A has a small enrichment of n-capture elements, which is as expected for pollution by an AGB star with an initial mass < 1.5 M⊙; CD-56 7708A is a previously unrecognized mild Ba-star, which is also expected due to pollution by an AGB star with an initial mass in the range of 1.5 − 3.0 M⊙; and HD 114174 has a very moderate excess of n-capture elements, which is in agreement with the expectation for a massive AGB star to have an initial mass > 3.0 M⊙.Conclusions. On the other hand, none of these stars show the excesses of C that are expected to go along with those of n-capture elements. This might be related to the fact that these stars are at the edges of the mass range where we expect nucleosynthesis related to thermal pulses. More work, both theoretical and observational, is required to better understand this issue.
|
|
| 3. |
- Bonavita, M., et al.
(författare)
-
New binaries from the SHINE survey
- 2022
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 663
-
Tidskriftsartikel (refereegranskat)abstract
- We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for sub-stellar companions to young stars using high contrast imaging. Although stars with known stellar companions within the SPHERE field of view (< 5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets observed so far. Twenty-seven per cent of the systems have three or more components. Given the heterogeneity of the sample in terms of observing conditions and strategy, tailored routines were used for data reduction and analysis, some of which were specifically designed for these datasets. We then combined SPHERE data with literature and archival data, TESS light curves, and Gaia parallaxes and proper motions for an accurate characterisation of the systems. Combining all data, we were able to constrain the orbits of 25 systems. We carefully assessed the completeness of our sample for separations between 50–500 mas (corresponding to periods of a few years to a few decades), taking into account the initial selection biases and recovering part of the systems excluded from the original list due to their multiplicity. This allowed us to compare the binary frequency for our sample with previous studies and highlight interesting trends in the mass ratio and period distribution. We also found that, when such an estimate was possible, the values of the masses derived from dynamical arguments were in good agreement with the model predictions. Stellar and orbital spins appear fairly well aligned for the 12 stars that have enough data, which favours a disk fragmentation origin. Our results highlight the importance of combining different techniques when tackling complex problems such as the formation of binaries and show how large samples can be useful for more than one purpose.
|
|
| 4. |
- Bonnefoy, M., et al.
(författare)
-
The GJ 504 system revisited Combining interferometric, radial velocity, and high contrast imaging data
- 2018
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 618
-
Tidskriftsartikel (refereegranskat)abstract
- Context. The G-type star GJ504A is known to host a 3-35 M-Jup companion whose temperature, mass, and projected separation all contribute to making it a test case for planet formation theories and atmospheric models of giant planets and light brown dwarfs. Aims. We aim at revisiting the system age, architecture, and companion physical and chemical properties using new complementary interferometric, radial-velocity, and high-contrast imaging data. Methods. We used the CHARA interferometer to measure GJ504A's angular diameter and obtained an estimation of its radius in combination with the HIPPARCOS parallax. The radius was compared to evolutionary tracks to infer a new independent age range for the system. We collected dual imaging data with IRDIS on VLT/SPHERE to sample the near-infrared (1.02-2.25 mu m) spectral energy distribution (SED) of the companion. The SED was compared to five independent grids of atmospheric models (petitCODE, Exo-REM, BT-SETTL, Morley et al., and ATMO) to infer the atmospheric parameters of GJ 504b and evaluate model-to-model systematic errors. In addition, we used a specific model grid exploring the effect of different C/O ratios. Contrast limits from 2011 to 2017 were combined with radial velocity data of the host star through the MESS2 tool to define upper limits on the mass of additional companions in the system from 0.01 to 100 au. We used an MCMC fitting tool to constrain the companion's orbital parameters based on the measured astrometry, and dedicated formation models to investigate its origin. Results. We report a radius of 1.35 +/- 0.04 R-circle dot for GJ504A. The radius yields isochronal ages of 21 +/- 2 Myr or 4.0 +/- 1.8 Gyr for the system and line-of-sight stellar rotation axis inclination of 162.4(-4.3)(+3.8) degrees or 18.6(-3.8)(+4.3) degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual-band images. The complete 1-4 mu m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages (<= 1.5Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All atmospheric models yield T-eff = 550 +/- 50 K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics; it is not degenerate with the C/O ratio. We derive log L/L-circle dot = 6.15 +/- 0.15 dex for the companion from the empirical analysis and spectral synthesis. The luminosity and T-eff yield masses of M = 1.3(-0.3)(+0.6) M-Jup and M = 23(-9)(+10) M-Jup for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity is lower than 0.55. The posterior on GJ 504b's orbital inclination suggests a misalignment with the rotation axis of GJ 504A. We exclude additional objects (90% prob.) more massive than 2.5 and 30 M-Jup with semi-major axes in the range 0.01-80 au for the young and old isochronal ages, respectively. Conclusions. The mass and semi-major axis of GJ 504b are marginally compatible with a formation by disk-instability if the system is 4 Gyr old. The companion is in the envelope of the population of planets synthesized with our core-accretion model. Additional deep imaging and spectroscopic data with SPHERE and JWST should help to confirm the possible spin-orbit misalignment and refine the estimates on the companion temperature, luminosity, and atmospheric composition.
|
|
| 5. |
- Lagrange, A. -M., et al.
(författare)
-
Post-conjunction detection of beta Pictoris b with VLT/SPHERE
- 2019
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 621
-
Tidskriftsartikel (refereegranskat)abstract
- Context. With an orbital distance comparable to that of Saturn in the solar system, beta Pictoris b is the closest (semi-major axis similar or equal to 9 au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to beta Pictoris have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters.Aims. We aimed at further constraining beta Pictoris b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit.Methods. We used SPHERE at the VLT to precisely monitor the orbital motion of beta beta Pictoris b since first light of the instrument in 2014.Results. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6 au) and prevented further detection. We redetected beta Pictoris b on the northeast side of the disk at a separation of 139 mas and a PA of 30 degrees in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of a = 9.0 +/- 0.5 au (1 sigma), it definitely excludes previously reported possible long orbital periods, and excludes beta Pictoris b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.
|
|
| 6. |
- Lazzoni, C., et al.
(författare)
-
Dynamical models to explain observations with SPHERE in planetary systems with double debris belts
- 2018
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 611
-
Tidskriftsartikel (refereegranskat)abstract
- Context. A large number of systems harboring a debris disk show evidence for a double belt architecture. One hypothesis for explaining the gap between the debris belts in these disks is the presence of one or more planets dynamically carving it. For this reason these disks represent prime targets for searching planets using direct imaging instruments, like the Spectro-Polarimetric High-constrast Exoplanet Research (SPHERE) at the Very Large Telescope. Aims. The goal of this work is to investigate this scenario in systems harboring debris disks divided into two components, placed, respectively, in the inner and outer parts of the system. All the targets in the sample were observed with the SPHERE instrument, which performs high-contrast direct imaging, during the SHINE guaranteed time observations. Positions of the inner and outer belts were estimated by spectral energy distribution fitting of the infrared excesses or, when available, from resolved images of the disk. Very few planets have been observed so far in debris disks gaps and we intended to test if such non-detections depend on the observational limits of the present instruments. This aim is achieved by deriving theoretical predictions of masses, eccentricities, and semi-major axes of planets able to open the observed gaps and comparing such parameters with detection limits obtained with SPHERE. Methods. The relation between the gap and the planet is due to the chaotic zone neighboring the orbit of the planet. The radial extent of this zone depends on the mass ratio between the planet and the star, on the semi-major axis, and on the eccentricity of the planet, and it can be estimated analytically. We first tested the different analytical predictions using a numerical tool for the detection of chaotic behavior and then selected the best formula for estimating a planet's physical and dynamical properties required to open the observed gap. We then apply the formalism to the case of one single planet on a circular or eccentric orbit. We then consider multi-planetary systems: two and three equal-mass planets on circular orbits and two equal-mass planets on eccentric orbits in a packed configuration. As a final step, we compare each couple of values (M-p, a(p)), derived from the dynamical analysis of single and multiple planetary models, with the detection limits obtained with SPHERE. Results. For one single planet on a circular orbit we obtain conclusive results that allow us to exclude such a hypothesis since in most cases this configuration requires massive planets which should have been detected by our observations. Unsatisfactory is also the case of one single planet on an eccentric orbit for which we obtained high masses and/or eccentricities which are still at odds with observations. Introducing multi planetary architectures is encouraging because for the case of three packed equal-mass planets on circular orbits we obtain quite low masses for the perturbing planets which would remain undetected by our SPHERE observations. The case of two equal-mass planets on eccentric orbits is also of interest since it suggests the possible presence of planets with masses lower than the detection limits and with moderate eccentricity. Our results show that the apparent lack of planets in gaps between double belts could be explained by the presence of a system of two or more planets possibly of low mass and on eccentric orbits whose sizes are below the present detection limits.
|
|
| 7. |
- Rodet, L., et al.
(författare)
-
Dynamical masses of M-dwarf binaries in young moving groups I. The case of TWA22 and GJ2060
- 2018
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 618
-
Tidskriftsartikel (refereegranskat)abstract
- Context. Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages (<200 Myr) and in the low-mass regime (< 1 M-circle dot). GJ 2060AB and TWA22AB are two rare astrometric M-dwarf binaries, respectively members of the AB Doradus (AB Dor) and Beta Pictoris (beta Pic) moving groups. As their dynamical mass can be measured to within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups.Aims. We provide the first dynamical mass measurement of GJ 2060 and a refined measurement of the total mass of TWA22. We also characterize the atmospheric properties of the individual components of GJ 2060 that can be used as inputs to the evolutionary models.Methods. We used NaCo and SPHERE observations at VLT and archival Keck/NIRC2 data to complement the astrometric monitoring of the binaries. We combined the astrometry with new HARPS radial velocities (RVs) and FEROS RVs of GJ 2060. We used a Markov chain MonteCarlo (MCMC) module to estimate posteriors on the orbital parameters and dynamical masses of GJ 2060AB and TWA22AB from the astrometry and RVs. Complementary data obtained with the integral field spectrograph VLT /SINFONI were gathered to extract the individual near-infrared (1.1-2.5 mu m) medium-resolution (R similar to 1500 2000) spectra of GJ 2060A and B. We compared the spectra to those of known objects and to grids of BT-SETTL model spectra to infer the spectral type, bolometric luminosities, and temperatures of those objects.Results. We find a total mass of 0 : 18 +/- 0 : 02 M-circle dot for TWA22, which is in good agreement with model predictions at the age of the fi Pic moving group. We obtain a total mass of 1 : 09 +/- 0 : 10 M-circle dot for GJ 2060. We estimate a spectral type of M1 +/- 0.5, L/L-circle dot = -1.20 +/- 0.05 dex, and T-eff = 3700 +/- 100 K for GJ 2060 A. The B component is a M3 +/- 0 : 5 dwarf with L/L-circle dot = 1.63 +/- 0.05 dex and T-eff = 3400 +/- 100 K. The dynamical mass of GJ 2060AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an AB Dor age in the range 50-150 Myr. It is 10%-20% (1-2 sigma, depending on the assumed age) above the model's predictions, corresponding to an underestimation of 0.10-0.20 M fi. Coevality suggests a young age for the system (similar to 50 Myr) according to most evolutionary models.Conclusions. TWA22 validates the predictions of recent evolutionary tracks at similar to 20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the stars hosts a tight companion. Alternatively, this would confirm the model's tendency to underestimate the mass of young low-mass stars.
|
|
| 8. |
- Zurlo, A., et al.
(författare)
-
Orbital and dynamical analysis of the system around HR 8799 New astrometric epochs from VLT/SPHERE and LBT/LUCI
- 2022
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 666
-
Tidskriftsartikel (refereegranskat)abstract
- Context. HR 8799 is a young planetary system composed of four planets and a double debris belt. Being the first multi-planetary system discovered with the direct imaging technique, it has been observed extensively since 1998. This wide baseline of astrometric measurements, counting over 50 observations in 20 years, permits a detailed orbital and dynamical analysis of the system.Aims. To explore the orbital parameters of the planets, their dynamical history, and the planet-to-disk interaction, we made follow-up observations of the system during the VLT/SPHERE guaranteed time observation program. We obtained 21 observations, most of them in favorable conditions. In addition, we observed HR 8799 with the instrument LUCI at the Large Binocular Telescope (LBT).Methods. All the observations were reduced with state-of-the-art algorithms implemented to apply the spectral and angular differential imaging method. We re-reduced the SPHERE data obtained during the commissioning of the instrument and in three open-time programs to have homogeneous astrometry. The precise position of the four planets with respect to the host star was calculated by exploiting the fake negative companions method. We obtained an astrometric precision of the order of 6 mas in the worst case and 1 mas in the best case. To improve the orbital fitting, we also took into account all of the astrometric data available in the literature. From the photometric measurements obtained in different wavelengths, we estimated the masses of the planets following the evolutionary models.Results. We obtained updated parameters for the orbits with the assumption of coplanarity, relatively small eccentricities, and periods very close to the 2:1 resonance. We also refined the dynamical mass of each planet and the parallax of the system (24.49 ± 0.07 mas), which overlap with the recent Gaia eDR3/DR3 estimate. Hydrodynamical simulations suggest that inward migration of the planets caused by the interaction with the disk might be responsible for the planets being locked in resonance. We also conducted detailed N-body simulations indicating possible positions of a putative fifth planet with a mass below the present detection limits of ≃3 MJup.
|
|
| 9. |
- Chauvin, G., et al.
(författare)
-
Investigating the young solar system analog HD 95086 A combined HARPS and SPHERE exploration
- 2018
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 617
-
Tidskriftsartikel (refereegranskat)abstract
- Context. HD 95086 (A8V, 17 Myr) hosts a rare planetary system for which a multi-belt debris disk and a giant planet of 4-5 Mjup have been directly imaged.Aims. Our study aims to characterize the global architecture of this young system using the combination of radial velocity and direct imaging observations. We want to characterize the physical and orbital properties of HD 95086 b, search for additional planets at short and wide orbits and image the cold outer debris belt in scattered light.Methods. We used HARPS at the ESO 3.6 m telescope to monitor the radial velocity of HD 95086 over two years and investigate the existence of giant planets at less than 3 au orbital distance. With the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE at VLT, we imaged the faint circumstellar environment beyond 10 au at six epochs between 2015 and 2017.Results. We do not detect additional giant planets around HD 95086. We identify the nature (bound companion or background contaminant) of all point-like sources detected in the IRDIS field of view. None of them correspond to the ones recently discovered near the edge of the cold outer belt by ALMA. HD 95086 b is resolved for the first time in J-band with IFS. Its near-infrared spectral energy distribution is well fitted by a few dusty and/or young L7-L9 dwarf spectral templates. The extremely red 1-4 mu m spectral distribution is typical of low-gravity objects at the L/T spectral type transition. The planet's orbital motion is resolved between January 2015 and May 2017. Together with past NaCo measurements properly re-calibrated, our orbital fitting solutions favor a retrograde low to moderate-eccentricity orbit e = 0.2(-0.2)(+0.3), with a semi-major axis similar to 52 au corresponding to orbital periods of similar to 288 yr and an inclination that peaks at i = 141 degrees, which is compatible with a planet-disk coplanar configuration. Finally, we report the detection in polarimetric differential imaging of the cold outer debris belt between 100 and 300 au, consistent in radial extent with recent ALMA 1.3 mm resolved observations.
|
|
| 10. |
- Lazzoni, C., et al.
(författare)
-
The search for disks or planetary objects around directly imaged companions : a candidate around DH Tauri B
- 2020
-
Ingår i: Astronomy and Astrophysics. - : EDP Sciences. - 0004-6361 .- 1432-0746. ; 641
-
Tidskriftsartikel (refereegranskat)abstract
- Context. In recent decades, thousands of substellar companions have been discovered with both indirect and direct methods of detection. While the majority of the sample is populated by objects discovered using radial velocity and transit techniques, an increasing number have been directly imaged. These planets and brown dwarfs are extraordinary sources of information that help in rounding out our understanding of planetary systems.Aims. In this paper, we focus our attention on substellar companions detected with the latter technique, with the primary goal of investigating their close surroundings and looking for additional companions and satellites, as well as disks and rings. Any such discovery would shed light on many unresolved questions, particularly with regard to their possible formation mechanisms.Methods. To reveal bound features of directly imaged companions, whether for point-like or extended sources, we need to suppress the contribution from the source itself. Therefore, we developed a method based on the negative fake companion technique that first estimates the position in the field of view (FoV) and the flux of the imaged companion with high precision, then subtracts a rescaled model point spread function (PSF) from the imaged companion, using either an image of the central star or another PSF in the FoV. Next it performs techniques, such as angular differential imaging, to further remove quasi-static patterns of the star (i.e., speckle contaminants) that affect the residuals of close-in companions.Results. After testing our tools on simulated companions and disks and on systems that were chosen ad hoc, we applied the method to the sample of substellar objects observed with SPHERE during the SHINE GTO survey. Among the 27 planets and brown dwarfs we analyzed, most objects did not show remarkable features, which was as expected, with the possible exception of a point source close to DH Tau B. This candidate companion was detected in four different SPHERE observations, with an estimated mass of ~1MJup, and a mass ratio with respect to the brown dwarf of 1∕10. This binary system, if confirmed, would be the first of its kind, opening up interesting questions for the formation mechanism, evolution, and frequency of such pairs. In order to address the latter, the residuals and contrasts reached for 25 companions in the sample of substellar objects observed with SPHERE were derived. If the DH Tau Bb companion is real, the binary fraction obtained is ~7%, which is in good agreement with the results obtained for field brown dwarfs.Conclusions. While there may currently be many limitations affecting the exploration of bound features to directly imaged exoplanets and brown dwarfs, next-generation instruments from the ground and space (i.e., JWST, ELT, and LUVOIR) will be able to image fainter objects and, thus, drive the application of this technique in upcoming searches for exo-moons and circumplanetary disks.
|
|
