Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random Effects”

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

Investigating Grain Structure and Microcracking in SiCf-SiCm Composites Using 4D-STEM

Silicon carbide (SiC) fiber-reinforced SiC matrix composites (SiCf-SiCm) are promising materials for accident-tolerant fuel (ATF) cladding in light water reactors. This paper, using four-dimensional scanning transmission electron microscopy (4D-STEM), studied the microstructure of SiCf-SiCm at nanoscale and provided grain statistics at the CVI/fiber region. Here, the results reveal that CVI region mainly contain irregular sub-micron grains (115 nm to 1800 nm) with a preferred orientation, while the fiber region has nano equiaxed grains (21 nm to 78 nm) with random orientations. Centered dark field imaging over the fiber regions verified the grain size measurement by 4D-STEM. Through 4D-STEM virtual dark field imaging and orientation mapping, intragranular propagation was identified as the fracture mechanism for a microcrack observed within the CVI region. The amorphous pyrolytic carbon layer was also shown to effectively arrest the microcrack.

4D-STEM↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Elucidating Primary Degradation Mechanisms in High-Cycling-Capacity, Compositionally Tunable High-Entropy Hydrides

The hydrogen sorption properties of single-phase bcc (TiVNb) 100–x Cr x alloys (x = 0–35) are reported. All alloys absorb hydrogen quickly at 25 °C, forming fcc hydrides with storage capacity depending on the Cr content. Here, a thermodynamic destabilization of the fcc hydride is observed with increasing Cr concentration, which agrees well with previous compositional machine learning models for metal hydride thermodynamics. The steric effect or repulsive interactions between Cr–H might be responsible for this behavior. The cycling performances of the TiVNbCr alloy show an initial decrease in capacity, which cannot be explained by a structural change. Pair distribution function analysis of the total X-ray scattering on the first and last cycled hydrides demonstrated an average random fcc structure without lattice distortion at short-range order. If the as-cast alloy contains a very low density of defects, the first hydrogen absorption introduces dislocations and vacancies that cumulate into small vacancy clusters, as revealed by positron annihilation spectroscopy. Finally, the main reason for the capacity drop seems to be due to dislocations formed during cycling, while the presence of vacancy clusters might be related to the lattice relaxation. Having identified the major contribution to the capacity loss, compositional modifications to the TiVNbCr system can now be explored that minimize defect formation and maximize material cycling performance.

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

Machine learning substitutional defect formation energies in ABO 3 perovskites

Perovskite oxides are a promising material platform for use in a wide range of technological applications including electronics, sensors, fuel cells, and catalysis. This is owing to the extraordinary tunability of their physical and chemical properties via defect engineering. The feasibility and the stability of a defect, such as a substitutional dopant, in the host lattice is usually obtained via experiments and/or through detailed quantum mechanical calculations. Both of these conventional routes are expensive and time consuming. An alternative is a data-driven machine learning (ML)-based approach. In this work, we have applied ML techniques to identify the factors that influence defect formation energy, which is an important measure of the stability of the defects, in perovskite oxides. Using 13 elemental properties as features and random forest regression, we demonstrate a systematic approach to down-selecting from the larger set of features to those that are important, establishing a framework for accurate predictions of the defect formation energy. We quantitatively show that the most important factors that control the dopant stability are the dopant ionic size, heat of formation, effective tolerance factor, and oxidation state. Our work reveals previously unknown correlations, chemical trends, and the interplay between stability and underlying chemistries. Hence, these results showcase the efficacy of ML tools in identifying and quantifying different feature-dependencies and provide a promising route toward dopant selection in the perovskites. Herein, we have developed a framework that itself is general and can be potentially applied to other material classes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Pressure Drop Correlation Improvement for the Near-Wall Region of Pebble-Bed Reactors

