White Dwarf Phase Research Articles

Main sequence dynamo magnetic fields emerging in the white dwarf phase

Main sequence dynamo magnetic fields emerging in the white dwarf phase

Tidal dissipation in evolved low- and intermediate-mass stars

As the observed occurrence for planets or stellar companions orbiting low- and intermediate-mass evolved stars is increasing, so is the importance of understanding and evaluating the strength of their interactions. This is important for the further evolution of both our own Earth-Sun system and most of the observed exoplanetary systems. One of the most fundamental mechanisms behind this interaction is the tidal dissipation in these stars, as it is one of the engines of the orbital and rotational evolution of star-planet and star-star systems. This article builds upon previous works that studied the evolution of the tidal dissipation along the pre-main sequence and the main sequence of low- and intermediate-mass stars and found a strong link between the structural and rotational evolution of stars and tidal dissipation. This article provides, for the first time, a complete picture of tidal dissipation along the entire evolution of low- and intermediate-mass stars, including the advanced phases of evolution. Using stellar evolutionary models, the internal structure of the star was computed from the pre-main sequence all the way up to the white dwarf phase for stars with initial masses between 1 and 4 M$_ Using this internal structure, the tidal dissipation was computed along the entire stellar evolution. Tidal dissipation was separated into two components: the dissipation of the equilibrium (non-wave-like) tide and the dissipation of the dynamical (wave-like) tide. For evolved stars, the dynamical tide is constituted by progressive internal gravity waves. The evolution of the tidal dissipation was investigated for both the equilibrium and dynamical tides, and the results were compared. The significance of both the equilibrium and dynamical tide dissipation becomes apparent within distinct domains of the parameter space. The dissipation of the equilibrium tide is dominant when the star is large or the companion is far from the star. Conversely, the dissipation of the dynamical tide is important when the star is small or the companion is close to the star. The size and location of these domains depend on the masses of both the star and the companion, as well as on the evolutionary phase. Both the equilibrium and the dynamical tides are important in evolved stars, and therefore both need to be taken into account when studying the tidal dissipation in evolved stars and the evolution of the planetary and/or stellar companions orbiting them.

Planetary Engulfment Prognosis within the ρ CrB System

Exoplanets have been detected around stars at various stages of their lives, ranging from young stars emerging from formation to the latter stages of evolution, including white dwarfs and neutron stars. Post-main-sequence stellar evolution can result in dramatic, and occasionally traumatic, alterations to the planetary system architecture, such as tidal disruption of planets and engulfment by the host star. The ρ CrB system is a particularly interesting case of advanced main-sequence evolution, due to the relative late age and brightness of the host star, its similarity to solar properties, and the harboring of four known planets. Here, we use stellar evolution models to estimate the expected trajectory of the stellar properties of ρ CrB, especially over the coming 1.0–1.5 billion yr as it evolves off the main sequence. We show that the inner three planets (e, b, and c) are engulfed during the red giant phase and asymptotic giant branch, likely destroying those planets via either evaporation or tidal disruption at the fluid-body Roche limit. The outer planet, planet d, is briefly engulfed by the star several times toward the end of the asymptotic giant branch, but the stellar mass loss and subsequent changing planetary orbit may allow the survival of the planet into the white dwarf phase of the stellar evolution. We discuss the implications of this outcome for similar systems and describe the consequences for planets that may lie within the habitable zone of the system.

Editorial: Planetary nebulae as tools for astrophysics

Planetary nebulae (PNe) represent the ejected outer layers of an evolved low-to-intermediate mass 11star, and consist of ionized, neutral, and molecular components, along with dust. In total, the nebular 12 material contains the stellar envelope that was ejected during the AGB phase or through binary 13 interactions with a companion. The PN central star is formed by the former core of the AGB or RGB 14 progenitor of the PNe, which is now contracting towards the white dwarf phase. The PN central star 15 is a post-AGB object which has evolved crossing the HR diagram, and has an effective temperature 16 of 28,000 K up to more than 250,000 K. It emits a large fraction of its light as UV photons that ionize 17 the envelope. PNe are mainly studied through imaging, photometry and low-and high-resolution 18 spectroscopy techniques which allow determination of the physical conditions and chemical 19 abundances in the nebula. Through this kind of analysis it is possible to study the physical processes 20 in ionized plasmas, the evolution and characteristics of the central star, along with the mechanisms of 21 ejection of the nebula and the evolution of the interstellar medium (ISM) -in external galaxies as 22 well as in the Milky Way -which will be enriched with elements produced in the stellar 23 nucleosynthesis and brought to the stellar surface via several dredge-up episodes. Further, an accurate 24 census of PNe provides an important constraint on a galaxy's star formation history.

