We address the turbulent fragmentation scenario for the origin of the stellar initial mass function (IMF), using a large set of numerical simulations of randomly driven supersonic MHD turbulence. The turbulent fragmentation model successfully predicts the main features of the observed stellar IMF assuming an isothermal equation of state without any stellar feedback. As a test of the model, we focus on the case of a magnetized isothermal gas, neglecting stellar feedback, while pursuing a large dynamic range in both space and timescales covering the full spectrum of stellar masses from brown dwarfs to massive stars. Our simulations represent a generic 4 pc region within a typical Galactic molecular cloud, with a mass of 3000 Msun and an rms velocity 10 times the isothermal sound speed and 5 times the average Alfven velocity, in agreement with observations. We achieve a maximum resolution of 50 au and a maximum duration of star formation of 4.0 Myr, forming up to a thousand sink particles whose mass distribution closely matches the observed stellar IMF. A large set of medium-size simulations is used to test the sink particle algorithm, while larger simulations are used to test the numerical convergence of the IMF and the dependence of the IMF turnover on physical parameters predicted by the turbulent fragmentation model. We find a clear trend toward numerical convergence and strong support for the model predictions, including the initial time evolution of the IMF. We conclude that the physics of isothermal MHD turbulence is sufficient to explain the origin of the IMF.

In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, 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). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with larger jet radius due to the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation study these kinetic instabilities lead to new types of instabilities in global jets. In the electron-proton jet simulation a modified recollimation occurs and jet electrons are strongly perturbed. In the electron-positron jet simulation mixed kinetic instabilities occur at early times followed by a turbulence-like structure. Simulations using much larger (and longer) systems are further required in order to thoroughly investigate the evolution of global 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.

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.

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.

Turbulence is ubiquitous in molecular clouds (MCs), but its origin is still unclear because MCs are usually assumed to live longer than the turbulence dissipation time. Interstellar medium (ISM) turbulence is likely driven by SN explosions, but it has never been demonstrated that SN explosions can establish and maintain a turbulent cascade inside MCs consistent with the observations. In this work, we carry out a simulation of SN-driven turbulence in a volume of (250 pc)$^3$, specifically designed to test if SN driving alone can be responsible for the observed turbulence inside MCs. We find that SN driving establishes a velocity scaling consistent with the usual scaling laws of supersonic turbulence, suggesting that previous idealized simulations of MC turbulence, driven with a random, large-scale volume force, were correctly adopted as appropriate models for MC turbulence, despite the artificial driving. We also find that the same scaling laws extend to the interior of MCs, and that the velocity-size relation of the MCs selected from our simulation is consistent with that of MCs from the Outer-Galaxy Survey, the largest MC sample available. The mass-size relation and the mass and size probability distributions also compare successfully with those of the Outer Galaxy Survey. Finally, we show that MC turbulence is super-Alfv\'{e}nic with respect to both the mean and rms magnetic-field strength. We conclude that MC structure and dynamics are the natural result of SN-driven turbulence.

One of the key open questions in the study of relativistic jets is their interaction with the environment. Here, we study the initial evolution of both electron-proton and electron-positron relativistic jets, focusing on their lateral interaction with the ambient plasma. We trace the generation and evolution of the toroidal magnetic field generated by both kinetic Kelvin-Helmholtz (kKH) and Mushroom instabilities (MI). This magnetic field collimates the jet. We show that in electron-proton jet, electrons are perpendicularly accelerated with jet collimation. The magnetic polarity switches from the clockwise to anti-clockwise in the middle of jet, as the instabilities weaken. For the electron-positron jet, we find strong mixture of electron-positron with the ambient plasma, that results in the creation of a bow shock. Merger of magnetic field current filaments generate density bumps which initiate a forward shock. The strong mixing between jet and ambient particles prevents full development of the jet on the studied scale. Our results therefore provide a direct evidence for both jet collimation and particle acceleration in the created bow shock. Differences in the magnetic field structures generated by electron-proton and electron-positron jets may contribute to observable differences in the polarized properties of emission by electrons.

We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and Lorentz factor. In the velocity shear between jet and ambient plasmas, strong AC ($e^{\pm}$ plasmas) or DC ($e^-$- $p^+$ plasmas) magnetic fields are generated via the kinetic Kelvin-Helmholtz instability (kKHI), and the magnetic field structure also depends on the jet Lorentz factor. We have calculated, self-consistently, the radiation from electrons accelerated in shock generated magnetic fields. The spectra depend on the jet's initial Lorentz factor and temperature via the resulting particle acceleration and magnetic field generation. Our ongoing "Global" jet simulations containing shocks and velocity shears will provide us with the ability to calculate and model the complex time evolution and/or spectral structure observed from gamma-ray bursts, AGN jets, and supernova remnants.

We use the adaptive mesh refinement code RAMSES to model the formation of protoplanetary disks in realistic star formation environments. The resolution scales over up to 29 powers of two ($\sim$ 9 orders of magnitude) covering a range from outer scales of 40 pc to inner scales of 0.015 AU. The accretion rate from a 1.5 solar mass envelope peaks near $10^{-4}$ $\mspy$ about 6 kyr after sink particle formation and then decays approximately exponentially, reaching $10^{-6}$ $\mspy$ in 100 kyr. The models suggest universal scalings of physical properties with radius during the main accretion phase, with kinetic and / or magnetic energy in approximate balance with gravitational energy. Efficient accretion is made possible by the braking action of the magnetic field, which nevertheless allows a near-Keplerian disk to grow to a 100 AU size. The magnetic field strength ranges from more than 10 G at 0.1 AU to less than 1 mG at 100 AU, and drives a time dependent bipolar outflow, with a collimated jet and a broader disk wind.

We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath configuration is slightly different even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field $E_{\rm z}$ and the magnetic field $B_{\rm y}$. After the $B_{\rm y}$ component is excited, an induced electric field $E_{\rm x}$ becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios $m_{\rm i}/m_{\rm e} = 1836$ and $m_{\rm i}/m_{\rm e} = 20$ are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case ($\gamma_{\rm j} = 1.5$) is larger than for a relativistic jet case ($\gamma_{\rm j} = 15$).

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 $\gamma$--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 $\gamma$--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.

We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants.

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.

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\'{e}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, $\epsilon_{\rm ff}$, on the above parameters. Our main results are: i) $\epsilon_{\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 $\epsilon_{\rm ff}$, but only to a point, beyond which $\epsilon_{\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, $\epsilon_{\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 $\epsilon_{\rm ff}$, and depending only on $t_{\rm ff}/t_{\rm dyn}$: $\epsilon_{\rm ff} \approx \epsilon_{\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.

lasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. Our initial results of a jet-ambient interaction with anti-parallel magnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to reconnection and associated particle acceleration. We will investigate the radiation in transient stage as a possible generation mechanism of precursors of prompt emission. In our simulations we calculate the radiation from electrons in the shock region. The detailed properties of this radiation are important for understanding the complex time evolution and spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

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.

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.

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.

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.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system. We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system generates a jitter/synchrotron spectrum.

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.

The explosive BN/KL outflow emerging from OMC1 behind the Orion Nebula may have been powered by the dynamical decay of a non-hierarchical multiple system $\sim$500 years ago that ejected the massive stars I, BN, and source n, with velocities of about 10 to 30 km s$^{-1}$. New proper motion measurements of H$_2$ features show that within the errors of measurement, the outflow originated from the site of stellar ejection. Combined with published data, these measurements indicate an outflow age of $\sim$500 years, similar to the time since stellar ejection. The total kinetic energy of the ejected stars and the outflow is about 2 to $6 \times 10^{47}$ ergs. It is proposed that the gravitational potential energy released by the formation of a short-period binary, most likely source I, resulted in stellar ejection and powered the outflow. A scenario is presented for the formation of a compact, non-hierarchical multiple star system, its decay into an ejected binary and two high-velocity stars, and launch of the outflow. Three mechanisms may have contributed to the explosion in the gas: (i) Unbinding of the circum-cluster envelope following stellar ejection, (ii) disruption of circumstellar disks and high-speed expulsion of the resulting debris during the final stellar encounter, and (iii) the release of stored magnetic energy. Plausible proto-stellar disk end envelope properties can produce the observed outflow mass, velocity, and kinetic energy distributions. The ejected stars may have acquired new disks by fall-back or Bondi-Hoyle accretion with axes roughly orthogonal to their velocities. The expulsion of gas and stars from OMC1 may have been driven by stellar interactions.

