Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stefan model”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Chapmanite [Fe 2 Sb(Si 2 O 5 )O 3 (OH)]: thermodynamic properties and formation in low-temperature environments

Abstract. of synthetic Sb 2 O 5 , MgSb 2 O 6 (analogue of the mineral byströmite), Mg[Sb(OH) 6 ] 2 ∙6H 2 O (brandholzite), and natural chapmanite [(Fe 1.88 Al 0.12 )Sb(Si 2 O 5 )O 3 (OH)]. Enthalpies of reactions, including formation enthalpies, were evaluated using reference compounds Sb, Sb 2 O 3 , Sb 2 O 5 , and other phases, with high-temperature oxide melt solution calorimetry in lead borate and sodium molybdate solvents. Heat capacity and entropy were determined by relaxation and differential scanning calorimetry. The best set of Δ f H o (kJ mol -1 ) and S o (J mol -1 K -1 ) is byströmite -1733.0±3.6, 139.3±1.0; brandholzite -5243.1±3.6, 571.0±4.0; and chapmanite -3164.9±4.7, 305.1±2.1. The data for chapmanite give Δ f G o of -2973.6±4.7 kJ mol -1 and log K=-17.10 for the dissolution reaction (Fe 1.88 Al 0.12 )Sb(Si 2 O 5 )O 3 (OH) + 6H + → 1.88Fe 3+ + 0.12Al 3+ + 2SiO$_2^0$ + Sb(OH)$_3^0$ + 2H 2 O. Analysis of the data showed that chapmanite is finely balanced in terms of its stability with schafarzikite (FeSb 2 O 4 ) and tripuhyite (FeSbO 4 ) under a specific, narrow range of conditions when both aqueous Fe(III) and Sb(III) are abundant. In such a model, chapmanite is metastable by a narrow margin but could be stabilized by high SiO$_2^0$(aq) activities. Natural assemblages of chapmanite commonly contain abundant amorphous silica, suggesting that this mechanism may be indeed responsible for the formation of chapmanite. Chapmanite probably forms during low-temperature hydrothermal overprint of pre-existing Sb ores under moderately reducing conditions; the slightly elevated temperatures may help to overcome the kinetic barrier for its crystallization. During weathering, sheet silicates may adsorb Sb 3+ in tridentate hexanuclear fashion, thus exposing their chapmanite-like surfaces to the surrounding aqueous environment. Formation of chapmanite, as many other sheet silicates, under ambient conditions, is unlikely.

58 GEOSCIENCES↗
DOE OSTI2024

A portable parton-level event generator for the high-luminosity LHC

The rapid deployment of computing hardware different from the traditional CPU+RAM model in data centers around the world mandates a change in the design of event generators for the Large Hadron Collider, in order to provide economically and ecologically sustainable simulations for the high-luminosity era of the LHC. Parton-level event generation is one of the most computationally demanding parts of the simulation and is therefore a prime target for improvements. We present a production-ready leading-order parton-level event generation framework capable of utilizing most modern hardware and discuss its performance in the standard candle processes of vector boson and top-quark pair production with up to five additional jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

PAMTRA 1.0: the Passive and Active Microwave radiative TRAnsfer tool for simulating radiometer and radar measurements of the cloudy atmosphere

Forward models are a key tool to generate synthetic observations given knowledge of the atmospheric state. In this way, they are an integral part of inversion algorithms that aim to retrieve geophysical variables from observations or in data assimilation. Their application for the exploitation of the full information content of remote sensing observations becomes increasingly important when these are used to evaluate the performance of cloud-resolving models (CRMs). Herein, CRM profiles or fields provide the input to the forward model whose simulation results are subsequently compared to the observations. This paper introduces the freely available comprehensive microwave forward model PAMTRA (Passive and Active Microwave TRAnsfer), demonstrates its capabilities to simulate passive and active measurements across the microwave spectral region for upward- and downward-looking geometries, and illustrates how the forward simulations can be used to evaluate CRMs and to interpret measurements to improve our understanding of cloud processes. PAMTRA is unique as it treats passive and active radiative transfer (RT) in a consistent way with the passive forward model providing upwelling and downwelling polarized brightness temperatures and radiances for arbitrary observation angles. The active part is capable of simulating the full radar Doppler spectrum and its moments. PAMTRA is designed to be flexible with respect to instrument specifications and interfaces to many different formats of input and output, especially CRMs, spanning the range from bin-resolved microphysical output to one- and two-moment schemes, and to in situ measured hydrometeor properties. A specific highlight is the incorporation of the self-similar Rayleigh–Gans approximation (SSRGA) for both active and passive applications, which becomes especially important for the investigation of frozen hydrometeors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Diazo Linker Ligand Promotes Flexibility and Induced Fit Binding in a Microporous Copper Coordination Network

