Search NASA⌕ Search

SEARCH · Search NASA

Results for “nuclear energy density functional”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Endo-Fullerene and Doped Diamond Nanocrystallite Based Models of Qubits for Solid-State Quantum Computers

Models of encapsulated 1/2 nuclear spin H-1 and P-31 atoms in fullerene and diamond nanocrystallite, respectively, are proposed and examined with ab-initio local density functional method for possible applications as single quantum bits (qubits) in solid-state quantum computers. A H-1 atom encapsulated in a fully deuterated fullerene, C(sub 20)D(sub 20), forms the first model system and ab-initio calculation shows that H-1 atom is stable in atomic state at the center of the fullerene with a barrier of about 1 eV to escape. A P-31 atom positioned at the center of a diamond nanocrystallite is the second model system, and 3 1P atom is found to be stable at the substitutional site relative to interstitial sites by 15 eV, Vacancy formation energy is 6 eV in diamond so that substitutional P-31 atom will be stable against diffusion during the formation mechanisms within the nanocrystallite. The coupling between the nuclear spin and weakly bound (valance) donor electron coupling in both systems is found to be suitable for single qubit applications, where as the spatial distributions of (valance) donor electron wave functions are found to be preferentially spread along certain lattice directions facilitating two or more qubit applications. The feasibility of the fabrication pathways for both model solid-state qubit systems within practical quantum computers is discussed with in the context of our proposed solid-state qubits.

Park, Seongjun↗
DOE OSTI2023

Modeling of fission gas diffusion and release for Gd 2 O 3 doped UO 2

Uranium dioxide (UO 2 ) is the primary nuclear fuel in light water reactors, and its excess neutronic reactivity can be controlled by adding burnable absorbers, such as Gd 2 O 3 . This burnable absorber has a large neutron absorption cross-section, lowering the high reactivity of the reactor's initial fuel load. However, there needs to be more understanding of how added Gd 2 O 3 influences the properties of UO 2 under irradiation. To understand the behavior of defects and fission gas in the UO 2 /Gd 2 O 3 system under irradiation, we use cluster dynamics modeling supported by density functional theory calculations. First, we calculate the formation energies of Gd point and cluster defects, and evaluate the temperature-dependent defect concentrations using the defect formation energies and entropies. We show that Gd is soluble in UO 2 , introducing a negative charge in the system. Using this information, we adapted the cluster dynamics code Centipede to model the influence of Gd on U self-diffusion and Xe diffusion in UO 2 with 10 wt% Gd 2 O 3 . Also, we analyzed the Xe diffusion as a function of Gd 2 O 3 concentration, showing that the Xe diffusivity is decreased, which means that the athermal diffusivity due to electronic stopping persists at higher temperatures. In conclusion, the decrease in Xe diffusion means that more Xe stays in the matrix, decreasing the Xe release, and lowering its influence of fission gas release on the thermomechanical properties of UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Machine learning assisted prediction of tungsten heavy alloy plasma facing component performance for fusion energy applications

Tungsten and tungsten heavy alloys (WHAs), known for their remarkably high hardness, durability, and corrosion resistance, play a critical role in the thriving development of nuclear fusion reactors in recent years. However, the exploration in tungsten alloys for the nuclear-related applications has been limited by the difficulty of manufacturing and the complexity of experiments to reproduce the environment of nuclear reaction. Therefore, this project aims to utilize nanoscale simulation methods such as density functional theory (DFT) and molecular dynamics (MD) with the help of machine learning techniques to not only understand the mechanisms of tungsten alloys but also allow us to computationally predict their mechanical behaviors under extreme environments. One critical problem of the application of WHAs in nuclear reactors is the surface melting. In the current design of the SPARC reactor, the WHA, W97Ni2.1Fe0.9 or W97NiFe, is chosen to be the first wall components to confine the plasma where the particles are fiercely moving and colliding into each other to create nuclear fusion reaction. This process will generate extremely high heat flux onto these WHA tiles, leaving high surface temperature that could possibly melt the surface of the WHA tiles, As illustrated in Fig. 1(a). a laser experiment previously done illustrates that a rough surface damage would be made after the surface melting where the matrix area mainly composed of nickel and iron as shown in Fig. 1(b), will first melt and then leave vacancies between these tungsten grains. Unfortunately, these kinds of roughness on the first-wall components could be deadly to the plasma inside a Tokmak reactor because the heat that is supposed to dissipate at a designed ratio through the tiles may in turn be excessively absorbed and accumulated on any uneven area of the surface, which will eventually make the whole nuclear reaction fail. In this project, we will introduce a machine learning potential, Allegro, based on DFT calculation and then build a MD model for W-Ni-Fe alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Stable, Low Permeable TEMPO Catholyte for Aqueous Total Organic Redox Flow Batteries