S-type stars are late-type giants whose atmosphere is enriched in carbon and s-process elements because of either extrinsic pollution by a binary companion or intrinsic nucleosynthesis and dredge-up on the thermally-pulsing AGB. A large grid of S-star model atmospheres has been computed covering the range 2700 < Teff < 4000 K with 0.5 < C/O < 0.99. ZrO and TiO band strength indices as well as VJHKL photometry are needed to disentangle Teff, C/O and [s/Fe]. A "best-model finding tool" was developed using a set of well-chosen indices and checked against photometry as well as low- and high-resolution spectroscopy. It is found that applying M-star model atmospheres (i.e., with a solar C/O ratio) to S stars can lead to errors on Teff up to 400K. We constrain the parameter space occupied by S stars of the vast sample of Henize stars in terms of Teff, [C/O] and [s/Fe].

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.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. We have also performed simulations with electron-ion jets. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability for electron-positron jets and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks for both cases. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value for pair plasmas. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system with two different cases for Lorentz factors (15 and 100). We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability.

We present a preliminary local thermodynamic equilibrium (LTE) abundance analysis of the template halo red giant HD122563 based on a realistic, three-dimensional (3D), time-dependent, hydrodynamical model atmosphere of the very metal-poor star. We compare the results of the 3D analysis with the abundances derived by means of a standard LTE analysis based on a classical, 1D, hydrostatic model atmosphere of the star. Due to the different upper photospheric temperature stratifications predicted by 1D and 3D models, we find large, negative, 3D-1D LTE abundance differences for low-excitation OH and Fe I lines. We also find trends with lower excitation potential in the derived Fe LTE abundances from Fe I lines, in both the 1D and 3D analyses. Such trends may be attributed to the neglected departures from LTE in the spectral line formation calculations.

Implications of the recently discovered systematic abundance difference between the Sun and two collections of `solar twins' are discussed. The differences can be understood as an imprint on the abundances of the solar convection zone caused by the lock-up of heavy elements in the planets. Such a scenario also leads naturally to possible solutions of two other abundance peculiarities; 1) the discrepancy between photospheric abundances derived from accurate 3-D models of the solar photosphere and the abundance of heavy elements in the solar interior deduced from helioseismology, and 2) the abundance pattern of Jupiter, which can either--with great difficulty--be interpreted as a general and similar overabundance of both common elements such as carbon, nitrogen and sulphur and rare inert gases such as argon, krypton and xenon, or--much more simply--as an under-abundance of hydrogen.

Numerical simulations have by now revealed that the fine scale structure of the penumbra in general and the Evershed effect in particular is due to overturning convection, mainly confined to gaps with strongly reduced magnetic field strength. The Evershed flow is the radial component of the overturning convective flow visible at the surface. It is directed outwards -- away from the umbra -- because of the broken symmetry due to the inclined magnetic field. The dark penumbral filament cores visible at high resolution are caused by the 'cusps' in the magnetic field that form above the gaps. Still remaining to be established are the details of what determines the average luminosity of penumbrae, the widths, lengths, and filling factors of penumbral filaments, and the amplitudes and filling factors of the Evershed flow. These are likely to depend at least partially also on numerical aspects such as limited resolution and model size, but mainly on physical properties that have not yet been adequately determined or calibrated, such as the plasma beta profile inside sunspots at depth and its horizontal profile, the entropy of ascending flows in the penumbra, etc.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique.

Results of realistic simulations of solar surface convection on the scale of supergranules (96 Mm wide by 20 Mm deep) are presented. The simulations cover only 10% of the geometric depth of the solar convection zone, but half its pressure scale heights. They include the hydrogen, first and most of the second helium ionization zones. The horizontal velocity spectrum is a power law and the horizontal size of the dominant convective cells increases with increasing depth. Convection is driven by buoyancy work which is largest close to the surface, but significant over the entire domain. Close to the surface buoyancy driving is balanced by the divergence of the kinetic energy flux, but deeper down it is balanced by dissipation. The damping length of the turbulent kinetic energy is 4 pressure scale heights. The mass mixing length is 1.8 scale heights. Two thirds of the area is upflowing fluid except very close to the surface. The internal (ionization) energy flux is the largest contributor to the convective flux for temperatures less than 40,000 K and the thermal energy flux is the largest contributor at higher temperatures. This data set is useful for validating local helioseismic inversion methods. Sixteen hours of data are available as four hour averages, with two hour cadence, at steinr.msu.edu/~bob/96averages, as idl save files. The variables stored are the density, temperature, sound speed, and three velocity components. In addition, the three velocity components at 200 km above mean continuum optical depth unity are available at 30 sec. cadence.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless, relativistic shocks, particle (electron, positron, and ion) acceleration is due to plasma waves and their associated instabilities (e.g., the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform, small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting "jitter" radiation from deflected electrons has different properties compared to synchrotron radiation, which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general, and supernova remnants.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the presence of relativistic jets, instabilities such as the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability create collisionless shocks, which are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons in small-scale magnetic fields has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation, a case of diffusive synchrotron radiation, may be important to understand the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Modeling plasma physical processes in astrophysical context demands for both detailed kinetics and large scale development of the electromagnetic field densities. We present a new framework for modeling plasma physics of hot tenuous plasmas by a two-split scheme, in which the large scale fields are modeled by means of a particle-in-cell (PIC) code, and in which binary collision processes and single-particle processes are modeled through a Monte-Carlo approach. Our novel simulation tool -- the PhotonPlasma code -- is a unique hybrid model; it combines a highly parallelized (Vlasov) particle-in-cell approach with continuous weighting of particles and a sub-Debye Monte-Carlo binary particle interaction framework. As an illustration of the capabilities we present results from a numerical study of Gamma-Ray Burst - Circumburst Medium interaction and plasma preconditioning via Compton scattering. We argue that important microphysical processes can only viably be investigated by means of hybrid codes such as the PhotonPlasma code. Our first results from 3D simulations with this new simulation tool suggest that magnetic fields and plasma filaments are created in the wakefield of prompt gamma-ray bursts. Furthermore, the photon flux density gradient impacts on particle acceleration in the burst head and wakefield. We discuss some possible implications of the circumburst medium being preconditioned for a trailing afterglow shock front. We also discuss important improvements for future studies of GRB wakefields processes, using the PhotonPlasma code.

We have constructed a grid of about 10,000 spherically symmetric and plane-parallel models with the MARCS program, and make it available for public use. Parameter ranges are: Teff=2500 to 8000 K, log g =log(GM/R2)= -1 to 5 (cgs) with various masses and radii, [Me/H]=-5 to +1, with [Alpha/Fe] = 0.0 and 0.4 and different choices of C and N abundances to also represent stars of types R, S and N, and with microturbulence parameters from 1 to 5 km/s. We also list fluxes in approximately 108,000 wavelength points. Underlying assumptions in addition to 1D stratification include hydrostatic equilibrium, MLT convection and LTE. A number of general properties of the models are discussed, in relation to the effects of changing blanketing and sphericity. Models are compared with other available grids and excellent agreement is found with plane-parallel models of Castelli and Kurucz within the overlapping parameter range. Although there are departures from the spherically symmetric NextGen models, the agreement with more recent PHOENIX models is gratifying. The models of the grid show regularities, but some interesting departures from general patterns occur for the coolest models due to the molecular opacities. We have tested rules of thumb concerning effects of blanketing and sphericity and found them to often be astonishingly accurate. Some interesting new phenomena have been discovered, such as the intricate coupling between blanketing and sphericity, and the strong effects of carbon enhancement on metal-poor models. We give further details of models and comparisons with observations in subsequent papers.

We discuss a numerical 3D radiation-MHD simulation of penumbral fine structure in a small sunspot. This simulation shows the development of short filamentary structures with horizontal flows, similar to observed Evershed flows, and an inward propagation of these structures at a speed compatible with observations. We conclude that the Evershed flow represents the horizontal flow component of overturning convection in gaps with strongly reduced field strength. The top of the flow is always directed outward--away from the umbra-- because of the broken symmetry due to the inclined magnetic field. Upflows occur in the inner parts of the gaps and most of the gas turns over radially (outwards and sideways), and descends back down again. The ascending, cooling and overturning flow tends to bend magnetic field lines down, forcing a weakening of the field that makes it easier for gas located in an adjacent layer--further in--to initiate a similar sequence of motion, aided by lateral heating, thus causing the inward propagation of the filament.

During the last two decades, the focus of star formation research has shifted from understanding the collapse of a single dense core into a star to studying the formation hundreds to thousands of stars in molecular clouds. In this chapter, we overview recent observational and theoretical progress toward understanding star formation on the scale of molecular clouds and complexes, i.e the macrophysics of star formation. We begin with an overview of recent surveys of young stellar objects (YSOs) in molecular clouds and embedded clusters, and we outline an emerging picture of cluster formation. We then discuss the role of turbulence to both support clouds and create dense, gravitationally unstable structures, with an emphasis on the role of magnetic fields (in the case of distributed stars) and feedback (in the case of clusters) to slow turbulent decay and mediate the rate and density of star formation. The discussion is followed by an overview of how gravity and turbulence may produce observed scaling laws for the properties of molecular clouds, stars and star clusters, and how the observed, low star formation rate may result from self regulated star formation. We end with some concluding remarks, including a number of questions to be addressed by future observations and simulations.

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.

Recent three-dimensional radiative hydrodynamics simulations of protoplanetary disks report disparate disk behaviors, and these differences involve the importance of convection to disk cooling, the dependence of disk cooling on metallicity, and the stability of disks against fragmentation and clump formation. To guarantee trustworthy results, a radiative physics algorithm must demonstrate the capability to handle both the high and low optical depth regimes. We develop a test suite that can be used to demonstrate an algorithm's ability to relax to known analytic flux and temperature distributions, to follow a contracting slab, and to inhibit or permit convection appropriately. We then show that the radiative algorithm employed by Meji\'a (2004) and Boley et al. (2006) and the algorithm employed by Cai et al. (2006) and Cai et al. (2007, in prep.) pass these tests with reasonable accuracy. In addition, we discuss a new algorithm that couples flux-limited diffusion with vertical rays, we apply the test suite, and we discuss the results of evolving the Boley et al. (2006) disk with this new routine. Although the outcome is significantly different in detail with the new algorithm, we obtain the same qualitative answers. Our disk does not cool fast due to convection, and it is stable to fragmentation. We find an effective $\alpha\approx 10^{-2}$. In addition, transport is dominated by low-order modes.

We present results of numerical 3D MHD simulations with radiative energy transfer of fine structure in a small sunspot of about 4 Mm width. The simulations show the development of filamentary structures and flow patterns that are, except for the lengths of the filaments, very similar to those observed. The filamentary structures consist of gaps with reduced field strength relative to their surroundings. Calculated synthetic images show dark cores like those seen in the observations; the dark cores are the result of a locally elevated $\tau=1$ surface. The magnetic field in these cores is weaker and more horizontal than for adjacent brighter structures, and the core support a systematic outflow. Movies show migration of the dark-cored structures towards the umbra, and fragments of magnetic flux that are carried away from the spot by a large scale `moat flow'. We conclude that the simulations are in qualitative agreement with observed penumbra filamentary structures, Evershed flows and moving magnetic features.

