Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cluster computing”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetic field induced phase transitions of the triangular spin model

In this work, motivated by a recent inelastic neutron scattering experiment on YbMgGaO 4 , I study the homogeneous spin model on the triangular lattice using the cluster mean-field theory and exact diagonalization. Phase diagrams with and without magnetic field are computed. The results show that the magnetic field can induce the phase transition from the spin liquid state (or the stripe state) to the 120 ° antiferromagnetic state. These phase transitions are suppressed by the next-nearest-neighbor exchange interaction J 2 / J 1 and vanish as J 2 / J 1 > 0.11 . I analyze a parameter space at J 2 / J 1 = 0.1 , where a field induced phase transition can occur, and find that the deviation of the theoretical spin excitation energy from the experimental data is only about 0.054. These results imply that an effective homogeneous spin model still applies to YbMgGaO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Pairing correlations in the cuprates: A numerical study of the three-band Hubbard model

We study the three-band Hubbard model for the copper oxide plane of the high-temperature superconducting cuprates using determinant quantum Monte Carlo and the dynamical cluster approximation (DCA) and provide a comprehensive view of the pairing correlations in this model using these methods. Specifically, we compute the pair-field susceptibility and study its dependence on temperature, doping, interaction strength, and charge-transfer energy. Using the DCA, we also solve the Bethe-Salpeter equation for the two-particle Green's function in the particle-particle channel to determine the transition temperature to the superconducting phase on smaller clusters. Our calculations reproduce many aspects of the cuprate phase diagram and indicate that there is an “optimal” value of the charge-transfer energy for the model where T c is largest. Finally, these results have implications for our understanding of superconductivity in both the cuprates and other doped charge-transfer insulators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Single-Walled Carbon Nanotubes with Confined Chalcogens as the Catalysts and Electrodes for Oxygen Reduction Reaction in Fuel Cells

The goal of this project is to synthesize and characterize a new non-metal electro-catalyst for oxygen reduction reaction (ORR) for fuel cell applications. The intended catalyst is a composite material composed of sulfur chains encapsulated in narrow diameter single-walled carbon nanotubes (S@SWNTs). S@SWNTs were successfully synthesized through sulfur vapor infusion method, and validated with Raman spectroscopy. However, the electrochemical analysis on the ORR catalytic activity of S@SWNTs indicated that it had low ORR catalytic activity. Our theoretical study based on density functional theory (DFT) revealed that the poor oxygen adsorption (low binding energy) on the surface of S@SWNTs was the bottleneck of the entire catalytic reaction. The focus of the project was subsequently pivoted to the development of a new non-noble metal ORR catalyst that could provide durability in acidic electrolyte. The hypothesis was to encapsulate small iron (Fe) clusters in SWNTs (Fe@SWNTs) can provide ORR electro-catalytic activity and long durability in acidic environment. DFT-based computational studies were carried out to explore the feasibility of the Fe@SWNTs catalyst. The theoretical study indicated that Fe@SWNTs indeed could catalyze the ORR with lower theoretical overpotential than platinum (Pt). However, its weaker bonding energy with oxygen was the bottleneck of the overall reaction. On the other hand, Fe clusters (composed of 7 Fe atoms) encapsulated in nitrogen-doped SWNTs (Fe 7 @N 4 WSNTs) showed proper oxygen adsorption by the Fe cluster and low theoretical overpotential of ORR. Comparing to the Fe single atom catalyst on N-doped SWNTs (Fe-N 4 SWNTs), which is one the best non-noble metal ORR catalysts reported in the literature, the Fe 7 @N 4 WSNTs showed lower overpotential and better resistance to acidic environment. Fe encapsulated N-doped SWNTs were synthesized with ferrocene as the Fe precursors through vapor infusion method, and experimental validate is underway. This study theoretically demonstrated the feasibility of a new type of non-noble metal electro-catalyst for ORR that could have high catalytic activity and long durability.

08 HYDROGEN↗
DOE OSTI2025

Understanding the Existence of a Na 2 Dimer in a High-Spin State

