Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure calculations”

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

Runtime performance of a GAMESS quantum chemistry application offloaded to GPUs

Summary Computational chemistry is at the forefront of solving urgent societal problems, such as polymer upcycling and carbon capture. The complexity of modeling these processes at appropriate length and time scales is mainly manifested in the number and types of chemical species involved in the reactions and may require models of several thousand atoms and large basis sets to accurately capture the chemical complexity and heterogeneity in the physical and chemical processes. The quantum chemistry package General Atomic and Molecular Electronic Structure System (GAMESS) has a wide array of methods that can efficiently and accurately treat complex chemical systems. In this work, we have used the GAMESS Effective Fragment Molecule Orbital (EFMO) method for electronic structure calculation of a challenging mesoporous silica nanoparticle (MSN) model surrounded by about 4700 water molecules to investigate the strong scaling and GPU offloading on hybrid CPU‐GPU nodes. Experiments were performed on the Perlmutter platform at the National Energy Research Scientific Computing Center. Good strong scaling and load balancing have been observed on up to 88 hybrid nodes for different settings of the execution parameters for the calculation considered here. When GPUs are oversubscribed by offloading work from multiple CPU processes, using the NVIDIA multi‐process service (MPS) has consistently reduced time to solution and energy consumed. Additionally, for some configuration parameter settings, oversubscription with MPS improved performance by up to 5.8% over the case without oversubscription.

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

Influence of Peripheral Substituents on Fe(II) Spin State in Complexes with Tridentate Schiff‐Base Ligands

Abstract Six Fe(II) complexes with substituted 6‐((quinoline‐8‐ylimino)methyl)phenolates (R 2 qsal) were synthesized. X‐ray crystal structure determination for [Fe(R 2 qsal) 2 ], where R=F ( 1 ), Cl ( 2 ), Br ( 3 ), I ( 4 ), CF 3 /H ( 5 ), or CN/H ( 6 ), revealed mononuclear complexes that crystallize either solvent‐free in case of 1 , 5 , and 6 or as CH 2 Cl 2 solvates in the case of 2 , 3 , and 4 . Complexes 1 – 4 , which contain π‐donating substituents, exhibit only high‐spin (HS) state, as shown by crystal structure analysis and magnetic measurements. In contrast, complexes 5 and 6 , which contain electron‐withdrawing substituents, exist in the low‐spin (LS) state at lower temperatures but exhibit crossover to the HS state above 300 K. Electronic structure calculations show a good correlation between the relative energies of the LS and HS states and the experimentally observed behavior.

Gakiya‐Teruya, Miguel↗
DOE OSTI2025

Driving rapid atomic order in MnAl via low-magnitude magnetic field annealing

Application of a mild (60 mT), uniform magnetic field during short-term thermal treatment of kinetically retained, atomically disordered (paramagnetic) ε-MnAl was found to deliver a significant ~50 % increase in the formation of L1 0 atomically ordered (ferromagnetic) τ-MnAl product phase, compared to that produced by conventional (i.e., zero-field) annealing under identical thermal conditions. The magnetic field, applied in a passive closed-circuit configuration during annealing, induced significant changes in the structural, magnetic, and phase evolution of the material. Computational results based on electronic structure calculations demonstrate that the effective magnetic susceptibility of τ-MnAl is sensitive to the orientation, rather than the magnitude, of an applied magnetic field in the vicinity of the Curie temperature. The uniaxial magnetocrystalline anisotropy of the L1 0 structure is proposed to act as a filter for selective propagation of the population of τ-MnAl variants that are favorably aligned with the applied field. In this manner, crystallographic “gridlock” is alleviated that would otherwise arise from the coexistence of multiple, energetically equivalent τ-phase variants within the parent ε-phase matrix. These results confirm that static, low-magnitude magnetic field annealing is able to accelerate L1 0 atomic ordering in the MnAl system and likely can exert similar influences in relevant magnetic systems, facilitating efficient tailoring of structure-sensitive magnetic properties for the manufacture of magnetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Selective phosphate removal with manganese oxide composite anion exchange membranes in membrane capacitive deionization