We extend semi-analytical computations of excitation rates for solar oscillation modes to those of other solar-like oscillating stars to compare them with recent observations. Numerical 3D simulations of surface convective zones of several solar-type oscillating stars are used to characterize the turbulent spectra as well as to constrain the convective velocities and turbulent entropy fluctuations in the uppermost part of the convective zone of such stars. These constraints, coupled with a theoretical model for stochastic excitation, provide the rate 'P' at which energy is injected into the p-modes by turbulent convection. These energy rates are compared with those derived directly from the 3D simulations. The excitation rates obtained from the 3D simulations are systematically lower than those computed from the semi-analytical excitation model. We find that Pmax, the excitation rate maximum, scales as (L/M)^s where s is the slope of the power law and L and M are the mass and luminosity of the 1D stellar model built consistently with the associated 3D simulation. The slope is found to depend significantly on the adopted form of the eddy time-correlation ; using a Lorentzian form results in s=2.6, whereas a Gaussian one gives s=3.1. Finally, values of Vmax, the maximum in the mode velocity, are estimated from the computed power laws for Pmax and we find that Vmax increases as (L/M)^sv. Comparisons with the currently available ground-based observations show that the computations assuming a Lorentzian eddy time-correlation yield a slope, sv, closer to the observed one than the slope obtained when assuming a Gaussian. We show that the spatial resolution of the 3D simulations must be high enough to obtain accurate computed energy rates.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.

We review five mechanisms for forming brown dwarfs: (i) turbulent fragmentation of molecular clouds, producing very low-mass prestellar cores by shock compression; (ii) collapse and fragmentation of more massive prestellar cores; (iii) disc fragmentation; (iv) premature ejection of protostellar embryos from their natal cores; and (v) photo-erosion of pre-existing cores overrun by HII regions. These mechanisms are not mutually exclusive. Their relative importance probably depends on environment, and should be judged by their ability to reproduce the brown-dwarf IMF, the distribution and kinematics of newly formed brown dwarfs, the binary statistics of brown dwarfs, the ability of brown dwarfs to retain discs, and hence their ability to sustain accretion and outflows. This will require more sophisticated numerical modelling than is presently possible, in particular more realistic initial conditions and more realistic treatments of radiation transport, angular momentum transport and magnetic fields. We discuss the minimum mass for brown dwarfs, and how brown dwarfs should be distinguished from planets.

An efficient algorithm for calculating radiative transfer on massively parallel computers using domain decomposition is presented. The integral formulation of the transfer equation is used to divide the problem into a local but compute-intensive part for calculating the intensity and optical depth integrals, and a nonlocal part for communicating the intensity between adjacent processors. The waiting time of idle processors during the nonlocal communication part does not have a severe impact on the scaling. The wall clock time thus scales nearly linearly with the inverse number of processors.

We investigate the 50-year old hypothesis that the magnetic fields of the Ap stars are stable equilibria that have survived in these stars since their formation. With numerical simulations we find that stable magnetic field configurations indeed appear to exist under the conditions in the radiative interior of a star. Confirming a hypothesis by Prendergast (1956), the configurations have roughly equal poloidal and toroidal field strengths. We find that tori of such twisted fields can form as remnants of the decay of an unstable random initial field. In agreement with observations, the appearance at the surface is an approximate dipole with smaller contributions from higher multipoles, and the surface field strength can increase with the age of the star. The results of this paper were summarised by Braithwaite & Spruit (2004).

Numerical simulations of turbulent stratified convection are used to study models with approximately the same convective flux, but different radiative fluxes. As the radiative flux is decreased, for constant convective flux: the entropy jump at the top of the convection zone becomes steeper, the temperature fluctuations increase and the velocity fluctuations decrease in magnitude, and the distance that low entropy fluid from the surface can penetrate increases. Velocity and temperature fluctuations follow mixing length scaling laws.

The excitation of internal gravity waves by penetrative convective plumes is investigated using 2-D direct simulations of compressible convection. The wave generation is quantitatively studied from the linear response of the radiative zone to the plumes penetration, using projections onto the g-modes solutions of the associated linear eigenvalue problem for the perturbations. This allows an accurate determination of both the spectrum and amplitudes of the stochastically excited modes. Using time-frequency diagrams of the mode amplitudes, we then show that the lifetime of a mode is around twice its period and that during times of significant excitation up to 40% of the total kinetic energy may be contained into g-modes.

A forward model is described in which we synthesize spectra from an ab-initio 3D MHD simulation of an outer stellar atmosphere, where the coronal heating is based on braiding of magnetic flux due to photospheric footpoint motions. We discuss the validity of assumptions such as ionization equilibrium and investigate the applicability of diagnostics like the differential emission measure inversion. We find that the general appearance of the synthesized corona is similar to the solar corona and that, on a statistical basis, integral quantities such as average Doppler shifts or differential emission measures are reproduced remarkably well. The persistent redshifts in the transition region, which have puzzled theorists since their discovery, are explained by this model as caused by the flows induced by the heating through braiding of magnetic flux. While the model corona is only slowly evolving in intensity, as is observed, the amount of structure and variability in Doppler shift is very large. This emphasizes the need for fast coronal spectroscopy, as the dynamical response of the corona to the heating process manifests itself in a comparably slow evolving coronal intensity but rapid changes in Doppler shift.

