Aims: We use new deep XMM-Newton observations of the northern disk of M 31 to trace the hot interstellar medium (ISM) in unprecedented detail and to characterise the physical properties of the X-ray emitting plasmas. Methods: We used all XMM-Newton data up to and including our new observations to produce the most detailed image yet of the hot ISM plasma in a grand design spiral galaxy such as our own. We compared the X-ray morphology to multi-wavelength studies in the literature to set it in the context of the multi-phase ISM. We performed spectral analyses on the extended emission using our new observations as they offer sufficient depth and count statistics to constrain the plasma properties. Data from the Panchromatic Hubble Andromeda Treasury were used to estimate the energy injected by massive stars and their supernovae. We compared these results to the hot gas properties. Results: The brightest emission regions were found to be correlated with populations of massive stars, notably in the 10 kpc star-forming ring. The plasma temperatures in the ring regions are ~0.2 keV up to ~0.6 keV. We suggest this emission is hot ISM heated in massive stellar clusters and superbubbles. We derived X-ray luminosities, densities, and pressures for the gas in each region. We also found large extended emission filling low density gaps in the dust morphology of the northern disk, notably between the 5 kpc and 10 kpc star-forming rings. We propose that the hot gas was heated and expelled into the gaps by the populations of massive stars in the rings. Conclusions. It is clear that the massive stellar populations are responsible for heating the ISM to X-ray emitting temperatures, filling their surroundings, and possibly driving the hot gas into the low density regions. Overall, the morphology and spectra of the hot gas in the northern disk of M 31 is similar to other galaxy disks.

Aims. Determination of K- and L-shell cross sections of the carbon atom and ions using the modified relativistic binary encounter Bethe (MRBEB) method, a simple analytical scheme based on one atomic parameter that allows determining electron-impact ionization cross sections. The quality of the cross sections calculated with the MRBEB method is shown through: (i) comparison with those obtained with the general ionization processes in the presence of electrons and radiation (GIPPER) code and the flexible atomic code (FAC), and (ii) determination of their effects on the ionic structure and cooling of an optically thin plasma. Results. The three sets of cross sections show deviations among each other in different energy regions. The largest deviations occur near and in the peak maximum. Ion fractions and plasma emissivities of an optically thin plasma that evolves under collisional ionization equilibrium, derived using each set of cross sections, show deviations that decrease with increase in temperature and ionization degree. In spite of these differences, the calculations using the three sets of cross sections agree overall. Conclusions. A simple model like the MRBEB is capable of providing cross sections similar to those calculated with more sophisticated quantum mechanical methods in the GIPPER and FAC codes.

We carried out new observations of two fields in the northern ring of M31 with XMM-Newton with two exposures of 100 ks each and obtained a complete list of X-ray sources down to a sensitivity limit of ~7 x 10^34 erg s^-1 (0.5 - 2.0 keV). The major objective of the observing programme was the study of the hot phase of the ISM in M31. The analysis of the diffuse emission and the study of the ISM is presented in a separate paper. We analysed the spectral properties of all detected sources using hardness ratios and spectra if the statistics were high enough. We also checked for variability. We cross-correlated the source list with the source catalogue of a new survey of the northern disc of M31 carried out with Chandra and Hubble (Panchromatic Hubble Andromeda Treasury, PHAT) as well as with other existing catalogues. We detected a total of 389 sources, including 43 foreground stars and candidates and 50 background sources. Based on the comparison to the Chandra/PHAT survey, we classify 24 hard X-ray sources as new candidates for X-ray binaries (XRBs). In total, we identified 34 XRBs and candidates and 18 supernova remnants (SNRs) and candidates. Three of the four brightest SNRs show emission mainly below 2 keV, consistent with shocked ISM. The spectra of two of them also require an additional component with a higher temperature. The SNR [SPH11] 1535 has a harder spectrum and might suggest that there is a pulsar-wind nebula inside the SNR. We find five new sources showing clear time variability. We also studied the spectral properties of the transient source SWIFT J004420.1+413702, which shows significant variation in flux over a period of seven months (June 2015 to January 2016) and associated change in absorption. Based on the likely optical counterpart detected in the Chandra/PHAT survey, the source is classified as a low-mass X-ray binary.

