Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FIT”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

An extension of high-order shock-fitted detonation propagation in explosives

An update to the numerical shock-fitting scheme from [1] is proposed. The updated scheme locally lowers the order of the scheme in the region of secondary non-fitted waves that interact with the fitted front. In addition, it utilizes a non-uniform mapping from the physical space to the computational space along the streamline direction to offer computational savings over the use of an uniform spacing. The utility of the scheme is first demonstrated on planar ideal gas detonating flows. It is then applied to detonation propagation in a two-dimensional axisymmetric geometry while utilizing an Arrhenius Wescott-Stewart-Davis model for PBX 9502.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

WeakIdent: Weak formulation for identifying differential equation using narrow-fit and trimming

Data-driven identification of differential equations is an interesting but challenging problem, especially when the given data are corrupted by noise. When the governing differential equation is a linear combination of various differential terms, the identification problem can be formulated as solving a linear system, with the feature matrix consisting of linear and nonlinear terms multiplied by a coefficient vector. This product is equal to the time derivative term, and thus generates dynamical behaviors. The goal is to identify the correct terms that form the equation to capture the dynamics of the given data. We propose a general and robust framework to recover differential equations using a weak formulation with two new mechanisms, narrow-fit and trimming, for both ordinary and partial differential equations (ODEs and PDEs). The weak formulation facilitates an efficient and robust way to handle noise, and two new mechanisms, narrow-fit and trimming, improve the coefficient support and value recoveries respectively. For each sparsity level, Subspace Pursuit is utilized to find an initial set of support from the large dictionary. Then, we focus on highly dynamic regions (rows of the feature matrix), and error normalize the feature matrix in the narrow-fit step. The support is further updated via trimming the terms that contribute the least. Finally, the support set of features with the smallest Cross-Validation error is chosen as the result. A comprehensive set of numerical experiments are presented for both systems of ODEs and PDEs with various noise levels. The proposed method gives a robust recovery of the coefficients, and a significant denoising effect which can handle up to 100% noise-to-signal ratio for some equations. We compare the proposed method with several state-of-the-art algorithms for the recovery of differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Model Fit to Experimental Data for Foam-Assisted Deep Vadose Zone Remediation

Foam has been regarded as a promising means of remeidal amendment delivery to overcome subsurface heterogeneity in subsurface remediation processes. This study investigates how a foam model, developed by Method of Characteristics and fractional flow analysis in the companion paper of Roostapour and Kam (2012), can be applied to make a fit to a set of existing laboratory flow experiments (Zhong et al., 2009) in an application relevant to deep vadose zone remediation. This study reveals a few important insights regarding foam-assisted deep vadose zone remediation: (i) the mathematical framework established for foam modeling can fit typical flow experiments matching wave velocities, saturation history , and pressure responses; (ii) the set of input parameters may not be unique for the fit, and therefore conducting experiments to measure basic model parameters related to relative permeability, initial and residual saturations, surfactant adsorption and so on should not be overlooked; and (iii) gas compressibility plays an important role for data analysis, thus should be handled carefully in laboratory flow experiments. Foam kinetics, causing foam texture to reach its steady-state value slowly, may impose additional complications.

Roostapour, A.↗
DOE OSTI2023

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.

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

Physics-Driven Construction of Compact Primitive Gaussian Density Fitting Basis Sets

We present a model-assisted density fitting (MADF) basis set generator, an algorithm for generating primitive atomic Gaussian density fitting (DF) basis sets (DFBSs) from a contracted Gaussian orbital basis set (OBS). The MADF algorithm produces DFBSs suitable for accurate robust DF approximation of 2-particle interactions in mean-field and correlated electronic structures. The algorithm is designed to (a) saturate the OBS product space by a large regularized set of primitive solid-harmonic Gaussian shells with nonuniform distribution of exponents, followed by (b) pruning of the shells according to their contributions to the 2- body energy of a correlated atomic ensemble. Building the DFBS generator model almost exclusively on mathematical and physical principles allows one to limit the number of parameters that control the density fitting error to three, with a single set of parameters sufficient for computations with all basis cardinal numbers, with and without correlation of core electrons, with and without scalar and spin-dependent relativistic effects, spanning almost all of the Periodic Table. Performance assessment included basis sets up to quadruple-ζ quality from several major basis set families, using molecules composed of main-group, d-block, and f-block elements. The resulting DF errors in Hartree−Fock and second-order MP2 energies (with relativistic all-electron treatments, when appropriate) were on the order of 20 and 10 μE h per electron, respectively.