We report aqueous organic redox flow batteries (AORFBs) are highly attractive for large-scale energy storage because of their nonflammability, low cost, and sustainability. (2,2,6,6-Tetramethylpiperidin-1-yl)oxyl (TEMPO) derivatives, a class of redox active molecules bearing air-stable free nitroxyl radicals and high redox potential (>0.8 V vs NHE), has been identified as promising catholytes for AORFBs. However, reported TEMPO based molecules are either permeable through ion exchange membranes or not chemically stable enough for long-term energy storage. Herein, a new TEMPO derivative functionalized with a dual-ammonium dicationic group, N 1 , N 1 , N 1 , N 3 , N 3 , 2, 2, 6, 6-nonamethyl-N 3 -(piperidinyloxy)propane-1,3-bis(ammonium) dichloride (N 2 -TEMPO) as a stable, low permeable catholyte for AORFBs is reported. Ultraviolet–visible (UV–vis) and proton nuclear magnetic resonance ( 1 H-NMR) spectroscopic studies reveal its exceptional stability and ultra-low permeability (1.49 × 10 -12 cm 2 s -1 ). Coupled with 1,1'-bis[3-(trimethylammonio)propyl]-4,4'-bipyridinium tetrachloride ((NPr) 2 V) as an anolyte, a 1.35 V N2-TEMPO/(NPr)2V AORFB with 0.5 m electrolytes (9.05 Wh L -1 ) delivers a high power density of 114 mW cm -2 and 100% capacity retention for 400 cycles at 60 mA cm-2. At 1.0 m electrolyte concentrations, the N 2 -TEMPO/(NPr) 2 V AORFB achieves an energy density of 18.1 Wh L -1 and capacity retention of 90% for 400 cycles at 60 mA cm -2 .

TEMPO↗
DOE OSTI2022

Elucidating the Synergic Effect in Nanoscale MoS 2 /TiO 2 Heterointerface for Na‐Ion Storage

Abstract Interface engineering in electrode materials is an attractive strategy for enhancing charge storage, enabling fast kinetics, and improving cycling stability for energy storage systems. Nevertheless, the performance improvement is usually ambiguously ascribed to the “synergetic effect”, the fundamental understanding toward the effect of the interface at molecular level in composite materials remains elusive. In this work, a well‐defined nanoscale MoS 2 /TiO 2 interface is rationally designed by immobilizing TiO 2 nanocrystals on MoS 2 nanosheets. The role of heterostructure interface between TiO 2 and MoS 2 by operando synchrotron X‐ray diffraction (sXRD), solid‐state nuclear magnetic resonance, and density functional theory calculations is investigated. It is found that the existence of a hetero‐interfacial electric field can promote charge transfer kinetics. Based on operando sXRD, it is revealed that the heterostructure follows a solid‐solution reaction mechanism with small volume changes during cycling. As such, the electrode demonstrates ultrafast Na + ions storage of 300 mAh g −1 at 10 A g −1 and excellent reversible capacity of 540 mAh g −1 at 0.2 A g −1 . This work provides significant insights into understanding of heterostructure interface at molecular level, which suggests new strategies for creating unconventional nanocomposite electrode materials for energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Time-Dependent Density Functional Theory Description of 238 U⁡(n,f), 240,242 Pu⁢(n,f), and 237 Np(n,f) Reactions

In nuclei with an odd nucleon number the nonvanishing spin number density is the source of a pseudomagnetic field, which favors the splitting of the nucleon Cooper pairs. Such a pseudomagnetic field is generated always in the dynamics of any nucleus, but its effects on Cooper pairs are significantly enhanced in the dynamic evolution of nuclei with an odd number of nucleons. We present for the first time a microscopic study of the induced fission of the odd neutron compound nuclei 239 U, 241,243 Pu, and the odd proton, odd neutron compound nucleus 238 Np, performed within the time-dependent density functional theory extended to superfluid fermion systems, without any simplifying assumptions, with controlled numerical approximations, and for a very large number of initial conditions. Because of the presence of the unpaired odd nucleon(s), the time-reversal symmetry of the fission compound nucleus is spontaneously broken, an aspect routinely neglected in the most advanced microscopic approaches of the past. The emerging fission fragment properties are quite similar to the properties of fission fragments of neighboring even-even nuclei. The time from saddle-to-scission is often significantly longer in odd-odd or odd-mass nuclei than for even-even nuclei, since systems with unpaired nucleons are easier to excite and the potential energy surfaces of these nuclei have more structure, often resembling a very complicated obstacle course, rather than a more direct evolution of the nuclear shape from the top of the outer fission barrier to the scission configuration. The Pauli blocking approximation, often invoked in the literature, expected to inhibit the fission of nuclei with unpaired nucleons, is surprisingly strongly violated during the fission dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Level density within a micro-macroscopic approach

