AMRVAC and relativistic hydrodynamic simulations for gamma-ray burst afterglow phases

We apply a novel adaptive mesh refinement (AMR) code, AMRVAC (Adaptive Mesh Refinement version of the Versatile Advection Code), to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold interstellar medium (ISM). We do this for both 1D isotropic and full 2D jet-like fireball models. This is relevant for gamma-ray bursts (GRBs), and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell–ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor  γ= 100  up to the almost Newtonian     phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jet-like GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.     

Metal enrichment history of the proto-galactic interstellar medium

To investigate the metal enrichment history of the primordial interstellar medium (ISM), we have studied the long-term evolution of supernova remnants (SNRs) and how SNRs distribute the heavy metals into the ISM when they explode. With the assumed IMF for massive stars, we have computed the multiple supernova explosions and evolution in an inhomogeneous ISM. We compare the predicted metallicity distribution of metal deficient halo stars with the observed one. This revised version was published online in September 2006 with corrections to the Cover Date.     

Chemodynamical Modelling of Galaxy Formation and Evolution

We present our recently developed 3-dimensional chemodynamical code for galaxy evolution. This code follows the evolution of different galactic components like stars, dark matter and different components of the interstellar medium (ISM), i.e. a diffuse gaseous phase and the molecular clouds. Stars and dark matter are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse gas and a sticky particle scheme for the molecular clouds. Additionally, the galactic components are coupled by several phase transitions like star formation, stellar death or condensation and evaporation processes within the ISM. As an example we show the dynamical and chemical evolution of a star forming dwarf galaxy with a total baryonic mass of 2 ċ 10⁹ M_⊙. After a moderate collapse phase the stars and the molecular clouds follow an exponential radial distribution, whereas the diffuse gas shows a central depression as a result of stellar feedback. The metallicities of the galactic components behave quite differently with respect to their temporal evolution as well as their radial distribution. Especially, the ISM is at no stage well mixed. This revised version was published online in September 2006 with corrections to the Cover Date.     

Coupled convection and tidal dissipation in Europa’s ice shell

We performed 2D numerical simulations of oscillatory tidal flexing to study the interrelationship between tidal dissipation (calculated using the Maxwell model) and a heterogeneous temperature structure in Europa’s ice shell. Our 2D simulations show that, if the temperature is spatially uniform, the tidal dissipation rate peaks when the Maxwell time is close to the tidal period, consistent with previous studies. The tidal dissipation rate in a convective plume encased in a different background temperature depends on both the plume and background temperature. At a fixed background temperature, the dissipation increases strongly with plume temperature at low temperatures, peaks, and then decreases with temperature near the melting point when a melting-temperature viscosity of 10¹³ Pa s is used; however, the peak occurs at significantly higher temperature in this heterogeneous case than in a homogeneous medium for equivalent rheology. For constant plume temperature, the dissipation rate in a plume decreases as the surrounding temperature increases; plumes that are warmer than their surroundings can exhibit enhanced heating not only relative to their surroundings but relative to the Maxwell-model prediction for a homogeneous medium at the plume temperature. These results have important implications for thermal feedbacks in Europa’s ice shell.To self-consistently determine how convection interacts with tidal heating that is correctly calculated from the time-evolving heterogeneous temperature field, we coupled viscoelastic simulations of oscillatory tidal flexing (using Tekton) to long-term simulations of the convective evolution (using ConMan). Our simulations show that the tidal dissipation rate resulting from heterogeneous temperature can have a strong impact on thermal convection in Europa’s ice shell. Temperatures within upwelling plumes are greatly enhanced and can reach the melting temperature under plausible tidal-flexing amplitude for Europa. A pre-existing fracture zone (at least 6 km deep) promotes the concentration of tidal dissipation (up to ∼20 times more than that in the surroundings), leading to lithospheric thinning. This supports the idea that spatially variable tidal dissipation could lead locally to high temperatures, partial melting, and play an important role in the formation of ridges, chaos, or other features.     

The interaction of planetary nebulae and their asymptotic giant branch progenitors with the interstellar medium

Interaction with the interstellar medium (ISM) cannot be ignored in understanding planetary nebula (PN) evolution and shaping. In an effort to understand the range of shapes observed in the outer envelopes of PNe, we have run a comprehensive set of three-dimensional hydrodynamic simulations, from the beginning of the asymptotic giant branch (AGB) superwind phase until the end of the post-AGB/PN phase. A 'triple-wind' model is used, including a slow AGB wind, fast post-AGB wind and third wind reflecting the linear movement through the ISM. A wide range of stellar velocities, mass-loss rates and ISM densities have been considered.
We find that ISM interaction strongly affects outer PN structures, with the dominant shaping occurring during the AGB phase. The simulations predict four stages of PN–ISM interaction whereby (i) the PN is initially unaffected, (ii) then limb-brightened in the direction of motion, (iii) then distorted with the star moving away from the geometric centre, and (iv) finally so distorted that the object is no longer recognizable as a PN and may not be classed as such. Parsec-size shells around PNe are predicted to be common. The structure and brightness of ancient PNe are largely determined by the ISM interaction, caused by rebrightening during the second stage; this effect may address the current discrepancies in Galactic PN abundance. The majority of PNe will have tail structures. Evidence for strong interaction is found for all known PNe in globular clusters.     

