Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameters”

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

Monte Carlo method for constructing confidence intervals with unconstrained and constrained nuisance parameters in the NOvA experiment

Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. Wilks' theorem provides a simple way to construct confidence intervals on model parameters, but it only applies under certain conditions. These conditions, such as nested hypotheses and unbounded parameters, are often violated in neutrino oscillation measurements and other experimental scenarios. Monte Carlo methods can address these issues, albeit at increased computational cost. In the presence of nuisance parameters, however, the best way to implement a Monte Carlo method is ambiguous. Furthermore, this paper documents the method selected by the NOvA experiment, the profile construction. It presents the toy studies that informed the choice of method, details of its implementation, and tests performed to validate it. It also includes some practical considerations which may be of use to others choosing to use the profile construction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Finding simplicity: unsupervised discovery of features, patterns, and order parameters via shift-invariant variational autoencoders *

Abstract Recent advances in scanning tunneling and transmission electron microscopies (STM and STEM) have allowed routine generation of large volumes of imaging data containing information on the structure and functionality of materials. The experimental data sets contain signatures of long-range phenomena such as physical order parameter fields, polarization, and strain gradients in STEM, or standing electronic waves and carrier-mediated exchange interactions in STM, all superimposed onto scanning system distortions and gradual changes of contrast due to drift and/or mis-tilt effects. Correspondingly, while the human eye can readily identify certain patterns in the images such as lattice periodicities, repeating structural elements, or microstructures, their automatic extraction and classification are highly non-trivial and universal pathways to accomplish such analyses are absent. We pose that the most distinctive elements of the patterns observed in STM and (S)TEM images are similarity and (almost-) periodicity, behaviors stemming directly from the parsimony of elementary atomic structures, superimposed on the gradual changes reflective of order parameter distributions. However, the discovery of these elements via global Fourier methods is non-trivial due to variability and lack of ideal discrete translation symmetry. To address this problem, we explore the shift-invariant variational autoencoders (shift-VAEs) that allow disentangling characteristic repeating features in the images, their variations, and shifts that inevitably occur when randomly sampling the image space. Shift-VAEs balance the uncertainty in the position of the object of interest with the uncertainty in shape reconstruction. This approach is illustrated for model 1D data, and further extended to synthetic and experimental STM and STEM 2D data. We further introduce an approach for training shift-VAEs that allows finding the latent variables that comport to known physical behavior. In this specific case, the condition is that the latent variable maps should be smooth on the length scale of the atomic lattice (as expected for physical order parameters), but other conditions can be imposed. The opportunities and limitations of the shift VAE analysis for pattern discovery are elucidated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Qualitative and quantitative enhancement of parameter estimation for model-based diagnostics using automatic differentiation with an application to inertial fusion

Parameter estimation using observables is a fundamental concept in the experimental sciences. Mathematical models that represent the physical processes can enable reconstructions of the experimental observables and greatly assist in parameter estimation by turning it into an optimization problem which can be solved by gradient-free or gradient-based methods. In this work, the recent rise in flexible frameworks for developing differentiable scientific computing programs is leveraged in order to dramatically accelerate data analysis of a common experimental diagnostic relevant to laser–plasma and inertial fusion experiments, Thomson scattering. A differentiable Thomson-scattering data analysis tool is developed that uses reverse-mode automatic differentiation (AD) to calculate gradients. By switching from finite differencing to reverse-mode AD, three distinct outcomes are achieved. First, gradient descent is accelerated dramatically to the extent that it enables near real-time usage in laser–plasma experiments. Second, qualitatively novel quantities which require $\mathcal{O}(10^3)$ parameters can now be included in the analysis of data which enables unprecedented measurements of small-scale laser–plasma phenomena. Third, uncertainty estimation approaches that leverage the value of the Hessian become accurate and efficient because reverse-mode AD can be used for calculating the Hessian.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Quantum chemically calculated Abraham parameters for quantifying and predicting polymer hydrophobicity

