Search NASA⌕ Search

Engineering topics

Evans, James W.

Publications and source records attributed to Evans, James W..

At least 37 records · Page 2

Reaction processes at step edges on S-decorated Cu(111) and Ag(111) surfaces: MD analysis utilizing machine learning derived potentials

We report a variety of complexation, reconstruction, and sulfide formation processes can occur at step edges on the {111} surfaces of coinage metals (M) in the presence of adsorbed S under ultra-high vacuum conditions. Given the cooperative many-atom nature of these reaction processes, Molecular Dynamics (MD) simulation of the associated dynamics is instructive. However, only quite restricted Density Functional Theory (DFT)-level ab initio MD is viable. Thus, for M = Ag and Cu, we instead utilize the DeePMD framework to develop machine-learning derived potentials, retaining near-DFT accuracy for the M–S systems, which should have broad applicability. These potentials are validated by comparison with DFT predictions for various key quantities related to the energetics of S on M(111) surfaces. The potentials are then utilized to perform extensive MD simulations elucidating the above diverse restructuring and reaction processes at step edges. Key observations from MD simulations include the formation of small metal–sulfur complexes, especially MS 2 ; development of a local reconstruction at A-steps featuring an S-decorated {100} motif; and 3D sulfide formation. Additional analysis yields further information on the kinetics for metal–sulfur complex formation, where these complexes can strongly enhance surface mass transport, and on the propensity for sulfide formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular interactions in diffusion-controlled aldol condensation with mesoporous silica nanoparticles

The aldol reaction of p-nitrobenzaldehyde in amino-catalyzed mesoporous silica nanoparticles (MSN) has revealed varying catalytic activity with the size of the pores of MSN. The pore size dependence related to the reactivity indicates that the diffusion process is important. A detailed molecular-level analysis for understanding diffusion requires assessment of the noncovalent interactions of the molecular species involved in the aldol reaction with each other, with the solvent, and with key functional groups on the pore surface. Such an analysis is presented here based upon the effective fragment potential (EFP). The EFP method can calculate the intermolecular interactions, decomposed into Coulomb, polarization, dispersion, exchange-repulsion, and charge-transfer interactions. In this study, the potential energy surfaces corresponding to each intermolecular interaction are analyzed for homo- and hetero-dimers with various configurations. The monomers that compose dimers are five molecules such as p-nitrobenzaldehyde, acetone, n-hexane, propylamine, and silanol. The results illustrate that the dispersion interaction is crucial in most dimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rotational and translational diffusion of liquid n-hexane: EFP-based molecular dynamics analysis

Molecular Dynamics (MD) simulations based on the Effective Fragment Potential (EFP) method are utilized to provide a comprehensive assessment of diffusion in liquid n-hexane. We decompose translational diffusion into components along and orthogonal to the long axis of the molecule. Rotational diffusion is decomposed into tumbling and spinning motions about this axis. Our analysis yields four corresponding diffusion coefficients which are related to diagonal entries in the complete 6 × 6 diffusion tensor accounting for the three rotational and three translational degrees of freedom and for the potential coupling between them. However, coupling between different degrees of freedom is expected to be minimal for a natural choice of the molecular body-fixed axis, so then off-diagonal entries in the tensor are negligible. This expectation is supported by a hydrodynamic analysis of the diffusion tensor which treats the liquid surrounding the molecule being tracked as a viscous continuum. Thus, the EFP MD analysis provides a comprehensive characterization of diffusion and also reveals expected shortcomings of the hydrodynamic treatment, particularly for rotational diffusion, when applied to neat liquids.

74 ATOMIC AND MOLECULAR PHYSICS↗

Preface for the special collection commemorating the career of Pat Thiel

This Special Collection in the Journal of Vacuum Science and Technology A (JVST A) has been compiled to commemorate the career of Professor Patricia (Pat) Thiel as a scientific leader in the international surface science community. Pat passed away in September 2020 due to breast cancer. This Collection includes 26 articles from colleagues, collaborators, and former Thiel group members. The topics overlap Pat’s scientific interests and advances in areas that include chemisorption, reconstruction, and complex formation on metal surfaces; submonolayer nucleation and growth and kinetic roughening of epitaxial thin films; coarsening of surface nanostructures; intercalation; surface phenomena for intermetallics including quasicrystals; and atomic-scale friction. We wish to thank the authors for their contributions towards producing this impressive Special Collection.

