SEARCH · Engineering Papers
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.
Analysis of multimodal characterization of Natrophosphate suggests that the crystalline structure is preserved across a range of synthesis conditions.
Bialkali photocathodes, such as cesium potassium antimonide (CsK 2 Sb), can generate a high-brightness electron beam using a high-power green laser. These photocathode materials have potential applications in advanced accelerators and electron microscopes. It is known that the quantum efficiency (QE) of these photocathodes is affected severely by their substrates; however, reusability of the substrates is not well known. Here, we use graphene, silicon (Si), and molybdenum (Mo) substrates to evaluate the effects of substrates on the QE of redeposited CsK 2 Sb photocathodes after thermal cleanings. We found the QE of CsK 2 Sb photocathodes redeposited on a graphene substrate after a thermal cleaning at 500 °C remained largely unchanged. On the other hand, the QE of redeposited photocathodes on Si and Mo substrates after thermal cleaning at the same temperature decreased drastically. We used X-ray photoelectron spectroscopy (XPS) to quantitatively evaluate the residues of photocathodes after thermal cleaning at 400 °C and 500 °C. We found that the Sb, K, and Cs are removed by thermal cleaning at 500 °C for the graphene substrate, but all or the majority of these elements remained on the Si and Mo substrates. The results were consistent with our density functional theory (DFT) calculations for the case of Si, which we investigated. Furthermore, our angle-resolved photoemission spectroscopy (ARPES) on graphene indicated that its intrinsic electronic structure is preserved after photocathode deposition and thermal cleaning at 500 °C. Hence, we attributed the difference in amount of photocathode residue to the unique dangling-bond- free surface of inert graphene. Our results provide a foundation for graphene-based reusable substrates for high-QE semiconductor photocathodes.
Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.
Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.
Being three-dimensional, stellarators have the advantage that plasma currents are not essential for creating rotational-transform; however, the external current-carrying coils in stellarators can have strong geometrical shaping, which can complicate the construction. Reducing the inter-coil electromagnetic forces acting on strongly shaped 3D coils and the stress on the support structure while preserving the favorable properties of the magnetic field is a design challenge. In this work, we recognize that the inter-coil forces are the gradient of the vacuum magnetic energy. We introduce an objective functional built on the usual quadratic flux on a prescribed target surface together with a weighed penalty on the vacuum energy. The Euler–Lagrange equation for stationary states is derived, and numerical illustrations are computed using a modern stellarator optimization framework. A study of the effect of the energy functional on the inter-coil forces is conducted and the energy is shown to be a promising quantity in producing coils with low forces.
Differential-algebraic equations (DAEs) integrate ordinary differential equations (ODEs) with algebraic constraints, providing a fundamental framework for developing models of dynamical systems characterized by time-scale separation, conservation laws and physical constraints. While sparse optimization has revolutionized model development by allowing data-driven discovery of parsimonious models from a library of possible equations, existing approaches for dynamical systems assume DAEs can be reduced to ODEs by eliminating variables before model discovery. This assumption limits the applicability of such methods for DAE systems with unknown constraints and time scales. We introduce sparse optimization for differential-algebraic systems (SODAs), a data-driven method for the identification of DAEs in their explicit form. By discovering the algebraic and dynamic components sequentially without prior identification of the algebraic variables, this approach leads to a sequence of convex optimization problems. It has the advantage of discovering interpretable models that preserve the structure of the underlying physical system. To this end, SODAs improves since SODAs is singular numerical stability when handling high correlations between library terms, caused by near-perfect algebraic relationships, by iteratively refining the conditioning of the candidate library. We demonstrate the performance of our method on biological, mechanical and electrical systems, showcasing its robustness to noise in both simulated time series and real-time experimental data.
The most important law of radioactivity is that of the exponential decay. In the realm of quantum mechanics, however, this decay law is neither rigorous nor fundamental. The deviations from the exponential decay have been observed experimentally at the early stage of a decay process, but there is little evidence for non-exponential behavior at long times. Yet such long-term non-exponentiality is expected theoretically to probe the non-resonant background components of the initial wave function which preserve the structural interference and the memory of how the state was created. In this paper, we propose new observables that can be used for experimental investigations of the post-exponential decay regime, including the decay of threshold resonances, particle correlations in three-body decays, and interference between near-lying resonances. While the specific examples presented in this work pertain to atomic nuclei, the properties of non-exponential decay are generic, i.e., they apply to other many-body open quantum systems, such as hadrons, atoms, molecules, and nanostructures.
Recent advancements in metal additive manufacturing (AM) show great potential to revolutionize the design and manufacturing of electromagnetic components used in the field of electrical engineering. Lattice structures directly printed by AM processes typically offer better structural performance with reduced weight, such as high stiffness, surface area, elongation, energy absorption, and porosity, than the solid counterpart. This paper aims to study the electromagnetic modeling of lattice structures used in additively manufactured magnetic cores or windings for electric machines. Three dimensional (3D) electromagnetic finite element (FE) analysis with high performance computing (HPC) shows the highest fidelity in predicting the electromagnetic performance of designs with lattice structures by preserving complex geometry details. FE-based homogenization methods have also been explored to potentially speed up concept design. A case study based on an additively manufactured axial-flux permanent magnet machine with a Hilbert pattern stator validates the discussed electromagnetic modeling approaches.
Ti(Si)N coatings were deposited on Si wafers via magnetron sputtering physical vapor deposition technique. Scanning electron microscopy and atomic force microscopy were used to study the morphologies, grain size, and thicknesses of the different films. X-ray diffraction confirmed the formation of the titanium nitride cubic phase. X-ray photoelectron spectroscopy (XPS) was used to study the silicon and oxygen levels in the films, and high resolution XPS was used to determine the surface oxide to nitride ratios as well as the silicon–nitrogen environment. The hardness of these films was then determined with nanoindentation. It was shown that grain size and film density could be controlled by the Si level and the ion bombardment to eliminate the oxygen level in the nitride coatings, which led to the harder films. It was suggested that due to the miscibility of Si, a ternary TiSiN solid solution was formed. Stronger (200) orientation of the TiN B1-NaCl crystals was associated with increased hardness of the coatings. In conclusion, this crystalline structure was preserved after incorporating Si.
This software showcases a latent space dynamics identification method, namely tLaSDI, that embeds the first and second principles of thermodynamics. The latent variables are learned through an autoencoder as a nonlinear dimension reduction model. The latent dynamics are constructed by a neural network-based model that precisely preserves certain structures for the thermodynamic laws through the GENERIC formalism. An abstract error estimate is established, which provides a new loss formulation involving the Jacobian computation of autoencoder. The autoencoder and the latent dynamics are simultaneously trained to minimize the new loss. Computational examples demonstrate the effectiveness of tLaSDI, which exhibits robust generalization ability, even in extrapolation. In addition, an intriguing correlation is empirically observed between a quantity from tLaSDI in the latent space and the behaviors of the full-state solution.
This code implements a Fully Homomorphic Encryption (FHE) system, enabling secure computation on encrypted data without requiring decryption. It supports encryption, decryption, and homomorphic operations like matrix multiplication and addition. This code is adaptable for integrating FHE into linear-time invariant (LTI) systems, including digital control and filtering. With proper configuration from subject matter expertise, encrypted system parameters and signals can be manipulated to perform tasks like state updates, output calculations, and convolution in the encrypted domain. By preserving the structure of LTI systems while ensuring privacy, the framework facilitates secure applications in areas such as autonomous systems, signal processing, and industrial automation. The code initializes the encryption system using parameters provided in the env dictionary. These parameters include the ciphertext modulus, key dimension, plaintext fixed-point scaling factor, and noise bound. During initialization, a secret key is generated, which is essential for encrypting and decrypting data securely. The modular design allows users to tailor these parameters to specific use cases or security requirements. The code implements multiple cryptographic schemes. The learning with errors (LWE) encryption method encodes cleartext message to their plaintext fixed-point representation then encrypted into ciphertext space with additive noise. This noise ensures the security of the scheme, relying on the computational hardness of the LWE problem. The code also includes the Gentry-Sahai-Waters (GSW) scheme based off the LWE problem. Homomorphic matrix multiplication is performed between the LWE and GSW to encrypted data. This is achieved using a decomposition function on the LWE ciphertext during the multiplication operation. For higher-dimensional data, the code includes a method to encrypt entire matrices (GSWMat) using GSW encryption. These encrypted matrices can then be used for homomorphic matrix multiplications (MatMult). The decryption function uses the secret key to recover the original plaintext, removing the added noise and scaling that was originally applied during encryption.
MOOSEenger is an open-source, terminal-first chat application for the MOOSE ecosystem that couples specialized parsing of MOOSE documentation and “.i” input files with retrieval-augmented generation to deliver grounded answers about multiphysics modeling and workflows. It includes dedicated readers for MOOSE-style HTML and a pyhit-based parser that uses the MOOSE syntax tree to preserve block structure and attach retrieval metadata. A data-ingestion pipeline performs semantic chunking into atomic facts and stores them hierarchically in a local Chroma vector database that maintains parent–child relationships across documents; the system can ingest directories, individual files, and single-page web content, and it provides CRUD operations (insert, update, delete) to manage the corpus. At query time, relevant chunks are embedded, retrieved, and fused into the model context, with interactive features such as token streaming, persistent chat history, and dynamic RAG (retrieval triggered by user input or intermediate model output). Deployment is flexible: MOOSEenger runs with local Ollama models or remote Hugging Face/OpenAI backends—typically coordinating generation, lightweight tagging/summarization, and embeddings across three models—and it also supports a server mode and integration with the VS Code Continue interface.
The goal of this IACMI Phase II Technical Collaboration Project was to establish the viability of producing affordable, recycled composite parts, which were produced from fibers with the properties in the range of those reclaimed through the controlled pyrolysis technology demonstrated in Phase I [1]. For the purpose of this report, we use the terms “recovered carbon fiber” (rCF) and “recovered glass fiber” (rGF) as the material reclaimed from controlled pyrolysis. This fiber recycling technology utilizes the inherent energy in composites for fuel and preserves the structural value of glass fiber (GF) and carbon fiber (CF) for reuse. It can also be used for recycling other waste streams at the same facility, thus spreading capital risk and recovery across multiple industries and achieving economy of scale. A second goal to be addressed was the viability of using recycled fibers as reinforcing materials in additive manufacturing applications.
The goal of this IACMI Phase II Technical Collaboration Project was to establish the viability of producing affordable, recycled composite parts, which were produced from fibers with the properties in the range of those reclaimed through the controlled pyrolysis technology demonstrated in Phase I. For the purpose of this report, we use the terms “recovered carbon fiber” (rCF) and “recovered glass fiber” (rGF) as the material reclaimed from controlled pyrolysis. This fiber recycling technology utilizes the inherent energy in composites for fuel and preserves the structural value of glass fiber (GF) and carbon fiber (CF) for reuse. It can also be used for recycling other waste streams at the same facility, thus spreading capital risk and recovery across multiple industries and achieving economy of scale. A second goal to be addressed was the viability of using recycled fibers as reinforcing materials in additive manufacturing applications.
The ICDC focused largely on heterogeneous catalysis of oxidative chemical transformations, especially those that came to prominence as natural gas (wet shale gas) supplanted crude oil as the primary feedstock in North America for carbon-based manufacturing of commodity chemicals. Driven by fracking-enabled access to enormous stores of shale-trapped wet gas, the feedstock change from oil to gas enabled a shift away from high-temperature catalytic breakdown of large hydrocarbons to catalytic build-up of C4, C3, C2, and C1 into desired intermediates, but now at much more modest temperatures. In turn, the much lower thermal force accompanying these kinds of reactions facilitates catalyst stabilization and characterization, including operando characterization, while slowing on-stream catalyst evolution (i.e., preserving catalyst structure), accommodating catalyst complexity, and enabling creation of uniform arrays of heterogeneous catalysts with single-atom or near-single-atom precision—at least for model systems. Together these features have proven to be a boon to hypothesis-driven catalysis science and to experimentally testable computation-based design advances and mechanism explication. The ICDC has greatly facilitated the integration of catalyst synthesis, advanced characterization, performance assessment, and computational elucidation—with advances in each of these realms spurring achievements in catalysis science and complementary realms.
We introduce and investigate matrix approximation by decomposition into a sum of radial basis function (RBF) components. An RBF component is a generalization of the outer product between a pair of vectors, where an RBF function replaces the scalar multiplication between individual vector elements. Even though the RBF functions are positive definite, the summation across components is not restricted to convex combinations and allows us to compute the decomposition for any real matrix that is not necessarily symmetric or positive definite. We formulate the problem of seeking such a decomposition as an optimization problem with a nonlinear and non-convex loss function. Several modern versions of the gradient descent method, including their scalable stochastic counterparts, are used to solve this problem. We provide extensive empirical evidence of the effectiveness of the RBF decomposition and that of the gradient-based fitting algorithm. While being conceptually motivated by singular value decomposition (SVD), our proposed nonlinear counterpart outperforms SVD by drastically reducing the memory required to approximate a data matrix with the same L2 error for a wide range of matrix types. For example, it leads to 2 to 6 times memory save for Gaussian noise, graph adjacency matrices, and kernel matrices. Moreover, this proximity-based decomposition can offer additional interpretability in applications that involve, e.g., capturing the inner low-dimensional structure of the data, retaining graph connectivity structure, and preserving the acutance of images.
With the development of nanometer accuracy stitching interferometry, ion beam figuring (IBF) of x-ray mirrors can now be achieved with unprecedented performance. However, the process of producing x-ray diffraction gratings on these surfaces may degrade the figure quality due to process errors introduced during the ruling of the grating grooves. To address this challenge, we have investigated the post-production correction of gratings using IBF, where stitching interferometry is used to provide in-process feedback. A concern with ion beam correction in this case is that ions will induce enough surface mobility of atoms to cause smoothing of the grating structure and degradation of diffraction efficiency. In this study we found however that it is possible to achieve a nanometer-level planarity of the global grating surface with IBF, while preserving the grating structure. The preservation was so good, that we could not detect a change in the diffraction efficiency after ion beam correction. This is of major importance in achieving ultra-high spectral resolution, and the preservation of brightness for coherent x-ray beams.
Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.