SEARCH · Search NASA

Results for “Time dependent 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2021

Physics and Chemistry of Te and HgTe-based Ternary Melts

Historically, the theoretical treatment of the liquid phases has always been more difficult and complicated than that for solid and gas phases. A liquid has no lattice structure as crystalline solids and the atoms/molecules in the liquid can migrate through it relatively rapidly. On the other hand, it is also interacting with many other atoms/molecules so that the simplifications of the kinetic theory of gases cannot be employed. For more complicated liquids, such as the liquids of high ionicity and those containing hydrogen bonds and electric dipoles, the understanding is far from complete. At the same time, accurate information on the physics and chemistry of semiconductor melts is needed for the quantitative descriptions of the process of crystal growth from melt. The pre-crystallization phenomena in the liquid phase are critical because the properties of the grown crystals depend on the state and structure of the melt as well as the thermal history of the melt during solidification process. However, the data on the liquid phase, such as thermophysical properties of semiconductor melts are scarce, especially for the HgTe-based II-VI ternary compound semiconductors because of their high vapor pressure and extreme toxicity. Analysis of the thermophysical properties of the melt can provide information about structural transitions of the melt during the solidification process. From a broader point of view, the structure of liquids is much more complicated than the crystalline solids, especially the relaxation behavior through different thermal histories. The theory of hetero-phase fluctuations of liquids is applicable to any many-body systems including condensed-matter physics, field theory, physics of nuclear-matter, cosmology, biology and even sociology. This book summarizes the physics and chemistry from the experimental measurements and theoretical analyses of phase diagram, thermodynamic properties, density, thermal conductivity, viscosity, and electrical conductivity on the binary, pseudo-binary and ternary melts of the most advanced IR-detector material systems of HgCdTe and HgZnTe as well as the analyses of these results. The main objectives of this study are: (1) to provide the phase diagrams and thermodynamic properties of Hg-Cd-Te and Hg-Zn-Te systems through quantitatively fitting the experimental data by assuming an associated solution model for the liquid phase, (2) to experimentally measure the thermophysical properties of the Hg-Cd-Te and Hg-Zn-Te melts, including density, viscosity, electrical conductivity and thermal conductivity as functions of temperature and composition and (3) to enhance the fundamental knowledge of hetero-phase fluctuations and relaxation phenomena in the melts and extend our understanding of the solidification process in order to interpret the experimental results of crystal growth so as to improve the melt growth processes of the compound semiconductor. The physics and chemistry of Te and HgTe-based ternary melts were explored through the studies of the structural transformation during melting, the supercooling during solidification, the relaxation phenomena after rapid cooling of the melts and the metal-semiconductor transition in the melts through the analyses of electrical conductivity and Lorenz number. An in-depth study on the thermophysical properties and their time-dependent structural dynamic processes taking place in the vicinity of the solid-liquid phase transition of the narrow homogeneity range HgTe-based ternary semiconductors as well as the analysis of the homogenization process in the melt will also be presented.

Thermophysical properties↗

DOE OSTI2025

Refining the two-band model for highly compensated semimetals using thermoelectric coefficients

