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 19 records

Broadband distortion modeling in Lyman-α forest BAO fitting

Recently, the Lyman-α absorption observed in the spectra of high-redshift quasars has been used as a tracer of large-scale structure by means of the three-dimensional Lyman-α forest auto-correlation function at redshift z≃ 2.3, but the need to fit the quasar continuum in every absorption spectrum introduces a broadband distortion that is difficult to correct and causes a systematic error for measuring any broadband properties. Here, we describe a k-space model for this broadband distortion based on a multiplicative correction to the power spectrum of the transmitted flux fraction that suppresses power on scales corresponding to the typical length of a Lyman-α forest spectrum. In implementing the distortion model in fits for the baryon acoustic oscillation (BAO) peak position in the Lyman-α forest auto-correlation, we find that the fitting method recovers the input values of the linear bias parameter b F and the redshift-space distortion parameter β F for mock data sets with a systematic error of less than 0.5%. Applied to the auto-correlation measured for BOSS Data Release 11, our method improves on the previous treatment of broadband distortions in BAO fitting by providing a better fit to the data using fewer parameters and reducing the statistical errors on βF and the combination b F (1+β F ) by more than a factor of seven. The measured values at redshift z=2.3 are βF=1.39 +0.11 +0.24 +0.38 -0.10 -0.19 -0.28 and bF(1+βF)=-0.374 +0.007 +0.013 +0.020 -0.007 -0.014 -0.022 (1σ, 2σ and 3σ statistical errors). Our fitting software and the input files needed to reproduce our main results are publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying and Reducing Curve-Fitting Uncertainty in Isc: Preprint

Current-voltage (I-V) curve measurements of photovoltaic (PV) devices are used to determine performance parameters and to establish traceable calibration chains. Measurement standards specify localized curve fitting methods, e.g., straight-line interpolation/extrapolation of the I-V curve points near short-circuit current, Isc. By considering such fits as statistical linear regressions, uncertainties in the performance parameters are readily quantified. However, the legitimacy of such a computed uncertainty requires that the model be a valid (local) representation of the I-V curve and that the noise be sufficiently well characterized. Using more data points often has the advantage of lowering the uncertainty. However, more data points can make the uncertainty in the fit arbitrarily small, and this fit uncertainty misses the dominant residual uncertainty due to so-called model discrepancy. Using objective Bayesian linear regression for straight-line fits for Isc, we investigate an evidence-based method to automatically choose data windows of I-V points with reduced model discrepancy. We also investigate noise effects. Uncertainties, aligned with the Guide to the Expression of Uncertainty in Measurement (GUM), are quantified throughout.

14 SOLAR ENERGY↗

Combined fit of spectrum and composition data as measured by the Pierre Auger Observatory

In this paper, we present a combined fit of a simple astrophysical model of UHECR sources to both the energy spectrum and mass composition data measured by the Pierre Auger Observatory. The fit has been performed for energies above 5 • 10 18 eV, i.e. the region of the all-particle spectrum above the so-called 'ankle' feature. The astrophysical model we adopted consists of identical sources uniformly distributed in a comoving volume, where nuclei are accelerated through a rigidity-dependent mechanism. The fit results suggest sources characterized by relatively low maximum injection energies, hard spectra and heavy chemical composition. We also show that uncertainties about physical quantities relevant to UHECR propagation and shower development have a non-negligible impact on the fit results.

79 ASTRONOMY AND ASTROPHYSICS↗

Two Outer Membrane Proteins Contribute to Caulobacter crescentus Cellular Fitness by Preventing Intracellular S-Layer Protein Accumulation

