Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gravitational collapse”

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.

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multimessenger signals of heavy axionlike particles in core-collapse supernovae: Two-dimensional simulations

Core-collapse supernovae are a useful laboratory to probe the nature of exotic particles. If axionlike particles (ALPs) are produced in supernovae, they can affect the transfer of energy and leave traces in observational signatures. In this work, we present results from two-dimensional supernova models including the effects of the production and the absorption of ALPs that couple with photons. It is found that the additional heating induced by ALPs can enhance the diagnostic energy of explosion, E diag . For example, for moderate ALP-photon coupling, we find explosion energies ~0.6×10 51 erg compared to our reference model without ALPs of ~0.4×10 51 erg in the first ~0.5 s postbounce explored in this work. Further, our findings indicate that when the coupling constant is sufficiently high, the neutrino luminosities and mean energies are decreased because of the additional cooling of the proto-neutron star via ALPs. The gravitational wave amplitude is also reduced because the mass accretion on the proto-neutron star is suppressed. Although the ALP-photon coupling can foster explodability, including enhancing the explosion energy closer to recent observations, more long-term simulations in three spatial dimensions are needed to draw robust conclusions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

The formation phase of the solar nebula

Hydrodynamical calculations of the collapse of an axisymmetric, rotating protostellar cloud, with radiation transport and without magnetic fields, are presented. The collapse is assumed to start from a centrally condensed sphere of radius 5 x 10 to the 15th cm, a mean density of 4 x 10 to the -15th g/cu cm, a total mass of 1 solar mass, and a total angular momentum of 10 to the 53rd g sq cm per sec. The numerical grid is chosen to resolve the region of disk formation between 1 and 60 AU from the center. Frequency-dependent radiative transfer calculations show how the emergent spectrum of the structure depends upon viewing angle with respect to the rotation axis and how the observed isophotal contours should depend on wavelength and viewing angle. The central part of the protostar, interior to 1 AU, is not resolved numerically but is modeled approximately. At the end of the calculation, this region is found to have a mass of 0.6 solar mass and a ratio of rotational to gravitational energy of about 0.4, sufficiently large to be unstable to nonaxisymmetric perturbations. Although the disk is gravitationally stable according to the local Toomre criterion, the nonaxisymmetric structure in the center is likely to lead to angular momentum transport.

Bodenheimer, Peter↗
DOE OSTI2020

BSE versus StarTrack: Implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences

As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events. In this work, we present new developments in the stellar-remnant formation and related schemes of the current N-body evolution program NBODY7. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or BSE of the NBODY7 code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised StarTrack population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. All the above newly implemented schemes are also shared with the standalone versions of SSE and BSE. Furthermore, all these demonstrations are performed with both the updated standalone BSE and the updated NBODY7/BSE. When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing. Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Scattering Amplitudes and QFT Insight (Final Scientific Report)

The calculation of scattering amplitudes provides an invariant window into the physical dynamic content of relativistic quantum field theory. Yet with traditional methods, for empirically relevant theories, these calculations scale with a factorial complexity in external particles and precision. If we aspire to collapse the theoretical uncertainty obscuring new physics hidden at all scales from the microscopic probed at high energy colliders to the largest effective field theory in the universe governing the evolution of large scale structure, this complexity challenge necessitates new ideas and methods in calculation. Novel approaches like the color-kinematics duality and the associated double-copy construction have drastically simplified the situation, relating both gauge and gravity theory predictions to a much smaller kernel of invariant kinematic data. This project supporting research by the PI's Amplitudes and Insights group at Northwestern University looked to push insight deep into both the IR and the UV by establishing how novel structures must constrain the predictions of counterterms in both gauge and gravity theories, setting the groundwork to exploring the high energy behavior of particular gravity theories via constituent gauge-theory calculations, and confronting fundamental challenges at the interface between QFT amplitudes analysis and next generation gravitational wave science as well as well as inflationary and large scale structure cosmology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Self-gravity and dissipation in polar rings