Gamma ray bursts are among the most energetic events in the known universe. A highly relativistic fireball is ejected. In most cases the burst itself is followed by an afterglow, emitted under deceleration as the fireball plunges through the circum-stellar media. To interpret the observations of the afterglow emission, two physical aspects need to be understood: 1) The origin and nature of the magnetic field in the fireball and 2) the particle velocity distribution function behind the shock. Both are necessary in existing afterglow models to account for what is believed to be synchrotron radiation. To answer these questions, we need to understand the microphysics at play in collisionless shocks. Using 3D particle-in-cell simulations we can gain insight in the microphysical processes that take place in such shocks. We discuss the results of such computer experiments. It is shown how a Weibel-like two-stream plasma instability is able to create a strong transverse intermittent magnetic field and points to a connected mechanism for in situ particle acceleration in the shock region.

Accretion rates of order 10^-8 M_\odot/yr are observed in young protostars of approximately a solar mass with evidence of circumstellar disks. The accretion rate is significantly lower for protostars of smaller mass, approximately proportional to the second power of the stellar mass, \dot{M}_accr\propto M^2. The traditional view is that the observed accretion is the consequence of the angular momentum transport in isolated protostellar disks, controlled by disk turbulence or self--gravity. However, these processes are not well understood and the observed protostellar accretion, a fundamental aspect of star formation, remains an unsolved problem. In this letter we propose the protostellar accretion rate is controlled by accretion from the large scale gas distribution in the parent cloud, not by the isolated disk evolution. Describing this process as Bondi--Hoyle accretion, we obtain accretion rates comparable to the observed ones. We also reproduce the observed dependence of the accretion rate on the protostellar mass. These results are based on realistic values of the ambient gas density and velocity, as inferred from numerical simulations of star formation in self--gravitating turbulent clouds.

We propose to interpret the stellar IMF as a property of the turbulence in the star--forming gas. Gravitationally unstable density enhancements in the turbulent flow collapse and form stars. Their mass distribution can be derived analytically from the power spectrum of the turbulent flow and the isothermal shock jump conditions in the 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 supersonic 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 Salpeter's 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 Log--Normal distribution of the gas density. The intermittent nature of this distribution is responsible for the generation of a significant number of collapsing cores of brown dwarf mass.

The origin of brown dwarfs (BDs) is an important component of the theory of star formation, because BDs are approximately as numerous as solar mass stars. It has been suggested that BDs originate from the gravitational fragmentation of protostellar disks, a very different mechanism from the formation of hydrogen burning stars. We propose that BDs are instead formed by the process of turbulent fragmentation, like more massive stars. In numerical simulations of turbulence and star formation we find that gravitationally unstable density peaks of BD mass are commonly formed by the turbulent flow. These density peaks collapse into BD mass objects with circumstellar disks, like more massive protostars. We rely on numerical experiments with very large resolution, achieved with adaptive mesh refinement (AMR). The turbulence simulation presented here is the first AMR turbulence experiment ever attempted and achieves an effective resolution of 1024^3 computational zones. The star formation simulation achieves an effective resolution of (10^6)^3 computational zones, from a cloud size of 5pc to protostellar disks resolved down to 1AU.

We compute the rates P at which energy is injected into the p modes for a set of 3D simulations of outer layers of stars. We found that Pmax - the maximum in P - scales as (L/M)^s where s is the slope of the power law, L and M are the luminosity and the mass of the 1D stellar models associated with the simulations. The slope is found to depend significantly on the adopted representation for the turbulent eddy-time correlation function, chi_k. According to the expected performances of COROT, it will likely be possible to measure Pmax as a function of L/M and to constrain the properties of stellar turbulence as the turbulent eddy time-correlation.

Cool stars like our Sun are surrounded by a million degree hot outer atmosphere, the corona. Since more than 60 years the physical nature of the processes heating the corona to temperatures well in excess of those on the stellar surface remain puzzling. Recent progress in observational techniques and numerical modeling now opens a new window to approach this problem. We present the first coronal emission line spectra synthesized from three-dimensional numerical models describing the evolution of the dynamics and energetics as well as of the magnetic field in the corona. In these models the corona is heated through motions on the stellar surface that lead to a braiding of magnetic field lines inducing currents which are finally dissipated. These forward models enable us to synthesize observed properties like (average) emission line Doppler shifts or emission measures in the outer atmosphere, which until now have not been understood theoretically, even though many suggestions have been made in the past. As our model passes these observational tests, we conclude that the flux braiding mechanism is a prime candidate for being the dominant heating process of the magnetically closed corona of the Sun and solar-like stars.

Hydrodynamical, 3D simulations of the outer layers of the Sun and Alpha Cen A are used to obtain constraints on the properties of turbulent convection in such stars. These constraints enable us to compute - on the base of a theoretical model of stochastic excitation - the rate P at which p modes are excited by turbulent convection in those two stars. Results are then compared with solar seismic observations and recent observations of Alpha Cen A. For the Sun, a good agreement between observations and computed P is obtained. For Alpha Cen A a large discrepancy is obtained which origin cannot be yet identified: it can either be caused by the present data quality which is not sufficient for our purpose or by the way the intrinsic amplitudes and the life-times of the modes are determined or finally attributed to our present modelling. Nevertheless, data with higher quality or/and more adapted data reductions will likely provide constraints on the p-mode excitation mechanism in Alpha Cen A.

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.

We present the first center-to-limb G-band images synthesized from high resolution simulations of solar magneto-convection. Towards the limb the simulations show "hilly" granulation with dark bands on the far side, bright granulation walls and striated faculae, similar to observations. At disk center G-band bright points are flanked by dark lanes. The increased brightness in magnetic elements is due to their lower density compared with the surrounding intergranular medium. One thus sees deeper layers where the temperature is higher. At a given geometric height, the magnetic elements are cooler than the surrounding medium. In the G-band, the contrast is further increased by the destruction of CH in the low density magnetic elements. The optical depth unity surface is very corrugated. Bright granules have their continuum optical depth unity 80 km above the mean surface, the magnetic elements 200-300 km below. The horizontal temperature gradient is especially large next to flux concentrations. When viewed at an angle, the deep magnetic elements optical surface is hidden by the granules and the bright points are no longer visible, except where the "magnetic valleys" are aligned with the line of sight. Towards the limb, the low density in the strong magnetic elements causes unit line-of-sight optical depth to occur deeper in the granule walls behind than for rays not going through magnetic elements and variations in the field strength produce a striated appearance in the bright granule walls.

Supersonic turbulent flows of magnetized gas are believed to play an important role in the dynamics of star-forming clouds in galaxies. Understanding statistical properties of such flows is crucial for developing a theory of star formation. In this letter we propose a unified approach for obtaining the velocity scaling in compressible and super--Alfv\'{e}nic turbulence, valid for arbitrary sonic Mach number, \ms. We demonstrate with numerical simulations that the scaling can be described with the She--L\'{e}v\^{e}que formalism, where only one parameter, interpreted as the Hausdorff dimension of the most intense dissipative structures, needs to be varied as a function of \ms. Our results thus provide a method for obtaining the velocity scaling in interstellar clouds once their Mach numbers have been inferred from observations.

We present results from three-dimensional particle simulations of collisionless shocks with relativistic counter-streaming ion-electron plasmas. Particles are followed over many skin depths downstream of the shock. Open boundaries allow the experiments to be continued for several particle crossing times. The experiments confirm the generation of strong magnetic and electric fields by a Weibel-like kinetic streaming instability, and demonstrate that the electromagnetic fields propagate far downstream of the shock. The magnetic fields are predominantly transversal, and are associated with merging ion current channels. The total magnetic energy grows as the ion channels merge, and as the magnetic field patterns propagate down stream. The electron populations are quickly thermalized, while the ion populations retain distinct bulk speeds in shielded ion channels and thermalize much more slowly. These results may help explain the origin of the magnetic fields responsible for afterglow synchrotron/jitter radiation from Gamma-Ray Bursts.

The excitation of gravity waves by penetrative convective plumes is investigated using 2D direct simulations of compressible convection. The oscillation field is measured by a new technique based on the projection of our simulation data onto the theoretical g-modes solutions of the associated linear eigenvalue problem. This allows us to determine both the excited modes and their corresponding amplitudes accurately.

