SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Although massive stars play a critical role in the production of turbulent energy in the ISM, in the formation and destruction of molecular clouds, and ultimately in the dynamical and chemodynamical evolution of galaxies, our understanding of the sequence of events which leads to their formation is still rather limited.
Cosmochemical analyses of meteorites imply that maximum temperatures in the inner solar nebula were on the order of 1300 K, yet standard viscous accretion disk models predict much lower midplane temperatures (approx. 300 K at 2 AU to 3 AU) in a minimum mass nebula. A second-order accurate radiative hydrodynamics code has been used to construct models of the late-phase solar nebula appropriate for low-mass star formation (M is approximately 10(exp -6) to 10(exp -5) solar-M yr(exp -1). For a minimum mass (0.02 solar-M) nebula and a solar-mass protostar, the new models show that compressional heating due to mass accretion onto the nebula and subsequent vertical contraction of the nebula are sufficient to lead to midplane temperatures T(sub m) greater than 1400 K at 1 AU and T(sub m) greater than 1000 K at 2.5 AU.
We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕ ) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to $M^{2/3}_{p}$. The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10 -5 M ⊕ yr -1 . In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙ , the solids’ accretion falls below 10 -6 M ⊕ yr -1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕ , predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10 -7 to ≈10 -6 L ⊙ .
Hydrodynamical calculations of the evolution of a collapsing, rotating axisymmetric 10 solar masses molecular clump, including the effects of radiative acceleration but without magnetic fields, are represented. The initial cloud is assumed to be uniformly rotating, centrally condensed sphere with rho is proportional to r(exp -2). Several cases are considered, in which both the overall clump size and the total amount of angular momentum are varied. The calculations show how a warm, quasi-hydrostatic disk surrounding a central unresolved core of only a few solar masses forms and grows in size and mass. The disk is encased in two distinct accretion shock fronts, both of which are several scale heights above the equatorial plane. At the end of the calculation of our standard case, the central unresolved region is found to have a mass of 2.7 solar masses and a ratio of rotational to gravitational energy of approximately 0.45, sufficiently large to be unstable to nonaxisymmetric perturbations. In addition, the inner portions of the disk containing most of the mass are unstable according to the local Toomre criterion, implying that also in this region nonaxisymmetric perturbations will lead to rapid evolution. Under the assumption that gravitational torques would transport angular momentum out of this region, a central core of less than or approximately 8 solar masses with a stable disk of greater than or approximately = 2 solar masses should result. Frequency-dependent radiative transfer calculations of the standard case at selected ages show how the continuum spectrum of the structure depends on the disk's orientation and age and how the observed isophotal contours vary with wavelength. Because of the strong dependence on viewing angle, continuum spectra alone should not be used to estimate the evolutionary stage of development of these objects. Comparable results were obtained for the other cases considered.
We calculate numerically the collapse of slowly rotating, nonmagnetic, massive molecular clumps of masses 30,60, and 120 Stellar Mass, which conceivably could lead to the formation of massive stars. Because radiative acceleration on dust grains plays a critical role in the clump's dynamical evolution, we have improved the module for continuum radiation transfer in an existing two-dimensional (axial symmetry assumed) radiation hydrodynamic code. In particular, rather than using "gray" dust opacities and "gray" radiation transfer, we calculate the dust's wavelength-dependent absorption and emission simultaneously with the radiation density at each wavelength and the equilibrium temperatures of three grain components: amorphous carbon particles. silicates, and " dirty ice " -coated silicates. Because our simulations cannot spatially resolve the innermost regions of the molecular clump, however, we cannot distinguish between the formation of a dense central cluster or a single massive object. Furthermore, we cannot exclude significant mass loss from the central object(s) that may interact with the inflow into the central grid cell. Thus, with our basic assumption that all material in the innermost grid cell accretes onto a single object. we are able to provide only an upper limit to the mass of stars that could possibly be formed. We introduce a semianalytical scheme for augmenting existing evolutionary tracks of pre-main-sequence protostars by including the effects of accretion. By considering an open outermost boundary, an arbitrary amount of material could, in principal, be accreted onto this central star. However, for the three cases considered (30, 60, and 120 Stellar Mass originally within the computation grid), radiation acceleration limited the final masses to 3 1.6, 33.6, and 42.9 Stellar Mass, respectively, for wavelength-dependent radiation transfer and to 19.1, 20.1, and 22.9 Stellar Mass. for the corresponding simulations with gray radiation transfer. Our calculations demonstrate that massive stars can in principle be formed via accretion through a disk. The accretion rate onto the central source increases rapidly after one initial free-fall time and decreases monotonically afterward. By enhancing the nonisotropic character of the radiation field, the accretion disk reduces the effects of radiative acceleration in the radial direction - a process we call the "flashlight effect." The flashlight effect is further amplified in our case by including the effects of frequency-dependent radiation transfer. We conclude with the warning that a careful treatment of radiation transfer is a mandatory requirement for realistic simulations of the formation of massive stars.
Shock breakout emission is among the first observable signals in a wide variety of astrophysical phenomena, including neutron star (NS) mergers, and it can be the dominant component in low-luminosity short γ-ray bursts, as exemplified by GRB 170817A. In this work, we investigate the cocoon shock breakout emission in NS mergers and how its signal depends on the outermost layers of the ejecta profile, which we derive from general relativistic radiation hydrodynamic simulations. We study the formation of the cocoon as a consequence of a relativistic jet propagating through the ejecta. To explore the influence of the outermost layers of the ejecta on the breakout emission, we explore cases where the ejecta has a sharp cutoff or an extended smooth tail. We find that the shock breakout emission is strongly influenced by the shape of the ejecta outer layers, with extended tails yielding results consistent with the observed properties of GRB 170817A, whereas sharp cutoffs overestimate the radiated energy. Using a Bayesian analysis, we estimate the best fit parameters of the central engine, considering both accreting black hole and magnetized neutron star scenarios. Our findings indicate a slight preference for the scenarios where the engine is a black hole. Furthermore, our work probes the nature of neutron star mergers and highlights the importance of the shape of the ejecta profile in modeling early electromagnetic counterparts to these mergers.