Search NASA⌕ Search

SEARCH · Search NASA

Results for “Debye 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
NASA NTRS1991

Wave mode identification of electrostatic noise observed with ISEE 3 in the deep tail boundary layer

The characteristics of the VLF electrostatic noise observed with ISEE 3 in the low-latitude boundary layer of distant geomagnetic tail are examined using a display format for the wave dynamic spectra different from that used by Scarf et al. (1984). It is shown that the observed noise is composed of impulsive bursts. The results of the detailed analysis of the noise parameters are used to develop a model of plasma wave behavior in the plasma rest frame. A hypothesis is proposed that the wide frequency extent of the noise spectra is composed of Doppler effects of waves propagating nearly omnidirectionally within the plasma rest frame, which is moving with the electron bulk speed. On the basis of this hypothesis, the wavelength of the observed waves were determined from the width of the frequency extent and the measured electron bulk speed. It is shown that the wavelength ranges from 2 to 8 times the plasma Debye length.

Tsutsui, M.↗
DOE OSTI2016

Effective electrostatic interactions among charged thermo-responsive microgels immersed in a simple electrolyte

This work explores the nature and thermodynamic behavior of the effective electrostatic interactions among charged microgels immersed in a simple electrolyte, taking special interest in the effects due to the thermally induced variation of the microgel size while the remaining parameters (microgel charge and concentration, plus the amount of added salt) are kept constant. To this end, the rigorous approach obtained from applying the precise methodology of the dressed ion theory to the proper definition of the effective direct correlation functions, which emerge from tracing-out the degrees of freedom of the microscopic ions, is employed to provide an exact description of the parameters characterizing such interactions: screening length, effective permittivity, and renormalized charges. A model solution with three components is assumed: large permeable anionic spheres for the microgels, plus small charged hard spheres of equal size for the monovalent cations and anions. The two-body correlations among the components of this model suspension, used as the input for the determination of the effective interaction parameters, are here calculated by using the hyper-netted chain approximation. It is then found that at finite microgel concentrations the values of these parameters change as the microgel size increases, even though the ionic strength of the supporting electrolyte and the bare charge of the microgels remain fixed during this process. The variation of the screening length, as well as that of the effective permittivity, is rather small, but still interesting in view of the fact that the corresponding Debye length stays constant. The renormalized charges, in contrast, increase markedly as the microgels swell. The ratio of the renormalized charge to the corresponding analytic result obtained in the context of an extended linear response theory allows us to introduce an effective charge that accounts for the non-linear effects induced by the short-ranged association of microions to the microgels. The behavior of these effective charges as a function of the amount of added salt and the macroion charge, size, and concentration reveals the interplay among all these system parameters.

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

Inherent issues regarding the use of in situ x-ray diffraction measurements to determine temperature in shock-compressed metals

Temperature determination in shock compressed solids constitutes an important and long-standing scientific need. Since the reduction of Bragg diffracted peaks due to temperature increase, using the Debye-Waller factor, is well established, we examined the use of this approach to determine temperatures in shock compressed gold and platinum by representing the shocked state as a superposition of density and temperature changes. Comparison of the calculated and measured diffraction peaks did not show good agreement, because x-ray diffraction (XRD) profiles in the shock compressed state are not governed solely by density and temperature changes. XRD results are also influenced measurably by shock wave induced microstructural changes. Our results demonstrate that contributions from microstructural changes need to be incorporated and modeled in the theoretical analysis to use XRD measurements for reliable temperature determination in shock compressed solids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Electron-lattice energy relaxation in laser-excited thin-film Au-insulator heterostructures studied by ultrafast MeV electron diffraction

We apply time-resolved MeV electron diffraction to study the electron-lattice energy relaxation in thin film Au-insulator heterostructures. Through precise measurements of the transient Debye-Waller-factor, the mean-square atomic displacement is directly determined, which allows to quantitatively follow the temporal evolution of the lattice temperature after short pulse laser excitation. Data obtained over an extended range of laser fluences reveal an increased relaxation rate when the film thickness is reduced or the Au-film is capped with an additional insulator top-layer. This behavior is attributed to a cross-interfacial coupling of excited electrons in the Au film to phonons in the adjacent insulator layer(s). Analysis of the data using the two-temperature-model taking explicitly into account the additional energy loss at the interface(s) allows to deduce the relative strength of the two relaxation channels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

A molecular Debye-Hückel theory and its applications to electrolyte solutions: The size asymmetric case