The leakage and accumulation of plastic in the environment is a significant and growing problem with numerous detrimental impacts and has led to a push toward the design and development of more environmentally benign materials. To this end, we have developed a quantum chemistry-based model for predicting the mobility of polymer materials from molecular structure. Hydrophobicity is used as a surrogate for mobility given that hydrophobic interactions drive much of the partitioning of contaminants in and out of various environmentally relevant compartments. To model polymer hydrophobicity, we adjusted a previously developed Quantum Chemically Calculated Abraham Parameter model to calculate Abraham parameters of small molecules from molecular structure information. The resulting model predicted the octanol-water partition coefficient (K OW ) of polymer repeating units with a root mean square error (RMSE) of 0.48 (log scale). Additionally, the hydrophobicity of high molecular weight polymer materials was captured through solubility parameters and Nile red staining experiments from the literature and predicted with RMSEs of 1.21 (J/cc) 0.5 and 3.42 nm, respectively. Finally, to test the environmental applicability of the model, the relative adsorption capacity of three polymers was predicted and used to unify sorption isotherms across multiple sorbates and polymer sorbents.

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

Generalization Across Experimental Parameters in Neural Network Analysis of High-Resolution Transmission Electron Microscopy Datasets

Neural networks are promising tools for high-throughput and accurate transmission electron microscopy (TEM) analysis of nanomaterials, but are known to generalize poorly on data that is “out-of-distribution” from their training data. Given the limited set of image features typically seen in high-resolution TEM imaging, it is unclear which images are considered out-of-distribution from others. Here, we investigate how the choice of metadata features in the training dataset influences neural network performance, focusing on the example task of nanoparticle segmentation. We train and validate neural networks across curated, experimentally collected high-resolution TEM image datasets of nanoparticles under various imaging and material parameters, including magnification, dosage, nanoparticle diameter, and nanoparticle material. Overall, we find that our neural networks are not robust across microscope parameters, but do generalize across certain sample parameters. Additionally, data preprocessing can have unintended consequences on neural network generalization. Our results highlight the need to understand how dataset features affect deployment of data-driven algorithms.

42 ENGINEERING↗
DOE OSTI2023

Locoregional Treatment of Glioblastoma With Targeted α Therapy: [213Bi]Bi-DOTA–Substance P Versus [225Ac]Ac-DOTA–Substance P—Analysis of Influence Parameters

Background Glioblastoma (GB) is the most malignant primary brain tumor. Therefore, introduction of new treatment options is critically important. The aim of this study was to assess local treatment with α emitters [ 213 Bi]Bi-DOTA–substance P (SP) and [ 225 Ac]Ac-DOTA-SP. Methods Treatment was performed as salvage therapy in patients with recurrent primary and secondary GB. [ 213 Bi]Bi-DOTA-SP with injected activity 1.85 GBq per cycle was used in 20 primary (48.2 ± 11.8 years old) and in 9 secondary (38.8 ± 10.8 years old) GB patients and [ 225 Ac]Ac-DOTA-SP in 15 primary (45.1 ± 9.9 years old) and in 6 secondary (37.8 ± 6.4 years old) GB patients with a dose escalation scheme (10, 20, and 30 MBq). Results Local treatment with [ 213 Bi]Bi-DOTA-SP and [ 225 Ac]Ac-DOTA-SP was well tolerated with only few adverse effects. There was no statistically significant difference between [ 213 Bi]Bi-DOTA-SP and [ 225 Ac]Ac-DOTA-SP groups in survival parameters. For primary GB, survival parameters of patients treated with [ 213 Bi]Bi-DOTA-SP and [ 225 Ac]Ac-DOTA-SP were as follows(in months): progression-free survival time, 2.7 versus 2.4; OS-d (overall survival from time of diagnosis to death from any cause), 23.6 versus 21.0; OS-t (overall survival from the start of treatment to death from any cause), 7.5 versus 5.0; and OS-r (overall survival from recurrence in primary tumors to death from any cause), 10.9 versus 12.0. Survival parameters of secondary GB patients treated with [ 213 Bi]Bi-DOTA-SP and [ 225 Ac]Ac-DOTA-SP were as follows (in months): progression-free survival time, 5.8 versus 2.4; OS-d, 52.3 versus 65.0; OS-t, 16.4 versus 16.0; and OS-c (overall survival from conversion into secondary GB multiforme to death from any cause), 18.4 versus 36.0. Conclusions The similarity results of 213 Bi or 225 Ac may suggest that the local treatment of brain tumors can be greatly simplified. The experience to date shows that local radioisotope treatment of brain tumors requires further dosimetry studies, taking into account the complexity of biological processes.

