Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

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 181 records · Page 10

The Convolutional Multiple Whole Profile (CMWP) Fitting Method, a Global Optimization Procedure for Microstructure Determination

The analysis of line broadening in X-ray and neutron diffraction patterns using profile functions constructed on the basis of well-established physical principles and TEM observations of lattice defects has proven to be a powerful tool for characterizing microstructures in crystalline materials. These principles are applied in the convolutional multiple-whole-profile (CMWP) procedure to determine dislocation densities, crystallite size, stacking fault and twin boundary densities, and intergranular strains. The different lattice defect contributions to line broadening are separated by considering the hkl dependence of strain anisotropy, planar defect broadening and peak shifts, and the defect dependent profile shapes. The Levenberg–Marquardt (LM) peak fitting procedure can be used successfully to determine crystal defect types and densities as long as the diffraction patterns are relatively simple. However, in more complicated cases like hexagonal materials or multiple-phase patterns, using the LM procedure alone may cause uncertainties. Here, we extended the CMWP procedure by including a Monte Carlo statistical method where the LM and a Monte Carlo algorithm were combined in an alternating manner. The updated CMWP procedure eliminated uncertainties and provided global optimized parameters of the microstructure in good correlation with electron microscopy methods.

36 MATERIALS SCIENCE↗

A Principled Framework to Assess the Information-Theoretic Fitness of Brain Functional Sub-Circuits

In systems and network neuroscience, many common practices in brain connectomic analysis are often not properly scrutinized. One such practice is mapping a predetermined set of sub-circuits, like functional networks (FNs), onto subjects’ functional connectomes (FCs) without adequately assessing the information-theoretic appropriateness of the partition. Another practice that goes unchallenged is thresholding weighted FCs to remove spurious connections without justifying the chosen threshold. This paper leverages recent theoretical advances in Stochastic Block Models (SBMs) to formally define and quantify the information-theoretic fitness (e.g., prominence) of a predetermined set of FNs when mapped to individual FCs under different fMRI task conditions. Our framework allows for evaluating any combination of FC granularity, FN partition, and thresholding strategy, thereby optimizing these choices to preserve the important topological features of the human brain connectomes. By applying to the Human Connectome Project with Schaefer parcellations at multiple levels of granularity, the framework showed that the common thresholding value of 0.25 was indeed information-theoretically valid for group-average FCs, despite its previous lack of justification. Our results pave the way for the proper use of FNs and thresholding methods, and provide insights for future research in individualized parcellations.

Duong-Tran, Duy (ORCID:0009000944967575)↗

Methods in PES-Learn: Direct-Fit Machine Learning of Born–Oppenheimer Potential Energy Surfaces

The release of PES-L EARN version 1.0 as an open-source software package for the automatic construction of machine learning models of semi-global molecular potential energy surfaces (PESs) is presented. Improvements to PES-L EARN ’s interoperability are stressed with new Python API that simplifies workflows for PES construction via interaction with QCSchema input and output infrastructure. In addition, a new machine learning method is introduced to PES-L EARN : kernel ridge regression (KRR). The capabilities of KRR are emphasized with examination of select semi-global PESs. All machine learning methods available in PES-L EARN are benchmarked with benzene and ethanol datasets from the rMD17 database to illustrate PES-L EARN ’s performance ability. Fitting performance and timings are assessed for both systems. Finally, the ability to predict gradients with neural network models is presented and benchmarked with ethanol and benzene. PES-L EARN is an active project and welcomes community suggestions and contributions.

kernel ridge regression↗

Precise Dynamical Masses and Orbital Fits for β Pic b and β Pic c

We present a comprehensive orbital analysis to the exoplanets β Pictoris b and c that resolves previously reported tensions between the dynamical and evolutionary mass constraints on β Pic b. We use the Markov Chain Monte Carlo orbit code orvara to fit 15 years of radial velocities and relative astrometry (including recent GRAVITY measurements), absolute astrometry from Hipparcos and Gaia, and a single relative radial velocity measurement between β Pic A and b. We measure model-independent masses of 9.3{sub −2.5}{sup +2.6} M {sub Jup} for β Pic b and 8.3 ± 1.0 M {sub Jup} for β Pic c. These masses are robust to modest changes to the input data selection. We find a well-constrained eccentricity of 0.119 ± 0.008 for β Pic b, and an eccentricity of 0.21{sub −0.09}{sup +0.16} for β Pic c, with the two orbital planes aligned to within ∼05. Both planets’ masses are within ∼1σ of the predictions of hot-start evolutionary models and exclude cold starts. We validate our approach on N-body synthetic data integrated using REBOUND. We show that orvara can account for three-body effects in the β Pic system down to a level ∼5 times smaller than the GRAVITY uncertainties. Systematics in the masses and orbital parameters from orvara’s approximate treatment of multiplanet orbits are a factor of ∼5 smaller than the uncertainties we derive here. Future GRAVITY observations will improve the constraints on β Pic c’s mass and (especially) eccentricity, but improved constraints on the mass of β Pic b will likely require years of additional radial velocity monitoring and improved precision from future Gaia data releases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dissecting Nearby Galaxies with piXedfit. I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