We developed a molecular Debye-Hückel theory for electrolyte solutions with size asymmetry, where the dielectric response of an electrolyte solution is described by a linear combination of Debye-Hückel-like response modes. Furthermore, as the size asymmetry of an electrolyte solution leads to a charge imbalanced border zone around a solute, the dielectric response to the solute is characterized by two types of charge sources, namely, a bare solute charge and a charge distribution due to size asymmetry. These two kinds of charge sources are screened by the solvent differently, our theory presents a method to calculate the mean electric potential as well as the electrostatic contributions to thermodynamic properties. Finally, the theory was successfully applied to binary as well as multi-component primitive models of electrolyte solutions.

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

Ultrafast non-radiative dynamics of atomically thin MoSe 2

Non-radiative energy dissipation in photoexcited materials and resulting atomic dynamics provide a promising pathway to induce structural phase transitions in two-dimensional materials. However, these dynamics have not been explored in detail thus far because of incomplete understanding of interaction between the electronic and atomic degrees of freedom, and a lack of direct experimental methods to quantify real-time atomic motion and lattice temperature. Here, we explore the ultrafast conversion of photoenergy to lattice vibrations in a model bi-layered semiconductor, molybdenum diselenide, MoSe 2 . Specifically, we characterize sub-picosecond lattice dynamics initiated by the optical excitation of electronic charge carriers in the high electron-hole plasma density regime. Our results focuses on the first ten picosecond dynamics subsequent to photoexcitation before the onset of heat transfer to the substrate, which occurs on a ~100 picosecond time scale. Photoinduced atomic motion is probed by measuring the time dependent Bragg diffraction of a delayed mega-electronvolt femtosecond electron beam. Transient lattice temperatures are characterized through measurement of Bragg peak intensities and calculation of the Debye-Waller factor (DWF). These measurements show a sub-picosecond decay of Bragg diffraction and a correspondingly rapid rise in lattice temperatures. We estimate a high quantum yield for the conversion of excited charge carrier energy to lattice motion under our experimental conditions, indicative of a strong electron-phonon interaction. First principles nonadiabatic quantum molecular dynamics simulations (NAQMD) on electronically excited MoSe 2 bilayers reproduce the observed picosecond-scale increase in lattice temperature and ultrafast conversion of photoenergy to lattice vibrations. Calculation of excited-state phonon dispersion curves suggests that softened vibrational modes in the excited state are involved in efficient and rapid energy transfer between the electronic system and the lattice.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Validation of NASCAP-2K Spacecraft-Environment Interactions Calculations

The recently released Nascap-2k, version 2.0, three-dimensional computer code models interactions between spacecraft surfaces and low-earth-orbit, geosynchronous, auroral, and interplanetary plasma environments. It replaces the earlier three-dimensional spacecraft interactions codes NASCAP/GEO, NASCAP/LEO, POLAR, and DynaPAC. Nascap-2k has improved numeric techniques, a modern user interface, and a simple, interactive satellite surface definition module (Object ToolKit). We establish the accuracy of Nascap-2k both by comparing computed currents and potentials with analytic results and by comparing Nascap-2k results with published calculations using the earlier codes. Nascap-2k predicts Langmuir-Blodgett or Parker-Murphy current collection for a nearly spherical (100 surfaces) satellite in a short Debye length plasma depending on the absence or presence of a magnetic field. A low fidelity (in geometry and time) Nascap-2k geosynchronous charging calculation gives the same results as the corresponding low fidelity NASCAP/GEO calculation. A high fidelity calculation (using the Nascap-2k improved geometry and time stepping capabilities) gives higher potentials, which are more consistent with typical observations. Nascap-2k predicts the same current as a function of applied potential as was observed and calculated by NASCAP/LEO for the SPEAR I rocket with a bipolar sheath. A Nascap-2k DMSP charging calculation gives results similar to those obtained using POLAR and consistent with observation.

Davis, V. A.↗
DOE OSTI2016

On the connection between microbursts and nonlinear electronic structures in planetary radiation belts

Using a dynamical-system approach, we have investigated the efficiency of large-amplitude whistler waves for causing microburst precipitation in planetary radiation belts by modeling the microburst energy and particle fluxes produced as a result of nonlinear wave–particle interactions. We show that wave parameters, consistent with large-amplitude oblique whistlers, can commonly generate microbursts of electrons with hundreds of keV-energies as a result of Landau trapping. Relativistic microbursts (>1 MeV) can also be generated by a similar mechanism, but require waves with large propagation angles θ{sub kB}>50{sup ∘} and phase-speeds v{sub Φ}⩾c/9. Using our result for precipitating density and energy fluxes, we argue that holes in the distribution function of electrons near the magnetic mirror point can result in the generation of double layers and electron solitary holes consistent in scales (of the order of Debye lengths) to nonlinear structures observed in the radiation belts by the Van Allen Probes. Our results indicate a relationship between nonlinear electrostatic and electromagnetic structures in the dynamics of planetary radiation belts and their role in the cyclical production of energetic electrons (E⩾100 keV) on kinetic timescales, which is much faster than previously inferred.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Determining hexavalent chromium transport properties in alkaline nuclear waste using nuclear magnetic resonance spectroscopy

This study focuses on the transport properties of hexavalent chromium, specifically the chromate anion, to improve predictive models and environmental remediation strategies for Cr(VI) migration. Using 53 Cr Nuclear Magnetic Resonance (NMR) spectroscopy, the research quantifies chromate in multicomponent electrolytes replicating nuclear waste conditions at the Hanford Site in Washington State. The consistency of the 53 Cr NMR signal integral with chromate concentration, despite varying matrix compositions, establishes it as a reliable concentration indicator. The transport properties of chromate in an alkaline solution were assessed using relaxation-based measurements via saturation recovery and Carr-Purcell-Meiboom-Gill experiments, determining spin-lattice and spin-spin relaxation times. These measurements, combined with the Bloembergen-Purcell-Pound equation, helped estimate the rotational correlation time and the 53 Cr self-diffusion coefficient using Stokes-Einstein-Debye and Stokes-Einstein equations. Direct measurements were obtained through pulsed field gradient stimulated echo 53 Cr NMR spectroscopy. Monte Carlo simulations further estimated uncertainty propagation. The results enhance comprehension of chromate transport and highlight prospects for identifying transport properties of NMR-active nuclei, traditionally considered unreachable.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

First-Principles Study of Phase Stability of bcc XZn (X = Cu, Ag, and Au) Alloys

First-principles density-functional theory is used to study the phase stability/instability and anomalies in the formation of the high-temperature bcc phases of X Zn (X = Cu, Ag, and Au) alloys. Although from perhaps a naive point of view, their properties are expected to monotonically depend on the noble-metal (X) column position in the periodic table, this is not the case. For example, the middle-column AgZn alloy has a lower bcc order-disorder (critical) temperature than the CuZn and AuZn alloys above and below in the column. It is shown that this and other nonmonotonic behaviors can be explained in terms of a competition between atomic-size effects and X-atom d-orbital spatial extent. For example, charge-density studies and pair-potential modeling of X Zn alloys show that the effective Ag-Zn bond is significantly weaker than either the Cu-Zn or Au-Zn bond at their respective equilibrium lattice constants. Furthermore, we find that an increased atomic-core size effect initially weakens the X-Zn bonding as one goes from CuZn to AgZn, but then the larger d-orbital spatial extent for higher principal quantum numbers becomes a more dominant effect and increases the bonding from AgZn to AuZn. In conclusion, this study is focused on the highly symmetric cubic high-temperature phases, where relative bond-strength magnitudes should be far more important than any bond-directionality effects; the lattice parameters, bulk moduli, elastic constants, Debye temperatures, heats of formation, and order-disorder temperatures for the bcc phases of the three X Zn alloys are calculated and compared with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2016

On the Connection Between Microbursts and Nonlinear Electronic Structures in Planetary Radiation Belts

Using a dynamical-system approach, we have investigated the efficiency of large-amplitude whistler waves for causing microburst precipitation in planetary radiation belts by modeling the microburst energy and particle fluxes produced as a result of nonlinear wave-particle interactions. We show that wave parameters, consistent with large amplitude oblique whistlers, can commonly generate microbursts of electrons with hundreds of keV-energies as a result of Landau trapping. Relativistic microbursts (greater than 1 MeV) can also be generated by a similar mechanism, but require waves with large propagation angles Theta (sub k)B greater than 50 degrees and phase-speeds v(sub phi) greater than or equal to c/9. Using our result for precipitating density and energy fluxes, we argue that holes in the distribution function of electrons near the magnetic mirror point can result in the generation of double layers and electron solitary holes consistent in scales (of the order of Debye lengths) to nonlinear structures observed in the radiation belts by the Van Allen Probes. Our results indicate a relationship between nonlinear electrostatic and electromagnetic structures in the dynamics of planetary radiation belts and their role in the cyclical production of energetic electrons (E greater than or equal to 100 keV) on kinetic timescales, which is much faster than previously inferred.

plasmas↗
NASA NTRS2005

NASA GRC and MSFC Space-Plasma Arc Testing Procedures

Tests of arcing and current collection in simulated space plasma conditions have been performed at the NASA Glenn Research Center (GRC) in Cleveland, Ohio, for over 30 years and at the Marshall Space flight Center (MSFC) for almost as long. During this period, proper test conditions for accurate and meaningful space simulation have been worked out, comparisons with actual space performance in spaceflight tests and with real operational satellites have been made, and NASA has achieved our own internal standards for test protocols. It is the purpose of this paper to communicate the test conditions, test procedures, and types of analysis used at NASA GRC and MSFC to the space environmental testing community at large, to help with international space-plasma arcing testing standardization. To be discussed are: 1. Neutral pressures, neutral gases, and vacuum chamber sizes. 2. Electron and ion densities, plasma uniformity, sample sizes, and Debye lengths. 3. Biasing samples versus self-generated voltages. Floating samples versus grounded. 4. Power supplies and current limits. Isolation of samples from power supplies during arcs. Arc circuits. Capacitance during biased arc-threshold tests. Capacitance during sustained arcing and damage tests. Arc detection. Preventing sustained discharges during testing. 5. Real array or structure samples versus idealized samples. 6. Validity of LEO tests for GEO samples. 7. Extracting arc threshold information from arc rate versus voltage tests. 8 . Snapover and current collection at positive sample bias. Glows at positive bias. Kapton pyrolization. 9. Trigger arc thresholds. Sustained arc thresholds. Paschen discharge during sustained arcing. 10. Testing for Paschen discharge thresholds. Testing for dielectric breakdown thresholds. Testing for tether arcing. 11. Testing in very dense plasmas (ie thruster plumes). 12. Arc mitigation strategies. Charging mitigation strategies. Models. 13. Analysis of test results. Finally, the necessity of testing will be emphasized, not to the exclusion of modeling, but as part of a complete strategy for determining when and if arcs will occur, and preventing them from occurring in space.

Ferguson, Dale C.a↗
NASA NTRS1983

The gas-grain interaction in the interstellar medium - Thermal accommodation and trapping

The paper develops a numerical model for calculating thermal accommodation coefficients alpha sub T and trapping functions f sub t for gases incident on solid surfaces. The method is especially designed for astrophysical applications in that it treats economically and with moderate accuracy (+ or - 20%) the dependences of alpha sub T and f sub t on finite and different surface and gas temperatures for a large number of gas-surface combinations. In particular, the method is applied to the astrophysical combinations of hydrogen and helium gases incident on graphite, silicon, and ice surfaces. Graphs are presented of the dependence of alpha sub T and f sub t on interstellar gas temperatures in the range of 10 to 10,000 K and grain temperatures in the range 10 to 1000 K, assuming the current estimates of the gas-surface physical parameters such as the composition and the Debye temperature of the grain material, the repulsive range of the surface potential, and the gas-grain adsorption energy.

Burke, J. R.↗
NASA NTRS1995

The history of re-connection and the concept of the solar wind plasma with relatively small electrical conductivity

Petschek's 're-connection' model, aspiring to be universal, treated as a boundary problem meets unresolvable difficulties connected with impossibility to specify correctly boundary and initial conditions. This problem was incorrectly formulated. Hence, ineradicable logarithmic singularities occurred on the boundary surfaces. Attempts to eliminate them by incorporating the finite electrical conductivity are incorrect. This should lead to the change in the equation type, boundary condition type and in consequence to the change in solutions. Besides, the slow mode shocks cannot be driven by small internal source. As an alternative a new plasma concept is suggested. The state of fully ionized plasma in space depends completely on the entropy of the plasma heating source and on the process in which plasma is involved. The presumptive source of the solar wind creation - the induction electric field of the solar origin - has very low entropy. The state of plasma should be very far from the thermodynamic equilibrium. Debye's screening is not complete. The excitation of the powerful resonant self-consistent electric fields in plasma provides low electric conductivity. The MHD problems should be treated in frameworks of dissipative theories.

Chertkov, A. D.↗
DOE OSTI2021

A Systematic Way to Extend the Debye–Hückel Theory beyond Dilute Electrolyte Solutions

An extended Debye–Hückel theory with fourth order gradient term is developed for electrolyte solutions; namely, the electric potential φ(r) of the bulk electrolyte solution can be described by ∇ 2 φ(r) = κ 2 φ(r) + L Q 2 ∇ 4 φ(r), where the parameters κ and L Q are chosen to reproduce the first two roots of the dielectric response function of the bulk solution. Three boundary conditions for solving the electric potential problem are proposed based upon the continuity conditions of involving functions at the dielectric boundary, with which a boundary element method for the electric potential of a solute with a general geometrical shape and charge distribution is derived. Solutions for the electric potential of a spherical ion and a diatomic molecule are found and used to calculate their electrostatic solvation energies. In this work, the validity of the theory is successfully demonstrated when applied to binary as well as multicomponent primitive models of electrolyte solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