Packed beds play an important role in several engineering fields, with their applications in nuclear energy being driven by the development of next-generation reactors utilizing pebble fuel. The random nature of a packed pebble bed creates a flow field that is complex and difficult to predict. Porous media models are an attractive option for modeling pebble-bed reactors (PBRs), as they provide intermediate fidelity results and are computationally efficient. Porous media models, however, rely on the use of correlations to estimate the effect of complicated flow features on the pressure drop and heat transfer in the system. Existing correlations were developed to predict the average behavior of the bed, but they are inaccurate in the near-wall region where the presence of the wall affects the pebble packing. This work aims to investigate the accuracy of a porous media model using the Kerntechnischer Ausschuss (KTA) correlation, the most common pressure drop correlation for PBRs compared to the high-fidelity large eddy simulation (LES). A bed of 1568 pebbles is investigated at Reynolds numbers from 625 to 10 000. The bed is divided into five concentric subdomains to compare the average velocity, friction losses, and form losses between the porous media and LES codes. The comparison between the LES simulation and the KTA correlation revealed that the KTA correlation largely underpredicts the form losses in the near-wall region, leading to an overprediction of the velocity near the wall by nearly 30%. An investigation of the form losses across the range of Reynolds numbers in the LES results provided additional insight into how the KTA correlation may be improved to better predict these spatial effects in a pebble bed. These data suggest that the form coefficient near the wall must be increased by 48% while decreasing the form coefficient of the inner bulk region of the bed by 15%. The implementation of these improvements to the KTA correlation in a porous media model produced a radial velocity profile that saw significantly improved agreement with the LES results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

The PSZ-MCMF catalogue of Planck clusters over the DES region

ABSTRACT We present the first systematic follow-up of Planck Sunyaev–Zeldovich effect (SZE) selected candidates down to signal-to-noise (S/N) of 3 over the 5000 deg2 covered by the Dark Energy Survey. Using the MCMF cluster confirmation algorithm, we identify optical counterparts, determine photometric redshifts, and richnesses and assign a parameter, fcont, that reflects the probability that each SZE-optical pairing represents a random superposition of physically unassociated systems rather than a real cluster. The new PSZ-MCMF cluster catalogue consists of 853 MCMF confirmed clusters and has a purity of 90 per cent. We present the properties of subsamples of the PSZ-MCMF catalogue that have purities ranging from 90 per cent to 97.5 per cent, depending on the adopted fcont threshold. Halo mass estimates M500, redshifts, richnesses, and optical centres are presented for all PSZ-MCMF clusters. The PSZ-MCMF catalogue adds 589 previously unknown Planck identified clusters over the DES footprint and provides redshifts for an additional 50 previously published Planck-selected clusters with S/N>4.5. Using the subsample with spectroscopic redshifts, we demonstrate excellent cluster photo-z performance with an RMS scatter in Δz/(1 + z) of 0.47 per cent. Our MCMF based analysis allows us to infer the contamination fraction of the initial S/N>3 Planck-selected candidate list, which is ∼50 per cent. We present a method of estimating the completeness of the PSZ-MCMF cluster sample. In comparison to the previously published Planck cluster catalogues, this new S/N>3 MCMF confirmed cluster catalogue populates the lower mass regime at all redshifts and includes clusters up to z∼1.3.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Density functionals combined with van der Waals corrections for graphene adsorbed on layered materials