The magnetic field strength in molecular clouds is a fundamental quantity for theories of star formation. It is estimated by Zeeman splitting measurements in a few dense molecular cores, but its volume--averaged value within large molecular clouds (over several parsecs) is still uncertain. In this work we provide a new method to constrain the average magnetic field strength in molecular clouds. We compare the power spectrum of gas density of molecular clouds with that of two $350^3$ numerical simulations of supersonic MHD turbulence. The numerical simulation with approximate equipartition of kinetic and magnetic energies (model A) yields the column density power spectrum $P(k)\propto k^{-2.25\pm 0.01}$, the super--Alfv\'{e}nic simulation (model B) $P(k)\propto k^{-2.71\pm 0.01}$. The column density power spectrum of the Perseus, Taurus and Rosetta molecular cloud complexes is found to be well approximated by a power law, $P_{\rm o}(k)\propto k^{-a}$, with $a=2.74\pm 0.07$, $2.74\pm 0.08$ and $2.76\pm 0.08$ respectively. We conclude that the observations are consistent with the presence of super--Alfv\'{e}nic turbulence in molecular clouds (model B) while model A is inconsistent (more than 99% confidence) with the observations.

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).

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.

A 3D simulation of the upper part of the solar convective zone is used to obtain information on the frequency component, chi_k, of the correlation product of the turbulent velocity field. This component plays an important role in the stochastic excitation of acoustic oscillations. A time analysis of the solar simulation shows that a gaussian function does not correctly reproduce the nu-dependency of chi_k inferred from the 3D simuation in the frequency range where the acoustic energy injected into the solar p modes is important (nu ~ 2 - 4 mHz). The nu-dependency of chi_k is fitted with different analytical functions which can then conveniently be used to compute the acoustic energy supply rate P injected into the solar radial oscillations. With constraints from a 3D simulation, adjustment of free parameters to solar data is no longer necessary and is not performed here. The result is compared with solar seismic data. Computed values of P obtained with the analytical function which fits best chi_k are found ~ 2.7 times larger than those obtained with the gaussian model and reproduce better the solar seismic observations. This non-gaussian description also leads to a Reynolds stress contribution of the same order as the one arising from the advection of the turbulent fluctuations of entropy by the turbulent motions. Some discrepancy between observed and computed P values still exist at high frequency and possible causes for this discrepancy are discussed.

We present numerical results from plasma particle simulations of collisionless shocks and ultra-relativistic counter-streaming plasmas. We demonstrate how the field-particle interactions lead to particle acceleration behind the shock-front. Further, we demonstrate how ultra relativistic counter-streaming plasmas create large scale patchy magnetic field structures and that these field structures propagate down-stream of the shock front. These results may help explain the origin of the magnetic fields and accelerated electrons responsible for afterglow synchrotron radiation from gamma ray bursts.

Analyses of a 3D simulation of the upper layers of a solar convective envelope provide constraints on the physical quantities which enter the theoretical formulation of a stochastic excitation model of solar p modes, for instance the convective velocities and the turbulent kinetic energy spectrum. These constraints are then used to compute the acoustic excitation rate for solar p modes, P. The resulting values are found ~5 times larger than the values resulting from a computation in which convective velocities and entropy fluctuations are obtained with a 1D solar envelope model built with the time-dependent, nonlocal Gough (1977) extension of the mixing length formulation for convection (GMLT). This difference is mainly due to the assumed mean anisotropy properties of the velocity field in the excitation region. The 3D simulation suggests much larger horizontal velocities compared to vertical ones than in the 1D GMLT solar model. The values of P obtained with the 3D simulation constraints however are still too small compared with the values inferred from solar observations. Improvements in the description of the turbulent kinetic energy spectrum and its depth dependence yield further increased theoretical values of P which bring them closer to the observations. It is also found that the source of excitation arising from the advection of the turbulent fluctuations of entropy by the turbulent movements contributes ~ 65-75 % to the excitation and therefore remains dominant over the Reynolds stress contribution. The derived theoretical values of P obtained with the 3D simulation constraints remain smaller by a factor ~3 compared with the solar observations. This shows that the stochastic excitation model still needs to be improved.

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.

We present here numerical simulations of surface solar convection which cover a box of 30$\times30\times$3.2 Mm$^3$ with a resolution of 315$\times315\times$82, which is used to investigate the dynamics of scales larger than granulation. No structure resembling supergranulation is present; possibly higher Reynolds numbers (i.e. higher numerical resolution), or magnetic fields, or greater depth are necessary. The results also show interesting aspects of granular dynamics which are briefly presented, like extensive p-mode ridges in the k-$\omega$ diagram and a ringlike distribution of horizontal vorticity around granules. At large scales, the horizontal velocity is much larger than the vertical velocity and the vertical motion is dominated by p-mode oscillations.

The excitation rate P of solar p-modes is computed with a model of stochatic excitation which involves constraints on the averaged properties of the solar turbulence. These constraints are obtained from a 3D simulation. Resulting values for P are found ~9 times larger than when the calculation assumes properties of turbulent convection which are derived from an 1D solar model based on Gough (1977)'s formulation of the mixing-length theory (GMLT). This difference is mainly due to the assumed values for the mean anisotropy of the velocity field in each case.Calculations based on 3D constraints bring the P maximum closer to the observational one.We also compute P for several models of intermediate mass stars (1Mo < M < 2Mo). Differences in the values of P_max between models computed with the classical mixing-length theory and GMLT models are found large enough for main sequence stars to suggest that measurements of P in this mass range will be able to discriminate between different models of turbulent convection.

Our recent study, based on a 3D numerical simulation, suggested that the stochastic excitation of solar p-modes has a non gaussian nature and a non gaussian model was proposed which results in a good agreement between the computed rate $P$ at which p-modes are excited and the solar seismic observational constraints. In the present work, we study some consequences of this non-gaussian model on the rate P at which p-modes are excited in intermediate mass stars (1 Mo < M < 2 Mo). The non gaussian model changes substantially the spectrum of P at high frequency compared with the commonly assumed gaussian model. The largest effects on P are found for stars in the mass range ~ 1.5-1.6 Mo. They are found large enough that observatins in solar-like oscillating stars observed with space based seismology missions (e.g. COROT) will be able to confirm or reject the gaussian nature of stochastic excitation.

We simulate convection near the solar surface, where the continuum optical depth is of order unity. Hence, to determine the radiative heating and cooling in the energy conservation equation, we must solve the radiative transfer equation (instead of using the diffusion or optically thin cooling approximations). A method efficient enough to calculate the radiation for thousands of time steps is needed. We assume LTE and a non-gray opacity grouped into 4 bins according to strength. We perform a formal solution of the Feautrier equation along a vertical and four straight, slanted, rays (at four azimuthal angles which are rotated 15 deg. every time step). We present details of our method. We also give some results: comparing simulated and observed line profiles for the Sun, showing the importance of 3D transfer for the structure of the mean atmosphere and the eigenfrequencies of p-modes, illustrating Stokes profiles for micropores, and analyzing the effect of radiation on p-mode asymmetries.

Magneto-convection simulations on meso-granule and granule scales near the solar surface are used to study small scale dynamo activity, the emergence and disappearance of magnetic flux tubes, and the formation and evolution of micropores. From weak seed fields, convective motions produce highly intermittent magnetic fields in the intergranular lanes which collect over the boundaries of the underlying meso-granular scale cells. Instances of both emerging magnetic flux loops and magnetic flux disappearing from the surface occur in the simulations. We show an example of a flux tube collapsing to kG field strength and discuss how the nature of flux disappearance can be investigated. Observed stokes profiles of small magnetic structures are severely distorted by telescope diffraction and seeing. Because of the strong stratification, there is little recycling of plasma and field in the surface layers. Recycling instead occurs by exchange with the deep layers of the convection zone. Plasma and field from the surface descend through the convection zone and rise again toward the surface. Because only a tiny fraction of plasma rising up from deep in the convection zone reaches the surface due to mass conservation, little of the magnetic energy resides in the near surface layers. Thus the dynamo acting on weak incoherent fields is global, rather than a local surface dynamo.

