Source author record

Jan-Pieter Paardekooper

Jan-Pieter Paardekooper appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

8 published item(s)

preprint2015arXiv

3D hydrodynamical and radiative transfer modeling of Eta Carinae's colliding winds

We present results of full 3D hydrodynamical and radiative transfer simulations of the colliding stellar winds in the massive binary system Eta Carinae. We accomplish this by applying the SimpleX algorithm for 3D radiative transfer on an unstructured Voronoi-Delaunay grid to recent 3D smoothed particle hydrodynamics (SPH) simulations of the binary colliding winds. We use SimpleX to obtain detailed ionization fractions of hydrogen and helium, in 3D, at the resolution of the original SPH simulations. We investigate several computational domain sizes and Luminous Blue Variable primary star mass-loss rates. We furthermore present new methods of visualizing and interacting with output from complex 3D numerical simulations, including 3D interactive graphics and 3D printing. While we initially focus on Eta Car, the methods employed can be applied to numerous other colliding wind (WR 140, WR 137, WR 19) and dusty 'pinwheel' (WR 104, WR 98a) binary systems. Coupled with 3D hydrodynamical simulations, SimpleX simulations have the potential to help determine the regions where various observed time-variable emission and absorption lines form in these unique objects.

preprint2015arXiv

3D Printing Meets Computational Astrophysics: Deciphering the Structure of Eta Carinae's Inner Colliding Winds

We present the first 3D prints of output from a supercomputer simulation of a complex astrophysical system, the colliding stellar winds in the massive (>120 M_Sun), highly eccentric (e ~ 0.9) binary star system Eta Carinae. We demonstrate the methodology used to incorporate 3D interactive figures into a PDF journal publication and the benefits of using 3D visualization and 3D printing as tools to analyze data from multidimensional numerical simulations. Using a consumer-grade 3D printer (MakerBot Replicator 2X), we successfully printed 3D smoothed particle hydrodynamics (SPH) simulations of Eta Carinae's inner (r ~ 110 au) wind-wind collision interface at multiple orbital phases. The 3D prints and visualizations reveal important, previously unknown 'finger-like' structures at orbital phases shortly after periastron (phi ~ 1.045) that protrude radially outward from the spiral wind-wind collision region. We speculate that these fingers are related to instabilities (e.g. thin-shell, Rayleigh-Taylor) that arise at the interface between the radiatively-cooled layer of dense post-shock primary-star wind and the fast (3000 km/s), adiabatic post-shock companion-star wind. The success of our work and easy identification of previously unrecognized physical features highlight the important role 3D printing and interactive graphics can play in the visualization and understanding of complex 3D time-dependent numerical simulations of astrophysical phenomena.

preprint2015arXiv

3D radiative transfer simulations of Eta Carinae's inner colliding winds - II: Ionization structure of helium at periastron

Spectral observations of the massive colliding wind binary Eta Carinae show phase-dependent variations, in intensity and velocity, of numerous helium emission and absorption lines throughout the entire 5.54-year orbit. Approaching periastron, the 3D structure of the wind-wind interaction region (WWIR) gets highly distorted due to the eccentric ($e \sim 0.9$) binary orbit. The secondary star ($η_{\mathrm{B}}$) at these phases is located deep within the primary's dense wind photosphere. The combination of these effects is thought to be the cause of the particularly interesting features observed in the helium lines at periastron. We perform 3D radiative transfer simulations of $η$ Car's interacting winds at periastron. Using the SimpleX radiative transfer algorithm, we post-process output from 3D smoothed particle hydrodynamic simulations of the inner 150 au of the $η$ Car system for two different primary star mass-loss rates ($\dot{M}_{η_{\mathrm{A}}}$). Using previous results from simulations at apastron as a guide for the initial conditions, we compute 3D helium ionization maps. We find that, for higher $\dot{M}_{η_{\mathrm{A}}}$, $η_{\mathrm{B}}$ He$^{0+}$-ionizing photons are not able to penetrate into the pre-shock primary wind. He$^{+}$ due to $η_{\mathrm{B}}$ is only present in a thin layer along the leading arm of the WWIR and in a small region close to the stars. Lowering $\dot{M}_{η_{\mathrm{A}}}$ allows $η_{\mathrm{B}}$'s ionizing photons to reach the expanding unshocked secondary wind on the apastron side of the system, and create a low fraction of He$^{+}$ in the pre-shock primary wind. With apastron on our side of the system, our results are qualitatively consistent with the observed variations in strength and radial velocity of $η$ Car's helium emission and absorption lines, which helps better constrain the regions where these lines arise.

preprint2015arXiv

The First Billion Years Project: The escape fraction of ionizing photons in the epoch of reionization

Proto-galaxies forming in low-mass dark matter haloes are thought to provide the majority of ionizing photons needed to reionize the Universe, due to their high escape fractions of ionizing photons. We study how the escape fraction in high-redshift galaxies relates to the physical properties of the halo in which the galaxies form, by computing escape fractions in more than 75000 haloes between redshifts 27 and 6 that were extracted from the First Billion Years project, high-resolution cosmological hydrodynamics simulations of galaxy formation. We find that the main constraint on the escape fraction is the gas column density in a radius of 10 pc around the stellar populations, causing a strong mass dependence of the escape fraction. The lower potential well in haloes with virial mass below 1e8 solar mass results in low column densities that can be penetrated by radiation from young stars (age < 5 Myr). In haloes with virial mass above 1e8 solar mass supernova feedback is important, but only 30% of the haloes in this mass range have an escape fraction higher than 1%. We find a large range of escape fractions in haloes with similar properties, caused by different distributions of the dense gas in the halo. This makes it very hard to predict the escape fraction on the basis of halo properties and results in a highly anisotropic escape fraction. The strong mass dependence, the large spread and the large anisotropy of the escape fraction may strongly affect the topology of reionization and is something current models of cosmic reionization should strive to take into account.