Abstract Flexible organic linkers represent an intuitive and effective strategy to design flexible metal–organic materials. We report herein a systematic study concerning the effect of varying the central bond of mixed pyridyl‐benzoate linkers, L, upon the flexibility of three isostructuralkddtopology microporous coordination networks (CNs) of formula ML 2 :X‐kdd‐1‐Cu,1= L = (E)‐4‐(pyridin‐4‐yldiazenyl)benzoate;X‐kdd‐2‐Cu,2= L = (E)‐4‐(2‐(pyridin‐4‐yl)vinyl)benzoate; the previously reportedX‐kdd‐3‐Cu,3= L = 4‐(pyridin‐4‐ylethynyl)benzoate. As revealed by single crystal x‐ray diffraction (SCXRD) and gas sorption studies,X‐kdd‐1‐Cu, exhibited gate‐opening during CO 2 and hydrocarbon (C2 and C8) sorption experiments whereas the other two CNs did not. Insight into these phase transformations was gained from in situ variable‐pressure and variable temperature powder X‐ray diffraction (PXRD), SCXRD, and modeling. Rotation of ligand1around the diazo bond, torsion angle changes between phenyl and carboxylate moieties, and deformation of the Cu‐based rod building blocks enabled activatedX‐kdd‐1‐Cuto form new phases with C8 isomers and CH 2 Cl 2 , CH 2 Cl 2 inducing contraction of the activated phase. Computational studies suggest that1enables flexibility thanks to its lower barrier of deformation versus2or3. This study teaches that diazo moieties could offer a general strategy to enhance the flexibility of CNs.

Chemistry↗
DOE OSTI2024

Sparse-grid discontinuous Galerkin methods for the Vlasov–Poisson–Lenard–Bernstein model

Sparse-grid methods have recently gained interest in reducing the computational cost of solving high-dimensional kinetic equations. In this paper, we construct adaptive and hybrid sparse-grid methods for the Vlasov–Poisson–Lenard–Bernstein (VPLB) model. This model has applications to plasma physics and is simulated in two reduced geometries: a 0x3v space homogeneous geometry and a 1x3v slab geometry. Here we use the discontinuous Galerkin (DG) method as a base discretization due to its high-order accuracy and ability to preserve important structural properties of partial differential equations. We utilize a multiwavelet basis expansion to determine the sparse-grid basis and the adaptive mesh criteria. We analyze the proposed sparse-grid methods on a suite of three test problems by computing the savings afforded by sparse-grids in comparison to standard solutions of the DG method. The results are obtained using the adaptive sparse-grid discretization library ASGarD.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Joint galaxy–galaxy lensing and clustering constraints on galaxy formation

ABSTRACT We compare predictions for galaxy–galaxy lensing profiles and clustering from the Henriques et al. public version of the Munich semi-analytical model (SAM) of galaxy formation and the IllustrisTNG suite, primarily TNG300, with observations from KiDS + GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, and stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast \gt 10^{11}\, \mathrm{ M}_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode active galactic nucleus accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high-mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only (χ2 = 1.81 compared to χ2 = 7.79 for the SAM) and locally brightest galaxy samples (χ2 = 3.80 compared to χ2 = 5.01). With added dust corrections to the colours it matches the SDSS clustering signal well for red low-mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\, {{\ \rm per\ cent}}$ excessive lensing signals for intermediate-mass red galaxies with 10.2 < log10M*[M⊙] < 11.2 at $r \approx 0.6\, h^{-1}\, \text{Mpc}$, which require further theoretical development.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Biochemical and structural insights into SARS-CoV-2 polyprotein processing by Mpro