The discharge of excessive phosphorous into water bodies can lead to serious eutrophication threatening aquatic ecosystem. Membrane capacitive deionization (MCDI) is an effective platform for deionizing aqueous streams; however, conventional MCDI is unable to selectively remove targeted ions from a liquid mixture. Here, in this work, we fabricated manganese oxide composite anion exchange membranes (AEMs) for MCDI to enhance phosphate removal selectivity from sodium chloride-sodium dihydrogen phosphate (10:1 M ratio) aqueous mixtures. We systematically investigated several critical factors, such as constant current or voltage operation, applied voltage amount, process stream pH, and manganese oxide (Mn 2 O 3 ) content in the AEM, on phosphate removal efficiency and phosphate selectivity. A trade-off was observed between phosphate removal and selectivity when increasing the cell voltage. Under the best conditions, a MCDI unit with a 20 wt% Mn 2 O 3 composite AEM and a bipolar membrane facilitated high phosphate removal efficiency of ≥ 31.8 % and a phosphate over chloride selectivity of 1.1 while showing stability for at least 30 cycles. To help understand how Mn 2 O 3 composite AEM boosts phosphate selectivity, static electronic structure calculations were performed, and they revelated that hydrogen phosphate absorption on Mn 2 O 3 composite AEM was 314 kcal/mol more exothermic than that on pristine AEM while chloride adsorption on Mn 2 O 3 composite AEM was 2.2 kcal/mol less exothermic than that on a pristine AEM. Overall, this work presents an effective strategy for selectively removing phosphate from model wastewater solutions and the mechanistic understanding that governs ion selectivity in composite ion-exchange membranes used in MCDI.

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

Absolute photoionization cross sections of brominated organics and determination of an effective C—Br bond cross section

For this work, absolute photoionization cross sections for five brominated species were measured using multiplexed synchrotron photoionization mass spectrometry at the Advanced Light Source (Lawrence Berkeley National Laboratory). These data enable reliable detection of bromine-containing compounds and estimation of C—Br bond photoionization cross sections. Complementary electronic-structure calculations provided adiabatic ionization energies in good agreement with experiment, while thermochemical appearance energies were measured directly. Bond dissociation energies for key fragmentation pathways of the neutral parent molecules were also analyzed. At 11 eV, cross sections were 53.8 ± 7.50 Mb (sp 2 C—Br) and 48.8 ± 7.70 Mb (sp 3 C—Br).

Brominated compounds↗
DOE OSTI2025

Siegert-pseudostate formulation with B-splines

Siegert states (SSs) serve as a useful basis for studying quantum scattering from finite-range potentials. Since they form a discrete instead of continuous set of eigen-solutions, SSs are convenient for performing electronic structure calculations in atoms, molecules, and plasmas. Numerical instabilities may arise, however, in the computation of SSs if the potential vanishes for some extended region, a situation commonly occurring in plasma calculations. Here, in this paper, we identify the cause of these instabilities as the use of non-localized radial basis functions. We thus advocate the use of localized radial basis functions, in particular B-splines, for more robust computations of SSs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Structure and disorder resulting in ultralow thermal conductivity in the defect chalcopyrite AgInSnSe 4

Quaternary chalcogenides are attractive for a variety of technological fields of interest due to the diverse physical properties these materials possess. Herein, the structural and thermal properties of the disordered chalcopyrite AgInSnSe 4 are reported revealing a distinct relation between its structure and physical properties. Modeling of the experimental temperature-dependent thermal properties revealed an ultralow thermal conductivity (0.47 Wm -1 K -1 at room temperature) due to lattice anharmonicity, a low speed of sound and a low Debye temperature. Moreover, first principles electronic structure calculations revealed that the Ag-Se tetrahedra within the crystal structure possess relatively weak bonds due to occupied antibonding states from p-d orbital hybridization, which suppress the thermal conductivity. Here, in order to further quantify our findings, we extended our analyses to include comparisons with other ternary and quaternary adamantine materials, all of which possess a fourfold tetrahedral coordination of atoms, in revealing the origin of the thermal properties in AgInSnSe 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Analysis of correlations between intrinsic ductility and electronic density of states in refractory alloys

