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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Distinguishing prompt-collapse binary neutron star mergers from binary black Holes: Tidal effects and remnant properties

We study the properties of remnants formed in prompt-collapse binary neutron star mergers. We consider nonspinning neutron star binaries over a range of total masses and mass ratios across a set of 22 equations of state, totaling 107 numerical relativity simulations. We report the final mass and spin of the systems (including the accretion disk and ejecta) to be constrained in a narrow range—0.98 ≲ 𝑀 𝑓 /𝑀 ≲ 0.99 for the mass and 0.85 ≲ 𝑎 𝑓 ≲ 0.95 for the dimensionless spin—regardless of the binary configuration and matter effects. This sets them apart from binary black hole merger remnants. We assess the detectability of the postmerger signal in a future 40 km Cosmic Explorer observatory and find that the signal-to-noise ratio in the postmerger of an optimally located and oriented binary at a distance of 100 Mpc can range from <1 to 8, depending on the binary configuration and equation of state, with a majority of them greater than 4 in the set of simulations that we consider. We also consider the distinguishability between prompt-collapse binary neutron star and binary black hole mergers with the same masses and spins. We find that Cosmic Explorer will be able to distinguish such systems primarily via the measurement of tidal effects in the late inspiral. Neutron star binaries with reduced tidal deformability $\tilde{Λ}$ as small as ∼ 3.5 can be identified up to a distance of 100 Mpc, while neutron star binaries with $\tilde{Λ}$ ∼ 22 can be identified to distances greater than 250 Mpc. This is larger than the distance up to which the postmerger will be visible. Finally, we discuss the possible implications of our findings for the equation of state of neutron stars from the gravitational wave event GW230529.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Revisiting constraints on asymmetric dark matter from collapse in white dwarf stars

The runaway collapse phase of a small dark matter cluster inside a white dwarf star encompasses a reversible stage, where heat can be transferred back and forth between nuclear and dark matter. Induced nuclear burning phases are stable and early carbon depletion undermines previous claims of type Ia supernova ignition. Instead, mini black holes are formed at the center of the star that either evaporate or accrete stellar material until a macroscopic sub-Chandrasekhar-mass black hole is formed. In the latter case, a 0.1 to 1 second lasting electromagnetic transient signal can be detected upon ejection of the white dwarf’s potential magnetic field. Binary systems that transmute to black holes and subsequently merge emit gravitational waves. Advanced LIGO should detect one such sub-Chandrasekhar binary black hole inspiral per year, while future Einstein telescope-like facilities will detect thousands per year. The effective spin parameter distribution is peaked at 0.2 and permits future studies to disentangle from primordial sub-Chandrasekhar black holes. Such signatures are compatible with current direct detection constraints, as well as with neutron star constraints in the case of bosonic dark matter, even though they remain in conflict with the fermionic case for part of the parameter space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Filling the black hole mass gap: Avoiding pair instability in massive stars through addition of nonnuclear energy

In standard stellar evolution, stars with masses ranging from approximately 150 to 240M ⊙ are expected to evolve to a pair instability supernova with no black hole (BH) remnant. This evolutionary behavior leads to a predicted gap in the black hole mass function from approximately 50 to 140M ⊙ . Yet the LIGO and Virgo Collaborations recently discovered black holes of masses 66M ⊙ and 85M ⊙ in the gravitational wave event GW190521. Here, we propose a new method to populate the BH mass gap. If an energy source is added throughout the star in addition to nuclear fusion, it is possible for the altered evolution to avoid the complete destruction of a pair instability supernova, and instead a BH remnant is left behind. An example of an extra energy source is dark matter annihilation within the star, but our results hold more generally. We show this phenomenon by exploring the effect of adding an energy source independent of temperature and density to a 180M ⊙ star, using the MESA one-dimensional stellar evolution software. If ~50% of the star’s energy is due to this new source, the star is capable of avoiding the pair instability entirely and evolving towards a core-collapse supernova and ultimately a BH remnant with mass ~ 120M ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Efficient method for estimating the time evolution of the proto-neutron star mass and radius from a supernova neutrino signal

In this paper, we present a novel method to estimate the time evolution of the proto-neutron star (PNS) structure from the neutrino signal in a core-collapse supernova (CCSN). Employing recent results from multidimensional CCSN simulations, we delve into a relation between the total emitted neutrino energy (TONE) and PNS mass/radius, and we find that they are strongly correlated with each other. We fit the relation by simple polynomial functions connecting the TONE to the mass and radius of the PNS as a function of time. By combining another fitting function representing the correlation between the TONE and the cumulative number of events at each neutrino observatory, the PNS mass and radius can be retrieved from purely observed neutrino data. We demonstrate retrievals of PNS mass and radius from mock data of the neutrino signal, and we assess the capability of our proposed method. While underlining the limitations of the method, we also discuss the importance of the joint analysis with the gravitational wave signal. This would reduce uncertainties of parameter estimations in our method, and may narrow down the possible neutrino oscillation model. The proposed method is a very easy and inexpensive computation, which will be useful in real data analysis of the CCSN neutrino signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Current Sheet Alignment in Oblique Black Hole Magnetospheres: A Black Hole Pulsar?

Abstract We study the magnetospheric evolution of a nonaccreting spinning black hole (BH) with an initially inclined split monopole magnetic field by means of 3D general relativistic magnetohydrodynamic simulations. This serves as a model for a neutron star (NS) collapse or a BH–NS merger remnant after the inherited magnetosphere has settled into a split monopole field creating a striped wind. We show that the initially inclined split monopolar current sheet aligns over time with the BH equatorial plane. The inclination angle evolves exponentially toward alignment, with an alignment timescale that is inversely proportional to the square of the BH angular velocity, where higher spin results in faster alignment. Furthermore, magnetic reconnection in the current sheet leads to exponential decay of event-horizon-penetrating magnetic flux with nearly the same timescale for all considered BH spins. In addition, we present relations for the BH mass and spin in terms of the period and alignment timescale of the striped wind. The explored scenario of a rotating, aligning, and reconnecting current sheet can potentially lead to multimessenger electromagnetic counterparts to a gravitational-wave event due to the acceleration of particles powering high-energy radiation, plasmoid mergers resulting in coherent radio signals, and pulsating emission due to the initial misalignment of the BH magnetosphere.

Astronomy & Astrophysics↗
DOE OSTI2020

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-14+21 M o˙ and 66-18+17 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-16+28 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.11+0.30 Gpc -3 yr -1 . 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.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
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↗
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↗
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.↗
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↗
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↗