Core-collapse supernova (SN) ejecta are probably structured on both small and large scales, with greater deviations from spherical symmetry nearer the explosion site. Here, we present 2D and 3D gray radiation-hydrodynamics simulations of type II SN light curves from red (RSG) and blue supergiant (BSG) star explosions to investigate the impact on SN observables of inhomogeneities in density or composition, with a characteristic scale set to a few percent of the local radius. Clumping is found to hasten the release of stored radiation, boosting the early time luminosity and shortening the photospheric phase. Around the photosphere, radiation leaks between the clumps where the photon mean free path is greater. Since radiation is stored uniformly in volume, a greater clumping can increase this leakage by storing more and more mass into smaller and denser clumps containing less and less radiation energy. An inhomogeneous medium in which different regions recombine at different temperatures can also impact the light curve. Clumping can thus be a source of diversity in SN brightness. Clumping may lead to a systematic underestimate of ejecta masses from light curve modeling, although a significant offset seems to require a large density contrast of a few tens between clumps and interclump medium.

Convection is an important physical process in astrophysics well-studied using numerical simulations under the Boussinesq and/or anelastic approximations. However these approaches reach their limits when compressible effects are important in the high Mach flow regime, e.g. in stellar atmospheres or in the presence of accretion shocks. In order to tackle these issues, we propose a new high performance and portable code, called "ARK" with a numerical solver well-suited for the stratified compressible Navier-Stokes equations. We take a finite volume approach with machine precision conservation of mass, transverse momentum and total energy. Based on previous works in applied mathematics we propose the use of a low Mach correction to achieve a good precision in both low and high Mach regimes. The gravity source term is discretized using a well-balanced scheme in order to reach machine precision hydrostatic balance. This new solver is implemented using the Kokkos library in order to achieve high performance computing and portability across different architectures (e.g. multi-core, many-core, and GP-GPU). We show that the low-Mach correction allows to reach the low-Mach regime with a much better accuracy than a standard Godunov-type approach. The combined well-balanced property and the low-Mach correction allowed us to trigger Rayleigh-B\'enard convective modes close to the critical Rayleigh number. Furthermore we present 3D turbulent Rayleigh-B\'enard convection with low diffusion using the low-Mach correction leading to a higher kinetic energy power spectrum. These results are very promising for future studies of high Mach and highly stratified convective problems in astrophysics.

By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth oceans, fingering convection in stellar atmospheres, and moist convection in Earth atmosphere are deriving from the same general diabatic convective instability. We show also that "radiative convection" triggered by CO/CH4 transition with radiative transfer in the atmospheres of brown dwarfs is analog to moist and thermohaline convection. We derive a generalization of the mixing length theory to include the effect of source terms in 1D codes. We show that CO/CH4 radiative convection could significantly reduce the temperature gradient in the atmospheres of brown dwarfs similarly to moist convection in Earth atmosphere thus possibly explaining the reddening in brown-dwarf spectra. By using idealized two-dimensional hydrodynamic simulations in the Ledoux unstable regime, we show that compositional source terms can indeed provoke a reduction of the temperature gradient. The L/T transition could be explained by a bifurcation between the adiabatic and diabatic convective transports and could be seen as a giant cooling crisis: an analog of the boiling crisis in liquid/steam-water convective flows. This mechanism with other chemical transitions could be present in many giant and earth-like exoplanets. The study of the impact of different parameters (effective temperature, compositional changes) on CO/CH4 radiative convection and the analogy with Earth moist and thermohaline convection is opening the possibility to use brown dwarfs to better understand some aspects of the physics at play in the climate of our own planet.

Magnetar power is believed to be at the origin of numerous super-luminous supernovae (SNe) of Type Ic, arising from compact, hydrogen-deficient, Wolf-Rayet type stars. Here, we investigate the properties that magnetar power would have on standard-energy SNe associated with 15-20Msun blue or red supergiant (BSG/RSG) stars. We use a combination of Eulerian grey radiation-hydrodynamics and non-LTE steady-state radiative transfer to study their dynamical, photometric, and spectroscopic properties. Adopting magnetar fields of 1, 3.5, 7 x 10^14G and rotational energies of 0.4, 1, and 3 x 10^51erg, we produce bolometric light curves with a broad maximum covering 50-150d and a magnitude of 10^43-10^44erg/s. The spectra at maximum light are analogous to those of standard SNe II-P but bluer. Although the magnetar energy is channelled roughly in equal proportion between SN kinetic energy and SN luminosity, the latter may be boosted by a factor 10-100 compared to a standard SN II. This influence breaks the observed relation between brightness and expansion rate of standard Type II SNe. Magnetar energy injection also delays recombination and may even cause re-ionization, with a reversal in photospheric temperature and velocity. Depositing the magnetar energy in a narrow mass shell at the ejecta base leads to the formation of a dense shell at a few 1000km/s, which causes a light-curve bump at the end of the photospheric phase. Depositing this energy over a broad range of mass in the inner ejecta, to mimic the effect of multi-dimensional fluid instabilities, prevents the formation of a dense shell and produces an earlier-rising and smoother light curve. The magnetar influence on the SN radiation is generally not visible prior to 20-30d, during which one may discern a BSG from a RSG progenitor. We propose a magnetar model for the super-luminous Type II SN OGLE-SN14-073.