Standard density functionals like Perdew-Burke-Ernzerhof (PBE) or Strongly Constrained and Appropriately Normed (SCAN) need a correction to account for long-range van der Waals (vdW) interaction. The damped Zaremba-Kohn model (dZK) starts from a formula for the vdW interaction of a distant atom with a solid surface, both with known dielectric properties, damps this formula at short range, and then treats an adsorbed molecule or atomic layer as a collection of renormalized atoms. Here, we extend this model to graphene adsorbed on semiconducting layered materials [bulk graphite and hBN, and multilayer transition metal dichalcogenides (TMD)] by including the C 4 asymptotic term and multiple electrostatic image effects due to the two surfaces of the substrate slabs in the vdW calculations. The resulting SCAN-vdW-dZK and PBE-vdW-dZK give approximately the same results for the systems considered here, in agreement with available reference values. The predicted binding energies are roughly 25% lower than those from SCAN+rVV10 (revised Vydrov and Van Voorhis 2010), and ~15% lower than those from PBE+rVV10. Since SCAN+rVV10 usually overbinds, the predicted binding energies by SCAN-vdW-dZK and PBE-vdW-dZK are expected to be closer to the true values. The predicted equilibrium binding distances from SCAN-vdW-dZK and PBE-vdW-dZK are slightly larger (~ 0.1 Å) than those from SCAN+rVV10, and close to those from SCAN. The binding energy depends upon the number of substrate layers more strongly in vdW-dZK than in rVV10. The C 4 -term contributions can be 40% of the total vdW interactions, and the C 5 term contributes about 10%. The effects of images and the back surfaces of slabs can contribute about 4–10%. The vdW interaction energy power laws from the vdW-dZK model for graphene adsorbed on multilayer MoS 2 show slowly varying decay to the pairwise exponent –4 with increasing separation $\textit{D}$, very similar to those obtained from the random-phase approximation and renormalization group approaches by Ambrosetti et al. Both the PBE-vdW-dZK and SCAN-vdW-dZK give a greater increase in interlayer binding energy when the TMD substrate changes from monolayer to four-layer for the graphene/TMD adsorption systems. This is consistent with the relevant electron-energy-loss spectroscopy experimental results showing increased dielectric response from the substrate with increasing substrate layer number.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Domain wall motion across microstructural features in polycrystalline ferroelectric films

This paper describes the effect of microstructural features such as grain boundaries and triple points on the pinning of domain wall motion in perovskite Pb(Zr 0.52 Ti 0.48 )O 3 (PZT) films on polycrystalline SrRuO 3 /SrTiO 3 substrates. Here, spatial variability in the collective domain wall dynamics was assessed using non-linearity mapping via Band Excitation Piezoresponse Force Microscopy (BE-PFM). Collocating the non-linearity maps with triple point locations (as visualized by EBSD) allowed for exploration of the effects that local microstructure (e.g., grain boundary) have on domain wall motion. It was found that the extrinsic behavior varied with both the misorientation angle and the proximity to the grain boundary. The width of influence of individual grain boundaries on the motion of domain walls was a function of the character of the grain boundary; random grain boundaries exhibit deeper minima in α d /d 33,initial and larger widths of influence (up to 905 nm) compared to coincident site lattice (CSL) boundaries (up to 572 nm). Additionally, triple points containing larger numbers of random boundaries exhibited non-Rayleigh behavior to greater distances, suggesting that the triple point provides either a deep potential minimum or a region where domain wall motion is unfavorable.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling and simulation of vascular tumors embedded in evolving capillary networks

Here in this work, we present a coupled 3D–1D model of solid tumor growth within a dynamically changing vascular network to facilitate realistic simulations of angiogenesis. Additionally, the model includes erosion of the extracellular matrix, interstitial flow, and coupled flow in blood vessels and tissue. We employ continuum mixture theory with stochastic Cahn–Hilliard type phase-field models of tumor growth. The interstitial flow is governed by a mesoscale version of Darcy’s law. The flow in the blood vessels is controlled by Poiseuille flow, and Starling’s law is applied to model the mass transfer in and out of blood vessels. The evolution of the network of blood vessels is orchestrated by the concentration of the tumor angiogenesis factors (TAFs); blood vessels grow towards the increasing TAFs concentrations. This process is not deterministic, allowing random growth of blood vessels and, therefore, due to the coupling of nutrients in tissue and vessels, makes the growth of tumors stochastic. We demonstrate the performance of the model by applying it to a variety of scenarios. Numerical experiments illustrate the flexibility of the model and its ability to generate satellite tumors. Simulations of the effects of angiogenesis on tumor growth are presented as well as sample-independent features of cancer.

