Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular geometry”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Validation of ablation model for polyethylene using pulsed x-ray and proton exposures

The surface erosion of polyethylene is interrogated using pulsed x rays at the Z Machine (Sandia National Laboratories) and with proton beams at the Gamble II generator (Naval Research Laboratory) to validate a coupled model for volumetric thermal ablation, photoionization, finite-rate decomposition, and molecular recombination of radicals. The intense radiation pulses (up to ∼1014W/m2 over tens of nanoseconds) are used to generate one-dimensional vapor flows with low ionization fractions and a simplified geometry compared to typical laser ablation, allowing for evaluation of the model under local thermal equilibrium conditions. Areal momentum carried by the ensuing uniaxial hydrodynamic shock is used to indicate the extent of ablation. The threshold fluence for ablation is found to be in close correspondence with the bulk melt transition, and reasonable agreement with the model is obtained for peak temperatures in polyethylene up to 5500 K and heating rates up to 1011K/s where thermal decomposition reactions are also active.

Moore, N. W. (ORCID:0000000272259674)↗
DOE OSTI2020

Materials Data on H7ClO3 by Materials Project

H5O2H7O3ClH2OCl crystallizes in the triclinic P1 space group. The structure is zero-dimensional and consists of one hydrochloric acid, monohydrate molecule; one molecular hydrogen;dihydrate molecule; and one H7O3Cl cluster. In the H7O3Cl cluster, there are seven inequivalent H1+ sites. In the first H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 1.00 Å. In the second H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the third H1+ site, H1+ is bonded in a distorted single-bond geometry to one O2- and one Cl1- atom. The H–O bond length is 1.03 Å. The H–Cl bond length is 1.88 Å. In the fourth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. In the fifth H1+ site, H1+ is bonded in a distorted linear geometry to two O2- atoms. There is one shorter (1.02 Å) and one longer (1.61 Å) H–O bond length. In the sixth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 1.00 Å. In the seventh H1+ site, H1+ is bonded in a linear geometry to two O2- atoms. There is one shorter (1.14 Å) and one longer (1.29 Å) H–O bond length. There are three inequivalent O2- sites. In the first O2- site, O2- is bonded in a trigonal non-coplanar geometry to three H1+ atoms. In the second O2- site, O2- is bonded in a distorted water-like geometry to three H1+ atoms. In the third O2- site, O2- is bonded in a trigonal non-coplanar geometry to three H1+ atoms. Cl1- is bonded in a single-bond geometry to one H1+ atom.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on H5ClO6 by Materials Project

H5O2ClO4 crystallizes in the orthorhombic P2_12_12_1 space group. The structure is zero-dimensional and consists of four molecular hydrogen;dihydrate molecules and four ClO4 clusters. In each ClO4 cluster, there are four inequivalent O sites. In the first O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. In the second O site, O is bonded in a distorted single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. In the third O site, O is bonded in a distorted single-bond geometry to one Cl atom. The O–Cl bond length is 1.49 Å. In the fourth O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.44 Å. Cl is bonded in a tetrahedral geometry to four O atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on H5ClO6 by Materials Project

H5O2ClO4 crystallizes in the orthorhombic P2_12_12_1 space group. The structure is zero-dimensional and consists of four molecular hydrogen;dihydrate molecules and four ClO4 clusters. In each ClO4 cluster, there are three inequivalent O sites. In the first O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.47 Å. In the second O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.45 Å. In the third O site, O is bonded in a single-bond geometry to one Cl atom. The O–Cl bond length is 1.46 Å. Cl is bonded in a tetrahedral geometry to four O atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

One-colour (~220 nm) resonance-enhanced (S 1 –S 0 ) multi-photon dissociation of acetylene: probe of the C 2 A 1 Π u – X 1 Σ + g band by frequency-modulation spectroscopy

In a recent paper, we demonstrated that one-colour (~220 nm), resonance-enhanced (S 1 –S 0 ), photodissociation of acetylene generates strong $C_2$ Swan band ($d^3Π_g–a^3Π_u$) and $C_2$ Deslandres-d'Azambuja band ($C^1Π_g–A^1Π_u$) fluorescence, and long-lived (>3 µs) fluorescence from an electronically-excited $C_2$H * species. It was not known whether the $C_2A^1Π_u$ and $X^1Σ^+_g$ states are also directly populated in this process. In this paper, multiple vibration-rotation transitions between the $C_2A$-state v = 2 and the X-state v = 0 level are examined by time-resolved frequency-modulation (FM) spectroscopy. The photolysis laser wavelength is tuned into resonance at the one-photon level with S 1 –S 0 transitions that populate individual rotational levels of the S 1 $_{trans}$-conformer $3^2,3^3,$ and $3^4$ vibrational states. By comparing the phase of the FM signals from the $C_2A–X$ transitions with that from the Rb D 1 -line absorption transition, we determine that, for all of the probed A–X transitions, the X-state level is more populated than the A-state level. We propose that the acetylene S 1 level is excited by the second photon to an acetylene dissociation precursor state, which undergoes sequential C-H bond-breaking to produce the $C_2X$ state. The dissociation precursor is assigned as the $1^1B_g(trans-bent)/1^1B_1(cis-bent)$ valence state, which correlates to a doubly-excited configuration, $(1π_u)^2(1π_g)^2$, at linear geometry. Based on the rotational distributions of the $C_2X$-state fragments, we believe that at least one of the transition states involved in the photolysis via S 1 $3^4$ has a larger CC-H bend-angle for the departing H-atom than that involved in the S 1 $3^2$ and $3^3$ photolysis.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Phosphatidylinositol phosphates modulate interactions between the StarD4 sterol trafficking protein and lipid membranes

