Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical methods”

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

Breaking the mold: Overcoming the time constraints of molecular dynamics on general-purpose hardware

The evolution of molecular dynamics (MD) simulations has been intimately linked to that of computing hardware. For decades following the creation of MD, simulations have improved with computing power along the three principal dimensions of accuracy, atom count (spatial scale), and duration (temporal scale). Since the mid-2000s, computer platforms have, however, failed to provide strong scaling for MD, as scale-out central processing unit (CPU) and graphics processing unit (GPU) platforms that provide substantial increases to spatial scale do not lead to proportional increases in temporal scale. Important scientific problems therefore remained inaccessible to direct simulation, prompting the development of increasingly sophisticated algorithms that present significant complexity, accuracy, and efficiency challenges. While bespoke MD-only hardware solutions have provided a path to longer timescales for specific physical systems, their impact on the broader community has been mitigated by their limited adaptability to new methods and potentials. In this work, we show that a novel computing architecture, the Cerebras wafer scale engine, completely alters the scaling path by delivering unprecedentedly high simulation rates up to 1.144 M steps/s for 200 000 atoms whose interactions are described by an embedded atom method potential. This enables direct simulations of the evolution of materials using general-purpose programmable hardware over millisecond timescales, dramatically increasing the space of direct MD simulations that can be carried out. In this paper, we provide an overview of advances in MD over the last 60 years and present our recent result in the context of historical MD performance trends.

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

Efficiently evaluating loop integrals in the EFTofLSS using QFT integrals with massive propagators

We develop a new way to analytically calculate loop integrals in the Effective Field Theory of Large Scale-Structure. Previous available methods show severe limitations beyond the one-loop power spectrum due to analytical challenges and computational and memory costs. Our new method is based on fitting the linear power spectrum with cosmology-independent functions that resemble integer powers of quantum field theory massive propagators with complex masses. A remarkably small number of them is sufficient to reach enough accuracy. Similarly to former approaches, the cosmology dependence is encoded in the coordinate vector of the expansion of the linear power spectrum in our basis. We first produce cosmology-independent tensors where each entry is the loop integral evaluated on a given combination of basis vectors. For each cosmology, the evaluation of a loop integral amounts to contracting this tensor with the coordinate vector of the linear power spectrum. The 3-dimensional loop integrals for our basis functions can be evaluated using techniques familiar to particle physics, such as recursion relations and Feynman parametrization. We apply our formalism to evaluate the one-loop bispectrum of galaxies in redshift space. The final analytical expressions are quite simple and can be evaluated with little computational and memory cost. We show that the same expressions resolve the integration of all one-loop N-point function in the EFTofLSS. This method, which is originally presented here, has already been applied in the first one-loop bispectrum analysis of the BOSS data to constraint ΛCDM parameters and primordial non-Gaussianities [1, 2].

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Enhancing the accuracy of XPS calculations: Exploring hybrid basis set schemes for CVS-EOMIP-CCSD calculations

Reliable computational methodologies and basis sets for modeling x-ray spectra are essential for extracting and interpreting electronic and structural information from experimental x-ray spectra. In particular, the trade-off between numerical accuracy and computational cost due to the size of the basis set is a major challenge, since molecular orbitals undergo extreme relaxation in the core-hole state. To gain clarity on the changes in electronic structure induced by the formation of a core-hole, the use of sufficiently flexible basis for expanding the orbitals, particularly for the core region, has been shown to be essential. This work focuses on the refinement of core-hole ionized state calculations using the equation-of-motion coupled cluster family of methods through an extensive analysis on the effectiveness of “hybrid” and mixed basis sets. In this investigation, we utilize the CVS-EOMIP-CCSD method in combination and construct hybrid basis sets piecewise from readily available Dunning’s correlation consistent basis sets in order to calculate x-ray ionization energies (IEs) for a set of small gas phase molecules. Our results provide insights into the impact of basis sets on the CVS-EOMIP-CCSD calculations of K-edge IEs of first-row p-block elements. Furthermore, these insights enable us to understand more about the basis set dependence of the core IEs computed and allow us to establish a protocol for deriving reliable and cost-effective theoretical estimates for computing IEs of small molecules containing such elements.

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

Conformal maps of viscous electron flow in the Gurzhi crossover