In studying compensated semimetals, the two-band model has proven extremely useful in capturing electrical conductivity under magnetic field, as a function of density and mobility of electronlike and holelike carriers. However, it rarely offers practical insight into magnetothermoelectric properties. Here, in this work, we report the field dependence of thermoelectric (TE) coefficients in a highly compensated semimetal NbSb 2 , where we find the Seebeck (𝑆 𝑥⁢𝑥 ) and Nernst (𝑆 𝑥⁢𝑦 ) coefficients increase quadratically and linearly with applied magnetic field, respectively. Such field dependence was predicted in previous work that studied a system of two parabolic bands, within semiclassical Boltzmann transport theory when the following two conditions are simultaneously met: 𝜔 𝑐 ⁢𝜏 ≫ 1 and tan ⁡𝜃 𝐻 ≪ 1, where 𝜔 𝑐 and 𝜏 refer to the cyclotron frequency and relaxation time, respectively, and 𝜃 𝐻 is the Hall angle. Under these conditions, we find the field dependence of the TE coefficients directly provides a relation between the electronlike (𝑛 𝑒 ) and holelike (𝑛 ℎ ) carrier densities, which in turn can be used to refine two-band model fitting. With this, we find the compensation factor ($\frac{|Δ⁢𝑛|}{𝑛_{𝑒}}$, where Δ⁢𝑛 = 𝑛 𝑒 −𝑛 ℎ ) of NbSb 2 is two orders of magnitude smaller than what was found in unrestricted fitting, resulting in a larger saturation-field scale for magnetoresistance. Within the same framework of the semiclassical theory, we also deduce that the thermoelectric Hall angle tan⁡ 𝜃 𝛾 = $\frac{𝑆_{𝑥⁢𝑦}}{𝑆_{𝑥⁢𝑥}}$ can be expressed as ($\frac{Δ⁢𝑛}{𝑛_{𝑒}}$ ×𝜔 𝑐 ⁢𝜏) −1 , which serves as a parameter to predict the degree of compensation. Our findings offer crucial insights into identifying empirical conditions for field-induced enhancement of TE performance and into engineering-efficient thermoelectric devices based on semimetallic materials.

electrical conductivity↗

DOE OSTI2025

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗

DOE OSTI2024

Electronic structure and optical properties of halide double perovskites from a Wannier-localized optimally-tuned screened range-separated hybrid functional

Halide double perovskites are a chemically diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and ab initio many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here, in this study, we apply the WOT-SRSH functional to five halide double perovskites and compare the results with those obtained from other known functionals and previous GW calculations. We also use the approach as a starting point for GW calculations and we compute the band structures and optical absorption spectrum for Cs 2 AgBiBr 6 , using both time-dependent DFT and the GW-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Thermodynamics of Liquid Uranium from Atomistic and Ab Initio Modeling

We present thermodynamic properties for liquid uranium obtained from classical molecular dynamics (MD) simulations and the first-principles theory. The coexisting phases method incorporated within MD modeling defines the melting temperature of uranium in good agreement with the experiment. The calculated melting enthalpy is in agreement with the experimental range. Classical MD simulations show that ionic contribution to the total specific heat of uranium does not depend on temperature. The density of states at the Fermi level, which is a crucial parameter in the determination of the electronic contribution to the total specific heat of liquid uranium, is calculated by ab initio all electron density functional theory (DFT) formalism applied to the atomic configurations generated by classical MD. The calculated specific heat of liquid uranium is compared with the previously calculated specific heat of solid γ-uranium at high temperatures. The liquid uranium cannot be supercooled below T sc ≈ 800 K or approximately about 645 K below the calculated melting point, although, the self-diffusion coefficient approaches zero at T D ≈ 700 K. Uranium metal can be supercooled about 1.5 times more than it can be overheated. The features of the temperature hysteresis are discussed.

36 MATERIALS SCIENCE↗

NASA NTRS2005

Shock Waves in a Bose-Einstein Condensate

A paper presents a theoretical study of shock waves in a trapped Bose-Einstein condensate (BEC). The mathematical model of the BEC in this study is a nonlinear Schroedinger equation (NLSE) in which (1) the role of the wave function of a single particle in the traditional Schroedinger equation is played by a space- and time-dependent complex order parameter (x,t) proportional to the square root of the density of atoms and (2) the atoms engage in a repulsive interaction characterized by a potential proportional to | (x,t)|2. Equations that describe macroscopic perturbations of the BEC at zero temperature are derived from the NLSE and simplifying assumptions are made, leading to equations for the propagation of sound waves and the transformation of sound waves into shock waves. Equations for the speeds of shock waves and the relationships between jumps of velocity and density across shock fronts are derived. Similarities and differences between this theory and the classical theory of sound waves and shocks in ordinary gases are noted. The present theory is illustrated by solving the equations for the example of a shock wave propagating in a cigar-shaped BEC.

Kulikov, Igor↗

DOE OSTI2025

Operando probing dynamic migration of copper carbonyl during electrocatalytic CO2 reduction

