Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice Integrable Models”

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 199 records · Page 11

Machine learning surrogate for charged particle beam dynamics with space charge based on a recurrent neural network with aleatoric uncertainty

In this work, we develop a machine learning (ML) model with aleatoric uncertainty for the low energy beam transport (LEBT) region of the LANSCE linear accelerator in which we model the transport of a space-charge-dominated 750 keV proton beam through a lattice of 22 quadrupole magnets. Our ML model is developed based on data generated by a Kapchinsky–Vladimirsky (KV) envelope model of beam transport. We show that a recurrent neural network can be used as a dynamical surrogate model for fast prediction of the LEBT beam envelope. Furthermore, we endow the model with the prediction of aleatoric uncertainty and compare three different approaches. We demonstrate that the ML-based uncertainty quantification models are well calibrated and produce good estimates of the regions where the model is less certain about its predictions. This ML framework is a necessary step in the development of a real-time virtual diagnostic tool with uncertainty quantification that can be integrated into more complex downstream tasks (e.g., adaptive control or learning flexible control policies via reinforcement learning) for improved efficiency in beam operations. In future work, we plan to expand on this preliminary study by considering more realistic envelope models that include longitudinal momentum spread and dispersive effects in bending magnets, as well as particle tracking codes with 3D space charge (such as and ). Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Spontaneous Lithiation of Binary Oxides during Epitaxial Growth on LiCoO 2

Epitaxial growth is a powerful tool for synthesizing heterostructures and integrating multiple functionalities. Interfacial mixing can readily occur at temperatures required for complete film crystallization and can significantly modify the properties of layered structures, particularly for those containing energy storage materials with smaller cations. Here we show a two-step sequence involving the growth of an epitaxial LiCoO 2 cathode layer followed by the deposition of a binary transition metal oxide (WO 3 , TiO 2 , and others) in which controlled lithiation of the binary oxide occurs. Orientation-controlled epitaxial synthesis of the model solid-state-electrolyte Li 2 WO 4 and model anode material Li 4 Ti 5 O 12 occurs as WO 3 and TiO 2 nucleate their respective host lattices and attract Li ions from the underlying cathode. The cathode layer readily provides a tunable amount of Li ions for the controlled lithiation of the subsequent layer. We demonstrate that this approach can be used for energy materials discovery and exploring different combinations of epitaxial interfaces that can serve as well-defined model systems for mechanistic studies of energy storage and conversion processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Artificial-Intelligence Aided Design and Synthesis of Novel Layered 2D Multi-Principal Element Materials for Energy Storage (Final Report)

This DOE-EPSCoR project aimed to predict, synthesize, and characterize novel layered two-dimensional (2D) high-entropy materials (HEMs). These 2D-HEMs, composed of multiple principal elements in nearly equal concentrations, are distinct from traditional 2D materials (typically containing two or three elements) and conventional alloys (dominated by a single primary element with minor secondary additions). Their unique structural and compositional features enable significant lattice strain accommodation, resulting in enhanced electrode performance and potential applications in catalysis, hydrogen storage, sensing, quantum information technologies, and flexible electronics. The research focused on addressing four fundamental questions: (i) What combinations of elements can form stable and synthesizable 2D-HEMs? (ii) What mechanisms drive the stability and synthesizability of crystalline single-phase 2D-HEMs? (iii) How do local chemical disorder and defects influence the macroscopic electronic and mechanical properties? and (iv) What charge storage mechanisms are active in selectively synthesized 2D-HEMs for battery and supercapacitor electrode applications? To achieve these goals, the project employed an integrated theory-experiment approach, incorporating high-throughput first-principles calculations, theoretical modeling, data mining, experimental synthesis, and advanced characterization techniques. The advanced computing resources and state-of-the-art experimental characterization facilities at Oak Ridge National Laboratory (ORNL) were leveraged through collaboration. Beyond scientific advancements, the project contributed to workforce development. Two postdoctoral researchers and three graduate students at the University of Maine were trained through co-advising by ORNL scientists and collaborative interactions, strengthening their expertise in cutting-edge materials science.

36 MATERIALS SCIENCE↗

Topological Green’s Function Zeros in an Exactly Solved Model and Beyond