Surface layers, or S-layers, are two-dimensional protein arrays that form the outermost layer of many bacteria and archaea. They serve several functions, including physical protection of the cell from environmental threats. The high abundance of S-layer proteins necessitates a highly efficient export mechanism to transport the S-layer protein from the cytoplasm to the cell exterior. Caulobacter crescentus is unique in that it has two homologous, seemingly redundant outer membrane proteins, RsaF a and RsaF b , which together with other components form a type I protein translocation pathway for S-layer export. These proteins have homology to Escherichia coli TolC, the outer membrane channel of multidrug efflux pumps. Here we provide evidence that, unlike TolC, RsaF a and RsaF b are not involved in either the maintenance of membrane stability or the active export of antimicrobial compounds. Rather, RsaF a and RsaF b are required to prevent intracellular accumulation and aggregation of the S-layer protein RsaA; deletion of RsaF a and RsaF b led to a general growth defect and lowered cellular fitness. Using Western blotting, transmission electron microscopy, and transcriptome sequencing (RNA-seq), we show that loss of both RsaF a and RsaF b led to accumulation of insoluble RsaA in the cytoplasm, which in turn caused upregulation of a number of genes involved in protein misfolding and degradation pathways. These findings provide new insight into the requirement for RsaF a and RsaF b in cellular fitness and tolerance to antimicrobial agents and further our understanding of the S-layer export mechanism on both the transcriptional and translational levels in C. crescentus . IMPORTANCE Decreased growth rate and reduced cell fitness are common side effects of protein production in overexpression systems. Inclusion bodies typically form inside the cell, largely due to a lack of sufficient export machinery to transport the overexpressed proteins to the extracellular environment. This phenomenon can conceivably also occur in natural systems. As one example of a system evolved to prevent intracellular protein accumulation, our study demonstrates that Caulobacter crescentus has two homologous outer membrane transporter proteins that are involved in S-layer export. This is an interesting case study that demonstrates how bacteria can evolve redundancy to ensure adequate protein export functionality and maintain high cellular fitness. Moreover, we provide evidence that these two outer membrane proteins, although being the closest C. crescentus homologs to TolC in E. coli , do not process TolC functionality in C. crescentus .

59 BASIC BIOLOGICAL SCIENCES↗

A Fast Universal Kinematic Fitting Code for Low-Energy Nuclear Physics: FUNKI_FIT

We present an open-source kinematic fitting routine designed for low-energy nuclear physics applications. Although kinematic fitting is commonly used in high-energy particle physics, it is rarely used in low-energy nuclear physics, despite its effectiveness. A FORTRAN and ROOT C++ version of the FUNKI_FIT kinematic fitting code have been developed and published open access. The FUNKI_FIT code is universal in the sense that the constraint equations can be easily modified to suit different experimental set-ups and reactions. Two case studies for the use of this code, utilising experimental and Monte–Carlo data, are presented: (1) charged-particle spectroscopy using silicon-strip detectors; (2) charged-particle spectroscopy using active target detectors. The kinematic fitting routine provides an improvement in resolution in both cases, demonstrating, for the first time, the applicability of kinematic fitting across a range of nuclear physics applications. The ROOT macro has been developed in order to easily apply this technique in standard data analysis routines used by the nuclear physics community.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Crop‐associated virus infection in a native perennial grass: reduction in plant fitness and dynamic patterns of virus detection

Summary To understand the eco‐evolutionary significance of plant viruses in nature, we must (i) quantify the effects of infection on plant fitness and (ii) recognize that native plants are increasingly exposed to crop‐associated viruses. Studies of perennials are particularly needed: most of our knowledge of plant‐virus interactions is from annuals, yet long‐lived species dominate landscapes. Here we used aster models for life‐history analysis and longitudinal measures of plant virus status to evaluate multi‐year consequences of crop virus infection in a native perennial. We used Barley yellow dwarf virus acquired from wheat to inoculate seedlings of P anicum virgatum L. (switchgrass), a North American prairie grass. We grew inoculated and mock‐inoculated individuals of two ecotypes for 3 years in the field. We measured plant size, infection status and fitness components. Aster modelling provided integrated multi‐year measures of fitness. Crop virus inoculation reduced multi‐year native plant fitness by 30% over 2 years despite generally asymptomatic infection and evidence of resistance. This reduction was greater than predicted from individual fitness components or most size measures. Ecotypes differed in response, with the lowland ecotype experiencing higher apparent recovery from infection. Virus treatment in the upland ecotype delayed flowering phenology and reduced seed filling. Synthesis . Our use of field experimentation, surveys of plant infection status and aster modelling demonstrates a rigorous and broadly applicable approach for quantifying the effects of viruses and other microbes on multi‐year plant fitness. We found that a crop virus had negative multi‐year effects on native plant fitness even after infection was no longer detected. Viruses may have substantial effects on native vegetation with domestication of landscapes and agricultural expansion.

