Search

We examine variations of the stellar initial mass function (IMF) in extreme environments within the formalism derived by Hennebelle \& Chabrier. We focus on conditions encountered in progenitors of massive early type galaxies and starburst regions. We show that, when applying the concept of turbulent Jeans mass as the characteristic mass for fragmentation in a turbulent medium, instead of the standard thermal Jeans mass for purely gravitational fragmentation, the peak of the IMF in such environments is shifted towards smaller masses, leading to a bottom-heavy IMF, as suggested by various observations. In very dense and turbulent environments, we predict that the high-mass tail of the IMF can become even steeper than the standard Salpeter IMF, with a limit for the power law exponent $\alpha\simeq -2.7$, in agreement with recent observational determinations. This steepening is a direct consequence of the high densities and Mach values in such regions but also of the time dependence of the fragmentation process, as incorporated in the Hennebelle-Chabrier theory. We provide analytical parametrizations of these IMFs in such environments, to be used in galaxy evolution calculations. We also calculate the star formation rates and the mass-to-light ratios expected under such extreme conditions and show that they agree well with the values inferred in starburst environments and massive high-redshift galaxies. This reinforces the paradigm of star formation as being a universal process, i.e. the direct outcome of gravitationally unstable fluctuations in a density field initially generated by large scale shock-dominated turbulence. This globally enables us to infer the variations of the stellar IMF and related properties for atypical galactic conditions., Comment: To appear in ApJ

There is observational evidence that central compact objects (CCOs) in supernova remnants have moderately strong magnetic fields $B\sim10^{11}$ G. Meanwhile, available models of partially ionized hydrogen atmospheres of neutron stars with strong magnetic fields are restricted to $B\gtrsim10^{12}$ G. We extend the equation of state and radiative opacities, presented in previous papers for $10^{12}\mbox{ G}\lesssim B \lesssim 10^{15}$ G, to weaker fields. An equation of state and radiative opacities for a partially ionized hydrogen plasma are obtained at magnetic fields $B$, temperatures $T$, and densities $\rho$ typical for atmospheres of CCOs and other isolated neutron stars with moderately strong magnetic fields. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for $3\times10^{10}\mbox{ G}\lesssim B\lesssim 10^{12}$ G, $10^5$ K $\lesssim T\lesssim 10^7$ K, and a wide range of densities. Atmosphere models and spectra are calculated to verify the applicability of the results and to determine the range of magnetic fields and effective temperatures where the incomplete ionization of the hydrogen plasma is important., Comment: 11 pages, 7 figures, accepted for publication in A&A

Understanding the dominant brown dwarf and giant planet formation processes, and finding out whether these processes rely on completely different mechanisms or share common channels represents one of the major challenges of astronomy and remains the subject of heated debates. It is the aim of this review to summarize the latest developments in this field and to address the issue of origin by confronting different brown dwarf and giant planet formation scenarios to presently available observational constraints. As examined in the review, if objects are classified as "Brown Dwarfs" or "Giant Planets" on the basis of their formation mechanism, it has now become clear that their mass domains overlap and that there is no mass limit between these two distinct populations. Furthermore, while there is increasing observational evidence for the existence of non-deuterium burning brown dwarfs, some giant planets, characterized by a significantly metal enriched composition, might be massive enough to ignite deuterium burning in their core. Deuterium burning (or lack of) thus plays no role in either brown dwarf or giant planet formation. Consequently, we argue that the IAU definition to distinguish these two populations has no physical justification and brings scientific confusion. In contrast, brown dwarfs and giant planets might bear some imprints of their formation mechanism, notably in their mean density and in the physical properties of their atmosphere. Future direct imaging surveys will undoubtedly provide crucial information and perhaps provide some clear observational diagnostics to unambiguously distinguish these different astrophysical objects., Comment: 24 pages, 2 figures. Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. S. Klessen, C. P. Dullemond, Th. Henning

This chapter reviews the most recent advancements on the topic of terrestrial and giant planet interiors, including Solar System and extrasolar objects. Starting from an observed mass-radius diagram for known planets in the Universe, we will discuss the various types of planets appearing in this diagram and describe internal structures for each type. The review will summarize the status of theoretical and experimental works performed in the field of equation of states (EOS) for materials relevant to planetary interiors and will address the main theoretical and experimental uncertainties and challenges. It will discuss the impact of new EOS on interior structures and bulk composition determination. We will discuss important dynamical processes which strongly impact the interior and evolutionary properties of planets (e.g plate tectonics, semiconvection) and describe non standard models recently suggested for our giant planets. We will address the case of short-period, strongly irradiated exoplanets and critically analyse some of the physical mechanisms which have been suggested to explain their anomalously large radius., Comment: 24 pages, 8 figures, Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning."

