Source author record

P. Lesaffre

P. Lesaffre appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

17 published item(s)

preprint2022arXiv

BxC: a swift generator for 3D magnetohydrodynamic turbulence

Magnetohydrodynamic turbulence is central to laboratory and astrophysical plasmas, and is invoked for interpreting many observed scalings. Verifying predicted scaling law behaviour requires extreme-resolution direct numerical simulations (DNS), with needed computing resources excluding systematic parameter surveys. We here present an analytic generator of realistically looking turbulent magnetic fields, that computes 3D ${\cal{O}}(1000^3)$ solenoidal vector fields in minutes to hours on desktops. Our model is inspired by recent developments in 3D incompressible fluid turbulence theory, where a Gaussian white noise vector subjected to a non-linear transformation results in an intermittent, multifractal random field. Our $B\times C$ model has only few parameters that have clear geometric interpretations. We directly compare a (costly) DNS with a swiftly $B\times C$-generated realization, in terms of its (i) characteristic sheet-like structures of current density, (ii) volume-filling aspects across current intensity, (iii) power-spectral behaviour, (iv) probability distribution functions of increments for magnetic field and current density, structure functions, spectra of exponents, and (v) partial variance of increments. The model even allows to mimic time-evolving magnetic and current density distributions and can be used for synthetic observations on 3D turbulent data cubes.

preprint2021arXiv

A method to statistically characterize turbulent data with physically motivated parameters, illustrated on a centroid velocity map

We investigate the potential of a recently proposed model for 3D compressible MHD turbulence (Chevillard et al. 2010; Durrive et al. 2021) to be used as a tool to characterize statistically 2D and 3D turbulent data. This model is parametrized by a dozen of free (intuitive, physically motivated) parameters, which control the statistics of the fields (density, velocity and magnetic fields). The present study is a proof of concept study: (i) we restrict ourselves to the incompressible hydrodynamical part of the model, (ii) we consider as data centroid velocity maps, and (iii) we let only three of the free parameters vary (namely the correlation length, the Hurst parameter and the intermittency parameter). Within this framework, we demonstrate that, given a centroid velocity map, we can find in an automated manner (i.e. by a Markov Chain Monte Carlo analysis) values of the parameters such that the model resembles the given map, i.e. which reproduces its statistics fairly well. Hence, thanks to this procedure, one may characterize statistically, and thus compare, various turbulent data. In other words, we show how this model may be used as a metric to compare observational or simulated data sets. In addition, because this model is numerically particularly fast (nearly 500 times faster than the numerical simulation we use to generate our reference data) it may be used as a surrogate model. Finally, by this process we also initiate the first systematic exploration of the parameter space of this model. Doing so, we show how the parameters impact the visual and the statistical properties of centroid velocity maps, and exhibit the correlations between the various parameters, providing new insight into the model.

preprint2021arXiv

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses.

preprint2021arXiv

An intuitive parametric model for 3D compressible hydrodynamical and MHD turbulence

An analytical model for three-dimensional incompressible turbulence was recently introduced in the hydrodynamics community which, with only a few parameters, shares many properties of experimental and numerical turbulence, notably intermittency (non-Gaussianity), the energy cascade (skewness), and vorticity alignment properties. In view of modeling astrophysical environments, we introduce a manner to extend to compressible fluids the three-dimensional turbulent velocity field model of Chevillard et al. (2010), as well as the three 3D turbulent magnetic field models of Durrive et al. (2020), following the same procedure based on the concept of multiplicative chaos. Our model provides a complementary tool to numerical simulations, as it enables us to generate very quickly fairly realistic velocity fields and magnetic fields, the statistics of which are controllable with intuitive parameters. Therefore our model will also provide a useful tool for observers in astrophysics. Finally, maybe even more than the model itself, it is the very procedure that matters the most: our model is modular, in the sense that it is constructed gradually, with intuitive and physically motivated steps, so that it is prone to many further improvements.

preprint2021arXiv

Molecules, shocks, and disk in the axi-symmetric wind of the MS-type AGB star RS Cancri

