Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PBE”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

First-principles study of the Stark shift effect on the zero-phonon line of the NV center in diamond

Point defects in semiconductors are attractive candidates for quantum information science applications owing to their ability to act as spin-photon interface or single-photon emitters. However, the coupling between the change of dipole moment upon electronic excitation and stray electric fields in the vicinity of the defect, an effect known as Stark shift, can cause significant spectral diffusion in the emitted photons. In this work, using first principles computations, we revisit the methodology to compute the Stark shift of point defects up to the second order. The approach consists of applying an electric field on a defect in a slab and monitoring the changes in the computed zero-phonon line (i.e., difference in energy between the ground and excited state) obtained from constraining the orbital occupations (constrained-DFT). Here, we study the Stark shift of the negatively charged nitrogen-vacancy (NV) center in diamond using this slab approach. We discuss and compare two approaches to ensure a negatively charged defect in a slab and we show that converged values of the Stark shift measured by the change in dipole moment between the ground and excited states (Δ⁢μ) can be obtained. We obtain a Stark shift of Δ⁢μ = 2.68⁢D using the semilocal GGA-PBE functional and of Δ⁢μ = 2.23⁢D using the HSE hybrid functional. These values are in good agreement with experimental results. We also show that modern theory of polarization can be used on constrained-DFT to obtain Stark shifts in very good agreement with the slab computations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Putting error bars on density functional theory dataset

This dataset contains submission files and raw output files from high-throughput DFT simulations to analyze the systemic errors in lattice constant, bulk moduli and formation energy predictions for a range of binary and ternary oxides using four exchange correlation functionals (LDA, PBE, PBEsol and vdW-DF-C09). This data was then used as the basis for employing materials informatics methods to predict the expected errors in the lattice constants of the studied compounds. Predicted errors were also used to better the DFT-predicted lattice parameters. Our results emphasize the link between the computed errors and the electron density and hybridization errors of a functional. In essence, these results provide “error bars” for choosing a functional for the creation of high-accuracy, high-throughput datasets as well as avenues for the development of XC functionals with enhanced performance, thereby enabling the accelerated discovery and design of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Lattice structure and dynamics of sparse molecular crystals: OsO4 and RuO4

This dataset contains input and output files from DFT simulations used to reproduce the electronic and phonon calculations of bulk and molecular OsO₄ and RuO₄. The files include data from initial electronic structure and phonon calculations performed using different functionals: PBE, PBEsol, and vdW-DF-optB86b. The computed phonon frequencies are compared with Raman crystal and gas-phase frequencies from existing literature, while the calculated phonon density of states (PhDOS) is compared with experimental PhDOS data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

MnRhBi3: A Cleavable Antiferromagnetic Metal

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (Chem. Mater. 2024, 36, 11306-11316). The calculations characterize MnRhBi3, an orthorhombic (Cmmm) van der Waals-layered intermetallic compound that cleaves easily between neighboring Bi layers. The dataset is organized into three calculation types: (i) Bulk: Structural relaxations of the periodic MnRhBi3 crystal in antiferromagnetic (AFM) and ferromagnetic (FM) configurations, using the vdW-DF-optB86b functional. These provide the equilibrium lattice constants, magnetic energy differences (AFM is 0.5 meV/f.u. lower than FM), and magnetic moments (4.4 µB/Mn, 0.17 µB/Rh, 0.18 µB/Bi) reported in Table 1 of the main text. (ii) Slab: Same magnetic configurations computed with an 18 Ang vacuum layer introduced between Bi layers, used to calculate the cleavage energy Ec = 0.56 J/m2 (AFM) and 0.57 J/m2 (FM), establishing MnRhBi3 as a van der Waals-layered material comparable to graphite, MoS2, and CrI3. (iii) ELF: Single-point calculation on the relaxed bulk AFM geometry with LELF=.TRUE., producing the ELFCAR file used to generate electron localization function isosurfaces and contour maps (Fig. 2, main text) showing Bi lone pairs directed into the van der Waals gaps. All folders contain CONTCAR, INCAR, KPOINTS, OUTCAR, and POSCAR. The ELF/ folder additionally contains ELFCAR. Calculations were performed using VASP 6.3.2 with PBE + vdW-DF-optB86b, PAW potentials, and an energy cutoff of 800 eV.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Programmable Phase Selection between Altermagnetic and Noncentrosymmetric Polymorphs of MnTe on InP via Molecular Beam Epitaxy