Hydrogen molecules (H2) come in two forms in the interstellar medium, ortho- and para-hydrogen, corresponding to the two different spin configurations of the two hydrogen atoms. The relative abundances of the two flavours in the interstellar medium are still very uncertain, and this abundance ratio has a significant impact on the thermal properties of the gas. In the context of star formation, theoretical studies have recently adopted two different strategies when considering the ortho:para ratio (OPR) of H2 molecules; the first considers the OPR to be frozen at 3:1 while the second assumes that the species are in thermal equilibrium. As the OPR potentially affects the protostellar cores which form as a result of the gravitational collapse of a dense molecular cloud, the aim of this paper is to quantify precisely what role the choice of OPR plays in the properties and evolution of the cores. We used two different ideal gas equations of state for a hydrogen and helium mix in a radiation hydrodynamics code to simulate the collapse of a dense cloud and the formation of the first and second Larson cores; the first equation of state uses a fixed OPR of 3:1 while the second assumes thermal equilibrium. Simulations using an equilibrium ratio collapse faster at early times and show noticeable oscillations around hydrostatic equilibrium, to the point where the core expands for a short time right after its formation before resuming its contraction. In the case of a fixed 3:1 OPR, the core's evolution is a lot smoother. The OPR was however found to have little impact on the size, mass and radius of the two Larson cores. We conclude that if one is solely interested in the final properties of the cores when they are formed, it does not matter which OPR is used. On the other hand, if one's focus lies primarily in the evolution of the first core, the choice of OPR becomes important., Comment: 9 pages, 5 figures. Accepted for publication in Astronomy & Astrophysics

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them., Comment: 24 pages, 5 figures, Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. S. Klessen, C. P. Dullemond, Th. Henning

We present a new model of warm dense matter that represents an intermediate approach between the relative simplicity of ''one-ion'' average atom models and the more realistic but computationally expensive ab initio simulation methods. Physical realism is achieved primarily by including the correlations in the plasma that surrounds a central ion. The plasma is described with the Ornstein-Zernike integral equations theory of fluids, which is coupled to an average atom model for the central ion. In this contribution we emphasize the key elements and approximations and how they relate to and expand upon a typical average atom model. Besides being relatively inexpensive computationally, this approach offers several advantages over ab initio simulations but also has a number of limitations. The model is validated by comparisons with numerical solutions for the pair distribution function of the ions from ab initio simulations for several elements and a wide range of plasma conditions. Simulations results are reproduced remarkably well and simpler limiting theories are recovered as well. This model has many potential applications to calculation of the properties of warm dense matter such as the equation of state and conductivities for a wide range of temperatures and densities., Comment: 26 pages, 10 figures, to appear in the proceedings of the Computational Challenges in Warm Dense Matter Workshop, (Springer, in press)

We study the impact of plasma correlation effects on nonresonant thermonuclear reactions for various stellar objects, namely in the liquid envelopes of neutron stars, and the interiors of white dwarfs, low-mass stars, and substellar objects. We examine in particular the effect of electron screening on the enhancement of thermonuclear reactions in dense plasmas within and beyond the linear mixing rule approximation as well as the corrections due to quantum effects at high density. In addition, we examine some recent unconventional (Yukawa-potential and "quantum-tail") theoretical results on stellar thermonuclear fusions and show that these scenarios do not apply to stellar conditions., Comment: 9 pages, 4 figures; Proceedings of the 14th International Conference on Physics of Nonideal Plasmas "PNP-14" (9-14 September 2012, Rostock, Germany)