{\bf Context.} Numerical models of the evolution of interstellar and integalactic plasmas often assume that the adiabatic parameter $\gamma$ (the ratio of the specific heats) is constant (5/3 in monoatomic plasmas). However, $\gamma$ is determined by the total internal energy of the plasma, which depends on the ionic and excitation state of the plasma. Hence, the adiabatic parameter may not be constant across the range of temperatures available in the interstellar medium. {\bf Aims.} We aim to carry out detailed simulations of the thermal evolution of plasmas with Maxwell-Boltzmann and non-thermal ($\kappa$ and $n$) electron distributions in order to determine the temperature variability of the total internal energy and of the adiabatic parameter. {\bf Methods.} The plasma, composed of H, He, C, N, O, Ne, Mg, Si, S, and Fe atoms and ions, evolves under collisional ionization equilibrium conditions, from an initial temperature of $10^9$ K. The calculations include electron impact ionization, radiative and dielectronic recombinations and line excitation. The ionization structure was calculated (...) using the Gauss elimination method with scaled partial pivoting. Numerical integrations (...) were carried out using the double-exponential over a semi-finite interval method. In both methods a precision of $10^{-15}$ is adopted. {\bf Results.} The total internal energy of the plasma is mainly dominated by the ionization energy for temperatures lower than $8\times 10^4$ K with the excitation energy having a contribution of less than one percent. In thermal and non-thermal plasmas composed of H, He, and metals, the adiabatic parameter evolution is determined by the H and He ionizations leading to a profile in general having three transitions. However, for $\kappa$ distributed plasmas these three transitions are not observed for $\kappa<15$ (...). In general, $\gamma$ varies from 1.01 to 5/3.

Tracking the thermal evolution of plasmas, characterized by an n-distribution, using numerical simulations, requires the determination of the emission spectra and of the radiative losses due to free-free emission from the correspond- ing temperature averaged and total Gaunt factors. Detailed calculations of the latter are presented, associated to n-distributed electrons with the parameter n ranging from 1 (corresponding to the Maxwell-Boltzmann distribu- tion) to 100. The temperature averaged and total Gaunt factors, with decreasing n tend to those obtained with the Maxwell-Boltzmann distribution. Radiative losses due to free-free emission in a plasma evolving under collisional ionization equilibrium conditions and composed by H, He, C, N, O, Ne, Mg, Si, S, and Fe ions, are presented. These losses decrease with the decrease in the parameter n reaching a minimum when n = 1, and, thus converging to the losses of a thermal plasma. Tables of the thermal averaged and total Gaunt factors calculated for n distributions and a wide range electron and photon energies are presented.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

Aims. Optically thin plasmas may deviate from thermal equilibrium and thus, electrons (and ions) are no longer described by the Maxwellian distribution. Instead they can be described by $\kappa$-distributions. The free-free spectrum and radiative losses depend on the temperature-averaged (over the electrons distribution) and total Gaunt factors, respectively. Thus, there is a need to calculate and make available these factors to be used by any software that deals with plasma emission. Methods. We recalculated the free-free Gaunt factor for a wide range of energies and frequencies using hypergeometric functions of complex arguments and the Clenshaw recurrence formula technique combined with approximations whenever the difference between the initial and final electron energies is smaller than $10^{-10}$ in units of $z^2Ry$. We used double and quadruple precisions. The temperature- averaged and total Gaunt factors calculations make use of the Gauss-Laguerre integration with 128 nodes. Results. The temperature-averaged and total Gaunt factors depend on the $\kappa$ parameter, which shows increasing deviations (with respect to the results obtained with the use of the Maxwellian distribution) with decreasing $\kappa$. Tables of these Gaunt factors are provided.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

Until recently the dynamical evolution of the interstellar medium (ISM) was simulated using collisional ionization equilibrium (CIE) conditions. However, the ISM is a dynamical system, in which the plasma is naturally driven out of equilibrium due to atomic and dynamic processes operating on different timescales. A step forward in the field comprises a multi-fluid approach taking into account the joint thermal and dynamical evolutions of the ISM gas.

Using three-dimensional non-equilibrium ionization (NEI) hydrodynamical simulation of the interstellar medium (ISM), we study the electron density, $n_{e}$, in the Galactic disk and compare it with the values derived from dispersion measures towards pulsars with known distances located up to 200 pc on either side of the Galactic midplane. The simulation results, consistent with observations, can be summarized as follows: (i) the DMs in the simulated disk lie between the maximum and minimum observed values, (ii) the log <n_e> derived from lines of sight crossing the simulated disk follows a Gaussian distribution centered at \mu=-1.4 with a dispersion \sigma=0.21, thus, the Galactic midplane <n_e>=0.04\pm 0.01$ cm$^{-3}$, (iii) the highest electron concentration by mass (up to 80%) is in the thermally unstable regime (200<T<10^{3.9} K), (iv) the volume occupation fraction of the warm ionized medium is 4.9-6%, and (v) the electrons have a clumpy distribution along the lines of sight.

