Search NASA⌕ Search

SEARCH · Search NASA

Results for “BCA”

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 19 records

Soot Aerosol In The Atmosphere: Pole-to-Pole Distribution And Contributions by Aircraft

Interest in the distribution of black carbon (soot) aerosol (BCA) in the atmosphere is warranted for the following reasons: (1) BCA has the highest absorption cross section of any compound known, thus it can absorb solar radiation to cause atmospheric warming; (2) BCA is a strong adsorber of gases, thus it can catalyze heterogeneous chemical reactions to modify the chemical composition of the atmosphere; (3) If aircraft emission is the major source of atmospheric BCA, it can serve as an atmospheric tracer of aircraft exhaust. We collect BCA particles greater than or equal to 0.02 micrometer diameter by wires mounted on both the DC-8 and ER-2 aircraft. After return to the laboratory, the wires are examined with a field emission scanning electron microscope to identify BCA particles by their characteristic morphology. Typically, BCA exists in the atmosphere as small particles of complex morphology. The particle sizes at the source are measured in tens of Angstrom units; after a short residence time in the atmosphere, individual particles coalesce to loosely packed agglomerates of typical dimensions 0.01 to 0.1 micrometer. We approximate the size of each BCA aggregate by that of a sphere of equivalent volume. This is done by computing the volume of a sphere whose diameter is the mean between averaged minimum and maximum dimensions of the BCA particle. While this procedure probably underestimates the actual surface area, it permits us to compare BCA size distributions among themselves and with other types of aerosols. When statistically justified, we fit lognormal distributions to the data points to determine number concentrations, geometric mean radii, standard deviations, BCA surface areas and volumes. Results to date permit the following conclusions: (1) BCA concentration in the northern stratosphere averages 0.6 ng per cubic meters. This amount is one part in 10(exp 4) after a volcanic eruption (e.g., Pinatubo) increasing to about one percent during volcanic quiescence. In the northern troposphere, BCA concentration averages 3.2 ng per cubic meters, or 0.3 percent of the background aerosol. (2) Applying an BCA emission index EI(BCA)=5 x 10(exp -5), measured in the exhaust wake of a Concorde supersonic jet aircraft, to realistic estimates of fuel burnt by the current and projected fleets permits us to conclude that: (i) Most BCA in the northern stratosphere results from aircraft emissions; (ii) Most BCA in the northern troposphere results from other sources than aircraft; (iii) A projected supersonic fleet will increase the northern stratospheric BCA concentration by one order of magnitude, unless the emission index is substantially reduced. (3) A strong gradient between the northern and southern hemispheres indicates that mixing across the equator is greatly inhibited in relation to atmospheric residence times of BCA. (4) The single scatter albedo of BCA/"background" aerosol mixtures suggests a cooling effect for most of the globe; an exemption is the Arctic because of the high surface albedo of the snow/ice covered earth's surface.

Pueschel, R. F.↗

Stratospheric Sulfuric Acid and Black Carbon Aerosol Measured During POLARIS and its Role in Ozone Chemistry

Stratospheric aerosol can affect the environment in three ways. Sulfuric acid aerosol have been shown to act as sites for the reduction of reactive nitrogen and chlorine and as condensation sites to form Polar Stratospheric Clouds, under very cold conditions, which facilitate ozone depletion. Recently, modeling studies have suggested a link between BCA (Black Carbon Aerosol) and ozone chemistry. These studies suggest that HNO3, NO2, and O3 may be reduced heterogeneously on BCA particles. The ozone reaction converts ozone to oxygen molecules, while HNO3 and NO2 react to form NOx. Finally, a buildup of BCA could reduce the single-scatter albedo of aerosol below a value of 0.98, a critical value that has been postulated to change the effect of stratospheric aerosol from cooling to warming. Correlations between measured BCA amounts and aircraft usage have been reported. Attempts to link BCA to ozone chemistry and other stratospheric processes have been hindered by questions concerning the amount of BCA that exists in the stratosphere, the magnitude of reaction probabilities, and the scarcity of BCA measurements. The Ames Wire Impactors (AWI) participated in POLARIS as part of the complement of experiments on the NASA ER-2. One of our main objectives was to determine the amount of aerosol surface area, particularly BCA, available for reaction with stratospheric constituents and assess if possible, the importance of these reactions. The AWI collects aerosol and BCA particles on thin Palladium wires that are exposed to the ambient air in a controlled manner. The samples are returned to the laboratory for subsequent analysis. The product of the AWI analysis is the size, surface area, and volume distributions, morphology and elemental composition of aerosol and BCA. This paper presents results from our experiments during POLARIS and puts these measurements in the context of POLARIS and other missions in which we have participated. It describes modifications to the AWI data analysis procedures in which the collection of BCA is modeled as a fractal aggregate. The new analysis results in an increase in BCA surface area of approximately 24 and an increase in mass of 7-10 from the previous method. For the current study, BCA surface area is used in computer models that attempt to predict measured NOx/NOy ratios and O3 depletion rates. Inclusion of the HNO3 reaction with BCA in one model tends to improve the agreement of calculated to measured NOx/NOy ratio. However, it was found that these trends are viable only if the reactions are catalytic.

