Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ROUGHNESS”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗
DOE OSTI2025

The glass-ceiling convective regime and the origin and diversity of coronae on Venus

Venus and Earth are rocky planets of roughly the same size and bulk density, yet their surface volcanic and tectonic features appear substantially different. On Venus, the coexistence of large volcanic highlands—interpreted as the surface expression of long-lived mantle plumes—alongside coronae, smaller features thought to be caused by transient thermal diapirs, remains enigmatic. Using two-dimensional numerical models of mantle convection with sharp and broad mineral phase transitions for pyrolite, we show that both scales of upwellings can be generated in a stagnant lid planet with an interior temperature 250 to 400 K warmer than Earth’s. The smaller plumes originate from a ~600 km deep internal layer that exists as a consequence of the different sequence of mineral phase transitions that occur in warmer mantles less processed and differentiated by partial melting and volcanism. Future models that include melting will provide further tests of our hypothesis.

Science & Technology - Other Topics↗
DOE OSTI2026

How earthquakes organize stress

Stress is not uniform in the Earth. Therefore, we must use natural experiments to measure the distribution of stresses and related quantities, rather than single values. For instance, dynamic triggering shows that faults are uniformly distributed over their loading cycles in Southern California. The probability that a fault ruptures across a barrier measures the in situ energy distribution. Fault roughness reflects the distribution of strength. These natural experiments produce observable distributions that are surprisingly consistent and suggest some degree of self-organization in the Earth’s crust. Once established, the functional form of the distributions can be used to track changes in response to earthquakes as well as to distinguish fundamentally different fault systems. Transient fault locking before stress release in laboratory experiments can be interpreted as a consequence of self-organization of fault stress. The robust self-organization of multiple variables in earthquake systems suggests that the most consequential mechanical outcome of earthquakes may be the redistribution of stress and the strain energy associated with it. The low friction on a fault during seismic slip as inferred by temperature measurements of the Tohoku earthquake is consistent with dissipation playing a secondary role to this redistribution process. Through stress redistribution and interaction, subduction zone faults tend to synchronize, perhaps due to their geometric simplicity, while the continental system of Southern California cannot synchronize, perhaps due to the complexity of the fault network. Earthquakes organize stress in the crust and produce a suite of well-defined, consistent distributions.

earthquakes↗
DOE OSTI2022

Accurate Compression of Tabulated Chemistry Models with Partition of Unity Networks

Tabulated chemistry models are widely used to simulate large-scale turbulent fires in applications including energy generation and fire safety. Tabulation via piecewise Cartesian interpolation suffers from the curse-of-dimensionality, leading to a prohibitive exponential growth in parameters and memory usage as more dimensions are considered. Artificial neural networks (ANNs) have attracted attention for constructing surrogates for chemistry models due to their ability to perform high-dimensional approximation. However, due to well-known pathologies regarding the realization of suboptimal local minima during training, in practice they do not converge and provide unreliable accuracy. Partition of unity networks (POUnets) are a recently introduced family of ANNs which preserve notions of convergence while performing high-dimensional approximation, discovering a mesh-free partition of space which may be used to perform optimal polynomial approximation. In this work, we assess their performance with respect to accuracy and model complexity in reconstructing unstructured flamelet data representative of nonadiabatic pool fire models. Our results show that POUnets can provide the desirable accuracy of classical spline-based interpolants with the low memory footprint of traditional ANNs while converging faster to significantly lower errors than ANNs. For example, we observe POUnets obtaining target accuracies in two dimensions with 40 to 50 times less memory and roughly double the compression in three dimensions. We also address the practical matter of efficiently training accurate POUnets by studying convergence over key hyperparameters, the impact of partition/basis formulation, and the sensitivity to initialization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The Impulse Approximation Scattering Function and Its Use in Monte Carlo Photon Transport Simulations

The viability of using the impulse approximation scattering function in Monte Carlo photon transport simulations is explored. This scattering function can be constructed from the double differential incoherent scattering cross section developed by Ribberfors and Berggren. A commonly used method for modeling photon Doppler broadening, which is referred to as the hybrid Doppler broadening method, can also be derived from this cross section. A new photon Doppler broadening method, called the consistent Doppler broadening method, is derived and discussed. This method eliminates some of the commonly employed approximations in the hybrid Doppler broadening method, in part, by using the impulse approximation scattering function. Integrated incoherent cross sections generated using the impulse approximation scattering function and the widely used Waller-Hartree scattering function are in good agreement above 20 keV. Below 20 keV, differences as high as 70% are observed, which differs from the roughly 5% differences observed by Ribberfors for some of the materials. Integral and spectral quantities for two problems are also generated using the Monte Carlo photon transport capabilities of the Framework for REsearch in Nuclear ScIence and Engineering (FRENSIE). Due to the small, relative result differences observed when using the impulse approximation scattering function, it is considered a viable alternative to the Waller-Hartree scattering function. In addition, some small, but expected, differences in spectral fluxes at low energies can be avoided by adopting the consistent Doppler broadening method.

