Search

Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs in the shocked regions. Simulations show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields and for particle acceleration. These magnetic fields contribute to the electron’s transverse eflection behind the shock. The “jitter” radiation from deflected electrons in turbulent magnetic fields has properties different from synchrotron radiation calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure of gamma-ray bursts, relativistic jets in general, and supernova remnants. In order to calculate radiation from first principles and go beyond the standard synchrotron model, we have used PIC simulations. We present synthetic spectra to compare with the spectra obtained from Fermi observations., Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe, 600

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., Accepted by ApJ

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.

The morphology and kinematics of molecular clouds (MCs) are best explained as the consequence of super--sonic turbulence. Super--sonic turbulence fragments MCs into dense sheets, filaments and cores and large low density ``voids'', via the action of highly radiative shocks. We refer to this process as "turbulent fragmentation". In this work we derive the mass distribution of gravitationally unstable cores generated by the process of turbulent fragmentation. The mass distribution above one solar mass depends primarily on the power spectrum of the turbulent flow and on the jump conditions for isothermal shocks in a magnetized gas. For a power spectrum index \beta=-1.74, consistent with Larson's velocity dispersion--size relation as well as with new numerical and analytic results on super--sonic turbulence, we obtain a power law mass distribution of dense cores with a slope equal to 3/(4-\beta) = 1.33, consistent with the slope of the stellar IMF. Below one solar mass, the mass distribution flattens and turns around at a fraction of a solar mass, as observed for the stellar IMF in a number of stellar clusters, because only the densest cores are gravitationally unstable. The mass distribution at low masses is determined by the probability distribution of the gas density, which is known to be approximately Log--Normal for an isothermal turbulent gas. The intermittent nature of the turbulent density distribution is thus responsible for the existence of a significant number of small collapsing cores, even of sub--stellar mass. Since turbulent fragmentation is unavoidable in super--sonically turbulent molecular clouds, and given the success of the present model in predicting the observed shape of the stellar IMF, we conclude that turbulent fragmentation is essential to the origin of the stellar IMF.

We present a direct numerical and analytical study of driven supersonic MHD turbulence that is believed to govern the dynamics of star-forming molecular clouds. We describe statistical properties of the turbulence by measuring the velocity difference structure functions up to the fifth order. In particular, the velocity power spectrum in the inertial range is found to be close to E(k) \~ k^{-1.74}, and the velocity difference scales as ~ L^{0.42}. The results agree well with the Kolmogorov--Burgers analytical model suggested for supersonic turbulence in [astro-ph/0108300]. We then generalize the model to more realistic, fractal structure of molecular clouds, and show that depending on the fractal dimension of a given molecular cloud, the theoretical value for the velocity spectrum spans the interval [-1.74 ... -1.89], while the corresponding window for the velocity difference scaling exponent is [0.42 ... 0.78]., Comment: 17 pages, 6 figures included

We report on realistic simulations of solar surface convection that are essentially parameter-free, but include detailed physics in the equation of state and radiative energy exchange. The simulation results are compared quantitatively with observations. Excellent agreement is obtained for the distribution of the emergent continuum intensity, the profiles of weak photospheric lines, the p-mode frequencies, the asymmetrical shape of the mode velocity and intensity spectra, the p-mode excitation rate, and the depth of the convection zone. We describe how solar convection is non-local. It is driven from a thin surface thermal boundary layer where radiative cooling produces low entropy gas which forms the cores of the downdrafts in which most of the buoyancy work occurs. Turbulence and vorticity are mostly confined to the intergranular lanes and underlying downdrafts. Finally, we present some preliminary results on magneto-convection.

Supersonic random motions are observed in dark clouds and are traditionally interpreted as Alfven waves, but the possibility that these motions are super-Alfvenic has not been ruled out. In this work we report the results of numerical experiments in two opposite regimes; M_a ~ 1 and M_a >> 1, where M_a is the initial Alfvenic Mach number --the ratio of the rms velocity to the Alfven speed. Our results show that models with M_a >> 1 are consistent with the observed properties of molecular clouds that we have tested --statistics of extinction measurements, Zeeman splitting measurements of magnetic field strength, line width versus integrated antenna temperature of molecular emission line spectra, statistical B-n relation, and scatter in that relation-- while models with M_a ~ 1 have properties that are in conflict with the observations. We find that both the density and the magnetic field in molecular clouds may be very intermittent. The statistical distributions of magnetic field and gas density are related by a power law, with an index that decreases with time in experiments with decaying turbulence. After about one dynamical time it stabilizes at B ~ n^{0.4}. Magnetically dominated cores form early in the evolution, while later on the intermittency in the density field wins out, and also cores with weak field can be generated, by mass accretion along magnetic field lines., 10 figures, 2 tables included

Lada et al. (1994) have described a method for studying the distribution of dust in dark clouds using infrared imaging surveys. In particular they show that the method provides some information about the structure of the gas (dust) on scales smaller than their resolution. In the present work we clarify the nature of the information provided by their method. We show that: 1) the 3D density field of the gas is well described by a Log-Normal distribution down to very small scales; 2) the power spectrum and the standard deviation of the 3D density field can be constrained; 3) the origin of such a structure of the density field is likely to be the supersonic turbulence in the gas. In fact we find a qualitative and quantitative agreement between the predictions based on recent numerical simulations of supersonic turbulence (Nordlund and Padoan 1996; Padoan, Nordlund and Jones 1996) and the constraints given by the infrared dust extinction measurements.

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., Comment: 15 pages, 13 figures accepted for publication in ApJ

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., Comment: 17 pages, 5 color figures, revised version to appear in ApJ, 735, July 2011