Supersonic turbulence fragments the interstellar medium into dense sheets, filaments, cores and large low density voids. The turbulence is driven on large scales, probably predominantly by supernovae. The scaling properties of supersonic turbulence are well described by a new analytical theory. The distribution of proto-stellar core masses depends primarily on the power spectrum of the turbulent flow, and on the jump conditions for near isothermal shocks in a magnetized gas. For the predicted velocity power spectrum index beta=1.74, consistent with results of numerical experiments of supersonic turbulence as well as with Larson's velocity-size relation, one obtains by scaling arguments a power law mass distribution of dense cores with a slope equal to 3/(4-beta) = 1.33, consistent with the slope of the Salpeter stellar initial mass function (IMF). Results from numerical simulations confirm this scaling. Both the analytical model for the stellar IMF and its numerical estimate show that turbulent fragmentation can also explain the origin of brown dwarfs.

A 3D simulation of the upper part of the solar convective zone is used to derive constraints about the averaged and dynamic properties of solar turbulent convection. Theses constraints are then used to compute the acoustic energy supply rate P(nu) injected into the solar radial oscillations according to the theoretical expression in Samadi & Goupil (2001). The result is compared with solar seismic data. Assuming, as it is usually done, a gaussian model for the frequency (nu) component chi_k(nu) of the model of turbulence, it is found that the computed P(nu) is underestimated compared with the solar seismic data by a factor ~ 2.5. A frequency analysis of the solar simulation shows that the gaussian model indeed does not correctly model chi_k(nu) in the frequency range where the acoustic energy injected into the solar p-modes is important (nu ~ 2 - 4 mHz). One must consider an additional non-gaussian component for chi_k(nu) to reproduce its behavior. Computed values of P obtained with this non-gaussian component reproduce better the solar seismic observations. This non-gaussian component leads to a Reynolds stress contribution of the same order than the one arising from the advection of the turbulent fluctuations of entropy by the turbulent motions.

We compute the structure function scaling of the integrated intensity images of two J=1-0 13CO maps of Taurus and Perseus. The scaling exponents of the structure functions follow the velocity scaling of supersonic turbulence, suggesting that turbulence plays an important role in the fragmentation of cold interstellar clouds. The data also allows to verify the validity of the two basic assumptions of the hierarchical symmetry model, originally proposed for the derivation of the velocity structure function scaling. This shows that the same hierarchical symmetry holds for the projected density field of cold interstellar clouds.

A magnetically dominated plasma driven by motions on boundaries at which magnetic field lines are anchored is forced to dissipate the work being done upon it, no matter how small the electrical resistivity. Numerical experiments have clarified that dissipation is achieved through the formation of a hierarchy of electrical current sheets. The probability distribution function of the local winding of magnetic field lines is nearly Gaussian, with a width of the order unity. The dissipation is highly irregular in space and time, but the average level of dissipation is well described by a scaling law that is independent of the electrical resistivity. If the boundary driving is suspended for a period of time the magnetic dissipation rapidly drops to insignificant levels, leaving the magnetic field in a nearly force-free state. Renewed boundary driving leads to a quick return to dissipation levels compatible with the rate of boundary work, with dissipation starting much more rapidly than when starting from idealized initial conditions with a uniform magnetic field. Application of these concepts to the solar corona lends credibility to realistic, three-dimensional numerical models that predict emission measures, coronal structures, and heating rates compatible with observations.

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 have performed 3-D numerical magnetohydrodynamic (MHD) jet experiments to study the instabilities associated with strongly toroidal magnetic fields and determine if such magnetic configurations in jets are as unstable as similar situations are found to be in the laboratory and in analytical estimates. A perfectly conducting Keplerian disc with fixed density, rotational velocity and pressure is used as a lower boundary for the jet. The energy equation is solved, with the inclusion of self-consistently computed heating by viscous and magnetic dissipation. The resulting jets evolve into time-dependent, non-axisymmetric configurations, but we find only minor disruption of the jets by for example the kink instability. We find that magnetic dissipation may have profound effects on the jet flow as: 1) it turns on in highly wound up magnetic field regions and helps to prevent critical kink situations; 2) it influences jet dynamics by re-organizing the magnetic field structure and increasing thermal pressure in the jet; and 3) it influences mass loading by increasing temperature and pressure at the base of the jet.

The solar acoustic p-mode line profiles are asymmetric. Velocity spectra have more power on the low-frequency sides, whereas intensity profiles show the opposite sense of asymmetry. Numerical simulations of the upper convection zone have resonant p-modes with the same asymmetries and asymmetry reversal as the observed modes. The temperature and velocity power spectra at optical depth $\tau_{\rm cont} = 1$ have the opposite asymmetry as is observed for the intensity and velocity spectra. At a fixed geometrical depth, corresponding to $<\tau_{\rm cont}>=1$, however, the temperature and velocity spectra have the same asymmetry. This indicates that the asymmetry reversal is produced by radiative transfer effects and not by correlated noise.

The heating of the solar corona and the puzzle of the slender high reaching magnetic loops seen in observations from the Transition Region And Coronal Explorer(TRACE) has been investigated through 3D numerical simulations, and found to be caused by the well observed plasma flows in the photosphere displacing the footpoints of magnetic loops in a nearly potential configuration. It is found that even the small convective displacements cause magnetic dissipation sufficient to heat the corona to temperatures of the order of a million Kelvin. The heating is intermittent in both space and time - at any one height and time it spans several orders of magnitude, and localized heating causes transonic flows along field lines, which explains the observed non-hydrostatic equilibrium of loops that are bright in emission measure.

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.

The interstellar medium (ISM) provides a unique laboratory for highly supersonic, driven hydrodynamics turbulence. We present a theory of such turbulence, confirm it by numerical simulations, and use the results to explain observational properties of interstellar molecular clouds, the regions where stars are born.

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.

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 <|\Delta u|>~ 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].

We present 3-d numerical magneto-hydrodynamic simulations of a buoyant, twisted magnetic flux rope embedded in a stratified, solar-like model convection zone. The flux rope is given an initial twist such that it neither kinks nor fragments during its ascent. Moreover, its magnetic energy content with respect to convection is chosen so that the flux rope retains its basic geometry while being deflected from a purely vertical ascent by convective flows. The simulations show that magnetic flux is advected away from the core of the flux rope as it interacts with the convection. The results thus support the idea that the amount of toroidal flux stored at or near the bottom of the solar convection zone may currently be underestimated.

We have computed line emission cooling rates for the main cooling species in models of interstellar molecular clouds. The models are based on numerical simulations of super-sonic magneto-hydrodynamic (MHD) turbulence. Non-LTE radiative transfer calculations have been performed to properly account for the complex density and velocity structures in the MHD simulations. Three models are used. Two of the models are based on MHD simulations with different magnetic field strength and the third includes the computation of self-gravity (in the super-Alfvenic regime of turbulence). The density and velocity fields in the simulations are determined self-consistently by the dynamics of super-sonic turbulence. The models are intended to represent molecular clouds with linear size L~6 pc and mean density <n>~300 cm^-3, with the density exceeding 10^4 cm^-3 in the densest cores. We present 12CO, 13CO, C18O, O2, OI, CI and H2O cooling rates in isothermal clouds with kinetic temperatures 10-80K. Analytical approximations are derived for the cooling rates. The inhomogeneity of the models reduces photon trapping and enhances the cooling in the densest parts of the clouds. Compared with earlier models the cooling rates are less affected by optical depth effects and are therefore higher. The main effects comes, however, from the density variation since cooling efficiency increases with density. This is very important for the cooling of the clouds as a whole since most cooling is provided by gas with density above the average.

Using a numerical simulation of compressible convection with radiative transfer mimicking the solar photosphere, we compare the velocity field derived from granule motions to the actual velocity field of the plasma. We thus test the idea that granules may be used to trace large-scale velocity fields at the sun's surface. Our results show that this is indeed the case provided the scale separation is sufficient. We thus estimate that neither velocity fields at scales less than 2500 km nor time evolution at scales shorter than 0.5 hr can be faithfully described by granules. At larger scales the granular motions correlate linearly with the underlying fluid motions with a slope of ~< 2 reaching correlation coefficients up to ~0.9.