Aims. We present the first high-resolution non-equilibrium ionization simulation of the joint evolution of the Local Bubble (LB) and Loop I superbubbles in the turbulent supernova-driven interstellar medium (ISM). The time variation and spatial distribution of the Li-like ions Civ, Nv, and Ovi inside the LB are studied in detail. Methods. This work uses the parallel adaptive mesh refinement code EAF-PAMR coupled to the newly developed atomic and molecular plasma emission module E(A+M)PEC, featuring the time-dependent calculation of the ionization structure of H through Fe, using the latest revision of solar abundances. The finest AMR resolution is 1 pc within a grid that covers a representative patch of the Galactic disk (with an area of 1 kpc^2 in the midplane) and halo (extending up to 10 kpc above and below the midplane). Results. The evolution age of the LB is derived by the match between the simulated and observed absorption features of the Li-like ions Civ, Nv, and Ovi . The modeled LB current evolution time is bracketed between 0.5 and 0.8 Myr since the last supernova reheated the cavity in order to have N(Ovi) < 8 \times 10^12 cm-2, log[N(Civ) /N(Ovi) ] < -0.9 and log[N(Nv) /N(Ovi) ] < -1 inside the simulated LB cavity, as found in Copernicus, IUE, GHRS-IST and FUSE observations.

E(A+M)PEC traces the ionization structure, cooling and emission spectra of plasmas. It is written in OpenCL, runs in NVIDIA Graphics Processor Units and can be coupled to any HD or MHD code to follow the dynamical and thermal evolution of any plasma in, e.g., the interstellar medium (ISM).

We present the first to date three-dimensional high-resolution hydrodynamical simulation tracing the non-equilibrium ionization evolution (using the Eborae Atomic and Molecular Plasma Emission Code - E(A+M)PEC) of the Local Bubble and Loop I superbubbles embedded in a turbulent supernova-driven interstellar medium.

The ISM, powered by SNe, is turbulent and permeated by a magnetic field (with a mean and a turbulent component). It constitutes a frothy medium that is mostly out of equilibrium and is ram pressure dominated on most of the temperature ranges, except for T< 200 K and T> 1E6 K, where magnetic and thermal pressures dominate, respectively. Such lack of equilibrium is also imposed by the feedback of the radiative processes into the ISM flow. Many models of the ISM or isolated phenomena, such as bubbles, superbubbles, clouds evolution, etc., take for granted that the flow is in the so-called collisional ionization equilibrium (CIE). However, recombination time scales of most of the ions below 1E6 K are longer than the cooling time scale. This implies that the recombination lags behind and the plasma is overionized while it cools. As a consequence cooling deviates from CIE. This has severe implications on the evolution of the ISM flow and its ionization structure. Here, besides reviewing several models of the ISM, including bubbles and superbubbles, the validity of the CIE approximation is discussed, and a presentation of recent developments in modeling the ISM by taking into account the time-dependent ionization structure of the flow in a full-blown numerical 3D high resolution simulation is presented.

High-resolution non-ideal magnetohydrodynamical simulations of the turbulent magnetized ISM, powered by supernovae types Ia and II at Galactic rate, including self-gravity and non-equilibriuim ionization (NEI), taking into account the time evolution of the ionization structure of H, He, C, N, O, Ne, Mg, Si, S and Fe, were carried out. These runs cover a wide range (from kpc to sub-parsec) of scales, providing resolution independent information on the injection scale, extended self-similarity and the fractal dmension of the most dissipative structures.

We review the evolution of the interstellar medium in disc galaxies, and show, both analytically and by numerical 3D hydrodynamic simulations, that the disc-halo connection is an essential ingredient in understanding the evolution of star forming galaxies. Depending on the star formation rate of the underlying gaseous disc, a galactic fountain is established. If the star formation rate is sufficiently high and/or cosmic rays are well coupled to the thermal plasma, a galactic wind will be formed and lead to a secular mass loss of the galaxy. Such a wind leaves a unique imprint on the soft X-ray spectra in edge-on galaxies, with delayed recombination being one of its distinctive features. We argue that synthetic spectra, obtained from self-consistent dynamical and thermal modelling of a galactic outflow, should be treated on an equal footing as observed spectra. We show that it is thus possible to successfully fit the spectrum of the starburst galaxy NGC 3079.

In this paper we review the modeling of the Local Bubble (LB) with special emphasis on the progress we have made since the last major conference "The Local Bubble and Beyond (I)" held in Garching in 1997. Since then new insight was gained into the possible origin of the LB, with a moving group crossing its volume during the last 10 - 15 Myr being most likely responsible for creating a local cavity filled with hot recombining gas. Numerical high resolution 3D simulations of a supernova driven inhomogeneous interstellar medium show that we can reproduce both the extension of the LB and the OVI column density in absorption measured with FUSE for a LB age of 13.5 - 14.5 Myr. We further demonstrate that the LB evolves like an ordinary superbubble expanding into a density stratified medium by comparing analytical 2D Kompaneets solutions to NaI contours, representing the extension of the local cavity. These results suggest that LB blow-out into the Milky Way halo has occurred roughly 5 Myr ago.