There is substantial evidence for extensive nonvesicular sterol transport in cells. For example, lipid transfer by the steroidogenic acute regulator-related proteins (StarD) containing a StarT domain has been shown to involve several pathways of nonvesicular trafficking. Among the soluble StarT domain–containing proteins, StarD4 is expressed in most tissues and has been shown to be an effective sterol transfer protein. However, it was unclear whether the lipid composition of donor or acceptor membranes played a role in modulating StarD4-mediated transport. Here, we used fluorescence-based assays to demonstrate a phosphatidylinositol phosphate (PIP)-selective mechanism by which StarD4 can preferentially extract sterol from liposome membranes containing certain PIPs (especially, PI(4,5)P 2 and to a lesser degree PI(3,5)P2). Monophosphorylated PIPs and other anionic lipids had a smaller effect on sterol transport. This enhancement of transport was less effective when the same PIPs were present in the acceptor membranes. Furthermore, using molecular dynamics (MD) simulations, we mapped the key interaction sites of StarD4 with PIP-containing membranes and identified residues that are important for this interaction and for accelerated sterol transport activity. We show that StarD4 recognizes membrane-specific PIPs through specific interaction with the geometry of the PIP headgroup as well as the surrounding membrane environment. Finally, we also observed that StarD4 can deform membranes upon longer incubations. Taken together, these results suggest a mechanism by which PIPs modulate cholesterol transfer activity via StarD4.

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

Designing molecular qubits: computational insights into first-row and group 6 transition metal complexes

In the realm of optically addressable qubits, a previously synthesized and characterized Cr( IV ) pseudo-tetrahedral complex, featuring four strongly donating ligands surrounding the chromium center, has demonstrated potential as a qubit candidate. This study proposes analogs of this complex through a metal substitution strategy, extending the investigation to different complexes based on metal centers selected from first-row and Group 6 transition metals. Computational modeling based on multiconfigurational methods CASPT2 and MC-PDFT was utilized to calculate energy gaps between ground and excited electronic spin states, and zero-field splitting parameters. Simulations were applied to each equilibrium geometry and related deformations based on vibrational modes. All results align with previous experimental findings, but also show that qubits based on V and Ti centers could be more electronically stable than the Cr one, suggesting a lower electronic features dependency from their related geometry. In some cases geometrical deformations provide changes in relative energy gaps between triplet and singlet excited state, that could potentially swap, offering a different initialization process, and some inspiration for ligand design based on such deformations. Additionally, this study identifies an unsynthesized Ti( II ) compound as a promising candidate for molecular qubits. This finding highlights the role of computational multireference methods in the rational design of qubit systems.