99 GENERAL AND MISCELLANEOUS↗

Sulfur-enhanced dynamics of coinage metal(111) surfaces: Step edges versus terraces as locations for metal-sulfur complex formation

The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.

36 MATERIALS SCIENCE↗

Thermodynamics and kinetics of H adsorption and intercalation for graphene on 6 H -SiC(0001) from first-principles calculations

Previous experimental observations for H intercalation under graphene on SiC surfaces motivate the clarification of configuration stabilities and kinetic processes related to intercalation. From first-principles density-functional-theory calculations, we analyze H adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that in the low-H coverage regime (described by a single H atom within a sufficiently large supercell), intercalation into the gallery between TLG and BLG or into the gallery underneath BLG is more favorable thermodynamically than adsorption on top of TLG. However, intercalation into the gallery between TLG and BLG is most favorable. We obtain energy barriers of about 1.3 and 2.3 eV for a H atom diffusing on and under TLG, respectively. From an additional analysis of the energy landscape in the vicinity of a step on the TLG, we assess how readily one guest H atom on the TLG terrace can directly penetrate the TLG into the gallery between TLG and BLG versus crossing a TLG step to access the gallery. Finally, we also perform density functional theory calculations for higher H coverages revealing a shift in favorability to intercalation of H underneath BLG and characterizing the variation with H coverage in interlayer spacings.

36 MATERIALS SCIENCE↗

Phase transitions in Schloegl's second model for autocatalysis on a Bethe lattice

Schloegl's second model (also known as the quadratic contact process) on a lattice involves spontaneous particle annihilation at rate p and autocatalytic particle creation at empty sites with n ≥ 2 occupied neighbors. The particle creation rate for exactly n occupied neighbors is selected here as n(n - 1)/[z(z - 1)] for lattice coordination number z. We analyze this model on a Bethe lattice. Precise behavior for stochastic models on regular periodic infinite lattices is usually surmised from kinetic Monte Carlo simulation on a finite lattice with periodic boundary conditions. However, the persistence of boundary effects for a Bethe lattice complicates this process, e.g., by inducing spatially heterogenous states. This motivates the exploration of various boundary conditions and unconventional simulation ensembles on the Bethe lattice to predict behavior for infinite size. Here, we focus on z = 3, and predict a discontinuous transition to the vacuum state on the infinite lattice when p exceeds a threshold value of around 0.053.

97 MATHEMATICS AND COMPUTING↗

Dy adsorption on and intercalation under graphene on 6 H -SiC(0001) surface from first-principles calculations

Previous experimental observations motivate clarification of configuration stabilities and kinetic processes for intercalation of guest atoms into a layered van der Waals material such as a graphene-SiC system. From our first-principles density functional theory (DFT) calculations, we analyze Dy adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that intercalation of a single Dy atom into the gallery between TLG and BLG is more favorable than adsorption on TLG but that intercalation into the gallery underneath BLG is highly unfavorable. We obtain diffusion barriers of ~0.45 and 0.54 eV for a Dy atom diffusing on and under TLG, respectively. We find that the direct penetration of a Dy atom from the graphene top into the gallery under TLG is almost inhibited below a temperature of ~1400 K due to a large global barrier of at least ~3.5 eV. Instead, we find that a single Dy atom on TLG can easily intercalate by crossing a TLG step (e.g., a zigzag step presaturated by a Dy chain or a reconstructed zigzag step zz57). Additionally, we also perform DFT calculations for different Dy coverages to demonstrate how the favorability of Dy intercalation, as well as the corresponding interlayer spacings, depend on the coverage. Consequently, we can provide general insight and guidance for extensively studied systems involving intercalation of foreign atoms into graphene on a SiC substrate.

36 MATERIALS SCIENCE↗

Energy barriers for Dy and H penetrating graphene on 6 H -SiC(0001) and freestanding bilayer graphene from first-principles calculations

