Search NASA⌕ Search

SEARCH · Search NASA

Results for “Classical density functional theory”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Many-body quantum dynamics slows down at low density

We study quantum many-body systems with a global U(1) conservation law, focusing on a theory of N interacting fermions with charge conservation, or N interacting spins with one conserved component of total spin. We define an effective operator size at finite chemical potential through suitably regularized out-of-time-ordered correlation functions. The growth rate of this density-dependent operator size vanishes algebraically with charge density; hence we obtain new bounds on Lyapunov exponents and butterfly velocities in charged systems at a given density, which are parametrically stronger than any Lieb-Robinson bound. We argue that the density dependence of our bound on the Lyapunov exponent is saturated in the charged Sachdev-Ye-Kitaev model. We also study random automaton quantum circuits and Brownian Sachdev-Ye-Kitaev models, each of which exhibit a different density dependence for the Lyapunov exponent, and explain the discrepancy. We propose that our results are a cartoon for understanding Planckian-limited energy-conserving dynamics at finite temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Charge transfer rate constants for the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran: The spin-boson model vs the linearized semiclassical approximation

Charge transfer rate constants were calculated for the carotenoid-porphyrin-C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran. The calculation was based on mapping the all-atom anharmonic Hamiltonian of this system onto the spin-boson Hamiltonian. The mapping was based on discretizing the spectral density from the time correlation function of the donor–acceptor potential energy gap, as obtained from all-atom molecular dynamics simulations. Different spin-boson Hamiltonians were constructed for each of the possible transitions between the three excited electronic states in two different triad conformations. The rate constants of three possible transitions were calculated via the quantum-mechanically exact Fermi’s golden rule (FGR), as well as a progression of more approximate expressions that lead to the classical Marcus expression. The advantage of the spin-boson approach is that once the mapping is established, the quantum-mechanically exact FGR and the hierarchy of approximations are known in closed form. The classical Marcus charge transfer rate constants obtained with the spin-boson Hamiltonians were found to reproduce those obtained from all-atom simulations with the linearized semiclassical approximation, thereby confirming the equivalence of the two approaches for this system. Furthermore, within the spin-boson Hamiltonian, we also found that the quantum-mechanically exact FGR rate constants were significantly enhanced compared to the classical Marcus theory rate constants for two out of three transitions in one of the two conformations under consideration. The results confirm that mapping to the spin-boson model can yield accurate predictions for charge transfer rate constants in a system as complex as CPC 60 dissolved in tetrahydrofuran.

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

DOE STTR Phase I Final Report Report: Machine-learning Based Prediction of Thermal Limits for Conjugated Organic Materials

The SCANN-DT Phase I STTR project led by NLM Photonics and partnering with the National Renewable Energy Laboratory (NREL) sought to apply machine learning techniques based on graph neural networks (GNNs) towards the prediction of decomposition temperatures (Td) of organic semiconductors, based on prior work on bond dissociation energy (BDE) prediction as implemented in NREL’s ALFABET prediction tool. Using a curated set of experimental decomposition energies, the project examined GNN-based, classical quantitative structure-property relationship (QSPR) based on DFT calculations, and combinations of both methods to predict Td. While the best MAEs in Td achieved were near 40°C, below project targets, the project led to improvements in the ALFABET model, improvements in cloud-based implementations of NWChem software, and a substantial dataset of calculations on medium-sized conjugated organic molecules.

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

The Role of Cation Coordination in the Electrical and Optical Properties of Amorphous Transparent Conducting Oxides

Amorphous oxide semiconductor materials have demonstrated numerous advantages without compromise of electrical properties as compared to their crystalline counterparts, yet understanding of the fundamental principles allowing this has remained elusive. To study the origins of enhanced optoelectronic properties, we apply high-throughput, combinatorial sputtering, structural and spectral mapping, and computationally intensive ab initio molecular dynamics simulations with density functional theory to a ternary, post-transition metal oxide system, namely, zinc tin oxide. The deposited thin films exhibit a high figure of merit, achieving carrier densities in the range of 1019 to 1020 cm–3 and carrier mobilities up to 35 cm2/Vs. These results highlight the role of local distortions and cation coordination in determining the microscopic origins of carrier generation and transport. In particular, we identify the strong likelihood of Sn undercoordination in both Zn-poor and Zn-rich phases leading to the high carrier concentrations observed. This not only diverges from the still widespread historical indictment of oxygen vacancies controlling carrier population in crystalline oxides but also provides a comprehensive framework to describe the unique structure–property relationships using specific structural and electronic descriptors in disordered phase materials.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2022