A New Third Planet and the Dynamical Architecture of the HD 33142 Planetary System* *Based on observations collected at the European Organization for Astronomical Research in the Southern Hemisphere under ESO programs 60.A-9700, 60.A-9036, 097.C-0090, 0100.C-0414, 0101.C-0232, 0102.C-0338, and MPG programs 088.C-0892, 099.A-9009, 0100.A-9006.

Based on recently-taken and archival HARPS, FEROS, and HIRES radial velocities (RVs), we present evidence for a new planet orbiting the first ascent red giant star HD 33142 (with an improved mass estimate of M1.52 ± 0.03 M ⊙), already known to host two planets. We confirm the Jovian-mass planets HD 33142b and c, with periods of P b = days and P c = days and minimum dynamical masses of = M Jup and = MJup, respectively. Furthermore, our periodogram analysis of the precise RVs shows strong evidence for a short-period Doppler signal in the residuals of a two-planet Keplerian fit, which we interpret as a third, Saturn-mass planet with = M Jup in a close-in orbit with an orbital period of P d = days. We study the dynamical behavior of the three-planet system configuration with an N-body integration scheme, finding it to be long-term stable with the planets alternating between low and moderate eccentricity episodes. We also perform N-body simulations, including stellar evolution and second-order dynamical effects such as planet–stellar tides and stellar mass loss on the way to the white dwarf phase. We find that planets HD 33142b, c, and d are likely to be engulfed near the tip of the red giant branch phase due to tidal migration. These results make the HD 33142 system an essential benchmark for planet population statistics of the multiple-planet systems found around evolved stars.

Mathematical encoding within multiresonant planetary systems as SETI beacons

ABSTRACT How might an advanced alien civilization manipulate the orbits within a planetary system to create a durable signpost that communicates its existence? While it is still debated whether such a purposeful advertisement would be prudent and wise, we propose that mean-motion resonances between neighbouring planets – with orbital periods that form integer ratios – could in principle be used to encode simple sequences that one would not expect to form in nature. In this letter, we build four multiresonant planetary systems and test their long-term orbital stability. The four systems each contain six or seven planets and consist of (i) consecutive integers from 1 to 6; (ii) prime numbers from 2 to 11; (iii) the Fibonacci sequence from 1 to 13; and (iv) the Lazy Caterer sequence from 1 to 16. We built each system using N-body simulations with artificial migration forces. We evaluated the stability of each system over the full 10 Gyr integration of the Sun’s main-sequence phase. We then tested the stability of these systems for an additional 10 Gyr, during and after post-main-sequence evolution of the central stars (assumed to be Sun-like) to their final, white dwarf phase. The only system that was destabilized was the consecutive integer sequence (system i). The other three sequences therefore represent potential Search for Extraterrestrial Intelligence (SETI) beacons.

Study of UV bright sources in globular cluster NGC 4590 using Ultraviolet Imaging Telescope (UVIT) observations