3D–1D coupled blood flow models↗
DOE OSTI2023

Capturing the softening in T91 steel annealed with molten LBE through interfacial gradient plasticity

In order to study the effect of annealing with molten lead–bismuth eutectic alloy (mLBE) on the mechanical properties of T91 steel, micropillar compression tests were analyzed through the consideration of mechanical interface energy terms within gradient plasticity. The stress–strain curves of the micropillars showed significant stochastic effects, characterized by a different elastic and plastic behavior response for each pillar. Among the various pillars, some exhibited a similar trilinear behavior. Scanning electron microscopy (SEM) images of these specimens demonstrated localized slip deformation and slip planes after compression, while transmission electron microscopy (TEM) images showed the presence of a severe slip zone along the random grain boundaries (RGBs), indicating that the GBs played a dominant role in the deformation/slip. By employing interfacial gradient plasticity that can explicitly account for the presence of GBs, the trilinear response was captured, allowing for the determination of the mechanical interface parameter. As anticipated the pillars which underwent severe slip at the GBs had a lower value for the mechanical interface parameter. Finally, annealing with mLBE can, therefore, in some cases result in softening in the overall stress–strain, as it lowers the mechanical interface energy of GBs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamic probabilistic risk assessment and game theory for cyber security risk analysis in nuclear power plants

Nuclear Power Plants and energy systems have become more prone to cyber-attacks with their digitalization and the increased use of smart equipment. Hence, it is important to quantify the risk associated with cyber-attacks in such systems. Dynamic Probabilistic Risk Assessment which involves studying the evolution of a system due to random events and operator and attacker actions during a cyber-attack by employing a physics-based model of the system is a suitable framework to quantify cybersecurity risk in nuclear power plants. In addition to the plant dynamics, it is also important to model the strategies of the attackers and plant operators for an effective cybersecurity risk assessment. Game theory provides a set of necessary tools to model such strategic interactions. In this research, a framework that integrates dynamic probabilistic risk assessment with game theory for cybersecurity risk analysis in nuclear power plants is presented. The mathematical formulation is derived based on the theory of continuous event trees. We propose a game theory based action model, that utilizes physics-based rewards to define the strategies of attackers and operators at every decision epoch. As a case study, the risk associated with cyber-attacks on the digital components in the secondary side of a pressurized water reactor is studied using a reduced order model. A set of attacker actions and a set of operator actions are defined for the system. The operator and attacker interactions were modelled using simultaneous game, their action policies were computed using the concept of mixed strategy Nash equilibrium and the evolution of the system was studied.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Entropy-Driven Structural Evolution in Ceramic Oxides

High-entropy ceramics, with five or more elements randomly occupying the same cation crystallographic sites, offer vast compositional diversity and unique properties for material design and applications. However, for many dissimilar elements, entropic stabilization cannot overcome the enthalpic barrier to cation substitution. As a result, most high-entropy ceramics incorporate only a few similar elements, limiting the in-depth exploration of the effect of entropy on ceramic properties. Here, we first use density functional theory to model fluorite crystal structures composed of 1-10 elements and then experimentally present practical fluorite oxide nanostructures containing 1, 3, 8, and 15 metals, as well as a record-breaking 25-element high-entropy ceramic incorporating a diverse palette of rare-earth, transition, alkaline, p-block, and noble metals. As entropy increases, structural and configurational disorder in the solid solution rises, altering structural features such as lattice distortion, crystallinity, homogeneity, defect density, and thermal stability. This research provides new insights and understanding of the role of entropy in stabilizing compositionally complex ceramics.

Liu, Shuo↗
DOE OSTI2023

Influences of cloud vertical overlapping on the calculated cloud albedo and their validation with satellite observations