Alexander, Helen M.↗

Search for Low-Mass Dark Matter with CDMSlite Using a Profile Likelihood Fit

The Cryogenic Dark Matter Search low ionization threshold experiment (CDMSlite) searches for interactions between dark matter particles and germanium nuclei in cryogenic detectors. The experiment has achieved a low energy threshold with improved sensitivity to low-mass (<10 GeV/c2) dark matter particles. We present an analysis of the final CDMSlite data set, taken with a different detector than was used for the two previous CDMSlite data sets. This analysis includes a data "salting" method to protect against bias, improved noise discrimination, background modeling, and the use of profile likelihood methods to search for a dark matter signal in the presence of backgrounds. We achieve an energy threshold of 70 eV and significantly improve the sensitivity for dark matter particles with masses between 2.5 and 10 GeV/c2 compared to previous analyses. We set an upper limit on the dark matter-nucleon scattering cross section in germanium of 5.4×10-42 cm2 at 5 GeV/c2, a factor of ~2.5 improvement over the previous CDMSlite result.

Agnese, R↗

Automatic Multiple Experiment Simulation and Fitting (Ames-Fit)

AMES-Fit is a program used to automatically fit the multi-field solid-state NMR spectra of half-integer quadrupolar nuclei. Due to the high dimensional space, gradient algorithms have failed to address the fitting of such data, which is at present done manually. AMES-Fit diverges from these approaches by using an adaptive step size random search algorithm to fit the NMR spectra to consistently find the global best fit parameters.

Perras, Frederic↗

Semi-exact concentric atomic density fitting: Reduced cost and increased accuracy compared to standard density fitting

A local density fitting scheme is considered in which atomic orbital (AO) products are approximated using only auxiliary AOs located on one of the nuclei in that product. The possibility of variational collapse to an unphysical “attractive electron” state that can affect such density fitting [P. Merlot, T. Kjærgaard, T. Helgaker, R. Lindh, F. Aquilante, S. Reine, and T. B. Pedersen, J. Comput. Chem. 34, 1486 (2013)] is alleviated by including atom-wise semidiagonal integrals exactly. Our approach leads to a significant decrease in the computational cost of density fitting for Hartree–Fock theory while still producing results with errors 2–5 times smaller than standard, nonlocal density fitting. Our method allows for large Hartree–Fock and density functional theory computations with exact exchange to be carried out efficiently on large molecules, which we demonstrate by benchmarking our method on 200 of the most widely used prescription drug molecules. Our new fitting scheme leads to smooth and artifact-free potential energy surfaces and the possibility of relatively simple analytic gradients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cryo_fit : Democratization of flexible fitting for cryo-EM

Cryo-electron microscopy (cryo-EM) is becoming a method of choice for describing native conformations of biomolecular complexes at high resolution. The rapid growth of cryo-EM in recent years has created a high demand for automated solutions, both in hardware and software. Flexible fitting of atomic models to three-dimensional (3D) cryo-EM reconstructions by molecular dynamics (MD) simulation is a popular technique but often requires technical expertise in computer simulation. This work introduces cryo_fit, a package for the automatic flexible fitting of atomic models in cryo-EM maps using MD simulation. Additionally, the package is integrated with the Phenix software suite. The module was designed to automate the multiple steps of MD simulation in a reproducible manner, as well as facilitate refinement and validation through Phenix. Through the use of cryo_fit, scientists with little experience in MD simulation can produce high quality atomic models automatically and better exploit the potential of cryo-EM.

59 BASIC BIOLOGICAL SCIENCES↗

Three dimensional equilibrium solutions for a current-carrying reversed-field pinch plasma with a close-fitting conducting shell

