The Use of Conditional Variational Autoencoders in Generating Stellar Spectra

This thesis studies the early phases of massive stars and their impact on the surrounding. The capabilities of continuum radiative transfer (RT) codes to interpret the observations are also investigated. The main results of this work are: 1) Two massive star-forming regions are observed in the infrared. The thermal emission from the ultra-compact HII regions is resolved and the spectral type of the ionizing stars is estimated. The hot cores are not detected thus implying line-of-sight extinction larger than 200 visual magnitude. 2) The first mid-infrared interferometric measurements towards a young massive star resolve thermal emission on scales of 30-50 AU probing the size of the predicted disk. The visibility curve differs from those of intermediate-mass stars. 3) The close vicinity of Theta1 C Ori are imaged using the NACO adaptive optics system. The binary proplyd Orion 168-326 and its interaction with the wind from Theta1 C Ori are resolved. A proplyd uniquely seen face-on is also identified. 4) Five RT codes are compared in a disk configuration. The solutions provide the first 2D benchmark and serve to test the reliability of other RT codes. The images/visibilities from two RT codes are compared for a distorted disk. The parameter range in which such a distortion is detectable with MIDI is explored.

We present an overview and several important upgrades to the Vector Discrete Ordinate Radiative Transfer (VDISORT) code. VDISORT is a polarized (vector) radiative transfer code that can be applied to a wide range of research problems including the Earth’s atmosphere and ocean system. First, a solution is developed to the complex algebraic eigenvalue problem resulting when the b2 component of the Stokes scattering matrix is non-zero. This solution is needed to compute the V component of the Stokes vector I=[I∥,I⊥,U,V]T. Second, a significant improvement in computational efficiency is obtained by reducing the dimension of the algebraic eigenvalue by a factor of 2 resulting in a speed increase of about 23 = 8. Third, an important upgrade of the VDISORT code is obtained by developing and implementing a method to enable output at arbitrary polar angles by the integration of the source function (ISF) method for partially reflecting Lambertian as well as general non-Lambertian surfaces. Fourth, a pseudo-spherical treatment has been implemented to provide important corrections for Earth curvature effects at near horizontal solar zenith and observation (viewing) polar angles. Fifth, a post-processing single-scattering correction procedure has been developed to enhance the accuracy and speed for strongly forward-peaked scattering. With these significant improvements the results from the upgraded version of the VDISORT code match published benchmark results for Rayleigh scattering, Mie scattering, and scattering by non-spherical cirrus particles. The performance of VDISORT for a polarized incident beam source is equally satisfactory. The VDISORT vector radiative transfer code is made public and freely available for use by the growing polarimetric research community including the space-borne polarimeters on the future NASA PACE and AOS missions.

Radiative interaction of atmosphere and surface: Write-up with elements of code

We study the formation of the Stokes profiles of the He I multiplet at 10 830 Å when relaxing two of the approximations that are typically considered in the modeling of this multiplet. Specifically, these are the lack of self-consistent radiation transfer and the assumption of equal illumination of the individual multiplet components. This He I multiplet is among the most important for the diagnostics of the outer solar atmosphere from spectropolarimetric observations, especially in prominences, filaments, and spicules. However, the aptness of these approximations is yet to be assessed, especially in situations where the optical thickness is on the order of one (or greater) and the radiation transfer has a significant impact in the local anisotropy as well as the ensuing spectral line polarization. This issue becomes particularly relevant in the ongoing development of new inversion tools that take into account multi-dimensional radiation transfer effects. To relax these approximations, we generalized the multi-term equations for the atomic statistical equilibrium to allow for a differential illumination of the multiplet components and implement them in a one-dimensional (1D) radiative transfer code. We find that even for this simple geometry and relatively small optical thickness, both radiation transfer and differential illumination effects have a significant impact on the emergent polarization profiles. These effects should be taken into account in order to avoid potentially significant errors when inferring the magnetic field vector.

In this study, Gxsview (Geometry and CROSS Section VIEWer), which is a stand-alone multi-platform integrated visualization tool for 3D geometry, 2D geometry, cross section, and input file, was developed. Its geometry engine covers almost all geometric features used in the major Monte Carlo radiation transport code, PHITS and MCNP, supporting nested universe, fill, lattice, like-but, and complement functions. The viewer is equipped with modern graphical user interface and an ability to pan, zoom, and rotate the geometry model in real time by using a mouse. The proposed tool is capable of exporting in 2D (png, jpg and xpm) or 3D (stl, vtk, vtp, and ply) formats, thereby realizing the collaboration with other software or 3D printers.