ABSTRACT We have studied ultraviolet (UV) bright sources in the Galactic globular cluster (GGC) NGC 4590 using Ultraviolet Imaging Telescope onboard the AstroSat satellite. Using UV-optical colour–magnitude diagrams (CMDs), we have identified and characterized the sources of different evolutionary stages i.e. blue horizontal branch stars (BHBs), extremely blue horizontal branch stars (EHBs), blue straggler stars (BSs), variable stars, etc. We estimated effective temperature (Teff), gravity (log (g)), luminosity (Lbol), and hence the radius (R) of these hot stars by fitting spectral energy distribution with the help of stellar atmosphere models. Two new far-UV bright cluster member stars situated near the core of the cluster have been detected; one of them is an EHB star and the other one is either in its post-blue hook evolutionary phase or in white dwarf phase. The evolutionary status of all the hot stars, identified in the cluster, has been investigated by using various evolutionary models. We find the massive and younger BSs are concentrated at the centre of the cluster whereas the older and less massive BSs are distributed throughout the cluster. The BSs normalized radial distribution seems to be bi-modal with a minimum located at rmin = 4.3 rc. We calculated A+ parameter of the cluster which is obtained using cumulative normalized radial distribution of horizontal branch stars (HBs) and BSs. We measured this value up to half-mass radius of the cluster to be +0.13, which indicates that NGC 4590 is one of the youngest clusters among dynamically intermediate age GGCs with a dynamical age of 0.423 ± 0.096 Gyr.

Planetary magnetosphere evolution around post-main-sequence stars

ABSTRACT Accompanying the mounting detections of planets orbiting white dwarfs and giant stars are questions about their physical history and evolution, particularly regarding detectability of their atmospheres and potential for habitability. Here we determine how the size of planetary magnetospheres evolves over time from the end of the main sequence through to the white dwarf phase due to the violent winds of red giant and asymptotic giant branch stars. By using a semi-analytic prescription, we investigate the entire relevant phase space of planet type, planet orbit, and stellar host mass ($1\!-\!7\, {\rm M}_{\odot }$). We find that a planetary magnetosphere will always be quashed at some point during the giant branch phases unless the planet’s magnetic field strength is at least two orders of magnitude higher than Jupiter’s current value. We also show that the time variation of the stellar wind and density generates a net increase in wind ram pressure and does not allow a magnetosphere to be maintained at any time for field strengths less than 10−5 T (0.1 G). This lack of protection hints that currently potentially habitable planets orbiting white dwarfs would have been previously inhospitable.

Host-star and exoplanet compositions: a pilot study using a wide binary with a polluted white dwarf

ABSTRACT Planets and stars ultimately form out of the collapse of the same cloud of gas. Whilst planets, and planetary bodies, readily loose volatiles, a common hypothesis is that they retain the same refractory composition as their host star. This is true within the Solar system. The refractory composition of chondritic meteorites, Earth, and other rocky planetary bodies are consistent with solar, within the observational errors. This work aims to investigate whether this hypothesis holds for exoplanetary systems. If true, the internal structure of observed rocky exoplanets can be better constrained using their host star abundances. In this paper, we analyse the abundances of the K-dwarf, G200-40, and compare them to its polluted white dwarf companion, WD 1425+540. The white dwarf has accreted planetary material, most probably a Kuiper belt-like object, from an outer planetary system surviving the star’s evolution to the white dwarf phase. Given that binary pairs are chemically homogeneous, we use the binary companion, G200-40, as a proxy for the composition of the progenitor to WD 1425+540. We show that the elemental abundances of the companion star and the planetary material accreted by WD 1425+540 are consistent with the hypothesis that planet and host-stars have the same true abundances, taking into account the observational errors.

Heavy elements nucleosynthesis on accreting white dwarfs: building seeds for the p-process

ABSTRACT The origin of the proton-rich trans-iron isotopes in the Solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor’s rapidly accreting white dwarf (RAWD) phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf (WD) models with initial masses ≥ 0.85 $\, \mathrm{M}_\odot$ using the stellar code mesa The near-surface layers of the 1, 1.26, 1.32 and 1.38 $\, \mathrm{M}_\odot$ models are most representative of the regions in which the bulk of the p nuclei are produced during SNIa explosions, and for these models we also calculate the neutron-capture nucleosynthesis in the external layers. Contrary to previous RAWD models at lower mass, we find that the H-shell flashes are the main site of n-capture nucleosynthesis. We find high neutron densities up to several 1015 cm−3 in the most massive WDs. Through the recurrence of the H-shell flashes these intermediate neutron densities can be sustained effectively for a long time leading to high-neutron exposures with a strong production up to Pb. Both the neutron density and the neutron exposure increase with increasing the mass of the accreting WD. Finally, the SNIa nucleosynthesis is calculated using the obtained abundances as seeds. We obtain solar to supersolar abundances for p-nuclei with A > 96. Our models show that SNIa are a viable p-process production site.