Strawa, Anthony W.↗

Pole-to-Pole Distribution of Stratospheric Black Carbon (Soot) Aerosol from Aircraft

The distribution of black carbon (soot) aerosol (BCA) in the atmosphere is of interest for several reasons: (1) Because BCA has the highest absorption cross section of any compound known, it can absorb solar radiation to cause atmospheric warming. (2) Because it is a strong adsorber of gases, it can catalyze heterogeneous reactions to change the chemical composition of the atmosphere.(3) If aircraft are a major source of BCA, it is an important tracer of aircraft emissions. Analysis for BCA of impactor samples from Arctic and Antarctic deployments, utilizing particle morphology of scanning electron microscopy images, permits the following conclusions: (1) The BCA concentration in the northern stratosphere varies between 0 and 2.6 ng m-3 averaging 0.6 ng/cu m. (2) This BCA loading is commensurate with estimated fuel consumptions in the stratosphere by the current commercial fleet and an emission index E=0.03 g BCA per kg fuel burnt which was measured in jet exhaust at al titude.Thus, most stratospheric BCA in the northern stratosphere results from aircraft emissions. The background BCA concentration in the southern stratosphere varies between 0 and 0.6 ng cu m averaging 0.1 ng/cu m. This strong meridional gradient implies that stratospheric BCA residence time- is shorter than are mixing times between hemispheres. Projected annual fuel consumption of a future supersonic commercial fleet is 7E13 g. This fleet would increase stratospheric BCA loadings by a factor of 2-3, because almost all fuel would be burnt above the tropopause. An improved EI(BCA) by a factor of ten would result in an increase of stratospheric BCA loadings by approximately 50 %.

Pueschel, R. F.↗

Soot Aerosols in the Atmosphere: Contributions by Aircraft

Interest in the distribution of black carbon (soot) aerosol (BCA) in the atmosphere is based on the following: (1) Because BCA has the highest absorption cross section of any compound know, it can absorb solar radiation to cause atmospheric warming; (2) Because BCA is a strong adsorber of gases, it can catalyze heterogeneous chemical reactions to modify the chemical composition of the atmosphere; (3) If aircraft emission is the major source of BCA, it can serve as an atmospheric tracer of aircraft exhaust. We collect BCA particles as small as 0.02 micrometers by wires mounted on both the DC-8 and ER-2 aircraft. After return to the laboratory, the wires are examined with a field emission scanning electron microscope to identify BCA particles by their characteristics morphology, Typically, BCA exists in the atmosphere as small particles of complex morphology. The particle sizes at the source are measured in tens of Angstrom units; after a short residence time in the atmosphere, individual particles coalesce to loosely packed agglomerates of typical dimensions 0.01 to 0.1 micrometer. We approximate the size of each BCA aggregate by that of a sphere of equivalent volume. This is done by computing the volume of a sphere whose diameter is the mean between averaged minimum and maximum dimensions of the BCA particle. While this procedure probably underestimates the actual surface area, it permits us to compare BCA size distributions among themselves and with other types of aerosols.

Pueschel, R. F.↗

Black carbon (soot) aerosol in the lower stratosphere and upper troposphere

As determined by impactor samplers flown on ER-2 and DC-8 aircraft, black carbon aerosol (BCA) mass loadings in the stratosphere average 0.6 nanograms per standard cubic meter, or 0.01 percent of the total aerosol. Upper tropospheric BCA increases to 0.3 percent. Low stratospheric concentration is commensurate with present commercial air traffic fuel consumption, given the following assumptions: the BCA emissions are 0.1 grams per kilogram of fuel consumed, 10 percent of route mileage is above the tropopause, and average BCA stratospheric residence time is about one year. Taking BCA into account, the stratospheric single scatter albedo is about 0.99. Using parameters for planned supersonic commercial aircraft, whose emissions will be predominantly in the stratosphere, it is shown that such traffic will double stratospheric BCA concentration. This would reduce the aerosol single scattering albedo by one percent, and double the BCA surface area that is available for heterogeneous chemistry.