Early-time observations of the Type II supernovae (SNe) 2013cu and 2013fs have revealed an interaction of ejecta with material near the star surface. Unlike the Type IIn SN2010jl, which interacts with a dense wind for ~1yr, the interaction ebbs after 2-3d, suggesting a dense and compact circumstellar envelope. Here, we use multi-group radiation-hydrodynamics and non-local-thermodynamic-equilibrium radiative transfer to explore the properties of red supergiant (RSG) star explosions embedded in a variety of dense envelopes. We consider the cases of an extended static atmosphere or a steady-state wind, adopting a range of mass loss rates. The shock-breakout signal, the SN radiation up to 10d, and the ejecta dynamics are strongly influenced by the properties of this nearby environment. This compromises the use of early-time observations to constrain Rstar. The presence of narrow lines for 2-3d in 2013fs and 2013cu require a cocoon of material of ~0.01Msun out to 5-10Rstar. Spectral lines evolve from electron-scattering to Doppler broadened, with a growing blueshift of their emission peaks. Recent studies propose a super-wind phase with a mass loss rate from 0.001 up to 1Msun/yr in the last months/years of the RSG life, although there is no observational constraint that this external material is a steady-state outflow. Alternatively, observations may be explained by the explosion of a RSG star inside its complex atmosphere. Indeed, spatially resolved observations reveal that RSG stars have extended atmospheres, with the presence of downflows and upflows out to several Rstar, even in a standard RSG like Betelgeuse. Mass loading in the region intermediate between star and wind can accommodate the 0.01Msun needed to explain the observations of 2013fs. A puzzling super-wind phase prior to core-collapse may be therefore superfluous [Abridged].

Some interacting supernovae (SNe) of type IIn show a sizeable continuum polarisation suggestive of a large scale asymmetry in the circumstellar medium (CSM) and/or the SN ejecta. Here, we extend the recent work of Dessart et al. on super-luminous SNe IIn and perform axially-symmetric (i.e., 2D) multi-group radiation hydrodynamics simulations to explore the impact of an imposed large scale density asymmetry. When the CSM is asymmetric, the latitudinal variation of the radial optical depth $\tau$ introduces a strong flux redistribution from the higher-density CSM regions, where the shock luminosity is larger, towards the lower-density CSM regions where photons escape more freely --- this redistribution ceases when $\tau<$1. Along directions where the CSM density is larger, the shock deceleration is stronger and its progression slower, producing a non-spherical cold-dense shell (CDS). For an oblate CSM density distribution, the photosphere (CDS) has an oblate (prolate) morphology when $\tau>$1. When the CSM is symmetric and the ejecta asymmetric, the flux redistribution within the CSM now tends to damp the latitudinal variation of the luminosity at the shock. It then requires a larger ejecta asymmetry to produce a sizeable latitudinal variation in the emergent flux. When the interaction is between a SN ejecta and a relic disk, the luminosity boost at early times scales with the disk opening angle -- forming a super-luminous SN IIn this way requires an unrealistically thick disk. In contrast, interaction with a disk of modest thickness/mass can yield a power that rivals radioactive decay of a standard SN II at nebular times.

Using radiation-hydrodynamics and radiative-transfer simulations, we explore the origin of the spectral diversity of interacting supernovae (SNe) of type IIn. We revisit SN1994W and investigate the dynamical configurations that can give rise to spectra with narrow lines at all times. We find that a standard ~10Msun 10^51erg SN ejecta ramming into a 0.4Msun dense CSM is inadequate for SN1994W, as it leads to the appearance of broad lines at late times. This structure, however, generates spectra that exhibit the key morphological changes seen in SN1998S. For SN1994W, we consider a completely different configuration, which involves the interaction at a large radius of a low mass inner shell with a high mass outer shell. Such a structure may arise in an 8-12Msun star from a nuclear flash (e.g., of Ne) followed within a few years by core collapse. Our simulations show that the large mass of the outer shell leads to the complete braking of the inner shell material, the formation of a slow dense shell, and the powering of a luminous SN IIn, even for a low inner shell energy. Early on, our model line profiles are typical of SNe IIn, exhibiting narrow cores and broad electron-scattering wings. As observed in SN1994W, they also remain narrow at late times. Our SN1994W model invokes two low energy ejections, both atypical of observed massive stars, and illustrates the diversity of configurations leading to SNe IIn. These results also highlight the importance of spectra to constrain the dynamical properties and understand the origin of SNe IIn.