We study the evaporation and condensation of CO and CO_2 during the embedded stages of low-mass star formation by using numerical simulations. We focus on the effect of luminosity bursts, similar in magnitude to FUors and EXors, on the gas-phase abundance of CO and CO_2 in the protostellar disk and infalling envelope. The evolution of a young protostar and its environment is followed based on hydrodynamical models using the thin-disk approximation, coupled with a stellar evolution code and phase transformations of CO and CO_2. The accretion and associated luminosity bursts in our model are caused by disk gravitational fragmentation followed by quick migration of the fragments onto the forming protostar. We found that bursts with luminosity on the order of 100-200 L_sun can evaporate CO ices in part of the envelope. The typical freeze-out time of the gas-phase CO onto dust grains in the envelope (a few kyr) is much longer than the burst duration (100-200 yr). This results in an increased abundance of the gas-phase CO in the envelope long after the system has returned into a quiescent stage. In contrast, luminosity bursts can evaporate CO_2 ices only in the disk, where the freeze-out time of the gas-phase CO_2 is comparable to the burst duration. We thus confirm that luminosity bursts can leave long-lasting traces in the abundance of gas-phase CO in the infalling envelope, enabling the detection of recent bursts as suggested by previous semi-analytical studies., Comment: 12 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses., Comment: 15 pages, 11 figures, accepted for publication in A&A

The present paper extends our previous theory of the stellar initial mass function (IMF) by including the time-dependence, and by including the impact of magnetic field. The predicted mass spectra are similar to the time independent ones with slightly shallower slopes at large masses and peak locations shifted toward smaller masses by a factor of a few. Assuming that star-forming clumps follow Larson type relations, we obtain core mass functions in good agreement with the observationally derived IMF, in particular when taking into account the thermodynamics of the gas. The time-dependent theory directly yields an analytical expression for the star formation rate (SFR) at cloud scales. The SFR values agree well with the observational determinations of various Galactic molecular clouds. Furthermore, we show that the SFR does not simply depend linearly on density, as sometimes claimed in the literature, but depends also strongly on the clump mass/size, which yields the observed scatter. We stress, however, that {\it any} SFR theory depends, explicitly or implicitly, on very uncertain assumptions like clump boundaries or the mass of the most massive stars that can form in a given clump, making the final determinations uncertain by a factor of a few. Finally, we derive a fully time-dependent model for the IMF by considering a clump, or a distribution of clumps accreting at a constant rate and thus whose physical properties evolve with time. In spite of its simplicity, this model reproduces reasonably well various features observed in numerical simulations of converging flows. Based on this general theory, we present a paradigm for star formation and the IMF., Comment: accepted for publication in ApJ

As they keep cooling and contracting, Solar System giant planets radiate more energy than they receive from the Sun. Applying the first and second principles of thermodynamics, one can determine their cooling rate, luminosity, and temperature at a given age. Measurements of Saturn's infrared intrinsic luminosity, however, reveal that this planet is significantly brighter than predicted for its age. This excess luminosity is usually attributed to the immiscibility of helium in the hydrogen-rich envelope, leading to "rains" of helium-rich droplets. Existing evolution calculations, however, suggest that the energy released by this sedimentation process may not be sufficient to resolve the puzzle. Here, we demonstrate using planetary evolution models that the presence of layered convection in Saturn's interior, generated, like in some parts of Earth oceans, by the presence of a compositional gradient, significantly reduces its cooling. It can explain the planet's present luminosity for a wide range of configurations without invoking any additional source of energy. This suggests a revision of the conventional homogeneous adiabatic interior paradigm for giant planets, and questions our ability to assess their heavy element content. This reinforces the possibility for layered convection to help explaining the anomalously large observed radii of extrasolar giant planets., Comment: Published in Nature Geoscience. Online publication date: April 21st, 2013. Accepted version before journal editing and with Supplementary Information

Radiative shocks are shocks in a gas where the radiative energy and flux coming from the very hot post-shock material are non-negligible in the shock's total energy budget, and are often large enough to heat the material ahead of the shock. Many simulations of radiative shocks, both in the contexts of astrophysics and laboratory experiments, use a grey treatment of radiative transfer coupled to the hydrodynamics. However, the opacities of the gas show large variations as a function of frequency and this needs to be taken into account if one wishes to reproduce the relevant physics. We have performed radiation hydrodynamics simulations of radiative shocks in Ar using multigroup (frequency dependent) radiative transfer with the HERACLES code. The opacities were taken from the ODALISC database. We show the influence of the number of frequency groups used on the dynamics and morphologies of subcritical and supercritical radiative shocks in Ar gas, and in particular on the extent of the radiative precursor. We find that simulations with even a low number of groups show significant differences compared to single-group (grey) simulations, and that in order to correctly model such shocks, a minimum number of groups is required. Results appear to eventually converge as the number of groups increases above 50. We were also able to resolve in our simulations of supercritical shocks the adaptation zones which connect the cooling layer to the final post-shock state and the precursor. Inside these adaptation zones, we find that the radiative flux just ahead of the shock in one or several high-opacity groups can heat the gas to a temperature higher than the post-shock temperature. Through the use of Hugoniot curves, we have checked the consistency of our radiation hydrodynamics scheme by showing that conservation of mass, momentum and energy (including radiative flux) holds. ABRIDGED., Comment: 19 pages, 12 figures, accepted for publication in JQSRT