Context. Models of active galactic nuclei (AGN) suggest that their circumnuclear media are complex with clumps and filaments, while recent observations hint towards polar extended structures of gas and dust, as opposed to the classical torus paradigm. The X-ray band could form an interesting observational window to study these circumnuclear media in great detail. Aims. We want to extend the radiative transfer code SKIRT with the X-ray processes that govern the broad-band X-ray spectra of obscured AGN, to study the structure of AGN circumnuclear media in full 3D, based on their reflected X-ray emission. Methods. We extended the SKIRT code with Compton scattering on free electrons, photo-absorption and fluorescence by cold atomic gas, scattering on bound electrons, and extinction by dust. This includes a novel treatment of extreme-forward scattering by dust, and a detailed description of anomalous Rayleigh scattering. To verify our X-ray implementation, we performed the first dedicated benchmark of X-ray torus models, comparing five X-ray radiative transfer codes. Results. The resulting radiative transfer code covers the X-ray to millimetre wavelength range self-consistently, has all the features of the established SKIRT framework, is publicly available, and is fully optimised to operate in arbitrary 3D geometries. In the X-ray regime, we find an excellent agreement with the simulation results of the MYTORUS and REFLEX codes, which validates our X-ray implementation. We find some discrepancies with other codes, which illustrates the complexity of X-ray radiative transfer and motivates the need for a robust framework that can handle non-linear 3D radiative transfer effects. We illustrate the 3D nature of the code by producing synthetic X-ray images and spectra of clumpy torus models. Conclusions. SKIRT forms a powerful new tool to model circumnuclear media in full 3D, and make predictions for the X-ray band in addition to the dust-dominated infrared-to-UV wavelength range. The new X-ray functionalities of the SKIRT code allow for uncomplicated access to a broad suite of 3D X-ray models for AGN that can easily be tested and modified. This will be particularly useful with the advent of X-ray microcalorimeter observations in the near future.

Context. Radiative transfer modelling of expanding stellar envelopes is an important task in their analysis. To account for inhomogeneities and deviations from spherical symmetry, it is necessary to develop a 3 D approach to radiative transfer modelling. Aims. We present a 3 D Monte Carlo code for radiative transfer modelling, which is aimed to calculate the plasma ionisation and excitation state with the statistical equilibrium equations, moreover, to implement photon-matter coupling. As a first step, we present our Monte Carlo radiation transfer routines developed and tested from scratch. Methods. The background model atmosphere (the temperature, density, and velocity structure) can use an arbitrary grid referred to as the modGrid. The radiative transfer was solved using the Monte Carlo method in a Cartesian grid, referred to as the propGrid. This Cartesian grid was created based on the structure of the modGrid; correspondence between these two grids was set at the beginning of the calculations and then kept fixed. The propGrid can be either regular or adaptive; two modes of adaptive grids were tested. The accuracy and calculation speed for different propGrids was analysed. Photon interaction with matter was handled using the Lucy’s macroatom approach. Test calculations using our code were compared with the results obtained by a different Monte Carlo radiative transfer code. Results. Our method and the related code for the 3 D radiative transfer using the Monte Carlo and macroatom methods offer an accurate and reliable solution for the radiative transfer problem, and are especially promising for the inclusion and treatment of 3 D inhomogeneities.

SKIRT: The design of a suite of input models for Monte Carlo radiative transfer simulations

Extension of radiative transfer code MOMO, matrix-operator model to the thermal infrared – Clear air validation by comparison to RTTOV and application to CALIPSO-IIR

This and companion papers by Harrington et al. and Blecic et al. present the Bayesian Atmospheric Radiative Transfer (bart) code, an open-source, open-development package to characterize extrasolar planet atmospheres. bart combines a thermochemical equilibrium abundance (tea), a radiative transfer (Transit), and a Bayesian statistical (mc3) module to constrain atmospheric temperatures and molecular abundances for given spectroscopic observations. Here we describe the Transit radiative transfer package, an efficient line-by-line radiative transfer C code for one-dimensional atmospheres, developed by P. Rojo and further modified by the UCF exoplanet group. This code produces transmission and hemisphere-integrated emission spectra. Transit handles line-by-line opacities from HITRAN, Partridge & Schwenke (H2O), Schwenke (TiO), and Plez (VO) and collision-induced absorption from Borysow, HITRAN, and ExoMol. Transit emission spectra models agree with models from C. Morley (private communication) within a few percent. We applied bart to the Spitzer and Hubble transit observations of the Neptune-sized planet HAT-P-11b. Our analysis of the combined HST and Spitzer data generally agrees with those from previous studies, finding atmospheric models with enhanced metallicity (≳100× solar) and high-altitude clouds (≲1 mbar level). When analyzing only the HST data, our models favor high-metallicity atmospheres, in contrast with the previous analysis by Chachan et al. We suspect that this discrepancy arises from the different choice of chemistry modeling (free constant-with-altitude versus thermochemical equilibrium) and the enhanced parameter correlations found when neglecting the Spitzer observations. The bart source code and documentation are available at https://github.com/exosports/BART.