Herein we investigate the impact of geometric constriction on the viscous flow of electron liquid through quantum point contacts. We provide analysis on the electric potential distribution given the setup of a slit configuration and use the method of conformal mapping to obtain analytical results. The potential profile can be tested and contrasted experimentally with the scanning tunneling potentiometry technique. We discuss intricate physics that underlies the Gurzhi effect, i.e., the enhancement of conductivity in the viscous flow, and compare results for different boundary conditions. In addition, we calculate the temperature dependence of the momentum relaxation time as a result of impurity assisted quasiballistic interference effects and discuss various correlational corrections that lead to the violation of Matthiessen's rule in the hydrodynamic regime. We caution that spatially inhomogeneous profiles of current in the Gurzhi crossover between Ohmic and Stokes flows might also appear in the nonhydrodynamic limit where nonlocality plays an important role. This conclusion is corroborated by calculation of dispersive conductivity in the weakly impure limit.

2-dimensional systems↗
DOE OSTI2020

Toward a systematic improvement of the fixed-node approximation in diffusion Monte Carlo for solids—A case study in diamond

While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for ab initio electronic structure calculations, in practice, the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions with approximate nodes (or zeros). This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced. This approach is illustrated in the case of the carbon diamond using Slater–Jastrow trial wavefunctions including up to one million Slater determinants. Fixed-node DMC energies obtained with such large expansions are much improved, and the fixed-node error is found to decrease monotonically and smoothly as a function of the number of determinants in the trial wavefunction, a property opening the way to a better control of this error. The cohesive energy extrapolated to the thermodynamic limit is in close agreement with the estimated experimental value. Interestingly, this is also the case at the single-determinant level, thus, indicating a very good error cancellation in carbon diamond between the bulk and atomic total fixed-node energies when using single-determinant nodes.

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

Integrated analysis of X-ray diffraction patterns and pair distribution functions for machine-learned phase identification

Abstract To bolster the accuracy of existing methods for automated phase identification from X-ray diffraction (XRD) patterns, we introduce a machine learning approach that uses a dual representation whereby XRD patterns are augmented with simulated pair distribution functions (PDFs). A convolutional neural network is trained directly on XRD patterns calculated using physics-informed data augmentation, which accounts for experimental artifacts such as lattice strain and crystallographic texture. A second network is trained on PDFs generated via Fourier transform of the augmented XRD patterns. At inference, these networks classify unknown samples by aggregating their predictions in a confidence-weighted sum. We show that such an integrated approach to phase identification provides enhanced accuracy by leveraging the benefits of each model’s input representation. Whereas networks trained on XRD patterns provide a reciprocal space representation and can effectively distinguish large diffraction peaks in multi-phase samples, networks trained on PDFs provide a real space representation and perform better when peaks with low intensity become important. These findings underscore the importance of using diverse input representations for machine learning models in materials science and point to new avenues for automating multi-modal characterization.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Crystal Growth of Quaternary RE 2 EuSi 2 S 8 ( RE = Ce–Nd, Sm, Gd, Tb) Using Flux-Assisted Boron Chalcogen Mixture (BCM) Method: Investigation of Magnetic and Luminescence Properties

A series of quaternary rare-earth containing thiosilicates with the general formula RE 2 EuSi 2 S 8 (RE = Ce–Nd, Sm, Gd, Tb) has been synthesized via the flux-assisted boron chalcogen mixture (BCM) crystal growth method. High-quality single crystals were obtained, and their crystal structures were determined by single-crystal X-ray diffraction. The RE 2 EuSi 2 S 8 series crystallizes in the trigonal system, adopting the space group R-3c. Polycrystalline samples were employed for physical property measurements, including magnetic susceptibility measurements, UV–visible diffuse reflectance, and photoluminescent response. Magnetic data of RE 2 EuSi 2 S 8 ( RE = Ce, Nd, and Gd) were collected over the 2–300 K temperature range. The samples were paramagnetic behavior with negative Weiss constants (θ W = −11.05, −10.55, and −1.35 K respectively). Their thermal stability was investigated using thermogravimetric analysis (TGA). Optical band gaps, estimated from diffuse reflectance spectra, were determined to be 2.2(1) eV for Ce 2 EuSi 2 S 8 , 1.8(1) eV for Nd 2 EuSi 2 S 8 , and 1.7(1) eV for Gd 2 EuSi 2 S 8 respectively. Finally, photoluminescence measurements were collected on Ce 2 EuSi 2 S 8 and Tb 2 EuSi 2 S 8 single crystals.

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

Structural Investigation of Six Quinary Sulfides Synthesized via the Flux-Assisted Boron-Chalcogen Mixture (BCM) Method: Eu 2+ Containing Members of the RE 3 MTQ 7 (M and T = Transition or Main Group Metals, Q = Chalcogens) Family