We present an overview of the Chandra ACIS Survey of M33 (ChASeM33): A Deep Survey of the Nearest Face-on Spiral Galaxy. The 1.4 Ms survey covers the galaxy out to $R \approx 18\arcmin (\approx 4$ kpc). These data provide the most intensive, high spatial resolution assessment of the X-ray source populations available for the confused inner regions of M33. Mosaic images of the ChASeM33 observations show several hundred individual X-ray sources as well as soft diffuse emission from the hot interstellar medium. Bright, extended emission surrounds the nucleus and is also seen from the giant \hii regions NGC 604 and IC 131. Fainter extended emission and numerous individual sources appear to trace the inner spiral structure. The initial source catalog, arising from $\sim$~2/3 of the expected survey data, includes 394 sources significant at the $3\sigma$ confidence level or greater, down to a limiting luminosity (absorbed) of $\sim$1.6\ergs{35} (0.35 -- 8.0 keV). The hardness ratios of the sources separate those with soft, thermal spectra such as supernova remnants from those with hard, non-thermal spectra such as X-ray binaries and background active galactic nuclei. Emission extended beyond the Chandra point spread function is evident in 23 of the 394 sources. Cross-correlation of the ChASeM33 sources against previous catalogs of X-ray sources in M33 results in matches for the vast majority of the brighter sources and shows 28 ChASeM33 sources within 10\arcsec of supernova remnants identified by prior optical and radio searches. This brings the total number of such associations to 31 out of 100 known supernova remnants in M33.

We study, by means of adaptive mesh refinement hydro- and magnetohydrodynamical simulations that cover a wide range of scales (from kpc to sub-parsec), the dimension of the most dissipative structures and the injection scale of the turbulent interstellar gas, which we find to be about 75 pc, in agreement with observations. This is however smaller than the average size of superbubbles, but consistent with significant density and pressure changes in the ISM, which leads to the break-up of bubbles locally and hence to injection of turbulence. The scalings of the structure functions are consistent with log-Poisson statistics of supersonic turbulence where energy is dissipated mainly through shocks. Our simulations are different from previous ones by other authors as (i) we do not assume an isothermal gas, but have temperature variations of several orders of magnitude and (ii) we have no artificial forcing of the fluid with some ad hoc Fourier spectrum, but drive turbulence by stellar explosions at the Galactic rate, self-regulated by density and temperature thresholds imposed on the ISM gas.

The supernova-driven interstellar medium in star-forming galaxies has Reynolds numbers of the order of $10^{6}$ or even larger. We study, by means of adaptive mesh refinement hydro- and magnetohydrodynamical simulations that cover the full available range (from 10 kpc to sub-parsec) scales, the statistical properties of the turbulent interstellar gas and the dimension of the most dissipative structures. The scalings of the structure functions are consistent with a log-Poisson statistics of supersonic turbulence where energy is dissipated mainly through shocks.

It is shown that a number of key observations of the Galactic ISM can be understood, if it is treated as a highly compressible and turbulent medium energized predominantly by supernova explosions (and stellar winds). We have performed extensive numerical high resolution 3D hydrodynamical and magnetohydrodynamical simulations with adaptive mesh refinement over sufficiently long time scales to erase memory effects of the initial setup. Our results show, in good agrement with observations, that (i) volume filling factors of the hot medium are modest (typically below 20%), (ii) global pressure is far from uniform due to supersonic (and to some extent superalfvenic) turbulence, (iii) a significant fraction of the mass (about 60%) in the warm neutral medium is in the thermally unstable regime (500 K < T < 5000 K), (iv) the average number density of OVI in absorption is 1.81 10^{-8} cm^{-3}, in excellent agreement with Copernicus and FUSE data, and its distribution is rather clumpy, consistent with its measured dispersion with distance.

We present a new unbiased search and analysis of all B stars in the solar neighbourhood (within a volume of 400 pc diameter) using the Arivel data base to track down the remains of the OB associations, which hosted the supernovae responsible for the Local Bubble in the interstellar gas. We find after careful dereddening and by comparison with theoretical isochrones, that besides the Upper Scorpius the Upper Centaurus Lupus and Lower Centaurus Crux subgroups are the youngest stellar associations in the solar neighbourhood with ages of 20 to 30 Myr, in agreement with previous work. In search for the ``smoking gun'' of the origin of the Local Bubble, we have traced the paths of the associations back into the past and found that they entered the present bubble region 10 to 15 Myr ago. We argue that the Local Bubble began to form then and estimate that 14 to 20 supernovae have gone off since. It is shown that the implied energy input is sufficient to excavate a bubble of the presently observed size.