The latest evolutionary phases of low- and intermediate mass stars are characterized by complex physical processes like turbulence, convection, stellar pulsations, magnetic fields, condensation of solid particles, and the formation of massive outflows that inject freshly produced heavy elements and dust particles into the interstellar medium. We use the Northern Extended Millimeter Array (NOEMA) to obtain spatially and spectrally resolved observations of the semi-regular Asymptotic Giant Branch star RS Cancri to shed light on the morpho-kinematic structure of its inner, wind forming environment by applying detailed 3-D reconstruction modeling and LTE radiative transfer calculations. We detect 32 lines of 13 molecules and isotopologs (CO, SiO, SO, SO$_2$, H$_2$O, HCN, PN), including several transitions from vibrationally excited states. HCN, H$^{13}$CN, millimeter vibrationally excited H$_2$O, SO, $^{34}$SO, SO$_2$, and PN are detected for the first time in RS Cnc. Evidence for rotation is seen in HCN, SO, SO$_2$, and SiO(v=1). From CO and SiO channel maps, we find an inner, equatorial density enhancement, and a bipolar outflow structure with a mass loss rate of $1 \times 10^{-7}M_\odot{\rm yr}^{-1}$ for the equatorial region and of $2 \times 10^{-7}M_\odot{\rm yr}^{-1}$ for the polar outflows. The $^{12}$CO/$^{13}$CO ratio is measured to be $\sim20$ on average, $24\pm2$ in the polar outflows and $19\pm3$ in the equatorial region. We do not find direct evidence of a companion that might explain this kind of kinematic structure, and explore the possibility that a magnetic field might be the cause of it. The innermost molecular gas is influenced by stellar pulsation and possibly by convective cells that leave their imprint on broad wings of certain molecular lines, such as SiO and SO.

preprint2020arXiv

Magnetic fields from Multiplicative Chaos

An analytical model for fully developed three-dimensional incompressible turbulence was recently proposed in the hydrodynamics community, based on the concept of multiplicative chaos. It consists of a random field represented by means of a stochastic integral, which, with only a few parameters, shares many properties with experimental and numerical turbulence, including in particular energy transfer through scales (the cascade) and intermittency (non-Gaussianity) which is most conveniently controlled with a single parameter. Here, we propose three models extending this approach to MHD turbulence. Our formulae provide physically motivated 3D models of a turbulent velocity field and magnetic field coupled together. Besides its theoretical value, this work is meant to provide a tool for observers: a dozen of physically meaningful free parameters enter the description, which is useful to characterize astrophysical data.

preprint2019arXiv

Probing the cold magnetized Universe with SPICA-POL (B-BOP)

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

preprint2016arXiv

Tracing extended low-velocity shocks through SiO emission - Case study of the W43-MM1 ridge

Previous literature suggests that the densest structures in the interstellar medium form through colliding flows, but patent evidence of this process is still missing. Recent literature proposes using SiO line emission to trace low-velocity shocks associated with cloud formation through collision. In this paper we investigate the bright and extended SiO(2-1) emission observed along the ~5 pc-long W43-MM1 ridge to determine its origin. We used high angular resolution images of the SiO(2-1) and HCN(1-0) emission lines obtained with the IRAM plateau de Bure (PdBI) interferometer and combined with data from the IRAM 30 m radiotelescope. These data were complemented by a Herschel column density map of the region. We performed spectral analysis of SiO and HCN emission line profiles to identify protostellar outflows and spatially disentangle two velocity components associated with low- and high-velocity shocks. Then, we compared the low-velocity shock component to a dedicated grid of one-dimensional (1D) radiative shock models. We find that the SiO emission originates from a mixture of high-velocity shocks caused by bipolar outflows and low-velocity shocks. Using SiO and HCN emission lines, we extract seven bipolar outflows associated with massive dense cores previously identified within the W43-MM1 mini-starburst cluster. Comparing observations with dedicated Paris-Durham shock models constrains the velocity of the low-velocity shock component from 7 to 12km/s. The SiO arising from low-velocity shocks spreads along the complete length of the ridge. Its contribution represents at least 45% and up to 100% of the total SiO emission depending on the area considered. The low-velocity component of SiO is most likely associated with the ridge formation through colliding flows or cloud-cloud collision.

preprint2015arXiv

Sulphur-bearing molecules in diffuse molecular clouds: new results from SOFIA/GREAT and the IRAM 30 m telescope

