Source author record

AAke Nordlund

AAke Nordlund 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

15works
11topics
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

15 published item(s)

preprint2016arXiv

Microscopic Processes in Global Relativistic Jets Containing Helical Magnetic Fields

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

preprint2016arXiv

Supernova Driving. II. Compressive Ratio in Molecular-Cloud Turbulence

The compressibility of molecular cloud (MC) turbulence plays a crucial role in star formation models, because it controls the amplitude and distribution of density fluctuations. The relation between the compressive ratio (the ratio of powers in compressive and solenoidal motions) and the statistics of turbulence has been previously studied systematically only in idealized simulations with random external forces. In this work, we analyze a simulation of large-scale turbulence (250 pc) driven by supernova (SN) explosions that has been shown to yield realistic MC properties. We demonstrate that SN driving results in MC turbulence with a broad lognormal distribution of the compressive ratio, with a mean value $\approx 0.3$, lower than the equilibrium value of $\approx 0.5$ found in the inertial range of isothermal simulations with random solenoidal driving. We also find that the compressibility of the turbulence is not noticeably affected by gravity, nor are the mean cloud radial (expansion or contraction) and solid-body rotation velocities. Furthermore, the clouds follow a general relation between the rms density and the rms Mach number similar to that of supersonic isothermal turbulence, though with a large scatter, and their average gas density PDF is described well by a lognormal distribution, with the addition of a high-density power-law tail when self-gravity is included.

preprint2016arXiv

Tracking the distribution of $^{26}$Al and $^{60}$Fe during the early phases of star and disk evolution

The short-lived $^{26}$Al and $^{60}$Fe radionuclides are synthesized and expelled in the interstellar medium by core-collapse supernova events. The solar system's first solids, calcium-aluminium refractory inclusions (CAIs), contain evidence for the former presence of the $^{26}$ Al nuclide defining the canonical $^{26}$Al/$^{27}$ Al ratio of $\sim5 \times10^{-5}$. A different class of objects temporally related to canonical CAIs are CAIs with fractionation and unidentified nuclear effects (FUN CAIs), which record a low initial $^{26}$Al/$^{27}$Al of $10^{-6}$. The contrasting level of $^{26}$Al between these objects is often interpreted as reflecting the admixing of the $^{26}$Al nuclide during the early formative phase of the Sun. We use giant molecular cloud (GMC) scale adaptive mesh-refinement numerical simulations to trace the abundance of $^{26}$Al and $^{60}$Fe in star-forming gas during the early stages of accretion of individual low mass protostars. We find that the $^{26}$Al/$^{27}$Al and $^{60}$Fe/$^{56}$Fe ratios of accreting gas within a vicinity of 1000 AU of the stars follow the predicted decay curves of the initial abundances at time of star formation without evidence of spatial or temporal heterogeneities for the first 100 kyr of star formation. Therefore, the observed differences in $^{26}$Al/$^{27}$Al ratios between FUN and canonical CAIs are likely not caused by admixing of supernova material during the early evolution of the proto-Sun. Selective thermal processing of dust grains is a more viable scenario to account for the heterogeneity in $^{26}$Al/$^{27}$Al ratios at the time of solar system formation.

preprint2013arXiv

Abundance of 26Al and 60Fe in evolving Giant Molecular Clouds

The nucleosynthesis and ejection of radioactive $^{26}$Al (t$_{1/2} \sim$ 0.72\,Myr) and $^{60}$Fe, (t$_{1/2} \sim$ 2.5\,Myr) into the interstellar medium is dominated by the stellar winds of massive stars and supernova type II explosions. Studies of meteorites and their components indicate that the initial abundances of these short-lived radionuclides in the solar protoplanetary disk were higher than the background levels of the galaxy inferred from $γ$--ray astronomy and models of the galactic chemical evolution. This observation has been used to argue for a late-stage addition of stellar debris to the Solar System's parental molecular cloud or, alternatively, the solar protoplanetary disk, thereby requiring a special scenario for the formation of our Solar System. Here, we use supercomputers to model---from first principles---the production, transport and admixing of freshly synthesized $^{26}$Al and $^{60}$Fe in star-forming regions within giant molecular clouds. Under typical star-formation conditions, the levels of $^{26}$Al in most star-forming regions are comparable to that deduced from meteorites, suggesting that the presence of short-lived radionuclides in the early Solar System is a generic feature of the chemical evolution of giant molecular clouds. The $^{60}$Fe/$^{26}$Al yield ratio of $\approx 0.2$ calculated from our simulations is consistent with the galactic value of 0.15 $\pm$ 0.06 inferred from $γ$--ray astronomy but is significantly higher than most current Solar System measurements indicate. We suggest that estimates based on differentiated meteorites and some chondritic components may not be representative of the initial $^{60}$Fe abundance of the bulk Solar System.