Pueschel, R. F.↗

Black Carbon (Soot) Aerosol in the Lower Stratosphere and Upper Troposphere

We have used two approaches to measure BCA in the stratosphere. The first method consists of collecting aerosols by impaction on quartz rods, upon which are mounted strips of polystyrene filter material. Each 25 mm by 0.5 mm filter strip is alternately masked and exposed in 5 mm long sections. After sampling, the optical density of the exposed sections is measured and compared with the optical density of the masked (unexposed) sections. The filter strip is a Lambertian collector, the scattering property of which is not affected by the deposition of liquid sulfuric acid droplets which dominate the stratospheric aerosol. Hence, all attenuation greater in the exposed than in the non-exposed sections is due to absorption by BCA or rare opaque materials such as meteoritic iron. In a second approach we expose carbon-coated gold or palladium wires to air outside the aircraft boundary layer [Farlow et al., 1979; Pueschel et al., 1989]. With a field emission scanning electron microscope we are able to characterize BCA particles in addition to the dominant sulfuric acid droplets. Typically, soot is comprised of 20 nm diameter spheres linked together as branching chains or loosely packed aggregates (Figure 1a). The microscope has a lateral image resolution of 1.5 nm and can clearly resolve individual spheres making up BCA aggregates. Their morphology is uniquely different from that of the liquid sulfuric acid-water aerosol (Figure 1b). We approximate the size of each BCA particle by that of a sphere of equivalent volume and fit a lognormal size distribution to both the BCA and H2SO4/H2O aerosols. Both sampler types have been mounted on NASA's ER-2 high-altitude research aircraft (altitude ceiling of 70,000 feet), and on NASA's DC-8 (40,000 feet ceiling) aircraft to sample aerosols from a significant portion of the atmosphere.

Pueschel, R. F.↗

Physical and Optical/Radiative Properties of Arctic Aerosols: Potential Effects on Arctic Climate

We have determined the abundance of light-scattering sulfuric acid (H2SO4/H2O) and light-absorbing black carbon aerosol (BCA) in Spring 1992 in the Arctic atmosphere by airborne in situ sampling with impactors, and measured particle sizes and morphologies by scanning electron microscopy. The mass of BCA in the Arctic troposphere is one percent of the total aerosol, reduced to one part in 104 in the stratosphere. A Mie algorithm permits the calculation of the optical properties of the various aerosol components, and an algorithm developed by Ackerman and Toon and modified to serve our needs lets us calculate the optical effects of the black carbon aerosol that is mixed internally with the sulfuric acid aerosol. It follows that the effect of internally-mixed BCA on the aerosol scattering and absorption properties depends on its location within the droplet. BCA concentrated near the droplet surface has a greater effect on absorption of solar radiation than does the same amount of BCA located near its center. Single scatter albedos of the combined system are omega(sub 0)=1.0 in the post-Pinatubo Arctic stratosphere, and as low as 0.94 in the troposphere. The aerosol has the potential to regionally warm the Arctic earth-atmosphere system, because of the high surface albedo of the snow-covered Arctic.

Pueschel, R. F.↗

Unraveling Vulnerabilities in Endocrine Therapy-Resistant HER2+/ER+ Breast Cancer

