SEARCH · Search NASA
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.
This is the poster I will present at the GRC Aqueous Corrosion meeting detailing our latest work on integrating Machine Learning into the Computational Calculations of Galvanic Corrosion
A computer-implemented method of preventing computationally explosive calculations. The method includes obtaining, by a processor of the computer, measured data of one of a physical process or a physical object; performing hierarchical numerical modeling of a physical process inclusive of an Earth model containing at least one of (a) infrastructure in the ground and (b) a formation feature in the ground, wherein predicted data is generated; comparing the measured data to the predicted data to calculate an estimated error; analyzing the estimated error via an inversion process to update the at least one of the Earth model and infrastructure model so as to reduce the estimated error and to determine a final composite Earth model of at least one of the infrastructure and the feature; and using the final composite Earth model to characterize at least one of the process and the physical object.
Calculating the effects of coherent synchrotron radiation (CSR) is one of the most computationally expensive tasks in accelerator physics. Here, we use convolutional neural networks (CNNs), along with a latent conditional diffusion (LCD) model, trained on physics-based simulations to speed up calculations. Specifically, we produce the 3D CSR wakefields generated by electron bunches in circular orbit in the steady-state condition. Two datasets are used for training and testing the models: wakefields generated by three-dimensional Gaussian electron distributions and wakefields from a sum of up to 25 three-dimensional Gaussian distributions. The CNNs are able to accurately produce the 3D wakefields ∼250–1000 times faster than the numerical calculations, while the LCD achieves a gain of a factor of ∼34. We also test the extrapolation and out-of-distribution generalization ability of the models. They generalize well on distributions with larger spreads than what they were trained on but struggle with smaller spreads.
Calculating the equilibrium properties of condensed-matter systems is one of the promising applications of near-term quantum computing. Recently, hybrid quantum-classical time-series algorithms have been proposed to efficiently extract these properties from a measurement of the Loschmidt amplitude ⟨ ψ | e − i H ^ t | ψ ⟩ from initial states | ψ ⟩ and a time evolution under the Hamiltonian H ^ up to short times t . In this work, we study the operation of this algorithm on a present-day quantum computer. Specifically, we measure the Loschmidt amplitude for the Fermi-Hubbard model on a 16 -site ladder geometry (32 orbitals) on the Quantinuum H2-1 trapped-ion device. We assess the effect of noise on the Loschmidt amplitude and implement algorithm-specific error-mitigation techniques. By using a thus-motivated error model, we numerically analyze the influence of noise on the full operation of the quantum-classical algorithm by measuring expectation values of local observables at finite energies. Finally, we estimate the resources needed for scaling up the algorithm. Published by the American Physical Society 2024
The objective is to evaluate the capabilities of multiple CSDs (provided by NDG Systems) using Hadoop Filesystem and Apache Spark. The independent variables are: number of CSDs, 0, 1, 2, 4, or 6; size of dataset, 1 GB, 5 GB, 10 GB; type of dataset, one large file with all of the data, 10 files, 100 files. The dependent variables are: job time; execution time. the constants are operations on the dataset.
Hydrogen storage materials form a crucial research topic for future energy utilization employing hydrogen and among those of interest magnesium diboride (MgB 2 ) has shown its prevalence. In this study, a first-principles analytical adsorption model of one hydrogen molecule in the vicinity of various magnesium diboride crystal surfaces was developed in order to obtain surface thermodynamic properties as a function of molecular and lattice properties. Henry's law constant (K H ) and isosteric heat of adsorption (ΔH ads ) indicators of the affinity between a gaseous molecule and a solid surface are thus calculated. The results in this paper not only address questions pertaining to the first stage of hydrogen storage processes but also advance the understanding of physisorption thermodynamics of a neutral molecule (H 2 ) coming in contact with a layered metallic-like surface (MgB 2 ). Although the model is built from a framework of classical calculations, quantum effects are incorporated as the fractional charge of the ions on the free surfaces, which is essential for the calculation of analytic thermodynamic values that approximate calculations from other methods. To benchmark our theoretical models, periodic density functional calculations were performed to determine the interactions between H 2 and different MgB 2 surfaces from first-principles. By considering both the top and sublayers of MgB 2 in calculating interaction energy, we have analytically and computationally calculated the interaction energies of H 2 molecules and MgB 2 's terminated planes, and witnessed the strong dependence of interaction energies on surface charges. We have also observed a dipole flipping phenomenon which explains the discontinuity seen in the interaction energy graph of Mg(0001). Both analytical and computational results showed heat of adsorption at zero coverage varying at a very low range (<7 kJ mol -1 ).
Herein we propose a quantum-classical hybrid variational algorithm, the quantum orbital minimization method (qOMM), for obtaining the ground state and low-lying excited states of a Hermitian operator. Given parametrized ansatz circuits representing eigenstates, qOMM implements quantum circuits to represent the objective function in the orbital minimization method and adopts a classical optimizer to minimize the objective function with respect to the parameters in ansatz circuits. The objective function has an orthogonality constraint implicitly embedded, which allows qOMM to apply a different ansatz circuit to each input reference state. We carry out numerical simulations that seek to find excited states of H 2 , LiH, and a toy model consisting of four hydrogen atoms arranged in a square lattice in the STO-3G basis with UCCSD ansatz circuits. Comparing the numerical results with existing excited states methods, qOMM is less prone to getting stuck in local minima and can achieve convergence with more shallow ansatz circuits.
The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.
Ligand selectivity to specific lanthanide (Ln) ions is key to the separation of rare earth elements from each other. Ligand selectivity can be quantified with relative stability constants (measured experimentally) or relative binding energies (calculated computationally). The relative stability constants of EDTA (ethylenediaminetetraacetic acid) with La 3+ , Eu 3+ , Gd 3+ , and Lu 3+ were predicted from relative binding energies, which were quantified using electronic structure calculations with relativistic effects and based on the molecular structures of Ln–EDTA complexes in solution from density functional theory molecular dynamics simulations. The protonation state of an EDTA amine group was varied to study pH ~7 and ~11 conditions. Further, simulations at 25 °C and 90 °C were performed to elucidate how structures of Ln–EDTA complexes varying with temperature are related to complex stabilities at different pH conditions. Relative stability trends are predicted from computation for varying Ln 3+ ions (La, Eu, Gd, Lu) with a single ligand (EDTA at pH ~11), as well as for a single Ln 3+ ion (La) with varying ligands (EDTA at pH ~7 and ~11). Changing the protonation state of an EDTA amine site significantly changes the solution structure of the Ln–EDTA complex resulting in a reduction of the complex stability. As a result, increased Ln–ligand complex stability is correlated to reduced structural variations in solution upon an increase in temperature.
This poster highlights the use of computational hydrogen electrode approach to explain the improved activity and selectivity of Ag nanoparticle catalysts for CO2 to CO conversion observed in the experiments. The calculations predict a charge transfer from the Ag nanoparticle to a defective carbon surface, stabilizing the *COOH intermediate through reduced antibonding orbital overlap, significantly reducing the *COOH formation energy barrier, and improving CO2-to-CO conversion selectivity compared with Ag nanocluster on defect-free carbon. These results provide new insights into carbon-supported electrocatalysts for CO2RR and introduce a new approach for creating active and selective nanocatalysts.
Nanostructures in the form of ellipsoids, prolate spheroids, rings, and cylinders are known to exhibit resonant surface and cavity modes with applications in nanophotonics and plasmonics and, more recently, in novel quantum experiments, in which control of plasmons and their interactions with plasmons, photons, phonons, excitons, and quantum emitters are desired. Nanorods and nanowires are examples of plasmonic structures with spectral properties of potential use as interconnects and circuit components. Additionally, estimates of the surface properties of these components are needed in circuit design and integrated systems. Here, we present a quantum Hamiltonian for the cylindrical surface charge density. We then study the photon excitation of plasmons on the cylindrical surface and calculate their scattering and radiative decay rate. Nonradiative decay of plasmons induces an efficient heating of the nanoparticle and can photoacoustically excite mechanical oscillations. Computational calculations are also presented for the plasmonic modes and the ensuing excitations of nanomechanical eigenmodes of nanoparticles with near-cylindrical symmetries.
In this study, the variational quantum eigensolver (VQE) on a quantum simulator is used in calculating ground state electronic structure properties of the LiH n , n = 1–3, complexes including their singly charged ions. Results calculated using classical electronic structure algorithms are also included. We investigate the use of the unitary coupled cluster with singles and doubles (UCCSD) Ansatz using VQE within Qiskit and compare results to full configuration interaction (FCI) calculations. Computed ground state energies, electron affinities, ionization potentials, and dipole moments are considered. We report the first-of-its-kind simulated quantum computing results of selected LiH n species and use the parity orbital to qubit mapping scheme. We find that VQE/UCCSD results are comparable to classical coupled clusters with singles and doubles for all considered systems with respect to FCI. A VQE calculation cost evaluation is included in which we evaluate performance using both Jordan–Wigner and parity orbital to qubit mapping schemes. We also discuss some of the current limitations of utilizing VQE for the study of chemical systems.
Ever since ferroelectricity was discovered in HfO2, the question of its origin remains controversial. Here, we probe this question using a combination of Landau theory of phase transitions and first-principles computations. In such an approach, the energy landscape associated with the phase transition between cubic and different experimentally demonstrated phases of HfO2 (tetragonal, monoclinic, orthorhombic Pbca, orthorhombic Pnma, and orthorhombic Pca21) is explored using density functional theory calculations. Computations revealed that stabilization of all but orthorhombic Pbca phase is driven by a single unstable zone-boundary antipolar mode X2−. When coupled with zone-center modes (Γ1+ and Γ3+), it stabilizes the tetragonal phase. Coupling with four additional modes (Γ5+, X3−, X5−, X5+) results in the monoclinic phase, which is the ground state of the material. If, however, Γ5+ mode is replaced with Γ4− mode, orthorhombic polar phase Pca21 is stabilized. The application of this framework to examine the effect of electric field on the ferroelectric phase of hafnia reveals that the field of 5 MV/cm is capable of stabilizing ferroelectric phase over the monoclinic one at 0 K.
New polymeric blend membranes for CO 2 separation were synthesized based on insights from molecular dynamics simulation. A molecular-level structure-property relationship in polymers of intrinsic microporosity (PIM) based blend membranes, was investigated in detail computationally. Calculated local density profiles and energy of interaction of the blend membranes, composed of PIM-1 and various polyphosphazenes, showed that using the polyphosphazene with a higher concentration of ether side chains can improve the compatibility with PIM-1. Furthermore, based on the findings of computational studies, blend membranes were experimentally fabricated from PIM-1 and polyphosphazenes with various polyether side chain concentrations. Polyether concentration in polyphosphazenes was correlated with the film properties and gas transport performance of the blend membranes. Blend membranes showed very high CO 2 permeability (3100-5300 barrer) and improved CO 2 /N 2 selectivity (24-28), outperforming all other PIM-based blend membranes reported to date. Moreover, the CO 2 permeability performance of the blend membranes was tested 566 hours under real post-combustion flue gas from a coal-fired power plant, including CO 2 , N 2 , H 2 O, O 2 , SO x and NO x .
The utilization of solar energy, an abundant and renewable resource, shows great promise. Various technologies have been employed to capture this vast energy potential. In this study, the effectiveness of natural dye-sensitized solar cells (DSSCs) was explored through experimental and computational methods. Organic dyes derived from specific plant species were examined, with a solvent mixture of ethanol, methanol, and tetrahydrofuran used for extraction. A UV–vis spectrophotometer was utilized to measure the sensitizers’ optical characteristics. The study employed density functional theory (DFT), implemented in Gaussian 09 W software, to perform computational calculations. The B3LYP method and a basis set of 6-31G++(d,p) were chosen for optimizing the geometries and energetics of dye molecules. The DFT results indicated that the dye molecules had a bandgap in the range of 2.16–2.38 eV. The photochemical performance of all four fabricated DSSC-based devices was evaluated to be between 0.33% and 1.04%, and the Euryops pectinatus sensitizer demonstrated the highest efficiency of 1.04% among the studied sensitizers.
Five bismuth(III)-organic phases that consist of supramolecular assemblies of Bi-2,6-pyridinedicarboxylate structural units and substituted 1,10-phenanthroline molecules (R-Phen; R = H, 5-methyl, 5-chloro, 2,9-dimethyl, and 2,9-dicholoro) were synthesized. All five compounds exhibited solid-state photoluminescence. Whereas the phases containing 2,9-dimethylphenanthroline (Me 2 Phen) and 2,9-dichlorophenanthroline (Cl 2 Phen) displayed solely phosphorescence, the structures built from 5-methylphenanthroline and 5-chlorophenanthroline showed exclusively fluorescence. The remaining phase, consisting of phenanthrolinium, exhibited both fluorescence and phosphoresence. It was determined that phosphorescence arises from triplet state emission (T 1 → S 0 ) of substituted R-Phen units while fluorescence originates from Bi(III) coordinated pyrdinedicarboxylate ligands. Bismuth induces spin–orbit coupling for triplet state population and additionally acts as a heavy metal attenuator for X-ray luminescence (radioluminescence). The electronic structure was mapped and excitation pathway investigated via density functional theory calculations. Computational findings indicate favorable conditions for triplet energy transfer from donor Bi(III)-organic units to acceptor R-Phen derivatives. It is proposed that for the phosphorescent compounds, strong π–π interactions promote electron transfer, whereas the compounds that exhibit purely fluorescence lack any such π–π interactions and undergo triplet energy transfer. In conclusion, these results provide a useful platform for probing structure–property relationships of luminescent bismuth-organic compounds, and specifically highlights the role of noncovalent interactions in achieving room temperature phosphorescence and radioluminescence.
Herein, the synthesis of Cu 3 (HAB) x (TATHB) 2-x (HAB: hexaaminobenzene, TATHB: triaminotrihydroxybenzene) is reported. Synthetic improvement of Cu 3 (TATHB) 2 leads to a more crystalline framework with higher electrical conductivity value than previously reported. The improved crystallinity and analogous structure between TATHB and HAB enable the synthesis of Cu 3 (HAB) x (TATHB) 2-x with ligand compositions precisely controlled by precursor ratios. The electrical conductivity is tuned from 4.2 × 10 -8 to 2.9 × 10 -5 S cm -1 by simply increasing the nitrogen content in the crystal lattice. Furthermore, computational calculation supports that the solid solution facilitates the band structure tuning. In conclusion, it is envisioned that the findings not only shed light on the ligand-dependent structure–property relationship but create new prospects in synthesizing multicomponent electrically conductive metal-organic frameworks (MOFs) for tailoring optoelectronic device applications.
Here, we predict two different universal scaling regimes for the quantum transmission of metallic nanodevices following the addition of a small amount of uncorrelated disorder. A nanodevice is connected to two thin semi-infinite uniform leads, and the Non-Equilibrium Green’s Function (NEGF) methodology yields the electron transmission $\mathscr{T}(E)$ as a function of the injected electron energy $E$. Ballistic nanodevices have no disorder and have $\mathscr{T}(E)=1$for all $E$ that allow electron propagation in the leads. Quantum dragon nanodevices can have extremely strong properly correlated disorder, and still have $\mathscr{T}(E)=1$for all $E$. Additional uncorrelated site disorder leads to Fano resonances in $\mathscr{T}(E)$. Averaging over the uncorrelated disorder we predict using perturbation theory two universal scaling regimes for $\mathscr{T}_{ave}(E)$. The functional form of both universal scaling regimes depend on the device length and width, energy, and variance of the uncorrelated disorder. The second scaling regime, valid for small but somewhat larger uncorrelated disorder than the first scaling regime, also has the form dependent on the density of states of the system. These two scaling regimes are demonstrated to be valid via large scale computer calculations.