@phdthesis{Feldmeier2001, author = {Feldmeier, Achim}, title = {Hydrodynamics of astrophysical winds driven by scattering in spectral lines}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-0000388}, school = {Universit{\"a}t Potsdam}, year = {2001}, abstract = {Liniengetriebene Winde werden durch Impuls{\"u}bertrag von Photonen auf ein Plasma bei Absorption oder Streuung in zahlreichen Spektrallinien beschleunigt. Dieser Prozess ist besonders effizient f{\"u}r ultraviolette Strahlung und Plasmatemperaturen zwischen 10^4 K und 10^5 K. Zu den astronomischen Objekten mit liniengetriebenen Winden geh{\"o}ren Sterne der Spektraltypen O, B und A, Wolf-Rayet-Sterne sowie Akkretionsscheiben verschiedenster Gr{\"o}ßenordnung, von Scheiben um junge Sterne und in kataklysmischen Ver{\"a}nderlichen bis zu Quasarscheiben. Es ist bislang nicht m{\"o}glich, das vollst{\"a}ndige Windproblem numerisch zu l{\"o}sen, also die Hydrodynamik, den Strahlungstransport und das statistische Gleichgewicht dieser Str{\"o}mungen gleichzeitig zu behandeln. Die Betonung liegt in dieser Arbeit auf der Windhydrodynamik, mit starken Vereinfachungen in den beiden anderen Gebieten. Wegen pers{\"o}nlicher Beteiligung betrachte ich drei Themen im Detail. 1. Windinstabilit{\"a}t durch Dopplerde-shadowing des Gases. Die Instabilit{\"a}t bewirkt, dass Windgas in dichte Schalen komprimiert wird, die von starken Stoßfronten begrenzt sind. Schnelle Wolken entstehen im Raum zwischen den Schalen und stoßen mit diesen zusammen. Dies erzeugt R{\"o}ntgenflashes, die die beobachtete R{\"o}ntgenstrahlung heißer Sterne erkl{\"a}ren k{\"o}nnen. 2. Wind runway durch radiative Wellen. Der runaway zeigt, warum beobachtete liniengetriebene Winde schnelle, kritische L{\"o}sungen anstelle von Brisenl{\"o}sungen (oder shallow solutions) annehmen. Unter bestimmten Bedingungen stabilisiert der Wind sich auf masse{\"u}berladenen L{\"o}sungen, mit einem breiten, abbremsenden Bereich und Knicken im Geschwindigkeitsfeld. 3. Magnetische Winde von Akkretionsscheiben um Sterne oder in aktiven Galaxienzentren. Die Linienbeschleunigung wird hier durch die Zentrifugalkraft entlang korotierender poloidaler Magnetfelder und die Lorentzkraft aufgrund von Gradienten im toroidalen Feld unterst{\"u}tzt. Ein Wirbelblatt, das am inneren Scheibenrand beginnt, kann zu stark erh{\"o}hten Massenverlustraten f{\"u}hren.}, language = {en} } @inproceedings{OPUS4-1574, title = {Clumping in hot-star winds : proceedings of an international workshop held in Potsdam, Germany, 18. - 22. June 2007}, editor = {Hamann, Wolf-Rainer and Feldmeier, Achim and Oskinova, Lida}, publisher = {Universit{\"a}tsverlag Potsdam}, address = {Potsdam}, isbn = {978-3-940793-33-1}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-13981}, pages = {254}, year = {2007}, abstract = {Stellar winds play an important role for the evolution of massive stars and their cosmic environment. Multiple lines of evidence, coming from spectroscopy, polarimetry, variability, stellar ejecta, and hydrodynamic modeling, suggest that stellar winds are non-stationary and inhomogeneous. This is referred to as 'wind clumping'. The urgent need to understand this phenomenon is boosted by its far-reaching implications. Most importantly, all techniques to derive empirical mass-loss rates are more or less corrupted by wind clumping. Consequently, mass-loss rates are extremely uncertain. Within their range of uncertainty, completely different scenarios for the evolution of massive stars are obtained. Settling these questions for Galactic OB, LBV and Wolf-Rayet stars is prerequisite to understanding stellar clusters and galaxies, or predicting the properties of first-generation stars. In order to develop a consistent picture and understanding of clumped stellar winds, an international workshop on 'Clumping in Hot Star Winds' was held in Potsdam, Germany, from 18. - 22. June 2007. About 60 participants, comprising almost all leading experts in the field, gathered for one week of extensive exchange and discussion. The Scientific Organizing Committee (SOC) included John Brown (Glasgow), Joseph Cassinelli (Madison), Paul Crowther (Sheffield), Alex Fullerton (Baltimore), Wolf-Rainer Hamann (Potsdam, chair), Anthony Moffat (Montreal), Stan Owocki (Newark), and Joachim Puls (Munich). These proceedings contain the invited and contributed talks presented at the workshop, and document the extensive discussions.}, language = {en} } @article{KustererNagelHartmannetal.2014, author = {Kusterer, D. -J. and Nagel, T. and Hartmann, S. and Werner, K. and Feldmeier, Achim}, title = {Monte Carlo radiation transfer in CV disk winds: application to the AM CVn prototype}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {561}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201321438}, pages = {10}, year = {2014}, abstract = {Context. AMCVn systems are ultracompact binaries in which a (semi-) degenerate star transfers helium-dominated matter onto a white dwarf. They are effective gravitational-wave emitters and potential progenitors of Type Ia supernovae. Aims. To understand the evolution of AMCVn systems it is necessary to determine their mass-loss rate through their radiation-driven accretion-disk wind. We constructed models to perform quantitative spectroscopy of P Cygni line profiles that were detected in UV spectra. Methods. We performed 2.5D Monte Carlo radiative transfer calculations in hydrodynamic wind structures by making use of realistic NLTE spectra from the accretion disk and by accounting for the white dwarf as an additional photon source. Results. We present first results from calculations in which LTE opacities are used in the wind model. A comparison with UV spectroscopy of the AMCVn prototype shows that the modeling procedure is potentially a good tool for determining mass-loss rates and abundances of trace metals in the helium-rich wind.}, language = {en} } @phdthesis{Hutter2014, author = {Hutter, Anne}, title = {Unveiling the epoch of reionization by simulations and high-redshift galaxies}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-76998}, school = {Universit{\"a}t Potsdam}, pages = {vi, 155}, year = {2014}, abstract = {The Epoch of Reionization marks after recombination the second major change in the ionization state of the universe, going from a neutral to an ionized state. It starts with the appearance of the first stars and galaxies; a fraction of high-energy photons emitted from galaxies permeate into the intergalactic medium (IGM) and gradually ionize the hydrogen, until the IGM is completely ionized at z~6 (Fan et al., 2006). While the progress of reionization is driven by galaxy evolution, it changes the ionization and thermal state of the IGM substantially and affects subsequent structure and galaxy formation by various feedback mechanisms. Understanding this interaction between reionization and galaxy formation is further impeded by a lack of understanding of the high-redshift galactic properties such as the dust distribution and the escape fraction of ionizing photons. Lyman Alpha Emitters (LAEs) represent a sample of high-redshift galaxies that are sensitive to all these galactic properties and the effects of reionization. In this thesis we aim to understand the progress of reionization by performing cosmological simulations, which allows us to investigate the limits of constraining reionization by high-redshift galaxies as LAEs, and examine how galactic properties and the ionization state of the IGM affect the visibility and observed quantities of LAEs and Lyman Break galaxies (LBGs). In the first part of this thesis we focus on performing radiative transfer calculations to simulate reionization. We have developed a mapping-sphere-scheme, which, starting from spherically averaged temperature and density fields, uses our 1D radiative transfer code and computes the effect of each source on the IGM temperature and ionization (HII, HeII, HeIII) profiles, which are subsequently mapped onto a grid. Furthermore we have updated the 3D Monte-Carlo radiative transfer pCRASH, enabling detailed reionization simulations which take individual source characteristics into account. In the second part of this thesis we perform a reionization simulation by post-processing a smoothed-particle hydrodynamical (SPH) simulation (GADGET-2) with 3D radiative transfer (pCRASH), where the ionizing sources are modelled according to the characteristics of the stellar populations in the hydrodynamical simulation. Following the ionization fractions of hydrogen (HI) and helium (HeII, HeIII), and temperature in our simulation, we find that reionization starts at z~11 and ends at z~6, and high density regions near sources are ionized earlier than low density regions far from sources. In the third part of this thesis we couple the cosmological SPH simulation and the radiative transfer simulations with a physically motivated, self-consistent model for LAEs, in order to understand the importance of the ionization state of the IGM, the escape fraction of ionizing photons from galaxies and dust in the interstellar medium (ISM) on the visibility of LAEs. Comparison of our models results with the LAE Lyman Alpha (Lya) and UV luminosity functions at z~6.6 reveals a three-dimensional degeneracy between the ionization state of the IGM, the ionizing photons escape fraction and the ISM dust distribution, which implies that LAEs act not only as tracers of reionization but also of the ionizing photon escape fraction and of the ISM dust distribution. This degeneracy does not even break down when we compare simulated with observed clustering of LAEs at z~6.6. However, our results show that reionization has the largest impact on the amplitude of the LAE angular correlation functions, and its imprints are clearly distinguishable from those of properties on galactic scales. These results show that reionization cannot be constrained tightly by exclusively using LAE observations. Further observational constraints, e.g. tomographies of the redshifted hydrogen 21cm line, are required. In addition we also use our LAE model to probe the question when a galaxy is visible as a LAE or a LBG. Within our model galaxies above a critical stellar mass can produce enough luminosity to be visible as a LBG and/or a LAE. By finding an increasing duty cycle of LBGs with Lya emission as the UV magnitude or stellar mass of the galaxy rises, our model reveals that the brightest (and most massive) LBGs most often show Lya emission. Predicting the Lya equivalent width (Lya EW) distribution and the fraction of LBGs showing Lya emission at z~6.6, we reproduce the observational trend of the Lya EWs with UV magnitude. However, the Lya EWs of the UV brightest LBGs exceed observations and can only be reconciled by accounting for an increased Lya attenuation of massive galaxies, which implies that the observed Lya brightest LAEs do not necessarily coincide with the UV brightest galaxies. We have analysed the dependencies of LAE observables on the properties of the galactic and intergalactic medium and the LAE-LBG connection, and this enhances our understanding of the nature of LAEs.}, language = {en} } @phdthesis{Kappel2015, author = {Kappel, David}, title = {Multi-spectrum retrieval of maps of Venus' surface emissivity in the infrared}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-85301}, school = {Universit{\"a}t Potsdam}, pages = {xix, 226}, year = {2015}, abstract = {The main goal of this cumulative thesis is the derivation of surface emissivity data in the infrared from radiance measurements of Venus. Since these data are diagnostic of the chemical composition and grain size of the surface material, they can help to improve knowledge of the planet's geology. Spectrally resolved images of nightside emissions in the range 1.0-5.1 μm were recently acquired by the InfraRed Mapping channel of the Visible and InfraRed Thermal Imaging Spectrometer (VIRTIS-M-IR) aboard ESA's Venus EXpress (VEX). Surface and deep atmospheric thermal emissions in this spectral range are strongly obscured by the extremely opaque atmosphere, but three narrow spectral windows at 1.02, 1.10, and 1.18 μm allow the sounding of the surface. Additional windows between 1.3 and 2.6 μm provide information on atmospheric parameters that is required to interpret the surface signals. Quantitative data on surface and atmosphere can be retrieved from the measured spectra by comparing them to simulated spectra. A numerical radiative transfer model is used in this work to simulate the observable radiation as a function of atmospheric, surface, and instrumental parameters. It is a line-by-line model taking into account thermal emissions by surface and atmosphere as well as absorption and multiple scattering by gases and clouds. The VIRTIS-M-IR measurements are first preprocessed to obtain an optimal data basis for the subsequent steps. In this process, a detailed detector responsivity analysis enables the optimization of the data consistency. The measurement data have a relatively low spectral information content, and different parameter vectors can describe the same measured spectrum equally well. A usual method to regularize the retrieval of the wanted parameters from a measured spectrum is to take into account a priori mean values and standard deviations of the parameters to be retrieved. This decreases the probability to obtain unreasonable parameter values. The multi-spectrum retrieval algorithm MSR is developed to additionally consider physically realistic spatial and temporal a priori correlations between retrieval parameters describing different measurements. Neglecting geologic activity, MSR also allows the retrieval of an emissivity map as a parameter vector that is common to several spectrally resolved images that cover the same surface target. Even applying MSR, it is difficult to obtain reliable emissivity maps in absolute values. A detailed retrieval error analysis based on synthetic spectra reveals that this is mainly due to interferences from parameters that cannot be derived from the spectra themselves, but that have to be set to assumed values to enable the radiative transfer simulations. The MSR retrieval of emissivity maps relative to a fixed emissivity is shown to effectively avoid most emissivity retrieval errors. Relative emissivity maps at 1.02, 1.10, and 1.18 μm are finally derived from many VIRTIS-M-IR measurements that cover a surface target at Themis Regio. They are interpreted as spatial variations relative to an assumed emissivity mean of the target. It is verified that the maps are largely independent of the choice of many interfering parameters as well as the utilized measurement data set. These are the first Venus IR emissivity data maps based on a consistent application of a full radiative transfer simulation and a retrieval algorithm that respects a priori information. The maps are sufficiently reliable for future geologic interpretations.}, language = {en} } @article{ThomasFeldmeier2016, author = {Thomas, Timon and Feldmeier, Achim}, title = {Radiative waves in stellar winds with line scattering}, series = {Monthly notices of the Royal Astronomical Society}, volume = {460}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stw1008}, pages = {1923 -- 1933}, year = {2016}, abstract = {Photospheric radiation can drive winds from hot, massive stars by direct momentum transfer through scattering in bound-bound transitions of atmospheric ions. The line radiation force should cause a new radiative wave mode. The dispersion relation from perturbations of the line force was analysed so far either in Sobolev approximation or for pure line absorption. The former does not include the line-driven instability, and the latter cannot account for upstream propagating, radiative waves. We consider a non-Sobolev line force that includes scattering in a simplified way, accounting however for the important line-drag effect. We derive a new dispersion relation for radiative waves, and analyse wave propagation using Fourier methods, and by numerical solution of an integro-differential equation. The existence of an upstream propagating, dispersive radiative wave mode is demonstrated.}, language = {en} } @article{MeyerKreplinKrausetal.2019, author = {Meyer, Dominique M.-A. and Kreplin, Alexander and Kraus, S. and Vorobyov, E. I. and Haemmerl{\´e}, Lionel and Eisl{\"o}ffel, Jochen}, title = {On the ALMA observability of nascent massive multiple systems formed by gravitational instability}, series = {Monthly notices of the Royal Astronomical Society}, volume = {487}, journal = {Monthly notices of the Royal Astronomical Society}, number = {4}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stz1585}, pages = {4473 -- 4491}, year = {2019}, abstract = {Massive young stellar objects (MYSOs) form during the collapse of high-mass pre-stellar cores, where infalling molecular material is accreted through a centrifugally balanced accretion disc that is subject to efficient gravitational instabilities. In the resulting fragmented accretion disc of the MYSO, gaseous clumps and low-mass stellar companions can form, which will influence the future evolution of massive protostars in the Hertzsprung-Russell diagram. We perform dust continuum radiative transfer calculations and compute synthetic images of disc structures modelled by the gravito-radiation-hydrodynamics simulation of a forming MYSO, in order to investigate the Atacama Large Millimeter/submillimeter Array (alma) observability of circumstellar gaseous clumps and forming multiple systems. Both spiral arms and gaseous clumps located at similar or equal to a few from the protostar can be resolved by interferometric alma Cycle 7 C43-8 and C43-10 observations at band 6 (), using a maximal 0.015 aracsec beam angular resolution and at least exposure time for sources at distances of . Our study shows that substructures are observable regardless of their viewing geometry or can be inferred in the case of an edge-viewed disc. The observation probability of the clumps increases with the gradually increasing efficiency of gravitational instability at work as the disc evolves. As a consequence, large discs around MYSOs close to the zero-age-main-sequence line exhibit more substructures than at the end of the gravitational collapse. Our results motivate further observational campaigns devoted to the close surroundings of the massive protostars S255IR-NIRS3 and NGC 6334I-MM1, whose recent outbursts are a probable signature of disc fragmentation and accretion variability.}, language = {en} } @article{ThomasPfrommer2019, author = {Thomas, T. and Pfrommer, Christoph}, title = {Cosmic-ray hydrodynamics}, series = {Monthly notices of the Royal Astronomical Society}, volume = {485}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stz263}, pages = {2977 -- 3008}, year = {2019}, abstract = {Star formation in galaxies appears to be self-regulated by energetic feedback processes. Among the most promising agents of feedback are cosmic rays (CRs), the relativistic ion population of interstellar and intergalactic plasmas. In these environments, energetic CRs are virtually collisionless and interact via collective phenomena mediated by kinetic-scale plasma waves and large-scale magnetic fields. The enormous separation of kinetic and global astrophysical scales requires a hydrodynamic description. Here, we develop a new macroscopic theory for CR transport in the self-confinement picture, which includes CR diffusion and streaming. The interaction between CRs and electromagnetic fields of Alfvenic turbulence provides the main source of CR scattering, and causes CRs to stream along the magnetic field with the Alfven velocity if resonant waves are sufficiently energetic. However, numerical simulations struggle to capture this effect with current transport formalisms and adopt regularization schemes to ensure numerical stability. We extent the theory by deriving an equation for the CRmomentum density along the mean magnetic field and include a transport equation for the Alfven-wave energy. We account for energy exchange of CRs and Alfven waves via the gyroresonant instability and include other wave damping mechanisms. Using numerical simulations, we demonstrate that our new theory enables stable, self-regulated CR transport. The theory is coupled to magnetohydrodynamics, conserves the total energy and momentum, and correctly recovers previous macroscopic CR transport formalisms in the steady-state flux limit. Because it is free of tunable parameters, it holds the promise to provide predictable simulations of CR feedback in galaxy formation.}, language = {en} } @article{KubatovaSzecsiSanderetal.2019, author = {Kubatova, Brankica and Szecsi, D. and Sander, Andreas Alexander Christoph and Kubat, Jiř{\´i} and Tramper, F. and Krticka, Jiri and Kehrig, C. and Hamann, Wolf-Rainer and Hainich, Rainer and Shenar, Tomer}, title = {Low-metallicity massive single stars with rotation}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {623}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201834360}, pages = {32}, year = {2019}, abstract = {Context. Metal-poor massive stars are assumed to be progenitors of certain supernovae, gamma-ray bursts, and compact object mergers that might contribute to the early epochs of the Universe with their strong ionizing radiation. However, this assumption remains mainly theoretical because individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the Small Magellanic Cloud. Aims. Here we explore the predictions of the state-of-the-art theories of stellar evolution combined with those of stellar atmospheres about a certain type of metal-poor (0.02 Z(circle dot)) hot massive stars, the chemically homogeneously evolving stars that we call Transparent Wind Ultraviolet INtense (TWUIN) stars. Methods. We computed synthetic spectra corresponding to a broad range in masses (20 130 M-circle dot) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage. We investigated the influence of mass loss and wind clumping on spectral appearance and classified the spectra according to the Morgan-Keenan (MK) system. Results. We find that TWUIN stars show almost no emission lines during most of their core hydrogen-burning lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission, and even some metal lines (oxygen and carbon, but still almost no nitrogen) are detected. Mass loss and clumping play a significant role in line formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early-O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early-O type star observed in a low-metallicity galaxy could be the result of chemically homogeneous evolution and might therefore be the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe because they are hot without showing prominent emission lines during most of their lifetime.}, language = {en} } @article{HainichRamachandranShenaretal.2019, author = {Hainich, Rainer and Ramachandran, Varsha and Shenar, Tomer and Sander, Andreas Alexander Christoph and Todt, Helge Tobias and Gruner, David and Oskinova, Lida and Hamann, Wolf-Rainer}, title = {PoWR grids of non-LTE model atmospheres for OB-type stars of various metallicities}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {621}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201833787}, pages = {12}, year = {2019}, abstract = {The study of massive stars in different metallicity environments is a central topic of current stellar research. The spectral analysis of massive stars requires adequate model atmospheres. The computation of such models is difficult and time-consuming. Therefore, spectral analyses are greatly facilitated if they can refer to existing grids of models. Here we provide grids of model atmospheres for OB-type stars at metallicities corresponding to the Small and Large Magellanic Clouds, as well as to solar metallicity. In total, the grids comprise 785 individual models. The models were calculated using the state-of-the-art Potsdam Wolf-Rayet (PoWR) model atmosphere code. The parameter domain of the grids was set up using stellar evolution tracks. For all these models, we provide normalized and flux-calibrated spectra, spectral energy distributions, feedback parameters such as ionizing photons, Zanstra temperatures, and photometric magnitudes. The atmospheric structures (the density and temperature stratification) are available as well. All these data are publicly accessible through the PoWR website.}, language = {en} } @article{Meyer2021, author = {Meyer, Dominique M.-A.}, title = {On the bipolarity of Wolf-Rayet nebulae}, series = {Monthly notices of the Royal Astronomical Society}, volume = {507}, journal = {Monthly notices of the Royal Astronomical Society}, number = {4}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stab2426}, pages = {4697 -- 4714}, year = {2021}, abstract = {Wolf-Rayet stars are amongst the rarest but also most intriguing massive stars. Their extreme stellar winds induce famous multiwavelength circumstellar gas nebulae of various morphologies, spanning from circles and rings to bipolar shapes. This study is devoted to the investigation of the formation of young, asymmetric Wolf-Rayet gas nebulae and we present a 2.5-dimensional magneto-hydrodynamical toy model for the simulation of Wolf-Rayet gas nebulae generated by wind-wind interaction. Our method accounts for stellar wind asymmetries, rotation, magnetization, evolution, and mixing of materials. It is found that the morphology of the Wolf-Rayet nebulae of blue supergiant ancestors is tightly related to the wind geometry and to the stellar phase transition time interval, generating either a broadened peanut-like or a collimated jet-like gas nebula. Radiative transfer calculations of our Wolf-Rayet nebulae for dust infrared emission at 24 mu m show that the projected diffuse emission can appear as oblate, bipolar, ellipsoidal, or ring structures. Important projection effects are at work in shaping observed Wolf-Rayet nebulae. This might call a revision of the various classifications of Wolf-Rayet shells, which are mostly based on their observed shape. Particularly, our models question the possibility of producing pre-Wolf-Rayet wind asymmetries, responsible for bipolar nebulae like NGC 6888, within the single red supergiant evolution channel scenario. We propose that bipolar Wolf-Rayet nebulae can only be formed within the red supergiant scenario by multiple/merged massive stellar systems, or by single high-mass stars undergoing additional, e.g. blue supergiant, evolutionary stages prior to the Wolf-Rayet phase.}, language = {en} } @phdthesis{Thomas2022, author = {Thomas, Timon}, title = {Cosmic-ray hydrodynamics: theory, numerics, applications}, doi = {10.25932/publishup-56384}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-563843}, school = {Universit{\"a}t Potsdam}, pages = {334}, year = {2022}, abstract = {Cosmic rays (CRs) are a ubiquitous and an important component of astrophysical environments such as the interstellar medium (ISM) and intracluster medium (ICM). Their plasma physical interactions with electromagnetic fields strongly influence their transport properties. Effective models which incorporate the microphysics of CR transport are needed to study the effects of CRs on their surrounding macrophysical media. Developing such models is challenging because of the conceptional, length-scale, and time-scale separation between the microscales of plasma physics and the macroscales of the environment. Hydrodynamical theories of CR transport achieve this by capturing the evolution of CR population in terms of statistical moments. In the well-established one-moment hydrodynamical model for CR transport, the dynamics of the entire CR population are described by a single statistical quantity such as the commonly used CR energy density. In this work, I develop a new hydrodynamical two-moment theory for CR transport that expands the well-established hydrodynamical model by including the CR energy flux as a second independent hydrodynamical quantity. I detail how this model accounts for the interaction between CRs and gyroresonant Alfv{\´e}n waves. The small-scale magnetic fields associated with these Alfv{\´e}n waves scatter CRs which fundamentally alters CR transport along large-scale magnetic field lines. This leads to the effects of CR streaming and diffusion which are both captured within the presented hydrodynamical theory. I use an Eddington-like approximation to close the hydrodynamical equations and investigate the accuracy of this closure-relation by comparing it to high-order approximations of CR transport. In addition, I develop a finite-volume scheme for the new hydrodynamical model and adapt it to the moving-mesh code Arepo. This scheme is applied using a simulation of a CR-driven galactic wind. I investigate how CRs launch the wind and perform a statistical analysis of CR transport properties inside the simulated circumgalactic medium (CGM). I show that the new hydrodynamical model can be used to explain the morphological appearance of a particular type of radio filamentary structures found inside the central molecular zone (CMZ). I argue that these harp-like features are synchrotron-radiating CRs which are injected into braided magnetic field lines by a point-like source such as a stellar wind of a massive star or a pulsar. Lastly, I present the finite-volume code Blinc that uses adaptive mesh refinement (AMR) techniques to perform simulations of radiation and magnetohydrodynamics (MHD). The mesh of Blinc is block-structured and represented in computer memory using a graph-based approach. I describe the implementation of the mesh graph and how a diffusion process is employed to achieve load balancing in parallel computing environments. Various test problems are used to verify the accuracy and robustness of the employed numerical algorithms.}, language = {en} }