Statistical level density $ρ$($E, A$) is derived for nucleonic system with a given energy $E$, particle number $A$ and other integrals of motion in the micro-macroscopic approximation beyond the standard saddle-point method of the Fermi gas model. This level density reaches the two limits; the well-known Fermi gas grand-canonical ensemble limit for a large entropy $S$ related to large excitation energies, and the finite micro-canonical limit for a small combinatorical entropy $S$ at low excitation energies. In conclusion, the inverse level density parameter $K$ as function of the particle number $A$ in the semiclassical periodic orbit theory, taking into account the extended Thomas-Fermi and Strutinsky shell corrections, is calculated and compared with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Computing chemical potentials with machine-learning-accelerated simulations to accurately predict thermodynamic properties of molten salts

The successful design and deployment of next-generation nuclear technologies heavily rely on thermodynamic data for relevant molten salt systems. However, the lack of accurate force fields and efficient methods has limited the quality of thermodynamic predictions from atomistic simulations. Here we propose an efficient free energy framework for computing chemical potentials, which is the central free energy quantity behind many thermodynamic properties. We accelerate our simulations without sacrificing accuracy by using machine learning interatomic potentials trained on density functional theory (DFT) data. Using lithium chloride as our model system, we compute chemical potentials with DFT-accuracy for solid and liquid phases by transmuting ions into noninteracting particles. Notably, in the liquid phase, we demonstrate consistency whether we transmute one ion pair or the entire system into ideal gas particles. By locating the temperature where the chemical potential of solid and liquid phases cross, we predict a melting point of 880 ± 18 K for lithium chloride, which is remarkably close to the experimental value of 883 K. With this successful demonstration, we lay the foundation for high-throughput thermodynamic predictions of many properties that can be derived from the chemical potentials of the minority and majority components in molten salts.