In order to characterize the Madison Symmetric Torus (MST) reversed-field pinch(RFP)plasmas that bifurcate to a helical equilibrium, the V3FIT equilibrium reconstruction code was modified to include a conducting boundary. RFPplasmas become helical at a high plasma current, which induces large eddy currents in MST's thick aluminum shell. The V3FIT conducting boundary accounts for the contribution from these eddy currents to external magnetic diagnostic coil signals. This implementation of V3FIT was benchmarked against MSTFit, a 2D Grad-Shafranov solver, for axisymmetric plasmas. The two codes both fit B measurement loops around the plasma minor diameter with qualitative agreement between each other and the measured field. Fits in the 3D case converge well, with q-profile and plasma shape agreement between two distinct toroidal locking phases. Greater than 60% of the measured n = 5 component of B at r = a is due to eddy currents in the shell, as calculated by the conducting boundary model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Alternative difference analysis scheme combining R -space EXAFS fit with global optimization XANES fit for X-ray transient absorption spectroscopy

Time-resolved X-ray absorption spectroscopy (TR-XAS), based on the laser-pump/X-ray-probe method, is powerful in capturing the change of the geometrical and electronic structure of the absorbing atom upon excitation. TR-XAS data analysis is generally performed on the laser-on minus laser-off difference spectrum. Here, a new analysis scheme is presented for the TR-XAS difference fitting in both the extended X-ray absorption fine-structure (EXAFS) and the X-ray absorption near-edge structure (XANES) regions. R -space EXAFS difference fitting could quickly provide the main quantitative structure change of the first shell. The XANES fitting part introduces a global non-derivative optimization algorithm and optimizes the local structure change in a flexible way where both the core XAS calculation package and the search method in the fitting shell are changeable. The scheme was applied to the TR-XAS difference analysis of Fe(phen) 3 spin crossover complex and yielded reliable distance change and excitation population.

Zhan, Fei↗

Pulmonary lobe segmentation based on ridge surface sampling and shape model fitting

Purpose: Performing lobe-based quantitative analysis of the lung in computed tomography (CT) scans can assist in efforts to better characterize complex diseases such as chronic obstructive pulmonary disease (COPD). While airways and vessels can help to indicate the location of lobe boundaries, segmentations of these structures are not always available, so methods to define the lobes in the absence of these structures are desirable. Methods: The authors present a fully automatic lung lobe segmentation algorithm that is effective in volumetric inspiratory and expiratory computed tomography (CT) datasets. The authors rely on ridge surface image features indicating fissure locations and a novel approach to modeling shape variation in the surfaces defining the lobe boundaries. The authors employ a particle system that efficiently samples ridge surfaces in the image domain and provides a set of candidate fissure locations based on the Hessian matrix. Following this, lobe boundary shape models generated from principal component analysis (PCA) are fit to the particles data to discriminate between fissure and nonfissure candidates. The resulting set of particle points are used to fit thin plate spline (TPS) interpolating surfaces to form the final boundaries between the lung lobes. Results: The authors tested algorithm performance on 50 inspiratory and 50 expiratory CT scans taken from the COPDGene study. Results indicate that the authors' algorithm performs comparably to pulmonologist-generated lung lobe segmentations and can produce good results in cases with accessory fissures, incomplete fissures, advanced emphysema, and low dose acquisition protocols. Dice scores indicate that only 29 out of 500 (5.85%) lobes showed Dice scores lower than 0.9. Two different approaches for evaluating lobe boundary surface discrepancies were applied and indicate that algorithm boundary identification is most accurate in the vicinity of fissures detectable on CT. Conclusions: The proposed algorithm is effective for lung lobe segmentation in absence of auxiliary structures such as vessels and airways. The most challenging cases are those with mostly incomplete, absent, or near-absent fissures and in cases with poorly revealed fissures due to high image noise. However, the authors observe good performance even in the majority of these cases.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

NPTFit: A Code Package for Non-Poissonian Template Fitting

We present NPTFit, an open-source code package, written in Python and Cython, for performing non-Poissonian template fits (NPTFs). The NPTF is a recently developed statistical procedure for characterizing the contribution of unresolved point sources (PSs) to astrophysical data sets. The NPTF was first applied to Fermi gamma-ray data to provide evidence that the excess of ∼GeV gamma-rays observed in the inner regions of the Milky Way likely arises from a population of sub-threshold point sources, and the NPTF has since found additional applications studying sub-threshold extragalactic sources at high Galactic latitudes. The NPTF generalizes traditional astrophysical template fits to allow for the ability to search for populations of unresolved PSs that may follow a given spatial distribution. NPTFit builds upon the framework of the fluctuation analyses developed in X-ray astronomy, thus it likely has applications beyond those demonstrated with gamma-ray data. The NPTFit package utilizes novel computational methods to perform the NPTF efficiently. The code is available at http://github.com/bsafdi/NPTFit and up-to-date and extensive documentation may be found at http://nptfit.readthedocs.io.