Abstract Breast tumors overexpressing human epidermal growth factor receptor (HER2) confer intrinsic resistance to endocrine therapy (ET), and patients with HER2/estrogen receptor–positive (HER2+/ER+) breast cancer (BCa) are less responsive to ET than HER2–/ER+. However, real-world evidence reveals that a large subset of patients with HER2+/ER+ receive ET as monotherapy, positioning this treatment pattern as a clinical challenge. In the present study, we developed and characterized 2 in vitro models of ET-resistant (ETR) HER2+/ER+ BCa to identify possible therapeutic vulnerabilities. To mimic ETR to aromatase inhibitors (AIs), we developed 2 long-term estrogen deprivation (LTED) cell lines from BT-474 (BT474) and MDA-MB-361 (MM361). Growth assays, PAM50 subtyping, and genomic and transcriptomic analyses, followed by validation and functional studies, were used to identify targetable differences between ET-responsive parental and ETR-LTED HER2+/ER+ cells. Compared to their parental cells, MM361 LTEDs grew faster, lost ER, and increased HER2 expression, whereas BT474 LTEDs grew slower and maintained ER and HER2 expression. Both LTED variants had reduced responsiveness to fulvestrant. Whole-genome sequencing of aggressive MM361 LTEDs identified mutations in genes encoding transcription factors and chromatin modifiers. Single-cell RNA sequencing demonstrated a shift towards non-luminal phenotypes, and revealed metabolic remodeling of MM361 LTEDs, with upregulated lipid metabolism and ferroptosis-associated antioxidant genes, including GPX4. Combining a GPX4 inhibitor with anti-HER2 agents induced significant cell death in both MM361 and BT474 LTEDs. The BT474 and MM361 AI-resistant models capture distinct phenotypes of HER2+/ER+ BCa and identify altered lipid metabolism and ferroptosis remodeling as vulnerabilities of this type of ETR BCa.

60 APPLIED LIFE SCIENCES↗

Investigation of Seal-to-Floor Effects on Semi-Span Transonic Models

In an effort to achieve the maximum possible Reynolds number (Re) when conducting production testing for flight loads aerodynamic databases, it has been the preferred practice of The Boeing Company / Commercial Airplanes (BCA) -- Loads and Dynamics Group since the early 1990's to test large scale semi-span models in the 11- By 11-Foot Transonic Wind Tunnel (TWT) leg of the Unitary Plan Wind Tunnel (UPWT) at the NASA Ames Research Center (ARC). There are many problems related to testing large scale semi-span models of high aspect ratio flexible transport wings, such as; floor boundary layer effects, wing spanwise wall effects, solid blockage buoyancy effects, floor mechanical interference effects, airflow under the model effects, or tunnel flow gradient effects. For most of these issues, BCA has developed and implemented either standard testing methods or numerical correction schemes and these will not be discussed in this document. Other researchers have reported on semi-span transonic testing correction issues, however most of the reported research has been for low Mach testing. Some of the reports for low Mach testing address the difficult problem of preventing undesirable airflow under a semi-span model while ensuring unrestricted main balance functionality, however, for transonic models this issue has gone unresolved. BCA has been cognizant for sometime that there are marked differences in wing pressure distributions from semi-span transonic model testing than from full model or flight testing. It has been suspected that these differences are at least in part due to airflow under the model. Previous efforts by BCA to address this issue have proven to be ineffective or inconclusive and in one situation resulted in broken hardware. This paper reports on a Boeing-NASA collaborative investigation based on a series of small tests conducted between June 2006 and November 2007 in the 11 by 11 foot Transonic Wind Tunnel at NASA Ames on three large commercial jet transport configurations to assess the effects of sealing a semi-span model to the floor and to investigate efficient sealing and testing techniques. This document will show how sealing the model to the floor has a small but remarkably far reaching spanwise effect on wing pressures, wing local section forces and wing force summations.

Sleppy, Mark A.↗

Mapping the Binary Covalent Alloy Space to Pursue Superior Nitrogen Reduction Reaction Catalysts

The electrochemical nitrogen reduction reaction (NRR) has the potential to decarbonize industrial ammonia production. However, NRR has poor activity and selectivity versus the competing hydrogen evolution reaction for catalysts that adhere to scaling relations. Overcoming the limitations imposed by scaling relations requires more complex catalyst materials, however, evaluating materials beyond simple metal systems is a large combinatorial problem that requires an improved understanding of the electrocatalyst surface to rationally guide the discovery of superior catalysts. The study uses grand canonical density functional theory to uncover NRR trends on a large and disparate set of binary covalent alloys (BCA) with variable compositions and active-site geometries. The studied BCAs generally follow scaling relations, albeit with larger variance and several systems that significantly break scaling. BCAs with early- to mid-transition metals tend to lie near the volcano peak and activate the N 2 triple bond via a side-on binding configuration. Trends in the BCA space cannot be readily predicted using simple electronic descriptors, which is ascribed to the large geometric variability of the BCA surfaces. Furthermore, it is anticipated that these findings will provide a foundation for the rational design of superior NRR electrocatalysts with increasing material complexity.

36 MATERIALS SCIENCE↗

Molecular dynamics simulations of reflection and sputtering behavior of boron under deuterium ion irradiation