Cloud albedo is expected to influence cloud radiative forcing in addition to cloud fraction, and inadequate description of the cloud overlapping effects on the cloud fraction may influence the simulated cloud fraction, and thus the relative cloud radiative forcing (RCRF) and cloud albedo. In this study, we first present a new formula by extending that presented previously in Liu et al. (2011) to consider multilayer clouds directly in the relationship between cloud albedo, cloud fraction and RCRF, and then quantitatively evaluate the effects of different cloud vertical overlapping structures, represented by the decorrelation length scales (L cf ), on the simulated cloud albedos. We use the BCC_AGCM2.0_CUACE/Aero model with simultaneous validation by observations from the Clouds and the Earth's Radiation Energy System (CERES) satellite. When L cf < 4 km (i.e., the cloud overlap is closer to the random overlap), the simulated cloud albedos are generally in good agreement with the satellite-based albedos for December– February and June–August; when L cf ≥ 4 km (i.e., the cloud vertical overlap is closer to the maximum overlap), the difference between simulated and observed cloud albedos became larger, due mainly to significant differences in cloud fractions and RCRF. Further quantitative analysis shows that the relative Euclidean distance, which represents the degree of overall model–observation disagreement, increases with the L cf for all three variables (cloud albedo, cloud fraction and RCRF), indicating the importance of cloud vertical overlapping in determining the accuracy of the calculated cloud albedo for multilayer clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sensitivity of fatigue reliability in wind turbines: effects of design turbulence and the Wöhler exponent

Fatigue assessment of wind turbines involves three main sources of uncertainty: material resistance, load, and the damage accumulation model. Many studies focus on increasing the accuracy of fatigue load assessment to improve the fatigue reliability. Probabilistic modeling of the wind's turbulence standard deviation is an example of an approach used for this purpose. Editions 3 and 4 of the IEC standard for the design of wind energy generation systems (IEC 61400-1) suggest different probability distributions as alternatives for the representative turbulence in the normal turbulence model (NTM) of edition 1. There are debates on whether the suggested distributions provide conservative reliability levels, as the established design safety factors are calibrated based on the representative turbulence approach. The current study addresses the debate by comparing annual reliability based on different scenarios of NTM using a probabilistic approach. More importantly, it elaborates on the relative importance of load assessment accuracy in defining the fatigue reliability. Using the DTU 10 MW reference wind turbine and the first-order reliability method (FORM), we study the changes in the annual reliability level and its sensitivity to the three main random inputs. We perform the study considering the blade root flapwise and the tower base fore–aft moments, assuming different fatigue exponents in each load channel. The results show that integration over distributions of turbulence in each mean wind speed results in less conservative annual reliability levels than representative turbulence. The difference in the reliability levels varies according to turbulence distribution and the fatigue exponent. In the case of the tower base, the difference in the annual reliability index after 20 years can be up to 50 %. However, the model and material uncertainty have much higher effects on the reliability levels compared to load uncertainty. Knowledge about such differences in the reliability levels due to the choice of turbulence distribution is especially important, as it impacts the extent of lifetime extension through reliability reassessments.

17 WIND ENERGY↗
DOE OSTI2020

Sub-micrometer random-pyramid texturing of silicon solar wafers with excellent surface passivation and low reflectance

Select emerging solar cell technologies, ranging from very narrow screen-printed fingers to perovskite/silicon tandems, would benefit from a reduction of the height of the standard pyramidal texture of monocrystalline silicon wafers to below 1 μm. However, fabricating such small surface features usually compromises the requisite low reflectance and high passivation quality. Through systematic design of experiments, we demonstrate that conventional alkaline-based wet-chemical etching processes can be tuned to produce dense, spatially homogeneous, and uniformly sized sub-micrometer pyramids. Precise texture size can be further controlled by adding potassium silicate to the texturing solution to effectively slow the pyramid growth rate and promote nucleation. With pyramid size distribution data—extracted from multiple atomic force microscopy images across each wafer—we demonstrate 16 nanotextures with average pyramid sizes ranging from 62 to 512 nm, with the largest having only 0.3% of its pyramids taller than 1 μm. Six of the nanotextures exhibit light trapping and surface defect densities comparable to state-of-the-art microtextured wafers. In conclusion, the best nanotexture achieved an AM1.5G-weighted reflectance of 11.8% and a minority-carrier lifetime of 3.4 ms after passivation with amorphous silicon, demonstrating that smaller textures need not deteriorate solar cell performance.