Studies of inclined rings inside galaxy potentials have mostly considered the influence of self-gravity and viscous dissipation separately. In this study, we construct models of highly inclined ('polar') rings in an external potential including both self-gravity and dissipation due to a drag force. We do not include pressure forces and thus ignore shock heating that dominates the evolution of gaseous rings inside strongly nonspherical potentials. We adopt an oblate spheroidal scale-free logarithmic potential with axis ratio q = 0.85 and an initial inclination of 80 deg for the self-gravitating rings. We find that stellar (dissipationless) rings suffer from mass loss during their evolution. Mass loss also drives a secular change of the mean inclination toward the poles of the potential. As much as half of the ring mass escapes in the process and forms an inner and an outer shell of precessing orbits. If the remaining mass is more than approximately 0.02 of the enclosed galaxy mass, rings remain bound and do not fall apart from differential precession. The rings precess at a constant rate for more than a precession period tau(sub p) finding the configuration predicted by Sparke in 1986 which warps at larger radii toward the poles of the potential. We model shear viscosity with a velocity-dependent drag force and find that nuclear inflow dominates over self-gravity if the characteristic viscous inflow time scale tau(sub vi) is shorter than approximately 25(tau(sub p)). Rings with (tau(sub vi))/(tau(sub p)) less than or approximately equal to 25 collapse toward the nucleus of the potential within one precession period independent of the amount of self-gravity. Our results imply that stars and gas in real polar rings exhibit markedly different dynamical evolutions.

Dubinski, John↗
NASA NTRS2005

Discrete element simulations of gravitational volcanic deformation: Deformation structures and geometries - 1

We have carried out two-dimensional particle dynamics simulations of granular piles subject to frictional Coulomb failure criteria to gain a first-order understanding of different modes of gravitational deformation within volcanoes. Under uniform basal and internal strength conditions, granular piles grow self-similarly, developing distinctive stratigraphies, morphologies, and structures. Piles constructed upon cohesive substrates exhibit particle avalanching, forming outward dipping strata and angle of repose slopes. Systematic decreases in basal strength lead to progressively deeper and steeper internal detachment faults and slip along a basal decollement; landslide forms grade from shallow slumps to deep-seated landslide and, finally, to axial subsidence and outward flank displacements, or volcanic spreading. Surface slopes decrease and develop concave up morphologies with decreasing decollement strength; depositional layers tilt progressively inward. Spatial variations in basal strength cause lateral transitions in pile structure, stratigraphy, and morphology. This approximation of volcanoes as Coulomb granular piles reproduces the richness of deformational structures and surface morphologies in many volcanic settings. The gentle slopes of Hawaiian volcanoes and Olympus Mons on Mars suggest weak basal decollements that enable volcanic spreading. High-angle normal faults, favored above weak decollements, are interpreted in both settings and may explain catastrophic sector collapse in Hawaii and broad aureole deposits surrounding Olympus Mons. In contrast, steeper slopes and shallow detachment faults predominate in the Canary Islands, thought to lack a weak decollement, favoring smaller, more frequent slope failures than predicted for Hawaii. The numerical results provide a useful predictive tool for interpreting dynamic behavior and associated geologic hazards of active volcanoes.

Morgan, Julia K.↗
DOE OSTI2025

Explaining Nonmerger Gamma-Ray Bursts and Broad-lined Supernovae with Close Binary Progenitors with Black Hole Central Engines

For over 25 yr, the origin of long-duration gamma-ray bursts (lGRBs) has been linked to the collapse of rotating massive stars. However, we have yet to pinpoint the stellar progenitor powering these transients. Moreover, the dominant engine powering the explosions remains open to debate. Observations of both lGRBs, supernovae associated with these GRBs, such as broad-line (BL) stripped-envelope (type Ic) supernovae (hereafter, Ic-BL), supernovae (SNe), and perhaps superluminous SNe, fast blue optical transients, and fast x-ray transients, may provide clues to both engines and progenitors. In this paper, we conduct a detailed study of the tight-binary formation scenario for lGRBs, comparing this scenario to other leading progenitor models. Combining this progenitor scenario with different lGRB engines, we can compare to existing data and make predictions for future observational tests. We find that the combination of the tight-binary progenitor scenario with the black hole accretion disk engine can explain lGRBs, low-luminosity GRBs, ultra-long GRBs, and Ic-BL. We discuss the various progenitor properties required for these different subclasses and note such systems would be future gravitational-wave merger sources. We show that the current literature on other progenitor-engine scenarios cannot explain all of these transient classes with a single origin, motivating additional work. We find that the tight-binary progenitor with a magnetar engine is excluded by existing observations. The observations can be used to constrain the properties of stellar evolution, the nature of the GRB, and the associated SN engines in lGRBs and Ic-BL. We discuss the future observations needed to constrain our understanding of these rare, but powerful, explosions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Insights into cosmological structure formation with machine learning.