Approximation↗
DOE OSTI2016

What is the best method to fit time-resolved data? A comparison of the residual minimization and the maximum likelihood techniques as applied to experimental time-correlated, single-photon counting data

The need for measuring fluorescence lifetimes of species in subdiffraction-limited volumes in, for example, stimulated emission depletion (STED) microscopy, entails the dual challenge of probing a small number of fluorophores and fitting the concomitant sparse data set to the appropriate excited-state decay function. This need has stimulated a further investigation into the relative merits of two fitting techniques commonly referred to as “residual minimization” (RM) and “maximum likelihood” (ML). Fluorescence decays of the well-characterized standard, rose bengal in methanol at room temperature (530 ± 10 ps), were acquired in a set of five experiments in which the total number of “photon counts” was approximately 20, 200, 1000, 3000, and 6000 and there were about 2–200 counts at the maxima of the respective decays. Each set of experiments was repeated 50 times to generate the appropriate statistics. Each of the 250 data sets was analyzed by ML and two different RM methods (differing in the weighting of residuals) using in-house routines and compared with a frequently used commercial RM routine. Convolution with a real instrument response function was always included in the fitting. While RM using Pearson’s weighting of residuals can recover the correct mean result with a total number of counts of 1000 or more, ML distinguishes itself by yielding, in all cases, the same mean lifetime within 2% of the accepted value. For 200 total counts and greater, ML always provides a standard deviation of <10% of the mean lifetime, and even at 20 total counts there is only 20% error in the mean lifetime. Here, the robustness of ML advocates its use for sparse data sets such as those acquired in some subdiffraction-limited microscopies, such as STED, and, more importantly, provides greater motivation for exploiting the time-resolved capacities of this technique to acquire and analyze fluorescence lifetime data.

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

Fitness Landscape-Guided Engineering of Locally Supercharged Virus-like Particles with Enhanced Cell Uptake Properties

Protein-based nanoparticles are useful models for the study of self-assembly and attractive candidates for drug delivery. Virus-like particles (VLPs) are especially promising platforms for expanding the repertoire of therapeutics that can be delivered effectively as they can deliver many copies of a molecule per particle for each delivery event. However, their use is often limited due to poor uptake of VLPs into mammalian cells. In this study, we use the fitness landscape of the bacteriophage MS2 VLP as a guide to engineer capsid variants with positively charged surface residues to enhance their uptake into mammalian cells. By combining mutations with positive fitness scores that were likely to produce assembled capsids, we identified two key double mutants with internalization efficiencies as much as 67-fold higher than that of wtMS2. Internalization of these variants with positively charged surface residues depends on interactions with cell surface sulfated proteoglycans, and yet, they are biophysically similar to wtMS2 with low cytotoxicity and an overall negative charge. Additionally, the best-performing engineered MS2 capsids can deliver a potent anticancer small-molecule therapeutic with efficacy levels similar to antibody-drug conjugates. Through this work, we were able to establish fitness landscape-based engineering as a successful method for designing VLPs with improved cell penetration. These findings suggest that VLPs with positive surface charge could be useful in improving the delivery of small-molecule- and nucleic acid-based therapeutics.

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

SDSS-IV MaStar: Stellar parameter determination with continuum-supplemented full-spectrum fitting