We present numerical simulations that include 1-D Eulerian multi-group radiation-hydrodynamics, 1-D non-LTE radiative transfer, and 2-D polarised radiative transfer for super-luminous interacting supernovae (SNe). Our reference model is a ~10Msun inner shell with 10^51erg ramming into a ~3Msun cold outer shell (the circumstellar-medium, or CSM) that extends from 10^15cm to 2x10^16cm and moves at 100km/s. We discuss the light curve evolution, which cannot be captured adequately with a grey approach. In these interactions, the shock-crossing time through the optically-thick CSM is much longer than the photon diffusion time. Radiation is thus continuously leaking from the shock through the CSM, in disagreement with the shell-shocked model that is often invoked. Our spectra redden with time, with a peak distribution in the near-UV during the first month gradually shifting to the optical range over the following year. Initially Balmer lines exhibit a narrow line core and the broad line wings that are characteristic of electron scattering in the SNe IIn atmospheres (CSM). At later times they also exhibit a broad blue shifted component which arises from the cold dense shell. Our model results are broadly consistent with the bolometric light curve and spectral evolution observed for SN2010jl. Invoking a prolate pole-to-equator density ratio in the CSM, we can also reproduce the ~2% continuum polarisation, and line depolarisation, observed in SN2010jl. By varying the inner shell kinetic energy and the mass and extent of the outer shell, a large range of peak luminosities and durations, broadly compatible with super-luminous SNe IIn like 2010jl or 2006gy, can be produced.

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses.

Radiative shocks are shocks in a gas where the radiative energy and flux coming from the very hot post-shock material are non-negligible in the shock's total energy budget, and are often large enough to heat the material ahead of the shock. Many simulations of radiative shocks, both in the contexts of astrophysics and laboratory experiments, use a grey treatment of radiative transfer coupled to the hydrodynamics. However, the opacities of the gas show large variations as a function of frequency and this needs to be taken into account if one wishes to reproduce the relevant physics. We have performed radiation hydrodynamics simulations of radiative shocks in Ar using multigroup (frequency dependent) radiative transfer with the HERACLES code. The opacities were taken from the ODALISC database. We show the influence of the number of frequency groups used on the dynamics and morphologies of subcritical and supercritical radiative shocks in Ar gas, and in particular on the extent of the radiative precursor. We find that simulations with even a low number of groups show significant differences compared to single-group (grey) simulations, and that in order to correctly model such shocks, a minimum number of groups is required. Results appear to eventually converge as the number of groups increases above 50. We were also able to resolve in our simulations of supercritical shocks the adaptation zones which connect the cooling layer to the final post-shock state and the precursor. Inside these adaptation zones, we find that the radiative flux just ahead of the shock in one or several high-opacity groups can heat the gas to a temperature higher than the post-shock temperature. Through the use of Hugoniot curves, we have checked the consistency of our radiation hydrodynamics scheme by showing that conservation of mass, momentum and energy (including radiative flux) holds. ABRIDGED.

This paper is a numerical study of the condensation of the warm neutral medium (WNM) into cold neutral medium (CNM) structures under the effect of turbulence and thermal instability. Using low resolution simulations we explored the impact of the WNM initial density and properties of the turbulence (stirring in Fourier with a varying mix of solenoidal and compressive modes) on the cold gas formation. Two sets of initial conditions which match the observations were selected to produce high resolution simulations (1024^3) allowing to study in details the properties of the produced dense structures. For typical values of the density, pressure and velocity dispersion of the WNM in the solar neighborhood, the turbulent motions of the HI can not provoque the phase transition from WNM to CNM, whatever their amplitude and their distribution in solenoidal and compressive modes. On the other hand we show that a quasi-isothermal increase in WNM density of a factor of 2 to 4 is enough to induce the phase transition, leading to the transition of about 40 percent of the gas to the cold phase within 1 Myr. Given the observed properties of the HI in the local ISM, the WNM and individual CNM structures in the local ISM are sub or transsonic and their dynamics are tightly interwoven. The velocity field bears the evidence of subsonic turbulence with a 2D power spectrum following the Kolmogorov law as P(k) \sim k^{-8/3} while the density is highly contrasted with a singificantly shallower power spectrum, reminiscent of what is observed in the cold ISM. Supra-thermal line width observed for CNM might be the result of relative velocity between cold structures. Finally, the cold structures denser than 5 cm^{-3} reproduce well the laws M \sim L^{2.25-2.28} and sigma(v) \sim 0.5-0.8L^{1/3} generally observed in molecular clouds.

