Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave function”

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

DFT Analysis of the Binding of Rare Earth Nitrates at Internal and External Surfaces of MCM-22

MCM-22-type zeolites constitute a well-characterized tunable class of aluminosilicates suitable for elucidating the fundamental aspects of the binding of rare earth elements (REEs) in layered materials. Starting from the pure silica version ITQ-1 with a unit cell of Si 72 O 144 , a model for periodic bulk crystalline MCM-22 with a finite Al concentration is provided by replacing a Si atom with an Al atom in the unit cell at suitable tetrahedral sites near an internal pore surface. Then, a H atom is added to an O atom bridging Si and Al atoms to create a Brønsted acid site (BAS). There are no internal silanol groups in this bulk model. To generate a model for an external surface, we adopt the fully hydroxylated surface structure of a layer within the ITQ-1 precursor with two silanols per lateral unit cell. A BAS on the external surface can be generated by replacing a near-surface Si atom with an Al atom and adding a H atom, as above. The strength of binding at a BAS of REE, X, taken to be present in the solution phase as nitrates, is determined from the energy change in the reaction X(NO 3 ) 3 + ≡Si–{OH}–Al≡ → ≡Si–{OX(NO 3 ) 2 }–Al≡+ HNO 3 . The strength of binding at the silanols is determined similarly. Binding energies are determined from two approaches. The first performs periodic plane-wave density functional theory (DFT) total energy analysis for an entire unit cell of MCM-22. The second utilizes cluster models capturing the local environment of REE binding sites and performs DFT analysis with localized basis sets. The two approaches yield consistent results for Nd, revealing similarly strong binding at either an external or internal BAS, but much weaker binding at a silanol site. This is consistent with the picture deduced from recent experiments. Here, we also comment on binding at Al-bridged siloxane sites, which have been suggested as alternative binding sites to BAS.

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

Reducing Qubit Requirements while Maintaining Numerical Precision for the Variational Quantum Eigensolver: A Basis-Set-Free Approach

We present a basis-set-free approach to the variational quantum eigensolver using an adaptive representation of the spatial part of molecular wave functions.Our approach directly determines system-specific representations of qubit Hamiltonians while fully omitting globally defined basis sets. In this work, we use directly determined pair-natural orbitals on the level of second-order perturbation theory. This results in compact qubit Hamiltonians with high numerical accuracy. We demonstrate initial applications with compact Hamiltonians on up to 22 qubits where conventional representation would for the same systems require 40–100 or more qubits. Here, we further demonstrate reductions in the quantum circuits through the structure ofthe pair-natural orbitals.

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

Rotational Modulation of à 2 A State Photodissociation of HCO via Renner-Teller Nonadiabatic Transitions