The recent observation of a high-spin Na 2 dimer formed on the surface of liquid helium nanodroplets raises some fundamental questions, as the ground state of Na 2 is known to have zero spin. Is it protected against spontaneous dissociation? What is its binding energy and interatomic distance? Is it stable at a higher temperature? Using calculations based on density functional theory (with and without long-range interaction) and coupled cluster methods, CCSD(T), we show that the bonding in the high-spin Na 2 dimer is governed by van der Waals interaction with binding energy (bond length) varying between −0.030 eV (5.108 Å) and −0.192 eV (4.231 Å), depending on the computational method used. Thus, the experimental method used by Kresin and coworkers can be very useful to study larger metastable high-spin clusters such as Li 4 which was predicted in 1985 to have a tetrahedral structure carrying a magnetic moment of 2 μ B , while its ground state is planar and nonmagnetic.

Basis sets↗
DOE OSTI2024

Membrane lipids drive formation of KRAS4b-RAF1 RBDCRD nanoclusters on the membrane

The oncogene RAS, extensively studied for decades, presents persistent gaps in understanding, hindering the development of effective therapeutic strategies due to a lack of precise details on how RAS initiates MAPK signaling with RAF effector proteins at the plasma membrane. Recent advances in X-ray crystallography, cryo-EM, and super-resolution fluorescence microscopy offer structural and spatial insights, yet the molecular mechanisms involving protein-protein and protein-lipid interactions in RAS-mediated signaling require further characterization. This study utilizes single-molecule experimental techniques, nuclear magnetic resonance spectroscopy, and the computational Machine-Learned Modeling Infrastructure (MuMMI) to examine KRAS4b and RAF1 on a biologically relevant lipid bilayer. MuMMI captures long-timescale events while preserving detailed atomic descriptions, providing testable models for experimental validation. Both in vitro and computational studies reveal that RBDCRD binding alters KRAS lateral diffusion on the lipid bilayer, increasing cluster size and decreasing diffusion. RAS and membrane binding cause hydrophobic residues in the CRD region to penetrate the bilayer, stabilizing complexes through β-strand elongation. These cooperative interactions among lipids, KRAS4b, and RAF1 are proposed as essential for forming nanoclusters, potentially a critical step in MAP kinase signal activation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Cluster perturbation theory X: A parallel implementation of Lagrangian perturbation series for the coupled cluster singles and doubles ground-state energy through fifth order

We describe an efficient implementation of cluster perturbation and Møller-Plesset Lagrangian energy series through fifth order that target the coupled cluster singles and doubles energy utilizing the resolution of the identity approximation. We illustrate the computational performance of the implementation by performing ground state energy calculations on systems with up to 1200 basis functions using a single node and by comparison to conventional CCSD calculations. We further show that our hybrid MPI/OMP parallel implementation that also utilizes graphical processing units can be used to obtain fifth order energies on systems with almost 1200 basis functions with a 90 minute "time to solution" running on Frontier at Oak Ridge National Laboratory.

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

Boosting H I -Galaxy Cross-Clustering Signal through Higher-Order Cross-Correlations

After reionization, neutral hydrogen (${\rm H\, \small {I}}$) traces the large-scale structure (LSS) of the Universe, enabling ${\rm H\, \small {I}}$ intensity mapping (IM) to capture the LSS in 3D and constrain key cosmological parameters. We present a new framework utilizing higher-order cross-correlations to study ${\rm H\, \small {I}}$ clustering around galaxies, tested using real-space data from the IllustrisTNG300 simulation. This approach computes the joint distributions of k-nearest neighbor (kNN) optical galaxies and the ${\rm H\, \small {I}}$ brightness temperature field smoothed at relevant scales (the kNN-field framework), providing sensitivity to all higher-order cross-correlations, unlike two-point statistics. To simulate ${\rm H\, \small {I}}$ data from actual surveys, we add random thermal noise and apply a simple foreground cleaning model, filtering out Fourier modes of the brightness temperature field with k ∥ < k min,∥ . Under current levels of thermal noise and foreground cleaning, typical of a Canadian Hydrogen Intensity Mapping Experiment (CHIME)-like survey, the ${\rm H\, \small {I}}$-galaxy cross-correlation signal in our simulations, using the kNN-field framework, is detectable at >30σ across r = [3, 12] h –1 Mpc. In contrast, the detectability of the standard two-point correlation function (2PCF) over the same scales depends strongly on the foreground filter: a sharp k ∥ filter can spuriously boost detection to 8σ due to position-space ringing, whereas a less sharp filter yields no detection. Nonetheless, we conclude that kNN-field cross-correlations are robustly detectable across a broad range of foreground filtering and thermal noise conditions, suggesting their potential for enhanced constraining power over 2PCFs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Ionization potential of radium monofluoride

