Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice physics”

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 703 records · Page 39

Predicting the Mineral Composition of Dust Aerosols: Representing Key Processes - Part 1

Soil dust aerosols created by wind erosion are typically assigned globally uniform physical and chemical properties within Earth system models, despite known regional variations in the mineral content of the parent soil. Mineral composition of the aerosol particles is important to their interaction with climate, including shortwave absorption and radiative forcing, nucleation of cloud droplets and ice crystals, heterogeneous formation of sulfates and nitrates, and atmospheric processing of iron into bioavailable forms that increase the productivity of marine phytoplankton. Here, aerosol mineral composition is derived by extending a method that provides the composition of a wet-sieved soil. The extension accounts for measurements showing significant differences between the mineral fractions of the wetsieved soil and the emitted aerosol concentration. For example, some phyllosilicate aerosols are more prevalent at silt sizes, even though they are nearly absent at these diameters in a soil whose aggregates are dispersed by wet sieving. We calculate the emitted mass of each mineral with respect to size by accounting for the disintegration of soil aggregates during wet sieving. These aggregates are emitted during mobilization and fragmentation of the original undispersed soil that is subject to wind erosion. The emitted aggregates are carried far downwind from their parent soil. The soil mineral fractions used to calculate the aggregates also include larger particles that are suspended only in the vicinity of the source. We calculate the emitted size distribution of these particles using a normalized distribution derived from aerosol measurements. In addition, a method is proposed for mixing minerals with small impurities composed of iron oxides. These mixtures are important for transporting iron far from the dust source, because pure iron oxides are more dense and vulnerable to gravitational removal than most minerals comprising dust aerosols. A limited comparison to measurements from North Africa shows that the model extensions result in better agreement, consistent with a more extensive comparison to global observations as well as measurements of elemental composition downwind of the Sahara, as described in companion articles.

aerosols↗

The Geographic Distribution of Boulder Halo Craters at Mid-to-High Latitudes on Mars

Extensive evidence exists for ground ice at mid-to-high latitudes on Mars, including results from neutron spectroscopy [1-3], thermal properties [4-5], geomorphology [e.g., 6-9], and the in situ observations of Mars Phoenix [10]. This ground ice has been hypothesized to be emplaced diffusively and fill pores [11], or to have accumulated by ice and dust deposition that draped or mantled the terrain [7, 12]. These two processes are not mutually exclusive; both potentially have occurred on Mars [5]. One of the landforms found in areas where ground ice is common on Mars are boulder halo craters [e.g., 13-15] (Figure 1), which are topographically muted impact craters that are filled by ice-rich regolith. They are outlined by boulders that trace a circular outline of the original crater rim. Boulder halos generally have distinctly higher boulder densities than the surrounding background plains and have few boulders in their interiors. The mechanism of boulder halo crater formation is somewhat uncertain. Our working model is that an impact event occurs with sufficient size to excavate to a depth greater than the boulder-poor, ice-rich soils. Excavated boulders are deposited around the crater's rim and in its proximal ejecta. Quite rapidly [14], the crater becomes infilled by icy soil. Rather than being buried, boulders in the halo remain at the surface, perhaps be-cause they 'float' relative to finer-grained materials [14, 16]. Regardless of the details of this process, the life-time of boulders at the surface is much greater than the timescale needed to remove most of the craters' topography. Physical weathering of rocks must be greatly out-paced by crater infilling (the opposite of what is typical, e.g., on the Moon [17]). The rapidity of this infilling is easiest to understand if icy mantling material is deposited and accumulates, rather than simply being added by pore filling of soils. If this model is correct, boulder halos only form when they excavate rock-producing materials from beneath the upper surface. Thus, the distribution and size of craters that result in boulders halos may provide in-sight into the thickness of the ice-rich surface layer in different locations. Note that this thickness is necessarily that of the ice-rich layer at the time of impact, not at present. This study is an initial survey of boulder halo crater locations in the 50deg to 80degN and 50deg to 80degS latitude bands on Mars.

Ice↗

Bounding the Role of Black Carbon in the Climate System: a Scientific Assessment