Dissolving salt is not equivalent to applying a pressure on water

Abstract Salt water is ubiquitous, playing crucial roles in geological and physiological processes. Despite centuries of investigations, whether or not water’s structure is drastically changed by dissolved ions is still debated. Based on density functional theory, we employ machine learning based molecular dynamics to model sodium chloride, potassium chloride, and sodium bromide solutions at different concentrations. The resulting reciprocal-space structure factors agree quantitatively with neutron diffraction data. Here we provide clear evidence that the ions in salt water do not distort the structure of water in the same way as neat water responds to elevated pressure. Rather, the computed structural changes are restricted to the ionic first solvation shells intruding into the hydrogen bond network, beyond which the oxygen radial-distribution function does not undergo major change relative to neat water. Our findings suggest that the widely cited pressure-like effect on the solvent in Hofmeister series ionic solutions should be carefully revisited.

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

An extended moments model of quantum efficiency for metals and semiconductors

The complexity of photocathode designs and detector materials, and the need to model their performance for short pulse durations, response to high frequency photons, presence of coatings and/or thinness of the absorptive layer, necessitates modifications to three-step and moments models of photoemission that are used in simulation codes. Here, methods to include input from computationally intensive approaches, such as density functional theory to model optical properties and transfer matrix approaches to treat emission from the surface or transport past coatings, by means of parametric models are demonstrated. First, a technique to accurately represent optical behavior so as to model reflectivity and penetration depth is given. Second, modifications to bulk models arising from the usage of thin film architectures, and a means to rapidly calculate them, is provided. Third, a parameterization to model the impact of wells associated with coatings and surface layers on the transmission probably is given. In all cases, the methods are computationally efficient and designed to allow for including input from numerically intensive approaches that would otherwise be unavailable for simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Magnetic and structural properties of Mn X NiSn ( X = Mn, Fe, Co)

Crystal structure and magnetic properties of Heusler alloys MnXNiSn (X = Mn, Fe, Co) are investigated using density functional theory and compared with experimental results. The parent alloy Mn 2 NiSn, which crystallizes in the inverse Heusler structure, is found to be ferrimagnetic, in agreement with previous experimental and theoretical work. The Fe and Co substitutions cause the alloys to assume a fairly well-ordered Y structure and enhance the magnetization substantially. We find that the strong nearest neighbour Mn-X exchange changes from antiferromagnetic (X = Mn) to ferromagnetic (X = Fe, Co), which explains and actually overestimates the experimental changes. A striking feature of the system is that Fe and Co have opposite effects on the Curie temperature T c : they reduce and enhance T c , respectively. We qualitatively explain this behaviour in terms of two-sublattice model based on the nearest-neighbour exchange.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Direct Measurement of Ice-Ablator Interface Motion for Instability Mitigation in Indirect Drive ICF Implosions

In indirect drive inertial confinement fusion (ICF) implosions hydrodynamic instability growth at the imploding capsule ablator-DT fuel interface can reduce fuel compressibility and inject ablator into the hot spot hence reducing hot spot pressure and temperature. As a mitigation strategy, a gentle acceleration of this interface is predicted by simulations and theory to significantly reduce this instability growth in the early stage of the implosion. We have performed high-contrast, time-resolved x-ray refraction enhanced radiography (RER) to accurately measure the level of acceleration as a function of the initial laser drive time history for indirect-drive implosions on the National Ignition Facility. Here, we demonstrate a transition from no acceleration to 20 ± 1.8μm ns -2 acceleration by tweaking the drive that should reduce the initial instabilities by an order of magnitude at high modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Controlled n –Doping of Naphthalene Diimide–Based Two–Dimensional Polymers

