Non-ideal magnetohydrodynamic simulations of subcritical pre-stellar cores with non-equilibrium chemistry

A. Tritsis*, C. Federrath, K. Willacy, K. Tassis

*Corresponding author for this work

    Research output: Contribution to journalArticlepeer-review

    7 Citations (Scopus)


    Non-ideal magnetohydrodynamic (MHD) effects are thought to be gravity's closest ally in overcoming the support of magnetic fields and in forming stars. Here, we modify the publicly available version of the adaptive mesh refinement code flash (Fryxell et al. 2000; Dubey et al. 2008) to include a detailed treatment of non-ideal MHD and study such effects in collapsing pre-stellar cores. We implement two very extended non-equilibrium chemical networks, the largest of which is comprised of ∼300 species and includes a detailed description of deuterium chemistry. The ambipolar diffusion, Ohmic and Hall resistivities are then self-consistently calculated from the abundances of charged species. We present a series of 2D axisymmetric simulations where we vary the chemical model, cosmic ray ionization rate, and grain distribution. We benchmark our implementation against ideal MHD simulations and previously published results. We show that, at high densities (n H2> ∼106∼ cm-3), the ion that carries most of the perpendicular and parallel conductivities is not H3+ as was previously thought, but is instead D 3+.

    Original languageEnglish
    Pages (from-to)4420-4435
    Number of pages16
    JournalMonthly Notices of the Royal Astronomical Society
    Issue number3
    Publication statusPublished - 1 Mar 2022


    Dive into the research topics of 'Non-ideal magnetohydrodynamic simulations of subcritical pre-stellar cores with non-equilibrium chemistry'. Together they form a unique fingerprint.

    Cite this