Search NASA⌕ Search

SEARCH · Search NASA

Results for “Point pattern analysis”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Measurement of Surface Deformation Related to the December 2018 Mt. Etna Eruption Using Time-Series Interferometry and Magma Modeling for Hazard Zone Mapping

Mount Etna has erupted several times since it was first formed. Recently, Mount Etna began erupting again over 24–27 December 2018. Because it erupts frequently, Mount Etna should be observed on a frequent basis. From June 2018 to October 2019, 34 and 56 interferometric synthetic aperture radar (InSAR) images were acquired from the ascending and descending tracks of the Sentinel-1 satellite, respectively. We employed the Stanford Method for Persistent Scatterers (StaMPS) and a refined small baseline subset (SBAS) InSAR method to produce a surface deformation time-series map. In the time-series analysis, the phase signal remained unaltered with time. The Okada model was then applied to the result to generate a modeled interferogram, and the Q-LavHA program was run to generate a lava flow prediction model. A direct comparison of the results showed that Persistent Scatterers Interferometry (PSI)-StaMPS and the refined SBAS technique were comparable in terms of the displacement pattern, with slightly different velocity values obtained for individual points. In particular, a velocity range of −25 to 21 cm/yr was obtained from PSI-StaMPS, whereas a range of −30 to 25 cm/yr was obtained from the refined SBAS method. Upon computation of the vertical and east-west displacement components based on ascending and descending track data using both methods, deformation velocities of 51.5 and 52.5 cm/yr in the westerly direction on the western flank of Mount Etna were obtained from PSI-StaMPS and the refined SBAS method, respectively, whereas on the eastern flank, deformation toward the east was estimated to occur at a velocity of 50.1 or 54.2 cm/yr, respectively. PSIStaMPS estimated a vertical deformation velocity of −5.3 to 18.3 cm/yr, whereas the refined SBAS method produced a velocity range of approximately −7 to 19 cm/yr. The interferogram obtained via Okada modeling showed two fault sources in the 2018 Mount Etna eruption and a total volume change of approximately 12.39 × 10 6 m 3 . From the modeling results, a lava flow prediction model was generated using the Q-LavHA program. The approaches described in this study can be used by government officials, authorities, and other decision-makers to monitor and assess the risk of volcanic activity in the region.

Suci Ramayanti↗
NASA NTRS2008

Analysis of a Dynamic Multi-Track Airway Concept for Air Traffic Management

The Dynamic Multi-track Airways (DMA) Concept for Air Traffic Management (ATM) proposes a network of high-altitude airways constructed of multiple, closely spaced, parallel tracks designed to increase en-route capacity in high-demand airspace corridors. Segregated from non-airway operations, these multi-track airways establish high-priority traffic flow corridors along optimal routes between major terminal areas throughout the National Airspace System (NAS). Air traffic controllers transition aircraft equipped for DMA operations to DMA entry points, the aircraft use autonomous control of airspeed to fly the continuous-airspace airway and achieve an economic benefit, and controllers then transition the aircraft from the DMA exit to the terminal area. Aircraft authority within the DMA includes responsibility for spacing and/or separation from other DMA aircraft. The DMA controller is responsible for coordinating the entry and exit of traffic to and from the DMA and for traffic flow management (TFM), including adjusting DMA routing on a daily basis to account for predicted weather and wind patterns and re-routing DMAs in real time to accommodate unpredicted weather changes. However, the DMA controller is not responsible for monitoring the DMA for traffic separation. This report defines the mature state concept, explores its feasibility and performance, and identifies potential benefits. The report also discusses (a) an analysis of a single DMA, which was modeled within the NAS to assess capacity and determine the impact of a single DMA on regional sector loads and conflict potential; (b) a demand analysis, which was conducted to determine likely city-pair candidates for a nationwide DMA network and to determine the expected demand fraction; (c) two track configurations, which were modeled and analyzed for their operational characteristic; (d) software-prototype airborne capabilities developed for DMA operations research; (e) a feasibility analysis of key attributes in the concept design; (f) a near-term, transitional application of the DMA concept as a proving ground for new airborne technologies; and (g) conclusions. The analysis indicates that the operational feasibility of a national DMA network faces significant challenges, especially for interactions between DMAs and between DMA and non-DMA traffic. Provided these issues are resolved, sectors near DMAs could experience significant local capacity benefits.