We present a stellar parameter catalog built to accompany the MaStar Stellar Library, which is a comprehensive collection of empirical, medium-resolution stellar spectra. We constructed this parameter catalog by using a multicomponent χ 2 fitting approach to match MaStar spectra to models generated by interpolating the ATLAS9-based BOSZ model spectra. The total χ 2 for a given model is defined as the sum of components constructed to characterize narrow-band features of observed spectra (e.g., absorption lines) and the broadband continuum shape separately. Extinction and systematics due to flux calibration were taken into account in the fitting. The χ 2 distribution for a given region of model space was sampled using a Markov chain Monte Carlo (MCMC) algorithm, the data from which were then used to extract atmospheric parameter estimates (T eff , log g, [Fe/H], and [α/Fe]), their corresponding uncertainties, and direct extinction measurements. Two methods were used to extract parameters and uncertainties: one that accepts the MCMC’s prescribed minimum-χ 2 result, and one that uses Bayesian inference to compute a likelihood-weighted mean from the χ 2 distribution sampled by the MCMC. Results were evaluated for internal consistency using repeat observations where available and by comparing them with external data sets (e.g., APOGEE-2 and Gaia DR2). Our spectral-fitting exercise reveals possible deficiencies in current theoretical model spectra, illustrating the potential power of MaStar spectra for helping to improve the models. This paper represents an update to the parameters that were originally presented with SDSS-IV DR17. The MaStar parameter catalog containing our BestFit results is available on the SDSS-IV DR17 website as part of version 2 of the MaStar stellar parameter value-added catalog.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Density fitting for three-electron integrals in explicitly correlated electronic structure theory

The principal challenge in using explicitly correlated wavefunctions for molecules is the evaluation of nonfactorizable integrals over the coordinates of three or more electrons. Immense progress was made in tackling this problem through the introduction of a single-particle resolution of the identity. Decompositions of sufficient accuracy can be achieved, but only with large auxiliary basis sets. Density fitting is an alternative integral approximation scheme, which has proven to be very reliable for two-electron integrals. Here, we extend density fitting to the treatment of all three-electron integrals that appear at the MP2-F12/3*A level of theory. We demonstrate that the convergence of energies with respect to auxiliary basis size is much more rapid with density fitting than with the traditional resolution-of-the-identity approach.

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

Magnetizability and rotational g tensors for density fitted local second-order Møller-Plesset perturbation theory using gauge-including atomic orbitals

In this paper, we present theory and implementation of an efficient program for calculating magnetizabilities and rotational g tensors of closed-shell molecules at the level of local second-order Møller-Plesset perturbation theory (MP2) using London orbitals. Density fitting is employed to factorize the electron repulsion integrals with ordinary Gaussians as fitting functions. The presented program for the calculation of magnetizabilities and rotational g tensors is based on a previous implementation of NMR shielding tensors reported by S. Loibl and M. Schütz [J. Chem. Phys. 137, 084107 (2012)]. Extensive test calculations show (i) that the errors introduced by density fitting are negligible, and (ii) that the errors of the local approximation are still rather small, although larger than for nuclear magnetic resonance (NMR) shielding tensors. Electron correlation effects for magnetizabilities are tiny for most of the molecules considered here. MP2 appears to overestimate the correlation contribution of magnetizabilities such that it does not constitute an improvement over Hartree-Fock (when comparing to higher-order methods like CCSD(T)). For rotational g tensors the situation is different and MP2 provides a significant improvement in accuracy over Hartree-Fock. The computational performance of the new program was tested for two extended systems, the larger comprising about 2200 basis functions. It turns out that a magnetizability (or rotational g tensor) calculation takes about 1.5 times longer than a corresponding NMR shielding tensor calculation.

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

Quantum chi-squared and goodness of fit testing

A quantum mechanical hypothesis test is presented for the hypothesis that a certain setup produces a given quantum state. Although the classical and the quantum problems are very much related to each other, the quantum problem is much richer due to the additional optimization over the measurement basis. A goodness of fit test for i.i.d quantum states is developed and a max-min characterization for the optimal measurement is introduced. We find the quantum measurement which leads both to the maximal Pitman and Bahadur efficiencies, and determine the associated divergence rates. We discuss the relationship of the quantum goodness of fit test to the problem of estimating multiple parameters from a density matrix. These problems are found to be closely related and we show that the largest error of an optimal strategy, determined by the smallest eigenvalue of the Fisher information matrix, is given by the divergence rate of the goodness of fit test.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Permutationally invariant fitting of intermolecular potential energy surfaces: A case study of the Ne-C 2 H 2 system