Radiology, Nuclear Medicine & Medical Imaging↗
DOE OSTI2026

Critical structural parameter determining magnetic phases in the F⁡e 2 ⁢M⁢o 3 ⁢O 8 altermagnet system

A systematic structural investigation of F⁡e 2 ⁢M⁢o 3⁢ O 8 as a function of pressure, temperature, and magnetic field reveals that the 𝑃⁢6 3 ⁢𝑚⁢𝑐 space group remains robust over a wide range of conditions. No changes in the long-range crystal structure are observed for pressures up to 10 GPa, temperatures between 11 K and 300 K, and magnetic fields up to 9 T. The magnetostructural response, quantified by Δ⁢𝑐/𝑐, is determined for magnetic fields applied transverse to the crystallographic 𝑐 axis, demonstrating strong magnetoelastic coupling. The well-known magnetic-field-induced transition is confirmed to be first order and isostructural, occurring between two distinct altermagnetic states. Importantly, this transition can also be accessed using magnetic fields applied within the 𝑎𝑏 plane. When expressed in terms of the 𝑐/𝑎 ratio, the structural evolution under pressure (0–10 GPa) maps continuously onto that observed across the full Zn substitution range in (F⁡e 1−𝑦⁢ Z⁢n 𝑦 )2 M⁢o 3 ⁢O 8 ⁢(0 ≤ 𝑦 ≤ 1), indicating a common underlying structural control parameter. Furthermore, these results demonstrate that the 𝑐-axis lattice parameter is the primary structural degree of freedom governing magnetic behavior under external tuning parameters. More broadly, they indicate that magnetic order in this family of complex oxides (𝐴 2 ⁢M⁢o 3 ⁢O 8 , 𝐴 = Co,Mn,Ni) can be systematically tuned via pressure. This tunability highlights these materials as promising platforms for controllable magnetic switching, particularly in thin films integrated with piezoelectric substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rethinking 𝛼−RuCl 3 : Parameters, models, and phase diagram

RuCl 3 was likely the first ever deliberately synthesized ruthenium compound, following the discovery of the 44 Ru element in 1844. For a long time it was known as an oxidation catalyst, with its physical properties being discrepant and confusing, until a decade ago when its allotropic form 𝛼−RuCl 3 rose to exceptional prominence. This “rediscovery” of 𝛼−RuCl 3 has not only reshaped the hunt for a material manifestation of the Kitaev spin liquid, but it has opened the floodgates of theoretical and experimental research in the many unusual phases and excitations that the anisotropic-exchange magnets as a class of compounds have to offer. Given its importance for the field of Kitaev materials, it is astonishing that the low-energy spin model that describes this compound and its possible proximity to the much-desired spin-liquid state is still a subject of significant debate ten years later. In the present study, we argue that the existing key phenomenological observations put strong natural constraints on the effective microscopic spin model of 𝛼−RuCl 3 , and specifically on its spin-orbit-induced anisotropic-exchange parameters that are responsible for the nontrivial physical properties of this material. These constraints allow one to focus on the relevant region of the multidimensional phase diagram of the 𝛼−RuCl 3 model, suggest an intuitive description of it via a different parametrization of the exchange matrix, offer a unifying view on the earlier assessments of its parameters, and bring closer together several approaches to the derivation of anisotropic-exchange models. We explore extended phase diagrams relevant to the 𝛼−RuCl 3 parameter space using quasiclassical, Luttinger-Tisza, exact diagonalization, and density-matrix renormalization-group methods, demonstrating a remarkably close quantitative accord between them on the general structure and hierarchy of the phases, with the zigzag, ferromagnetic, and incommensurate phases that are proximate to each other. As a result, one of the highlights is the detailed agreement on the nature of the incommensurate phases that realize two distinct counterrotating helical states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Suppression of superconducting parameters by correlated quasi-two-dimensional magnetic fluctuations