The radiative transfer equation (RTE) models the transport of light inside photonic scattering samples such as paint, foam and tissue. Analytic approximations to solve the RTE fail for samples with strong absorption and dominant anisotropic scattering and predict unphysical negative energy densities and the diffuse flux in the wrong direction. Here we fully characterize the unphysical regions of three popular approximations to the RTE for a slab, namely the P1 approximation (or diffusion approximation), the P3 approximation, and a popular modification to P3 that corrects the forward scattering in the approximation. We find that the delta function correction to P3 eliminates the unphysical range in the forward scattering. In addition, we compare the predictions of these analytical methods to exact Monte Carlo simulations for the physical and unphysical regions. We present maps of relative errors for the albedo and the anisotropy of the scatterers for a realistic index contrast typical of a polymer slab in air and optical thickness.

Radiative transfer analyses of Titan’s tropical atmosphere

Context. The analysis of spectral energy distributions (SEDs) of protoplanetary disks to determine their physical properties is known to be highly degenerate. Hence, a full Bayesian analysis is required to obtain parameter uncertainties and degeneracies. The main challenge here is computational speed, as one proper full radiative transfer model requires at least a couple of CPU minutes to compute. Aims. We performed a full Bayesian analysis for 30 well-known protoplanetary disks to determine their physical disk properties, including uncertainties and degeneracies. To circumvent the computational cost problem, we created neural networks (NNs) to emulate the SED generation process. Methods. We created two sets of MCFOST Monte Carlo radiative transfer disk models to train and test two NNs that predict SEDs for continuous and discontinuous disks, with 18 and 26 free model parameters, respectively. A Bayesian analysis was then performed on 30 protoplanetary disks with SED data collected by the FP7-Space DIANA project to determine the posterior distributions of all parameters. We ran this analysis twice, (i) with old distances and additional parameter constraints as used in a previous study, to compare results, and (ii) with updated distances and free choice of parameters to obtain homogeneous and unbiased model parameters. We evaluated the uncertainties in the determination of physical disk parameters from SED analysis, and detected and quantified the strongest degeneracies. Results. The NNs are able to predict SEDs within ~1 ms with uncertainties of about 5% compared to the true SEDs obtained by the radiative transfer code. We find parameter values and uncertainties that are significantly different from previous values obtained by χ2 fitting. Comparing the global evidence for continuous and discontinuous disks, we find that 26 out of 30 objects are better described by disks that have two distinct radial zones. The analysed sample shows a significant trend for massive disks to have small scale heights, which is consistent with lower midplane temperatures in massive disks. We find that the frequently used analytic relationship between disk dust mass and millimetre-flux systematically underestimates the dust mass for high-mass disks (dust mass ≥10−4 M⊙). We determine how well the dust mass can be determined with our method for different numbers of flux measurements. As a byproduct, we created an interactive graphical tool that instantly returns the SED predicted by our NNs for any parameter combination.

The authors have investigated the possibility of elaborating a new generation of radiative transfer models for climate studies based on the neural network technique. The authors show that their neural network–based model, NeuroFlux, can be used successfully for accurately deriving the longwave radiative budget from the top of the atmosphere to the surface. The reliable sampling of the earth’s atmospheric situations in the new version of the TIGR (Thermodynamic Initial Guess Retrieval) dataset, developed at the Laboratoire de Meteorologie Dynamique, allows for an efficient learning of the neural networks. Two radiative transfer models are applied to the computation of the radiative part of the dataset: a line-by-line model and a band model. These results have been used to infer the parameters of two neural network–based radiative transfer codes. Both of them achieve an accuracy comparable to, if not better than, the current general circulation model radiative transfer codes, and they are much fast...

Abstract. Aerosol optical depth retrieval over land surface using remote sensing employs the use of radiative transfer simulations and/or simultaneous measurements of atmospheric parameters at the time of satellite pass. Also, an accurate estimate of land surface parameters is also required in order to separate the atmospheric component from the land surface reflectance reaching at-sensor. In addition to empirical and semi-empirical approaches, amongst the most widely used methods to retrieve the aerosol properties from satellite measurements are radiative transfer codes used in either forward or inverse modes. As most of them are computationally complex, henceforth, efforts are made to formulate approximate models. In this study, we have tried to estimate aerosol optical depth using one such established physically based model, namely, SMART (Simple Model for Atmospheric Radiative Transfer) code in multiple scattering approximation for aerosols over first band (0.52–0.59 μm) of RESOURCESAT-AWiFS sensor. The aim of the analysis was to find out an approach to decouple aerosol effects from Top of atmosphere signals recorded by AWiFS sensor using multiple scattering approximations for aerosols. The model is first calibrated for aerosol asymmetry parameter for one dataset each of summer and winter seasons respectively and subsequently validated for 4 different datasets (2 summer and 2 winter) against the MODIS atmosphere product for aerosol optical depth. The results show that the difference between simulated vs. MODIS AOD fall within MODIS expected errors for the aerosol product.