Wing, David J.↗
NASA NTRS2000

Interannual Atmospheric Oscillations and Their Gravitational Effects

Past studies have examined the effects of the interannual atmospheric oscillations (IAO), or often called teleconnection patterns, on the Earth's rotation, such as El Nino/Southern Oscillation (ENSO), the quasi-biennial oscillation, and to a lesser extent the North Atlantic Oscillation (NAO). The present study focuses on IAO's effects on the gravitational field and geocenter; the purpose is to be able to identify and isolate the contributions of each prominent IAO in relation to the total contribution of the atmosphere, in terms of their magnitudes, their geographical patterns, and their interannual time history. We use the 40-year NCEP reanalysis of the monthly, global atmospheric surface pressure field as our basic data set. The method we apply to isolate the IAOs is the empirical orthogonal function (EOF) decomposition which is widely used in meteorological investigations. We do the EOF analysis first on IAOs' seasonal signals (by "collapsing" the 40-year series into 12 mean-months for each grid point) and obtain estimates for their respective contributions. Then we remove these seasonal signals from the data to focus on the (broad-band) interannual EOFs. We examine ENSO, NAO, North Pacific Oscillation, and other less prominent IAOs that can be identified from our data set, and compute their respective contributions to the variation of global gravitational field and geocenter motion. Finally, we compare the results with the available observational data, and discuss the implications w.r.t. the upcoming space missions such as GRACE.

Chao, Benjamin F.↗
DOE OSTI2020

A spatial-statistical investigation of surface expressions associated with cyclic steaming in the Midway-Sunset Oil Field, California

In the Midway Sunset Oil Field in Central California, operators inject steam into the shallow diatomite formation to enhance heavy oil recovery through imbibition, wettability alteration, and viscosity reduction, among other mechanisms. The injected steam, however, does not always remain in the reservoir or return through the wells. In two zones in the study area, the steam comes out at the surface, creating sinkholes, seeps, and steam outlets. These phenomena, called “surface expressions,” pose safety and environmental hazards. Even though these surface expressions are a widespread problem in Central California, they are not well documented and understood. Possible causes of the surface expressions include: high injection pressure, structurally controlled flow patterns, leakage of steam through old improperly abandoned wells, high injection volumes, or flow along naturally occurring faults, among other possible factors. This work examines attributes of the zones with surface expressions in order to determine factors that may contribute to their occurrence. Spatial statistical analysis using logistic regression, random forests, and classification trees is used to explore the relationship between the surface expressions and geological and production-related attributes. The results point to a significant spatial correlation between the surface expressions and two predictors: concentration of plugged wells and geologic seal thickness. The results guide follow-up studies to further investigate the role of well abandonment and seal thickness in the occurrence of surface expressions.

02 PETROLEUM↗
DOE OSTI2024

Extraction of Vibration Data with Imaging

To date, the primary sensing technology used to measure the vibration response has been accelerometers and strain gages mounted directly to the structure and using either wired or, more recently, wireless telemetry. Cost issues with these sensors and the associated data acquisition systems typically limit the numbers that are deployed on in situ structures. Although there are a few structures with larger sensing counts that in some cases exceed over 1000 sensors, more typical numbers range from ten to one hundred sensors resulting in low spatial resolution when they are applied to physically large systems. When one considers that nuclear power plant structures usually have complex geometries, material properties, connectivity and boundary conditions, it is clear these current approaches to vibration measurements can only provide limited information about a system’s dynamics response characteristics. As an alternative, many non-contact measurement technologies have emerged, including point wise measurement methods such as Global Positioning System (GPS), microwave interferometry, and laser Doppler vibrometry (LDV), as well as simultaneous full-field measurement methods such as electronic speckle pattern interferometry, holography interferometry, and muon tomography, some of which can provide high spatial resolution measurements. Among these methods, digital video imaging techniques have emerged as a feasible solution for full-field vibration measurements that provide significantly more detailed dynamic response information because every pixel becomes a measurement point. Furthermore, recent advances in image processing and computer vision algorithms have been successfully used to process video data for experimental and operational modal analysis. Such full-field measurements have the potential to significantly improve many current structural assessment procedures including system identification (modal parameter estimation), structural health monitoring, load reconstruction, model validation, and model updating. Furthermore, more recent full-field imaging techniques can be accomplished with relatively low-cost, commercially-available off-the-shelf cameras. However, these measurement procedures have other limitations that must be considered such as the ability to only measure visibly accessible points on a structure and a more limited dynamic range and bandwidth than can be achieved with accelerometers or strain gages.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1976