For this work, a series of six quinary rare-earth sulfides Ce 4+ 1.85 Eu 2+ 1.15 Na 0.30 SiS 7 , Ce 4+ 1.91 Eu 2+ 1.09 K 0.18 SiS 7 , Ce 4+ 1.96 Eu 2+ 1.04 Rb 0.08 SiS 7 , Ce 4+ 1.98 Eu 2+ 1.02 Cs 0.05 SiS 7 , Ce 4+ 1.97 Eu 2+ 1.03 Ag 0.06 SiS 7 , and Ce 4+ 1.50 Eu 2+ 1.50 CuSiS 7 were obtained in an alkali iodide flux using the boron-chalcogen mixture (BCM) method. Single crystal X-ray diffraction was used to determine the structures of the high quality single crystals that were grown; their elemental compositions were confirmed by energy-dispersive spectroscopy (EDS). The compounds crystallize in the hexagonal crystal system in the noncentrosymmetric space group P63. The crystal structure consists of a three-dimensional network composed of mixed cerium and europium bicapped trigonal prisms, isolated SiS4 tetrahedra, and monovalent metals (Na, K, Rb, Cs, Ag, and Cu) located in cavities created by linked Ce/EuS 8 polyhedra. The structures are charge-balanced when Ce and Eu are in their +4 and +2 oxidation states, respectively. The effective magnetic moment of Ce 1.50 4+ Eu 1.50 2+ CuSiS 7 determined from the temperature dependence of the magnetic susceptibility data is consistent with the presence of Ce 4+ and Eu 2+ . Clear correlations between the alkali ion site occupancy, the ionic radius of the alkali cations, and the average bond length of Ce 4+ /Eu 2+ –S, were established. UV–vis diffuse reflectance data were collected for Ce 1.50 4+ Eu 1.50 2+ CuSiS 7 and a band gap of 1.9(1) eV was established.

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

Crystal Growth of Quaternary Rare Earth Selenosilicates by Using the Flux-Assisted Boron Chalcogen Mixture Method: Investigation of Their Magnetic and Optical Properties

A series of rare earth magnesium selenosilicates, RE3Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd, Tb, Dy) were obtained as single crystal using the flux assisted boron chalcogen mixture (BCM) method. The structures of the crystals were determined by single-crystal X-ray diffraction. The RE3Mg 0.5 SiSe 7 series crystallizes in the hexagonal crystal system in the space group P6 3 . Polycrystalline powders were synthesized to perform physical property measurements. Magnetic measurements over the 2–300 K temperature range reveal that Ce 3 Mg 0.5 SiSe 7 and Gd 3 Mg 0.5 SiSe 7 exhibit paramagnetic behavior with negative Weiss constants (θ W = −14.50, θ W = −6.13 K, respectively). The optical properties of RE 3 Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd) were measured by ultraviolet–visible (UV–vis) diffuse reflectance. Density functional theory (DFT) electronic structure calculations were performed. Furthermore, a second harmonic generation measurement was performed on a polycrystalline powder of Ce 3 Mg 0.5 SiSe 7 and was found to be SHG active with an efficiency of 0.11 times the standard potassium dihydrogen phosphate (KDP).

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

Synthesis of a series of rare-earth-based multi-anion chalcogenide iodides RE 3 Si 2 Se x S 8− x I (RE = La, Ce, Pr, and Nd) using the flux-assisted boron–chalcogen mixture method

Single crystals and polycrystalline powders of rare earth mixed chalcogenide iodides La 3 Si 2 Se 1.21 S 6.79 I, Ce 3 Si 2 Se 1.39 S 6.61 I, Pr 3 Si 2 Se 1.22 S 6.78 I, and Nd 3 Si 2 Se 1.18 S 6.82 I were prepared using the reactive flux-assisted boron–chalcogen mixture (BCM) method at 850 °C. All compounds crystallized in the monoclinic crystal system, space group C2/c (space group number 15). The series adopts the La 3 Si 2 O 8 Cl structure type, containing isolated SiQ 4 tetrahedra connected by REQ 8 (RE = La, Ce, Pr and Nd) polyhedra; this arrangement creates tunnels that are filled by I atoms. The partial substitution of S by Se was carried out to modulate the optical properties. Phase pure samples and uniform solid solutions were obtained for all compositions as determined using powder X-ray diffraction patterns. Polycrystalline powders were used for physical property measurements, including magnetic susceptibility and UV-Vis diffuse reflectance. The solid-state UV-Vis data for the polycrystalline La 3 Si 2 Se 1.21 S 6.79 I, Ce 3 Si 2 Se 1.39 S 6.61 I, and Pr 3 Si 2 Se 1.22 S 6.78 I samples revealed band gaps of E g = 2.5(1), 2.2(1), and 2.3(1) eV, typical of semiconductors. Magnetic measurements indicated that Ce 3 Si 2 Se 1.39 S 6.61 I and Nd 3 Si 2 Se 1.18 S 6.82 I exhibit paramagnetic behavior with slightly negative Weiss constants θ = −25 and −38. The photoluminescence spectrum of Ce 3 Si 2 Se 1.39 S 6.61 I exhibits a broad emission band around ∼493 nm.

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