High entropy alloys (HEAs) correspond to a new and emerging class of materials that allows us to explore a large composition space to tune mechanical strength and thermal stability. Therefore, to design better alloys, it is important to scan the high-dimensional space of chemistry, composition and temperature. Here, to facilitate this search, we present a method to screen intrinsically ductile body centered cubic (BCC) refractory alloys from electronic structure calculations by using the density of states (DOS) at the Fermi level, g(μ F ). This correlation between intrinsic ductility and g(μ F ) is tested by analyzing group V (V, Nb, Ta) and VI (Mo, W) refractory metals, binary alloys, such as W-Nb, W-V, Mo-Nb and Mo-V, and refractory alloys for which experimental stress-strain measurements are available. In addition, we perform a high-throughput exploration of the entire composition space of a recently proposed alloy system, CrMoNbV, and identify compositions that exhibit high intrinsic ductility.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular Mass Growth Processes to Polycyclic Aromatic Hydrocarbons through Radical–Radical Reactions Exploiting Photoionization Reflectron Time-of-Flight Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) represent critical building blocks in molecular mass growth processes to carbonaceous nanoparticles, referred to as interstellar and circumstellar grains along with soot particles in astrophysical environments and combustion systems, respectively. Recent advancements on elucidating elementary steps to PAHs have utilized reactions of aromatic radicals, resonantly stabilized free radicals, and aliphatic radicals with closed shell hydrocarbons. However, the role of radical–radical reactions (RRRs) leading to PAHs has remained largely unexplored on the molecular level due to preceding experimental challenges in producing sufficiently high number densities of radical reactants for isomer-selective detection of products from bimolecular and termolecular reactions. This Account offers the latest developments in our knowledge on the mechanisms and pathways to PAHs via RRRs probed in a chemical microreactor at temperatures as high as 1600 K. Product preservation in a molecular beam coupled with synchrotron vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry and photoelectron photoion coincidence spectroscopy enabled isomer-selective detection of PAHs of up to three rings by their photoionization efficiency curves, which were fit with a linear combination of reference curves for identification. Experiments were combined with computational fluid dynamics modeling of the physicochemical processes in the microreactor, as well as high-level electronic structure calculations to reveal the reaction pathways of each system. Six distinct reaction mechanisms were discovered in this work: propargyl addition─benzannulation (PABA), methyl addition─ring expansion (MARE), cyclopentadienyl addition─naphthylization (CPAN), fulvenallenyl addition─cyclization─aromatization (FACA), benzyl addition─aromatization (BAA), and phenyl addition─pentacyclization (PAP). By systematically varying the number of carbon atoms in the radical reactants, molecular mass growth processes involving reactions between radicals with odd numbers of carbon atoms access aromatics carrying one, two, or three six-membered rings, whereas reactions between even- and odd-carbon-numbered radicals produce aromatics combining five- and six-membered rings. Our investigations reveal unconventional cycloadditions on excited state triplet surfaces, additions of radicals to low spin density carbon-centered radicals, spiroaromatic and fulvene-type intermediates, and highly strained bicyclic reaction intermediates, challenging current perceptions of PAH molecular mass growth processes. All of the listed mechanisms, except for FACA, feature endoergic reactions or barriers which lie above the separated reactants and therefore might be central to circumstellar environments of carbon-rich stars and planetary nebulae as their descendants, but they play no role in the gas phase of cold molecular clouds where temperatures as low as 10 K dominate. Altogether, this work provides detailed reaction mechanisms of PAH growth processes, advancing our knowledge of the chemistry of carbonaceous matter in the universe.

Addition reactions↗
DOE OSTI2025

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

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

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

Solid-State Hydroxide Ion Conductivity in Silver(I) Oxide, Ag 2 O

Silver (I) oxide, Ag 2 O, precipitated as microcrystals by combining aqueous silver(I) nitrate and KOH solutions, was found to be a solid-state hydroxide ion conductor with ionic conductivity on the order of 10 -3 S/cm. The proton chemical shifts at 4.87 and -7.35 ppm measured by solid-state 1H-NMR experiments are attributed to water molecules and in-lattice OH - coordinated to silver, respectively. The lack of spinning sidebands around the 4.87 ppm peak indicates rapid reorientation on the NMR timescale, suggesting that the water molecules are adsorbed to the surface of the Ag 2 O crystals. Pulsed field gradient measurements gave similar diffusion coefficients (2 x 10 -7 cm 2 /s at 298K) for all three proton environments indicating chemical exchange between sites on the millisecond timescale. The activation energy for OH - diffusion measured by NMR (0.18 eV) was comparable to that obtained by conductivity measurements and density functional theory (DFT) electronic structure calculations. Furthermore, the calculated Pourbaix diagram of Ag 2 O is consistent with the slightly lower sample density observed in He pycnometry and thermogravimetric measurements.

