Computational Design of High-Entropy Rare Earth Disilicates as Next-Generation Thermal/Environmental Barrier Coatings
For presentation at Materials Science Technology 2023, Columbus, OH, October 1-4, 2023.
Hao, Shiqiang↗
SEARCH · Search NASA
Results for “DEBYE FUNCTION”
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.
Giant Dipole Moments: Remarkable Effects Mono‐, Di‐, and Tri‐ Hydrated 5,6‐Diaminobenzene‐1,2,3,4‐Tetracarbonnitrile
The molecule 5,6-diaminobenzene-1,2,3,4-tetracarbonnitrile (MOI) was first synthesized by Müllen and coworkers in 2016 and boasts an ultrastrong dipole moment of $14.1\pm 0.7$ Debye in THF. Gas phase DFT computations do not fully reflect this ultrastrong dipole moment, demonstrating the role of solvent in increasing this dipole moment. Here, we investigate the effect of solvent molecule position on the dipole moment of this species, computationally examining systems with giant dipole moments. These systems are optimized in the gas phase with the B3LYP functional, employing the aug-cc-pVTZ and def2-TZVP basis sets, as well as the B3LYP-D3BJ/aug-cc-pVTZ functional in Orca. Single point DLPNO-CCSD/aug-cc-pVDZ results were obtained from Orca and Psi4, as well as DLPNO-CCSD(T)/CBS information from Psi4. Additionally, these are compared to the dipole moments of di- and tri-hydrated systems, and the SMD models for THF and water at the B3LYP/aug-cc-pVTZ level of theory. The dissociation energies, HOMO-LUMO energy gaps, and dipole moments are presented. These metrics show the nh1nh1′ THF system boasts the largest dissociation energy and dipole moment of the singly solvated systems, due to its strong hydrogen bonding. The importance of solvent placement is highlighted and may guide the synthesis of macromolecules or organic frameworks incorporating the MOI or MOI-like subunits. Remarkably, a single solvent molecule provides a good model for the difference between the gas phase and solvated species. The predicted gas phase dipole moments computed with B3LYP/aug-cc-pVTZ for the MOI, its monohydrated complex, dihydrated complex, and its trihydrated complex are 9.6, 14.2, 16.0, and 16.8 Debye, respectively.
dipole↗
Theoretical study of the freezing of polystyrene sphere suspensions
A theoretical study of melting in aqueous suspensions of polystyrene spheres is presented. The Helmholtz free energies of the liquid and solid phase of the suspensions are calculated as functions of sphere number density and electrolyte concentrations. The results tend to show that the freezing curve of such suspensions is that of a conventional classical liquid with repulsive short-range interactions. The interactions can be treated satisfactorily within a slightly modified Debye-Hueckel approximation. The modifications include size correction and a correction for nonlinear screening. The results are confirmed by analogous calculations for the solid phase, and for the line along which liquid and solid free energies are equal.
Shih, W.-H.↗
Artifact identification in X-ray diffraction data using machine learning methods
In situ synchrotron high-energy X-ray powder diffraction (XRD) is highly utilized by researchers to analyze the crystallographic structures of materials in functional devices ( e.g. battery materials) or in complex sample environments ( e.g. diamond anvil cells or syntheses reactors). An atomic structure of a material can be identified by its diffraction pattern along with a detailed analysis of the Rietveld refinement which yields rich information on the structure and the material, such as crystallite size, microstrain and defects. For in situ experiments, a series of XRD images is usually collected on the same sample under different conditions ( e.g. adiabatic conditions) yielding different states of matter, or is simply collected continuously as a function of time to track the change of a sample during a chemical or physical process. In situ experiments are usually performed with area detectors and collect images composed of diffraction patterns. For an ideal powder, the diffraction pattern should be a series of concentric Debye–Scherrer rings with evenly distributed intensities in each ring. For a realistic sample, one may observe different characteristics other than the typical ring pattern, such as textures or preferred orientations and single-crystal diffraction spots. Textures or preferred orientations usually have several parts of a ring that are more intense than the rest, whereas single-crystal diffraction spots are localized intense spots owing to diffraction of large crystals, typically >10 µm. In this work, an investigation of machine learning methods is presented for fast and reliable identification and separation of the single-crystal diffraction spots in XRD images. The exclusion of artifacts during an XRD image integration process allows a precise analysis of the powder diffraction rings of interest. When it is trained with small subsets of highly diverse datasets, the gradient boosting method can consistently produce high-accuracy results. The method dramatically decreases the amount of time spent identifying and separating single-crystal diffraction spots in comparison with the conventional method.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Carbon-nitrogen complex as an alternative to the 𝑇 center in Si
A first-principles study of a carbon-nitrogen (CN) impurity complex in silicon is presented as an isoelectronic alternative to the T center [(CCH) Si ]. The latter has been pursued for applications in quantum information science, yet its sensitivity to the presence of hydrogen is still problematic. Our proposed complex has no hydrogen, thereby eliminating this issue. First, we show that the CN complex is stable against decomposition into substitutional and interstitial defects. Next, we show that due to being isoelectronic to the T center, the CN complex has a similar electronic structure, and therefore could be used in similar applications. We assess several low-energy configurations of the CN complex, finding (CN) Si to be stable and have the largest Debye-Waller factor. We predict a zero-phonon line (ZPL) of 828 meV (in the telecom S-band) and a radiative lifetime of 4.2 μs, comparable to the T center. Due to the presence of a bound exciton, choice of the exchange-correlation functional and also supercell-size scaling of the ZPL and transition dipole moment require special scrutiny; we rigorously justify our extrapolation schemes that allow computing values in the dilute limit.
36 MATERIALS SCIENCE↗
Simple phenomenological modeling of transition-region capacitance of forward-biased p-n junction diodes and transistor diodes
The derivation of a simple expression for the capacitance C(V) associated with the transition region of a p-n junction under a forward bias is derived by phenomenological reasoning. The treatment of C(V) is based on the conventional Shockley equations, and simpler expressions for C(V) result that are in general accord with the previous analytical and numerical results. C(V) consists of two components resulting from changes in majority carrier concentration and from free hole and electron accumulation in the space-charge region. The space-charge region is conceived as the intrinsic region of an n-i-p structure for a space-charge region markedly wider than the extrinsic Debye lengths at its edges. This region is excited in the sense that the forward bias creates hole and electron densities orders of magnitude larger than those in equilibrium. The recent Shirts-Gordon (1979) modeling of the space-charge region using a dielectric response function is contrasted with the more conventional Schottky-Shockley modeling.
Lindholm, F. A.↗
Low frequency dielectric relaxation in boracites
In order to elucidate the nature of the imperfections which adversely affect pyroelectric processes in boracites, the thermal and dispersive characteristics of the low frequency dielectric response in iron-iodide boracite (Fe3B7O13I) and copper-chloride boracite (Cu3B7O13Cl) have been investigated. These characteristics were measured as a function of crystallographic orientation and applied field in both the ferroelectric and paraelectric states. The low frequency dielectric relaxation of 100 line oriented multi-domain copper-chloride boracite clearly indicates the dipole nature of the lattice imperfections. The activation energies calculated from a noninteracting Debye model, are 0.53 eV in the ferroelectric phase and 0.10 eV in the paraelectric phase.
Kim, Q.↗
Formation of the Electric Double Layer and its Effects on Moving Bodies in a Space Plasma Environment
In this paper we solve the self-consistent Vlasov and Poisson equations by a numerical method to determine the local distribution function of the ion and the electron, within a thin layer near the moving body, respectively. Using these ion and electron distributions, the number density for the ions and electrons are determined, such that, the electric potential is obtained within this thin layer (i.e., measured by Debye length). Numerical results are presented for temporal evolution of the electron and ion density and its corresponding electric potential within the layer which shows the formation of electric double layer and its structures. From these numerical results, we are able to determine the maximum conditions of the electric potential, it may create satellite anomaly.
Yang, Qianli↗
Thermal Conductivity of β-Phase Ga 2 O 3 and (Al x Ga 1– x ) 2 O 3 Heteroepitaxial Thin Films
Heteroepitaxy of β-phase gallium oxide (β-Ga 2 O 3 ) thin films on foreign substrates shows promise for the development of next-generation deep ultraviolet solar blind photodetectors and power electronic devices. In this work, the influences of the film thickness and crystallinity on the thermal conductivity of ($\bar{2}01$)-oriented β-Ga 2 O 3 heteroepitaxial thin films were investigated. Unintentionally doped β-Ga 2 O 3 thin films were grown on c-plane sapphire substrates with off-axis angles of 0° and 6° toward $\langle$$11\bar{2}0$$\rangle$ via metal–organic vapor phase epitaxy (MOVPE) and low-pressure chemical vapor deposition. The surface morphology and crystal quality of the β-Ga 2 O 3 thin films were characterized using scanning electron microscopy, X-ray diffraction, and Raman spectroscopy. The thermal conductivities of the β-Ga 2 O 3 films were measured via time-domain thermoreflectance. The interface quality was studied using scanning transmission electron microscopy. The measured thermal conductivities of the submicron-thick β-Ga 2 O 3 thin films were relatively low as compared to the intrinsic bulk value. The measured thin film thermal conductivities were compared with the Debye–Callaway model incorporating phononic parameters derived from first-principles calculations. The comparison suggests that the reduction in the thin film thermal conductivity can be partially attributed to the enhanced phonon-boundary scattering when the film thickness decreases. They were found to be a strong function of not only the layer thickness but also the film quality, resulting from growth on substrates with different offcut angles. Growth of β-Ga 2 O 3 films on 6° offcut sapphire substrates was found to result in higher crystallinity and thermal conductivity than films grown on on-axis c-plane sapphire. However, the β-Ga 2 O 3 films grown on 6° offcut sapphire exhibit a lower thermal boundary conductance at the β-Ga 2 O 3 /sapphire heterointerface. In addition, the thermal conductivity of MOVPE-grown ($\bar{2}01$)-oriented β-(Al x Ga 1–x ) 2 O 3 thin films with Al compositions ranging from 2% to 43% was characterized. Because of phonon-alloy disorder scattering, the β-(Al x Ga 1–x ) 2 O 3 films exhibit lower thermal conductivities (2.8–4.7 W/m∙K) than the β-Ga 2 O 3 thin films. The dominance of the alloy disorder scattering in β-(Al x Ga 1–x ) 2 O 3 is further evidenced by the weak temperature dependence of the thermal conductivity. This work provides fundamental insight into the physical interactions that govern phonon transport within heteroepitaxially grown β-phase Ga 2 O 3 and (Al x Ga 1–x ) 2 O 3 thin films and lays the groundwork for the thermal modeling and design of β-Ga 2 O 3 electronic and optoelectronic devices.
36 MATERIALS SCIENCE↗
A prediction of the thermodynamic, thermophysical, and mechanical properties of CrTaO 4 from first principles
It is reported that the self-forming CrTaO 4 oxide scale can protect refractory high-entropy alloys from oxidation, superior to Cr 2 O 3 . In this paper, the phase stability, mechanical, and thermal properties of three polymorphous phases of CrTaO 4 are systematically investigated from first-principles density functional theory calculations. The mechanical properties predicted using the strain–energy methods indicated that all three phases are mechanically stable. The temperature dependence of elastic constants and polycrystalline moduli of three phases demonstrated the thermal softening as temperature increase. The Helmholtz-free energies as a function of volume and temperature are derived from phonon dispersions within the quasi-harmonic approximation at six strained volumes. The calculated apparent bulk coefficients of thermal expansion of these three phases are evaluated, the highest value approximately 13.4× 10 –6 K –1 within a temperature range of 500–2000 K for the rutile I4 1 md phase. The lattice thermal conductivity calculated by the Debye–Callaway model suggested that the rutile type I4 1 md phase has the lowest value of approximately 2.1 W/m/K at 1800 K. The other two phases, C2/m and P2/c, exhibit higher values due to relatively lower Grüneisen parameters and larger phonon velocities. The melting point of CrTaO 4 is predicted to be between 1975 and 2449 K using ab initio molecular dynamics simulations. Finally, this work provides a comprehensive theoretical understanding of the thermodynamic, mechanical, and thermal properties for the new material CrTaO 4 and serves as an example of a viable computational design strategy for improved oxidation resistance of refractory alloys at high temperatures.
36 MATERIALS SCIENCE↗
Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field
A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.
Physics↗
Rainbow Scattering by a Coated Sphere
We examine the behavior of the first-order rainbow for a coated sphere by using both ray theory and Aden-Kerker wave theory as the radius of the core alpha(sub 12) and the thickness of the coating beta are varied. As the ratio beta/alpha(sub 12) increases from 10(sup -4) to 0.33, we find three classes of rainbow phenomena that cannot occur for a homogeneous-sphere rainbow. For beta/alpha(sub 12) approx less than 10(sup -3), the rainbow intensity is an oscillatory function of the coating thickness, for beta/alpha(sub 12) approx. 10(sup -2), the first-order rainbow breaks into a pair of twin rainbows, and for beta/alpha(sub 12) approx. 0.33, various rainbow-extinction transitions occur. Each of these effects is analyzed, and their physical interpretations are given. A Debye series decomposition of coated-sphere partial-wave scattering amplitudes is also performed and aids in the analysis.
Lock, James A.↗
Superconducting compounds and alloys research
Resistivity measurements as a function of temperature were performed on alloys of the binary material system In sub(1-x) Bi sub x for x varying between 0 and 1. It was found that for all single-phase alloys (the pure elements, alpha-In, and the three intermetallic compounds) at temperatures sufficiently above the Debye-temperature, the resistivity p can be expressed as p = a sub o T(n), where a sub o and n are composition-dependent constants. The same exponential relationship can also be applied for the sub-system In-In2Bi, when the two phases are in compositional equilibrium. Superconductivity measurements on single and two-phase alloys can be explained with respect to the phase diagram. There occur three superconducting phases (alpha-In, In2Bi, and In5Bi3) with different transition temperatures in the alloying system. The magnitude of the transition temperatures for the various intermetallic phases of In-Bi is such that the disappearance or occurrence of a phase in two component alloys can be demonstrated easily by means of superconductivity measurements.
Otto, G.↗
Complete Equations of State for Cyclotetramethylene Tetranitramine
Complete equations of state for the volumetric deformation of the β-polymorph of cyclotetramethylene tetranitramine (β-HMX) have been derived from its Helmholtz free energy. Dispersion-corrected density functional theory calculations were used to compute the dependence of the frequencies of the vibrational normal modes on specific volume. Here, the normal mode frequencies were used directly to generate both a tabular equation of state and an approximate model where the sum over normal mode frequencies is replaced by a set of five Debye models.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
Lattice dynamics and thermodynamics for δ -plutonium from density functional theory
Here, we present results from density functional theory (DFT) calculations of the lattice dynamics (phonons) and thermodynamics for δ-phase plutonium. The fully relativistic electronic structure is calculated assuming a three-dimensional noncollinear magnetic structure in conjunction with DFT and the general gradient approximation for the electron exchange and correlation interactions. The electronic-structure model is further enhanced by addressing strong orbital-orbital coupling via the conventional orbital-polarization (OP) scheme as has been successfully done for plutonium. The temperature dependence of the phonons is calculated within the self-consistent ab initio lattice dynamics approach. The obtained phonons compare very well with measurements although a modest overestimation of the transverse L-point [ξξξ] phonon is acknowledged. Calculated thermal vibration amplitudes and the associated Debye-Waller temperatures are close to experiments. Lattice, electronic, and magnetic contributions to the heat capacity are predicted and consistent to a few percent with that deduced from experimental data. Good agreement is only achieved when a magnetic contribution to the specific heat is recognized. The parameter-free DFT+OP electronic model is thus capable of predicting phonon properties and thermodynamic behavior of δ-phase plutonium rather accurately.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Low-temperature thermostatics of face-centered-cubic metallic hydrogen
The thermostatic properties of a high-symmetry phase of metallic hydrogen with atomic sphere radius between 0.1 and 1.5 bohr are studied, with special emphasis accorded to electronic screening and quantum proton motion. The electron-proton and proton-proton interactions receive a perturbation treatment based on the Singwi dielectric function, while the proton motion is handled by self-consistent harmonic approximation. Quantum behavior is found to be less pronounced than expected, and nuclear magnetism is absent. The phonon spectrum is, however, affected by screening and large proton motion. The zero-point vibrational energy and the superconducting critical temperature are below previous estimates. The crystalline-defect formation energies are a few times the Debye energy, which implies that defects contribute significantly to melting at the lower particle densities.
Caron, L. G.↗
Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate
We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTiO3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