SARS-CoV-2, a human coronavirus, is the causative agent of the COVID-19 pandemic. Its genome is translated into two large polyproteins subsequently cleaved by viral papain-like protease and main protease (Mpro). Polyprotein processing is essential yet incompletely understood. We studied Mpro-mediated processing of the nsp7-11 polyprotein, whose mature products include cofactors of the viral replicase, and identified the order of cleavages. Integrative modeling based on mass spectrometry (including hydrogen-deuterium exchange and cross-linking) and x-ray scattering yielded a nsp7-11 structural ensemble, demonstrating shared secondary structural elements with individual nsps. The pattern of cross-links and HDX footprint of the C145A Mpro and nsp7-11 complex demonstrate preferential binding of the enzyme active site to the polyprotein junction sites and additional transient contacts to help orient the enzyme on its substrate for cleavage. Last, proteolysis assays were used to characterize the effect of inhibitors/binders on Mpro processing/inhibition using the nsp7-11 polyprotein as substrate.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗
DOE OSTI2025

Higgs production via vector-boson fusion at the LHC

In this article, we summarise the recent experimental measurements and theoretical work on Higgs boson production via vector-boson fusion at the LHC. Along with this, we provide state-of-the-art predictions at fixed order as well as with parton-shower corrections within the Standard Model at 13.6 TeV. The results are presented in the form of multi-differential distributions as well as in the Simplified Template Cross Section bins. All materials and outputs of this study are available on public repositories. Finally, following findings in the literature, recommendations are made to estimate theoretical uncertainties related to parton-shower corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Machine learning the microscopic form of nematic order in twisted double-bilayer graphene

Modern scanning probe techniques, such as scanning tunneling microscopy, provide access to a large amount of data encoding the underlying physics of quantum matter. In this work, we show how convolutional neural networks can be used to learn effective theoretical models from scanning tunneling microscopy data on correlated moiré superlattices. Moiré systems are particularly well suited for this task as their increased lattice constant provides access to intra-unit-cell physics, while their tunability allows for the collection of high-dimensional data sets from a single sample. Using electronic nematic order in twisted double-bilayer graphene as an example, we show that incorporating correlations between the local density of states at different energies allows convolutional neural networks not only to learn the microscopic nematic order parameter, but also to distinguish it from heterostrain. These results demonstrate that neural networks are a powerful method for investigating the microscopic details of correlated phenomena in moiré systems and beyond.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Hyperglycemia alters mitochondrial respiration efficiency and mitophagy in human podocytes

Highlights: • Impairment of mitochondrial dynamics underlies the pathogenesis of type 2 diabetes. • Mitochondrial respiration is altered in glomeruli isolated from diabetic STZ-rats. • High glucose induces mitochondrial fission and fragmentation in human podocytes. • Mitochondrial fusion, biogenesis and mitophagy are reduced in HG-treated podocytes. • HG can injure the kidney by alteration of podocytes and glomerular bioenergetics. Podocytes constitute the outer layer of the renal glomerular filtration barrier. Their energy requirements strongly depend on efficient oxidative respiration, which is tightly connected with mitochondrial dynamics. We hypothesized that hyperglycemia modulates energy metabolism in glomeruli and podocytes and contributes to the development of diabetic kidney disease. We found that oxygen consumption rates were severely reduced in glomeruli from diabetic rats and in human podocytes that were cultured in high glucose concentration (30 mM; HG). In these models, all of the mitochondrial respiratory parameters, including basal and maximal respiration, ATP production, and spare respiratory capacity, were significantly decreased. Podocytes that were treated with HG showed a fragmented mitochondrial network, together with a decrease in expression of the mitochondrial fusion markers MFN1, MFN2, and OPA1, and an increase in the activity of the fission marker DRP1. We showed that markers of mitochondrial biogenesis, such as PGC-1α and TFAM, decreased in HG-treated podocytes. Moreover, PINK1/parkin-dependent mitophagy was inhibited in these cells. These results provide evidence that hyperglycemia impairs mitochondrial dynamics and turnover, which may underlie the remarkable deterioration of mitochondrial respiration parameters in glomeruli and podocytes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Understanding the structure and mechanism of Na + diffusion in NASICON solid-state electrolytes and the effect of Sc- and Al/Y-substitution