We study spatially resolved properties (on spatial scales of ~1–2 kpc out to at least 3 effective radii) of the stars, dust, and gas in 10 nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point-spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy-balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX–βrelation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of the star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kiloparsec scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), and gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of the molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center.

79 ASTRONOMY AND ASTROPHYSICS↗

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

Particle transport constraints via Bayesian spectral fitting of multiple atomic lines

Optimized operation of fusion devices demands detailed understanding of plasma transport, a problem that must be addressed with advances in both measurement and data analysis techniques. In this work, we adopt Bayesian inference methods to determine experimental particle transport, leveraging opportunities from high-resolution He-like ion spectra in a tokamak plasma. The Bayesian spectral fitting code is used to analyze resonance (w), forbidden (z), intercombination (x, y), and satellite (k, j) lines of He-like Ca following laser blow-off injections on Alcator C-Mod. This offers powerful transport constraints since these lines depend differently on electron temperature and density, but also differ in their relation to Li-like, He-like, and H-like ion densities, often the dominant Ca charge states over most of the C-Mod plasma radius. Using synthetic diagnostics based on the AURORA package, we demonstrate improved effectiveness of impurity transport inferences when spectroscopic data from a progressively larger number of lines are included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSUF University Research Reactor Fitness Workshop

University research reactors have been a cornerstone of nuclear engineering research and education since the first reactor was deployed at North Carolina State University in the 1950s. The population of university reactors grew to a high of almost 80 in 1970 but has dropped to 25 operating today, also shrinking the community, availability of vendors and the talent pipeline. The nation’s fleet of TRTRs serves as a vital resource for educating the next generation of nuclear engineers and as a flexible testbed for innovative reactor technologies. Typically, the reactors have a very small staff with responsibilities ranging through research and operations, isotope production, regulatory compliance, facility maintenance, and university teaching obligations. Knowledge transfer and staffing were identified as critical areas in a recent study of US university research reactor fitness. This panel discusses the unique challenges and potential solutions to these issues.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Composite to metal end fitting joint

A composite torque pipe to metal fitting joint including concentric inner and outer shells formed with respective smooth conical surfaces projecting distally and respectively tapered radially inwardly and radially outwardly to form an annulus expanding in thickness for receipt of a tapered connector ring complementally shaped at one end of a composite pipe and bonded along the interface. The present invention includes selecting the sleeves with the corresponding shells and, forming the torque pipe with the tapered ring and bonding the tapered ring in the annulus.

Leslie, James C.↗

Fast Automatic Knot Placement Method for Accurate B-spline Curve Fitting

The choice of knot vector has immense influence on the resulting accuracy of a B-spline approximation of a curve. However, despite the significance of this problem and the various solutions that were proposed in the literature, optimizing the number and placement of knots remains a difficult task. This paper presents a novel method for the approximation of a curve by a B-spline of arbitrary order, which automatically determines a knot vector that achieves high approximation quality. Additionally, at the core of our approach is a feature function that characterizes the amount and spatial distribution of geometric details in the input curve by estimating its derivatives. Knots are then selected in such a way as to evenly distribute the feature contents across their intervals. A comparison to the state of the art for a wide variety of curves shows that our method is faster and achieves more accurate reconstruction results, while typically reducing the number of necessary knots.

97 MATHEMATICS AND COMPUTING↗

Tracker Fits for shot 1s1674

Gas gun shock-to-detonation experiments provide data for calibrating HE burn models. At LANL the principal diagnostic is a magnetic gauge package that provides data on the lead shock trajectory, which is used to determine the run-distance to detonation, and the lead shock speed with run distance, which together with jump-off velocity of the velocity gauges determines Hugoniot locus points for the reactants.

42 ENGINEERING↗

A Diazo Linker Ligand Promotes Flexibility and Induced Fit Binding in a Microporous Copper Coordination Network