The interplay of topological electronic band structures and strong interparticle interactions provides a promising path towards the constructive design of robust, long-range entangled many-body systems. As a prototype for such systems, we here study an exactly integrable, local model for a fractionalized topological insulator. Using a controlled perturbation theory about this limit, we demonstrate the existence of topological bands of zeros in the exact fermionic Green’s function and show that in this model they do affect the topological invariant of the system, but not the quantized transport response. Close to (but prior to) the Higgs transition signaling the breakdown of fractionalization, the topological bands of zeros acquire a finite “lifetime.” We also discuss the appearance of edge states and edge zeros at real space domain walls separating different phases of the system. This model provides a fertile ground for controlled studies of the phenomenology of Green’s function zeros and the underlying exactly solvable lattice gauge theory illustrates the synergetic cross pollination between solid-state theory, high-energy physics, and quantum information science.

Bollmann, Steffen↗

Adaptive tau-leaping methods for microscopic-lattice kinetic Monte Carlo simulations

Traditional Kinetic Monte Carlo (KMC) approaches, rooted in Gillespie’s stochastic simulation algorithm, become computationally demanding in systems with a large range of timescales. The goal of this work is to propose and study new adaptive lattice-KMC time integration strategies for spatially non-uniform systems. To that end, two novel adaptive tau-leaping methods and their corresponding time integration strategies are developed based on the idea of the “n-fold” direct KMC method. These strategies allow for the simultaneous execution of multiple reactions, advancing time by adaptively selected coarse increments. We present numerical experiments comparing the proposed methods with existing approaches in a catalytic surface kinetics application involving ammonia decomposition.

Bimolecular reactions↗

Finite temperature auxiliary field quantum Monte Carlo in the canonical ensemble

Finite temperature auxiliary field-based quantum Monte Carlo methods, including determinant quantum Monte Carlo and Auxiliary Field Quantum Monte Carlo (AFQMC), have historically assumed pivotal roles in the investigation of the finite temperature phase diagrams of a wide variety of multidimensional lattice models and materials. Despite their utility, however, these techniques are typically formulated in the grand canonical ensemble, which makes them difficult to apply to condensates such as superfluids and difficult to benchmark against alternative methods that are formulated in the canonical ensemble. Working in the grand canonical ensemble is furthermore accompanied by the increased overhead associated with having to determine the chemical potentials that produce desired fillings. Given this backdrop, in this work, we present a new recursive approach for performing AFQMC simulations in the canonical ensemble that does not require knowledge of chemical potentials. To derive this approach, we exploit the convenient fact that AFQMC solves the many-body problem by decoupling many-body propagators into integrals over one-body problems to which non-interacting theories can be applied. We benchmark the accuracy of our technique on illustrative Bose and Fermi–Hubbard models and demonstrate that it can converge more quickly to the ground state than grand canonical AFQMC simulations. We believe that our novel use of HS-transformed operators to implement algorithms originally derived for non-interacting systems will motivate the development of a variety of other methods and anticipate that our technique will enable direct performance comparisons against other many-body approaches formulated in the canonical ensemble.

Shen, Tong (ORCID:0000000206292468)↗

Neural network interatomic potential-driven analysis of phase stability in Ti–V alloys at the atomistic scale

The evolution of the ω phase in titanium–vanadium (Ti–V) alloys is critical for their mechanical properties, particularly in aerospace and biomedical applications. Here, this study employs a Rapid Artificial Neural Network (RANN) potential to model the ω phase evolution at the atomistic level, demonstrating a high degree of consistency with experimental observations, unlike the Modified Embedded Atom Method (MEAM), which fails to capture this phase transformation accurately. RANN simulations replicate key phenomena such as the nucleation of α precipitates at ω/β interfaces and accurate lattice orientations, enhancing our understanding of phase stability and transformation kinetics. The findings affirm that RANN potentials can significantly improve the prediction accuracy of complex material behaviors, offering a powerful tool for designing advanced materials with tailored properties such as solute effect in various stacking fault energies. This approach not only bridges the gap between theoretical predictions and empirical data but also sets a new direction for future research in materials science, emphasizing the integration of machine learning techniques in the development and optimization of new alloys.

36 MATERIALS SCIENCE↗

Efficient lattice QCD computation of radiative-leptonic-decay form factors at multiple positive and negative photon virtualities