In this work, we consider a clean layered magnetic superconductor in which a continuous magnetic transition takes place inside a superconducting state. We assume that the exchange interaction between superconducting and magnetic subsystems is weak so that superconductivity is not destroyed at the magnetic transition. A representative example of such material is Rb Eu Fe 4 As 4 . We investigate the suppression of the superconducting gap and superfluid density by correlated magnetic fluctuations in the vicinity of the magnetic transition. The influence of nonuniform exchange field on superconducting parameters is very sensitive to the relation between the magnetic correlation length ξ h and superconducting coherence length ξ s defining the ‘scattering’ ( ξ h < ξ s ) and ‘smooth’ ( ξ h > ξ s ) regimes. As a small uniform exchange field does not affect the superconducting gap and superfluid density at zero temperature, smoothening of the spatial variations of the exchange field reduces its effects on these parameters. We develop a quantitative description of this ‘scattering-to-smooth’ crossover for the case of quasi-two-dimensional magnetic fluctuations realized in Rb Eu Fe 4 As 4 . Since the magnetic-scattering energy scale is comparable with the gap in the crossover region, the standard quasiclassical approximation is not applicable and full microscopic treatment is required. We find that the corrections to both the gap and superfluid density increase proportionally to ξ h until it remains much smaller than ξ s . In the opposite limit, when the correlation length exceeds the coherence length both parameters have much weaker dependence on ξ h . Moreover, the gap correction may decrease with increasing of ξ h in the immediate vicinity of the magnetic transition if it is located at temperature much lower than the superconducting transition. We also find that the crossover between the two regimes is unexpectedly broad: The standard scattering approximation becomes sufficient only when ξ h is substantially smaller than ξ s .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measurement of the sound velocity and Grüneisen parameter of polystyrene at inertial confinement fusion conditions

The principal Hugoniot, sound velocity, and Grüneisen parameter of polystyrene were measured at conditions relevant to shocks in inertial confinement fusion implosions, from 100 to 1000 GPa. The sound velocity is in good agreement with quantum molecular dynamics calculations and all tabular equation of state models at pressures below 200 GPa. Above 200 GPa, the experimental results agree with two of the examined tables, but do not agree with the most recent table developed for design of ICF experiments. The Grüneisen parameter increases with density below ~3.1 g/cm 3 and approaches the asymptotic value for an ideal gas after complete dissociation. This behavior is in good agreement with quantum molecular dynamics results and previous work but is not represented by any of the tabular models. Finally, the discrepancy between tabular models and experimental measurement of the sound velocity and Grüneisen parameter is sufficient to impact simulations of ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Non-Gaussian transverse momentum fluctuations from impact parameter fluctuations

The transverse momentum per particle, [p t ], fluctuates event by event in ultrarelativistic nucleus-nucleus collisions, for a given multiplicity. These fluctuations are small and approximately Gaussian, but a nonzero skewness has been predicted on the basis of hydrodynamic calculations, and seen experimentally. We argue that the mechanism driving the skewness is that, if the system thermalizes, the mean transverse momentum increases with impact parameter for a fixed collision multiplicity. We postulate that fluctuations are Gaussian at fixed impact parameter, and that non-Gaussianities solely result from impact parameter fluctuations. Using recent data on the variance of [p t ] fluctuations, we make quantitative predictions for their skewness and kurtosis as a function of the collision multiplicity. Here, we predict, in particular, a spectacular increase of the skewness below the knee of the multiplicity distribution, followed by a fast decrease.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Impact of cross-section uncertainties on supernova neutrino spectral parameter fitting in the Deep Underground Neutrino Experiment