This dataset contains DFT input and output files supporting the theoretical modeling in the associated publication (ACS Appl. Mater. Interfaces 2026, 18, 15654-15664). The calculations model the interfacial energetics of two MnTe polymorphs — NiAs-MnTe (hexagonal, alpha phase) and ZnS-MnTe (cubic, gamma phase) — on InP(111) substrates with two surface terminations: In-terminated InP(111)A and P-terminated InP(111)B. This gives four interface configurations: NiAs on In-terminated (experimentally observed), NiAs on P-terminated (computed for comparison), ZnS on In-terminated (computed for comparison), and ZnS on P-terminated (experimentally observed). The dataset is organized into four calculation types, each covering all four polymorph/termination combinations: (i) Slabs: Pristine MnTe/InP heterostructure slabs used to compute total energies and interface energy densities (Eint) for all four configurations, as reported in Fig. 6 of the main text. (ii) Disorder: Same slab geometries with a P_Te + Te_P antisite defect pair introduced near the interface, used to assess chemical intermixing effects on interface stability (Fig. S8, SI). (iii) Strain: Pristine slab calculations with in-plane lattice parameters strained by -1% and +1% relative to the InP lattice constant, used to evaluate strain-dependent interface energetics (Fig. S9, SI). (iv) Charge_Density: Single-point calculations on the full heterostructure, the isolated InP slab, and the isolated MnTe slab at fixed geometry, used to compute differential charge density plots showing interfacial charge accumulation and depletion as a function of surface termination (Fig. S10, SI). Each calculation folder contains INCAR, KPOINTS, POSCAR, CONTCAR, OUTCAR, and POTCAR_info.txt (PAW potential information, excluding the full POTCAR due to VASP licensing restrictions). The calculations were performed using VASP 6.4.3 with PBE exchange-correlation, PAW potentials, a Hubbard correction of Ueff = 5 eV on Mn d-states, and A-type AFM spin initialization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Further development of finite-temperature density functional theory (Technical Report on FG02-08ER46496)

The final year of funding was spent supporting research into extending conditional probability theory to generate the thermal dependence of PBE. The necessary background work led to a paper being published on the uniform electron gas (Dennis Perchak, Ryan J. McCarty, and Kieron Burke, Phys. Rev. B 105, 165143 (2022).). The graduate student John Kozlowski also contributed to a mathematical paper about DFT (Steven Crisostomo, Ryan Pederson, John Kozlowski, Bhupalee Kalita, Antonio C. Cancio, Kiril Datchev, Adam Wasserman, Suhwan Song, and Kieron Burke, Letters in Mathematical Physics 113, 42 (2023).).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Powering the Blue Economy: A Survey of Station-Keeping Methods for Mooringless Platforms