Numerical simulations of the multi-phase interstellar medium have been carried out, using a 3D, nonlinear, magnetohydrodynamic, shearing-box model, with random motions driven by supernova explosions. These calculations incorporate the effects of magnetic fields and rotation in 3D; these play important dynamical roles in the galaxy, but are neglected in many other simulations. The supernovae driving the motions are not arbitrarily imposed, but occur where gas accumulates into cold, dense clouds; their implementation uses a physically motivated model for the evolution of such clouds. The process is self-regulating, and produces mean supernova rates as part of the solution. Simulations with differing mean density show a power law relation between the supernova rate and density, with exponent 1.7; this value is within the range suggested from observations (taking star formation rate as a proxy for supernova rate). The global structure of the supernova driven medium is strongly affected by the presence of magnetic fields; e.g. for one solution the filling factor of hot gas is found to vary from 0.19 (with no field) to 0.12 (with initial mid-plane field B_0=6 muG)., 6 pages, 4 figures, v2: minor revisions, now accepted for publication in: Magnetic fields and star formation: theory versus observations, eds: A. I. Gomez de Castro et al., Kluwer Acad. Publ., Dordrecht

Statistical parameters of the ISM driven by thermal energy injections from supernova explosions have been obtained from 3D, nonlinear, magnetohydrodynamic, shearing-box simulations for spiral arm and interarm regions. The density scale height obtained for the interarm regions is 50% larger than within the spiral arms because of the higher gas temperature. The filling factor of the hot gas is also significantly larger between the arms and depends sensitively on magnetic field strength., 6 pages, 4 figures, submitted to proceedings issue of Astrophysics & Space Science following the JENAM 2002 Workshop: "From observations to self-consistent modeling of the ISM in galaxies"

Hundreds of brown dwarfs (BDs) have been discovered in the last few years in stellar clusters and among field stars. BDs are almost as numerous as hydrogen burning stars and so a theory of star formation should also explain their origin. The ``mystery'' of the origin of BDs is that their mass is two orders of magnitude smaller than the average Jeans' mass in star--forming clouds, and yet they are so common. In this work we investigate the possibility that gravitationally unstable protostellar cores of BD mass are formed directly by the process of turbulent fragmentation. Supersonic turbulence in molecular clouds generates a complex density field with a very large density contrast. As a result, a fraction of BD mass cores formed by the turbulent flow are dense enough to be gravitationally unstable. We find that with density, temperature and rms Mach number typical of cluster--forming regions, turbulent fragmentation can account for the observed BD abundance., 11 pages, 3 figures, ApJ submitted Error in equation 1 has been corrected. Improved figures

Numerical simulations of kinematic dynamo action in steady and 3-d ABC flows are presented with special focus on growth rates and multiple periods of the prescribed velocity field. It is found that the difference in growth rate is due to differences in the recycling of the weakest part of the magnetic field. Differences in the topology in cases with and without stagnation points in the imposed velocity field are also investigated, and it is found that the cigar-like structures that develop in the classical A=B=C dynamos, are replaced by ribbon structures in cases where the flow is without stagnation points., 7 pages, 12 PS figures, submitted to A&A (+ slight changes from previous version)

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.

This work presents a new physical model of the star formation rate (SFR), which is 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, αvir, the rms sonic Mach number, , and the ratio of mean gas pressure to mean magnetic pressure, β0. The SFR is predicted to decrease with increasing αvir (stronger turbulence relative to gravity), to increase with increasing (for constant values of αvir), and to depend weakly on β0 for values typical of star forming regions (-20 and β0 1-20). In the unrealistic limit of β0 → ∞, 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). The star-formation simulations used to test the model result in an approximately constant SFR, after an initial transient phase. The dependence of the SFR on the virial parameter is shown to agree very well with the theoretical predictions.

Recent PIC simulations of relativistic electron-positron (electron-ion) jets injected into a stationary medium show that particle acceleration occurs in the shocked regions. Simulations show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields and for particle acceleration. These magnetic fields contribute to the electron's transverse deflection behind the shock. The “jitter” radiation from deflected electrons in turbulent magnetic fields has different properties from synchrotron radiation calculated in a uniform magnetic field. This jitter radiation may be important for understanding the complex time evolution and/or spectral structure of gamma-ray bursts, relativistic jets in general, and supernova remnants. In order to calculate radiation from first principles and go beyond the standard synchrotron model, we have used PIC simulations. We will present detailed spectra for conditions relevant to various astrophysical sites of collisionless shock formation. In particular we will discuss application to GRBs and SNRs.

We present results from numerical simulations of nonlinear MHD dynamo action produced by three-dimensional flows that become turbulent for high values of the fluid Reynolds number. The magnitude of the forcing function driving the flow is allowed to evolve with time in such way as to maintain an approximately constant velocity amplitude (and average kinetic energy) when the flow becomes hydrodynamically unstable. It is found that the saturation level of the dynamo increases with the fluid Reynolds number (at constant magnetic Prandtl number), and that the average growth rate approaches an asymptotic value for high fluid Reynolds number. The generation and destruction of magnetic field is examined during the laminar and turbulent phase of the flow and it is found that in the neighborhood of strong magnetic "flux cigars" Joule dissipation is balanced by the work done against the Lorentz force, while the steady increase of magnetic energy occurs mainly through work done in the weak part of the magnetic field., Comment: 9 pages, 12 figures, accepted by A&A (frontpage illustration): Corrected typos