A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the 𝒪⁡(10) MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the 𝜈 𝑒 component of the supernova flux, enabling a wide variety of physics and astrophysics measurements. A key requirement for a correct interpretation of these measurements is a good understanding of the energy-dependent total cross section 𝜎⁡(𝐸 𝜈 ) for charged-current 𝜈 𝑒 absorption on argon. In the context of a simulated extraction of supernova 𝜈 𝑒 spectral parameters from a toy analysis, we investigate the impact of 𝜎⁡(𝐸 𝜈 ) modeling uncertainties on DUNE’s supernova neutrino physics sensitivity for the first time. We find that the currently large theoretical uncertainties on 𝜎⁡(𝐸 𝜈 ) must be substantially reduced before the 𝜈𝑒 flux parameters can be extracted reliably; in the absence of external constraints, a measurement of the integrated neutrino luminosity with less than 10% bias with DUNE requires 𝜎⁡(𝐸 𝜈 ) to be known to about 5%. The neutrino spectral shape parameters can be known to better than 10% for a 20% uncertainty on the cross-section scale, although they will be sensitive to uncertainties on the shape of 𝜎⁡(𝐸 𝜈 ). A direct measurement of low-energy 𝜈𝑒-argon scattering would be invaluable for improving the theoretical precision to the needed level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Impact of recent updates to neutrino oscillation parameters on the effective Majorana neutrino mass in 0 ν β β decay

We investigate how recent updates to neutrino oscillation parameters and the sum of neutrino masses influence the sensitivity of neutrinoless double-beta ( 0 ν β β ) decay experiments. Incorporating the latest cosmological constraints on the sum of neutrino masses and laboratory measurements on oscillations, we determine the sum of neutrino masses for both the normal hierarchy (NH) and the inverted hierarchy (IH). Our analysis reveals a narrow range for the sum of neutrino masses, approximately 0.06 eV / c 2 for NH and 0.102 eV / c 2 for IH. Utilizing these constraints, we calculate the effective Majorana masses for both NH and IH scenarios, establishing the corresponding allowed regions. Importantly, we find that the minimum neutrino mass is nonzero, as constrained by the current oscillation parameters. Additionally, we estimate the half-life of 0 ν β β decay using these effective Majorana masses for both NH and IH. Our results suggest that upcoming ton-scale experiments will comprehensively explore the IH scenario, while 100-ton-scale experiments will effectively probe the parameter space for the NH scenario, provided the background index can achieve 1 event/kton-year in the region of interest. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Bendability parameter for twisted ribbons to describe longitudinal wrinkling and delineate the near-threshold regime

We propose a dimensionless bendability parameter, ε -1 =[(h/W) 2 ⁢T -1 ] -1 , for wrinkling of thin, twisted ribbons with thickness h, width W, and tensional strain T. Bendability permits efficient collapse of data for wrinkle onset, wavelength, critical stress, and residual stress, demonstrating longitudinal wrinkling's primary dependence on this parameter. This parameter also allows us to distinguish the highly bendable range (ε -1 >20) from moderately bendable samples (ε -1 ϵ(0,20]). We identify scaling relations to describe longitudinal wrinkles that are valid across our entire set of simulated ribbons. When restricted to the highly bendable regime, simulations confirm theoretical near-threshold (NT) predictions for wrinkle onset and wavelength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Domain Wall Reactions in Multiple-Order Parameter Ferroelectrics

The motion of domain walls is crucial for ferroelectric switching. Conventionally, the switching dynamics is believed to be determined by the motion of one or a few low-energy domain wall types of dominant population. Here, we challenge this conventional idea in multiple-order-parameter ferroelectrics. Using hafnia as example, we show that the multiple-order-parameter nature not only provides various mobile domain walls and defectlike immobile domain walls, but also enables the domain wall reactions. In analogy with chemical reactions where substances react to form new substances, domain walls could also react to form other domain walls during switching. We identify several elementary domain wall reaction types including synthesis, decomposition, and exchange reactions. Domain walls are continually changed by these reactions during switching so that the switching behavior reflects the statistical average of many domain wall types with distinct mobility and stability. These reactions also lead to phenomenon like remanent nuclei and defect-site nucleation that facilitate switching and lower coercive field. Finally, the concept of domain wall reaction is not limited to hafnia but can be generalized to any multiple-order-parameter ferroelectric. As a result, this Letter conceptualizes domain wall reaction, expands the theory of ferroelectric switching, and suggests a practical way for defect engineering to control switching behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Automated prediction of lattice parameters from X-ray powder diffraction patterns