We have observed five sulphur-bearing molecules in foreground diffuse molecular clouds lying along the sight-lines to five bright continuum sources. We have used the GREAT instrument on SOFIA to observe the 1383 GHz $^2Π_{3/2} J=5/2-3/2$ transitions of SH towards the star-forming regions W31C, G29.96-0.02, G34.3+0.1, W49N and W51, detecting foreground absorption towards all five sources; and the EMIR receivers on the IRAM 30m telescope at Pico Veleta to detect the H$_2$S 1(10)-1(01), CS J=2-1 and SO 3(2)-2(1) transitions. In nine foreground absorption components detected towards these sources, the inferred column densities of the four detected molecules showed relatively constant ratios, with N(SH)/N(H$_2$S) in the range 1.1 - 3.0, N(CS)/N(H$_2$S) in the range 0.32 - 0.61, and N(SO)/N(H$_2$S) in the range 0.08 - 0.30. The observed SH/H$_2$ ratios - in the range (0.5-2.6) $\times 10^{-8}$ - indicate that SH (and other sulphur-bearing molecules) account for << 1% of the gas-phase sulphur nuclei. The observed abundances of sulphur-bearing molecules, however, greatly exceed those predicted by standard models of cold diffuse molecular clouds, providing further evidence for the enhancement of endothermic reaction rates by elevated temperatures or ion-neutral drift. We have considered the observed abundance ratios in the context of shock and turbulent dissipation region (TDR) models. Using the TDR model, we find that the turbulent energy available at large scale in the diffuse ISM is sufficient to explain the observed column densities of SH and CS. Standard shock and TDR models, however, fail to reproduce the column densities of H$_2$S and SO by a factor of about 10; more elaborate shock models - in which account is taken of the velocity drift, relative to H$_2$, of SH molecules produced by the dissociative recombination of H$_3$S$^+$ - reduce this discrepancy to a factor ~ 3.

preprint2014arXiv

Dense molecular globulettes and the dust arc towards the runaway O star AE Aur (HD 34078)

Some runaway stars are known to display IR arc-like structures around them, resulting from their interaction with surrounding interstellar material. The properties of these features as well as the processes involved in their formation are still poorly understood. We aim at understanding the physical mechanisms that shapes the dust arc observed near the runaway O star AEAur (HD34078). We obtained and analyzed a high spatial resolution map of the CO(1-0) emission that is centered on HD34078, and that combines data from both the IRAM interferometer and 30m single-dish antenna. The line of sight towards HD34078 intersects the outer part of one of the detected globulettes, which accounts for both the properties of diffuse UV light observed in the field and the numerous molecular absorption lines detected in HD34078's spectra, including those from highly excited H2 . Their modeled distance from the star is compatible with the fact that they lie on the 3D paraboloid which fits the arc detected in the 24 μm Spitzer image. Four other compact CO globulettes are detected in the mapped area. These globulettes have a high density and linewidth, and are strongly pressure-confined or transient. The good spatial correlation between the CO globulettes and the IR arc suggests that they result from the interaction of the radiation and wind emitted by HD 34078 with the ambient gas. However, the details of this interaction remain unclear. A wind mass loss rate significantly larger than the value inferred from UV lines is favored by the large IR arc size, but does not easily explain the low velocity of the CO globulettes. The effect of radiation pressure on dust grains also meets several issues in explaining the observations. Further observational and theoretical work is needed to fully elucidate the processes shaping the gas and dust in bow shocks around runaway O stars. (Abridged)

preprint2013arXiv

A two-dimensional mixing length theory of convective transport

The helioseismic observations of the internal rotation profile of the Sun raise questions about the two-dimensional (2D) nature of the transport of angular momentum in stars. Here we derive a convective prescription for axisymmetric (2D) stellar evolution models. We describe the small scale motions by a spectrum of unstable linear modes in a Boussinesq fluid. Our saturation prescription makes use of the angular dependence of the linear dispersion relation to estimate the anisotropy of convective velocities. We are then able to provide closed form expressions for the thermal and angular momentum fluxes with only one free parameter, the mixing length. We illustrate our prescription for slow rotation, to first order in the rotation rate. In this limit, the thermodynamical variables are spherically symetric, while the angular momentum depends both on radius and latitude. We obtain a closed set of equations for stellar evolution, with a self-consistent description for the transport of angular momentum in convective regions. We derive the linear coefficients which link the angular momentum flux to the rotation rate ($Λ$- effect) and its gradient ($α$-effect). We compare our results to former relevant numerical work.