In previous work [D. Giusti, Methods for high-precision determinations of radiative-leptonic decay form factors using lattice QCD, Phys. Rev. D 107, 074507 (2023)], we showed that form factors for radiative leptonic decays of pseudoscalar mesons can be determined efficiently and with high precision from lattice QCD using the “three-dimensional (3D) method,” in which three-point functions are computed for all values of the current insertion time and the time integral is performed at the data-analysis stage. Here, we demonstrate another benefit of the 3D method: the form factors can be extracted for any number of nonzero photon virtualities from the same three-point functions at no extra cost. We present results for the $D_s → ℓνγ*$ vector form factor as a function of photon energy and photon virtuality, for both positive and negative virtuality, for a single ensemble with 340 MeV pion mass and 0.11 fm lattice spacing. In our analysis, we separately consider the two different time orderings and the different quark flavors in the electromagnetic current. We discuss in detail the behavior of the unwanted exponentials contributing to the three-point functions, as well as the choice of fit models and fit ranges used to remove them for various values of the virtuality. While positive photon virtuality is relevant for decays to multiple charged leptons, negative photon virtuality suppresses soft contributions and is of interest in QCD-factorization studies of the form factors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Machine-learning prediction for quasiparton distribution function matrix elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-x dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken x) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small- x reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_za)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements, especially in finer lattice units, increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning algorithms can make predictions of correlators to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different prediction accuracy and systematic errors. The predictions from smaller displacement z to larger ones work better than those for momentum p due to the higher correlation among the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

LatticeAnalytics: Strut-Level Visualization and Inspection of Additively Manufactured Lattice Structures

Additive manufacturing (AM) is revolutionizing the production of custom components with complex internal geometries, essential for high-performance applications in diverse fields such as medicine and defense. These AM parts optimize strength while minimizing weight by utilizing internal lattice structures consisting of large quantities of small interconnected struts. However, the complexity of these structures, combined with the challenges of using X-ray Computed Tomography (XCT) data, makes validation of part reliability difficult. This ultimately inhibits the development of novel parts for our collaborating material scientists. Here, we introduce LatticeAnalytics, a novel framework specifically designed for visual inspection of defects in these lattice structures. Our framework offers an end-to-end solution that includes the data management of XCT scans, enables remote access for geographically dispersed teams through a web-based dashboard, and incorporates novel visualizations. Our analysis is facilitated by a coarse alignment between the lattice’s nominal model, a spatial graph, and the XCT data. We employ a simple VR-based approach for fast and rough alignment, followed by an offline registration and identification of the struts. With the nodes and struts aligned and identified in the volume, our framework allows querying of subvolumes containing a single strut at multiple resolutions. This avoids computation over the entire lattice and also allow for easy parallelization of down-stream computations, such as strut-specific metrics. To depict a fast overview of the strut quality, we introduce two innovative visual encodings, crucial for our collaborators’ research in creating novel AM parts: the Contour View and the Roughness Map, which depict critical geometrical and surface features of individual struts in standardized two 2D views. We evaluated the integrated system through expert interviews. The feedback confirms the framework’s practicality and its effectiveness in enhancing current inspection workflows. It solves major bottlenecks for our collaborators, ultimately helping them create novel parts with advanced properties.