We develop an analytical equation of state (EOS) for magnetized fully ionized plasmas, which cover a wide range of temperatures and densities, from low-density classical plasmas to relativistic, quantum plasma conditions. This EOS directly applies to calculations of structure and evolution of strongly magnetized white dwarfs and neutron stars. We review available analytical and numerical results for thermodynamic functions of the nonmagnetized and magnetized Coulomb gases, liquids, and solids. We propose a new analytical expression for the free energy of solid Coulomb mixtures. Based on recent numerical results, we construct analytical approximations for the thermodynamic functions of harmonic Coulomb crystals in quantizing magnetic fields. The analytical description ensures a consistent evaluation of all astrophysically important thermodynamic functions based on the first, second, and mixed derivatives of the free energy. Our numerical code for calculation of thermodynamic functions based on these approximations is made publicly available. Using this code, we calculate and discuss the effects of electron screening and magnetic quantization on the position of the melting point in a range of densities and magnetic fields relevant to white dwarfs and outer envelopes of neutron stars. We consider also the thermal and mechanical structure of a magnetar envelope and argue that it can have a frozen surface which covers the liquid ocean above the solid crust., Comment: 16 pages, 8 figures. In v.2, a typo in Eq.(53) is fixed. In v.3, numerical factors in two inline formulae are corrected (thanks to Dr. Emil Truhlik)

We have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how we have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics.

The idea that non steady accretion during the embedded phase of protostar evolution can produce the observed luminosity spread in the Herzsprung-Russell diagram (HRD) of young clusters has recently been called into question. Observations of Fu Ori, for instance, suggest an expansion of the star during strong accretion events whereas the luminosity spread implies a contraction of the accreting objects, decreasing their radiating surface. In this paper, we present a global scenario based on calculations coupling episodic accretion histories derived from numerical simulations of collapsing cloud prestellar cores of various masses and subsequent protostar evolution. Our calculations show that, assuming an initial protostar mass $\mi \sim 1\,\mjup$, typical of the second Larson's core, both the luminosity spread in the HRD and the inferred properties of Fu Ori events (mass, radius, accretion rate) can be explained by this scenario, providing two conditions. First, there must be some variation within the fraction of accretion energy absorbed by the protostar during the accretion process. Second the range of this variation should increase with increasing accretion burst intensity, and thus with the initial core mass and final star mass. The numerical hydrodynamics simulations of collapsing cloud prestellar cores indeed show that the intensity of the accretion bursts correlates with the mass and initial angular momentum of the prestellar core. Massive prestellar cores with high initial angular momentum are found to produce intense bursts characteristic of Fu Ori like events. Our results thus suggest a link between the burst intensities and the fraction of accretion energy absorbed by the protostar, with some threshold in the accretion rate, of the order of $10^{-5}\msolyr$, delimitating the transition from "cold" to "hot" accretion. [Abridged], Comment: 23 pages, 5 figures, ApJ accepted

Radiative transfer plays a major role in the process of star formation. Many simulations of gravitational collapse of a cold gas cloud followed by the formation of a protostellar core use a grey treatment of radiative transfer coupled to the hydrodynamics. However, dust opacities which dominate extinction show large variations as a function of frequency. In this paper, we used frequency-dependent radiative transfer to investigate the influence of the opacity variations on the properties of Larson's first core. We used a multigroup M1 moment model in a 1D radiation hydrodynamics code to simulate the spherically symmetric collapse of a 1 solar mass cloud core. Monochromatic dust opacities for five different temperature ranges were used to compute Planck and Rosseland means inside each frequency group. The results are very consistent with previous studies and only small differences were observed between the grey and multigroup simulations. For a same central density, the multigroup simulations tend to produce first cores with a slightly higher radius and central temperature. We also performed simulations of the collapse of a 10 and 0.1 solar mass cloud, which showed the properties of the first core to be independent of the initial cloud mass, with again no major differences between grey and multigroup models. For Larson's first collapse, where temperatures remain below 2000 K, the vast majority of the radiation energy lies in the IR regime and the system is optically thick. In this regime, the grey approximation does a good job reproducing the correct opacities, as long as there are no large opacity variations on scales much smaller than the width of the Planck function. The multigroup method is however expected to yield more important differences in the later stages of the collapse when high energy (UV and X-ray) radiation is present and matter and radiation are strongly decoupled., Comment: 9 pages, 5 figures, accepted for publication in A&A