preprint2012arXiv

A Simple Law of Star Formation

We show that supersonic MHD turbulence yields a star formation rate (SFR) as low as observed in molecular clouds (MCs), for characteristic values of the free-fall time divided by the dynamical time, $t_{\rm ff}/t_{\rm dyn}$, the alfvénic Mach number, ${\cal M}_{\rm a}$, and the sonic Mach number, ${\cal M}_{\rm s}$. Using a very large set of deep adaptive-mesh-refinement simulations, we quantify the dependence of the SFR per free-fall time, $ε_{\rm ff}$, on the above parameters. Our main results are: i) $ε_{\rm ff}$ decreases exponentially with increasing $t_{\rm ff}/t_{\rm dyn}$, but is insensitive to changes in ${\cal M}_{\rm s}$, for constant values of $t_{\rm ff}/t_{\rm dyn}$ and ${\cal M}_{\rm a}$. ii) Decreasing values of ${\cal M}_{\rm a}$ (stronger magnetic fields) reduce $ε_{\rm ff}$, but only to a point, beyond which $ε_{\rm ff}$ increases with a further decrease of ${\cal M}_{\rm a}$. iii) For values of ${\cal M}_{\rm a}$ characteristic of star-forming regions, $ε_{\rm ff}$ varies with ${\cal M}_{\rm a}$ by less than a factor of two. We propose a simple star-formation law, based on the empirical fit to the minimum $ε_{\rm ff}$, and depending only on $t_{\rm ff}/t_{\rm dyn}$: $ε_{\rm ff} \approx ε_{\rm wind} \exp(-1.6 \,t_{\rm ff}/t_{\rm dyn})$. Because it only depends on the mean gas density and rms velocity, this law is straightforward to implement in simulations and analytical models of galaxy formation and evolution.

preprint2012arXiv

Photon-Plasma: a modern high-order particle-in-cell code

We present the Photon-Plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks.

preprint2011arXiv

Astrophysical turbulence modeling

The role of turbulence in various astrophysical settings is reviewed. Among the differences to laboratory and atmospheric turbulence we highlight the ubiquitous presence of magnetic fields that are generally produced and maintained by dynamo action. The extreme temperature and density contrasts and stratifications are emphasized in connection with turbulence in the interstellar medium and in stars with outer convection zones, respectively. In many cases turbulence plays an essential role in facilitating enhanced transport of mass, momentum, energy, and magnetic fields in terms of the corresponding coarse-grained mean fields. Those transport properties are usually strongly modified by anisotropies and often completely new effects emerge in such a description that have no correspondence in terms of the original (non coarse-grained) fields.

preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

preprint2011arXiv

Radiation from sub-Larmor scale magnetic fields

Spontaneous rapid growth of strong magnetic fields is rather ubiquitous in high-energy density environments ranging from astrophysical sources (e.g., gamma-ray bursts and relativistic shocks), to reconnection, to laser-plasma interaction laboratory experiments, where they are produced by kinetic streaming instabilities of the Weibel type. Relativistic electrons propagating through these sub-Larmor-scale magnetic fields radiate in the jitter regime, in which the anisotropy of the magnetic fields and the particle distribution have a strong effect on the produced radiation. Here we develop the general theory of jitter radiation, which includes (i) anisotropic magnetic fields and electron velocity distributions, (ii) the effects of trapped electrons and (iii) extends the description to large deflection angles of radiating particles thus establishing a cross-over between the classical jitter and synchrotron regimes. Our results are in remarkable agreement with the radiation spectra obtained from particle-in-cell simulations of the classical Weibel instability. Particularly interesting is the onset of the field growth, when the transient hard synchrotron-violating spectra are common as a result of the dominant role of the trapped population. This effect can serve as a distinct observational signature of the violent field growth in astrophysical sources and lab experiments. It is also interesting that a system with small-scale fields tends to evolve toward the small-angle jitter regime, which can, under certain conditions, dominate the overall emission of a source.