We model the polarized thermal dust emission from protostellar cores that are assembled by super-sonic turbulent flows in molecular clouds. Self-gravitating cores are selected from a three dimensional simulation of super-sonic and super-Alfvenic magneto-hydrodynamic (MHD) turbulence. The polarization is computed in two ways. In model A it is assumed that dust properties and grain alignment efficiency are uniform; in model B it is assumed that grains are not aligned at visual extinction larger than 3 mag. The main results of this work are: i) Values of the degree of polarization P between 1 and 10% are typical, despite the super-Alfvenic nature of the turbulence; ii) A steep decrease of P with increasing values of the sub-mm dust continuum intensity I is always found in self--gravitating cores selected from the MHD simulations, if grains are not aligned above a certain value of visual extinction (model B); iii) The same behavior is hard to reproduce if grains are aligned independently of visual extinction (model A); iv) The Chandrasekhar-Fermi formula, corrected by a factor f=0.4, provides an approximate estimate of the average magnetic field strength in the cores. Sub-mm dust continuum polarization maps of quiescent protostellar cores and Bok globules always show a decrease in P with increasing value of I consistent with the predictions of our model B. We therefore conclude that sub-mm polarization maps of quiescent cores do not map the magnetic field inside the cores at visual extinction larger than approximately 3 mag. There is no inconsistency between the results from optical and near-IR polarized absorption of background stars, and the observed polarization of sub-mm dust continuum from quiescent cores. In both cases, grains at large visual extinction appear to be virtually unaligned.

Warped, precessing accretion discs appear in a range of astrophysical systems, for instance the X-ray binary Her X-1 and in the active nucleus of NGC4258. In a warped accretion disc there are horizontal pressure gradients that drive an epicyclic motion. We have studied the interaction of this epicyclic motion with the magnetohydrodynamic turbulence in numerical simulations. We find that the turbulent stress acting on the epicyclic motion is comparable in size to the stress that drives the accretion, however an important ingredient in the damping of the epicyclic motion is its parametric decay into inertial waves.

Super-sonic turbulence fragments molecular clouds (MC) into a very complex density field with density contrasts of several orders of magnitude. A fraction of the gas is locked into dense and gravitationally bound cores, which collapse as proto-stars. This process can be studied with numerical simulations of super-sonic self-gravitating turbulence. In this work, we use numerical simulations of magneto-hydrodynamic (MHD), super-sonic, super-Alfv\'{e}nic and self-gravitating turbulence to compute the mass distribution of collapsing proto-stellar cores, which are selected as local density maxima. We find that the mass distribution of collapsing cores is consistent with the stellar initial mass function (IMF), suggesting that super-sonic turbulence may be responsible for the generation of the IMF. To support this conclusion we also show that the physical properties of the numerically selected cores are in agreement with the properties of observed NH$_3$ cores and that their magnetic field strength is consistent with Zeeman splitting measurements. In turbulent MCs, star formation occurs via the gravitational collapse of super-critical cores, formed by the turbulent flow, sub-critical cores being irrelevant for the process of star formation.

The fragmentation of molecular clouds (MC) into protostellar cores is a central aspect of the process of star formation. Because of the turbulent nature of super-sonic motions in MCs, it has been suggested that dense structures such as filaments and clumps are formed by shocks in a turbulent flow. In this work we present strong evidence in favor of the turbulent origin of the fragmentation of MCs. The most generic result of turbulent fragmentation is that dense post shock gas traces a gas component with a smaller velocity dispersion than lower density gas, since shocks correspond to regions of converging flows, where the kinetic energy of the turbulent motion is dissipated. Using synthetic maps of spectra of molecular transitions, computed from the results of numerical simulations of super--sonic turbulence, we show that the dependence of velocity dispersion on gas density generates an observable relation between the rms velocity centroid and the integrated intensity (column density), Sigma(V_0)-I, which is indeed found in the observational data. The comparison between the theoretical model (maps of synthetic 13CO spectra), with 13CO maps from the Perseus, Rosette and Taurus MC complexes, shows excellent agreement in the Sigma(V_0)-I relation. The Sigma(V_0)-I relation of different observational maps with the same total rms velocity are remarkably similar, which is a strong indication of their origin from a very general property of the fluid equations, such as the turbulent fragmentation process.

We have studied the fate of initial magnetic fields in the hot halo gas out of which the visible parts of galaxies form, using three-dimensional numerical MHD-experiments. The halo gas undergoes compression by several orders of magnitude in the subsonic cooling flow that forms the cold disk. The magnetic field is carried along and is amplified considerably in the process, reaching muG levels for reasonable values of the initial ratio of magnetic to thermal energy density.

The interaction of magnetic fields and stratified convection was studied in the context of the solar and late type stellar dynamos by using numerical 3D MHD simulations. The topology of stratified asymmetric and over-turning convection enables a pumping mechanism that may render the magnetic flux storage problem obsolete. The inclusion of open boundary conditions leads to a considerable flux loss unless the open boundary is placed close to the physical boundary. Simulations including solar-like latitudinal shear indicates that a toroidal field of several tens of kilo-Gauss may be held down by the pumping mechanism.

The computationally demanding nature of radiative-hydrodynamical simulations of stellar surface convection warrants an investigation of the sensitivity of the convective structure and spectral synthesis to the numerical resolution and dimension of the simulations, which is presented here. With too coarse a resolution the predicted spectral lines tend to be too narrow, reflecting insufficient Doppler broadening from the convective motions, while at the currently highest affordable resolution the line shapes have converged essentially perfectly to the observed profiles. Similar conclusions are drawn from the line asymmetries and shifts. In terms of abundances, weak FeI and FeII lines show a very small dependence (~0.02 dex) while for intermediate strong lines with significant non-thermal broadening the sensitivity increases (~0.10 dex). Problems arise when using 2D convection simulations to describe an inherent 3D phenomenon, which translates to inaccurate atmospheric velocity fields and temperature and pressure structures. In 2D the theoretical line profiles tend to be too shallow and broad compared with the 3D calculations and observations, in particular for intermediate strong lines. In terms of abundances, the 2D results are systematically about 0.1 dex lower than for the 3D case for FeI lines. Furthermore, the predicted line asymmetries and shifts are much inferior in 2D. Given these shortcomings and computing time considerations it is better to use 3D simulations of even modest resolution than high-resolution 2D simulations.

Realistic ab-initio 3D, radiative-hydrodynamical convection simulations of the solar granulation have been applied to FeI and FeII line formation. In contrast to classical analyses based on hydrostatic 1D model atmospheres the procedure contains no adjustable free parameters but the treatment of the numerical viscosity in the construction of the 3D, time-dependent, inhomogeneous model atmosphere and the elemental abundance in the 3D spectral synthesis. However, the numerical viscosity is introduced purely for numerical stability purposes and is determined from standard hydrodynamical test cases with no adjustments allowed to improve the agreement with the observational constraints from the solar granulation. The non-thermal line broadening is mainly provided by the Doppler shifts arising from the convective flows in the solar photosphere and the solar oscillations. The almost perfect agreement between the predicted temporally and spatially averaged line profiles for weak Fe lines with the observed profiles and the absence of trends in derived abundances with line strengths, seem to imply that the micro- and macroturbulence concepts are obsolete in these 3D analyses. Furthermore, the theoretical line asymmetries and shifts show a very satisfactory agreement with observations with an accuracy of typically 50-100 m/s on an absolute velocity scale. The remaining minor discrepancies point to how the convection simulations can be refined further.

The solar photospheric Fe abundance has been determined using realistic ab initio 3D, time-dependent, hydrodynamical model atmospheres. The study is based on the excellent agreement between the predicted and observed line profiles directly rather than equivalent width, since the intrinsic Doppler broadening from the convective motions and oscillations provide the necessary non-thermal broadening. Thus, three of the four hotly debated parameters (equivalent widths, microturbulence and damping enhancement factors) in the center of the recent solar Fe abundance dispute regarding FeI lines no longer enter the analysis, leaving the transition probabilities as the main uncertainty. Both FeI (using the samples of lines of both the Oxford and Kiel studies) and FeII lines have been investigated, which give consistent results: log FeI = 7.44 +- 0.05 and log FeII = 7.45 +- 0.10. Also the wings of strong FeI lines return consistent abundances, log FeII = 7.42 +- 0.03, but due to the uncertainties inherent in analyses of strong lines we give this determination lower weight than the results from weak and intermediate strong lines. In view of the recent slight downward revision of the meteoritic Fe abundance log Fe = 7.46 +- 0.01, the agreement between the meteoritic and photospheric values is very good, thus appearingly settling the debate over the photospheric Fe abundance from FeI lines.