Black carbon aerosol plays a unique and important role in Earth's climate system. Black carbon is a type of carbonaceous material with a unique combination of physical properties. This assessment provides an evaluation of black-carbon climate forcing that is comprehensive in its inclusion of all known and relevant processes and that is quantitative in providing best estimates and uncertainties of the main forcing terms: direct solar absorption; influence on liquid, mixed phase, and ice clouds; and deposition on snow and ice. These effects are calculated with climate models, but when possible, they are evaluated with both microphysical measurements and field observations. Predominant sources are combustion related, namely, fossil fuels for transportation, solid fuels for industrial and residential uses, and open burning of biomass. Total global emissions of black carbon using bottom-up inventory methods are 7500 Gg/yr in the year 2000 with an uncertainty range of 2000 to 29000. However, global atmospheric absorption attributable to black carbon is too low in many models and should be increased by a factor of almost 3. After this scaling, the best estimate for the industrial-era (1750 to 2005) direct radiative forcing of atmospheric black carbon is +0.71 W/sq m with 90% uncertainty bounds of (+0.08, +1.27)W/sq m. Total direct forcing by all black carbon sources, without subtracting the preindustrial background, is estimated as +0.88 (+0.17, +1.48) W/sq m. Direct radiative forcing alone does not capture important rapid adjustment mechanisms. A framework is described and used for quantifying climate forcings, including rapid adjustments. The best estimate of industrial-era climate forcing of black carbon through all forcing mechanisms, including clouds and cryosphere forcing, is +1.1 W/sq m with 90% uncertainty bounds of +0.17 to +2.1 W/sq m. Thus, there is a very high probability that black carbon emissions, independent of co-emitted species, have a positive forcing and warm the climate. We estimate that black carbon, with a total climate forcing of +1.1 W/sq m, is the second most important human emission in terms of its climate forcing in the present-day atmosphere; only carbon dioxide is estimated to have a greater forcing. Sources that emit black carbon also emit other short-lived species that may either cool or warm climate. Climate forcings from co-emitted species are estimated and used in the framework described herein. When the principal effects of short-lived co-emissions, including cooling agents such as sulfur dioxide, are included in net forcing, energy-related sources (fossil fuel and biofuel) have an industrial-era climate forcing of +0.22 (0.50 to +1.08) W/sq m during the first year after emission. For a few of these sources, such as diesel engines and possibly residential biofuels, warming is strong enough that eliminating all short-lived emissions from these sources would reduce net climate forcing (i.e., produce cooling). When open burning emissions, which emit high levels of organic matter, are included in the total, the best estimate of net industrial-era climate forcing by all short-lived species from black-carbon-rich sources becomes slightly negative (0.06 W/sq m with 90% uncertainty bounds of 1.45 to +1.29 W/sq m). The uncertainties in net climate forcing from black-carbon-rich sources are substantial, largely due to lack of knowledge about cloud interactions with both black carbon and co-emitted organic carbon. In prioritizing potential black-carbon mitigation actions, non-science factors, such as technical feasibility, costs, policy design, and implementation feasibility play important roles. The major sources of black carbon are presently in different stages with regard to the feasibility for near-term mitigation. This assessment, by evaluating the large number and complexity of the associated physical and radiative processes in black-carbon climate forcing, sets a baseline from which to improve future climate forcing estimates.

black carbon↗

Evolutionary characteristics of a bimodal grain model

A bimodal model of interstellar grains consisting of silicate cores of size approximately 0.06 micron with modified ice mantles of approximately 0.1 micron and very small particles of silicate and/or graphite of size approximately 0.005 micron is followed through various stages of evolution. Starting with a distribution of core-mantle cylinders producing average wavelength dependence of polarization and extinction, changes in physical and optical characteristics of the model are studied going first into regions of dense cloud condensation and coming out at the other end of star formation into young H2 regions. Physical justification is presented to show that not only in dense clouds must the core mantle grains be larger than normal, but also in young H2 regions.

Greenberg, J. M.↗

Navier-Stokes computations for exotic airfoils

