Source author record

Neil Vaytet

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

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph.IM

Catalog footprint

What is connected

7works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

Discs and outflows in the early phases of massive star formation: influence of magnetic fields and ambipolar diffusion

We study mass accretion and ejection in the vicinity of massive star forming cores using high-resolution (5 au) 3D AMR numerical simulations. We investigate the mechanisms at the origin of outflows and characterise the properties of the disc forming around massive protostars. We include both protostellar radiative feedback via PMS evolutionary tracks and magnetic ambipolar diffusion. We studied 3 different cases: purely hydrodynamical, ideal MHD, and ambipolar diffusion. In the resistive models, we investigate the effects the initial amplitude of both magnetic field and rotation have on the properties of the massive protostellar system. We use simple criteria to identify the outflow and disc material and follow their evolution as the central star accretes mass up to 20 solar mass. The outflow is completely different when magnetic fields are introduced, so that magnetic processes are the main driver of the outflow up to stellar masses of ~20 solar mass. The disc properties depend on the physics included. The disc formed in the ideal and resistive runs show opposite properties in terms of plasma beta and of magnetic fields topology. While the disc in the ideal case is dominated by the magnetic pressure and the toroidal magnetic fields, the one formed in the resistive runs is dominated by the thermal pressure and has essentially vertical magnetic fields in the inner regions (R<200 au). We find that magnetic processes dominate the early evolution of massive protostellar systems (<20 solar mass) and shapes the accretion/ejection as well as the disc formation. Ambipolar diffusion is mainly at work at disc scales and regulates its properties. Our finding for the outflow and disc properties are reminiscent of low-mass star formation, suggesting that accretion and ejection in young massive and low-mass protostars are regulated by the same physical processes at the early stages.

preprint2016arXiv

A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

preprint2015arXiv

Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).

preprint2015arXiv

Multigroup radiation hydrodynamics with flux-limited diffusion and adaptive mesh refinement

Radiative transfer plays a key role in the star formation process. Due to a high computational cost, radiation-hydrodynamics simulations performed up to now have mainly been carried out in the grey approximation. In recent years, multi-frequency radiation-hydrodynamics models have started to emerge, in an attempt to better account for the large variations of opacities as a function of frequency. We wish to develop an efficient multigroup algorithm for the adaptive mesh refinement code RAMSES which is suited to heavy proto-stellar collapse calculations. Due to prohibitive timestep constraints of an explicit radiative transfer method, we constructed a time-implicit solver based on a stabilised bi-conjugate gradient algorithm, and implemented it in RAMSES under the flux-limited diffusion approximation. We present a series of tests which demonstrate the high performance of our scheme in dealing with frequency-dependent radiation-hydrodynamic flows. We also present a preliminary simulation of a three-dimensional proto-stellar collapse using 20 frequency groups. Differences between grey and multigroup results are briefly discussed, and the large amount of information this new method brings us is also illustrated. We have implemented a multigroup flux-limited diffusion algorithm in the RAMSES code. The method performed well against standard radiation-hydrodynamics tests, and was also shown to be ripe for exploitation in the computational star formation context.

preprint2014arXiv

On the role of the H2 ortho:para ratio in gravitational collapse during star formation

Hydrogen molecules (H2) come in two forms in the interstellar medium, ortho- and para-hydrogen, corresponding to the two different spin configurations of the two hydrogen atoms. The relative abundances of the two flavours in the interstellar medium are still very uncertain, and this abundance ratio has a significant impact on the thermal properties of the gas. In the context of star formation, theoretical studies have recently adopted two different strategies when considering the ortho:para ratio (OPR) of H2 molecules; the first considers the OPR to be frozen at 3:1 while the second assumes that the species are in thermal equilibrium. As the OPR potentially affects the protostellar cores which form as a result of the gravitational collapse of a dense molecular cloud, the aim of this paper is to quantify precisely what role the choice of OPR plays in the properties and evolution of the cores. We used two different ideal gas equations of state for a hydrogen and helium mix in a radiation hydrodynamics code to simulate the collapse of a dense cloud and the formation of the first and second Larson cores; the first equation of state uses a fixed OPR of 3:1 while the second assumes thermal equilibrium. Simulations using an equilibrium ratio collapse faster at early times and show noticeable oscillations around hydrostatic equilibrium, to the point where the core expands for a short time right after its formation before resuming its contraction. In the case of a fixed 3:1 OPR, the core's evolution is a lot smoother. The OPR was however found to have little impact on the size, mass and radius of the two Larson cores. We conclude that if one is solely interested in the final properties of the cores when they are formed, it does not matter which OPR is used. On the other hand, if one's focus lies primarily in the evolution of the first core, the choice of OPR becomes important.

preprint2013arXiv

Spatio-kinematic modelling of Abell 65, a double-shelled planetary nebula with a binary central star

We present the first detailed spatio-kinematical analysis and modelling of the planetary nebula Abell 65, which is known to host a post-common envelope, binary, central star system. As such, this object is of great interest in studying the link between nebular morphology and central star binarity. [OIII]5007A and H-alpha+[NII]6584A longslit spectra and imagery of Abell 65 were obtained with the Manchester Echelle Spectrometer on the 2.1-m telescope at the San Pedro Martir Observatory (MES-SPM). Further [OIII]5007A longslit spectra were obtained with the Ultraviolet and Visual Echelle Spectrograph on the Very Large Telescope (VLT-UVES). These data were used to develop a spatio-kinematical model for the [OIII]5007A emission from Abell 65. A "best-fit" model was found by comparing synthetic spectra and images rendered from the model to the data. The model comprises an outer shell and an inner shell, with kinematical ages of 15000 +/- 5000 yr kpc^-1 and 8000 +/- 3000 yr kpc^-1, respectively. Both shells have peanut-shaped bipolar structures with symmetry axes at inclinations of 55 +/- 10 deg (to the line-of-sight) for the outer shell and 68 +/- 10 deg for the inner shell. The near-alignment between the nebular shells and the binary orbital inclination (of 68 +/- 2 deg) is strongly indicative that the binary is responsible for shaping the nebula. Abell 65 is one of a growing number of planetary nebulae (seven to date, including Abell 65 itself) for which observations and modelling support the shaping influence of a central binary.

preprint2012arXiv

Simulations of protostellar collapse using multigroup radiation hydrodynamics. I. The first collapse

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.