NASICON (sodium superionic conductor) based ceramics are one of the most promising classes of solid-state electrolytes for all-solid-state batteries. However, the mechanism of sodium ion diffusion is not understood in great detail since there is still a discrepancy between reported average structure models, local structures, and the number and position of sodium sites. To close this gap, we investigate the underlying diffusion mechanism and structural changes governing the Na + transport in Na 3.4 Zr 2 Si 2.4 P 0.6 O 12 using quasielastic neutron scattering (QENS) and powder X-ray diffraction (XRD). In the temperature range from 298 K to 640 K, the correlations between structural changes of a monoclinic C2/c to rhombohedral R $\bar{3}$c phase transition and the result of ion diffusion are investigated. The analysis of the quasielastic neutron scattering data reveals two quasielastic components corresponding to the Chudley-Elliott jump-diffusion model. It clearly shows two different Na + diffusion processes, local and long-range, on two different time and length scales and allows calculations of their corresponding activation energies. Additionally, the effects of Sc 3+ and Al 3+ /Y 3+ aliovalent substitution of Zr 4+ ions on the crystal structure and Na + diffusion are also studied. We can distinguish a local, chain, and cross-chain diffusion mechanism based on correlated QENS and XRD comparison of relevant nearest crystallographic Na–Na distances. The results reveal that the Na + diffusion in these NASICONs is three-dimensional and can provide guidelines on how dopants and changes in the crystal structure can affect the Na + conductivity.

25 ENERGY STORAGE↗
DOE OSTI2022

Boron nitride on SiC(0001)

Here, in the field of van der Waals heterostructures, the twist angle between stacked two-dimensional layers has been identified to be of utmost importance for the properties of the heterostructures. In this context, we previously reported the growth of a single layer of unconventionally oriented epitaxial graphene that forms in a surfactant atmosphere [F. C. Bocquet et al., Phys. Rev. Lett. 125, 106102 (2020)]. The resulting G-R0° layer is aligned with the SiC lattice, and hence represents an important milestone towards high-quality twisted bilayer graphene, a frequently investigated model system in this field. Here, we focus on the surface structures obtained in the same surfactant atmosphere, but at lower preparation temperatures at which a boron nitride template layer forms on SiC(0001). In a comprehensive study based on complementary experimental and theoretical techniques, we find—in contrast to the literature—that this template layer is a hexagonal B x N y layer, but not high-quality hBN. It is aligned with the SiC lattice and gradually replaced by low-quality graphene in the 0° orientation of the B x N y template layer upon annealing.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Snow-eater heat waves of the western United States

Abrupt snowmelt, triggered by rain-on-snow events or "snow-eater heat waves," can cause flooding, initiate or accelerate snow drought, and affect water availability. However, the characteristics (e.g., area, duration, and frequency), impacts, and trends of snow-eater heat waves have received little attention. To address this gap, we developed a method to identify snow-eater heat waves and estimate their melt potential using 20th Century Reanalysis version 3 air temperature data, the TempestExtremes algorithm, and an operational snowmelt model (SNOW-17) across 1850-2015. Melt season snow-eater heat waves typically last 3 to 5 days, with three to five events, doubling snowmelt rates. Seven of 11 spring superfloods are shown to coincide with snow-eater heat waves. Since the 1850s, snow-eater heat waves have increased in area and frequency, decreased in duration, and shifted earlier in the melt season. Incorporating snow-eater heat-wave impacts into SNOW-17 enhances extreme melt estimates, improving water management support tools.

