Observational Insights into White Dwarf Planetary Systems Laura Kathryn Rogers Institute of Astronomy University of Cambridge This dissertation is submitted for the degree of Doctor of Philosophy Jesus College August 2022 To my astronomically amazing and supportive parents, Sue and Paul. Declaration I hereby declare that except where specific reference is made to the work of others, the contents of this dissertation are original and have not been submitted in whole or in part for consideration for any other degree or qualification in this, or any other university. I state that except where specific reference is made in the text to the work of others, including any collaborators, the contents of this dissertation are the result of my own work. This dissertation contains fewer than 60,000 words including summary/abstract, tables, footnotes and appendices. Chapter 2 has been adapted from the publication: ‘Near-infrared variability in dusty white dwarfs: tracing the accretion of planetary material’ L. K. Rogers, et al., MNRAS, Volume 494, Issue 2, May 2020, Pages 2861–2874. Chapter 4 has been adapted from work submitted to MNRAS: ‘Tracing the composition of exoplanetary building blocks using seven polluted white dwarfs’ L. K. Rogers, et al. (2022, submitted). Laura Kathryn Rogers August 2022 Abstract Observational Insights into White Dwarf Planetary Systems Laura Kathryn Rogers In recent decades the number of known planets has escalated from the eight solar system planets to over 5000 exoplanets; the focus has now shifted to their characterisation. Spec- troscopic studies of white dwarfs ‘polluted’ by planetary material are a unique laboratory allowing measurements of their bulk composition. This thesis focuses on these polluted white dwarfs to observationally investigate how the material ultimately accretes onto the white dwarfs, and the inferences made about the composition of exoplanetary bodies. These white dwarfs become polluted due to the scattering of exoplanetary bodies on star grazing orbits where they tidally disrupt, producing dusty debris detectable as excess infrared emission, and then subsequently accrete onto the white dwarf. The scattering and accretion are expected to be stochastic processes with variability predicted on human time-scales. Chapter 2 reports near-infrared (JHK) monitoring campaigns of the dust emission with the UKIRT/WFCAM and NTT/SOFI. Over timescales of hours, days, months, and years no statistically significant variation is found. Chapter 3 reports spectroscopic monitoring campaigns of the metal features in the white dwarfs using SALT/HRS and Magellan/MIKE. Across more than 10 years and thousands of sinking timescales, no unambiguous statistically significant variability in the amount of material in the photospheres of the white dwarfs is found. Both results agree that the processes driving the accretion do so at a constant rate. Polluted white dwarfs with infrared emission tend to be the most heavily polluted systems. Chapter 4 presents the first composition studies from a novel programme which identifies new heavily polluted white dwarfs discovered from their infrared excess. The planetary material that polluted these seven white dwarfs are broadly consistent with rocky material, but show some compositional and geological diversity. Some of the white dwarfs appear to have accreted a fragment of a larger core-mantle differentiated body. Using oxygen budgeting, two white dwarfs are discovered to have accreted oxygen-rich material which may imply water-rich bodies. Also, evidence points towards one of the white dwarfs accreting material from two distinct planetary bodies, this challenges the commonly used assumption that there is one body in the white dwarf’s atmosphere at once. Acknowledgements I would firstly like to thank my wonderful supervisors, Amy Bonsor and Simon Hodgkin, without which I would have never made it to the end of the PhD. Thank you for the endless support and for always making time to help and encourage me during my PhD (and for putting up with the endless silly questions). I will forever be grateful for Amy as my primary advisor; she has been such a supportive, encouraging, and knowledgeable supervisor relating to both research and life. I wish I could thank Amy a million times over, she has been such an amazing role model and mentor and I will always look up to and aspire to be like her. I am incredibly grateful for funding from the Science and Technology Facilities Coun- cil which allowed me to pursue my dream of a PhD in astronomy. I would also like to acknowledge funding from the University of Cambridge and Jesus College, Cambridge. I am so grateful to my ‘unofficial’ supervisor, Siyi Xu; Siyi has supported me throughout my PhD, her willingness to teach me and help me develop as a researcher has meant the world to me and I am so grateful for her time, knowledge and patience, especially considering the 11 hour time difference! I am so thankful to numerous members of the astronomical community who have helped me throughout my PhD, but would especially like to thank: Beth Klein, Ted von Hippel, and Patrick Dufour for their help and guidance. I have loved learning about observing and have been very fortunate to gain lots of observational experience during my PhD. I am incredibly grateful to Mike Shara and Ted von Hippel for enabling me to obtain my very first data set using the HRS/SALT (reported in Chapter 3). Thank you so much for your generosity, knowledge, and advice as I have learnt infinite amounts from this project. I was also incredibly fortunate to go on an observing trip to Chile with Anais Gonneau as part of my first accepted PI proposal, not only did I learn lots from her, but we had such a wonderful time. I would also like to thank Claudia Aguilera Gómez for teaching me all she knows about observations and for letting me watch her remotely to observe. I have thoroughly enjoyed my time as a PhD student at the IoA, and have especially enjoyed the friendly and collaborative feel of the department. Thank you to everyone both past and present for making the IoA such an amazing place to study for my PhD, even if some of it was during a global pandemic! I have really appreciated the wealth of knowledge xat the IoA and want to thank Mark Wyatt, Mike Irwin, and Matt Auger for always being willing to answer my questions. I would distinctly like to thank Paul Hewett and Debbie Peterson, I will be eternally grateful for their time, help and guidance. They both made it possible and easy for me to continue and finish my PhD under difficult circumstances. My PhD has been very unconventional, beginning full time and transitioning to part time after a period of absence due to a head injury. However, the positive side of doing a PhD over 6 years meant that I got to overlap with 10 cohorts of amazing and inspiring PhD and Masters students. It has been awesome sharing my time at the IoA with each and every one of them, and I am so thankful for their support, friendships, and laughs. I would especially like to thank all those who started their PhD at the same time as me, with special thanks to my fellow polluted white dwarf enthusiast John, I learnt so much from John including the best goalkeepers to pick for fantasy league. I would also like to thank the wonderful humans from the years below who welcomed me into their years with open arms. I would especially like to thank: Amy, Sophie, Anjali, Aoife, Chiara, Matthew, and Luis, your friendship and support has been incredible. I also want to shout out to my numerous amazing office mates, the white dwarf group, and the board games crews for always keeping me sane and happy. I want to thank all my family and friends for their non-stop encouragement, care, and love throughout this journey. I have chosen to dedicate my thesis to my amazing parents, Paul and Sue, who have always shown me unconditional love and support. Thank you so much for always being there for me, for encouraging me when I didn’t believe in myself, and looking after me when I needed it most. To my amazing siblings Sarah, Matthew and Sophie, thanks for always being there for me and I am so sorry for all the times you have had to listen to me go on about astronomy! To my niece, Ella, thank you for always making me smile, I love that you love the moon more than I do! To my friends outside of my PhD, thank you so much for keeping me sane over these years and for the never ending love, laughs, and support: Rosie, Kris, Thomas, Kitty, Georgia, and George, and to all my Jesus College and Cambridge friends, thank you for making this such a memorable experience. Finally, I want to thank my incredible partner, Dave. Life has been a roller coaster with the PhD, head injury, and pandemic, but you have always been my calm constant that keeps me grounded. You make me a better person and your love and support throughout the PhD is something I am thankful for, far beyond words. Table of contents List of figures xv List of tables xxv 1 Introduction 1 1.1 Solar System . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 1.1.1 Rocky bodies in the solar system . . . . . . . . . . . . . . . . . . . 2 1.2 Exoplanets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 1.2.1 Discovery . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 1.2.2 Characterisation . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 1.3 White dwarfs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 1.3.1 Formation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11 1.3.2 White dwarf properties . . . . . . . . . . . . . . . . . . . . . . . . 13 1.3.3 Spectral classification . . . . . . . . . . . . . . . . . . . . . . . . . 14 1.4 Polluted White Dwarfs . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 1.4.1 How do the white dwarfs become polluted? . . . . . . . . . . . . . 17 1.4.2 Composition of polluting bodies . . . . . . . . . . . . . . . . . . . 18 1.4.3 Circumstellar Discs . . . . . . . . . . . . . . . . . . . . . . . . . . 21 1.4.4 The Link with the Main Sequence . . . . . . . . . . . . . . . . . . 33 1.5 This Thesis in Context . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 2 Dust Variability 37 2.1 Dust variability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38 2.1.1 UKIRT Observations . . . . . . . . . . . . . . . . . . . . . . . . . 38 2.1.2 Variability Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . 40 2.1.3 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 2.1.4 Conclusions from the UKIRT survey . . . . . . . . . . . . . . . . . 58 2.2 Dust variability . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 xii Table of contents 2.2.1 NTT Observations and data reduction . . . . . . . . . . . . . . . . 60 2.2.2 Variability Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . 60 2.2.3 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64 2.2.4 Conclusions from the NTT/SOFI survey . . . . . . . . . . . . . . . 68 3 Metal line variability 71 3.1 Observations and Data Reduction . . . . . . . . . . . . . . . . . . . . . . . 72 3.1.1 Target Selection . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72 3.1.2 Spectroscopic observations . . . . . . . . . . . . . . . . . . . . . . 72 3.2 Variability Analysis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78 3.3 Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79 3.3.1 WD 1929+011 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80 3.3.2 G29-38 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89 3.3.3 WD 0106−328 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 3.3.4 WD 0408−041 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 3.3.5 WD 1457−086 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93 3.4 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 93 3.4.1 Discussion of methods . . . . . . . . . . . . . . . . . . . . . . . . 93 3.4.2 G29-38 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95 3.4.3 WD 0106−328 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95 3.4.4 Accretion rates . . . . . . . . . . . . . . . . . . . . . . . . . . . . 96 3.4.5 White dwarf models . . . . . . . . . . . . . . . . . . . . . . . . . 98 3.4.6 The link with circumstellar variability . . . . . . . . . . . . . . . . 99 3.5 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100 4 Tracing the composition of exoplanetary building blocks 103 4.1 Observations and Data Reduction . . . . . . . . . . . . . . . . . . . . . . . 104 4.1.1 Targets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104 4.1.2 X-shooter . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107 4.1.3 HIRES . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107 4.1.4 MIKE . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108 4.2 Modelling Methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108 4.2.1 White Dwarf Parameters . . . . . . . . . . . . . . . . . . . . . . . 108 4.2.2 Line Identification . . . . . . . . . . . . . . . . . . . . . . . . . . 110 4.2.3 Abundance of planetary material . . . . . . . . . . . . . . . . . . . 110 4.2.4 Interpreting the observed abundances . . . . . . . . . . . . . . . . 114 4.3 Results from individual systems . . . . . . . . . . . . . . . . . . . . . . . 118 Table of contents xiii 4.3.1 Gaia J0006+2858 . . . . . . . . . . . . . . . . . . . . . . . . . . . 119 4.3.2 Gaia J0347+1624 . . . . . . . . . . . . . . . . . . . . . . . . . . . 119 4.3.3 Gaia J0510+2315 . . . . . . . . . . . . . . . . . . . . . . . . . . . 120 4.3.4 Gaia J0611−6931 . . . . . . . . . . . . . . . . . . . . . . . . . . . 122 4.3.5 Gaia J0644−0352 . . . . . . . . . . . . . . . . . . . . . . . . . . . 122 4.3.6 WD 1622+587 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124 4.3.7 Gaia J2100+2122 . . . . . . . . . . . . . . . . . . . . . . . . . . . 124 4.4 Comparing pollution in white dwarfs with and without detectable circumstel- lar gas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 126 4.4.1 Accretion rates . . . . . . . . . . . . . . . . . . . . . . . . . . . . 128 4.4.2 Metal abundances . . . . . . . . . . . . . . . . . . . . . . . . . . . 131 4.5 Discussion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 131 4.5.1 Effects of uncertainties on elemental abundances . . . . . . . . . . 133 4.5.2 Volatiles in exoplanetary systems . . . . . . . . . . . . . . . . . . 134 4.5.3 Evidence for the accretion of multiple bodies? . . . . . . . . . . . . 135 4.5.4 Are high Ca/Mg ratios explained by the accretion of a high proportion of refractory material? . . . . . . . . . . . . . . . . . . . . . . . . 136 4.5.5 Evidence for the formation of iron cores in exoplanetary systems . . 137 4.5.6 Magnetic field of WD 1622+587 . . . . . . . . . . . . . . . . . . . 138 4.6 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139 5 Conclusions and Future Work 141 5.1 Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141 5.2 Future Work . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 144 References 147 6 Appendix 1: Additional tables and figures from Chapter 2 159 6.1 Observation Information . . . . . . . . . . . . . . . . . . . . . . . . . . . 159 6.2 SEDs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 161 7 Appendix 2: Additional tables and figures from Chapter 4 169 7.1 Spectral Energy Distributions . . . . . . . . . . . . . . . . . . . . . . . . . 169 7.2 Spectral Lines and abundances . . . . . . . . . . . . . . . . . . . . . . . . 169 7.3 Equivalent Width Upper Limits . . . . . . . . . . . . . . . . . . . . . . . . 169 7.4 Model fits to spectral lines . . . . . . . . . . . . . . . . . . . . . . . . . . 176 7.5 Bayesian output . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 181 List of figures 1.1 Abundance ratio of element X to Mg, normalised to solar for compositions of: CI Chondrites, bulk Earth, Earth’s mantle, and Earth’s core. Data are from: McDonough (2003b); Grevesse et al. (2007). ∗ The core composition is: 2/3 core and 1/3 mantle. . . . . . . . . . . . . . . . . . . . . . . . . . . 4 1.2 The relative abundance of X/Mg for bulk Earth (open circles) and classes of carbonaceous chondrites (CM class is shown as filled circles, CV as filled squares, and CO as open triangles), normalised to CI chondrites, plotted against the log temperature at which 50% of the element would have con- densed from a gas at 10−4 bar. The shaded regions defines the area covered by the three carbonaceous chondrites classes. The volatility depletion trend in bulk Earth and carbonaceous chondrites is clear. Figure from McDonough (2003a). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 1.3 Orbital period versus planet mass for exoplanets with both measured/estimated orbital period and mass. The colours highlight the different discovery meth- ods. Figure using data from the NASA exoplanet archive1. . . . . . . . . . 8 1.4 A diagram highlighting the evolution of a solar mass star from the zero age main sequence (ZAMS) to the white dwarf phase. Figure from Carroll & Ostlie (2017). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 1.5 Example spectra of white dwarf spectral types from LAMOST. Figure from Gentile Fusillo et al. (2015). . . . . . . . . . . . . . . . . . . . . . . . . . 16 1.6 A diagram to demonstrate the polluted white dwarf systems, taken from Jura & Young (2014). The composition of the pollution in the atmosphere will be discussed in Section 1.4.2 and the planetary debris will be discussed in Section 1.4.3. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 1.7 Bar charts showing the pollutant mass fractions of: O, Mg, Si, Fe, C, for 19 polluted white dwarfs. The bulk compositions of the Earth and comet Halley are shown on the left for comparison. Figure from Xu & Bonsor (2021). . . 22 xvi List of figures 1.8 The spectral energy distribution of G29-38. It shows the first infrared excess discovered around a polluted white dwarf (Zuckerman & Becklin, 1987). The data was taken using the NASA Infrared Telescope Facility. An excess at near to mid-infrared wavelengths is clear. . . . . . . . . . . . . . . . . . 23 1.9 Spitzer/IRS spectrum of the dust around G29-38. The observed spectrum was divided by a blackbody of temperature 930 K to give the black data points. The best fitting linear combination of 12 minerals is plotted as the red dashed line, to show the potential composition of the dust. The contribution of each mineral is shown by the coloured lines (Reach et al., 2009). The composition of this dusty material should match the composition of the material polluting the star. Water ice, amorphous silicates and forsterites all have distinguishable features in the mid-IR. . . . . . . . . . . . . . . . . . 24 1.10 The Ca II emission line profiles for WD 1226+110 from Gänsicke et al. (2006). The double peaked profiles are clear, arising due to Keplerian rotation. There are clear asymmetries in the red and blue Doppler shifted lines, indicating the disc is eccentric. . . . . . . . . . . . . . . . . . . . . 26 1.11 Figure from Wyatt et al. (2014) showing model predictions for accretion rate over time for 20 stars modelled to be similar to G29-38. Changes in the accretion rate can be up to an order of magnitude over timescales of less than a year. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 1.12 Figure from Swan et al. (2020) showing the mid-infrared variability observed in a sample of ∼ 40 white dwarfs using the Spitzer space telescope. The largest flux changes seem to occur over longer baselines with the largest changes being those with circumstellar gaseous discs. . . . . . . . . . . . 29 1.13 A diagram from Rafikov (2011a) illustrating an interacting dust disc (black) and gaseous disc (grey). The radial extent of the dust disc is defined on the inner edge by the sublimation radius of solids and the outer edge by the tidal disruption radius of the white dwarf. The gas disc has the same radial extent as the dust, with a reservoir close to the white dwarf, feeding the atmosphere with metallic gas. The movement of high metallicity materials is demonstrated by the arrows. . . . . . . . . . . . . . . . . . . . . . . . . . 32 1.14 This figure is taken from Bonsor et al. (2017). It demonstrates the accretion rate plotted against the temperature of the white dwarfs. The plot has been subdivided into four regions to explain the lack of observations of infrared excesses around polluted white dwarfs. . . . . . . . . . . . . . . . . . . . . 33 List of figures xvii 2.1 A normalised histogram showing the distribution of photometrically cor- rected magnitudes from the LIGHTCURVES software using all dithered stacked frames for WD 1145+017. Each dithered stacked frame contains 5×10 second exposures. This white dwarf has 102 measurements in the K band over 782 days. The Gaussian curve shows the median best-fitting Gaus- sian, with a median magnitude of 17.467, and a median standard deviation of 0.068 mags. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42 2.2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45 2.2 The median observed magnitudes in the (a) J, (b) H and (c) K bands plotted against the median standard deviation of this magnitude for all white dwarfs and field stars. The two different observing modes, with 5 and 10 second frame times, are analysed separately. The density plots on the left hand side show the parameters for all field stellar objects in the field of view, there are of the order 104 field stars in each plot. The blue points show the white dwarfs, with the objects discussed in Section 2.1.2 labelled. The right hand plots show the median magnitude and standard deviation with 16th and 84th percentile confidence levels for field stars, in bins of width 0.5 mag. This forms a distribution about that median which should not comprise contaminant variable stars. The white dwarfs are over-plotted in black with the same values as in the left hand plots, but also including the 16th and 84th percentiles. This plot can be used to distinguish white dwarfs which are variable from those that follow the field star distribution. . . . . . . . . . . 46 2.3 A lightcurve in the J, H and K bands for WD 1145+017. At each observation date a number of observations were taken in the J, H and K bands, totalling 28 minutes. The dotted lines show the mean magnitude for each band calculated from the first observation date, 02/05/15. This highlights how there was an apparent dip on 09/05/17 in the flux in the J and H bands attributed to dust transit features. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47 2.4 A cumulative distribution showing the level of K band variability that can be ruled out for the white dwarfs. The blue line shows the median cumulative distribution, with the contours as the cumulative distribution on the upper and lower limits (16th and 84th percentiles). For 90% of the sample, the K band flux varied by less than 10%. . . . . . . . . . . . . . . . . . . . . . . 50 xviii List of figures 2.5 (a) The maximum K band median variability of each white dwarf that would not have been detected over the time-scale of the survey. (b) The maximum dust variability of each white dwarf that would not have been detected over the time-scale of the survey. The variability constraint of the dust component to the K band flux was calculated using equation 2.1. . . . . . . . . . . . . 53 2.6 Predictions for the maximum median variability at 4.5 µm that would have escaped detection based on the maximum median K band variability not detected with this UKIRT survey. The arrows represent that each axis is a constraint on the level of variability. Most objects had Spitzer 4.5 µm data, but for four objects WISE data was used as no Spitzer data was available. Swan et al. (2019b, 2020) assessed the level of Spitzer and WISE mid infrared variability at levels above these limits, future analyses that combine the results from these studies are crucial for understanding these dusty white dwarfs. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54 2.7 A lightcurve of WD 0408−041 showing the UKIRT K band data, and the Spitzer/WISE data at which Farihi et al. (2018) observed a drop in the mid- infrared flux. The grey regions represent the 1 and 3σ median variability constraint for the K band, and then scaled up to Spitzer wavelengths using equations 2.1 and 2.2 (assuming no colour change). The blue region shows the period of time with UKIRT observations. Figure adapted from Farihi et al. (2018). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55 2.8 A histogram showing the characteristic survey sampling to search for vari- ability. The plot shows the distribution of all time-scales between frames of a given white dwarf, for all white dwarfs. This demonstrates that for all white dwarfs the survey was sensitive to minute/hour time-scales and year time-scales. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57 2.9 A toy model in which every white dwarf in the sample is considered to have a K band flux that varies by ∆% every t days. The density plot indicates the probability of detecting that combination of ∆ and t for white dwarfs in the UKIRT survey. The smallest grid scale is 1 day and 0.1%. . . . . . . . . . 57 2.10 KS band measurements for the seven white dwarfs in the sample, all offset by 0.2 mags. Observations were taken separated by hours and days. ∆K is the output KS band magnitude subtracted from the median magnitude of the lightcurve, such that > 0 is brighter than the median magnitude and < 0 is fainter. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63 List of figures xix 2.11 The lightcurves of the four stars in the field of view of WD 1226+110 that were used to perform differential photometry, each offset from one another. ∆K is the output KS band magnitude subtracted from the median magnitude of the lightcurve, such that > 0 is brighter than the median magnitude and < 0 is fainter. The median KS band magnitude of each comparison star is labelled. The lightcurve labelled ‘offset’ is the calculated offset from the four stars that was applied to the WD 1226+110 lightcurve to perform differential photometry. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65 2.12 K band lightcurve spanning ∼ 6 hours for WD 1226+110. ∆K is the output KS band magnitude subtracted from the median magnitude of the lightcurve, such that > 0 means WD 1226+110 is brighter than the median magnitude and < 0 means it is fainter. The errors are large due to non-photometric conditions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66 2.13 Lomb-Scargle periodogram for WD 1226+110 using the lightcurve in Figure 2.12. The vertical dashed line shows the 123.4 minute period observed by Manser et al. (2019) in the gaseous component of the disc for WD 1226+110. The periodogram reveals two peaks at 0.011 days and 0.039 days but have false alarm probabilities of 0.62 and 0.60 respectively and are unlikely to be real. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66 3.1 The Mg II 4481.125 Å spectral line from the MIKE/Magellan (blue) and HRS/SALT (orange) spectra spanning more than a decade for WD 0106−328. The spectra have been continuum normalised and are offset by 0.5. Spectra with SNR < 10 are not plotted. . . . . . . . . . . . . . . . . . . . . . . . . 81 3.2 The Mg II 4481.125 Å spectral line from the MIKE/Magellan (blue) and HRS/SALT (orange) spectra spanning more than a decade for WD 0408−041. The spectra have been continuum normalised and are offset by 0.5. Spectra with SNR < 10 are not plotted. . . . . . . . . . . . . . . . . . . . . . . . . 82 3.3 The Mg II 4481.125 Å spectral line from the MIKE/Magellan (blue) and HRS/SALT (orange) spectra spanning more than a decade for WD 1457−086. The spectra have been continuum normalised and are offset by 0.5. . . . . . 83 3.4 The Mg II 4481.125 Å spectral line from the MIKE/Magellan (blue) and HRS/SALT (orange) spectra spanning more than a decade for WD 1929+011. The spectra have been continuum normalised and are offset by 0.5. . . . . 84 xx List of figures 3.5 The Mg II 4481.125 Å spectral line from the MIKE/Magellan (blue) and HRS/SALT (orange) spectra spanning more than a decade for G29-38. The spectra have been continuum normalised and are offset by 0.15. Spectra with SNR < 10 are not plotted. . . . . . . . . . . . . . . . . . . . . . . . . . . 85 3.6 The Ca II 3933 Å spectral line from the MIKE/Magellan (blue) for G29-38. The spectra have been continuum normalised and are offset by 0.5. Spectra with SNR < 10 are not plotted. . . . . . . . . . . . . . . . . . . . . . . . . 86 3.7 The Mg II (top), Ca II (middle), and Si II (bottom) equivalent width (EW) measurements over time for WD 1929+011. Measurements from MIKE/Magellan are reported as blue data points between 5000–6000 days, and measurements from HRS/SALT are reported as orange data points after 7000 days. . . . . 87 3.8 Comparison of the equivalent width (EW) measurements from the MCMC analysis on the MIKE/Magellan (M) and HRS/SALT (S) stacked spectra for WD 1929+011. Markers represent different lines of Mg, Si, Fe and O with their transition wavelength along the x axis. . . . . . . . . . . . . . . . . . 88 3.9 The Mg II doublet (top), Ca II K line (middle), and Mg I line, 5183.602 Å (bottom) changes of equivalent width (EW) over time for G29-38. Mea- surements from MIKE/Magellan are reported as blue data points between 5000–7000 days, and measurements from HRS/SALT are reported as orange data points after 7000 days. . . . . . . . . . . . . . . . . . . . . . . . . . . 90 3.10 The stacked MIKE/Magellan and SALT/HRS data for the Mg II doublet for WD 0106−328. The darker blue line is the MIKE/Magellan spectra interpolated and convolved to the resolution of HRS/SALT. . . . . . . . . 92 3.11 The equivalent width of the Mg II doublet over time for WD 0106−328, WD 0408−041, and WD 1457−086. . . . . . . . . . . . . . . . . . . . . . 94 3.12 (a) The best fitting model to the 2020-08-24 HRS/SALT spectrum of the Mg- II doublet WD 1929+011. The best fitting Mg abundance is [Mg/H] = −4.25. Models are plotted every 0.2 dex spreading to +/- 1 dex in [Mg/H]. Models kindly provided by P. Dufour (Dufour et al., 2012). (b) The equivalent width of the Mg II doublet over time for WD 1929+011. The equivalent width for each model in (a) is plotted, showing the expected equivalent width values for each [Mg/H] abundance. . . . . . . . . . . . . . . . . . . . . . . . . . 97 3.13 The updated sinking timescales using 3D models including convection. The 1D sinking timescales reported in Table 3.2 are plotted for the three white dwarfs in the temperature regime of interest. Figure adapted from Cunning- ham et al. (2019) Figure 22. . . . . . . . . . . . . . . . . . . . . . . . . . . 99 List of figures xxi 3.14 The infrared fluxes of the white dwarfs with confirmed dust emission over time. S1 and S2 are the Spitzer channels 1 and 2 with data from Swan et al. (2020) (private communications). The K band data are from Rogers et al. (2020), as reported in Chapter 2. The horizontal dashed lines are the median values for each photometric band. . . . . . . . . . . . . . . . . . . . . . . 101 4.1 Model fits (red lines) to the Hydrogen Balmer line profiles from the X-shooter spectra for Gaia J0006+2858, Gaia J0611−6931, and Gaia J2100+2122. The best-fitting model parameters are labelled on each panel. Doubled peaked emission lines from the circumstellar gas discs can be seen in the wings of the hydrogen lines. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110 4.2 Model fits (red lines) to the X-shooter Helium lines for Gaia J0644−0352. The best-fitting white dwarf parameters and hydrogen abundance are labelled. The green lines show the points used to re-normalise the spectrum to the models. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 111 4.3 Model fits to the HIRES data using the average abundances for WD 1622+587. Dashed orange line shows the model for a magnetic field of 0kG, and the blue solid line shows the model for a magnetic field of 10kG. . . . . . . . 112 4.4 (a) [Mg/Ca] vs [Si/Ca] abundance ratios of the pollutant, normalised to solar, derived from the spectroscopically determined white dwarf parameters. The ellipses show a 1σ error ellipse on the abundance ratio. The Hypatia catalogue (Hinkel et al., 2014) of nearby main sequence stars is plotted as scatter points with a density histogram in the densest regions (around solar values - dashed lines at 0, 0). The red arrows show how sinking in the white dwarf atmosphere and heating affects the abundance ratios. Limits are shown with arrows. (b) [O/Ca] vs [Fe/Ca] abundance ratios normalised to solar, for reference the value of [O/Ca] normalised to solar for Bulk Earth is -0.7. An additional red arrow describes how the [Fe/Ca] ratio changes as the core-mass-fraction (CMF) of the planetary body is increased. . . . . . . . 115 4.5 The likely division of the observed oxygen abundance between metal oxides (MgO, SiO2, Al2O3, CaO, and FeO) and water-ice (yellow), for Gaia J0611−6931 and Gaia 0644−0352. The approach is described in Section 4.2.4, using abundances adjusted for sinking in steady-state. For Gaia J0611−6931 the fraction of excess oxygen is 0.48, this excess oxygen is likely carried in water-ice. For Gaia 0644−0352, the oxygen can all be accounted for in metal oxides, and the pollutant was most likely dry and rocky. . . . . . . . . . . . 121 xxii List of figures 4.6 The number ratio of elemental abundances in the parent body relative to Mg, normalised to solar (grey dash dot line). Abundances for the white dwarfs, Gaia J0611−6931 and Gaia J2100+2122, are shown as green diamond and black hexagon data points respectively, with 1σ errors, derived using the spectroscopic temperature and log(g). Abundances have been adjusted for sinking assuming steady state. The Fe/Mg for Gaia J0611−6931 lies 5.8σ above solar, and for Gaia J2100+2122 the Fe/Mg value lies 3.2σ above solar; both systems are enhanced in Fe compared to Mg. . . . . . . . . . . . . . . 123 4.7 The number ratio of elemental abundances in the pollutant material relative to Mg, normalised to solar (dashed line). Observed abundances for the white dwarf, Gaia J0644−0352, are shown as black points, with 1σ errors, derived using the spectroscopic temperature and log(g). Sinking effects are included in the model. The blue line represents the Bayesian model with the maximum likelihood (max Z). The black line represents the maximum likelihood median model with differentiation included and a 1σ error as the shaded dark grey region. The body that was accreted by Gaia J0644−0352 may be mantle rich, but there is insufficient statistical evidence to include an additional parameter in the model fit. The shaded light grey region indicates the range of abundances seen in nearby stars using the catalogue from Brewer & Fischer (2016). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125 4.8 The accretion rates of polluted white dwarfs with and without detectable gaseous discs in emission as a function of their effective temperatures. The data for the accretion rates of polluted white dwarfs without circumstellar gaseous discs are from Figure 3 in Xu et al. (2019), see reference for details. The accretion rates for gaseous disc systems are from Table 4.4 including the seven reported in this work. . . . . . . . . . . . . . . . . . . . . . . . 126 4.9 Single blackbody dust temperature for the dust emission from the literature versus the [Ca/Mg] ratio, normalised to solar, of the pollutant body. Those in this work are black squares, and those from the literature are grey. Horizon- tal dashed grey lines show apparent condensation temperatures for Al, Ca, Mg, Si and Fe (Lodders, 2003). The material polluting Gaia J0347+1624 is depleted in moderate volatiles and has a composition most like Ca-Al–rich inclusions, it is hypothesised that the material underwent incomplete conden- sation during formation explaining the extreme [Ca/Mg]. Therefore the dust surrounding Gaia J0347+1624 can survive to higher temperatures, and so the dust emission has a hotter overall inferred blackbody temperature. . . . . . 127 List of figures xxiii 4.10 (a) Comparison of the [Si/Ca] against [Mg/Ca] composition of the material that has polluted white dwarfs which have a detectable circumstellar gaseous disc compared to those without. (b) Comparison of the [Fe/Ca] against [O/Ca] composition. The grey ellipses show a 1σ , 2σ and 3σ error ellipse based on solar abundances with Ca, Mg and Si errors of 0.1 dex. If a white dwarf accreted material of solar abundance with typical errors of 0.1 dex, the contours show the distribution of expected abundances. The arrows show how the abundances are affected by: sinking in the white dwarfs atmosphere, heating of the pollutant and increasing the core mass fraction (CMF). . . . 132 7.1 Spectral energy distributions for the seven white dwarfs. The X-shooter spectra for the four white dwarfs observed are shown in grey, this is scaled to fit the spectroscopic models as the spectra suffer from flux loss due to non-ideal weather conditions and slit losses. There is a gap in the data between 6360 Å– 6375 Å and telluric absorption features are present in the reddest parts of the spectra. The photometric data points are given in black, errors are plotted but are smaller than the data points. The photometric and spectroscopic model fits are shown in orange (solid line) and blue (dashed line) where available. Missing GALEX FUV or NUV fluxes implies either a non-detection in that band, or the flux was flagged as it contained an artefact. 170 7.2 Model fits (red lines) to the HIRES data using the average abundances for Gaia J0006+2858. There is additional non-photospheric absorption blue- wards of the Ca II K line. . . . . . . . . . . . . . . . . . . . . . . . . . . . 176 7.3 Model fits (red lines) to the HIRES data using the average abundances for Gaia J0347+1624. Silicon is shown for reference, the absorption is not significant to 3σ and is therefore taken as an upper limit. . . . . . . . . . . 177 7.4 Model fits (red lines) to the HIRES data using the average abundances for Gaia J0510+2315. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 177 7.5 Model fits (red lines) to the MIKE data using the average abundances for Gaia J0611−6931. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 178 7.6 Model fits (red lines) to the HIRES data using the average abundances for Gaia J0644−0352. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 179 7.7 Model fits (red lines) to the HIRES data using the average abundances for Gaia J2100−2122. There is an additional non-photospheric absorption bluewards of the Ca II K line. . . . . . . . . . . . . . . . . . . . . . . . . . 180 List of tables 2.1 The UKIRT sample of dusty polluted white dwarfs and their properties. H/He column highlights the dominant element in the atmosphere. . . . . . . . . . 39 2.2 The resulting median magnitude and standard deviation found from the Gaussian fits to the photometry from the UKIRT observations in the J, H and K bands. The errors are quoted as the 16th and 84th percentiles from the posterior distributions. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44 2.3 The level of variability that can be detected at each magnitude in the K band, calculated from the field stars. The objects which were observed with 5 and 10 second frames were analysed separately. The median percentage variability is quoted for the field stars in bins with a width of one mag, along with the 16th and 84th percentiles. . . . . . . . . . . . . . . . . . . . . . . 49 2.4 Observations of the white dwarfs listing the number of sets of JHK photome- try obtained and the baseline for these observations. . . . . . . . . . . . . . 61 2.5 Table listing dates of observations for the white dwarfs in the sample, and number of exposures. N comp is the number of comparison stellar objects in the frame. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62 3.1 The sample of polluted white dwarfs stating their effective temperature and log(g) derived from a spectroscopic fit in previous works. The dust column highlights whether they are confirmed to have an infrared excess from circumstellar dust. . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73 3.2 The pollution levels and sinking timescales for the sample of white dwarfs. Sinking timescales interpolated from Koester et al. (2020). . . . . . . . . . 73 3.3 WD 0106−328 observation details. SNR determined from the continuum around the Mg II line at 4481.125 Å as defined in Section 3.1.2. . . . . . . 75 3.4 WD 0408−041 observation details. SNR determined from the continuum around the Mg II line at 4481.125 Å as defined in Section 3.1.2. . . . . . . 75 xxvi List of tables 3.5 WD 1457−086 observation details. SNR determined from the continuum around the Mg II line at 4481.125 Å as defined in Section 3.1.2. . . . . . . 76 3.6 WD 1929+011 observation details. SNR determined from the continuum around the Mg II line at 4481.125 Å as defined in Section 3.1.2. . . . . . . 76 3.7 G29-38 observation details. SNR determined from the continuum around the Mg II line at 4481.125 Å as defined in Section 3.1.2. . . . . . . . . . . . . 77 3.8 The median equivalent width (mÅ) and error for the HRS/SALT (S) data and the MIKE/Magellan (M) data obtained from the stacked spectra for the Mg II doublet at 4481.125 Å. . . . . . . . . . . . . . . . . . . . . . . . . . . . 87 3.9 The median equivalent width and error for the stacked HRS/SALT (S) data and the MIKE/Magellan (M) data for weaker lines for WD 0106−328, WD 1929+011, and G29-38. . . . . . . . . . . . . . . . . . . . . . . . . . 88 3.10 The inferred accretion rates assuming a bulk Earth composition for the best fitting [Mg/H] abundance for the 2020-08-24 data set (−4.25, in italics), and for abundance +0.2, +0.4, −0.2, and −0.4 dex from the measure [Mg/H] value. Macc, 10 yrs reports the mass of material that would accrete over a 10 year baseline. Assuming the abundance in the photosphere changed from [Mg/H] = −4.25, ∆Macc, 10 yrs reports the change in mass between the new [Mg/H] value and [Mg/H] = −4.25 that would accrue over 10 years, reported as a fraction of the mass of Ceres, a dwarf planet in the asteroid belt, and the fraction of the mass of Comet Halley. . . . . . . . . . . . . . . . . . . . . 96 3.11 Infrared variability from previous studies. The K band data is from Rogers et al. (2020) where the variability measure is an upper limit that corresponds to the measurement limit of the observations. The Spitzer data are from Swan et al. (2020) showing the percentage change and significance of the data.100 4.1 Stellar parameters derived from the spectroscopic (spec) and photometric (phot) fitting methods. Distances (D) are inferred from Gaia parallaxes. . . 104 4.2 Observations of the seven white dwarfs listing dates of observations, exposure times in seconds, and SNR as defined in Section 4.1.2. The SNR for X- shooter UVB, HIRESb and MIKE-blue were calculated from the continuum around the Ca II K line (3933.7Å). The SNR for X-shooter VIS, HIRESr and MIKE-red were calculated from the continuum around 6600Å. . . . . . . . 105 List of tables xxvii 4.3 Number abundances (log n(Z)/n(H(e)) of the pollutant material using the photometric and spectroscopic white dwarf parameters separately. For the derivation of the upper limits and the analysis of the abundances the spectro- scopic solutions are used. Photometric abundances for Gaia J0510+2315 are less reliable and are quoted in brackets. . . . . . . . . . . . . . . . . . . . 106 4.4 Total accretion rates based on Mg abundances for white dwarfs with an observable gaseous disc with emission features. This is calculated us- ing M˙ = (100/15.8)×MWD×10q×10[Mg/H(e)]×AMg/H(e)/τMg, where q= log10(MCVZ/MWD), AMg/H(e) is the atomic mass of Mg divided by the atomic mass of H or He, depending on the dominant atmospheric (atm) constituent, and τMg is the sinking time of Mg. . . . . . . . . . . . . . . . . . . . . . . 130 6.1 All the polluted white dwarfs with an infrared excess observed with the UKIRT and the dates they were observed. . . . . . . . . . . . . . . . . . . 160 6.2 All the polluted white dwarfs with an infrared excess observed with the UKIRT survey and the total number of frames. Each frame consists of a dithered stack of 5 exposures, each with the corresponding frame time of 5 or 10 seconds. FT = frame time for each individual exposure in the dithered stack, FN = total number of dithered stacks over all observations. . . . . . . 167 7.1 The absorption lines in Gaia J0006+2858. If two lines lie within ∼ 10 Å, equivalent widths and abundances were fitted simultaneously. . . . . . . . 171 7.2 The absorption lines measured in the HIRES spectrum of Gaia J0347+1624. There are also Na D lines at both 5890 Angstrom and 5896 Angstrom. However, in checking the ISM model, this line of sight intercepts LIC (RV of 22.36 km/s) and Hyades (RV of 14.16 km/s). This is close to the measured RV of 15 km/s so it is deduced that these lines are not photospheric. . . . . 171 7.3 The absorption lines in the HIRES spectra of Gaia J0510+2315. If two lines lie within ∼ 10 Å, equivalent widths and abundances were fitted simultane- ously. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 172 7.4 The absorption lines in Gaia J0611−6931. If two lines lie within ∼ 10 Å, equivalent widths and abundances were fitted simultaneously. . . . . . . . 173 7.5 The absorption lines in Gaia J0644−0352. If two lines lie within ∼ 10 Å, equivalent widths and abundances were fitted simultaneously. . . . . . . . 174 7.6 The absorption lines measured in the WD 1622+587 HIRES spectra. If two lines lie within ∼ 10 Å, equivalent widths and abundances were fitted simultaneously. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 175 xxviii List of tables 7.7 The absorption lines measured in Gaia J2100+2122. If two lines lie within ∼ 10 Å, equivalent widths and abundances were fitted simultaneously. . . . . 175 7.8 The upper limit equivalent widths in mÅ of the material polluting the seven white dwarfs in this study. ∗ denotes when gaseous emission is present at this wavelength. When calculating abundance upper limit, add 12 percent onto the equivalent width quoted to be conservative. . . . . . . . . . . . . 176 7.9 The resulting median parameters from the Bayesian models, errors are given as 16th and 84th percentiles. A description of the parameters is presented in Section 4.2.4. The basic model uses a best-fitting initial stellar composition, ‘Stellar metallicity’; this is omitted from this table as most white dwarfs are equally likely to have accreted material that started with a wide range of initial compositions. The pollution fraction and the most likely accretion phase: build up (BU), steady state (SS) or declining phase (DP) are output for the basic primitive model. Italics represent optional complex parameters that are only invoked when the Bayes factor is > 1. The fragment core fraction (fcf) is quoted based on number fraction, and in brackets mass fraction, which is calculated assuming the output core and mantle composition of the median Bayesian model. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 182 Chapter 1 Introduction Astronomy has entered the era of large observational surveys, providing us with unprece- dented information about the Universe, and our place within in. Exoplanets are now known to be ubiquitous with planets predicted to outnumber stars in the Milky Way (Cassan et al., 2012). Over the last few decades the field of exoplanetary science has rapidly developed with the focus shifting away from discovery and towards the characterisation of exoplanets and exoplanetary populations. Key questions driving the field are: what are the compositions of exoplanets? How similar are exoplanet compositions and system architectures to that of the solar system? Do other planetary systems have the conditions to develop extrater- restrial life? White dwarfs are unique laboratories that can contribute to studying these questions. Between 25–50% of white dwarfs have accreted planetary material from their remnant planetary systems (Zuckerman et al., 2003, 2010; Koester et al., 2014; Wilson et al., 2019). Spectroscopic observations of these white dwarfs reveal the bulk composition of the exoplanetary material accreted. This is the only observational method to directly infer exoplanetary interior composition. These polluted white dwarf systems and what they reveal about the composition of exoplanetary material is the focus of this thesis. This chapter will first introduce the knowledge gained from studying planetary composition and architecture of the solar system, then will discuss this in relation to exoplanets, and finally will combine these topics to introduce white dwarfs and their remnant planetary systems. 1.1 Solar System Before the 1990s, the solar system was the only planetary system known. Therefore, research on planets concentrated on the solar system planets and small bodies consisting of the four terrestrial bodies: Mercury, Venus, Earth and Mars, the four gas giants: Jupiter, Saturn, Uranus and Neptune; and the small bodies: asteroid belt, Kuiper-belt, Oort cloud and their 2 Introduction various constituents. The solar system formed from a collapsing cloud of dust and gas around 4.567 billion years ago forming the proto-Sun with an orbiting disc of gas and dust (Connelly et al., 2012). The first solids condensed out of the protoplanetary nebula which constituted the building blocks for the planets. The process of planetesimal formation is an extensive field of scientific research, for further details refer to the review by Raymond & Morbidelli (2022). Although the solar system is the only planetary system where planetary material can be directly sampled, numerous open questions regarding the formation and evolution of the solar system remain. 1.1.1 Rocky bodies in the solar system Small bodies Meteorites are debris from outer space that survive the journey through the Earth’s at- mosphere to land on the surface. Laboratory measurements reveal their composition and structure and thus gives insight into solar system formation, composition, and evolution. The main classifications of meteorites are chondritic and achondritic. Chondritic meteorites are primitive asteroids that contain the earliest solids and consequently define the age of the solar system as 4.567 billion years (Connelly et al., 2012). They provide the abundances of the material that condensed to form the solar system constituents, which provides powerful information about the formation of planetary bodies. Achondritic meteorites have undergone geological and collisional processing (melting and core–mantle differentiation, which are discussed in more detail in the following section). These meteorites include those from the asteroid belt, Mars, and the Moon, and can consequently reveal the formation and internal structure of the terrestrial planets (Mittlefehldt et al., 1998; Weisberg et al., 2006). Chondritic meteorites are composed of: chondrules, refractory inclusions (0.01—10 % volume), metallic Fe and Ni (<0.1—70% volume), and are surrounded by a fine-grain matrix (1—80% volume). Chondrules are millimeter sized particles that form by the melting of material in the solar nebula, which are subsequently deposited onto the mid-plane of the protoplanetary disc with other particles. Chondrules are mostly composed of olivine and pyroxene minerals which have crystallised rapidly between temperatures of ∼ 1800 and 1300 K (Scott, 2007). Refractory inclusions are crystalline silicates and oxides which form at high-temperatures (> 1300 K) from processes such as melting, evaporation and condensation. The inclusions consist of calcium-aluminium rich inclusions (CAI’s) and amoeboid olivine aggregates (AOA). Fe and Ni metallic grains can be present both within and outside of chondrules. Similar to chondrules and inclusions they also form at the highest temperatures. 1.1 Solar System 3 The fine-grain matrix material contains volatile-rich minerals, which fills the space between the other constituents of the chondritic meteorite (Alexander et al., 2018). The three major classes of chondrites are: carbonaceous, ordinary, and enstatite. These classes are sub-divided into 15 groups with variations that are categorised by chondrule sizes, the quantity and composition of the refractory inclusions, the amount of metallic iron and nickel, and the relative fraction of mineral matrix (Bizzarro et al., 2017). Carbonaceous chondrites contain a significant fraction of volatiles, including carbon, and are thought to originate furthest from the Sun, they account to 5% of known meteorites. The CI group of carbonaceous chondrites have a composition most similar to the solar photosphere, as seen in Figure 1.1, and chemically are the most primitive group of meteorites. Ordinary chondrites are most common, consisting predominantly of chondrules, with small components of refractory inclusions, Fe, Ni metals, and matrix. The rarest are enstatite which account to 2% of known meteorites, which contain a significant fraction of the mineral enstatite with Fe taking the form of metallic-Fe rather than oxides. Enstatites likely formed in reduced environments close to the proto-Sun. Meteorites are thought to be the result of the fragmentation of asteroids that fall to Earth. Understanding how the classes of meteorites link to classes of asteroids has been subject to intensive study as their connection and properties can inform about the history and evolution of the solar system. Asteroids have experienced significant dynamical and collisional evolution which affects their composition, size, shape, rotation and location; making it difficult to model the evolution of the small bodies. Asteroids are defined based on their albedo and spectral features. C-class (carbonaceous) asteroids are thought to be the precursors to carbonaceous chondrite meteorites and are typically associated with the outer asteroid belt, they have low albedos and can contain up to 10% water by mass. Whereas, the S-type (stony) asteroids have much higher albedos, have spectral signatures of silicates and are dry with less than 0.1% water by mass, these asteroids are believed to be associated with the ordinary chondrite class of meteorites. (Michel et al., 2015; Alexander et al., 2018). The asteroid belt contains a mixture of both volatile-poor asteroids (S-type) in the inner regions and volatile-rich asteroids (C-type) in the outer regions. Originally it was hypothesised that this was due to the ice lines and formation conditions, but abundance and isotopic studies revealed these asteroids formed in different reservoirs (Warren, 2011; Budde et al., 2016; Kleine et al., 2020). Significant progress has been made modelling the solar system formation and evolution, but no proposed model is able to explain all observations in the solar system, including how the asteroid belt became a mix of two distinct reservoirs of small bodies (e.g. Raymond & Morbidelli, 2022). 4 Introduction Al Ti Ca Ni Fe Cr Si Na O −1.0 −0.5 0.0 0.5 1.0 [X /M g] - [X /M g] ¯ Solar Bulk Earth CI Chondrites Earth Mantle Earth Core* Fig. 1.1 Abundance ratio of element X to Mg, normalised to solar for compositions of: CI Chondrites, bulk Earth, Earth’s mantle, and Earth’s core. Data are from: McDonough (2003b); Grevesse et al. (2007). ∗ The core composition is: 2/3 core and 1/3 mantle. 1.1 Solar System 5 Core–mantle differentiation The Earth has an iron rich core and is surrounded by an iron poor mantle. These layers are key to the Earth being habitable and maintaining habitable conditions over time. The molten Earth’s core generates a magnetic field which is crucial for shielding the Earth from harmful solar winds. The mantle and crust result in tectonic activity which regulates the temperature and atmospheric conditions of the planet, for example, by regulating the quantities of carbon dioxide in the atmosphere via the carbon cycle, thus keeping the temperature stable. Core–mantle differentiation is the process by which melting causes siderophilic (‘iron loving’) elements to separate and sink towards the centre of the body, and lithophilic (‘rock loving’) elements to go towards the surface of the body where they can combine readily with oxygen. In order to melt the material such that core–mantle separation can occur, a heat source is required. For larger bodies greater than ∼ 1500 km, the release of gravitational potential energy during formation can result in large scale melting and differentiation into a core and mantle. However, in the solar system the decay of 26Al is required to form cores in asteroids as they are not massive enough to generate heat via gravitational potential energy. Modelling indicates that for planetesimals greater than ∼ 20 km in diameter, 26Al heating could cause melting and core formation (Hevey & Sanders, 2006; Elkins-Tanton, 2012). Terrestrial planets The terrestrial planets in the inner solar system are all composed of silicate material and metals (e.g. iron). All four terrestrial planets are differentiated and have a core component and a mantle/crust component. Differentiation into core and mantle is expected to have occurred early in the formation of terrestrial bodies as the gravitational potential energy released upon accretion was significant enough to melt their interior. The mass fraction of the metallic core in Mercury, Venus, Earth, and Mars decreases as a function of distance from the Sun from approximately 3/4, to 1/3 (Venus and Earth), to 1/5, respectively (Sohl & Schubert, 2007). The structure and composition of the Earth’s interior has been subject to extensive studies. For example, measuring the velocities of seismic waves propagating through the Earth reveals the density profile and therefore layered structure of the Earth; and studying rock samples from the crust and upper mantle reveals the composition of the outer layers of the Earth. The best estimate of the bulk constituents of the Earth by mass is: 32% iron, 29.7% oxygen, 16.1% silicon, 15.4% magnesium, 1.82% nickel, 1.71% calcium, 1.59% aluminium, plus others in smaller amounts (McDonough, 2003b). The core of the Earth is enhanced in siderophilic elements, such as iron and nickel, and the mantle of the Earth is depleted in siderophilic 6 Introduction Fig. 1.2 The relative abundance of X/Mg for bulk Earth (open circles) and classes of carbona- ceous chondrites (CM class is shown as filled circles, CV as filled squares, and CO as open triangles), normalised to CI chondrites, plotted against the log temperature at which 50% of the element would have condensed from a gas at 10−4 bar. The shaded regions defines the area covered by the three carbonaceous chondrites classes. The volatility depletion trend in bulk Earth and carbonaceous chondrites is clear. Figure from McDonough (2003a). 1.2 Exoplanets 7 elements, as can be seen in Figure 1.1. The bulk composition of the Earth, disregarding volatiles, matches that of the most primitive meteorites, CI chondrites, as seen in Figure 1.1. These meteorites are assumed to be the building blocks of the Earth, although exactly which is debated. The volatiles, especially water, significantly contributed to the Earth becoming a habitable planet. Figure 1.2 shows the volatility depletion trend in bulk Earth. The chemical composi- tion of the Earth was largely inherited from the condensation and accretion of nearby material in the protoplanetary disc. The standard model of the disc which the solar system formed has the snow line at 2.7 au (Hayashi, 1981), the Earth formed within this line so appears to have formed dry. However, water is present on the Earth with a mass fraction of ∼ 0.1%, greater than would be expected based on its formation location (Van Dishoeck et al., 2014). It is actively debated as to whether the volatiles were already present on Earth or whether volatile-rich bodies were later accreted. If they were later accreted, it is thought that volatile delivery from icy planetesimals is crucial for the present day levels of water on Earth (Wänke, 1981; McDonough, 2003a; Ballhaus et al., 2013; Yoshizaki & McDonough, 2021). The ratio of deuterium, a stable isotope of hydrogen, to hydrogen ([D/H]) is a useful diagnostic of where in the solar system objects formed as the [D/H] ratio changes depending on the location in the protoplanetary disc the bodies originated from (Hallis, 2017). Therefore, the [D/H] ratio can be used to investigate what bodies contributed to seeding the Earth with water and where in the solar system they originated from. The [D/H] ratio of the outer asteroid belt and the Kuiper belt have similar values to the Earth’s oceans and as they also have high fractions of water, it is hypothesised that these reservoirs significantly contributed to water delivery on Earth and the terrestrial planets (Van Dishoeck et al., 2014). Instabilities which cause mixing of the different reservoirs are pivotal for understanding of the solar system evolution and volatile delivery, for further discussions see (Raymond & Morbidelli, 2022). 1.2 Exoplanets 1.2.1 Discovery It wasn’t until the 1990s that the first exoplanets were discovered (e.g. Wolszczan & Frail, 1992). The confirmation of exoplanets gauged much interest and excitement worldwide as people began to wonder what these other worlds are like, could they potentially be habitable and host life, and ultimately whether we are alone in the universe. The advancement of modern technology has expanded our knowledge immensely; to date, over 5,000 exoplanets have been discovered, Figure 1.3 shows the exoplanets discovered 8 Introduction 100 102 104 106 108 Period (days) 10 1 101 103 M as s ( Ea rth M as se s) Transit Radial Velocity Imaging Microlensing Pulsar Timing Pulsar Timing Fig. 1.3 Orbital period versus planet mass for exoplanets with both measured/estimated orbital period and mass. The colours highlight the different discovery methods. Figure using data from the NASA exoplanet archive1. 1.2 Exoplanets 9 thus far1. Exoplanets appear to be ubiquitous and outnumber stars in the Milky Way (Cassan et al., 2012). The first exoplanets discovered were around a millisecond pulsar, found using precise timing measurements of the pulses (Wolszczan & Frail, 1992). Following this, the first exoplanet discovered around a Sun-like star was 51 Pegasi b (Mayor & Queloz, 1995), a Jupiter mass planet discovered with the radial velocity technique. The radial velocity technique detects the presence of an orbiting companion from the Doppler shift in the stellar spectral lines. The star and planet orbit the common centre of mass; if the stellar mass is known, the minimum mass of the planet can be inferred (Fischer et al., 2015). Despite the successes of these two detection methods delivering the first exoplanet discoveries, more than three quarters of exoplanets were discovered with the transit technique. The transit technique photometrically monitors a star; as a planet passes in front of the stellar disc a flux drop proportional to the relative sizes of the star and the planet is observed. This method provides information on the radius of the planet, assuming the radius of the star is known, and constrains the inclination of the orbit. The first exoplanet transit detection was around HD 209458, where radial velocity measurements had previously confirmed an orbiting planet (Charbonneau et al., 1999). Other techniques have been employed to detect exoplanets, such as, direct imaging, microlensing, and astrometry. The current detection statistics are low, but with technology developments, current (e.g. Gaia, Collaboration et al., 2016) and future space missions (e.g. LIFE, Quanz et al., 2021), as well as the upcoming ground based extremely large telescopes, many more planets are expected to be discovered with these methods. Figure 1.3 shows the exoplanets discovered so far with the regions in parameter space that these techniques are most sensitive to clearly standing out. The success of the Kepler and TESS space missions significantly contributed to the discovery of transiting exoplanets (e.g. Borucki et al., 2010; Guerrero et al., 2021), which has allowed the field of exoplanet demographics to flourish (e.g. Batalha, 2014). The transit technique is sensitive to larger planets, on close in orbits as this increases the transit depth and the probability of the transit (as the star and planet must be aligned). The radial velocity technique is sensitive to more massive planets on shorter orbits, as these produce larger amplitude signals and multiple orbital periods can be measured. With current and future generations of higher-resolution and higher-stability instruments, radial velocity surveys can push to discover Earth mass planets around Sun like stars. However, stellar activity signals now dominate the noise so in order to discover planets that are increasingly lower in mass and on longer orbital periods, stellar variability needs to be better understood and quantified (e.g. Rajpaul et al., 2015; Vanderburg et al., 2016; Nava et al., 2019). Direct imaging is most sensitive to planets on 1Data from: https://exoplanets.nasa.gov/discovery/discoveries-dashboard/ on 26/07/2022 10 Introduction wide (> 5 au) orbits, the best targets tend to be those that are young and massive as these are self-luminous due to heat from formation. The photons from the planet are directly measured demonstrating the power of this technique as these planets can be studied both photometrically and spectroscopically. Exoplanets detected via microlensing events happen when a background object becomes aligned with an intermediate star with an orbiting planet. The star and planet act as a lens and cause brightening of the background source. This technique does not depend as strongly on planetary mass so has a wide discovery space which includes free-floating planets (Gaudi, 2012). 1.2.2 Characterisation With the number of exoplanet discoveries increasing everyday, the research focus has shifted towards their characterisation. The field of characterisation is focused on the composition of their interior and atmosphere, and the potential of these worlds to be habitable. To constrain a planet’s bulk composition, measurements of its mass and radius are required. The first exoplanet characterised in this way was HD 209458, radial velocity measurements provided a minimum planetary mass (M sin i), and transit photometry provided the inclination (i), resulting in the determination of the true planetary mass. The combined mass estimate from the radial velocity technique with the radius estimate from the transit technique provides an estimate for the average density of the exoplanet, the average density was low and due to its size it was designated a gas giant (Charbonneau et al., 1999), therefore beginning the era of exoplanet characterisation. Theoretical mass-radius relationships derived from interior models can be compared to derived mass and radii of exoplanets to infer the interior structure of the exoplanets. Various compositions and interior structures can be invoked in the models, for example, the models can be in the form of a single composition, e.g. 100% Fe, or 100% H/He; or multi-layered composition with combinations of: Fe core, MgSiO3 mantle, H2O layer; or more exotic compositions (e.g. Seager et al., 2007; Dorn et al., 2015). However, degeneracies arise as different compositions can produce similar mass-radius curves, thus introducing uncertainties in the determination of bulk compositions. In order to more accurately determine the bulk composition of these bodies, additional input information or alternate measurement techniques are required. The study of exoplanet atmospheres reveals information about the outer layers of an exoplanet (e.g. Crossfield, 2015; Madhusudhan, 2019). Over the coming decades the study of exoplanetary atmospheres will be a key focus of the field. For transiting planets, the opacity of the atmosphere varies with wavelength due to its composition, and measuring the transit depth as a function of wavelength, the ‘transmission spectrum’, can reveal information 1.3 White dwarfs 11 about its atmospheric composition. Transiting planets are also eclipsed by their host star, subtracting the spectrum of the star alone from the combined star and planet spectrum, reveals the thermal ‘emission spectrum’ of the planet. Current space missions such as Hubble and JWST have been revealing the spectrum of increasingly smaller planets. Directly imaged planets allow direct spectroscopy without contamination from the stellar spectrum to reveal their thermal emission spectra. With plans for future large telescopes e.g. LUVOIR, LIFE missions, this technique will enable characterisation of non-transiting planets, not possible via other methods and on longer period orbits. With increasing amounts of data and computing power, the field of exoplanet characteri- sation is moving towards a combined approach, where models and results from sub-fields are combined to gain a global understanding of individual exoplanets, for example, interior and atmospheric coupling (e.g. Madhusudhan et al., 2020). The study of exoplanet atmospheres reveals the bulk properties of exoplanets (e.g. planetary mass, radius, atmospheric proper- ties), these can be input into exoplanet interior models to constrain the internal structure and thermodynamic conditions. Continued efforts to combine observations and computational techniques from sub-fields will help to piece together what these other worlds are like. 1.3 White dwarfs 1.3.1 Formation White dwarfs are the final product of stellar evolution for stars with an initial main sequence mass in the range 0.07 to 8 M⊙. It is predicted that as the Milky Way ages more than 97% of the stars will eventually evolve to white dwarfs (Fontaine et al., 2001). Figure 1.4 shows how a solar mass star evolves from the zero age main sequence (ZAMS) to the white dwarf phase. On the main sequence, stars burn hydrogen to helium via the proton-proton chain, building up helium in their core. Stars spend ∼ 90% of their lifetimes on the main sequence. Once hydrogen fusion in the core ceases, the central temperatures are not sufficient to ignite helium fusion, however, the hydrogen shell surrounding the core reaches the temperature and pressure required for proton-proton chain fusion, and hydrogen shell burning commences. The stellar envelope expands, cooling the star and increasing the luminosity, positioning the star on the sub giant branch (SGB). As the shell burns hydrogen to helium, helium is dumped onto the core, increasing its mass and eventually the helium core becomes degenerate. The stellar envelope continues expanding causing a rapid increase in luminosity and it travels up the red giant branch (RGB). Due to the larger temperature gradients, energy transport is via convection rather than radiation. This causes large convection zones which dredges 12 Introduction Fig. 1.4 A diagram highlighting the evolution of a solar mass star from the zero age main sequence (ZAMS) to the white dwarf phase. Figure from Carroll & Ostlie (2017). 1.3 White dwarfs 13 up material that has been modified by nuclear processes - the first convective dredge up. Helium continues to be produced and dumped onto the core, and as the core is degenerate, the increased mass causes core contraction. Eventually the conditions in the core are sufficient for helium to fuse via the triple alpha process; this causes the helium flash. The core expands and is no longer degenerate, so helium core burning and hydrogen shell burning occurs as the star settles onto the horizontal branch (HB). The HB is the equivalent of the hydrogen burning main sequence for helium burning, except the lifetime is much shorter due to the higher luminosity. Once the star has exhausted helium in the core, the core contraction results in a degenerate core consisting of C and O. The temperature is not sufficient for C and O to burn in the core, but the core contraction generates sufficient energy that helium and hydrogen shell burning can occur. The outer layers of the star expands and luminosity increases moving it up the asymptotic giant branch (AGB) to transition to a red super-giant. Instabilities in the burning shells result in thermal runaways driving thermal pulses which cause the outer layers of the star to be ejected producing a planetary nebula. The remnant stellar core has insufficient gravity to fuse so becomes a carbon-oxygen white dwarf supported entirely by electron degeneracy pressure. The average mass of a white dwarf is ∼ 0.6 M⊙ with radii similar to that of the Earth. They form extremely hot (>100,000 K) and radiate as they fade over long cooling times, τcool: τcool ∝ L−5/7, (1.1) where L is their luminosity (Mestel, 1952; Van Horn, 1971). This demonstrates the connection between a white dwarfs age and its temperature, with the coolest systems being the oldest. Consequently, white dwarfs can be used to constrain the ages of stellar populations in the Milky Way (e.g. Fontaine et al., 2001). It is expected that white dwarfs will eventually cool and fade into a black dwarf (Adams & Laughlin, 1997). 1.3.2 White dwarf properties White dwarfs are powerful laboratories to study numerous aspects of physics and astrophysics, including: degenerate states of matter, age history of stellar populations in the Milky Way, and the bulk composition of exoplanets. However, without knowledge of the stellar properties, this work is not possible. The following discussion focuses on the determination of the effective temperature and log(g). There are two main methods to observationally determine these parameters: the spectroscopic method and the photometric method. The spectroscopic method fits white dwarf models to the H/He lines in their spectrum. The cores of white dwarfs are made of degenerate matter due to the high pressures (∼ 106 g cm−3), this electron degeneracy pressure prevents the further collapse of the core. These degenerate 14 Introduction C-O cores are surrounded by a layer of non-degenerate gas, most often He, with a thin H layer on top. Therefore, the spectra of most white dwarfs exhibit only absorption lines from hydrogen and/or helium, and these absorption lines are pressure broadened by a variety of mechanisms. For example, the Stark Broadening Effect; due to the extreme temperatures, the envelopes of the white dwarfs are often ionised, these moving charges give rise to local electric fields. As a consequence of the extreme densities, the neutral emitters feel the presence of the electric field which results in the shifting of energy levels via the Stark Broadening Effect (Stark, 1913), one of the main causes of the broadening of the spectral lines in the temperature range of interest in this thesis (> 15,000 K). The profiles of the Balmer lines, and optical helium lines are very sensitive to temperature, therefore, comparing the observed H/He absorption profiles with those predicted from theoretical models, the white dwarf parameters can be determined, this is the spectroscopic method (e.g. Holberg et al., 1985; Bergeron et al., 1992, 2011). The photometric method fits white dwarf models to broad band photometry. For this method, the white dwarf models are convolved with the band-passes of the photometric filters and are fitted to the broad band photometry. Combining the parallax to give the distance, theoretical mass-radius relations, and the synthetic photometry gives the best-fitting effective temperature and surface gravity, (log(g) = log(M)−2log(R)+4.437 in solar units). Reddening becomes important for objects > 100 pc, and so it is crucial to de-redden the photometry before fitting for accurate parameter determination. Genest-Beaulieu & Bergeron (2019) highlights the importance of including the Sloan Digital Sky Survey (SDSS) ‘u’ band for accurate and consistent parameter determination as often the optical photometry lies on the Rayleigh Jeans slope, but the ‘u’ band samples the peak of the white dwarf’s spectral energy distribution (SED), therefore, including ‘u’ band photometry allows for improved estimates of the effective temperature. 1.3.3 Spectral classification White dwarfs tend to be classified based on their spectra. Spectral classifications begin with the letter ‘D’ which highlights that it is a degenerate object. The main spectral classifications of white dwarfs are: DA, DB, DC, DO, DQ, and DZ, Figure 1.5 shows an example spectrum from each sub-class. Most white dwarfs are DAs, where strong hydrogen absorption lines are present in their atmospheric spectra due to an outer hydrogen envelope on top of a helium layer. DBs, however, show strong helium absorption lines (McCook & Sion, 1999), it is hypothesised that they fused their outer hydrogen envelope during the post main sequence and are left with a pure helium envelope. The pressure broadening of the spectral lines (discussed above), makes the hydrogen and helium lines strong and observable. The strength of the DA 1.4 Polluted White Dwarfs 15 and DB spectral features change by significant amounts between stars due to the white dwarf parameters and differing histories, for example, due to accretion events, convective mixing or rotation (Fontaine & Michaud, 1979; Chayer et al., 1995a,b). DC white dwarfs contain no hydrogen or helium lines, this most often occurs when DA white dwarfs cool below 6,000 K, or when DB white dwarfs cool below 11,000 K. Therefore, below 6,000 K it is not possible to differentiate between hydrogen or helium dominated white dwarf atmospheres. DO white dwarfs contain absorption lines of ionised helium (He II), it is thought that around 45,000 K these DO white dwarfs transition to DB white dwarfs. DQ white dwarfs have spectra that show strong carbon absorption features. Cool DQs are thought to have deep convection zones which can result in carbon dredge up from the core to the atmosphere. Hot DQs are thought to arise due to binary mergers. DZs are white dwarfs which contain absorption lines from metals. They are cool DC white dwarfs which have accreted metals into their atmosphere. It is possible for a white dwarf to have multiple classifications and for those the letters are listed in order of which absorption features dominate, for example, DAZ and DBZ are DA and DB white dwarfs which also show strong metal absorption features. The combination of strong gravitational and electric fields in the outer stellar layers causes elements to become separated by diffusion. This leaves only the lightest elements, hydrogen or helium, in the observable atmosphere. The separation of the elements is expected to occur quickly as diffusion time scales are much shorter than white dwarf cooling times (∼ 1 Gyr for an effective temperature of 10,000 K). This picture is complicated by stellar winds, convection, accretion and radiative levitation (Vauclair et al., 1979; Koester & Chanmugam, 1990) which can disrupt the layers. 1.4 Polluted White Dwarfs In recent years evidence for planetary systems around white dwarf stars has grown, from the discovery of planets orbiting white dwarfs (e.g. Vanderburg et al., 2020; Blackman et al., 2021), to approximately 25–50% of white dwarfs being observed with planetary material polluting their atmospheres (e.g. Zuckerman et al., 2003, 2010; Koester et al., 2014; Hollands et al., 2017). The latter are titled polluted white dwarfs. The rapid gravitational settling times in comparison to the white dwarf’s cooling age implies that for metals to remain in the white dwarfs’ atmosphere ongoing accretion from an external reservoir is required (Koester, 2009). Figure 1.6 shows a schematic of these polluted white dwarf systems including the settling of the metals. vMa2 was the first single white dwarf discovered, but spectroscopic follow up revealed it had strong metal lines and was therefore designated as a ‘F type’ star (Van Maanen, 16 Introduction Fig. 1.5 Example spectra of white dwarf spectral types from LAMOST. Figure from Gen- tile Fusillo et al. (2015). 1.4 Polluted White Dwarfs 17 Fig. 1.6 A diagram to demonstrate the polluted white dwarf systems, taken from Jura & Young (2014). The composition of the pollution in the atmosphere will be discussed in Section 1.4.2 and the planetary debris will be discussed in Section 1.4.3. 1917). Almost a century later it was determined that this was a white dwarf which had been polluted with planetary material. Therefore, exoplanetary material has been observed for over a century. The first unambiguous detection of metal pollution that was attributed to a white dwarf was G74−7, this was hypothesised to be accretion from the interstellar medium (ISM) (Lacombe et al., 1983; Dupuis et al., 1993). This hypothesis involves a white dwarf transitioning a dense region of the ISM; dust and gas become gravitationally attracted to the star and therefore polluting its atmosphere. Much evidence invalidated this claim with the main argument being that the abundances of the accreted material did not correlate with the position and kinematics of the white dwarf in comparison to the ISM (Hansen & Liebert, 2003; Farihi et al., 2010a). It is now unambiguously known that these polluted white dwarfs are accreting material from their remnant planetary systems. 1.4.1 How do the white dwarfs become polluted? The favoured theory regarding the pollution source is that it originates from planetesimals scattered on eccentric orbits towards the white dwarf where they tidally disrupt, migrate inwards, sublimate, and finally accrete onto the atmosphere of the white dwarf (as seen in Figure 1.6). These planetesimals have masses comparable to solar system asteroids (Debes 18 Introduction & Sigurdsson, 2002; Jura, 2003; Farihi et al., 2010a; Jura & Young, 2014; Veras et al., 2014). Additional support for the asteroid tidal disruption model is the discovery of disintegrating planetesimals transiting the polluted and dusty white dwarf WD 1145+017 (Vanderburg et al., 2015). The deepest transit of this star has a period of 4.5 hours and blocks 60% of the flux from the star in the optical (Rappaport et al., 2016; Gary et al., 2017). Following this discovery, transiting debris was found around ZTF J032833.52−121945.27 with periods of 9.9 and 11.2 hrs (Vanderbosch et al., 2021). For the asteroid tidal disruption model, planetary bodies need to be scattered inwards from the surviving outer planetary system of the white dwarf. There are many dynamical mechanisms which can perturb planetesimals onto star-grazing orbits, the favoured theories cite planets as the source. Stellar mass loss can induce instabilities whereby a planet or planets perturb planetesimals onto star-grazing orbits (Debes & Sigurdsson, 2002; Bonsor et al., 2011; Veras et al., 2011; Debes et al., 2012a; Mustill et al., 2018). There is now evidence that planets can survive to the white dwarf phase (Gänsicke et al., 2019; Vanderburg et al., 2020), supporting these models. Alternative mechanisms to scatter planetesimals include: perturbations due to wide binary companions (Hamers & Portegies Zwart, 2016; Veras et al., 2016; Petrovich & Muñoz, 2017; Stephan et al., 2017; Smallwood et al., 2018), and the liberation of exo-moons (Payne et al., 2016; Trierweiler et al., 2022). 1.4.2 Composition of polluting bodies When assessing the composition of pollutant bodies, it must be acknowledged that certain phenomena occurring inside white dwarfs can contribute to metals polluting their atmo- spheres. For white dwarfs with an effective temperature greater than ∼ 25,000 K, radiative levitation can lead to an observable amount of heavy elements in the atmosphere (Vauclair et al., 1979; Chayer et al., 1995a,b). Convection zones that extend deep into the white dwarf interior can also cause the dredge up of heavy elements such that they appear in the outer layers of the convective zone of white dwarfs (Pelletier et al., 1986). This causes difficulty in determining the origin of polluting heavy elements for these temperature ranges. However, for white dwarfs with effective temperatures between 10,000 K and 25,000 K, the pollutant must be external, and therefore, the bulk composition of the pollutant can be directly measured using spectroscopic techniques. Spectroscopic studies of polluted white dwarfs combined with white dwarf atmospheric models allow the measurement of the bulk composition of extrasolar planetesimals making them unique laboratories. So far, 23 heavy elements have been discovered across all polluted white dwarfs (see references in Table 1. of Klein et al., 2021). The polluted white dwarf GD 362 has absorption features from 19 elements, the most elements detected for a given 1.4 Polluted White Dwarfs 19 white dwarf (e.g. Zuckerman et al., 2007; Xu et al., 2013; Melis & Dufour, 2016; Xu et al., 2017). For the majority of systems, the polluting materials’ abundance resembles rocky material similar to asteroids and the terrestrial planets in the inner solar system as seen in Figure 1.7 (Jura & Young, 2014). The inferred accreted mass is consistent with the masses of solar system asteroids (∼ 1020 g). The polluted white dwarf SDSS J073842.56+183509.06 has the most mass inferred in its photosphere, 7×1023 g, similar to the mass of the dwarf planet Ceres in the solar system (Dufour et al., 2012). Accretion phase: In order to reconstruct the abundances of the parent body that was accreted by the white dwarf, it is crucial to consider the differential sinking of elements and the affects this can have on the observed abundances. There are three accretion and diffusion scenarios that could occur. Build up stage - when accretion has just started so the observed abundances correspond closely to the actual abundances of the parent body, nX(A) par nX(B) par = nX(A)WD nX(B)WD , (1.2) where nX(A) par and nX(B) par are the abundances for element A and B respectively of the parent body before being accreted by the white dwarf, and nX(A)WD and nX(B)WD are the derived number abundances for element A and B from the observations of the photosphere of the white dwarf (Koester, 2009; Harrison et al., 2018). Steady state – when accretion and diffusion processes have reached an equilibrium. For this phase the abundances derived for the accreted material need to be modified to consider gravitational settling for the elements in the white dwarf photosphere, nX(A) par nX(B) par = τB τA X(A)WD X(B)WD , (1.3) where τA and τB are the diffusion time-scales of element A and B respectively through the white dwarfs photosphere. Declining phase – after the parent body has been fully accreted by the white dwarf, the abundances decrease exponentially with the decay factors depending on the time since accretion ceased, t, nX(A) par nX(B) par = nX(A)WD nX(B)WD e−t/τB e−t/τA . (1.4) 20 Introduction The sinking timescales of the heavy metals in the photospheres of white dwarfs varies depending both on the element sinking and the white dwarf properties. Warm DAZ white dwarfs have sinking timescales on the order of days–years, whereas DBZ white dwarfs have much longer sinking timescales on the order of 102 to 106 years. The sinking timescales are most often based on 1D models (e.g. Koester, 2009). However, recent works have been incorporating various mixing mechanisms (e.g. convective instability, thermohaline instability2) at the surfaces of white dwarfs in 3D atmosphere models, this has the effect of increasing the sinking timescale of metals in white dwarfs (e.g. Cunningham et al., 2019). As inferences about the pollutants relies on accurate knowledge of the sinking timescales, further work is required to confirm the sinking timescales. Differentiation: Fragments of core-mantle differentiated bodies have been observed in the atmospheres of polluted white dwarfs, where observations show enhancement or depletion in siderophilic elements (core-rich or mantle-rich respectively) (e.g. Zuckerman et al., 2011; Hollands et al., 2018; Harrison et al., 2018). Figure 1.7 shows that some white dwarfs are en- hanced in elements commonly found in the mantle (Mg, Si) compared to core (Fe), e.g. SDSS J1043+08556 (S1043), and others are enhanced in core elements, e.g. PG 0843+516 (PG0843). If the white dwarfs are accreting the collisional fragments of parent bodies that formed an iron core, a large (as much as two thirds) fraction of white dwarfs may have accreted fragments of larger core-mantle differentiated bodies (Bonsor et al., 2020). This implies that the process of (iron) core formation in exo-asteroids is likely ubiquitous across exoplanetary systems. Water: Water is widely recognised as one of the key components required for life, and therefore, there is much astrophysical interest in its detection. Water-rich exoplanetary material which pollutes white dwarfs cannot be directly detected in the white dwarf photospheres as it dissociates into oxygen and hydrogen. However, by calculating how much oxygen would be expected to be in the form of metal oxides (MgO, SiO2, Al2O3, CaO, and FeO), the oxygen remainder can be attributed to water (Klein et al., 2010). A handful of objects have accreted bodies rich in water-ice (Farihi et al., 2011, 2013; Raddi et al., 2015; Hoskin et al., 2020; 2The thermohaline instability arises due to unstable gradients in composition which become stabilised by a temperature gradient. Therefore, in polluted white dwarfs the thermohaline instability could lead to additional mixing as the significant mean molecular weight of the material that has accreted compared to the atmosphere of hydrogen below can lead to this instability arising. 1.4 Polluted White Dwarfs 21 Klein et al., 2021), Figure 1.7 shows WD 1425+540 (WD 1425), WD 1232+563 (WD 1232), SDSS 1242+5226 (S 1242), WD 0738 (WD J0738+1835), and GD 61 have excess oxygen and are inferred to be ice-rich objects. WD 1425+540 has accreted a fragment of a Kuiper-belt like icy body rich in both water and nitrogen ices (Xu et al., 2017). Trace hydrogen is detected in up to 75% of DB white dwarfs (Koester & Kepler, 2015). Gentile Fusillo et al. (2017) found a link between trace hydrogen and metal polluted white dwarfs, demonstrating another way to infer the accretion of water-rich planetesimals onto white dwarfs. Carbon: Carbon is a key species for the understanding of life. It can be in the form of graphite, organics or ices and is a volatile element, hence the amount of carbon in solid species in the solar system is depleted in comparison to the solar photosphere. There is only one report of carbon rich material, in Ton 345; however, the abundance is debated. Ground based spectroscopy revealed that it was carbon rich with a composition similar to Kuiper Belt objects with small amounts of ice due to the low oxygen abundance (Jura et al., 2015), whilst space based UV observations revealed a much lower abundance of carbon (Wilson et al., 2015). Therefore, there is no definitive evidence for any white dwarf accreting carbon-rich exoplanetary bodies, and all but one demonstrate a carbon depletion (Jura, 2006; Wilson et al., 2016). The body that does not show a carbon depletion is that polluting WD 1425+540, discovered by Xu et al. (2017), as mentioned in the previous section. This body was also the first exo-planetary debris to be detected with Nitrogen; this insinuates that this body formed in a Kuiper belt-like region. UV spectroscopy is most sensitive to detecting carbon, therefore with more UV spectra of polluted white dwarfs providing detections and stringent limits on the carbon abundance, further understanding about the form of carbon in exo-planetary worlds will be revealed. 1.4.3 Circumstellar Discs Observations of dust discs: Dust debris from tidally disrupted planetesimals has been discovered via excess infrared emission above that expected from the white dwarf around 1.5–4% of white dwarfs (e.g. Becklin et al., 2005; Kilic et al., 2006; Jura et al., 2007b; Rebassa-Mansergas et al., 2019; Wilson et al., 2019; Xu et al., 2020). G29-38, one of the brightest and most heavily polluted white dwarfs, was the first polluted white dwarf to be discovered with excess emission at infrared wavelengths, as seen in Figure 1.8. It was initially hypothesised that the excess was caused by a brown dwarf companion (Zuckerman & Becklin, 1987). However, pulsation 22 Introduction Fig. 1.7 Bar charts showing the pollutant mass fractions of: O, Mg, Si, Fe, C, for 19 polluted white dwarfs. The bulk compositions of the Earth and comet Halley are shown on the left for comparison. Figure from Xu & Bonsor (2021). timing in the optical and infrared found echoing signals which were difficult to explain with the presence of a brown dwarf companion. Also, infrared spectra demonstrated an excess emission at 10 µm, triple the intensity expected from a brown dwarf (Graham et al., 1990). These arguments favour the infrared excess resulted from a dust disc close to the white dwarf. The number of infrared excess detections remained at one for nearly two decades. Ground based observations discovered GD 362 as the second white dwarf with an infrared excess associated with a dust disc, with GD 40 closely following this (Becklin et al., 2005; Kilic et al., 2005, 2006). The launch of the Spitzer Space Telescope in 2003 exponentially increased the number of detected discs. To date, tens of white dwarfs have been found with an infrared excess associated with a dust disc, and all are found to also be polluted, implying that the circumstellar material is accreting onto the atmosphere (e.g. Jura et al., 2007b). Infrared spectra of the dust can inform about the composition of the dust grains. Low signal-to-noise Spitzer spectra from the infrared spectrograph (IRS) of a number of polluted white dwarfs demonstrates strong 10 µm silicate emission features, originating from the orbiting dust disc (Jura et al., 2007a, 2009a). All spectra show the 10 µm feature, with the red wing extending to 12 µm. This feature is characteristic of micron-sized olivine glasses, associated with a tidally disrupted asteroid. G29-38 is the only dusty white dwarf bright enough to have a sufficiently high signal-to-noise infrared spectrum that an in depth mineralogy study could be conducted. Reach et al. (2009) reported a 1–35 µm spectrum, 1.4 Polluted White Dwarfs 23 Fig. 1.8 The spectral energy distribution of G29-38. It shows the first infrared excess discovered around a polluted white dwarf (Zuckerman & Becklin, 1987). The data was taken using the NASA Infrared Telescope Facility. An excess at near to mid-infrared wavelengths is clear. as shown in Figure 1.9, and determined the possible mineral compositions of the dust disc. Disc modelling found it likely that the disc is formed of amorphous carbon, amorphous and crystalline silicates, water ice and metal sulphides. Xu et al. (2014) compared the disc abundances with those in the white dwarf atmosphere and found that the overall agreement is consistent, except for a few elements. JWST will provide high quality mid-infrared (mid-IR) spectra that covers strong features from different minerals for discs fainter than G29-38. This will allow further investigation into the correlation between the atmospheric pollution and the disc material, and reveal both the chemical and mineral compositions of extrasolar planetesimals. It has been suggested that there may be observable signatures of dust beyond the tidal radius, in the form of an asteroid belt or Kuiper belt analogue. This could supply the planetesimals that eventually accrete onto the white dwarfs. Xu et al. (2013) used the Herschel Space Observatory to search for a parent population of dusty bodies around GD 362, the data resulted in null detections and stringent upper limits were placed. Farihi et al. (2014) used the Herschel Space Observatory and the Atacama Large Millimetre/submillimetre Array (ALMA) to search for this parent population around the closest and brightest dusty white dwarf, G29-38. Both Herschel Space Observatory and ALMA data resulted in null detections, 24 Introduction Em iss iv ity = O bs er ve d / B ν( 93 0 K) Wavelength (μm) am. carbon carbonates PAHwater vap. water ice sulfides phyllosilicates cryst. pyr. am. silicates cryst. forsterite 5 10 30 3515 20 25 0.0 0.5 1.0 Fig. 1.9 Spitzer/IRS spectrum of the dust around G29-38. The observed spectrum was divided by a blackbody of temperature 930 K to give the black data points. The best fitting linear combination of 12 minerals is plotted as the red dashed line, to show the potential composition of the dust. The contribution of each mineral is shown by the coloured lines (Reach et al., 2009). The composition of this dusty material should match the composition of the material polluting the star. Water ice, amorphous silicates and forsterites all have distinguishable features in the mid-IR. 1.4 Polluted White Dwarfs 25 putting strong upper limits on the amount of cold dust in the systems. Further observations of dusty white dwarfs using the full capabilities of ALMA could be beneficial. Contrary to this, there are detections of cold dust, however, they are discovered around hot white dwarfs ∼ 100,000 K. These hot systems are different to the cool metal rich white dwarfs. Su et al. (2007) discovered excess infrared emission from the central region of the Helix nebula. The central star is a hot white dwarf, WD 2226-210, with an effective temperature of ∼ 110,000 K. Spitzer detected excess emission at 8, 24 and 70 µm, which was hypothesised to be from a dust disc, with a blackbody temperature in the range 90–130 K. It is thought the dust originates from Kuiper belt-like objects, which have survived the violent evolution of the host star and collided to form a disc. This demonstrates it is possible for cold dust to survive and exist around white dwarf systems. Motivated by this work, Chu et al. (2011) used the Spitzer Multi-band Imaging Photometer to investigate a sample of 71 hot white dwarfs, or pre-white dwarfs to search for excess emission at 24 µm. Excess emission was discovered for a number of the sources so it can be hypothesised that dust at ∼ 100 K should be found around ≳ 15 % of white dwarfs and pre-white dwarfs. However, these cool infrared excesses are widely debated as a number of these objects have been discovered to host companions. Observations of gas discs: 21 of the white dwarfs with dust debris also show evidence of circumstellar gas in emission near the same radius as the dust (Gänsicke et al., 2006, 2007, 2008; Melis et al., 2010; Farihi et al., 2012a; Melis et al., 2012; Brinkworth et al., 2012; Debes et al., 2012b; Dennihy et al., 2020; Melis et al., 2020; Gentile Fusillo et al., 2020). These systems are identified by their double peaked emission features, usually strongest at the Ca II infrared triplet (8498 Å, 8542 Å, and 8662 Å). WD 1226+110 was the first of such detections, where Ca II emission lines in the optical and near infrared (NIR) were observed. The Ca II emission line profiles for WD 1226+110 are shown in Figure 1.10. The double peaked spectral lines are characteristic of Doppler broadening due to a Keplerian rotating gaseous disc. This has a rotational velocity law of ΣKαr−3/2 where the inner disc rotates faster than the outer disc. Gaia J0611−6931 has the most elements detected in emission, with observations of Ca, O, Si, Mg, Na and Fe (Dennihy et al., 2020; Melis et al., 2020). As there are a lack of hydrogen and helium spectral features in the gaseous discs, this implies the gas discs are metallic. The only gas disc found to consist of predominantly volatiles (H, O, S) is that around WD J091405.30+191412.25 as reported in Gänsicke et al. (2019), where the chemical composition matches that found deep in the atmospheric layers of icy giant planets. 26 Introduction Fig. 1.10 The Ca II emission line profiles for WD 1226+110 from Gänsicke et al. (2006). The double peaked profiles are clear, arising due to Keplerian rotation. There are clear asymmetries in the red and blue Doppler shifted lines, indicating the disc is eccentric. Disc Variability: The scattering of planetary bodies that lead to pollution is expected to be a stochastic process, with the potential for variability on human time scales as seen in Figure 1.11 by the predicted order of magnitude changes in the accretion rate over timescales of less than a year (Wyatt et al., 2014). Such dust variability has previously been identified. The first infrared variability associated with the dust was for white dwarf WD J0959-0200. The fluxes in both Spitzer (Infrared Array Camera) channels, 3.6 µm and 4.5 µm, dropped by approximately 35 % within a year baseline. The NIR K band flux also dropped by 20 %, emphasising the dust is close to the star and dissipating on short time scales (Xu & Jura, 2014). Several scenarios could explain the large drop in flux for WD J0959-0200. One scenario involves planetesimals being scattered inward by a large outer planet and impacting a pre-existing disc. Scattering from a size distribution of small bodies can be approximated as continuous accretion, but infrequent scattering of larger bodies can affect the disc and accretion signatures, with the frequency of impact being be as often as once every 30 years or less (Jura, 2008; Wyatt et al., 2014). Mutual collisions between fragments can cause the inner disc to evaporate. If this hypothesis is supported, Xu & Jura (2014) suggested that WD J0959-0200 lost approximately 3 % of the disc mass within 300 days, implying that white dwarf accretion is intrinsically variable. 1.4 Polluted White Dwarfs 27 Fig. 1.11 Figure from Wyatt et al. (2014) showing model predictions for accretion rate over time for 20 stars modelled to be similar to G29-38. Changes in the accretion rate can be up to an order of magnitude over timescales of less than a year. 28 Introduction Investigations into whether the accretion rate of white dwarfs varies as a function of time have been conducted previously for one white dwarf, G29-38 (von Hippel & Thompson, 2007; Debes & López-Morales, 2008). Further observations of accretion rate over time will allow insight to be gained into whether this model is correct. Another scenario could explain the variability if the white dwarfs are observed with both circumstellar gas and dust. The dust and gas may interact and cause the inner disc to evaporate, resulting in non-steady state accretion. WD J0959-0200 has both gas and dust discs, emphasising this assumption is reasonable. However, variability has now been observed to be ubiquitous in the mid-infrared even for those systems without observable gas discs, with no preference for increases or decreases in dust flux as shown in Figure 1.12 (Swan et al., 2019a, 2020). Variability in the circumstellar gas has also been observed for white dwarfs with cir- cumstellar gaseous discs (Wilson et al., 2014, 2015; Manser et al., 2015, 2016b; Dennihy et al., 2018). Variability tends to be morphological variations in the line profile shape, or increases/decreases in the equivalent widths of the spectral features. The variability is mostly attributed to the precession of eccentric rings, collisions, or planetesimals orbiting within the debris discs (Manser et al., 2019). Jura (2003) Flat and Opaque Disc Model: In 2003 it was well known that metals polluted the atmospheres of some white dwarfs, however, the only known infrared excess was that of G29-38. Jura (2003) developed a model to explain both the atmospheric pollutants and the infrared excess for G29-38. Scattered planetesimals perturbed onto highly eccentric orbits can fall within the Roche radius of the white dwarf. This causes them to become tidally disrupted, and over time the system relaxes into a flat orbiting dust disc. It is assumed that the white dwarf illuminates the thin, flat, optically thick disc and the incident flux is re-radiated in the infrared. The thermal emission of the discs are of the order 103 K, demonstrating these discs lie close to their host star. The closest solar system analogue are the rings of Saturn; the physics of these rings can mostly be applied to white dwarfs (Brahic, 1977; Dones, 1991). The temperature of the white dwarf disc, Tdisc, as a function of distance from the star, R, which has a stellar radius of R∗ and a temperature of T∗, is given by: Tdisc ≈ ( 2 3π )1/4 (R∗ R )3/4 T∗. (1.5) 1.4 Polluted White Dwarfs 29 Fig. 1.12 Figure from Swan et al. (2020) showing the mid-infrared variability observed in a sample of ∼ 40 white dwarfs using the Spitzer space telescope. The largest flux changes seem to occur over longer baselines with the largest changes being those with circumstellar gaseous discs. 30 Introduction The flux from the disc, is predicted to be: Fdisc = 2π cos i D2∗ ∫ Rout Rin Bν (Tdisc)RdR, (1.6) where D∗ is the distance to the white dwarf, i is the inclination angle of the disc, Rin and Rout are the inner and outer radii of the dust disc respectively and Bν is the Planck function. In this model, the inner radius is defined by the sublimation radius of solids. The sublimation radius, RS, is given by: RS = √ ε R∗ 2 ( T∗ TS )2 , (1.7) where ε is the ratio of the emissivity of the particle for the white dwarf radiation and for its own thermal radiation, and TS is the sublimation temperature for the particles in the disc (Rafikov, 2011a). The dust ring is expected to form from the tidal disruption of an asteroid, therefore the outer radius is defined by the tidal radius of the white dwarf. This radius, Rtide, is given by: Rtide =Ctide ( ρ∗ ρa )1/3 R∗, (1.8) where ρ∗ is the density of the white dwarf, ρa is the density of the asteroid and Ctide is a constant, O(1), dependent on the composition, rotation and orbit of the asteroid (Davidsson, 1999). The process of accretion from the disc onto the stellar atmosphere is complicated and refuted. After the asteroid is tidally disrupted, collisions result in the dissipation of energy, such that the dust settles onto a flat plane. Veras et al. (2015) established that the radiation from the white dwarf can circularise and compress the orbiting dust. The radiation from the white dwarf drives matter towards the sublimation radius via Poynting-Robertson (PR) drag (Robertson, 1937; Burns et al., 1979) shown in Figure 1.13. This delivers material to a reservoir of metallic gas close to the white dwarf. Viscous torques provide the means for metals to accrete onto the white dwarf surface (Rafikov, 2011a). Rafikov (2011b) demonstrated that PR drag is sufficient to achieve metal accretion rates of at least 108 g s−1. But this is not sufficient to explain the accretion rates of the highest accretors (∼1011 g s−1). Paradigm Shift: The model of the flat opaque disc as developed in Jura (2003) was revolutionary as it provided a framework to understand these enigmatic polluted white dwarf systems. However, there is increasing evidence for the requirement of an optically thin component, moving the field away from the flat opaque disc model. The 10 µm silicate emission features detected for 1.4 Polluted White Dwarfs 31 a number of white dwarfs with infrared excesses must result from an optically thin region (Jura et al., 2009a). In order to fit the strong 10 µm features and the shape of the infrared excesses, GD 362 (Jura et al., 2007a) and G29-38 (Reach et al., 2009) were fitted with models which invoked warped/flared discs. Large infrared variability is observed in a number of systems with both increases and decreases observed. In order to explain this collisional models are proposed where collisions cause infrared brightening as dust is produced, and then dimming occurs as the cascade grinds the material down such that it no longer emits at micron wavelengths (Farihi et al., 2018; Swan et al., 2020, 2021). The tidal disruption model implies that observable dust discs should be more frequent than is detected. Bonsor et al. (2017) found that no more than 3.3% of white dwarfs can have a wide, flat, opaque dust disc and suggest four reasons to explain the lack of discs observed. These reasons are extracted from Figure 1.14 where the accretion rate is plotted against the white dwarf temperature. The green region (A) may represent white dwarfs which have already accreted their disc. This is possible if a white dwarf has a long sinking time scale such that pollutants exist in the atmosphere after the full dissipation of the disc. The blue region (B) may represent white dwarfs which are surrounded by undetectable optically thin dust. The brown region (C) may explain the lack of infrared excess for hot white dwarfs where the sublimation radius of solids lies outside of the tidal disruption radius. This means that gas is accreted onto the stellar photosphere without a disc forming (Steckloff et al., 2021). The red region (D) represents an area of phase space where there is enhanced accretion, although it is difficult to explain the lack of infrared detections here. Although further work is required to confirm these hypotheses, this work emphasises the need for a paradigm shift away from the optically thick, razor thin model. As previously stated, 21 systems have been observed with both circumstellar dust and gas discs. These gas discs extend to a similar radii as the dust discs (e.g. Melis et al., 2012). This is surprising as gas is observed at radii larger than the sublimation radius for refractory elements. The vertical extent of the gas disc appears to be much greater than the flat, opaque dust disc. Different theories have been proposed to explain the production of the gas - sublimation of solids, and collisions. Metallic gas inwards of the sublimation radius feeds the stellar atmosphere via viscous torques. To conserve angular momentum a proportion of the gas produced at the sublimation radius viscously spreads outwards causing an overlap in the location of the dust and gas (Rafikov, 2011a; Metzger et al., 2012). This gas disc has a higher temperature than the dust disc enabling the gas disc to exist outside of the sublimation radius. This outwardly spreading gas causes drag on the dust particles and thus accelerates their accretion onto the white dwarf creating a runaway effect; this might explain the highest accretion rates observed in polluted white dwarfs (Rafikov, 2011a). To 32 Introduction Fig. 1.13 A diagram from Rafikov (2011a) illustrating an interacting dust disc (black) and gaseous disc (grey). The radial extent of the dust disc is defined on the inner edge by the sublimation radius of solids and the outer edge by the tidal disruption radius of the white dwarf. The gas disc has the same radial extent as the dust, with a reservoir close to the white dwarf, feeding the atmosphere with metallic gas. The movement of high metallicity materials is demonstrated by the arrows. investigate this, the accretion rates of white dwarfs with circumstellar gas needs to be studied in comparison to the population. A schematic of this theory is shown in Figure 1.13. An alternative explanation for the gas emission is collisional cascades of planetesimals within the Roche radius of the white dwarf, these collisions produce gaseous material (Kenyon & Bromley, 2017a,b). Observations of infrared variability in WD 0145+234 appear consistent with simple collisional cascade models (Swan et al., 2021). It is likely numerous routes of planetesimals arriving in the white dwarf atmosphere are active to explain the diversity in the architecture of the systems and accretion rates, a more thorough discussion appears in Brouwers et al. (2022). Dust provides information regarding how planetary material arrives into the atmosphere of these polluted white dwarfs. No theoretical models can encompass all dust observations, including the shapes of the SEDs, the silicate emission features, and the dust variability. Dust variability studies are crucial as they can reveal information about the optical depth and geometry of the circumstellar dust. However, as of 2016 (when this thesis work begun), studies of dust variability only focused on individual systems rather than the population, so creating theoretical models that can explain the population was impossible. Therefore, it is crucial to determine how frequent variability is within the population and at what wavelengths this variability is dominant. 1.4 Polluted White Dwarfs 33 Fig. 1.14 This figure is taken from Bonsor et al. (2017). It demonstrates the accretion rate plotted against the temperature of the white dwarfs. The plot has been subdivided into four regions to explain the lack of observations of infrared excesses around polluted white dwarfs. 1.4.4 The Link with the Main Sequence The asteroid destruction model relies on the survival of at least one large planet as well as planetesimals from the main sequence to the white dwarf phase. As nearly all stars are expected to end their lives as white dwarfs, looking at the planetary systems around main sequence stars can reveal what could survive. There are over 5,000 exoplanets confirmed to date, with a large range of radius, mass, semi-major axis and orbital period. Due to this diversity, it is likely that some exoplanets have the necessary location and properties to survive to the white dwarf phase. The post main sequence evolution of a star is a violent process, however dynamicists have hypothesised that certain planets can survive. Mustill & Villaver (2012) argued that terrestrial planets need an initial semi-major axis of at least 2.8 AU and gas giants need an initial semi-major axis of 5 AU to avoid being engulfed by the expanding envelope of the giant star. During the asymptotic giant branch phase, the star loses a large fraction of its initial mass. This mass loss occurs over such long time scales that it is possible for objects to have stable expanding orbits (Debes & Sigurdsson, 2002). Therefore, if planets have a sufficiently large semi-major axis, they could survive to the white dwarf phase and perturb asteroids into highly eccentric orbits. The fate of this asteroid could be: a collision with another surviving body, to be ejected from the planetary system, to directly hit the white dwarf, or to fall within the tidal disruption radius of the white dwarf (Villaver & Livio, 2007; 34 Introduction Veras et al., 2011, 2013; McDonald & Veras, 2021). A fraction of main sequence planets have been observed with these large semi-major axis, so the idea that they could survive to the white dwarf phase is reasonable. Bonsor & Wyatt (2010) found that debris belts with sufficiently massive planetesimals that lie beyond a few AU can survive to the white dwarf phase. However, the masses of the discs are reduced significantly; this is mainly due to the loss of small planetesimals. The small planetesimals are subsequently replenished via collisions once the system has settled into the white dwarf phase. It has been shown that 25–50% of white dwarfs are polluted (e.g. Zuckerman et al., 2003). Therefore, the study of these objects could provide information about main sequence exo-planetary systems, if we assume that this observed pollution is due to scattered planetesimals. This implies that at least 25–50% of low mass main sequence systems must have a massive planetesimal belt with enough debris to survive to the white dwarf phase and pollute the atmosphere of the white dwarf. Cold outer belts around main sequence stars have detection rates of 30% around A stars and about 20% for FGK stars (Marino, 2022). This is broadly consistent with the fraction of white dwarfs found to be polluted. Sub-stellar bodies have been discovered around white dwarfs. Numerous brown dwarf companions to white dwarfs have been discovered (Longstaff et al., 2019; van Roestel et al., 2021), but a limited number of planets have been found. Most notably, the giant transiting planet candidate, WD 1856+534b (Vanderburg et al., 2020). Although it is expected that planets around white dwarfs should be commonplace in order for 25–50% of white dwarfs to be polluted, few other planet candidates have been observed (Xu et al., 2015a; Belardi et al., 2016). This is unsurprising as the probability of observing a transiting planet is low due to small radii of white dwarfs and surviving planets are predicted to be far from the star (Debes & Sigurdsson, 2002). Gaia will enable planets to be discovered around white dwarfs astrometrically. Gaia eDR3 resulted in a candidate super-Jupiter around WD 0141−675, if confirmed this would be the first polluted white dwarf with an exoplanet detected (Collaboration et al., 2022). With the release of Gaia DR5 it is expected that Gaia will astrometrically discover 6±1 planets between 1.6–3.91 au with masses between 0.03 – 13 MJ (Sanderson et al., 2022) allowing investigations into the evolution of planetary systems from the main sequence to white dwarf phase. 1.5 This Thesis in Context Polluted white dwarfs are unique laboratories that can be used to directly sample the bulk composition of exoplanets; this is not possible via any other method. We are transitioning 1.5 This Thesis in Context 35 into the era of big data. For white dwarfs, Gaia has aided the discovery of > 350,000 white dwarfs (Gentile Fusillo et al., 2021), and with large spectroscopic surveys upcoming such as 4MOST, DESI, WEAVE, and SDSS-V, thousands of new polluted white dwarfs will be discovered. Each white dwarf tells its own unique story about the exoplanetary material surrounding it, therefore, it will be crucial to follow up with high resolution instruments across large wavelength ranges to ensure a thorough characterisation and understanding of this exoplanetary material. These individual studies have led to ground-breaking discoveries regarding the composition and geology of planetary material around other stars, such as the discovery of water-rich exoplanetary bodies (e.g. Farihi et al., 2011), which in turn has prompted further investigations into the possibility of life in exoplanetary systems. For exoplanets, over the next decades a myriad of rocky planets will be discovered and characterised with current missions such as TESS, CHEOPS, and JWST; and future missions, such as: PLATO, LUVOIR, and LIFE. With this wealth of new data, collaboration will be crucial, information about exoplanetary material from polluted white dwarf systems can help to answer big questions about the solar system, exoplanetary systems, and galactic populations, and similarly these subjects can help inform inferences from white dwarf planetary systems. In order to fully understand and characterise exoplanetary rocky worlds, as well as determine whether they could be habitable, it is crucial to understand what they are made of. Current models of the interiors of exoplanets suffer from degeneracies; however, the population of the bulk composition of exoplanetary material, as inferred from polluted white dwarfs, can be used as inputs into planetary interior modelling to help to break these degeneracies. With the wealth of data coming it is crucial for sub-fields to collaborate such that a new era of Earth 2.0 discovery can begin. Our place in the solar system and galaxy appears to be key to life developing on Earth, but it is unknown whether the solar system is common, or how planet formation is affected by our place in the galaxy. Polluted white dwarfs can reveal the population of exoplanetary material which can be used to answer these questions. Combining the inferred exoplanetary composition with Gaia astrometric information reveals whether different galactic populations have different planet compositions and the affects this has on the resultant planetary systems. This thesis focuses on white dwarf planetary systems, with an emphasis on using obser- vations to inform models of the disruption and accretion processes, and using large ground based spectrographs to make inferences about exoplanetary composition of individual sys- tems. Chapter 2 reports large ground-based near-infrared photometric monitoring campaigns of the circumstellar dust around ∼ 40 white dwarfs. There is increasing evidence for a paradigm shift away from the model where a planetesimal tidally disrupts and forms an opaque flat disc (Jura, 2003). Dust variability sheds light on the dynamics of the planetary 36 Introduction bodies, the distribution of the dust, and the optical depth and properties of the dust. Chapter 3 reports an optical spectroscopic monitoring campaign of 5 DAZ white dwarfs over decade long timescales. Stochastic models of the accretion of planetesimals predicts variability in the accretion rate, and therefore the amount of material in the atmosphere. Warm DAZ white dwarfs have accretion rates on the order of days, therefore enabling the search for variability in the accretion rate using spectroscopic monitoring. Constraining accretion rates provides valuable information to inform models of accretion and disc lifetimes. Chapter 4 reports the abundances for the material accreting onto seven white dwarf stars and the composition and geological history of these exoplanetary bodies. These findings can be used to aid interior structure models of exoplanets, investigate whether the composition of solar system bodies are unique, and study volatiles and their delivery in exoplanetary systems. Finally, Chapter 5 presents the conclusions from this thesis and provides suggestions for future work to contribute to questions raised. Chapter 2 Dust Variability The favoured theory regarding the source of white dwarf pollution is that it originates from planetesimals scattered onto eccentric orbits towards the white dwarf where they tidally disrupt, producing dusty debris detectable as excess infrared emission, and subsequently accrete onto the atmosphere of the white dwarf. This inwards scattering of planetary bodies that leads to pollution is expected to be a stochastic process (Wyatt et al., 2014), with variability predicted on human time-scales. Variations could come in the form of gaseous emission, dusty emission, or total flux caused by transiting debris. Studying the circumstellar dust via time series photometry of the infrared emission can reveal how this planetary material arrives in their atmospheres. The first white dwarf discovered to show dust emission variability was WD J0959−0200. The 3–5 µm flux dropped by 35% between observations separated by a year and appeared to be stable afterwards, and the near-infrared K band flux dropped by 18.5% between observations separated by 8 years (Xu & Jura, 2014). The large K band excess for this white dwarf implies that the dust is hot and close-in. These findings imply that the variability event was large and became stable quickly. Motivated by the large variability in WD J0959−0200, ground based near-infrared monitoring campaigns were conducted to search for variability in the dust for white dwarfs with infrared excesses. This chapter focuses on these monitoring campaigns and characterising how white dwarfs vary in the near-infrared in order to better understand how dusty material arrives in the atmospheres of the white dwarfs. Section 2.1 presents a near-infrared monitoring campaign of 34 white dwarfs with circumstellar dust emission using the UKIRT/WFCAM to search for variability on timescales of months – years. Section 2.2 presents a near-infrared monitoring campaign of 8 white dwarfs with dust emission using the NTT/SOFI to search for variability over timescales of hours – days. 38 Dust Variability 2.1 The UKIRT/WFCAM near-infrared monitoring A near-infrared monitoring campaign of 34 white dwarfs with infrared excesses was con- ducted with the aim to search for variability in the dust emission. Time series photometry of these white dwarfs from the United Kingdom Infrared Telescope (Wide Field Camera) in the J, H and K bands were obtained over baselines of up to three years. No statistically significant variation is found in the dust emission in all three near-infrared bands. Specifically, variability at ∼ 1.3% can be ruled out for the 13 white dwarfs brighter than 16th mag in K band, and at ∼ 10% for the 32 white dwarfs brighter than 18th mag over time-scales of three years. Although to date two white dwarfs, SDSS J095904.69−020047.6 (Xu & Jura, 2014) and WD 1226+110 (Xu et al., 2018b), have shown K band variability, in the sample no evidence of new K band variability at these levels is seen. One interpretation is that the tidal disruption events which lead to large variabilities are rare, occur on short time-scales, and after a few years the white dwarfs return to being stable in the near-infrared. 2.1.1 UKIRT Observations An infrared monitoring campaign was completed using the Wide Field Camera (WFCAM) on the UKIRT on Maunakea, Hawaii (PI: Xu). WFCAM operates in the near-infrared between 0.83 and 2.37 µm ; this includes the J, H and K filters (Hewett et al., 2006). WFCAM has four detectors composed of 2048 by 2048 18 micron pixels, with a pixel scale of 0.4′′. The sample consists of 34 white dwarfs with infrared excesses. Table 2.1 shows the sample, giving the white dwarf parameters and the nature of the infrared excess. In the sample, 30 white dwarfs have dust emission confirmed with Spitzer observations, the 4 remaining objects are identified by WISE and are labelled as dust candidates. JHK broadband photometry was obtained for the 34 targets with the UKIRT. 48 hours of UKIRT WFCAM time between semesters 2014B and 2017A were used to study long-term variability in the near-infrared. The dates of the observations for the 34 white dwarfs are shown in Table 6.1. Each white dwarf was observed up to 6 times across the 3 year campaign. The UKIRT observations probed short time-scales (minutes) between individual frames taken on a given night, and long time-scales (years) between observation dates. Further discussions about time-scales are given in Section 2.1.3. The observations were designed such that a signal- to-noise ratio (SNR) of 30 was obtained for the J and H bands, and 40 for the K band. To achieve the required SNR, the brightest white dwarfs in the sample required 25 s of exposure time for the J and H band, and 75 s for the K band. For the faintest white dwarfs, this required 150 s for the J band, 450 s for the H band, and 1250 s for the K band. As the infrared sky is bright, each frame consisted of a dithered stack of five 5 s or 10 s exposures. Table 6.2 2.1 Dust variability 39 Table 2.1 The UKIRT sample of dusty polluted white dwarfs and their properties. H/He column highlights the dominant element in the atmosphere. WD Name Other Name H/He Te f f ⋆ log(g)⋆ Nature of Excess † WD 0010+280 PG 0010+281 H 23700 7.731 Dustα WD 0106−328 HE 0106−3253 H 16200 8.031 Dustβ WD 0146+187 GD 16 He 11500 8.002 Dustγ WD 0300−013 GD 40 He 13600 8.023 Dustδ WD 0307+077 HS 0307+0746 H 10100 7.991 Dustβ WD 0408−041 GD 56 H 14200 7.961 Dustδ WD 0435+410 GD 61 He 15700 8.043 Dustε WD J0738+1835∗ SDSS J073842.57+183509.6 He 14000 8.404 Dustζ WD J0959−0200∗ SDSS J095904.69−020047.6 H 13300 8.065 Dustη WD 1015+161 PG 1015+161 H 19200 8.031 Dustδ WD 1018+410 PG 1018+411 H 21700 8.051 Dustθ WD 1041+092 Gaia DR2 3869060540584643328 H 17600 8.121 Dustζ WD 1116+026 GD 133 H 12600 8.046 Dustδ WD 1145+017 HE 1145+0145 He 15900 8.007 Dust ι WD 1145+288 Gaia DR2 4019789359821201536 H 12400 7.961 Candκ WD 1150−153 EC 11507−1519 H 11400 7.981 Dustλ WD J1221+1245∗ SDSS J122150.81+124513.3 H 12300 8.205 Dustη WD 1225−079 PG 1225−079 He 10800 8.008 Dustβ WD 1226+110 Gaia DR2 3904415787947492096 H 20900 8.111 Dustµ WD 1232+563 Gaia DR2 1571584539980588544 He 11800 8.303 Candκ WD 1349−230 HE 1349−2305 He 18200 8.139 Dustν WD 1456+298 G166-58 H 7390 8.001 Dustξ WD 1504+329 Gaia DR2 1288812212565231232 H 7180 8.0610 Candκ WD 1536+520 Gaia DR2 1595298501827000960 H 16700 7.706 Candκ WD 1551+175 Gaia DR2 1196531988354226560 He 15600 7.9611 Dustπ WD 1554+094 KUV 15519+1730 H 21800 7.6312 BD+Dustρ WD J1617+1620∗ SDSS J161717.04+162022.4 H 13200 8.041 Dustζ WD 1729+371 GD 362 H 10300 8.136 Dustξ WD 1929+011 Gaia DR2 4287654959563143168 H 25400 8.171 Dustθ WD 2132+096 HS 2132+0941 H 13100 7.9613 Dustπ WD 2207+121 Gaia DR2 2727904257071365760 He 14800 7.973 Dustσ WD 2221−165 HE 2221−1630 H 9900 8.121 Dustβ WD 2326+049 G29-38 H 11200 8.001 Dustτ WD 2328+107 PG 2328+108 H 17900 7.741 Dustθ Notes: ⋆ Temperature and log(g): Temperature and log(g) come from photometric fits, where possible, otherwise a new fit considering the Gaia distances was performed (see Section 2.1.3). (1) This work, (2) Koester et al. (2005), (3) Coutu et al. (2019), (4) Dufour et al. (2012), (5) Farihi et al. (2012a), (6) Dufour et al. (2017), (7) Vanderburg et al. (2015), (8) Klein et al. (2011), (9) Voss et al. (2007), (10) Barber et al. (2014), (11) Bergeron et al. (2011), (12) Farihi et al. (2017), (13) Gentile Fusillo et al. (2018). † Nature of infrared excess references: ‘dust’ refers to white dwarfs with confirmed dust detection with Spitzer observations and the reference is that of the Spitzer observations, otherwise it is labelled ‘cand’ (candidate), and the reference is the paper that first refers to dust emission. WD 1554+094 has a brown dwarf companion and dust. (α) Xu et al. (2015b), (β ) Farihi et al. (2010b), (γ) Farihi et al. (2009), (δ ) Jura et al. (2007b), (ε) Farihi et al. (2011), (ζ ) Brinkworth et al. (2012), (η) Farihi et al. (2012a), (θ ) Rocchetto et al. (2015), (ι) Xu et al. (2018a), (κ) Debes et al. (2011), (λ ) Jura et al. (2009a), (µ) Brinkworth et al. (2009), (ν) Girven et al. (2012), (ξ ) Farihi et al. (2008), (π) Bergfors et al. (2014), (ρ) Farihi et al. (2017), (σ ) Xu & Jura (2012), (τ) Reach et al. (2005). 40 Dust Variability (Appendix 1: Chapter 6) shows the exposure time for each of the five individual exposures making up the dithered stack. All data obtained with WFCAM were pipeline-processed by the Cambridge Astronomical Survey Unit using standard infrared photometry data reduction steps (CASU, Irwin et al., 2004; Dye et al., 2006). Further processing using the LIGHTCURVES software1 was executed to improve the precision of the photometry (Irwin et al., 2007). List driven photometry was performed using a master frame to force the centroid positions for the objects. Here, the master frame is the stacked image of all frames; this increases the SNR of the master frame and reduces errors associated with inaccurately placed centroids. Each frame used in the construction of the lightcurves was re-aligned to the same set of astrometric calibrators (from 2MASS) as the master frame, which gives significantly improved precision in aperture location when compared to centroiding for faint sources (see Irwin et al., 2007). For a particular object, all frames had an aperture between 3–5 pixels (1.2–2′′), depending on which aperture gave the lowest root mean square (RMS) value. Using a large number of non-variable stars, the zero-point shift for each frame was calculated to reduce atmospheric effects, and a 2D polynomial was fitted to the magnitude residuals to minimise the error associated with the position on the detector. LIGHTCURVES outputs a robustly calibrated lightcurve for each filter and for each object in the catalogue list. For variability searches this is essential as it allows us to probe changes down to the percent level. 2.1.2 Variability Analysis Analysis Method This work aimed to obtain high precision JHK photometry of a sample of dusty white dwarfs and to search for and constrain the level of variability in the photometry over the length of the survey. Throughout the analysis the Vega magnitude system was adopted. In order to robustly study the statistics of the photometry, the same analysis method was applied to all white dwarf observations regardless of the magnitude of the white dwarf and the number of measurements in the J, H and K bands. Gaussians were fitted to the distribution of photometrically corrected magnitudes from the LIGHTCURVES software using all dithered stacked frames, as demonstrated in Fig. 2.1. The total number of dithered stacked frames for each filter is shown in Table 6.2 (Appendix 1: Chapter 6), this ranges from 3 to 125, with the median number of frames for the J, H and K bands being 9, 12 and 20 respectively. A Markov Chain Monte Carlo (MCMC) approach was implemented to model the magnitude distribution of each star with a Gaussian profile, yielding posterior distributions for the 1https://github.com/mdwarfgeek/lightcurves 2.1 Dust variability 41 mean and standard deviation thereof. In the analysis, the median values from the posteriors are adopted. Uniform priors were used for the magnitude between 0 and 20, and standard deviation between 0 and 2. Larger upper values for the assumed prior range are used to test the sensitivity on the parameters. The results remained consistent independent of the size of the prior. A python package, PYMC2 (Patil et al., 2010), was used for MCMC parameter estimation. Five walkers were used, each resulting in 100,000 posterior samples, with the first 40% of the chain being discarded in the burn-in phase. The Gelman-Rubin statistic (Gelman et al., 1992) was implemented to check for chain convergence. For the two parameters (magnitude and standard deviation), all white dwarfs had values 1.000