preprint2011arXiv

The observable prestellar phase of the IMF

The observed similarities between the mass function of prestellar cores (CMF) and the stellar initial mass function (IMF) have led to the suggestion that the IMF is already largely determined in the gas phase. However, theoretical arguments show that the CMF may differ significantly from the IMF. In this Letter, we study the relation between the CMF and the IMF, as predicted by the IMF model of Padoan and Nordlund. We show that 1) the observed mass of prestellar cores is on average a few times smaller than that of the stellar systems they generate; 2) the CMF rises monotonically with decreasing mass, with a noticeable change in slope at approximately 3-5 solar masses, depending on mean density; 3) the selection of cores with masses larger than half their Bonnor-Ebert mass yields a CMF approximately consistent with the system IMF, rescaled in mass by the same factor as our model IMF, and therefore suitable to estimate the local efficiency of star formation, and to study the dependence of the IMF peak on cloud properties; 4) only one in five pre-brown-dwarf core candidates is a true progenitor to a brown dwarf.

preprint2010arXiv

Radiation Spectral Synthesis of Relativistic Filamentation

Radiation from many astrophysical sources, e.g. gamma-ray bursts and active galactic nuclei, is believed to arise from relativistically shocked collisionless plasmas. Such sources often exhibit highly transient spectra evolving rapidly, compared with source lifetimes. Radiation emitted from these sources is typically associated with non-linear plasma physics, complex field topologies and non-thermal particle distributions. In such circumstances a standard synchrotron paradigm may fail to produce accurate conclusions regarding the underlying physics. Simulating spectral emission and spectral evolution numerically in various relativistic shock scenarios is then the only viable method to determine the detailed physical origin of the emitted spectra. In this Letter we present synthetic radiation spectra representing the early stage development of the filamentation (streaming) instability of an initially unmagnetized plasma, which is relevant for both collisionless shock formation and reconnection dynamics in relativistic astrophysical outflows, as well as for laboratory astrophysics experiments. Results were obtained using a highly efficient "in situ" diagnostics method, based on detailed particle-in-cell modeling of collisionless plasmas. The synthetic spectra obtained here are compared with those predicted by a semi-analytical model for jitter radiation from the filamentation instability, the latter including self-consistent generated field topologies and particle distributions obtained from the simulations reported upon here. Spectra exhibit dependence on the presence - or absence - of an inert plasma constituent, when comparing baryonic plasmas (i.e. containing protons) with pair plasmas. The results also illustrate that considerable care should be taken when using lower-dimensional models to obtain information about the astrophysical phenomena generating observed spectra.

preprint2010arXiv

The Star Formation Rate of Supersonic MHD Turbulence

This work presents a new physical model of the star formation rate (SFR), verified with an unprecedented set of large numerical simulations of driven, supersonic, self-gravitating, magneto-hydrodynamic (MHD) turbulence, where collapsing cores are captured with accreting sink particles. The model depends on the relative importance of gravitational, turbulent, magnetic, and thermal energies, expressed through the virial parameter, alpha_vir, the rms sonic Mach number, M_S,0, and the ratio of mean gas pressure to mean magnetic pressure, beta_0. The SFR is predicted to decrease with increasing alpha_vir (stronger turbulence relative to gravity), to increase with increasing M_S,0 (for constant values of alpha_vir), and to depend weakly on beta_0 for values typical of star forming regions (M_S,0 ~ 4-20 and beta_0 ~ 1-20). In the unrealistic limit of beta_0 -> infinity, that is in the complete absence of a magnetic field, the SFR increases approximately by a factor of three, which shows the importance of magnetic fields in the star formation process, even when they are relatively weak (super-Alfvenic turbulence). In this non-magnetized limit, our definition of the critical density for star formation has the same dependence on alpha_vir, and almost the same dependence on M_S,0, as in the model of Krumholz and McKee, although our physical derivation does not rely on the concepts of local turbulent pressure and sonic scale. However, our model predicts a different dependence of the SFR on alpha_vir and M_S,0 than the model of Krumholz and McKee. The star-formation simulations used to test the model result in an approximately constant SFR, after an initial transient phase. Both the value of the SFR and its dependence on the virial parameter found in the simulations are shown to agree very well with the theoretical predictions.