Rhoades, Alan M↗
DOE OSTI2025

When can we detect lianas from space? Toward a mechanistic understanding of liana‐infested forest optics

Abstract Lianas, woody vines acting as structural parasites of trees, have profound effects on the composition and structure of tropical forests, impacting tree growth, mortality, and forest succession. Remote sensing could offer a powerful tool for quantifying the scale of liana infestation, provided the availability of robust detection methods. We analyze the consistency and global geographic specificity of spectral signals—reflectance across wavelengths—from liana‐infested tree crowns and forest stands, examining the underlying mechanisms of these signals. We compiled a uniquely comprehensive database, including leaf reflectance spectra from 5424 leaves, fine‐scale airborne reflectance data from 999 liana‐infested canopies, and coarse‐scale satellite reflectance data covering 775 ha of liana‐infested forest stands. To unravel the mechanisms of the liana spectral signal, we applied mechanistic radiative transfer models across scales, establishing a synthesis of the relative importance of different mechanisms, which we corroborate with field data on liana leaf chemistry and canopy structure. We find a consistent liana spectral signal at canopy and stand scales across globally distributed sites. This signature mainly arises at the canopy level due to direct effects of more horizontal leaf angles, resulting in a larger projected leaf area, and indirect effects from increased light scattering in the near and short‐wave infrared regions, linked to lianas' less costly leaf construction compared with trees on average. The existence of a consistent global spectral signal for lianas suggests that large‐scale quantification of liana infestation is feasible. However, because the traits responsible for the liana canopy‐reflectance signal are not exclusive to lianas, accurate large‐scale detection requires rigorously validated remote sensing methods. Our models highlight challenges in automated detection, such as potential misidentification due to leaf phenology, tree life history, topography, and climate, especially where the scale of liana infestation is less than a single remote sensing pixel. The observed cross‐site patterns also prompt ecological questions about lianas' adaptive similarities in optical traits across environments, indicating possible convergent evolution due to shared constraints on leaf biochemical and structural traits.

Environmental Sciences & Ecology↗
DOE OSTI2021

Kinetics of Hydride Transfer from Catalytic Metal-Free Hydride Donors to CO2

Selective reduction of CO 2 to formate represents an ongoing challenge in photoelectrocatalysis. To provide mechanistic insights, we investigate the kinetics of hydride transfer (HT) from a series of metal-free hydride donors to CO 2 . The observed dependence of experimental and calculated HT barriers on the thermodynamic driving force was modeled by using the Marcus hydride transfer formalism to obtain the insights into the effect of reorganization energies on the reaction kinetics. Our results indicate that even if the most ideal hydride donor were discovered, the HT to CO 2 would exhibit sluggish kinetics (<100 turnovers per second at -0.1 eV driving force), indicating that the conventional HT may not be an appropriate mechanism for solar conversion of CO 2 to formate. We propose that the conventional HT mechanism should not be considered for CO 2 reduction catalysis and argue that the orthogonal HT mechanism, previously proposed to address thermodynamic limitations of this reaction, may also lead to lower kinetic barriers for CO 2 reduction to formate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Fast and flexible neutrino decoupling. Part I. The Standard Model

Cosmological determinations of the number of relativistic neutrino species, N eff , are becoming increasingly accurate, and further improvements are expected both from CMB and BBN data. Given this context, we update the evaluation of N eff and the current entropy density via the momentum-averaged approach. This allows for a numerically fast description of neutrino decoupling, easily portable to an array of new physics scenarios. We revisit all aspects of this approach, including collision terms with full electron mass dependence, finite temperature QED corrections to the equation of state, neutrino oscillations, and the modelling of neutrino ensembles with effective chemical potentials. For integrated observables, our results differ by less than 0.04% from the solution of the momentum-dependent evolution equation. We outline how to extend the approach to BSM settings, and will highlight its power in Part II. To facilitate the practical implementation, we release a Mathematica and python code within nudec_BSM_v2, easily linkable to BBN codes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