Single crystals and shape-controlled nanocrystals are well known to exhibit facet-dependent catalytic properties. However, few studies have investigated how those nanocrystals evolve and (de)activate during reactions, calling for the development of nanoscale time-resolved operando methods. In this context, we have designed Cu nanocubes as a model system to elucidate the underlying driving force of dynamic nanocatalyst reconstruction during the CO2 reduction reaction (CO2RR). Operando electrochemical liquid-cell scanning transmission electron microscopy (EC-STEM) and synchrotron-based X-ray spectroscopy reveal the size- and potential-dependent complete transformation from (100)-oriented Cu@Cu2O nanocubes to polycrystalline metallic Cu nanograins under CO2RR conditions. In addition, machine learning-assisted operando four-dimensional STEM reveals that large Cu nanograins derived from nanocubes form mainly crystalline domains, while their smaller counterparts are more amorphous due to faster evolution kinetics. In situ Raman spectroscopy and density functional theory calculations suggest that CO drives the ejection of single Cu atoms, resulting in few-nanometre Cu clusters and the surface migration of highly mobile copper carbonyl (Cu–CO) species. Combined, these multimodal operando methods and theoretical approaches pave the way for understanding the complex structural evolution of energy-related nanocatalysts under electrochemical conditions.

Yang, Yao↗

NASA NTRS1983

Ultrasonic wave propagation in two-phase media: Spherical inclusions

The scattering theory, recently developed via the extended method of equivalent inclusion, is used to study the propagation of time-harmonic waves in two-phase media of elastic matrix with randomly distributed elastic spherical inclusion materials. The elastic moduli and mass density of the composite medium are determined as functions of frequencies when given properties and concentration of the spheres and the matrix. Velocity and attenuation of ultrasonic waves in two-phase media are determined for cases of distributed spheres and localized damage. An averaging theorem that requires the equivalence of the strain energy and the kinetic energy between the effective medium and the original matrix with spherical inhomogeneities is employed to derive the effective moduli and mass density. The functional dependency of these quantities upon frequencies and concentration provides a method of data analysis in ultrasonic evaluation of material properties. Numerical results or moduli, velocity and/or attenuation as functions of concentration of inclusion material, or porosity, are graphically displayed.

Fu, L. S.↗

NASA NTRS1984

Ultrasonic wave propagation in two-phase media - Spherical inclusions

Fu, L. S.↗

DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗

DOE OSTI2024

An Impedance Humidity Sensor Based on CVD Grown WSe 2 2D Films

2D materials-based planar devices have been sought after for gas and humidity sensor applications owing to their high sensitivity due to their ultrathin nature and increased surface area. Herein, 2D WSe 2 films-based planner devices were fabricated to evaluate their performance for potential application in humidity sensors. The structure and morphology of the WSe 2 sensing films were studied by Raman spectroscopy, photoluminescence, optical microscopy, and scanning electron microscopy (SEM) techniques. The relative humidity-dependent sensing performances of the 2D WSe 2 films were evaluated using current–voltage (I–V), relative impedance, and complex total impedance (electrochemical impedance spectroscopy, EIS) techniques. This sensor showed a very stable and repeatable performance investigated over the period of 18 months. The response and recovery times of the WSe 2 -based impedance sensors were found to be 40 and 30 s, respectively. The WSe 2 sensor exhibits a 2 ± 1% minimum hysteresis at a lower humidity level (50% RH) and around 8 ± 2% maximum hysteresis at relative humidity level (70% RH). This value of hysteresis is in the acceptable range during the initial investigation of any material to explore its potential for application in humidity sensors. The density functional theory (DFT) analysis was performed with the interaction from water molecules of the WSe 2 sample for a better understanding of the sensing mechanism of devices.

2D- WSe2↗

NASA NTRS1986

Chemical evolution in spiral and irregular galaxies

A brief review of models of chemical evolution of the interstellar medium in our galaxy and other galaxies is presented. These models predict the time variation and radial dependence of chemical composition in the gas as function of the input parameters; initial mass function, stellar birth rate, chemical composition of mass lost by stars during their evolution (yields), and the existence of large scale mass flows, like infall from the halo, outflow to the intergalactic medium or radial flows within a galaxy. At present there is a considerable wealth of observational data on the composition of HII regions in spiral and irregular galaxies to constrain the models. Comparisons are made between theory and the observed physical conditions. In particular, studies of helium, carbon, nitrogen and oxygen abundances are reviewed. In many molecular clouds the information we have on the amount of H2 is derived from the observed CO column density, and a standard CO/H2 ratio derived for the solar neighborhood. Chemical evolution models and the observed variations in O/H and N/O values, point out the need to include these results in a CO/H2 relation that should be, at least, a function of the O/H ratio. This aspect is also discussed.