The Local and Loop I superbubbles are the closest and best investigated supernova (SN) generated bubbles and serve as test laboratories for observations and theories of the interstellar medium. Since the morphology and dynamical evolution of bubbles depend on the ambient density and pressure distributions, a realistic modelling of the galactic environment is crucial for a detailed comparison with observations. We have performed 3D high resolution (down to 1.25 pc on a kpc-scale grid) hydrodynamic simulations of the Local Bubble (LB) and the neighbouring Loop I (L1) superbubble in a realistically evolving inhomogeneous background ISM, disturbed already by SN explosions at the Galactic rate for 200 Myr before the LB and L1 are generated. The LB is the result of 19 SNe occurring in a moving group, which passed through the present day local HI cavity. We can reproduce (i) the OVI column density in absorption within the LB in agreement with COPERNICUS and recent FUSE observations, giving N(OVI) <2 10^{13} cm^-2 and N(OVI)<7 10^{12} cm^-2, respectively, (ii) the observed sizes of the Local and Loop I superbubbles, (iii) the interaction shell between LB and L1, discovered with ROSAT, (iv) constrain the age of the LB to be 14.5+0.7/-0.4 Myr, (v) predict the merging of the two bubbles in about 3 Myr, when the interaction shell starts to fragment, (vi) the generation of blobs like the Local Cloud as a consequence of a dynamical instability. We find that evolving superbubbles strongly deviate from idealised self-similar solutions due to ambient pressure and density gradients, as well as due to turbulent mixing and mass loading. Hence, at later times the hot interior can break through the surrounding shell, which may also help to explain the puzzling energy "deficit" observed in LMC bubbles.

The statistical properties of interstellar turbulence are studied by means of three-dimensional high-resolution HD and MHD simulations of a SN-driven ISM. It is found that the longitudinal and transverse turbulent length scales have time averaged (over a period of 50 Myr) ratios of 0.5-0.6, almost similar to the one expected for isotropic homogeneous turbulence. The mean characteristic size of the larger eddies is found to be ~75 pc. Furthermore, the scalings of the structure functions measured in the simulated disk show unambiguous departure from the Kolmogorv (1941) model being consistent with the latest intermittency studies of supersonic turbulence (Politano & Pouquet 1995; Boldyrev 2002). Our results are independent of the resolution, indicating that convergence has been reached, and that the unresolved smaller dissipative scales do not feed back on the larger ones.

We study the ovi distribution in space and time in a representative section of the Galactic disk by 3D adaptive mesh refinement HD and MHD simulations of the ISM, including the disk-halo-disk circulation. The simulations describe a supernova driven ISM on large (10 kpc) and small (1.25 pc) scales over a sufficiently large timescale (400 Myrs) in order to establish a global dynamical equilibrium. The Ovi column density, N(Ovi), is monitored through lines of sight measurements at different locations in the simulated disk. One has been deliberately chosen to be inside of a hot bubble, like our own Local Bubble, while the other locations are random. We obtain a correlation between N(Ovi) and distance, which is independent of the observer's vantage point in the disk. In particular, the location of the observer inside a hot bubble does not have any influence on the correlation, because the contribution of an individual bubble (with a typical extension of 100 pc) is negligibly small. We find a remarkable agreement between the Ovi column densities (as a function of distance) and the averaged Ovi density (~1.8x10^{-8}$ cm$^{-3}$) in the disk from our simulations and the values observed with Copernicus, and FUSE. Our results strongly support the important r\^ole of turbulent mixing in the distribution of Ovi clumps in the ISM. Supernova induced turbulence is quite strong and unavoidable due to shearing motions in the ISM and operates on a large range of scales.

Dealing numerically with the turbulent nature and non-linearity of the physical processes involved in the ISM requires the use of sophisticated numerical schemes coupled to HD and MHD mathematical models. SNe are the main drivers of the interstellar turbulence by transferring kinetic energy into the system. This energy is dissipated by shocks (which is more efficient) and by molecular viscosity. We carried out adaptive mesh refinement simulations (with a finest resolution of 0.625 pc) of the turbulent ISM embedded in a magnetic field with mean field components of 2 and 3 $\mu$G. The time scale of our run was 400 Myr, sufficiently long to avoid memory effects of the initial setup, and to allow for a global dynamical equilibrium to be reached in case of a constant energy input rate. It is found that the longitudinal and transverse turbulent length scales have a time averaged (over a period of 50 Myr) ratio of 0.52-0.6, almost similar to the one expected for isotropic homogeneous turbulence. The mean characteristic size of the larger eddies is found to be \~75 pc in both runs. In order to check the simulations against observations, we monitored the OVI and HI column densities within a superbubble created by the explosions of 19 SNe having masses and velocities of the stars that exploded in vicinity of the Sun generating the Local Bubble. The model reproduces the FUSE absorption measurements towards 25 white dwarfs of the \ovi column density as function of distance and of \nhi. In particular for lines of sight with lengths smaller than 120 pc it is found that there is no correlation between N(OVI) and N(HI).

We review recent 3D high resolution numerical HD and MHD studies of a multi-component and multi-phase ISM. The computational grid was chosen large enough in the disk to represent typical ISM gas patterns, and in the halo to include the full extent of the Galactic fountain. The evolution of the ISM is driven by supernovae of all types and followed for 400 Myrs, long enough to get rid of memory effects of the initial distributions. The results are substantially different from previous classical analytical models and entail: (i) a very inhomogeneous distribution of ISM gas, (ii) a high level of supernova driven turbulence, (iii) an active Galactic fountain flow, reaching out to about 5 kpc on either side of the disk, (iv) a low volume filling factor for the hot gas of about 20%, (v) a very weak correlation between magnetic field and gas density, n, for n < 100 cm^-3, and, (vi) about 50% of the gas disk mass with n < 100 cm^-3, residing in classically thermally unstable temperature regimes.