14 SOLAR ENERGY↗
DOE OSTI2022

The nanocaterpillar's random walk: diffusion with ligand–receptor contacts

Particles with ligand–receptor contacts bind and unbind fluctuating “legs” to surfaces, whose fluctuations cause the particle to diffuse. Quantifying the diffusion of such “nanoscale caterpillars” is a challenge, since binding events often occur on very short time and length scales. Here, in this work, we derive an analytical formula, validated by simulations, for the long time translational diffusion coefficient of an overdamped nanocaterpillar, under a range of modeling assumptions. We demonstrate that the effective diffusion coefficient, which depends on the microscopic parameters governing the legs, can be orders of magnitude smaller than the background diffusion coefficient. Furthermore it varies rapidly with temperature, and reproduces the striking variations seen in existing data and our own measurements of the diffusion of DNA-coated colloids. Our model gives insight into the mechanism of motion, and allows us to ask: when does a nanocaterpillar prefer to move by sliding, where one leg is always linked to the surface, and when does it prefer to move by hopping, which requires all legs to unbind simultaneously? We compare a range of systems (viruses, molecular motors, white blood cells, protein cargos in the nuclear pore complex, bacteria such as Escherichia coli, and DNA-coated colloids) and present guidelines to control the mode of motion for materials design.

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

Investigating Particle Size‐Dependent Redox Kinetics and Charge Distribution in Disordered Rocksalt Cathodes

Abstract Understanding how various redox activities evolve and distribute in disordered rocksalt oxides (DRX) can advance insights into manipulating materials properties for achieving stable, high‐energy batteries. Herein, the authors present how the reaction kinetics and spatial distribution of redox activities are governed by the particle size of DRX materials. The size‐dependent electrochemical performance is attributed to the distinct cationic and anionic reaction kinetics at different sizes, which can be tailored to achieve optimal capacity and stability. Overall, the local charged domains in DRX particles display random heterogeneity caused by the isotropic delithiation pathways. Owing to the kinetic limitation, the micron‐sized particles exhibit a holistic “core‐shell” charge distribution, whereas sub‐micron particles show more uniform redox reactions throughout the particles and ensembles. Sub‐micron DRX particles exhibit increasing anionic redox activities yet inferior cycling stability. In summary, engineering particle size can effectively modulate how cationic and anionic redox activities evolve and distribute in DRX materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Constraining primordial non-Gaussianity from the large scale structure two-point and three-point correlation functions

Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by f NL ⁠, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spectrum. We propose a framework to measure f NL using the configuration space two-point correlation function (2pcf) monopole and three-point correlation function (3pcf) monopole of survey tracers. Our model estimates the effect of the scale-dependent bias induced by the presence of PNG on the 2pcf and 3pcf from the clustering of simulated dark matter haloes. We describe how this effect may be scaled to an arbitrary tracer of the cosmological matter density. The 2pcf and 3pcf of this tracer are measured to constrain the value of f NL ⁠. In LSS surveys, the effect of imaging systematics on two-point statistics is often degenerate with the PNG signal. Our proposed model employs three-point statistics primarily to break this degeneracy. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our method. Our forecast indicates the ability to constrain f NL to a precision of σf NL ≈ 22 with one year of DESI survey data, as well as the ability to constrain the imaging systematic weights in situ.

early Universe↗