We investigate the interplay between the ionization radiation from massive stars and the turbulence inside the surrounding molecular gas thanks to 3D numerical simulations. We used the 3D hydrodynamical code HERACLES to model an initial turbulent medium that is ionized and heated by an ionizing source. Three different simulations are performed with different mean Mach numbers (1, 2 and 4). A non-equilibrium model for the ionization and the associated thermal processes was used. This revealed to be crucial when turbulent ram pressure is of the same order as the ionized-gas pressure. The density structures initiated by the turbulence cause local curvatures of the dense shell formed by the ionization compression. When the curvature of the shell is sufficient, the shell collapse on itself to form a pillar while a smaller curvature leads to the formation of dense clumps that are accelerated with the shell and therefore remain in the shell during the simulation. When the turbulent ram pressure of the cold gas is sufficient to balance the ionized-gas pressure, some dense-gas bubbles have enough kinetic energy to penetrate inside the ionized medium, forming cometary globules. This suggests a direct relation in the observations between the presence of globules and the relative importance of the turbulence compared to the ionized-gas pressure. The probability density functions present a double peak structure when the turbulence is low relative to the ionized-gas pressure. This could be used in observations as an indication of the turbulence inside molecular clouds.

Radiative transfer plays a major role in the process of star formation. Many simulations of gravitational collapse of a cold gas cloud followed by the formation of a protostellar core use a grey treatment of radiative transfer coupled to the hydrodynamics. However, dust opacities which dominate extinction show large variations as a function of frequency. In this paper, we used frequency-dependent radiative transfer to investigate the influence of the opacity variations on the properties of Larson's first core. We used a multigroup M1 moment model in a 1D radiation hydrodynamics code to simulate the spherically symmetric collapse of a 1 solar mass cloud core. Monochromatic dust opacities for five different temperature ranges were used to compute Planck and Rosseland means inside each frequency group. The results are very consistent with previous studies and only small differences were observed between the grey and multigroup simulations. For a same central density, the multigroup simulations tend to produce first cores with a slightly higher radius and central temperature. We also performed simulations of the collapse of a 10 and 0.1 solar mass cloud, which showed the properties of the first core to be independent of the initial cloud mass, with again no major differences between grey and multigroup models. For Larson's first collapse, where temperatures remain below 2000 K, the vast majority of the radiation energy lies in the IR regime and the system is optically thick. In this regime, the grey approximation does a good job reproducing the correct opacities, as long as there are no large opacity variations on scales much smaller than the width of the Planck function. The multigroup method is however expected to yield more important differences in the later stages of the collapse when high energy (UV and X-ray) radiation is present and matter and radiation are strongly decoupled.

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

Radiative feedback from massive stars is a key process to understand how HII regions may enhance or inhibit star formation in pillars and globules at the interface with molecular clouds. We aim to contribute to model the interactions between ionization and gas clouds to better understand the processes at work. We study in detail the impact of modulations on the cloud-HII region interface and density modulations inside the cloud. We run three-dimensional hydrodynamical simulations based on Euler equations coupled with gravity using the HERACLES code. We implement a method to solve ionization/recombination equations and we take into account typical heating and cooling processes at work in the interstellar medium and due to ionization/recombination physics. UV radiation creates a dense shell compressed between an ionization front and a shock ahead. Interface modulations produce a curved shock that collapses on itself leading to stable growing pillar-like structures. The narrower the initial interface modulation, the longer the resulting pillar. We interpret pillars resulting from density modulations in terms of the ability of these density modula- tions to curve the shock ahead the ionization front. The shock curvature is a key process to understand the formation of structures at the edge of HII regions. Interface and density modulations at the edge of the cloud have a direct impact on the morphology of the dense shell during its formation. Deeper in the cloud, structures have less influence due to the high densities reached by the shell during its expansion.

Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper.

Radiative shocks play a dominant role in star formation. The accretion shocks on the first and second Larson's cores involve radiative processes and are thus characteristic of radiative shocks. In this study, we explore the formation of the first Larson's core and characterize the radiative and dynamical properties of the accretion shock, using both analytical and numerical approaches. We develop both numerical RHD calculations and a semi-analytical model that characterize radiative shocks in various physical conditions, for radiating or barotropic fluids. Then, we perform 1D spherical collapse calculations of the first Larson's core, using a grey approximation for the opacity of the material. We consider three different models for radiative transfer, namely: the barotropic approximation, the FLD approximation and the more complete M1 model. We investigate the characteristic properties of the collapse and of the first core formation. Comparison between the numerical results and our semi-analytical model shows that this latter reproduces quite well the core properties obtained with the numerical calculations. The accretion shock on the first Larson core is found to be supercritical, implying that all the accretion shock energy on the core is radiated away. The FLD approximation is found to agree quite well with the results based on the M1 model, and is thus appropriate to study the star formation process. In contrast, the barotropic approximation does not correctly describe the thermal properties of the gas during the collapse. We show that a consistent treatment of radiation and hydrodynamics is mandatory to correctly handle the cooling of the gas during the core formation and thus to obtain the correct mechanical and thermal properties for this latter.