79 ASTRONOMY AND ASTROPHYSICS↗

GADRAS-DRF Enhancements for Safeguards – Custom Peak Fitting to Enhance Model Fitting and Isotopics

One major software update was accomplished within the Gamma Detector Response and Analysis Software-Detector Response Function (GADRAS-DRF) package. This update allows users to adjust individual peak fits for use in subsequent analysis processes within GADRAS-DRF. A graphical user interface (GUI) was implemented so users can see the effect of their fit adjustments on the spectrum. This new feature will enhance the capability of the previously funded auto-enrichment and peak-based 1D model fit feature that was implemented in FY23 (Fiscal Year). Isotopic analysis was performed using M400 data obtained from uranium standards and the isotopic assessment is given with and without manual peak adjustments.

07 ISOTOPE AND RADIATION SOURCES↗

Derivation of general analytic gradient expressions for density-fitted post-Hartree-Fock methods: An efficient implementation for the density-fitted second-order Møller–Plesset perturbation theory

General analytic gradient expressions (with the frozen-core approximation) are presented for density-fitted post-HF methods. An efficient implementation of frozen-core analytic gradients for the second-order Møller–Plesset perturbation theory (MP2) with the density-fitting (DF) approximation (applying to both reference and correlation energies), which is denoted as DF-MP2, is reported. The DF-MP2 method is applied to a set of alkanes, conjugated dienes, and noncovalent interaction complexes to compare the computational cost of single point analytic gradients with MP2 with the resolution of the identity approach (RI-MP2) [F. Weigend and M. Häser, Theor. Chem. Acc. 97, 331 (1997); R. A. Distasio, R. P. Steele, Y. M. Rhee, Y. Shao, and M. Head-Gordon, J. Comput. Chem. 28, 839 (2007)]. In the RI-MP2 method, the DF approach is used only for the correlation energy. Our results demonstrate that the DF-MP2 method substantially accelerate the RI-MP2 method for analytic gradient computations due to the reduced input/output (I/O) time. Because in the DF-MP2 method the DF approach is used for both reference and correlation energies, the storage of 4-index electron repulsion integrals (ERIs) are avoided, 3-index ERI tensors are employed instead. Further, as in case of integrals, our gradient equation is completely avoid construction or storage of the 4-index two-particle density matrix (TPDM), instead we use 2- and 3-index TPDMs. Hence, the I/O bottleneck of a gradient computation is significantly overcome. Therefore, the cost of the generalized-Fock matrix (GFM), TPDM, solution of Z-vector equations, the back transformation of TPDM, and integral derivatives are substantially reduced when the DF approach is used for the entire energy expression. Further application results show that the DF approach introduce negligible errors for closed-shell reaction energies and equilibrium bond lengths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MICROCALORIMETER SPECTROSCOPY AT HIGH PULSE RATES: A MULTI-PULSE FITTING TECHNIQUE

Transition Edge Sensor microcalorimeters can measure X-ray and gamma-ray energies with very high energy resolution and high photon-collection efficiency. For this technology to reach its full potential in future X-ray observatories, each sensor must be able to measure hundreds or even thousands of photon energies per second. Current “optimal filtering” approaches to achieve the best possible energy resolution work only for photons that are well isolated in time, a requirement which is in direct conflict with the need for high-rate measurements. We describe a new analysis procedure to allow fitting for the pulse height of all photons even in the presence of heavy pulse pile-up. In the limit of isolated pulses, the technique reduces to standard optimal filtering with long records. We employ reasonable approximations to the noise covariance function in order to render this procedure computationally viable even for very long data records. The technique is employed to analyze X-ray emission spectra at 600 eV and 6 keV at rates up to 250 counts s{sup −1} in microcalorimeters having exponential signal decay times of approximately 1.2 ms.

79 ASTRONOMY AND ASTROPHYSICS↗