Two-dimensional polymers (2DPs) are promising as structurally well-defined, permanently porous, organic semiconductors. However, 2DPs are nearly always isolated as closed shell organic species with limited charge carriers, which leads to low bulk conductivities. Here, we enhance the bulk conductivity of two naphthalene diimide (NDI)-containing 2DP semiconductors by controllably n-doping the NDI units using cobaltocene (CoCp 2 ). Optical and transient microwave spectroscopy reveals that both as-prepared NDI-containing 2DPs are semiconducting with sub-2 eV optical bandgaps and photoexcited charge-carrier lifetimes of tens of nanoseconds. Following reduction with CoCp 2 , both 2DPs largely retain their periodic structures and exhibit optical and electron-spin resonance spectroscopic features consistent with the presence of NDI-radical anions. Here, while the native NDI-based 2DPs are electronically insulating, maximum bulk conductivities of >10 –4 S cm –1 are achieved by substoichiometric levels of n-doping. Density functional theory calculations show that the strongest electronic couplings in these 2DPs exist in the out-of-plane (π-stacking) crystallographic directions, which indicates that cross-plane electronic transport through NDI stacks is primarily responsible for the observed electronic conductivity. Taken together, this study underlines that controlled molecular doping is a useful approach to access structurally well-defined, paramagnetic, 2DP n-type semiconductors with measurable bulk electronic conductivities of interest for electronic or spintronic devices.

2D polymers↗
DOE OSTI2022

Bandgap analysis and carrier localization in cation-disordered ZnGeN 2

The bandgap of ZnGeN 2 changes with the degree of cation site disorder and is sought in light emitting diodes for emission at green to amber wavelengths. By combining the perspectives of carrier localization and defect states, we analyze the impact of different degrees of disorder on electronic properties in ZnGeN 2 , addressing a gap in current studies, which largely focus on dilute or fully disordered systems. The present study demonstrates changes in the density of states and localization of carriers in ZnGeN 2 calculated using bandgap-corrected density functional theory and hybrid calculations on partially disordered supercells generated using the Monte Carlo method. We use localization and density of states to discuss the ill-defined nature of a bandgap in a disordered material and identify site disorder and its impact on the structure as a mechanism controlling electronic properties and potential device performance. Decreasing the order parameter results in a large reduction of the bandgap. The reduction in bandgap is due, in part, to isolated, localized states that form above the valence band continuum associated with nitrogen coordinated by more zinc than germanium. The prevalence of defect states in all but the perfectly ordered structure creates challenges for incorporating disordered ZnGeN 2 into optical devices, but the localization associated with these defects provides insight into the mechanisms of electron/hole recombination in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

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

Chemistry-mediated Ostwald ripening in carbon-rich C/O systems at extreme conditions

Abstract There is significant interest in establishing a capability for tailored synthesis of next-generation carbon-based nanomaterials due to their broad range of applications and high degree of tunability. High pressure (e.g., shockwave-driven) synthesis holds promise as an effective discovery method, but experimental challenges preclude elucidating the processes governing nanocarbon production from carbon-rich precursors that could otherwise guide efforts through the prohibitively expansive design space. Here we report findings from large scale atomistically-resolved simulations of carbon condensation from C/O mixtures subjected to extreme pressures and temperatures, made possible by machine-learned reactive interatomic potentials. We find that liquid nanocarbon formation follows classical growth kinetics driven by Ostwald ripening (i.e., growth of large clusters at the expense of shrinking small ones) and obeys dynamical scaling in a process mediated by carbon chemistry in the surrounding reactive fluid. The results provide direct insight into carbon condensation in a representative system and pave the way for its exploration in higher complexity organic materials. They also suggest that simulations using machine-learned interatomic potentials could eventually be employed as in-silico design tools for new nanomaterials.

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

Reimagining the $e_g$ 1 Electronic State in Oxygen Evolution Catalysis: Oxidation-State-Modulated Superlattices as a New Type of Heterostructure for Maximizing Catalysis

We report that the discovery of solid-phase, inexpensive transition-metal-based water oxidation catalysts is a central goal for renewable energy, and has led to a general consensus that a partially populated metal e g d-electronic state is desirable, leading to favorable catalysis for certain elements in specific oxidation states. In manganese systems, the key species is manganese(III), whose high-spin d 4 electronic configuration places an unpaired electron in the e g orbital, which is postulated to contribute to electronic and structural features that support catalysis. Based on density functional theory calculations, it is predicted that electron transfer would be facilitated by a catalyst with alternating low- and high-Mn III -content sheets, which positions neighboring band edges in closer energetic proximity. The preparation of such catalysts is demonstrated for the first time and it is shown that the catalytic activity is maximized in these systems over more uniform, but more Mn III -rich systems. The best catalyst possesses alternating high-and low-average oxidation state sheets with interlayer Cs + ions, and has an overpotential of 450 mV at 10 mA, which represents an improvement of 250 mV over the best unmodified synthetic potassium birnessites. Using scanning tunneling spectroscopy, bandgap modulations consistent with the theoretically predicted band edge shifts are detected.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular Modes of Attosecond Charge Migration