Boronization is a commonly used method of wall conditioning in fusion reactors. The application of boron films to the plasma-facing materials results in enhanced plasma performance due to the reduction of intrinsic impurities. This is primarily driven by a reduction in oxygen content that is chemically trapped in the boron film. The reactive nature of these boron films also raises questions concerning interactions with hydrogen isotopes. In this work, boron-deuterium interactions were studied using molecular dynamics (MD). Reactive force field potentials were used to model the chemical interactions between B and D. An amorphous boron substrate was irradiated by D atoms at varying incident energies, 10 eV < E i < 150 eV and angles, 0° < α < 85°. The reflection probability was calculated and compared to results from the commonly-used binary collision approximation (BCA) method. This comparison found that the BCA underestimated the reflection probability at E i < 35 eV and α > 45 ∘ . The source of this discrepancy was found to be the surface binding energy model. The BCA calculation with an isotropic surface binding energy model was more closely aligned to the MD result. This, in combination with a correction function based on the MD results allows for corrections to the reflection probability of deuterium impinging on boron surfaces. The sputtering of the substrate material was also studied. While this study did not contain sufficient events to quantitatively describe the sputtering behavior, some qualitative results emerged: namely, chemical sputtering of B and D-containing molecules (BD, BD 2 , BD 3 ) at low ( < 20 eV) incident deuterium energies. This result suggests that chemical sputtering could be a significant factor in limiting boron coating lifetime when exposed to lower ion energies, such as those in detached plasmas. The results show that chemical interactions should be taken into account when modeling the interactions between D ions and B surfaces.

Boron↗

Crystallographic and temperature effects in low-energy collisions for plasma–material interactions

The interaction of plasma with materials is critical for the fundamental understanding of non-equilibrium processes and their wide application. Recent achievements in fusion energy research emphasize the importance of this problem. Because modelling and predicting plasma–material interactions (PMI) require considering a tremendous number of single plasma particle–surface interaction events, this type of modelling was, until recently, possible only within a binary collision approximation (BCA). The BCA approach considers materials as uniform temperature-insensitive media. The research presented here addresses the PMI problem within an atomistic-based approach for the first time at a statistically significant level. Approximately 10 5 molecular dynamics trajectories were generated for 100 eV deuterium ions interacting with a tungsten surface. The research demonstrated the critical importance of incorporating the discrete lattice structure of matter into the model. Deuterium penetration depth and fraction of backscattered deuterium ions strongly depend on the surface orientation, the impact incident ion directions, and the ions’ initial positions. On average, the BCA-based calculations underestimate the penetration depth by a factor of two or more, and the fraction of backscattered atoms is overestimated by a factor of two or more. As a result, increasing the material temperature from 500 to 3000 K reduces penetration depth by about 30% and increases the fraction of backscattered atoms by 250%.

36 MATERIALS SCIENCE↗

Single-Particle Measurements of Midlatitude Black Carbon and Light-Scattering Aerosols from the Boundary Layer to the Lower Stratosphere

A single-particle soot photometer (SP2) was flown on a NASA WB-57F high-altitude research aircraft in November 2004 from Houston, Texas. The SP2 uses laser-induced incandescence to detect individual black carbon (BC) particles in an air sample in the mass range of approx.3-300 fg (approx.0.15-0.7 microns volume equivalent diameter). Scattered light is used to size the remaining non-BC aerosols in the range of approx.0.17-0.7 microns diameter. We present profiles of both aerosol types from the boundary layer to the lower stratosphere from two midlatitude flights. Results for total aerosol amounts in the size range detected by the SP2 are in good agreement with typical particle spectrometer measurements in the same region. All ambient incandescing particles were identified as BC because their incandescence properties matched those of laboratory-generated BC aerosol. Approximately 40% of these BC particles showed evidence of internal mixing (e.g., coating). Throughout profiles between 5 and 18.7 km, BC particles were less than a few percent of total aerosol number, and black carbon aerosol (BCA) mass mixing ratio showed a constant gradient with altitude above 5 km. SP2 data was compared to results from the ECHAM4/MADE and LmDzT-INCA global aerosol models. The comparison will help resolve the important systematic differences in model aerosol processes that determine BCA loadings. Further intercomparisons of models and measurements as presented here will improve the accuracy of the radiative forcing contribution from BCA.

Schwartz, J. P.↗

Solar Wind Ion Sputtering of Sodium from Silicates Using Molecular Dynamics Calculations of Surface Binding Energies