Our modern understanding of cosmological structure formation posits that small matter density fluctuations present in the early Universe, as traced by the cosmic microwave background, grow via gravitational instability to form extended haloes of dark matter. A theoretical understanding of the structure, evolution and formation of dark matter haloes is an essential step towards unravelling the intricate connection between halo and galaxy formation, needed to test our cosmological model against data from upcoming galaxy surveys. Physical understanding of the process of dark matter halo formation is made difficult by the highly non-linear nature of the haloes' evolution. I describe a new approach to gain physical insight into cosmological structure formation based on machine learning. This approach combines the ability of machine learning algorithms to learn non-linear relationships, with techniques that enable us to physically interpret the learnt mapping. I describe applicatio ns of th e method, with the aim of investigating which aspects of the early universe density field impact the later formation of dark matter haloes. First I present a case where the process of halo formation is turned into a binary classification problem; the algorithm predicts whether or not dark matter `particles' in the initial conditions of a simulation will collapse into haloes of a given mass range. Second, I present its generalization to regression, where the algorithm infers the final mass of the halo to which each particle will later belong. I show that the initial tidal shear does not play a significant role compared to the initial density field in establishing final halo masses. Finally, I demonstrate that extending the framework to deep learning algorithms such as convolutional neural networks allows us to explore connections between the early universe and late time haloes beyond those studied by existing analytic approximations of halo collapse.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Paramagnetic Liquid Bridge in a Gravity-Compensating Magnetic Field

Magnetic levitation was used to stabilize cylindrical columns of a paramagnetic liquid in air between two solid supports. The maximum achievable length to diameter ratio R(sub max) was approx. (3.10 +/- 0.07), very close to the Rayleigh-Plateau limit of pi. For smaller R, the stability of the column was measured as a function of the Bond number, which could be continuously varied by adjusting the strength of the magnetic field. Liquid bridges supported by two solid surfaces have been attracting scientific attention since the time of Rayleigh and Plateau. For a cylindrical bridge of length L and diameter d, it was shown theoretically that in zero gravity the maximum slenderness ratio R (identically = L/d) is pi. The stability and ultimate collapse of such bridges is of interest because of their importance in a number of industrial processes and their potential for low gravity applications. In the presence of gravity, however, the cylindrical shape of an axisymmetric bridge tends to deform, limiting its stability and decreasing the maximum achievable value of R. Theoretical studies have discussed the stability and possible shapes of axisymmetric bridges. Experiments typically are performed in either a Plateau tank, in which the bridge is surrounded by a density-matched immiscible fluid, or in a space-borne microgravity environment. It has been shown, for example, that the stability limit R can be pushed beyond pi by using flow stabilization, by acoustic radiation pressure, or by forming columns in the presence of an axial electric field. In this work, magnetic levitation was used to simulate a low gravity environment and create quasi-cylindrical liquid columns in air. Use of a magnetic field permits us to continuously vary the Bond number B identically equal to (g)(rho)d(exp 2)/4(sigma), where g is the gravitational acceleration, rho is the density of the liquid, and sigma is the surface tension of the liquid in air. The dimensionless Bond number represents the relative importance of external forces acting on the liquid column to those due to surface tension. Our central result is that in a large magnetic field gradient we could create and stabilize columns of mixtures of water and paramagnetic manganese chloride tetrahydrate (MnCl2.4H2O), achieving a length to diameter ratio very close to pi.

Mahajan, Milind P.↗
NASA NTRS2013

Hubble Space Telescope Hx Imaging of Star-forming Galaxies at z approximately equal to 1-1.5: Evolution in the Size and Luminosity of Giant H II Regions