Currently, intercalation of foreign guest atoms into two-dimensional (2D) layered van der Waals materials is an active research area motivated in part by the development of next-generation energy-storage technologies and optoelectronic devices. One such extensively studied 2D material is the graphene-on-SiC system. To realize and control the desired intercalated structures, it is fundamentally important to understand the kinetic process of intercalation. For the intercalation of a guest atom into graphene layers on SiC substrate, a critical kinetic parameter is the energy barrier of a guest atom penetrating the perfect graphene top layer into the gallery under it. However, accurate theoretical calculations for such penetration barriers are unavailable in literature. From our first-principles density functional theory calculations, we obtain the global energy barriers of 3.47 and 1.80 eV for single Dy and H atoms penetrating the graphene top layer on a graphene buffer layer supported by a Si-terminated 6H-SiC(0001) substrate, respectively. Furthermore, for comparison as well as for examining the lateral strain effects, we also obtain the global barriers of 5.05 and 1.50 eV for single Dy and H atoms penetrating freestanding bilayer graphene with a tensile strain of about 8.8% to match our model for supported graphene, as well as the global barriers of 7.21 and 4.18 eV for penetrating unstrained freestanding bilayer graphene, respectively. From corresponding minimum energy paths with multiple energy minima and saddle points, we can also obtain various local energy barriers and the global backward barrier from the graphene gallery back to the top surface.

2D materials↗

Shape Stability of Truncated Octahedral fcc Metal Nanocrystals

Metallic nanocrystals (NCs) can be synthesized with tailored nonequilibrium shapes to enhance desired properties, e.g., octahedral fcc metal NCs optimize catalytic activity associated with {111} facets. However, maintenance of optimized properties requires stability against thermal reshaping. Thus, we analyze the reshaping of truncated fcc metal octahedra mediated by surface diffusion using a stochastic atomistic-level model with energetic input parameters for Pd. Furthermore, the model describes NC thermodynamics by an effective nearest-neighbor interaction and includes a realistic treatment of diffusive hopping for undercoordinated surface atoms. Kinetic Monte Carlo simulation reveals that the effective barrier, E eff , for the initial stage of reshaping is strongly tied to the degree of truncation of the vertices in the synthesized initial octahedral shapes. This feature is elucidated via exact analytic determination of the energy variation along the optimal reshaping pathway at low-temperature (T), which involves transfer of atoms from truncated {100} vertex facets to form new layers on {111} side facets. Deviations from predictions of the low-T analysis due to entropic effects are more prominent for higher T and larger NC sizes.

36 MATERIALS SCIENCE↗

Sierpinski Structure and Electronic Topology in Bi Thin Films on InSb(111)B Surfaces

Deposition of Bi on InSb(111)B reveals a striking Sierpinski-triangle (ST)-like structure in Bi thin films. Such a fractal geometric topology is further shown to turn off the intrinsic electronic topology in a thin film. Relaxation of a huge misfit strain of about 30 to 40% between Bi adlayer and substrate is revealed to drive the ST-like island formation. A Frenkel-Kontrova model is developed to illustrate the enhanced strain relief in the ST islands offsetting the additional step energy cost. Besides a sufficiently large tensile strain, forming ST-like structures also requires larger adlayer-substrate and intra-adlayer elastic stiffnesses, and weaker intra-adlayer interatomic interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane

To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. Finally, the results reported here are in reasonable agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanism of Metal Intercalation under Graphene through Small Vacancy Defects

Metal intercalation under graphene has attracted extensive experimental and theoretical research because of its capability to manipulate the electronic structure and properties of graphene. However, the pathways and mechanisms of intercalation are still not well understood. Here, we systematically investigate the intercalation process of metal atoms through graphene vacancies using first-principles calculations. We show that the energy barrier for metal atom penetration through the vacancies in graphene is small as long as the size of the vacancy is larger than a mono-vacancy. However, metal atoms are strongly bonded to the vacancy so that the detachment energy of a metal atom from the vacancy is extremely high. This inhibits the diffusion of the metal atom into the gallery beneath the surface to complete the intercalation process. On the other hand, our calculation results show that the detachment energy of a metal atom from a metal dimer at small vacancy defects is significantly reduced, making intercalation much easier. Therefore, the key step limiting the intercalation process is the detachment of the metal atoms from vacancy defects. This finding from our study provides useful insight into the defect-assisted intercalation mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Competitive formation of intercalated versus supported metal nanoclusters during deposition on layered materials with surface point defects