In star forming disk galaxies, matter circulation between stars and the interstellar gas, and, in particular the energy input by random and clustered supernova explosions, determine the dynamical and chemical evolution of the ISM, and hence of the galaxy as a whole. Using a 3D MHD code with adaptive mesh refinement developed for this purpose, we have investigated the r\^ole of magnetized matter circulation between the gaseous disk and the surrounding galactic halo. Special emphasis has been put on the effect of the magnetic field with respect to the volume and mass fractions of the different ISM ``phases'', the relative importance of ram, thermal and magnetic pressures, and whether the field can prevent matter transport from the disk into the halo. The simulations were performed on a grid with an area of 1 kpc$^{2}$, centered on the solar circle, extending $\pm 10$ kpc perpendicular to the galactic disk with a resolution as high as 1.25 pc. The simulations were run for a time scale of 400 Myr, sufficiently long to avoid memory effects of the initial setup, and to allow for a global dynamical equilibrium to be reached in case of a constant energy input rate. (...) We find that in general gas transport into the halo in 3D is not prevented by an initial disk parallel magnetic field, but only delayed initially, for as long as it is needed to punch holes into the thick magnetized gas disk. The mean volume filling factor of the hot phase in the disk is similar in HD and MHD (the latter with a total field strength of 4.4 $\mu$G) runs, amounting to $\sim 17-21%$ for the Galactic supernova rate.

State of the art models of the ISM use adaptive mesh refinement to capture small scale structures, by refining on the fly those regions of the grid where density and pressure gradients occur, keeping at the same time the existing resolution in the other regions. With this technique it became possible to study the ISM in star-forming galaxies in a global way by following matter circulation between stars and the interstellar gas, and, in particular the energy input by random and clustered supernova explosions, which determine the dynamical and chemical evolution of the ISM, and hence of the galaxy as a whole. In this paper we review the conditions for a self-consistent modelling of the ISM and present the results from the latest developments in the 3D HD/MHD global models of the ISM. Special emphasis is put on the effects of the magnetic field with respect to the volume and mass fractions of the different ISM ``phases'', the relative importance of ram, thermal and magnetic pressures, and whether the field can prevent matter transport from the disk into the halo. The simulations were performed on a grid with a square area of 1 kpc$^{2}$, centered on the solar circle, extending $\pm 10$ kpc perpendicular to the galactic disk with a resolution as high as 1.25 pc. The run time scale was 400 Myr, sufficiently long to avoid memory effects of the initial setup, and to allow for a global dynamical equilibrium to be reached in case of a constant energy input rate.

The Local Bubble (LB) is an X-ray emitting region extending 100 pc in radius in the Galactic plane and 300 pc perpendicular to it, and it is embedded in a somewhat larger H{\sc i} deficient cavity. Its origin and spectral properties in UV, EUV and X-rays are still poorly understood. We have performed 3D high resolution (down to 1.25 pc) hydrodynamic superbubble simulations of the LB and Loop I superbubble in a realistic \emph{inhomogeneous} background ISM, disturbed by supernova (SN) explosions at the Galactic rate. We can reproduce (i) the size of the bubbles (in contrast to similarity solutions), (ii) the interaction shell with Loop I, discovered with ROSAT, (iii) predict the merging of the two bubbles in about 3 Myr, when the interaction shell starts to fragment, and, (iv) the generation of blobs like the Local Cloud as a consequence of a dynamical instability. The \ion{O}{6} column densities are monitored and found to be in excellent agreement with \textsc{Copernicus} and FUSE absorption line data, showing LB column densities $<1.7 \times 10^{13} {\rm cm}^{-2}$, in contrast to other existing models.

The observational and theoretical state of Galactic and extragalactic bubbles are reviewed. Observations of superbubbles are discussed, with some emphasis on nearby bubbles such as the Local Bubble (LB) and the Loop I superbubble (LI). Analytical bubble theory is revisited, and similarity solutions, including the time-dependent energy input by supernova explosions according to a Galactic initial mass function (IMF), are studied. Since the agreement with observations is not convincing in case of the LB, we present high resolution 3D AMR simulations of the LB and LI in an inhomogeneous background medium. It is demonstrated that both the morphology and recently published FUSE data on OVI absorption line column densities can be well understood, if the LB is the result of about 20 supernova explosions from a moving group, and the LB age is about 14.7 Myrs.