We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M1) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well established kinetic model through Marshak wave tests with frequency dependent opacities. Very good agreement between the multigroup M1 and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.

Radiative feedback and magnetic field are understood to have a strong impact on the protostellar collapse. We present high resolution numerical calculations of the collapse of a 1 solar mass dense core in solid body rotation, including both radiative transfer and magnetic field. Using typical parameters for low-mass cores, we study thoroughly the effect of radiative transfer and magnetic field on the first core formation and fragmentation. We show that including the two aforementioned physical processes does not correspond to the simple picture of adding them separately. The interplay between the two is extremely strong, via the magnetic braking and the radiation from the accretion shock.

We study the rotational properties of magnetized and self-gravitating molecular cloud cores formed in 2 very high resolution 3D molecular cloud simulations.The simulations have been performed using the code RAMSES at an effective resolution of 4096^3.One simulation represents a mildly magnetically-supercritical cloud and the other a strongly magnetically-supercritical cloud.A noticeable difference between the 2 simulations is the core formation efficiency (CFE) of the high density cores.In the strongly supercritical simulations the CFE is ~33 % per free-fall time of the cloud tff,cl, whereas in the mildly supercritical simulations this value goes down to ~6%/tff,cl. A comparison of the intrinsic specific angular momentum j3D distributions of the cores with the distribitions of j2D derived using synthetic 2D velocity maps of the cores,shows that the synthetic observations tend to overestimate the true value of j by a factor of ~10.The origin of this discrepancy lies in the fact that contrary to the intrinsic determination which sums up the individual gas parcels contributions to j, the determination of j using the observational procedure which is based on a measurement on the global velocity gradient under the hypothesis of uniform rotation smoothes out the complex fluctuations present in the 3D velocity field. Our results provide a natural explanation for the discrepancy by a factor ~10 observed between the intrinsic 3D distributions of j and the corresponding distributions derived in real observations.We suggest that measurements of j which are based on the measurement of the observed global velocity gradients may need to be reduced by a factor of ~10 in order to derive a more accurate estimate of j in the cores.

We have studied the dynamics of strong radiative shocks generated with the high-energy subnanosecond iodine laser at Prague Asterix Laser System facility

We report on the investigation of strong radiative shocks generated with the high energy, sub-nanosecond iodine laser at PALS. These shock waves are characterized by a developed radiative precursor and their dynamics is analyzed over long time scales (~50 ns), approaching a quasi-stationary limit. We present the first preliminary results on the rear side XUV spectroscopy. These studies are relevant to the understanding of the spectroscopic signatures of accretion shocks in Classical T Tauri Stars.

It is established that both radiative transfer and magnetic field have a strong impact on the collapse and the fragmentation of prestellar dense cores, but no consistent calculation exists yet at such scales. We present original AMR calculations including magnetic field (in the ideal MHD limit) and radiative transfer, within the Flux Limited Diffusion approximation, of the collapse of a 1 solar mass dense core. We compare the results with calculations performed with a barotropic EOS. We show that radiative transfer has an important impact on the collapse and the fragmentation, through the cooling or heating of the gas, and is complementary of the magnetic field. A larger field yields a stronger magnetic braking, increasing the accretion rate and thus the effect of the radiative feedback. Even for a strongly magnetized core, where the dynamics of the collapse is dominated by the magnetic field, radiative transfer is crucial to determine the temperature and optical depth distributions, two potentially accessible observational diagnostics. A barotropic EOS cannot account for realistic fragmentation. The diffusivity of the numerical scheme, however, is found to strongly affect the output of the collapse, leading eventually to spurious fragmentation. Both radiative transfer and magnetic field must be included in numerical calculations of star formation to obtain realistic collapse configurations and observable signatures. Nevertheless, the numerical resolution and the robustness of the solver are of prime importance to obtain reliable results. When using an accurate solver, the fragmentation is found to always remain inhibited by the magnetic field, at least in the ideal MHD limit, even when radiative transfer is included.

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