42 ENGINEERING↗
DOE OSTI2025

Evaluation of PBR Spent Fuel Criticality and Dose Rate Compliance for Storage and Transportation

Spent tri-structural isotropic (TRISO)–based fuels have a strong track record in storage and transportation without documented incidents. This work seeks to reduce uncertainty to aid in more informed spent fuel management of TRISO-based fuels by modeling both fresh and spent pebble bed reactor (PBR) fuel and comparing the results to the regulatory standards from 10 CFR 71. SCALE was used for all modeling due to it having fast and accurate methods for handling PBR fuel modeling, as well as having an efficient method for shielding calculations in monaco with automated variance reduction using importance calculations (MAVRIC), which utilizes the consistent adjoint-driven importance sampling (CADIS) and the forward-weighted consistent adjoint-driven importance sampling (FW-CADIS) methods. KENO-VI was used for all criticality calculations, TSUNAMI was used for uncertainty quantification on k-effective, TRITON and the Oak Ridge isotope generation code (ORIGEN) were both used for depletion of the fuel, and MAVRIC was used for shielding calculations. For criticality assessments, this study focused on the requirement that the value of the neutron multiplication factor, k-effective (k-eff), would not exceed a peak value of 0.95, including uncertainty, with 95% confidence. Criticality was initially examined by modeling fresh fuel from three different designs—HTR-10 fuel, PBMR-400 fuel, and demonstration fuel representative of a TRISO-fueled modern high-temperature gas reactor (HTGR) design, henceforth referred to as Demo HTGR—and placing them into various sized containers with conditions described in 10 CFR 71 to quantify the peak k-eff state. When the peak value of 0.95 k-eff was exceeded, mitigation methods were examined in those scenarios. Burnup credit, pebble displacement in areas of strong neutron multiplication, and random pebble replacement using pebbles of various compositions and replacement fractions were examined. In summary, the criticality of PBR fuels can be well accounted for by restricting container size, taking credit for burnup, or by displacing/replacing pebbles. Uncertainty of the k-eff due to nuclear data uncertainties was recorded at ~0.6644%Δk/k, or roughly 664% mil (pcm). The nuclear data–induced uncertainty was relatively small and should not require significant modification in the design to be accounted for. Revisions to the evaluated nuclear data file values have been shown to have a larger impact than nuclear data–induced uncertainty. For dose rate aspects, U.S. Nuclear Regulatory Commission regulations require a maximum dose rate of 10 millirem per hour (mrem/h) at 2 meters. In examining the dose rate behavior of spent PBR fuel, the representative Demo HTGR fuel was modeled exclusively due to it possessing the highest target burnup of the examined fuels. Equilibrium cycle modeling methods were used to produce a higher-fidelity discharge isotopic composition than simple assumptions, such as reflected pebbles. The discharge composition was used as a source term in the fixed-source transport shielding calculations, and dose rates were calculated at 2 m for the shortest possible cooling time. The low concentration of fuel material led to dose rates that were in line with regulatory limits, despite the high burnup when compared to traditional light water reactor fuels. In conclusion, the methods employed in this study would require more work to further verify and validate and are limited to the criticality and dose rate analyses performed.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Fast matrix algebra for Bayesian model calibration

In Bayesian model calibration, evaluation of the likelihood function usually involves finding the inverse and determinant of a covariance matrix. When Markov Chain Monte Carlo (MCMC) methods are used to sample from the posterior, hundreds of thousands of likelihood evaluations may be required. In this paper, we demonstrate that the structure of the covariance matrix can be exploited, leading to substantial time savings in practice. Here, we also derive two simple equations for approximating the inverse of the covariance matrix in this setting, which can be computed in near-quadratic time. The practical implications of these strategies are demonstrated using a simple numerical case study and the "quack" R package. For a covariance matrix with 1000 rows, application of these strategies for a million likelihood evaluations leads to a speedup of roughly 4000 compared to the naive implementation

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The FRIB Decay Station: New Horizons with Rare Isotopes