By examining the product-state distribution of a prototypical non- adiabatic predissociation system, HCO($Ã^2A''–\tilde X^2A'$), we demonstrate here that the dissociation dynamics is strongly modulated by parent rotational quantum numbers. The predissociation of the nominal ($ν_{\text{C–H}} = 0, ν_{\text{bend}}, ν_{\text{C–O}} = 1$) vibronic levels of the $AÃ^2A''$ state surprisingly gives rise to both vibrational ground and excited states of the CO product, despite the assumed spectator nature of the CO moiety. This anomaly is attributed to the dependence of the lifetime of the vibronic resonance facilitated by the Renner–Teller interaction on the parent rotational angular momentum quantum numbers coupled with transient intensity borrowing from nearby vibronic resonances with $ν_{\text{C–O}} = 0$. This unique phenomenon is a purely quantum mechanical behavior that has no classical analogue.

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

Is the Matrix Completion of Reduced Density Matrices Unique?

Reduced density matrices are central to describing observables in many-body quantum systems. In electronic structure theory, the two-particle reduced density matrix (2-RDM) suffices to determine the energy and other key properties. Recent work has used matrix completion, leveraging the low-rank structure of RDMs and approximate theoretical models, to reconstruct the 2-RDM from partial data and thus reduce the computational cost. However, matrix completion is, in general, an under-determined problem. Revisiting Rosina’s theorem (Rosina, M. Queen’s Papers on Pure and Applied Mathematics , 1968, No. 11, 369), we here show that the matrix completion is unique under certain conditions, identifying the subset of 2-RDM elements that enables its exact reconstruction from incomplete information. Building on this, we introduce a hybrid quantum–stochastic algorithm that achieves exact matrix completion, demonstrated through applications to the Fermi–Hubbard model.

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

Pt-Assisted Carbon Remediation of Mo2C Materials for CO Disproportionation

Using the CO disproportionation (Boudouard) reaction as a probe reaction, an in-depth analysis of temperature-programmed pulse response data shows that the addition of Pt to Mo2C mitigates deactivation of Mo active sites by acting as a carbon collector. CO2 production on Mo2C and Pt/Mo2C materials is dependent on both the activation energy and the CO surface concentration. Detailed plane-wave density functional theory calculations of the CO adsorption and disproportion reactions on Mo2C-supported Pt nanoparticles (NPs) are reported. The Mo2C was modeled by the ß-Mo2C (100) surface, and the Pt/Mo2C interface was modeled by the addition of 12 Pt atoms to the Mo2C (100) surface (12Pt@Mo2C). The potential energy surfaces of the Boudouard reaction were calculated on pure Mo2C, 12Pt@Mo2C, and Pt (111) surfaces. CO dissociation readily occurs on the Mo2C (100) surface, but not on the Pt (111) surface, with the former being exothermic and the latter being endothermic. At the Pt/Mo2C interface, CO dissociation is still exothermic, but with a larger energy barrier. The Boudouard reaction takes place on the Mo2C region, where CO2 is formed from a surface O atom dissociated from one CO molecule in reaction with another CO molecule, leaving one C atom on the surface. C adsorption is preferential on the Pt site in comparison to the Mo site. The supported Pt domains can collect the remaining C atoms, facilitating further CO2 formation on the active Mo sites. A Bader charge analysis shows that the surface metal-carbon bond is a mixture of covalent and ionic bonds, whereas the surface metal-oxygen bond is ionic. Electron localization function (ELF) and partial charge density calculations agree well with the Bader charge analysis. These computational results are consistent with experimental observations of the interaction of CO with Mo2C nanotube supported Pt domains in the transient regime under far from equilibrium conditions. The Boudouard reaction is an important side reaction, and the unexpected role found for Pt as a carbon collector, with Mo serving as a disproportionation site, provides a unique vantage point for understanding carbon and coke formation on catalytic materials.

Pt-assisted carbon remediation, Pt/Mo2C interface,↗
DOE OSTI2023

Synthesis of Quantum-Confined Borophene Nanoribbons

Here, borophene nanoribbons (BNRs) are one-dimensional strips of atomically thin boron expected to exhibit quantum-confined electronic properties that are not present in extended two-dimensional borophene. While the parent material borophene has been experimentally shown to possess anisotropic metallicity and diverse polymorphic structures, the atomically precise synthesis of nanometer-wide BNRs has not yet been achieved. Here, we demonstrate the synthesis of multiple BNR polymorphs with well-defined edge configurations within the nanometer-scale terraces of vicinal Ag(977). Through atomic-scale imaging, spectroscopy, and first-principles calculations, the synthesized BNR polymorphs are characterized and found to possess distinct edge structures and electronic properties. For single-phase BNRs, v 1/6 -BNRs and v 1/5 -BNRs adopt reconstructed armchair edges and sawtooth edges, respectively. In addition, the electronic properties of single-phase v 1/6 -BNRs and v 1/5 -BNRs are dominated by Friedel oscillations and striped moiré patterns, respectively. On the other hand, mixed-phase BNRs possess quantum-confined states with increasing nodes in the electronic density of states at elevated biases. Overall, the high degree of polymorphism and diverse edge topologies in borophene nanoribbons provide a rich quantum platform for studying one-dimensional electronic states.

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

Coupling of Nondegenerate Topological Modes in Nitrogen Core-Doped Graphene Nanoribbons

Nitrogen core-doping of graphene nanoribbons (GNRs) allows trigonal planar carbon atoms along the backbone of GNRs to be substituted by higher-valency nitrogen atoms. The excess valence electrons are injected into the π-orbital system of the GNR, thereby changing not only its electronic occupation but also its topological properties. We have observed this topological change by synthesizing dilute nitrogen core-doped armchair GNRs with a width of five atoms (N 2 -5-AGNRs). The incorporation of pairs of trigonal planar nitrogen atoms results in the emergence of topological boundary states at the interface between doped and undoped segments of the GNR. These topological boundary states are offset in energy by approximately ΔE = 300 meV relative to the topological end states at the termini of finite 5-AGNRs. Scanning tunneling microscopy (STM) and spectroscopy (STS) reveal that for finite GNRs the two types of topological states can interact through a linear combination of orbitals, resulting in a pair of asymmetric hybridized states. This behavior is captured by an effective Hamiltonian of nondegenerate diatomic molecules, where the analogous interatomic hybridization interaction strength is tuned by the distance between GNR topological modes.

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

Toward an Understanding of Linear Scaling Relations through Energy Decomposition Analysis

The discovery of linear scaling relations has fundamentally changed the field of heterogeneous catalysis. The scaling relations have been rationalized based on the d-band theory, specifically a separation of sp and d electron contributions to adsorption energies. Within the framework of energy decomposition analysis, a full understanding of such a separation would require one to further break down the adsorption energy into distinct energy components such as electrostatics, polarization, charge transfer, and van der Waals interactions, and to examine the sp and d contributions to each of them. As a step in this direction, we analyzed the interaction energy between CH x (x = 1–4) adsorbates and fcc(100) transition metal surfaces (M = Cu, Ag, Au, Rh, and Pt), with the surfaces represented both as slabs in plane-wave density functional theory (pw-DFT) calculations and as atomic clusters in atomic-orbital basis density functional theory (ao-DFT) calculations. Through an absolutely localized molecular orbital (ALMO) based energy decomposition analysis of the ao-DFT adsorption energy, each of the interaction energy components (electrostatics, polarization, van der Waals, and charge transfer) was found to follow its own scaling relations, with an intricate interplay among these energy components yielding the overall scaling relations for the total adsorption energies. Using the recently introduced ALMO-based polarization and charge-transfer analysis schemes, we further dissected polarization into metal surface and adsorbate contributions, and charge transfer into metal → adsorbate and adsorbate → metal contributions. The contributions from the sp and d electrons of the metal to these terms were further quantified, and the dominant role of the metal d electrons was reaffirmed. These results shed light on how CHx adsorbates interact with metal surfaces and further reveal the physical origin of the scaling relations.

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

In Situ Velocity‐Strain Sensitivity Near the San Jacinto Fault Zone Analyzed Through Train Tremors

Abstract We utilize train tremors as P‐wave seismic sources to investigate velocity‐strain sensitivity near the San Jacinto Fault Zone. A dense nodal array deployed at the Piñon Flat Observatory is used to detect and identify repeating train energy emitted from a railway in the Coachella valley. We construct P‐wave correlation functions across the fault zone and estimate the spatially averaged dt / t versus strain sensitivity to be 6.25 × 10 4 . Through numerical simulations, we explore how the sensitivity decays exponentially with depth. The optimal solution reveals a subsurface sensitivity of 1.2 × 10 5 and a depth decay rate of 0.05 km −1 . This sensitivity aligns with previous findings but is toward the higher end, likely due to the fractured fault‐zone rocks. The depth decay rate, previously unreported, is notably smaller than assumed in empirical models. This raises the necessity of further investigations of this parameter, which is crucial to study stress and velocity variations at seismogenic depth.

Geology↗
DOE OSTI2023

An open database of computed bulk ternary transition metal dichalcogenides

Abstract We present a dataset of structural relaxations of bulk ternary transition metal dichalcogenides (TMDs) computed via plane-wave density functional theory (DFT). We examined combinations of up to two chalcogenides with seven transition metals from groups 4–6 in octahedral (1T) or trigonal prismatic (2H) coordination. The full dataset consists of 672 unique stoichiometries, with a total of 50,337 individual configurations generated during structural relaxation. Our motivations for building this dataset are (1) to develop a training set for the generation of machine and deep learning models and (2) to obtain structural minima over a range of stoichiometries to support future electronic analyses. We provide the dataset as individual VASP xml files as well as all configurations encountered during relaxations collated into an ASE database with the corresponding total energy and atomic forces. In this report, we discuss the dataset in more detail and highlight interesting structural and electronic features of the relaxed structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

CASPT2 molecular geometries of Fe( II ) spin-crossover complexes

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. Here we recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe–ligand bond distances.

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

Probing growth of metal–organic frameworks with X-ray scattering and vibrational spectroscopy

Nucleation and crystallization arising from liquid to solid phase are involved in a multitude of processes in fields ranging from materials science to biology. Controlling the thermodynamics and kinetics of growth is advantageous to help tune the formation of complex morphologies. In this report, we harness wide-angle X-ray scattering and vibrational spectroscopy to elucidate the mechanism for crystallization and growth of the metal–organic framework Co-MOF-74 within microscopic volumes enclosed in a capillary and an attenuated total reflection microchip reactor. The experiments reveal molecular and structural details of the growth processes, while the results of plane wave density functional calculations allow identification of lattice and linker modes in the formed crystals. Synthesis of the metal–organic framework with microscopic volumes leads to monodisperse and micron-sized crystals, in contrast to those typically observed under bulk reaction conditions. Reduction in the volume of reagents within the microchip reactor was found to accelerate the reaction rate. The coupling of spectroscopy with scattering to probe reactions in microscopic volumes promises to be a useful tool in the synthetic chemist's kit to understand chemical bonding and has potential in designing complex materials.

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

Electronic structure of NdO via slow photoelectron velocity-map imaging spectroscopy of NdO –

Electronically excited NdO is a possible product of the chemistry associated with the release of Nd into the ionosphere, and emission from these states may contribute to the observations following such experiments. To better characterize the energetics and spectroscopy of NdO, we report a combined experimental and theoretical study using slow photoelectron velocity-map imaging spectroscopy of cryogenically cooled NdO – anions (cryo-SEVI) supplemented by wave function-based quantum-chemical calculations. Using cryo-SEVI, we measure the electron affinity of NdO to be 1.0091(7) eV and resolve numerous transitions to low-lying electronic and vibrational states of NdO that are assigned with the aid of the electronic structure calculations. Additionally, temperature-dependent data suggest contributions from the (2)4.5 state of NdO – residing 2350 cm –1 above the ground anion state. Finally, photodetachment to higher-lying excited states of NdO is also reported, which may help to clarify observations from prior release experiments.

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

A non-hierarchical correlation discrete variable representation

The correlation discrete variable representation (CDVR) facilitates (multi-layer) multi-configurational time-dependent Hartree (MCTDH) calculations with general potentials. It employs a layered grid representation to efficiently evaluate all potential matrix elements appearing in the MCTDH equations of motion. The original CDVR approach and its multi-layer extension show a hierarchical structure: the size of the grids employed at the different layers increases when moving from an upper layer to a lower one. In this work, a non-hierarchical CDVR approach, which uses identically structured quadratures at all layers of the MCTDH wavefunction representation, is introduced. Here, the non-hierarchical CDVR approach crucially reduces the number of grid points required, compared to the hierarchical CDVR, shows superior scaling properties, and yields identical results for all three representations showing the same topology. Numerical tests studying the photodissociation of NOCl and the vibrational states of CH 3 demonstrate the accuracy of the non-hierarchical CDVR approach.

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

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

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

Physical origin of enhanced electrical conduction in aluminum-graphene composites

In this study, the electronic and transport properties of aluminum-graphene composite materials were investigated using the ab initio plane wave density functional theory. The interfacial structure is reported for several configurations. In some cases, the face-centered aluminum (111) surface relaxes in a nearly ideal registry with graphene, resulting in a remarkably continuous interface structure. The Kubo–Greenwood formula and space-projected conductivity were employed to study electronic conduction in aluminum single- and double-layer graphene-aluminum composite models. The electronic density of states at the Fermi level is enhanced by the graphene for certain aluminum–graphene interfaces, thus improving electronic conductivity. In double-layer graphene composites, conductivity varies non-monotonically with temperature, showing an increase between 300 and 400 K at short aluminum-graphene distances, unlike the consistent decrease in single-layer composites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