Intercalated metal nanoclusters (NCs) can be formed under the surface of graphite after sputtering to generate surface “portal” defects that allow deposited atoms to reach the subsurface gallery. Although, there is a competition between formation of supported NCs on top of the surface and intercalated NCs under the surface, the latter only dominating at sufficiently high temperature. A stochastic model incorporating appropriate system thermodynamics and kinetics is developed to capture this complex and competitive nucleation and growth process. Kinetic Monte Carlo simulation shows that the model captures experimental trends observed for Cu and other metals and reveals that higher temperatures are needed to facilitate detachment of atoms from supported NCs enabling them to reach the gallery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Encapsulation of metal nanoparticles at the surface of a prototypical layered material

Encapsulation of metal nanoparticles just below the surface of a prototypical layered material, graphite, is a recently discovered phenomenon. These encapsulation architectures have potential for tuning the properties of two-dimensional or layered materials, and additional applications might exploit the properties of the encapsulated metal nanoclusters themselves. The encapsulation process produces novel surface nanostructures and can be achieved for a variety of metals. Because these studies of near-surface intercalation are in their infancy, these systems provide a rich area for future studies. This Review presents the current progress on the encapsulation, including experimental strategies and characterization, as well as theoretical understanding which leads to the development of predictive capability. The Review closes with future opportunities where further understanding of the encapsulation is desired to exploit its applications.

encapsulation↗

Enhanced Nanostructure Dynamics on Au(111) with Adsorbed Sulfur due to Au–S Complex Formation

Chemisorbed species can enhance the fluxional dynamics of nanostructured metal surfaces which has implications for applications such as catalysis. Scanning tunneling microscopy studies at room temperature reveal that the presence of adsorbed sulfur (S) greatly enhances the decay rate of 2D Au islands in the vicinity of extended step edges on Au(111). This enhancement is already significant at S coverages, θS, of a few hundredths of a monolayer (ML), and is most pronounced for 0.1–0.3 ML where the decay rate is increased by a factor of around 30. For θ S close to saturation at about 0.6 ML, sulfur induces pitting and reconstruction of the entire surface, and Au islands are stabilized. Enhanced coarsening at lower θS is attributed to the formation and diffusion across terraces of Au–S complexes, particularly AuS 2 and Au4S 4 , with some lesser contribution from Au3S4. This picture is supported by density functional theory analysis of complex formation energies and diffusion barriers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic Preference for Atom Adsorption on versus Intercalation into Multilayer Graphene

The thermodynamic preference of a foreign atom for adsorption on versus intercalation into a graphitic surface is of fundamental and widespread interest. From an exhaustive first-principles density functional theory investigation for 38 typical elements over the periodic table, we reveal a quasilinear correlation between the Shannon effective ionic radius and the chemical-potential difference for a single atom from adsorption to intercalation at multilayer graphene surfaces. A critical Shannon radius is found to be around 0.10 nm, below (above) which intercalation (adsorption) is more favorable for elements with ionic-like bonding after intercalation. Single atoms with van der Waals-biased bonding show some deviation from the linear relationship, while single atoms for the elements with covalent-like bonding do not favor intercalation relative to adsorption. An energy decomposition analysis indicates that the chemical-potential difference determining the thermodynamic preference of a foreign atom for adsorption versus intercalation results from the competition between the electronic and elastic strain effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface structure of linear nanopores in amorphous silica: Comparison of properties for different pore generation algorithms

In this work, we compare the surface structure of linear nanopores in amorphous silica (a-SiO 2 ) for different versions of “pore drilling” algorithms (where the pores are generated by the removal of atoms from the preformed bulk a-SiO 2 ) and for “cylindrical resist” algorithms (where a-SiO 2 is formed around a cylindrical exclusion region). After adding H to non-bridging O, the former often results in a moderate to high density of surface silanol groups, whereas the latter produces a low density. The silanol surface density for pore drilling can be lowered by a final dehydroxylation step, and that for the cylindrical resist approach can be increased by a final hydroxylation step. In this respect, the two classes of algorithms are complementary. We focus on the characterization of the chemical structure of the pore surface, decomposing the total silanol density into components corresponding to isolated and vicinal mono silanols and geminal silanols. The final dehyroxylation and hydroxylation steps can also be tuned to better align some of these populations with the target experimental values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