Achieving Product Control in Furfural Hydrogenation Using Intermetallic Catalysts
Explore the source record for details and available documents.
Ward, Charles J.↗
Engineering topics
Liu, Da-Jiang
Publications and source records attributed to Liu, Da-Jiang.
Molecular simulation using transfer-learned potentials for the disordered nanoscale structure of nitrogen-doped nanoporous carbons
Machine learning (ML)-based molecular dynamics (MD) simulations of the formation of a class of N-doped nanoporous carbons are performed to assess their disordered partially graphitized nanoscale structure. The study is motivated by the effectiveness of so-called nitrogen assembly carbons (NACs) for catalysis applications. Benchmark simulations for pure-C disordered graphitic systems reveal the importance of reliably capturing the vdW component of the potentials in order to accurately describe the tendency for layering of disordered graphene-like sheets. In our modeling, this is achieved by a transfer learning strategy incorporating features of the energetics from the optB88-vdW DFT functional into potentials initially trained with a less expensive functional, thereby providing a superior description of the pure-C systems. Generation from MD simulations of realistic partially graphitized structures is significantly more challenging for N-doped versus for pure C systems. However, such structures are achieved by a tailored MD simulation protocol mimicking the experimental synthesis process and in particular incorporating an annealing and subsequent quenching stages. Simulated PXRD patterns effectively reproduce the features of experimental observations for NACs, including the appearance of a prominent but broad (002) peak at around 25, and the development of another weaker feature associated with in-layer ordering of mixed C-N graphene-like sheets.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Formation and coarsening of epitaxially-supported metal nanoclusters
This mini-review describes developments over the last ~30 years in characterizing the nucleation & growth of epitaxially-supported metal nanoclusters (NCs) or islands during vapor deposition, as well as their post-deposition coarsening. A beyond-mean-field treatment for homogeneous nucleation & growth corrects the deficiencies of traditional treatments in describing, e.g., the island size distribution, but also necessitates consideration of the spatial distribution of islands and their capture zones. We discuss advances in modeling capabilities, including those based upon on an ab-initio level treatment of periphery diffusion kinetics, for description of the non-equilibrium growth shapes of these NCs, focusing on 2D NCs. For post-deposition coarsening of arrays of NCs, there is generally a competition between Ostwald Ripening (OR) and Smoluchowski Ripening (SR). SR is also known as Particle Migration & Coalescence. For 2D NCs in homoepitaxial systems, conventional OR is observed on pristine fcc(111) surfaces, dramatically enhanced OR in the presence of even trace amounts of chalcogens for Cu(111) and Ag(111), and anomalous OR on anisotropic fcc(110) surfaces. The unexpected discovery of SR for fcc(100) homoepitaxial systems prompted extensive analysis of the underlying diffusivities of 2D NCs as a function of size, as well as of NC coalescence dynamics. A comprehensive understanding of these processes is now available. Self-assembly of 3D NCs during deposition, issues related to heterogeneous nucleation, directed assembly, NC growth structure selection, and coarsening are addressed. Finally, for SR of 3D epitaxial NCs, recent insights into the size-dependence of diffusivity are described.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Enhanced Activity from Coordinatively Unsaturated and Dynamic Zeolite-Bound Organoscandium Species
Scandium borohydride grafted into the micropores of the faujasite zeolite HY30 catalyzes the C–H borylation of benzene, whereas silica-grafted species are inactive. This catalytic activity may originate from grafting at a Brønsted acid site leading to a more electron-deficient rare earth center. Herein, we apply multinuclear double-resonance nuclear magnetic resonance (NMR) experiments to probe the structure and dynamics of zeolite- and silica-bound scandium borohydride complexes. The experiments reveal that scandium centers located within the zeolite micropores, in proximity to Al-created Brønsted sites, are more dynamic than rigid scandium sites grafted on silanols. Through a combination of NMR and molecular dynamics simulations, we show that the coordination of the scandium in the zeolite is labile, with the metal exchanging between two binding sites. As a result, the weak electron donation from the support that enables the movement of the Sc center leads to the formation of an undercoordinated metal center that cannot exist on silica, ultimately leading to the new catalytic activity of the species.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Thermal Shape Stability of fcc Metal Nanocrystals Synthesized with Faceted Nonequilibrium Shapes
Highly refined capabilities of the shape-controlled solution-phase synthesis of metal nanocrystals (NCs) allow the generation of NCs with faceted nonequilibrium shapes, which optimize properties for target applications such as catalysis and plasmonics. Often, for such applications and also for TEM analysis, the NCs are removed from the solution-phase environment. We explore the postsynthesis evolution of these metastable NCs in a high-vacuum TEM environment. Specifically, here we analyze their reshaping toward the equilibrium Wulff shapes mediated by surface diffusion, where such reshaping degrades the above-mentioned optimized properties. Typical sizes for these NCs range from 5 to 30 nm or 10 3 –10 6 atoms, and reshaping often occurs on the time scale of minutes for temperatures around, say, 400 °C. We discuss the development of predictive stochastic atomistic-level models for NC evolution with a realistic description of surface diffusion. These models, in contrast to Molecular Dynamics, can naturally address the relevant time and length scales for these systems. KMC simulation results for the stochastic models are described, focusing on the reshaping of slightly elongated nanorods and of mildly truncated octahedra and nanocubes. In addition, we review appropriate theoretical formulations for reshaping, which involves the nucleation and growth on 2D islands or layers on outer facets of the NC. We note the limitations of classical nucleation theory in some scenarios and demonstrate the successes of a more fundamental and general master equation-based analysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Fluorine spillover for ceria- vs silica-supported palladium nanoparticles: A MD study using machine learning potentials
Supported metallic nanoparticles play a central role in catalysis. However, predictive modeling is particularly challenging due to the structural and dynamic complexity of the nanoparticle and its interface with the support, given that the sizes of interest are often well beyond those accessible via traditional ab initio methods. With recent advances in machine learning, it is now feasible to perform MD simulations with potentials retaining near-density-functional theory (DFT) accuracy, which can elucidate the growth and relaxation of supported metal nanoparticles, as well as reactions on those catalysts, at temperatures and time scales approaching those relevant to experiments. Furthermore, the surfaces of the support materials can also be modeled realistically through simulated annealing to include effects such as defects and amorphous structures. We study the adsorption of fluorine atoms on ceria and silica supported palladium nanoparticles using machine learning potential trained by DFT data using the DeePMD framework. We show defects on ceria and Pd/ceria interfaces are crucial for the initial adsorption of fluorine, while the interplay between Pd and ceria and the reverse oxygen migration from ceria to Pd control spillover of fluorine from Pd to ceria at later stages. In contrast, silica supports do not induce fluorine spillover from Pd particles.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations
Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.
97 MATHEMATICS AND COMPUTING↗
Nucleation-mediated reshaping of facetted metallic nanocrystals: Breakdown of the classical free energy picture
Shape stability is key to avoiding degradation of performance for metallic nanocrystals synthesized with facetted non-equilibrium shapes to optimize properties for catalysis, plasmonics, and so on. Reshaping of facetted nanocrystals is controlled by the surface diffusion-mediated nucleation and growth of new outer layers of atoms. Here, Kinetic Monte Carlo (KMC) simulation of a realistic stochastic atomistic-level model is applied to precisely track the reshaping of Pd octahedra and nanocubes. Unexpectedly, separate constrained equilibrium Monte Carlo analysis of the free energy profile during reshaping reveals a fundamental failure of the classical nucleation theory (CNT) prediction for the reshaping barrier and rate. Why? Nucleation barriers can be relatively low for these processes, so the system is not locally equilibrated before crossing the barrier, as assumed in CNT. This claim is supported by an analysis of a first-passage problem for reshaping within a master equation framework for the model that reasonably captures the behavior in KMC simulations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Efficient transfer hydrodehalogenation of halophenols catalyzed by Pd supported on ceria
We report the transfer hydrodehalogenation (THD) of halophenols is efficiently catalyzed by palladium supported on high surface area ceria (Pd/CeO 2 ) under mild conditions (65 °C) using isopropanol (iPrOH) as hydrogen source. The reactivity of 4-halophenols (4-X-PhOH) varies in the order 4-F-PhOH > 4-Cl-PhOH > 4-Br-PhOH > > 4-I-PhOH and appears to be controlled by the desorption of halides from the catalyst surface. Kinetic analysis of the reactions and temperature programmed surface reaction (TPSR) experiments indicate that oxidative addition of C-X bonds and H-abstraction from isopropoxide compete for the same active sites on Pd. The catalyst was able to conduct the THD of various hazardous pollutants and emerging contaminants (dichlorodiphenyltrichloroethane (DDT), pentachlorophenol, pentafluorophenol and triclosan).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Spatial arrangement of dynamic surface species from solid-state NMR and machine learning-accelerated MD simulations
Here, the surface arrangement of motional organic functionalities is explored by experimental dipolar coupling measurements and the prediction of motionally-averaged coupling constant from molecular dynamics simulations. The use of machine learning potentials was key to reaching the timescale required. The distance between dynamic surface species are important in cooperative heterogeneous catalysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
1 H chemical shift anisotropy: a high sensitivity solid-state NMR dynamics probe for surface studies?
We report dynamics play significant roles in chemistry and biochemistry-molecular motions impact both large- and small-scale chemical reactions in addition to biochemical processes. In many systems, including heterogeneous catalysts, the characterization of dynamics remains a challenge. The most common approaches involve the solid-state NMR measurement of anisotropic interactions, in particular 2 H quadrupolar coupling and 1 H–X dipolar coupling, which generally require isotope enrichment. Due to the high sensitivity of 1 H NMR, 1 H chemical shift anisotropy (CSA) is a particularly enticing, and underexplored, dynamics probe. We carried out 1 H CSA and 1 H– 13 C dipolar coupling measurements in a series of model supported complexes to understand how 1 H CSA can be leveraged to gain dynamic information for heterogeneous catalysts. Mathematical descriptions are given for the dynamic averaging of the CSA tensor, and its dependence on orientation and asymmetry. The variability of the orientation of the tensor in the molecular frame, in addition to its magnitude and asymmetry, negatively impacts attempts to extract quantitative dynamic information. Nevertheless, 1 H CSA measurements can reveal useful qualitative insights into the motions of a particularly dilute site, such as from a surface species.
1H chemical shift anisotropy↗
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↗
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↗
Azo(xy) vs Aniline Selectivity in Catalytic Nitroarene Reduction by Intermetallics: Experiments and Simulations
Intermetallic nanoparticles are promising catalysts in hydrogenation and fuel cell technologies. Much is known about the ability of intermetallic nanoparticles to selectively reduce nitro vs alkene, alcohol, or halide functional groups; less is known about their selectivity toward aniline vs azo or azoxy condensation products that result from the reduction of a nitro group alone. Because azo(xy)arenes bear promise as dyes, chemical stabilizers, and building blocks to functional materials but can be difficult to isolate, developing high surface area nanoparticle catalysts that display azo(xy) selectivity is desirable. To address this question, we studied a family of nanocrystalline group 10 metal (Pd, Pt)- and group 14 metal (Ge, Sn, Pb)-containing intermetallics-Pd 2 Ge, Pd 2 Sn, Pd 3 Sn 2 , Pd 3 Pb, and PtSn-in the catalytic reduction of nitroarenes. In contrast to monometallic Au, Pt, and Pd nanoparticles and "random" PdxSn1 – x nanoalloys, which are selective for aniline, nanoparticles of atomically precise intermetallic Pd 2 Ge, Pd 2 Sn, Pd 3 Sn 2 , and PtSn prefer an indirect condensation pathway and have a high selectivity for the azo(xy) products. The only exception is Pd 3 Pb, the most active among the intermetallic nanoparticles studied here, which is instead selective for aniline. Employing a novel application of molecular dynamics-based on machine learned potentials within a DeePMD framework-to heterogeneous catalysis, we are able to identify key reaction species on the different types of catalysts employed, furthering our understanding of the unique selectivity of these materials. By demonstrating how intermetallic nanoparticles can be as active yet more selective than other more traditional catalysts, this work provides new physical insights and opens new opportunities in the use of these materials in other important chemical transformations and applications.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
Observing the three-dimensional dynamics of supported metal complexes
Dynamics are intricately linked with activity and selectivity when it comes to catalysis, as noted for instance in the enzymatic principles of induced fit and allostery, and yet the range of motions heterogeneous catalytic sites are able to undergo is poorly understood. Solid-state nuclear magnetic resonance (NMR) spectroscopy is perhaps the only tool capable of probing the rapid conformational dynamics found in heterogeneous catalysts but has historically been restricted by its low sensitivity, limiting the detail with which structures can be resolved. In this work, we apply solid-state NMR and dynamic nuclear polarization, in combination with density functional theory modeling, to reveal the high-resolution structure and motional freedom of a scandium supported complex in three dimensions. The results are contrasted with the study of the analogous homoleptic complex in the crystalline state, highlighting the impacts that surface structure may have on the dynamics of supported complexes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling
Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Extended families of critical and stationary droplets for nonequilibrium phase transitions in spatially discrete bistable systems
Bistable nonequilibrium systems are realized in catalytic reaction-diffusion processes, biological transport and regulation, spatial epidemics, etc. Behavior in spatially continuous formulations, described at the mean-field level by reaction-diffusion type equations (RDEs), often mimics that of classic equilibrium van der Waals type systems. When accounting for noise, similarities include a discontinuous phase transition at some value, p eq , of a control parameter, p , with metastability and hysteresis around p eq . For each p , there is a unique critical droplet of the more stable phase embedded in the less stable or metastable phase which is stationary (neither shrinking nor growing), and with size diverging as p → p eq . Spatially discrete analogs of these mean-field formulations, described by lattice differential equations (LDEs), are more appropriate for some applications, but have received less attention. It is recognized that LDEs can exhibit richer behavior than RDEs, specifically propagation failure for planar interphases separating distinct phases. Herein, we show that this feature, together with an orientation dependence of planar interface propagation also deriving from spatial discreteness, results in the occurrence of entire families of stationary droplets. The extent of these families increases approaching the transition and can be infinite if propagation failure is realized. In addition, there can exist a regime of generic two-phase coexistence where arbitrarily large droplets of either phase always shrink. Such rich behavior is qualitatively distinct from that for classic nucleation in equilibrium and spatially continuous nonequilibrium systems.
97 MATHEMATICS AND COMPUTING↗