An efficient hyperbolic grid generator with improvements for handling sharp corners and concave surfaces is combined with an efficient and accurate Navier-Stokes flow solver. This combination is applied to some rather complex two-dimensional airfoil configurations. Steady separated flow about an iced leading edge of an airfoil is presented. Unsteady viscous separated flows past an airfoil at two angles of attack with a spoiler deployed at 60 degrees are presented and compared with experiment. The spoiler computations are performed with two different topological maps of the physical domain to the computational domain. Innovative graphical techniques for both the static and unsteady display of the flow fields are presented and discussed.

Barth, T. J.↗

A Review and Reflections on the Sun-Climate Connection

The field of sun-climate is beset with an extraordinary number of numerical correlations attempting to relate various periodicities of solar activity with changes in the Earth's weather and climate. Signatures representing climatological variability have been sought for cycles as short as the solar 28-day rotational period up to Milankovich periods of thousands of years, although a majority of correlations have concentrated on the 11-year sunspot and 22-year Hale double sunspot cycles. For the shorter term, parameters including temperature, pressure, winds storm tracks, rainfall, and water levels in rivers and lakes, etc. have been correlated with solar variability. For longer periods, it has been necessary to seek more indirect evidence in ice cores, tree rings, and geologic deep sea cores. Other atmospheric parameters relating to atmospheric electricity and the global electric circuit have also been correlated in similar fashion. Unfortunately, few, if any, of this wide spectrum of numerical correlations have been associated with any viable physical explanation, making most studies in the field an exercise in numerical statistics. More recently, a few suggestions for plausible coupling processes have begun to appear. These, coupled with new and stronger correlations involving selective binning of climatological data sets have injected new life and hope to this field. An overview is given of the historical past and current perspectives, to evaluate possible avenues for defining physical linking processes in the future.

Goldberg, Richard A.↗

Seasonality of Biological and Physical Systems as Indicators of Climate Variation and Change

Evidence-based responses to climate change by society require operational and sustained information including biophysical indicator systems that provide up-to-date measures of trends and patterns against historical baselines. Two key components linking anthropogenic climate change to impacts on socio-ecological systems are the periodic inter- and intra-annual variations in physical climate systems (seasonality) and in plant and animal life cycles (phenology). We describe a set of national indicators that reflect sub-seasonal to seasonal drivers and responses of terrestrial physical and biological systems to climate change and variability at the national scale. Proposed indicators and metrics include seasonality of surface climate conditions (e.g., frost and freeze dates and durations), seasonality of freeze/thaw in freshwater systems (e.g., timing of stream runoff and durations of lake/river ice), seasonality in ecosystem disturbances (e.g., wildfire season timing and duration), seasonality in vegetated land surfaces (e.g., green-up and brown-down of landscapes), and seasonality of organismal life-history stages (e.g., timings of bird migration). Recommended indicators have strong linkages to variable and changing climates, include abiotic and biotic responses and feedback mechanisms, and are sufficiently simple to facilitate communication to broad audiences and stakeholders interested in understanding and adapting to climate change.

Seasonality↗

VIPER Rover: Flight Build and Environmental Test Status

The NASA Artemis Program plans to return humans to the Moon to stay. Extended human stays on the Moon will require substantial resources to sustain human presence, requiring continuous supplies delivered from the Earth. However, if some of the resources were indigenously available, Earth logistical requirements could be substantially reduced by “living off the land” with in-situ lunar resources. Local volatiles could be processed into propellants and human life-supporting needs, reducing risk of maintaining a permanent human presence on the Moon. LCROSS, LRO and other missions have confirmed the presence of lunar volatiles resources in polar regions, so the next step is to understand the physical distribution of those resources, as well as the scientific basis for how water got there, and why it is still there. The Volatiles Investigating Polar Exploration Resource (VIPER) is a surface mobility scientific platform, designed to spend ~100 days mapping and surveying four different Ice Stability Regions to understand the scientific nature and distribution of water and other volatiles. VIPER will also provide scientific mineralogical context of the lunar regolith, such as the presence of silicon and light metals in lunar regolith, providing a composite picture of resource availability and sustainment. This paper will discuss VIPER’s completion of the flight rover build, as well as current progress in environmental testing, preparedness for mission operations, and overall readying for launch integration with our CLPS partner. VIPER is managed within NASA’s Science Mission Directorate (SMD), utilizing the Commercial Lunar Payload Services (CLPS) lunar delivery model.