Here, the permutation invariant polynomial-neural network (PIP-NN) approach is extended to fit intermolecular potential energy surfaces (PESs). Specifically, three PESs were constructed for the Ne-C 2 H 2 system. PES1 is a full nine-dimensional PIP-NN PES directly fitted to ~42 000 ab initio points calculated at the level of CCSD(T)-F12a/cc-pCVTZ-F12, while the other two consist of the six-dimensional PES for C 2 H 2 [H. Han, A. Li, and H. Guo, J. Chem. Phys. 141, 244312 (2014)] and an intermolecular PES represented in either the PIP (PES2) or PIP-NN (PES3) form. The comparison of fitting errors and their distributions, one-dimensional cuts and two-dimensional contour plots of the PESs, as well as classical trajectory collisional energy transfer dynamics calculations shows that the three PESs are very similar. We conclude that full-dimensional PESs for non-covalent interacting molecular systems can be constructed efficiently and accurately by the PIP-NN approach for both the constituent molecules and intermolecular parts.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

MB-Fit: Software infrastructure for data-driven many-body potential energy functions

Many-body potential energy functions (MB-PEFs), which integrate data-driven representations of many-body short-range quantum mechanical interactions with physics-based representations of many-body polarization and long-range interactions, have recently been shown to provide high accuracy in the description of molecular interactions from the gas to the condensed phase. Here, we present MB-Fit, a software infrastructure for the automated development of MB-PEFs for generic molecules within the TTM-nrg (Thole-type model energy) and MB-nrg (many-body energy) theoretical frameworks. Besides providing all the necessary computational tools for generating TTM-nrg and MB-nrg PEFs, MB-Fit provides a seamless interface with the MBX software, a many-body energy and force calculator for computer simulations. Given the demonstrated accuracy of the MB-PEFs, particularly within the MB-nrg framework, we believe that MB-Fit will enable routine predictive computer simulations of generic (small) molecules in the gas, liquid, and solid phases, including, but not limited to, the modeling of quantum isomeric equilibria in molecular clusters, solvation processes, molecular crystals, and phase diagrams.

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

Quantitative assessment of fitting errors associated with streak camera noise in Thomson scattering data analysis

Thomson scattering measurements in high energy density experiments are often recorded using optical streak cameras. In the low-signal regime, noise introduced by the streak camera can become an important and sometimes the dominant source of measurement uncertainty. In this paper, we present a formal method of accounting for the presence of streak camera noise in our measurements. We present a phenomenological description of the noise generation mechanisms and present a statistical model that may be used to construct the covariance matrix associated with a given measurement. This model is benchmarked against simulations of streak camera images. We demonstrate how this covariance may then be used to weight fitting of the data and provide quantitative assessments of the uncertainty in the fitting parameters determined by the best fit to the data and build confidence in the ability to make statistically significant measurements in the low-signal regime, where spatial correlations in the noise become apparent. These methods will have general applicability to other measurements made using optical streak cameras.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Fit-free optical determination of electronic thermalization time in nematic iron-based superconductors

We present a nematic response function model (NRFM) for fit-free direct extraction of the characteristic time of ultrafast electronic thermalization in iron-based superconductors, materials with electronic nematicity. By combining the NRFM for polarization-dependent pump–probe measurements of electronic nematic response with the two-temperature model (TTM) for sub-picosecond quasiparticle relaxation, we quantify the electronic thermalization timescales and their anisotropy. The nematic response function is modeled as the difference in normalized reflectivity signals, revealing a pronounced sub-picosecond extremum in signal evolution that directly yields the characteristic electronic thermalization time. This method demonstrates that the NRFM is consistent with TTM fits of transient optical response, yielding electronic thermalization time constants on the order of 110–230 fs for the FeSe 1−x Te x and Ba(Fe 0.92 Co 0.08 ) 2 As 2 thin films. Here, the proposed approach can be applied to any material that exhibits electronic nematicity, providing a powerful tool for direct mapping of the relaxation time in nematic materials, avoiding complex experimental data-fitting procedures.

