ترغب بنشر مسار تعليمي؟ اضغط هنا

Powderday: Dust Radiative Transfer for Galaxy Simulations

70   0   0.0 ( 0 )
 نشر من قبل Desika Narayanan
 تاريخ النشر 2020
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.


قيم البحث

اقرأ أيضاً

74 - M. Juvela 2018
Thermal dust emission carries information on physical conditions and dust properties in many astronomical sources. Because observations represent a sum of emission along the line of sight, their interpretation often requires radiative transfer modell ing. We describe a new radiative transfer program SOC for computations of dust emission and examine its performance in simulations of interstellar clouds with external and internal heating. SOC implements the Monte Carlo radiative transfer method as a parallel program for shared memory computers. It can be used to study dust extinction, scattering, and emission. We tested SOC with realistic cloud models and examined the convergence and noise of the dust temperature estimates and of the resulting surface brightness maps. SOC has been demonstrated to produces accurate estimates for dust scattering and for thermal dust emission. It performs well with both with CPUs and with GPUs, the latter providing up to an order of magnitude speed-up. In the test cases, ALI improved the convergence rates but also was sensitive to Monte Carlo noise. Run-time refinement of the hierarchical-grid models did not help in reducing the run times required for a given accuracy of solution. The use of a reference field, without ALI, works more robustly. It also allows the run time to be optimised if the number of photon packages is increased only as the iterations progress. The use of GPUs in radiative transfer computations should be investigated further.
We present the first radiative transfer (RT) model of a non-edge-on disk galaxy in which the large-scale geometry of stars and dust is self-consistently derived through fitting of multiwavelength imaging observations from the UV to the submm. To this end we used the axi-symmetric RT model of Popescu et al. and a new methodology for deriving geometrical parameters, and applied this to decode the{spectral energy distribution (SED) of M33. We successfully account for both the spatial and spectral energy distribution, with residuals typically within $7%$ in the profiles of surface brightness and within $8%$ in the spatially-integrated SED. We predict well the energy balance between absorption and re-emission by dust, with no need to invoke modified grain properties, and we find no submm emission that is in excess of our model predictions. We calculate that $80pm8%$ of the dust heating is powered by the young stellar populations. We identify several morphological components in M33, a nuclear, an inner, a main and an outer disc, showing a monotonic trend in decreasing star-formation surface-density ($Sigma_{rm SFR}$) from the nuclear to the outer disc. In relation to surface density of stellar mass, the $Sigma_{rm SFR}$ of these components define a steeper relation than the main sequence of star-forming galaxies, which we call a structurally resolved main sequence. Either environmental or stellar feedback mechanisms could explain the slope of the newly defined sequence. We find the star-formation rate to be ${rm SFR}=0.28^{+0.02}_{-0.01}{rm M}_{odot}{rm yr}^{-1}$.
The radiative transport of photons through arbitrary three-dimensional (3D) structures of dust is a challenging problem due to the anisotropic scattering of dust grains and strong coupling between different spatial regions. The radiative transfer pro blem in 3D is solved using Monte Carlo or Ray Tracing techniques as no full analytic solution exists for the true 3D structures. We provide the first 3D dust radiative transfer benchmark composed of a slab of dust with uniform density externally illuminated by a star. This simple 3D benchmark is explicitly formulated to provide tests of the different components of the radiative transfer problem including dust absorption, scattering, and emission. This benchmark includes models with a range of dust optical depths fully probing cases that are optically thin at all wavelengths to optically thick at most wavelengths. This benchmark includes solutions for the full dust emission including single photon (stochastic) heating as well as two simplifying approximations: One where all grains are considered in equilibrium with the radiation field and one where the emission is from a single effective grain with size-distribution-averaged properties. A total of six Monte Carlo codes and one Ray Tracing code provide solutions to this benchmark. Comparison of the results revealed that the global SEDs are consistent on average to a few percent for all but the scattered stellar flux at very high optical depths. The image results are consistent within 10%, again except for the stellar scattered flux at very high optical depths. The lack of agreement between different codes of the scattered flux at high optical depths is quantified for the first time. We provide the first 3D dust radiative transfer benchmark and validate the accuracy of this benchmark through comparisons between multiple independent codes and detailed convergence tests.
We implement a state-of-the-art treatment of the processes affecting the production and Interstellar Medium (ISM) evolution of carbonaceous and silicate dust grains within SPH simulations. We trace the dust grain size distribution by means of a two-s ize approximation. We test our method on zoom-in simulations of four massive ($M_{200} geq 3 times 10^{14} M_{odot}$) galaxy clusters. We predict that during the early stages of assembly of the cluster at $z gtrsim 3$, where the star formation activity is at its maximum in our simulations, the proto-cluster regions are rich of dusty gas. Compared to the case in which only dust production in stellar ejecta is active, if we include processes occurring in the cold ISM,the dust content is enhanced by a factor $2-3$. However, the dust properties in this stage turn out to be significantly different than those observationally derived for the {it average} Milky Way dust, and commonly adopted in calculations of dust reprocessing. We show that these differences may have a strong impact on the predicted spectral energy distributions. At low redshift in star forming regions our model reproduces reasonably well the trend of dust abundances over metallicity as observed in local galaxies. However we under-produce by a factor of 2 to 3 the total dust content of clusters estimated observationally at low redshift, $z lesssim 0.5$ using IRAS, Planck and Herschel satellites data. This discrepancy can be solved by decreasing the efficiency of sputtering which erodes dust grains in the hot Intracluster Medium (ICM).
A large fraction of early-type galaxies (ETGs) host prominent dust features, and central dust rings are arguably the most interesting among them. We present here `Lord Of The Rings (LOTR), a new methodology which allows to integrate the extinction by dust rings in a 2D fitting modelling of the surface brightness distribution. Our pipeline acts in two steps, first using the surface fitting software GALFIT to determine the unabsorbed stellar emission, and then adopting the radiative transfer code SKIRT to apply dust extinction. We apply our technique to NGC 4552 and NGC 4494, two nearby ETGs. We show that the extinction by a dust ring can mimic, in a surface brightness profile, a central point source (e.g. an unresolved nuclear stellar cluster or an active galactic nucleus; AGN) superimposed to a `core (i.e. a central flattening of the stellar light commonly observed in massive ETGs). We discuss how properly accounting for dust features is of paramount importance to derive correct fluxes especially for low luminosity AGNs (LLAGNs). We suggest that the geometries of dust features are strictly connected with how relaxed is the gravitational potential, i.e. with the evolutionary stage of the host galaxy. Additionally, we find hints that the dust mass contained in the ring relates to the AGN activity.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا