Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “THOMAS-FERMI MODEL”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Correlation energy of the uniform electron gas determined by ground-state conditional probability density functional theory

Conditional-probability density functional theory (CP-DFT) is a formally exact method for finding correlation energies from Kohn-Sham DFT without evaluating an explicit energy functional. We present details on how to generate accurate exchange-correlation energies for the ground-state uniform gas. We also use the exchange hole in a CP antiparallel spin calculation to extract the high-density limit. We give a highly accurate analytic solution to the Thomas-Fermi model for this problem, showing its performance relative to Kohn-Sham and may be useful at high temperatures. We explore several approximations to the CP potential. Furthermore, results are compared to accurate parameterizations for both exchange-correlation energies and holes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Spontaneous Localization at a Potential Saddle Point from Edge State Reconstruction in a Quantum Hall Point Contact

Quantum point contacts (QPCs) are an essential component in mesoscopic devices. Here, in this work, we study the transmission of quantum Hall edge modes through a gate-defined QPC in monolayer graphene. We observe resonant tunneling peaks and a nonlinear conductance pattern characteristic of Coulomb-blockaded localized states. The in-plane electric polarizability reveals the states are localized at a classically unstable electrostatic saddle point. We explain this unexpected finding within a self-consistent Thomas-Fermi model, finding that localization of a zero-dimensional state at the saddle point is favored whenever the applied confinement potential is sufficiently soft compared to the Coulomb energy. Our results provide a direct demonstration of Coulomb-driven reconstruction at the boundary of a quantum Hall system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Wide-voltage-window reversible control of electronic transport in electrolyte-gated epitaxial BaSnO 3

The wide gap perovskite semiconductor BaSnO 3 has attracted much interest since the discovery of room temperature electron mobility up to 320 cm 2 V -1 s -1 in bulk crystals. Motivated by applications in oxide heterostructures, rapid progress has been made with BaSnO 3 films, although questions remain regarding transport mechanisms and mobility optimization. Here we report on a detailed study of epitaxial BaSnO 3 electric double layer transistors based on ion gel electrolytes, enabling wide-doping-range studies of transport in single films. The work spans an order of magnitude in initial n doping (similar to 10 19 to 10 20 cm -3 , with both oxygen vacancies and La), film thicknesses from 10-50 nm, and measurements of resistance, Hall effect, mobility, and magnetoresistance. In contrast with many oxides, electrolyte gating of BaSnO 3 is found to be essentially reversible over an exceptional gate voltage window (approaching +/- 4 V), even at 300 K, supported by negligible structural modification in operando synchrotron x-ray diffraction. We propose that this occurs due to a special situation in BaSnO3, where electrochemical gating via oxygen vacancies is severely limited by their low diffusivity. Wide-range reversible modulation of transport is thus achieved (in both electron accumulation and depletion modes), spanning strongly localized, weakly localized, and metallic regimes. Two-channel conduction analysis is then combined with self-consistent Schrodinger-Poisson and Thomas-Fermi modeling to extract accumulation layer electron densities and mobilities. Electrostatic electron densities approaching 10 14 cm -2 are shown to increase room temperature mobility by up to a factor of similar to 24. These results lay the groundwork for future studies of electron-density-dependent phenomena in high mobility BaSnO 3 , and significantly elucidate oxide electrolyte gating mechanisms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Electrostatic modulation of the lateral carrier density profile in field effect devices with nonlinear dielectrics

The properties of two-dimensional (2D) electronic systems are often effectively controlled using electrostatic gating. The geometry of such field effect devices influences the effectiveness of the gate and the carrier density profile in the 2D device. Here, in this work, we analyze the gate-induced spatial variations in the lateral carrier density in patterned LaAlO 3 /SrTiO 3 devices. We model the electrostatics of the 2D interface using the Thomas-Fermi approximation and compute the gate-induced charge distribution at the interface. We show that the electric field lines generated by the gate are focused at the edges of the device, causing an increased depletion near its edges. This effect is accentuated in LaAlO 3 /SrTiO 3 due to the large, nonlinear dielectric constant of the substrate, and the large distance between the gate electrode and the interface. We experimentally demonstrate one consequence of this effect by directly imaging current distributions in gated heterostructures, finding that insulating regions nucleate at the edges of the device due to the gate. Our results suggest that device geometry and choice of dielectric materials control the charge distribution in 2D systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Sesame-style decomposition of KS-DFT molecular dynamics for direct interrogation of nuclear models

A common paradigm used in the construction of equations of state is to decompose the thermodynamics into a superposition of three terms: a static-lattice cold curve, a contribution from the thermal motion of the nuclei, and a contribution from the thermal excitation of the electrons. While statistical mechanical models for crystals provide tractable framework for the nuclear contribution in the solid phase, much less is understood about the nuclear contribution above the % melt temperature $(C^{(nuc)}_ν≈3R)$ and how it should transition to the high-temperature limit $(C^{(nuc)}_ν\sim\frac{3}{2}R)$. Here, we describe an algorithm for extracting both the thermal nuclear and thermal electronic contributions from quantum molecular dynamics (QMD) simulations. We then use the VASP QMD package to probe thermal nuclear behavior of liquid aluminum at normal density to compare the results to semi-empirical models – the Johnson generic model, the Chisolm high-temperature liquid model, and the CRIS model. The thermal electronic energy is also calculated from the VASP results and compared to the Kohn-Sham calculations from QE-DFT simulations, Green’s function average atom calculations, and the Thomas-Fermi atomic model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

On the Question of the Existence of One-Dimensional Multielectron Atoms

The semiclassical Thomas–Fermi method of calculating the electrostatic field of a multielectron atom is modified for effectively one-dimensional atoms. The main Thomas–Fermi equation and the dimensionless, one-dimensional Thomas–Fermi equation are obtained, and it is shown that within the framework of the method “one-dimensional” neutral multielectron atoms cannot exist.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

On the Question of the Existence of Two-Dimensional Multi-Electron Atoms

The semiclassical Thomas–Fermi method is extended to the case of two-dimensional (planar) multi-electron atoms. It is shown that within the framework of this model a two-dimensional multi-electron atom cannot exist, due to the impossibility of fulfilling the physical asymptotic limit φ(r)|{sub r→0} → (Ze)/r of the solution of the corresponding two-dimensional Thomas–Fermi equation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Electric fields at finite temperature

