Engineering topics
Musgrave, Charles B.
Publications and source records attributed to Musgrave, Charles B..
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.
Electrochemical control over stoichiometry via cation intercalation into Chevrel-phase sulphides (Cu x Mo 6 S 8 , x = 1-3)
Herein this work presents a systematic investigation of the electrochemical intercalation of aqueous copper cations into the Chevrel phase (CP) Mo 6 S 8 and its effect on the host's electronic and structural characteristics as a function of stoichiometry. The electrochemical potentials at which copper intercalates into the CP Mo 6 S 8 host in 1 M CuSO 4 were identified for Cu x Mo 6 S 8 (x = 1–3) as 0.23 V, 0.11 V, and 0.055 V (vs. Cu/Cu 2+ ) for Cu 1 Mo 6 S 8 , Cu 2 Mo 6 S 8 , and Cu 3 Mo 6 S 8 , respectively, and correlated via ex situ PXRD to stoichiometric content. Potentials were also identified that yield non-stoichiometric phases with traceable Cu content-dependent structural changes. Ab initio DFT calculations were performed to probe thermodynamic stability of Cu x Mo 6 S 8 CP compositions and the metal-to-semiconductor transitions as a function of copper content. Lastly, we demonstrate that the SK-pre-edge intensity is a successful probe for intercalant content: as copper equivalents increase, so does charge transfer to the partially occupied S 3p orbitals, resulting in a decrease in the pre-edge intensity corresponding to the excitation of a S 1s-to-3p/Mo 4d hybridized orbital.
Dual atom catalysts for rapid electrochemical reduction of CO to ethylene
Strong CO adsorption and facile CO dimerization are the key challenges in electrochemical CO2 reduction towards multi-carbon (C2+) products. We recently showed that CoPc immobilized on a single-walled carbon nanotube can selectively reduce CO2 to methanol. This is enabled through molecular strain, which dramatically improves the CO adsorption energy to CoPc, which in turn facilitates methanol formation. We now examine the extended Phthalocyanine (PcEx) dual atom catalyst (DAC), which is intrinsically strained and contains two catalyst centers, making it a candidate for reducing CO to C2+ products. Using Quantum Mechanics (QM), we screened 20 elements embedded in the PcEx, seeking catalysts with weak hydrogen binding, strong CO binding, and facile CO dimerization. We identi>ied Fe, Ru, Co, and Ir as the best performers and subsequently evaluated the entire CO to C2H4 mechanism (9 steps) using each of these elements as catalysts. In terms of limiting potential and overall exergonicity, we identi>ied CoPcEx as the best catalyst, followed by IrPcEx. We then examined the full CO to C2H4 mechanism on the bimetallic IrCoPcEx catalyst using grand canonical QM to obtain the reaction energetics as a function of applied potential. We conclude that the bimetallic IrCoPcEx is most promising for ef>iciently converting CO to ethylene.
Tailoring the Near-Surface Environment of Rh Single-Atom Catalysts for Selective CO 2 Hydrogenation
For this study, we used a combination of experimental spectroscopies, density functional theory calculations, and CO 2 hydrogenation studies to investigate the effects of modifying single-atom Rh 1 /TiO 2 catalysts with functionalized phosphonic acid monolayers. We found that the deposition of specific amine-functionalized ligands resulted in an ~8× increase in site-specific CO 2 reduction turnover frequency at 150 °C and a ~2× increase at 250 °C. On-stream stability also improved following ligand deposition. The effect of the modifier on reactivity was highly sensitive to the proximity of the amine functional group to the surface, which was controlled by adjusting the length of the phosphonic acid tail. Furthermore, deposition of alkyl phosphonic acids without an amine functional group resulted in blocked CO 2 adsorption and a near-complete loss of catalytic activity. Infrared spectroscopy studies suggested that the amine group provided binding sites for CO 2 that enabled hydrogenation when the amine was positioned near a Rh1 site. Phosphonic acid-modified catalysts also exhibited high selectivity to CO over the series product methane; the selectivity effect was traced to modification of the Rh1 sites to favor CO desorption. Phosphonic acid deposition resulted in 80–90% loss of accessible Rh1 sites, likely due to blocking by tail groups. However, even with the loss of sites, under low-temperature reaction conditions, the rates of CO 2 hydrogenation were improved with the coatings, indicating that the remaining sites are highly efficient. Organic functionalization of the supports for atomically dispersed catalysts offers the opportunity to precisely control the positioning of functional groups in the vicinity of a well-defined active site, potentially enabling an additional level of control over active site design.
Revised Nitrogen Reduction Scaling Relations from Potential-Dependent Modeling of Chemical and Electrochemical Steps
The electrochemical nitrogen reduction reaction (NRR) is a promising route to enable carbon-free ammonia production. However, this reaction is limited by the poor activity and selectivity of current catalysts. The rational design of superior NRR electrocatalysts requires a detailed mechanistic understanding of current material limitations to inform how these might be overcome. The current understanding of how scaling limits NRR on metal catalysts is predicated on a simplified reaction pathway that considers only proton-coupled electron transfer (PCET) steps. Here, we apply grand-canonical density functional theory to investigate a more comprehensive NRR mechanism that includes both electrochemical and chemical steps on 30 metal surfaces in solvent under an applied potential. We applied ..phi..max, a grand-canonical adaptation of the Gmax thermodynamic descriptor, to evaluate trends in catalyst activity. This approach produces a ..phi..max "volcano" diagram for NRR activity scaling on metals that qualitatively differs from the scaling relations identified when only PCET steps are considered. NH3* desorption was found to limit the NRR activity for materials at the top of the volcano and truncate the volcano's peak at increasingly reducing potentials. These revised scaling relations may inform the rational design of superior NRR electrocatalysts. This approach is transferable to study other materials and reaction chemistries where both electrochemical and chemical steps are modeled under an applied potential.
Strain enhances the activity of molecular electrocatalysts via carbon nanotube supports
Abstract Support-induced strain engineering is useful for modulating the properties of two-dimensional materials. However, controlling strain of planar molecules is technically challenging due to their sub-2 nm lateral size. Additionally, the effect of strain on molecular properties remains poorly understood. Here we show that carbon nanotubes (CNTs) are ideal substrates for inducing optimum properties through molecular curvature. In a tandem-flow electrolyser with monodispersed cobalt phthalocyanine (CoPc) on single-walled CNTs (CoPc/SWCNTs) for CO 2 reduction, we achieve a methanol partial current density of >90 mA cm −2 with >60% selectivity, surpassing wide multiwalled CNTs at 16.6%. We report vibronic and X-ray spectroscopies to unravel the distinct local geometries and electronic structures induced by the strong molecule–support interactions. Grand canonical density functional theory confirms that curved CoPc/SWCNTs improve *CO binding to enable subsequent reduction, whereas wide multiwalled CNTs favour CO desorption. Our results show the important role of SWCNTs beyond catalyst dispersion and electron conduction.
Pd(II) and Rh(I) Catalytic Precursors for Arene Alkenylation: Comparative Evaluation of Reactivity and Mechanism Based on Experimental and Computational Studies
We combine experimental and computational investigations to compare and understand catalytic arene alkenylation using the Pd(II) and Rh(I) precursors Pd(OAc) 2 and [(η 2 -C 2 H 4 ) 2 Rh(µ-OAc)] 2 with arene, olefin and Cu(II) carboxylate at elevated temperature (> 120 °C). Under specific conditions, previous computational and experimental efforts have identified heterotrimetallic cyclic PdCu 2 (η 2 -C 2 H 4 ) 3 (μ-OPiv) 6 and [(η 2 -C 2 H 4 ) 2 Rh(µ-OPiv) 2 ] 2 (µ-Cu) (OPiv = pivalate) species as likely active catalysts for these processes. Further studies of catalyst speciation suggest a complicated equilibrium between Cu(II)-containing complexes containing one Rh or Pd atom with complexes containing two Rh or Pd atoms. At 120 °C, Rh catalysis produces styrene > 20-fold more rapidly than Pd. Also, at 120 °C, Rh is ~98% selective for styrene formation while Pd is ~82% selective. Furthermore, our studies indicate that Pd catalysis has a higher predilection toward olefin functionalization to form undesired vinyl ester, while Rh catalysis is more selective for arene/olefin coupling. However, at elevated temperatures, Pd converts vinyl ester and arene to vinyl arene, which is proposed to occur through low valent Pd(0) clusters that are formed in situ. Regardless of arene functionality, the regioselectivity for alkenylation of mono-substituted arenes with the Rh catalyst gives an approximate 2:1 meta:para ratio with minimal ortho C–H activation. In contrast, Pd selectivity is significantly influenced by arene electronics with electron-rich arenes giving an approximate 1:2:2 ortho:meta:para ratio while the electron deficient (α,α,α)-trifluorotoluene gives a 3:1 meta:para ratio with minimal ortho functionalization. Kinetic intermolecular arene ethenylation competition experiments find that Rh reacts most rapidly with benzene, and the rate of mono-substituted arene alkenylation does not correlate with arene electronics. In contrast, with Pd catalysis, electron-rich arenes react more rapidly than benzene while electron-deficient arenes react less rapidly than benzene. These experimental findings, in combination with computational results, are consistent with the arene C–H activation step for Pd catalysis involving significant η 1 -arenium character due to Pd-mediated electrophilic aromatic substitution character. In contrast, the mechanism for Rh catalysis is not sensitive to arene substituent electronics, which we propose indicates less electrophilic aromatic substitution character for the Rh-mediated arene C–H activation.
Rhodium-Catalyzed Alkenylation of Arenes with Multi-Substituted Olefins: Comparison of Selectivity and Reaction Rate as a Function of Olefin Identity
Rhodium-catalyzed arene alkenylation using Cu(II) carboxylates as the in situ oxidant and mono-substituted olefins has been previously reported (e.g., J. Am. Chem. Soc. 2019, 139, 5474; J. Am. Chem. Soc. 2018, 140, 17007; Organometallics 2019, 38, 3860; J. Am. Chem. Soc. 2020, 142, 10534). Herein, studies are extended to multi-substituted olefins with the goal of evaluating the effect of olefin substitution pattern and substituent identity on selectivity and turnover frequency. The influence of olefin substitution is probed by comparing the conversion of benzene to alkenyl arenes with ethylene, propylene, 1-butene, cis-2-butene, trans-2-butene, isobutene, 2-methyl-2-butene, and tetramethylethylene as well as the phenyl-substituted olefins and isomers of propenylbenzene. The rate of oxidative hydrophenylation for multi-substituted olefins follows the trend monosubstituted > disubstituted > trisubstituted, and tetrasubstituted olefins are unreactive. To probe the effect of substituent size on Markovnikov/ anti-Markovnikov regioselectivity, cyclohexyl, tert-butyl, isopropyl, ethyl, and methyl substituted α-olefins are compared. Selectivity for anti-Markovnikov products generally increases as substituent steric bulk is increased. Tolerance for some functionalized olefins is demonstrated. Here, the ortho/meta/para regioselectivity with mono-substituted arenes reveals that arene and olefin identity influences selectivity. Further mechanistic studies provide evidence for Curtin–Hammet control of ortho/meta/para regioselectivity with monosubstituted arenes.
Dataset of theoretical multinary perovskite oxides
Perovskite oxides (ternary chemical formula ABO 3 ) are a diverse class of materials with applications including heterogeneous catalysis, solid-oxide fuel cells, thermochemical conversion, and oxygen transport membranes. However, their multicomponent (chemical formula $A_xA^{'}_{1-x}B_yB^{'}_{1-y}O_3$) chemical space is underexplored due to the immense number of possible compositions. To expand the number of computed $A_xA^{'}_{1-x}B_yB^{'}_{1-y}O_3$ compounds we report a dataset of 66,516 theoretical multinary oxides, 59,708 of which are perovskites. First, 69,407 $A_{0.5}A^{'}_{0.5}B_{0.5}B^{'}_{0.5}O_3$ compositions were generated in the a - b + a - Glazer tilting mode using the computationally-inexpensive Structure Prediction and Diagnostic Software (SPuDS) program. Next, we optimized these structures with density functional theory (DFT) using parameters compatible with the Materials Project (MP) database. Our dataset contains these optimized structures and their formation (ΔH f ) and decomposition enthalpies (ΔH d ) computed relative to MP tabulated elemental references and competing phases, respectively. This dataset can be mined, used to train machine learning models, and rapidly and systematically expanded by optimizing more SPuDS-generated $A_{0.5}A^{'}_{0.5}B_{0.5}B^{'}_{0.5}O_3$ perovskite structures using MP-compatible DFT calculations.
In-operando synchrotron experiments of flash sintering carried out in current rate mode
Here, we report first-time results for in-operando flash sintering synchrotron experiments carried out in current rate mode where the specimen, held at a constant temperature, is fed current that is increased at a constant rate. These experiments are unique because the time dependence of the sintering behavior can be stretched out over a longer period (by changing the current rate) than in voltage-to-current experiments in which sintering occurs in a burst at the onset of the flash. Two results are presented: (i) A comparison of temperatures measured with the platinum standard to those predicted by the black body radiation model leading to estimates of the emissivity as a function of porosity whereby emissivity increases from 0.65 to 0.9 as the sample sinters from its green state to full density, and (ii) measurements of the excess lattice expansion as a function of density as the sample sinters continuously while the current is increased. The present work highlights the promise of current rate experiments to obtain results while the sample sinters gradually from its green density to full density (somewhat akin to conventional sintering) for gaining further insights into the mechanisms of flash sintering.
Chapter 2: Evaluating a Concentrating Solar Power Plantas an Extended-Duration Peaking Resource
We explore the ability of a concentrating solar power (CSP) plant with thermal energy storage (TES) to provide peaking capacity. We focus on future power systems, wherein net load patterns may be significantly different than they are today (e.g., due to higher renewable-energy penetrations). We examine 28 locations in the southwestern United States over an 18-year period. The hourly operation of the CSP plants are simulated to determine their potential to provide energy during an eight-hour peak-load window for each day up to 365 days per year. Our result shows that for the large majority of locations and years, CSP plants with certain configurations (i.e., in terms of solar field and TES sizes) can provide nearly 100% peak-load capacity. We examine also the amount of supplemental energy (e.g.,by using natural gas as a supplemental thermal-energy source) that would be required to ensure that a CSP plant could serve the eight highest-load hours of every day of the year. We find that in most cases, a CSP plant supplemented with natural gas would require less than 5% of the fuel that is used by a natural-gas fired power plant providing the same level of reliable capacity. A series of sensitivity analyses show that these results are robust to the number of peak-load hours and days that are considered and the configuration of the CSP plant.
Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution
Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.
Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory
Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.
How the Bioinspired Fe 2 Mo 6 S 8 Chevrel Breaks Electrocatalytic Nitrogen Reduction Scaling Relations
The nitrogen reduction reaction (NRR) is a renewable alternative to the energy- and CO 2 -intensive Haber–Bosch NH 3 synthesis process but is severely limited by the low activity and selectivity of studied electrocatalysts. The Chevrel phase Fe 2 Mo 6 S 8 has a surface Fe–S–Mo coordination environment that mimics the nitrogenase FeMo-cofactor and was recently shown to provide state-of-the-art activity and selectivity for NRR. Here, we elucidate the previously unknown NRR mechanism on Fe 2 Mo 6 S 8 via grand-canonical density functional theory (GC-DFT) that realistically models solvated and biased surfaces. Fe sites of Fe 2 Mo 6 S 8 selectively stabilize the key *NNH intermediate via a narrow band of free-atom-like surface d-states that selectively hybridize with p-states of *NNH, which results in Fe sites breaking NRR scaling relationships. These sharp d-states arise from an Fe–S bond dissociation during N 2 adsorption that mimics the mechanism of the nitrogenase FeMo-cofactor. Furthermore, we developed a new GC-DFT-based approach for calculating transition states as a function of bias (GC-NEB) and applied it to produce a microkinetic model for NRR at Fe 2 Mo 6 S 8 that predicts high activity and selectivity, in close agreement with experiments. Furthermore, our results suggest new design principles that may identify effective NRR electrocatalysts that minimize the barriers for *N 2 protonation and *NH 3 desorption and that may be broadly applied to the rational discovery of stable, multinary electrocatalysts for other reactions where narrow bands of surface d-states can be tuned to selectively stabilize key reaction intermediates and guide selectivity toward a target product. Furthermore, our results highlight the importance of using GC-DFT and GC-NEB to accurately model electrocatalytic reactions.