The computer simulation of automobile use patterns for defining battery requirements for electric cars

A Monte Carlo simulation process was used to develop the U.S. daily range requirements for an electric vehicle from probability distributions of trip lengths and frequencies and average annual mileage data. The analysis shows that a car in the U.S. with a practical daily range of 82 miles (132 km) can meet the needs of the owner on 95% of the days of the year, or at all times other than his long vacation trips. Increasing the range of the vehicle beyond this point will not make it more useful to the owner because it will still not provide intercity transportation. A daily range of 82 miles can be provided by an intermediate battery technology level characterized by an energy density of 30 to 50 watt-hours per pound (66 to 110 W-hr/kg). Candidate batteries in this class are nickel-zinc, nickel-iron, and iron-air. The implication of these results for the research goals of far-term battery systems suggests a shift in emphasis toward lower cost and greater life and away from high energy density.

Schwartz, H. J.↗
NASA NTRS2022

Delaware Basin Health and Air Quality: Spatiotemporal Analysis of Air Pollutants Collected from Ground and Space Instruments Around the Guadalupe Mountains and Carlsbad Caverns National Parks

Nitrogen dioxide (NO2) is a precursor for secondary air pollutants, which are associated with decreased visibility, and decreased ecosystem and respiratory health. NO2 is a growing threat to national parks within the Delaware Basin where nearby oil and gas activity contributes to deteriorating park conditions, implying adverse effects on the local tourism economy and public health. To demonstrate spatial and temporal patterns of air pollution in the parks, we examined average monthly, seasonal, and annual tropospheric column concentrations of NO2 in Carlsbad Caverns (CAVE) and Guadalupe Mountain (GUMO) National Parks. We used both the NASA Ozone Monitoring Instrument (OMI) and the European Space Agency (ESA) Tropospheric Monitoring Instrument (TROPOMI) to map NO2 tropospheric column densities. Using ground-based emissions values from the Environmental Protection Agency National Emissions Inventory (NEI) for two nearby natural gas processing plants (Indian Basin Gas Plant, South Carlsbad Plant), we extrapolated monthly trends from these point sources and compared seasonal emissions levels with the measurements recorded by OMI and TROPOMI. The NEI data show an 8% increase in flaring from 2013–2021. OMI measured a 38.3% NO2 increase over the Delaware Basin, 15.29% increase over CAVE, and 4.26% decrease over GUMO from 2011–2018. TROPOMI measured a -1% NO2 change over the Delaware Basin, 3% over CAVE, and 7% over GUMO from 2018–2020. The analysis indicates a positive correlation between emissions from fossil fuel exploration and NO2 concentrations above CAVE and GUMO. This information will inform National Park Service air quality monitoring and policy efforts to ensure compliance with the Clean Air Act.

Sean Cusick↗
NASA NTRS2004

Subsonic Aircraft With Regression and Neural-Network Approximators Designed

