-
Moving-mesh non-ideal magnetohydrodynamical simulations of the collapse of cloud cores to protostars
Authors:
Alexander C. Mayer,
Oliver Zier,
Thorsten Naab,
Rüdiger Pakmor,
Paola Caselli,
Alexei V. Ivlev,
Volker Springel,
Stefanie Walch
Abstract:
Magnetic fields have been shown both observationally and through theoretical work to be an important factor in the formation of protostars and their accretion disks. Accurate modelling of the evolution of the magnetic field in low-ionization molecular cloud cores requires the inclusion of non-ideal magnetohydrodynamics (MHD) processes, specifically Ohmic and ambipolar diffusion and the Hall effect…
▽ More
Magnetic fields have been shown both observationally and through theoretical work to be an important factor in the formation of protostars and their accretion disks. Accurate modelling of the evolution of the magnetic field in low-ionization molecular cloud cores requires the inclusion of non-ideal magnetohydrodynamics (MHD) processes, specifically Ohmic and ambipolar diffusion and the Hall effect. These have a profound influence on the efficiency of magnetic removal of angular momentum from protostellar disks and simulations that include them can avoid the `magnetic-braking catastrophe' in which disks are not able to form. However, the impact of the Hall effect, in particular, is complex and remains poorly studied. In this work, we perform a large suite of simulations of the collapse of cloud cores to protostars with several non-ideal MHD chemistry models and initial core geometries using the moving-mesh code {\small AREPO}. We find that the efficiency of angular momentum removal is significantly reduced with respect to ideal MHD, in line with previous results. The Hall effect has a varied influence on the evolution of the disk which depends on the initial orientation of the magnetic field. This extends to the outflows seen in a subset of the models, where this effect can act to enhance or suppress them and open up new outflow channels. We conclude, in agreement with a subset of the previous literature, that the Hall effect is the dominant non-ideal MHD process in some collapse scenarios and thus should be included in simulations of protostellar disk formation.
△ Less
Submitted 9 October, 2024;
originally announced October 2024.
-
Non-ideal magnetohydrodynamics on a moving mesh II: Hall effect
Authors:
Oliver Zier,
Alexander C. Mayer,
Volker Springel
Abstract:
In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square reconstruction on the unstructured mesh, which we also used in the companion paper for Ohmic and ambipolar diffusion. In an extensive study of simulations of a magnetic…
▽ More
In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square reconstruction on the unstructured mesh, which we also used in the companion paper for Ohmic and ambipolar diffusion. In an extensive study of simulations of a magnetic shock, we show that without additional magnetic diffusion our algorithm for the Hall effect becomes unstable at high resolution. We can however stabilise it by artificially increasing the Ohmic resistivity, $η_{\rm OR}$, so that it satisfies the condition $η_{\rm OR} \geq η_{\rm H} /5$, where $η_{\rm H}$ is the Hall diffusion coefficient. Adopting this solution we find second order convergence for the C-shock and are also able to accurately reproduce the dispersion relation of the whistler waves. As a first application of the new scheme, we simulate the collapse of a magnetised cloud with constant Hall parameter $η_{\rm H}$ and show that, depending on the sign of $η_{\rm H}$, the magnetic braking can either be weakened or strengthened by the Hall effect. The quasi-Lagrangian nature of the moving mesh method used here automatically increases the resolution in the forming core, making it well suited for more realistic studies with non-constant magnetic diffusivities in the future.
△ Less
Submitted 27 September, 2023;
originally announced September 2023.
-
Non-ideal magnetohydrodynamics on a moving mesh I: Ohmic and ambipolar diffusion
Authors:
Oliver Zier,
Volker Springel,
Alexander C. Mayer
Abstract:
Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we present a new numerical scheme for the first two resistive terms, which we implement in the moving-mesh code AREPO using the single-fluid approximation combine…
▽ More
Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we present a new numerical scheme for the first two resistive terms, which we implement in the moving-mesh code AREPO using the single-fluid approximation combined with a new gradient estimation technique based on a least-squares fit per interface. Through various test calculations including the diffusion of a magnetic peak, the structure of a magnetic C-shock, and the damping of an Alfvén wave, we show that we can achieve an accuracy comparable to the state-of-the-art code ATHENA++. We apply the scheme to the linear growth of the magnetorotational instability and find good agreement with the analytical growth rates. By simulating the collapse of a magnetised cloud with constant magnetic diffusion, we show that the new scheme is stable even for large density contrasts. Thanks to the Lagrangian nature of the moving mesh method the new scheme is thus well suited for intended future applications where a high resolution in the dense cores of collapsing protostellar clouds needs to be achieved. In a forthcoming work we will extend the scheme to the Hall effect.
△ Less
Submitted 21 July, 2023;
originally announced July 2023.
-
ΛCDM with baryons vs. MOND: the time evolution of the universal acceleration scale in the Magneticum simulations
Authors:
Alexander C. Mayer,
Adelheid F. Teklu,
Klaus Dolag,
Rhea-Silvia Remus
Abstract:
MOdified Newtonian Dynamics (MOND) is an alternative to the standard Cold Dark Matter (CDM) paradigm which proposes an alteration of Newton's laws of motion at low accelerations, characterized by a universal acceleration scale a_0. It attempts to explain observations of galactic rotation curves and predicts a specific scaling relation of the baryonic and total acceleration in galaxies, referred to…
▽ More
MOdified Newtonian Dynamics (MOND) is an alternative to the standard Cold Dark Matter (CDM) paradigm which proposes an alteration of Newton's laws of motion at low accelerations, characterized by a universal acceleration scale a_0. It attempts to explain observations of galactic rotation curves and predicts a specific scaling relation of the baryonic and total acceleration in galaxies, referred to as the Rotational Acceleration Relation (RAR), which can be equivalently formulated as a Mass Discrepancy Acceleration Relation (MDAR). The appearance of these relations in observational data such as SPARC has lead to investigations into the existence of similar relations in cosmological simulations using the standard ΛCDM model. Here, we report the existence of an RAR and MDAR similar to that predicted by MOND in ΛCDM using a large sample of galaxies extracted from a cosmological, hydrodynamical simulation (Magneticum). Furthermore, by using galaxies in Magneticum at different redshifts, a prediction for the evolution of the inferred acceleration parameter a_0 with cosmic time is derived by fitting a MOND force law to these galaxies. In Magneticum, the best fit for a_0 is found to increase by a factor of approximately 3 from redshift z = 0 to z = 2. This offers a powerful test from cosmological simulations to distinguish between MOND and ΛCDM observationally.
△ Less
Submitted 17 October, 2022; v1 submitted 9 June, 2022;
originally announced June 2022.
-
Growth of macroscopic-area single crystal polyacene thin films on arbitrary substrates
Authors:
Randall L. Headrick,
Hua Zhou,
Binran Wang,
Yiping Wang,
Greggory P. Carpenter,
Alex C. Mayer,
Matthew Lloyd,
George G. Malliaras,
Alexander Kazimirov,
John E. Anthony
Abstract:
Organic electronic materials have potential applications in a number of low-cost, large area electronic devices such as flat panel displays and inexpensive solar panels. Small molecules in the series Anthracene, Tetracene, Pentacene, are model molecules for organic semiconductor thin films to be used as the active layers in such devices. This has motivated a number of studies of polyacene thin fil…
▽ More
Organic electronic materials have potential applications in a number of low-cost, large area electronic devices such as flat panel displays and inexpensive solar panels. Small molecules in the series Anthracene, Tetracene, Pentacene, are model molecules for organic semiconductor thin films to be used as the active layers in such devices. This has motivated a number of studies of polyacene thin film growth and structure. Although the majority of these studies rely on vapor-deposited films, solvent-based deposition of films with improved properties onto non-crystalline substrates is desired for industrial production of devices. Improved ordering in thin films will have a large impact on their electronic properties, since grain boundaries and other defects are detrimental to carrier mobilities and lifetimes. Thus, a long-standing challenge in this field is to prepare largearea single crystal films on arbitrary substrates. Here we demonstrate a solvent-based method to deposit thin films of organic semiconductors, which is very general. Anthracene thin films with single-crystal domain sizes exceeding 1x1 cm2 can be prepared on various substrates by the technique. Films of 6,13-bis(triisopropylsilylethynyl)pentacene are also demonstrated to have grain sizes larger than 2x2 mm2. In contrast, films produced by conventional means such as vapor deposition or spin coating are polycrystalline with micron-scale grain sizes. The general propensity of these small molecules towards crystalline order in an optimized solvent deposition process shows that there is great potential for thin films with improved properties. Films prepared by these methods will also be useful in exploring the limits of performance in organic thin film devices.
△ Less
Submitted 17 July, 2011;
originally announced July 2011.