preprint2007arXiv

Two regimes of Turbulent Fragmentation and the stellar IMF from Primordial to Present Day Star Formation

The Padoan and Nordlund model of the stellar initial mass function (IMF) is derived from low order statistics of supersonic turbulence, neglecting gravity (e.g. gravitational fragmentation, accretion and merging). In this work the predictions of that model are tested using the largest numerical experiments of supersonic hydrodynamic (HD) and magneto-hydrodynamic (MHD) turbulence to date (~1000^3 computational zones) and three different codes (Enzo, Zeus and the Stagger Code). The model predicts a power law distribution for large masses, related to the turbulence energy power spectrum slope, and the shock jump conditions. This power law mass distribution is confirmed by the numerical experiments. The model also predicts a sharp difference between the HD and MHD regimes, which is recovered in the experiments as well, implying that the magnetic field, even below energy equipartition on the large scale, is a crucial component of the process of turbulent fragmentation. These results suggest that the stellar IMF of primordial stars may differ from that in later epochs of star formation, due to differences in both gas temperature and magnetic field strength. In particular, we find that the IMF of primordial stars born in turbulent clouds may be narrowly peaked around a mass of order 10 solar masses, as long as the column density of such clouds is not much in excess of 10^22 cm^-2.

preprint2001arXiv

Solar Convection and Magneto-Convection Simulations

Magneto-convection simulations with two scenarios have been performed: In one, horizontal magnetic field is advected into the computational domain by fluid entering at the bottom. In the other, an initially uniform vertical magnetic field is imposed on a snapshot of non-magnetic convection and allowed to evolve. In both cases, the field is swept into the intergranular lanes and the boundaries of the underlying mesogranules. The largest field concentrations at the surface reach pressure balance with the surrounding gas. They suppress both horizontal and vertical flows, which reduces the heat transport. They cool, become evacuated and their optical depth unity surface is depressed by several hundred kilometers. Micropores form, typically where a small granule disappears and surrounding flux tubes squeeze into its previous location.

preprint1999arXiv

Super-Sonic Turbulence in the Perseus Molecular Cloud

We compare the statistical properties of J=1-0 13CO spectra observed in the Perseus Molecular Cloud with synthetic J=1-0 13CO spectra, computed solving the non-LTE radiative transfer problem for a model cloud obtained as solutions of the three dimensional magneto-hydrodynamic (MHD) equations. The model cloud is a randomly forced super-Alfvenic and highly super-sonic turbulent isothermal flow. The purpose of the present work is to test if idealized turbulent flows, without self-gravity, stellar radiation, stellar outflows, or any other effect of star formation, are inconsistent or not with statistical properties of star forming molecular clouds. We present several statistical results that demonstrate remarkable similarity between real data and the synthetic cloud. Statistical properties of molecular clouds like Perseus are appropriately described by random super-sonic and super-Alfvenic MHD flows. Although the description of gravity and stellar radiation are essential to understand the formation of single protostars and the effects of star formation in the cloud dynamics, the overall description of the cloud and of the initial conditions for star formation can apparently be provided on intermediate scales without accounting for gravity, stellar radiation, and a detailed modeling of stellar outflows. We also show that the relation between equivalent line width and integrated antenna temperature indicates the presence of a relatively strong magnetic field in the core B1, in agreement with Zeeman splitting measurements.