Spatially quasi-periodic bifurcations from periodic traveling water waves and a method for detecting bifurcations using signed singular values

We present a method of detecting bifurcations by locating zeros of a signed version of the smallest singular value of the Jacobian. This enables the use of quadratically convergent root-bracketing techniques or Chebyshev interpolation to locate bifurcation points. Only positive singular values have to be computed, though the method relies on the existence of an analytic or smooth singular value decomposition (SVD). The sign of the determinant of the Jacobian, computed as part of the bidiagonal reduction in the SVD algorithm, eliminates slope discontinuities at the zeros of the smallest singular value. We use the method to search for spatially quasi-periodic traveling water waves that bifurcate from large-amplitude periodic waves. The water wave equations are formulated in a conformal mapping framework to facilitate the computation of the quasi-periodic Dirichlet-Neumann operator. We find examples of pure gravity waves with zero surface tension and overhanging gravity-capillary waves. In both cases, the waves have two spatial quasi-periods whose ratio is irrational. We follow the secondary branches via numerical continuation beyond the realm of linearization about solutions on the primary branch to obtain traveling water waves that extend over the real line with no two crests or troughs of exactly the same shape. The pure gravity wave problem is of relevance to ocean waves, where capillary effects can be neglected. Such waves can only exist through secondary bifurcation as they do not persist to zero amplitude. The gravity-capillary wave problem demonstrates the effectiveness of using the signed smallest singular value as a test function for multi-parameter bifurcation problems. This test function becomes mesh independent once the mesh is fine enough.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enhanced liquid metal wetting on oxide surfaces via patterned particles

Liquid metals typically do not intrinsically wet oxide surfaces due to chemical dissimilarity, preventing the formation of well-bonded interfaces in many applications. This paper investigates how the wetting and spreading of liquid silver (Ag) droplets on oxide surfaces can be enhanced via patterned nickel (Ni) particles using multiscale modeling. First, a force field (FF) for an Ag/yttria-stabilized-zirconia system was developed by fitting the interfacial binding energy to density-functional theory (DFT) calculations. Then, molecular dynamics (MD) simulations using these FF parameters revealed that the good wetting of Ag on Ni caused Ag droplets to be attracted to the Ni particles, promoting spreading. Thus, the optimal hexagonal Ni pattern was shown to enhance the wetting area by 224% compared to a bare oxide surface. Second, an analytical model predicting the maximum possible enhancement in wetting area as a function of the intrinsic wetting angle, size of the Ag droplet, and size of the Ni particles was derived and validated with MD simulations, to generalize this method to other substrates. Here, both the analytical model and MD results suggested that the ability of Ni patterns to enhance the Ag wetting area becomes more significant as the wettability of substrate becomes worse.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Minimum Detectable Intakes and Doses for Uranium Bioassays—Comparison between Alpha Spectrometry and ICP-MS

Naturally occurring uranium complicates monitoring for occupational exposures. There are several retroactive methods that can be used to monitor for occupational exposures, with benefits and drawbacks to each. Analysis of uranium in urine by mass spectrometry and alpha spectrometry is compared, and methods of determining an occupational exposure are presented. Furthermore, the minimum detectable concentrations from each analysis and a method for intake determination based on the analytical results are compared for various solubility types and mixtures. Mass spectrometry with radiochemical separation was found to be the most sensitive analysis for detecting occupational exposures to anthropogenic mixtures based on minimum detectable doses calculated from the proposed method for intake determination.

bioassay↗
DOE OSTI2026

A critical review on additive manufacturing of refractory alloys from a data analytics perspective- beyond nickel-based superalloys

Refractory alloys (RAs) are promising materials due to their exceptional physicochemical properties, but most research remains at the laboratory scale. For broader adoption, advancements in manufacturing are essential. Because their high stability makes conventional methods like machining and casting difficult, additive manufacturing (AM) is emerging as an effective approach for fabricating refractory alloy components. However, AM's repeated non-equilibrium thermal cycles introduce undesired features (e.g. defects, anisotropic microstructures, and residual stresses), which are magnified due to RAs’ unique properties. This paper comprehensively reviews the state-of-the-art methods of AM for refractory alloys. It explores data analytics techniques to establish design rules based on multi-fidelity experimental and computational methods. Furthermore, it investigates integrated, collaborative efforts to harmonise standalone databases, information, knowledge, and predictive models at multi-physics, multi-stage, and multi-scale. Unlike the existing literature that focuses primarily on material systems or process fundamentals, this work provides an integrated perspective on AM of refractory alloys from a data analytics standpoint, highlighting the roles of integrated computational materials engineering (ICME), verification, validation, and uncertainty quantification (VV&UQ), and digital twin-driven qualification in overcoming data scarcity and accelerating rapid qualification.