The role of matter circulation between the disk and halo in establishing the volume filling factors of the different ISM phases in the Galactic disk is investigated, using a modified version of the three-dimensional supernova-driven ISM model of Avillez (2000). We carried out adaptive mesh refinement simulations of the ISM with five supernova rates (in units of the Galactic value), $\sigma/\sigma_{Gal}=1$, 2, 4, 8 and 16 (corresponding to starburst conditions) using three finer level resolutions of 2.5, 1.25 and 0.625 pc, allowing us to understand how resolution would affect the volumes of gas phases in pressure equilibrium. We find that the volume filling factors of the different ISM phases depend sensitively on the existence of a duty cycle between the disk and halo acting as a pressure release valve for the hot (T> 10^5.5 K) phase in the disk.

We have performed high resolution 3D simulations of the Local Bubble (with 1.25 pc finest resolution) in a realistic background ISM, jointly with the dynamical evolution of the neighbouring Loop I superbubble. We can reproduce (i) the size of the bubbles (in contrast to similarity solutions), (ii) the interaction shell with Loop I, discovered with ROSAT, (iii) predict the merging of the two bubbles in about 3 Myr, when the interaction shell starts to fragment, and, (iv) the generation of blobs like the Local Cloud as a consequence of a dynamical instability.

We have performed high resolution 3D simulations with adaptive mesh refinement, following the ISM evolution in a star forming galaxy both on small (<1 pc) and large (>10 kpc) scales, enabling us to track structures in cooling shock compressed regions as well as the entire Galactic fountain flow. It is shown in an MHD run that the latter one is not inhibited by a large scale disk parallel magnetic field. The fountain plays a vital role in limiting the volume filling factor of the hot gas. Contrary to classical models most of the gas between 100K and 8000 K is found to be thermally unstable. On scales of superbubbles we find that the internal temperature structure is rather inhomogeneous for an old object like our Local Bubble, leading to low OVI column densities, consistent with observations.

We carried out 1.25 pc resolution MHD simulations of the ISM -- including the large scale galactic fountain -- by fully tracking the time-dependent evolution of the magnetic field and the formation of shock compressed regions in a supernova-driven ISM. The simulations show that large scale gas streams emerge, driven by SN explosions, which are responsible for the formation and destruction of shocked compressed layers with lifetimes up to some 15 Myr. The T\leq 10^{3} K gas is distributed in filaments which tend to show a preferred orientation due to the anisotropy of the flow induced by the galactic magnetic field. The simulations also show that the magnetic field has a high variability, it is largely uncorrelated with the density and it is driven by inertial motions. The latter is consistent with the fact that ram pressure dominates the flow for 10^{2}<T\leq 10^{6} K. For T> 10^{6} K thermal pressure dominates, while for T\leq 10^{2} K (stable branch) magnetic pressure takes over.

We have carried out 1.25 pc resolution MHD simulations of the ISM, on a Cartesian grid of $0 \leq (x,y) \leq 1$ kpc size in the galactic plane and $-10 \leq z \leq 10$ kpc into the halo, thus being able to fully trace the time-dependent evolution of the galactic fountain. The simulations show that large scale gas streams emerge, driven by SN explosions, which are responsible for the formation and destruction of shocked compressed layers. The shocked gas can have densities as high as 800 cm$^{-3}$ and lifetimes up to 15 Myr. The cold gas is distributed into filaments which tend to show a preferred orientation due to the anisotropy of the flow induced by the galactic magnetic field. Ram pressure dominates the flow in the unstable branch $10^{2}<$T$\leq 10^{3.9}$ K, while for T$\leq 100$ K (stable branch) magnetic pressure takes over. Near supernovae thermal and ram pressures determine the dynamics of the flow. Up to 80% of the mass in the disk is concentrated in the thermally unstable regime $10^{2}<$T$\leq 10^{3.9}$ K with $\sim30%$ of the disk mass enclosed in the T$\leq 10^{3}$ K gas. The hot gas in contrast is controlled by the thermal pressure, since magnetic field lines are swept towards the dense compressed walls.

"The purpose of numerical models is not numbers but insight." (Hamming) In the spirit of this adage, and of Don Cox's approach to scientific speaking, we discuss the questions that the latest generation of numerical models of the interstellar medium raise, at least for us. The energy source for the interstellar turbulence is still under discussion. We review the argument for supernovae dominating in star forming regions. Magnetorotational instability has been suggested as a way of coupling disk shear to the turbulent flow. Models make evident that the unstable wavelengths are very long compared to thermally unstable wavelengths, with implications for star formation in the outer galaxy and low surface brightness disks. The perennial question of the factors determining the hot gas filling factor in a SN-driven medium remains open, in particular because of the unexpectedly strong turbulent mixing at the boundaries of hot cavities seen in the models. The formation of molecular clouds in the turbulent flow is also poorly understood. Dense regions suitable for cloud formation clearly form even in the absence of self-gravity, although their ultimate evolution remains to be computed.