We present Hubble Space Telescope/Wide Field Camera 3 narrow-band imaging of the Hα emission in a sample of eight gravitationally lensed galaxies at z = 1-1.5. The magnification caused by the foreground clusters enables us to obtain a median source plane spatial resolution of 360 pc, as well as providing magnifications in flux ranging from approximately 10× to approximately 50×. This enables us to identify resolved star-forming HII regions at this epoch and therefore study their Hα luminosity distributions for comparisons with equivalent samples at z approximately 2 and in the local Universe. We find evolution in the both luminosity and surface brightness of HII regions with redshift. The distribution of clump properties can be quantified with an HII region luminosity function, which can be fit by a power law with an exponential break at some cut-off, and we find that the cut-off evolves with redshift. We therefore conclude that 'clumpy' galaxies are seen at high redshift because of the evolution of the cut-off mass; the galaxies themselves follow similar scaling relations to those at z = 0, but their HII regions are larger and brighter and thus appear as clumps which dominate the morphology of the galaxy. A simple theoretical argument based on gas collapsing on scales of the Jeans mass in a marginally unstable disc shows that the clumpy morphologies of high-z galaxies are driven by the competing effects of higher gas fractions causing perturbations on larger scales, partially compensated by higher epicyclic frequencies which stabilize the disc.

H imaging↗
DOE OSTI2020

Newton versus the machine: solving the chaotic three-body problem using deep neural networks

ABSTRACT Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system’s chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.

Breen, Philip G.↗
NASA NTRS2020

Properties and Astrophysical Implications of the 150 M Binary Black Hole Merger GW190521

The gravitational-wave signal GW190521 is consistent with a binary black hole (BBH) merger source at redshift 0.8 with unusually high component masses, 85 (+21,-14) Mꙩ and 66 (+17,-18) Mꙩ, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range 65–120 Mꙩ. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger (142 (+28,-16) Mꙩ) classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular BBH coalescence, we detail the physical properties of GW190521’s source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of the coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be 0.13 (+0.30,-0.11)Gp/cu.c yr. We discuss the astrophysical implications of GW190521 for stellar collapse and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescences, or via hierarchical mergers of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

R. Abbott↗
DOE OSTI2023

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO’s and Advanced Virgo’s third observing run