Abstract Flexible organic linkers represent an intuitive and effective strategy to design flexible metal–organic materials. We report herein a systematic study concerning the effect of varying the central bond of mixed pyridyl‐benzoate linkers, L, upon the flexibility of three isostructuralkddtopology microporous coordination networks (CNs) of formula ML 2 :X‐kdd‐1‐Cu,1= L = (E)‐4‐(pyridin‐4‐yldiazenyl)benzoate;X‐kdd‐2‐Cu,2= L = (E)‐4‐(2‐(pyridin‐4‐yl)vinyl)benzoate; the previously reportedX‐kdd‐3‐Cu,3= L = 4‐(pyridin‐4‐ylethynyl)benzoate. As revealed by single crystal x‐ray diffraction (SCXRD) and gas sorption studies,X‐kdd‐1‐Cu, exhibited gate‐opening during CO 2 and hydrocarbon (C2 and C8) sorption experiments whereas the other two CNs did not. Insight into these phase transformations was gained from in situ variable‐pressure and variable temperature powder X‐ray diffraction (PXRD), SCXRD, and modeling. Rotation of ligand1around the diazo bond, torsion angle changes between phenyl and carboxylate moieties, and deformation of the Cu‐based rod building blocks enabled activatedX‐kdd‐1‐Cuto form new phases with C8 isomers and CH 2 Cl 2 , CH 2 Cl 2 inducing contraction of the activated phase. Computational studies suggest that1enables flexibility thanks to its lower barrier of deformation versus2or3. This study teaches that diazo moieties could offer a general strategy to enhance the flexibility of CNs.

Chemistry↗

Flavor dependence of unpolarized quark transverse momentum distributions from a global fit

We present an extraction of the unpolarized transverse-momentum-dependent parton distribution and fragmentation functions that takes into account possible differences between quark flavors and final-state hadrons. The extraction is based on experimental measurements from Drell-Yan processes and semi-inclusive deep-inelastic scattering, whose combination is essential to distinguish flavor differences. The analysis is carried out at N 3 LL accuracy. The extracted flavor-dependent distributions give a very good description of the data (χ 2 /N dat = 1.08). The resulting error bands take fully into account also the uncertainties in the determination of the corresponding collinear distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Unpolarized transverse momentum distributions from a global fit of Drell-Yan and semi-inclusive deep-inelastic scattering data

We present an extraction of unpolarized transverse-momentum-dependent parton distribution and fragmentation functions based on more than two thousand data points from several experiments for two different processes: semi-inclusive deep-inelastic scattering and Drell-Yan production. The baseline analysis is performed using the Monte Carlo replica method and resumming large logarithms at N 3 LL accuracy. The resulting description of the data is very good (χ 2 /N dat = 1.06). For semi-inclusive deep-inelastic scattering, predictions for multiplicities are normalized by factors that cure the discrepancy with data introduced by higher-order perturbative corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fitting Matérn smoothness parameters using automatic differentiation

The Mat$\acute{e}$rn covariance function is ubiquitous in the application of Gaussian processes to spatial statistics and beyond. Perhaps the most important reason for this is that the smoothness parameter $\nu$ gives complete control over the mean-square differentiability of the process, which has significant implications for the behavior of estimated quantities such as interpolants and forecasts. Unfortunately, derivatives of the Mat$\acute{e}$rn covariance function with respect to $\nu$ require derivatives of the modified second-kind Bessel function $K$ $\nu$ with respect to $\nu$. While closed form expressions of these derivatives do exist, they are prohibitively difficult and expensive to compute. For this reason, many software packages require fixing $\nu$ as opposed to estimating it, and all existing software packages that attempt to offer the functionality of estimating $\nu$ use finite difference estimates for $\partial$ $\nu$ $K$ $\nu$ . In this work, we introduce a new implementation of $K$$\nu$ that has been designed to provide derivatives via automatic differentiation (AD), and whose resulting derivatives are significantly faster and more accurate than those computed using finite differences. Here, we provide comprehensive testing for both speed and accuracy and show that our AD solution can be used to build accurate Hessian matrices for second-order maximum likelihood estimation in settings where Hessians built with finite difference approximations completely fail.

97 MATHEMATICS AND COMPUTING↗

Combining PCA and nonlinear fitting of peak models to re-evaluate C 1s XPS spectrum of cellulose

Cellulose is an example of a material that responds to XPS by the creation of new chemistry not present in the as-received sample. While improvements in instrumentation may be seen in general as beneficial to surface science, recent studies have shown that the consequences for some materials are detrimental. Here, in this work, these problems are illustrated through an analysis of cellulose spectra obtained during a degradation study. C 1s spectra are decomposed into two well-formed component curves that are open to chemical interpretation. In particular, a component-curve representative of pure cellulose is obtained as well as a second component curve that implies cellulose is degraded through the creation of carbon chemistry involving C—O, C$=$O and O—C$=$O. Since cellulose is a crystalline material, formed through the alignment of molecules under the influence of hydrogen bonds, the analysis and findings presented in this paper are relevant to any material analyzed by XPS whose properties are dependent on hydrogen bonds. The analysis techniques are based on an informed vectorial approach, which extracts directly from data spectral shapes that are used to monitor sample degradation via linear least squares optimization. Related mathematics of Principal Component Analysis and linear analysis are presented.

36 MATERIALS SCIENCE↗