25 ENERGY STORAGE↗
DOE OSTI2026

Mo 4 FeGa 17.25– x Ge x : Complementary Point Substitutions, Buffering Frameworks, and Merging of the 18-n and Octet Bonding Schemes

We present the discovery of Mo 4 FeGa 17.25–x Ge x (x ∼ 2.1(4), based on determination of Ge content), a complex gallide featuring a variety of point substitution phenomena. Its crystal structure is derived from the Ti 2 Ni type, in which a diamond network of face-sharing octahedra is interpenetrated by a second diamond network of vertex-sharing stella quadrangula. However, while in the ideal Ti 2 Ni type these two frameworks run uninterrupted through the crystal, Mo 4 FeGa 17.25–x Ge x shows three variations. First, the inner tetrahedron of every other stella quadrangula is replaced with a main group atom, creating tetrahedra reminiscent of the Zintl phase NaTl. Next, a selection of Ga 6 octahedra are filled with Fe atoms, in a manner analogous to the stuffed AuCu 3 -type phases. Finally, the refined crystal structure shows that in each unit cell a single Ga2 atom in the octahedral network is substituted with a dumbbell of Ga/Ge atoms. Electronic structure calculations on ordered models of Mo 4 FeGa 17.25–x Ge x reveal a narrow band near the Fermi energy, which is explained in terms of the 18-n and octet bonding schemes using reversed approximation Molecular Orbital analysis. A DFT-chemical pressure analysis connects the tetrahedron/atom and atom/dumbbell substitutions to the relief of atomic packing tensions and highlights soft atomic motions within the octahedral framework and driving forces for atom/dumbbell substitution. The combination of soft vibrational modes, disorder, and a narrow band gap could make this phase of interest for potential thermoelectric properties.

Chemical structure↗
DOE OSTI2025

Affinity for OH – Produces Four-Coordinated Zn 2+ Impurities in Hydrated Amorphous Calcium Carbonate

Using ab initio based molecular dynamics and electronic structure calculations, we show that Zn impurities in hydrated amorphous calcium carbonate (ACC) have a much lower coordination number than other divalent impurities due to covalent interactions between the 3d Zn shell and the oxygen atoms of the carbonate and water groups. Further, the local structure around Zn in ACC, including the predicted low coordination number, is confirmed by X-ray absorption spectroscopy of synthetic Zn-bearing ACC. The strong Zn–O chemical interaction leads to substantial water dissociation and slightly disrupts the hydrogen bonding network. Implications of Zn 2+ incorporation for ACC stability are discussed.

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

Understanding the Hydronium Cation in the Solid State: A Study in Synthetic Hydronium Uranyl Phosphate and Arsenate Mineral Systems and Their Irradiation Stability

Here, we present neutron diffraction, electronic structure calculations, and optical vibrational spectroscopic characterization of hydronium uranyl phosphate (HUP), the synthetic analog of chernikovite ((H 3 O)(UO 2 )(PO 4 )·3H 2 O) to gain insight into structural features of the hydronium cation in the solid state. HUP crystallizes in P2 1 /c instead of previously reported P4/ncc 1 at room temperature and Pccn below 302 K. Lower symmetry was required due to hydronium and interstitial water positions within the structure. Crystallographic positions of hydronium cations determined from powder neutron diffraction data are bolstered by Raman and attenuated total reflectance infrared spectroscopic measurements and are further informed by density functional theory with phonon eigenvector analysis for spectral assignments. Finally, HUP and its arsenate analog (trogerite, (H 3 O)(UO 2 )(AsO 4 )·3H 2 O), (HUAs) were studied using He 2+ irradiation as an analog for α radiolysis to investigate the irradiation stability of these phases and the stability of hydronium cations in the solid state. Dose studies were employed wherein each sample was irradiated to 5, 10, 15, 25, and 50 MGy. Structural insight regarding irradiated materials is gained using optical vibrational spectroscopy and powder X-ray diffraction. Surprising irradiation stability of HUAs was found up to 50 MGy of dose, which may have important implications for understanding and modeling the geologic stability of legacy U.