In May 2022, the Facility for Rare Isotope Beams (FRIB), located on the campus of Michigan State University (MSU), began delivering exotic isotopes to an international community of scientists. New discoveries are now being reported from radioactive decay of neutron-rich nuclei near N = 20 and N = 28.FRIB is expected to produce roughly 80% of the unstable or radioactive isotopes predicted to exist up to uranium (Z = 92). The new user facility is supported by the U.S. Department of Energy, and it is operated by MSU. A high-power superconducting linear accelerator, shaped like a paper-clip, drives the production of these rare isotopes by colliding stable nuclei moving at half the speed of light with a rotating, water-cooled graphite tar-get. These collisions cause the primary stable beam to fragment into a wide variety of unstable nuclei, which can be subsequently filtered through a multistage magnetic separator, the Advanced Rare Isotope Separator, and transported to one of several experimental stations. The FRIB Decay Station initiator (FDSi) (see Figure 1) was developed to enable comprehensive radio-active decay studies of the exotic nuclei produced by FRIB and it was used in the first two experiments in 2022. Further, the FDSi is a highly reconfigurable multidetector system with two focal planes (FP1 for discrete spectroscopy and FP2 for total absorption spectroscopy) that can be optimized for the specific science goals of each experiment. It is designed, built, and operated by a community of users with the sup-port of U.S. funding agencies, including the Department of Energy and National Science Foundation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Realistic operation of two residential cordwood-fired outdoor hydronic heater appliances—Part 3: Optical properties of black and brown carbon emissions

Residential biomass combustion is a source of carbonaceous aerosol. Inefficient combustion, particularly of solid fuels produces large quantities of black and brown carbon (BC and BrC). These particle types are important as they have noted effects on climate forcing and human health. One method of measuring these quantities is by measurement of aerosol light-absorption and scattering, which can be performed using an aethalometer and nephelometer, respectively. These instruments are widely deployed in the study of ambient air and are frequently used in air quality modeling and source apportionment studies. In this study, we will describe (1) a method for measuring primary BC and BrC emissions from two residential log-fired wood hydronic heaters and (2) the BC and BrC emission from these devices over a wide range of operating conditions, such as cold-starts, warm-starts, four different levels of output ranging from 15% to 100% maximum rated output, and periods of repeated cycling. The range in flue-gas BC concentrations, measured using an aethalometer at the 880 nanometer (nm) wavelength, were between 5.09 × 10 2 and 2.24 × 10 4 micrograms per cubic meter (µg/m 3 ) while the scattering coefficient of the flue-gas, measured by a nephelometer at 880 nm, ranged between 2.20 × 10 3 and 8.56 × 10 5 inverse megameters (Mm –1 ). The BrC concentrations, measured using the 370 nm wavelength of an aethalometer, were between 9.10 × 10 1 and 3.56 × 10 4 µg/m 3 . The calculated Angstrom Absorption Exponent (AAE) of the flue-gas aerosol ranged between 1.54 and 3.63. Performing a comparison between the measured BC concentration and an external particulate matter (PM) concentration showed that overall BC makes up roughly a quarter of the PM emitted by either of the two appliances. Further for both appliances, the cold-start and the test phase immediately following it had the highest BC and BrC concentrations, the highest measured scattering coefficient, as well as a low AAE.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink.A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m 2 . It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m 2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m 2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

CFD Simulation of Helium Flow Loop Test Section

A helium flow loop is being assembled at Oak Ridge National Laboratory to analyze heat transfer enhancement for systems such as blanket and divertor components. To efficiently identify optimum geometries for heat transfer enhancement in these applications, simulation work is performed to optimize test section designs that are built and tested in the helium flow loop that operates at 4 MPa and a mass flow rate of 100 g/s. Different ribbed geometries that examine rib shape, rib height, rib orientation, rib spacing, and three dimensional orientation are modeled and simulated in STAR-CCM+ to compare their ability to remove heat and mitigate pressure drop. Following the simulations, models are selected and manufactured for the helium flow loop tests. Simulations initially focus on a hydrodynamic study to determine the appropriate mesh and physics models and then add a heat flux to analyze the heat transfer abilities of the models. The simulations are run in steady state and use a Reynolds-averaged Navier-Stokes k-ε turbulence model. The helium is modeled as an ideal gas. The simulation explores models of geometries that enhance the heat transfer and decrease pressure drop with an overall goal of increasing fluid collision with the wall. Enhanced geometries are simulated to select appropriate designs for manufacturing, and preliminary experimental results are used to validate the simulations. Furthermore, the factors that are being analyzed in the comparison between the experimental and the simulated results include matching thermocouple temperatures, pressure drop, roughness, and fluid velocity.