Partial differential equations for the electric potential at finite temperature, taking into account the thermal Euler–Heisenberg contribution to the electromagnetic Lagrangian are derived. This complete temperature dependence introduces quantum corrections to several well known equations such as the Thomas–Fermi and the Poisson–Boltzmann equation. Our unified approach allows at the same time to derive other similar equations which take into account the effect of the surrounding heat bath on electric fields. We vary our approach by considering a neutral plasma as well as the screening caused by electrons only. The effects of changing the statistics from Fermi–Dirac to the Tsallis statistics and including the presence of a magnetic field are also investigated. Some useful applications of the above formalism are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to leading order (–0.221 Z 5/3 and –0.021 Z ln Z respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to E x (N,Z) ≈ –0.354 N 2/3 Z (as known before only for Z >> N >> 1) and E c ≈ –0.02 N ln N . These asymptotic estimates are most correct for atomic ions with large N and Z >> N, but we find that they are qualitatively and semi- quantitatively correct even for small N and for N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. Here, it is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Can “Two-” and “One-Dimensional” Multielectron Atoms Exist?

The quasi-classical Thomas–Fermi method is applied to 2D and 1D multielectron atoms. In terms of this method, such atoms are shown not to exist because of the fact that the physical boundary conditions that are analogous to the 3D version of the theory, where boundary conditions are met, cannot be fulfilled. Our theoretical results can be experimentally tested. Atomic number Z{sub 1,} {sub 2max} (~10{sup 2}?) is assumed to exist in terms of this method. At Z > Z{sub 1,2max}, low-dimensional multielectron atoms cannot exist, in contrast to oneor two-electron atoms and, e.g., an experimentally detected Bose condensate of low-dimensional atoms with Z ~ 10 (Na).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Fermi-Teller Theory of Low Velocity Ionization Losses Applied to Monopoles

Recently there has developed a renewed interest in the ionization losses of slow monopoles. An excellent summary of the literature is contained in a paper of Ahlen and Kinoshital (hereafter referred to as A.K.). A.K. have derived lower bounds for the ionization losses of a slow monopole passing through matter. However, they neglected the effects of an energy gap, in insulating materials, on the excitation of atomic electrons. Using the Thomas-Fermi statistical model of atomic wave functions, we have calculated lower bounds for ionization losses in "insulators" including effects due to an energy gap. The Fermi-Teller theory 2 (hereafter referred-to as F.T.) was originally used to predict the stopping, by ionization losses, of slow charged particles in materials. The theory is based on a calculation of the energy lost to a uniform "Fermi sea". However, the authors note that the actual "Fermi velocity" cancels out in the derivation. Therefore their calculated results also apply to a "Thomas-Fermi" atom, in which each volume element is considered to be a Fermi sea filled to the top of the potential well with atomic electrons.

Ritson, D.↗
DOE OSTI2017

A Model for the Characterization of the Polarizability of Thin Films Independently of the Thickness of the Film

It is known that the dielectric properties of thin films can be modified relative to the bulk material because the interaction between film and substrate influences the mobility of the atoms or molecules in the first layers. Here we show that a strong scale effect occurs in nanometer size octadecylammine thin films. Here, this effect is attributed to the different distribution of molecules depending on the size of the film. To accurately describe this effect, we have developed a model which is a reinterpretation of the linearized Thomas-Fermi approximation. Within this model, we have been able to characterize the polarizability of thin films independently of the thickness of the film.

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

Extending and Smoothing Two-dimension Equation of State Simulation Data

The Multiphase Equation of State (MEOS) project works to produce high quality equation of state tables which are used in computationally intensive simulations of materials in different conditions. An equation of state table describes the state of matter under certain physical conditions such as pressure, volume, temperature, or internal energy. MEOS uses many different models and combinations thereof to produce accurate tables consisting of continuous smooth data, derivatives, and higher-order derivatives. Accurate and smooth data are important factors in producing precise simulations. This report focuses on improving the electron tables produced from Purgatorio and Thomas-Fermi data. Purgatorio data is extremely accurate but often jagged and discontinuous in the lower temperature and density region. Thomas-Fermi is a model that provides smooth data throughout. Together, this produces an EOS data table that is both accurate and smooth. This report describes the implementation of a new feature that allows users to define more detailed regions in the Purgatorio table to be replaced with Thomas-Fermi. As a result, MEOS can generate a table that is both more smooth and still accurate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Path integral Monte Carlo simulations of warm dense aluminum

Here, we perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of $0.1-32.4$~g$\,$cm$^{-3}$ and $10^4-10^8$~K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multi-shock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multi-shock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and Z to reveal the evolution of the plasma structure and ionization behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Equation of state of boron nitride combining computation, modeling, and experiment

The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (actex), and all-electron Green's function Korringa-Kohn-Rostoker (mecca), and compute the pressure and internal energy of BN over a broad range of densities and temperatures. Our experiments were conducted at the Omega laser facility and the Hugoniot response of BN to unprecedented pressures (1200–2650 GPa). The EOSs computed using different methods cross validate one another in the warm-dense matter regime, and the experimental Hugoniot data are in good agreement with our theoretical predictions. By comparing the EOS results from different methods, we assess that the largest discrepancies between theoretical predictions are ≲ 4% in pressure and ≲ 3% in energy and occur at 10 6 K, slightly below the peak compression that corresponds to the K -shell ionization regime. At these conditions, we find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the high-pressure, high-temperature regime. Moreover, the recently developed SQ and FOE methods produce EOS data that have significantly smaller statistical error bars than PIMC, and so represent significant advances for efficient computation at high temperatures. The shock Hugoniot predicted by PIMC, actex, and mecca shows a maximum compression ratio of 4.55 ± 0.05 for an initial density of 2.26 g/cm 3 , higher than the Thomas-Fermi predictions by about 5%. In addition, we construct tabular EOS models that are consistent with the first-principles simulations and the experimental data. Our findings clarify the ionic and electronic structure of BN over a broad range of temperatures and densities and quantify their roles in the EOS and properties of this material. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum Scales of Galaxies from Self-interacting Ultralight Dark Matter

We derive the characteristic scales for physical quantities of dwarf galaxies, such as mass, size, acceleration, and angular momentum, within the self-interacting ultralight dark matter (ULDM) model. Due to the small mass of ULDM, even minor self-interactions can drastically alter these scales in the Thomas-Fermi limit. We suggest that these characteristic scales are connected to mysteries of observed galaxies. Oscillation of ULDM field can explain the current cosmological density of dark matter. Many cosmological constraints suggest that the energy scale $\tilde{m}$ for self-interacting ULDM is typically of the order $10~eV$, whereas the mass $m$ for the non-interacting case is around $10^{-21}~eV$. Self-interacting ULDM provides the better explanation for cosmological observations than the non-interacting case.

Astrophysics of Galaxies (astro-ph.GA)↗
DOE OSTI2022

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Atom-in-jellium equations of state in the high-energy-density regime

Recent path-integral Monte Carlo and quantum molecular dynamics simulations have shown that computationally efficient average-atom models can predict thermodynamic states in warm dense matter to within a few percent. One such atom-in-jellium model has typically been used to predict the electron-thermal behavior only, although it was previously developed to predict the entire equation of state (EOS). In this work, we report completely atom-in-jellium EOS calculations for Be, Al, Si, Fe, and Mo, as elements representative of a range of atomic number and low-pressure electronic structure. Comparing the more recent method of pseudoatom molecular dynamics, atom-in-jellium results were similar: sometimes less accurate, sometimes more. All these techniques exhibited pronounced effects of electronic shell structure in the shock Hugoniot which are not captured by Thomas-Fermi based EOS. These results demonstrate the value of a hierarchical approach to EOS construction, using average-atom techniques with shell structure to populate a wide-range EOS surface efficiently, complemented by more rigorous three-dimensional multiatom calculations to validate and adjust the EOS.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