Researcher profile

Daniel Kitzmann

Daniel Kitzmann contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2022arXiv

3D radiative-transfer for exoplanet atmospheres. gCMCRT: a GPU accelerated MCRT code

Radiative-transfer (RT) is a key component for investigating atmospheres of planetary bodies. With the 3D nature of exoplanet atmospheres being important in giving rise to their observable properties, accurate and fast 3D methods are required to be developed to meet future multi-dimensional and temporal data sets. We develop an open source GPU RT code, gCMCRT, a Monte Carlo RT forward model for general use in planetary atmosphere RT problems. We aim to automate the post-processing pipeline, starting from direct global circulation model (GCM) output to synthetic spectra. We develop albedo, emission and transmission spectra modes for 3D and 1D input structures. We include capability to use correlated-k and high-resolution opacity tables, the latter of which can be Doppler shifted inside the model. We post-process results from several GCM groups including ExoRad, SPARC/MITgcm THOR, UK Met Office UM, Exo-FMS and the Rauscher model. Users can therefore take advantage of desktop and HPC GPU computing solutions. gCMCRT is well suited for post-processing large GCM model grids produced by members of the community and for high-resolution 3D investigations.

preprint2022arXiv

Astrometric Accelerations as Dynamical Beacons: Discovery and Characterization of HIP 21152 B, the First T-Dwarf Companion in the Hyades

Benchmark brown dwarf companions with well-determined ages and model-independent masses are powerful tools to test substellar evolutionary models and probe the formation of giant planets and brown dwarfs. Here, we report the independent discovery of HIP~21152~B, the first imaged brown dwarf companion in the Hyades, and conduct a comprehensive orbital and atmospheric characterization of the system. HIP~21152 was targeted in an ongoing high-contrast imaging campaign of stars exhibiting proper motion changes between Hipparcos and Gaia, and was also recently identified by Bonavita et al. (2022) and Kuzuhara et al. (2022). Our Keck/NIRC2 and SCExAO/CHARIS imaging of HIP~21152 revealed a comoving companion at a separation of $0.37^{\prime\prime}$ (16 au). We perform a joint orbit fit of all available relative astrometry and radial velocities together with the Hipparcos-Gaia proper motions, yielding a dynamical mass of $24^{+6}_{-4}\,\mathrm{M_{Jup}}$, which is $1{-}2σ$ lower than evolutionary model predictions. Hybrid grids that include the evolution of cloud properties best reproduce the dynamical mass. We also identify a comoving wide-separation ($1837^{\prime\prime}$ or $7.9 \times 10^4 \, \mathrm{au}$) early-L dwarf with an inferred mass near the hydrogen-burning limit. Finally, we analyze the spectra and photometry of HIP~21152~B using the Saumon & Marley (2008) atmospheric models and a suite of retrievals. The best-fit grid-based models have $f_{\mathrm{sed}}=2$, indicating the presence of clouds, $T_{\mathrm{eff}}=1400 \, \mathrm{K}$, and $\log{g}=4.5 \, \mathrm{dex}$. These results are consistent with the object's spectral type of $\mathrm{T0\pm1}$. As the first benchmark brown dwarf companion in the Hyades, HIP~21152~B joins the small but growing number of substellar companions with well-determined ages and dynamical masses.

preprint2022arXiv

Retrieval Study of Brown Dwarfs Across the L-T Sequence

A large suite of 228 atmospheric retrievals is performed on a curated sample of 19 brown dwarfs spanning the L0 to T8 spectral types using the open-source Helios-r2 retrieval code, which implements the method of short characteristics for radiative transfer and a finite-element description of the temperature-pressure profile. Surprisingly, we find that cloud-free and cloudy (both gray and non-gray) models are equally consistent with the archival SpeX data from the perspective of Bayesian model comparison. Only upper limits for cloud properties are inferred if log-uniform priors are assumed, but the cloud optical depth becomes constrained if a uniform prior is used. Water is detected in all 19 objects and methane is detected in all of the T dwarfs, but no obvious trend exists across effective temperature. As carbon monoxide is only detected in a handful of objects, the inferred carbon-to-oxygen ratios are unreliable. The retrieved radius generally decreases with effective temperature, but the values inferred for some T dwarfs are implausibly low and may indicate missing physics or chemistry in the models. For the early L dwarfs, the retrieved surface gravity depends on whether the gray or non-gray cloud model is preferred. Future data are necessary for constraining cloud properties and the vertical variation of chemical abundances, the latter of which is needed for distinguishing between the chemical instability versus traditional cloud interpretation of the L-T transition.