{It is well established that the atomic interstellar hydrogen is filling the galaxies and constitutes the building blocks of molecular clouds.} {To understand the formation and the evolution of molecular clouds, it is necessary to investigate the dynamics of turbulent and thermally bistable as well as barotropic flows.} {We perform high resolution 3-dimensional hydrodynamical simulations of 2-phase, isothermal and polytropic flows.} {We compare the density PDF and Mach number density relation in the various simulations and conclude that 2-phase flows behave rather differently than polytropic flows. We also extract the clumps and study their statistical properties such as mass spectrum, mass-size relation and internal velocity dispersion. In each case, it is found that the behaviour is well represented by a simple powerlaw. While the various exponents inferred are very similar for the 2-phase, isothermal and polytropic flows, we nevertheless find significant differences, such as for example the internal velocity dispersion which is smaller for 2-phase flows.} {The structures statistics are very similar to what has been inferred from observations, in particular the mass spectrum, the mass-size relation and the velocity dispersion-size relation are all powerlaws whose indices well agree with the observed values. Our results suggest that in spite of various statistics being similar for 2-phase and polytropic flows, they nevertheless present significant differences, stressing the necessity to consider the proper thermal structure of the interstellar atomic hydrogen for computing its dynamics as well as the formation of molecular clouds.}

Supernova explosions inject a considerable amount of energy into the interstellar medium (ISM) in regions with high to moderate star formation rates. In order to assess whether the driving of turbulence by supernovae is also important in the outer Galactic disk, where the star formation rates are lower, we study the spatial distribution of molecular cloud (MC) inclinations with respect to the Galactic plane. The latter contains important information on the nature of the mechanism of energy injection into the ISM. We analyze the spatial correlations between the position angles (PAs) of a selected sample of MCs (the largest clouds in the catalogue of the outer Galaxy published by Heyer et al. 2001). Our results show that when the PAs of the clouds are all mapped to values into the [0,90]degrees interval, there is a significant degree of spatial correlation between the $PA$s on spatial scales in the range of 100-800 pc. These scales are of the order of the sizes of individual SN shells in low density environments such as those prevailing in the outer Galaxy and where the metallicity of the ambient gas is of the order of the solar value or smaller. These findings suggest that individual SN explosions, occurring in the outer regions of the Galaxy and in likewise spiral galaxies, albeit at lower rates, continue to play an important role in shaping the structure and dynamics of the ISM in those regions. The SN explosions we postulate here are likely associated with the existence of young stellar clusters in the far outer regions of the Galaxy and the UV emission and low levels of star formation observed with the GALEX satellite in the outer regions of local galaxies.

Radiative shocks are found in various astrophysical objects and particularly at different stages of stellar evolution. Studying radiative shocks, their topology, and thermodynamical properties is therefore a starting point to understanding their physical properties. This study has become possible with the development of large laser facilities, which has provided fresh impulse to laboratory astrophysics. We present the main characteristics of radiative shocks modeled using cylindrical simulations. We focus our discussion on the importance of multi-dimensional radiative-transfer effects on the shock topology and dynamics. We present results obtained with our code HERACLES for conditions corresponding to experiments already performed on laser installations. The multi-dimensional hydrodynamic code HERACLES is specially adapted to laboratory astrophysics experiments and to astrophysical situations where radiation and hydrodynamics are coupled. The importance of the ratio of the photon mean free path to the transverse extension of the shock is emphasized. We present how it is possible to achieve the stationary limit of these shocks in the laboratory and analyze the angular distribution of the radiative flux that may emerge from the walls of the shock tube. Implications of these studies for stellar accretion shocks are presented.

{It has recently been proposed that giant molecular complexes form at the sites where streams of diffuse warm atomic gas collide at transonic velocities.} {We study the global statistics of molecular clouds formed by large scale colliding flows of warm neutral atomic interstellar gas under ideal MHD conditions. The flows deliver material as well as kinetic energy and trigger thermal instability leading eventually to gravitational collapse.} {We perform adaptive mesh refinement MHD simulations which, for the first time in this context, treat self-consistently cooling and self-gravity.} {The clouds formed in the simulations develop a highly inhomogeneous density and temperature structure, with cold dense filaments and clumps condensing from converging flows of warm atomic gas. In the clouds, the column density probability density distribution (PDF) peaks at $\sim 2 \times 10^{21} \psc$ and decays rapidly at higher values; the magnetic intensity correlates weakly with density from $n \sim 0.1$ to $10^4 \pcc$, and then varies roughly as $n^{1/2}$ for higher densities.} {The global statistical properties of such molecular clouds are reasonably consistent with observational determinations. Our numerical simulations suggest that molecular clouds formed by the moderately supersonic collision of warm atomic gas streams.}

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare thoroughly these two methods within their standard framework. We use the AMR code RAMSES and the SPH code DRAGON. Our physical model is as simple as possible and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled. Higher resolution figures can be downloaded at http://www-dapnia.cea.fr/Projets/COAST/paper_amrvssph.pdf

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare here thoroughly these two methods with numerical resolution requirements deduced from previous studies. Our physical model is as simple as possible, and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled.