preprint2015arXiv

The impact of reionization on the formation of supermassive black hole seeds

Direct collapse black holes (DCBHs) formed from the collapse of atomically-cooled primordial gas in the early Universe are strong candidates for the seeds of supermassive BHs. DCBHs are thought to form in atomic cooling haloes in the presence of a strong molecule-dissociating, Lyman-Werner (LW) radiation field. Given that star forming galaxies are likely to be the source of the LW radiation in this scenario, ionizing radiation from these galaxies may accompany the LW radiation. We present cosmological simulations resolving the collapse of primordial gas into an atomic cooling halo, including the effects of both LW and ionizing radiation. We find that in cases where the gas is not self-shielded from the ionizing radiation, the collapse can be delayed by ~ 25 Myr. When the ionized gas does collapse, the free electrons that are present catalyze H2 formation. In turn, H2 cooling becomes efficient in the center of the halo, and DCBH formation is prevented. We emphasize, however, that in many cases the gas collapsing into atomic cooling haloes at high redshift is self-shielding to ionizing radiation. Therefore, it is only in a fraction of such haloes in which DCBH formation is prevented due to reionization.

preprint2014arXiv

3D radiative transfer simulations of Eta Carinae's inner colliding winds - I. Ionization structure of helium at apastron

The highly eccentric binary system Eta Carinae shows numerous time-variable emission and absorption features. These observational signatures are the result of interactions between the complex three-dimensional (3D) wind-wind collision regions and photoionization by the luminous stars. Specifically, helium presents several interesting spectral features that provide important clues on the geometry and physical properties of the system and the individual stars. We use the SimpleX algorithm to post-process 3D smoothed particle hydrodynamics simulation output of the interacting winds in Eta Car in order to obtain the fractions of ionized helium assuming three different primary star mass-loss rates. The resultant ionization maps constrain the regions where helium is singly- and doubly-ionized. We find that reducing the primary's mass-loss rate increases the volume of He+. Lowering the primary mass-loss rate produces large variations in the volume of He+ in the pre-shock primary wind on the periastron side of the system. Our results show that binary orientations in which apastron is on our side of the system are more consistent with available observations. We suggest that small variations in the primary's mass-loss rate might explain the observed increase in HeI absorption in recent decades, although numerous questions regarding this scenario remain open. We also propose that the absence of broad HeI lines in the spectra of Eta Car between its 1890's eruption and 1944 might be explained by the companion star's He0+ ionizing photons not being able to penetrate the wind-wind interaction region, due to a higher primary mass-loss rate at that time (by a factor >2, compared to the present value).

preprint2012arXiv

The First Billion Years project - IV: Proto-galaxies reionising the Universe

The contribution of stars in galaxies to cosmic reionisation depends on the star formation history in the Universe, the abundance of galaxies during reionisation, the escape fraction of ionising photons and the clumping factor of the inter-galactic medium (IGM). We compute the star formation rate and clumping factor during reionisation in a cosmological volume using a high-resolution hydrodynamical simulation. We post-process the output with detailed radiative transfer simulations to compute the escape fraction of ionising photons. Together, this gives us the opportunity to assess the contribution of galaxies to reionisation self-consistently. The strong mass and redshift dependence of the escape fraction indicates that reionisation occurred between z=15 and z=10 and was mainly driven by proto-galaxies forming in dark-matter haloes with masses between 1e7 and 1e8 solar mass. More massive galaxies that are rare at these redshifts and have significantly lower escape fractions contribute less photons to the reionisation process than the more-abundant low-mass galaxies. Star formation in the low-mass haloes is suppressed by radiative feedback from reionisation, therefore these proto-galaxies only contribute when the part of the Universe they live in is still neutral. After z~10, massive galaxies become more abundant and provide most of the ionising photons. In addition, we find that Population (Pop) III stars are too short-lived and not frequent enough to have a major contribution to reionisation. Although the stellar component of the proto-galaxies that produce the bulk of ionising photons during reionisation is too faint to be detected by the James Webb Space Telescope (JWST), these sources are brightest in the H-alpha and Ly-alpha recombination lines, which will likely be detected by JWST in deep surveys.

preprint2011arXiv

The escape of ionising radiation from high-redshift dwarf galaxies

The UV escape fraction from high-redshift galaxies plays a key role in models of cosmic reionisation. Because it is currently not possible to deduce the escape fractions during the epoch of reionisation from observations, we have to rely on numerical simulations. Our aim is to better constrain the escape fraction from high-redshift dwarf galaxies, as these are the most likely sources responsible for reionising the Universe. We employ a N-body/SPH method that includes realistic prescriptions for the physical processes that are important for the evolution of dwarf galaxies. These models are post-processed with radiative transfer to determine the escape fraction of ionising radiation. We perform a parameter study to assess the influence of the spin parameter, gas fraction and formation redshift of the galaxy and study the importance of numerical parameters as resolution, source distribution and local gas clearing. We find that the UV escape fraction from high-redshift dwarf galaxies that have formed a rotationally supported disc lie between 1e-5 and 0.1. The mass and angular momentum of the galaxy are the most important parameters that determine the escape fraction. We compare our results to previous work and discuss the uncertainties of our models. The low escape fraction we find for high-redshift dwarf galaxies is balanced by their high stellar content, resulting in an efficiency parameter for stars that is only marginally lower than the values found by semi-analytic models of reionisation. We therefore conclude that dwarf galaxies play an important role in cosmic reionisation also after the initial starburst phase, when the gas has settled into a disc.