We review large scale modelling of the ISM with emphasis on the importance to include the disk-halo-disk duty cycle and to use a dynamical refinement of the grid (in regions where steep variations of density and pressure occur) for a realistic modelling of the ISM. We also discuss the necessity of convergence of the simulation results by comparing 0.625, 1.25 and 2.5 pc resolution simulations and show that a minimum grid resolution of 1.25 pc is required for quantitatively reliable results, as there is a rapid convergence for Delta x \leq 1.1 pc.

We present the first to date high resolution calculations of the ISM down to scales of 0.625 pc of the global and local ISM. The simulations show the morphology and structure of the different ISM phases and reproduce many of the features that have been observed in the Milky Way and other galaxies. In particular, they show that the hot gas has a moderately low volume filling factor (~20%) even in the absence of magnetic fields. Also, cold gas is mainly concentrated in filamentary structures running perpendicular to the midplane forming and dissipating within 10-12 Myr. Compression is the dominant process for their formation, but thermal instability also plays a role. Also the evolution of the Local Bubble is simulated by multi-supernova explosions; calculated extensions after ~13 Myr match observations.

We study the mixing of chemical species in the interstellar medium (ISM). Recent observations suggest that the distribution of species such as deuterium in the ISM may be far from homogeneous. This raises the question of how long it takes for inhomogeneities to be erased in the ISM, and how this depends on the length scale of the inhomogeneities. We added a tracer field to the three-dimensional, supernova-driven ISM model of Avillez (2000) to study mixing and dispersal in kiloparsec-scale simulations of the ISM with different supernova (SN) rates and different inhomogeneity length scales. We find several surprising results. Classical mixing length theory fails to predict the very weak dependence of mixing time on length scale that we find on scales of 25--500 pc. Derived diffusion coefficients increase exponentially with time, rather than remaining constant. The variance of composition declines exponentially, with a time constant of tens of Myr, so that large differences fade faster than small ones. The time constant depends on the inverse square root of the supernova rate. One major reason for these results is that even with numerical diffusion exceeding physical values, gas does not mix quickly between hot and cold regions.

In order to understand the time scales for complete mixing in a supernova-driven ISM we used the model of Avillez (2000) to follow the mixing and dispersal of inhomogeneities with different length scales. We included tracer fields to follow the inhomogeneities in kpc-scale simulations with different supernova rates. We find that an initially inhomogeneous ISM still shows local variations and a clumpy distribution on time scales of several hundred Myr. This clumpy structure does not disapear over a time of a few tens of Myr, even if the rate of SNe is increased to ten times the Galactic rate. These local variations occur because, even in the presence of numerical diffusion, gas does not mix quickly between hot and cold regions. The lower limit that we are able to place on the mixing time at the kpc scale is already longer than the time for chemistry to reach local equilibrium in the gas. The simulations also show that the time scale to erase inhomogeneities in the ISM is nearly independent of their length scale over the range 25-500 pc, contrary to what is expected from classical mixing length theory.

Recent HI emission observations of the Southern Galactic hemisphere have revealed a mushroom-like structure extending from z=-70 to -450 pc, composed of a stem and a cap. Similar structures occur in three-dimensional simulations of a dynamic galactic disk driven by isolated and clustered supernovae. Using these simulations, we show that hot gas in the Galactic disk that is not evacuated through chimneys expands into the cooler gas of the thick disk, forming mushroom-shaped structures. This new class of objects traces buoyant flow of hot gas into the thick disk.

The results of a three-dimensional model for disk-halo interaction are presented here. The model considers explicitly the input of energy and mass by isolated a nd correlated supernovae in the disk. Once disrupted by the explosions, the disk never returns to its initial state. Instead it approaches a state where a thin HI disk is formed in the Galactic plane overlayed by thick HI and HII gas disk w ith scale heights of 500 pc and of 1 to 1.5 kpc, respectively. The upper parts o f the thick HII disk (the diffuse ionized medium) act as a disk-halo interface a nd its formation and stability are directly correlated to the supernova rate per unit area in the simulated disk.

Three-dimensional simulations of the disk-halo interaction show the formation of a thick HI and HII gas disk with different scale heights. The thick HI disk prevents the disk gas from expanding freely upwards, unless some highly energetic event such as chimneys occurs, whereas the thick HII disk acts as a disk-halo interaction region from where the hot ionized gas flows freely into the halo. The upflowing gas reaches the maximum height at $z\sim 9.3\pm 1$ kpc becoming thermally unstable due to radiative losses, and condenses into HI clouds. Because the major fraction of the gas is gravitationally bound to the Galaxy, the cold gas returns to the disk. The descending clouds will have at some height high velocities. In a period of 200 Myr of fountain evolution, some 10 percent of the total number of clouds are HVCs.

Gas that escapes from the Galactic disk, rises into the halo, cools and falls back, constitutes a "Galactic Fountain". Three-Dimensional simulations show that such a fountain model reproduces many of the features that have been observed in the Galaxy and other galaxies such as M31 and M33. Here, these results are reported.