Preliminary results of the multisite time-series campaign on the sdB pulsator PG 1325+101

I present a model for the formation and evolution of a massive disk galaxy, within a growing dark halo whose mass evolves according to cosmological simulations of structure formation. The galactic evolution is simulated with a new 3D chemo-dynamical code, including dark matter, stars and a multi-phase ISM. We follow the evolution from redshift z = 4.85 until the present epoch. The energy release by massive stars and supernovae prevents a rapid collapse of the baryonic matter and delays the maximum star formation until redshift z ≈ 1. The galaxy forms radially from inside-out and vertically from top-to-bottom. The feedback of stars leads to turbulent motions and large-scale flows in the ISM. As one result the galactic disk is significantly enriched by chemical elements synthesized in bulge stars.     

Mechanism of thermonuclear burning propagation in a helium layer on a neutron star surface: A simplified adiabatic model

Some thermonuclear X-ray bursters exhibit a high-frequency (about 300 Hz or more) brightness modulation at the rising phase of some bursts. These oscillations are explained by inhomogeneous heating of the surface layer on a rapidly rotating neutron star due to the finite propagation speed of thermonuclear burning. We suggest and substantiate a mechanism of this propagation that is consistent with experimental data. Initially, thermonuclear ignition occurs in a small region of the neutron star surface layer. The burning products rapidly rise and spread in the upper atmospheric layers due to turbulent convection. The accumulation of additional matter leads to matter compression and ignition at the bottom of the layer. This determines the propagation of the burning front. To substantiate this mechanism, we use the simplifying assumptions about a helium composition of the neutron star atmosphere and its initial adiabatic structure with a density of 1.75 × 10⁸ g cm⁻³ at the bottom. 2D numerical simulations have been performed using a modified particle method in the adiabatic approximation.     

Star formation in a multi-phase ISM

We present a 3d code for the dynamical evolution of a multi-phase interstellar medium (ISM) coupled to stars via star formation (SF) and feedback processes. The multi-phase ISM consists of clouds (sticky particles) and diffuse gas (SPH): exchange of matter, energy and momentum is achieved by drag (due to ram pressure) and condensation or evaporation processes. The cycle of matter is completed by SF and feedback by SNe and PNe. A SF scheme based on a variable SF efficiency as proposed by Elmegreen and Efremov (1997) is presented. For a Milky Way type galaxy we get a SF rate of ∼1 M_⊙ yr^-1 with an average SF efficiency of ∼5%. This revised version was published online in August 2006 with corrections to the Cover Date.     

The influence of gravitational acceleration on the supernova-driven Parker instability

Within a framework of 2D magnetohydrodynamic (MHD) simulations, we explore the dynamical regimes initiated by a supernova explosion in a magnetized stratified interstellar medium (ISM). We concentrate on the formation of large-scale magnetic structures and outflows connected with the Parker instability. For the sake of simplicity we only show models with a fixed explosion energy corresponding to a single supernova (SN) occurring in host galaxies with different fixed values of the gravitational acceleration g and different ratios of specific heats. We show that in general, depending on these two parameters, three different regimes are possible: a slowly growing Parker instability on time-scales much longer than the galactic rotation period for small g; the Parker instability growing at roughly the rotation period, which for ratios of specific heats larger than one is accompanied by an outflow resulting from the explosion for intermediate g; and a rapidly growing instability and a strong blow-out flow for large g . By means of numerical simulations and analytical estimates we show that the explosion energy and gravitational acceleration which separate the three regimes scale as Eg ²∼constant in the 2D case. We expect that in the 3D case this scaling law is Eg ³∼constant . Our simulations demonstrate furthermore that a single SN explosion can lead to the growth of multiple Parker loops in the disc and large-scale magnetic field loops in the halo, extending over 2–3 kpc horizontally and up to 3 kpc vertically above the mid-plane of the disc.     

G-dwarf metallicity distribution curves assuming an inhomogeneous interstellar medium