Additive manufacturing↗
DOE OSTI2020

The Lorentzian inversion formula and the spectrum of the 3d O(2) CFT

We study the spectrum and OPE coefficients of the three-dimensional critical O(2) model, using four-point functions of the leading scalars with charges 0, 1, and 2 ( s , $\phi$, and t ). We obtain numerical predictions for low-twist OPE data in several charge sectors using the extremal functional method. We compare the results to analytical estimates using the Lorentzian inversion formula and a small amount of numerical input. We find agreement between the analytic and numerical predictions. We also give evidence that certain scalar operators lie on double-twist Regge trajectories and obtain estimates for the leading Regge intercepts of the O(2) model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Gravitational production of heavy particles during and after inflation

Abstract We investigate the gravitational production of a scalar fieldχwith a mass exceeding the Hubble scale during inflationm χ ≳H I , employing both analytical and numerical approaches. We demonstrate that the steepest descent method effectively captures the epochs and yields of gravitational production in a compact and simple analytical framework. These analytical results align with the numerical solutions of the field equation. Our study covers three spacetime backgrounds: de Sitter, power-law inflation, and the Starobinsky inflation model. Within these models, we identify two distinct phases of particle production: during and after inflation. During inflation, we derive an accurate analytic expression for the particle production rate, accounting for a varying Hubble rate. After inflation, the additional burst of particle production depends on the inflaton mass around its minimum. When this mass is smaller than the Hubble scale during inflation,H I , there is no significant extra production. However, if the inflaton mass is larger, post-inflation production becomes the dominant contribution. Furthermore, we explore the implications of gravitationally produced heavy fields for dark matter abundance, assuming their cosmological stability.

Physics↗
DOE OSTI2021

Quantitation of amiodarone and N-desethylamiodarone in single HepG2 cells by single-cell printing-liquid vortex capture-mass spectrometry

Quantitative measure of a drug and its associated metabolite(s) with single-cell resolution is often limited by sampling throughput or other compromises that limit broad use. Here, we demonstrate the use of single-cell printing-liquid vortex capture-mass spectrometry (SCP-LVC-MS) to quantitatively measure the intracellular concentrations of amiodarone (AMIO) and its metabolite, N-desethylamiodarone (NDEA), from thousands of single cells across several AMIO incubation concentrations ranging from 0 to 10 μM. Concentrations obtained by SCP-LVC-MS were validated through comparison with average assays and traditional measurement of cells in bulk. Average of SCP-LVC-MS measurements and aggregate vial collection assay the concentrations differed by < 5%. Both AMIO and NDEA had clear log-normal distributions with similar standard deviation of concentrations in the cell population. The mean of both AMIO and NDEA intracellular concentrations were positively correlated with AMIO incubation concentration, increasing from 0.026 to 0.520 and 0.0055 to 0.048 mM for AMIO and NDEA, respectively. The standard deviation of AMIO and NDEA log-normal distribution fits were relatively similar in value across incubation concentrations, 0.15–0.19 log 10 (mM), and exhibited a linear trend with respect to each other. The single cell-resolved conversion ratio of AMIO to NDEA increased with decreasing incubation concentration, 7 ± 2%, 18 ± 3%, and 20 ± 7% for 10.0, 1.0, and 0.1 μM AMIO incubation concentrations, respectively. Finally, association with simultaneously measured lipids had several ions with statistically significant difference in intensity but no clear correlations with AMIO intracellular content was observed.

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

Standard operating procedures for the comprehensive and reliable analysis of cannabis terpenes

Terpenes are the primary determinants of cannabis flower aroma, and ongoing research tests their potential for impacting the overall experience. Frustratingly, despite the importance of terpenes in cannabis physiology and commercial uses, literature reports vary widely regarding the major constituents of volatile blends and the concentrations of individual terpenes. Here, in this article, we provide detailed descriptions of complementary approaches that will allow researchers to determine the identity and quantity of cannabis terpenes unequivocally and reliably. These standard operating procedures will guide decisions about which method to employ to address specific analytical goals. We are including two application examples to illustrate the utility of different approaches for tackling the analysis of terpenes in cannabis flower samples.

Polito, Joshua↗