The term “ocean platform” is used to reference everything from stationary, typically moored, buoys to mobile water vehicles, whether they operate on the ocean’s surface or underwater. For certain applications for which relatively stationary station-keeping conditions are desired, the use of mooring systems is not always a viable alternative either for economic, environmental, regulatory, or otherwise practical reasons, or a combination thereof, (e.g., short deployments, sensitive ecosystems, very deep project sites). Maintaining a platform at a single waypoint or reference location without being moored would require additional control systems and a power source to counteract the drift forces that would naturally displace it. Mobile platforms, which are usually untethered except for remotely operated vehicles, typically require energy input to power their station-keeping capabilities so that they hold or control their location in the ocean. Currently, most of these platforms use combustion engines or batteries for this purpose, which, depending on the specific systems, may be costly, pollute the environment, or create limitations on the length of the deployment. However, powering this kind of platforms with surrounding renewable resources (waves, currents, winds, or sun) has been identified as a promising solution to expand their application. The intent of this report is to investigate station-keeping methods for various ocean platforms that are not moored or otherwise anchored to the ocean floor, or another platform or vessel, paying particular interest to technologies that use marine renewable resources to power their operation, because that is of particular interest to the U.S. Department of Energy’s Powering the Blue Economy (PBE) initiative. As a first step, 72 articles and technical reports related to mooringless station-keeping methods were collected for review. The preliminary literature review provided a broad overview of common themes across the literature from which a descriptive methodology for analyzing various platforms was developed. That is, station-keeping methods were categorized based on their predominant energy source and consumption (renewable, nonrenewable, or hybrid if the platform uses renewable and nonrenewable resources equally), and their localization strategy (drift reduction, “path-planning or “waypoint-holding”). In addition, platform types were segregated into the following groups: buoys, surface drifters, and unoccupied surface vehicles (USVs); offshore renewable energy systems; and unoccupied underwater vehicles (UUVs). The main types of station-keeping methods encountered in this report achieve their intended localization strategy by means of drift mitigation, steering, and/or propulsion. Drift mitigation is commonly accomplished via drogues and sea anchors. Stand-along steering subsystems use control surfaces (e.g., ship rudder, wing sail, etc.) that react to ocean currents, waves, or winds to provide varying-degrees of course adjustments. Combined steering and propulsion subsystems include differential thrusters, directional thrusters separate from a primary thruster that cause the platform to pitch up/down or yaw clockwise/counterclockwise, or vectored thrusters that direct the propulsion in a range of directions relative to the platform’s local coordinate system. Propulsion is often achieved by running a motor and applying active control strategies but can also involve buoyancy shifts and using sails to generate lifting forces that propel a platform in a desired direction. Future research is primarily expected to take place in the form of a technoeconomic analysis that would aim to determine the technological viability, cost, and added value of mooringless station-keeping use cases identified through this research, including docking for UUV recharging or for georeferencing drifter buoys, deep-sea floating wind farms, U.S. Navy sonar arrays, and a Pacific Ocean wave buoy network.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

KCl-UCl 3 molten salts investigated by Ab Initio Molecular Dynamics (AIMD) simulations: A comparative study with three dispersion models

Ab Initio Molecular Dynamics (AIMD) simulations are performed on molten KCl-UCl 3 salt mixtures to determine energies, heat capacities, and densities. The density-dependent energy correction (DFT-dDsC), Grimme et al.’s DFT-D3, and Langreth & Lundqvist (vdW-cx) models are used for dispersion forces and combined with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation potential with a Hubbard U param eter for the 5$f$ electrons of uranium. After validating predictions for the end-member systems to literature data, KCl-UCl 3 mixtures are studied at select temperatures. Densities and energies both deviate from ideal solution behavior, with the maximum deviation occurring around 36% UCl 3 for mixing energies and slightly lower (29% UCl 3 ) for densities. Compared to the NaCl-UCl 3 system, which was previously investigated using the same simulation methodologies, the KCl-UCl 3 density and mixing energy deviations from ideal solution behavior are larger by almost a factor of two. No deviation from ideal solution behavior for heat capacity was observed. The AIMD predictions for mixing energies and densities agree qualitatively with experimental data, though the spread in data obtained from the various dispersion force models utilized, measurements, and empirical estimates makes strong conclusions difficult. The dependence of thermodynamic and thermophysical properties on composition is correlated with the local chemistry of the solution phase, in particular, the tendency of UCl 3 to form network structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗
DOE OSTI2026

A Play-Based Exploration of CO 2 Storage in the Illinois Basin: Development and Application of the Play-Based Exploration Workflow

This report documents the workflow developed and applied for a play-based exploration (PBE) of carbon capture and storage (CCS) potential in the Illinois Basin, conducted under U.S. Department of Energy Cooperative Agreement DE-FE0032366. The project was led by the Illinois State Geological Survey (ISGS) at the University of Illinois Urbana-Champaign in collaboration with Visage Energy. The goal of this workflow is to identify locations where subsurface, surface, and societal conditions align to support safe and effective geologic CO 2 storage. To achieve this, the project team developed a systematic, seven-step workflow for conducting play-based exploration: (1) defining play elements, (2) collecting data, (3) constructing a geodatabase, (4) defining suitability criteria, (5) conducting suitability mapping, (6) generating Common Risk Segment (CRS) maps, and (7) creating composite play maps. The workflow is designed to be repeatable and adaptable for CCS assessments in other regions or for alternative subsurface energy applications such as hydrogen or natural gas storage.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

