Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear inverse modeling”

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

Integrated modeling of RF-induced tungsten erosion at ICRH antenna structures in the WEST tokamak *

This paper introduces STRIPE (Simulated Transport of RF Impurity Production and Emission), an advanced modeling framework developed to analyze material erosion and the global transport of eroded impurities originating from radio-frequency (RF) antenna structures in magnetic confinement fusion devices. STRIPE integrates multiple physics modules: SolEdge3x for scrape-off-layer plasma profiles, COMSOL for 3D RF rectified sheath potentials, RustBCA for erosion yields and surface interactions, and global impurity transport for 3D ion energy-angle distributions and impurity transport. The framework is applied to an ion cyclotron RF-heated L-mode discharge (#57877) in the WEST tokamak, where it predicts a thirty-fold increase in gross tungsten erosion at antenna limiters during the transition from ohmic to ICRH operation. Additionally, under ICRH conditions, a tenfold enhancement in erosion is observed when comparing RF sheath effects to purely thermal sheath conditions. High-charge-state oxygen ions ($\mathrm{O}$ 6+ and above) are identified as the dominant contributors to tungsten sputtering. To validate the model, a synthetic diagnostic tool based on inverse photon efficiency (S/XB coefficients) from the ColRadPy collisional-radiative model enables direct comparison with spectroscopic measurements. Model predictions using a plasma composition of 1% oxygen and 99% deuterium show good agreement with observed W − I (400.9 nm) emission for discharge #57877, supporting the accuracy of the STRIPE framework. This study focuses specifically on gross erosion calculations to demonstrate STRIPE’s capabilities. Future extensions of this work will incorporate net erosion, re-deposition, self-sputtering effects, and whole-device modeling of sputtered tungsten impurity transport. STRIPE is also being applied to other RF-heated linear and toroidal devices, offering valuable insights for antenna design, impurity control, and performance optimization in next-generation fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Three Hundred Project: Modeling baryon and hot-gas fraction evolution in simulated clusters

The baryon fraction of galaxy clusters, expressed as the ratio between the mass in baryons (including both stars and cold or hot gas) and the total mass, is a powerful tool to provide information on the cosmological parameters, while the hot-gas fraction provides indications on the physics of the intracluster plasma and its interplay with the processes that drive galaxy formation. Using cosmological hydrodynamical simulations of about 300 simulated massive galaxy clusters with a median mass M 500 ≈ 7 × 10 14 M ⊙ at z = 0, we model the relations between total mass and either baryon fraction or the hot gas fractions at overdensities Δ = 2500, 500, and 200 with respect to the cosmic critical density, and their evolution from z ∼ 0 to z ∼ 1.3. We utilized the simulated galaxy clusters from the Three Hundred project, which include star formation and feedback from both supernovae and active galactic nuclei. We fit the simulation results for such scaling relations against three analytic forms (linear, quadratic, and logarithmic in a logarithmic plane) and three forms for the redshift dependence, and we considered as a variable both the inverse of the cosmic scale factor, (1 + z), and the Hubble expansion rate, E(z). We show that power-law dependencies on cluster mass poorly describe the investigated relations. A power law fails to simultaneously capture the flattening of the total baryon and gas fractions at high masses, their drop at low masses, and the transition between these two regimes. The other two functional forms provide a more accurate description of the curvature in mass scaling. The fractions measured within smaller radii exhibit a stronger evolution than those measured within larger radii. From the analysis of these simulations, we evince that as long as we include systems in the mass range herein investigated, the baryon or gas fraction can be accurately related to the total mass through either a parabola or a logarithm in the logarithmic plane. The trends are common to all modern hydro simulations, although the amplitude of the drop at low masses might differ. Being able to observationally determine the gas fraction in groups will thus provide constraints on the baryonic physics.

galaxy clusters↗
NASA NTRS2012

Application of Artificial Neural Networks to the Development of Improved Multi-Sensor Retrievals of Near-Surface Air Temperature and Humidity Over Ocean

Improved estimates of near-surface air temperature and air humidity are critical to the development of more accurate turbulent surface heat fluxes over the ocean. Recent progress in retrieving these parameters has been made through the application of artificial neural networks (ANN) and the use of multi-sensor passive microwave observations. Details are provided on the development of an improved retrieval algorithm that applies the nonlinear statistical ANN methodology to a set of observations from the Advanced Microwave Scanning Radiometer (AMSR-E) and the Advanced Microwave Sounding Unit (AMSU-A) that are currently available from the NASA AQUA satellite platform. Statistical inversion techniques require an adequate training dataset to properly capture embedded physical relationships. The development of multiple training datasets containing only in-situ observations, only synthetic observations produced using the Community Radiative Transfer Model (CRTM), or a mixture of each is discussed. An intercomparison of results using each training dataset is provided to highlight the relative advantages and disadvantages of each methodology. Particular emphasis will be placed on the development of retrievals in cloudy versus clear-sky conditions. Near-surface air temperature and humidity retrievals using the multi-sensor ANN algorithms are compared to previous linear and non-linear retrieval schemes.

Roberts, J. Brent↗
NASA NTRS1978

Evidence for a non-random magnetization of the moon

An approximate inversion, using a specific model, of Apollo subsatellite vector magnetometer data has provided estimates for the locations, scale sizes, and magnetization characteristics of strongly magnetized lunar crustal materials across a portion of the lunar far side. The nearly model-independent local bulk directions of magnetization, while otherwise random, exhibit a tendency to lie in or near planes perpendicular to the present lunar spin axis. The apparent existence of such a preferred plane of magnetization must be regarded as unfavorable for proposed local, presumably random sources of large-scale, coherent magnetization and implies that any assumed internal magnetizing dipole or external linear magnetizing field must also have been oriented mainly perpendicular to the present spin axis. The otherwise random orientations of these vectors when projected onto the equatorial plane suggests a simple explanation for the negligible lunar global dipole moment.

Hood, L.↗
DOE OSTI2021

The anion exchange properties of AFm (hydrocalumite-group) minerals defined from solubility experiments and crystallographic information

Solubility experiments with binary series of OH-, Cl-, SO{sub 4}- and CO{sub 3}-AFm were re-examined with a non-ideal, multicomponent solid solution model. The electrical double layer on the mineral's outer surface, which neutralizes about 5% of the anionic charge, can be included. The choice of the end-members in the solid solution and the reaction formulation can help to eliminate the ‘non-ideality’ for OH-SO{sub 4} and Cl-OH exchange. CO{sub 3}-Cl-OH exchange is not ideal in the solubility experiments, nor can it be, because the mineral volume is not a linear combination of the end-members. The solubilities of the end-members are related to the anion charge and inversely to the anion-layer thickness. OH is part of the anion-charge, even when a pure, other anion form is the target of the synthesis. Pure OH-AFm is unstable and converts into hydrogarnet and portlandite, depending on the precipitation rate of hydrogarnet.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
DOE OSTI2022

DIAS: A Data-Informed Active Subspace Regularization Framework for Inverse Problems

This paper presents a regularization framework that aims to improve the fidelity of Tikhonov inverse solutions. At the heart of the framework is the data-informed regularization idea that only data-uninformed parameters need to be regularized, while the data-informed parameters, on which data and forward model are integrated, should remain untouched. We propose to employ the active subspace method to determine the data-informativeness of a parameter. The resulting framework is thus called a data-informed (DI) active subspace (DIAS) regularization. Four proposed DIAS variants are rigorously analyzed, shown to be robust with the regularization parameter and capable of avoiding polluting solution features informed by the data. They are thus well suited for problems with small or reasonably small noise corruptions in the data. Furthermore, the DIAS approaches can effectively reuse any Tikhonov regularization codes/libraries. Though they are readily applicable for nonlinear inverse problems, we focus on linear problems in this paper in order to gain insights into the framework. Various numerical results for linear inverse problems are presented to verify theoretical findings and to demonstrate advantages of the DIAS framework over the Tikhonov, truncated SVD, and the TSVD-based DI approaches.

Nguyen, Hai↗
DOE OSTI2022

New insights into Mn 2+ and Mg 2+ inhibition of calcite growth

Impurity ion and isotope partitioning into carbonate minerals provide a window into the molecular processes occurring at the fluid-mineral interface during crystal growth. Here, we employ calcium isotope fractionation together with process-based modeling to elucidate the mechanisms by which two divalent cations with starkly contrasting compatibility, magnesium and manganese, inhibit calcite growth and incorporate into the mineral lattice. Calcite growth inhibition by Mg 2+ is log-linear and K Mg is on the order of 0.02–0.03 throughout the range of {Mg 2+ }/{Ca 2+ } studied here (0.01–2.6). Mn 2+ exhibits much stronger log-linear growth rate inhibition at low Mn 2+ concentrations (fluid {Mn 2+ }/{Ca 2+ } = 0.001–0.02). Mn 2+ is readily incorporated into the calcite lattice to form a calcite-rhodochrosite solid solution, with large partition coefficients (K Mn 4.6–15.6) inversely correlated to growth rate. For both Mn 2+ and Mg 2+ , calcium isotope fractionation is found to be invariant with {Me 2+ }/{Ca 2+ } despite more than an order of magnitude decline in growth rate. This invariant Δ 44/40 Ca suggests that the presence of Mn 2 + or Mg 2+ does not significantly change the relative rates of Ca 2+ attachment and detachment at kink sites during growth, indicative of a dominantly kink blocking inhibition mechanism. Because the partitioning behavior dictates that Mn 2+ must attach to the surface significantly faster than Ca 2+ , attachment of Mn 2+ is likely to be as a non-monomer species such as an ion pair or possibly a larger polynuclear cluster. We propose that calcite growth rate inhibition by Mn is determined by the kinetics of carbonate attachment at Mn-occupied kink sites, potentially due to slow re-orientation kinetics of carbonate ions that have formed an inner-sphere complex with Mn 2+ at the surface but must reorient to incorporate into the lattice. We demonstrate that patterns in Mg 2+ partitioning and inhibition behavior are broadly consistent with growth inhibition driven by slow Mg 2+ -aquo complex dehydration relative to Ca 2+ but argue that this mechanism likely represents one endmember scenario, seen in Mg-calcite growth at low supersaturations and net precipitation rates. During growth at faster net precipitation rates, some portion of Mg 2+ is likely incorporated as a partially hydrated or otherwise complexed species, but calcite growth remains significantly inhibited by the kinetics of CO 3 2- attachment at Mg 2+ kink sites. These findings suggest a hybrid classical/nonclassical growth mechanism whereby Ca 2+ incorporates largely as a free ion at kink sites while Mn 2+ and some portion of Mg 2+ are incorporated via non-monomer attachment. This pattern may be generalizable; trace constituent cations with aquo-complex desolvation rates significantly slower than the mineral growth rate preferentially incorporate as a non-monomer species during otherwise classical crystal growth.

58 GEOSCIENCES↗
DOE OSTI2022

Interpreting radial correlation Doppler reflectometry using gyrokinetic simulations

A linear response, local model for the DBS amplitude applied to gyrokinetic simulations shows that radial correlation Doppler reflectometry measurements (RCDR, Schirmer et al 2007 Plasma Phys. Control. Fusion 49 1019) are not sensitive to the average turbulence radial correlation length, but to a correlation length that depends on the binormal wavenumber $k_{\perp}$ selected by the Doppler backscattering (DBS) signal. Nonlinear gyrokinetic simulations show that the turbulence naturally exhibits a nonseparable power law spectrum in wavenumber space, leading to a power law dependence of the radial correlation length with binormal wavenumber $l_r \sim C k_{\perp}^{-\alpha} (\alpha \approx 1)$ which agrees with the inverse proportionality relationship between the measured $l_r$ and $k_{\perp} $ observed in experiments (Fernández-Marina et al 2014 Nucl. Fusion 54 072001). This new insight indicates that RCDR characterizes the eddy aspect ratio in the perpendicular plane to the magnetic field. It also motivates future use of a nonseparable turbulent spectrum to quantitatively interpret RCDR and potentially other turbulence diagnostics. The radial correlation length is only measurable when the radial resolution at the cutoff location $W_\text n$ satisfies $W_\text n \ll l_r$, while the measurement becomes dominated by $W_\text n$ for $W_\text n \gg l_r$. This suggests that $l_r$ is likely to be inaccessible for electron-scale DBS measurements (${k_{\perp}}{\rho_{\text s}} \gt 1$). The effect of $W_\text n$ on ion-scale radial correlation lengths could be nonnegligible.

radial correlation Doppler reflectometry↗
DOE OSTI2022

Skyrmions and hopfions in three-dimensional frustrated magnets

Here, a model of an inversion-symmetric frustrated spin system is introduced which hosts three-dimensional extensions of magnetic skyrmions. In the continuum approximation, this model reduces to a nonlinear sigma model on a squashed sphere that has a natural interpolating parameter. At one limit of the parameter, the model reduces to a frustrated magnetic system earlier considered by Sutcliffe as a host to hopfions, and in the other limit, it becomes very similar to the 3D Skyrme model. To better understand the relation between hopfions and 3D skyrmions, a model interpolating between the Faddeev-Niemi and the Skyrme models is reconsidered and it is shown that energies of the solitons obey a linear Bogomol'nyi-Prasad-Sommerfeld bound. The 3D skyrmions in the frustrated magnetic model are found and compared to the rational map ansatz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2005

Adjoint Sensitivity Analysis of Orbital Mechanics: Application to Computations of Observables' Partials with Respect to Harmonics of the Planetary Gravity Fields

An approach is presented to the inversion of gravity fields based on evaluation of partials of observables with respect to gravity harmonics using the solution of adjoint problem of orbital dynamics of the spacecraft. Corresponding adjoint operator is derived directly from the linear operator of the linearized forward problem of orbital dynamics. The resulting adjoint problem is similar to the forward problem and can be solved by the same methods. For given highest degree N of gravity harmonics desired, this method involves integration of N adjoint solutions as compared to integration of N2 partials of the forward solution with respect to gravity harmonics in the conventional approach. Thus, for higher resolution gravity models, this approach becomes increasingly more effective in terms of computer resources as compared to the approach based on the solution of the forward problem of orbital dynamics.

gravity fields↗
DOE OSTI2023

Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys

Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.

density functional theory↗
NASA NTRS1972

The general linear inverse problem - Implication of surface waves and free oscillations for earth structure.

The discrete general linear inverse problem reduces to a set of m equations in n unknowns. There is generally no unique solution, but we can find k linear combinations of parameters for which restraints are determined. The parameter combinations are given by the eigenvectors of the coefficient matrix. The number k is determined by the ratio of the standard deviations of the observations to the allowable standard deviations in the resulting solution. Various linear combinations of the eigenvectors can be used to determine parameter resolution and information distribution among the observations. Thus we can determine where information comes from among the observations and exactly how it constraints the set of possible models. The application of such analyses to surface-wave and free-oscillation observations indicates that (1) phase, group, and amplitude observations for any particular mode provide basically the same type of information about the model; (2) observations of overtones can enhance the resolution considerably; and (3) the degree of resolution has generally been overestimated for many model determinations made from surface waves.

Wiggins, R. A.↗
NASA NTRS1996

Possible effects of protracted exposure on the additivity of risks from space radiations

Conventional radiation risk assessments are presently based on the additivity assumption. This assumption states that risks from individual components of a complex radiation field involving many different types of radiation can be added to yield the total risk of the complex radiation field. If the assumption is not correct, the summations and integrations performed to obtain the presently quoted risk estimates are not appropriate. This problem is particularly important in the area of space radiation risk evaluation because of the many different types of high- and low-LET radiation present in the galactic cosmic ray environment. For both low- and high-LET radiations at low enough dose rates, the present convention is that the addivity assumption holds. Mathematically, the total risk, Rtot is assumed to be Rtot = summation (i) Ri where the summation runs over the different types of radiation present. If the total dose (or fluence) from each component is such that the interaction between biological lesions caused by separate single track traversals is negligible within a given cell, it is presently considered to be reasonable to accept the additivity assumption. However, when the exposure is protracted over many cell doubling times (as will be the case for extended missions to the moon or Mars), the possibility exists that radiation effects that depend on multiple cellular events over a long time period, such as is probably the case in radiation-induced carcinogenesis, may not be additive in the above sense and the exposure interval may have to be included in the evaluation procedure. It is shown, however, that "inverse" dose-rate effects are not expected from intermediate LET radiations arising from the galactic cosmic ray environment due to the "sensitive-window-in-the-cell-cycle" hypothesis.

NASA Discipline Number 45-10↗
NASA NTRS2010

Tool for Forecasting Cool-Season Peak Winds Across Kennedy Space Center and Cape Canaveral Air Force Station

The expected peak wind speed for the day is an important element in the daily morning forecast for ground and space launch operations at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45th Weather Squadron (45 WS) must issue forecast advisories for KSC/CCAFS when they expect peak gusts for >= 25, >= 35, and >= 50 kt thresholds at any level from the surface to 300 ft. In Phase I of this task, the 45 WS tasked the Applied Meteorology Unit (AMU) to develop a cool-season (October - April) tool to help forecast the non-convective peak wind from the surface to 300 ft at KSC/CCAFS. During the warm season, these wind speeds are rarely exceeded except during convective winds or under the influence of tropical cyclones, for which other techniques are already in use. The tool used single and multiple linear regression equations to predict the peak wind from the morning sounding. The forecaster manually entered several observed sounding parameters into a Microsoft Excel graphical user interface (GUI), and then the tool displayed the forecast peak wind speed, average wind speed at the time of the peak wind, the timing of the peak wind and the probability the peak wind will meet or exceed 35, 50 and 60 kt. The 45 WS customers later dropped the requirement for >= 60 kt wind warnings. During Phase II of this task, the AMU expanded the period of record (POR) by six years to increase the number of observations used to create the forecast equations. A large number of possible predictors were evaluated from archived soundings, including inversion depth and strength, low-level wind shear, mixing height, temperature lapse rate and winds from the surface to 3000 ft. Each day in the POR was stratified in a number of ways, such as by low-level wind direction, synoptic weather pattern, precipitation and Bulk Richardson number. The most accurate Phase II equations were then selected for an independent verification. The Phase I and II forecast methods were compared using an independent verification data set. The two methods were compared to climatology, wind warnings and advisories issued by the 45 WS, and North American Mesoscale (NAM) model (MesoNAM) forecast winds. The performance of the Phase I and II methods were similar with respect to mean absolute error. Since the Phase I data were not stratified by precipitation, this method's peak wind forecasts had a large negative bias on days with precipitation and a small positive bias on days with no precipitation. Overall, the climatology methods performed the worst while the MesoNAM performed the best. Since the MesoNAM winds were the most accurate in the comparison, the final version of the tool was based on the MesoNAM winds. The probability the peak wind will meet or exceed the warning thresholds were based on the one standard deviation error bars from the linear regression. For example, the linear regression might forecast the most likely peak speed to be 35 kt and the error bars used to calculate that the probability of >= 25 kt = 76%, the probability of >= 35 kt = 50%, and the probability of >= 50 kt = 19%. The authors have not seen this application of linear regression error bars in any other meteorological applications. Although probability forecast tools should usually be developed with logistic regression, this technique could be easily generalized to any linear regression forecast tool to estimate the probability of exceeding any desired threshold . This could be useful for previously developed linear regression forecast tools or new forecast applications where statistical analysis software to perform logistic regression is not available. The tool was delivered in two formats - a Microsoft Excel GUI and a Tool Command Language/Tool Kit (Tcl/Tk) GUI in the Meteorological Interactive Data Display System (MIDDS). The Microsoft Excel GUI reads a MesoNAM text file containing hourly forecasts from 0 to 84 hours, from one model run (00 or 12 UTC). The GUI then displays e peak wind speed, average wind speed, and the probability the peak wind will meet or exceed the 25-, 35- and 50-kt thresholds. The user can display the Day-1 through Day-3 peak wind forecasts, and separate forecasts are made for precipitation and non-precipitation days. The MIDDS GUI uses data from the NAM and Global Forecast System (GFS), instead of the MesoNAM. It can display Day-1 and Day-2 forecasts using NAM data, and Day-1 through Day-5 forecasts using GFS data. The timing of the peak wind is not displayed, since the independent verification showed that none of the forecast methods performed significantly better than climatology. The forecaster should use the climatological timing of the peak wind (2248 UTC) as a first guess and then adjust it based on the movement of weather features.

Barrett, Joe H., III↗
DOE OSTI2022

LaSDI: Parametric Latent Space Dynamics Identification

Enabling fast and accurate physical simulations with data has become an important area of computational physics to aid in inverse problems, design-optimization, uncertainty quantification, and other various decision-making applications. This paper presents a data-driven framework for parametric latent space dynamics identification procedure that enables fast and accurate simulations. The parametric model is achieved by building a set of local latent space model and designing an interaction among them. An individual local latent space dynamics model achieves accurate solution in a trust region. By letting the set of trust region to cover the whole parameter space, our model shows an increase in accuracy with an increase in training data. Herein we introduce two different types of interaction mechanisms, i.e., point-wise and region-based approach. Both linear and nonlinear data compression techniques are used. We illustrate the framework of Latent Space Dynamics Identification (LaSDI) enable a fast and accurate solution process on various partial differential equations, i.e., Burgers’ equations, radial advection problem, and nonlinear heat conduction problem, achieving 0 (100)x speed-up and 0 (1)% relative error with respect to the corresponding full order models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Distribution of centrality measures on undirected random networks via the cavity method

The Katz centrality of a node in a complex network is a measure of the node’s importance as far as the flow of information across the network is concerned. For ensembles of locally tree-like undirected random graphs, this observable is a random variable. Its full probability distribution is of interest but difficult to handle analytically because of its “global” character and its definition in terms of a matrix inverse. Leveraging a fast Gaussian Belief Propagation-Cavity algorithm to solve linear systems on tree-like structures, we show that i) the Katz centrality of a single instance can be computed recursively in a very fast way, and ii) the probability P ( K ) that a random node in the ensemble of undirected random graphs has centrality K satisfies a set of recursive distributional equations, which can be analytically characterized and efficiently solved using a population dynamics algorithm. We test our solution on ensembles of Erdős-Rényi and Scale Free networks in the locally tree-like regime, with excellent agreement. The analytical distribution of centrality for the configuration model conditioned on the degree of each node can be employed as a benchmark to identify nodes of empirical networks with over- and underexpressed centrality relative to a null baseline. We also provide an approximate formula based on a rank- 1 projection that works well if the network is not too sparse, and we argue that an extension of our method could be efficiently extended to tackle analytical distributions of other centrality measures such as PageRank for directed networks in a transparent and user-friendly way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