We excite an epicyclic motion, whose amplitude depends on the vertical position, $z$, in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays we can obtain information about the interaction between the warp and the disc turbulence. A high amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, $\alpha_{\rm v}$, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, $\alpha_{\rm h}$, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, $\alpha_{\rm v}$ is approximately equal to $\alpha_{\rm h}$ and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness ($\sim 0.01$) of $\alpha_{\rm v}$ and $\alpha_{\rm h}$ we conclude that for $\beta = p_{\rm gas}/p_{\rm mag} \sim 10$ the timescale for diffusion or damping of a warp is much shorter than the usual viscous timescale. Finally, we review the astrophysical implications.

We use the three-dimensional hydrodynamic code of Stein and Nordlund to realistically simulate the upper layers of the solar convection zone in order to study physical characteristics of solar oscillations. Our first result is that the properties of oscillation modes in the simulation closely match the observed properties. Recent observations from SOHO/MDI and GONG have confirmed the asymmetry of solar oscillation line profiles, initially discovered by Duvall et al. In this paper we compare the line profiles in the power spectra of the Doppler velocity and continuum intensity oscillations from the SOHO/MDI observations with the simulation. We also compare the phase differences between the velocity and intensity data. We have found that the simulated line profiles are asymmetric and have the same asymmetry reversal between velocity and intensity as observed. The phase difference between the velocity and intensity signals is negative at low frequencies and jumps in the vicinity of modes as is also observed. Thus, our numerical model reproduces the basic observed properties of solar oscillations, and allows us to study the physical properties which are not observed.

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.

Realistic 3-dimensional (3D), radiative hydrodynamical surface convection simulations of the metal-poor halo stars HD 140283 and HD 84937 have been performed. Due to the dominance of adiabatic cooling over radiative heating very low atmospheric temperatures are encountered. The lack of spectral lines in these metal-poor stars thus causes much steeper temperature gradients than in classical 1D hydrostatic model atmospheres where the temperature of the optically thin layers is determined by radiative equilibrium. The modified atmospheric structures cause changes in the emergent stellar spectra. In particular, the primordial Li abundances may have been overestimated by 0.2-0.35 dex with 1D model atmospheres. However, we caution that our result assumes local thermodynamic equilibrium (LTE), while the steep temperature gradients may be prone to e.g. over-ionization.

Differences between observed and theoretical eigenfrequencies of the Sun have characteristics which identify them as arising predominantly from properties of the oscillations in the vicinity of the solar surface: in the super-adiabatic, convective boundary layer and above. These frequency differences may therefore provide useful information about the structure of these regions, precisely where the theory of solar structure is most uncertain. In the present work we use numerical simulations of the outer part of the Sun to quantify the influence of turbulent convection on solar oscillation frequencies. Separating the influence into effects on the mean model and effects on the physics of the modes, we find that the main model effects are due to the turbulent pressure that provides additional support against gravity, and thermal differences between average 3-D models and 1-D models. Surfaces of constant pressure in the visible photosphere are elevated by about 150 km, relative to a standard envelope model. As a result, the turning points of high-frequency modes are raised, while those of the low-frequency modes remain essentially unaffected. The corresponding gradual lowering of the mode frequencies accounts for most of the frequency difference between observations and standard solar models. Additional effects are expected to come primarily from changes in the physics of the modes, in particular from the modulation of the turbulent pressure by the oscillations.

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.

The question of the shape of the density PDF for supersonic turbulence is addressed, using both analytical and numerical methods. For isothermal supersonic turbulence, the PDF is Log-Normal, with a width that scales approximately linearly with the Mach number. For a polytropic equation of state, with an effective gamma smaller than one, the PDF becomes skewed and reminiscent of (but not equal to) a power law on the high density side.

The interaction of magnetic fields and convection is investigated in the context of the coronal heating problem. We study the motions of photospheric magnetic elements using filtergrams obtained at the Swedish Vacuum Solar Telescope at La Palma. We use potential-field modeling to extrapolate the magnetic and velocity fields to larger height. We find that the velocity in the chromosphere can be locally enhanced at the separatrix surfaces between neighboring flux tubes. The predicted velocities are several km/s, significantly larger than those of the photospheric flux tubes, which may have important implications for coronal heating. s

The dynamics of molecular clouds is characterized by supersonic random motions in the presence of a magnetic field. We study this situation using numerical solutions of the three-dimensional compressible magneto-hydrodynamic (MHD) equations in a regime of highly supersonic random motions. The non-LTE radiative transfer calculations are performed through the complex density and velocity fields obtained as solutions of the MHD equations, and more than 5x10^5 synthetic molecular spectra are obtained. We use a numerical flow without gravity or external forcing. The flow is super-Alfvenic and corresponds to model A of Padoan and Nordlund (1997). Synthetic data consist of sets of 90x90 synthetic spectra with 60 velocity channels, in five molecular transitions: J=1-0 and J=2-1 for 12CO and 13CO, and J=1-0 for CS. Though we do not consider the effects of stellar radiation, gravity, or mechanical energy input from discrete sources, our models do contain the basic physics of magneto-fluid dynamics and non-LTE radiation transfer and are therefore more realistic than previous calculations. As a result, these synthetic maps and spectra bear a remarkable resemblance to the corresponding observations of real clouds.

We compare the statistical properties of J=1-0 13CO spectra observed in the Perseus Molecular Cloud with synthetic J=1-0 13CO spectra from a 5 pc model cloud. The synthetic spectra are computed solving the non-LTE radiative transfer problem for a model cloud obtained as solutions of the 3-D magneto-hydrodynamic (MHD) equations in both the highly supersonic and super-Alfvenic regimes of random flows. We present several statistical results that demonstrate remarkable similarity between real data and the synthetic cloud. The three-dimensional structure and dynamics of molecular clouds like Perseus are appropriately described by random supersonic 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 described on intermediate scales without accounting for gravity and stellar radiation.

In this work we use models of molecular clouds (MC), and non-LTE radiative transfer calculations, to obtain a theoretical calibration of the relation between LTE 13CO column density and true column density in MCs. The cloud models consist of 3 dimensional grids of density and velocity fields obtained as solutions of the compressible magneto-hydrodynamic equations in a 128x128x128 periodic grid in both the supersonic and super-Alfvenic regimes. Due to the random nature of the velocity field and the presence of shocks, the density spans a continuous range of values covering over 5-6 orders of magnitude (from ~0.1 to ~10^5 cm^-3). As a result, the LTE column density can be calibrated over 3 orders of magnitude. We find that LTE column density of molecular clouds typically underestimates the mean 13CO true column density by a factor ranging from 1.3 to 7. These results imply that the standard LTE methods for the derivation of column densities from CO data systematically underestimate the true values independent of other major sources of uncertainty such as the relative abundance of CO.

Supersonic random motions are observed in dark clouds and are traditionally interpreted as Alfv\'{e}n 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: beta ~1 and beta << 1, where beta is the ratio of gas pressure and magnetic pressure: beta=P_g/P_m. Our results, combined with observational tests, show that the model with beta ~ 1 is consistent with the observed properties of molecular clouds, while the model with beta << 1 has several properties that are in conflict with the observations. We also 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. Magnetically dominated cores form early in the evolution, while later on the intermittency in the density field wins out.

We propose that the stellar initial mass function (IMF) is universal in the sense that its functional form arises as a consequence of the statistics of random supersonic flows. A model is developed for the origin of the stellar IMF, that contains a dependence on the average physical parameters (temperature, density, velocity dispersion) of the large scale site of star formation. The model is based on recent numerical experiments of highly supersonic random flows that have a strong observational counterpart. It is shown that a Miller-Scalo like IMF is naturally produced by the model for the typical physical conditions in molecular clouds. A more ``massive'' IMF in star bursts is also predicted.

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.

Following the report of solar-like oscillations in the G0 V star eta Boo (Kjeldsen et al. 1995, AJ 109, 1313), a first attempt to model the observed frequencies was made by Christensen-Dalsgaard et al. (1995, ApJ Letters, in press). This attempt succeeded in reproducing the observed frequency separations, although there remained a difference of about 10 microHz between observed and computed frequencies. In those models, the near-surface region of the star was treated rather crudely. Here we consider more sophisticated models that include non-local mixing-length theory, turbulent pressure and nonadiabatic oscillations.