A Play-Based Exploration of CO 2 Storage in the Illinois Basin (Final Technical Report)

This report documents the results of "A Play-Based Exploration of CO₂ Storage in the Illinois Basin" (DE-FE0032366), a project funded by the U.S. Department of Energy Office of Fossil Energy and Carbon Management and conducted by the Illinois State Geological Survey (ISGS) at the University of Illinois Urbana-Champaign in partnership with Visage Energy. The project adapted play-based exploration (PBE), a systematic basin-scale evaluation methodology from the petroleum industry, to screen areas of Illinois for commercial geologic carbon storage (GCS) in Cambro-Ordovician strata. The traditional play concept was expanded to encompass three play element groups, subsurface geologic factors, surface features, and societal factors, yielding 24 play elements with defensible suitability criteria applied through a five-tier classification scheme. An integrated geospatial database was assembled from ISGS, MGSC, MRCI, NATCARB, and public data sources, supported by significant data-improvement work including correction of legacy well locations, digitization of more than 1,300 well construction records using the DOE CATALOG team's OGRRE tool, compilation of a statewide 2D seismic database, and production of a refined fault and fold geodatabase.

58 GEOSCIENCES↗
DOE OSTI2025

Datasets for Custom-trained Machine-learning Interatomic Potentials: Nitric Acid Aqueous Solution

This dataset was generated using an iterative active learning strategy with the ArcaNN software package (https://github.com/arcann-chem/arcann_training) to train machine-learning interatomic potentials (MLIPs) for aqueous nitric acid. Each active-learning cycle consisted of three stages: (1) training, (2) exploration, and (3) labeling. The initial training set comprised approximately 800 randomly selected configurations from a previous study by Lewis et al. (https://doi.org/10.1021/jp205510q), which investigated nitric acid solutions at 2, 3, 4, and 5 mol/L. For all configurations, single-point calculations of atomic forces and total energies were performed at the quantum density functional theory BLYP-D2 and PBE-D3 levels of theory using the CP2K Quickstep module. Valence electrons were treated explicitly, while core electrons on all atoms were represented by norm-conserving Goedecker–Teter–Hutter (GTH) pseudopotentials. Long-range dispersion interactions were accounted for using Grimme dispersion corrections. Wave functions were expanded in a mixed Gaussian-and-plane-wave scheme using TZV2P-MOLOPT basis sets for all elements and an 800 Ry auxiliary plane-wave cutoff for the electron density. Self-consistent field convergence was accelerated using orbital transformation and Direct Inversion in the Iterative Subspace, with a convergence threshold of 10^{-6}. All single-point calculations were carried out in periodic orthorhombic cells whose dimensions match those of the molecular configurations sampled from earlier trajectories. The CELL_REF keyword in CP2K was used to define a fixed reference cell, ensuring consistency in the reference data used for MLIP training, particularly when cell fluctuations are present in NpT simulations. The resulting high-fidelity energies and forces constitute the ground-truth labels used to train the MLIPs contained in this dataset.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2023

A density functional theory analysis of the adsorption and surface chemistry of inorganic iodine species on graphitea

In the event of a nuclear accident, fission products may be released into the environment. The release of 131 I is of particular concern to human health. Iodine can be captured using a number of materials and frequently, this is accomplished with activated carbon impregnated with organic bases. Previous studies have used DFT and the graphite (0001) surface as a surrogate for adsorption, those studies focus on the species I • , I 2 , and CH 3 I. In this work we perform an ab initio study of the adsorption onto the surface of a graphite sheet of I 2 , CH 3 I, and inorganic acidic iodine species (HI, HOI, HIO 2 , and HIO 3 ), which were selected to examine the possible effect of oxidation state on adsorption. The PBE exchange-correlation functional with D3 dispersion was employed. It was found that for molecular iodine, the iodine atoms tended to either situate above the center of a hexagonal site on the graphite or directly atop a carbon atom with the lighter components resting closer to the graphite. For each species the relative binding energies spanned the range of 21–33 kJ mol -1 and graphite-iodine distance was in the range of 3.52–3.93 Å. In all cases we found no significant charge transfer between the iodine species and the graphite, thus we conclude that all the iodine species studied undergo strong physisorption to the graphite.

Ritzmann, Andrew M.↗
DOE OSTI2022

Powering the Blue Economy: Exploring Opportunities for Marine Renewable Energy in Maritime Markets

To spur economic growth and revitalize the ocean the U.S. Department of Energy’s (DOE’s) Water Power Technologies Offce (WPTO) launched the Powering the Blue Economy™ (PBE) initiative, which aims to foster long-term, sustainable growth of the blue economy by: • Protecting the ocean and understanding and leveraging its immense power • Learning the power needs of emerging coastal and maritime markets • Advancing marine renewable energy technologies. Remote and island communities, areas hit hard by natural disasters, and ocean researchers can also beneft from renewable energy found in the ocean.

Powering the blue economy, marine, marine energy↗
DOE OSTI2022

Powering the Blue Economy: Foundational Research and Development

To spur economic growth and revitalize the ocean, the U.S. Department of Energy's (DOE's) Water Power Technologies Office (WPTO) launched the Powering the Blue Economy™ (PBE) initiative, which aims to foster long-term, sustainable growth of the blue economy by protecting the ocean and understanding and leveraging its immense power, learning the power needs of emerging coastal and maritime markets, and advancing marine renewable energy technologies.