The ionization potential (IP) of radium monofluoride (RaF) was measured to be 4.969(2)[10] eV, revealing a relativistic enhancement in the series of alkaline earth monofluorides. The results are in agreement with a relativistic coupled-cluster prediction of 4.981(7) eV, incorporating up to quantum electrodynamics corrections. Using the same computational methodology, an improved calculation for the dissociation energy (𝐷 0 ) of 5.54(5) eV is presented. This confirms that RaF joins the group of diatomic molecules for which 𝐷 0 > IP, paving the way for precision control and interrogation of its Rydberg states.

Wilkins, S. G. [Massachusetts Inst. of Technology ↗
DOE OSTI2021

Abelian SU(N) 1 chiral spin liquids on the square lattice

In the physics of the Fractional Quantum Hall (FQH) effect, a zoo of Abelian topological phases can be obtained by varying the magnetic field. Aiming to reach the same phenomenology in spin-like systems, in this work we propose a family of SU(N)-symmetric models in the fundamental representation, on the square lattice with short-range interactions restricted to triangular units, a natural generalization for arbitrary N of an SU(3) model studied previously where time-reversal symmetry is broken explicitly. Guided by the recent discovery of SU(2) 1 and SU(3) 1 chiral spin liquids (CSL) on similar models we search for topological SU(N) 1 CSL in some range of the Hamiltonian parameters via a combination of complementary numerical methods such as exact diagonaliza- tions (ED), infinite density matrix renormalization group (iDMRG) and infinite Projected Entangled Pair State (iPEPS). Extensive ED on small (periodic and open) clusters up to N = 10 and an innovative SU(N)-symmetric version of iDMRG to compute entanglement spectra on (infinitely-long) cylinders in all topological sectors pro- vide unambiguous signatures of the SU(N) 1 character of the chiral liquids. An SU(4)-symmetric chiral PEPS, constructed in a manner similar to its SU(2) and SU(3) analogs, is shown to give a good variational ansatz of the N = 4 ground state, with chiral edge modes originating from the PEPS holographic bulk-edge correspondence. Finally, we discuss the possible observation of such Abelian CSL in ultracold atom setups where the possibility of varying N provides a tuning parameter similar to the magnetic field in the physics of the FQH effect.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Identifying Entangled Physics Relationships Through Sparse Matrix Decomposition to Inform Plasma Fusion Design

We report a sustainable burn platform through inertial confinement fusion (ICF) has been an ongoing challenge for over 50 years. Mitigating engineering limitations and improving the current design involves an understanding of the complex coupling of physical processes. While sophisticated simulation codes are used to model ICF implosions, these tools contain necessary numerical approximation but miss physical processes that limit predictive capability. Identification of relationships between controllable design inputs to ICF experiments and measurable outcomes (e.g., neutron yield, neutron velocity, areal density) from performed experiments can help guide the future design of experiments and development of simulation codes, to potentially improve the accuracy of the computational models used to simulate ICF experiments. We use sparse matrix decomposition methods to identify clusters of a few related design variables. Sparse principal component analysis (SPCA) identifies groupings that are related to the physical origin of the variables (laser, hohlraum, and capsule). A variable importance analysis finds that in addition to variables highly correlated with neutron yield, such as picket power and laser energy, variables that represent a dramatic change of the ICF design, such as number of pulse steps, are also very important. The obtained sparse components are then used to train a random forest (RF) regression surrogate for predicting total yield. The RF performance on the training and testing data compares with the performance of the RF trained using all the design variables considered. This work is intended to inform design changes in future ICF experiments by augmenting the expert intuition and simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Synthase-selected sorting approach identifies a beta-lactone synthase in a nudibranch symbiotic bacterium

Background - Nudibranchs comprise a group of > 6000 marine soft-bodied mollusk species known to use secondary metabolites (natural products) for chemical defense. The full diversity of these metabolites and whether symbiotic microbes are responsible for their synthesis remains unexplored. Another issue in searching for undiscovered natural products is that computational analysis of genomes of uncultured microbes can result in detection of novel biosynthetic gene clusters; however, their in vivo functionality is not guaranteed which limits further exploration of their pharmaceutical or industrial potential. To overcome these challenges, we used a fluorescent pantetheine probe, which produces a fluorescent CoA-analog employed in biosynthesis of secondary metabolites, to label and capture bacterial symbionts actively producing these compounds in the mantle of the nudibranch Doriopsilla fulva. Results - We recovered the genome of Candidatus Doriopsillibacter californiensis from the Ca. Tethybacterales order, an uncultured lineage of sponge symbionts not found in nudibranchs previously. It forms part of the core skin microbiome of D. fulva and is nearly absent in its internal organs. We showed that crude extracts of D. fulva contained secondary metabolites that were consistent with the presence of a beta-lactone encoded in Ca. D. californiensis genome. Beta-lactones represent an underexplored group of secondary metabolites with pharmaceutical potential that have not been reported in nudibranchs previously. Conclusions - Altogether, this study shows how probe-based, targeted sorting approaches can capture bacterial symbionts producing secondary metabolites in vivo.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Using Apptainer in a Pilot-based Distributed Workload

GlideinWMS is a pilot and pressure-based workload manager for distributed scientific computing. Many experiments like CMS and Fermilab’s Neutrino experiments use it to provision elastic clusters for their analysis and simulations, split into close to a million concurrent jobs. Most user jobs require containers, and the pilots use Apptainer to set up the desired platform. For the pilots that run as regular batch jobs, Apptainer is safer, lighter, and easier to use than other containerization solutions. Many images used by the pilots are expanded SIF images distributed via the CernVM-FS: this combination is very efficient. At Fermilab, for example, we store on GitHub Dockerfiles that mimic the platform in the worker nodes of local clusters. GitHub workflows build and push the images to Docker Hub, and a service periodically pulls and converts them to the expanded SIF images in the CernVM-FS, so the scientists can find a familiar environment everywhere. Apptainer has also been used to run services inside the pilot jobs, like benchmarks that characterize the worker node being used, or a Triton Inference Server that allows sharing a GPU with all the jobs that run in parallel on a node.

Mambelli, Marco [Fermilab] (ORCID:0000000294892681↗
DOE OSTI2023

Localized High-Concentration Electrolytes Get More Localized

Liquid electrolytes in batteries are typically treated as macroscopically homogeneous ionic transport media, despite having complex chemical composition and atomistic solvation structures, thus leaving a knowledge gap of microstructural characteristics. A promising electrolyte system design has recently been developed with the title of a localized high-concentration electrolyte (LHCE). These electrolytes are a summation of three components – a charge-carrying salt, a solvent that solvates the salt, and a diluent that is sparingly interactive with the salt in terms of ion solvation. This creates the perception of a high-concentration electrolyte (HCE) with high anion-cation interactions that supports a salt-derived, stable solid-electrolyte interphase, while an LHCE improves bulk electrolyte properties, such as ionic conductivity and viscosity, when compared to an HCE. Much trial by error design of this LHCE concept has been supported by simulations to confirm ionic interactions, cycling data to confirm cell stability, and post-mortem analysis to provide evidence of mechanisms. However, a systematic design to extend the capabilities of LHCEs is still to be desired. Here, we analyzed a specific LHCE mixture to dictate the best viable mixture that extends the localization of solvation structures while maintaining viable cell-level capabilities. We seek to observe and understand solvent-to-salt interactions that impact electrolyte performance. These interactions range from single-ion solvation, multidentate interactions (i.e., solvation-driven clustering), to more widespread network-level interactions. With the support of both experimental and computational analysis of different electrolyte features, we have begun a methodology of formulating the most ideal LHCE mixtures for improving the cyclability of high-capacity battery chemistries.

25 ENERGY STORAGE↗
DOE OSTI2022

Scalable probabilistic estimates of electric vehicle charging given observed driver behavior

To prepare for rapid growth in global electric vehicle adoption, grid and policy planners depend on detailed forecasts of future charging demand. In this paper we propose a novel holistic, scalable, probabilistic framework to produce large-scale estimates of electric vehicle charging load for long-term planning that capture real drivers’ charging patterns. Our framework captures the uncertainty and stochasticity in charging demand by taking a graphical modeling approach. It has three core elements: driver groups, charging segment choices, and charging session time and energy requirements. The framework uses hierarchical clustering to group drivers by their charging histories, capturing their heterogeneous behaviors and preferences across different segments or types of charging. The framework uses probabilistic mixture models for each driver group’s sessions to identify the unique charging behaviors observed within each segment. We illustrate its application with a large data set from California, profiling the charging patterns and unique driver clusters it identifies. Using the model knobs representing drivers’ battery capacities, behavior, and segment access we present scenarios for California’s charging demand in 2030 with 8 million passenger electric vehicles. Peak charging demand ranged from 3.3 to 8.7 GW across scenarios. Furthermore, each was calculated in under 45 s on a laptop computer.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Dispersion-enhanced sequential batch sampling for adaptive contour estimation

In computer simulation and optimal design, sequential batch sampling offers an appealing way to iteratively stipulate optimal sampling points based upon existing selections and efficiently construct surrogate modeling. Nonetheless, the issue of near duplicates poses tremendous quandary for sequential learning. It refers to the situation that selected critical points cluster together in each sampling batch, which are individually but not collectively informative towards the optimal design. Near duplicates severely diminish the computational efficiency as they barely contribute extra information towards update of the surrogate. To address this issue, we impose a dispersion criterion on concurrent selection of sampling points, which essentially forces a sparse distribution of critical points in each batch, and demonstrate the effectiveness of this approach in adaptive contour estimation. Specifically, we adopt Gaussian process surrogate to emulate the simulator, acquire variance reduction of the critical region from new sampling points as a dispersion criterion, and combine it with the modified expected improvement (EI) function for critical batch selection. The critical region here is the proximity of the contour of interest. This proposed approach is vindicated in numerical examples of a two-dimensional four-branch function, a four-dimensional function with a disjoint contour of interest and a time-delay dynamic system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Arrested-motility states in populations of shape-anisotropic active Janus particles

The emergence of large-scale collective phenomena from simple interactions between individual units is a hallmark of active matter systems. Active colloids with alignment-dominated interparticle interactions tend to develop orientational order and form motile coherent states, such as flocks and swarms. Alternatively, a combination of self-propulsion and excluded-volume interactions results in self-trapping and active phase separation into dense clusters. Here, we reveal unconventional arrested-motility states in ensembles of active discoidal particles powered by induced-charge electrophoresis. Combining experiments and computational modeling, we demonstrate that the shape asymmetry of the particles promotes the hydrodynamically assisted formation of active particles’ bound states in a certain range of excitation parameters, ultimately leading to a spontaneous collective state with arrested motility. Unlike the jammed clusters obtained through self-trapping, the arrested-motility phase remains sparse, dynamic, and reconfigurable. The demonstrated mechanism of phase separation seeded by bound state formation in ensembles of oblate active particles is generic and should be applicable to other active colloidal systems.

42 ENGINEERING↗
DOE OSTI2024

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗
DOE OSTI2020

Time-sliced quantum circuit partitioning for modular architectures

Current quantum computer designs will not scale. To scale beyond small prototypes, quantum architectures will likely adopt a modular approach with clusters of tightly connected quantum bits and sparser connections between clusters. We exploit this clustering and the statically-known control flow of quantum programs to create tractable partitioning heuristics which map quantum circuits to modular physical machines one time slice at a time. Specifically, we create optimized mappings for each time slice, accounting for the cost to move data from the previous time slice and using a tunable lookahead scheme to reduce the cost to move to future time slices. We compare our approach to a traditional statically-mapped, owner-computes model. Our results show strict improvement over the static mapping baseline. Here, we reduce the non-local communication overhead by 89.8% in the best case and by 60.9% on average. Our techniques, unlike many exact solver methods, are computationally tractable.

Baker, Jonathan M.↗