Here, first-principles calculations are employed to elucidate the modes of attosecond charge migration (CM) in halogenated hydrocarbon chains. We use constrained density functional theory (DFT) to emulate the creation of a localized hole on the halogen and follow the subsequent dynamics via time-dependent DFT. We find low-frequency CM modes (~ 1 eV) that propagate across the molecule and study their dependence on length, bond order, and halogenation. We observe that the CM speed (~ 4Å/ fs) is largely independent of molecule length, but is lower for triple-bonded versus double-bonded molecules. Additionally, as the halogen mass increases, the hole travels in a more particlelike manner as it moves across the molecule. These heuristics will be useful in identifying molecules and optimal CM detection methods for future experiments, especially for halogenated hydrocarbons which are promising targets for ionization-triggered CM.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Predicted thermophysical properties of UN, PuN and (U,Pu)N

Molecular dynamics and density functional theory simulations are used to predict the lattice and electronic contributions of thermophysical properties for UN, PuN, and mixed (U,Pu)N systems. The properties predicted include the lattice parameter, linear thermal expansion, enthalpy, and specific heat capacity, as a function of temperature. The simulation predictions for high temperature specific heat capacity are compared against experimental measurements to understand the behavior, and why differences in the experimental measurements are observed. The influence of adding U vacancies, N interstitials, and Pu to UN is also examined. For this, a new PuN potential parameter set is developed and used with the Kocevski UN potential, enabling the dynamics of mixed (U,Pu)N systems to be studied. How defects impact the thermophysical properties is important for understanding fuel behavior under different reactor conditions, and these mechanistic predictions can be used to support fuel performance codes where data is scarce.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol

In this work, synchrotron X-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol and D-sorbitol) from above the glass transition temperatures $T_g$ to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently reflecting the different natures of bonding. Upon cooling from 400 K, the O$\cdots$O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrest at $T_g$; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above $T_g$, glycerol loses the fewest HBs per OH, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C$\cdots$C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above $T_g$ causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

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

Stoichiometry dependent properties of cerium hydride: An active learning developed interatomic potential study

Cerium hydride has a variety of interesting properties, including a known lattice contraction and densification with increasing hydrogen content. However, precise stoichiometric control is not experimentally straightforward and ab initio approaches are not computationally feasible for many properties such as melting and low temperature diffusion. Therefore, we develop a machine-learned interatomic potential for cerium hydride that is valid for H to Ce ratios from 2.0 to 3.0. A query-by-committee active learning approach is used to develop the training set. Leveraging classical molecular dynamics simulations, we assess a range of properties and provide fundamental mechanisms for the trends with stoichiometry. Finally, a majority of the properties follow the trend of lattice contraction, being governed by the stronger lattice binding induced by adding octahedral atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum real-time evolution of entanglement and hadronization in jet production: Lessons from the massive Schwinger model

The possible link between entanglement and thermalization, and the dynamics of hadronization are addressed by studying the real-time response of the massive Schwinger model coupled to external sources. This setup mimics the production and fragmentation of quark jets, as the Schwinger model and quantum chromodynamics (QCD) share the properties of confinement and chiral symmetry breaking. By using simulations of quantum dynamics on classical hardware, we study the entanglement between the produced jets, and observe the growth of the corresponding entanglement entropy in time. This growth arises from the increased number of contributing eigenstates of the reduced density matrix with sufficiently large and close eigenvalues. We also investigate the physical nature of these eigenstates, and find that at early times they correspond to fermionic Fock states. We then observe the transition from these fermionic Fock states to mesonlike bound states as a function of time. In other words, we observe how hadronization develops in real time. At late times, the local observables at midrapidity (such as the fermion density and the electric field) approach approximately constant values, suggesting the onset of equilibrium and approach to thermalization. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