preprint2022arXiv

Sunbathing under white light -- 3D modelling of brown dwarf - white dwarf atmospheres with strong UV irradiation

The atmospheres of brown dwarfs orbiting in close proximity to their parent white dwarf represent some of the most extreme irradiated environments known. Understanding their complex dynamical mechanisms pushes the limits of theoretical and modelling efforts, making them valuable objets to study to test contemporary understanding of irradiated atmospheres. We use the Exo-FMS GCM to simulate the brown dwarfs WD0137-349B, SDSS J141126.20+200911.1B and EPIC212235321B, first coupled to a multi-banded grey radiative-transfer scheme then a spectral correlated-k scheme with high temperature opacity tables. We then post-process the GCM results using gCMCRT to compare to available observational data. Our GCM models predict strongly temperature inverted atmospheres, spanning many decades in pressure due to impact of UV band heating. Post-processing of our models suggest that the day-night contrast is too small in the GCM results. We therefore suggest that the formation of cloud particles as well as atmospheric drag effects such as magnetic drag are important considerations in setting the day-night temperature contrast for these objects.

preprint2022arXiv

The PEPSI Exoplanet Transit Survey (PETS) I: Investigating the presence of a silicate atmosphere on the super-Earth 55 Cnc e

The study of exoplanets and especially their atmospheres can reveal key insights on their evolution by identifying specific atmospheric species. For such atmospheric investigations, high-resolution transmission spectroscopy has shown great success, especially for Jupiter-type planets. Towards the atmospheric characterization of smaller planets, the super-Earth exoplanet 55 Cnc e is one of the most promising terrestrial exoplanets studied to date. Here, we present a high-resolution spectroscopic transit observation of this planet, acquired with the PEPSI instrument at the Large Binocular Telescope. Assuming the presence of Earth-like crust species on the surface of 55 Cnc e, from which a possible silicate-vapor atmosphere could have originated, we search in its transmission spectrum for absorption of various atomic and ionized species such as Fe , Fe+, Ca , Ca+, Mg and K , among others. Not finding absorption for any of the investigated species, we are able to set absorption limits with a median value of 1.9 x RP. In conclusion, we do not find evidence of a widely extended silicate envelope on this super-Earth reaching several planetary radii.

preprint2022arXiv

Titanium oxide and chemical inhomogeneity in the atmosphere of the exoplanet WASP-189b