VIPER↗

Viper Rover: Flight Build and Environmental Test Status

The NASA Artemis Program plans to return humans to the Moon to stay. Extended human stays on the Moon will require substantial resources to sustain human presence, requiring continuous supplies delivered from the Earth. However, if some of the resources were indigenously available, Earth logistical requirements could be substantially reduced by “living off the land” with in-situ lunar resources. Local volatiles could be processed into propellants and human life-supporting needs, reducing risk of maintaining a permanent human presence on the Moon. LCROSS, LRO and other missions have confirmed the presence of lunar volatiles resources in polar regions, so the next step is to understand the physical distribution of those resources, as well as the scientific basis for how water got there, and why it is still there. The Volatiles Investigating Polar Exploration Resource (VIPER) is a surface mobility scientific platform, designed to spend ~100 days mapping and surveying four different Ice Stability Regions to understand the scientific nature and distribution of water and other volatiles. VIPER will also provide scientific mineralogical context of the lunar regolith, such as the presence of silicon and light metals in lunar regolith, providing a composite picture of resource availability and sustainment. This paper will discuss VIPER’s completion of the flight rover build, as well as current progress in environmental testing, preparedness for mission operations, and overall readying for launch integration with our CLPS partner. VIPER is managed within NASA’s Science Mission Directorate (SMD), utilizing the Commercial Lunar Payload Services (CLPS) lunar delivery model.

rover↗

VIPER Rover Flight Build and Environmental Test Status

The NASA Artemis Program plans to return humans to the Moon to stay. Extended human stays on the Moon will require substantial resources to sustain human presence, requiring continuous supplies delivered from the Earth. However, if some of the resources were indigenously available, Earth logistical requirements could be substantially reduced by “living off the land” with in-situ lunar resources. Local volatiles could be processed into propellants and human life-supporting needs, reducing risk of maintaining a permanent human presence on the Moon. LCROSS, LRO and other missions have confirmed the presence of lunar volatiles resources in polar regions, so the next step is to understand the physical distribution of those resources, as well as the scientific basis for how water got there, and why it is still there. The Volatiles Investigating Polar Exploration Resource (VIPER) is a surface mobility scientific platform, designed to spend ~100 days mapping and surveying four different Ice Stability Regions to understand the scientific nature and distribution of water and other volatiles. VIPER will also provide scientific mineralogical context of the lunar regolith, such as the presence of silicon and light metals in lunar regolith, providing a composite picture of resource availability and sustainment. This paper will discuss VIPER’s completion of the flight rover build, as well as current progress in environmental testing, preparedness for mission operations, and overall readying for launch integration with our CLPS partner. VIPER is managed within NASA’s Science Mission Directorate (SMD), utilizing the Commercial Lunar Payload Services (CLPS) lunar delivery model.

rover↗

The Lunar Atmosphere: History, Status, Current Problems, and Context

After decades of speculation and fruitless searches, the lunar atmosphere was first observed by Apollo surface and orbital instruments between 1970 and 1972. With the demise of Apollo in 1972, and the termination of funding for Apollo lunar ground station studies in 1977, the field withered for many years, but has recently enjoyed a renaissance. This reflowering has been driven by the discovery and exploration of sodium and potassium in the lunar exosphere by groundbased observers, the detection of metal ions derived from the Moon in interplanetary space, the possible discoveries of H2O ice at the poles of the Moon and Mercury, and the detections of tenuous atmospheres around more remote sites in the solar system, including Mercury and the Galilean satellites. In this review we summarize the present state of knowledge about the lunar atmosphere, describe the important physical processes taking place within it, and then discuss related topics including a comparison of the lunar atmosphere to other surface boundary exospheres in the solar system.

Stern, S. Alan .↗

An Efficient Monte Carlo Method for Modeling Radiative Transfer in Protoplanetary Disks