Miao, Haichao [Lawrence Livermore National Laborat↗

Actinides and Correlated Electron Materials

The Actinides and Correlated Electron Materials area of leadership spans Los Alamos National Laboratory competency in actinide materials research dating to the Manhattan Project as articulated in the Integrated Plutonium Science and Research Strategy and competency in strongly correlated electron systems dating back to at least the early 1980s. This area of leadership focuses on the goals of discovering, understanding, and controlling emergent electronic states and predictive performance of actinide materials. They are quintessentially linked by the fact that the physics of actinides—and plutonium in particular—are governed by strong electronic correlations. Not only is the electronic structure of actinides dictated by fine details of electron correlations, but chemical bonding and physical structure are as well. Hence, by addressing the first goal of this leadership area we can significantly accelerate progress on the second goal. To understand such matter requires probing the intertwined spin, charge, orbital, and lattice degrees of freedom with greater precision and developing models that accurately predict the consequences of these coupled degrees of freedom, on multiple length and time scales and including acute reactivity and effects of self-irradiation phenomena in these materials.

36 MATERIALS SCIENCE↗

CFD and Experimental Data Comparisons for Conventional and AFC-Enabled CRM High-Lift Configurations

Numerical simulations have been performed for a simplified high-lift (SHL) version of the Common Research Model (CRM) configuration, where the Fowler flaps of a representative conventional high-lift (CRM-HL) configuration are replaced by a set of simple hinged flaps. These hinged flaps are equipped with integrated modular active flow control (AFC) cartridges on the suction surface of the flap shoulder, and the resulting geometry is known as the CRM-SHL-AFC configuration. The main objective is to make use of AFC devices on the CRM-SHL-Configuration to produce the aerodynamic performance (lift) comparable to that of the CRM-HL configuration over a large angle of attack range encompassing maximum lift conditions. For comparison purposes, computations are also per-formed for the CRM-HL configuration. In the current paper, PowerFLOWR©, a CFD code based on the Lattice Boltzmann method (LBM), is used to simulate the entire flow field associated with the CRM-SHL-AFC configuration equipped with a combination of two different types of AFC actuators. The transonic version of the PowerFLOWR© code that has been validated for high-speed flows is used to simulate the flow field generated by the high-momentum actuators required to mitigate reversed flow regions on the suction surfaces of the main wing and the flap. This study is focused on the AFC systems and actuator arrangements that had emerged based on the parametric studies conducted previously at the nominal landing condition. Comparisons of the numerical solutions for lift and surface pressures are presented here with the experimental data, demonstrating the usefulness of CFD for predicting the flow field and lift characteristics of AFC-enabled high-lift configurations over a broad angle of attack range. The numerical solutions predict the expected trends in aerodynamic forces with angle of attack variation.

Veer N Vatsa↗

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗

Influence of misfit dislocations on ionic conductivity at oxide interfaces

Mismatched complex oxide thin films and heterostructures have gained significant traction for use as electrolytes in intermediate temperature solid oxide fuel cells, wherein interfaces exhibit variation in ionic conductivity as compared to the bulk. Although misfit dislocations present at interfaces in these structures impact ionic conductivity, the fundamental mechanisms responsible for this effect are not well understood. To this end, a kinetic lattice Monte Carlo (KLMC) model was developed to trace oxygen vacancy diffusion at misfit dislocations in SrTiO 3 /BaZrO 3 heterostructures and elucidate the atomistic mechanisms governing ionic diffusion at oxide interfaces. The KLMC model utilized oxygen vacancy migration energy barriers computed using molecular statics. While some interfaces promote oxygen vacancy diffusion, others impede their transport. Fundamental factors such as interface layer chemistry, misfit dislocation structure, and starting and ending sites of migrating ions play a crucial role in oxygen diffusivity. Molecular dynamics (MD) simulations were further performed to support qualitative trends for oxygen vacancy diffusion. Overall, the agreement between KLMC and MD is quite good, though MD tends to predict slightly higher conductivities, perhaps a reflection of nuanced structural relaxations that are not captured by KLMC. The current framework comprising KLMC modeling integrated with molecular statics offers a powerful tool to perform mechanistic studies focused on ionic transport in thin film oxide electrolytes and facilitate their rational design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of thermal properties of magnetic materials with unknown entropy

Theoretical approaches to study thermal properties of magnetic materials typically require accurate models of magnetic interactions in order to define the entropy. Here we introduce a complementary approach for examining thermal properties in magnetic systems where an accepted model for such interactions does not exist. In place of a specific model for magnetic interactions, the approach integrates measurements of temperature dependent magnetization of the studied material into a first principles computational scheme. The approach calculates system pressure from thermally disordered microstates that properly incorporate vibrational and spin subsystems at each temperature as well as the coupling between these subsystems. We apply the approach to calculate phonon modes and to investigate the anomalously low thermal expansion of the classical Invar alloy, Fe 0.65 Ni 0.35 . Here, the calculated phonon dispersions for Invar are in excellent agreement with measured data. The Invar thermal expansion is shown to remain small between 50 K and room temperature, consistent with the experimentally observed low thermal expansion value in this same temperature range. This anomalously small thermal expansion is directly connected to a small positive contribution from lattice thermal disorder that is nearly canceled by a large negative magnetic disorder contribution. By contrast, calculations for bcc Fe show a much larger thermal expansion, consistent with experiment, which is dominated by a large contribution from lattice thermal disorder that is reduced only slightly by a small negative contribution from that of magnetism. These findings give insights into the unusual nature of magnetism and spin-lattice coupling in Invar and Fe. In addition, they give promising preliminary support to the presented new methodology as a complementary way to investigate thermal properties of magnetic materials. The success achieved on Invar and Fe motivates future testing of the approach on other magnetic materials.

36 MATERIALS SCIENCE↗

Constraints on the two-dimensional pseudospin-$\frac{1}{2}$Mott insulator description of $\mathrm{Sr_2IrO_4}$

Sr 2 IrO 4 has often been described via a simple, one-band pseudospin-$\frac{1}{2}$ model subject to electron-electron interactions on a square lattice, fostering analogies with cuprate superconductors believed to be well described by a similar model. In this work we argue - based on a detailed study of the low-energy electronic structure by circularly polarized spin and angle-resolved photoemission spectroscopy combined with dynamical mean-field theory calculations - that a pseudospin-$\frac{1}{2}$ model fails to capture the full complexity of the system. We show instead that a realistic multiband Hubbard Hamiltonian, accounting for the full correlated t 2g manifold, provides a detailed description of the interplay between spin-orbital entanglement and electron-electron interactions and yields quantitative agreement with experiments. Our analysis establishes that the j 3/2 states make up a substantial percentage of the low-energy spectral weight, i.e., approximately 74% as determined from the integration of the j-resolved spectral function in the 0 to -1.64eV energy range. Therefore, the results in our work are of relevance not only to Ir-based materials but also more generally to multiorbital materials with closely spaced energy scales.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parametric Deformation of Discrete Geometry for Aerodynamic Shape Design

We present a versatile discrete geometry manipulation platform for aerospace vehicle shape optimization. The platform is based on the geometry kernel of an open-source modeling tool called Blender and offers access to four parametric deformation techniques: lattice, cage-based, skeletal, and direct manipulation. Custom deformation methods are implemented as plugins, and the kernel is controlled through a scripting interface. Surface sensitivities are provided to support gradient-based optimization. The platform architecture allows the use of geometry pipelines, where multiple modelers are used in sequence, enabling manipulation difficult or impossible to achieve with a constructive modeler or deformer alone. We implement an intuitive custom deformation method in which a set of surface points serve as the design variables and user-specified constraints are intrinsically satisfied. We test our geometry platform on several design examples using an aerodynamic design framework based on Cartesian grids. We examine inverse airfoil design and shape matching and perform lift-constrained drag minimization on an airfoil with thickness constraints. A transport wing-fuselage integration problem demonstrates the approach in 3D. In a final example, our platform is pipelined with a constructive modeler to parabolically sweep a wingtip while applying a 1-G loading deformation across the wingspan. This work is an important first step towards the larger goal of leveraging the investment of the graphics industry to improve the state-of-the-art in aerospace geometry tools.

Anderson, George R.↗

Airframe Noise Simulations of a Full-Scale Large Civil Transport in Landing Configuration

This paper summarizes the results obtained from an extensive computational campaign to accurately predict full-scale landing gear noise for large civil transports. A highly accurate digital model of a full-scale Boeing 777-300ER aircraft with as-flown nose and main landing gear components was developed for use in the simulations. Two aircraft configurations were selected: nose and main landing gear deployed with wing high-lift devices retracted, and nose and main landing gear deployed with wing high-lift devices deflected. The two configurations were simulated without and with toboggan fairings installed on the main gear to represent the principal configurations evaluated during the 2005 QTD2 flight test. All simulations were performed with the lattice Boltzmann solver PowerFLOW® to resolve and capture the highly complex, unsteady flow field in the immediate vicinity of the aircraft. The far-field noise sig-nature of the aircraft was computed via a Ffowcs-Williams and Hawkings integral approach, with flow quantities on a permeable surface enclosing the source regions used as input. Synthetic pressure records at ground array microphone locations used during the QDT2 test were employed to generate narrowband acoustic maps and integrated far-field noise spectra. With high-lift devices retracted, the predicted spectra showed that landing gear noise is equivalent to total airframe noise, with no other airframe sources appearing within 10 dB of gear peak levels. Application of a toboggan fairing to the main gear produced modest noise reductions of 1-2 dB across the resolved frequency range. With high-lift devices deflected, undercarriage noise was within 3-4 dB of the total airframe noise, thus comprising nearly half of the total airframe noise. For this configuration, the toboggan fairing did not produce a reduction in noise, corroborating trends previously observed in QTD2 flight test data.

airframe noise↗