While conventional interior models for Jupiter and Saturn are based on the simplistic assumption of a solid core surrounded by a homogeneous gaseous envelope, we derive new models with an inhomogeneous distribution of heavy elements, i.e. a gradient of composition, within these planets. Such a compositional stratification hampers large scale convection which turns into double-diffusive convection, yielding an inner thermal profile which departs from the traditionally assumed adiabatic interior, affecting these planet heat content and cooling history. To address this problem, we develop an analytical approach of layered double-diffusive convection and apply this formalism to Solar System gaseous giant planet interiors. These models satisfy all observational constraints and yield a metal enrichment for our gaseous giants up to 30 to 60% larger than previously thought. The models also constrain the size of the convective layers within the planets. As the heavy elements tend to be redistributed within the gaseous envelope, the models predict smaller than usual central cores inside Saturn and Jupiter, with possibly no core for this latter. These models open a new window and raise new challenges on our understanding of the internal structure of giant (solar and extrasolar) planets, in particular on the determination of their heavy material content, a key diagnostic for planet formation theories., Comment: Accepted for publication in Astronomy and Astrophysics

We study the plasma correlation effects on nonresonant thermonuclear reactions of carbon and oxygen in the interiors of white dwarfs and liquid envelopes of neutron stars. We examine the effects of electron screening on thermodynamic enhancement of thermonuclear reactions in dense plasmas beyond the linear mixing rule. Using these improved enhancement factors, we calculate carbon and oxygen ignition curves in white dwarfs and neutron stars. The energy balance and ignition conditions in neutron star envelopes are evaluated, taking their detailed thermal structure into account. The result is compared to the simplified "one-zone model," which is routinely used in the literature. We also consider the effect of strong magnetic fields on the ignition curves in the ocean of magnetars., Comment: 9 pages, 5 figures. In v2 and v3, typos are corrected in Eqs.(22) and (3), respectively

The equation of state (EOS) for partially ionized carbon, oxygen, and carbon-oxygen mixtures at temperatures 3\times10^5 K <~ T <~ 3\times10^6 K is calculated over a wide range of densities, using the method of free energy minimization in the framework of the chemical picture of plasmas. The free energy model is an improved extension of our model previously developed for pure carbon (Phys. Rev. E, 72, 046402; arXiv:physics/0510006). The internal partition functions of bound species are calculated by a self-consistent treatment of each ionization stage in the plasma environment taking into account pressure ionization. The long-range Coulomb interactions between ions and screening of the ions by free electrons are included using our previously published analytical model, recently improved, in particular for the case of mixtures. We also propose a simple but accurate method of calculation of the EOS of partially ionized binary mixtures based on detailed ionization balance calculations for pure substances., Comment: 10 pages, 10 figures, accepted for publication in Phys. Rev. E

We present an analytical determination of the star formation rate (SFR) in molecular clouds, based on a time-dependent extension of our analytical theory of the stellar initial mass function (IMF). The theory yields SFR's in good agreement with observations, suggesting that turbulence {\it is} the dominant, initial process responsible for star formation. In contrast to previous SFR theories, the present one does not invoke an ad-hoc density threshold for star formation; instead, the SFR {\it continuously} increases with gas density, naturally yielding two different characteristic regimes, thus two different slopes in the SFR vs gas density relationship, in agreement with observational determinations. Besides the complete SFR derivation, we also provide a simplified expression, which reproduces reasonably well the complete calculations and can easily be used for quick determinations of SFR's in cloud environments. A key property at the heart of both our complete and simplified theory is that the SFR involves a {\it density-dependent dynamical time}, characteristic of each collapsing (prestellar) overdense region in the cloud, instead of one single mean or critical freefall timescale. Unfortunately, the SFR also depends on some ill determined parameters, such as the core-to-star mass conversion efficiency and the crossing timescale. Although we provide estimates for these parameters, their uncertainty hampers a precise quantitative determination of the SFR, within less than a factor of a few., Comment: accepted for publication in ApJL