At the NASA Glenn Research Center, NASA Langley Research Center's Flight Optimization System (FLOPS) and the design optimization testbed COMETBOARDS with regression and neural-network-analysis approximators have been coupled to obtain a preliminary aircraft design methodology. For a subsonic aircraft, the optimal design, that is the airframe-engine combination, is obtained by the simulation. The aircraft is powered by two high-bypass-ratio engines with a nominal thrust of about 35,000 lbf. It is to carry 150 passengers at a cruise speed of Mach 0.8 over a range of 3000 n mi and to operate on a 6000-ft runway. The aircraft design utilized a neural network and a regression-approximations-based analysis tool, along with a multioptimizer cascade algorithm that uses sequential linear programming, sequential quadratic programming, the method of feasible directions, and then sequential quadratic programming again. Optimal aircraft weight versus the number of design iterations is shown. The central processing unit (CPU) time to solution is given. It is shown that the regression-method-based analyzer exhibited a smoother convergence pattern than the FLOPS code. The optimum weight obtained by the approximation technique and the FLOPS code differed by 1.3 percent. Prediction by the approximation technique exhibited no error for the aircraft wing area and turbine entry temperature, whereas it was within 2 percent for most other parameters. Cascade strategy was required by FLOPS as well as the approximators. The regression method had a tendency to hug the data points, whereas the neural network exhibited a propensity to follow a mean path. The performance of the neural network and regression methods was considered adequate. It was at about the same level for small, standard, and large models with redundancy ratios (defined as the number of input-output pairs to the number of unknown coefficients) of 14, 28, and 57, respectively. In an SGI octane workstation (Silicon Graphics, Inc., Mountainview, CA), the regression training required a fraction of a CPU second, whereas neural network training was between 1 and 9 min, as given. For a single analysis cycle, the 3-sec CPU time required by the FLOPS code was reduced to milliseconds by the approximators. For design calculations, the time with the FLOPS code was 34 min. It was reduced to 2 sec with the regression method and to 4 min by the neural network technique. The performance of the regression and neural network methods was found to be satisfactory for the analysis and design optimization of the subsonic aircraft.

Patnaik, Surya N.↗
DOE OSTI2023

Crowd cluster data in the USA for analysis of human response to COVID-19 events and policies

We provide data on daily social contact intensity of clusters of people at different types of Points of Interest (POI) by zip code in Florida and California. This data is obtained by aggregating fine-scaled details of interactions of people at the spatial resolution of 10 m, which is then normalized as a social contact index. We also provide the distribution of cluster sizes and average time spent in a cluster by POI type. This data will help researchers perform fine-scaled, privacy-preserving analysis of human interaction patterns to understand the drivers of the COVID-19 epidemic spread and mitigation. Current mobility datasets either provide coarse-level metrics of social distancing, such as radius of gyration at the county or province level, or traffic at a finer scale, neither of which is a direct measure of contacts between people. We use anonymized, de-identified, and privacy-enhanced location-based services (LBS) data from opted-in cell phone apps, suitably reweighted to correct for geographic heterogeneities, and identify clusters of people at non-sensitive public areas to estimate fine-scaled contacts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

The chemical composition of globular clusters in the Local Group

We present detailed chemical abundance measurements for 45 globular clusters (GCs) associated with galaxies in (and, in one case, beyond) the Local Group. The measurements are based on new high-resolution integrated-light spectra of GCs in the galaxies NGC 185, NGC 205, M 31, M 33, and NGC 2403, combined with reanalysis of previously published observations of GCs in the Fornax dSph, WLM, NGC 147, NGC 6822, and the Milky Way. The GCs cover the range -2.8 < [Fe/H] < -0.1 and we determined abundances for Fe, Na, Mg, Si, Ca, Sc, Ti, Cr, Mn, Ni, Cu, Zn, Zr, Ba, and Eu. Corrections for non local thermodynamic equilibrium effects are included for Na, Mg, Ca, Ti, Mn, Fe, Ni, and Ba, building on a recently developed procedure. For several of the galaxies, our measurements provide the first quantitative constraints on the detailed composition of their metal-poor stellar populations. Overall, the GCs in different galaxies exhibit remarkably uniform abundance patterns of the α, iron-peak, and neutron-capture elements, with a dispersion of less than 0.1 dex in [α/Fe] for the full sample. There is a hint that GCs in dwarf galaxies are slightly less α-enhanced (by ~0.04 dex on average) than those in larger galaxies. One GC in M 33 (HM33-B) resembles the most metal-rich GCs in the Fornax dSph (Fornax 4) and NGC 6822 (SC7) by having α-element abundances closer to scaled-solar values, possibly hinting at an accretion origin. A principal components analysis shows that the α-element abundances strongly correlate with those of Na, Sc, Ni, and Zn. Several GCs with [Fe/H] < -1.5 are deficient in Mg compared to other α-elements. We find no GCs with strongly enhanced r-process abundances as reported for metal-poor stars in some ultra-faint dwarfs and the Magellanic Clouds. The similarity of the abundance patterns for metal-poor GCs in different environments points to similar early enrichment histories and only allow for minor variations in the initial mass function.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