The interstellar atomic hydrogen is known to be a 2-phase medium in which turbul ence plays an important r\^ole. Here we present high resolution numerical simulations describing the gas from tens of parsec down to hundreds of AU. This high resolut ion allows to probe numerically, the small scale structures which naturally arises from the turbulence and the 2-phase physics.

{We study in some details the statistical properties of the turbulent 2-phase interstellar atomic gas.{We present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen which describe it over 3 to 4 orders of magnitude in spatial scales.}{The simulations produce naturally small scale structures having either large or small column density. It is tempting to propose that the former are connected to the tiny small scale structures observed in the ISM. We compute the mass spectrum of CNM structures and find that ${\cal N}(M) dM \propto M ^{-1.7} dM$, which is remarkably similar to the mass spectrum inferred for the CO clumps. We propose a theoretical explanation based on a formalism inspired from the Press & Schecter (1974) approach and used the fact that the turbulence within WNM is subsonic. This theory predicts ${\cal N}(M) \propto M ^{-5/3}$ in 2D and ${\cal N}(M) \propto M ^{-16/9}$ in 3D. We compute the velocity and the density power-spectra and conclude that, although the latter is rather flat, as observed in supersonic isothermal simulations, the former follows the Kolmogorov prediction and is dominated by its solenoidal component. This is due to the bistable nature of the flow which produces large density fluctuations even when the rms Mach number (of WNM) is not large. We also find that, whereas the energy at large scales is mainly in the WNM, at smaller scales, it is dominated by the kinetic energy of the CNM fragments.}

{It is necessary to understand the dynamics of the atomic gas to use complex modeling and to carry out detailed comparisons between theoretical models and observations.}{In a companion paper, we present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen. Here, we further characterize these simulations and we compare our results with various observations.}{We give statistics of the column density and velocity along the line of sight and show that they compare favorably with observations of high-latitude lines of sight. We compute synthetic HI spectra and qualitatively discuss the information that could be inferred if these spectra were observed. Finally, we extract CNM clouds and study their physical properties finding strong similarities with real clouds. In particular, we find that the clouds follow Larson-type relations, i.e. $M \propto L^\gamma$, where $\gamma \simeq 1.7$ (we speculate that in 3D, $\gamma \simeq 2.5$) and $\sqrt{<\delta v^2>} \propto L^{0.4}$. We also find that the distribution, ${\cal N}(N)$, of the column density, $N$, of the CNM structures formed in the simulation follows ${\cal N}(N) \propto N^{-1.2}$ which is marginally compatible with the observational result obtained by Heiles & Troland (2005). From the mass-size relation and the mass spectrum, we derive an exponent for the column density distribution close to the value obtained in the numerical simulation.} {We conclude that the simulations reproduce various observational features reasonably well....}

We present a numerical and analytical study of the thermal fragmentation of a turbulent flow of interstellar hydrogen. We first present the different dynamical processes and the large range of spatial (and temporal) scales that need to be adequately represented in numerical simulations. Next, we present bidimensionnal simulations of turbulent converging flows which induce the dynamical condensation of the warm neutral phase into the cold phase. We then analyse the cold structures and the fraction of unstable gas in each simulation, paying particular attention to the influence of the degree of turbulence. When the flow is very turbulent a large fraction of the gas remains in the thermally unstable domain. This unstable gas forms a filamentary network. We show that the fraction of thermally unstable gas is strongly correlated with the level of turbulence of the flow. We then develop a semi-analytical model to explain the origin of this unstable gas. This simple model is able to reproduce quantitatively the fraction of unstable gas observed in the simulations and its correlation with turbulence. Finally, we stress the fact that even when the flow is very turbulent and in spite of the fact that a large fraction of the gas is maintained dynamically in the thermally unstable domain, the classical picture of a 2-phase medium with stiff thermal fronts and local pressure equilibrium turns out to be still relevant in the vicinity of the cold structures.

We present in this paper the numerical treatment of the coupling between hydrodynamics and radiative transfer. The fluid is modeled by classical conservation laws (mass, momentum and energy) and the radiation by the grey moment $M_1$ system. The scheme introduced is able to compute accurate numerical solution over a broad class of regimes from the transport to the diffusive limits. We propose an asymptotic preserving modification of the HLLE scheme in order to treat correctly the diffusion limit. Several numerical results are presented, which show that this approach is robust and have the correct behavior in both the diffusive and free-streaming limits. In the last numerical example we test this approach on a complex physical case by considering the collapse of a gas cloud leading to a proto-stellar structure which, among other features, exhibits very steep opacity gradients.

We address the extended problem of component separation for CMB applications when a mixture of both astrophysical and instrumental components are present in the observations, and show how standard methods can be adapted to handle this more complex inversion problem.

Scattering of starlight by dust, molecules and electrons in spiral galaxies will produce a modification of the direct intensity and a polarization in the observed light. We treat the case where the distribution of scatterers can be considered to be optically thin, and derive semi-analytic expressions for the resolved intensity and polarized intensity for Thomson, Rayleigh, and more general scattering mechanisms. These expressions are applied to a parametric model spiral galaxies. It is further shown that in the case of Thomson and Rayleigh scattering, and when scatterers and stars are distributed with rotational symmetry, the total polarized flux depends on the inclination, $i$, of the galactic axis to the line of sight according to a simple $\sin ^2 i$ law. This generalises the well known result for pointlike and spherical light sources. By using a method based on spherical harmonics, we generalise this law for more general mechanisms, and show that to good approximation, the $\sin ^2 i$ law still holds for the class of models considered.

The polarization of the CMBR scattered by galaxy clusters in the kinematic Sunyaev Zeldovich effect depends on the transverse velocity of the cluster. This polarizing effect is proportional to the transverse velocity squared, and so weaker that the change in intensity due to the radial motion in the kinematic effect. The value given by Sunyaev and Zeldovich, and which is frequently cited, underestimates the polarizing effect by a factor of ten. We show furthermore that the polarization has a strong frequency dependence. This means that the polarization should be detectable with the new generation of CMBR probes, at least for some clusters. Thus this effect offers, almost uniquely, a method of obtaining the vectorial velocity of clusters.

The measurement of the integrated optical polarization of weakly gravitationally lensed galaxies can provide considerable constraints on lens models. The method outlined depends on fact that the orientation of the direction of optical polarization is not affected by weak gravitational lensing. The angle between the semi-major axis of the imaged galaxy and the direction of integrated optical polarization thus informs one of the distortion produced by the gravitational lensing. Although the method depends on the polarimetric measurement of faint galaxies, large telescopes and improved techniques should make such measurements possible in the near future.

We compare the mass functions obtained analytically, in the framework of an extended Press & Schechter (PS) formalism, in a previous paper to the numerical mass functions obtained in N -- body simulations, using different algorithms to define objects in the density field. After discussing the properties of the algorithms, we show that the mass function obtained using the friend -- of -- friend algorithm reproduces best the scaling behaviors predicted in the extended PS formalism. Following this statistical analysis, we show that it is possible in the framework of our extended PS formalism to reproduce the mass function but also, and for the first time, the initial statistical properties of structures and their collapse time. This allow to present a ``coherent'' picture of structure formation which can account for the initial, final and dynamical properties of structures.

We investigate some modifications to the Press & Schechter (1974) (PS) prescription resulting from shear and tidal effects. These modifications rely on more realistic treatments of the collapse process than the standard approach based on the spherical model. First, we show that the mass function resulting from a new approximate Lagrangian dynamic (Audit & Alimi 96), contains more objects at high mass, than the classical PS mass function and is well fitted by a PS-like function with a threshold density of $\delta_c \simeq 1.4$. However, such a Lagrangian description can underestimate the epoch of structure formation since it defines it as the collapse of the first principal axis. We therefore suggest some analytical prescriptions, for computing the collapse time along the second and third principal axes, and we deduce the corresponding mass functions. The collapse along the third axis is delayed by the shear and the number of objects of high mass then decreases. Finally, we show that the shear also strongly affects the formation of low-mass halos. This dynamical effect implies a modification of the low-mass slope of the mass function and allows the reproduction of the observed luminosity function of field galaxies. In a companion paper, we present results of numerical simulations which complete this work.

In this paper we present a new local Lagrangian approximation to the gravitational dynamics of cold matter. We describe the dynamics of a Lagrangian fluid element through only one quantity, the deformation tensor. We show that this tensor is clearly suited to the study of gravitational dynamics and, moreover, that knowing its evolution is enough to completely describe a fluid element. Our approximate dynamics treat exactly the conservation of mass, of the velocity divergence and of the shear and is exact in the case of planar, cylindrical and spherical collapses. It also reproduces very accurately the evolution of all dynamical quantities for a very wide class of initial conditions as illustrated by a detailed comparison with the homogeneous ellipsoid model. Beside, we highlight, for the first time, the important dynamical role played by the Newtonian counterpart of the magnetic part of the Weyl tensor.