Bartenev, Alexander [University of Puerto Rico, Ma↗
DOE OSTI2013

ON THE ROBUSTNESS OF z = 0-1 GALAXY SIZE MEASUREMENTS THROUGH MODEL AND NON-PARAMETRIC FITS

We present the size-stellar mass relations of nearby (z = 0.01-0.02) Sloan Digital Sky Survey galaxies, for samples selected by color, morphology, Sérsic index n, and specific star formation rate. Several commonly employed size measurement techniques are used, including single Sérsic fits, two-component Sérsic models, and a non-parametric method. Through simple simulations, we show that the non-parametric and two-component Sérsic methods provide the most robust effective radius measurements, while those based on single Sérsic profiles are often overestimates, especially for massive red/early-type galaxies. Using our robust sizes, we show for all sub-samples that the mass-size relations are shallow at low stellar masses and steepen above ∼3-4 × 10{sup 10} M{sub ☉}. The mass-size relations for galaxies classified as late-type, low-n, and star-forming are consistent with each other, while blue galaxies follow a somewhat steeper relation. The mass-size relations of early-type, high-n, red, and quiescent galaxies all agree with each other but are somewhat steeper at the high-mass end than previous results. To test potential systematics at high redshift, we artificially redshifted our sample (including surface brightness dimming and degraded resolution) to z = 1 and re-fit the galaxies using single Sérsic profiles. The sizes of these galaxies before and after redshifting are consistent and we conclude that systematic effects in sizes and the size-mass relation at z ∼ 1 are negligible. Interestingly, since the poorer physical resolution at high redshift washes out bright galaxy substructures, single Sérsic fitting appears to provide more reliable and unbiased effective radius measurements at high z than for nearby, well-resolved galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Determining neutron star properties by fitting oblate-star waveform models to X-ray burst oscillations

We describe sophisticated new Bayesian analysis methods that make it possible to estimate quickly the masses and radii of rapidly rotating, oblate neutron stars by fitting oblate-star waveform models to energy-resolved observations of the waveforms of X-ray burst oscillations produced by such stars. We find that a 25% variation of the temperature of the hot spot in the north-south direction does not significantly bias estimates of the mass M and equatorial radius R{sub eq} derived by fitting a model that assumes a uniform-temperature spot. Our results show that fits of oblate-star waveform models to waveform data can simultaneously determine M and R{sub eq} with 1σ uncertainties ≲7% if (1) the star's rotation rate is ≳ 600 Hz; (2) the spot center and observer's sightline are both within 30{sup ∘} of the star's rotational equator; and (3) ∼10{sup 7} counts are collected during burst oscillations that have a fractional rms amplitude ≳10%. A fractional amplitude of ∼10% is realistic, and the accepted NICER and proposed LOFT and AXTAR space missions could collect 10{sup 7} counts from a single star by combining data from many X-ray bursts from the star. Uncertainties ≲7% are small enough to improve substantially our understanding of cold, ultradense matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The PDF4LHC21 combination of global PDF fits for the LHC Run III

A precise knowledge of the quark and gluon structure of the proton, encoded by the parton distribution functions (PDFs), is of paramount importance for the interpretation of high-energy processes at present and future lepton–hadron and hadron–hadron colliders. Motivated by recent progress in the PDF determinations carried out by the CT, MSHT, and NNPDF groups, we present an updated combination of global PDF fits: PDF4LHC21. It is based on the Monte Carlo combination of the CT18, MSHT20, and NNPDF3.1 sets followed by either its Hessian reduction or its replica compression. Extensive benchmark studies are carried out in order to disentangle the origin of the differences between the three global PDF sets. In particular, dedicated fits based on almost identical theory settings and input datasets are performed by the three groups, highlighting the role played by the respective fitting methodologies. We compare the new PDF4LHC21 combination with its predecessor, PDF4LHC15, demonstrating their good overall consistency and a modest reduction of PDF uncertainties for key LHC processes such as electroweak gauge boson production and Higgs boson production in gluon fusion. We study the phenomenological implications of PDF4LHC21 for a representative selection of inclusive, fiducial, and differential cross sections at the LHC. The PDF4LHC21 combination is made available via the LHAPDF library and provides a robust, user-friendly, and efficient method to estimate the PDF uncertainties associated to theoretical calculations for the upcoming Run III of the LHC and beyond.

quantum chromodynamics↗