A comparative study of Conroy and Monte Carlo methods applied to multiple quadratures and multiple scattering

An efficient numerical method of multiple quadratures, the Conroy method, is applied to the problem of computing multiple scattering contributions in the radiative transfer through realistic planetary atmospheres. A brief error analysis of the method is given and comparisons are drawn with the more familiar Monte Carlo method. Both methods are stochastic problem-solving models of a physical or mathematical process and utilize the sampling scheme for points distributed over a definite region. In the Monte Carlo scheme the sample points are distributed randomly over the integration region. In the Conroy method, the sample points are distributed systematically, such that the point distribution forms a unique, closed, symmetrical pattern which effectively fills the region of the multidimensional integration. The methods are illustrated by two simple examples: one, of multidimensional integration involving two independent variables, and the other, of computing the second order scattering contribution to the sky radiance.

Deepak, A.↗
NASA NTRS1985

Vegetation analysis for arid lands geobotany

Three primary study sites were selected for measurement of plant phenological properties and spectral analysis. The sites selected represented typical sagebrush, creosote bush, and saltbush communities in Owens Valley, CA. Community composition was studied at these three sites plus five burned sites. Ten 50 m transects at each locality were measured for percent cover (over 10 cm) by a given species. On each transect two point quarter and five nearest neighbor analyses were conducted. These data provided percent cover, cover by area, plant size, tendency for association, and recolonization patterns after a disturbance. Plots representing percentage plant cover for six sites are included.

Barbour, M. G.↗
NASA NTRSDate not supplied

Regional Variations in the Diurnal Cycle of Tropical Precipitation as Represented by IMERG, ERA5, and Spaceborne Ku Radar

The diurnal cycle of precipitation is highly regional and is typically a product of multiple competing effects that can be highly localized. The diurnal cycle in high precipitation regions such as the Amazon and the Maritime Continent are of particular interest, especially due to the complex coastal effects which take place over the Maritime Continent. The high spatial and temporal resolution provided by the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (GPM) mission (IMERG) dataset, is used in this study to examine the fine-scale features of the diurnal cycle in these regions. Using an 18-year (2000 – 2018) record of IMERG precipitation observations, diurnal and semidiurnal phase and amplitude are calculated using a fast Fourier transform (FFT) method on precipitation averaged for each half-hour of the day at 0.1°x0.1° spatial resolution. We first introduce an objective method of identifying locations where the diurnal signals are robust and strong. Clear patterns of precipitation phase propagation with distance from shore are shown over both regions, with the diurnal phase and amplitude exhibiting a strong dependence on the distance from the coastline. Semidiurnal cycles are generally weaker than the diurnal cycle except in some isolated locations. Similar analysis is also conducted on the ERA5 reanalysis data in order to evaluate the model’s representation of the precipitation diurnal cycle. The model captures the broad scale patterns of diurnal variability but does not capture all the fine scale patterns nor the exact timing that is observed by IMERG. Comparisons are also made to a long record Ku radar dataset created by combining Tropical Rainfall Measuring Mission (TRMM) and GPM observations, thus providing an additional point of comparison for the timing of the ERA5 precipitation peak, since the timing precipitation can be different, even in between observational datasets.

L J Hayden↗
DOE OSTI2026

A Bayesian desmearing algorithm for Bonse–Hart USANS with anisotropic scattering

Ultra-small-angle neutron scattering (USANS) using Bonse–Hart optics provides micrometer-scale structural insights but suffers from severe slit-geometry smearing. While well-established for isotropic systems, quantitative desmearing of anisotropic data remains a challenge because conventional corrections break down for non-radial scattering. In this work, we address this by developing a resolution-aware Bayesian framework that explicitly incorporates anisotropy via an affine deformation to the scattering pattern, guided by the principle of parsimony. This results in orientation-resolved point-spread functions that enable a self-consistent determination of both the resolution and deformation parameters. Using Gaussian process regression with uncertainty quantification and a probabilistic correction for multiple scattering, we demonstrate the framework’s effectiveness through numerical benchmarks and experimental studies of a stretched polymer melt. Our approach enables the seamless integration of SANS and USANS data, facilitating quantitative structural analysis of deformed materials at nanometer to micrometer scales.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electric current paths in a Si:P delta-doped device imaged by nitrogen-vacancy diamond magnetic microscopy