42 ENGINEERING↗
DOE OSTI2023

Wall diffuser velocity effects on American Shad ( Alosa sapidissima ) inside a fishway entrance channel

Attraction water for fishways is typically introduced through a diffuser inside the entrance channel, often through the floor or wall. In the spring of 2019, this laboratory study examined how 151 adult American Shad (Alosa sapidissima) responded to different gross velocities through a wall diffuser inside a full-scale fishway entrance channel. Two velocity conditions were studied, 0.152 m/s and 0.305 m/s, both without turning vanes inside the auxiliary water channel. The fish were tracked using the passive integrated transponder telemetry technique. The results of the experiments showed no difference in American Shad behavior when exhibited to the low and high velocity treatments. Moreover, shad passed the diffuser in roughly 3 out of every 4 attempts, regardless of the treatment. However, the similarity in shad behavior and passage performance is believed to be more of a result of the similarity in flow fields that resulted from the lack of flow guidance devices inside the auxiliary water channel. Furthermore, these findings therefore highlight the importance of properly maintained flow guidance devices, an often-overlooked component of an auxiliary water system.

13 HYDRO ENERGY↗
DOE OSTI2021

EVAP: A two-photon imaging tool to study conformational changes in endogenous Kv2 channels in live tissues

A primary goal of molecular physiology is to understand how conformational changes of proteins affect the function of cells, tissues, and organisms. Here, we describe an imaging method for measuring the conformational changes of the voltage sensors of endogenous ion channel proteins within live tissue, without genetic modification. We synthesized GxTX-594, a variant of the peptidyl tarantula toxin guangxitoxin-1E, conjugated to a fluorophore optimal for two-photon excitation imaging through light-scattering tissue. We term this tool EVAP (Endogenous Voltage-sensor Activity Probe). GxTX-594 targets the voltage sensors of Kv2 proteins, which form potassium channels and plasma membrane–endoplasmic reticulum junctions. GxTX-594 dynamically labels Kv2 proteins on cell surfaces in response to voltage stimulation. To interpret dynamic changes in fluorescence intensity, we developed a statistical thermodynamic model that relates the conformational changes of Kv2 voltage sensors to degree of labeling. We used two-photon excitation imaging of rat brain slices to image Kv2 proteins in neurons. We found puncta of GxTX-594 on hippocampal CA1 neurons that responded to voltage stimulation and retain a voltage response roughly similar to heterologously expressed Kv2.1 protein. Our findings show that EVAP imaging methods enable the identification of conformational changes of endogenous Kv2 voltage sensors in tissue.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Characterizing gravitational wave detector networks: from A # to cosmic explorer

Gravitational-wave observations by the laser interferometer gravitational-wave observatory (LIGO) and Virgo have provided us a new tool to explore the Universe on all scales from nuclear physics to the cosmos and have the massive potential to further impact fundamental physics, astrophysics, and cosmology for decades to come. In this paper we have studied the science capabilities of a network of LIGO detectors when they reach their best possible sensitivity, called A # , given the infrastructure in which they exist and a new generation of observatories that are factor of 10 to 100 times more sensitive (depending on the frequency), in particular a pair of L-shaped cosmic explorer (CE) observatories (one 40 km and one 20 km arm length) in the US and the triangular Einstein telescope with 10 km arms in Europe. We use a set of science metrics derived from the top priorities of several funding agencies to characterize the science capabilities of different networks. The presence of one or two A # observatories in a network containing two or one next generation observatories, respectively, will provide good localization capabilities for facilitating multimessenger astronomy (MMA) and precision measurement of the Hubble parameter. Two CE observatories are indispensable for achieving precise localization of binary neutron star events, facilitating detection of electromagnetic counterparts and transforming MMA. Their combined operation is even more important in the detection and localization of high-redshift sources, such as binary neutron stars, beyond the star-formation peak, and primordial black hole mergers, which may occur roughly 100 million years after the Big Bang. The addition of the Einstein Telescope to a network of two CE observatories is critical for accomplishing all the identified science metrics including the nuclear equation of state, cosmological parameters, the growth of black holes through cosmic history, but also make new discoveries such as the presence of dark matter within or around neutron stars and black holes, continuous gravitational waves from rotating neutron stars, transient signals from supernovae, and the production of stellar-mass black holes in the early Universe. For most metrics the triple network of next generation terrestrial observatories are a factor 100 better than what can be accomplished by a network of three A # observatories.