For nearly 40 yr, studies of exosphere formation on airless bodies have been hindered by uncertainties in our understanding of the underlying ion collisional sputtering by the solar wind (SW). These ion impacts on airless bodies play an important role in altering their surface properties and surrounding environment. Much of the collisional sputtering data needed for exosphere studies come from binary collision approximation (BCA) sputtering models. These depend on the surface binding energy (SBE) for the atoms sputtered from the impacted material. However, the SBE is not reliably known for many materials important for planetary science, such as plagioclase feldspars and sodium pyroxenes. BCA models typically approximate the SBE using the cohesive energy for a monoelemental solid. We use molecular dynamics (MD) to provide the first accurate SBE data we are aware of for Na sputtered from the above silicate minerals, which are expected to be important for exospheric formation at Mercury and the Moon. The MD SBE values are ∼8 times larger than the Na monoelemental cohesive energy. This has a significant effect on the predicted SW ion sputtering yield and energy distribution of Na and the formation of the corresponding Na exosphere. We also find that the SBE is correlated with the coordination number of the Na atoms within the substrate and with the cohesive energy of the Na-bearing silicate. Our MD SBE results will enable more accurate BCA predictions for the SW ion sputtering contribution to the Na exosphere of Mercury and the Moon.

Lunar surface↗

Effect of cellulose reducing ends and primary hydroxyl groups modifications on cellulose-cellulase interactions and cellulose hydrolysis

Cellulose reducing ends are believed to play a vital role in the cellulose recalcitrance to enzymatic conversion. However, their role in insoluble cellulose accessibility and hydrolysis is not clear. Thus, in this study, reducing ends of insoluble cellulose derived from various sources were modified by applying reducing and/or oxidizing agents. Here, the effects of cellulose reducing ends modification on cellulose reducing ends, cellulose structure, and cellulose accessibility to cellulase were evaluated along with the impact on cellulose hydrolysis with complete as well purified cellulase components. Sodium borohydride (NaBH 4 ) reduction and sodium chlorite-acetic acid (SC/AA) oxidation were able to modify more than 90% and 60% of the reducing ends, respectively, while the bicinchoninic acid (BCA) reagent applied for various cycles oxidized cellulose reducing ends to various extents. X-ray diffractograms of the treated solids showed that these treatments did not change the cellulose crystalline structure and the change in crystallinity index was insignificant. Surprisingly, it was found that the cellulose reducing ends modification, either through selective NaBH 4 reduction or BCA oxidation, had a negligible impact on cellulose accessibility as well on cellulose hydrolysis rates or final conversions with complete cellulase at loadings as low as 0.5 mg protein/g cellulose. In fact, in contrast to what is traditionally believed, modifications of cellulose reducing ends by these two methods had no apparent impact on cellulose conversion with purified cellulase components and their synergy. However, SC/AA oxidation resulted in significant drop in cellulose conversion (10%–50%) with complete as well purified cellulase components. Nonetheless, further research revealed that the cause for drop in cellulose conversion for the SC/AA oxidation case was due to primary hydroxyl groups (PHGs) oxidation and not the oxidation of reducing ends. Furthermore, it was found that the PHGs modification affects cellulose accessibility and slows the cellulase uptake as well resulting in significant drop in cellulose conversions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗

Sputtering of Surfaces by Ion Irradiation: A Comparison of Molecular Dynamics and Binary Collision Approximation Models to Laboratory Measurements

We compare various sputtering simulation methods to experimental results in both the low energy (<1 keV) and high energy (≥1 keV) impact regimes for argon ions impacting a pure copper substrate at normal incidence. Our results indicate that for high energy impacts, both binary collision approximation (BCA) and molecular dynamics methods can be used to generate reasonable predictions for the yield and energy distribution of the sputtered atoms. We also find reasonable agreement between the theoretical and experimental results down to impact energies of 600 eV. However, at 200 eV impact energies, significant discrepancies appear between the experimental and theoretical ejecta energy distributions in the peak position, the width of the energy distribution, and the magnitude of the high energy tail. These discrepancies appear to arise from the experimental results being only for atoms sputtered normal to the substrate surface, whereas the theoretical results are integrated over all 2π solid angles above the surface. Using the BCA code SDTrimSP and limiting the results to only atoms emitted within ±15° of the surface normal brings theory and experiment into reasonable agreement. These results suggest that for low energy impacts, the energy distribution of sputtered atoms is highly dependent on the emission angle of the ejecta.

Solar wind↗