Monte Carlo methods have been shown to be effective and versatile in modeling radiative transfer processes to calculate model temperature profiles for protoplanetary disks. Temperatures profiles are important for connecting physical structure to observation and for understanding the conditions for planet formation and migration. However, certain areas of the disk such as the optically thick disk interior are under-sampled, or are of particular interest such as the snow line (where water vapor condenses into ice) and the area surrounding a protoplanet. To improve the sampling, photon packets can be preferentially scattered and reemitted toward the preferred locations at the cost of weighting packet energies to conserve the average energy flux. Here I report on the weighting schemes developed, how they can be applied to various models, and how they affect simulation mechanics and results. We find that improvements in sampling do not always imply similar improvements in temperature accuracies and calculation speeds.

radiative transfer↗

Organics and Ices in the Outer Solar System: Connections to the Interstellar Medium

The solar nebula, that aggregate of gas and dust that formed the birthplace of the Sun, planets and plethora of small bodies comprising the Solar System, originated in a molecular cloud that is thought to have spawned numerous additional stars, some with their own planets and attendant small bodies. The question of the chemical and physical reprocessing of the original interstellar materials in the solar nebula has challenged both theory and observations. The acquisition and analysis of samples of comet and asteroid solids, and a growing suite of in-situ and close-up analyses of relatively unaltered small Solar System bodies now adds critical new dimensions to the study of the origin and evolution of the early solar nebula. Better understanding the original composition of the material from which our solar nebula formed, and the processing that material experienced, will aid in formulations of chemistry that might occur in other solar systems. While we seek to understand the compositional history of planetary bodies in our own Solar System, we will inevitably learn more about the materials that comprise exoplanets and their surrounding systems.

SOLAR NEBULA↗

Evolution of large-scale plasma structures in comets: Kinematics and physics