ABSTRACT We describe a search for gravitational waves from compact binaries with at least one component with mass $0.2$–$1.0 \, \mathrm{M}_\odot$ and mass ratio q ≥ 0.1 in Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and Advanced Virgo data collected between 2019 November 1, 15:00 utc and 2020 March 27, 17:00 utc. No signals were detected. The most significant candidate has a false alarm rate of $0.2 \, \rm {yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO’s and Advanced Virgo’s third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim \, 0.6$ (at 90 per cent confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions, we are unable to rule out fPBH = 1. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound fDBH < 10−5 on the fraction of atomic dark matter collapsed into black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nucleosynthesis in outflows of compact objects and detection prospects of associated kilonovae

ABSTRACT We perform a comparative analysis of nucleosynthesis yields from binary neutron star (BNS) mergers, black hole-neutron star (BHNS) mergers, and core-collapse supernovae (CCSNe) with the goal of determining which are the most dominant sources of r-process enrichment observed in stars. We find that BNS and BHNS binaries may eject similar mass distributions of robust r-process nuclei post-merger (up to third peak and actinides, A ∼ 200−240), after accounting for the volumetric event rates. Magnetorotational (MR) CCSNe likely undergo a weak r-process (up to A ∼ 140) and contribute to the production of light element primary process (LEPP) nuclei, whereas typical thermal, neutrino-driven CCSNe only synthesize up to first r-process peak nuclei (A ∼ 80−90). We also find that the upper limit to the rate of MR CCSNe is $\lesssim 1~{{\ \rm per\ cent}}$ the rate of typical thermal CCSNe; if the rate was higher, then weak r-process nuclei would be overproduced. Although the largest uncertainty is from the volumetric event rate, the prospects are encouraging for confirming these rates in the next few years with upcoming surveys. Using a simple model to estimate the resulting kilonova light curve from mergers and our set of fiducial merger parameters, we predict that ∼7 BNS and ∼2 BHNS events will be detectable per year by the Vera C. Rubin Observatory (LSST), with prior gravitational wave (GW) triggers.

Astronomy & Astrophysics↗
DOE OSTI2022

Evidence for X-Ray Emission in Excess to the Jet-afterglow Decay 3.5 yr after the Binary Neutron Star Merger GW 170817: A New Emission Component

For the first ~3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm -3 . We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 10 38 erg s -1 at 1234 days, the recently detected X-ray emission represents a ≥ 3.2σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥ 3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥ 3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Dark matter halo mass functions and density profiles from mass and energy cascade

Abstract Halo abundance and structure play a central role for modeling structure formation and evolution. Without relying on a spherical or ellipsoidal collapse model, we analytically derive the halo mass function and cuspy halo density (inner slope of −4/3) based on the mass and energy cascade theory in dark matter flow. The hierarchical halo structure formation leads to halo or particle random walk with a position-dependent waiting time $$\tau _g$$ τ g . First, the inverse mass cascade from small to large scales leads to the halo random walk in mass space with $$\tau _g\propto m_h^{-\lambda }$$ τ g ∝ m h - λ , where $$m_h$$ m h is the halo mass and $$\lambda$$ λ is a halo geometry parameter with predicted value of 2/3. The corresponding Fokker-Planck solution for halo random walk in mass space gives rise to the halo mass function with a power-law behavior on small scale and exponential decay on large scale. This can be further improved by considering two different $$\lambda$$ λ for haloes below and above a critical mass scale $$m_h^*$$ m h ∗ , i.e. a double- $$\lambda$$ λ halo mass function. Second, a double- $$\gamma$$ γ density profile can be derived based on the particle random walk in 3D space with a position-dependent waiting time $$\tau _g \propto \Phi (r)^{-1} \propto r^{-\gamma }$$ τ g ∝ Φ ( r ) - 1 ∝ r - γ , where $$\Phi$$ Φ is the gravitational potential and r is the particle distance to halo center. Theory predicts $$\gamma =2/3$$ γ = 2 / 3 that leads to a cuspy density profile with an inner slope of −4/3, consistent with the predicted scaling laws from energy cascade. The Press-Schechter mass function and Einasto density profile are just special cases of proposed models. The small scale permanence can be identified due to the scale-independent rate of mass and energy cascade, where density profiles of different halo masses and redshifts converge to the $$-4/3$$ - 4 / 3 scaling law ( $$\rho _h \propto r^{-4/3}$$ ρ h ∝ r - 4 / 3 ) on small scales. Theory predicts the halo number density scales with halo mass as $$\propto m_h^{-1.9}$$ ∝ m h - 1.9 , while the halo mass density scales as $$\propto m_h^{4/9}$$ ∝ m h 4 / 9 . Results were compared against the Illustris simulations. This new perspective provides a theory for nearly universal halo mass functions and density profiles.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

The early evolution of protostellar disks

We consider the origin and intital growth of the disks that form around protostars during the collapse of rotating molecular cloud cores. These disks are assumed to be inviscid and pressure free, and to have masses small compared to those of their central stars. We find that there exist three distinct components-an outer disk, in which shocked gas moves with comparable azimuthal and radical velocities; and inner disk, where material follows nearly circular orbits, but spirals slowly toward the star because of the drag exerted by adjacent onfalling matter, and a turbulent ring adjoining the first two regions. Early in the evolution, i.e., soon after infalling matter begins to miss the star, only the outer disk is present, and the total mass acceration rate onto the protostar is undiminished. Once the outer disk boundary grows to more than 2.9 times the stellar radius, first the ring, and then the inner disk appear. Thereafter, the radii of all three components expand as t(exp 3). The mass of the ring increase with time and is always 13% of the total mass that has fallen from the cloud. Concurrently with the buildup of the inner disk and ring, the accretion rate onto the star falls off. However, the protostellar mass continue to rise, asymptotically as t(exp 1/4). We calculated the radiated flux from the inner and outer disk components due to the release of gravitational potential energy. The flux from the inner disk is dominant and rises steeply toward the stellar surface. We also determine the surface temperature of the inner disk as a function of radius. The total disk luminosity decreases slowly with time, while the contributions from the ring and inner disk both fall as t(exp -2).

Stahler, Steven W.↗
DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