Einstein telescope↗
DOE OSTI2020

Mechanisms and dynamics of the NH + 2 + H + and NH + + H + + H fragmentation channels upon single-photon double ionization of NH 3

We present state-selective measurements on the NH$_2^{+}$ + H$^{+}$ and NH$^{+}$ + H$^{+}$ + H dissociation channels following single-photon double ionization at 61.5 eV of neutral NH$_{3}$, where the two photoelectrons and two cations are measured in coincidence using 3-D momentum imaging. Three dication electronic states are identified to contribute to the NH$_2^{+}$ + H$^{+}$ dissociation channel, where the excitation in one of the three states undergoes intersystem crossing prior to dissociation, producing a cold NH$_2^+$ fragment. In contrast, the other two states directly dissociate, producing a ro-vibrationally excited NH$_2^+$ fragment with roughly 1 eV of internal energy. The NH$^{+}$ + H$^{+}$ + H channel is fed by direct dissociation from three intermediate dication states, one of which is shared with the NH$_2^{+}$ + H$^{+}$ channel. We find evidence of autoionization contributing to each of the double ionization channels. Finally, the distributions of the relative emission angle between the two photoelectrons, as well as the relative angle between the recoil axis of the molecular breakup and the polarization vector of the ionizing field, are also presented to provide insight on the ionization mechanisms for the different dication states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Superconducting films of MgB 2 via ion beam mixing of Mg/B multilayers

The growth of smooth superconducting MgB 2 films in a non-epitaxial regime is challenging. Here, we study the formation of superconducting MgB 2 films by solid-phase reactive inter-diffusion of sputter-deposited Mg/B multilayers, employing ion beam mixing to disperse the multilayers prior to thermal annealing. The multilayers are intermixed by room-temperature bombardment with 500 keV Xe ions to doses up to 4 × 10 16 cm –2 , followed by thermal annealing to form MgB 2 . Results show that such an intermixing step leads to a dramatic reduction in surface roughness of superconducting films. However, lattice defects produced by ion bombardment reduce the critical superconducting transition temperature, an effect which scales monotonically with ion dose. The critical temperature can be recovered by an additional defect annealing step at $700\ ^{\circ}{\rm C}$.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Single-beam ion source enhanced growth of transparent conductive thin films

A single-beam ion source was developed and used in combination with magnetron sputtering to modulate the film microstructure. The ion source emits a single beam of ions that interact with the deposited film and simultaneously enhances the magnetron discharge. The magnetron voltage can be adjusted over a wide range, from approximately 240 to 130 V, as the voltage of the ion source varies from 0 to 150 V, while the magnetron current increases accordingly. The low-voltage high-current magnetron discharge enables a 'soft sputtering mode', which is beneficial for thin-film growth. Indium tin oxide (ITO) thin films were deposited at room temperature using a combined single-beam ion source and magnetron sputtering. Furthermore, the ion beam resulted in the formation of polycrystalline ITO thin films with significantly reduced resistivity and surface roughness. Single-beam ion-source-enhanced magnetron sputtering has many potential applications in which low-temperature growth of thin films is required, such as coatings for organic solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomic scale etching of diamond: insights from molecular dynamics simulations

Diamond is a promising material for multiple applications in quantum information processing and sensing as well as applications in microelectronics. However, diamond devices can be limited by surface defects that compromise charge stability and spin coherence, among others. Improved strategies in plasma etching of diamond could play an important role in minimizing or eliminating these defects. In this work, we explore plasma-assisted atomic scale etching of diamond using argon ions (Ar + ), hydrogen ions (H + ) and hydrogen atoms (H). We employ classical molecular dynamics (MD) simulations and test several interatomic potentials based on the Reactive Empirical Bond Order (REBO) form with comparisons to a variety of published experimental results. We performed MD simulations of low-energy hydrogen ($\leqslant$50 eV) and argon ( $\leqslant$200 eV) ion bombardment of diamond surfaces. Ar + bombardment can be used to locally smooth initially rough diamond surfaces via the formation of an amorphous C layer, the thickness of which increases with argon ion energy. Subsequent exposure with hydrogen ions (or fast neutrals) will selectively etch this amorphous C layer, leaving the underlying diamond layer mostly intact if the H energy is maintained below about 10 eV. The simulations suggest that combining Ar + smoothing with selective, near threshold energy H removal of amorphous C can be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices.

74 ATOMIC AND MOLECULAR PHYSICS↗