The temperature of an atmosphere decreases with increasing altitude, unless a shortwave absorber exists that causes a temperature inversion. Ozone plays this role in the Earth`s atmosphere. In the atmospheres of highly irradiated exoplanets, shortwave absorbers are predicted to be titanium oxide (TiO) and vanadium oxide (VO). Detections of TiO and VO have been claimed using both low and high spectral resolution observations, but later observations have failed to confirm these claims or overturned them. Here we report the unambiguous detection of TiO in the ultra-hot Jupiter WASP-189b using high-resolution transmission spectroscopy. This detection is based on applying the cross-correlation technique to many spectral lines of TiO from 460 to 690 nm. Moreover, we report detections of metals, including neutral and singly ionised iron and titanium, as well as chromium, magnesium, vanadium and manganese (Fe, Fe+, Ti, Ti+, Cr, Mg, V, Mn). The line positions of the detected species differ, which we interpret as a consequence of spatial gradients in their chemical abundances, such that they exist in different regions or dynamical regimes. This is direct observational evidence for the three-dimensional thermo-chemical stratification of an exoplanet atmosphere derived from high-resolution ground-based spectroscopy.

preprint2021arXiv

A Comparative Study of Atmospheric Chemistry with VULCAN

We present an update of the open-source photochemical kinetics code VULCAN (Tsai et al. 2017; https://github.com/exoclime/VULCAN) to include C-H-N-O-S networks and photochemistry. Additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross-sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an intercomparison of HD 189733b models between Moses et al. (2011), Venot et al. (2012), and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep troposphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and ammonia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP-33b, HD 189733b, GJ 436b, and 51 Eridani b. We look into the effects of the C/O ratio and chemistry of titanium/vanadium species for WASP-33b; we revisit HD 189733b for the effects of sulfur and carbon condensation; the effects of internal heating and vertical mixing ($K_{\textrm{zz}}$) are explored for GJ 436b; we test updated planetary properties for 51 Eridani b with S$_8$ condensates. We find sulfur can couple to carbon or nitrogen and impact other species such as hydrogen, methane, and ammonia. The observable features of the synthetic spectra and trends in the photochemical haze precursors are discussed for each case.

preprint2021arXiv

Physically-motivated basis functions for temperature maps of exoplanets

Thermal phase curves of exoplanet atmospheres have revealed temperature maps as a function of planetary longitude, often by sinusoidal decomposition of the phase curve. We construct a framework for describing two-dimensional temperature maps of exoplanets with mathematical basis functions derived for a fluid layer on a rotating, heated sphere with drag/friction, which are generalizations of spherical harmonics. These basis functions naturally produce physically-motivated temperature maps for exoplanets with few free parameters. We investigate best practices for applying this framework to temperature maps of hot Jupiters by splitting the problem into two parts: (1) we constrain the temperature map as a function of latitude by tuning the basis functions to reproduce general circulation model (GCM) outputs, since disk-integrated phase curve observations do not constrain this dimension; and (2) we infer the temperature maps of real hot Jupiters using original reductions of several Spitzer phase curves, which directly constrain the temperature variations with longitude. The resulting phase curves can be described with only three free parameters per bandpass -- an efficiency improvement over the usual five or so used to describe sinusoidal decompositions of phase curves. Upon obtaining the hemispherically averaged dayside and nightside temperatures, the standard approach would be to use zero-dimensional box models to infer the Bond albedo and redistribution efficiency. We elucidate the limitation of these box models by demonstrating that negative Bond albedos may be obtained due to a choice of boundary condition on the nightside temperature. We propose generalized definitions for the Bond albedo and heat redistribution efficiency for use with two-dimensional (2D) temperature maps. Open-source software called kelp is provided to efficiently compute these phase curves.

preprint2020arXiv

Atmospheric Characterization via Broadband Color Filters on the PLAnetary Transits and Oscillations of stars (PLATO) Mission

We assess broadband color filters for the two fast cameras on the PLAnetary Transits and Oscillations (PLATO) of stars space mission with respect to exoplanetary atmospheric characterization. We focus on Ultra Hot Jupiters and Hot Jupiters placed 25pc and 100pc away from the Earth and low mass low density planets placed 10pc and 25pc away. Our analysis takes as input literature values for the difference in transit depth between the broadband lower (500 to 675nm) wavelength interval (hereafter referred to as blue) and the upper (675-1125nm) broadband wavelength interval (hereafter referred to as red) for transmission, occultation and phase curve analyses. Planets orbiting main sequence central stars with stellar classes F, G, K and M are investigated. We calculate the signal-to-noise ratio with respect to photon and instrument noise for detecting the difference in transit depth between the two spectral intervals. Results suggest that bulk atmospheric composition and planetary geometric albedos could be detected for (Ultra) Hot Jupiters up to about 100pc (about 25pc) with strong (moderate) Rayleigh extinction. Phase curve information could be extracted for Ultra Hot Jupiters orbiting K and G dwarf stars up to 25pc away. For low mass low density planets, basic atmospheric types (primary and water-dominated) and the presence of sub-micron hazes in the upper atmosphere could be distinguished for up to a handful of cases up to about 10pc.

preprint2020arXiv

Exploring the atmospheric dynamics of the extreme ultra-hot Jupiter KELT-9b using TESS photometry

We carry out a phase-curve analysis of the KELT-9 system using photometric observations from NASA's Transiting Exoplanet Survey Satellite (TESS). The measured secondary eclipse depth and peak-to-peak atmospheric brightness modulation are $650^{+14}_{-15}$ ppm and $566\pm16$ ppm, respectively. The planet's brightness variation reaches maximum $31\pm5$ minutes before the midpoint of the secondary eclipse, indicating a $5\overset{\circ}{.}2\pm0\overset{\circ}{.}9$ eastward shift in the dayside hot spot from the substellar point. We also detect stellar pulsations on KELT-9 with a period of $7.58695\pm0.00091$ hours. The dayside emission of KELT-9b in the TESS bandpass is consistent with a blackbody brightness temperature of $4600\pm100$ K. The corresponding nightside brightness temperature is $3040\pm100$ K, comparable to the dayside temperatures of the hottest known exoplanets. In addition, we detect a significant phase-curve signal at the first harmonic of the orbital frequency and a marginal signal at the second harmonic. While the amplitude of the first harmonic component is consistent with the predicted ellipsoidal distortion modulation assuming equilibrium tides, the phase of this photometric variation is shifted relative to the expectation. Placing KELT-9b in the context of other exoplanets with phase-curve observations, we find that the elevated nightside temperature and relatively low day-night temperature contrast agree with the predictions of atmospheric models that include H$_{2}$ dissociation and recombination. The nightside temperature of KELT-9b implies an atmospheric composition containing about 50% molecular and 50% atomic hydrogen at 0.1 bar, a nightside emission spectrum that deviates significantly from a blackbody, and a 0.5-2.0 $μ$m transmission spectrum that is featureless at low resolution.

preprint2020arXiv

Helios-r2 -- A new Bayesian, open-source retrieval model for brown dwarfs and exoplanet atmospheres

We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2, which is applicable to medium-resolution emission spectra of brown dwarfs, in preparation for precision atmospheric spectroscopy in the era of the James Webb Space Telescope. The model is available as open-source code on the Exoclimes Simulation Platform. We subject Helios-r2 to a battery of tests of varying difficulty. The simplest test involves a mock retrieval on a forward model generated using the same radiative transfer technique, the same implementation of opacities, and the same chemistry model. The least trivial test involves a mock retrieval on synthetic spectra from the Sonora model grid, which uses a different radiative transfer technique, a different implementation of opacities, and a different chemistry model. A calibration factor, which is included to capture uncertainties in the brown dwarf radius, distance to the brown dwarf and flux calibration of the spectrum, may compensate, sometimes erroneously, for discrepancies in modeling choices and implementation. We analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown dwarf companions in the Epsilon Indi system. The retrieved surface gravities are consistent with previous studies and/or values inferred from dynamical masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on how to reject unphysical values of the retrieved brown dwarf radii. The inferred radii and corresponding masses should be taken with great caution. The retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical equilibrium is assumed.

preprint2020arXiv

Hunt for Starspots in HARPS Spectra of G and K Stars

We present a method for detecting starspots on cool stars using the cross-correlation function (CCF) of high resolution molecular spectral templates applied to archival high-resolution spectra of G and K stars observed with HARPS/HARPS-N. We report non-detections of starspots on the Sun even when the Sun was spotted, the solar twin 18 Scorpii, and the very spotted Sun-like star HAT-P-11, suggesting that Sun-like starspot distributions will be invisible to the CCF technique, and should not produce molecular absorption signals which might be confused for signatures of exoplanet atmospheres. We detect strong TiO absorption in the T Tauri K-dwarfs LkCa 4 and AA Tau, consistent with significant coverage by cool regions. We show that despite the non-detections, the technique is sensitive to relatively small spot coverages on M dwarfs and large starspot areas on Sun-like stars.

preprint2020arXiv

Information content of JWST-NIRSPEC transmission spectra of warm Neptunes

Warm Neptunes offer a rich opportunity for understanding exo-atmospheric chemistry. With the upcoming James Webb Space Telescope (JWST), there is a need to elucidate the balance between investments in telescope time versus scientific yield. We use the supervised machine learning method of the random forest to perform an information content analysis on a 11-parameter model of transmission spectra from the various NIRSpec modes. The three bluest medium-resolution NIRSpec modes (0.7 - 1.27 microns, 0.97 - 1.84 microns, 1.66 - 3.07 microns) are insensitive to the presence of CO. The reddest medium-resolution mode (2.87 - 5.10 microns) is sensitive to all of the molecules assumed in our model: CO, CO2, CH4, C2H2, H2O, HCN and NH3. It competes effectively with the three bluest modes on the information encoded on cloud abundance and particle size. It is also competitive with the low-resolution prism mode (0.6 - 5.3 microns) on the inference of every parameter except for the temperature and ammonia abundance. We recommend astronomers to use the reddest medium-resolution NIRSpec mode for studying the atmospheric chemistry of 800-1200 K warm Neptunes; its corresponding high-resolution counterpart offers diminishing returns. We compare our findings to previous JWST information content analyses that favor the blue orders, and suggest that the reliance on chemical equilibrium could lead to biased outcomes if this assumption does not apply. A simple, pressure-independent diagnostic for identifying chemical disequilibrium is proposed based on measuring the abundances of H2O, CO and CO2.

preprint2020arXiv

Interpreting High-Resolution Spectroscopy of Exoplanets Using Cross-Correlations and Supervised Machine Learning

We present a new method for performing atmospheric retrieval on ground-based, high-resolution data of exoplanets. Our method combines cross-correlation functions with a random forest, a supervised machine learning technique, to overcome challenges associated with high-resolution data. A series of cross-correlation functions are concatenated to give a "CCF-sequence" for each model atmosphere, which reduces the dimensionality by a factor of ~100. The random forest, trained on our grid of ~65,000 models, provides a likelihood-free method of retrieval. The pre-computed grid spans 31 values of both temperature and metallicity, and incorporates a realistic noise model. We apply our method to HARPS-N observations of the ultra-hot Jupiter KELT-9b, and obtain a metallicity consistent with solar (logM = $-0.2\pm0.2$). Our retrieved transit chord temperature (T = $6000^{+0}_{-200}$K) is unreliable as the ion cross-correlations lie outside of the training set, which we interpret as being indicative of missing physics in our atmospheric model. We compare our method to traditional nested-sampling, as well as other machine learning techniques, such as Bayesian neural networks. We demonstrate that the likelihood-free aspect of the random forest makes it more robust than nested-sampling to different error distributions, and that the Bayesian neural network we tested is unable to reproduce complex posteriors. We also address the claim in Cobb et al. (2019) that our random forest retrieval technique can be over-confident but incorrect. We show that this is an artefact of the training set, rather than the machine learning method, and that the posteriors agree with those obtained using nested-sampling.

preprint2019arXiv

Supervised Machine Learning for Intercomparison of Model Grids of Brown Dwarfs: Application to GJ 570D and the Epsilon Indi B Binary System

Self-consistent model grids of brown dwarfs involve complex physics and chemistry, and are often computed using proprietary computer codes, making it challenging to identify the reasons for discrepancies between model and data as well as between the models produced by different research groups. In the current study, we demonstrate a novel method for analyzing brown dwarf spectra, which combines the use of the Sonora, AMES-Cond and HELIOS model grids with the supervised machine learning method of the random forest. Besides performing atmospheric retrieval, the random forest enables information content analysis of the three model grids as a natural outcome of the method, both individually on each grid and by comparing the grids against one another, via computing large suites of mock retrievals. Our analysis reveals that the different choices made in modelling the alkali line shapes hinder the use of the alkali lines as gravity indicators. Nevertheless, the spectrum longward of 1.2 micron encodes enough information on the surface gravity to allow its inference from retrieval. Temperature may be accurately and precisely inferred independent of the choice of model grid, but not the surface gravity. We apply random forest retrieval to three objects: the benchmark T7.5 brown dwarf GJ 570D; and Epsilon Indi Ba (T1.5 brown dwarf) and Bb (T6 brown dwarf), which are part of a binary system and have measured dynamical masses. For GJ 570D, the inferred effective temperature and surface gravity are consistent with previous studies. For Epsilon Indi Ba and Bb, the inferred surface gravities are broadly consistent with the values informed by the dynamical masses.