Torres-Peimbert, S.↗

DOE OSTI2026

Identifying Strain Stacking Boundaries between Multiphase Domains in Atomically Thin Two-Dimensional Magnets

Stacking engineering of van der Waals materials is an important strategy to control the materials’ properties, such as electronic correlations, ferroelectricity, and layer-dependent two-dimensional magnetism. A timely testbed for the study of the latter is atomically thin chromium trihalides (CrX 3 , X = Cl, Br, I). Notably, by understanding the sliding mechanism between different stacking sequences, control of the stacking arrangement, and thus magnetic properties in CrX 3 , can be achieved. Such insight, however, is currently lacking. Here, in this study, advanced electron microscopy methods are used to identify multiple stacking sequences corresponding to different bulk phases in atomically thin CrX 3 (X = Cl and Br) down to bilayer thickness and with lateral domain sizes as small as tens of nanometers. Indications of nanometer scale transitions and interactions at the stacking boundaries are found, including a universally preferred sliding direction that is consistent with density functional theory calculations and the strain fields at lateral heterostructure boundaries. This study demonstrates the necessity to consider local stacking structures when interpreting averaged magnetic properties measured with macroscale probes. Additionally, the preferred sliding direction insight from this study provides a strategy to control stacking sequence in atomically thin CrX 3 samples during the exfoliation and sample fabrication process.

DFT calculations↗

DOE OSTI2025

Multireference diffusion Monte Carlo reaches 2D materials

Abstract Quantum confinement in 2D materials strongly enhances electronic correlation effects. Therefore, predicting the properties of these unique materials, with both a high level of accuracy and computational efficiency, without relying on adjustable parameters or functionals, remains an outstanding theoretical challenge. The majority of theoretical studies are based on the approximations of density functional theory (DFT). The reliability of DFT predictions are heavily dependent on the choice of an approximated exchange-correlation functional. Here, we estimate the magnitude of impact of correlation on the total energy for the quintessential 2D material, graphene, by performing and comparing state-of-the-art selected CI and quantum Monte Carlo extrapolated calculations for a single unit cell at the$$\Gamma$$point. We demonstrate that Self-Healing Diffusion Monte Carlo (SHDMC) obtains a very compact, but high-quality wavefunction for this system that lacks the strong basis set dependence displayed by state of the art quantum chemistry methods. The SHDMC wavefunction is of higher quality compared to that obtained from sCI, in the same orbital basis, while being$$\sim$$ 1000 times smaller in terms of determinant count compared to sCI. We also demonstrate that extrapolating SHDMC results to the infinite determinant limit compares extremely well with complete basis set extrapolated sCI. Our work paves the way for future validation of SHDMC applied to challenging 2D materials.

Science & Technology - Other Topics↗

DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

DOE OSTI2025

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1978

A magnetohydrodynamic theory for the lunar response to time variations in a spatially uniform ambient magnetic field

An analytic theory based on principles of magnetohydrodynamics is derived for interpreting the lunar magnetic response to magnetic field fluctuations in the lobes of the geomagnetic tail. Only the case of magnetic perturbations parallel to a spatially uniform ambient magnetic field is treated. The theoretical frequency-dependent MHD transfer function depends on both the lunar electrical conductivity profile and the Alfven speed of the tail lobe plasma. For circular frequencies much less than the Alfven speed/lunar radius, the MHD transfer function reduces to the vacuum transfer function, while at higher frequencies the finite propagation speed of magnetic disturbances affects the response. These characteristics of the MHD transfer function are in qualitative agreement with observation. A carefully selected set of magnetic field transient events could be used together with this theory to infer lunar electrical conductivity profiles and the tail lobe plasma density (via its dependence on the Alfven speed).

Hood, L. L.↗