marine energy↗
DOE OSTI2022

Soft and transferable pseudopotentials from multi-objective optimization

Ab initio pseudopotentials are a linchpin of modern molecular and condensed matter electronic structure calculations. In this work, we employ multi-objective optimization to maximize pseudopotential softness while maintaining high accuracy and transferability. To accomplish this, we develop a formulation in which softness and accuracy are simultaneously maximized, with accuracy determined by the ability to reproduce all-electron energy differences between Bravais lattice structures, whereupon the resulting Pareto frontier is scanned for the softest pseudopotential that provides the desired accuracy in established transferability tests. We employ an evolutionary algorithm to solve the multi-objective optimization problem and apply it to generate a comprehensive table of optimized norm-conserving Vanderbilt (ONCV) pseudopotentials (https://github.com/SPARC-X/SPMS-psps). Here, we show that the resulting table is softer than existing tables of comparable accuracy, while more accurate than tables of comparable softness. The potentials thus afford the possibility to speed up calculations in a broad range of applications areas while maintaining high accuracy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Bridging the gap between simulated and experimental ionic conductivities in lithium superionic conductors

Lithium superionic conductors (LSCs) are of major importance as solid electrolytes for next-generation all-solid-state lithium-ion batteries. While ab initio molecular dynamics have been extensively applied to study these materials, there are often large discrepancies between predicted and experimentally measured ionic conductivities and activation energies due to the high temperatures and short time scales of such simulations. Here, we present a strategy to bridge this gap using moment tensor potentials (MTPs). We show that MTPs trained on energies and forces computed using the van der Waals optB88 functional yield much more accurate lattice parameters, which in turn leads to accurate prediction of ionic conductivities and activation energies for the Li 0·33 La 0·56 TiO 3 , Li 3 YCl 6 and Li 7 P 3 S 11 LSCs. NPT MD simulations using the optB88 MTPs also reveal that all three LSCs undergo a transition between two quasi-linear Arrhenius regimes at relatively low temperatures. This transition can be traced to an increase in the number and diversity of diffusion pathways, in some cases with a change in the dimensionality of diffusion. Furthermore, this work presents not only an approach to develop high accuracy MTPs, but also outlines the diffusion characteristics for LSCs which is otherwise inaccessible through ab initio computation.

25 ENERGY STORAGE↗