A key step in the analysis of powder X-ray diffraction (PXRD) data is the accurate determination of unit-cell lattice parameters. This step often requires significant human intervention and is a bottleneck that hinders efforts towards automated analysis. This work develops a series of one-dimensional convolutional neural networks (1D-CNNs) trained to provide lattice parameter estimates for each crystal system. A mean absolute percentage error of approximately 10% is achieved for each crystal system, which corresponds to a 100- to 1000-fold reduction in lattice parameter search space volume. The models learn from nearly one million crystal structures contained within the Inorganic Crystal Structure Database and the Cambridge Structural Database and, due to the nature of these two complimentary databases, the models generalize well across chemistries. A key component of this work is a systematic analysis of the effect of different realistic experimental non-idealities on model performance. It is found that the addition of impurity phases, baseline noise and peak broadening present the greatest challenges to learning, while zero-offset error and random intensity modulations have little effect. However, appropriate data modification schemes can be used to bolster model performance and yield reasonable predictions, even for data which simulate realistic experimental non-idealities. In order to obtain accurate results, a new approach is introduced which uses the initial machine learning estimates with existing iterative whole-pattern refinement schemes to tackle automated unit-cell solution.

42 ENGINEERING↗
DOE OSTI2026

A Machine Learning based Approach of Estimating Equivalent Circuit Model Parameters at Different SoCs of Li-ion Batteries from Voltage Relaxation

Abstract: In this study, an approach of estimating the equivalent circuit model (ECM) parameters for Li-ion batteries (LIBs) is proposed based on the voltage value at different intervals while relaxing the LIB after discharge. The typical approach for estimating ECM parameters of a LIB is to conduct electrochemical impedance spectroscopy (EIS) measurements at different frequencies and fit them to a predefined circuit model, which requires additional measuring arrangements and specialized devices. The proposed methodology utilizes four different voltages at 0s, 60s, 360s, and 1800s alongside the specific state of charge (SoC) value for a specific constant discharge current value of ~1C until the relaxation stage to train and evaluate three regression-based machine learning models— Support Vector Regression (SVR), Extreme Gradient Boosting (XGBoost), and Gaussian Process Regression (GPR)—for estimating the ECM parameters of the selected model. Bayesian optimization is employed for hyperparameter tuning to achieve optimal performance for all the regressor models, among which, the GPR provided the best performance with the root-mean-squared error (RMSE) of less than 4x10-4 on average for the resistive components and less than 0.27 for capacitive components with excellent R2 scores. The simplicity of the approach enables it to eliminate the need for sophisticated measuring equipment and computation power.

Sagar, Md. Samiul [The University of Alabama (UA)]↗
DOE OSTI2022

Two-Stage Optimization Framework for Detecting and Correcting Parameter Cyber-Attacks in Power System State Estimation

One major tool of Energy Management Systems for monitoring the status of the power grid is State Estimation. Since the results of state estimation are used within the energy management system, the security of the state estimation process is most important. The focus research in this area is on detecting False Data Injection attacks on measurements. While this is important, State Estimation also rely on database that are used to describe the relationship between measurements and systems' states. This paper presents a two-stage programming framework to detect and correct attacks in the parameters of the measurement model used by the state estimation process in the Energy Management System. In the first stage, an estimate of the line parameters ratios are obtained. In the second stage, the estimated ratios from stage I are used in a Bi-Level model for obtaining a final estimate of the measurements' model parameters. Hence, the presented framework does not only unify the detection and correction in a single optimization run, but also provide a monitoring scheme for the SE database that is typically considered static. In addition, in the two stages, linear programming framework is preserved. For validation, the IEEE 118 bus system is used for implementation. The results of this paper illustrate the effectiveness of the proposed model for detecting attacks in the database used in the state estimation process.

state estimation, two-stage optimization, cyber-ph↗