preprint2013arXiv

Low-velocity shocks traced by extended SiO emission along the W43 ridges: witnessing the formation of young massive clusters

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

preprint2010arXiv

Uncertainties and robustness of the ignition process in type Ia supernovae

It is widely accepted that the onset of the explosive carbon burning in the core of a CO WD triggers the ignition of a SN Ia. The features of the ignition are among the few free parameters of the SN Ia explosion theory. We explore the role for the ignition process of two different issues: firstly, the ignition is studied in WD models coming from different accretion histories. Secondly, we estimate how a different reaction rate for C-burning can affect the ignition. Two-dimensional hydrodynamical simulations of temperature perturbations in the WD core ("bubbles") are performed with the FLASH code. In order to evaluate the impact of the C-burning reaction rate on the WD model, the evolution code FLASH_THE_TORTOISE from Lesaffre et al. (2006) is used. In different WD models a key role is played by the different gravitational acceleration in the progenitor's core. As a consequence, the ignition is disfavored at a large distance from the WD center in models with a larger central density, resulting from the evolution of initially more massive progenitors. Changes in the C reaction rate at T < 5e8 K slightly influence the ignition density in the WD core, while the ignition temperature is almost unaffected. Recent measurements of new resonances in the C-burning reaction rate (Spillane et al. 2007) do not affect the core conditions of the WD significantly. This simple analysis, performed on the features of the temperature perturbations in the WD core, should be extended in the framework of the state-of-the-art numerical tools for studying the turbulent convection and ignition in the WD core. Future measurements of the C-burning reactions cross section at low energy, though certainly useful, are not expected to affect dramatically our current understanding of the ignition process.

preprint2007arXiv

Time-dependent analytic solutions of quasi-steady shocks with cooling

I present time-dependent analytical solutions of quasi-steady shocks with cooling, where quasi-steady shocks are objects composed of truncated steady-state models of shocks at any intermediate time. I compare these solutions to simulations with a hydrodynamical code and finally discuss quasi-steady shocks as approximations to time-dependent shocks. Large departure of both the adiabatic and steady-state approximations from the quasi-steady solution emphasise the importance of the cooling history in determining the trajectory of a shock.

preprint2006arXiv

The C-flash and the ignition conditions of type Ia supernovae

Thanks to a stellar evolution code able to compute through the C-flash we link the binary population synthesis of single degenerate progenitors of type Ia supernovae (SNe Ia) to their physical condition at the time of ignition. We show that there is a large range of possible ignition densities and we detail how their probability distribution depends on the accretion properties. The low density peak of this distribution qualitatively reminds of the clustering of the luminosities of Branch-normal SNe Ia. We tighten the possible range of initial physical conditions for explosion models: they form a one-parameter family, independent of the metallicity. We discuss how these results may be modified if we were to relax our hypothesis of a permanent Hachisu wind or if we were to include electron captures.

preprint2004arXiv

The convective Urca process

One possible fate of an accreting white dwarf is explosion in a type Ia supernova. However, the route to the thermonuclear runaway has always been uncertain owing to the lack of a convective model consistent with the Urca process. We derive a formalism for convective motions involving two radial flows. This formalism provides a framework for convective models that guarantees self-consistency for chemistry and energy budget, allows time-dependence and describes the interaction of convective motions with the global contraction or expansion of the star. In the one-stream limit, we reproduce several already existing convective models and allow them to treat chemistry. We also suggest as a model easy to implement in a stellar evolution code. We apply this formalism to convective Urca cores in Chandrasekhar mass white dwarfs. We stress that in degenerate matter, nuclear reactions that change the number of electrons strongly influence the convective velocities. We point out the sensitivity of the energy budget on the mixing. We illustrate our model by computing {\it stationary} convective cores with Urca nuclei. We show that even a very small mass fraction of Urca nuclei ($10^{-8}$) strongly influences the convective velocities. Finally, we present preliminary computations of the late evolution of a close to Chandrasekhar mass C+O white dwarf including the convective Urca process.