Search

We present a geometric multigrid solver for the M1 model of radiative transfer without source terms. In radiative hydrodynamics applications, the radiative transfer needs to be solved implicitly because of the fast propagation speed of photons relative to the fluid velocity. The M1 model is hyperbolic and can be discretized with an HLL solver, whose time implicit integration can be done using a nonlinear Jacobi method. One can show that this iterative method always preserves the admissible states, such as positive radiative energy and reduced flux less than 1. To decrease the number of iterations required for the solver to converge, and therefore to decrease the computational cost, we propose a geometric multigrid algorithm. Unfortunately, this method is not able to preserve the admissible states. In order to preserve the admissible state states, we introduce a pseudo-time such that the solution of the problem on the coarse grid is the steady state of a differential equation in pseudo-time. We present preliminary results showing the decrease of the number of iterations and computational cost as a function of the number of multigrid levels used in the method. These results suggest that nonlinear multigrid methods can be used as a robust implicit solver for hyperbolic systems such as the M1 model.

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., Comment: accepted for publication in A&A

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., Comment: in press (accepted for publication in A&A on 18 Jan 2018)

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]., Comment: submitted to A&A

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., Comment: Accepted to MNRAS

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., Comment: Accepted to MNRAS

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., Comment: paper accepted to MNRAS

Nowadays, climate models rely on couplers. Each complete climate model is broken into different sub-models (oceanic, atmospheric,...), each one working on a different grid. The coupler brings these models together and interpolates the physical quantities between the grids. However, neither the coupler nor sometimes the sub-models themselves know precisely the grid cell edges. They only know the grid cell corners (vertices) and the true grid cell areas. Thus, the coupler has to make assumptions about the grid cell edges in order to compute the grid cell intersections. For first-order schemes, the most straightforward way to interpolate scalar quantities is to directly use these approximate grid cell intersections, that don' take the true grid cell areas into account. It is the method used in the "conservative" regridding option implemented in the widely used spherical interpolation package SCRIP. We show that is doesn't preserve integrals in the general case, whether the coupler using the SCRIP-generated weights transmits directly the intensive quantity, or its extensive counterpart (that is the same quantity multiplied by the true area of the individual source grid cell). We show how to modify the interpolation scheme to preserve integrals in the general case., Comment: 16 pages, 1 figure

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., Comment: 9 pages, 5 figures, accepted for publication in A&A

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., Comment: 15 pages, 11 figures, accepted for publication in A&A

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., Comment: 16 pages, 9 figures, A&A accepted

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., Comment: 4 pages, 2 figures ; to appear in "IAU Symposium 270: Computational Star formation", Eds. J. Alves, B. Elmegreen, J. Girart, V. Trimble

We have studied the dynamics of strong radiative shocks generated with the high-energy subnanosecond iodine laser at Prague Asterix Laser System facility, Comment: with small correction in Fig.12

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., Comment: accepted to ApJ

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., Comment: 21 pages, 1 table, 7 figures

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., Comment: Accepted to MNRAS main journal. Additional discussion and and one figure with error estimates is added. 7 pages, 7 figures. Main conclusions unchanged

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., Comment: 8 pages, 7 figures, accepted for publication in A&A

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, Comment: 16 pages, 16 figures, accepted for publication in A&A

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., Comment: Submitted to MNRAS letter. 5 pages using mnras file style. email: edouard.audit@obspm.fr

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., Comment: 13 pages, 11 figures, uses mnras style file. Accepted for publication in MNRAS

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., Comment: A&A in press, 13 pages, 8 Postscript figures

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., Comment: 12 pages, latex 2e, accepted for publication in Astronomy & Astrophysics

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., Comment: 12 pages,l-aa.sty macro. latex with postscript figures. Accepted in A&A. Send Preprint request to edouard@gin.obspm.fr