Linking the formation and fate of exo-Kuiper belts within Solar system analogues

ABSTRACT Escalating observations of exo-minor planets and their destroyed remnants both passing through the Solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here, we explore how the structure of a 40–1000 au annulus of planetesimals orbiting inside of a Solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt’s evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals’ gravitational responses to analogues of the four Solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of ≳100 km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close (≪103 au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.

Constraining planet formation around 6–8 M⊙ stars

ABSTRACT Identifying planets around O-type and B-type stars is inherently difficult; the most massive known planet host has a mass of only about $3\, \mathrm{M}_{\odot }$. However, planetary systems which survive the transformation of their host stars into white dwarfs can be detected via photospheric trace metals, circumstellar dusty and gaseous discs, and transits of planetary debris crossing our line of sight. These signatures offer the potential to explore the efficiency of planet formation for host stars with masses up to the core-collapse boundary at $\approx 8\, \mathrm{M}_{\odot }$, a mass regime rarely investigated in planet formation theory. Here, we establish limits on where both major and minor planets must reside around $\approx 6\rm {-}8\, \mathrm{M}_{\odot }$ stars in order to survive into the white dwarf phase. For this mass range, we find that intact terrestrial or giant planets need to leave the main sequence beyond approximate minimum star–planet separations of, respectively, about 3 and 6 au. In these systems, rubble pile minor planets of radii 10, 1.0, and 0.1 km would have been shorn apart by giant branch radiative YORP spin-up if they formed and remained within, respectively, tens, hundreds, and thousands of au. These boundary values would help distinguish the nature of the progenitor of metal pollution in white dwarf atmospheres. We find that planet formation around the highest mass white dwarf progenitors may be feasible, and hence encourage both dedicated planet formation investigations for these systems and spectroscopic analyses of the highest mass white dwarfs.

Asteroseismology of evolved stars to constrain the internal transport of angular momentum

Context. The internal characteristics of stars, such as their core rotation rates, are obtained via asteroseismic observations. A comparison of core rotation rates found in this way with core rotation rates as predicted by stellar evolution models demonstrate a large discrepancy. This means that there must be a process of angular momentum transport missing in the current theory of stellar evolution. A new formalism was recently proposed to fill in for this missing process, which has the Tayler instability as its starting point (by Fuller et al. 2019, MNRAS, 485, 3661, hereafter referred to as “Fuller-formalism”). Aims. We investigate the effect of the Fuller-formalism on the internal rotation of stellar models with an initial mass of 2.5 M⊙. Methods. Stellar evolution models, including the Fuller-formalism, of intermediate-mass stars were calculated to make a comparison between asteroseismically obtained core rotation rates in the core He burning phase and in the white dwarf phase. Results. Our main results show that models including the Fuller-formalism can match the core rotation rates obtained for the core He burning phases. However, these models are unable to match the rotation rates obtained for white dwarfs. When we exclude the Fuller-formalism at the end of the core He burning phase, the white dwarf rotation rates of the models match the observed rates. Conclusions. We conclude that in the present form, the Fuller-formalism cannot be the sole solution for the missing process of angular momentum transport in intermediate-mass stars.

The s process in rotating low-mass AGB stars

Aims. We investigate the s-process during the AGB phase of stellar models whose cores are enforced to rotate at rates consistent with asteroseismology observations of their progenitors and successors. Methods. We calculated new 2 M⊙, Z = 0.01 models, rotating at 0, 125, and 250 km s−1 at the start of main sequence. An artificial, additional viscosity was added to enhance the transport of angular momentum in order to reduce the core rotation rates to be in agreement with asteroseismology observations. We compared rotation rates of our models with observed rotation rates during the MS up to the end of core He burning, and the white dwarf phase. Results. We present nucleosynthesis calculations for these rotating AGB models that were enforced to match the asteroseismic constraints on rotation rates of MS, RGB, He-burning, and WD stars. In particular, we calculated one model that matches the upper limit of observed rotation rates of core He-burning stars and we also included a model that rotates one order of magnitude faster than the upper limit of the observations. The s-process production in both of these models is comparable to that of non-rotating models. Conclusions. Slowing down the core rotation rate in stars to match the above mentioned asteroseismic constraints reduces the rotationally induced mixing processes to the point that they have no effect on the s-process nucleosynthesis. This result is independent of the initial rotation rate of the stellar evolution model. However, there are uncertainties remaining in the treatment of rotation in stellar evolution, which need to be reduced in order to confirm our conclusions, including the physical nature of our approach to reduce the core rotation rates of our models, and magnetic processes.

Photoionization modelling of the Wolf-Rayet planetary nebulae IC 4663 based on multiwavelength observations

Observations of the central star planetary nebulae (CSPN) IC 4663 revealed its spectral features similar to WN sequence of massive Wolf-Rayet (WR) stars. Because CSPN IC 4663 spectrum is dominated by broad He II and N V emission lines, it is classified as [WN3] spectral type. The similarity on the central stars composition of O(He) and IC 4663 provides evidence for evolutionary sequence [WN] → O(He). Evolution of CSPN IC 4663 followed evolution low-mass star, but experienced born-again scenario when it enters white dwarf phase. So it goes back to Asymptotic Giant Branch (AGB) phase. To investigate central star properties, the nebula and its evolutionary status, we constructed the spectral energy distribution (SED) model by using photoionization code CLOUDY based on data set from optic to far-IR wavelength (∼3 – 100 μm) and include the Postdam Wolf-Rayet (PoWR) stellar model atmosphere. The comparison between the nebular abundances from photoionization model result and AGB model gives initial 3M⊙ star with Z = 0.02. The nebula ionized mass derived from the photoionization model also consistent with the prediction from AGB model.

Constraining transport of angular momentum in stars

Context. Transport of angular momentum has been a challenging topic within the stellar evolution community, even more since the recent asteroseismic surveys. All published studies on rotation using asteroseismic observations show a discrepancy between the observed and calculated rotation rates, indicating there is an undetermined process of angular momentum transport active in these stars. Aims. We aim to constrain the efficiency of this process by investigating rotation rates of 2.5 M⊙ stars. Methods. First, we investigated whether the Tayler-Spruit dynamo could be responsible for the extra transport of angular momentum for stars with an initial mass of 2.5 M⊙. Then, by computing rotating models including a constant additional artificial viscosity, we determined the efficiency of the missing process of angular momentum transport by comparing the models to the asteroseismic observations of core helium burning stars. Parameter studies were performed to investigate the effect of the stellar evolution code used, initial mass, and evolutionary stage. We evolved our models into the white dwarf phase, and provide a comparison to white dwarf rotation rates. Results. The Tayler-Spruit dynamo is unable to provide enough transport of angular momentum to reach the observed values of the core helium burning stars investigated in this paper. We find that a value for the additional artificial viscosity νadd around 107 cm2 s−1 provides enough transport of angular momentum. However, the rotational period of these models is too high in the white dwarf phase to match the white dwarf observations. From this comparison we infer that the efficiency of the missing process must decrease during the core helium burning phase. When excluding the νadd during core helium burning phase, we can match the rotational periods of both the core helium burning stars and white dwarfs.

Improving the distances of post-AGB objects in the Milky Way

Post-AGB (PAGB) stars are short-lived, low-intermediate mass objects transitioning from the asymptotic giant branch (AGB) to the white dwarf (WD) phase. These objects are characterised by a constant, core-mass dependent luminosity and a large infrared excess from the dusty envelope ejected at the top of the AGB. PAGB stars provide insights into the evolution of their direct descendants, planetary nebulae (PNe). Calculation of physical characteristics of PAGB are dependent on accurately determined distances scarcely available in the literature. Using the Torun catalogue for PAGB objects, supplemented with archival data, we have determined distances to the known population of Galactic PAGB stars. This is by modelling their spectral energy distributions (SED) with black bodies and numerically integrating over the entire wavelength range to determine the total integrated object flux. For most PAGB stars we assumed their luminosities are based on their positional characteristics and stellar evolution models. RV Tauri stars however are known to follow a period-luminosity relation (PLR) reminiscent of type-2 Cepheids. For these variable PAGB stars we determined their luminosities via the PLR and hence their distances. This allows us to overcome the biggest obstacle to characterising these poorly understood objects that play a vital part in Galactic chemical enrichment.

Full-lifetime simulations of multiple unequal-mass planets across all phases of stellar evolution

We know that planetary systems are just as common around white dwarfs as around main sequence stars. However, self-consistently linking a planetary system across these two phases of stellar evolution through the violent giant branch poses computational challenges, and previous studies restricted architectures to equal-mass planets. Here, we remove this constraint and perform over 450 numerical integrations over a Hubble time (14 Gyr) of packed planetary systems with unequal-mass planets. We characterize the resulting trends as a function of planet order and mass. We find that intrusive radial incursions in the vicinity of the white dwarf become less likely as the dispersion amongst planet masses increases. The orbital meandering which may sustain a sufficiently dynamic environment around a white dwarf to explain observations is more dependent on the presence of terrestrial-mass planets than any variation in planetary mass. Triggering unpacking or instability during the white dwarf phase is comparably easy for systems of unequal-mass planets and systems of equal-mass planets; instabilities during the giant branch phase remain rare and require fine-tuning of initial conditions. We list the key dynamical features of each simulation individually as a potential guide for upcoming discoveries.

ANGULAR MOMENTUM TRANSPORT WITHIN EVOLVED LOW-MASS STARS

Asteroseismology of 1.0-2.0 Msun red giants by the Kepler satellite has enabled the first definitive measurements of interior rotation in both first ascent red giant branch (RGB) stars and those on the Helium burning clump. The inferred rotation rates are 10-30 days for the ~0.2Msun He degenerate cores on the RGB and 30-100 days for the He burning core in a clump star. Using the MESA code we calculate state-of-the-art stellar evolution models of low mass rotating stars from the zero-age main sequence to the cooling white dwarf (WD) stage. We include transport of angular momentum due to rotationally induced instabilities and circulations, as well as magnetic fields in radiative zones (generated by the Tayler-Spruit dynamo). We find that all models fail to predict core rotation as slow as observed on the RGB and during core He burning, implying that an unmodeled angular momentum transport process must be operating on the early RGB of low mass stars. Later evolution of the star from the He burning clump to the cooling WD phase appears to be at nearly constant core angular momentum. We also incorporate the adiabatic pulsation code, ADIPLS, to explicitly highlight this shortfall when applied to a specific Kepler asteroseismic target, KIC8366239. The MESA inlist adopted to calculate the models in this paper can be found at \url{https://authorea.com/1608/} (bottom of the document).

Long-term evolution of three-planet systems to the post-main sequence and beyond

We study the stability of systems of three giant planets orbiting 3 - 8 M_Sol stars at orbital distances of >10 au as the host star ages through the Main Sequence (MS) and well into the White Dwarf (WD) stage. Systems are stable on the MS if the planets are separated by more than ~9 Hill radii. Most systems surviving the MS will remain stable until the WD phase, although planets scattered onto small pericentres in unstable systems can be swallowed by the expanding stellar envelope when the star ascends the giant branches. Mass loss at the end of the asymptotic giant branch triggers delayed instability in many systems, leading to instabilities typically occurring at WD cooling ages of a few 100 Myr. This instability occurs both in systems that survived the star's previous evolution unscathed, and in systems that previously underwent scattering instabilities. The outcome of such instability around WDs is overwhelmingly the ejection of one of the planets from the system, with several times more ejections occurring during the WD phase than during the MS. Furthermore, few planets are scattered close to the WD, just outside the Roche limit, where they can be tidally circularised. Hence, we predict that planets in WD systems rarely dynamically evolve to become "hot Jupiters". Nor does it appear that the observed frequency of metal pollution in WD atmospheres can be entirely explained by planetesimals being destabilised following instability in systems of multiple giant planets, although further work incorporating low-mass planets and planetesimals is needed.