Sauza-de la Vega, Arturo [Univ. of Chicago, IL (Un↗
DOE OSTI2025

Molecular Design of Al(II) Intermediates for Small Molecule Activation

Promoting societally important small molecule activation processes with earth-abundant metals is foundational for a sustainable chemistry future. In this context, mapping new reaction pathways that would enable abundant main-group elements to mimic the behaviors of d- and f-block elements is facilitated by exploring unusual oxidation states. The most abundant metal on earth, aluminum, has been well studied in the Lewis acidic +III and Lewis basic +I oxidation states but rarely in the potentially biphilic +II oxidation state until recently, when a renaissance of Al(II) chemistry emerged from a range of research groups. In this Perspective, we review the chemistry of mononuclear Al radicals, including both Al-centered radicals (i.e., Al(II) compounds) and redox non-innocent systems (i.e., formally Al(II) species that are physically Al(III) with ligand-centered radicals), with an emphasis on small molecule reactivity. We also provide a meta-analysis of the Al(II) literature to summarize how different design strategies (e.g., redox non-innocence, strained coordination geometries) have been shown to impart biphilic character to Al radicals and tune their behavior, thus allowing Al radicals to mimic the chemistry of certain d- and f-block metal ions such as Ti(III) and Sm(II). We expect these molecular design concepts to inform future Al(II) studies as the chemistry of this unusual oxidation state of Al continues to grow.

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

Identifying crack tip position and stress intensity factors from displacement data

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, in this work, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lowering the Activation Barriers for Lithium-Ion Conductivity through Orientational Disorder in the Cyanide Argyrodite Li6PS5CN

Rapid advancements in safe and high-energy-density energy storage are predicated on identifying new solid-state ion conductors with low activation energies and high ionic conductivities for all-solid-state battery technologies. Halide argyrodites are among some of the top candidates for solid-state electrolytes, as they can achieve ionic conductivities that approach liquid electrolytes. Incorporating dynamic pseudohalide species in argyrodite solid electrolytes presents an exciting opportunity to exploit lattice dynamics as a design principle to modulate the ion conduction properties of solid-state ion conductors. In the present study, we have prepared the new argyrodite Li 6 PS 5 CN containing orientationally disordered cyanide ions. The new cyanide argyrodite Li 6 PS 5 CN exhibits an activation barrier to Li-ion transport of 471 ± 25 meV and a room-temperature ionic conductivity of 6(2) × 10 –5 S cm –1 in comparison to the activation barrier of 502 ± 16 meV and an ionic conductivity of 2.3(1) × 10 –4 S cm –1 measured for the bromide analogue Li 6 PS 5 Br. Structural studies of both compounds by high-resolution X-ray diffraction indicate that Li 6 PS 5 CN and Li 6 PS 5 Br adopt nearly identical crystal structures with similar lattice parameters, which indicates that lower activation barriers in Li 6 PS 5 CN arise due to the cyanide ion itself rather than due to changes in the geometry of conduction pathways in the local lithium environment. The orientational disorder of the quadrupolar cyanide ion in Li 6 PS 5 CN points to a complex interplay of lattice polarizability and molecular dynamics that lower the activation barrier for lithium-ion conductivity in the cyanide argyrodite.

25 ENERGY STORAGE↗
DOE OSTI2024

Ligand Many-Body Expansion as a General Approach for Accelerating Transition Metal Complex Discovery

Methods that accelerate the evaluation of molecular properties are essential for chemical discovery. While some degree of ligand additivity has been established for transition metal complexes, it is underutilized in asymmetric complexes, such as the square pyramidal coordination geometries highly relevant to catalysis. To develop predictive methods beyond simple additivity, we apply a many-body expansion to octahedral and square pyramidal complexes and introduce a correction based on adjacent ligands (i.e., the cis interaction model). We first test the cis interaction model on adiabatic spin-splitting energies of octahedral Fe(II) complexes, predicting DFT-calculated values of unseen binary complexes to within an average of 1.4 kcal/mol. Uncertainty analysis reveals the optimal basis, comprising the homoleptic and mer symmetric complexes. We next show that the cis model (i.e., the cis interaction model solved for the optimal basis) infers both DFT- and CCSD(T)-calculated model catalytic reaction energies to within 1 kcal/mol on average. The cis model predicts low-symmetry complexes with reaction energies outside the range of binary complex reaction energies. We observe that trans interactions are unnecessary for most monodentate systems but can be important for some combinations of ligands, such as complexes containing a mixture of bidentate and monodentate ligands. Lastly, we demonstrate that the cis model may be combined with Δ-learning to predict CCSD(T) reaction energies from exhaustively calculated DFT reaction energies and the same fraction of CCSD(T) reaction energies needed for the cis model, achieving around 30% of the error from using the CCSD(T) reaction energies in the cis model alone.

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

Analysis of the bonding in tetrahedrane and phosphorus-substituted tetrahedranes

The bonding structures of tetrahedrane, phosphatetrahedrane, diphosphatetrahedrane and triphosphatetrahedrane are studied by employing an intrinsic quasi-atomic orbital analysis. Ethane, cyclopropane and tetrahedral P 4 are employed as reference systems. The orbital analysis is paired with the computation of strain energies via isodesmic reactions. The results show that the increase in geometric strain upon transition from ethane to cyclopropane to tetrahedrane weakens the CC bonds, despite leading to shorter C–C interatomic distances. With the increase in strain, the orbitals centered on C and involved in the bonding of the cage structure are observed to have elevated p-character, and the orbital structure of C deviates from sp 3 hybridization. The systematic substitution of CH groups by P atoms in the cage structure of tetrahedrane leads to stronger CC bonds, larger angles in the cage structures of the resulting phosphatetrahedranes, lower p-character in the orbitals involved in the bonding of the cages, and lower strain energies. It is found that P is more amenable to strained molecular arrangements than is C, and that the propensity of a given atom to hybridize s and p functions, or the lack thereof, has implications in the stability of molecules with strained geometries. The combination of the calculations presented here with the existing literature provides insight into the apparent propensity of tetrahedrane and P 4 to ‘break’ their tetrahedral structures. Trends in the bonding interactions, such as bond strengths, s- and p-orbital characters and charge transfer are identified and related to the strain energy observed in each of the analyzed systems.

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

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA’s diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Materials Data on Zn2Hg2H3(SO5)2 by Materials Project

(HgZnSO5)4(H2)3 crystallizes in the orthorhombic Pbcm space group. The structure is three-dimensional and consists of four molecular hydrogen molecules and one HgZnSO5 framework. In the HgZnSO5 framework, there are two inequivalent Hg+1.50+ sites. In the first Hg+1.50+ site, Hg+1.50+ is bonded in a 2-coordinate geometry to six O2- atoms. There are a spread of Hg–O bond distances ranging from 2.10–2.97 Å. In the second Hg+1.50+ site, Hg+1.50+ is bonded in a distorted linear geometry to two O2- atoms. There are one shorter (2.07 Å) and one longer (2.10 Å) Hg–O bond lengths. Zn2+ is bonded to five O2- atoms to form ZnO5 trigonal bipyramids that share corners with three equivalent SO4 tetrahedra and an edgeedge with one ZnO5 trigonal bipyramid. There are a spread of Zn–O bond distances ranging from 2.01–2.21 Å. S5+ is bonded to four O2- atoms to form SO4 tetrahedra that share corners with three equivalent ZnO5 trigonal bipyramids. There are a spread of S–O bond distances ranging from 1.46–1.52 Å. There are six inequivalent O2- sites. In the first O2- site, O2- is bonded in a 2-coordinate geometry to one Hg+1.50+, one Zn2+, and one S5+ atom. In the second O2- site, O2- is bonded to two Hg+1.50+ and two equivalent Zn2+ atoms to form a mixture of corner and edge-sharing OZn2Hg2 tetrahedra. In the third O2- site, O2- is bonded in a bent 150 degrees geometry to one Zn2+ and one S5+ atom. In the fourth O2- site, O2- is bonded to two Hg+1.50+ and two equivalent Zn2+ atoms to form a mixture of corner and edge-sharing OZn2Hg2 tetrahedra. In the fifth O2- site, O2- is bonded in a distorted bent 150 degrees geometry to one Zn2+ and one S5+ atom. In the sixth O2- site, O2- is bonded in a single-bond geometry to one Hg+1.50+ and one S5+ atom.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Real-Time Monitoring of Strain Accumulation and Relief during Epitaxy of Ultrathin Co Ferrite Films with Varied Co Content

Ultrathin Co x Fe 3-x O 4 films of high structural quality and with different Co content (x = 0.6–1.2) were prepared by reactive molecular beam epitaxy on MgO(001) substrates. Epitaxy of these ferrite films is extensively monitored by means of time-resolved (operando) X-ray diffraction recorded in out-of-plane geometry to characterize the temporal evolution of the film structure. The Co ferrite films show high crystalline ordering and smooth film interfaces independent of their Co content. All Co x Fe 3-x O 4 films exhibit enhanced compressive out-of-plane strain during the early stages of growth, which partly releases with increasing film thickness. When the Co content of the ferrite films increases, the vertical-layer distances increase, accompanied by slightly increasing film roughnesses. The latter result is supported by surface-sensitive low-energy electron diffraction as well as X-ray reflectivity measurements on the final films. In contrast, the substrate–film interface roughness decreases with increasing Co content, which is confirmed with X-ray reflectivity measurements. In addition, the composition and electronic structure of the ferrite films is characterized by means of hard X-ray photoelectron spectroscopy performed after film growth. The experiments reveal the expected increasing Fe 3+ /Fe 2+ cation ratios for a higher Co content.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab-Initio Simulation Studies of Cr Solvation in Fluoride Molten Salts

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF 2 , and 3LiF-AlF3. We find that Cr 0 ,Cr 2+ ,Cr 3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF 2 and 3LiF-AlF 3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent.

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

Defect-Free Growth of Decagonal Quasicrystals around Obstacles

Quasicrystals are aperiodic solids with exotic physical properties. The mechanisms driving their growth are far from understood, particularly in the presence of rigid obstacles. Here, using in situ synchrotron x-ray tomography and molecular dynamics simulations, we investigate the interaction between decagonal quasicrystals and obstacles (shrinkage pores) during solidification of an Al 79 ⁢Co 6 ⁢Ni 15 alloy. The results reveal defect-free quasicrystal growth around the pores, independent of pore size and geometry. Our findings point to a universal feature of aperiodic solids: their ability to accommodate structural disruptions due to additional (phasonic) degrees of freedom. This Letter opens new possibilities for synthesizing large-scale, single quasicrystals from a liquid metal precursor.

alloys↗
DOE OSTI2026

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