The recently-developed ability to control phosphorous-doping of silicon at an atomic level using scanning tunneling microscopy, a technique known as atomic precision advanced manufacturing (APAM), has allowed us to tailor electronic devices with atomic precision, and thus has emerged as a way to explore new possibilities in Si electronics. In these applications, critical questions include where current flow is actually occurring in or near APAM structures as well as whether leakage currents are present. In general, detection and mapping of current flow in APAM structures are valuable diagnostic tools to obtain reliable devices in digital-enhanced applications. In this report, we used nitrogen-vacancy (NV) centers in diamond for wide-field magnetic imaging (with a few-mm field of view and micron-scale resolution) of magnetic fields from surface currents flowing in an APAM test device made of a P delta-doped layer on a Si substrate, a standard APAM witness material. We integrated a diamond having a surface NV ensemble with the device (patterned in two parallel mm-sized ribbons), then mapped the magnetic field from the DC current injected in the APAM device in a home-built NV wide-field microscope. The 2D magnetic field maps were used to reconstruct the surface current densities, allowing us to obtain information on current paths, device failures such as choke points where current flow is impeded, and current leakages outside the APAM-defined P-doped regions. Analysis on the current density reconstructed map showed a projected sensitivity of ~0.03 A m -1 , corresponding to a smallest-detectable current in the 200 μm wide APAM ribbon of ~6 μA. These results demonstrate the failure analysis capability of NV wide-field magnetometry for APAM materials, opening the possibility to investigate other cutting-edge microelectronic devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1994

Global-scale modes of surface temperature variability on interannual to century timescales

Using 100 years of global temperature anomaly data, we have performed a singluar value decomposition of temperature variations in narrow frequency bands to isolate coherent spatio-temporal modes of global climate variability. Statistical significance is determined from confidence limits obtained by Monte Carlo simulations. Secular variance is dominated by a globally coherent trend; with nearly all grid points warming in phase at varying amplitude. A smaller, but significant, share of the secular variance corresponds to a pattern dominated by warming and subsequent cooling in the high latitude North Atlantic with a roughly centennial timescale. Spatial patterns associated with significant peaks in variance within a broad period range from 2.8 to 5.7 years exhibit characteristic El Nino-Southern Oscillation (ENSO) patterns. A recent transition to a regime of higher ENSO frequency is suggested by our analysis. An interdecadal mode in the 15-to-18 years period and a mode centered at 7-to-8 years period both exhibit predominantly a North Atlantic Oscillation (NAO) temperature pattern. A potentially significant decadal mode centered on 11-to-12 years period also exhibits an NAO temperature pattern and may be modulated by the century-scale North Atlantic variability.

Mann, Michael E.↗
DOE OSTI2026

Revealing the Defect‐Driven Ferroelectric Mechanisms of Aluminum Nitride

Wurtzite III-nitride compounds are CMOS-compatible with widespread industrial interest to exercise ferroelectricity, despite their polar structure being highly resistant to polarization reversal. Here, we induce and tune ferroelectric properties in w-AlN via direct-write ion-beam processing, using nanoscale patterned defect engineering as a post-growth alternative to conventional cation substitution. Nanometric piezoresponse spectroscopy of the focused He + beam patterned defect concentrations in ferroelectric Al 0.92 B 0.08 N measures a localized 10x enhancement in effective piezoresponse and 40% reduction in switching barrier. The irradiation-induced point defects convert piezoelectric AlN into a ferroelectric system with site-saturated nucleation and raise the dielectric susceptibility, switched polarization, and effective piezoelectric coefficient. Enhanced defect-lattice interactions in AlN increase carrier conduction and phonon scattering loss but preserve long-range crystallinity. Here, based on atomistic analysis of nudged elastic band density functional theory calculations and reactive force field simulations, both nitrogen vacancies and defect complexes disrupt bond ordering, facilitating a line-by-line low-barrier switching of pristine AlN.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