We show that the explicit assumption of a chemically inhomogeneous interstellar medium allows a better reproduction of the metallicity distribution of G-dwarfs in the solar neighbourhood. The inhomogeneity is considered by assuming that at any time stars are born with a spread in their metallicities, the spread being a Gaussian in the logarithm of the metallicity around the mean metallicity of that epoch. We show that for various simple models of chemical evolution, the fit to the G-dwarf metallicity curve improves considerably once the above assumption is applied. We show that the parameters obtained from the fitting also give acceptable predictions for the age-metallicity relation. We also find that if we use a G-dwarf metallicity function corrected for the scale height inflation of stars, the conventional models of chemical evolution cannot match the shape of the curve, at least under the instantaneous recycling approximation applied to a chemically homogeneous ISM. Under the inhomogeneous ISM approximation, the predicted shapes are found to be better, though not totally satisfactory.     

Gas phase processes affecting galactic evolution

Gas processes affecting star formation are reviewed with an emphasis on gravitational and magnetic instabilities as a source of turbulence. Gravitational instabilities are pervasive in a multi-phase medium, even for sub-threshold column densities, suggesting that only an ISM with a pure-warm phase can stop star formation. The instabilities generate turbulence, and this turbulence influences the structure and timing of star formation through its effect on the gas distribution and density. The final trigger for star formation is usually direct compression by another star or cluster. The star formation rate is apparently independent of the detailed mechanisms for star formation, and determined primarily by the total mass of gas in a dense form. If the density distribution function is a log-normal, as suggested by turbulence simulations, then this dense gas mass can be calculated and the star formation rate determined from first principles. The results suggest that only 10^-4 of the ISM mass actively participates in the star formation process and that this fraction does so because its density is larger than 10⁵ cm^-3, at which point several key processes affecting dynamical equilibrium begin to break down. This revised version was published online in August 2006 with corrections to the Cover Date.     

Evolution of Galaxies with Triaxial Haloes

We use particle simulations to study the motion of gas in galaxy models with mildly non-axisymmetric dark matter haloes with nearly constant density cores. In particular the effect of varying the dissipation rate is studied. We find that even very weak dissipation may cause inflow of material within the core radius towards the centre, and thus lead to the formation of a central mass concentration. Typically, a total of 10⁸ M _⊙ solar masses are accreted inside the central 100 pc in a few Gyr. This, in turn, destabilizes the trajectories in the central region. It is suggested that these processes may lead to the formation of bulge-like structures from discs, the extent of which will depend on the halo core radius and initial asymmetry. This and other possible consequences are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Coupled convection and tidal dissipation in Europa’s ice shell using non-Newtonian grain-size-sensitive (GSS) creep rheology

We present self-consistent, fully coupled two-dimensional (2D) numerical models of thermal evolution and tidal heating to investigate how convection interacts with tidal dissipation under the influence of non-Newtonian grain-size-sensitive creep rheology (plausibly resulting from grain boundary sliding) in Europa’s ice shell. To determine the thermal evolution, we solved the convection equations (using finite-element code ConMan) with the tidal dissipation as a heat source. For a given heterogeneous temperature field at a given time, we determined the tidal dissipation rate throughout the ice shell by solving for the tidal stresses and strains subject to Maxwell viscoelastic rheology (using finite-element code Tekton). In this way, the convection and tidal heating are fully coupled and evolve together. Our simulations show that the tidal dissipation rate can have a strong impact on the onset of thermal convection in Europa’s ice shell under non-Newtonian GSS rheology. By varying the ice grain size (1-10 mm), ice-shell thickness (20-120 km), and tidal-strain amplitude (0-4 × 10⁻⁵), we study the interrelationship of convection and conduction regimes in Europa’s ice shell. Under non-Newtonian grain-size-sensitive creep rheology and ice grain size larger than 1 mm, no thermal convection can initiate in Europa’s ice shell (for thicknesses <100 km) without tidal dissipation. However, thermal convection can start in thinner ice shells under the influence of tidal dissipation. The required tidal-strain amplitude for convection to occur decreases as the ice-shell thickness increases. For grain sizes of 1-10 mm, convection can occur in ice shells as thin as 20-40 km with the estimated tidal-strain amplitude of 2 × 10⁻⁵ on Europa.     

Statistical technique to study the connection between large‐scale structures in the ISM and low‐mass star formation

Applying a color index selection the Point Source Catalog of the Two Micron All Sky Survey (2MASS PSC) has been searched for Classical T Tauri (CTT) stars in the 2nd and 3rd Galactic quadrant based on their apparent K_S excess. The selection resulted in 3872 reliable CTT candidates. The obtained CTT sample is extended enough for statistical examination of the inhomogeneities in their distribution due to correlation with structures in the ISM, like infrared loops. A correlation was observed between the presence of dust loops and the CTT density. The latter shows an excess on loops with respect to that expected from random fluctuation in a homogeneous distribution matching with the observed overall distribution. Monte Carlo simulations were used to show the significance of the excess. The results imply that the formation of a fraction of CTTs was triggered during the loop formation. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Sulphur abundances in halo stars from multiplet 3 at 1045 nm