Cometary and solar wind data from December 1985 through April 1986 are presented for the purpose of determining the solar wind conditions associated with comet plasma tail disconnection events (DE's). The cometary data are from The International Halley Watch Atlas of Large-Scale Phenomena (Brandt, Niedner, and Rahe, 1992). In addition, we present the kinematic analysis of 4 DE's, those of Dec. 13.5 and 31.2, 1985, and Feb. 21.7 and 28.7, 1986. The circumstances of these DE's clearly illustrate the need to analyze DE's in groups. In situ solar wind measurements from IMP-8, ICE, and PVO were used to construct the variation of solar wind speed, density, and dynamic pressure during this interval. Data from these same spacecraft plus Vega-1 were used to determine the time of 48 current sheet crossings. These data were fitted to heliospheric current sheet curves extrapolated from the corona into the heliosphere in order to determine the best-fit source surface radius for each Carrington rotation. Comparison of the solar wind conditions and 16 DE's in Halley's comet (the four DE's discussed in this paper and 12 DE's in the literature) leaves little doubt that DE's are associated primarily with crossings of the heliospheric current sheet and apparently not with any other property of the solar wind. If we assume that there is a single or primary physical mechanism and that Halley's DE's are representative, efforts at simulation should concentrate on conditions at current sheet crossings. The mechanisms consistent with this result are sunward magnetic reconnection and tailward magnetic reconnection, if tailward reconnection can be triggered by the sector boundary crossing.

Brandt, John C.↗

The Disruption and Demise of Periodic Comet Shoemaker-Levy 9

The impact of the fragmented comet Shoemaker-Levy 9 (SL9) into Jupiter this July promises to change our understanding of the outer solar system. More than twenty mountain-sized conglomerates of ice and rock will hit the atmosphere at approx. 50 km/s over the course of a week beginning July 16, releasing approx. 10(exp 4) to 10(exp8) megatons of energy per burst, and providing unique and perhaps pivotal clues to the properties of comets and the physics of massive atmospheres. Because the fragments will strike the far side of Jupiter, data acquisition, analysis and interpretation will be quite sensitive to the actual size and energy of the fragments. We therefore examine an event which took place two summers ago, unnoticed and unobserved: the disruption of SL9 into a "string of pearls' as it passed within the Roche limit at perijove. We first demonstrate, on the basis of timescales of tidal interaction, that the comet could not have broken into 20+ fragments through a hierarchy of brittle fracture events. Next, noting that the tidal stress was too weak to have even fragmented an uncompressed mass of freshly fallen snow, we run models for a strengthless comet held together only by self-gravity. We explore the initial size, density, and rotation. We conclude that a 4 km diameter comet (smaller if a prograde rotator) of density approx. 0.5 g/cu cm disrupts and disperses into a chain of fragments similar to Shoemaker-Levy 9, whether we begin with 21, 85, 169, 700 or 2000 sub-grains. Gravitational reaccumulation is evidently the answer, and there is no need to invoke the presence of 21 "cometesimals" as the subscale of the comet. To explain how a comet can be weaker than uncompacted snow, we show that the ring-plane crossing prior to perijove could have caused total damage. Finally, we compute the tidal stress on impactors as they approach Jupiter this July. Objects of various density are moderately distorted but not disrupted by the time they strike the planet.

Asphaug, Erik↗

An Automated, Low Mass, Low Power Drill for Acquiring Subsurface Samples of Ground Ice for Astrobiology Studies on Earth and on Mars

As a project that is part of NASA s Astrobiology Technology & Instrument Development Program (ASTID), we are developing a low mass (approx.20kg) drill that will be operated without drilling fluids and at very low power levels (approx.60 watts electrical) to access and retrieve samples from permafrost regions of Earth and Mars. The drill, designed and built as a joint effort by NASA JSC and Baker-Hughes International, takes the form of a down-hole unit attached to a cable so that it can, in principle, be scaled easily to reach significant depths. A parallel laboratory effort is being carried out at UC Berkeley to characterize the physics of dry drilling under martian conditions of pressure, temperature and atmospheric composition. Data from the UCB and JSC laboratory experiments are being used as input to a drill simulation program which is under development to provide autonomous control of the drill. The first Arctic field test of the unit is planned for May 2004. A field expedition to Eureka on Ellesmere Island in Spring 2003 provided an introduction for several team members to the practical aspects of drilling under Arctic conditions. The field effort was organized by Wayne Pollard of McGill University and Christopher McKay of NASA ARC. A conventional science drill provided by New Zealand colleagues was used to recover ground ice cores for analysis of their microbial content and also to develop techniques using tracers to track the depth of penetration of contamination from the core surface into the interior of the samples.

Briggs, G. A.↗

Evaluation of Higher-order Quadrature Schemes in Improving Computational Efficiency for Orientation-averaged Single-Scattering Properties of Nonspherical Ice Particles

We evaluate several high-order quadrature schemes for accuracy and efficacy in obtaining orientation-averaged single-scattering properties (SSPs). We use the highly efficient MIDAS to perform electromagnetic scattering calculations to evaluate the gain in efficiency from these schemes. MIDAS is shown to be superior to DDSCAT, a popular discrete dipole approximation (DDA) method. This study is motivated by the fact that quality physical precipitation retrievals rely on using accurate orientation-averaged SSPs derived from realistic hydrometeors as input to radiative transfer Models (RTMs). The DDA has been a popular choice for single-scattering calculations, due to its versatility with respect to target geometry. However, being iterative-solver-based (ISB), the most used DDA codes, e.g. DDSCAT and ADDA, must solve the scattering problem for each orientation of the target separately. As the size parameter and geometric anisotropy of the hydrometeor increase, the number of orientations needed to obtain accurate orientation-averages can increase drastically and so does the computation cost incurred by the ISB-DDA methods. MIDAS is a Direct-Solver-Based (DSB) code, its decomposition of the original large matrix with a high rank into multiple more manageable smaller matrices of lower ranks makes it much more computationally efficient while maintaining excellent accuracy. In addition, direct solvers consider all requested orientations at once, giving MIDAS further advantage over popular ISB-DDA methods. MIDAS, when combined with high-order quadrature for orientation averaging, can be greater than three orders of magnitude more efficient in obtaining RTM-ready SSPs of complex-shaped hydrometeors than existing ISB-DDA methods, with the native quadrature schemes they offer.

Ines Fenni↗

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.↗