actinides↗
DOE OSTI2025

Beyond CCSD(T) Accuracy at Lower Scaling with Auxiliary Field Quantum Monte Carlo

We introduce a black-box auxiliary field quantum Monte Carlo (AFQMC) approach to perform highly accurate electronic structure calculations using configuration interaction singles and doubles (CISD) trial states. This method consistently provides more accurate energy estimates than coupled cluster singles and doubles with perturbative triples (CCSD(T)), often regarded as the gold standard in quantum chemistry. This level of precision is achieved at a lower asymptotic computational cost, scaling as O(N 6 ) compared to the O(N 7 ) scaling of CCSD(T). Furthermore, we provide numerical evidence supporting these findings through results for challenging main group and transition metal-containing molecules.

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

Accelerating Instanton Theory with the Line Integral Nudged Elastic Band Method and Gaussian Process Regression

Quantum tunneling plays a fundamental role in many chemical reactions, particularly proton transfer processes. Ring polymer instanton theory offers a practical framework for computing tunneling rates in complex molecular systems. However, applying the ring polymer instanton method with a potential energy surface generated on-the-fly using electronic structure calculations can be computationally demanding. Here, in this work, we present a new efficient implementation of the ring polymer instanton method by combining the Line Integral Nudged Elastic Band (LI-NEB) approach with Gaussian Process Regression (GPR). We benchmarked this method on prototypical ground-state proton transfer systems, including the benchmark gas-phase hydrogen abstraction reaction H + CH 4 → H 2 + CH 3 , malonaldehyde, and Z-3-amino-propenal (aminopropenal). Our results show that this approach is an order of magnitude faster than traditional instanton algorithms while maintaining excellent agreement with their tunneling rates. This development opens the door to studying proton transfer in larger systems with improved efficiency.

chemical physics↗
DOE OSTI2025

Energetic and Electronic Properties of AcX and LaX (X = O and F)

The bonding and spectroscopic properties of LaX and AcX (X = O and F) diatomic molecules were studied by highlevel ab initio CCSD(T) and SO-CASPT2 electronic structure calculations. Bond dissociation energies (BDEs) were calculated at the Feller−Peterson−Dixon (FPD) level. Potential energy curves and spectroscopic constants for the lowest-lying spin−orbit Ω states were obtained at the SO-CASPT2/aQ-DK level. A dense manifold of excited states was described for the monofluorides with the ground states well separated from the excited states. The spectroscopic parameters were in good agreement with those reported experimentally for LaO and LaF. For the diatomic molecules containing actinides, no experimental data of these parameters was found, but the results were consistent with other high-level calculations. The BDEs calculated at the FPD level were 791.3 (LaO), 705.2 (AcO), 650.0 (LaF), and 678.6 (AcF) kJ/ mol. The NBO analysis showed that the monofluorides are essentially ionic, which explains why the BDE(AcF) is higher than BDE(LaF); for the monoxides, covalent contributions involving the d orbitals of the metal and the p orbitals of the oxygen are stronger for LaO than AcO, which explains the higher BDE for LaO. Finally, the bond orders are predicted to be 2 for LaF and AcF, 3 for AcO, and higher than 3 for LaO.

Basis sets↗
DOE OSTI2025

Formation of Water Networks on Anionic Perylene

We present infrared photodissociation spectra of hydrated perylene anion clusters with up to four water molecules, as well as electronic structure calculations based on density functional theory. Water molecules form weak hydrogen bonds to the π system of the perylene anion. For clusters with more than one water molecule, water−water hydrogen bonds are formed, which generally appear to be stronger than water−π hydrogen bonds, especially for the trihydrate and tetrahydrate. The resulting water networks exist as water subclusters on the surface of the carbon frame of perylene. Here, we observe temperature-dependent dynamic effects, which highlight large amplitude motions of the water network and the shallowness of the potential energy surfaces governing the structures of these clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