Sulphur is a volatile α ‐element which is not locked into dust grains in the interstellar medium (ISM). Hence, its abundance does not need to be corrected for dust depletion when comparing the ISM to the stellar atmospheres. The abundance of sulphur in the photosphere of metal‐poor stars is a matter of debate: according to some authors, [S/Fe] versus [Fe/H] forms a plateau at low metallicity, while, according to other studies, there is a large scatter or perhaps a bimodal distribution. In metal‐poor stars sulphur is detectable by its lines of multiplet 1 at 920 nm, but this range is heavily contaminated by telluric absorptions, and one line of the multiplet is blended by the hydrogen Paschen ζ line. We study the possibility of using multiplet 3 (at 1045 nm) for deriving the sulphur abundance because this range, now observable at the VLT with the infra‐red spectrograph CRIRES, is little contaminated by telluric absorption and not affected by blends at least in metal‐poor stars. We compare the abundances derived from multiplets 1 and 3, taking into account NLTE corrections and 3D effects. Here we present the results for a sample of four stars, although the scatter is less pronounced than in previous analysis, we cannot find a plateau in [S/Fe], and confirm the scatter of the sulphur abundance at low metallicity (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the importance of high-redshift intergalactic voids

We investigate the properties of 1D flux 'voids' (connected regions in the flux distribution above the mean-flux level) by comparing hydrodynamical simulations of large cosmological volumes with a set of observed high-resolution spectra at z ∼ 2. After addressing the effects of box size and resolution, we study how the void distribution changes when the most significant cosmological and astrophysical parameters are varied. We find that the void distribution in the flux is in excellent agreement with predictions of the standard Λcold dark matter (ΛCDM) cosmology, which also fits other flux statistics remarkably well. We then model the relation between flux voids and the corresponding 1D gas-density field along the line of sight and make a preliminary attempt to connect the 1D properties of the gas-density field to the 3D dark matter distribution at the same redshift. This provides a framework that allows statistical interpretations of the void population at high redshift using observed quasar spectra, and eventually it will enable linking the void properties of the high-redshift universe with those at lower redshifts, which are better known.     

Star Formation and Cosmological Simulations

The inclusion of a detailed modeling of the short-scale baryonic physics in a large-scale cosmological simulation is crucial for a better comparison between observations and predictions from cosmological models. From a set of 3D hydrodynamical simulations which include a chemical model to account for the complex physics of the ISM at a sub-grid scale, we have been able to obtain a statistically significant sample of galaxy-type halos with observational properties, like colors and luminosities for different cosmological scenarios. From this data base, we have studied a number of different things, like Tully-Fisher relations, luminosity functions and environmental effects. Despite the progress made during the last few years in the modeling of the physics of ISM and star formation, more work is clearly needed. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the thermal instability in numerical models of the interstellar medium

We investigate with 3D hydrodynamical simulations the role played by thermal processes in the dynamical evolution of the interstellarmedium (ISM). A parametric approach of the coolingprocess shows that the observed mass fraction of the cold (< 300 K)and unstable gas (300K < T < 6000K) can not be produced by turbulentcompression or background heating of the medium alone. An analysis of theproperties of the clouds that are formed by the combined effect of the thermal and gravitational instability shows that the cloud’s scaling relations imprinted by the thermal instability (TI) are in good agreementwith observational values.     

Approaching a homogeneous galaxy distribution: results from the Stromlo–APM redshift survey

Recent results from a number of redshift surveys suggest that the Universe is well described by an inhomogeneous, fractal distribution on the largest scales probed. This distribution has been found to have fractal dimension, D , approximately equal to 2.1, in contrast to a homogeneous distribution in which the dimension should approach the value 3 as the scale is increased. In this paper we demonstrate that estimates of D , based on the conditional density of galaxies, are prone to bias from several sources. These biases generally result in a smaller measured fractal dimension than the true dimension of the sample. We illustrate this behaviour in application to the Stromlo–APM redshift survey, showing that this data set in fact provides evidence for fractal dimension increasing with survey depth. On the largest scale probed, r ≈60  h ⁻¹ Mpc, we find evidence for a distribution with dimension D =2.76±0.10. A comparison between this sample and mock Stromlo–APM catalogues taken from N -body simulations (which assume a CDM cosmology) reveals a striking similarity in the behaviour of the fractal dimension. Thus we find no evidence for inhomogeneity in excess of that expected from conventional cosmological theory. We consider biases affecting future large surveys and demonstrate, using mock SDSS catalogues, that this survey will be able to measure the fractal dimension on scales at which we expect to see full turn-over to homogeneity, in an accurate and unbiased way.