We present a simple dimensional argument to illustrate the impact of nonthermal support from turbulent velocity dispersion on the shape of the prestellar core mass function (CMF), precursor of the stellar initial mass function (IMF). The argument demonstrates the need to invoke such support to recover the Salpeter slope in the high-mass part of the CMF/IMF, whereas pure thermal support leads to a much steeper slope. This simple dimensional argument clearly highlights the results obtained in the complete Hennebelle-Chabrier theory of the IMF., Comment: Research Note. To appear in A&A

We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field geometries, amplification rates, and saturation levels, decreasing strongly at the transition from subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity., Comment: 4 pages, 3 figures, published in Physical Review Letters (Phys. Rev. Lett. 107, 114504 (2011)), for more information and movies go to http://www.ita.uni-heidelberg.de/~chfeder/pubs/dynamo_prl/dynamo_prl.shtml?lang=en

Radiative shocks play a dominant role in star formation. The accretion shocks on the first and second Larson's cores involve radiative processes and are thus characteristic of radiative shocks. In this study, we explore the formation of the first Larson's core and characterize the radiative and dynamical properties of the accretion shock, using both analytical and numerical approaches. We develop both numerical RHD calculations and a semi-analytical model that characterize radiative shocks in various physical conditions, for radiating or barotropic fluids. Then, we perform 1D spherical collapse calculations of the first Larson's core, using a grey approximation for the opacity of the material. We consider three different models for radiative transfer, namely: the barotropic approximation, the FLD approximation and the more complete M1 model. We investigate the characteristic properties of the collapse and of the first core formation. Comparison between the numerical results and our semi-analytical model shows that this latter reproduces quite well the core properties obtained with the numerical calculations. The accretion shock on the first Larson core is found to be supercritical, implying that all the accretion shock energy on the core is radiated away. The FLD approximation is found to agree quite well with the results based on the M1 model, and is thus appropriate to study the star formation process. In contrast, the barotropic approximation does not correctly describe the thermal properties of the gas during the collapse. We show that a consistent treatment of radiation and hydrodynamics is mandatory to correctly handle the cooling of the gas during the core formation and thus to obtain the correct mechanical and thermal properties for this latter., Comment: 15 pages, 11 figures, accepted for publication in A&A

Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper., Comment: 16 pages, 9 figures, A&A accepted

We quantify the systematic impact of the non-spherical shape of transiting planets and brown dwarfs, due to tidal forces and rotation, on the observed transit depth. Such a departure from sphericity leads to a bias in the derivation of the transit radius from the light curve and affects the comparison with planet structure and evolution models which assume spherical symmetry. As the tidally deformed planet projects its smallest cross section area during the transit, the measured effective radius is smaller than the one of the unperturbed spherical planet. This effect can be corrected by calculating the theoretical shape of the observed planet. We derive simple analytical expressions for the ellipsoidal shape of a fluid object (star or planet) accounting for both tidal and rotational deformations and calibratre it with fully numerical evolution models in the 0.3Mjup-75Mjup mass range. Our calculations yield a 20% effect on the transit depth, i.e. a 10% decrease of the measured radius, for the extreme case of a 1Mjup planet orbiting a Sun-like star at 0.01AU. For the closest planets detected so far (< 0.05 AU), the effect on the radius is of the order of 1 to 10%, by no means a negligible effect, enhancing the puzzling problem of the anomalously large bloated planets. These corrections must thus be taken into account for a correct determination of the radius from the transit light curve. Our analytical expressions can be easily used to calculate these corrections, due to the non-spherical shape of the planet, on the observed transit depth and thus to derive the planet's real equilibrium radius. They can also be used to model ellipsoidal variations of the stellar flux now detected in the CoRoT and Kepler light curves. We also derive directly usable analytical expressions for the moment of inertia, oblateness and Love number (k_2) of a fluid planet as a function of its mass., Comment: 19 pages, 6 figures, 5 tables. Published in A&A. Correction of minor errors in Appendix B. An electronic version of the grids of planetary models is available at http://perso.ens-lyon.fr/jeremy.leconte/JLSite/JLsite/Exoplanets_Simulations.html