Gibson, Luke D. [Oak Ridge National Laboratory (OR↗
DOE OSTI2020

In-Situ Helium Implantation and TEM Investigation of Radiation Tolerance to Helium Bubble Damage in Equiaxed Nanocrystalline Tungsten and Ultrafine Tungsten-TiC Alloy

The use of ultrafine and nanocrystalline materials is a proposed pathway to mitigate irradiation damage in nuclear fusion components. Here, we examine the radiation tolerance of helium bubble formation in 85 nm (average grain size) nanocrystalline-equiaxed-grained tungsten and an ultrafine tungsten-TiC alloy under extreme low energy helium implantation at 1223 K via in-situ transmission electron microscope (TEM). Helium bubble damage evolution in terms of number density, size, and total volume contribution to grain matrices has been determined as a function of He+ implantation fluence. The outputs were compared to previously published results on severe plastically deformed (SPD) tungsten implanted under the same conditions. Large helium bubbles were formed on the grain boundaries and helium bubble damage evolution profiles are shown to differ among the different materials with less overall damage in the nanocrystalline tungsten. Compared to previous works, the results in this work indicate that the nanocrystalline tungsten should possess a fuzz formation threshold more than one order of magnitude higher than coarse-grained tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fission and the r -process nucleosynthesis of translead nuclei in neutron star mergers

In this report we study the impact of fission on the production and destruction of translead nuclei during the r-process nucleosynthesis occurring in neutron-star mergers. Abundance patterns and rates of nuclear energy production are obtained for different ejecta conditions using three sets of stellar reaction rates, one of which is based on microscopic and consistent calculations of nuclear masses, fission barriers, and collective inertias. We show that the accumulation of fissioning material during the r process can strongly affect the free neutron abundance after the r-process freeze-out. This leads to a significant impact on the abundances of heavy nuclei that undergo α decay or spontaneous fission, affecting the radioactive energy production by the ejecta at timescales relevant for kilonova emission.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Radiation induced athermal diffusivity in uranium mononitride

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuel considered for light water reactors and advanced reactor designs. Properties like self- and fission gas diffusivity need to be better understood, given that they influence key fuel performance phenomena such as fission gas swelling and release. In particular, the radiation induced athermal (D 3 ) diffusivity remains challenging to accurately predict and has only been sparsely characterized in UN, despite its importance as it likely governs diffusion at the low temperatures this high-thermal-conductivity fuel form may operate. Molecular Dynamics simulations are used to estimate the mean square displacement induced by a primary knock-on atom (PKA) with a given kinetic energy. These results are combined with the PKA energy distributions obtained from binary collision approximation calculations to obtain the displacement due to a particular fission fragment. Finally, this is combined with experimental fission fragment yields to determine the displacement due to an average fission event and, thus, express the athermal diffusivity as a function of the fission rate density. These results are in excellent agreement with available experimental data. In conclusion, a particular importance is given to the understanding and the quantification of the variability of these results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Measurements of dihadron correlations relative to the event plane in Au+Au collisions at $\sqrt{s_{NN}} =200$ GeV*

Dihadron azimuthal correlations containing a high transverse momentum ($ p_{T} $) trigger particle are sensitive to the properties of the nuclear medium created at RHIC through the strong interactions occurring between the traversing parton and the medium, i.e. jet-quenching. Previous measurements revealed a strong modification to dihadron azimuthal correlations in Au+Au collisions with respect to p+p and d+Au collisions. The modification increases with the collision centrality, suggesting a path-length or energy density dependence to the jet-quenching effect. This paper reports STAR measurements of dihadron azimuthal correlations in mid-central (20%-60%) Au+Au collisions at $ \sqrt{s_{\rm{NN}}} = 200 $ GeV as a function of the trigger particle's azimuthal angle relative to the event plane, $ \phi_{s} = | \phi_{t}- \psi_{{\rm{EP}}}| $. The azimuthal correlation is studied as a function of both the trigger and associated particle $ p_{T} $. The subtractions of the combinatorial background and anisotropic flow, assuming Zero Yield At Minimum (ZYAM), are described. The correlation results are first discussed with subtraction of the even harmonic (elliptic and quadrangular) flow backgrounds. The away-side correlation is strongly modified, and the modification varies with $ \phi_{s} $, with a double-peak structure for out-of-plane trigger particles. The near-side ridge (long range pseudo-rapidity $ \Delta\eta $ correlation) appears to drop with increasing $ \phi_{s} $ while the jet-like component remains approximately constant. The correlation functions are further studied with the subtraction of odd harmonic triangular flow background arising from fluctuations. It is found that the triangular flow, while responsible for the majority of the amplitudes, is not sufficient to explain the $ \phi_{s} $-dependence of the ridge or the away-side double-peak structure. The dropping ridge with $ \phi_{s} $ could be attributed to a $ \phi_{s} $-dependent elliptic anisotropy; however, the physics mechanism of the ridge remains an open question. Even with a $ \phi_{s} $-dependent elliptic flow, the away-side correlation structure is robust. These results, with extensive systematic studies of the dihadron correlations as a function of $ \phi_{s} $, trigger and associated particle $ p_{T} $, and the pseudo-rapidity range $ \Delta\eta $, should provide stringent inputs to help understand the underlying physics mechanisms of jet-medium interactions in high energy nuclear collisions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu 12-x MxSb 4 S 13 , where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m -1 K -1 ) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive. In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu 10.5 NiZn 0.5 Sb 4 S 13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu 14 Sb 4 S 13 is revealed by Cu self-diffusivity, nuclear density map and “nudged elastic band” (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu1 0 Zn 2 Sb 4 S 13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu 12 Sb 4 S 13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li 3 PO 4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability.

Li, Junchao↗
DOE OSTI2019

Nonadiabatic excited-state intramolecular proton transfer in 3-hydroxyflavone: S2 state involvement via multi-mode effect

The excited-state intramolecular proton transfer process in 3-hydroxyflavone is investigated based on the computed structural parameters and energetics of stationary points of vibronically coupled S 1 -S 2 potential energy surfaces. A conical intersection close to the Franck-Condon point on S1 is identified. The minimum energy of the conical intersection is found to be near-degenerate with the equilibrium minimum of S 2 . Quantum nuclear wavepacket simulations revealed a small population transfer from the “bright” S 1 to “dark” S 2 on a time scale shorter than the O—H stretching vibrational period. Such a nonadiabatic transition opens up the possibility of new photophysical and photochemical pathways, including the proton transfer via S 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